Genetic variation exists for telomeric array organization within and among the genomes of normal, immortalized, and transformed chicken systems

Hi-C analysis of genomic contacts revealed karyotype abnormalities in chicken HD3 cell line

BACKGROUND

Karyotype abnormalities are frequent in immortalized continuous cell lines either transformed or derived from primary tumors. Chromosomal rearrangements can cause dramatic changes in gene expression and affect cellular phenotype and behavior during in vitro culture. Structural variations of chromosomes in many continuous mammalian cell lines are well documented, but chromosome aberrations in cell lines from other vertebrate models often remain understudied. The chicken LSCC-HD3 cell line (HD3), generated from erythroid precursors, was used as an avian model for erythroid differentiation and lineage-specific gene expression. However, karyotype abnormalities in the HD3 cell line were not assessed. In the present study, we applied high-throughput chromosome conformation capture to analyze 3D genome organization and to detect chromosome rearrangements in the HD3 cell line.

RESULTS

We obtained Hi-C maps of genomic interactions for the HD3 cell line and compared A/B compartments and topologically associating domains between HD3 and several other cell types. By analysis of contact patterns in the Hi-C maps of HD3 cells, we identified more than 25 interchromosomal translocations of regions ≥ 200 kb on both micro- and macrochromosomes. We classified most of the observed translocations as unbalanced, leading to the formation of heteromorphic chromosomes. In many cases of microchromosome rearrangements, an entire microchromosome together with other macro- and microchromosomes participated in the emergence of a derivative chromosome, resembling "chromosomal fusions'' between acrocentric microchromosomes. Intrachromosomal inversions, deletions and duplications were also detected in HD3 cells. Several of the identified simple and complex chromosomal rearrangements, such as between GGA2 and GGA1qter; GGA5, GGA4p and GGA7p; GGA4q, GGA6 and GGA19; and duplication of the sex chromosome GGAW, were confirmed by FISH.

CONCLUSIONS

In the erythroid progenitor HD3 cell line, in contrast to mature and immature erythrocytes, the genome is organized into distinct topologically associating domains. The HD3 cell line has a severely rearranged karyotype with most of the chromosomes engaged in translocations and can be used in studies of genome structure-function relationships. Hi-C proved to be a reliable tool for simultaneous assessment of the spatial genome organization and chromosomal aberrations in karyotypes of birds with a large number of microchromosomes.

Inbreeding is associated with shorter early-life telomere length in a wild passerine

Inbreeding can have negative effects on survival and reproduction, which may be of conservation concern in small and isolated populations. However, the physiological mechanisms underlying inbreeding depression are not well-known. The length of telomeres, the DNA sequences protecting chromosome ends, has been associated with health or fitness in several species. We investigated effects of inbreeding on early-life telomere length in two small island populations of wild house sparrows (Passer domesticus) known to be affected by inbreeding depression. Using genomic measures of inbreeding we found that inbred nestling house sparrows (n = 371) have significantly shorter telomeres. Using pedigree-based estimates of inbreeding we found a tendency for inbred nestling house sparrows to have shorter telomeres (n = 1195). This negative effect of inbreeding on telomere length may have been complemented by a heterosis effect resulting in longer telomeres in individuals that were less inbred than the population average. Furthermore, we found some evidence of stronger effects of inbreeding on telomere length in males than females. Thus, telomere length may reveal subtle costs of inbreeding in the wild and demonstrate a route by which inbreeding negatively impacts the physiological state of an organism already at early life-history stages.

Marek's Disease Virus Telomeric Integration Profiles of Neoplastic Host Tissues Reveal Unbiased Chromosomal Selection and Loss of Cellular Diversity during Tumorigenesis

