The first complete T2T Assemblies of Cattle and Sheep Y-Chromosomes uncover remarkable divergence in structure and gene content.

Reference genomes of cattle and sheep have lacked contiguous assemblies of the sex-determining Y chromosome. We assembled complete and gapless telomere to telomere (T2T) Y chromosomes for these species. The pseudo-autosomal regions were similar in length, but the total chromosome size was substantially different, with the cattle Y more than twice the length of the sheep Y. The length disparity was accounted for by expanded ampliconic region in cattle. The genic amplification in cattle contrasts with pseudogenization in sheep suggesting opposite evolutionary mechanisms since their divergence 18MYA. The centromeres also differed dramatically despite the close relationship between these species at the overall genome sequence level. These Y chromosome have been added to the current reference assemblies in GenBank opening new opportunities for the study of evolution and variation while supporting efforts to improve sustainability in these important livestock species that generally use sire-driven genetic improvement strategies.

Enhanced bovine genome annotation through integration of transcriptomics and epi-transcriptomics datasets facilitates genomic biology.

BACKGROUND: The accurate identification of the functional elements in the bovine genome is a fundamental requirement for high-quality analysis of data informing both genome biology and genomic selection. Functional annotation of the bovine genome was performed to identify a more complete catalog of transcript isoforms across bovine tissues. RESULTS: A total of 160,820 unique transcripts (50% protein coding) representing 34,882 unique genes (60% protein coding) were identified across tissues. Among them, 118,563 transcripts (73% of the total) were structurally validated by independent datasets (PacBio isoform sequencing data, Oxford Nanopore Technologies sequencing data, de novo assembled transcripts from RNA sequencing data) and comparison with Ensembl and NCBI gene sets. In addition, all transcripts were supported by extensive data from different technologies such as whole transcriptome termini site sequencing, RNA Annotation and Mapping of Promoters for the Analysis of Gene Expression, chromatin immunoprecipitation sequencing, and assay for transposase-accessible chromatin using sequencing. A large proportion of identified transcripts (69%) were unannotated, of which 86% were produced by annotated genes and 14% by unannotated genes. A median of two 5' untranslated regions were expressed per gene. Around 50% of protein-coding genes in each tissue were bifunctional and transcribed both coding and noncoding isoforms. Furthermore, we identified 3,744 genes that functioned as noncoding genes in fetal tissues but as protein-coding genes in adult tissues. Our new bovine genome annotation extended more than 11,000 annotated gene borders compared to Ensembl or NCBI annotations. The resulting bovine transcriptome was integrated with publicly available quantitative trait loci data to study tissue-tissue interconnection involved in different traits and construct the first bovine trait similarity network. CONCLUSIONS: These validated results show significant improvement over current bovine genome annotations.

Single-cell insights into development of the bovine placenta†.

A central determinant of pregnancy success is proper development of the conceptus (embryo/fetus and associated extraembryonic membranes including the placenta). Although the gross morphology and histology of the bovine placenta have been well studied, the cellular and molecular mechanisms regulating placenta development and trophoblast differentiation and function remain essentially undefined. Here, single-cell transcriptome (scRNA-seq) analysis was performed on the day 17 bovine conceptus and chorion of day 24, 30, and 50 conceptuses (n = 3-4 samples per day) using the 10X Genomics platform. Bioinformatic analyses identified cell types and their ontogeny including trophoblast, mesenchyme, and immune cells. Loss of interferon tau-expressing trophoblast uninucleate cells occurred between days 17 and 30, whereas binucleate cells, identified based on expression of placental lactogen (CSH2) and specific pregnancy-associated glycoprotein genes (PAGs), first appeared on day 24. Several different types of uninucleate cells were present in day 24, 30, and 50 samples, but only one (day 24) or two types of binucleate cells (days 30 and 50). Cell trajectory analyses provided a conceptual framework for uninucleate cell development and binucleate cell differentiation, and bioinformatic analyses identified candidate transcription factors governing differentiation and function of the trophoblasts. The digital atlas of cell types in the developing bovine conceptus reported here serves as a resource to discover key genes and biological pathways regulating its development during the critical periods of implantation and placentation.

