Genetic clonal diversity predicts progression to esophageal adenocarcinoma

Neoplasms are thought to progress to cancer through genetic instability generating cellular diversity and clonal expansions driven by selection for mutations in cancer genes. Despite advances in the study of molecular biology of cancer genes, relatively little is known about evolutionary mechanisms that drive neoplastic progression. It is unknown, for example, which may be more predictive of future progression of a neoplasm: genetic homogenization of the neoplasm, possibly caused by a clonal expansion, or the accumulation of clonal diversity. Here, in a prospective study, we show that clonal diversity measures adapted from ecology and evolution can predict progression to adenocarcinoma in the premalignant condition known as Barrett's esophagus, even when controlling for established genetic risk factors, including lesions in TP53 (p53; ref. 6) and ploidy abnormalities. Progression to cancer through accumulation of clonal diversity, on which natural selection acts, may be a fundamental principle of neoplasia with important clinical implications.     

Dichotomy of single-nucleotide polymorphism haplotypes in olfactory receptor genes and pseudogenes

Substantial efforts are focused on identifying single-nucleotide polymorphisms (SNPs) throughout the human genome, particularly in coding regions (cSNPs), for both linkage disequilibrium and association studies. Less attention, however, has been directed to the clarification of evolutionary processes that are responsible for the variability in nucleotide diversity among different regions of the genome. We report here the population sequence diversity of genomic segments within a 450-kb cluster of olfactory receptor (OR) genes on human chromosome 17. We found a dichotomy in the pattern of nucleotide diversity between OR pseudogenes and introns on the one hand and the closely interspersed intact genes on the other. We suggest that weak positive selection is responsible for the observed patterns of genetic variation. This is inferred from a lower ratio of polymorphism to divergence in genes compared with pseudogenes or introns, high non-synonymous substitution rates in OR genes, and a small but significant overall reduction in variability in the entire OR gene cluster compared with other genomic regions. The dichotomy among functionally different segments within a short genomic distance requires high recombination rates within this OR cluster. Our work demonstrates the impact of weak positive selection on human nucleotide diversity, and has implications for the evolution of the olfactory repertoire.     

Dopa decarboxylase (Ddc) affects variation in Drosophila longevity

Mutational analyses in model organisms have shown that genes affecting metabolism and stress resistance regulate life span, but the genes responsible for variation in longevity in natural populations are largely unidentified. Previously, we mapped quantitative trait loci (QTLs) affecting variation in longevity between two Drosophila melanogaster strains. Here, we show that the longevity QTL in the 36E;38B cytogenetic interval on chromosome 2 contains multiple closely linked QTLs, including the Dopa decarboxylase (Ddc) locus. Complementation tests to mutations show that Ddc is a positional candidate gene for life span in these strains. Linkage disequilibrium (LD) mapping in a sample of 173 alleles from a single population shows that three common molecular polymorphisms in Ddc account for 15.5% of the genetic contribution to variance in life span from chromosome 2. The polymorphisms are in strong LD, and the effects of the haplotypes on longevity suggest that the polymorphisms are maintained by balancing selection. DDC catalyzes the final step in the synthesis of the neurotransmitters, dopamine and serotonin. Thus, these data implicate variation in the synthesis of bioamines as a factor contributing to natural variation in individual life span.     

Adaptive evolution of the tumour suppressor BRCA1 in humans and chimpanzees. Australian Breast Cancer Family Study

Mutations in BRCA1 (ref. 1) confer an increased risk of female breast cancer. In a genome-wide scan of linkage disequilibrium (LD), a high level of LD was detected among microsatellite markers flanking BRCA1 (ref. 3), raising the prospect that positive natural selection may have acted on this gene. We have used the predictions of evolutionary genetic theory to investigate this further. Using phylogeny-based maximum likelihood analysis of the BRCA1 sequences from primates and other mammals, we found that the ratios of replacement to silent nucleotide substitutions on the human and chimpanzee lineages were not different from one another (P=0.8), were different from those of other primate lineages (P=0.004) and were greater than 1 (P=0.04). This is consistent with the historic occurrence of positive darwinian selection pressure on the BRCA1 protein in the human and chimpanzee lineages. Analysis of genetic variation in a sample of female Australians of Northern European origin showed evidence for Hardy-Weinberg (HW) disequilibrium at polymorphic sites in BRCA1, consistent with the possibility that natural selection is affecting genotype frequencies in modern Europeans. The clustering of between-species variation in the region of the gene encoding the RAD51-interaction domain of BRCA1 suggests the maintenance of genomic integrity as a possible target of selection.     

