Epigenomics and genotype-phenotype association analyses reveal conserved genetic architecture of complex traits in cattle and human.

BACKGROUND: Lack of comprehensive functional annotations across a wide range of tissues and cell types severely hinders the biological interpretations of phenotypic variation, adaptive evolution, and domestication in livestock. Here we used a combination of comparative epigenomics, genome-wide association study (GWAS), and selection signature analysis, to shed light on potential adaptive evolution in cattle. RESULTS: We cross-mapped 8 histone marks of 1300 samples from human to cattle, covering 178 unique tissues/cell types. By uniformly analyzing 723 RNA-seq and 40 whole genome bisulfite sequencing (WGBS) datasets in cattle, we validated that cross-mapped histone marks captured tissue-specific expression and methylation, reflecting tissue-relevant biology. Through integrating cross-mapped tissue-specific histone marks with large-scale GWAS and selection signature results, we for the first time detected relevant tissues and cell types for 45 economically important traits and artificial selection in cattle. For instance, immune tissues are significantly associated with health and reproduction traits, multiple tissues for milk production and body conformation traits (reflecting their highly polygenic architecture), and thyroid for the different selection between beef and dairy cattle. Similarly, we detected relevant tissues for 58 complex traits and diseases in humans and observed that immune and fertility traits in humans significantly correlated with those in cattle in terms of relevant tissues, which facilitated the identification of causal genes for such traits. For instance, PIK3CG, a gene highly specifically expressed in mononuclear cells, was significantly associated with both age-at-menopause in human and daughter-still-birth in cattle. ICAM, a T cell-specific gene, was significantly associated with both allergic diseases in human and metritis in cattle. CONCLUSION: Collectively, our results highlighted that comparative epigenomics in conjunction with GWAS and selection signature analyses could provide biological insights into the phenotypic variation and adaptive evolution. Cattle may serve as a model for human complex traits, by providing additional information beyond laboratory model organisms, particularly when more novel phenotypes become available in the near future.

Longitudinal study of humoral immunity to bovine coronavirus, virus shedding, and treatment for bovine respiratory disease in pre-weaned beef calves.

BACKGROUND: Bovine coronavirus (BCV) is associated with respiratory infections in cattle of all ages; however, a temporal study to evaluate the effect of BCV immunity on virus shedding and bovine respiratory disease (BRD) incidence in pre-weaned beef calves has not been reported. Thus, we report here a prospective study in three herds of crossbred beef calves (n = 817) with endemic BCV. Serial blood samples for measurement of serum anti-BCV antibody titers and nasal swabs for detection of BCV and other common viral and bacterial BRD pathogens were collected from all calves or subsets of calves at predetermined times from birth through weaning. The calves were monitored for BRD and those that developed signs of respiratory disease were sampled for diagnostic testing. To discover additional risk factors that could have influenced BRD development, sequence analysis of the BCV strain(s) circulating in each herd, and the prevalence of common opportunistic bacterial pathogens in the upper respiratory tract of sick and apparently healthy cattle were also evaluated. RESULTS: Two hundred forty-eight of the 817 study calves (30.4%) were treated for BRD prior to weaning; 246 of those were from a single herd involved in two outbreaks of BRD leading to mass treatment of all calves in that group. Molecular diagnostic testing found BCV and Histophilus somni in nasal swabs taken at the time of BRD treatment. Between herd analyses revealed anti-BCV serum antibody abundance did not associate with the incidence of BRD or BCV shedding, though these measurements may have been hindered by the long periods between sample collections. Analysis of the BCV spike gene hypervariable region revealed four polymorphisms in 15 isolates from the three herds, making strain variation unlikely to account for differences in treatment rates between herds. Persistent or recurrent shedding episodes of BCV occurred in some animals treated for BRD. CONCLUSION: Co-detection of BCV and H. somni at the time of the disease outbreak suggests that these pathogens contributed to disease pathogenesis. Developing appropriate control measures for respiratory BCV infections may help decrease the incidence of pre-weaning BRD. The role of antibodies in protection must still be further defined.

Association of Circulating Transfer RNA fragments with antibody response to Mycoplasma bovis in beef cattle.

