An atlas of regulatory elements in chicken: A resource for chicken genetics and genomics.

A comprehensive characterization of regulatory elements in the chicken genome across tissues will have substantial impacts on both fundamental and applied research. Here, we systematically identified and characterized regulatory elements in the chicken genome by integrating 377 genome-wide sequencing datasets from 23 adult tissues. In total, we annotated 1.57 million regulatory elements, representing 15 distinct chromatin states, and predicted about 1.2 million enhancer-gene pairs and 7662 super-enhancers. This functional annotation of the chicken genome should have wide utility on identifying regulatory elements accounting for gene regulation underlying domestication, selection, and complex trait regulation, which we explored. In short, this comprehensive atlas of regulatory elements provides the scientific community with a valuable resource for chicken genetics and genomics.

Histone remodeling reflects conserved mechanisms of bovine and human preimplantation development.

How histone modifications regulate changes in gene expression during preimplantation development in any species remains poorly understood. Using CUT&Tag to overcome limiting amounts of biological material, we profiled two activating (H3K4me3 and H3K27ac) and two repressive (H3K9me3 and H3K27me3) marks in bovine oocytes, 2-, 4-, and 8-cell embryos, morula, blastocysts, inner cell mass, and trophectoderm. In oocytes, broad bivalent domains mark developmental genes, and prior to embryonic genome activation (EGA), H3K9me3 and H3K27me3 co-occupy gene bodies, suggesting a global mechanism for transcription repression. During EGA, chromatin accessibility is established before canonical H3K4me3 and H3K27ac signatures. Embryonic transcription is required for this remodeling, indicating that maternally provided products alone are insufficient for reprogramming. Last, H3K27me3 plays a major role in restriction of cellular potency, as blastocyst lineages are defined by differential polycomb repression and transcription factor activity. Notably, inferred regulators of EGA and blastocyst formation strongly resemble those described in humans, as opposed to mice. These similarities suggest that cattle are a better model than rodents to investigate the molecular basis of human preimplantation development.

Transcriptome annotation of 17 porcine tissues using nanopore sequencing technology.

The annotation of animal genomes plays an important role in elucidating molecular mechanisms behind the genetic control of economically important traits. Here, we employed long-read sequencing technology, Oxford Nanopore Technology, to annotate the pig transcriptome across 17 tissues from two Yorkshire littermate pigs. More than 9.8 million reads were obtained from a single flow cell, and 69 781 unique transcripts at 50 108 loci were identified. Of these transcripts, 16 255 were found to be novel isoforms, and 22 344 were found at loci that were novel and unannotated in the Ensembl (release 102) and NCBI (release 106) annotations. Novel transcripts were mostly expressed in cerebellum, followed by lung, liver, spleen, and hypothalamus. By comparing the unannotated transcripts to existing databases, there were 21 285 (95.3%) transcripts matched to the NT database (v5) and 13 676 (61.2%) matched to the NR database (v5). Moreover, there were 4324 (19.4%) transcripts matched to the SwissProt database (v5), corresponding to 11 356 proteins. Tissue-specific gene expression analyses showed that 9749 transcripts were highly tissue-specific, and cerebellum contained the most tissue-specific transcripts. As the same samples were used for the annotation of cis-regulatory elements in the pig genome, the transcriptome annotation generated by this study provides an additional and complementary annotation resource for the Functional Annotation of Animal Genomes effort to comprehensively annotate the pig genome.

Prediction of transcript isoforms in 19 chicken tissues by Oxford Nanopore long-read sequencing.

