Insights into the evolution of Darwin's finches from comparative analysis of the Geospiza magnirostris genome sequence

Background

A classical example of repeated speciation coupled with ecological diversification is the evolution of 14 closely related species of Darwin’s (Galápagos) finches (Thraupidae, Passeriformes). Their adaptive radiation in the Galápagos archipelago took place in the last 2–3 million years and some of the molecular mechanisms that led to their diversification are now being elucidated. Here we report evolutionary analyses of genome of the large ground finch, Geospiza magnirostris.

Results

13,291 protein-coding genes were predicted from a 991.0 Mb G. magnirostris genome assembly. We then defined gene orthology relationships and constructed whole genome alignments between the G. magnirostris and other vertebrate genomes. We estimate that 15% of genomic sequence is functionally constrained between G. magnirostris and zebra finch. Genic evolutionary rate comparisons indicate that similar selective pressures acted along the G. magnirostris and zebra finch lineages suggesting that historical effective population size values have been similar in both lineages. 21 otherwise highly conserved genes were identified that each show evidence for positive selection on amino acid changes in the Darwin's finch lineage. Two of these genes (Igf2r and Pou1f1) have been implicated in beak morphology changes in Darwin’s finches. Five of 47 genes showing evidence of positive selection in early passerine evolution have cilia related functions, and may be examples of adaptively evolving reproductive proteins.

Conclusions

These results provide insights into past evolutionary processes that have shaped G. magnirostris genes and its genome, and provide the necessary foundation upon which to build population genomics resources that will shed light on more contemporaneous adaptive and non-adaptive processes that have contributed to the evolution of the Darwin’s finches.

A draft sequence of the Neandertal genome

Neandertals, the closest evolutionary relatives of present-day humans, lived in large parts of Europe and western Asia before disappearing 30,000 years ago. We present a draft sequence of the Neandertal genome composed of more than 4 billion nucleotides from three individuals. Comparisons of the Neandertal genome to the genomes of five present-day humans from different parts of the world identify a number of genomic regions that may have been affected by positive selection in ancestral modern humans, including genes involved in metabolism and in cognitive and skeletal development. We show that Neandertals shared more genetic variants with present-day humans in Eurasia than with present-day humans in sub-Saharan Africa, suggesting that gene flow from Neandertals into the ancestors of non-Africans occurred before the divergence of Eurasian groups from each other.

Mutation discovery in the mouse using genetically guided array capture and resequencing

Forward genetics (phenotype-driven approaches) remain the primary source for allelic variants in the mouse. Unfortunately, the gap between observable phenotype and causative genotype limits the widespread use of spontaneous and induced mouse mutants. As alternatives to traditional positional cloning and mutation detection approaches, sequence capture and next-generation sequencing technologies can be used to rapidly sequence subsets of the genome. Application of these technologies to mutation detection efforts in the mouse has the potential to significantly reduce the time and resources required for mutation identification by abrogating the need for high-resolution genetic mapping, long-range PCR, and sequencing of individual PCR amplimers. As proof of principle, we used array-based sequence capture and pyrosequencing to sequence an allelic series from the classically defined Kit locus (approximately 200 kb) from each of five noncomplementing Kit mutants (one known allele and four unknown alleles) and have successfully identified and validated a nonsynonymous coding mutation for each allele. These data represent the first documentation and validation that these new technologies can be used to efficiently discover causative mutations. Importantly, these data also provide a specific methodological foundation for the development of large-scale mutation detection efforts in the laboratory mouse.

Analysis of copy number variants and segmental duplications in the human genome: Evidence for a change in the process of formation in recent evolutionary history

Segmental duplications (SDs) are operationally defined as >1 kb stretches of duplicated DNA with high sequence identity. They arise from copy number variants (CNVs) fixed in the population. To investigate the formation of SDs and CNVs, we examine their large-scale patterns of co-occurrence with different repeats. Alu elements, a major class of genomic repeats, had previously been identified as prime drivers of SD formation. We also observe this association; however, we find that it sharply decreases for younger SDs. Continuing this trend, we find only weak associations of CNVs with Alus. Similarly, we find an association of SDs with processed pseudogenes, which is decreasing for younger SDs and absent entirely for CNVs. Next, we find that SDs are significantly co-localized with each other, resulting in a highly skewed "power-law" distribution and chromosomal hotspots. We also observe a significant association of CNVs with SDs, but find that an SD-mediated mechanism only accounts for some CNVs (<28%). Overall, our results imply that a shift in predominant formation mechanism occurred in recent history: approximately 40 million years ago, during the "Alu burst" in retrotransposition activity, non-allelic homologous recombination, first mediated by Alus and then the by newly formed CNVs themselves, was the main driver of genome rearrangements; however, its relative importance has decreased markedly since then, with proportionally more events now stemming from other repeats and from non-homologous end-joining. In addition to a coarse-grained analysis, we performed targeted sequencing of 67 CNVs and then analyzed a combined set of 270 CNVs (540 breakpoints) to verify our conclusions.

The nonhomologous end joining factor Artemis suppresses multi-tissue tumor formation and prevents loss of heterozygosity

