Chromosomal mapping of chicken mega-telomere arrays to GGA9, 16, 28 and W using a cytogenomic approach

The evolution of early-life telomere length, pace-of-life and telomere-chromosome length dynamics in birds

Telomeres, the short DNA sequences that protect chromosome ends, are an ancient molecular structure, which is highly conserved across most eukaryotes. Species differ in their telomere lengths, but the causes of this variation are not well understood. Here, we demonstrate that mean early-life telomere length is an evolutionary labile trait across 57 bird species (representing 35 families in 12 orders) with the greatest trait diversity found among passerines. Among these species, telomeres are significantly shorter in fast-lived than in slow-lived species, suggesting that telomere length may have evolved to mediate trade-offs between physiological requirements underlying the diversity of pace-of-life strategies in birds. This association was attenuated when excluding studies that may include interstitial telomeres in the estimation of mean telomere length. Curiously, within some species, larger individual chromosome size predicts longer telomere lengths on that chromosome, leading to the hypothesis that telomere length also covaries with chromosome length across species. We show that longer mean chromosome length or genome size tends to be associated with longer mean early-life telomere length (measured across all chromosomes) within a phylogenetic framework constituting up to 31 bird species. These associations were strengthened when excluding highly influential outliers. However, sensitivity analyses suggested that they were susceptible to sample size effects and not robust to the exclusion of studies that may include interstitial telomeres. Combined, our analyses generalize patterns previously found within a few species and provide potential adaptive explanations for the 10-fold variation in telomere lengths observed among birds.

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.

Telomeres, species differences, and unusual telomeres in vertebrates: presenting challenges and opportunities to understanding telomere dynamics

There has been increasing interest in the use of telomeres as biomarkers of stress, cellular ageing and life-histories. However, the telomere landscape is a diverse feature, with noticeable differences between species, a fact which is highlighted by the unusual telomeres of various vertebrate organisms. We broadly review differences in telomere dynamics among vertebrates, and emphasize the need to understand more about telomere processes and trends across species. As part of these species differences, we review unusual telomeres in vertebrates. This includes mega-telomeres, which are present across a diverse set of organisms, but also focusing on the unusual telomeres traits of marsupials and monotremes, which have seen little to no prior discussion, yet uniquely stand out from other unusual telomere features discovered thus far. Due to the presence of at least two unique telomere features in the marsupial family Dasyuridae, as well as to the presence of physiological strategies semelparity and torpor, which have implications for telomere life-histories in these species, we suggest that this family has a very large potential to uncover novel information on telomere evolution and dynamics.

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.

Evolutionary rates of mammalian telomere-stability genes correlate with karyotype features and female germline expression

Telomeres protect the ends of eukaryotic chromosomes and are essential for cell viability. In mammals, telomere dynamics vary with life history traits (e.g. body mass and longevity), suggesting differential selection depending on physiological characteristics. Telomeres, in analogy to centromeric regions, also represent candidate meiotic drivers and subtelomeric DNA evolves rapidly. We analyzed the evolutionary history of mammalian genes implicated in telomere homeostasis (TEL genes). We detected widespread positive selection and we tested two alternative hypotheses: (i) fast evolution is driven by changes in life history traits; (ii) a conflict with selfish DNA elements at the female meiosis represents the underlying selective pressure. By accounting for the phylogenetic relationships among mammalian species, we show that life history traits do not contribute to shape diversity of TEL genes. Conversely, the evolutionary rate of TEL genes correlates with expression levels during meiosis and episodes of positive selection across mammalian species are associated with karyotype features (number of chromosome arms). We thus propose a telomere drive hypothesis, whereby (sub)telomeres and telomere-binding proteins are engaged in an intra-genomic conflict similar to the one described for centromeres.

Sex chromosome repeats tip the balance towards speciation

Because sex chromosomes, by definition, carry genes that determine sex, mutations that alter their structural and functional stability can have immediate consequences for the individual by reducing fertility, but also for a species by altering the sex ratio. Moreover, the sex-specific segregation patterns of heteromorphic sex chromosomes make them havens for selfish genetic elements that not only create suboptimal sex ratios but can also foster sexual antagonism. Compensatory mutations to mitigate antagonism or return sex ratios to a Fisherian optimum can create hybrid incompatibility and establish reproductive barriers leading to species divergence. The destabilizing influence of these selfish elements is often manifest within populations as copy number variants (CNVs) in satellite repeats and transposable elements (TE) or as CNVs involving sex-determining genes, or genes essential to fertility and sex chromosome dosage compensation. This review catalogs several examples of well-studied sex chromosome CNVs in Drosophilids and mammals that underlie instances of meiotic drive, hybrid incompatibility and disruptions to sex differentiation and sex chromosome dosage compensation. While it is difficult to pinpoint a direct cause/effect relationship between these sex chromosome CNVs and speciation, it is easy to see how their effects in creating imbalances between the sexes, and the compensatory mutations to restore balance, can lead to lineage splitting and species formation.

