Transmission ratio distortion regions in the context of genomic evaluation and their effects on reproductive traits in cattle.

Transmission ratio distortion (TRD), which is a deviation from Mendelian expectations, has been associated with basic mechanisms of life such as sperm and ova fertility and viability at developmental stages of the reproductive cycle. In this study different models including TRD regions were tested for different reproductive traits [days from first service to conception (FSTC), number of services, first service nonreturn rate (NRR), and stillbirth (SB)]. Thus, in addition to a basic model with systematic and random effects, including genetic effects modeled through a genomic relationship matrix, we developed 2 additional models, including a second genomic relationship matrix based on TRD regions, and TRD regions as a random effect assuming heterogeneous variances. The analyses were performed with 10,623 cows and 1,520 bulls genotyped for 47,910 SNPs, 590 TRD regions, and several records ranging from 9,587 (FSTC) to 19,667 (SB). The results of this study showed the ability of TRD regions to capture some additional genetic variance for some traits; however, this did not translate into higher accuracy for genomic prediction. This could be explained by the nature of TRD itself, which may arise in different stages of the reproductive cycle. Nevertheless, important effects of TRD regions were found on SB (31 regions) and NRR (18 regions) when comparing at-risk versus control matings, especially for regions with allelic TRD pattern. Particularly for NRR, the probability of observing nonpregnant cow increases by up to 27% for specific TRD regions, and the probability of observing stillbirth increased by up to 254%. These results support the relevance of several TRD regions on some reproductive traits, especially those with allelic patterns that have not received as much attention as recessive TRD patterns.

Identification of novel alternative splicing associated with mastitis disease in Holstein dairy cows using large gap read mapping.

BACKGROUND: Mastitis is a very common disease in the dairy industry that producers encounter daily. Transcriptomics, using RNA-Sequencing (RNA - Seq) technology, can be used to study the functional aspect of mastitis resistance to identify animals that have a better immune response to mastitis. When the cow has mastitis, not only genes but also specific mRNA isoforms generated via alternative splicing (AS) could be differentially expressed (DE), leading to the phenotypic variation observed. Therefore, the objective of this study was to use large gap read mapping to identify mRNA isoforms DE between healthy and mastitic milk somatic cell samples (N = 12). These mRNA isoforms were then categorized based on being 1) annotated mRNA isoforms for gene name and length, 2) annotated mRNA isoforms with different transcript length and 3) novel mRNA isoforms of non - annotated genes. RESULTS: Analysis identified 333 DE transcripts (with at least 2 mRNA isoforms annotated, with at least one being DE) between healthy and mastitic samples corresponding to 303 unique genes. Of these 333 DE transcripts between healthy and mastitic samples, 68 mRNA isoforms are annotated in the bovine genome reference (ARS.UCD.1.2), 249 mRNA isoforms had novel transcript lengths of known genes and 16 were novel transcript lengths of non - annotated genes in the bovine genome reference (ARS.UCD.1.2). Functional analysis including gene ontology, gene network and metabolic pathway analysis was performed on the list of 288 annotated and unique DE mRNA isoforms. In total, 67 significant metabolic pathways were identified including positive regulation of cytokine secretion and immune response. Additionally, numerous DE novel mRNA isoforms showed potential involvement with the immune system or mastitis. Lastly, QTL annotation analysis was performed on coding regions of the DE mRNA isoforms, identifying overlapping QTLs associated with clinical mastitis and somatic cell score. CONCLUSION: This study identified novel mRNA isoforms generated via AS that could lead to differences in the immune response of Holstein dairy cows and be potentially implemented in future breeding programs.

Genome-wide association study using haplotype libraries and repeated-measures model to identify candidate genomic regions for stillbirth in Holstein cattle.

Reduced fertility is one of the main causes of economic losses on dairy farms, resulting in economic losses estimated at $938 per stillbirth case in Holstein herds. The identification of genomic regions associated with stillbirth could help to develop better management and breeding strategies aimed to reduce the frequency of undesirable gestation outcomes. Here, 10,570 cows and 50,541 birth records were used to perform a haplotype-based GWAS. A total of 41 significantly associated pseudo-SNPs (haplotypes within haplotype blocks converted to a binary classification) were identified after Bonferroni adjustment for multiple tests. A total of 117 positional candidate genes were annotated within or close (in a 200-kb interval) to significant pseudo-SNPs (haplotype blocks). The guilt-by-association functional prioritization identified 31 potential functional candidate genes for reproductive performance out of the 117 positional candidate genes annotated. These genes play crucial roles in biological processes associated with pregnancy persistence, fetus development, immune response, among others. These results helped us to better understand the genetic basis of stillbirth in dairy cattle and may be useful for the prediction of stillbirth in Holstein cattle, helping to reduce the related economic losses caused by this phenotype.

