Third Report on Chicken Genes and Chromosomes 2015

The molecule events expression of TGF-β/Smad signaling pathway in morphological and structural developmental characteristics of gonads in goose embryos

China is the largest producer and consumer of geese with significant social and economic value in agriculture. The Jilin White Goose, known for its excellent egg-laying and reproductive characteristics, is a prominent breeding breed in the northeast of China widely used for cross-breeding.Gonad development is a complex process, which will differentiate into testes or ovaries, thus affecting sex determination. Since the level of gonadal development in geese of different sexes is also related to their productivity, the study of gonadal function in geese is of great theoretical and practical importance. In this study, we used the embryonic gonads of Jilin white geese as the research object to reveal the morphologic and histologic developmental characteristics and process of the gonads of goose embryos from the 9th day to the 28th day and to investigate the mechanisms regulating the differentiation of the gonads. The results of ALP-PAS staining by histomorphological observations showed that the primordial germ cells in the female gonads were fully differentiated at embryonic day 19 (E19) and in the male gonads at embryonic day 18 (E18).Previous studies indicate that the TGF-β /Smad signaling pathway is crucial for sex differentiation and gamete formation in animals. In our study, we quantitatively analyzed key genes and proteins in this pathway within goose embryos. We found that TGF-βRI/II receptors and Smad2/3 are involved in gonadal development early in differentiation. Smad2 protein levels peaked at E18 in male gonads and E19 in female gonads, while Smad3 peaked at E10 for both. It has been suggested that there is a novel interaction between TGF-β and Wnt/β-catenin signaling pathway. In this study, the expression of LEF1, a protein downstream of the Wnt/β-catenin signaling pathway, peaked at E19 in the female gonads and at E18 in the male gonads, and the protein expression level of β-catenin was significantly increased in the female gonads at E19 and E28 and peaked at E28 in the male gonads. β-catenin expression levels were significantly higher in the male gonads than in the other periods and reached a minimum at E28. We injected the TMAO drug, which activates the TGF-β /Smad pathway, into goose embryos using the embryo injection technique to study how it affects gene and protein expression in this signaling pathway. We found that in female gonads, the levels of Smad2/3 proteins were lower than in the control group at E18, but higher than the control group at E19. In male gonads, the results for Smad3 proteins matched those of the females, while the results for Smad2 proteins were the opposite. This shows that TMAO activation of the TGF-β /Smad signaling pathway may impact gonadal development in goose embryos. In conclusion, this study suggest that the TGF-β /Smad and Wnt signaling pathways may enhance gonadal tissue development in goose embryos. The study identified key time points in gonadogenesis, providing a new theoretical basis for further research on the molecular mechanisms of gonadal development.

Near telomere-to-telomere assembly of the Tarim pigeon (Columba livia) genome

Pigeons serve as important model animals and commercial poultry. The Tarim pigeon, as a breed of Columba livia, is a locally indigenous breed unique to China. While the genome of C. livia was published in 2013, its assembly was fragmented and incomplete. In this study, we generated a near telomere-to-telomere assembly of the pigeon genome using the sequencing platform of PacBio HiFi, Nanopore long reads and Hi-C. The assembled genome spans 1295.8 Mb, with a contig N50 size of 49 Mb and a scaffold N50 size of 85 Mb. Approximately 98.4% of the assembly is anchored onto 41 chromosomes, with a BUSCO assessment indicating a completeness of 97.2%. And telomeres were identified at both ends of the four chromosomes. A total of 21,450 genes were annotated. The genome assembly of C. livia lays the foundation for understanding their genetic composition and evolutionary history and contributes to the pigeon breeding industry. Additionally, it will provide a basis for further management and conservation of pigeon breed diversity.

Insights into avian molecular cytogenetics-with reptilian comparisons

In last 100 years or so, much information has been accumulated on avian karyology, genetics, physiology, biochemistry and evolution. The chicken genome project generated genomic resources used in comparative studies, elucidating fundamental evolutionary processes, much of it funded by the economic importance of domestic fowl (which are also excellent model species in many areas). Studying karyotypes and whole genome sequences revealed population processes, evolutionary biology, and genome function, uncovering the role of repetitive sequences, transposable elements and gene family expansion. Knowledge of the function of many genes and non-expressed or identified regulatory components is however still lacking. Birds (Aves) are diverse, have striking adaptations for flight, migration and survival and inhabit all continents most islands. They also have a unique karyotype with ~ 10 macrochromosomes and ~ 30 microchromosomes that are smaller than other reptiles. Classified into Palaeognathae and Neognathae they are evolutionarily close, and a subset of reptiles. Here we overview avian molecular cytogenetics with reptilian comparisons, shedding light on their karyotypes and genome structure features. We consider avian evolution, then avian (followed by reptilian) karyotypes and genomic features. We consider synteny disruptions, centromere repositioning, and repetitive elements before turning to comparative avian and reptilian genomics. In this context, we review comparative cytogenetics and genome mapping in birds as well as Z- and W-chromosomes and sex determination. Finally, we give examples of pivotal research areas in avian and reptilian cytogenomics, particularly physical mapping and map integration of sex chromosomal genes, comparative genomics of chicken, turkey and zebra finch, California condor cytogenomics as well as some peculiar cytogenetic and evolutionary examples. We conclude that comparative molecular studies and improving resources continually contribute to new approaches in population biology, developmental biology, physiology, disease ecology, systematics, evolution and phylogenetic systematics orientation. This also produces genetic mapping information for chromosomes active in rearrangements during the course of evolution. Further insights into mutation, selection and adaptation of vertebrate genomes will benefit from these studies including physical and online resources for the further elaboration of comparative genomics approaches for many fundamental biological questions.

Annotations of four high-quality indigenous chicken genomes identify more than one thousand missing genes in subtelomeric regions and micro-chromosomes with high G/C contents

BACKGROUND

Although multiple chicken genomes have been assembled and annotated, the numbers of protein-coding genes in chicken genomes and their variation among breeds are still uncertain due to the low quality of these genome assemblies and limited resources used in their gene annotations. To fill these gaps, we recently assembled genomes of four indigenous chicken breeds with distinct traits at chromosome-level. In this study, we annotated genes in each of these assembled genomes using a combination of RNA-seq- and homology-based approaches.

RESULTS

We identified varying numbers (17,497-17,718) of protein-coding genes in the four indigenous chicken genomes, while recovering 51 of the 274 "missing" genes in birds in general, and 36 of the 174 "missing" genes in chickens in particular. Intriguingly, based on deeply sequenced RNA-seq data collected in multiple tissues in the four breeds, we found 571 ~ 627 protein-coding genes in each genome, which were missing in the annotations of the reference chicken genomes (GRCg6a and GRCg7b/w). After removing redundancy, we ended up with a total of 1,420 newly annotated genes (NAGs). The NAGs tend to be found in subtelomeric regions of macro-chromosomes (chr1 to chr5, plus chrZ) and middle chromosomes (chr6 to chr13, plus chrW), as well as in micro-chromosomes (chr14 to chr39) and unplaced contigs, where G/C contents are high. Moreover, the NAGs have elevated quadruplexes G frequencies, while both G/C contents and quadruplexes G frequencies in their surrounding regions are also high. The NAGs showed tissue-specific expression, and we were able to verify 39 (92.9%) of 42 randomly selected ones in various tissues of the four chicken breeds using RT-qPCR experiments. Most of the NAGs were also encoded in the reference chicken genomes, thus, these genomes might harbor more genes than previously thought.

CONCLUSION

The NAGs are widely distributed in wild, indigenous and commercial chickens, and they might play critical roles in chicken physiology. Counting these new genes, chicken genomes harbor more genes than originally thought.

High quality assemblies of four indigenous chicken genomes and related functional data resources

