Identification of a two-marker-haplotype on Bos taurus autosome 18 associated with somatic cell score in German Holstein cattle

<p>Abstract</p> <p>Background</p> <p>The somatic cell score (SCS) is implemented in routine sire evaluations in many countries as an indicator trait for udder health. Somatic cell score is highly correlated with clinical mastitis, and in the German Holstein population quantitative trait loci (QTL) for SCS have been repeatedly mapped on <it>Bos taurus </it>autosome 18 (BTA18). In the present study, we report a refined analysis of previously detected QTL regions on BTA18 with the aim of identifying marker and marker haplotypes in linkage disequilibrium with SCS. A combined linkage and linkage disequilibrium approach was implemented, and association analyses of marker genotypes and maternally inherited two-marker-haplotypes were conducted to identify marker and haplotypes in linkage disequilibrium with a locus affecting SCS in the German Holstein population.</p> <p>Results</p> <p>We detected a genome-wide significant QTL within marker interval 9 (<it>HAMP_c.366+109G>A </it>- <it>BMS833</it>) in the middle to telomeric region on BTA18 and a second putative QTL in marker interval 12-13 (<it>BB710 </it>- <it>PVRL2_c.392G>A</it>). Association analyses with genotypes of markers flanking the most likely QTL positions revealed the microsatellite marker <it>BMS833 </it>(interval 9) to be associated with a locus affecting SCS within the families investigated. A further analysis of maternally inherited two-marker haplotypes and effects of maternally inherited two-marker-interval gametes indicated haplotype <it>249-G </it>in marker interval 12-13 (<it>BB710 </it>- <it>PVRL2_c.392G>A</it>) to be associated with SCS in the German Holstein population.</p> <p>Conclusion</p> <p>Our results confirmed previous QTL mapping results for SCS and support the hypothesis that more than one locus presumably affects udder health in the middle to telomeric region of BTA18. However, a subsequent investigation of the reported QTL regions is necessary to verify the two-QTL hypothesis and confirm the association of two-marker-haplotype <it>249-G </it>in marker interval 12-13 (<it>BB710 </it>- <it>PVRL2_c.392G>A</it>) with SCS. For this purpose, higher marker density and multiple-trait and multiple-QTL models are required to narrow down the position of the causal mutation or mutations affecting SCS in German Holstein cattle.</p

Cryptosporidium parvum infection alters the intestinal mucosa transcriptome in neonatal calves: implications for immune function

One of the leading causes of infectious diarrhea in newborn calves is the apicomplexan protozoan Cryptosporidium parvum (C. parvum). However, little is known about its immunopathogenesis. Using next generation sequencing, this study investigated the immune transcriptional response to C. parvum infection in neonatal calves. Neonatal male Holstein-Friesian calves were either orally infected (N = 5) or not (CTRL group, N = 5) with C. parvum oocysts (gp60 subtype IIaA15G2R1) at day 1 of life and slaughtered on day 7 after infection. Total RNA was extracted from the jejunal mucosa for short read. Differentially expressed genes (DEGs) between infected and CTRL groups were assessed using DESeq2 at a false discovery rate &lt; 0.05. Infection did not affect plasma immunohematological parameters, including neutrophil, lymphocyte, monocyte, leucocyte, thrombocyte, and erythrocyte counts as well as hematocrit and hemoglobin concentration on day 7 post infection. The immune-related DEGs were selected according to the UniProt immune system process database and were used for gene ontology (GO) and pathway enrichment analysis using Cytoscape (v3.9.1). Based on GO analysis, DEGs annotated to mucosal immunity, recognizing and presenting antigens, chemotaxis of neutrophils, eosinophils, natural killer cells, B and T cells mediated by signaling pathways including toll like receptors, interleukins, tumor necrosis factor, T cell receptor, and NF-KB were upregulated, while markers of macrophages chemotaxis and cytosolic pattern recognition were downregulated. This study provides a holistic snapshot of immune-related pathways induced by C. parvum in calves, including novel and detailed feedback and feedforward regulatory mechanisms establishing the crosstalk between innate and adaptive immune response in neonate calves, which could be utilized further to develop new therapeutic strategies

