Poultry genome sequences: progress and outstanding challenges

Genome-wide association study of copy number variations with shank traits in a F2 crossbred chicken population

Copy number variations (CNVs) are large-scale changes in the DNA sequence that can affect the genetic structure and phenotype of an organism. The purpose of this study was to investigate the existing CNVs and their associations with the shank diameter (ShD) and shank length (ShL) traits using data from an F2 crossbred chicken population. To carry out the study, 312 chickens were genotyped using the Illumina 60k SNP Beadchip. The shank traits of the birds were measured from day 1 to 12 weeks of age. penncnv and cnvruler tools were used to find copy numbers and regions with copy number changes (CNVR), respectively. The CNVRanger package was used to perform a genome-wide association study between shank traits and CNVs. Gene ontology research in CNVRs was carried out using the david database. In this investigation, 966 CNVs and 606 regions with copy number changes were discovered. The copy number states and variations were randomly distributed along the length of the autosomal chromosomes. Weeks 1-4, 9 and 12 of growth revealed a significant association of copy number variations with shank traits, false discovery rate (FDR-corrected p-value < 0.01), and the majority of CNVs that were statistically significant were found on chromosomes 1-3. These CNV segments are nearby genes such as KCNJ12, FGF6 and MYF5, which are fundamental to growth and development. In addition, gene set analyses revealed terms related to muscle physiology, regulation of cellular processes and potassium channels.

A Bird's-Eye View of Chromosomic Evolution in the Class Aves

Birds (Aves) are the most speciose of terrestrial vertebrates, displaying Class-specific characteristics yet incredible external phenotypic diversity. Critical to agriculture and as model organisms, birds have adapted to many habitats. The only extant examples of dinosaurs, birds emerged ~150 mya and >10% are currently threatened with extinction. This review is a comprehensive overview of avian genome ("chromosomic") organization research based mostly on chromosome painting and BAC-based studies. We discuss traditional and contemporary tools for reliably generating chromosome-level assemblies and analyzing multiple species at a higher resolution and wider phylogenetic distance than previously possible. These results permit more detailed investigations into inter- and intrachromosomal rearrangements, providing unique insights into evolution and speciation mechanisms. The 'signature' avian karyotype likely arose ~250 mya and remained largely unchanged in most groups including extinct dinosaurs. Exceptions include Psittaciformes, Falconiformes, Caprimulgiformes, Cuculiformes, Suliformes, occasional Passeriformes, Ciconiiformes, and Pelecaniformes. The reasons for this remarkable conservation may be the greater diploid chromosome number generating variation (the driver of natural selection) through a greater possible combination of gametes and/or an increase in recombination rate. A deeper understanding of avian genomic structure permits the exploration of fundamental biological questions pertaining to the role of evolutionary breakpoint regions and homologous synteny blocks.

Identification of the genetic characteristics of copy number variations in experimental specific pathogen-free ducks using whole-genome resequencing

BACKGROUND

Specific pathogen-free ducks are a valuable laboratory resource for waterfowl disease research and poultry vaccine development. High throughput sequencing allows the systematic identification of structural variants in genomes. Copy number variation (CNV) can explain the variation of important duck genetic traits. Herein, the genome-wide CNVs of the three experimental duck species in China (Jinding ducks (JD), Shaoxing ducks (SX), and Fujian Shanma ducks (SM)) were characterized using resequencing to determine their genetic characteristics and selection signatures.

RESULTS

We obtained 4,810 CNV regions (CNVRs) by merging 73,012 CNVs, covering 4.2% of the duck genome. Functional analysis revealed that the shared CNVR-harbored genes were significantly enriched for 31 gene ontology terms and 16 Kyoto Encyclopedia of Genes and Genomes pathways (e.g., olfactory transduction and immune system). Based on the genome-wide fixation index for each CNVR, growth (SPAG17 and PTH1R), disease resistance (CATHL3 and DMBT1), and thermoregulation (TRPC4 and SLIT3) candidate genes were identified in strongly selected signatures specific to JD, SM, and SX, respectively.

CONCLUSIONS

In conclusion, we investigated the genome-wide distribution of experimental duck CNVs, providing a reference to establish the genetic basis of different phenotypic traits, thus contributing to the management of experimental animal genetic resources.

Whole-genome resequencing to unveil genetic characteristics and selection signatures of specific pathogen-free ducks