Many lines of evidence indicate that red jungle fowl (RJF) is the primary ancestor of domestic chickens. Although multiple versions of RJF (galgal2-galgal5 and GRCg6a) and commercial chickens (GRCg7b/w and Huxu) genomes have been assembled since 2004, no high-quality indigenous chicken genomes have been assembled, hampering the understanding of chicken domestication and evolution. To fill the gap, we sequenced the genomes of four indigenous chickens with distinct morphological traits in southwest China, using a combination of short, long and Hi-C reads. We assembled each genome (~1.0 Gb) into 42 chromosomes with chromosome N50 90.5-90.9 Mb, amongst the highest quality of chicken genome assemblies. To provide resources for gene annotation and functional analysis, we also sequenced transcriptomes of 10 tissues for each of the four chickens. Moreover, we corrected many mis-assemblies and assembled missing micro-chromosomes 29 and 34-39 for GRCg6a. Our assemblies, sequencing data and the correction of GRCg6a can be valuable resources for studying chicken domestication and evolution.

Identification, biogenesis, and function prediction of a novel circRNA_3238 of chicken

With the development of high-throughput sequencing, circular RNA has come into people's vision and attracted more and more attention. Many studies have found that circular RNA plays an important role in a variety of biological processes and the occurrence and development of diseases. According to the previous sequencing results, circRNA_3238 was differentially expressed in ALV-J infected group and the non-infected group was selected for subsequent verification and analysis. We found that circRNA_3238 is a stable, circular transcript, which mainly exists in the cytoplasm. And it is widely expressed in various tissues of chickens, and highly expressed in lung, lymph, and bursa of fabricius. Bioinformatics results show that circRNA_3238 and the predicted target genes enriched MAPK signaling pathway, Notch signaling pathway, and other pathways related to disease or immune, revealing circRNA_3238 may indirectly regulate the process of ALV-J infection by regulating target genes.

Integrative comparative analysis of avian chromosome evolution by in-silico mapping of the gene ontology of homologous synteny blocks and evolutionary breakpoint regions

Avian chromosomes undergo more intra- than interchromosomal rearrangements, which either induce or are associated with genome variations among birds. Evolving from a common ancestor with a karyotype not dissimilar from modern chicken, two evolutionary elements characterize evolutionary change: homologous synteny blocks (HSBs) constitute common conserved parts at the sequence level, while evolutionary breakpoint regions (EBRs) occur between HSBs, defining the points where rearrangement occurred. Understanding the link between the structural organization and functionality of HSBs and EBRs provides insight into the mechanistic basis of chromosomal change. Previously, we identified gene ontology (GO) terms associated with both; however, here we revisit our analyses in light of newly developed bioinformatic algorithms and the chicken genome assembly galGal6. We aligned genomes available for six birds and one lizard species, identifying 630 HSBs and 19 EBRs. We demonstrate that HSBs hold vast functionality expressed by GO terms that have been largely conserved through evolution. Particularly, we found that genes within microchromosomal HSBs had specific functionalities relevant to neurons, RNA, cellular transport and embryonic development, and other associations. Our findings suggest that microchromosomes may have conserved throughout evolution due to the specificity of GO terms within their HSBs. The detected EBRs included those found in the genome of the anole lizard, meaning they were shared by all saurian descendants, with others being unique to avian lineages. Our estimate of gene richness in HSBs supported the fact that microchromosomes contain twice as many genes as macrochromosomes.

Evolutionary analysis of a complete chicken genome

Microchromosomes are prevalent in nonmammalian vertebrates [P. D. Waters et al., Proc. Natl. Acad. Sci. U.S.A. 118 (2021)], but a few of them are missing in bird genome assemblies. Here, we present a new chicken reference genome containing all autosomes, a Z and a W chromosome, with all gaps closed except for the W. We identified ten small microchromosomes (termed dot chromosomes) with distinct sequence and epigenetic features, among which six were newly assembled. Those dot chromosomes exhibit extremely high GC content and a high level of DNA methylation and are enriched for housekeeping genes. The pericentromeric heterochromatin of dot chromosomes is disproportionately large and continues to expand with the proliferation of satellite DNA and testis-expressed genes. Our analyses revealed that the 41-bp CNM repeat frequently forms higher-order repeats (HORs) at the centromeres of acrocentric chromosomes. The centromere core regions where the kinetochore attaches often encompass telomeric sequence (TTAGGG)n, and in a one of the dot chromosomes, the centromere core recruits an endogenous retrovirus (ERV). We further demonstrate that the W chromosome shares some common features with dot chromosomes, having large arrays of hypermethylated tandem repeats. Finally, using the complete chicken chromosome models, we reconstructed a fine picture of chordate karyotype evolution, revealing frequent chromosomal fusions before and after vertebrate whole-genome duplications. Our sequence and epigenetic characterization of chicken chromosomes shed insights into the understanding of vertebrate genome evolution and chromosome biology.

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).

Chromatin Interaction Responds to Breast Muscle Development and Intramuscular Fat Deposition Between Chinese Indigenous Chicken and Fast-Growing Broiler

Skeletal muscle development and intramuscular fat (IMF) content, which positively contribute to meat production and quality, are regulated by precisely orchestrated processes. However, changes in three-dimensional chromatin structure and interaction, a newly emerged mediator of gene expression, during the skeletal muscle development and IMF deposition have remained unclear. In the present study, we analyzed the differences in muscle development and IMF content between one-day-old commercial Arbor Acres broiler (AA) and Chinese indigenous Lushi blue-shelled-egg chicken (LS) and performed Hi-C analysis on their breast muscles. Our results indicated that significantly higher IMF content, however remarkably lower muscle fiber diameter was detected in breast muscle of LS chicken compared to that of AA broiler. The chromatin intra-interaction was prior to inter-interaction in both AA and LS chicken, and chromatin inter-interaction was heavily focused on the small and gene-rich chromosomes. For genomic compartmentalization, no significant difference in the number of B type compartments was found, but AA had more A type compartments versus LS. The A/B compartment switching of AA versus LS showed more A to B switching than B to A switching. There were no significant differences in the average sizes and distributions of topologically associating domains (TAD). Additionally, approximately 50% of TAD boundaries were overlapping. The reforming and disappearing events of TAD boundaries were identified between AA and LS chicken breast muscles. Among these, the HMGCR gene was located in the TAD-boundary regions in AA broilers, but in TAD-interior regions in LS chickens, and the IGF2BP3 gene was located in the AA-unique TAD boundaries. Both HMGCR and IGF2BP3 genes exhibited increased mRNA expression in one-day-old AA broiler breast muscles. It was demonstrated that the IGF2BP3 and HMGCR genes regulated by TAD boundary sliding were potential biomarkers for chicken breast muscle development and IMF deposition. Our data not only provide a valuable understanding of higher-order chromatin dynamics during muscle development and lipid accumulation but also reveal new insights into the regulatory mechanisms of muscle development and IMF deposition in chicken.

The homeobox gene TGIF1 is required for chicken ovarian cortical development and generation of the juxtacortical medulla

During early embryogenesis in amniotic vertebrates, the gonads differentiate into either ovaries or testes. The first cell lineage to differentiate gives rise to the supporting cells: Sertoli cells in males and pre-granulosa cells in females. These key cell types direct the differentiation of the other cell types in the gonad, including steroidogenic cells. The gonadal surface epithelium and the interstitial cell populations are less well studied, and little is known about their sexual differentiation programs. Here, we show the requirement of the homeobox transcription factor gene TGIF1 for ovarian development in the chicken embryo. TGIF1 is expressed in the two principal ovarian somatic cell populations: the cortex and the pre-granulosa cells of the medulla. TGIF1 expression is associated with an ovarian phenotype in estrogen-mediated sex reversal experiments. Targeted misexpression and gene knockdown indicate that TGIF1 is required, but not sufficient, for proper ovarian cortex formation. In addition, TGIF1 is identified as the first known regulator of juxtacortical medulla development. These findings provide new insights into chicken ovarian differentiation and development, specifically cortical and juxtacortical medulla formation.

Assessing evolutionary and developmental transcriptome dynamics in homologous cell types

Background

During development, complex organ patterns emerge through the precise temporal and spatial specification of different cell types. On an evolutionary timescale, these patterns can change, resulting in morphological diversification. It is generally believed that homologous anatomical structures are built—largely—by homologous cell types. However, whether a common evolutionary origin of such cell types is always reflected in the conservation of their intrinsic transcriptional specification programs is less clear.