Differential gene expression in dairy cows under negative energy balance and ketosis: A systematic review and meta-analysis.

Development of ketosis in high-producing dairy cows contributes to several animal health issues and highlights the need for a better understanding of the genetic basis of metabolic diseases. To evaluate the pattern of differential gene expression in the liver of cows under negative energy balance (NEB), and under subclinical and clinical ketosis, a meta-analysis of gene expression and genome-wide association studies results was performed. An initial systematic review identified 118 articles based on the key words "cow," "liver," "negative energy balance," "ketosis," "expression," "qPCR," "microarray," "proteomic," "RNA-Seq," and "GWAS." After further screening for only peer-reviewed and pertinent articles for gene expression during NEB and clinical and subclinical ketosis (considering plasma levels of ß-hydroxybutyrate), 20 articles were included in the analysis. From the systematic review, 430 significant SNPs identified by genome-wide association studies (GWAS) were assigned to genes reported in gene expression studies by considering chromosome and base pair positions in the ARS-UCD 1.2 bovine assembly. Venn diagrams were created to integrate the data obtained in the systematic review, and Gene Ontology enrichment analysis was carried out using official gene names. A QTL enrichment analysis was also performed to identify potential positional candidate loci. Twenty-four significant SNPs were located within the coordinates of differentially expressed genes located on chromosomes 2, 3, 6, 9, 11, 14, 27, and 29. Three significant metabolic pathways were associated with NEB and subclinical and clinical ketosis. In addition, 2 important genes, PPARA (peroxisome proliferator activated receptor alpha) and ACACA (acetyl-coenzyme A carboxylase α), were identified, which were differentially expressed in the 3 metabolic conditions. The PPARA gene is involved in the regulation of lipid metabolism and fatty liver disease and the ACACA gene encodes an enzyme that catalyzes the carboxylation of acetyl-coenzyme A to malonyl-coenzyme A, which is a rate-limiting step in fatty acid synthesis. Gene network analysis revealed co-expression interactions among 34 genes associated with functions involving fatty acid transport and fatty acid metabolism. For the annotated QTL, 9 QTL were identified for ketosis. The genes FN1 (fibronectin 1) and PTK2 (protein tyrosine kinase 2), which are mainly involved in cell adhesion and formation of extracellular matrix constituents, were enriched for QTL previously associated with the trait "ketosis" on chromosome 2 and for the trait "milk iron content" on chromosome 14, respectively. This integration of gene expression and GWAS data provides an additional understanding of the genetic background of NEB and subclinical and clinical ketosis in dairy cattle. Thus, it is a useful approach to identify biological mechanisms underlying these metabolic conditions in dairy cattle.

Identification of functional candidate variants and genes for feed efficiency in Holstein and Jersey cattle breeds using RNA-sequencing.

The identification of functional genetic variants and associated candidate genes linked to feed efficiency may help improve selection for feed efficiency in dairy cattle, providing economic and environmental benefits for the dairy industry. This study used RNA-sequencing data obtained from liver tissue from 9 Holstein cows [n = 5 low residual feed intake (RFI), n = 4 high RFI] and 10 Jersey cows (n = 5 low RFI, n = 5 high RFI), which were selected from a single population of 200 animals. Using RNA-sequencing, 3 analyses were performed to identify: (1) variants within low or high RFI Holstein cattle; (2) variants within low or high RFI Jersey cattle; and (3) variants within low or high RFI groups, which are common across both Holstein and Jersey cattle breeds. From each analysis, all variants were filtered for moderate, modifier, or high functional effect, and co-localized quantitative trait loci (QTL) classes, enriched biological processes, and co-localized genes related to these variants, were identified. The overlapping of the resulting genes co-localized with functional SNP from each analysis in both breeds for low or high RFI groups were compared. For the first two analyses, the total number of candidate genes associated with moderate, modifier, or high functional effect variants fixed within low or high RFI groups were 2,810 and 3,390 for Holstein and Jersey breeds, respectively. The major QTL classes co-localized with these variants included milk and reproduction QTL for the Holstein breed, and milk, production, and reproduction QTL for the Jersey breed. For the third analysis, the common variants across both Holstein and Jersey breeds, uniquely fixed within low or high RFI groups were identified, revealing a total of 86,209 and 111,126 functional variants in low and high RFI groups, respectively. Across all 3 analyses for low and high RFI cattle, 12 and 31 co-localized genes were overlapping, respectively. Among the overlapping genes across breeds, 9 were commonly detected in both the low and high RFI groups (INSRR, CSK, DYNC1H1, GAB1, KAT2B, RXRA, SHC1, TRRAP, PIK3CB), which are known to play a key role in the regulation of biological processes that have high metabolic demand and are related to cell growth and regeneration, metabolism, and immune function. The genes identified and their associated functional variants may serve as candidate genetic markers and can be implemented into breeding programs to help improve the selection for feed efficiency in dairy cattle.