Single-cell insights into epithelial morphogenesis in the neonatal mouse uterus.

The uterus is vital for successful reproduction in mammals, and two different types of epithelia (luminal and glandular) are essential for embryo implantation and pregnancy establishment. However, the essential cellular and molecular factors and pathways governing postnatal epithelium maturation, determination, and differentiation in developing uterus are yet to be elucidated. Here, the epithelium of the neonatal mouse uterus was isolated and subjected to single-cell transcriptome (scRNA-seq) analysis. Both the undifferentiated epithelium and determined luminal epithelium were heterogeneous and contained several different cell clusters based on single-cell transcription profiles. Substantial gene expression differences were evident as the epithelium matured and differentiated between postnatal days 1 to 15. Two new glandular epithelium-expressed genes (Gas6 and Cited4) were identified and validated by in situ hybridization. Trajectory analyses provided a framework for understanding epithelium maturation, lineage bifurcation, and differentiation. A candidate set of transcription factors and gene regulatory networks were identified that potentially direct epithelium lineage specification and morphogenesis. This atlas provides a foundation important to discover intrinsic cellular and molecular mechanisms directing uterine epithelium morphogenesis during a critical window of postnatal development.

Advanced Skeletal Ossification Is Associated with Genetic Variants in Chronologically Young Beef Heifers.

Osteogenesis is a developmental process critical for structural support and the establishment of a dynamic reservoir for calcium and phosphorus. Changes in livestock breeding over the past 100 years have resulted in earlier bone development and increased physical size of cattle. Advanced skeletal maturity is now commonly observed at harvest, with heifers displaying more mature bone than is expected at 30 months of age (MOA). We surmise that selection for growth traits and earlier reproductive maturity resulted in co-selection for accelerated skeletal ossification. This study examines the relationship of single nucleotide polymorphisms (SNPs) in 793 beef heifers under 30 MOA with USDA-graded skeletal maturity phenotypes (A-, B-, C- skeletal maturity). Further, the estrogen content of FDA-approved hormonal implants provided to heifers prior to harvest was evaluated in association with the identified SNPs and maturities. Association tests were performed, and the impact of the implants were evaluated as covariates against genotypes using a logistic regression model. SNPs from the ESR1, ALPL, PPARGC1B, SORCS1 genes, and SNPs near KLF14, ANKRD61, USP42, H1C1, OVCA2, microRNA mir-29a were determined to be associated with the advanced skeletal ossification phenotype in heifers. Higher dosage estrogen implants increased skeletal maturity in heifers with certain SNP genotypes.

Estrogen receptor alpha regulates uterine epithelial lineage specification and homeostasis.

Postnatal development of the uterus involves specification of undifferentiated epithelium into uterine-type epithelium. That specification is regulated by stromal-epithelial interactions as well as intrinsic cell-specific transcription factors and gene regulatory networks. This study utilized mouse genetic models of Esr1 deletion, endometrial epithelial organoids (EEO), and organoid-stromal co-cultures to decipher the role of Esr1 in uterine epithelial development. Organoids derived from wild-type (WT) mice developed a normal single layer of columnar epithelium. In contrast, EEO from Esr1 null mice developed a multilayered stratified squamous type of epithelium with basal cells. Co-culturing Esr1 null epithelium with WT uterine stromal fibroblasts inhibited basal cell development. Of note, estrogen treatment of EEO-stromal co-cultures and Esr1 conditional knockout mice increased basal epithelial cell markers. Collectively, these findings suggest that Esr1 regulates uterine epithelium lineage plasticity and homeostasis and loss of ESR1 promotes altered luminal-to-basal differentiation driven by ESR1-mediated paracrine factors from the stroma.