Results

Here, we developed a user‐friendly bioinformatics workflow to detect gene co‐expression modules and test for their conservation across developmental stages and species boundaries. Using a paradigm of morphological diversification, the tetrapod limb, and single‐cell RNA‐sequencing data from two distantly related species, chicken and mouse, we assessed the transcriptional dynamics of homologous cell types during embryonic patterning. With mouse limb data as reference, we identified 19 gene co‐expression modules with varying tissue or cell type‐restricted activities. Testing for co‐expression conservation revealed modules with high evolutionary turnover, while others seemed maintained—to different degrees, in module make‐up, density or connectivity—over developmental and evolutionary timescales.

Conclusions

We present an approach to identify evolutionary and developmental dynamics in gene co‐expression modules during patterning‐relevant stages of homologous cell type specification using single‐cell RNA‐sequencing data.

We present an approach to identify evolutionary and developmental dynamics in gene co‐expression modules during patterning‐relevant stages of homologous cell type specification using single‐cell RNA‐sequencing data.

Domestic chicken diversity: Origin, distribution, and adaptation

Chicken is the most numerous among the domesticated livestock species. Across cultures, religions, and societies, chicken is widely accepted with little or no taboo compared to other domestic animals. Its adaptability to diverse environmental conditions and demonstrated potential for breeding improvement provide a unique genetic resource for addressing the challenges of food security in a world impacted by climatic change and human population growth. Recent studies, shedding new knowledge on the chicken genomes, have helped reconstruct its past evolutionary history. Here, we review the literature concerning the origin, dispersion, and adaptation of domestic chicken. We highlight the role of human and natural selection in shaping the diversity of the species and provide a few examples of knowledge gaps that may be the focus of future research.

Sustainable intensification of indigenous village chicken production system: matching the genotype with the environment

The multi-purpose indigenous village chickens (IVCs) are kept in low- and middle-income countries. IVCs are hardy and are resilient to disease, hostile environment, global warming, and climate change. The IVCs are a little impacted by anthropogenic effects; consequently, they possess high genetic and phenotypic diversity. Likewise, the genetic structure of IVCs is principally shaped by natural selection, which enables them to accumulate high genetic polymorphism and to adaptively radiate. Regardless of this, the genetic wealth of IVCs has been eroded by indiscriminate crossbreeding. Emerging infectious and non-infectious diseases, flawed assumptions, predation, inadequate nutrition, poorly maintained night enclosures, and underdeveloped market infrastructure, and the overlooked multiple-use values and unique attributes of IVCs have threatened their mere survival. The IVCs lay a few eggs and produce less meat, which cannot meet the growing (existing) demand. However, the demand for IVC products is growing attributable to the flavor and texture of eggs and meat, and they are well-aligned with the subsistence farming system. The several use values and ecosystem services provided by IVCs have been increasingly realized. Enhanced production can be attained through sustainable use of local (genetic) resources and by scaling up and out best practices. Genetic improvement needs to mainly rely upon IVC genetic resources and should have to match the genotype with the environment. Moreover, it has to maintain the genetic polymorphism that has been accumulated for time immemorial to respond to unanticipated changes in the production system and consumers' demand. In this review, enhanced management, selection strategies, and genetic crosses including the crossing of commercial chickens with red junglefowl have been proposed to sustainably intensify the IVC production system.

When and how do sex-linked regions become sex chromosomes?

The attention given to heteromorphism and genetic degeneration of "classical sex chromosomes" (Y chromosomes in XY systems, and the W in ZW systems that were studied first and are best described) has perhaps created the impression that the absence of recombination between sex chromosomes is inevitable. I here argue that continued recombination is often to be expected, that absence of recombination is surprising and demands further study, and that the involvement of selection in reduced recombination is not yet well understood. Despite a long history of investigations of sex chromosome pairs, there is a need for more quantitative approaches to studying sex-linked regions. I describe a scheme to help understand the relationships between different properties of sex-linked regions. Specifically, I focus on their sizes (differentiating between small regions and extensive fully sex-linked ones), the times when they evolved, and their differentiation, and review studies using DNA sequencing in nonmodel organisms that are providing information about the processes causing these properties.

Comparative Mapping of the Macrochromosomes of Eight Avian Species Provides Further Insight into Their Phylogenetic Relationships and Avian Karyotype Evolution

Avian genomes typically consist of ~10 pairs of macro- and ~30 pairs of microchromosomes. While inter-chromosomally, a pattern emerges of very little change (with notable exceptions) throughout evolution, intrachromosomal changes remain relatively poorly studied. To rectify this, here we use a pan-avian universally hybridising set of 74 chicken bacterial artificial chromosome (BAC) probes on the macrochromosomes of eight bird species: common blackbird, Atlantic canary, Eurasian woodcock, helmeted guinea fowl, houbara bustard, mallard duck, and rock dove. A combination of molecular cytogenetic, bioinformatics, and mathematical analyses allowed the building of comparative cytogenetic maps, reconstruction of a putative Neognathae ancestor, and assessment of chromosome rearrangement patterns and phylogenetic relationships in the studied neognath lineages. We observe that, as with our previous studies, chicken appears to have the karyotype most similar to the ancestor; however, previous reports of an increased rate of intrachromosomal change in Passeriformes (songbirds) appear not to be the case in our dataset. The use of this universally hybridizing probe set is applicable not only for the re-tracing of avian karyotype evolution but, potentially, for reconstructing genome assemblies.

Identification, biogenesis and function prediction of novel circRNA during the chicken ALV-J infection

Circular RNA (circRNA) is a new non-coding RNA with a highly conserved and stable covalently closed loop structure, and it plays an important role in a variety of biological processes and the occurrence of diseases. Based on the sequencing results, circRNA_3079 had the most significant difference between the infected group and normal group, up to about 8 times. It has attracted our attention and was selected for further verification and analysis. Though the characteristics analysis of circRNA_3079 in chicken, we found circRNA_3079 is a stable, circular transcript, which mainly exists in the cytoplasm. And it is widely expressed in various tissues of chickens, and highly expressed in lung, spleen, lymph and bursa of fabricius. Bioinformatics analysis results showed that circRNA_3079 and the predicted target genes are enriched in many pathways related to immunity or tumors, such as p53 signaling pathway, Jak-STAT signaling pathway and NOD-like receptor signaling pathway, which revealed that circRNA_3079 may indirectly regulate the ALV-J infection process through the regulation of target genes.HIGHLIGHTSCircRNA_3079 is an abundant and stable circular RNA expressed in different tissues and cells in chicken.The circularization of circRNA_3079 does not depend on the reverse complementary repetitive sequence of the flanking sequence.CircRNA_3079 may indirectly regulate the ALV-J infection process.

Mapping epigenetic modifications on chicken lampbrush chromosomes

Background

The epigenetic regulation of genome is crucial for implementation of the genetic program of ontogenesis through establishing and maintaining differential gene expression. Thus mapping of various epigenetic modifications to the genome is relevant for studying the regulation of gene expression. Giant transcriptionally active lampbrush chromosomes are an established tool for high resolution physical mapping of the genome and its epigenetic modifications. This study is aimed at characterizing the epigenetic status of compact chromatin domains (chromomeres) of chicken lampbrush macrochromosomes.

Results

Distribution of three epigenetic modifications – 5-methylcytosine, histone H3 trimethylated at lysine 9 and hyperacetylated histone H4 – along the axes of chicken lampbrush chromosomes 1–4, Z and W was analyzed in details. Enrichment of chromatin domains with the investigated epigenetic modifications was indicated on the cytological chromomere-loop maps for corresponding chicken lampbrush chromosomes. Heterogeneity in the distribution of 5-methylcytosine and histone H3 trimethylated at lysine 9 along the chromosome axes was revealed.

Conclusions

On examples of certain chromomeres of chicken lampbrush chromosomes 1, 3, 4 and W we demonstrated that a combination of immunofluorescent staining and fluorescence in situ hybridization allows to relate the epigenetic status and a DNA sequence context of individual chromomeres.