Associations between feed efficiency and aspects of lactation curves in primiparous Holstein dairy cattle.

Genetic selection for improved feed efficiency in dairy cattle has received renewed attention over the last decade to address the needs of a growing global population. As milk yield is a critical component of feed efficiency metrics in dairy animals, our objective was to evaluate the associations between feed efficiency in primiparous Holstein cattle and parameters of a mathematical model describing individual lactation curves. The Dijkstra lactation curve model was fit to individual lactation records from 34 Holstein heifers with previously estimated measures of feed efficiency. We found that the optimal fit of the Dijkstra model was achieved using daily milk yield records up to 21 d in milk to capture the rise to peak milk yield and using monthly dairy herd improvement records for the remainder of lactation to accurately characterize lactation persistency. In the period of lactation before peak milk yield, improved feed efficiency was associated with a faster increase in daily milk yield over a shorter period of time at the expense of increased mobilization of body reserves; this serves to reinforce the concept that dairy cattle are primarily capital breeders versus income breeders. Feed efficiency in the period following peak lactation, as measured by gross feed efficiency, return over feed costs, and net energy efficiency of lactation, was positively associated with higher peak milk yield. The findings in early lactation suggest that estimates of feed efficiency could be improved by evaluating feed efficiency relative to conception, rather than parturition and lactation, to better account for the energy stored and released from body reserves in capital breeding.

Development and comparison of RNA-sequencing pipelines for more accurate SNP identification: practical example of functional SNP detection associated with feed efficiency in Nellore beef cattle.

BACKGROUND: Optimization of an RNA-Sequencing (RNA-Seq) pipeline is critical to maximize power and accuracy to identify genetic variants, including SNPs, which may serve as genetic markers to select for feed efficiency, leading to economic benefits for beef production. This study used RNA-Seq data (GEO Accession ID: PRJEB7696 and PRJEB15314) from muscle and liver tissue, respectively, from 12 Nellore beef steers selected from 585 steers with residual feed intake measures (RFI; n = 6 low-RFI, n = 6 high-RFI). Three RNA-Seq pipelines were compared including multi-sample calling from i) non-merged samples; ii) merged samples by RFI group, iii) merged samples by RFI and tissue group. The RNA-Seq reads were aligned against the UMD3.1 bovine reference genome (release 94) assembly using STAR aligner. Variants were called using BCFtools and variant effect prediction (VeP) and functional annotation (ToppGene) analyses were performed. RESULTS: On average, total reads detected for Approach i) non-merged samples for liver and muscle, were 18,362,086.3 and 35,645,898.7, respectively. For Approach ii), merging samples by RFI group, total reads detected for each merged group was 162,030,705, and for Approach iii), merging samples by RFI group and tissues, was 324,061,410, revealing the highest read depth for Approach iii). Additionally, Approach iii) merging samples by RFI group and tissues, revealed the highest read depth per variant coverage (572.59 ± 3993.11) and encompassed the majority of localized positional genes detected by each approach. This suggests Approach iii) had optimized detection power, read depth, and accuracy of SNP calling, therefore increasing confidence of variant detection and reducing false positive detection. Approach iii) was then used to detect unique SNPs fixed within low- (12,145) and high-RFI (14,663) groups. Functional annotation of SNPs revealed positional candidate genes, for each RFI group (2886 for low-RFI, 3075 for high-RFI), which were significantly (P < 0.05) associated with immune and metabolic pathways. CONCLUSION: The most optimized RNA-Seq pipeline allowed for more accurate identification of SNPs, associated positional candidate genes, and significantly associated metabolic pathways in muscle and liver tissues, providing insight on the underlying genetic architecture of feed efficiency in beef cattle.