The avian α-herpesvirus known as Marek's disease virus (MDV) linearly integrates its genomic DNA into host telomeres during infection. The resulting disease, Marek's disease (MD), is characterized by virally-induced lymphomas with high mortality. The temporal dynamics of MDV-positive (MDV+) transformed cells and expansion of MD lymphomas remain targets for further understanding. It also remains to be determined whether specific host chromosomal sites of MDV telomere integration confer an advantage to MDV-transformed cells during tumorigenesis. We applied MDV-specific fluorescence in situ hybridization (MDV FISH) to investigate virus-host cytogenomic interactions within and among a total of 37 gonad lymphomas and neoplastic splenic samples in birds infected with virulent MDV. We also determined single-cell, chromosome-specific MDV integration profiles within and among transformed tissue samples, including multiple samples from the same bird. Most mitotically-dividing cells within neoplastic samples had the cytogenomic phenotype of 'MDV telomere-integrated only', and tissue-specific, temporal changes in phenotype frequencies were detected. Transformed cell populations composing gonad lymphomas exhibited significantly lower diversity, in terms of heterogeneity of MDV integration profiles, at the latest stages of tumorigenesis (>50 days post-infection (dpi)). We further report high interindividual and lower intraindividual variation in MDV integration profiles of lymphoma cells. There was no evidence of integration hotspots into a specific host chromosome(s). Collectively, our data suggests that very few transformed MDV+ T cell populations present earlier in MDV-induced lymphomas (32-50 dpi), survive, and expand to become the dominant clonal population in more advanced MD lymphomas (51-62 dpi) and establish metastatic lymphomas.

One-step generation of a targeted knock-in calf using the CRISPR-Cas9 system in bovine zygotes

Background

The homologous recombination (HR) pathway is largely inactive in early embryos prior to the first cell division, making it difficult to achieve targeted gene knock-ins. The homology-mediated end joining (HMEJ)-based strategy has been shown to increase knock-in efficiency relative to HR, non-homologous end joining (NHEJ), and microhomology-mediated end joining (MMEJ) strategies in non-dividing cells.

Results

By introducing gRNA/Cas9 ribonucleoprotein complex and a HMEJ-based donor template with 1 kb homology arms flanked by the H11 safe harbor locus gRNA target site, knock-in rates of 40% of a 5.1 kb bovine sex-determining region Y (SRY)-green fluorescent protein (GFP) template were achieved in Bos taurus zygotes. Embryos that developed to the blastocyst stage were screened for GFP, and nine were transferred to recipient cows resulting in a live phenotypically normal bull calf. Genomic analyses revealed no wildtype sequence at the H11 target site, but rather a 26 bp insertion allele, and a complex 38 kb knock-in allele with seven copies of the SRY-GFP template and a single copy of the donor plasmid backbone. An additional minor 18 kb allele was detected that looks to be a derivative of the 38 kb allele resulting from the deletion of an inverted repeat of four copies of the SRY-GFP template.

Conclusion

The allelic heterogeneity in this biallelic knock-in calf appears to have resulted from a combination of homology directed repair, homology independent targeted insertion by blunt-end ligation, NHEJ, and rearrangement following editing of the gRNA target site in the donor template.

This study illustrates the potential to produce targeted gene knock-in animals by direct cytoplasmic injection of bovine embryos with gRNA/Cas9, although further optimization is required to ensure a precise single-copy gene integration event.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07418-3.

Avian Expression Patterns and Genomic Mapping Implicate Leptin in Digestion and TNF in Immunity, Suggesting That Their Interacting Adipokine Role Has Been Acquired Only in Mammals

In mammals, leptin and tumor-necrosis factor (TNF) are prominent interacting adipokines mediating appetite control and insulin sensitivity. While TNF pleiotropically functions in immune defense and cell survival, leptin is largely confined to signaling energy stores in adipocytes. Knowledge about the function of avian leptin and TNF is limited and they are absent or lowly expressed in adipose, respectively. Employing radiation-hybrid mapping and FISH-TSA, we mapped TNF and its syntenic genes to chicken chromosome 16 within the major histocompatibility complex (MHC) region. This mapping position suggests that avian TNF has a role in regulating immune response. To test its possible interaction with leptin within the immune system and beyond, we compared the transcription patterns of TNF, leptin and their cognate receptors obtained by meta-analysis of GenBank RNA-seq data. While expression of leptin and its receptor (LEPR) were detected in the brain and digestive tract, TNF and its receptor mRNAs were primarily found in viral-infected and LPS-treated leukocytes. We confirmed leptin expression in the duodenum by immunohistochemistry staining. Altogether, we suggest that whereas leptin and TNF interact as adipokines in mammals, in birds, they have distinct roles. Thus, the interaction between leptin and TNF may be unique to mammals.