Nonhomologous end joining (NHEJ) is a critical DNA repair pathway, with proposed tumor suppression functions in many tissues. Mutations in the NHEJ factor ARTEMIS cause radiation-sensitive severe combined immunodeficiency in humans and may increase susceptibility to lymphoma in some settings. We now report that deficiency for Artemis (encoded by Dclre1c/Art in mouse) accelerates tumorigenesis in several tissues in a Trp53 heterozygous setting, revealing tumor suppression roles for NHEJ in lymphoid and non-lymphoid cells. We also show that B-lineage lymphomas in these mice undergo loss of Trp53 heterozygosity by allele replacement, but arise by mechanisms distinct from those in Art Trp53 double null mice. These findings demonstrate a general tumor suppression function for NHEJ, and reveal that interplay between NHEJ and Trp53 loss of heterozygosity influences the sequence of multi-hit oncogenesis. We present a model where p53 status at the time of tumor initiation is a key determinant of subsequent oncogenic mechanisms. Because Art deficient mice represent a model for radiation-sensitive severe combined immunodeficiency, our findings suggest that these patients may be at risk for both lymphoid and non-lymphoid cancers.

Selective degradation of transcripts during meiotic maturation of mouse oocytes

There is massive destruction of transcripts during the maturation of mouse oocytes. The objective of this project was to identify and characterize the transcripts that are degraded versus those that are stable during the transcriptionally silent germinal vesicle (GV)-stage to metaphase II (MII)-stage transition using a microarray approach. A system for oocyte transcript amplification using both internal and 3'-poly(A) priming was utilized to minimize the impact of complex variations in transcript polyadenylation prevalent during this transition. Transcripts were identified and quantified using the Affymetrix Mouse Genome 430 v2.0 GeneChip. The significantly changed and stable transcripts were analyzed using Ingenuity Pathways Analysis and GenMAPP/MAPPFinder to characterize the biological themes underlying global changes in oocyte transcripts during maturation. It was concluded that the destruction of transcripts during the GV to MII transition is a selective rather than promiscuous process in mouse oocytes. In general, transcripts involved in processes that are associated with meiotic arrest at the GV-stage and the progression of oocyte maturation, such as oxidative phosphorylation, energy production, and protein synthesis and metabolism, were dramatically degraded. In contrast, transcripts encoding participants in signaling pathways essential for maintaining the unique characteristics of the MII-arrested oocyte, such as those involved in protein kinase pathways, were the most prominent among the stable transcripts.

Inheritance patterns of transcript levels in F1 hybrid mice

Genetic analysis of transcriptional regulation is a rapidly emerging field of investigation that promises to shed light on the regulatory networks that control gene expression. Although a number of such studies have been carried out, the nature and extent of the heritability of gene expression traits have not been well established. We describe the inheritance of transcript levels in liver tissue in the first filial (F1) generation of mice obtained from reciprocal crosses between the common inbred strains A/J and C57BL/6J. We obtain estimates of genetic and technical variance components from these data and demonstrate that shrinkage estimators can increase detectable heritability. Estimates of heritability vary widely from transcript to transcript, with one-third of transcripts showing essentially no heritability (<0.01) and one-quarter showing very high heritability (>0.50). Roughly half of all transcripts are differentially expressed between the two parental strains. Most transcripts show an additive pattern of inheritance. Dominance effects were observed for 20% of transcripts and a small number of transcripts were identified as showing an overdominance mode of inheritance. In addition, we identified 314 transcripts with expression levels that differ between the reciprocal F1 animals. These genes may be related to maternal effect.

A comparison of cDNA, oligonucleotide, and Affymetrix GeneChip gene expression microarray platforms

We have conducted a study to compare the variability in measured gene expression levels associated with three types of microarray platforms. Total RNA samples were obtained from liver tissue of four male mice, two each from inbred strains A/J and C57BL/6J. The same four samples were assayed on Affymetrix Mouse Genome Expression Set 430 GeneChips (MOE430A and MOE430B), spotted cDNA microarrays, and spotted oligonucleotide microarrays using eight arrays of each type. Variances associated with measurement error were observed to be comparable across all microarray platforms. The MOE430A GeneChips and cDNA arrays had higher precision across technical replicates than the MOE430B GeneChips and oligonucleotide arrays. The Affymetrix platform showed the greatest range in the magnitude of expression levels followed by the oligonucleotide arrays. We observed good concordance in both estimated expression level and statistical significance of common genes between the Affymetrix MOE430A GeneChip and the oligonucleotide arrays. Despite their apparently high precision, cDNA arrays showed poor concordance with other platforms.

Distribution of Nitrogen-Fixing Microorganisms along the Neuse River Estuary, North Carolina

Nitrogen fixation genes (nifH) were amplified and sequenced from DNA extracted from surface water samples collected from six stations along the length of the Neuse River Estuary, North Carolina, in order to determine the distribution of nitrogen-fixing organisms in the transition from fresh- to saltwater. Nitrogenase genes were detected in all samples by a nested polymerase chain reaction method, and the amplification products from the upriver, midriver, and downriver stations were cloned, sequenced, and used for phylogenetic analysis. The composition of nifH clone libraries from upriver, midriver, and downriver stations (each composed of 14 randomly selected clones) were very diverse (samples from upriver and midriver stations were composed of 14 unique sequences, downriver station composed of 7 unique sequences) and differed among the stations. Some phylotypes were found at more than one station, but were usually found in the upriver and midriver stations or in the midriver and downriver stations, indicating that the phylotypes were probably transported along the river. Cyanobacterial nifH were not found at the most upriver site, but were a large fraction of sequences (50%) recovered from the downriver station, where nitrate concentration was an order of magnitude lower and salinity was higher. In contrast, g proteobacteria nifH sequences were much more common at the midriver and upriver sites (58% and 64%, respectively), compared to the downriver site (14%). Results indicate that substantially different nitrogen-fixing assemblages are present along the river, reflecting differential watershed hydrological inputs, sedimentation, and environmental selection pressures, along the salinity gradient.