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.

Comparison of the somatic TADs and lampbrush chromomere-loop complexes in transcriptionally active prophase I oocytes.

In diplotene oocyte nuclei of all vertebrate species, except mammals, chromosomes lack interchromosomal contacts and chromatin is linearly compartmentalized into distinct chromomere-loop complexes forming lampbrush chromosomes. However, the mechanisms underlying the formation of chromomere-loop complexes remain unexplored. Here we aimed to compare somatic topologically associating domains (TADs), recently identified in chicken embryonic fibroblasts, with chromomere-loop complexes in lampbrush meiotic chromosomes. By measuring 3D-distances and colocalization between linear equidistantly located genomic loci, positioned within one TAD or separated by a TAD border, we confirmed the presence of predicted TADs in chicken embryonic fibroblast nuclei. Using three-colored FISH with BAC probes, we mapped equidistant genomic regions included in several sequential somatic TADs on isolated chicken lampbrush chromosomes. Eight genomic regions, each comprising two or three somatic TADs, were mapped to non-overlapping neighboring lampbrush chromatin domains - lateral loops, chromomeres, or chromomere-loop complexes. Genomic loci from the neighboring somatic TADs could localize in one lampbrush chromomere-loop complex, while genomic loci belonging to the same somatic TAD could be localized in neighboring lampbrush chromomere-loop domains. In addition, FISH-mapping of BAC probes to the nascent transcripts on the lateral loops indicates transcription of at least 17 protein-coding genes and 2 non-coding RNA genes during the lampbrush stage of chicken oogenesis, including genes involved in oocyte maturation and early embryo development.

Integrating Omics and CRISPR Technology for Identification and Verification of Genomic Safe Harbor Loci in the Chicken Genome.

BACKGROUND: One of the most prominent questions in the field of transgenesis is 'Where in the genome to integrate a transgene?'. Escape from epigenetic silencing and promoter shutdown of the transgene needs reliable genomic safe harbor (GSH) loci. Advances in genome engineering technologies combined with multi-omics bioinformatics data have enabled rational evaluation of GSH loci in the host genome. Currently, no validated GSH loci have been evaluated in the chicken genome. RESULTS: Here, we analyzed and experimentally examined two GSH loci in the genome of chicken cells. To this end, putative GSH loci including chicken HIPP-like (cHIPP; between DRG1 and EIF4ENIF1 genes) and chicken ROSA-like (cROSA; upstream of the THUMPD3 gene) were predicted using multi-omics bioinformatics data. Then, the durable expression of the transgene was validated by experimental characterization of continuously-cultured isogenous cell clones harboring DsRed2-ΔCMV-EGFP cassette in the predicted loci. The weakened form of the CMV promoter (ΔCMV) allowed the precise evaluation of GSH loci in a locus-dependent manner compared to the full-length CMV promoter. CONCLUSIONS: cHIPP and cROSA loci introduced in this study can be reliably exploited for consistent bio-manufacturing of recombinant proteins in the genetically-engineered chickens. Also, results showed that the genomic context dictates the expression of transgene controlled by ΔCMV in GSH loci.

Expanding duplication of the testis PHD Finger Protein 7 (PHF7) gene in the chicken genome.

Gene duplications increase genetic and phenotypic diversity and occur in complex genomic regions that are still difficult to sequence and assemble. PHD Finger Protein 7 (PHF7) acts during spermiogenesis for histone-to-histone protamine exchange and is a determinant of male fertility in Drosophila and the mouse. We aimed to explore and characterise in the chicken genome the expanding family of the numerous orthologues of the unique mouse Phf7 gene (highly expressed in the testis), observing the fact that this information is unclear and/or variable according to the versions of databases. We validated nine primer pairs by in silico PCR for their use in screening the chicken bacterial artificial chromosome (BAC) library to produce BAC-derived probes to detect and localise PHF7-like loci by fluorescence in situ hybridisation (FISH). We selected nine BAC that highlighted nine chromosomal regions for a total of 10 distinct PHF7-like loci on five Gallus gallus chromosomes: Chr1 (three loci), Chr2 (two loci), Chr12 (one locus), Chr19 (one locus) and ChrZ (three loci). We sequenced the corresponding BAC by using high-performance PacBio technology. After assembly, we performed annotation with the FGENESH program: there were a total of 116 peptides, including 39 PHF7-like proteins identified by BLASTP. These proteins share a common exon-intron core structure of 8-11 exons. Phylogeny revealed that the duplications occurred first between chromosomal regions and then inside each region. There are other duplicated genes in the identified BAC sequences, suggesting that these genomic regions exhibit a high rate of tandem duplication. We showed that the PHF7 gene, which is highly expressed in the rooster testis, is a highly duplicated gene family in the chicken genome, and this phenomenon probably concerns other bird species.

Single-Cell DNA Methylation Analysis of Chicken Lampbrush Chromosomes.

