Genome-wide association analysis of heifer livability and early first calving in Holstein cattle.

BACKGROUND: The survival and fertility of heifers are critical factors for the success of dairy farms. The mortality of heifers poses a significant challenge to the management and profitability of the dairy industry. In dairy farming, achieving early first calving of heifers is also essential for optimal productivity and sustainability. Recently, Council on Dairy Cattle Breeding (CDCB) and USDA have developed new evaluations of heifer health and fertility traits. However, the genetic basis of these traits has yet to be thoroughly studied. RESULTS: Leveraging the extensive U.S dairy genomic database maintained at CDCB, we conducted large-scale GWAS analyses of two heifer traits, livability and early first calving. Despite the large sample size, we found no major QTL for heifer livability. However, we identified a major QTL in the bovine MHC region associated with early first calving. Our GO analysis based on nearby genes detected 91 significant GO terms with a large proportion related to the immune system. This QTL in the MHC region was also confirmed in the analysis of 27 K bull with imputed sequence variants. Since these traits have few major QTL, we evaluated the genome-wide distribution of GWAS signals across different functional genomics categories. For heifer livability, we observed significant enrichment in promotor and enhancer-related regions. For early calving, we found more associations in active TSS, active Elements, and Insulator. We also identified significant enrichment of CDS and conserved variants in the GWAS results of both traits. By linking GWAS results and transcriptome data from the CattleGTEx project via TWAS, we detected four and 23 significant gene-trait association pairs for heifer livability and early calving, respectively. Interestingly, we discovered six genes for early calving in the Bovine MHC region, including two genes in lymph node tissue and one gene each in blood, adipose, hypothalamus, and leukocyte. CONCLUSION: Our large-scale GWAS analyses of two heifer traits identified a major QTL in the bovine MHC region for early first calving. Additional functional enrichment and TWAS analyses confirmed the MHC QTL with relevant biological evidence. Our results revealed the complex genetic basis of heifer health and fertility traits and indicated a potential connection between the immune system and reproduction in cattle.

Investigation of rumen long noncoding RNA before and after weaning in cattle.

BACKGROUND: This study aimed to identify long non-coding RNA (lncRNA) from the rumen tissue in dairy cattle, explore their features including expression and conservation levels, and reveal potential links between lncRNA and complex traits that may indicate important functional impacts of rumen lncRNA during the transition to the weaning period. RESULTS: A total of six cattle rumen samples were taken with three replicates from before and after weaning periods, respectively. Total RNAs were extracted and sequenced with lncRNA discovered based on size, coding potential, sequence homology, and known protein domains. As a result, 404 and 234 rumen lncRNAs were identified before and after weaning, respectively. However, only nine of them were shared under two conditions, with 395 lncRNAs found only in pre-weaning tissues and 225 only in post-weaning samples. Interestingly, none of the nine common lncRNAs were differentially expressed between the two weaning conditions. LncRNA averaged shorter length, lower expression, and lower conservation scores than the genome overall, which is consistent with general lncRNA characteristics. By integrating rumen lncRNA before and after weaning with large-scale GWAS results in cattle, we reported significant enrichment of both pre- and after-weaning lncRNA with traits of economic importance including production, reproduction, health, and body conformation phenotypes. CONCLUSIONS: The majority of rumen lncRNAs are uniquely expressed in one of the two weaning conditions, indicating a functional role of lncRNA in rumen development and transition of weaning. Notably, both pre- and post-weaning lncRNA showed significant enrichment with a variety of complex traits in dairy cattle, suggesting the importance of rumen lncRNA for cattle performance in the adult stage. These relationships should be further investigated to better understand the specific roles lncRNAs are playing in rumen development and cow performance.

Can routinely recorded reproductive events be used as indicators of disease emergence in dairy cattle? An evaluation of 5 indicators during the emergence of bluetongue virus in France in 2007 and 2008.