Disruption of CSF2RA in the bovine preimplantation embryo reduces development and affects embryonic gene expression in utero.

The hypothesis that CSF2 plays a role in the preimplantation development of the bovine embryo was tested by evaluating consequences of inactivation of CSF2RA (the functional receptor in the embryo) for development of embryos in utero. CRISPR/Cas9 was used to alter sequences on exon 5 and intron 5 of CSF2RA, Control embryos were injected with Cas9 mRNA only. Embryos > 16 cells at day 5 after insemination were transferred to synchronized recipient females in groups of 7 to 24. Embryos were flushed from the uterus two days later. The proportion of recovered embryos that developed to the blastocyst stage was lower for knockout embryos (39%) than for control embryos (63%). RNA sequencing of individual morulae and blastocysts indicated a total of 27 (morula) or 15 (blastocyst) differentially-expressed genes (false discovery rate <0.05). Gene set enrichment analysis indicated that the knockout affected genes playing roles in several functions including cell signaling and glycosylation. It was concluded that signaling through CSF2RA is not obligatory for development of the bovine preimplantation embryo to the blastocyst stage but that CSF2 signaling does enhance the likelihood that the embryo can become a blastocyst and result in specific changes in gene expression.

Single-nuclei RNA sequencing (snRNA-seq) uncovers trophoblast cell types and lineages in the mature bovine placenta.

Ruminants have a semi-invasive placenta, which possess highly vascularized placentomes formed by maternal endometrial caruncles and fetal placental cotyledons and required for fetal development to term. The synepitheliochorial placenta of cattle contains at least two trophoblast cell populations, including uninucleate (UNC) and binucleate (BNC) cells that are most abundant in the cotyledonary chorion of the placentomes. The interplacentomal placenta is more epitheliochorial in nature with the chorion developing specialized areolae over the openings of uterine glands. Of note, the cell types in the placenta and cellular and molecular mechanisms governing trophoblast differentiation and function are little understood in ruminants. To fill this knowledge gap, the cotyledonary and intercotyledonary areas of the mature day 195 bovine placenta were analyzed by single nuclei analysis. Single-nuclei RNA-seq analysis found substantial differences in cell type composition and transcriptional profiles between the two distinct regions of the placenta. Based on clustering and cell marker gene expression, five different trophoblast cell types were identified in the chorion, including proliferating and differentiating UNC and two different types of BNC in the cotyledon. Cell trajectory analyses provided a framework for understanding the differentiation of trophoblast UNC into BNC. The upstream transcription factor binding analysis of differentially expressed genes identified a candidate set of regulator factors and genes regulating trophoblast differentiation. This foundational information is useful to discover essential biological pathways underpinning the development and function of the bovine placenta.

A high-density genome-wide association with absolute blood monocyte count in domestic sheep identifies novel loci.

Monocytes are a core component of the immune system that arise from bone marrow and differentiate into cells responsible for phagocytosis and antigen presentation. Their derivatives are often responsible for the initiation of the adaptive immune response. Monocytes and macrophages are central in both controlling and propagating infectious diseases such as infection by Coxiella burnetii and small ruminant lentivirus in sheep. Genotypes from 513 Rambouillet, Polypay, and Columbia sheep (Ovis aries) were generated using the Ovine SNP50 BeadChip. Of these sheep, 222 animals were subsequently genotyped with the Ovine Infinium® HD SNP BeadChip to increase SNP coverage. Data from the 222 HD genotyped sheep were combined with the data from an additional 258 unique sheep to form a 480-sheep reference panel; this panel was used to impute the low-density genotypes to the HD genotyping density. Then, a genome-wide association analysis was conducted to identify loci associated with absolute monocyte counts from blood. The analysis used a single-locus mixed linear model implementing EMMAX with age and ten principal components as fixed effects. Two genome-wide significant peaks (p < 5x10-7) were identified on chromosomes 9 and 1, and ten genome-wide suggestive peaks (p < 1x10-5) were identified on chromosomes 1, 2, 3, 4, 9, 10, 15, and 16. The identified loci were within or near genes including KCNK9, involved into cytokine production, LY6D, a member of a superfamily of genes, some of which subset monocyte lineages, and HMGN1, which encodes a chromatin regulator associated with myeloid cell differentiation. Further investigation of these loci is being conducted to understand their contributions to monocyte counts. Investigating the genetic basis of monocyte lineages and numbers may in turn provide information about pathogens of veterinary importance and elucidate fundamental immunology.