Associations between IGF1, IGFBP2 and TGFß3 Genes Polymorphisms and Growth Performance of Broiler Chicken Lines

Marker-assisted selection based on fast and accurate molecular analysis of individual genes is considered an acceptable tool in the speed-up of the genetic improvement of production performance in chickens. The objective of this study was to detect the single nucleotide polymorphisms (SNPs) in the IGF1, IGFBP2 and TGFß3 genes, and to investigate their associations with growth performance (body weight (BW) and average daily gain (ADG) at 14, 21, 28, 35 and 42 days of age) and carcass traits in broilers. Performance (carcass) data (weight before slaughter; weights of the trunk, giblets, abdominal fat, breast muscle and thigh muscle; slaughter value and slaughter percentage), as well as blood samples for DNA extraction and SNP analysis, were obtained from 97 chickens belonging to two different lines (Hubbard F15 and Cobb E) equally divided between the two sexes. The genotypes were detected using polymerase chain reaction- restriction fragment length polymorphism (PCR-RFLP) methods with specific primers and restrictase for each gene. The statistical analysis discovered significant associations (p < 0.05) between the TGFβ3 SNP and the following parameters: BW at 21, 28 and 35 days, trunk weight and slaughter value. Association analysis of BWs (at 21, 28 and 35 days) and SNPs was always significant for codominant, dominant and overdominant genetic models, showing a possible path for genomic selection in these chicken lines. Slaughter value was significant for codominant, recessive and overdominant patterns, whereas other carcass traits were not influenced by SNPs. Based on the results of this study, we suggested that the TGFβ3 gene could be used as a candidate gene marker for chicken growth traits in the Hubbard F15 and Cobb E population selection programs, whereas for carcass traits further investigation is needed.

An eQTL in the cystathionine beta synthase gene is linked to osteoporosis in laying hens

BACKGROUND

Skeletal damage is a challenge for laying hens because the physiological adaptations required for egg laying make them susceptible to osteoporosis. Previously, we showed that genetic factors explain 40% of the variation in end of lay bone quality and we detected a quantitative trait locus (QTL) of large effect on chicken chromosome 1. The aim of this study was to combine data from the commercial founder White Leghorn population and the F2 mapping population to fine-map this QTL and understand its function in terms of gene expression and physiology.

RESULTS

Several single nucleotide polymorphisms on chromosome 1 between 104 and 110 Mb (galGal6) had highly significant associations with tibial breaking strength. The alternative genotypes of markers of large effect that flanked the region had tibial breaking strengths of 200.4 vs. 218.1 Newton (P < 0.002) and, in a subsequent founder generation, the higher breaking strength genotype was again associated with higher breaking strength. In a subsequent generation, cortical bone density and volume were increased in individuals with the better bone genotype but with significantly reduced medullary bone quality. The effects on cortical bone density were confirmed in a further generation and was accompanied by increased mineral maturity of the cortical bone as measured by infrared spectrometry and there was evidence of better collagen cross-linking in the cortical bone. Comparing the transcriptome of the tibia from individuals with good or poor bone quality genotypes indicated four differentially-expressed genes at the locus, one gene, cystathionine beta synthase (CBS), having a nine-fold higher expression in the genotype for low bone quality. The mechanism was cis-acting and although there was an amino-acid difference in the CBS protein between the genotypes, there was no difference in the activity of the enzyme. Plasma homocysteine concentration, the substrate of CBS, was higher in the poor bone quality genotype.

CONCLUSIONS

Validated markers that predict bone strength have been defined for selective breeding and a gene was identified that may suggest alternative ways to improve bone health in addition to genetic selection. The identification of how genetic variants affect different aspects of bone turnover shows potential for translational medicine.

Chickspress: a resource for chicken gene expression

High-throughput sequencing and proteomics technologies are markedly increasing the amount of RNA and peptide data that are available to researchers, which are typically made publicly available via data repositories such as the NCBI Sequence Read Archive and proteome archives, respectively. These data sets contain valuable information about when and where gene products are expressed, but this information is not readily obtainable from archived data sets. Here we report Chickspress (http://geneatlas.arl.arizona.edu), the first publicly available gene expression resource for chicken tissues. Since there is no single source of chicken gene models, Chickspress incorporates both NCBI and Ensembl gene models and links these gene sets with experimental gene expression data and QTL information. By linking gene models from both NCBI and Ensembl gene prediction pipelines, researchers can, for the first time, easily compare gene models from each of these prediction workflows to available experimental data for these products. We use Chickspress data to show the differences between these gene annotation pipelines. Chickspress also provides rapid search, visualization and download capacity for chicken gene sets based upon tissue type, developmental stage and experiment type. This first Chickspress release contains 161 gene expression data sets, including expression of mRNAs, miRNAs, proteins and peptides. We provide several examples demonstrating how researchers may use this resource.

Transcriptomic analysis of early B-cell development in the chicken embryo

The chicken bursa of Fabricius is a primary lymphoid tissue important for B-cell development. Our long-term goal is to understand the role of bursal microenvironment in an early B-cell differentiation event initiating repertoire development through immunoglobulin gene conversion in the chick embryo. We hypothesize that early bursal B-cell differentiation is guided by signals through cytokine receptors. Our theory is based on previous evidence for expression of the receptor tyrosine kinase superfamily members and interleukin receptors in unseparated populations of bursal B-cells and bursal tissue. Knowledge of the expressed genes that are responsible for B-cell differentiation is a prerequisite for understanding the bursal microenvironment's function. This project uses transcriptomic analysis to evaluate gene expression across early B-cell development. RNA-seq was performed with total RNA isolated from bursal B-cells at embryonic day (ED) 16 and ED 19 (n = 3). Approximately 90 million high-quality clean reads were obtained from the cDNA libraries. The analysis revealed differentially expressed genes involved in the Jak-STAT pathway, Wnt signaling pathway, MAPK signaling pathway, metabolic pathways including tyrosine metabolism, Toll-like receptor signaling pathway, and cell-adhesion molecules. The genes predicted to encode surface receptors, signal transduction proteins, and transcription factors identified in this study represent gene candidates for controlling B-cell development in response to differentiation factors in the bursal microenvironment.

Together But Different: The Subgenomes of the Bimodal Eleutherine Karyotypes Are Differentially Organized

Bimodal karyotypes are characterized by the presence of two sets of chromosomes of contrasting size. Eleutherine bulbosa (2n = 12) presents a bimodal karyotype with a large chromosome pair, which has a pericentric inversion in permanent heterozygosity with suppressed recombination, and five pairs of three to four times smaller chromosomes. Aiming to understand whether high copy number sequence composition differs between both chromosome sets, we investigated the repetitive DNA fraction of E. bulbosa and compared it to the chromosomal organization of the related Eleutherine latifolia species, not containing the pericentric inversion. We also compared the repetitive sequence proportions between the heteromorphic large chromosomes of E. bulbosa and between E. bulbosa and E. latifolia to understand the influence of the chromosome inversion on the dynamics of repetitive sequences. The most abundant repetitive families of the genome showed a similar chromosomal distribution in both homologs of the large pair and in both species, apparently not influenced by the species-specific inversions. The repeat families Ebusat1 and Ebusat4 are localized interstitially only on the large chromosome pair, while Ebusat2 is located in the centromeric region of all chromosomes. The four most abundant retrotransposon lineages are accumulated in the large chromosome pair. Replication timing and distribution of epigenetic and transcriptional marks differ between large and small chromosomes. The differential distribution of retroelements appears to be related to the bimodal condition and is not influenced by the nonrecombining chromosome inversions in these species. Thus, the large and small chromosome subgenomes of the bimodal Eleutherine karyotype are differentially organized and probably evolved by repetitive sequences accumulation on the large chromosome set.

Characterization of the endogenous retrovirus insertion in CYP19A1 associated with henny feathering in chicken

Background

Henny feathering in chickens is determined by a dominant mutation that transforms male-specific plumage to female-like plumage. Previous studies indicated that this phenotype is caused by ectopic expression in skin of CYP19A1 encoding aromatase that converts androgens to estrogen and thereby inhibits the development of male-specific plumage. A long terminal repeat (LTR) from an uncharacterized endogenous retrovirus (ERV) insertion was found in an isoform of the CYP19A1 transcript from henny feathering chicken. However, the complete sequence and the genomic position of the insertion were not determined.

Results

We used publicly available whole genome sequence data to determine the flanking sequences of the ERV, and then PCR amplified the entire insertion and sequenced it using Nanopore long reads and Sanger sequencing. The 7524 bp insertion contains an intact endogenous retrovirus that was not found in chickens representing 31 different breeds not showing henny feathering or in samples of the ancestral red junglefowl. The sequence shows over 99% sequence identity to the avian leukosis virus ev-1 and ev-21 strains, suggesting a recent integration. The ERV 3’LTR, containing a powerful transcriptional enhancer and core promoter with TATA box together with binding sites for EFIII and Ig/EBP inside the CYP19A1 5′ untranslated region, was detected partially in an aromatase transcript, which present a plausible explanation for ectopic expression of aromatase in non-ovarian tissues underlying the henny feathering phenotype.

Conclusions

We demonstrate that the henny feathering allele harbors an insertion of an intact avian leukosis virus at the 5’end of CYP19A1. The presence of this ERV showed complete concordance with the henny feathering phenotype both within a pedigree segregating for this phenotype and across breeds.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0181-4) contains supplementary material, which is available to authorized users.