Comparative expression profiling of E. coli and S. aureus inoculated primary mammary gland cells sampled from cows with different genetic predispositions for somatic cell score

BACKGROUND: During the past ten years many quantitative trait loci (QTL) affecting mastitis incidence and mastitis related traits like somatic cell score (SCS) were identified in cattle. However, little is known about the molecular architecture of QTL affecting mastitis susceptibility and the underlying physiological mechanisms and genes causing mastitis susceptibility. Here, a genome-wide expression analysis was conducted to analyze molecular mechanisms of mastitis susceptibility that are affected by a specific QTL for SCS on Bos taurus autosome 18 (BTA18). Thereby, some first insights were sought into the genetically determined mechanisms of mammary gland epithelial cells influencing the course of infection. METHODS: Primary bovine mammary gland epithelial cells (pbMEC) were sampled from the udder parenchyma of cows selected for high and low mastitis susceptibility by applying a marker-assisted selection strategy considering QTL and molecular marker information of a confirmed QTL for SCS in the telomeric region of BTA18. The cells were cultured and subsequently inoculated with heat-inactivated mastitis pathogens Escherichia coli and Staphylococcus aureus, respectively. After 1, 6 and 24 h, the cells were harvested and analyzed using the microarray expression chip technology to identify differences in mRNA expression profiles attributed to genetic predisposition, inoculation and cell culture. RESULTS: Comparative analysis of co-expression profiles clearly showed a faster and stronger response after pathogen challenge in pbMEC from less susceptible animals that inherited the favorable QTL allele 'Q' than in pbMEC from more susceptible animals that inherited the unfavorable QTL allele 'q'. Furthermore, the results highlighted RELB as a functional and positional candidate gene and related non-canonical Nf-kappaB signaling as a functional mechanism affected by the QTL. However, in both groups, inoculation resulted in up-regulation of genes associated with the Ingenuity pathways 'dendritic cell maturation' and 'acute phase response signaling', whereas cell culture affected biological processes involved in 'cellular development'. CONCLUSIONS: The results indicate that the complex expression profiling of pathogen challenged pbMEC sampled from cows inheriting alternative QTL alleles is suitable to study genetically determined molecular mechanisms of mastitis susceptibility in mammary epithelial cells in vitro and to highlight the most likely functional pathways and candidate genes underlying the QTL effect

From FAANG to Fork: Application of Highly Annotated Genomes to Improve Farmed Animal Production

Here we review and describe a set of research priorities to meet present and future challenges posed to farmed animal production that build on progress, successes and resources from the Functional Annotation of ANimal Genomes (FAANG) project

Combined analysis of data from two granddaughter designs: A simple strategy for QTL confirmation and increasing experimental power in dairy cattle

A joint analysis of five paternal half-sib Holstein families that were part of two different granddaughter designs (ADR- or Inra-design) was carried out for five milk production traits and somatic cell score in order to conduct a QTL confirmation study and to increase the experimental power. Data were exchanged in a coded and standardised form. The combined data set (JOINT-design) consisted of on average 231 sires per grandsire. Genetic maps were calculated for 133 markers distributed over nine chromosomes. QTL analyses were performed separately for each design and each trait. The results revealed QTL for milk production on chromosome 14, for milk yield on chromosome 5, and for fat content on chromosome 19 in both the ADR- and the Inra-design (confirmed within this study). Some QTL could only be mapped in either the ADR- or in the Inra-design (not confirmed within this study). Additional QTL previously undetected in the single designs were mapped in the JOINT-design for fat yield (chromosome 19 and 26), protein yield (chromosome 26), protein content (chromosome 5), and somatic cell score (chromosome 2 and 19) with genomewide significance. This study demonstrated the potential benefits of a combined analysis of data from different granddaughter designs

Dominance and parent-of-origin effects of coding and non-coding alleles at the acylCoA-diacylglycerol-acyltransferase (DGAT1) gene on milk production traits in German Holstein cows