Specific pathogen-free ducks are important high-grade laboratory animals, with a key role in research related to poultry biosecurity, production, and breeding. However, the genetic characteristics of experimental duck varieties remain poorly explored. Herein we performed whole-genome resequencing to construct a single nucleotide polymorphism genetic map of the genomes of 3 experimental duck varieties [Jinding ducks (JD), Shaoxing ducks (SX), and Fujian Shanma ducks (SM)] to determine their genetic characteristics and identify selection signatures. Subsequent analyses of population structure and genetic diversity revealed that each duck variety formed a monophyletic group, with SM showing richer genetic diversity than JD and SX. Further, on exploring shared selection signatures, we found 2 overlapping genomic regions on chromosome Z of all experimental ducks, which comprised immune response-related genes (IL7R and IL6ST). Moreover, growth and skeletal development (IGF1R and GDF5), meat quality (FoxO1), and stress resistance (HSP90B1 and Gpx8-b) candidate gene loci were identified in strongly selected signatures specific to JD, SM, and SX, respectively. Our results identified the population genetic basis of experimental ducks at the whole-genome level, providing a framework for future molecular investigations of genetic variations and phenotypic changes. We believe that such studies will eventually contribute to the management of experimental animal resources.

The science of genetically modified poultry

The exuberant development of targeted genome editing has revolutionized research on the chicken genome, generating chickens with beneficial parameters. The chicken model is a crucial experimental tool that can be utilized for drug manufacture, preclinical research, pathological observation, and other applications. In essence, tweaking the chicken’s genome has enabled the poultry industry to get more done with less, generating genetically modified chickens that lay eggs containing large amounts of lifesaving humanized drugs. The transition of gene editing from concept to practical application has been dramatically hastened by the development of programmable nucleases, bringing scientists closer than ever to the efficient producers of tomorrow’s medicines. Combining the developmental and physiological characteristics of the chicken with cutting-edge genome editing, the chicken furnishes a potent frontier that is foreseen to be actively pursued in the future. Herein we review the current and future prospects of gene editing in chickens and the contributions to the development of humanized pharmaceuticals.

The Innovative Informatics Approaches of High-Throughput Technologies in Livestock: Spearheading the Sustainability and Resiliency of Agrigenomics Research

For more than a decade, next-generation sequencing (NGS) has been emerging as the mainstay of agrigenomics research. High-throughput technologies have made it feasible to facilitate research at the scale and cost required for using this data in livestock research. Scale frameworks of sequencing for agricultural and livestock improvement, management, and conservation are partly attributable to innovative informatics methodologies and advancements in sequencing practices. Genome-wide sequence-based investigations are often conducted worldwide, and several databases have been created to discover the connections between worldwide scientific accomplishments. Such studies are beginning to provide revolutionary insights into a new era of genomic prediction and selection capabilities of various domesticated livestock species. In this concise review, we provide selected examples of the current state of sequencing methods, many of which are already being used in animal genomic studies, and summarize the state of the positive attributes of genome-based research for cattle (Bos taurus), sheep (Ovis aries), pigs (Sus scrofa domesticus), horses (Equus caballus), chickens (Gallus gallus domesticus), and ducks (Anas platyrhyncos). This review also emphasizes the advantageous features of sequencing technologies in monitoring and detecting infectious zoonotic diseases. In the coming years, the continued advancement of sequencing technologies in livestock agrigenomics will significantly influence the sustained momentum toward regulatory approaches that encourage innovation to ensure continued access to a safe, abundant, and affordable food supplies for future generations.

Assemblies of the genomes of parasitic wasps using meta-assembly and scaffolding with genetic linkage

Safe, effective biological-control introductions against invasive pests depend on narrowly host-specific natural enemies with the ability to adapt to a changing environment. As part of a project on the genetic architectures of these traits, we assembled and annotated the genomes of two aphid parasitoids, Aphelinus atriplicis and Aphelinus certus. We report here several assemblies of A. atriplicis made with Illumina and PacBio data, which we combined into a meta-assembly. We scaffolded the meta-assembly with markers from a genetic map of hybrids between A. atriplicis and A. certus. We used this genetic-linkage scaffolded (GLS) assembly of A. atriplicis to scaffold a de novo assembly of A. certus. The de novo assemblies of A. atriplicis differed in contiguity, and the meta-assembly of these assemblies was more contiguous than the best de novo assembly. Scaffolding with genetic-linkage data allowed chromosomal-level assembly of the A. atriplicis genome and scaffolding a de novo assembly of A. certus with this GLS assembly, greatly increased the contiguity of the A. certus assembly to the point where it was also at the chromosomal-level. However, completeness of the A. atriplicis assembly, as measured by percent complete, single-copy BUSCO hymenopteran genes, varied little among de novo assemblies and was not increased by meta-assembly or genetic scaffolding. Furthermore, the greater contiguity of the meta-assembly and GLS assembly had little or no effect on the numbers of genes identified, the proportions with homologs or functional annotations. Increased contiguity of the A. certus assembly provided modest improvement in assembly completeness, as measured by percent complete, single-copy BUSCO hymenopteran genes. The total genic sequence increased, and while the number of genes declined, gene length increased, which together suggest greater accuracy of gene models. More contiguous assemblies provide uses other than gene annotation, for example, identifying the genes associated with quantitative trait loci and understanding of chromosomal rearrangements associated with speciation.