Acute-onset smell and taste disorders in the context of COVID-19: a pilot multicentre polymerase chain reaction based case-control study.

BACKGROUND AND PURPOSE: Specific respiratory tract infections, including COVID-19, may cause smell and/or taste disorders (STDs) with increased frequency. The aim was to determine whether new-onset STDs are more frequent amongst COVID-19 patients than influenza patients. METHOD: This was a case-control study including hospitalized patients of two tertiary care centres. Consecutive patients positive for COVID-19 polymerase chain reaction (cases) and patients positive for influenza polymerase chain reaction (historical control sample) were assessed during specific periods, employing a self-reported STD questionnaire. RESULTS: Seventy-nine cases and 40 controls were included. No significant differences were found in basal features between the two groups. New-onset STDs were significantly more frequent amongst cases (31, 39.2%) than in the control group (5, 12.5 %) [adjusted odds ratio 21.4 (2.77-165.4, P = 0.003)]. COVID-19 patients with new-onset STDs were significantly younger than COVID-19 patients without STDs (52.6 ± 17.2 vs. 67.4 ± 15.1, P < 0.001). Amongst COVID-19 patients who presented STDs, 22 (70.9%) recalled an acute onset and it was an initial manifestation in 11 (35.5%). Twenty-five (80.6%) presented smell disorders (mostly anosmia, 14, 45.2%) and 28 (90.3%) taste disorders (mostly ageusia, 14, 45.2%). Only four (12.9 %) reported concomitant nasal obstruction. The mean duration of STD was 7.5 ± 3.2 days and 12 patients (40%) manifested complete recovery after 7.4 ± 2.3 days of onset. CONCLUSION: New-onset STDs were significantly more frequent amongst COVID-19 patients than influenza patients; they usually had an acute onset and were commonly an initial manifestation. The use of STD assessment in anamnesis as a hint for COVID-19 and to support individuals' self-isolation in the current epidemic context is suggested.

Discriminating between allele- and genotype-specific transmission ratio distortion.

Transmission ratio distortion (TRD) is defined as the observed deviation from the expected Mendelian inheritance of alleles from heterozygous parents. This phenomenon is attributed to various biological mechanisms acting on germ cells, embryos or fetuses, or even in early postnatal life. Current statistical approaches typically use two independent parametrizations assuming that TRD relies on allele- or genotype-related mechanisms, although they have never been tested and compared. This study compared allele- and genotype-related TRD models on simulated datasets with 1000 genotyped offspring and real data from 168 sire-dam-offspring beef cattle trios. The analysis of simulated datasets favored the true model of analysis in most cases (>93%), and a low percentage of missidentification occurred under (almost) null dominance (genotype-related model) or similar and moderate-to-low sire- and dam-specific TRD parameters (allele-related model). Moreover, the correlation between simulated and predicted distortion parameters was high (>0.97) under the true model. The comparison of allele- and genotype-related TRD models is an appealing tool to infer the biological source of TRD (i.e. haploid vs. diploid cells) when screening the whole genome. The analysis of beef cattle data corroborated a TRD region previously reported in chromosome 4, although discarding allele-related mechanisms and favoring the genotype-related model as the more reliable one. The results of this study highlight the relevance of implementing and comparing different parametrizations to capture all kinds of TRD, and to compare them using appropriate statistical methods.

Discovering lethal alleles across the turkey genome using a transmission ratio distortion approach.