BACKGROUND: High throughput sequencing allows identification of small non-coding RNAs. Transfer RNA Fragments are a class of small non-coding RNAs, and have been identified as being involved in inhibition of gene expression. Given their role, it is possible they may be involved in mediating the infection-induced defense response in the host. Therefore, the objective of this study was to identify 5' transfer RNA fragments (tRF5s) associated with a serum antibody response to M. bovis in beef cattle. RESULTS: The tRF5s encoding alanine, glutamic acid, glycine, lysine, proline, selenocysteine, threonine, and valine were associated (P < 0.05) with antibody response against M. bovis. tRF5s encoding alanine, glutamine, glutamic acid, glycine, histidine, lysine, proline, selenocysteine, threonine, and valine were associated (P < 0.05) with season, which could be attributed to calf growth. There were interactions (P < 0.05) between antibody response to M. bovis and season for tRF5 encoding selenocysteine (anticodon UGA), proline (anticodon CGG), and glutamine (anticodon TTG). Selenocysteine is a rarely used amino acid that is incorporated into proteins by the opal stop codon (UGA), and its function is not well understood. CONCLUSIONS: Differential expression of tRF5s was identified between ELISA-positive and negative animals. Production of tRF5s may be associated with a host defense mechanism triggered by bacterial infection, or it may provide some advantage to a pathogen during infection of a host. Further studies are needed to establish if tRF5s could be used as a diagnostic marker of chronic exposure.

Incorporation of genetic technologies associated with applied reproductive technologies to enhance world food production.

Animal breeding and reproductive physiology have been closely related throughout the history of animal production science. Artificial insemination provides the best method of increasing the influence of sires with superior genetics to improve production traits. Multiple ovulation embryo transfer (MOET) provides some ability to increase the genetic influence of the maternal line as well. The addition of genetic technologies to this paradigm allows for improved methods of selecting sires and dams carrying the best genes for production and yield of edible products and resistance to diseases and parasites. However, decreasing the number of influential parents within a population also increases the risk of propagating a recessive gene that could negatively impact the species (Reprod Domest Anim 44:792-796, 2009; BMC Genomics 11:337, 2010). Furthermore, antagonistic genotypic relationships between production traits and fertility (Anim Prod Sci 49:399-412, 2009; Anim Genet 43:442-446, 2012) suggest that care must be taken to ensure that increasing the frequency of genes with a positive influence on production does not negatively impact the fertility of the replacement females entering the herd.

Phylogeny of Shiga toxin-producing Escherichia coli O157 isolated from cattle and clinically ill humans.

Cattle are a major reservoir for Shiga toxin-producing Escherichia coli O157 (STEC O157) and harbor multiple genetic subtypes that do not all associate with human disease. STEC O157 evolved from an E. coli O55:H7 progenitor; however, a lack of genome sequence has hindered investigations on the divergence of human- and/or cattle-associated subtypes. Our goals were to 1) identify nucleotide polymorphisms for STEC O157 genetic subtype detection, 2) determine the phylogeny of STEC O157 genetic subtypes using polymorphism-derived genotypes and a phage insertion typing system, and 3) compare polymorphism-derived genotypes identified in this study with pulsed field gel electrophoresis (PFGE), the current gold standard for evaluating STEC O157 diversity. Using 762 nucleotide polymorphisms that were originally identified through whole-genome sequencing of 189 STEC O157 human- and cattle-isolated strains, we genotyped a collection of 426 STEC O157 strains. Concatenated polymorphism alleles defined 175 genotypes that were tagged by a minimal set of 138 polymorphisms. Eight major lineages of STEC O157 were identified, of which cattle are a reservoir for seven. Two lineages regularly harbored by cattle accounted for the majority of human disease in this study, whereas another was rarely represented in humans and may have evolved toward reduced human virulence. Notably, cattle are not a known reservoir for E. coli O55:H7 or STEC O157:H(-) (the first lineage to diverge within the STEC O157 serogroup), which both cause human disease. This result calls into question how cattle may have originally acquired STEC O157. The polymorphism-derived genotypes identified in this study did not surpass PFGE diversity assessed by BlnI and XbaI digestions in a subset of 93 strains. However, our results show that they are highly effective in assessing the evolutionary relatedness of epidemiologically unrelated STEC O157 genetic subtypes, including those associated with the cattle reservoir and human disease.