Marek's disease herpesvirus vaccines integrate into chicken host chromosomes yet lack a virus-host phenotype associated with oncogenic transformation

Marek's disease (MD) is a lymphotropic and oncogenic disease of chickens that can lead to death in susceptible and unvaccinated host birds. The causative pathogen, MD virus (MDV), a highly oncogenic alphaherpesvirus, integrates into host genome near the telomeres. MD occurrence is controlled across the globe by biosecurity, selective breeding for enhanced MD genetic resistance, and widespread vaccination of flocks using attenuated serotype 1 MDV or other serotypes. Despite over 40 years of usage, the specific mechanism(s) of MD vaccine-related immunity and anti-tumor effects are not known. Here we investigated the cytogenetic interactions of commonly used MD vaccine strains of all three serotypes (HVT, SB-1, and Rispens) with the host to determine if all were equally capable of host genome integration. We also studied the dynamic profiles of chromosomal association and integration of the three vaccine strains, a first for MD vaccine research. Our cytogenetic data provide evidence that all three MD vaccine strains tested integrate in the chicken host genome as early as 1 day after vaccination similar to oncogenic strains. However, a specific, transformation-associated virus-host phenotype observed for oncogenic viruses is not established. Our results collectively provide an updated model of MD vaccine-host genome interaction and an improved understanding of the possible mechanisms of vaccinal immunity. Physical integration of the oncogenic MDV genome into host chromosomes along with cessation of viral replication appears to have joint signification in MDV's ability to induce oncogenic transformation. Whereas for MD vaccine serotypes, a sustained viral replication stage and lack of the chromosome-integrated only stage were shared traits during early infection.

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.

Chicken rRNA Gene Cluster Structure

Ribosomal RNA (rRNA) genes, whose activity results in nucleolus formation, constitute an extremely important part of genome. Despite the extensive exploration into avian genomes, no complete description of avian rRNA gene primary structure has been offered so far. We publish a complete chicken rRNA gene cluster sequence here, including 5'ETS (1836 bp), 18S rRNA gene (1823 bp), ITS1 (2530 bp), 5.8S rRNA gene (157 bp), ITS2 (733 bp), 28S rRNA gene (4441 bp) and 3'ETS (343 bp). The rRNA gene cluster sequence of 11863 bp was assembled from raw reads and deposited to GenBank under KT445934 accession number. The assembly was validated through in situ fluorescent hybridization analysis on chicken metaphase chromosomes using computed and synthesized specific probes, as well as through the reference assembly against de novo assembled rRNA gene cluster sequence using sequenced fragments of BAC-clone containing chicken NOR (nucleolus organizer region). The results have confirmed the chicken rRNA gene cluster validity.

Comparative Cytogenetics of the Congo African Grey Parrot (Psittacus erithacus)

The Congo African grey parrot (Psittacus erithacus, PER) is an endemic species of Central Africa, valued for its intelligence and listed as vulnerable due to poaching and habitat destruction. Improved knowledge about the P. erithacus genome is needed to address key biological questions and conservation of this species. The P. erithacus genome was studied using conventional and molecular cytogenetic approaches including Zoo-FISH. P. erithacus has a 'typical' parrot karyotype with 2n = 62-64 and 8 pairs of macrochromosomes. A distinct feature was a sharp macro-microchromosome boundary. Telomeric sequences were present at all chromosome ends and interstitially in PER2q, the latter coinciding with a C-band. NORs mapped to 4 pairs of microchromosomes which is in contrast to a single NOR in ancestral type avian karyotypes. Zoo-FISH with chicken macrochromosomes GGA1-9 and Z revealed patterns of conserved synteny similar to many other avian groups, though neighboring synteny combinations of GGA6/7, 8/9, and 1/4 were distinctive only to parrots. Overall, P. erithacus shared more Zoo-FISH patterns with neotropical macaws than Australian species such as cockatiel and budgerigar. The observations suggest that Psittaciformes karyotypes have undergone more extensive evolutionary rearrangements compared to the majority of other avian genomes.

Virus and host genomic, molecular, and cellular interactions during Marek's disease pathogenesis and oncogenesis