Deviation from Mendelian inheritance expectations (transmission ratio distortion, TRD) has been observed in several species, including the mouse and humans. In this study, TRD was characterized in the turkey genome using both allelic (specific- and unspecific-parent TRD) and genotypic (additive- and dominance-TRD) parameterizations within a Bayesian framework. In this study, we evaluated TRD for 23 243 genotyped Turkeys across 56 393 autosomal SNPs. The analyses included 500 sires, 2013 dams and 11 047 offspring (trios). Three different haplotype sliding windows of 4, 10 and 20 SNPs were used across the autosomal chromosomes. Based on the genotypic parameterizations, 14 haplotypes showed additive and dominance TRD effects highlighting regions with a recessive TRD pattern. In contrast, the allelic model uncovered 12 haplotype alleles with the allelic TRD pattern which showed an underrepresentation of heterozygous offspring in addition to the absence of homozygous animals. For regions with the allelic pattern, only one particular region showed a parent-specific TRD where the penetrance was high via the dam, but low via the sire. The gene set analysis uncovered several gene ontology functional terms, Reactome pathways and several Medical Subject Headings that showed significant enrichment of genes associated with TRD. Many of these gene ontology functional terms (e.g. mitotic spindle assembly checkpoint, DRM complex and Aneuploidy), Reactome pathways (e.g. Mismatch repair) and Medical Subject Headings (e.g. Adenosine monophosphate) are known to be related to fertility, embryo development and lethality. The results of this study revealed potential novel candidate lethal haplotypes, functional terms and pathways that may enhance breeding programs in Turkeys through reducing mortality and improving reproduction rate.

Short communication: Characterizing ovine serum stress biomarkers during endotoxemia.

Breeding stress-resilient livestock is a potential strategy to help mitigate the negative effect of environmental and pathogenic stressors. The hypothalamic-pituitary-adrenal axis and immune system are activated during stress events and release mediators into the circulation that help restore physiological homeostasis. The purpose of this study was to assess a comprehensive set of circulatory mediators released in response to an acute immune stress challenge to identify candidate biomarkers that can be used for the selection of stress-resilient animals. Fifteen female lambs were stress challenged with an intravenous bolus of lipopolysaccharide (LPS; 400 ng/kg), and blood was collected from the jugular vein at 0, 2, 4, and 6 h after LPS challenge to identify and monitor candidate stress biomarkers; temperature was also recorded over time. Biomarker responses were evaluated with a repeated-measures model to compare time points with baseline values. As expected, all sheep had a monophasic febrile response to LPS challenge, and cortisol increased and returned to baseline by 6 h. The cytokines tumor necrosis factor-α, IL-6, IFN-γ (proinflammatory), and IL-10 (anti-inflammatory) increased, but only tumor necrosis factor-α returned to baseline during the monitoring period. The cytokines IL-1α, IL-1ß, IL-17α (proinflammatory), and IL-4 (anti-inflammatory) did not respond to LPS challenge. All chemokines (CCL2, CCL3, CCL4, CXCL10, and IL-8) responded to LPS challenge; however, only CCL2, CCL3, CCL4, and CXCL10 increased over time, and only CCL3, CCL4, and CXCL10 returned to baseline during the monitoring period. MicroRNA (miR-145, miR-233, and miR-1246) also increased and remained elevated during the study. In summary, the LPS challenge induced a strong stress response in Rideau-Dorset sheep that could be monitored with a distinct profile of circulatory biomarkers.

The dynamic behavior of feed efficiency in primiparous dairy cattle.

A variety of measures are currently used to evaluate feed efficiency in livestock, with some of the most popular being the linear regression of nutrient intake on the various energy sinks of the cohort being evaluated. The adoption of novel on-farm technologies has greatly increased the breadth and depth of data being collected on a daily basis, which creates the opportunity for the implementation of more advanced feed efficiency models. The objectives of the current study were to leverage data collected using high-throughput biometric technologies to evaluate changes of early-lactation feed efficiency on a daily basis and to develop a novel measure of feed efficiency based on the NRC net energy model. Daily records of nutrient intake, fat- and protein-corrected milk yield, body weight, and body condition score were used to model 7 different indices of feed efficiency (gross feed efficiency, return over feed costs, residual feed intake, residual net energy intake, and net energy efficiencies of lactation, maintenance, and body reserve flux) in primiparous Holstein heifers (n = 40) from 0 to 150 d in milk. We observed that feed efficiency was dynamic and fluctuated over the course of the experimental period regardless of the efficiency measure used, but most measures followed similar trends in changes of efficiency on a day-to-day basis. Of the measures evaluated, return over feed costs would be most suited for management purposes due to practicality and simplicity, whereas other measures, such as the efficiency of net energy flux to and from body reserves, present potential targets for inclusion in breeding goals. Net energy efficiency of lactation was estimated to have higher repeatability compared with residual feed or energy intake, which may reduce the time required to estimate whole-lactation efficiency. Because decreased net energy balance was strongly correlated with increased efficiency for all indices evaluated, future work investigating the associations of feed efficiency with the magnitude and duration of negative energy balance is warranted.