To identify and annotate transcript isoforms in the chicken genome, we generated Nanopore long-read sequencing data from 68 samples that encompassed 19 diverse tissues collected from experimental adult male and female White Leghorn chickens. More than 23.8 million reads with mean read length of 790 bases and average quality of 18.2 were generated. The annotation and subsequent filtering resulted in the identification of 55,382 transcripts at 40,547 loci with mean length of 1,700 bases. We predicted 30,967 coding transcripts at 19,461 loci, and 16,495 lncRNA transcripts at 15,512 loci. Compared to existing reference annotations, we found ∼52% of annotated transcripts could be partially or fully matched while ∼47% were novel. Seventy percent of novel transcripts were potentially transcribed from lncRNA loci. Based on our annotation, we quantified transcript expression across tissues and found two brain tissues (i.e., cerebellum and cortex) expressed the highest number of transcripts and loci. Furthermore, ∼22% of the transcripts displayed tissue specificity with the reproductive tissues (i.e., testis and ovary) exhibiting the most tissue-specific transcripts. Despite our wide sampling, ∼20% of Ensembl reference loci were not detected. This suggests that deeper sequencing and additional samples that include different breeds, cell types, developmental stages, and physiological conditions, are needed to fully annotate the chicken genome. The application of Nanopore sequencing in this study demonstrates the usefulness of long-read data in discovering additional novel loci (e.g., lncRNA loci) and resolving complex transcripts (e.g., the longest transcript for the TTN locus).

Profiling of open chromatin in developing pig (Sus scrofa) muscle to identify regulatory regions.

There is very little information about how the genome is regulated in domestic pigs (Sus scrofa). This lack of knowledge hinders efforts to define and predict the effects of genetic variants in pig breeding programs. To address this knowledge gap, we need to identify regulatory sequences in the pig genome starting with regions of open chromatin. We used the "Improved Protocol for the Assay for Transposase-Accessible Chromatin (Omni-ATAC-Seq)" to identify putative regulatory regions in flash-frozen semitendinosus muscle from 24 male piglets. We collected samples from the smallest-, average-, and largest-sized male piglets from each litter through five developmental time points. Of the 4661 ATAC-Seq peaks identified that represent regions of open chromatin, >50% were within 1 kb of known transcription start sites. Differential read count analysis revealed 377 ATAC-Seq defined genomic regions where chromatin accessibility differed significantly across developmental time points. We found regions of open chromatin associated with downregulation of genes involved in muscle development that were present in small-sized fetal piglets but absent in large-sized fetal piglets at day 90 of gestation. The dataset that we have generated provides a resource for studies of genome regulation in pigs and contributes valuable functional annotation information to filter genetic variants for use in genomic selection in pig breeding programs.

Pig genome functional annotation enhances the biological interpretation of complex traits and human disease.

The functional annotation of livestock genomes is crucial for understanding the molecular mechanisms that underpin complex traits of economic importance, adaptive evolution and comparative genomics. Here, we provide the most comprehensive catalogue to date of regulatory elements in the pig (Sus scrofa) by integrating 223 epigenomic and transcriptomic data sets, representing 14 biologically important tissues. We systematically describe the dynamic epigenetic landscape across tissues by functionally annotating 15 different chromatin states and defining their tissue-specific regulatory activities. We demonstrate that genomic variants associated with complex traits and adaptive evolution in pig are significantly enriched in active promoters and enhancers. Furthermore, we reveal distinct tissue-specific regulatory selection between Asian and European pig domestication processes. Compared with human and mouse epigenomes, we show that porcine regulatory elements are more conserved in DNA sequence, under both rapid and slow evolution, than those under neutral evolution across pig, mouse, and human. Finally, we provide biological insights on tissue-specific regulatory conservation, and by integrating 47 human genome-wide association studies, we demonstrate that, depending on the traits, mouse or pig might be more appropriate biomedical models for different complex traits and diseases.

Large-Scale Multiplexing Permits Full-Length Transcriptome Annotation of 32 Bovine Tissues From a Single Nanopore Flow Cell.