3D organization of chicken genome demonstrates evolutionary conservation of topologically associated domains and highlights unique architecture of erythrocytes' chromatin

How chromosomes are folded, spatially organized and regulated in three dimensions inside the cell nucleus are among the longest standing questions in cell biology. Genome-wide chromosome conformation capture (Hi-C) technique allowed identifying and characterizing spatial chromatin compartments in several mammalian species. Here, we present the first genome-wide analysis of chromatin interactions in chicken embryonic fibroblasts (CEF) and adult erythrocytes. We showed that genome of CEF is partitioned into topologically associated domains (TADs), distributed in accordance with gene density, transcriptional activity and CTCF-binding sites. In contrast to mammals, where all examined somatic cell types display relatively similar spatial organization of genome, chicken erythrocytes strongly differ from fibroblasts, showing pronounced A- and B- compartments, absence of typical TADs and formation of long-range chromatin interactions previously observed on mitotic chromosomes. Comparing mammalian and chicken genome architectures, we provide evidence highlighting evolutionary role of chicken TADs and their significance in genome activity and regulation.

Immune genes are hotspots of shared positive selection across birds and mammals

Consistent patterns of positive selection in functionally similar genes can suggest a common selective pressure across a group of species. We use alignments of orthologous protein-coding genes from 39 species of birds to estimate parameters related to positive selection for 11,000 genes conserved across birds. We show that functional pathways related to the immune system, recombination, lipid metabolism, and phototransduction are enriched for positively selected genes. By comparing our results with mammalian data, we find a significant enrichment for positively selected genes shared between taxa, and that these shared selected genes are enriched for viral immune pathways. Using pathogen-challenge transcriptome data, we show that genes up-regulated in response to pathogens are also enriched for positively selected genes. Together, our results suggest that pathogens, particularly viruses, consistently target the same genes across divergent clades, and that these genes are hotspots of host-pathogen conflict over deep evolutionary time.

Chromosome-level assembly reveals extensive rearrangement in saker falcon and budgerigar, but not ostrich, genomes

BACKGROUND

The number of de novo genome sequence assemblies is increasing exponentially; however, relatively few contain one scaffold/contig per chromosome. Such assemblies are essential for studies of genotype-to-phenotype association, gross genomic evolution, and speciation. Inter-species differences can arise from chromosomal changes fixed during evolution, and we previously hypothesized that a higher fraction of elements under negative selection contributed to avian-specific phenotypes and avian genome organization stability. The objective of this study is to generate chromosome-level assemblies of three avian species (saker falcon, budgerigar, and ostrich) previously reported as karyotypically rearranged compared to most birds. We also test the hypothesis that the density of conserved non-coding elements is associated with the positions of evolutionary breakpoint regions.

RESULTS

We used reference-assisted chromosome assembly, PCR, and lab-based molecular approaches, to generate chromosome-level assemblies of the three species. We mapped inter- and intrachromosomal changes from the avian ancestor, finding no interchromosomal rearrangements in the ostrich genome, despite it being previously described as chromosomally rearranged. We found that the average density of conserved non-coding elements in evolutionary breakpoint regions is significantly reduced. Fission evolutionary breakpoint regions have the lowest conserved non-coding element density, and intrachromomosomal evolutionary breakpoint regions have the highest.

CONCLUSIONS

The tools used here can generate inexpensive, efficient chromosome-level assemblies, with > 80% assigned to chromosomes, which is comparable to genomes assembled using high-density physical or genetic mapping. Moreover, conserved non-coding elements are important factors in defining where rearrangements, especially interchromosomal, are fixed during evolution without deleterious effects.

Banding cytogenetics of chimeric hybrids Coturnixcoturnix × Coturnixjaponica and comparative analysis with the domestic fowl

The Common quail Coturnixcoturnix Linnaeus, 1758 is a wild migratory bird which is distributed in Eurasia and North Africa, everywhere with an accelerating decline in population size. This species is protected by the Bonn and Berne conventions (1979) and by annex II/1 of the Birds Directive (2009). In Algeria, its breeding took place at the hunting centre in the west of the country. Breeding errors caused uncontrolled crosses between the Common quail and Japanese quail Coturnixjaponica Temminck & Schlegel, 1849. In order to help to preserve the natural genetic heritage of the Common quail and to lift the ambiguity among the populations of quail raised in Algeria, it seemed essential to begin to describe the chromosomes of this species in the country since no cytogenetic study has been reported to date. Fibroblast cultures from embryo and adult animal were initiated. Double synchronization with excess thymidine allowed us to obtain high resolution chromosomes blocked at prometaphase stage. The karyotype and the idiogram in GTG morphological banding (G-bands obtained with trypsin and Giemsa) corresponding to larger chromosomes 1-12 and ZW pair were thus established. The diploid set of chromosomes was estimated as 2N=78. Cytogenetic analysis of expected hybrid animals revealed the presence of a genetic introgression and cellular chimerism. This technique is effective in distinguishing the two quail taxa. Furthermore, the comparative chromosomal analysis of the two quails and domestic chicken Gallusgallusdomesticus Linnaeus, 1758 has been conducted. Differences in morphology and/or GTG band motifs were observed on 1, 2, 4, 7, 8 and W chromosomes. Neocentromere occurrence was suggested for Common quail chromosome 1 and Chicken chromosomes 4 and W. Double pericentric inversion was observed on the Common quail chromosome 2 while pericentric inversion hypothesis was proposed for Chicken chromosome 8. A deletion on the short arm of the Common quail chromosome 7 was also found. These results suggest that Common quail would be a chromosomally intermediate species between Chicken and Japanese quail. The appearance of only a few intrachromosomal rearrangements that occurred during evolution suggests that the organization of the genome is highly conserved between these three galliform species.

Combination of novel and public RNA-seq datasets to generate an mRNA expression atlas for the domestic chicken

Background

The domestic chicken (Gallus gallus) is widely used as a model in developmental biology and is also an important livestock species. We describe a novel approach to data integration to generate an mRNA expression atlas for the chicken spanning major tissue types and developmental stages, using a diverse range of publicly-archived RNA-seq datasets and new data derived from immune cells and tissues.

Results

Randomly down-sampling RNA-seq datasets to a common depth and quantifying expression against a reference transcriptome using the mRNA quantitation tool Kallisto ensured that disparate datasets explored comparable transcriptomic space. The network analysis tool Graphia was used to extract clusters of co-expressed genes from the resulting expression atlas, many of which were tissue or cell-type restricted, contained transcription factors that have previously been implicated in their regulation, or were otherwise associated with biological processes, such as the cell cycle. The atlas provides a resource for the functional annotation of genes that currently have only a locus ID. We cross-referenced the RNA-seq atlas to a publicly available embryonic Cap Analysis of Gene Expression (CAGE) dataset to infer the developmental time course of organ systems, and to identify a signature of the expansion of tissue macrophage populations during development.

