Use of locally weighted scatterplot smoothing (LOWESS) regression to study selection signatures in Piedmontese and Italian Brown cattle breeds

Selection is the major force affecting local levels of genetic variation in species. The availability of dense marker maps offers new opportunities for a detailed understanding of genetic diversity distribution across the animal genome. Over the last 50 years, cattle breeds have been subjected to intense artificial selection. Consequently, regions controlling traits of economic importance are expected to exhibit selection signatures. The fixation index (F_st) is an estimate of population differentiation, based on genetic polymorphism data, and it is calculated using the relationship between inbreeding and heterozygosity. In the present study, locally weighted scatterplot smoothing (LOWESS) regression and a control chart approach were used to investigate selection signatures in two cattle breeds with different production aptitudes (dairy and beef). F_st was calculated for 42 514 SNP marker loci distributed across the genome in 749 Italian Brown and 364 Piedmontese bulls. The statistical significance of F_st values was assessed using a control chart. The LOWESS technique was efficient in removing noise from the raw data and was able to highlight selection signatures in chromosomes known to harbour genes affecting dairy and beef traits. Examples include the peaks detected for BTA2 in the region where the myostatin gene is located and for BTA6 in the region harbouring the ABCG2 locus. Moreover, several loci not previously reported in cattle studies were detected.     

Predicting haplotype carriers from SNP genotypes in Bos taurus through linear discriminant analysis

Background

SNP (single nucleotide polymorphisms) genotype data are increasingly available in cattle populations and, among other things, can be used to predict carriers of specific haplotypes. It is therefore convenient to have a practical statistical method for the accurate classification of individuals into carriers and non-carriers. In this paper, we present a procedure combining variable selection (i.e. the selection of predictive SNPs) and linear discriminant analysis for the identification of carriers of a haplotype on BTA19 (Bos taurus autosome 19) known to be associated with reduced cow fertility. A population of 3645 Brown Swiss cows and bulls genotyped with the 54K SNP-chip was available for the analysis.

Results

The overall error rate for the prediction of haplotype carriers was on average very low (∼≤1%). The error rate was found to depend on the number of SNPs in the model and their density around the region of the haplotype on BTA19. The minimum set of SNPs to still achieve accurate predictions was 5, with a total test error rate of 1.59.

Conclusions

The paper describes a procedure to accurately identify haplotype carriers from SNP genotypes in cattle populations. Very few misclassifications were observed, which indicates that this is a very reliable approach for potential applications in cattle breeding.     

Use of partial least squares regression to predict single nucleotide polymorphism marker genotypes when some animals are genotyped with a low-density panel

High-density single nucleotide polymorphism (SNP) platforms are currently used in genomic selection (GS) programs to enhance the selection response. However, the genotyping of a large number of animals with high-throughput platforms is rather expensive and may represent a constraint for a large-scale implementation of GS. The use of low-density marker (LDM) platforms could overcome this problem, but different SNP chips may be required for each trait and/or breed. In this study, a strategy of imputation independent from trait and breed is proposed. A simulated population of 5865 individuals with a genome of 6000 SNP equally distributed on six chromosomes was considered. First, reference and prediction populations were generated by mimicking high- and low-density SNP platforms, respectively. Then, the partial least squares regression (PLSR) technique was applied to reconstruct the missing SNP in the low-density chip. The proportion of SNP correctly reconstructed by the PLSR method ranged from 0.78 to 0.97 when 90% and 50%, respectively, of genotypes were predicted. Moreover, data sets consisting of a mixture of actual and PLSR-predicted SNP or only actual SNP were used to predict genomic breeding values (GEBVs). Correlations between GEBV and true breeding values varied from 0.74 to 0.76, respectively. The results of the study indicate that the PLSR technique can be considered a reliable computational strategy for predicting SNP genotypes in an LDM platform with reasonable accuracy.     

Effects of double ligation of Stensen's duct on the rabbit parotid gland

