Multi-platform next-generation sequencing of the domestic turkey (Meleagris gallopavo): genome assembly and analysis

Ensembl 2012

The Ensembl project (http://www.ensembl.org) provides genome resources for chordate genomes with a particular focus on human genome data as well as data for key model organisms such as mouse, rat and zebrafish. Five additional species were added in the last year including gibbon (Nomascus leucogenys) and Tasmanian devil (Sarcophilus harrisii) bringing the total number of supported species to 61 as of Ensembl release 64 (September 2011). Of these, 55 species appear on the main Ensembl website and six species are provided on the Ensembl preview site (Pre!Ensembl; http://pre.ensembl.org) with preliminary support. The past year has also seen improvements across the project.

Emergence and evolution of the glycoprotein hormone and neurotrophin gene families in vertebrates

Background

The three vertebrate pituitary glycoprotein hormones (GPH) are heterodimers of a common α and a specific β subunit. In human, they are located on different chromosomes but in a similar genomic environment. We took advantage of the availability of genomic and EST data from two cartilaginous fish species as well as from two lamprey species to identify their repertoire of neurotrophin, lin7 and KCNA gene family members which are in the close environment of gphβ. Gphα and gphβ are absent outside vertebrates but are related to two genes present in both protostomes and deuterostomes that were named gpa2 and gpb5. Genomic organization and functional characteristics of their protein products suggested that gphα and gphβ might have been generated concomitantly by a duplication of gpa2 and gpb5 just prior to the radiation of vertebrates. To have a better insight into this process we used new genomic resources and tools to characterize the ancestral environment before the duplication occurred.

Results

An almost similar repertoire of genes was characterized in cartilaginous fishes as in tetrapods. Data in lampreys are either incomplete or the result of specific duplications and/or deletions but a scenario for the evolution of this genomic environment in vertebrates could be proposed. A number of genes were identified in the amphioxus genome that helped in reconstructing the ancestral environment of gpa2 and gpb5 and in describing the evolution of this environment in vertebrates.

Conclusion

Our model suggests that vertebrate gphα and gphβ were generated by a specific local duplication of the ancestral forms of gpa2 and gpb5, followed by a translocation of gphβ to a new environment whereas gphα was retained in the gpa2-gpb5 locus. The two rounds of whole genome duplication that occurred early in the evolution of vertebrates generated four paralogues of each gene but secondary gene losses or lineage specific duplications together with genomic rearrangements have resulted in the present organization of these genes, which differs between vertebrate lineages.

Intrachromosomal rearrangements in avian genome evolution: evidence for regions prone to breakpoints

It is generally believed that the organization of avian genomes remains highly conserved in evolution as chromosome number is constant and comparative chromosome painting demonstrated there to be very few interchromosomal rearrangements. The recent sequencing of the zebra finch (Taeniopygia guttata) genome allowed an assessment of the number of intrachromosomal rearrangements between it and the chicken (Gallus gallus) genome, revealing a surprisingly high number of intrachromosomal rearrangements. With the publication of the turkey (Meleagris gallopavo) genome it has become possible to describe intrachromosomal rearrangements between these three important avian species, gain insight into the direction of evolutionary change and assess whether breakpoint regions are reused in birds. To this end, we aligned entire chromosomes between chicken, turkey and zebra finch, identifying syntenic blocks of at least 250 kb. Potential optimal pathways of rearrangements between each of the three genomes were determined, as was a potential Galliform ancestral organization. From this, our data suggest that around one-third of chromosomal breakpoint regions may recur during avian evolution, with 10% of breakpoints apparently recurring in different lineages. This agrees with our previous hypothesis that mechanisms of genome evolution are driven by hotspots of non-allelic homologous recombination.

DNA ligase III: a spotty presence in eukaryotes, but an essential function where tested

DNA ligases are crucial for most DNA transactions, including DNA replication, repair, and recombination. Recently, DNA ligase III (Lig3) has been demonstrated to be crucial for cell survival due to its catalytic function in mitochondria. This review summarizes these recent results and reports on a hitherto unappreciated widespread phylogenetic presence of Lig3 in eukaryotes, including in some organisms before the divergence of metazoa. Analysis of these putative Lig3 homologs suggests that many of them are likely to be found in mitochondria and to be critical for mitochondrial function.

Rapid in vivo testing of drug response in multiple myeloma made possible by xenograft to turkey embryos

Background:

The best current xenograft model of multiple myeloma (MM) in immune-deficient non-obese diabetic/severe-combined immunodeficient mice is costly, animal maintenance is complex and several weeks are required to establish engraftment and study drug efficacy. More practical in vivo models may reduce time and drug development cost. We recently described a rapid low-cost xenograft model of human blood malignancies in pre-immune turkey. Here, we report application of this system for studying MM growth and the preclinical assessment of anticancer therapies.

Methods:

Cell lines and MM patient cells were injected intravenously into embryonic veins on embryonic day 11 (E11). Engraftment of human cells in haematopoietic organs was detected by quantitative real-time polymerase chain reaction, immunohistochemistry, flow cytometry and circulating free light chain.

Results:

Engraftment was detected after 1 week in all embryos injected with cell lines and in 50% of those injected with patient cells. Injection of bortezomib or lenalinomide 48 h after cell injection at therapeutic levels that were not toxic to the bone marrow dramatically reduced MM engraftment.

Conclusion:

The turkey embryo provides a practical, xenograft system to study MM and demonstrates the utility of this model for rapid and affordable testing therapeutics in vivo. With further development, this model may enable rapid, inexpensive personalised drug screening.

Sequencing technologies and genome sequencing

The high-throughput - next generation sequencing (HT-NGS) technologies are currently the hottest topic in the field of human and animals genomics researches, which can produce over 100 times more data compared to the most sophisticated capillary sequencers based on the Sanger method. With the ongoing developments of high throughput sequencing machines and advancement of modern bioinformatics tools at unprecedented pace, the target goal of sequencing individual genomes of living organism at a cost of $1,000 each is seemed to be realistically feasible in the near future. In the relatively short time frame since 2005, the HT-NGS technologies are revolutionizing the human and animal genome researches by analysis of chromatin immunoprecipitation coupled to DNA microarray (ChIP-chip) or sequencing (ChIP-seq), RNA sequencing (RNA-seq), whole genome genotyping, genome wide structural variation, de novo assembling and re-assembling of genome, mutation detection and carrier screening, detection of inherited disorders and complex human diseases, DNA library preparation, paired ends and genomic captures, sequencing of mitochondrial genome and personal genomics. In this review, we addressed the important features of HT-NGS like, first generation DNA sequencers, birth of HT-NGS, second generation HT-NGS platforms, third generation HT-NGS platforms: including single molecule Heliscope™, SMRT™ and RNAP sequencers, Nanopore, Archon Genomics X PRIZE foundation, comparison of second and third HT-NGS platforms, applications, advances and future perspectives of sequencing technologies on human and animal genome research.

Expression profiling without genome sequence information in a non-model species, Pandalid shrimp (Pandalus latirostris), by next-generation sequencing

While the study of phenotypic variation is a central theme in evolutionary biology, the genetic approaches available to understanding this variation are usually limited because of a lack of genomic information in non-model organisms. This study explored the utility of next-generation sequencing (NGS) technologies for studying phenotypic variations between 2 populations of a non-model species, the Hokkai shrimp (Pandalus latirostris; Decapoda, Pandalidae). Before we performed transcriptome analyses using NGS, we examined the genetic and phenotypic differentiation between the populations. Analyses using microsatellite DNA markers suggested that these populations genetically differed from one another and that gene flow is restricted between them. Moreover, the results of our 4-year field observations indicated that the egg traits varied genetically between the populations. Using mRNA extracted from the ovaries of 5 females in each population of Hokkai shrimp, we then performed a transcriptome analysis of the 2 populations. A total of 13.66 gigabases (Gb) of 75-bp reads was obtained. Further, 58,804 and 33,548 contigs for the first and second population, respectively, and 47,467 contigs for both populations were produced by de novo assembly. We detected 552 sequences with the former approach and 702 sequences with the later one; both sets of sequences showed greater than twofold differences in the expression levels between the 2 populations. Twenty-nine sequences were found in both approaches and were considered to be differentially expressed genes. Among them, 9 sequences showed significant similarity to functional genes. The present study showed a de novo assembly approach for the transcriptome of a non-model species using only short-read sequence data, and provides a strategy for identifying sequences showing significantly different expression levels between populations.

Short reads, circular genome: skimming solid sequence to construct the bighorn sheep mitochondrial genome

As sequencing technology improves, an increasing number of projects aim to generate full genome sequence, even for nonmodel taxa. These projects may be feasibly conducted at lower read depths if the alignment can be aided by previously developed genomic resources from a closely related species. We investigated the feasibility of constructing a complete mitochondrial (mt) genome without preamplification or other targeting of the sequence. Here we present a full mt genome sequence (16,463 nucleotides) for the bighorn sheep (Ovis canadensis) generated though alignment of SOLiD short-read sequences to a reference genome. Average read depth was 1240, and each base was covered by at least 36 reads. We then conducted a phylogenomic analysis with 27 other bovid mitogenomes, which placed bighorn sheep firmly in the Ovis clade. These results show that it is possible to generate a complete mitogenome by skimming a low-coverage genomic sequencing library. This technique will become increasingly applicable as the number of taxa with some level of genome sequence rises.

Mosquito genomics: progress and challenges

The whole-genome sequencing of mosquitoes has facilitated our understanding of fundamental biological processes at their basic molecular levels and holds potential for application to mosquito control and prevention of mosquito-borne disease transmission. Draft genome sequences are available for Anopheles gambiae, Aedes aegypti, and Culex quinquefasciatus. Collectively, these represent the major vectors of African malaria, dengue fever and yellow fever viruses, and lymphatic filariasis, respectively. Rapid advances in genome technologies have revealed detailed information on genome architecture as well as phenotype-specific transcriptomics and proteomics. These resources allow for detailed comparative analyses within and across populations as well as species. Next-generation sequencing technologies will likely promote a proliferation of genome sequences for additional mosquito species as well as for individual insects. Here we review the current status of genome research in mosquitoes and identify potential areas for further investigations.

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.

The genome of the green anole lizard and a comparative analysis with birds and mammals

The evolution of the amniotic egg was one of the great evolutionary innovations in the history of life, freeing vertebrates from an obligatory connection to water and thus permitting the conquest of terrestrial environments¹. Among amniotes, genome sequences are available for mammals² and birds^3–5, but not for non-avian reptiles. Here we report the genome sequence of the North American green anole lizard, Anolis carolinensis. We find that A. carolinensis microchromosomes are highly syntenic with chicken microchromosomes, yet do not exhibit the high GC and low repeat content that are characteristic of avian microchromosomes³. Also, A. carolinensis mobile elements are very young and diverse – more so than in any other sequenced amniote genome. This lizard genome’s GC content is also unusual in its homogeneity, unlike the regionally variable GC content found in mammals and birds⁶. We describe and assign sequence to the previously unknown A. carolinensis X chromosome. Comparative gene analysis shows that amniote egg proteins have evolved significantly more rapidly than other proteins. An anole phylogeny resolves basal branches to illuminate the history of their repeated adaptive radiations.

Emergence of male-biased genes on the chicken Z-chromosome: sex-chromosome contrasts between male and female heterogametic systems

There has been extensive traffic of male-biased genes out of the mammalian and Drosophila X-chromosomes, and there are also reports of an under-representation of male-biased genes on the X. This may reflect an adaptive process driven by natural selection where an autosomal location of male-biased genes is favored since male genes are only exposed to selection one-third of the time when X-linked. However, there are several alternative explanations to "out-of-the-X" gene movement, including mutational bias and a means for X-linked genes to escape meiotic sex chromosome inactivation (MSCI) during spermatogenesis. As a critical test of the hypothesis that genomic relocation of sex-biased genes is an adaptive process, I examined the emergence, and loss, of genes on the chicken Z-chromosome, i.e., a female heterogametic system (males ZZ, females ZW). Here, the analogous prediction would be an emergence of male-biased genes onto, not a loss from, the Z-chromosome because Z is found more often in males than autosomes are. I found that genes expressed in testis but not in ovary are highly over-represented among genes that have emerged on the Z-chromosome during avian evolution. Moreover, genes with male-biased expression are similarly over-represented among new Z-chromosomal genes. Interestingly, genes with female-biased expression have more often moved from than to the Z-chromosome. These observations show that male and female heterogametic organisms display opposing directionalities in the emergence and loss of sex-biased genes on sex chromosomes. This is consistent with theoretical models on the evolution of sexually antagonistic genes in which new mutations are at least partly dominant.

Gallus gallus NEU3 sialidase as model to study protein evolution mechanism based on rapid evolving loops

BACKGROUND

Large surface loops contained within compact protein structures and not involved in catalytic process have been proposed as preferred regions for protein family evolution. These loops are subjected to lower sequence constraints and can evolve rapidly in novel structural variants. A good model to study this hypothesis is represented by sialidase enzymes. Indeed, the structure of sialidases is a β-propeller composed by anti-parallel β-sheets connected by loops that suit well with the rapid evolving loop hypothesis. These features prompted us to extend our studies on this protein family in birds, to get insights on the evolution of this class of glycohydrolases.

RESULTS

Gallus gallus (Gg) genome contains one NEU3 gene encoding a protein with a unique 188 amino acid sequence mainly constituted by a peptide motif repeated six times in tandem with no homology with any other known protein sequence. The repeat region is located at the same position as the roughly 80 amino acid loop characteristic of mammalian NEU4. Based on molecular modeling, all these sequences represent a connecting loop between the first two highly conserved β-strands of the fifth blade of the sialidase β-propeller. Moreover this loop is highly variable in sequence and size in NEU3 sialidases from other vertebrates. Finally, we found that the general enzymatic properties and subcellular localization of Gg NEU3 are not influenced by the deletion of the repeat sequence.

CONCLUSION

In this study we demonstrated that sialidase protein structure contains a surface loop, highly variable both in sequence and size, connecting two conserved β-sheets and emerging on the opposite site of the catalytic crevice. These data confirm that sialidase family can serve as suitable model for the study of the evolutionary process based on rapid evolving loops, which may had occurred in sialidases. Giving the peculiar organization of the loop region identified in Gg NEU3, this protein can be considered of particular interest in such evolutionary studies and to get deeper insights in sialidase evolution.

Spliceosomal intron insertions in genome compacted ray-finned fishes as evident from phylogeny of MC receptors, also supported by a few other GPCRs

BACKGROUND

Insertions of spliceosomal introns are very rare events during evolution of vertebrates and the mechanisms governing creation of novel intron(s) remain obscure. Largely, gene structures of melanocortin (MC) receptors are characterized by intron-less architecture. However, recently a few exceptions have been reported in some fishes. This warrants a systematic survey of MC receptors for understanding intron insertion events during vertebrate evolution.

METHODOLOGY/PRINCIPAL FINDINGS

We have compiled an extended list of MC receptors from different vertebrate genomes with variations in fishes. Notably, the closely linked MC2Rs and MC5Rs from a group of ray-finned fishes have three and one intron insertion(s), respectively, with conserved positions and intron phase. In both genes, one novel insertion was in the highly conserved DRY motif at the end of helix TM3. Further, the proto-splice site MAG↑R is maintained at intron insertion sites in these two genes. However, the orthologs of these receptors from zebrafish and tetrapods are intron-less, suggesting these introns are simultaneously created in selected fishes. Surprisingly, these novel introns are traceable only in four fish genomes. We found that these fish genomes are severely compacted after the separation from zebrafish. Furthermore, we also report novel intron insertions in P2Y receptors and in CHRM3. Finally, we report ultrasmall introns in MC2R genes from selected fishes.

CONCLUSIONS/SIGNIFICANCE

The current repository of MC receptors illustrates that fishes have no MC3R ortholog. MC2R, MC5R, P2Y receptors and CHRM3 have novel intron insertions only in ray-finned fishes that underwent genome compaction. These receptors share one intron at an identical position suggestive of being inserted contemporaneously. In addition to repetitive elements, genome compaction is now believed to be a new hallmark that promotes intron insertions, as it requires rapid DNA breakage and subsequent repair processes to gain back normal functionality.

The genome sequence of the North-European cucumber (Cucumis sativus L.) unravels evolutionary adaptation mechanisms in plants

Cucumber (Cucumis sativus L.), a widely cultivated crop, has originated from Eastern Himalayas and secondary domestication regions includes highly divergent climate conditions e.g. temperate and subtropical. We wanted to uncover adaptive genome differences between the cucumber cultivars and what sort of evolutionary molecular mechanisms regulate genetic adaptation of plants to different ecosystems and organism biodiversity. Here we present the draft genome sequence of the Cucumis sativus genome of the North-European Borszczagowski cultivar (line B10) and comparative genomics studies with the known genomes of: C. sativus (Chinese cultivar – Chinese Long (line 9930)), Arabidopsis thaliana, Populus trichocarpa and Oryza sativa. Cucumber genomes show extensive chromosomal rearrangements, distinct differences in quantity of the particular genes (e.g. involved in photosynthesis, respiration, sugar metabolism, chlorophyll degradation, regulation of gene expression, photooxidative stress tolerance, higher non-optimal temperatures tolerance and ammonium ion assimilation) as well as in distributions of abscisic acid-, dehydration- and ethylene-responsive cis-regulatory elements (CREs) in promoters of orthologous group of genes, which lead to the specific adaptation features. Abscisic acid treatment of non-acclimated Arabidopsis and C. sativus seedlings induced moderate freezing tolerance in Arabidopsis but not in C. sativus. This experiment together with analysis of abscisic acid-specific CRE distributions give a clue why C. sativus is much more susceptible to moderate freezing stresses than A. thaliana. Comparative analysis of all the five genomes showed that, each species and/or cultivars has a specific profile of CRE content in promoters of orthologous genes. Our results constitute the substantial and original resource for the basic and applied research on environmental adaptations of plants, which could facilitate creation of new crops with improved growth and yield in divergent conditions.

Rampant horizontal transfer of SPIN transposons in squamate reptiles

Transposable elements (TEs) are highly abundant in the genome and capable of mobility, two properties that make them particularly prone to transfer horizontally between organisms. Although the impact of horizontal transfer (HT) of TEs is well recognized in prokaryotes, the frequency of this phenomenon and its contribution to genome evolution in eukaryotes remain poorly appreciated. Here, we provide evidence that a DNA transposon called SPIN has colonized the genome of 17 species of reptiles representing nearly every major lineage of squamates, including 14 families of lizards, snakes, and amphisbaenians. Slot blot analyses indicate that SPIN has amplified to high copy numbers in most of these species, ranging from 2,000-28,000 copies per haploid genome. In contrast, we could not detect the presence of SPIN in any of the turtles (seven species from seven families) and crocodiles (four species) examined. Genetic distances between SPIN sequences from species belonging to different squamate families are consistently very low (average = 0.1), considering the deep evolutionary divergence of the families investigated (most are >100 My diverged). Furthermore, these distances fall below interfamilial distances calculated for two genes known to have evolved under strong functional constraint in vertebrates (RAG1, average = 0.24 and C-mos, average = 0.27). These data, combined with phylogenetic analyses, indicate that the widespread distribution of SPIN among squamates is the result of at least 13 independent events of HTs. Molecular dating and paleobiogeographical data suggest that these transfers took place during the last 50 My on at least three different continents (North America, South America and, Africa). Together, these results triple the number of known SPIN transfer events among tetrapods, provide evidence for a previously hypothesized transoceanic movement of SPIN transposons during the Cenozoic, and further underscore the role of HT in the evolution of vertebrate genomes.

Whole genome QTL mapping for growth, meat quality and breast meat yield traits in turkey

Background

The turkey (Meleagris gallopavo) is an important agricultural species and is the second largest contributor to the world's poultry meat production. Demand of turkey meat is increasing very rapidly. Genetic markers linked to genes affecting quantitative traits can increase the selection response of animal breeding programs. The use of these molecular markers for the identification of quantitative trait loci, and subsequently fine-mapping of quantitative trait loci regions, allows for pinpointing of genes that underlie such economically important traits.

Results

The quantitative trait loci analyses of the growth curve, body weight, breast yield and the meat quality traits showed putative quantitative trait loci on 21 of the 27 turkey chromosomes covered by the linkage map. Forty-five quantitative trait loci were detected across all traits and these were found in 29 different regions on 21 chromosomes. Out of the 45 quantitative trait loci, twelve showed significant (p < 0.01) evidence of linkage while the remaining 33 showed suggestive evidence (p < 0.05) of linkage with different growth, growth curve, meat quality and breast yield traits.

Conclusion

A large number of quantitative trait loci were detected across the turkey genome, which affected growth, breast yield and meat quality traits. Pleiotropic effects or close linkages between quantitative trait loci were suggested for several of the chromosomal regions. The comparative analysis regarding the location of quantitative trait loci on different turkey, and on the syntenic chicken chromosomes, along with their phenotypic associations, revealed signs of functional conservation between these species.

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

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

Comparative BAC-based mapping in the white-throated sparrow, a novel behavioral genomics model, using interspecies overgo hybridization

Background

The genomics era has produced an arsenal of resources from sequenced organisms allowing researchers to target species that do not have comparable mapping and sequence information. These new "non-model" organisms offer unique opportunities to examine environmental effects on genomic patterns and processes. Here we use comparative mapping as a first step in characterizing the genome organization of a novel animal model, the white-throated sparrow (Zonotrichia albicollis), which occurs as white or tan morphs that exhibit alternative behaviors and physiology. Morph is determined by the presence or absence of a complex chromosomal rearrangement. This species is an ideal model for behavioral genomics because the association between genotype and phenotype is absolute, making it possible to identify the genomic bases of phenotypic variation.

Findings

We initiated a genomic study in this species by characterizing the white-throated sparrow BAC library via filter hybridization with overgo probes designed for the chicken, turkey, and zebra finch. Cross-species hybridization resulted in 640 positive sparrow BACs assigned to 77 chicken loci across almost all macro-and microchromosomes, with a focus on the chromosomes associated with morph. Out of 216 overgos, 36% of the probes hybridized successfully, with an average number of 3.0 positive sparrow BACs per overgo.

Conclusions

These data will be utilized for determining chromosomal architecture and for fine-scale mapping of candidate genes associated with phenotypic differences. Our research confirms the utility of interspecies hybridization for developing comparative maps in other non-model organisms.

Reference-guided assembly of four diverse Arabidopsis thaliana genomes

We present whole-genome assemblies of four divergent Arabidopsis thaliana strains that complement the 125-Mb reference genome sequence released a decade ago. Using a newly developed reference-guided approach, we assembled large contigs from 9 to 42 Gb of Illumina short-read data from the Landsberg erecta (Ler-1), C24, Bur-0, and Kro-0 strains, which have been sequenced as part of the 1,001 Genomes Project for this species. Using alignments against the reference sequence, we first reduced the complexity of the de novo assembly and later integrated reads without similarity to the reference sequence. As an example, half of the noncentromeric C24 genome was covered by scaffolds that are longer than 260 kb, with a maximum of 2.2 Mb. Moreover, over 96% of the reference genome was covered by the reference-guided assembly, compared with only 87% with a complete de novo assembly. Comparisons with 2 Mb of dideoxy sequence reveal that the per-base error rate of the reference-guided assemblies was below 1 in 10,000. Our assemblies provide a detailed, genomewide picture of large-scale differences between A. thaliana individuals, most of which are difficult to access with alignment-consensus methods only. We demonstrate their practical relevance in studying the expression differences of polymorphic genes and show how the analysis of sRNA sequencing data can lead to erroneous conclusions if aligned against the reference genome alone. Genome assemblies, raw reads, and further information are accessible through http://1001genomes.org/projects/assemblies.html.

The development and characterization of a 60K SNP chip for chicken

BACKGROUND

In livestock species like the chicken, high throughput single nucleotide polymorphism (SNP) genotyping assays are increasingly being used for whole genome association studies and as a tool in breeding (referred to as genomic selection). To be of value in a wide variety of breeds and populations, the success rate of the SNP genotyping assay, the distribution of the SNP across the genome and the minor allele frequencies (MAF) of the SNPs used are extremely important.

RESULTS

We describe the design of a moderate density (60k) Illumina SNP BeadChip in chicken consisting of SNPs known to be segregating at high to medium minor allele frequencies (MAF) in the two major types of commercial chicken (broilers and layers). This was achieved by the identification of 352,303 SNPs with moderate to high MAF in 2 broilers and 2 layer lines using Illumina sequencing on reduced representation libraries. To further increase the utility of the chip, we also identified SNPs on sequences currently not covered by the chicken genome assembly (Gallus_gallus-2.1). This was achieved by 454 sequencing of the chicken genome at a depth of 12x and the identification of SNPs on 454-derived contigs not covered by the current chicken genome assembly. In total we added 790 SNPs that mapped to 454-derived contigs as well as 421 SNPs with a position on Chr_random of the current assembly. The SNP chip contains 57,636 SNPs of which 54,293 could be genotyped and were shown to be segregating in chicken populations. Our SNP identification procedure appeared to be highly reliable and the overall validation rate of the SNPs on the chip was 94%. We were able to map 328 SNPs derived from the 454 sequence contigs on the chicken genome. The majority of these SNPs map to chromosomes that are already represented in genome build Gallus_gallus-2.1.0. Twenty-eight SNPs were used to construct two new linkage groups most likely representing two micro-chromosomes not covered by the current genome assembly.

CONCLUSIONS

The high success rate of the SNPs on the Illumina chicken 60K Beadchip emphasizes the power of Next generation sequence (NGS) technology for the SNP identification and selection step. The identification of SNPs from sequence contigs derived from NGS sequencing resulted in improved coverage of the chicken genome and the construction of two new linkage groups most likely representing two chicken micro-chromosomes.

Oxygen supply from the bird's eye perspective: globin E is a respiratory protein in the chicken retina

The visual process in the vertebrate eye requires high amounts of metabolic energy and thus oxygen. Oxygen supply of the avian retina is a challenging task because birds have large eyes, thick retinae, and high metabolic rates but neither deep retinal nor superficial capillaries. Respiratory proteins such as myoglobin may enhance oxygen supply to certain tissues, and thus the mammalian retina harbors high amounts of neuroglobin. Globin E (GbE) was recently identified as an eye-specific globin of chicken (Gallus gallus). Orthologous GbE genes were found in zebra finch and turkey genomes but appear to be absent in non-avian vertebrate classes. Analyses of globin phylogeny and gene synteny showed an ancient origin of GbE but did not help to assign it to any specific globin type. We show that the photoreceptor cells of the chicken retina have a high level of GbE protein, which accumulates to ∼10 μM in the total eye. Quantitative real-time RT-PCR revealed an ∼50,000-fold higher level of GbE mRNA in the eye than in the brain. Spectroscopic analysis and ligand binding kinetics of recombinant chicken GbE reveal a penta-coordinated globin with an oxygen affinity of P(50) = 5.8 torrs at 25 °C and 15 torrs at 41 °C. Together these data suggest that GbE helps to sustain oxygen supply to the avian retina.

Field guide to next-generation DNA sequencers

The diversity of available 2(nd) and 3(rd) generation DNA sequencing platforms is increasing rapidly. Costs for these systems range from < $100,000 to more than $1,000,000, with instrument run times ranging from minutes to weeks. Extensive trade-offs exist among these platforms. I summarize the major characteristics of each commercially available platform to enable direct comparisons. In terms of cost per megabase (Mb) of sequence, the Illumina and SOLiD platforms are clearly superior (≤ $0.10/Mb vs. > $10/Mb for 454 and some Ion Torrent chips). In terms of cost per nonmultiplexed sample and instrument run time, the Pacific Biosciences and Ion Torrent platforms excel, with the 454 GS Junior and Illumina MiSeq also notable in this regard. All platforms allow multiplexing of samples, but details of library preparation, experimental design and data analysis can constrain the options. The wide range of characteristics among available platforms provides opportunities both to conduct groundbreaking studies and to waste money on scales that were previously infeasible. Thus, careful thought about the desired characteristics of these systems is warranted before purchasing or using any of them. Updated information from this guide will be maintained at: http://dna.uga.edu/ and http://tomato.biol.trinity.edu/blog/.

Sim4db and Leaff: utilities for fast batch spliced alignment and sequence indexing

The large number of genomes that will be sequenced will need to be annotated with genes and other functional features. Aligning gene sequences from a related species to the target genome is an economical and highly reliable method to identify genes; unfortunately, existing tools have been lacking in sensitivity and speed. A program we reported, sim4cc, was shown to be highly accurate but is limited to comparing one cDNA with one genomic sequence. We present here an optimization of the tool, implemented in the packages sim4db and leaff. The new tool performs batch alignments of cDNA and genomic sequences in a fraction of the time required by its predecessor, and thus is very well suited for genome-wide analyses.

AVAILABILITY

Sim4db and leaff are written in C, C++ and Perl for Linux and other Unix platforms. Source code is distributed free of charge from http://sourceforge.net/projects/kmer/.

CONTACT

florea@umiacs.umd.edu

A macroevolutionary perspective on multiple sexual traits in the phasianidae (galliformes)

Traits involved in sexual signaling are ubiquitous among animals. Although a single trait appears sufficient to convey information, many sexually dimorphic species exhibit multiple sexual signals, which may be costly to signalers and receivers. Given that one signal may be enough, there are many microevolutionary hypotheses to explain the evolution of multiple signals. Here we extend these hypotheses to a macroevolutionary scale and compare those predictions to the patterns of gains and losses of sexual dimorphism in pheasants and partridges. Among nine dimorphic characters, including six intersexual signals and three indicators of competitive ability, all exhibited both gains and losses of dimorphism within the group. Although theories of intersexual selection emphasize gain and elaboration, those six characters exhibited greater rates of loss than gain; in contrast, the competitive traits showed a slight bias towards gains. The available models, when examined in a macroevolutionary framework, did not yield unique predictions, making it difficult to distinguish among them. Even with this limitation, when the predictions of these alternative models were compared with the heterogeneous patterns of evolution of dimorphism in phasianids, it is clear that many different selective processes have been involved in the evolution of sexual signals in this group.

Genomic insight into the common carp (Cyprinus carpio) genome by sequencing analysis of BAC-end sequences

Background

Common carp is one of the most important aquaculture teleost fish in the world. Common carp and other closely related Cyprinidae species provide over 30% aquaculture production in the world. However, common carp genomic resources are still relatively underdeveloped. BAC end sequences (BES) are important resources for genome research on BAC-anchored genetic marker development, linkage map and physical map integration, and whole genome sequence assembling and scaffolding.

Result

To develop such valuable resources in common carp (Cyprinus carpio), a total of 40,224 BAC clones were sequenced on both ends, generating 65,720 clean BES with an average read length of 647 bp after sequence processing, representing 42,522,168 bp or 2.5% of common carp genome. The first survey of common carp genome was conducted with various bioinformatics tools. The common carp genome contains over 17.3% of repetitive elements with GC content of 36.8% and 518 transposon ORFs. To identify and develop BAC-anchored microsatellite markers, a total of 13,581 microsatellites were detected from 10,355 BES. The coding region of 7,127 genes were recognized from 9,443 BES on 7,453 BACs, with 1,990 BACs have genes on both ends. To evaluate the similarity to the genome of closely related zebrafish, BES of common carp were aligned against zebrafish genome. A total of 39,335 BES of common carp have conserved homologs on zebrafish genome which demonstrated the high similarity between zebrafish and common carp genomes, indicating the feasibility of comparative mapping between zebrafish and common carp once we have physical map of common carp.

Conclusion

BAC end sequences are great resources for the first genome wide survey of common carp. The repetitive DNA was estimated to be approximate 28% of common carp genome, indicating the higher complexity of the genome. Comparative analysis had mapped around 40,000 BES to zebrafish genome and established over 3,100 microsyntenies, covering over 50% of the zebrafish genome. BES of common carp are tremendous tools for comparative mapping between the two closely related species, zebrafish and common carp, which should facilitate both structural and functional genome analysis in common carp.

Origin of amphibian and avian chromosomes by fission, fusion, and retention of ancestral chromosomes

Amphibian genomes differ greatly in DNA content and chromosome size, morphology, and number. Investigations of this diversity are needed to identify mechanisms that have shaped the evolution of vertebrate genomes. We used comparative mapping to investigate the organization of genes in the Mexican axolotl (Ambystoma mexicanum), a species that presents relatively few chromosomes (n = 14) and a gigantic genome (>20 pg/N). We show extensive conservation of synteny between Ambystoma, chicken, and human, and a positive correlation between the length of conserved segments and genome size. Ambystoma segments are estimated to be four to 51 times longer than homologous human and chicken segments. Strikingly, genes demarking the structures of 28 chicken chromosomes are ordered among linkage groups defining the Ambystoma genome, and we show that these same chromosomal segments are also conserved in a distantly related anuran amphibian (Xenopus tropicalis). Using linkage relationships from the amphibian maps, we predict that three chicken chromosomes originated by fusion, nine to 14 originated by fission, and 12-17 evolved directly from ancestral tetrapod chromosomes. We further show that some ancestral segments were fused prior to the divergence of salamanders and anurans, while others fused independently and randomly as chromosome numbers were reduced in lineages leading to Ambystoma and Xenopus. The maintenance of gene order relationships between chromosomal segments that have greatly expanded and contracted in salamander and chicken genomes, respectively, suggests selection to maintain synteny relationships and/or extremely low rates of chromosomal rearrangement. Overall, the results demonstrate the value of data from diverse, amphibian genomes in studies of vertebrate genome evolution.

Selective loss of glycogen synthase kinase-3α in birds reveals distinct roles for GSK-3 isozymes in tau phosphorylation

Mammalian glycogen synthase kinase-3 (GSK-3), a critical regulator in neuronal signaling, cognition, and behavior, exists as two isozymes GSK-3α and GSK-3β. Their distinct biological functions remains largely unknown. Here, we examined the evolutionary significance of each of these isozymes. Surprisingly, we found that unlike other vertebrates that harbor both GSK-3 genes, the GSK-3α gene is missing in birds. GSK-3-mediated tau phosphorylation was significantly lower in adult bird brains than in mouse brains, a phenomenon that was reproduced in GSK-3α knockout mouse brains. Tau phosphorylation was detected in brains from bird embryos suggesting that GSK-3 isozymes play distinct roles in tau phosphorylation during development. Birds are natural GSK-3α knockout organisms and may serve as a novel model to study the distinct functions of GSK-3 isozymes.

A fast, lock-free approach for efficient parallel counting of occurrences of k-mers

MOTIVATION

Counting the number of occurrences of every k-mer (substring of length k) in a long string is a central subproblem in many applications, including genome assembly, error correction of sequencing reads, fast multiple sequence alignment and repeat detection. Recently, the deep sequence coverage generated by next-generation sequencing technologies has caused the amount of sequence to be processed during a genome project to grow rapidly, and has rendered current k-mer counting tools too slow and memory intensive. At the same time, large multicore computers have become commonplace in research facilities allowing for a new parallel computational paradigm.

RESULTS

We propose a new k-mer counting algorithm and associated implementation, called Jellyfish, which is fast and memory efficient. It is based on a multithreaded, lock-free hash table optimized for counting k-mers up to 31 bases in length. Due to their flexibility, suffix arrays have been the data structure of choice for solving many string problems. For the task of k-mer counting, important in many biological applications, Jellyfish offers a much faster and more memory-efficient solution.

AVAILABILITY

The Jellyfish software is written in C++ and is GPL licensed. It is available for download at http://www.cbcb.umd.edu/software/jellyfish.

Next generation sequence assembly with AMOS

A Modular Open-Source Assembler (AMOS) was designed to offer a modular approach to genome assembly. AMOS includes a wide range of tools for assembly, including the lightweight de novo assemblers Minimus and Minimo, and Bambus 2, a robust scaffolder able to handle metagenomic and polymorphic data. This protocol describes how to configure and use AMOS for the assembly of Next Generation sequence data. Additionally, we provide three tutorial examples that include bacterial, viral, and metagenomic datasets with specific tips for improving assembly quality.

Poultry genome sequences: progress and outstanding challenges

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.

A fast, lock-free approach for efficient parallel counting of occurrences of k-mers

MOTIVATION

Counting the number of occurrences of every k-mer (substring of length k) in a long string is a central subproblem in many applications, including genome assembly, error correction of sequencing reads, fast multiple sequence alignment and repeat detection. Recently, the deep sequence coverage generated by next-generation sequencing technologies has caused the amount of sequence to be processed during a genome project to grow rapidly, and has rendered current k-mer counting tools too slow and memory intensive. At the same time, large multicore computers have become commonplace in research facilities allowing for a new parallel computational paradigm.

RESULTS

We propose a new k-mer counting algorithm and associated implementation, called Jellyfish, which is fast and memory efficient. It is based on a multithreaded, lock-free hash table optimized for counting k-mers up to 31 bases in length. Due to their flexibility, suffix arrays have been the data structure of choice for solving many string problems. For the task of k-mer counting, important in many biological applications, Jellyfish offers a much faster and more memory-efficient solution.

AVAILABILITY

The Jellyfish software is written in C++ and is GPL licensed. It is available for download at http://www.cbcb.umd.edu/software/jellyfish.

Ensembl 2011

The Ensembl project (http://www.ensembl.org) seeks to enable genomic science by providing high quality, integrated annotation on chordate and selected eukaryotic genomes within a consistent and accessible infrastructure. All supported species include comprehensive, evidence-based gene annotations and a selected set of genomes includes additional data focused on variation, comparative, evolutionary, functional and regulatory annotation. The most advanced resources are provided for key species including human, mouse, rat and zebrafish reflecting the popularity and importance of these species in biomedical research. As of Ensembl release 59 (August 2010), 56 species are supported of which 5 have been added in the past year. Since our previous report, we have substantially improved the presentation and integration of both data of disease relevance and the regulatory state of different cell types.

Biomarkers as management tools for industries in the pork production chain

In livestock production, traits such as meat quantity and quality are mainly determined by biological processes. Knowing the biological background of traits makes it possible to identify relevant processes, genes or molecules that may be used as indicators to rapidly and easily detect the status and phase of biological processes. Such indicators are called biomarkers. Thus biomarkers provide information about the status and phase of biological processes underlying particular traits. This information can be used to modulate these processes and thereby optimize the economic value of these traits. To produce products from livestock, pork industries are organized in a production chain – from the breeding industry to the retail industry. Industries in a chain interact with each other to optimize economic value. (Automated) detection of biomarkers could activate processes throughout the chain that can be robotized and computerized, thereby creating economic benefit for all participants in the chain. Biomarkers may be used on individual animals and on pooled samples representing a whole herd for herd management. Here we describe how biomarkers can be developed and how they may be used to enhance value for the pork production chain industries.

Paleogenomic data suggest mammal-like genome size in the ancestral amniote and derived large genome size in amphibians

An unsolved question in evolutionary genomics is whether amniote genomes have been expanding or contracting since the common ancestor of this diverse group. Here, we report on the polarity of amniote genome size evolution using genome size estimates for 14 extinct tetrapod genera from the Paleozoic and early Mesozoic Eras using osteocyte lacunae size as a correlate. We find substantial support for a phylogenetically controlled regression model relating genome size to osteocyte lacunae size (P of slopes <0.01, r²=0.65, phylogenetic signal λ=0.83). Genome size appears to have been homogeneous across Paleozoic crown-tetrapod lineages (average haploid genome size 2.9-3.7 pg) with values similar to those of extant mammals. The differentiation in genome size and underlying architecture among extant tetrapod lineages likely evolved in the Mesozoic and Cenozoic Eras, with expansion in amphibians, contractions along the diapsid lineage, and no directional change within the synapsid lineage leading to mammals.

A SNP based linkage map of the turkey genome reveals multiple intrachromosomal rearrangements between the turkey and chicken genomes

BACKGROUND

The turkey (Meleagris gallopavo) is an important agricultural species that is the second largest contributor to the world's poultry meat production. The genomic resources of turkey provide turkey breeders with tools needed for the genetic improvement of commercial breeds of turkey for economically important traits. A linkage map of turkey is essential not only for the mapping of quantitative trait loci, but also as a framework to enable the assignment of sequence contigs to specific chromosomes. Comparative genomics with chicken provides insight into mechanisms of genome evolution and helps in identifying rare genomic events such as genomic rearrangements and duplications/deletions.

RESULTS

Eighteen full sib families, comprising 1008 (35 F1 and 973 F2) birds, were genotyped for 775 single nucleotide polymorphisms (SNPs). Of the 775 SNPs, 570 were informative and used to construct a linkage map in turkey. The final map contains 531 markers in 28 linkage groups. The total genetic distance covered by these linkage groups is 2,324 centimorgans (cM) with the largest linkage group (81 loci) measuring 326 cM. Average marker interval for all markers across the 28 linkage groups is 4.6 cM. Comparative mapping of turkey and chicken revealed two inter-, and 57 intrachromosomal rearrangements between these two species.

CONCLUSION

Our turkey genetic map of 531 markers reveals a genome length of 2,324 cM. Our linkage map provides an improvement of previously published maps because of the more even distribution of the markers and because the map is completely based on SNP markers enabling easier and faster genotyping assays than the microsatellitemarkers used in previous linkage maps. Turkey and chicken are shown to have a highly conserved genomic structure with a relatively low number of inter-, and intrachromosomal rearrangements.