Genome-wide association analysis for ß-hydroxybutyrate concentration in Milk in Holstein dairy cattle.

BACKGROUND: Ketosis in dairy cattle has been shown to cause a high morbidity in the farm and substantial financial losses to dairy farmers. Ketosis symptoms, however, are difficult to identify, therefore, the amount of ketone bodies (mainly ß-hydroxybutyric acid, BHB) is used as an indicator of subclinical ketosis in cows. It has also been shown that milk BHB concentrations have a strong correlation with ketosis in dairy cattle. Mid-infrared spectroscopy (MIR) has recently became a fast, cheap and high-throughput method for analyzing milk components. The aim of this study was to perform a genome-wide association study (GWAS) on the MIR-predicted milk BHB to identify genomic regions, genes and pathways potentially affecting subclinical ketosis in North American Holstein dairy cattle. RESULTS: Several significant regions were identified associated with MIR-predicted milk BHB concentrations (indicator of subclinical ketosis) in the first lactation (SCK1) and second and later lactations (SCK2) in Holstein dairy cows. The strongest association was located on BTA6 for SCK1 and BTA14 on SCK2. Several SNPs on BTA6 were identified in regions and variants reported previously to be associated with susceptibility to ketosis and clinical mastitis in Jersey and Holstein dairy cattle, respectively. One highly significant SNP on BTA14 was found within the DGAT1 gene with known functions on fat metabolism and inflammatory response in dairy cattle. A region on BTA6 and three SNPs on BTA20 were found to overlap between SCK1 and SCK2. However, a novel region on BTA20 (55-63 Mb) for SCK2 was also identified, which was not reported in previous association studies. Enrichment analysis of the list of candidate genes within the identified regions for MIR-predicted milk BHB concentrations yielded molecular functions and biological processes that may be involved in the inflammatory response and lipid metabolism in dairy cattle. CONCLUSIONS: The results of this study confirmed several SNPs and genes identified in previous studies as associated with ketosis susceptibility and immune response, and also found a novel region that can be used for further analysis to identify causal variations and key regulatory genes that affect clinical/ subclinical ketosis.

Genetic and genomic analyses of embryo production in dairy cattle.

The Canadian dairy industry has been using invivo and invitro assisted reproductive technologies to produce embryos. Technological improvements have helped increase the number and quality of embryos produced, but genetic and genomic tools for improving these traits have yet to be assessed for the Canadian Holstein population. Genetic parameters and a genome-wide association study were performed in Canadian Holstein for the total number of embryos (NE) and the number of viable embryos (VE). Results showed potential for genetic selection for both NE and VE, with heritability estimates (± s.e.) of approximately 0.15±0.01. Genetic correlations between the number of embryos produced using different procedures (invivo and invitro) suggested that a similar number of embryos should be expected from a donor regardless of the procedure used. A region on chromosome 11 of the bovine genome was found to be significantly associated with the number of embryos, indicating a potential regulatory role of this region on embryo production. Overall, these findings are of interest for the Canadian dairy industry because they provide useful information for breeders that are interested in producing embryos from the elite donors in their herds or in the population using assisted reproductive technologies.

Genetic mechanisms regulating the host response during mastitis.

Mastitis is a very costly and common disease in the dairy industry. The study of the transcriptome from healthy and mastitic milk somatic cell samples using RNA-Sequencing technology can provide measurements of transcript levels associated with the immune response to the infection. The objective of this study was to characterize the Holstein milk somatic cell transcriptome from 6 cows to determine host response to intramammary infections. RNA-Sequencing was performed on 2 samples from each cow from 2 separate quarters, one classified as healthy (n = 6) and one as mastitic (n = 6). In total, 449 genes were differentially expressed between the healthy and mastitic quarters (false discovery rate <0.05, fold change >±2). Among the differentially expressed genes, the most expressed genes based on reads per kilobase per million mapped reads (RPKM) in the healthy group were associated with milk components (CSN2 and CSN3), and in the mastitic group they were associated with immunity (B2M and CD74). In silico functional analysis was performed using the list of 449 differentially expressed genes, which identified 36 significantly enriched metabolic pathways (false discovery rate <0.01), some of which were associated with the immune system, such as cytokine-cytokine interaction and cell adhesion molecules. Seven functional candidate genes were selected, based on the criteria of being highly differentially expressed between healthy and mastitic groups and significantly enriched in metabolic pathways that are relevant to the inflammatory process (GLYCAM1, B2M, CD74, BoLA-DRA, FCER1G, SDS, and NFKBIA). Last, we identified the differentially expressed genes that are located in quantitative trait locus regions previously known to be associated with mastitis, specifically clinical mastitis, somatic cell count, and somatic cell score. It was concluded that multiple genes within quantitative trait locus regions could potentially affect host response to mastitis-causing agents, making some cows more susceptible to intramammary infections. The identification of potential candidate genes with functional, statistical, biological, and positional relevance associated with host defense to infection will contribute to a better understanding of the underlying genetic architecture associated with mastitis. This in turn will improve the sustainability of agricultural practices by facilitating the selection of cows with improved host defense leading to increased resistance to mastitis.