Analysis of copy number variations among diverse cattle breeds.

Genomic structural variation is an important and abundant source of genetic and phenotypic variation. Here, we describe the first systematic and genome-wide analysis of copy number variations (CNVs) in modern domesticated cattle using array comparative genomic hybridization (array CGH), quantitative PCR (qPCR), and fluorescent in situ hybridization (FISH). The array CGH panel included 90 animals from 11 Bos taurus, three Bos indicus, and three composite breeds for beef, dairy, or dual purpose. We identified over 200 candidate CNV regions (CNVRs) in total and 177 within known chromosomes, which harbor or are adjacent to gains or losses. These 177 high-confidence CNVRs cover 28.1 megabases or approximately 1.07% of the genome. Over 50% of the CNVRs (89/177) were found in multiple animals or breeds and analysis revealed breed-specific frequency differences and reflected aspects of the known ancestry of these cattle breeds. Selected CNVs were further validated by independent methods using qPCR and FISH. Approximately 67% of the CNVRs (119/177) completely or partially span cattle genes and 61% of the CNVRs (108/177) directly overlap with segmental duplications. The CNVRs span about 400 annotated cattle genes that are significantly enriched for specific biological functions, such as immunity, lactation, reproduction, and rumination. Multiple gene families, including ULBP, have gone through ruminant lineage-specific gene amplification. We detected and confirmed marked differences in their CNV frequencies across diverse breeds, indicating that some cattle CNVs are likely to arise independently in breeds and contribute to breed differences. Our results provide a valuable resource beyond microsatellites and single nucleotide polymorphisms to explore the full dimension of genetic variability for future cattle genomic research.

Use of overlapping DNA pools to discern genetic differences despite pooling error.

Genotyping pools of commercial cattle and individual seedstock animals may reveal hidden relationships between sectors enabling use of commercial data for genetic evaluation. However, commercial data capture may be compromised by inexact pool formation. We aimed to estimate the concordance between distances or genomic covariance among pooling allele frequencies (PAFs) of DNA pools comprised of 100 animals with 0% or 50% overlap of animals in common between pools. Cattle lung samples were collected from a commercial beef processing plant on a single day. Six pools of 100 animals each were constructed so that overlap between pools was 0% or 50%. Two pools of all 200 animals were constructed to estimate PAFs for all 200 animals. Frozen lung tissue (0.01 g) from each animal was weighed into a tube containing a pool; there were two pools of 200 animals each and six pools of 100 animals each. Every contribution of an individual animal was an independent measurement to insure independence of pooling errors. Lung samples were kept on dried ice during the pooling process to keep them from thawing. The eight pools were then assayed for approximately 100,000 single nucleotide polymorphisms (SNP). PAF for each SNP and pool was based on the relative intensity of the two dyes used to detect the alleles rather than genotype calls which are not tractable from pooling data. Euclidean distances and genomic relationships among the PAFs for the eight pools were estimated and distances were tested for concordance with pool overlap using permutation-based analysis of distance. Distances among pools were concordant with the planned overlap of animals shared between pools (P = 0.0024); pool overlap accounted for 70% of the variation and pooling error accounted for 30%. Pools containing 100 animals with no overlap were the most distant from one another and pools with 50% overlap were the least distant. This work shows that we can discern differences in distance between pairs of overlapping DNA pools sharing 0% and 50% of the animals. Genomic correlations among nonoverlapping pools indicated that nonoverlapping pool pairs did not share many related animals because genomic correlations were near zero for these pairs. On the other hand, one pair of nonoverlapping pools likely contained related animals between pools because the correlation was 0.21. Pools sharing 50% overlap ranged in genomic relationship between 0.21 and 0.39 (N = 12).