A Missense Mutation in the MYBPH Gene Is Associated With Abdominal Fat Traits in Meat-Type Chickens

Chicken is an important source of protein for human nutrition and a model system for growth and developmental biology. Although the genetic architecture of quantitative traits in meat-type chickens has been the subject of ongoing investigation, the identification of mutations associated with carcass traits of economic interest remains challenging. Therefore, our aim was to identify predicted deleterious mutation, which potentially affects protein function, and test if they were associated with carcass traits in chickens. For that, we performed a genome-wide association analysis (GWAS) for breast, thigh and drumstick traits in meat-type chickens and detected 19 unique quantitative trait loci (QTL). We then used: (1) the identified windows; (2) QTL for abdominal fat detected in a previous study with the same population and (3) previously obtained whole genome sequence data, to identify 18 predicted deleterious single nucleotide polymorphisms (SNPs) in those QTL for further association with breast, thigh, drumstick and abdominal fat traits. Using the additive model, a predicted deleterious SNP c.482C > T (SIFT score of 0.4) was associated (p-value < 0.05) with abdominal fat weight and percentage. This SNP is in the second exon of the MYBPH gene, and its allele frequency deviates from Hardy–Weinberg equilibrium. In conclusion, our study provides evidence that the c.482C > T SNP in the MYBPH gene is a putative causal mutation for fat deposition in meat-type chickens.

Genome-wide detection of CNVs and their association with performance traits in broilers

Background

Copy number variations (CNVs) are a major type of structural genomic variants that underlie genetic architecture and phenotypic variation of complex traits, not only in humans, but also in livestock animals. We identified CNVs along the chicken genome and analyzed their association with performance traits. Genome-wide CNVs were inferred from Affymetrix® high density SNP-chip data for a broiler population. CNVs were concatenated into segments and association analyses were performed with linear mixed models considering a genomic relationship matrix, for birth weight, body weight at 21, 35, 41 and 42 days, feed intake from 35 to 41 days, feed conversion ratio from 35 to 41 days and, body weight gain from 35 to 41 days of age.

Results

We identified 23,214 autosomal CNVs, merged into 5042 distinct CNV regions (CNVRs), covering 12.84% of the chicken autosomal genome. One significant CNV segment was associated with BWG on GGA3 (q-value = 0.00443); one significant CNV segment was associated with BW35 (q-value = 0.00571), BW41 (q-value = 0.00180) and BW42 (q-value = 0.00130) on GGA3, and one significant CNV segment was associated with BW on GGA5 (q-value = 0.00432). All significant CNV segments were verified by qPCR, and a validation rate of 92.59% was observed. These CNV segments are located nearby genes, such as KCNJ11, MyoD1 and SOX6, known to underlie growth and development. Moreover, gene-set analyses revealed terms linked with muscle physiology, cellular processes regulation and potassium channels.

Conclusions

Overall, this CNV-based GWAS study unravels potential candidate genes that may regulate performance traits in chickens. Our findings provide a foundation for future functional studies on the role of specific genes in regulating performance in chickens.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07676-1.

High-Density Linkage Maps Based on Genotyping-by-Sequencing (GBS) Confirm a Chromosome-Level Genome Assembly and Reveal Variation in Recombination Rate for the Pacific Oyster Crassostrea gigas