<p>Abstract</p> <p>Background</p> <p>Substantial gene substitution effects on milk production traits have formerly been reported for alleles at the K232A and the promoter VNTR loci in the bovine acylCoA-diacylglycerol-acyltransferase 1 (<it>DGAT1</it>) gene by using data sets including sires with accumulated phenotypic observations of daughters (breeding values, daughter yield deviations). However, these data sets prevented analyses with respect to dominance or parent-of-origin effects, although an increasing number of reports in the literature outlined the relevance of non-additive gene effects on quantitative traits.</p> <p>Results</p> <p>Based on a data set comprising German Holstein cows with direct trait measurements, we first confirmed the previously reported association of <it>DGAT1 </it>promoter VNTR alleles with milk production traits. We detected a dominant mode of effects for the <it>DGAT1 </it>K232A and promoter VNTR alleles. Namely, the contrasts between the effects of heterozygous individuals at the <it>DGAT1 </it>loci differed significantly from the midpoint between the effects for the two homozygous genotypes for several milk production traits, thus indicating the presence of dominance. Furthermore, we identified differences in the magnitude of effects between paternally and maternally inherited <it>DGAT1 </it>promoter VNTR – K232A haplotypes indicating parent-of-origin effects on milk production traits.</p> <p>Conclusion</p> <p>Non-additive effects like those identified at the bovine <it>DGAT1 </it>locus have to be accounted for in more specific QTL detection models as well as in marker assisted selection schemes. The <it>DGAT1 </it>alleles in cattle will be a useful model for further investigations on the biological background of non-additive effects in mammals due to the magnitude and consistency of their effects on milk production traits.</p

Risk thresholds for alcohol consumption : combined analysis of individual-participant data for 599 912 current drinkers in 83 prospective studies

Background Low-risk limits recommended for alcohol consumption vary substantially across different national guidelines. To define thresholds associated with lowest risk for all-cause mortality and cardiovascular disease, we studied individual-participant data from 599 912 current drinkers without previous cardiovascular disease. Methods We did a combined analysis of individual-participant data from three large-scale data sources in 19 high-income countries (the Emerging Risk Factors Collaboration, EPIC-CVD, and the UK Biobank). We characterised dose-response associations and calculated hazard ratios (HRs) per 100 g per week of alcohol (12.5 units per week) across 83 prospective studies, adjusting at least for study or centre, age, sex, smoking, and diabetes. To be eligible for the analysis, participants had to have information recorded about their alcohol consumption amount and status (ie, non-drinker vs current drinker), plus age, sex, history of diabetes and smoking status, at least 1 year of follow-up after baseline, and no baseline history of cardiovascular disease. The main analyses focused on current drinkers, whose baseline alcohol consumption was categorised into eight predefined groups according to the amount in grams consumed per week. We assessed alcohol consumption in relation to all-cause mortality, total cardiovascular disease, and several cardiovascular disease subtypes. We corrected HRs for estimated long-term variability in alcohol consumption using 152 640 serial alcohol assessments obtained some years apart (median interval 5.6 years [5th-95th percentile 1.04-13.5]) from 71 011 participants from 37 studies. Findings In the 599 912 current drinkers included in the analysis, we recorded 40 310 deaths and 39 018 incident cardiovascular disease events during 5.4 million person-years of follow-up. For all-cause mortality, we recorded a positive and curvilinear association with the level of alcohol consumption, with the minimum mortality risk around or below 100 g per week. Alcohol consumption was roughly linearly associated with a higher risk of stroke (HR per 100 g per week higher consumption 1.14, 95% CI, 1.10-1.17), coronary disease excluding myocardial infarction (1.06, 1.00-1.11), heart failure (1.09, 1.03-1.15), fatal hypertensive disease (1.24, 1.15-1.33); and fatal aortic aneurysm (1.15, 1.03-1.28). By contrast, increased alcohol consumption was loglinearly associated with a lower risk of myocardial infarction (HR 0.94, 0.91-0.97). In comparison to those who reported drinking >0-100-200-350 g per week had lower life expectancy at age 40 years of approximately 6 months, 1-2 years, or 4-5 years, respectively. Interpretation In current drinkers of alcohol in high-income countries, the threshold for lowest risk of all-cause mortality was about 100 g/week. For cardiovascular disease subtypes other than myocardial infarction, there were no clear risk thresholds below which lower alcohol consumption stopped being associated with lower disease risk. These data support limits for alcohol consumption that are lower than those recommended in most current guidelines. Copyright (C) The Author(s). Published by Elsevier Ltd.Peer reviewe