Genetic evaluation of seedstock cattle could benefit from commercial data. There are hidden relationships between commercial and seedstock sectors because many commercial producers buy bulls from the seedstock sector. Relationships are hidden because pedigree is not tracked in commercial populations. Single nucleotide polymorphism genotypes could reveal these hidden relationships; however, genotyping can be cost prohibitive. Cost of commercial data capture could be decreased by pooling DNA which is a method to genotype groups of animals to use their data in genetic evaluation; however, error from inexact pool formation can complicate interpretation. Results from pools of overlapping random unrelated animals mimic the results from pools sharing relatives with the same degree of shared genomes. For example, a pool of progeny and a pool of the dams of the pooled progeny would produce the same result as two pools sharing 50% overlap of random unrelated animals. We can estimate the relatedness between unknown pools even in the presence of pooling error if an unknown pool comparison is similar to an overlapping pool comparison. Knowing the relationship between seedstock cattle and pools of commercial cattle may allow commercial data to enhance genetic evaluation of seedstock animals.

Probiotics in milk replacer affect the microbiome of the lung in neonatal dairy calves.

Introduction: Probiotics have been investigated for their many health benefits and impact on the microbiota of the gut. Recent data have also supported a gut-lung axis regarding the bacterial populations (microbiomes) of the two locations; however, little research has been performed to determine the effects of oral probiotics on the microbiome of the bovine respiratory tract. We hypothesized that probiotic treatment would result in changes in the lung microbiome as measured in lung lavage fluid. Our overall goal was to characterize bacterial populations in the lungs of calves fed probiotics in milk replacer and dry rations from birth to weaning. Methods: A group of 20 dairy calves was split into two treatment groups: probiotic (TRT; N = 10, milk replacer +5 g/d probiotics; Bovamine Dairy, Chr. Hansen, Inc., Milwaukee, WI) and control (CON; N = 10, milk replacer only). On day 0, birth weight was obtained, and calves were provided colostrum as per the dairy SOP. On day 2, probiotics were added to the milk replacer of the treated group and then included in their dry ration. Lung lavages were performed on day 52 on five random calves selected from each treatment group. DNA was extracted from lavage fluid, and 16S ribosomal RNA (rRNA) gene hypervariable regions 1-3 were amplified by PCR and sequenced using next-generation sequencing (Illumina MiSeq) for the identification of the bacterial taxa present. Taxa were classified into both operational taxonomic units (OTUs) and amplicon sequence variants (ASVs). Results: Overall, the evaluation of these samples revealed that the bacterial genera identified in the lung lavage samples of probiotic-fed calves as compared to the control calves were significantly different based on the OTU dataset (p < 0.05) and approached significance for the ASV dataset (p < 0.06). Additionally, when comparing the diversity of taxa in lung lavage samples to nasal and tonsil samples, taxa diversity of lung samples was significantly lower (p < 0.05). Discussion: In conclusion, analysis of the respiratory microbiome in lung lavage samples after probiotic treatment provides insight into the distribution of bacterial populations in response to oral probiotics and demonstrates that oral probiotics affect more than the gut microbiome.

Polymorphism of the follicle stimulating hormone receptor does not impact reproductive performance or in-vitro embryo production in beef heifers.

Assisted reproductive technologies are used to propagate desirable genetics in a shortened timeframe. Selected females undergo ovarian stimulation with the use of follicle stimulating hormone (FSH) to increase embryo recovery for subsequent transfer programs. The FSH receptor (FSHR) c.337 C > G variant was reported to have a reduction in viable embryo numbers in an ovarian stimulation protocol. We, therefore, hypothesized that FSHR c.337 C > G would result in reduced in-vitro blastocyst development. Beef heifers were genotyped and selected based on the c.337 C > G FSHR genotype (CC, CG, GG; n = 15-16/genotype). Estrus was synchronized with a Select Synch protocol and heifers were slaughtered 5 days after induced ovulation. Anterior pituitaries, serum and reproductive tracts were collected at slaughter for analysis. Cumulus oocyte complexes (COCs) were collected and pooled within genotype for in-vitro fertilization (IVF) and subsequent blastocyst development. No differences were observed in carcass weights, anterior pituitary weights, serum progesterone, corpus lutea weight, surface follicle counts, histological follicle counts or follicular fluid estradiol concentration (P > 0.1) due to FSHR genotype. Differences were observed for ovulation rates in the GG FSHR genotype group (P < 0.01). However, cleavage and blastocyst rates were not affected due to FSHR genotype (P > 0.1), following standard IVF protocols. The FSHR variant does not influence antral follicle counts, estradiol production, or in-vitro blastocyst development in beef heifers. The GG FSHR genotype had an increased ovulation rate, which may indicate a greater potential for twinning, but research with a larger population is warranted to support this hypothesis.