Studies of linkage and linkage mapping have advanced genetic and biological knowledge for over 100 years. In addition to their growing role, today, in mapping phenotypes to genotypes, dense linkage maps can help to validate genome assemblies. Previously, we showed that 40% of scaffolds in the first genome assembly for the Pacific oyster Crassostrea gigas were chimeric, containing single nucleotide polymorphisms (SNPs) mapping to different linkage groups. Here, we merge 14 linkage maps constructed of SNPs generated from genotyping-by-sequencing (GBS) methods with five, previously constructed linkage maps, to create a compendium of nearly 69 thousand SNPs mapped with high confidence. We use this compendium to assess a recently available, chromosome-level assembly of the C. gigas genome, mapping SNPs in 275 of 301 contigs and comparing the ordering of these contigs, by linkage, to their assembly by Hi-C sequencing methods. We find that, while 26% of contigs contain chimeric blocks of SNPs, i.e., adjacent SNPs mapping to different linkage groups than the majority of SNPs in their contig, these apparent misassemblies amount to only 0.08% of the genome sequence. Furthermore, nearly 90% of 275 contigs mapped by linkage and sequencing are assembled identically; inconsistencies between the two assemblies for the remaining 10% of contigs appear to result from insufficient linkage information. Thus, our compilation of linkage maps strongly supports this chromosome-level assembly of the oyster genome. Finally, we use this assembly to estimate, for the first time in a Lophotrochozoan, genome-wide recombination rates and causes of variation in this fundamental process.

Discovery of high affinity, pan-allelic, and pan-mammalian reactive antibodies against the myeloid checkpoint receptor SIRPα

Targeting the CD47-signal-regulatory protein α (SIRPα) pathway represents a novel therapeutic approach to enhance anti-cancer immunity by promoting both innate and adaptive immune responses. Unlike CD47, which is expressed ubiquitously, SIRPα expression is mainly restricted to myeloid cells and neurons. Therefore, compared to CD47-targeted therapies, targeting SIRPα may result in differential safety and efficacy profiles, potentially enabling lower effective doses and improved pharmacokinetics and pharmacodynamics. The development of effective SIRPα antagonists is restricted by polymorphisms within the CD47-binding domain of SIRPα, necessitating pan-allele reactive anti-SIRPα antibodies for therapeutic intervention in diverse patient populations. We immunized wild-type and human antibody transgenic chickens with a multi-allele and multi-species SIRPα regimen in order to discover pan-allelic and pan-mammalian reactive anti-SIRPα antibodies suitable for clinical translation. A total of 200 antibodies were isolated and screened for SIRPα reactivity from which approximately 70 antibodies with diverse SIRPα binding profiles, sequence families, and epitopes were selected for further characterization. A subset of anti-SIRPα antibodies bound to both human SIRPα v1 and v2 alleles with high affinity ranging from low nanomolar to picomolar, potently antagonized the CD47/SIRPα interaction, and potentiated macrophage-mediated antibody-dependent cellular phagocytosis in vitro. X-ray crystal structures of five anti-SIRPα antigen-binding fragments, each with unique epitopes, in complex with SIRPα (PDB codes 6NMV, 6NMU, 6NMT, 6NMS, and 6NMR) are reported. Furthermore, some of the anti-SIRPα antibodies cross-react with cynomolgus SIRPα and various mouse SIRPα alleles (BALB/c, NOD, BL/6), which can facilitate preclinical to clinical development. These properties provide an attractive rationale to advance the development of these anti-SIRPα antibodies as a novel therapy for advanced malignancies. Abbreviations: ADCC: antibody-dependent cellular cytotoxicity; ADCP: antibody-dependent cellular phagocytosis; CFSE: carboxyfluorescein succinimidyl ester; Fab: fragment antigen binding; Fc: fragment crystallizable; FcγR: Fcγ receptor; Ig: immunoglobulin; IND: investigational new drug; MDM⊘: monocyte-derived macrophage; NOD: non-obese diabetic; scFv: single chain fragment variable; SCID: severe combined immunodeficiency; SIRP: signal-regulatory protein.

Chickens with humanized immunoglobulin genes generate antibodies with high affinity and broad epitope coverage to conserved targets

Transgenic animal platforms for the discovery of human monoclonal antibodies have been developed in mice, rats, rabbits and cows. The immune response to human proteins is limited in these animals by their tolerance to mammalian-conserved epitopes. To expand the range of epitopes that are accessible, we have chosen an animal host that is less phylogenetically related to humans. Specifically, we generated transgenic chickens expressing antibodies from immunoglobulin heavy and light chain loci containing human variable regions and chicken constant regions. From these birds, paired human light and heavy chain variable regions are recovered and cloned as fully human recombinant antibodies. The human antibody-expressing chickens exhibit normal B cell development and raise immune responses to conserved human proteins that are not immunogenic in mice. Fully human monoclonal antibodies can be recovered with sub-nanomolar affinities. Binning data of antibodies to a human protein show epitope coverage similar to wild type chickens, which we previously showed is broader than that produced from rodent immunizations.