DNA methylation is an essential epigenetic regulation mechanism implicated in transcription and replication control, developmental reprogramming, retroelements silencing and other genomic processes. During mammalian development, a specific DNA methylation pattern should be established in germ cells to allow embryonic development. Less is known about germ cell DNA methylation in other species. To close this gap, we performed a single-cell methylome analysis of chicken diplotene oocytes. We comprehensively characterized methylation patterns in these cells, obtained methylation-based chicken genome segmentation and identified oocyte-specific methylated gene promoters. Our data show that despite the formation of specific transcriptionally hyperactive genome architecture in chicken diplotene oocytes, methylation patterns in these cells closely resemble genomic distribution observed in somatic tissues.

Vaccination and Host Marek's Disease-Resistance Genotype Significantly Reduce Oncogenic Gallid alphaherpesvirus 2 Telomere Integration in Host Birds.

Marek's disease (MD) is an infectious disease characterized by lymphomas and high mortality in susceptible chickens. The causative and ubiquitous alpha-herpesvirus known as MD virus (MDV) integrates into host telomeres during early infection through latency, known to be an important phase for oncogenic transformation. Herein, we sought to determine the influence of vaccination and host genetics on the temporal dynamics of MDV-host genome interactions. We studied integration profiles using 2 MD vaccines that vary in protective efficacy in 2 genetic lines that differ in MD resistance/susceptibility. Virus integration of both oncogenic MDV and vaccine strains was observed in both MD susceptible and resistant birds, however, the lines differed in their dynamic telomere-integration profiles. Notably, the resistant host genotype exhibited a smaller percentage of replicating cells with the virus telomere-integrated only phenotype as compared to the susceptible genotype. Vaccination with Rispens, the most protective MD vaccine, also reduced the establishment of the virus telomere-integrated only phenotype, suggesting a significant role of the phenotype in MD lymphoma development. The effect of Rispens vaccination was most dramatic in the susceptible genotype. These results suggest important connections between vaccinal immunity, MDV telomere integration, virus-induced oncogenesis, and virus-host genome interactions in the context of host genetics and disease susceptibility.

Three-dimensional architecture of tandem repeats in chicken interphase nucleus.

Tandem repeats belong to a class of genomic repetitive elements that form arrays of head-to-tail monomers. Due to technical difficulties in sequencing and assembly of large tandem repeat arrays, it remains largely unknown by which mechanisms tandem-repeat-containing regions aid in maintenance of ordered radial genome organization during interphase. Here we analyzed spatial distribution of several types of tandem repeats in interphase nuclei of chicken MDCC-MSB1 cells and somatic tissues relative to heterochromatin compartments and nuclear center. We showed that telomere and subtelomere repeats generally localize at the nuclear or chromocenters periphery. A tandem repeat known as CNM, typical for centromere regions of gene-dense microchromosomes, forms interchromosome clusters and occupies DAPI-positive chromocenters that appear predominantly within the nuclear interior. In contrast, centromere-specific tandem repeats of the majority of gene-poor macrochromosomes are embedded into the peripheral layer of heterochromatin. Chicken chromocenters rarely comprise centromere sequences of both macro- and microchromosomes, whose territories localize in different radial nuclear zones. Possible mechanisms of observed tandem repeats positioning and its implication in highly ordered arrangement of chromosome territories in chicken interphase nucleus are discussed.

Sixteen additional enhancers associated with the chicken Sox2 locus outside the central 50-kb region.

The transcription factor Sox2 plays a central role in the regulation of neuro-sensory development, and many other developmental processes. To gain an in depth understanding of the Sox2 gene regulation, we previously investigated the Sox2-proximal 50-kb region of the chicken genome to determine enhancers based on functional assays using chicken embryo electroporation. We identified 11 enhancers with specificity for neuro-sensory tissues. In this study, we extended the analysis of Sox2 locus-associated enhancers to a 200-kb region and identified 16 additional enhancers with functions in neuro-sensory development. These enhancers roughly correspond to a fraction of the sequence blocks that are highly conserved between chicken and mammalian genomes. The neural enhancers were activated in sequence, thereby creating a complex pattern of functional overlaps in the developing central nervous system (CNS). The variations in the specificities of the sensory enhancers also reflected the intermediate steps of sensory tissue development. This study provides an example where a single transcription factor gene has numerous regulatory elements that allow it to fulfill many functional roles in different biological contexts.

Large-scale genome-wide SNP analysis reveals the rugged (and ragged) landscape of global ancestry, phylogeny, and demographic history in chicken breeds. / å¤§è§„æ¨¡å ¨åŸºå› ç»„SNPåˆ†æžæ­ç¤ºäº†é¸¡å“ç§çš„å ¨çƒç¥–å ˆã€ç§ç¾¤å‘å±•å’Œç§ç¾¤åŽ†å²çš„å¤æ‚(å’Œå¤šæ ·)çš„é—ä¼ å›¾è°±.