A comprehensive annotation of transcript isoforms in domesticated species is lacking. Especially considering that transcriptome complexity and splicing patterns are not well-conserved between species, this presents a substantial obstacle to genomic selection programs that seek to improve production, disease resistance, and reproduction. Recent advances in long-read sequencing technology have made it possible to directly extrapolate the structure of full-length transcripts without the need for transcript reconstruction. In this study, we demonstrate the power of long-read sequencing for transcriptome annotation by coupling Oxford Nanopore Technology (ONT) with large-scale multiplexing of 93 samples, comprising 32 tissues collected from adult male and female Hereford cattle. More than 30 million uniquely mapping full-length reads were obtained from a single ONT flow cell, and used to identify and characterize the expression dynamics of 99,044 transcript isoforms at 31,824 loci. Of these predicted transcripts, 21% exactly matched a reference transcript, and 61% were novel isoforms of reference genes, substantially increasing the ratio of transcript variants per gene, and suggesting that the complexity of the bovine transcriptome is comparable to that in humans. Over 7,000 transcript isoforms were extremely tissue-specific, and 61% of these were attributed to testis, which exhibited the most complex transcriptome of all interrogated tissues. Despite profiling over 30 tissues, transcription was only detected at about 60% of reference loci. Consequently, additional studies will be necessary to continue characterizing the bovine transcriptome in additional cell types, developmental stages, and physiological conditions. However, by here demonstrating the power of ONT sequencing coupled with large-scale multiplexing, the task of exhaustively annotating the bovine transcriptome - or any mammalian transcriptome - appears significantly more feasible.

Simplification of culture conditions and feeder-free expansion of bovine embryonic stem cells.

Bovine embryonic stem cells (bESCs) extend the lifespan of the transient pluripotent bovine inner cell mass in vitro. After years of research, derivation of stable bESCs was only recently reported. Although successful, bESC culture relies on complex culture conditions that require a custom-made base medium and mouse embryonic fibroblasts (MEF) feeders, limiting the widespread use of bESCs. We report here simplified bESC culture conditions based on replacing custom base medium with a commercially available alternative and eliminating the need for MEF feeders by using a chemically-defined substrate. bESC lines were cultured and derived using a base medium consisting of N2B27 supplements and 1% BSA (NBFR-bESCs). Newly derived bESC lines were easy to establish, simple to propagate and stable after long-term culture. These cells expressed pluripotency markers and actively proliferated for more than 35 passages while maintaining normal karyotype and the ability to differentiate into derivatives of all three germ lineages in embryoid bodies and teratomas. In addition, NBFR-bESCs grew for multiple passages in a feeder-free culture system based on vitronectin and Activin A medium supplementation while maintaining pluripotency. Simplified conditions will facilitate the use of bESCs for gene editing applications and pluripotency and lineage commitment studies.

Transcription initiation mapping in 31 bovine tissues reveals complex promoter activity, pervasive transcription, and tissue-specific promoter usage.

Characterizing transcription start sites is essential for understanding the regulatory mechanisms that control gene expression. Recently, a new bovine genome assembly (ARS-UCD1.2) with high continuity, accuracy, and completeness was released; however, the functional annotation of the bovine genome lacks precise transcription start sites and contains a low number of transcripts in comparison to human and mouse. By using the RAMPAGE approach, this study identified transcription start sites at high resolution in a large collection of bovine tissues. We found several known and novel transcription start sites attributed to promoters of protein-coding and lncRNA genes that were validated through experimental and in silico evidence. With these findings, the annotation of transcription start sites in cattle reached a level comparable to the mouse and human genome annotations. In addition, we identified and characterized transcription start sites for antisense transcripts derived from bidirectional promoters, potential lncRNAs, mRNAs, and pre-miRNAs. We also analyzed the quantitative aspects of RAMPAGE to produce a promoter activity atlas, reaching highly reproducible results comparable to traditional RNA-seq. Coexpression networks revealed considerable use of tissue-specific promoters, especially between brain and testicle, which expressed several genes in common from alternate loci. Furthermore, regions surrounding coexpressed modules were enriched in binding factor motifs representative of each tissue. The comprehensive annotation of promoters in such a large collection of tissues will substantially contribute to our understanding of gene expression in cattle and other mammalian species, shortening the gap between genotypes and phenotypes.

