Search

A comprehensive characterization of regulatory elements in the chicken genome across tissues will have substantial impacts on both fundamental and applied research. Here, we systematically identified and characterized regulatory elements in the chicken genome by integrating 377 genome-wide sequencing datasets from 23 adult tissues. In total, we annotated 1.57 million regulatory elements, representing 15 distinct chromatin states, and predicted about 1.2 million enhancer-gene pairs and 7662 super-enhancers. This functional annotation of the chicken genome should have wide utility on identifying regulatory elements accounting for gene regulation underlying domestication, selection, and complex trait regulation, which we explored. In short, this comprehensive atlas of regulatory elements provides the scientific community with a valuable resource for chicken genetics and genomics.

MHCY is a second major histocompatibility complex-like gene region in chickens originally identified by the presence of major histocompatibility complex class I-like and class II-like gene sequences. Up to now, the MHCY gene region has been poorly represented in genomic sequence data. A high density of repetitive sequence and multiple members of several gene families prevented the accurate assembly of short-read sequence data for MHCY. Identified here by single-molecule real-time sequencing sequencing of BAC clones for the Gallus gallus Red Jungle Fowl reference genome are 107 MHCY region genes (45 major histocompatibility complex class I-like, 41 c-type-lectin-like, 8 major histocompatibility complex class IIβ, 8 LENG9-like, 4 zinc finger protein loci, and a single only zinc finger-like locus) located amid hundreds of retroelements within 4 contigs representing the region. Sequences obtained for nearby ribosomal RNA genes have allowed MHCY to be precisely mapped with respect to the nucleolar organizer region. Gene sequences provide insights into the unusual structure of the MHCY class I molecules. The MHCY class I loci are polymorphic and group into 22 types based on predicted amino acid sequences. Some MHCY class I loci are full-length major histocompatibility complex class I genes. Others with altered gene structure are considered gene candidates. The amino acid side chains at many of the polymorphic positions in MHCY class I are directed away rather than into the antigen-binding groove as is typical of peptide-binding major histocompatibility complex class I molecules. Identical and nearly identical blocks of genomic sequence contribute to the observed multiplicity of identical MHCY genes and the large size (>639 kb) of the Red Jungle Fowl MHCY haplotype. Multiple points of hybridization observed in fluorescence in situ hybridization suggest that the Red Jungle Fowl MHCY haplotype is made up of linked, but physically separated genomic segments. The unusual gene content, the evidence of highly similar duplicated segments, and additional evidence of variation in haplotype size distinguish polymorphic MHCY from classical polymorphic major histocompatibility complex regions.

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 (&gt, 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.

Gene regulatory elements are central drivers of phenotypic variation and thus of critical importance towards understanding the genetics of complex traits. The Functional Annotation of Animal Genomes consortium was formed to collaboratively annotate the functional elements in animal genomes, starting with domesticated animals. Here we present an expansive collection of datasets from eight diverse tissues in three important agricultural species: chicken (Gallus gallus), pig (Sus scrofa), and cattle (Bos taurus). Comparative analysis of these datasets and those from the human and mouse Encyclopedia of DNA Elements projects reveal that a core set of regulatory elements are functionally conserved independent of divergence between species, and that tissue-specific transcription factor occupancy at regulatory elements and their predicted target genes are also conserved. These datasets represent a unique opportunity for the emerging field of comparative epigenomics, as well as the agricultural research community, including species that are globally important food resources., In order to interpret non-coding variants, information about regulatory elements in the genome is essential. Here, the authors annotate regulatory elements in chicken, pig and cattle, and characterize conservation of these elements between species.

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.

Summary An embryonic lethal mutation in chicken named cleft primary palate (cpp) is inherited in an autosomal recessive mode and results in a severely truncated upper beak. In this study, genotyping and sequencing techniques were employed to advance our genetic and genomic knowledge of the mutation’s chromosomal location, candidate region and possible causative element using a congenic inbred line. Herein, the candidate region for the cpp developmental mutation was established as a ca. 5.1 Mb region of chicken chromosome 11 (GGA 11) through the use of a 600K Affymetrix SNP array. The SNPs identified from this array linked to cpp were used to genotype individuals from the congenic inbred line over several generations and thereby fine‐map the causative region resulting in an approximately 200 kb size reduction. This candidate region (4.9 Mb) was sequenced via capture array in a cohort of 24 individuals, including carriers, mutants and their wild type (wt) siblings. Interestingly, the GGA 11 region for cpp encompasses the predicted centromere location and is thus unlikely to be highly disrupted by further recombination. Here we report on the variation unique to the cpp mutation, i.e. single‐nucleotide variants and insertions or deletions. Although the candidate region contains several genes of interest with regard to the cpp phenotype, only one cpp‐linked variant was predicted to have a significant physiological effect by causing a frameshift mutation in ESRP2, which has a role in tissue‐specific splicing during development.

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.

The chicken wingless-2 (wg-2) mutation is inherited in an autosomal recessive fashion, and the resulting phenotype in mutant (wg-2/wg-2) individuals is a developmental syndrome characterized by absent wings, truncated legs, craniofacial as well as skin and feather defects, and kidney malformations. Mapping and genotyping established that the mutation resides within 227 kilobases (kb) of chromosome 12 in a wg-2 congenic inbred line. A capture array was designed to target and sequence the candidate region along with flanking DNA in 24 birds from the line. Many point mutations and insertions or deletions were identified, and analysis of the linked variants indicated a point mutation predicted to cause a premature stop codon in the RAF1 gene. Expression studies were conducted inclusive of all genes in the candidate region. Interestingly, RAF1 transcription was elevated, yet the protein was absent in the mutants relative to normal individuals. RAF1 encodes a protein integral to the Ras/Raf/MAPK signaling pathway controlling cellular proliferation, and notably, human RASopathies are developmental syndromes caused by germline mutations in genes of this pathway. Our work indicates RAF1 as the priority candidate causative gene for wg-2 and provides a new animal model to study an important signaling pathway implicated in limb development, as well as RASopathies.

Chromosome complements in avian species include both constant and variable features that must be considered in studies to map genes, identify sites of integration of foreign DNA, and track donor cells to tissue sites in chimeras. Detailed cytogenetic studies in the domestic chicken are restricted primarily to the first 10 chromosome pairs, each of which is morphologically distinct. High-resolution banding permits analyses of yet smaller chromosomes (e.g., chromosomes 1 – 15). The primary constitutive features of the avian genome have been revealed by G banding for structural landmarks, C banding for constitutive heterochromatin, RBG banding for DNA replication patterns, and in situ hybridization for localizing highly repeated sequences to centromeric and telomeric regions. Variable features of the chromosome complement and of gene expression in selected repeated gene families provide convenient metaphase and interphase phenotypes for studies in vitro and in vivo. A significant amount of genomic variability has been observed among early chick embryos (1 – 15% aberration frequency among genetic lines), including haploidy, triploidy, trisomy, and mosaicism. The inadvertent selection of an aberrant embryo as a source of cells for chimeric production could yield surprising or confusing results if cells are not karyotyped. Variability at the ribosomal RNA gene cluster generates polymorphic nucleolar (PNU) patterns in interphase cells. Polymorphic cells have one macro- and one micronucleolus per cell (Pp). This phenotype is easily diagnosed in cytological preparations. We have developed chicken genetic strains with defined PNU patterns. A particular nucleolar variant is expressed in the embryo and at all other stages including adult. We have detected the PNU phenotype in stage X embryos, in embryonic tissues representing ectoderm, mesoderm, and endoderm, and in feather pulp cells from chickens. Heterozygous PNU cells should be easy to track to tissue sites when transferred to nonpolymorphic recipient embryos. While development of homozygous PNU embryos (pp) is arrested at early primitive streak formation, heterozygotes (Pp) develop normally. Studies with rDNA variants should be useful for investigating the relative contributions of maternal ribosomes versus embryo-derived rRNA in early developmental processes.

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.

The chicken has been a particularly useful model for the study of craniofacial development and disease for over a century due to their relatively large size, accessibility, and amenability for classical bead implantation and transplant experiments. Several naturally occurring mutant lines with craniofacial anomalies also exist and have been heavily utilized by developmental biologist for several decades. Two of the most well known lines, talpid(2) (ta(2)) and talpid(3) (ta(3)), represent the first spontaneous mutants to have the causative genes identified. Despite having distinct genetic causes, both mutants have recently been identified as ciliopathic. Excitingly, both of these mutants have been classified as models for human craniofacial ciliopathies: Oral-facial-digital syndrome (ta(2)) and Joubert syndrome (ta(3)). Herein, we review and compare these two models of craniofacial disease and highlight what they have revealed about the molecular and cellular etiology of ciliopathies. Furthermore, we outline how applying classical avian experiments and new technological advances (transgenics and genome editing) with naturally occurring avian mutants can add a tremendous amount to what we currently know about craniofacial ciliopathies.

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.

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

Next-generation sequencing (NGS) and expression technologies were utilized to investigate the genes and sequence elements in a 586 kb region of chicken chromosome 1 associated with the autosomal recessive diplopodia-1 (dp-1) mutation. This mutation shows a syndromic phenotype similar to known human developmental abnormalities (e.g., cleft palate, polydactyly, omphalocele [exposed viscera]). Toward our goal to ascertain the variant responsible, the entire 586 kb region was sequenced following utilization of a specifically designed capture array and to confirm/validate fine-mapping results. Bioinformatic analyses identified a total of 6142 sequence variants, which included SNPs, indels, and gaps. Of these, 778 SNPs, 146 micro-indels, and 581 gaps were unique to the UCD-Dp-1.003 inbred congenic line; those found within exons and splice sites were studied for contribution to the mutant phenotype. Upon further validation with additional mutant samples, a smaller subset (of variants [51]) remains linked to the mutation. Additionally, utilization of specific samples in the NGS technology was advantageous in that fine-mapping methodologies eliminated an additional 326 kb of sequence information on chromosome 1. Predicted and confirmed protein-coding genes within the smaller 260 kb region were assessed for their developmental expression patterns over several stages of early embryogenesis in regions/tissues of interest (e.g., digits, craniofacial region). Based on these results and known function in other vertebrates, 2 genes within 5 kb of each other, MRE11 and GPR83, are proposed as high-priority candidates for the dp-1 mutation.

It is now over 10 years since the first avian genome [International Chicken Genome Sequencing Consortium, 2004] and the first complete avian karyotype [Masabanda et al., 2004] were both published; however, until 2014, avian cytogenetics has focused heavily on descriptive studies [e.g. Griffin et al., 2007, 2008; Skinner et al., 2009; Volker et al., 2010] with less attention to its functional relevance. Last year, however, saw 2 landmark efforts in the chromosomal studies of birds: a special issue of Chromosome Research in April and the announcement of recently completed sequences of multiple new avian genomes in Science and the BMC journals (taking the total number sequenced to over 50) in December. Studying the chromosomes of birds is, perhaps for the first time, telling us more about avian biology, function and evolution than it ever has... Conclusions. The most recent advances in avian cytogenetics have culminated in great promise not only for the study of bird karyotypes, but also for providing insight into the mechanisms of chromosome evolution in general. New avenues for investigation include gene regulation; for instance, it will become necessary to map accurately the physical location of polyadenylation and transcription start sites, important reference points that define promoters and post-transcriptional regulation. It will also become possible to sequence full-length transcripts, to allow accurate identification of alternate splicing events and their controlling elements. The ENCODE (Encyclopedia of DNA Elements) project has helped to define functional elements of the human genome, including those aforementioned as well as other chromatin signals, e.g. active chromatin, enhancers, insulators, methylation domains, etc. An effort of agENCODE is underway to include agriculturally important birds such as chicken, turkey, duck, quail, and perhaps ostrich. The study of cytogenetics will be essential here in helping to define higher-order structures in nuclear organization that show regulatory interactions within and between chromosomes. Finally, reconstruction of evolutionary events allows us to study genome organization and function not only in extant but, by extrapolation, in extinct species also. Reconstruction of avian-reptilian ancestral karyotypes will allow us to define chromosomal rearrangements in long-dead species that have captured the public imagination. Here be dragons!

Wingless-2 (wg-2) is an autosomal recessive mutation in chicken that results in an embryonic lethal condition. Affected individuals exhibit a multisystem syndrome characterized by absent wings, truncated legs, and craniofacial, kidney, and feather malformations. Previously, work focused on phenotype description, establishing the autosomal recessive pattern of Mendelian inheritance and placing the mutation on an inbred genetic background to create the congenic line UCD Wingless-2.331. The research described in this paper employed the complementary tools of breeding, genetics, and genomics to map the chromosomal location of the mutation and successively narrow the size of the region for analysis of the causative element. Specifically, the wg-2 mutation was initially mapped to a 7 Mb region of chromosome 12 using an Illumina 3 K SNP array. Subsequent SNP genotyping and exon sequencing combined with analysis from improved genome assemblies narrowed the region of interest to a maximum size of 227 kb. Within this region, 3 validated and 3 predicted candidate genes are found, and these are described. The wg-2 mutation is a valuable resource to contribute to an improved understanding of the developmental pathways involved in chicken and avian limb development as well as serving as a model for human development, as the resulting syndrome shares features with human congenital disorders.

talpid(2) is an avian autosomal recessive mutant with a myriad of congenital malformations, including polydactyly and facial clefting. Although phenotypically similar to talpid(3), talpid(2) has a distinct facial phenotype and an unknown cellular, molecular and genetic basis. We set out to determine the etiology of the craniofacial phenotype of this mutant. We confirmed that primary cilia were disrupted in talpid(2) mutants. Molecularly, we found disruptions in Hedgehog signaling. Post-translational processing of GLI2 and GLI3 was aberrant in the developing facial prominences. Although both GLI2 and GLI3 processing were disrupted in talpid(2) mutants, only GLI3 activator levels were significantly altered in the nucleus. Through additional fine mapping and whole-genome sequencing, we determined that the talpid(2) phenotype was linked to a 1.4 Mb region on GGA1q that contained the gene encoding the ciliary protein C2CD3. We cloned the avian ortholog of C2CD3 and found its expression was ubiquitous, but most robust in the developing limbs and facial prominences. Furthermore, we found that C2CD3 is localized proximal to the ciliary axoneme and is important for docking the mother centriole to the ciliary vesicle and cell membrane. Finally, we identified a 19 bp deletion in talpid(2) C2CD3 that produces a premature stop codon, and thus a truncated protein, as the likely causal allele for the phenotype. Together, these data provide insight into the cellular, molecular and genetic etiology of the talpid(2) phenotype. Our data suggest that, although the talpid(2) and talpid(3) mutations affect a common ciliogenesis pathway, they are caused by mutations in different ciliary proteins that result in differences in craniofacial phenotype.

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.

Marek's disease virus (MDV) is an oncogenic α-herpesvirus that induces Marek's disease characterized by fatal lymphomas in chickens. Here, we explored the timing during pathogenesis when the virus integrates into the host genome, the cell type involved, the role of viral integration on cellular transformation, and tumor clonality. Three immune organs of chicken (thymus, bursa, and spleen) were extracted following infection with either an oncogenic or a non-oncogenic strain of MDV. Using molecular cytogenetics, cells were investigated for viral integration at key time points throughout pathogenesis. Integration profiling of tumors (early to late stage) was conducted. Virus integration was widespread in B and T lymphocytes based on their abundance in bursa and thymus, respectively. Viral replication was detected early after infection as was viral integration into the host genome. Integration is a natural part of the MDV herpesvirus life cycle. In addition, our data using a non-oncogenic virus establish that although integration is a hallmark of tumor cell populations, integration alone is not sufficient for cellular transformation. Our results provide evidence for progression of lineage clonality within tumors. Understanding the features of integration provides insight into the mechanisms of herpesvirus pathology which could lead to disease mitigation strategies.

Trisomy mapping is a powerful method for assigning genes to chicken microchromosome 16 (GGA 16). The single chicken nucleolar organizer region (NOR), the 2 major histocompatibility complex regions (MHC-Y and MHC-B), and CD1 genes were all previously assigned to GGA 16 using trisomy mapping. Here, we combined array comparative genomic hybridization with trisomy mapping to screen unassigned genomic scaffolds (consigned temporarily to chrUn_random) for sequences originating from GGA 16. A number of scaffolds mapped to GGA 16. Among these were scaffolds that contain genes for olfactory (OR) and cysteine-rich domain scavenger (SRCR) receptors, along with a number of genes that encode putative immunoglobulin-like receptors and other molecules. We used high-resolution cytogenomic analyses to confirm assignment of OR and SRCR genes to GGA 16 and to pinpoint members of these gene families to the q-arm in partially overlapping regions between the centromere and the NOR. Southern blots revealed sequence polymorphism within the OR/SRCR region and linkage with the MHC-Y region, thereby providing evidence for conserved linkage between OR genes and the MHC within birds. This work localizes OR genes to the vicinity of the chicken MHC and assigns additional genes, including immune defense genes, to GGA 16.

The first build of the chicken genome sequence appeared in March, 2004 – the first genome sequence of any animal agriculture species. That sequence was done primarily by whole genome shotgun Sanger sequencing, along with the use of an extensive BAC contig-based physical map to assemble the sequence contigs and scaffolds and align them to the known chicken chromosomes and linkage groups. Subsequent sequencing and mapping efforts have improved upon that first build, and efforts continue in search of missing and/or unassembled sequence, primarily on the smaller microchromosomes and the sex chromosomes. In the past year, a draft turkey genome sequence of similar quality has been obtained at a much lower cost primarily due to the development of ‘next-generation’ sequencing techniques. However, assembly and alignment of the sequence contigs and scaffolds still depended on a detailed BAC contig map of the turkey genome that also utilized comparison to the existing chicken sequence. These 2 land fowl (Galliformes) genomes show a remarkable level of similarity, despite an estimated 30–40 million years of separate evolution since their last common ancestor. Among the advantages offered by these sequences are routine re-sequencing of commercial and research lines to identify the genetic correlates of phenotypic change (for example, selective sweeps), a much improved understanding of poultry diversity and linkage disequilibrium, and access to high-density SNP typing and association analysis, detailed transcriptomic and proteomic studies, and the use of genome-wide marker- assisted selection to enhance genetic gain in commercial stocks.

Division-dependent telomere shortening correlating with age triggers senescence on a cellular level and telomere dysfunction can facilitate oncogenesis. Therefore, the study of telomere biology is critical to the understanding of aging and cancer. The domestic chicken, a classic model for the study of developmental biology, possesses a telomere genome with highly conserved aspects and distinctive features which make it uniquely suited for the study of telomere maintenance mechanisms, their function and dysfunction. The purpose of this review is to highlight the chicken as a model for aging research, specifically as a model for telomere and telomerase research, and to increase its utility as such by describing developments in the study of chicken telomeres and telomerase in the context of related research in human and mouse.

Vertebrate immune responses to ectoparasites influence pathogen transmission and host fitness costs. Few studies have characterized natural immune responses to ectoparasites and resultant fitness effects on the ectoparasite. These are critical gaps in understanding vertebrate–ectoparasite interaction, disease ecology and host-parasite co-adaptation. This study focused on an ectoparasite of birds—the northern fowl mite (NFM) (Ornithonyssus sylviarum). Based on prior evidence that chickens develop resistance to NFM, these experiments tested two hypotheses: (i) skin inflammation blocks mite access to blood, impairing development, reproduction and survival; and (ii) host immunogenetic variation influences the inflammatory response and subsequent effects on the ectoparasite. On infested hosts, histology of skin inflammation revealed increased epidermal cell number and size, immigration of leukocytes and deposition of serous exudates on the skin surface. Survival of adult mites and their offspring decreased as the area of skin inflammation increased during an infestation. Inflammation increased the distance to blood vessels beyond the length of mite mouthparts (100–160 μm) and prevented protonymphs and adults from reaching a blood source. Consequently, protonymphs could not complete development, evidenced by a significant inverse relationship between inflammation and protonymph feeding success, as well as an increasing protonymph/adult ratio. Adult females were unable to feed and reproduce, indicated by an inverse relationship between inflammation and egg production, and decreasing female/juvenile ratio. These combined impacts of host inflammation reversed NFM population growth. Intensity of inflammation was influenced by the genotype of the major histocompatibility complex (MHC), supporting previous research that linked these immunological loci with NFM resistance. Overall, these data provide a model for a mechanism of avian resistance to an ectoparasitic arthropod and the fitness costs to the parasite of that host defense.

Four mega-telomere loci were mapped to chicken chromosomes 9, 16, 28, and the W sex chromosome by dual-color fluorescence in situ hybridization using a telomeric sequence probe and BAC clones previously assigned to chicken chromosomes. The in-common features of the mega-telomere chromosomes are that microchromosomes are involved rather than macrochromosomes; in three cases (9, 16, 28) acrocentrics are involved with the mega-telomeres mapping to the p arms. Three of the four chromosomes (9, 16, W) encode tandem repeats which in two cases (9 and 16) involve the ribosomal DNA arrays (the 5S and 18S-5.8S-28S gene repeats, respectively). All involved chromosomes have a typical-sized telomere on the opposite terminus. Intra- and interindividual variation for mega-telomere distribution are discussed in terms of karyotype abnormalities and the potential for mitotic instability of some telomeres. The diversity and distribution of telomere array quantity in the chicken genome should be useful in contributing to research related to telomere length regulation – how and by what mechanism genomes and individual chromosomes establish and maintain distinct sets of telomere array sizes, as well as for future studies related to stability of the chicken genome affecting development, growth, cellular lifespan and disease. An additional impact of this study includes the listing of BAC clones (26 autosomal and six W BACs tested) that were cytogenetically verified; this set of BACs provide a useful tool for future cytogenetic analyses of the microchromosomes.

Oral-facial-digital syndrome (OFD) is a ciliopathy that is characterized by oral-facial abnormalities, including cleft lip and/or palate, broad nasal root, dental anomalies, micrognathia and glossal defects. In addition, these individuals have several other characteristic abnormalities that are typical of a ciliopathy, including polysyndactyly, polycystic kidneys and hypoplasia of the cerebellum. Recently, a subset of OFD cases in humans has been linked to mutations in the centriolar protein C2 Ca2+-dependent domain-containing 3 (C2CD3). Our previous work identified mutations in C2CD3 as the causal genetic lesion for the avian talpid2 mutant. Based on this common genetic etiology, we re-examined the talpid2 mutant biochemically and phenotypically for characteristics of OFD. We found that, as in OFD-affected individuals, protein-protein interactions between C2CD3 and oral-facial-digital syndrome 1 protein (OFD1) are reduced in talpid2 cells. Furthermore, we found that all common phenotypes were conserved between OFD-affected individuals and avian talpid2 mutants. In light of these findings, we utilized the talpid2 model to examine the cellular basis for the oral-facial phenotypes present in OFD. Specifically, we examined the development and differentiation of cranial neural crest cells (CNCCs) when C2CD3-dependent ciliogenesis was impaired. Our studies suggest that although disruptions of C2CD3-dependent ciliogenesis do not affect CNCC specification or proliferation, CNCC migration and differentiation are disrupted. Loss of C2CD3-dependent ciliogenesis affects the dispersion and directional persistence of migratory CNCCs. Furthermore, loss of C2CD3-dependent ciliogenesis results in dysmorphic and enlarged CNCC-derived facial cartilages. Thus, these findings suggest that aberrant CNCC migration and differentiation could contribute to the pathology of oral-facial defects in OFD., Summary: This study utilizes a naturally occurring avian mutant known as talpid2 to determine the cellular basis for the oral-facial phenotypes present in oral-facial-digital syndrome.

Primordial germ cells (PGCs) are the precursors of sperm and eggs. In most animals, segregation of the germ line from the somatic lineages is one of the earliest events in development; in avian embryos, PGCs are first identified in an extra-embryonic region, the germinal crescent, after approximately 18 h of incubation. After 50-55 h of development, PGCs migrate to the gonad and subsequently produce functional sperm and oocytes. So far, cultures of PGCs that remain restricted to the germ line have not been reported in any species. Here we show that chicken PGCs can be isolated, cultured and genetically modified while maintaining their commitment to the germ line. Furthermore, we show that chicken PGCs can be induced in vitro to differentiate into embryonic germ cells that contribute to somatic tissues. Retention of the commitment of PGCs to the germ line after extended periods in culture and after genetic modification combined with their capacity to acquire somatic competence in vitro provides a new model for developmental biology. The utility of the model is enhanced by the accessibility of the avian embryo, which facilitates access to the earliest stages of development and supplies a facile route for the reintroduction of PGCs into the embryonic vasculature. In addition, these attributes create new opportunities to manipulate the genome of chickens for agricultural and pharmaceutical applications.

Poultry genetic resources, which are valuable for research, span an impressive gamut from breeds to highly specialized inbred lines. The community of scientists utilizing specialized lines is broad, including researchers in medicine, basic biology, and agricultural science. The majority of specialized research lines used by such scientists are held at land grant universities. Over the prior 2 decades, hundreds of lines were eliminated. This pattern continues today with no evidence of abatement. Awareness and visibility of the causes and ongoing problems have been highlighted via a number of high-profile forums. Given the large community of scientists and the negative impact on future advances in biological, medical, and agricultural research as these genetic resources dwindle, the issue is of national interest and warrants federal funding to support a united network of avian and poultry stocks centers.

Male and female embryonic stem (ES) cell lines were derived from the area pellucidae of Stage X (EG&K) chicken embryos. These ES cell lines were grown in culture for extended periods of time and the majority of the cells retained a diploid karyotype. When reintroduced into Stage VI-X (EG&K) recipient embryos, the cES cells were able to contribute to all somatic tissues. By combining irradiation of the recipient embryo with exposure of the cES cells to the embryonic environment in diapause, a high frequency and extent of chimerism was obtained. High-grade chimeras, indistinguishable from the donor phenotype by feather pigmentation, were produced. A transgene encoding GFP was incorporated into the genome of cES cells under control of the ubiquitous promoter CX and GFP was widely expressed in somatic tissues. Although cES cells made extensive contributions to the somatic tissues, contribution to the germline was not observed.

Telomere-binding proteins, their interaction partners and transcription factors play a prominent role in telomere maintenance and telomerase activation. We examined mRNA expression levels of tankyrase 1 and 2, TRF1 and 2, c-myc, TERT and TR in Gallus domesticus, the domestic chicken, by quantitative real-time PCR, establishing expression profiles for three contrasting cell systems: the pluripotent gastrula, differentiated embryo fibroblasts and transformed DT40 cells. All seven genes were up-regulated in DT40 cells compared to telomerase-negative CEFs and a majority of the genes were also up-regulated in the gastrula relative to CEFs. Surprisingly, we found TERT and TR transcripts in CEFs, albeit at low levels. TRF1 was down-regulated in the six CEF cultures by the time of culture growth arrest. A marked increase in the TRF2:TRF1 ratio occurred at or near senescence in all of the CEF cultures studied, with the most elevated ratio found in a short-lived culture in which TRF1 mRNA levels decreased two-fold and TRF2 levels increased 21-fold. This culture also showed highly reduced, degraded telomeres by Southern blot analysis. These data suggest that genes involved in telomere maintenance and telomerase induction are expressed differentially in pluripotent, differentiated and transformed cell systems.

The chicken genome, like those of most avian species, contains numerous microchromosomes that cannot be distinguished by size alone. Unique properties attributed to the microchromosomes include high GC content and gene density, and an enhanced recombination rate. Previously, microchromosome GGA 17 was shown to align with the consensus genetic linkage group E41W17, and bacterial artificial chromosome (BAC) clones containing E41W17 markers were isolated and assigned on the physical BAC map as well as the recently assembled draft chicken genome sequence. For this study, these same BACS were utilized as probes for fluorescence in-situ hybridization (FISH) to develop the GGA 17 cytogenetic map. Here we detail the chromosome order of ten BAC DNAs, thereby deriving a cytogenetic map of GGA 17 that is simultaneously integrated with both the linkage map and genome sequence. The location of the FISH probes together with the morphological appearance of the chromosome suggested that GGA 17 is an acrocentric chromosome whose cytogenetic map orientation is reversed from that currently indicated by the linkage map and draft genome sequence. The reversed orientation and the centromere location of GGA 17 were confirmed experimentally by dual-colour FISH hybridization using terminal BACs and the centromere-specific CNM oligonucleotide as probes. An advantage of this cyto-genomic approach is the improved alignment of the sequence and linkage maps with cytogenetic features such as the centromere, telomeres, p and q arms, and staining patterns indicating GC versus AT content.

Telomerase activity is essential for maintaining the termini of linear chromosomes. Telomerase consists of both a RNA and a specialized reverse transcriptase. Our objective for this study was to determine the molecular and cytogenetic features of the chicken telomerase reverse transcriptase (chTERT) gene and protein. The TERT mRNA from gastrula stage embryos was found to be 4497 bp in length, translating into a protein of 1346 amino acids (aa). The chTERT protein shares 45% aa identity with human TERT (hTERT). A distinctive feature of chTERT, as compared to human and other vertebrate TERTs, is the larger size of the protein due mainly to a considerably longer N-terminal flexible linker region (144 aa longer than in human). Chicken TERT was mapped to chromosome 2q21 near an interstitial telomere site. Several transcription factor binding motifs in the 5Vflanking/promoter region of chTERT were similar to those found associated with hTERT (E-box, Ik1, MAZ, Sp1 sites), whereas several c-Myb sites were found associated with chTERT only and c-Ets-2 and WT1 were associated with hTERT only. Results presented here should promote structure–function studies of chTERT, as well as contribute to the comparative analysis of TERT regulation and function in vertebrates utilizing the telomere clock mechanism to different degrees. D 2004 Elsevier B.V. All rights reserved.

Telomeres are the complex nucleoprotein structures at the termini of linear chromosomes. Telomeric DNA consists of a highly conserved hexanucleotide arranged in tandem repeats. Telomerase, a ribonucleoprotein of the reverse transcriptase family, specifies the sequence of telomeric DNA and maintains telomere array length. Numerous studies in model organisms established the significance of telomere structure and function in regulating genome stability, cellular aging, and oncogenesis. Our overall research objectives are to understand the organization of the telomere arrays in chicken in the context of the unusual organization and specialized features of this higher vertebrate genome (which include a compact genome, numerous microchromosomes, and high recombination rate) and to elucidate the role telomeres play in genome stability impacting cell function and life span. Recent studies found that the chicken genome contains three overlapping size classes of telomere arrays that differ in location and age-related stability: Class I 0.5 to 10 kb, Class II 10 to 40 kb, and Class III 40 kb to 2 Mb. Some notable features of chicken telomere biology are that the chicken genome contains ten times more telomeric DNA than the human genome and the Class III telomere arrays are the largest described for any vertebrate species. In vivo, chicken telomeres (Class II) shorten in an age-related fashion and telomerase activity is high in early stage embryos and developing organs but down-regulates during late embryogenesis or postnatally in most somatic tissues. In vitro, chicken cells down-regulate telomerase activity unless transformed. Knowledge of chicken telomere biology contributes information relevant to present and future biotechnology applications of chickens in vivo and chicken cells in vitro.

Although vertebrate telomeres are highly conserved, telomere dynamics and telomerase profiles vary among species. The objective of the present study was to examine telomerase activity and telomere length profiles of transformed and non-transformed avian cells in vitro. Non-transformed chicken embryo fibroblasts (CEFs) showed little or no telomerase activity from the earliest passages through senescence. Unexpectedly, a single culture of particularly long-lived senescent CEFs showed telomerase activity after over 250 days in culture. Transformed avian lines (six chicken, two quail and one turkey) and tumor samples (two chicken) exhibited telomerase activity. Telomere length profiles of non-transformed CEF cultures derived from individual embryos of an inbred line (UCD 003) exhibited cycles of shortening and lengthening with a substantial net loss of telomeric DNA by senescence. The telomere length profiles of several transformed cell lines resembled telomere length profiles of senescent CEFs in that they exhibited little of the typical smear of terminal restriction fragments (TRFs) suggesting that these transformed cells may possess a reduced amount of telomeric DNA. These results show that avian telomerase activity profiles are consistent with the telomerase activity profiles of human primary and transformed cells. Further, monitoring of telomere lengths of primary cells provides evidence for a dynamic series of changes over the lifespan of any specific cell culture ultimately resulting in net telomeric DNA loss by senescence.

The 5S ribosomal (r) RNA genes encode a small (∼120-bp) highly-conserved component of the large ribosomal subunit. The objective of the present research was to study the molecular and cytogenetic organization of the chicken 5S rDNA. A predominant 2.2-kb gene (5Sα) consisting of a coding and intergenic spacer (IGS) region was identified in ten research and commercial populations. A variant gene repeat of 0.6 kb (5Sβ) was observed in some of the populations. Genetic linkage analysis and cytogenetic localization by fluorescence in-situ hybridization assigned the 5S rDNA to chromosome 9. The 5S rDNA array was determined to be 80.2 ± 7.0 kb upon electrophoretic sizing following EcoRV digestion. Sequence analysis of 5Sα IGS regions revealed considerable conservation between chicken subspecies (98.4% identity) as well as homology with vertebrate Pol III promoter and regulatory sequence motifs. Minor intraindividual sequence variation within 1000 bp of IGS was observed in four cloned Red Jungle Fowl (Gallus gallus gallus) 5Sα repeats (95.5% identity in this region). Sequence comparisons between IGS regions of 5Sα and 5Sβ genes indicated two short continuous (>20 bp) and many short non-continuous homologous regions as well as other conserved features such as promoter and termination motifs.

Telomerase RNA (TR) is essential for telomerase activity and the maintenance of telomere length in proliferating cell populations. The objective of the present research was to define the cytogenetic and molecular genomic organization of chicken TR (chTR). The chTR exists as a single copy gene (TERC, alias TR), mapping to chromosome 9 (GGA9). The loci on the q arm of GGA9 map to three chromosomes in human with five of the nine GGA9q loci mapping to HSA3q. Sequencing of the chTERC locus (3,763 bp) from the UCD 001 genome (Red Jungle Fowl) included: 604 bp 5′, 465 coding, and 2,694 bp 3′ (from –604 to +3159). Sequence analysis included homology searches conducted on several levels including comparisons among different chicken genotypes, Marek’s disease virus (MDV) sequences, plus human and murine. We provide evidence for distal 5′ and 3′ sequence homology between chTERC and the MDV genome among other known regions of homology (promoter and coding), elaborate on 5′ transcription factor binding motifs among the various genomes as well as show type and number of TERT-related motifs 3′ of chicken TR (e.g., Sp1, c-Myb, c-Myc, AP2, among others). Surrounding the gene are more than 25 Sp1 sites, over 20 oncogene transcription factor binding motifs and numerous hormonal and other specialized binding motifs. Knowledge of 5′ and 3′ chTERC regulatory elements will be useful for investigating normal control mechanisms during growth and development as well as investigating the potential for dysregulation of this important gene during oncogenesis, especially among different genotypes.

The nucleolus organizer region (NOR) encodes the tandemly repeated 18S, 5 ·8S and 28S ribosomal (r) RNA genes. The NORs of broiler and layer commercial chicken pure lines were studied to establish the type and extent of genetic variation at this important locus. The parameters studied were gene copy number, repeat size, and diversity of NOR-types. The populations were organized into three groups for analysis including brown-egg broiler (13 lines), brown-egg layer (six lines), and white-egg layer (eight lines). The ribosomal gene copy number average of the white-egg layer populations was significantly lower (329 genes) than that of the brown-egg layers (372 genes); the brown-egg broiler ribosomal gene average was intermediate (350 genes). The white-egg layer populations exhibited a ribosomal repeat unit average size of 36 kb, significantly different from the brown-egg layer and brown-egg broiler average repeat unit size of 32 ·5 and 33 ·9 kb, respectively. NOR array size was similar among the three groups (6 mb). The brown-egg broiler populations exhibited polymorphic NOR patterns, intra-and interline, whereas the white-egg layer populations were essentially monomorphic for NOR-type; brown-egg layers exhibited an intermediate level of NOR diversity. Some NOR array characteristics may be a function of breed origin as brown-egg commercial populations, both broilers and layers, have similar breed origins and exhibited similarities for predominant repeat unit size as compared with white-egg layer populations. However, the finding that brown-egg broiler lines typically exhibit a greater number of segregating NOR-types than brown-egg layer lines suggests that the selection schemes of broiler vs. layer pure line populations may also have influenced the degree of variation at this gene complex.

The molecular organization of the 18S, 5.8S, and 28S ribosomal RNA gene repeat units, located at the single nucleolus organizer region (NOR) locus in the chicken, was investigated in genetically distinct populations of research and commercial chickens. Substantial gene repeat variation within and among NORs was documented. Intact ribosomal gene repeat size ranged from 11 kb to over 50 kb. Unique combinations of ribosomal genes, of different size, were specific to particular populations. It was determined that the basis for the ribosomal gene repeat size variation was intergenic spacer (IGS) length heterogeneity. Interestingly, in different populations, the location of the variation that contributes to length heterogeneity was specific to particular IGS subregions. In addition to IGS variation, an inbred line of Red Jungle Fowl exhibited coding region variation. Ribosomal gene copy number variation was also studied, and line averages ranged from 279 to 368. Average rDNA array size (a function of copy number and gene repeat length) was calculated for each of the populations and found to vary over a range of two megabases, from 5 to 7 Mb.Key words: rDNA, NOR, IGS, genetic variation, chicken.

Ribosomal (r) DNA genotypes (rRNA gene copy number) and nucleolar phenotypes (nucleoli number and size) were studied in dam and sire commercial broiler pure lines from three primary breeder sources. Thirteen lines were studied to determine whether directionally selected broiler pure lines contain higher numbers of rRNA genes than a control line unselected for performance traits. Eight of the 13 lines exhibited rRNA gene copy averages between 261 and 331 copies, three lines had averages between 365 and 380, and two lines had average copy numbers equal to or greater than 450 rRNA genes. The overall source copy number average from one breeder company exhibited a value (402 rRNA genes) significantly different from the control value (300 rRNA genes). Nucleoli number and relative-size were examined in 9 of the 13 lines to establish ploidy and determine the population incidence of nucleolar size polymorphisms. All of the individuals examined for nucleolar phenotype expressed two nucleoli, indicating that gene copy number variation in those lines was generally unrelated to haploidy, aneuploidy, or polyploidy. A high frequency of individuals exhibited nucleolar size polymorphisms (line values of 57 to 87%). The results suggest that multiple nucleolus organizer region (NOR) types are segregating within and among broiler pure lines and that these NOR types contain variable numbers of rRNA genes that differ in nucleogenesis capacity.

A new strain of chickens (mPNU) that segregates a severely deleted rDNA cluster was studied. Individuals heterozygous (+/p2) and homozygous (p2/p2) for the deletion were found to have 56 and 27%, respectively, of the normal complement of rRNA genes (290 copies/cell). Morphogenesis, cellular rRNA levels, and nucleolar sizes, were investigated and compared in normal +/+, +/p2, and p2/p2 embryos. Cellular rRNA contents were similar among the three genotypes at stage X, but subsequently during gastrulation, p2/p2 levels were reduced to 56% and eventually to 43% of +/+. Viability and morphogenesis were normal in p2/p2 embryos until the initial primitive streak stage of gastrulation. However, further development was abnormal and characterized by disrupted axis formation. In +/+ and +/p2 embryos, rRNA levels and nucleolar sizes increased during early development; however, the profile of these increases differed temporally and quantitatively between the genotypes. The +/p2 embryos, at the full streak stage of gastrulation, exhibited reduced rRNA levels and nucleolar sizes (80% of +/+), yet the +/p2 embryos developed normally. These studies establish a minimum copy number requirement lower than previously demonstrated, that is, a rDNA genotype with only 56% of the normal gene complement (∼160 genes) is compatible with early embryonic viability. Also, a rRNA threshold was detected: rRNA levels that were 56% of +/+ failed to support normal gastrulation; however, even under the circumstance of reduced rRNA levels (43% of control), some aspects of gastrulation apparently continue (cell migration and invagination). The teratogenic development of p2/p2 embryos is a biological consequence unique from that found in other metazoan models of rDNA-deficiency, and will be useful as a model to investigate mechanisms of vertebrate gastrulation and axis formation.

Chicken developmental mutants are valuable for discovering sequences and pathways controlling amniote development. Herein we applied the advanced technologies of targeted sequence genomic capture enrichment and next-generation sequencing to discover the causative element for three inherited mutations affecting craniofacial, limb and/or organ development. Since the mutations (coloboma, diplopodia-1 and wingless-2) were bred into a congenic line series and previously mapped to different chromosomes, each targeted mutant causative region could be compared to that of the other two congenic partners, thereby providing internal controls on a single array. Of the ~73 million 50-bp sequence reads, ~76% were specific to the enriched targeted regions with an average target coverage of 132-fold. Analysis of the three targeted regions (2.06 Mb combined) identified line-specific single nucleotide polymorphism (SNPs) and micro (1–3 nt) indels. Sequence content for regions indicated as gaps in the reference genome was generated, thus contributing to its refinement. Additionally, Mauve alignments were constructed and indicated putative chromosomal rearrangements. This is the first report of targeted capture array technology in an avian species, the chicken, an important vertebrate model, the work highlights the utility of employing advanced technologies in an organism with only a “draft stage” reference genome sequence.