Salivary gland duct ligation is an alternative to gland excision for treating sialorrhea or reducing salivary gland size prior to tumor excision. Duct ligation also is used as an approach to study salivary gland aging, regeneration, radiotherapy, sialolithiasis and sialadenitis. Reports conflict about the contribution of each salivary cell population to gland size reduction after ductal ligation. Certain cell populations, especially acini, reportedly undergo atrophy, apoptosis and proliferation during reduction of gland size. Acini also have been reported to de-differentiate into ducts. These contradictory results have been attributed to different animal or salivary gland models, or to methods of ligation. We report here a bilateral double ligature technique for rabbit parotid glands with histologic observations at 1, 7, 14, 30, 60 days after ligation. A large battery of special stains and immunohistochemical procedures was employed to define the cell populations. Four stages with overlapping features were observed that led to progressive shutdown of gland activities: 1) marked atrophy of the acinar cells occurred by 14 days, 2) response to and removal of the secretory material trapped in the acinar and ductal lumens mainly between 30 and 60 days, 3) reduction in the number of parenchymal (mostly acinar) cells by apoptosis that occurred mainly between 14–30 days, and 4) maintenance of steady-state at 60 days with a low rate of fluid, protein, and glycoprotein secretion, which greatly decreased the number of leukocytes engaged in the removal of the luminal contents. The main post- ligation characteristics were dilation of ductal and acinar lumens, massive transient infiltration of mostly heterophils (rabbit polymorphonuclear leukocytes), acinar atrophy, and apoptosis of both acinar and ductal cells. Proliferation was uncommon except in the larger ducts. By 30 days, the distribution of myoepithelial cells had spread from exclusively investing the intercalated ducts pre-ligation to surrounding a majority of the residual duct-like structures, many of which clearly were atrophic acini. Thus, both atrophy and apoptosis made major contributions to the post-ligation reduction in gland size. Structures also occurred with both ductal and acinar markers that suggested acini differentiating into ducts. Overall, the reaction to duct ligation proceeded at a considerably slower pace in the rabbit parotid glands than has been reported for the salivary glands of the rat.     

Sequence of a cDNA for mouse thymidylate synthase reveals striking similarity with the prokaryotic enzyme

We report the nucleotide sequence of a cloned cDNA, pMTS-3, that contains a 1-kb insert corresponding to mouse thymidylate synthase (E.C. 2.1.1.45). The open reading frame of 921 nucleotides from the first AUG to the termination codon specifies a protein with a molecular mass of 34,962 daltons. The predicted amino acid sequence is 90% identical with that of the human enzyme. The mouse sequence also has an extremely high degree of similarity (as much as 55% identity) with prokaryotic thymidylate synthase sequences, indicating that thymidylate synthase is among the most highly conserved proteins studied to date. The similarity is especially pronounced (as much as 80% identity) in the 44-amino-acid region encompassing the binding site for deoxyuridylic acid. The cDNA sequence also suggests that mouse thymidylate synthase mRNA lacks a 3' untranslated region, since the termination codon, UAA, is followed immediately by a poly(A) segment.      

Inhibition of a nutrient-dependent pinocytosis in dictyostelium discoideum by the amino acid analogue hadacidin

In the present study we examine the effects of the drug hadacidin (N-formyl-N- hydroxyglycine) on pinocytosis in the eukaryotic microorganism dictyostelium discoideum. At concentrations of up to approximately 8 mg/ml, hadacidin inhibited the rate of pinocytosis of fluorescein isothiocyanate (FITC) dextran in cells in growth medium in a concentration-dependent manner but had no effect on cells in starvation medium. Because hadacidin also inhibits cellular proliferation at this concentration, the relationship between growth rate and pinocytosis was studied further using another drug, cerulenin, to produce growth-arrest. These experiments showed no changes in the rate pinocytosis even after complete cessation of cellular proliferation. Other studies showed that the transfer of cells from growth to starvation medium reduced the rate of pinocytosis by approximately 50 percent. A reduction of similar magnitude occurred if cells were transferred from growth to starvation medium containing hadacidin. Also, no additional reduction in pinocytosis occurred when cells that had been treated with hadacidin were transferred to starvation medium containing hadacidin. These cells were able to take up [(14)C]hadacidin in the starvation medium. In contrast to the results with hadacidin-treated cells, cells in a cerulenin-induced state of growth-arrest when transferred to starvation medium exhibited the same 50 percent reduction in pinocytosis observed in cells not previously exposed to either drug. Cells treated with azide, in either growth or starvation medium, exhibited an immediate inhibition of all pinocytotic activity. After the transfer of log-phase cells to starvation medium supplemented with glucose, the reduction in rate was only approximately 10-15 percent. In contrast, a 50 percent reduction was observed after supplementation of starvation medium with sucrose, KCl, or concanavalin A. Maintaining the cells in growth medium containing hadacidin for as long as 16 h had no effect on the rate at which cells aggregated. These results are consistent with the conclusion that D. discoideum exhibits two types of pinocytotic activity: one that is nutrient dependent and the other independent of nutrients. This latter activity persists in starvation medium and is unaffected by hadacidin, whereas the nutrient-dependent activity is present in growth medium and is inhibited by hadacidin.     

Exploring the relationship between bacterial genera and lipid metabolism in bovine rumen