Discovery of common and rare genetic risk variants for colorectal cancer.

To further dissect the genetic architecture of colorectal cancer (CRC), we performed whole-genome sequencing of 1,439 cases and 720 controls, imputed discovered sequence variants and Haplotype Reference Consortium panel variants into genome-wide association study data, and tested for association in 34,869 cases and 29,051 controls. Findings were followed up in an additional 23,262 cases and 38,296 controls. We discovered a strongly protective 0.3% frequency variant signal at CHD1. In a combined meta-analysis of 125,478 individuals, we identified 40 new independent signals at P < 5 × 10-8, bringing the number of known independent signals for CRC to ~100. New signals implicate lower-frequency variants, Krüppel-like factors, Hedgehog signaling, Hippo-YAP signaling, long noncoding RNAs and somatic drivers, and support a role for immune function. Heritability analyses suggest that CRC risk is highly polygenic, and larger, more comprehensive studies enabling rare variant analysis will improve understanding of biology underlying this risk and influence personalized screening strategies and drug development.Goncalo R Abecasis has received compensation from 23andMe and Helix. He is currently an employee of Regeneron Pharmaceuticals. Heather Hampel performs collaborative research with Ambry Genetics, InVitae Genetics, and Myriad Genetic Laboratories, Inc., is on the scientific advisory board for InVitae Genetics and Genome Medical, and has stock in Genome Medical. Rachel Pearlman has participated in collaborative funded research with Myriad Genetics Laboratories and Invitae Genetics but has no financial competitive interest

Beyond Back Splicing, a Still Poorly Explored World: Non-Canonical Circular RNAs

Most of the circRNAs reported to date originate from back splicing of a pre-mRNA, and these exonic circRNAs are termed canonical circRNAs. Our objective was to provide an overview of all other (non-canonical) circRNAs that do not originate from the junction of two exons and to characterize their common properties. Those generated through a failure of intron lariat debranching are the best known, even though studies on them are rare. These circRNAs retain the 2&prime;&ndash;5&prime; bond derived from the intron lariat, and this feature probably explains the difficulties in obtaining efficient reverse transcription through the circular junction. Here, we provide an unprecedented overview of non-canonical circRNAs (lariat-derived intronic circRNAs, sub-exonic circRNAs, intron circles, tricRNAs), which all derive from non-coding sequences. As there are few data suggesting their involvement in cellular regulatory processes, we believe that it is early to propose a general function for circRNAs, even for lariat-derived circRNAs. We suggest that their small size and probably strong secondary structures could be major obstacles to their reliable detection. Nevertheless, we believe there are still several possible ways to advance our knowledge of this class of non-coding RNA

In-Depth Analysis Reveals Production of Circular RNAs from Non-Coding Sequences

International audienceThe sequencing of total RNA depleted for ribosomal sequences remains the method of choice for the study of circRNAs. Our objective was to characterize non-canonical circRNAs, namely not originating from back splicing and circRNA produced by non-coding genes. To this end, we analyzed a dataset from porcine testis known to contain about 100 intron-derived circRNAs. Labelling reads containing a circular junction and originating from back splicing provided information on the very small contribution of long non-coding genes to the production of canonical circRNAs. Analyses of the other reads revealed two origins for non-canonical circRNAs: (1) Intronic sequences for lariat-derived intronic circRNAs and intron circles, (2) Mono-exonic genes (mostly non-coding) for either a new type of circRNA (including only part of the exon: sub-exonic circRNAs) or, even more rarely, mono-exonic canonical circRNAs. The most complex set of sub-exonic circRNAs was produced byRNase_MRP(ribozyme RNA). We specifically investigated the intronic circRNA ofATXN2L, which is probably an independently transcribed sisRNA (stable intronic sequence RNA). We may be witnessing the emergence of a new non-coding gene in the porcine genome. Our results are evidence that most non-canonical circRNAs originate from non-coding sequences