Functional annotations of three domestic animal genomes provide vital resources for comparative and agricultural research.

Gene regulatory elements are central drivers of phenotypic variation and thus of critical importance towards understanding the genetics of complex traits. The Functional Annotation of Animal Genomes consortium was formed to collaboratively annotate the functional elements in animal genomes, starting with domesticated animals. Here we present an expansive collection of datasets from eight diverse tissues in three important agricultural species: chicken (Gallus gallus), pig (Sus scrofa), and cattle (Bos taurus). Comparative analysis of these datasets and those from the human and mouse Encyclopedia of DNA Elements projects reveal that a core set of regulatory elements are functionally conserved independent of divergence between species, and that tissue-specific transcription factor occupancy at regulatory elements and their predicted target genes are also conserved. These datasets represent a unique opportunity for the emerging field of comparative epigenomics, as well as the agricultural research community, including species that are globally important food resources.

A comparative analysis of chromatin accessibility in cattle, pig, and mouse tissues.

BACKGROUND: Although considerable progress has been made towards annotating the noncoding portion of the human and mouse genomes, regulatory elements in other species, such as livestock, remain poorly characterized. This lack of functional annotation poses a substantial roadblock to agricultural research and diminishes the value of these species as model organisms. As active regulatory elements are typically characterized by chromatin accessibility, we implemented the Assay for Transposase Accessible Chromatin (ATAC-seq) to annotate and characterize regulatory elements in pigs and cattle, given a set of eight adult tissues. RESULTS: Overall, 306,304 and 273,594 active regulatory elements were identified in pig and cattle, respectively. 71,478 porcine and 47,454 bovine regulatory elements were highly tissue-specific and were correspondingly enriched for binding motifs of known tissue-specific transcription factors. However, in every tissue the most prevalent accessible motif corresponded to the insulator CTCF, suggesting pervasive involvement in 3-D chromatin organization. Taking advantage of a similar dataset in mouse, open chromatin in pig, cattle, and mice were compared, revealing that the conservation of regulatory elements, in terms of sequence identity and accessibility, was consistent with evolutionary distance; whereas pig and cattle shared about 20% of accessible sites, mice and ungulates only had about 10% of accessible sites in common. Furthermore, conservation of accessibility was more prevalent at promoters than at intergenic regions. CONCLUSIONS: The lack of conserved accessibility at distal elements is consistent with rapid evolution of enhancers, and further emphasizes the need to annotate regulatory elements in individual species, rather than inferring elements based on homology. This atlas of chromatin accessibility in cattle and pig constitutes a substantial step towards annotating livestock genomes and dissecting the regulatory link between genome and phenome.

Chromatin remodeling in bovine embryos indicates species-specific regulation of genome activation.

The shift from maternal to embryonic control is a critical developmental milestone in preimplantation development. Widespread transcriptomic and epigenetic remodeling facilitate this transition from terminally differentiated gametes to totipotent blastomeres, but the identity of transcription factors (TF) and genomic elements regulating embryonic genome activation (EGA) are poorly defined. The timing of EGA is species-specific, e.g., the timing of murine and human EGA differ significantly. To deepen our understanding of mammalian EGA, here we profile changes in open chromatin during bovine preimplantation development. Before EGA, open chromatin is enriched for maternal TF binding, similar to that observed in humans and mice. During EGA, homeobox factor binding becomes more prevalent and requires embryonic transcription. A cross-species comparison of open chromatin during preimplantation development reveals strong similarity in the regulatory circuitry underlying bovine and human EGA compared to mouse. Moreover, TFs associated with murine EGA are not enriched in cattle or humans, indicating that cattle may be a more informative model for human preimplantation development than mice.