Utilizing the chicken as an animal model for human craniofacial ciliopathies

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.

Second-Generation Linkage Maps for the Pacific Oyster Crassostrea gigas Reveal Errors in Assembly of Genome Scaffolds

The Pacific oyster Crassostrea gigas, a widely cultivated marine bivalve mollusc, is becoming a genetically and genomically enabled model for highly fecund marine metazoans with complex life-histories. A genome sequence is available for the Pacific oyster, as are first-generation, low-density, linkage and gene-centromere maps mostly constructed from microsatellite DNA makers. Here, higher density, second-generation, linkage maps are constructed from more than 1100 coding (exonic) single-nucleotide polymorphisms (SNPs), as well as 66 previously mapped microsatellite DNA markers, all typed in five families of Pacific oysters (nearly 172,000 genotypes). The map comprises 10 linkage groups, as expected, has an average total length of 588 cM, an average marker-spacing of 1.0 cM, and covers 86% of a genome estimated to be 616 cM. All but seven of the mapped SNPs map to 618 genome scaffolds; 260 scaffolds contain two or more mapped SNPs, but for 100 of these scaffolds (38.5%), the contained SNPs map to different linkage groups, suggesting widespread errors in scaffold assemblies. The 100 misassembled scaffolds are significantly longer than those that map to a single linkage group. On the genetic maps, marker orders and intermarker distances vary across families and mapping methods, owing to an abundance of markers segregating from only one parent, to widespread distortions of segregation ratios caused by early mortality, as previously observed for oysters, and to genotyping errors. Maps made from framework markers provide stronger support for marker orders and reasonable map lengths and are used to produce a consensus high-density linkage map containing 656 markers.

Cloning of a functional 25-hydroxyvitamin D-1α-hydroxylase in zebrafish (Danio rerio)

Activation of precursor 25-hydroxyvitamin D3 (25D) to hormonal 1,25-dihydroxyvitamin D3 (1,25D) is a pivotal step in vitamin D physiology, catalysed by the enzyme 25-hydroxyvitamin D-1α-hydroxylase (1α-hydroxylase). To establish new models for assessing the physiological importance of the 1α-hydroxylase-25D-axis, we used Danio rerio (zebrafish) to characterize expression and biological activity of the gene for 1α-hydroxylase (cyp27b1). Treatment of day 5 zebrafish larvae with inactive 25D (5-150 nM) or active 1,25D (0.1-10 nM) induced dose responsive expression (15-95-fold) of the vitamin D-target gene cyp24a1 relative to larvae treated with vehicle, suggesting the presence of Cyp27b1 activity. A full-length zebrafish cyp27b1 cDNA was then generated using RACE and RT-PCR methods. Sequencing of the resulting clone revealed an open reading frame encoding a protein of 505 amino acids with 54% identity to human CYP27B1. Transfection of a cyp27b1 expression vector into HKC-8, a human kidney proximal tubular epithelial cell line, enhanced intracrine metabolism of 25D to 1,25D resulting in greater than twofold induction of CYP24A1 mRNA expression and a 25-fold increase in 1,25D production compared to empty vector. These data indicate that we have cloned a functional zebrafish CYP27B1, representing a phylogenetically distant branch from mammals of this key enzyme in vitamin D metabolism. Further analysis of cyp27b1 expression and activity in zebrafish may provide new perspectives on the biological importance of 25D metabolism.

Applied animal genomics: results from the field

Genomic selection (GS) is the use of statistical methods to estimate the genetic merit of a genotyped animal based on prediction equations derived from large ancestral populations with both phenotypes and genotypes. It has revolutionized the dairy cattle breeding industry and has been implemented with varying degrees of success in other animal breeding programs, including swine, poultry, and beef cattle. The findings of empirical field studies applying GS to the breeding sectors of these main animal protein industries are reviewed. Several translational considerations must be addressed before implementing GS in genetic improvement programs. These include determining and obtaining economically relevant phenotypes and determining the optimal size of the training population, cost-effective genotyping strategies, the practicality of field implementation, and the relative costs versus the benefits of the realized rate of genetic gain. GS may additionally change the optimal breeding scheme design, and studies that address this consideration are also reviewed briefly.

The proteome of the calcified layer organic matrix of turkey (Meleagris gallopavo) eggshell