Gene number determination and genetic polymorphism of the gamma delta T cell co-receptor WC1 genes.

BACKGROUND: WC1 co-receptors belong to the scavenger receptor cysteine-rich (SRCR) superfamily and are encoded by a multi-gene family. Expression of particular WC1 genes defines functional subpopulations of WC1(+) γδ T cells. We have previously identified partial or complete genomic sequences for thirteen different WC1 genes through annotation of the bovine genome Btau_3.1 build. We also identified two WC1 cDNA sequences from other cattle that did not correspond to sequences in the Btau_3.1 build. Their absence in the Btau_3.1 build may have reflected gaps in the genome assembly or polymorphisms among animals. Since the response of γδ T cells to bacterial challenge is determined by WC1 gene expression, it was critical to understand whether individual cattle or breeds differ in the number of WC1 genes or display polymorphisms. RESULTS: Real-time quantitative PCR using DNA from the animal whose genome was sequenced ("Dominette") and sixteen other animals representing ten breeds of cattle, showed that the number of genes coding for WC1 co-receptors is thirteen. The complete coding sequences of those thirteen WC1 genes is presented, including the correction of an error in the WC1-2 gene due to mis-assembly in the Btau_3.1 build. All other cDNA sequences were found to agree with the previous annotation of complete or partial WC1 genes. PCR amplification and sequencing of the most variable N-terminal SRCR domain (domain 1 which has the SRCR "a" pattern) of each of the thirteen WC1 genes showed that the sequences are highly conserved among individuals and breeds. Of 160 sequences of domain 1 from three breeds of cattle, no additional sequences beyond the thirteen described WC1 genes were found. Analysis of the complete WC1 cDNA sequences indicated that the thirteen WC1 genes code for three distinct WC1 molecular forms. CONCLUSION: The bovine WC1 multi-gene family is composed of thirteen genes coding for three structural forms whose sequences are highly conserved among individual cattle and breeds. The sequence diversity necessary for WC1 genes to function as a multi-genic pattern recognition receptor array is encoded in the genome, rather than generated by recombinatorial diversity or hypermutation.

Recent Emergence of Bovine Coronavirus Variants with Mutations in the Hemagglutinin-Esterase Receptor Binding Domain in U.S. Cattle.

Bovine coronavirus (BCoV) has spilled over to many species, including humans, where the host range variant coronavirus OC43 is endemic. The balance of the opposing activities of the surface spike (S) and hemagglutinin-esterase (HE) glycoproteins controls BCoV avidity, which is critical for interspecies transmission and host adaptation. Here, 78 genomes were sequenced directly from clinical samples collected between 2013 and 2022 from cattle in 12 states, primarily in the Midwestern U.S. Relatively little genetic diversity was observed, with genomes having >98% nucleotide identity. Eleven isolates collected between 2020 and 2022 from four states (Nebraska, Colorado, California, and Wisconsin) contained a 12 nucleotide insertion in the receptor-binding domain (RBD) of the HE gene similar to one recently reported in China, and a single genome from Nebraska collected in 2020 contained a novel 12 nucleotide deletion in the HE gene RBD. Isogenic HE proteins containing either the insertion or deletion in the HE RBD maintained esterase activity and could bind bovine submaxillary mucin, a substrate enriched in the receptor 9-O-acetylated-sialic acid, despite modeling that predicted structural changes in the HE R3 loop critical for receptor binding. The emergence of BCoV with structural variants in the RBD raises the possibility of further interspecies transmission.

Evaluating Accuracy of DNA Pool Construction Based on White Blood Cell Counts.