The worldwide chicken gene pool encompasses a remarkable, but shrinking, number of divergently selected breeds of diverse origin. This study was a large-scale genome-wide analysis of the landscape of the complex molecular architecture, genetic variability, and detailed structure among 49 populations. These populations represent a significant sample of the world's chicken breeds from Europe (Russia, Czech Republic, France, Spain, UK, etc.), Asia (China), North America (USA), and Oceania (Australia). Based on the results of breed genotyping using the Illumina 60K single nucleotide polymorphism (SNP) chip, a bioinformatic analysis was carried out. This included the calculation of heterozygosity/homozygosity statistics, inbreeding coefficients, and effective population size. It also included assessment of linkage disequilibrium and construction of phylogenetic trees. Using multidimensional scaling, principal component analysis, and ADMIXTURE-assisted global ancestry analysis, we explored the genetic structure of populations and subpopulations in each breed. An overall 49-population phylogeny analysis was also performed, and a refined evolutionary model of chicken breed formation was proposed, which included egg, meat, dual-purpose types, and ambiguous breeds. Such a large-scale survey of genetic resources in poultry farming using modern genomic methods is of great interest both from the viewpoint of a general understanding of the genetics of the domestic chicken and for the further development of genomic technologies and approaches in poultry breeding. In general, whole genome SNP genotyping of promising chicken breeds from the worldwide gene pool will promote the further development of modern genomic science as applied to poultry.

Assignment of the somatic A/B compartments to chromatin domains in giant transcriptionally active lampbrush chromosomes.

BACKGROUND: The three-dimensional configuration of the eukaryotic genome is an emerging area of research. Chromosome conformation capture outlined genome segregation into large scale A and B compartments corresponding mainly to transcriptionally active and repressive chromatin. It remains unknown how the compartmentalization of the genome changes in growing oocytes of animals with hypertranscriptional type of oogenesis. Such oocytes are characterized by highly elongated chromosomes, called lampbrush chromosomes, which acquire a typical chromomere-loop appearance, representing one of the classical model systems for exploring the structural and functional organization of chromatin domains. RESULTS: Here, we compared the distribution of A/B compartments in chicken somatic cells with chromatin domains in lampbrush chromosomes. We found that in lampbrush chromosomes, the extended chromatin domains, restricted by compartment boundaries in somatic cells, disintegrate into individual chromomeres. Next, we performed FISH-mapping of the genomic loci, which belong to A or B chromatin compartments as well as to A/B compartment transition regions in embryonic fibroblasts on isolated lampbrush chromosomes. We found, that in chicken lampbrush chromosomes, clusters of dense compact chromomeres bearing short lateral loops and enriched with repressive epigenetic modifications generally correspond to constitutive B compartments in somatic cells. A compartments align with lampbrush chromosome segments with smaller, less compact chromomeres, longer lateral loops, and a higher transcriptional status. Clusters of small loose chromomeres with relatively long lateral loops show no obvious correspondence with either A or B compartment identity. Some genes belonging to facultative B (sub-) compartments can be tissue-specifically transcribed during oogenesis, forming distinct lateral loops. CONCLUSIONS: Here, we established a correspondence between the A/B compartments in somatic interphase nucleus and chromatin segments in giant lampbrush chromosomes from diplotene stage oocytes. The chromomere-loop structure of the genomic regions corresponding to interphase A and B compartments reveals the difference in how they are organized at the level of chromatin domains. The results obtained also suggest that gene-poor regions tend to be packed into chromomeres.

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.

The genome of the chicken DT40 bursal lymphoma cell line.

The chicken DT40 cell line is a widely used model system in the study of multiple cellular processes due to the efficiency of homologous gene targeting. The cell line was derived from a bursal lymphoma induced by avian leukosis virus infection. In this study we characterized the genome of the cell line using whole genome shotgun sequencing and single nucleotide polymorphism array hybridization. The results indicate that wild-type DT40 has a relatively normal karyotype, except for whole chromosome copy number gains, and no karyotype variability within stocks. In a comparison to two domestic chicken genomes and the Gallus gallus reference genome, we found no unique mutational processes shaping the DT40 genome except for a mild increase in insertion and deletion events, particularly deletions at tandem repeats. We mapped coding sequence mutations that are unique to the DT40 genome; mutations inactivating the PIK3R1 and ATRX genes likely contributed to the oncogenic transformation. In addition to a known avian leukosis virus integration in the MYC gene, we detected further integration sites that are likely to de-regulate gene expression. The new findings support the hypothesis that DT40 is a typical transformed cell line with a relatively intact genome; therefore, it is well-suited to the role of a model system for DNA repair and related processes. The sequence data generated by this study, including a searchable de novo genome assembly and annotated lists of mutated genes, will support future research using this cell line.

Calibration of mutation rates reveals diverse subfamily structure of galliform CR1 repeats.

Chicken Repeat 1 (CR1) repeats are the most abundant family of repeats in the chicken genome, with more than 200,000 copies accounting for approximately 80% of the chicken interspersed repeats. CR1 repeats are believed to have arisen from the retrotransposition of a small number of master elements, which gave rise to the 22 CR1 subfamilies as previously reported in Repbase. We performed a global assessment of the divergence distributions, phylogenies, and consensus sequences of CR1 repeats in the chicken genome. We identified and validated 57 chicken CR1 subfamilies and further analyzed the correlation between these subfamilies and their regional GC contents. We also discovered one novel lineage-specific CR1 subfamilies in turkeys when compared with chickens. We built an evolutionary tree of these subfamilies and concluded that CR1 repeats may play an important role in reshaping the structure of bird genomes.