Natural selection has driven population differentiation in modern humans

The considerable range of observed phenotypic variation in human populations may reflect, in part, distinctive processes of natural selection and adaptation to variable environmental conditions. Although recent genome-wide studies have identified candidate regions under selection, it is not yet clear how natural selection has shaped population differentiation. Here, we have analyzed the degree of population differentiation at 2.8 million Phase II HapMap single-nucleotide polymorphisms. We find that negative selection has globally reduced population differentiation at amino acid-altering mutations, particularly in disease-related genes. Conversely, positive selection has ensured the regional adaptation of human populations by increasing population differentiation in gene regions, primarily at nonsynonymous and 5'-UTR variants. Our analyses identify a fraction of loci that have contributed, and probably still contribute, to the morphological and disease-related phenotypic diversity of current human populations.     

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.     

Dynamic evolution of the innate immune system in Drosophila

The availability of complete genome sequence from 12 Drosophila species presents the opportunity to examine how natural selection has affected patterns of gene family evolution and sequence divergence among different components of the innate immune system. We have identified orthologs and paralogs of 245 Drosophila melanogaster immune-related genes in these recently sequenced genomes. Genes encoding effector proteins, and to a lesser extent genes encoding recognition proteins, are much more likely to vary in copy number across species than genes encoding signaling proteins. Furthermore, we can trace the apparent recent origination of several evolutionarily novel immune-related genes and gene families. Using codon-based likelihood methods, we show that immune-system genes, and especially those encoding recognition proteins, evolve under positive darwinian selection. Positively selected sites within recognition proteins cluster in domains involved in recognition of microorganisms, suggesting that molecular interactions between hosts and pathogens may drive adaptive evolution in the Drosophila immune system.     

Identification of genes involved in Drosophila melanogaster geotaxis,a complex behavioral trait

Identifying the genes involved in polygenic traits has been difficult. In the 1950s and 1960s, laboratory selection experiments for extreme geotaxic behavior in fruit flies established for the first time that a complex behavioral trait has a genetic basis. But the specific genes responsible for the behavior have never been identified using this classical model. To identify the individual genes involved in geotaxic response, we used cDNA microarrays to identify candidate genes and assessed fly lines mutant in these genes for behavioral confirmation. We have thus determined the identities of several genes that contribute to the complex, polygenic behavior of geotaxis.     

Skint1, the prototype of a newly identified immunoglobulin superfamily gene cluster, positively selects epidermal gammadelta T cells

B cells, alphabeta T cells and gammadelta T cells are conserved lymphocyte subtypes encoding their antigen receptors from somatically rearranged genes. alphabeta T cells undergo positive selection in the thymus by engagement of their T cell receptors (TCRs) with self-peptides presented by major histocompatibility complex molecules. The molecules that select gammadelta T cells are unknown. Vgamma5+Vdelta1+ cells comprise 90% of mouse epidermal gammadelta T cells. By mapping and genetic complementation using a strain showing loss of Vgamma5+Vdelta1+ cells due to a failure of thymic selection, we show that this defect is caused by mutation in Skint1, a newly identified gene expressed in thymus and skin that encodes a protein with immunoglobulin-like and transmembrane domains. Skint1 is the prototypic member of a rapidly evolving family of at least 11 genes in mouse, with greatest similarity to the butyrophilin genes. These findings define a new family of proteins mediating key epithelial-immune interactions.     

Evidence for co-evolution of gene order and recombination rate

There is increasing evidence in eukaryotic genomes that gene order is not random, even allowing for tandem duplication. Notably, in numerous genomes, genes of similar expression tend to be clustered. Are there other reasons for clustering of functionally similar genes? If genes are linked to enable genetic, rather than physical clustering, then we also expect that clusters of certain genes might be associated with blocks of reduced recombination rates. Here we show that, in yeast, essential genes are highly clustered and this clustering is independent of clustering of co-expressed genes and of tandem duplications. Adjacent pairs of essential genes are preferentially conserved through evolution. Notably, we also find that clusters of essential genes are in regions of low recombination and that larger clusters have lower recombination rates. These results suggest that selection acts to modify both the fine-scale intragenomic variation in the recombination rate and the distribution of genes and provide evidence for co-evolution of gene order and recombination rate.     

Robustness against mutations in genetic networks of yeast