Marek's Disease Virus (MDV) is a chicken alphaherpesvirus that causes paralysis, chronic wasting, blindness, and fatal lymphoma development in infected, susceptible host birds. This disease and its protective vaccines are highly relevant research targets, given their enormous impact within the poultry industry. Further, Marek's disease (MD) serves as a valuable model for the investigation of oncogenic viruses and herpesvirus patterns of viral latency and persistence--as pertinent to human health as to poultry health. The objectives of this article are to review MDV interactions with its host from a variety of genomic, molecular, and cellular perspectives. In particular, we focus on cytogenetic studies, which precisely assess the physical status of the MDV genome in the context of the chicken host genome. Combined, the cytogenetic and genomic research indicates that MDV-host genome interactions, specifically integration of the virus into the host telomeres, is a key feature of the virus life cycle, contributing to the viral achievement of latency, transformation, and reactivation of lytic replication. We present a model that outlines the variety of virus-host interactions, at the multiple levels, and with regard to the disease states.

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.

Repetitive Sequence and Sex Chromosome Evolution in Vertebrates

Sex chromosomes are the most dynamic entity in any genome having unique morphology, gene content, and evolution. They have evolved multiple times and independently throughout vertebrate evolution. One of the major genomic changes that pertain to sex chromosomes involves the amplification of common repeats. It is hypothesized that such amplification of repeats facilitates the suppression of recombination, leading to the evolution of heteromorphic sex chromosomes through genetic degradation of Y or W chromosomes. Although contrasting evidence is available, it is clear that amplification of simple repetitive sequences played a major role in the evolution of Y and W chromosomes in vertebrates. In this review, we present a brief overview of the repetitive DNA classes that accumulated during sex chromosome evolution, mainly focusing on vertebrates, and discuss their possible role and potential function in this process.

Genome-wide characterization of insertion and deletion variation in chicken using next generation sequencing

Insertion and deletion (INDEL) is one of the main events contributing to genetic and phenotypic diversity, which receives less attention than SNP and large structural variation. To gain a better knowledge of INDEL variation in chicken genome, we applied next generation sequencing on 12 diverse chicken breeds at an average effective depth of 8.6. Over 1.3 million non-redundant short INDELs (1-49 bp) were obtained, the vast majority (92.48%) of which were novel. Follow-up validation assays confirmed that most (88.00%) of the randomly selected INDELs represent true variations. The majority (95.76%) of INDELs were less than 10 bp. Both the detected number and affected bases were larger for deletions than insertions. In total, INDELs covered 3.8 Mbp, corresponding to 0.36% of the chicken genome. The average genomic INDEL density was estimated as 0.49 per kb. INDELs were ubiquitous and distributed in a non-uniform fashion across chromosomes, with lower INDEL density in micro-chromosomes than in others, and some functional regions like exons and UTRs were prone to less INDELs than introns and intergenic regions. Nearly 620,253 INDELs fell in genic regions, 1,765 (0.28%) of which located in exons, spanning 1,358 (7.56%) unique Ensembl genes. Many of them are associated with economically important traits and some are the homologues of human disease-related genes. We demonstrate that sequencing multiple individuals at a medium depth offers a promising way for reliable identification of INDELs. The coding INDELs are valuable candidates for further elucidation of the association between genotypes and phenotypes. The chicken INDELs revealed by our study can be useful for future studies, including development of INDEL markers, construction of high density linkage map, INDEL arrays design, and hopefully, molecular breeding programs in chicken.

The genomic landscape of species divergence in Ficedula flycatchers

Unravelling the genomic landscape of divergence between lineages is key to understanding speciation. The naturally hybridizing collared flycatcher and pied flycatcher are important avian speciation models that show pre- as well as postzygotic isolation. We sequenced and assembled the 1.1-Gb flycatcher genome, physically mapped the assembly to chromosomes using a low-density linkage map and re-sequenced population samples of each species. Here we show that the genomic landscape of species differentiation is highly heterogeneous with approximately 50 'divergence islands' showing up to 50-fold higher sequence divergence than the genomic background. These non-randomly distributed islands, with between one and three regions of elevated divergence per chromosome irrespective of chromosome size, are characterized by reduced levels of nucleotide diversity, skewed allele-frequency spectra, elevated levels of linkage disequilibrium and reduced proportions of shared polymorphisms in both species, indicative of parallel episodes of selection. Proximity of divergence peaks to genomic regions resistant to sequence assembly, potentially including centromeres and telomeres, indicate that complex repeat structures may drive species divergence. A much higher background level of species divergence of the Z chromosome, and a lower proportion of shared polymorphisms, indicate that sex chromosomes and autosomes are at different stages of speciation. This study provides a roadmap to the emerging field of speciation genomics.

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

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 mapping and candidate gene discovery of chicken developmental mutants and genome-wide variation analysis of MHC congenics