In response to increasing risks of emerging infectious diseases, syndromic surveillance can be a suitable approach to detect outbreaks of such diseases across a large territory in an early phase. To implement a syndromic surveillance system, the primary challenge is to find appropriate health-related data. The objective of this study was to evaluate whether routinely collected dates of reproductive events in dairy cattle could be used to build indicators of health anomalies for syndromic surveillance. The evaluation was performed on data collected in France between 2003 and 2009. First, a set of 5 indicators was proposed to assess several types of reproductive disorders. For each indicator, the demographic coverage over the total number of cattle at risk was analyzed in time and space. Second, the ability to detect an emerging disease in an early phase was retrospectively evaluated during epidemics of bluetongue serotypes 1 and 8 (BTV-1, BTV-8) in France in 2007 and 2008. Reproductive indicators were analyzed weekly during these epidemics for each indicator in each infected French district (16 in 2007 and 50 in 2008 out of 94 districts). The indicators were able to detect the BTV epidemics despite their low demographic coverage on a weekly basis relatively to total number of cattle (median=1.21%; range=0-11.7%). Four indicators related to abortions, late embryonic death, and short gestations were abnormally elevated during both BTV epidemics. Median times to abnormal elevations in these indicators were 20 to 71 d after the first notification of clinical signs of BTV by veterinarians. These results demonstrate that reproduction data can be used as indicators of disease emergences, whereas in the specific case of these BTV epidemics, detection via these indicators was later than clinical detection by veterinarians. The emergence of bluetongue in 2007 in France was associated with gestations that were a few days shorter than expected. A short gestation indicator underwent high elevations relative to prior random fluctuations and was the earliest (out of the 4 indicators) to show abnormal elevations, making it possible to detect this emergence.

Investigation of lncRNA in Bos taurus Mammary Tissue during Dry and Lactation Periods.

This study aims to collect RNA-Seq data from Bos taurus samples representing dry and lactating mammary tissue, identify lncRNA transcripts, and analyze findings for their features and functional annotation. This allows for connections to be drawn between lncRNA and the lactation process. RNA-Seq data from 103 samples of Bos taurus mammary tissue were gathered from publicly available databases (60 dry, 43 lactating). The samples were filtered to reveal 214 dry mammary lncRNA transcripts and 517 lactating mammary lncRNA transcripts. The lncRNAs met common lncRNA characteristics such as shorter length, fewer exons, lower expression levels, and less sequence conservation when compared to the genome. Interestingly, several lncRNAs showed sequence similarity to genes associated with strong hair keratin intermediate filaments. Human breast cancer research has associated strong hair keratin filaments with mammary tissue cellular resilience. The lncRNAs were also associated with several genes/proteins that linked to pregnancy using expression correlation and gene ontology. Such findings indicate that there are crucial relationships between the lncRNAs found in mammary tissue and the development of the tissue, to meet both the animal's needs and our own production needs; these relationships should be further investigated to ensure that we continue to breed the most resilient, efficient dairy cattle.

Application of syndromic surveillance on routinely collected cattle reproduction and milk production data for the early detection of outbreaks of Bluetongue and Schmallenberg viruses.

This study aimed to evaluate the use of routinely collected reproductive and milk production data for the early detection of emerging vector-borne diseases in cattle in the Netherlands and the Flanders region of Belgium (i.e., the northern part of Belgium). Prospective space-time cluster analyses on residuals from a model on milk production were carried out to detect clusters of reduced milk yield. A CUSUM algorithm was used to detect temporal aberrations in model residuals of reproductive performance models on two indicators of gestation length. The Bluetongue serotype-8 (BTV-8) epidemics of 2006 and 2007 and the Schmallenberg virus (SBV) epidemic of 2011 were used as case studies to evaluate the sensitivity and timeliness of these methods. The methods investigated in this study did not result in a more timely detection of BTV-8 and SBV in the Netherlands and BTV-8 in Belgium given the surveillance systems in place when these viruses emerged. This could be due to (i) the large geographical units used in the analyses (country, region and province level), and (ii) the high level of sensitivity of the surveillance systems in place when these viruses emerged. Nevertheless, it might be worthwhile to use a syndromic surveillance system based on non-specific animal health data in real-time alongside regular surveillance, to increase the sense of urgency and to provide valuable quantitative information for decision makers in the initial phase of an emerging disease outbreak.