The rumen is characterised by a complex microbial ecosystem, which is particularly active in lipid metabolism. Several studies demonstrated a role of diet and breed on bacterial community profile, with the effect on metabolic pathways. Despite the knowledge achieved on metabolism and the bacterial profile, little is known about the relationship between individual bacteria and metabolic pathways. Therefore, a multivariate approach was used to search for possible relationships between bacteria and products of several pathways. The correlation between rumen bacterial community composition and rumen lipid metabolism was assessed in 40 beef steers (20 Maremmana and 20 Aubrac) reared with the same system and fed the same diet. A canonical discriminant analysis combined with a canonical correlation analysis (CCA) was performed to explore this correlation. The variables showing a Pearson correlation higher than 0.6 as absolute value and significant were retained for CCA considering the relationship of bacterial composition with several metabolic pathways. The results indicated that some bacterial genera could have significant impacts on the presence of several fatty acids. However, the relationship between genera and fatty acid changes according to the breed, demonstrating that the metabolic pathways change according to the host genetic background, related to breed evolution, although there is also an intra-breed genetic background which should not be ignored. In Maremmana, Succiniclasticum and Rikenellaceae_RC9_gut_group showed a high positive correlation with dimethylacetals (DMAs) DMA_C13:0, DMA_C14:0, DMA_C14:0iso, DMA_C15:0, DMA_C15:0iso, and DMA_C18:0. Prevotellaceae_UCG-003 correlates with C18:3c9c12c15 and C18:1t11, while Fibrobacter and Succiniclasticum correlate with C18:2c9t11 and Lachnospiraceae_NK3A20_group correlates with C18:1c12. Prevotellaceae_UCG-003, Ruminococcaceae_UCG-010, and Oribacterium showed a positive correlation with C13:0iso, and C17:0. Conversely, in Aubrac, Treponema_2 and Rikenellaceae_RC9_gut_group correlated with DMA_C14:0iso, DMA_C16:0iso, DMA_C17:0iso, while Ruminococcaceae_UCG-010, Christensenellaceae_R-7_group and Ruminococcaceae_NK4A214_group correlated with DMA_C18:1t11, DMA_C14:0, DMA_C18:1c12. Acetitomaculum correlated with C18:2c9c12, C18:1c12, C18:1c13, C18:1t12 and Lachnospiraceae_NK3A20_group with C18:1t6-8 and C18:1t9. Saccharofermentas, Ruminococcaceae_UCG-010 and Rikenellaceae_RC9_gut_group correlated with C18:2c9t11 while, Prevotellaceae_UCG-001 and Ruminococcus_1 correlated with C14:0iso, C15:0, C15:0iso, C17:0. Saccharofermentans, Rikenellaceae_RC9_gut_group, Ruminococcaceae_NK4A214_group, and Ruminococcaceae_UCG-010 correlated with C13:1c12 and C16:0iso. These results lead to hypothesise a possible association between several metabolic pathways and one or a few bacterial genera. If these associations are confirmed by further investigations that verify the causality of a bacterial genus with a particular metabolic process, it will be possible to deepen the knowledge on the activity of the rumen population in lipid metabolism. This approach appears to be a promising tool for uncovering the correlation between bacterial genera and products of rumen lipid metabolism.     

The discovery of aminopyrazines as novel, potent Na(v)1.7 antagonists: hit-to-lead identification and SAR

Herein the discovery of a novel class of aminoheterocyclic Na(v)1.7 antagonists is reported. Hit compound 1 was potent but suffered from poor pharmacokinetics and selectivity. The compact structure of 1 offered a modular synthetic strategy towards a broad structure-activity relationship analysis. This analysis led to the identification of aminopyrazine 41, which had vastly improved hERG selectivity and pharmacokinetic properties.     

Factors that drive the increasing use of FFPE tissue in basic and translational cancer research

The decision to use 10% neutral buffered formalin fixed, paraffin embedded (FFPE) archival pathology material may be dictated by the cancer research question or analytical technique, or may be governed by national ethical, legal and social implications (ELSI), biobank, and sample availability and access policy. Biobanked samples of common tumors are likely to be available, but not all samples will be annotated with treatment and outcomes data and this may limit their application. Tumors that are rare or very small exist mostly in FFPE pathology archives. Pathology departments worldwide contain millions of FFPE archival samples, but there are challenges to availability. Pathology departments lack resources for retrieving materials for research or for having pathologists select precise areas in paraffin blocks, a critical quality control step. When samples must be sourced from several pathology departments, different fixation and tissue processing approaches create variability in quality. Researchers must decide what sample quality and quality tolerance fit their specific purpose and whether sample enrichment is required. Recent publications report variable success with techniques modified to examine all common species of molecular targets in FFPE samples. Rigorous quality management may be particularly important in sample preparation for next generation sequencing and for optimizing the quality of extracted proteins for proteomics studies. Unpredictable failures, including unpublished ones, likely are related to pre-analytical factors, unstable molecular targets, biological and clinical sampling factors associated with specific tissue types or suboptimal quality management of pathology archives. Reproducible results depend on adherence to pre-analytical phase standards for molecular in vitro diagnostic analyses for DNA, RNA and in particular, extracted proteins. With continuing adaptations of techniques for application to FFPE, the potential to acquire much larger numbers of FFPE samples and the greater convenience of using FFPE in assays for precision medicine, the choice of material in the future will become increasingly biased toward FFPE samples from pathology archives. Recognition that FFPE samples may harbor greater variation in quality than frozen samples for several reasons, including variations in fixation and tissue processing, requires that FFPE results be validated provided a cohort of frozen tissue samples is available.