There are two principal mechanisms that are responsible for the ability of an organism's physiological and developmental processes to compensate for mutations. In the first, genes have overlapping functions, and loss-of-function mutations in one gene will have little phenotypic effect if there are one or more additional genes with similar functions. The second mechanism has its origin in interactions between genes with unrelated functions, and has been documented in metabolic and regulatory gene networks. Here I analyse, on a genome-wide scale, which of these mechanisms of robustness against mutations is more prevalent. I used functional genomics data from the yeast Saccharomyces cerevisiae to test hypotheses related to the following: if gene duplications are mostly responsible for robustness, then a correlation is expected between the similarity of two duplicated genes and the effect of mutations in one of these genes. My results demonstrate that interactions among unrelated genes are the major cause of robustness against mutations. This type of robustness is probably an evolved response of genetic networks to stabilizing selection.     

Genome-wide association study of leaf architecture in the maize nested association mapping population

US maize yield has increased eight-fold in the past 80 years, with half of the gain attributed to selection by breeders. During this time, changes in maize leaf angle and size have altered plant architecture, allowing more efficient light capture as planting density has increased. Through a genome-wide association study (GWAS) of the maize nested association mapping panel, we determined the genetic basis of important leaf architecture traits and identified some of the key genes. Overall, we demonstrate that the genetic architecture of the leaf traits is dominated by small effects, with little epistasis, environmental interaction or pleiotropy. In particular, GWAS results show that variations at the liguleless genes have contributed to more upright leaves. These results demonstrate that the use of GWAS with specially designed mapping populations is effective in uncovering the basis of key agronomic traits.     

Recombination and linkage disequilibrium in Arabidopsis thaliana

Linkage disequilibrium (LD) is a major aspect of the organization of genetic variation in natural populations. Here we describe the genome-wide pattern of LD in a sample of 19 Arabidopsis thaliana accessions using 341,602 non-singleton SNPs. LD decays within 10 kb on average, considerably faster than previously estimated. Tag SNP selection algorithms and 'hide-the-SNP' simulations suggest that genome-wide association mapping will require only 40%-50% of the observed SNPs, a reduction similar to estimates in a sample of African Americans. An Affymetrix genotyping array containing 250,000 SNPs has been designed based on these results; we demonstrate that it should have more than adequate coverage for genome-wide association mapping. The extent of LD is highly variable, and we find clear evidence of recombination hotspots, which seem to occur preferentially in intergenic regions. LD also reflects the action of selection, and it is more extensive between nonsynonymous polymorphisms than between synonymous polymorphisms.     

Global diversity and evidence for coevolution of KIR and HLA

The killer immunoglobulin-like receptor (KIR) gene cluster shows extensive genetic diversity, as do the HLA class I loci, which encode ligands for KIR molecules. We genotyped 1,642 individuals from 30 geographically distinct populations to examine population-level evidence for coevolution of these two functionally related but unlinked gene clusters. We observed strong negative correlations between the presence of activating KIR genes and their corresponding HLA ligand groups across populations, especially KIR3DS1 and its putative HLA-B Bw4-80I ligands (r = -0.66, P = 0.038). In contrast, we observed weak positive relationships between the various inhibitory KIR genes and their ligands. We observed a negative correlation between distance from East Africa and frequency of activating KIR genes and their corresponding ligands, suggesting a balance between selection on HLA and KIR loci. Most KIR-HLA genetic association studies indicate a primary influence of activating KIR-HLA genotypes in disease risk; concomitantly, activating receptor-ligand pairs in this study show the strongest signature of coevolution of these two complex genetic systems as compared with inhibitory receptor-ligand pairs.     

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.     

Evidence that positive selection drives Y-chromosome degeneration in Drosophila miranda

Why does the Y chromosome harbor so few functional loci? Evolutionary theory predicts that Y chromosomes degenerate because they lack genetic recombination. Both positive and negative selection models have been invoked to explain this degeneration, as both can result in the recurrent fixation of linked deleterious mutations on a nonrecombining Y chromosome. To distinguish between these models, I investigated patterns of nucleotide variability along 37 kb of the recently formed neo-Y chromosome in Drosophila miranda. Levels of nucleotide variability on this chromosome are 30 times lower than in highly recombining portions of the genome. Both positive and negative selection models can result in reduced variability levels, but their effects on the frequency spectrum of mutations differ. Using coalescent simulations, I show that the patterns of nucleotide variability on the neo-Y chromosome are unlikely under deleterious mutation models (including background selection and Muller's ratchet) but are expected under recent positive selection. These results implicate positive selection as an important force driving the degeneration of Y chromosomes; adaptation at a few loci, possibly increasing male fitness, occurs at the cost of most other genes on this chromosome.