Implementation of Bayesian methods to identify SNP and haplotype regions with transmission ratio distortion across the whole genome: TRDscan v.1.0.

Realized deviations from the expected Mendelian inheritance of alleles from heterozygous parents have been previously reported in a broad range of organisms (i.e., transmission ratio distortion; TRD). Various biological mechanisms affecting gametes, embryos, fetuses, or even postnatal offspring can produce patterns of TRD. However, knowledge about its prevalence and potential causes in livestock species is still scarce. Specific Bayesian models have been recently developed for the analyses of TRD for biallelic loci, which accommodated a wide range of population structures, enabling TRD investigation in livestock populations. The parameterization of these models is flexible and allows the study of overall (parent-unspecific) TRD and sire- and dam-specific TRD. This research aimed at deriving Bayesian models for fitting TRD on the basis of haplotypes, testing the models for both haplotype- and SNP-based methods in simulated data and actual Holstein genotypes, and developing a specific software for TRD analyses. Results obtained on simulated data sets showed that the statistical power of the analysis increased with sample size of trios (n), proportion of heterozygous parents, and the magnitude of the TRD. On the other hand, the statistical power to detect TRD decreased with the number of alleles at each loci. Bayesian analyses showed a strong Pearson correlation coefficient (≥0.97) between simulated and estimated TRD that reached the significance level of Bayes factor ≥10 for both single-marker and haplotype analyses when n ≥ 25. Moreover, the accuracy in terms of the mean absolute error decreased with the increase of the sample size and increased with the number of alleles at each loci. Using real data (55,732 genotypes of Holstein trios), SNP- and haplotype-based distortions were detected with overall TRD, sire-TRD, or dam-TRD, showing different magnitudes of TRD and statistical relevance. Additionally, the haplotype-based method showed more ability to capture TRD compared with individual SNP. To discard possible random TRD in real data, an approximate empirical null distribution of TRD was developed. The program TRDscan v.1.0 was written in Fortran 2008 language and provides a powerful statistical tool to scan for TRD regions across the whole genome. This developed program is freely available at http://www.casellas.info/files/TRDscan.zip.

Invited review: Determination of large-scale individual dry matter intake phenotypes in dairy cattle.

Feed efficiency has been widely studied in many areas of dairy science and is currently seeing renewed interest in the field of breeding and genetics. A critical part of determining how efficiently an animal utilizes feed is accurately measuring individual dry matter (DM) intake. Currently, multiple methods are used to measure feed intake or determine the DM content of that feed, resulting in different levels of accuracy of measurement. Furthermore, the scale at which data need to be collected for use in genetic analyses makes some methodologies impractical. This systematic review aims to provide an overview of the current methodologies used to measure both feed intake in ruminants and DM content of feedstuffs, current methods to predict individual DM intake, and applications of large-scale intake measurements. Overall, advances in milk spectral data analysis present a promising method of estimating individual DM intake on a herd scale with further validation of prediction models. Although measurements of individual feed intake rely on the same underlying principle, the methods selected are largely dictated by the costs of capital, labor, and necessary analyses. Finally, DM methodologies were synthesized into a comprehensive protocol for use in a variety of feedstuffs.

Pre- and post-puberty expression of genes and proteins in the uterus of Bos indicus heifers: the luteal phase effect post-puberty.