Pooling individual samples prior to DNA extraction can mitigate the cost of DNA extraction and genotyping; however, these methods need to accurately generate equal representation of individuals within pools. The objective of this study was to determine accuracy of pool construction of blood samples based on white blood cell counts compared to two common DNA quantification methods. Fifty individual bovine blood samples were collected, and then pooled with all individuals represented in each pool. Pools were constructed with the target of equal representation of each individual animal based on number of white blood cells, spectrophotometric readings, spectrofluorometric readings, and whole blood volume with 9 pools per method and a total of 36 pools. Pools and individual samples that comprised the pools were genotyped using a commercially available genotyping array. ASReml was used to estimate variance components for individual animal contribution to pools. The correlation between animal contributions between two pools was estimated using bivariate analysis with starting values set to the result of a univariate analysis. Adonis test on distance matrix from the animal correlation showed clustering with method, and higher correlations between methods than within (P < 1 × 10-6). White blood cell count was predictive of sample representation when compared to pooling based on DNA concentration. Therefore, constructing pools using white blood cell counts prior to DNA extraction may reduce cost associated with DNA extraction and genotyping and improve representation of individuals in a pool.

A deletion mutation in bovine SLC4A2 is associated with osteopetrosis in Red Angus cattle.

BACKGROUND: Osteopetrosis is a skeletal disorder of humans and animals characterized by the formation of overly dense bones, resulting from a deficiency in the number and/or function of bone-resorbing osteoclast cells. In cattle, osteopetrosis can either be induced during gestation by viral infection of the dam, or inherited as a recessive defect. Genetically affected calves are typically aborted late in gestation, display skull deformities and exhibit a marked reduction of osteoclasts. Although mutations in several genes are associated with osteopetrosis in humans and mice, the genetic basis of the cattle disorder was previously unknown. RESULTS: We have conducted a whole-genome association analysis to identify the mutation responsible for inherited osteopetrosis in Red Angus cattle. Analysis of >54,000 SNP genotypes for each of seven affected calves and nine control animals localized the defective gene to the telomeric end of bovine chromosome 4 (BTA4). Homozygosity analysis refined the interval to a 3.4-Mb region containing the SLC4A2 gene, encoding an anion exchanger protein necessary for proper osteoclast function. Examination of SLC4A2 from normal and affected animals revealed a approximately 2.8-kb deletion mutation in affected calves that encompasses exon 2 and nearly half of exon 3, predicted to prevent normal protein function. Analysis of RNA from a proven heterozygous individual confirmed the presence of transcripts lacking exons 2 and 3, in addition to normal transcripts. Genotyping of additional animals demonstrated complete concordance of the homozygous deletion genotype with the osteopetrosis phenotype. Histological examination of affected tissues revealed scarce, morphologically abnormal osteoclasts displaying evidence of apoptosis. CONCLUSIONS: These results indicate that a deletion mutation within bovine SLC4A2 is associated with osteopetrosis in Red Angus cattle. Loss of SLC4A2 function appears to induce premature cell death, and likely results in cytoplasmic alkalinization of osteoclasts which, in turn, may disrupt acidification of resorption lacunae.

Differentiation of Mannheimia haemolytica genotype 1 and 2 strains by visible phenotypic characteristics on solid media.

Genotype 2 Mannheimia haemolytica associate with the lungs of cattle with bovine respiratory disease more frequently than genotype 1 strains. Different colony colors and morphologies were identified between genotype 1 and 2 solid media cultures. Genotype of strains, and frequency differences between them in mixed cultures are discernable by visual inspection.

Classification of 16S rRNA reads is improved using a niche-specific database constructed by near-full length sequencing.

Surveys of microbial populations in environmental niches of interest often utilize sequence variation in the gene encoding the ribosomal small subunit (the 16S rRNA gene). Generally, these surveys target the 16S genes using semi-degenerate primers to amplify portions of a subset of bacterial species, sequence the amplicons in bulk, and assign to putative taxonomic categories by comparison to databases purporting to connect specific sequences in the main variable regions of the gene to specific organisms. Due to sequence length constraints of the most popular bulk sequencing platforms, the primers selected amplify one to three of the nine variable regions, and taxonomic assignment is based on relatively short stretches of sequence (150-500 bases). We demonstrate that taxonomic assignment is improved through reduced unassigned reads by including a survey of near-full-length sequences specific to the target environment, using a niche of interest represented by the upper respiratory tract (URT) of cattle. We created a custom Bovine URT database from these longer sequences for assignment of shorter, less expensive reads in comparisons of the upper respiratory tract among individual animals. This process improves the ability to detect changes in the microbial populations of a given environment, and the accuracy of defining the content of that environment at increasingly higher taxonomic resolution.