Deep landscape update of dispersed and tandem repeats in the genome model of the red jungle fowl, Gallus gallus, using a series of de novo investigating tools

BACKGROUND

The program RepeatMasker and the database Repbase-ISB are part of the most widely used strategy for annotating repeats in animal genomes. They have been used to show that avian genomes have a lower repeat content (8-12 %) than the sequenced genomes of many vertebrate species (30-55 %). However, the efficiency of such a library-based strategies is dependent on the quality and completeness of the sequences in the database that is used. An alternative to these library based methods are methods that identify repeats de novo. These alternative methods have existed for a least a decade and may be more powerful than the library based methods. We have used an annotation strategy involving several complementary de novo tools to determine the repeat content of the model genome galGal4 (1.04 Gbp), including identifying simple sequence repeats (SSRs), tandem repeats and transposable elements (TEs).

RESULTS

We annotated over one Gbp. of the galGal4 genome and showed that it is composed of approximately 19 % SSRs and TEs repeats. Furthermore, we estimate that the actual genome of the red jungle fowl contains about 31-35 % repeats. We find that library-based methods tend to overestimate TE diversity. These results have a major impact on the current understanding of repeats distributions throughout chromosomes in the red jungle fowl.

CONCLUSIONS

Our results are a proof of concept of the reliability of using de novo tools to annotate repeats in large animal genomes. They have also revealed issues that will need to be resolved in order to develop gold-standard methodologies for annotating repeats in eukaryote genomes.

Chromosomal-Level Assembly of the Asian Seabass Genome Using Long Sequence Reads and Multi-layered Scaffolding

We report here the ~670 Mb genome assembly of the Asian seabass (Lates calcarifer), a tropical marine teleost. We used long-read sequencing augmented by transcriptomics, optical and genetic mapping along with shared synteny from closely related fish species to derive a chromosome-level assembly with a contig N50 size over 1 Mb and scaffold N50 size over 25 Mb that span ~90% of the genome. The population structure of L. calcarifer species complex was analyzed by re-sequencing 61 individuals representing various regions across the species' native range. SNP analyses identified high levels of genetic diversity and confirmed earlier indications of a population stratification comprising three clades with signs of admixture apparent in the South-East Asian population. The quality of the Asian seabass genome assembly far exceeds that of any other fish species, and will serve as a new standard for fish genomics.

Mother-offspring and nest-mate resemblance but no heritability in early-life telomere length in white-throated dippers

Telomeres are protective DNA-protein complexes located at the ends of eukaryotic chromosomes, whose length has been shown to predict life-history parameters in various species. Although this suggests that telomere length is subject to natural selection, its evolutionary dynamics crucially depends on its heritability. Using pedigree data for a population of white-throated dippers (Cinclus cinclus), we test whether and how variation in early-life relative telomere length (RTL, measured as the number of telomeric repeats relative to a control gene using qPCR) is transmitted across generations. We disentangle the relative effects of genes and environment and test for sex-specific patterns of inheritance. There was strong and significant resemblance among offspring sharing the same nest and offspring of the same cohort. Furthermore, although offspring resemble their mother, and there is some indication for an effect of inbreeding, additive genetic variance and heritability are close to zero. We find no evidence for a role of either maternal imprinting or Z-linked inheritance in generating these patterns, suggesting they are due to non-genetic maternal and common environment effects instead. We conclude that in this wild bird population, environmental factors are the main drivers of variation in early-life RTL, which will severely bias estimates of heritability when not modelled explicitly.