An improved ovine reference genome assembly to facilitate in-depth functional annotation of the sheep genome.

BACKGROUND: The domestic sheep (Ovis aries) is an important agricultural species raised for meat, wool, and milk across the world. A high-quality reference genome for this species enhances the ability to discover genetic mechanisms influencing biological traits. Furthermore, a high-quality reference genome allows for precise functional annotation of gene regulatory elements. The rapid advances in genome assembly algorithms and emergence of sequencing technologies with increasingly long reads provide the opportunity for an improved de novo assembly of the sheep reference genome. FINDINGS: Short-read Illumina (55× coverage), long-read Pacific Biosciences (75× coverage), and Hi-C data from this ewe retrieved from public databases were combined with an additional 50× coverage of Oxford Nanopore data and assembled with canu v1.9. The assembled contigs were scaffolded using Hi-C data with Salsa v2.2, gaps filled with PBsuitev15.8.24, and polished with Nanopolish v0.12.5. After duplicate contig removal with PurgeDups v1.0.1, chromosomes were oriented and polished with 2 rounds of a pipeline that consisted of freebayes v1.3.1 to call variants, Merfin to validate them, and BCFtools to generate the consensus fasta. The ARS-UI_Ramb_v2.0 assembly is 2.63 Gb in length and has improved continuity (contig NG50 of 43.18 Mb), with a 19- and 38-fold decrease in the number of scaffolds compared with Oar_rambouillet_v1.0 and Oar_v4.0. ARS-UI_Ramb_v2.0 has greater per-base accuracy and fewer insertions and deletions identified from mapped RNA sequence than previous assemblies. CONCLUSIONS: The ARS-UI_Ramb_v2.0 assembly is a substantial improvement in contiguity that will optimize the functional annotation of the sheep genome and facilitate improved mapping accuracy of genetic variant and expression data for traits in sheep.

Variation in type two taste receptor genes is associated with bitter tasting phenylthiocarbamide consumption in mature Targhee and Rambouillet rams.

Bitter taste perception in sheep can lead to avoidance of specific types of forage, such as sagebrush, which is present on many rangeland grazing systems in the Intermountain West. In humans, bitter taste perception is influenced by variation in several TAS2R genes, including more extensively studied TAS2R38 and TAS2R16. We hypothesize that variation in taste receptor genes in sheep is associated with bitter taste. Therefore, the objective of this study was to examine variation in TAS2R genes in relation to consumption of a bitter tasting compound phenylthiocarbamide (PTC) which determines bitter "taster" and "non-taster" status in humans. Rambouillet and Targhee rams (n = 26) were offered various concentrations of PTC solution (0.2-12.29 mM) and water in a side-by-side presentation during two experiments. Blood was collected for DNA isolation and sequencing. Nineteen TAS2R genes were amplified and sequenced with long read Oxford Nanopore MinION technology. A total of 1,049 single nucleotide polymorphisms (SNPs) and 26 haplotypes were identified in these genes. Of these, 24 SNPs and 11 haplotypes were significantly (P < 0.05) associated with PTC consumption in TAS2R3, TAS2R5, TAS2R8, TAS2R9, TAS2R16, TAS2R31-like, TAS2R38, TAS2R39, and TAS2R42-like. Over 50% of the SNPs resulted in a change in amino acid sequence and several resided in potential regulatory regions, which could have downstream functional consequences and influence bitter taste perception in sheep. Further research is needed to validate these associations and elucidate the mechanisms that link variation in TAS2R genes to bitter taste perception in sheep. This may enable producers to select sheep more likely to consume bitter forage such as sagebrush as a flock and rangeland management strategy.