Conclusion

Expression profiles obtained from public RNA-seq datasets – despite being generated by different laboratories using different methodologies – can be made comparable to each other. This meta-analytic approach to RNA-seq can be extended with new datasets from novel tissues, and is applicable to any species.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4972-7) contains supplementary material, which is available to authorized users.

Whole-Genome Resequencing of Red Junglefowl and Indigenous Village Chicken Reveal New Insights on the Genome Dynamics of the Species

The red junglefowl Gallus gallus is the main progenitor of domestic chicken, the commonest livestock species, outnumbering humans by an approximate ratio of six to one. The genetic control for production traits have been well studied in commercial chicken, but the selection pressures underlying unique adaptation and production to local environments remain largely unknown in indigenous village chicken. Likewise, the genome regions under positive selection in the wild red junglefowl remain untapped. Here, using the pool heterozygosity approach, we analyzed indigenous village chicken populations from Ethiopia, Saudi Arabia, and Sri Lanka, alongside six red junglefowl, for signatures of positive selection across the autosomes. Two red junglefowl candidate selected regions were shared with all domestic chicken populations. Four candidates sweep regions, unique to and shared among all indigenous domestic chicken, were detected. Only one region includes annotated genes (TSHR and GTF2A1). Candidate regions that were unique to each domestic chicken population with functions relating to adaptation to temperature gradient, production, reproduction and immunity were identified. Our results provide new insights on the consequence of the selection pressures that followed domestication on the genome landscape of the domestic village chicken.

Banding cytogenetics of the Barbary partridge Alectoris barbara and the Chukar partridge Alectoris chukar (Phasianidae): a large conservation with Domestic fowl Gallus domesticus revealed by high resolution chromosomes

The development of avian cytogenetics is significantly behind that of mammals. In fact, since the advent of cytogenetic techniques, fewer than 1500 karyotypes have been established. The Barbary partridge Alectoris barbara Bonnaterre, 1790 is a bird of economic interest but its genome has not been studied so far. This species is endemic to North Africa and globally declining. The Chukar partridge Alectoris chukar Gray, 1830 is an introduced species which shares the same habitat area as the Barbary partridge and so there could be introgressive hybridisation. A cytogenetic study has been initiated in order to contribute to the Barbary partridge and the Chukar partridge genome analyses. The GTG, RBG and RHG-banded karyotypes of these species have been described. Primary fibroblast cell lines obtained from embryos were harvested after simple and double thymidine synchronisation. The first eight autosomal pairs and Z sex chromosome have been described at high resolution and compared to those of the domestic fowl Gallus domesticus Linnaeus, 1758. The diploid number was established as 2n = 78 for both partridges, as well as for most species belonging to the Galliformes order, underlying the stability of chromosome number in avian karyotypes. Wide homologies were observed for macrochromosomes and gonosome except for chromosome 4, 7, 8 and Z which present differences in morphology and/or banding pattern. Neocentromere occurrence was suggested for both partridges chromosome 4 with an assumed paracentric inversion in the Chukar partridge chromosome 4. Terminal inversion in the long arm of the Barbary partridge chromosome Z was also found. These rearrangements confirm that the avian karyotypes structure is conserved interchromosomally, but not at the intrachromosomal scale.

Chicken Interferon-induced Protein with Tetratricopeptide Repeats 5 Antagonizes Replication of RNA Viruses

The intracellular actions of interferon (IFN)-regulated proteins, including IFN-induced proteins with tetratricopeptide repeats (IFITs), attribute a major component of the protective antiviral host defense. Here we applied genomics approaches to annotate the chicken IFIT locus and currently identified a single IFIT (chIFIT5) gene. The profound transcriptional level of this effector of innate immunity was mapped within its unique cis-acting elements. This highly virus- and IFN-responsive chIFIT5 protein interacted with negative sense viral RNA structures that carried a triphosphate group on its 5' terminus (ppp-RNA). This interaction reduced the replication of RNA viruses in lentivirus-mediated IFIT5-stable chicken fibroblasts whereas CRISPR/Cas9-edited chIFIT5 gene knockout fibroblasts supported the replication of RNA viruses. Finally, we generated mosaic transgenic chicken embryos stably expressing chIFIT5 protein or knocked-down for endogenous chIFIT5 gene. Replication kinetics of RNA viruses in these transgenic chicken embryos demonstrated the antiviral potential of chIFIT5 in ovo. Taken together, these findings propose that IFIT5 specifically antagonize RNA viruses by sequestering viral nucleic acids in chickens, which are unique in innate immune sensing and responses to viruses of both poultry and human health significance.

A survey of functional genomic variation in domesticated chickens

BACKGROUND

Deleterious genetic variation can increase in frequency as a result of mutations, genetic drift, and genetic hitchhiking. Although individual effects are often small, the cumulative effect of deleterious genetic variation can impact population fitness substantially. In this study, we examined the genome of commercial purebred chicken lines for deleterious and functional variations, combining genotype and whole-genome sequence data.

RESULTS

We analysed over 22,000 animals that were genotyped on a 60 K SNP chip from four purebred lines (two white egg and two brown egg layer lines) and two crossbred lines. We identified 79 haplotypes that showed a significant deficit in homozygous carriers. This deficit was assumed to stem from haplotypes that potentially harbour lethal recessive variations. To identify potentially deleterious mutations, a catalogue of over 10 million variants was derived from 250 whole-genome sequenced animals from three purebred white-egg layer lines. Out of 4219 putative deleterious variants, 152 mutations were identified that likely induce embryonic lethality in the homozygous state. Inferred deleterious variation showed evidence of purifying selection and deleterious alleles were generally overrepresented in regions of low recombination. Finally, we found evidence that mutations, which were inferred to be evolutionally intolerant, likely have positive effects in commercial chicken populations.

CONCLUSIONS

We present a comprehensive genomic perspective on deleterious and functional genetic variation in egg layer breeding lines, which are under intensive selection and characterized by a small effective population size. We show that deleterious variation is subject to purifying selection and that there is a positive relationship between recombination rate and purging efficiency. In addition, multiple putative functional coding variants were discovered in selective sweep regions, which are likely under positive selection. Together, this study provides a unique molecular perspective on functional and deleterious variation in commercial egg-laying chickens, which can enhance current genomic breeding practices to lower the frequency of undesirable variants in the population.

Transformation of animal genomics by next-generation sequencing technologies: a decade of challenges and their impact on genetic architecture

For more than a quarter of a century, sequencing technologies from Sanger's method to next-generation high-throughput techniques have provided fascinating opportunities in the life sciences. The continuing upward trajectory of sequencing technologies will improve livestock research and expedite the development of various new genomic and technological studies with farm animals. The use of high-throughput technologies in livestock research has increased interest in metagenomics, epigenetics, genome-wide association studies, and identification of single nucleotide polymorphisms and copy number variations. Such studies are beginning to provide revolutionary insights into biological and evolutionary processes. Farm animals, such as cattle, swine, and horses, have played a dual role as economically and agriculturally important animals as well as biomedical research models. The first part of this study explores the current state of sequencing methods, many of which are already used in animal genomic studies, and the second part summarizes the state of cattle, swine, horse, and chicken genome sequencing and illustrates its achievements during the last few years. Finally, we describe several high-throughput sequencing approaches for the improved detection of known, unknown, and emerging infectious agents, leading to better diagnosis of infectious diseases. The insights from viral metagenomics and the advancement of next-generation sequencing will strongly support specific and efficient vaccine development and provide strategies for controlling infectious disease transmission among animal populations and/or between animals and humans. However, prospective sequencing technologies will require further research and in-field testing before reaching the marketplace.

Genome-wide characterization of genetic variants and putative regions under selection in meat and egg-type chicken lines

Background