Use of monthly collected milk yields for the detection of the emergence of the 2007 French BTV epizootic.

Two culicoides-borne diseases, Bluetongue (BTV) and Schmallenberg, have emerged in the European cattle population since 2006. Other diseases transmitted by these vectors could emerge. This justifies the development of syndromic surveillance programs whereby one or several indicators would be routinely monitored for the early detection of emerging diseases. The aim of this study was to evaluate milk yield from milk recording in dairy cattle as an indicator to be included in an emerging disease surveillance system. It was hypothesized that emergences would result in episodes of low milk production clustered in space and time. The 2007 BTV epizootic in France was used as a case study. Because it had already emerged in neighbouring countries, the disease emergence was expected and notification was mandatory. Herd-test-day milk productions were predicted for the entire country for 2006 and 2007 from herd historical data using linear mixed models. The differences between observed and predicted milk productions were averaged per week and per municipality and used as input for a space-time prospective scan statistic. Log likelihood ratios (LLR) associated with clusters were used to define alarms. The threshold chosen was a trade-off between detection timeliness and the number of false alarms per week. The first four BTV notifications occurred on the 12th (two notifications), 13th and 27th of July 2007. The 12th of July was considered to be the date of emergence. Alarms occurring before the 1st of March 2007 were considered to be false alarms. Using an LLR of 50, there were an average of 1.7 false alarms per week and the BTV emergence was detected seven weeks after emergence. Using an LLR of 100, there were an average of 0.8 false alarms per week and the BTV emergence was detected 9 weeks after emergence. Detection may have been delayed because of a discontinuation of milk recording between mid-July and mid-August. The first cluster with an LLR>100 located in the emergence area was further investigated. A difference between observed and predicted production of >1 kg/cow/day was observed around the time of emergence. However, a difference of equal magnitude was observed during the year preceding the outbreak. Milk production predicted from herd history alone did not allow the detection of the 2007 BTV emergence in France. Further research should be conducted on improving the prediction of test-day milk yield and on combining it with other indicators based on routinely collected data.

Evaluation of a continuous indicator for syndromic surveillance through simulation. application to vector borne disease emergence detection in cattle using milk yield.

Two vector borne diseases, caused by the Bluetongue and Schmallenberg viruses respectively, have emerged in the European ruminant populations since 2006. Several diseases are transmitted by the same vectors and could emerge in the future. Syndromic surveillance, which consists in the routine monitoring of indicators for the detection of adverse health events, may allow an early detection. Milk yield is routinely measured in a large proportion of dairy herds and could be incorporated as an indicator in a surveillance system. However, few studies have evaluated continuous indicators for syndromic surveillance. The aim of this study was to develop a framework for the quantification of both disease characteristics and model predictive abilities that are important for a continuous indicator to be sensitive, timely and specific for the detection of a vector-borne disease emergence. Emergences with a range of spread characteristics and effects on milk production were simulated. Milk yields collected monthly in 48 713 French dairy herds were used to simulate 576 disease emergence scenarios. First, the effect of disease characteristics on the sensitivity and timeliness of detection were assessed: Spatio-temporal clusters of low milk production were detected with a scan statistic using the difference between observed and simulated milk yields as input. In a second step, the system specificity was evaluated by running the scan statistic on the difference between observed and predicted milk yields, in the absence of simulated emergence. The timeliness of detection depended mostly on how easily the disease spread between and within herds. The time and location of the emergence or adding random noise to the simulated effects had a limited impact on the timeliness of detection. The main limitation of the system was the low specificity i.e. the high number of clusters detected from the difference between observed and predicted productions, in the absence of disease.