Self-fertilization sweeps up variation in the worm genome

A new study reports a comprehensive survey of genetic diversity in natural populations of the nematode Caenorhabditis elegans. Their analyses suggest that recent chromosome-scale selective sweeps have reduced C. elegans genetic diversity worldwide and strongly structured genetic variation across its genome.     

A genome-wide map of diversity in Plasmodium falciparum

Genetic variation allows the malaria parasite Plasmodium falciparum to overcome chemotherapeutic agents, vaccines and vector control strategies and remain a leading cause of global morbidity and mortality. Here we describe an initial survey of genetic variation across the P. falciparum genome. We performed extensive sequencing of 16 geographically diverse parasites and identified 46,937 SNPs, demonstrating rich diversity among P. falciparum parasites (pi = 1.16 x 10(-3)) and strong correlation with gene function. We identified multiple regions with signatures of selective sweeps in drug-resistant parasites, including a previously unidentified 160-kb region with extremely low polymorphism in pyrimethamine-resistant parasites. We further characterized 54 worldwide isolates by genotyping SNPs across 20 genomic regions. These data begin to define population structure among African, Asian and American groups and illustrate the degree of linkage disequilibrium, which extends over relatively short distances in African parasites but over longer distances in Asian parasites. We provide an initial map of genetic diversity in P. falciparum and demonstrate its potential utility in identifying genes subject to recent natural selection and in understanding the population genetics of this parasite.     

Whole-genome sequencing of multiple Arabidopsis thaliana populations

The plant Arabidopsis thaliana occurs naturally in many different habitats throughout Eurasia. As a foundation for identifying genetic variation contributing to adaptation to diverse environments, a 1001 Genomes Project to sequence geographically diverse A. thaliana strains has been initiated. Here we present the first phase of this project, based on population-scale sequencing of 80 strains drawn from eight regions throughout the species' native range. We describe the majority of common small-scale polymorphisms as well as many larger insertions and deletions in the A. thaliana pan-genome, their effects on gene function, and the patterns of local and global linkage among these variants. The action of processes other than spontaneous mutation is identified by comparing the spectrum of mutations that have accumulated since A. thaliana diverged from its closest relative 10 million years ago with the spectrum observed in the laboratory. Recent species-wide selective sweeps are rare, and potentially deleterious mutations are more common in marginal populations.     

Conservation and diversification of Wnt signaling function during the evolution of nematode vulva development

Cell-fate specification and cell-cell signaling have been well studied during vulva development in Caenorhabditis elegans and provide a paradigm in evolutionary developmental biology. Pristionchus pacificus has been developed as a 'satellite' organism with an integrated physical and genetic map that allows detailed comparisons to C. elegans. A common aspect of vulva formation in both species is the polarization of the P7.p lineage, which is responsible for vulval symmetry. In C. elegans, Wnt signaling is crucial for P7.p cell-fate patterning; nothing is known about vulval symmetry in P. pacificus. We isolated mutations that disrupt polarization of the P7.p lineage in P. pacificus and found that the corresponding gene encodes a Frizzled-like molecule. In addition, mutations in Ppa-lin-17 (encoding Frizzled) and morpholino knock-down of Ppa-lin-44 (encoding Wnt), Ppa-egl-20 (encoding Wnt), Ppa-mig-5 (encoding Dsh), Ppa-apr-1 (encoding APC) and Ppa-bar-1 (encoding beta-catenin) results in gonad-independent vulva differentiation, indicating that these genes have a role in a negative signaling process. In contrast, in C. elegans, Wnt signaling has a positive role in vulva induction, and mutations in bar-1 result in a hypoinduced phenotype. Therefore, whereas the molecular mechanisms that generate vulval symmetry are conserved, the genetic control of vulva induction diversified during evolution.     

Structural diversity and African origin of the 17q21.31 inversion polymorphism