Background

Chicken eggshell mineralization is a prominent model for biomineralization not only because of its importance for avian reproduction but also because of the commercial interest associated with eggshell quality. An analysis and comparison of the protein constituents of eggshells of several species would contribute to a better understanding of the shell mineralization process. The recent publication of the turkey genome sequence now provides a basis for the in-depth analysis of the turkey eggshell proteome.

Results

Proteomic analysis of turkey acid-soluble and acid-insoluble organic eggshell matrix yielded 697 identified proteins/protein groups. However, intensity-based absolute quantification (iBAQ) results indicated that the 47 most abundant identified proteins already constituted 95% of the total turkey eggshell matrix proteome. Forty-four of these proteins were also identified in chicken eggshell matrix previously. Despite these similarities there were important and unexpected differences. While ovocleidin-116 and ovocalyxin-36 were major proteins constituting approximately 37% of the identified proteome, other members of the group of so-called eggshell-specific proteins were not identified. Thus ovocalyxin-21 and ovocalyxin-32 were missing among matrix proteins. Conversely, major turkey eggshell proteins were not detected in chicken, such as the bone protein periostin, the mammalian counterpart of which is involved in many aspects of bone metabolism and which represented 10-11% of the total identified proteome.

Conclusions

Even members of the same avian family show important differences in eggshell matrix composition and more studies on the proteome and the transcriptome level will be necessary to identify a common toolkit of eggshell mineralization and to work out species differences among functional eggshell protein sets and their role in eggshell production.

Evolution of the let-7 microRNA family

The increase of bodyplan complexity in early bilaterian evolution is correlates with the advent and diversification of microRNAs. These small RNAs guide animal development by regulating temporal transitions in gene expression involved in cell fate choices and transitions between pluripotency and differentiation. One of the two known microRNAs whose origins date back before the bilaterian ancestor is mir-100. In Bilateria, it appears stably associated in polycistronic transcripts with let-7 and mir-125, two key regulators of development. In vertebrates, these three microRNA families have expanded to form a complex system of developmental regulators. In this contribution, we disentangle the evolutionary history of the let-7 locus, which was restructured independently in nematodes, platyhelminths, and deuterostomes. The foundation of a second let-7 locus in the common ancestor of vertebrates and urochordates predates the vertebrate-specific genome duplications, which then caused a rapid expansion of the let-7 family.

Accelerated evolution of 3'avian FOXE1 genes, and thyroid and feather specific expression of chicken FoxE1

Background

The forkhead transcription factor gene E1 (FOXE1) plays an important role in regulation of thyroid development, palate formation and hair morphogenesis in mammals. However, avian FOXE1 genes have not been characterized and as such, codon evolution of FOXE1 orthologs in a broader evolutionary context of mammals and birds is not known.

Results

In this study we identified the avian FOXE1 gene in chicken, turkey and zebra finch, all of which consist of a single exon. Chicken and zebra finch FOXE1 are uniquely located on the sex-determining Z chromosome. In situ hybridization shows that chicken FOXE1 is specifically expressed in the developing thyroid. Its expression is initiated at the placode stage and is maintained during the stages of vesicle formation and follicle primordia. Based on this expression pattern, we propose that avian FOXE1 may be involved in regulating the evagination and morphogenesis of thyroid. Chicken FOXE1 is also expressed in growing feathers. Sequence analysis identified two microdeletions in the avian FOXE1 genes, corresponding to the loss of a transferable repression domain and an engrailed homology motif 1 (Eh1) C-terminal to the forkhead domain. The avian FOXE1 proteins exhibit a significant sequence divergence of the C-terminus compared to those of amphibian and mammalian FOXE1. The codon evolution analysis (dN/dS) of FOXE1 shows a significantly increased dN/dS ratio in the avian lineages, consistent with either a relaxed purifying selection or positive selection on a few residues in avian FOXE1 evolution. Further site specific analysis indicates that while relaxed purifying selection is likely to be a predominant cause of accelerated evolution at the 3'-region of avian FOXE1, a few residues might have evolved under positive selection.

Conclusions

We have identified three avian FOXE1 genes based on synteny and sequence similarity as well as characterized the expression pattern of the chicken FOXE1 gene during development. Our evolutionary analyses suggest that while a relaxed purifying selection is likely to be the dominant force driving accelerated evolution of avian FOXE1 genes, a few residues may have evolved adaptively. This study provides a basis for future genetic and comparative biochemical studies of FOXE1.

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.