Comparative cytogenomics of poultry: mapping of single gene and repeat loci in the Japanese quail (Coturnix japonica)

Well-characterized molecular and cytogenetic maps are yet to be established in Japanese quail (Coturnix japonica). The aim of the current study was to cytogenetically map and determine linkage of specific genes and gene complexes in Japanese quail through the use of chicken (Gallus gallus) and turkey (Meleagris gallopavo) genomic DNA probes and conduct a comparative study among the three genomes. Chicken and turkey clones were used as probes on mitotic metaphase and meiotic pachytene stage chromosomes of the three species for the purpose of high-resolution fluorescence in situ hybridization (FISH). The genes and complexes studied included telomerase RNA (TR), telomerase reverse transcriptase (TERT), 5S rDNA, 18S-5.8S-28S rDNA (i.e., nucleolus organizer region (NOR)), and the major histocompatibility complex (MHC). The telomeric profile of Japanese quail was investigated through the use of FISH with a TTAGGG-PNA probe. A range of telomeric array sizes were confirmed as found for the other poultry species. Three NOR loci were identified in Japanese quail, and single loci each for TR, TERT, 5S rDNA and the MHC-B. The MHC-B and one NOR locus were linked on a microchromosome in Japanese quail. We confirmed physical linkage of 5S rDNA and the TR gene on an intermediate-sized chromosome in quail, similar to both chicken and turkey. TERT localized to CJA 2 in quail and the orthologous chromosome region in chicken (GGA 2) and in turkey (MGA 3). The cytogenetic profile of Japanese quail was further developed by this study and synteny was identified among the three poultry species.

[Telomere biology of the chicken]

Telomeres, the nucleoprotein "caps" protecting the ends of linear chromosomes, are maintained by telomerase. Telomeres have important roles in maintaining genomic stability and preventing senescence or oncogenesis. Chicken is a classical model animal for genetic and developmental studies. With further development of chicken genomics, great progress has been made in research of chicken telomere and telomerase. This review describes recent advances and future research directions in chicken telomere biology.

The expression of preaxial polydactyly is influenced by modifying genetic elements and is not maintained by chromosomal inversion in an avian biomedical model

Polydactyly (Po) is a common mutation found in many vertebrates. The UCD-Po.003 congenic chicken line was previously characterized for Po inheritance (autosomal dominant) and the mutation was mapped to chromosome 2p. Here, we describe for the first time the range and variability of the phenotype in this congenic line. Further, we studied the hypothesis that a chromosomal inversion was responsible for the maintenance of a large (6.3 Mb) candidate gene region. Fluorescence in situ hybridization employing BACs encompassing a 10.7-Mb region of GGA2p showed that the Po chromosome was normal, i.e. exhibits the wild-type BAC order. Continued fine-mapping along with a change in breeding strategy reduced the size of the causative region to 1.43 Mb. Recent research indicates that the cause of preaxial Po resides within a 794-bp highly conserved zone of polarizing activity regulatory sequence (ZRS) element located in intron 5 of the LMBR1 gene; however, the ZRS polymorphism of interest is found in some but not all breeds of polydactylous chicken. Therefore, we sequenced the ZRS in 101 heterozygous and 30 unaffected (wild-type) individuals to establish the relevance of this region to the Po condition in the UCD-Po.003 congenic line. A single point mutation (C/A at coordinate GGA2p: 8,414,121) within the ZRS segregated with carrier status. The polydactylous UCD-Silkie line also maintains this SNP in addition to a single base deletion. An inheritance analysis of the phenotypic variation in UCD-Po.003 suggests recessive epistasis as the mode of inheritance for the additional modifying genetic elements, residing outside the ZRS, to impact the preaxial polydactyl phenotype. These results contribute to our understanding of the cause of Po in an important vertebrate model.

Genome-wide differential gene expression in immortalized DF-1 chicken embryo fibroblast cell line