Transcript profiling of bovine embryos implicates specific transcription factors in the maternal-to-embryo transition.

Full-grown oocytes are transcriptionally quiescent. Following maturation and fertilization, the early stages of embryonic development occur in the absence (or low levels) of transcription that results in a period of development relying on maternally derived products (e.g., mRNAs and proteins). Two critical steps occur during the transition from maternal to embryo control of development: maternal mRNA clearance and embryonic genome activation with an associated dramatic reprogramming of gene expression required for further development. By combining an RNA polymerase II inhibitor with RNA sequencing, we were able not only to distinguish maternally derived from embryonic transcripts in bovine preimplantation embryos but also to establish that embryonic gene activation is required for clearance of maternal mRNAs as well as to identify putative transcription factors that are likely critical for early bovine development.

Transcriptome profiling of individual rhesus macaque oocytes and preimplantation embryos.

Early mammalian embryonic transcriptomes are dynamic throughout the process of preimplantation development. Cataloging of primate transcriptomics during early development has been accomplished in humans, but global characterization of transcripts is lacking in the rhesus macaque: a key model for human reproductive processes. We report here the systematic classification of individual macaque transcriptomes using RNA-Seq technology from the germinal vesicle stage oocyte through the blastocyst stage embryo. Major differences in gene expression were found between sequential stages, with the 4- to 8-cell stages showing the highest level of differential gene expression. Analysis of putative transcription factor binding sites also revealed a striking increase in key regulatory factors in 8-cell embryos, indicating a strong likelihood of embryonic genome activation occurring at this stage. Furthermore, clustering analyses of gene co-expression throughout this period resulted in distinct groups of transcripts significantly associated to the different embryo stages assayed. The sequence data provided here along with characterizations of major regulatory transcript groups present a comprehensive atlas of polyadenylated transcripts that serves as a useful resource for comparative studies of preimplantation development in humans and other species.

Salmonella enterica Serovars Enteritidis Infection Alters the Indigenous Microbiota Diversity in Young Layer Chicks.

Avian gastrointestinal (GI) tracts are highly populated with a diverse array of microorganisms that share a symbiotic relationship with their hosts and contribute to the overall health and disease state of the intestinal tract. The microbiome of the young chick is easily prone to alteration in its composition by both exogenous and endogenous factors, especially during the early posthatch period. The genetic background of the host and exposure to pathogens can impact the diversity of the microbial profile that consequently contributes to the disease progression in the host. The objective of this study was to profile the composition and structure of the gut microbiota in young chickens from two genetically distinct highly inbred lines. Furthermore, the effect of the Salmonella Enteritidis infection on altering the composition makeup of the chicken microbiome was evaluated through the 16S rRNA gene sequencing analysis. One-day-old layer chicks were challenged with S. Enteritidis and the host cecal microbiota profile as well as the degree of susceptibility to Salmonella infection was examined at 2 and 7 days post infection. Our result indicated that host genotype had a limited effect on resistance to S. Enteritidis infection. Alpha diversity, beta diversity, and overall microbiota composition were analyzed for four factors: host genotype, age, treatment, and postinfection time points. S. Enteritidis infection in young chicks was found to significantly reduce the overall diversity of the microbiota population with expansion of Enterobacteriaceae family. These changes indicated that Salmonella colonization in the GI tract of the chickens has a direct effect on altering the natural development of the GI microbiota. The impact of S. Enteritidis infection on microbial communities was also more substantial in the late stage of infection. Significant inverse correlation between Enterobacteriaceae and Lachnospiraceae family in both non-infected and infected groups, suggested possible antagonistic interaction between members of these two taxa, which could potentially influences the overall microbial population in the gut. Our results also revealed that genetic difference between two lines had minimal effect on the establishment of microbiota population. Overall, this study provided preliminary insights into the contributing role of S. Enteritidis in influencing the overall makeup of chicken's gut microbiota.