Characterizing Genetic Regulatory Elements in Ovine Tissues.

The Ovine Functional Annotation of Animal Genomes (FAANG) project, part of the broader livestock species FAANG initiative, aims to identify and characterize gene regulatory elements in domestic sheep. Regulatory element annotation is essential for identifying genetic variants that affect health and production traits in this important agricultural species, as greater than 90% of variants underlying genetic effects are estimated to lie outside of transcribed regions. Histone modifications that distinguish active or repressed chromatin states, CTCF binding, and DNA methylation were used to characterize regulatory elements in liver, spleen, and cerebellum tissues from four yearling sheep. Chromatin immunoprecipitation with sequencing (ChIP-seq) was performed for H3K4me3, H3K27ac, H3K4me1, H3K27me3, and CTCF. Nine chromatin states including active promoters, active enhancers, poised enhancers, repressed enhancers, and insulators were characterized in each tissue using ChromHMM. Whole-genome bisulfite sequencing (WGBS) was performed to determine the complement of whole-genome DNA methylation with the ChIP-seq data. Hypermethylated and hypomethylated regions were identified across tissues, and these locations were compared with chromatin states to better distinguish and validate regulatory elements in these tissues. Interestingly, chromatin states with the poised enhancer mark H3K4me1 in the spleen and cerebellum and CTCF in the liver displayed the greatest number of hypermethylated sites. Not surprisingly, active enhancers in the liver and spleen, and promoters in the cerebellum, displayed the greatest number of hypomethylated sites. Overall, chromatin states defined by histone marks and CTCF occupied approximately 22% of the genome in all three tissues. Furthermore, the liver and spleen displayed in common the greatest percent of active promoter (65%) and active enhancer (81%) states, and the liver and cerebellum displayed in common the greatest percent of poised enhancer (53%), repressed enhancer (68%), hypermethylated sites (75%), and hypomethylated sites (73%). In addition, both known and de novo CTCF-binding motifs were identified in all three tissues, with the highest number of unique motifs identified in the cerebellum. In summary, this study has identified the regulatory regions of genes in three tissues that play key roles in defining health and economically important traits and has set the precedent for the characterization of regulatory elements in ovine tissues using the Rambouillet reference genome.

Genome-Wide Association Analyses of Fertility Traits in Beef Heifers.

The ability of livestock to reproduce efficiently is critical to the sustainability of animal agriculture. Antral follicle count (AFC) and reproductive tract scores (RTS) can be used to estimate fertility in beef heifers, but the genetic mechanisms influencing variation in these measures are not well understood. Two genome-wide association studies (GWAS) were conducted to identify the significant loci associated with these traits. In total, 293 crossbred beef heifers were genotyped on the Bovine GGP 50K chip and genotypes were imputed to 836,121 markers. A GWAS was performed with the AFC phenotype for 217 heifers with a multi-locus mixed model, conducted using the year, age at time of sampling and principal component analysis groupings as the covariates. The RTS GWAS was performed with 289 heifers using an additive correlation/trend test comparing prepubertal to pubertal heifers. The loci on chromosomes 2, 3 and 23 were significant in the AFC GWAS and the loci on chromosomes 2, 8, 10 and 11 were significant in the RTS GWAS. The significant region on chromosome 2 was similar between both analyses. These regions contained genes associated with cell proliferation, transcription, apoptosis and development. This study proposes candidate genes for beef cattle fertility, although future research is needed to elucidate the precise mechanisms.

A Reference Genome Assembly of American Bison, Bison bison bison.