BACKGROUND

When compared to primary chicken embryo fibroblast (CEF) cells, the immortal DF-1 CEF line exhibits enhanced growth rates and susceptibility to oxidative stress. Although genes responsible for cell cycle regulation and antioxidant functions have been identified, the genome-wide transcription profile of immortal DF-1 CEF cells has not been previously reported. Global gene expression in primary CEF and DF-1 cells was performed using a 4X44K chicken oligo microarray.

RESULTS

A total of 3876 differentially expressed genes were identified with a 2 fold level cutoff that included 1706 up-regulated and 2170 down-regulated genes in DF-1 cells. Network and functional analyses using Ingenuity Pathways Analysis (IPA, Ingenuity® Systems, http://www.ingenuity.com) revealed that 902 of 3876 differentially expressed genes were classified into a number of functional groups including cellular growth and proliferation, cell cycle, cellular movement, cancer, genetic disorders, and cell death. Also, the top 5 gene networks with intermolecular connections were identified. Bioinformatic analyses suggested that DF-1 cells were characterized by enhanced molecular mechanisms for cell cycle progression and proliferation, suppressing cell death pathways, altered cellular morphogenesis, and accelerated capacity for molecule transport. Key molecules for these functions include E2F1, BRCA1, SRC, CASP3, and the peroxidases.

CONCLUSIONS

The global gene expression profiles provide insight into the cellular mechanisms that regulate the unique characteristics observed in immortal DF-1 CEF cells.

A comparative physical map reveals the pattern of chromosomal evolution between the turkey (Meleagris gallopavo) and chicken (Gallus gallus) genomes

Background

A robust bacterial artificial chromosome (BAC)-based physical map is essential for many aspects of genomics research, including an understanding of chromosome evolution, high-resolution genome mapping, marker-assisted breeding, positional cloning of genes, and quantitative trait analysis. To facilitate turkey genetics research and better understand avian genome evolution, a BAC-based integrated physical, genetic, and comparative map was developed for this important agricultural species.

Results

The turkey genome physical map was constructed based on 74,013 BAC fingerprints (11.9 × coverage) from two independent libraries, and it was integrated with the turkey genetic map and chicken genome sequence using over 41,400 BAC assignments identified by 3,499 overgo hybridization probes along with > 43,000 BAC end sequences. The physical-comparative map consists of 74 BAC contigs, with an average contig size of 13.6 Mb. All but four of the turkey chromosomes were spanned on this map by three or fewer contigs, with 14 chromosomes spanned by a single contig and nine chromosomes spanned by two contigs. This map predicts 20 to 27 major rearrangements distinguishing turkey and chicken chromosomes, despite up to 40 million years of separate evolution between the two species. These data elucidate the chromosomal evolutionary pattern within the Phasianidae that led to the modern turkey and chicken karyotypes. The predominant rearrangement mode involves intra-chromosomal inversions, and there is a clear bias for these to result in centromere locations at or near telomeres in turkey chromosomes, in comparison to interstitial centromeres in the orthologous chicken chromosomes.

Conclusion

The BAC-based turkey-chicken comparative map provides novel insights into the evolution of avian genomes, a framework for assembly of turkey whole genome shotgun sequencing data, and tools for enhanced genetic improvement of these important agricultural and model species.

Molecular and cellular evidence for the alternative lengthening of telomeres (ALT) mechanism in chicken

Telomere maintenance is an important genetic mechanism controlling cellular proliferation. Normally, telomeres are maintained by telomerase which is downregulated upon cellular differentiation in most somatic cell lineages. Telomerase activity is upregulated in immortalized cells and cancers to support an infinite lifespan and uncontrolled cell growth; however, some immortalized and transformed cells lack telomerase activity. Telomerase-negative tumors and immortalized cells utilize an alternative mechanism for maintaining telomeres termed alternative lengthening of telomeres (ALT). This research explored evidence for the ALT pathway in chicken cell lines by studying nontransformed immortalized cell lines (DF-1 and OU2) and comparing them to a normal (mortal) cell line and a transformed cell line (DT40). The research consisted of molecular and cellular analyses including profiling of telomeric DNA (array sizing and total content), telomerase activity, and expression of genes involved in the telomerase, recombination, and ALT pathways. In addition, an immunofluorescence analysis for an ALT marker, i.e. ALT-associated promyelocytic leukemia bodies (APBs), was conducted. Evidence for ALT was observed in the telomerase-negative immortalized cell lines. Additionally, the APB marker was also found in the other cell systems. The attributes of the chicken provide an additional vertebrate model for investigation of the ALT pathway.

Defining the turkey MHC: identification of expressed class I- and class IIB-like genes independent of the MHC-B

The MHC of the turkey (Meleagris gallopavo) is divided into two genetically unlinked regions; the MHC-B and MHC-Y. Although previous studies found the turkey MHC-B to be highly similar to that of the chicken, little is known of the gene content and extent of the MHC-Y. This study describes two partially overlapping large-insert BAC clones that genetically and physically map to the turkey MHC chromosome (MGA18) but to a region that assorts independently of MHC-B. Within the sequence assembly, 14 genes were predicted including new class I- and class IIB-like loci. Additional unassembled sequences corresponded to multiple copies of the ribosomal RNA repeat unit (18S-5.8S-28S). Thus, this newly identified MHC region appears to represent a physical boundary of the turkey MHC-Y. High-resolution multi-color fluorescence in situ hybridization studies confirm rearrangement of MGA18 relative to the orthologous chicken chromosome (GGA16) in regard to chromosome architecture, but not gene order. The difference in centromere position between the species is indicative of multiple chromosome rearrangements or alternate events such as neocentromere formation/centromere inactivation in the evolution of the MHC chromosome. Comparative sequencing of commercial turkeys (six amplicons totaling 7.6 kb) identified 68 single nucleotide variants defining nine MHC-Y haplotypes. Sequences of the new class I- and class IIB-like genes are most similar to MHC-Y genes in the chicken. All three loci are expressed in the spleen. Differential transcription of the MHC-Y class IIB-like loci was evident as one class IIB-like locus was only expressed in some individuals.

Chromosomal integration of an avian oncogenic herpesvirus reveals telomeric preferences and evidence for lymphoma clonality

BACKGROUND

Herpesviruses are a major health concern for numerous organisms, including humans, causing both acute and chronic infections recurrent over an individual's lifespan. Marek's disease virus (MDV) is a highly contagious herpesvirus which causes a neoplastic condition in chicken populations. Several vertebrate-infecting herpesviruses have been shown to exist in an integrated state during latent periods of infection. However the status of MDV during latency has been a topic of debate.

RESULTS

Here we employed high-resolution multi-color fluorescence in situ hybridization (FISH) to show integration of MDV at the telomeres of chicken chromosomes. Cytogenomic mapping of the chromosomal integrations allowed us to examine the clonal relationships among lymphomas within individuals, whereas analysis of tumors from multiple individuals indicated the potential for chromosomal preferences.

CONCLUSIONS

Our data highlight that substantive genome-level interactions between the virus and host exist, and merit consideration for their potential impact and role in key aspects of herpesvirus pathobiology including infection, latency, cellular transformation, latency-breaks and viral evolution.

Telomere biology in Metazoa

In this review we present critical overview of some of the available literature on the fundamental biology of telomeres and telomerase in Metazoan. With the exception of Nematodes and Arthropods, the (TTAGGG)(n) sequence is conserved in most Metazoa. Available data show that telomerase-based end maintenance is a very ancient mechanism in unicellular and multicellular organisms. In invertebrates, fish, amphibian, and reptiles persistent telomerase activity in somatic tissues might allow the maintenance of the extensive regenerative potentials of these species. Telomerase repression among birds and many mammals suggests that, as humans, they may use replicative aging as a tumor protection mechanism.