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.

Numerous long non-coding RNAs (lncRNAs) have been identified and their roles in gene regulation in humans, mice, and other model organisms studied; however, far less research has been focused on lncRNAs in farm animal species. While previous studies in chickens, cattle, and pigs identified lncRNAs in specific developmental stages or differentially expressed under specific conditions in a limited number of tissues, more comprehensive identification of lncRNAs in these species is needed. The goal of the FAANG Consortium (Functional Annotation of Animal Genomes) is to functionally annotate animal genomes, including the annotation of lncRNAs. As one of the FAANG pilot projects, lncRNAs were identified across eight tissues in two adult male biological replicates from chickens, cattle, and pigs. Comprehensive lncRNA annotations for the chicken, cattle, and pig genomes were generated by utilizing RNA-seq from eight tissue types from two biological replicates per species at the adult developmental stage. A total of 9393 lncRNAs in chickens, 7235 lncRNAs in cattle, and 14,429 lncRNAs in pigs were identified. Including novel isoforms and lncRNAs from novel loci, 5288 novel lncRNAs were identified in chickens, 3732 in cattle, and 4870 in pigs. These transcripts match previously known patterns of lncRNAs, such as generally lower expression levels than mRNAs and higher tissue specificity. An analysis of lncRNA conservation across species identified a set of conserved lncRNAs with potential functions associated with chromatin structure and gene regulation. Tissue-specific lncRNAs were identified. Genes proximal to tissue-specific lncRNAs were enriched for GO terms associated with the tissue of origin, such as leukocyte activation in spleen. LncRNAs were identified in three important farm animal species using eight tissues from adult individuals. About half of the identified lncRNAs were not previously reported in the NCBI annotations for these species. While lncRNAs are less conserved than protein-coding genes, a set of positionally conserved lncRNAs were identified among chickens, cattle, and pigs with potential functions related to chromatin structure and gene regulation. Tissue-specific lncRNAs have potential regulatory functions on genes enriched for tissue-specific GO terms. Future work will include epigenetic data from ChIP-seq experiments to further refine these annotations.

Telomeres are the specialized ends of chromosomes consisting of highly conserved repeat (5′-TTAGGG-3′)_n sequences. Lack of information regarding the existence of an in vivo telomere clock function in birds, conflicting data regarding telomere array length in the chicken model, and the paucity of molecular telomere information for other avian species led us to study telomere array organization within and among 18 species and subspecies of birds. Most of the species contained between 2% and 4% telomere sequence per diploid genome. Arrays spanning 0.5–10 kb (Class I) and 10–40 kb (Class II) were observed in all of the species studied. Extremely long arrays, ranging from hundreds of kilobases to 1–2 Mb (Class III) were observed in all except two raptor species, the northern goshawk and American bald eagle. In chicken, there was evidence for shortening of the Class II arrays in vivo, based on intraindividual comparisons of somatic versus germline tissues in birds of different ages; terminally differentiated erythrocyte arrays were, on average, 2.3 kb shorter than sperm (germline) arrays. This study provides the first evidence for the existence of telomere arrays significantly larger than have been described for any vertebrate species to date and for developmentally programmed in vivo telomere shortening in the Aves taxa. The novel finding of megabase-sized telomere arrays may be an important feature of avian karyotypes that contain a large number of very small genetic units, the microchromosomes.

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.

Here we present a high-resolution cytogenomic analysis of chicken microchromosome 16. We established the location of the major histocompatibility complex (MHC)-B and -Y subregions relative to each other and to the nucleolus organizer region (NOR) encoding the 18S–5.8S–28S ribosomal DNA. To do so, we employed multicolor fluorescence in situ hybridization using large-insert bacterial artificial chromosome clones with fully sequenced inserts or repetitive sequence probes specific for the subregion of interest. We show that the MHC-Y and -B regions are located on the same side of the NOR, rather than opposite ends, as previously proposed. On the q arm, the MHC-Y is closely adjacent to the NOR, whereas the MHC-B is distal near the q-terminus. A relatively large GC-rich region separates the 2 MHC subregions and includes a specialized structure, a secondary constriction. We propose that the GC-rich large physical distance is the basis for the lack of genetic linkage between the NOR and MHC-B and between the MHC-Y and -B. An integrated model for GGA 16 is presented that incorporates gene complex order in the context of key architectural features including p and q arms, primary (centromere) and secondary constrictions, telomeres, as well as AT- and GC-rich regions.

Specialized populations and mutant genetic lines of poultry have contributed to the biological, biomedical, and the agricultural sciences for over 100 years. Despite extensive loss of these unique genetic resources, there is no funding mechanism or long-term strategy to conserve and develop these important research resources. As the draft sequence of the chicken genome is completed during the coming year, the authors of this Policy Forum urge national initiatives to maintain and use specialized avian genetic variants.