De novo assembly of the cattle reference genome with single-molecule sequencing.

BACKGROUND: Major advances in selection progress for cattle have been made following the introduction of genomic tools over the past 10-12 years. These tools depend upon the Bos taurus reference genome (UMD3.1.1), which was created using now-outdated technologies and is hindered by a variety of deficiencies and inaccuracies. RESULTS: We present the new reference genome for cattle, ARS-UCD1.2, based on the same animal as the original to facilitate transfer and interpretation of results obtained from the earlier version, but applying a combination of modern technologies in a de novo assembly to increase continuity, accuracy, and completeness. The assembly includes 2.7 Gb and is >250× more continuous than the original assembly, with contig N50 >25 Mb and L50 of 32. We also greatly expanded supporting RNA-based data for annotation that identifies 30,396 total genes (21,039 protein coding). The new reference assembly is accessible in annotated form for public use. CONCLUSIONS: We demonstrate that improved continuity of assembled sequence warrants the adoption of ARS-UCD1.2 as the new cattle reference genome and that increased assembly accuracy will benefit future research on this species.

MicroRNA transcriptome profiles during swine skeletal muscle development.

BACKGROUND: MicroRNA (miR) are a class of small RNAs that regulate gene expression by inhibiting translation of protein encoding transcripts. To evaluate the role of miR in skeletal muscle of swine, global microRNA abundance was measured at specific developmental stages including proliferating satellite cells, three stages of fetal growth, day-old neonate, and the adult. RESULTS: Twelve potential novel miR were detected that did not match previously reported sequences. In addition, a number of miR previously reported to be expressed in mammalian muscle were detected, having a variety of abundance patterns through muscle development. Muscle-specific miR-206 was nearly absent in proliferating satellite cells in culture, but was the highest abundant miR at other time points evaluated. In addition, miR-1 was moderately abundant throughout developmental stages with highest abundance in the adult. In contrast, miR-133 was moderately abundant in adult muscle and either not detectable or lowly abundant throughout fetal and neonate development. Changes in abundance of ubiquitously expressed miR were also observed. MiR-432 abundance was highest at the earliest stage of fetal development tested (60 day-old fetus) and decreased throughout development to the adult. Conversely, miR-24 and miR-27 exhibited greatest abundance in proliferating satellite cells and the adult, while abundance of miR-368, miR-376, and miR-423-5p was greatest in the neonate. CONCLUSION: These data present a complete set of transcriptome profiles to evaluate miR abundance at specific stages of skeletal muscle growth in swine. Identification of these miR provides an initial group of miR that may play a vital role in muscle development and growth.

Microbiome of the upper nasal cavity of beef calves prior to weaning12.

Disease incidence is intimately associated with an animal's commensal bacteria populations (microbiome), as microbes that are involved with morbidity and mortality are commonly found in animals with no sign of disease. An understanding of the animal's resident respiratory pathogens, in the upper nasal cavity prior to weaning, may help us to understand the impact of these pathogens on incidence of respiratory disease. For this research, the overall goal was to characterize bacterial populations associated with calves at an early age and through time periods prior to weaning in 3 herds at the U.S. Meat Animal Research Center. Nasal swabs from the upper nasal cavity were collected at initial vaccination (approximately 40 d of age), preconditioning (approximately 130 d of age), and weaning (approximately 150 d of age) in 2015 and 2016. DNA was extracted from nasal swabs and combined into 2 pools of 10 animals for each sampling time point, in each herd, for a total of 6 pools at each sampling time point and 18 pools for all sampling time points within each year. To evaluate and compare the microbiome of each pooled sample, hypervariable regions 1 through 3 along the 16S ribosomal RNA (rRNA) gene were amplified by PCR and sequenced using next-generation sequencing (Illumina MiSeq) for identification of the bacterial taxa present. Alpha and beta diversity were also measured. Overall, microbial communities were different between combinations of sampling year, herd location, and sampling time prior to weaning as shown by beta diversity. Analysis of these specific respiratory pathogens prior to weaning will present a clearer picture of the distribution of microbial populations in animals prior to weaning and not exhibiting clinical signs of respiratory disease. Therefore, evaluation of the animal's resident bacterial populations in the upper nasal cavity during different phases of the beef production system may help us to understand the impact of the microbiome on incidence of respiratory disease in cattle.