The chicken has been widely used in experimental research given its importance to agriculture and its utility as a model for vertebrate biology and biomedical pursuits for over 100 years. Herein we used advanced technologies to investigate the genomic characteristics of specialized chicken congenic genetic resources developed on a highly inbred background. An Illumina 3K chicken single nucleotide polymorphism (SNP) array was utilized to study variation within and among major histocompatibility complex (MHC)-congenic lines as well as investigate line-specific genomic diversity, inbreeding coefficients, and MHC B haplotype-specific GGA 16 SNP profiles. We also investigated developmental mutant-congenic lines to map a number of single-gene mutations using both the Illumina 3K array and a recently developed Illumina 60K chicken SNP array. In addition to identifying the chromosomes and specific subregions, the mapping results affirmed prior analyses indicating recessive or dominant and autosomal or sex chromosome modes of inheritance. Priority candidate genes are described for each mutation based on association with similar phenotypes in other vertebrates. These single-gene mutations provide a means of studying amniote development and in particular serve as invaluable biomedical models for similar malformations found in human.

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.

Integrative mapping analysis of chicken microchromosome 16 organization

BACKGROUND

The chicken karyotype is composed of 39 chromosome pairs, of which 9 still remain totally absent from the current genome sequence assembly, despite international efforts towards complete coverage. Some others are only very partially sequenced, amongst which microchromosome 16 (GGA16), particularly under-represented, with only 433 kb assembled for a full estimated size of 9 to 11 Mb. Besides the obvious need of full genome coverage with genetic markers for QTL (Quantitative Trait Loci) mapping and major genes identification studies, there is a major interest in the detailed study of this chromosome because it carries the two genetically independent MHC complexes B and Y. In addition, GGA16 carries the ribosomal RNA (rRNA) genes cluster, also known as the NOR (nucleolus organizer region). The purpose of the present study is to construct and present high resolution integrated maps of GGA16 to refine its organization and improve its coverage with genetic markers.

RESULTS

We developed 79 STS (Sequence Tagged Site) markers to build a physical RH (radiation hybrid) map and 34 genetic markers to extend the genetic map of GGA16. We screened a BAC (Bacterial Artificial Chromosome) library with markers for the MHC-B, MHC-Y and rRNA complexes. Selected clones were used to perform high resolution FISH (Fluorescent In Situ Hybridization) mapping on giant meiotic lampbrush chromosomes, allowing meiotic mapping in addition to the confirmation of the order of the three clusters along the chromosome. A region with high recombination rates and containing PO41 repeated elements separates the two MHC complexes.

CONCLUSIONS

The three complementary mapping strategies used refine greatly our knowledge of chicken microchromosome 16 organisation. The characterisation of the recombination hotspots separating the two MHC complexes demonstrates the presence of PO41 repetitive sequences both in tandem and inverted orientation. However, this region still needs to be studied in more detail.

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

This study investigated telomeric array organization of diverse chicken genotypes utilizing in vivo and in vitro cells having phenotypes with different proliferation potencies. Our experimental objective was to characterize the extent and nature of array variation present to explore the hypothesis that mega-telomeres are a universal and fixed feature of chicken genotypes. Four different genotypes were studied including normal (UCD 001, USDA-ADOL Line 0), immortalized (DF-1), and transformed (DT40) cells. Both cytogenetic and molecular approaches were utilized to develop an integrated view of telomeric array organization. It was determined that significant variation exists within and among chicken genotypes for chromosome-specific telomeric array organization and total genomic-telomeric sequence content. Although there was variation for mega-telomere number and distribution, two mega-telomere loci were in common among chicken genetic lines (GGA 9 and GGA W). The DF-1 cell line was discovered to maintain a complex derivative karyotype involving chromosome fusions in the homozygous and heterozygous condition. Also, the DF-1 cell line was found to contain the greatest amount of telomeric sequence per genome (17%) as compared to UCD 001 (5%) and DT40 (1.2%). The chicken is an excellent model for studying unique and universal features of vertebrate telomere biology, and characterization of the telomere length variation among genotypes will be useful in the exploration of mechanisms controlling telomere length maintenance in different cell types having unique phenotypes.

Architecture and organization of chicken microchromosome 16: order of the NOR, MHC-Y, and MHC-B subregions

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