Progesterone signaling and uterine function are crucial in terms of pregnancy establishment. To investigate how the uterine tissue and its secretion changes in relation to puberty, we sampled tissue and uterine fluid from six pre- and six post-pubertal Brahman heifers. Post-pubertal heifers were sampled in the luteal phase. Gene expression of the uterine tissue was investigated with RNA-sequencing, whereas the uterine fluid was used for protein profiling with mass spectrometry. A total of 4034 genes were differentially expressed (DE) at a nominal P-value of 0.05, and 26 genes were significantly DE after Bonferroni correction (P < 3.1 × 10-6 ). We also identified 79 proteins (out of 230 proteins) that were DE (P < 1 × 10-5 ) in the uterine fluid. When we compared proteomics and transcriptome results, four DE proteins were identified as being encoded by DE genes: OVGP1, GRP, CAP1 and HBA. Except for CAP1, the other three had lower expression post-puberty. The function of these four genes hypothetically related to preparation of the uterus for a potential pregnancy is discussed in the context of puberty. All DE genes and proteins were also used in pathway and ontology enrichment analyses to investigate overall function. The DE genes were enriched for terms related to ribosomal activity. Transcription factors that were deemed key regulators of DE genes are also reported. Transcription factors ZNF567, ZNF775, RELA, PIAS2, LHX4, SOX2, MEF2C, ZNF354C, HMG20A, TCF7L2, ZNF420, HIC1, GTF3A and two novel genes had the highest regulatory impact factor scores. These data can help to understand how puberty influences uterine function.

SNP detection using RNA-sequences of candidate genes associated with puberty in cattle.

Fertility traits, such as heifer pregnancy, are economically important in cattle production systems, and are therefore, used in genetic selection programs. The aim of this study was to identify single nucleotide polymorphisms (SNPs) using RNA-sequencing (RNA-Seq) data from ovary, uterus, endometrium, pituitary gland, hypothalamus, liver, longissimus dorsi muscle, and adipose tissue in 62 candidate genes associated with heifer puberty in cattle. RNA-Seq reads were assembled to the bovine reference genome (UMD 3.1.1) and analyzed in five cattle breeds; Brangus, Brahman, Nellore, Angus, and Holstein. Two approaches used the Brangus data for SNP discovery 1) pooling all samples, and 2) within each individual sample. These approaches revealed 1157 SNPs. These were compared with those identified in the pooled samples of the other breeds. Overall, 172 SNPs within 13 genes (CPNE5, FAM19A4, FOXN4, KLF1, LOC777593, MGC157266, NEBL, NRXN3, PEPT-1, PPP3CA, SCG5, TSG101, and TSHR) were concordant in the five breeds. Using Ensembl's Variant Effector Predictor, we determined that 12% of SNPs were in exons (71% synonymous, 29% nonsynonymous), 1% were in untranslated regions (UTRs), 86% were in introns, and 1% were in intergenic regions. Since these SNPs were discovered in RNA, the variants were predicted to be within exons or UTRs. Overall, 160 novel transcripts in 42 candidate genes and five novel genes overlapping five candidate genes were observed. In conclusion, 1157 SNPs were identified in 62 candidate genes associated with puberty in Brangus cattle, of which, 172 were concordant in the five cattle breeds. Novel transcripts and genes were also identified.

Reliability of the minimum basic dataset for diagnoses of cerebrovascular disease. / Fiabilidad del conjunto mínimo básico de datos en el diagnóstico de la enfermedad cerebrovascular.

INTRODUCTION: The minimum basic dataset is the largest available hospital care administrative database that is used in clinical studies and hospital management in association with diagnosis-related groups (DRGs). In 2011, the quality of the national MBDS in hospital discharges was audited, in order to assess its reliability. This paper presents a sub-analysis of the results from that analysis which are referred to cerebrovascular disease (CVD). METHODS: Using all discharge reports from the Spanish MBDS in 2009, a representative sample was obtained by stratified sampling and 11 209 records were evaluated. Outcome indicators were obtained to measure any differences observed between the national MBDS being evaluated and the hospital's original MBDS. Analysis of codes for CVD as a primary diagnosis was performed for ICD-9-CM diagnostic categories 430 through 438. We evaluated error rates in the selection and classification of main diagnoses, as well as in DRG assignment. RESULTS: There were 397 discharges of cases of CVD which included 21 different DRGs. Diagnostic coding showed a concordance rate of 81.87%; the selection error rate was 2.26% and the classification error rate was 15.87%. The error rate in the DRG was 16.12% and associated with the greatest impact on the mortality risk level. CONCLUSIONS: While the errors we observed must be taken into account, data suggest that the quality of the MBDS for CVD is sufficient to ensure delivery of valid information. The hospital discharge registry serves as a valuable tool for use in studies of this disease.