Bison are an icon of the American West and an ecologically, commercially, and culturally important species. Despite numbering in the hundreds of thousands today, conservation concerns remain for the species, including the impact on genetic diversity of a severe bottleneck around the turn of the 20th century and genetic introgression from domestic cattle. Genetic diversity and admixture are best evaluated at genome-wide scale, for which a high-quality reference is necessary. Here, we use trio binning of long reads from a bison-Simmental cattle (Bos taurus taurus) male F1 hybrid to sequence and assemble the genome of the American plains bison (Bison bison bison). The male haplotype genome is chromosome-scale, with a total length of 2.65 Gb across 775 scaffolds (839 contigs) and a scaffold N50 of 87.8 Mb. Our bison genome is ~13× more contiguous overall and ~3400× more contiguous at the contig level than the current bison reference genome. The bison genome sequence presented here (ARS-UCSC_bison1.0) will enable new research into the evolutionary history of this iconic megafauna species and provide a new tool for the management of bison populations in federal and commercial herds.

A Reference Genome Assembly of Simmental Cattle, Bos taurus taurus.

Genomics research has relied principally on the establishment and curation of a reference genome for the species. However, it is increasingly recognized that a single reference genome cannot fully describe the extent of genetic variation within many widely distributed species. Pangenome representations are based on high-quality genome assemblies of multiple individuals and intended to represent the broadest possible diversity within a species. A Bovine Pangenome Consortium (BPC) has recently been established to begin assembling genomes from more than 600 recognized breeds of cattle, together with other related species to provide information on ancestral alleles and haplotypes. Previously reported de novo genome assemblies for Angus, Brahman, Hereford, and Highland breeds of cattle are part of the initial BPC effort. The present report describes a complete single haplotype assembly at chromosome-scale for a fullblood Simmental cow from an F1 bison-cattle hybrid fetus by trio binning. Simmental cattle, also known as Fleckvieh due to their red and white spots, originated in central Europe in the 1830s as a triple-purpose breed selected for draught, meat, and dairy production. There are over 50 million Simmental cattle in the world, known today for their fast growth and beef yields. This assembly (ARS_Simm1.0) is similar in length to the other bovine assemblies at 2.86 Gb, with a scaffold N50 of 102 Mb (max scaffold 156.8 Mb) and meets or exceeds the continuity of the best Bos taurus reference assemblies to date.

Meiotic Recombination Differences in Rams from Three Breeds of Sheep in the US.

Meiotic recombination is an important contributor to genetic variation and ensures proper chromosome segregation during gametogenesis. Previous studies suggest that at least 1 crossover (CO) per chromosome arm is important to avoid mis-segregation. While the total number of COs per spermatocyte is known to differ in mice, this is only beginning to be evaluated in sheep. This study used a cytogenetic approach to quantify and compare the number of COs per spermatocyte in rams from 3 breeds of sheep: Suffolk, Icelandic, and Targhee. In total, 2,758 spermatocytes and over 170,000 COs were examined. Suffolk rams exhibited the lowest mean number of COs (61.1 ± 0.15) compared to Icelandic (63.5 ± 0.27) and Targhee (65.9 ± 0.26) rams. Significant differences in the number of COs per spermatocyte were observed between Suffolk, Icelandic, and Targhee breeds as well as within each breed. Additionally, the number and location of COs were characterized for homologous chromosomes in a subset of spermatocytes for each ram. A positive correlation was identified between the number of COs and the length of the homologous chromosome pair. Suffolk and Icelandic rams exhibited up to 7 COs per chromosome, while Targhee rams exhibited up to 9. Further, distinct CO location preferences on homologous chromosome pairs with 1, 2, 3, and 4 COs were observed in all 3 breeds. These data in sheep will aid in elucidating the mechanism of mammalian meiotic recombination, an important contributor to genetic diversity.

Complete Mitochondrial Genome Sequence of Bighorn Sheep.

We report here the complete mitochondrial genome sequence of a Rocky Mountain bighorn sheep (Ovis canadensis) in the United States. The circular genome has a size of 16,466 bp and contains 13 protein-coding genes, 22 tRNA genes, and 2 rRNA genes.