Meat and egg-type chickens have been selected for several generations for different traits. Artificial and natural selection for different phenotypes can change frequency of genetic variants, leaving particular genomic footprints throghtout the genome. Thus, the aims of this study were to sequence 28 chickens from two Brazilian lines (meat and white egg-type) and use this information to characterize genome-wide genetic variations, identify putative regions under selection using Fst method, and find putative pathways under selection.

Results

A total of 13.93 million SNPs and 1.36 million INDELs were identified, with more variants detected from the broiler (meat-type) line. Although most were located in non-coding regions, we identified 7255 intolerant non-synonymous SNPs, 512 stopgain/loss SNPs, 1381 frameshift and 1094 non-frameshift INDELs that may alter protein functions. Genes harboring intolerant non-synonymous SNPs affected metabolic pathways related mainly to reproduction and endocrine systems in the white-egg layer line, and lipid metabolism and metabolic diseases in the broiler line. Fst analysis in sliding windows, using SNPs and INDELs separately, identified over 300 putative regions of selection overlapping with more than 250 genes. For the first time in chicken, INDEL variants were considered for selection signature analysis, showing high level of correlation in results between SNP and INDEL data. The putative regions of selection signatures revealed interesting candidate genes and pathways related to important phenotypic traits in chicken, such as lipid metabolism, growth, reproduction, and cardiac development.

Conclusions

In this study, Fst method was applied to identify high confidence putative regions under selection, providing novel insights into selection footprints that can help elucidate the functional mechanisms underlying different phenotypic traits relevant to meat and egg-type chicken lines. In addition, we generated a large catalog of line-specific and common genetic variants from a Brazilian broiler and a white egg layer line that can be used for genomic studies involving association analysis with phenotypes of economic interest to the poultry industry.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4444-0) contains supplementary material, which is available to authorized users.

A dual transcript-discovery approach to improve the delimitation of gene features from RNA-seq data in the chicken model

The sequence of the chicken genome, like several other draft genome sequences, is presently not fully covered. Gaps, contigs assigned with low confidence and uncharacterized chromosomes result in gene fragmentation and imprecise gene annotation. Transcript abundance estimation from RNA sequencing (RNA-seq) data relies on read quality, library complexity and expression normalization. In addition, the quality of the genome sequence used to map sequencing reads, and the gene annotation that defines gene features, must also be taken into account. A partially covered genome sequence causes the loss of sequencing reads from the mapping step, while an inaccurate definition of gene features induces imprecise read counts from the assignment step. Both steps can significantly bias interpretation of RNA-seq data. Here, we describe a dual transcript-discovery approach combining a genome-guided gene prediction and a de novo transcriptome assembly. This dual approach enabled us to increase the assignment rate of RNA-seq data by nearly 20% as compared to when using only the chicken reference annotation, contributing therefore to a more accurate estimation of transcript abundance. More generally, this strategy could be applied to any organism with partial genome sequence and/or lacking a manually-curated reference annotation in order to improve the accuracy of gene expression studies.

Thyroid transcriptome analysis reveals different adaptive responses to cold environmental conditions between two chicken breeds

Selection for cold tolerance in chickens is important for improving production performance and animal welfare. The identification of chicken breeds with higher cold tolerance and production performance will help to target candidates for the selection. The thyroid gland plays important roles in thermal adaptation, and its function is influenced by breed differences and transcriptional plasticity, both of which remain largely unknown in the chicken thyroid transcriptome. In this study, we subjected Bashang Long-tail (BS) and Rhode Island Red (RIR) chickens to either cold or warm environments for 21 weeks and investigated egg production performance, body weight changes, serum thyroid hormone concentrations, and thyroid gland transcriptome profiles. RIR chickens had higher egg production than BS chickens under warm conditions, but BS chickens produced more eggs than RIRs under cold conditions. Furthermore, BS chickens showed stable body weight gain under cold conditions while RIRs did not. These results suggested that BS breed is a preferable candidate for cold-tolerance selection and that the cold adaptability of RIRs should be improved in the future. BS chickens had higher serum thyroid hormone concentrations than RIRs under both environments. RNA-Seq generated 344.3 million paired-end reads from 16 sequencing libraries, and about 90% of the processed reads were concordantly mapped to the chicken reference genome. Differential expression analysis identified 46-1,211 genes in the respective comparisons. With regard to breed differences in the thyroid transcriptome, BS chickens showed higher cell replication and development, and immune response-related activity, while RIR chickens showed higher carbohydrate and protein metabolism activity. The cold environment reduced breed differences in the thyroid transcriptome compared with the warm environment. Transcriptional plasticity analysis revealed different adaptive responses in BS and RIR chickens to cope with the cold, and showed higher responsiveness in BS compared with RIR chickens, suggesting greater adaptability of the thyroid in BS chickens. Moreover, 10,053 differential splicing events were revealed among the groups, with RNA splicing and processing, gene expression, transport, and metabolism being the main affected biological processes, identifying a valuable alternative splicing repertoire for the chicken thyroid. A short isoform of TPO (encoding thyroid peroxidase) containing multiple open reading frames was generated in both breeds by skipping exons 4 and 5 in the cold environment. These findings provide novel clues for future studies of the molecular mechanisms underlying cold adaptation and/or acclimation in chickens.

Genomic Region Containing Toll-Like Receptor Genes Has a Major Impact on Total IgM Antibodies Including KLH-Binding IgM Natural Antibodies in Chickens

Natural antibodies (NAb) are antigen binding antibodies present in individuals without a previous exposure to this antigen. Keyhole limpet hemocyanin (KLH)-binding NAb levels were previously associated with survival in chickens. This suggests that selective breeding for KLH-binding NAb may increase survival by means of improved general disease resistance. Genome-wide association studies (GWAS) were performed to identify genes underlying genetic variation in NAb levels. The studied population consisted of 1,628 adolescent layer chickens with observations for titers of KLH-binding NAb of the isotypes IgM, IgA, IgG, the total KLH-binding (IgT) NAb titers, total antibody concentrations of the isotypes IgM, IgA, IgG, and the total antibodies concentration in plasma. GWAS were performed using 57,636 single-nucleotide polymorphisms (SNP). One chromosomal region on chromosome 4 was associated with KLH-binding IgT NAb, and total IgM concentration, and especially with KLH-binding IgM NAb. The region of interest was fine mapped by imputing the region of the study population to whole genome sequence, and subsequently performing an association study using the imputed sequence variants. 16 candidate genes were identified, of which FAM114A1, Toll-like receptor 1 family member B (TLR1B), TLR1A, Krüppel-like factor 3 (KLF3) showed the strongest associations. SNP located in coding regions of the candidate genes were checked for predicted changes in protein functioning. One SNP (at 69,965,939 base pairs) received the maximum impact score from two independent prediction tools, which makes this SNP the most likely causal variant. This SNP is located in TLR1A, which suggests a fundamental role of TLR1A on regulation of IgM levels (i.e., KLH-binding IgM NAb, and total IgM concentration), or B cells biology, or both. This study contributes to increased understanding of (genetic) regulation of KLH-binding NAb levels, and total antibody concentrations.

Amphibian and Avian Karyotype Evolution: Insights from Lampbrush Chromosome Studies

Amphibian and bird karyotypes typically have a complex organization, which makes them difficult for standard cytogenetic analysis. That is, amphibian chromosomes are generally large, enriched with repetitive elements, and characterized by the absence of informative banding patterns. The majority of avian karyotypes comprise a small number of relatively large macrochromosomes and numerous tiny morphologically undistinguishable microchromosomes. A good progress in investigation of amphibian and avian chromosome evolution became possible with the usage of giant lampbrush chromosomes typical for growing oocytes. Due to the giant size, peculiarities of organization and enrichment with cytological markers, lampbrush chromosomes can serve as an opportune model for comprehensive high-resolution cytogenetic and cytological investigations. Here, we review the main findings on chromosome evolution in amphibians and birds that were obtained using lampbrush chromosomes. In particular, we discuss the data on evolutionary chromosomal rearrangements, accumulation of polymorphisms, evolution of sex chromosomes as well as chromosomal changes during clonal reproduction of interspecies hybrids.