Evolution of "determinants" in sex-determination: a novel hypothesis for the origin of environmental contingencies in avian sex-bias

Sex-determination is commonly categorized as either "genetic" or "environmental"-a classification that obscures the origin of this dichotomy and the evolution of sex-determining factors. The current focus on static outcomes of sex-determination provides little insight into the dynamic developmental processes by which some mechanisms acquire the role of sex determinants. Systems that combine "genetic" pathways of sex-determination (i.e., sex chromosomes) with "environmental" pathways (e.g., epigenetically induced segregation distortion) provide an opportunity to examine the evolutionary relationships between the two classes of processes and, ultimately, illuminate the evolution of sex-determining systems. Taxa with sex chromosomes typically undergo an evolutionary reduction in size of one of the sex chromosomes due to suppressed recombination, resulting in pronounced dimorphism of the sex chromosomes, and setting the stage for emergence of epigenetic compensatory mechanisms regulating meiotic segregation of heteromorphic sex chromosomes. Here we propose that these dispersed and redundant regulatory mechanisms enable environmental contingency in genetic sex-determination in birds and account for frequently documented context-dependence in avian sex-determination. We examine the evolution of directionality in such sex-determination as a result of exposure of epigenetic regulators of meiosis to natural selection and identify a central role of hormones in integrating female reproductive homeostasis, resource allocation to oocytes, and offspring sex. This approach clarifies the evolutionary relationship between sex-specific molecular genetic mechanisms of sex-determination and non-sex-specific epigenetic regulators of meiosis and demonstrates that both can determine sex. Our perspective shows how non-sex-specific mechanisms can acquire sex-determining function and, by establishing the explicit link between physiological integration of oogenesis and sex-determination, opens new avenues to the studies of adaptive sex-bias and sex-specific resource allocation in species with genetic sex-determination.

Addition of the microchromosome GGA25 to the chicken genome sequence assembly through radiation hybrid and genetic mapping

Background

The publication of the first draft chicken sequence assembly became available in 2004 and was updated in 2006. However, this does not constitute a definitive and complete sequence of the chicken genome, since the microchromosomes are notably under-represented. In an effort to develop maps for the microchromosomes absent from the chicken genome assembly, we developed radiation hybrid (RH) and genetic maps with markers isolated from sequence currently assigned to "chromosome Unknown" (chrUn). The chrUn is composed of sequence contigs not assigned to named chromosomes. To identify and map sequence belonging to the microchromosomes we used a comparative mapping strategy, and we focused on the small linkage group E26C13.

Results

In total, 139 markers were analysed with the chickRH6 panel, of which 120 were effectively assigned to the E26C13 linkage group, the remainder mapping elsewhere in the genome. The final RH map is composed of 22 framework markers extending over a 245.6 cR distance. A corresponding genetic map was developed, whose length is 103 cM in the East Lansing reference population. The E26C13 group was assigned to GGA25 (Gallus gallus chromosome 25) by FISH (fluorescence in situ hybridisation) mapping.

Conclusion

The high-resolution RH framework map obtained here covers the entire chicken chromosome 25 and reveals the existence of a high number of intrachromosomal rearrangements when compared to the human genome. The strategy used here for the characterization of GGA25 could be used to improve knowledge on the other uncharacterized small, yet gene-rich microchromosomes.

Characterization of the turkey MHC chromosome through genetic and physical mapping

Previous studies in the chicken have identified a single microchromosome (GGA16) containing the ribosomal DNA (rDNA) and two genetically unlinked MHC regions, MHC-B and MHC-Y. Chicken DNA sequence from these loci was used to develop PCR primers for amplification of homologous fragments from the turkey (Meleagris gallopavo). PCR products were sequenced and overgo probes were designed to screen the CHORI 260 turkey BAC library. BAC clones corresponding to the turkey rDNA, MHC-B and MHC-Y were identified. BAC end and subclone sequencing confirmed identity and homology of the turkey BAC clones to the respective chicken loci. Based on subclone sequences, single-nucleotide polymorphisms (SNPs) segregating within the UMN/NTBF mapping population were identified and genotyped. Analysis of SNP genotypes found the B and Y to be genetically unlinked in the turkey. Silver staining of metaphase chromosomes identified a single pair of microchromosomes with nucleolar organizer regions (NORs). Physical locations of the rDNA and MHC loci were determined by fluorescence in situ hybridization (FISH) of the BAC clones to metaphase chromosomes. FISH clearly positioned the rDNA distal to the Y locus on the q-arm of the MHC chromosome and the MHC-B on the p-arm. An internal telomere array on the MHC chromosome separates the B and Y loci.