Evaluating the microbiome of two sampling locations in the nasal cavity of cattle with bovine respiratory disease complex (BRDC).

Bovine respiratory disease complex (BRDC) is a multifactor disease, and disease incidence may be associated with an animal's commensal bacterial populations (microbiome) in the upper nasal cavity. Identifying these commensal bacterial populations in the upper nasal cavity may help us to understand the impact of the microbiome on incidence of BRDC in cattle. Various sampling techniques have previously been utilized to evaluate the microbiome of different locations of the upper nasal cavity in cattle. Therefore, our objective was to determine whether bacterial populations of the nasal cavity vary based on these sampling locations. Two common sampling techniques were evaluated, including 6-inch nasal swabs and deep nasopharyngeal swabs. Nasal swabs from calves were collected when the animal was diagnosed with BRDC after weaning in the feedlot in addition to collection of samples from asymptomatic cohorts. Samples were pooled in groups based on year the animal was in the feedlot (2015 or 2016), when the animal was diagnosed with BRDC (1 to 5 weeks after weaning), type of sample (6-inch nasal swab or deep nasopharyngeal swab), and health status (diagnosis with BRDC or control). Variable regions 1 through 3 along the 16S rRNA gene were amplified by PCR and sequenced using next-generation sequencing (Illumina MiSeq) for identification of the bacterial taxa present. Overall, sampling site did not consistently influence diversity of the bacterial populations of the upper nasal cavity. However, the effect of disease incidence on the microbiome was depended on sampling time after weaning (P = 0.0462) for 2015, while the main effects of sampling time after weaning (P = 0.00992) and disease phenotype (P = 0.012) were significant for 2016. These data for 2016 demonstrate that in addition to bacterial profiles changing throughout weaning, calves diagnosed with BRDC have different bacterial profiles compared to their control cohorts. In addition, evaluation of the microbiome identified predominant bacteria genera in the upper nasal cavity included those previously reported to be associated with cattle diagnosed with BRDC including Mycoplasma sp., Psychrobacter sp., and Mannheimia sp. In summary, these results demonstrate that shorter, less invasive 6-inch nasal swabs produce similar results to deep nasopharyngeal swabs.

Complete blood count data and leukocyte expression of cytokine genes and cytokine receptor genes associated with bovine respiratory disease in calves.

OBJECTIVE: The purpose of this study was to evaluate potential relationships between cytokine gene expression, complete blood counts (CBC) and animals that were sick or would become sick. The CBC and the transcript abundance of cytokines and their receptors expressed in leukocytes were measured from calves at two early timepoints, and again after diagnosis with bovine respiratory disease (BRD). RESULTS: Blood was collected from calves at pre-conditioning (n = 796) and weaning (n = 791) for CBC. Blood counts were also measured for the calves with BRD (n = 13), and asymptomatic calves (n = 75) after weaning. The CBC were compared for these animals at 3 time points. At diagnosis, neutrophils were higher and basophils lower in sick animals (P < 0.05). To further characterize BRD responses, transcript abundance of 84 cytokine genes were evaluated in 5 calves with BRD and 9 asymptomatic animals at all time points. There was more data for CBC than transcript abundance; hence, animal and temporary environmental correlations between CBC and transcript abundance were exploited to improve the power of the transcript abundance data. Expression of CCL16, CXCR1, CCR1 was increased in BRD positive animals compared to controls (P-corrected < 0.1). Cytokine expression data may help to provide insight into an animal's health.