The 17q21.31 inversion polymorphism exists either as direct (H1) or inverted (H2) haplotypes with differential predispositions to disease and selection. We investigated its genetic diversity in 2,700 individuals, with an emphasis on African populations. We characterize eight structural haplotypes due to complex rearrangements that vary in size from 1.08-1.49 Mb and provide evidence for a 30-kb H1-H2 double recombination event. We show that recurrent partial duplications of the KANSL1 gene have occurred on both the H1 and H2 haplotypes and have risen to high frequency in European populations. We identify a likely ancestral H2 haplotype (H2') lacking these duplications that is enriched among African hunter-gatherer groups yet essentially absent from West African populations. Whereas H1 and H2 segmental duplications arose independently and before human migration out of Africa, they have reached high frequencies recently among Europeans, either because of extraordinary genetic drift or selective sweeps.     

Systematic mapping of genetic interactions in Caenorhabditis elegans identifies common modifiers of diverse signaling pathways

Most heritable traits, including disease susceptibility, are affected by interactions between multiple genes. However, we understand little about how genes interact because very few possible genetic interactions have been explored experimentally. We have used RNA interference in Caenorhabditis elegans to systematically test approximately 65,000 pairs of genes for their ability to interact genetically. We identify approximately 350 genetic interactions between genes functioning in signaling pathways that are mutated in human diseases, including components of the EGF/Ras, Notch and Wnt pathways. Most notably, we identify a class of highly connected 'hub' genes: inactivation of these genes can enhance the phenotypic consequences of mutation of many different genes. These hub genes all encode chromatin regulators, and their activity as genetic hubs seems to be conserved across animals. We propose that these genes function as general buffers of genetic variation and that these hub genes may act as modifier genes in multiple, mechanistically unrelated genetic diseases in humans.     

A model-based approach for analysis of spatial structure in genetic data

Characterizing genetic diversity within and between populations has broad applications in studies of human disease and evolution. We propose a new approach, spatial ancestry analysis, for the modeling of genotypes in two- or three-dimensional space. In spatial ancestry analysis (SPA), we explicitly model the spatial distribution of each SNP by assigning an allele frequency as a continuous function in geographic space. We show that the explicit modeling of the allele frequency allows individuals to be localized on the map on the basis of their genetic information alone. We apply our SPA method to a European and a worldwide population genetic variation data set and identify SNPs showing large gradients in allele frequency, and we suggest these as candidate regions under selection. These regions include SNPs in the well-characterized LCT region, as well as at loci including FOXP2, OCA2 and LRP1B.     

Open-reading-frame sequence tags (OSTs) support the existence of at least 17,300 genes in C. elegans

The genome sequences of Caenorhabditis elegans, Drosophila melanogaster and Arabidopsis thaliana have been predicted to contain 19,000, 13,600 and 25,500 genes, respectively. Before this information can be fully used for evolutionary and functional studies, several issues need to be addressed. First, the gene number estimates obtained in silico and not yet supported by any experimental data need to be verified. For example, it seems biologically paradoxical that C. elegans would have 50% more genes than Drosophilia. Second, intron/exon predictions need to be tested experimentally. Third, complete sets of open reading frames (ORFs), or "ORFeomes," need to be cloned into various expression vectors. To address these issues simultaneously, we have designed and applied to C. elegans the following strategy. Predicted ORFs are amplified by PCR from a highly representative cDNA library using ORF-specific primers, cloned by Gateway recombination cloning and then sequenced to generate ORF sequence tags (OSTs) as a way to verify identity and splicing. In a sample (n=1,222) of the nearly 10,000 genes predicted ab initio (that is, for which no expressed sequence tag (EST) is available so far), at least 70% were verified by OSTs. We also observed that 27% of these experimentally confirmed genes have a structure different from that predicted by GeneFinder. We now have experimental evidence that supports the existence of at least 17,300 genes in C. elegans. Hence we suggest that gene counts based primarily on ESTs may underestimate the number of genes in human and in other organisms.     

A systematic RNAi screen identifies a critical role for mitochondria in C. elegans longevity

We report a systematic RNA interference (RNAi) screen of 5,690 Caenorhabditis elegans genes for gene inactivations that increase lifespan. We found that genes important for mitochondrial function stand out as a principal group of genes affecting C. elegans lifespan. A classical genetic screen identified a mutation in the mitochondrial leucyl-tRNA synthetase gene (lrs-2) that impaired mitochondrial function and was associated with longer-lifespan. The long-lived worms with impaired mitochondria had lower ATP content and oxygen consumption, but differential responses to free-radical and other stresses. These data suggest that the longer lifespan of C. elegans with compromised mitochrondria cannot simply be assigned to lower free radical production and suggest a more complex coupling of metabolism and longevity.     

Polygenic susceptibility to breast cancer and implications for prevention

The knowledge of human genetic variation that will come from the human genome sequence makes feasible a polygenic approach to disease prevention, in which it will be possible to identify individuals as susceptible by their genotype profile and to prevent disease by targeting interventions to those at risk. There is doubt, however, regarding the magnitude of these genetic effects and thus the potential to apply them to either individuals or populations. We have therefore examined the potential for prediction of risk based on common genetic variation using data from a population-based series of individuals with breast cancer. The data are compatible with a log-normal distribution of genetic risk in the population that is sufficiently wide to provide useful discrimination of high- and low-risk groups. Assuming all of the susceptibility genes could be identified, the half of the population at highest risk would account for 88% of all affected individuals. By contrast, if currently identified risk factors for breast cancer were used to stratify the population, the half of the population at highest risk would account for only 62% of all cases. These results suggest that the construction and use of genetic-risk profiles may provide significant improvements in the efficacy of population-based programs of intervention for cancers and other diseases.     

Measurement of the human allele frequency spectrum demonstrates greater genetic drift in East Asians than in Europeans

Large data sets on human genetic variation have been collected recently, but their usefulness for learning about history and natural selection has been limited by biases in the ways polymorphisms were chosen. We report large subsets of SNPs from the International HapMap Project that allow us to overcome these biases and to provide accurate measurement of a quantity of crucial importance for understanding genetic variation: the allele frequency spectrum. Our analysis shows that East Asian and northern European ancestors shared the same population bottleneck expanding out of Africa but that both also experienced more recent genetic drift, which was greater in East Asians.     

Caenorhabditis elegans expressed sequence tags identify gene families and potential disease gene homologues.

A database containing mapped partial cDNA sequences from Caenorhabditis elegans will provide a ready starting point for identifying nematode homologues of important human genes and determining their functions in C. elegans. A total of 720 expressed sequence tags (ESTs) have been generated from 585 clones randomly selected from a mixed-stage C. elegans cDNA library. Comparison of these ESTs with sequence databases identified 422 new C. elegans genes, of which 317 are not similar to any sequences in the database. Twenty-six new genes have been mapped by YAC clone hybridization. Members of several gene families, including cuticle collagens, GTP-binding proteins, and RNA helicases were discovered. Many of the new genes are similar to known or potential human disease genes, including CFTR and the LDL receptor.     

Evidence for substantial fine-scale variation in recombination rates across the human genome

Characterizing fine-scale variation in human recombination rates is important, both to deepen understanding of the recombination process and to aid the design of disease association studies. Current genetic maps show that rates vary on a megabase scale, but studying finer-scale variation using pedigrees is difficult. Sperm-typing experiments have characterized regions where crossovers cluster into 1-2-kb hot spots, but technical difficulties limit the number of studies. An alternative is to use population variation to infer fine-scale characteristics of the recombination process. Several surveys reported 'block-like' patterns of diversity, which may reflect fine-scale recombination rate variation, but limitations of available methods made this impossible to assess. Here, we applied a new statistical method, which overcomes these limitations, to infer patterns of fine-scale recombination rate variation in 74 genes. We found extensive rate variation both within and among genes. In particular, recombination hot spots are a common feature of the human genome: 47% (35 of 74) of genes showed substantive evidence for a hot spot, and many more showed evidence for some rate variation. No primary sequence characteristics are consistently associated with precise hot-spot location, although G+C content and nucleotide diversity are correlated with local recombination rate.     

A stress-sensitive reporter predicts longevity in isogenic populations of Caenorhabditis elegans

When both genotype and environment are held constant, 'chance' variation in the lifespan of individuals in a population is still quite large. Using isogenic populations of the nematode Caenorhabditis elegans, we show that, on the first day of adult life, chance variation in the level of induction of a green fluorescent protein (GFP) reporter coupled to a promoter from the gene hsp-16.2 predicts as much as a fourfold variation in subsequent survival. The same reporter is also a predictor of ability to withstand a subsequent lethal thermal stress. The level of induction of GFP is not heritable, and GFP expression levels in other reporter constructs are not associated with differences in longevity. HSP-16.2 itself is probably not responsible for the observed differences in survival but instead probably reflects a hidden, heterogeneous, but now quantifiable, physiological state that dictates the ability of an organism to deal with the rigors of living.     

Genetic variation in laboratory mice

Characterizing the patterns of genetic variation in an organism provides fundamental insight into the evolutionary history of the organism and defines the scope and nature of studies that must be designed to correlate genotype to phenotype. Given the pre-eminent role of the inbred mouse in biomedical research, considerable effort has been undertaken in recent years to describe more fully the nature and amount of genetic variation among the numerous strains of mice that are in widest use. Here, we discuss recent studies that have contributed to an emerging understanding of the unique variation patterns found in inbred strains of mice and how they have arisen through a combination of natural evolution and human-directed breeding. These preliminary results have ramifications for genetic research into complex biomedical traits and are the basis for the development of future variation resources.     

Genome-wide association study of flowering time and grain yield traits in a worldwide collection of rice germplasm

A high-density haplotype map recently enabled a genome-wide association study (GWAS) in a population of indica subspecies of Chinese rice landraces. Here we extend this methodology to a larger and more diverse sample of 950 worldwide rice varieties, including the Oryza sativa indica and Oryza sativa japonica subspecies, to perform an additional GWAS. We identified a total of 32 new loci associated with flowering time and with ten grain-related traits, indicating that the larger sample increased the power to detect trait-associated variants using GWAS. To characterize various alleles and complex genetic variation, we developed an analytical framework for haplotype-based de novo assembly of the low-coverage sequencing data in rice. We identified candidate genes for 18 associated loci through detailed annotation. This study shows that the integrated approach of sequence-based GWAS and functional genome annotation has the potential to match complex traits to their causal polymorphisms in rice.     

Population genetic structure of variable drug response.

Geographic patterns of genetic variation, including variation at drug metabolizing enzyme (DME) loci and drug targets, indicate that geographic structuring of inter-individual variation in drug response may occur frequently. This raises two questions: how to represent human population genetic structure in the evaluation of drug safety and efficacy, and how to relate this structure to drug response. We address these by (i) inferring the genetic structure present in a heterogeneous sample and (ii) comparing the distribution of DME variants across the inferred genetic clusters of individuals. We find that commonly used ethnic labels are both insufficient and inaccurate representations of the inferred genetic clusters, and that drug-metabolizing profiles, defined by the distribution of DME variants, differ significantly among the clusters. We note, however, that the complexity of human demographic history means that there is no obvious natural clustering scheme, nor an obvious appropriate degree of resolution. Our comparison of drug-metabolizing profiles across the inferred clusters establishes a framework for assessing the appropriate level of resolution in relating genetic structure to drug response.