New high copy tandem repeat in the content of the chicken W chromosome

The content of repetitive DNA in avian genomes is considerably less than in other investigated vertebrates. The first descriptions of tandem repeats were based on the results of routine biochemical and molecular biological experiments. Both satellite DNA and interspersed repetitive elements were annotated using library-based approach and de novo repeat identification in assembled genome. The development of deep-sequencing methods provides datasets of high quality without preassembly allowing one to annotate repetitive elements from unassembled part of genomes. In this work, we search the chicken assembly and annotate high copy number tandem repeats from unassembled short raw reads. Tandem repeat (GGAAA)_n has been identified and found to be the second after telomeric repeat (TTAGGG)n most abundant in the chicken genome. Furthermore, (GGAAA)_n repeat forms expanded arrays on the both arms of the chicken W chromosome. Our results highlight the complexity of repetitive sequences and update data about organization of sex W chromosome in chicken.

Systematic analysis of transcription start sites in avian development

Cap Analysis of Gene Expression (CAGE) in combination with single-molecule sequencing technology allows precision mapping of transcription start sites (TSSs) and genome-wide capture of promoter activities in differentiated and steady state cell populations. Much less is known about whether TSS profiling can characterize diverse and non-steady state cell populations, such as the approximately 400 transitory and heterogeneous cell types that arise during ontogeny of vertebrate animals. To gain such insight, we used the chick model and performed CAGE-based TSS analysis on embryonic samples covering the full 3-week developmental period. In total, 31,863 robust TSS peaks (>1 tag per million [TPM]) were mapped to the latest chicken genome assembly, of which 34% to 46% were active in any given developmental stage. ZENBU, a web-based, open-source platform, was used for interactive data exploration. TSSs of genes critical for lineage differentiation could be precisely mapped and their activities tracked throughout development, suggesting that non-steady state and heterogeneous cell populations are amenable to CAGE-based transcriptional analysis. Our study also uncovered a large set of extremely stable housekeeping TSSs and many novel stage-specific ones. We furthermore demonstrated that TSS mapping could expedite motif-based promoter analysis for regulatory modules associated with stage-specific and housekeeping genes. Finally, using Brachyury as an example, we provide evidence that precise TSS mapping in combination with Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR)-on technology enables us, for the first time, to efficiently target endogenous avian genes for transcriptional activation. Taken together, our results represent the first report of genome-wide TSS mapping in birds and the first systematic developmental TSS analysis in any amniote species (birds and mammals). By facilitating promoter-based molecular analysis and genetic manipulation, our work also underscores the value of avian models in unravelling the complex regulatory mechanism of cell lineage specification during amniote development.

The Agassiz's desert tortoise genome provides a resource for the conservation of a threatened species

Agassiz's desert tortoise (Gopherus agassizii) is a long-lived species native to the Mojave Desert and is listed as threatened under the US Endangered Species Act. To aid conservation efforts for preserving the genetic diversity of this species, we generated a whole genome reference sequence with an annotation based on deep transcriptome sequences of adult skeletal muscle, lung, brain, and blood. The draft genome assembly for G. agassizii has a scaffold N50 length of 252 kbp and a total length of 2.4 Gbp. Genome annotation reveals 20,172 protein-coding genes in the G. agassizii assembly, and that gene structure is more similar to chicken than other turtles. We provide a series of comparative analyses demonstrating (1) that turtles are among the slowest-evolving genome-enabled reptiles, (2) amino acid changes in genes controlling desert tortoise traits such as shell development, longevity and osmoregulation, and (3) fixed variants across the Gopherus species complex in genes related to desert adaptations, including circadian rhythm and innate immune response. This G. agassizii genome reference and annotation is the first such resource for any tortoise, and will serve as a foundation for future analysis of the genetic basis of adaptations to the desert environment, allow for investigation into genomic factors affecting tortoise health, disease and longevity, and serve as a valuable resource for additional studies in this species complex.

Validation of Alternative Transcript Splicing in Chicken Lines that Differ in Genetic Resistance to Marek's Disease

Utilizing RNA-seq data, 1,574 candidate genes with alternative splicing were previously identified between two chicken lines that differ in Marek's disease (MD) genetic resistance under control and Marek's disease virus infection conditions. After filtering out 1,530 genes with splice variants in the first or last exon, 44 genes were screened for possible exon loss or gain using PCR and gel electrophoresis. Consequently, 7 genes exhibited visually detectable differential expression of splice variants between lines 6 (MD resistant) and 7 (MD susceptible), and the resultant PCR products verified by DNA sequencing. Birds from inbred line 6 have transcripts that preferentially retain an exon compared to line 7 chickens for ITGB2, SGPL1, and COMMD5. Birds from inbred line 7 have alleles that preferentially retain an exon compared to line 6 for MOCS2. CCBL2 exon 1a is absent and ATAD1 exon 2 is truncated by 87 nucleotides in transcripts expressed by line 7 compared to those from line 6. For CHTF18, line 6 transcripts have an indel mutation with 7 additional nucleotides in exon 21 compared to line 7. The current study validates 7 genes with alternatively spliced isomers between the two chicken lines, which helps provide potential underlying mechanisms for the phenotypic differences.

Genetic Mapping of Quantitative Trait Loci for Egg Production and Egg Quality Traits in Chickens: a Review

Chickens display a wide spectrum of phenotypic variations in quantitative traits such as egg-related traits. Quantitative trait locus (QTL) analysis is a statistical method used to understand the relationship between phenotypic (trait measurements) and genotypic data (molecular markers). We have performed QTL analyses for egg-related traits using an original resource population based on the Japanese Large Game (Oh-Shamo) and the White Leghorn breeds of chickens. In this article, we summarize the results of our extensive QTL analyses for 11 and 66 traits for egg production and egg quality, respectively. We reveal that at least 30 QTL regions on 17 different chromosomes affect phenotypic variation in egg-related traits. Each locus had an age-specific effect on traits, and a variety in effects was also apparent, such as additive, dominance, and epistatic-interaction effects. Although genome-wide association study (GWAS) is suitable for gene-level resolution mapping of GWAS loci with additive effects, QTL mapping studies enable us to comprehensively understand genetic control, such as chromosomal regions, genetic contribution to phenotypic variance, mode of inheritance, and age-specificity of both common and rare alleles. QTL analyses also describe the relationship between genotypes and phenotypes in experimental populations. Accumulation of QTL information, including GWAS loci, is also useful for studies of population genomics approached without phenotypic data in order to validate the identified genomic signatures of positive selection. The combination of QTL studies and next-generation sequencing techniques with uncharacterized genetic resources will enhance current understanding of the relationship between genotypes and phenotypes in livestock animals.

A New Chicken Genome Assembly Provides Insight into Avian Genome Structure

The importance of the Gallus gallus (chicken) as a model organism and agricultural animal merits a continuation of sequence assembly improvement efforts. We present a new version of the chicken genome assembly (Gallus_gallus-5.0; GCA_000002315.3), built from combined long single molecule sequencing technology, finished BACs, and improved physical maps. In overall assembled bases, we see a gain of 183 Mb, including 16.4 Mb in placed chromosomes with a corresponding gain in the percentage of intact repeat elements characterized. Of the 1.21 Gb genome, we include three previously missing autosomes, GGA30, 31, and 33, and improve sequence contig length 10-fold over the previous Gallus_gallus-4.0. Despite the significant base representation improvements made, 138 Mb of sequence is not yet located to chromosomes. When annotated for gene content, Gallus_gallus-5.0 shows an increase of 4679 annotated genes (2768 noncoding and 1911 protein-coding) over those in Gallus_gallus-4.0. We also revisited the question of what genes are missing in the avian lineage, as assessed by the highest quality avian genome assembly to date, and found that a large fraction of the original set of missing genes are still absent in sequenced bird species. Finally, our new data support a detailed map of MHC-B, encompassing two segments: one with a highly stable gene copy number and another in which the gene copy number is highly variable. The chicken model has been a critical resource for many other fields of study, and this new reference assembly will substantially further these efforts.