Common genetic variants account for differences in gene expression among ethnic groups

Variation in DNA sequence contributes to individual differences in quantitative traits, but in humans the specific sequence variants are known for very few traits. We characterized variation in gene expression in cells from individuals belonging to three major population groups. This quantitative phenotype differs significantly between European-derived and Asian-derived populations for 1,097 of 4,197 genes tested. For the phenotypes with the strongest evidence of cis determinants, most of the variation is due to allele frequency differences at cis-linked regulators. The results show that specific genetic variation among populations contributes appreciably to differences in gene expression phenotypes. Populations differ in prevalence of many complex genetic diseases, such as diabetes and cardiovascular disease. As some of these are probably influenced by the level of gene expression, our results suggest that allele frequency differences at regulatory polymorphisms also account for some population differences in prevalence of complex diseases.     

The 'evolvability' of promiscuous protein functions

How proteins with new functions (e.g., drug or antibiotic resistance or degradation of man-made chemicals) evolve in a matter of months or years is still unclear. This ability is dependent on the induction of new phenotypic traits by a small number of mutations (plasticity). But mutations often have deleterious effects on functions that are essential for survival. How are these seemingly conflicting demands met at the single-protein level? Results from directed laboratory evolution experiments indicate that the evolution of a new function is driven by mutations that have little effect on the native function but large effects on the promiscuous functions that serve as starting point. Thus, an evolving protein can initially acquire increased fitness for a new function without losing its original function. Gene duplication and the divergence of a completely new protein may then follow.     

Ancestry and pharmacogenomics of relapse in acute lymphoblastic leukemia

Although five-year survival rates for childhood acute lymphoblastic leukemia (ALL) are now over 80% in most industrialized countries, not all children have benefited equally from this progress. Ethnic differences in survival after childhood ALL have been reported in many clinical studies, with poorer survival observed among African Americans or those with Hispanic ethnicity when compared with European Americans or Asians. The causes of ethnic differences remain uncertain, although both genetic and non-genetic factors are likely important. Interrogating genome-wide germline SNP genotypes in an unselected large cohort of children with ALL, we observed that the component of genomic variation that co-segregated with Native American ancestry was associated with risk of relapse (P = 0.0029) even after adjusting for known prognostic factors (P = 0.017). Ancestry-related differences in relapse risk were abrogated by the addition of a single extra phase of chemotherapy, indicating that modifications to therapy can mitigate the ancestry-related risk of relapse.     

Detection of regulatory variation in mouse genes

Functional polymorphism in genes can be classified as coding variation, altering the amino-acid sequence of the encoded protein, or regulatory variation, affecting the level or pattern of expression of the gene. Coding variation can be recognized directly from DNA sequence, and consequently its frequency and characteristics have been extensively described. By contrast, virtually nothing is known about the extent to which gene regulation varies in populations. Yet it is likely that regulatory variants are important in modulating gene function: alterations in gene regulation have been proposed to influence disease susceptibility and to have been the primary substrate for the evolution of species. Here, we report a systematic study to assess the extent of cis-acting regulatory variation in 69 genes across four inbred mouse strains. We find that at least four of these genes show allelic differences in expression level of 1.5-fold or greater, and that some of these differences are tissue specific. The results show that the impact of regulatory variants can be detected at a significant frequency in a genomic survey and suggest that such variation may have important consequences for organismal phenotype and evolution. The results indicate that larger-scale surveys in both mouse and human could identify a substantial number of genes with common regulatory variation.     

Evidence of widespread selection on standing variation in Europe at height-associated SNPs

Strong signatures of positive selection at newly arising genetic variants are well documented in humans, but this form of selection may not be widespread in recent human evolution. Because many human traits are highly polygenic and partly determined by common, ancient genetic variation, an alternative model for rapid genetic adaptation has been proposed: weak selection acting on many pre-existing (standing) genetic variants, or polygenic adaptation. By studying height, a classic polygenic trait, we demonstrate the first human signature of widespread selection on standing variation. We show that frequencies of alleles associated with increased height, both at known loci and genome wide, are systematically elevated in Northern Europeans compared with Southern Europeans (P < 4.3 × 10(-4)). This pattern mirrors intra-European height differences and is not confounded by ancestry or other ascertainment biases. The systematic frequency differences are consistent with the presence of widespread weak selection (selection coefficients ～10(-3)-10(-5) per allele) rather than genetic drift alone (P < 10(-15)).     

Flowering-time genes modulate meristem determinacy and growth form in Arabidopsis thaliana

Plants have evolved annual and perennial life forms as alternative strategies to adapt reproduction and survival to environmental constraints. In isolated situations, such as islands, woody perennials have evolved repeatedly from annual ancestors. Although the molecular basis of the rapid evolution of insular woodiness is unknown, the molecular difference between perennials and annuals might be rather small, and a change between these life strategies might not require major genetic innovations. Developmental regulators can strongly affect evolutionary variation and genes involved in meristem transitions are good candidates for a switch in growth habit. We found that the MADS box proteins SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) and FRUITFULL (FUL) not only control flowering time, but also affect determinacy of all meristems. In addition, downregulation of both proteins established phenotypes common to the lifestyle of perennial plants, suggesting their involvement in the prevention of secondary growth and longevity in annual life forms.     

Formation of stable epialleles and their paramutation-like interaction in tetraploid Arabidopsis thaliana

Polyploidization is found frequently in plants, and species previously considered to be diploid may show remnants of earlier polyploidization events on closer inspection of their genomes. The success of polyploids may lie in increased genetic redundancy supporting subsequent genetic diversification. Although doubling the genome does not generate diversity per se, recent studies show that rapid genomic rearrangements and changes in DNA modification and gene expression patterns are associated with polyploid formation. But recessive modifications will not become phenotypically apparent in early polyploid generations. Here we show that epialleles in tetraploid plants (but not in diploids) interact in trans and lead to heritable gene silencing persisting after segregation from the inactivating allele. This mechanism, resembling paramutation, leads to the establishment of functional epigenetic homozygosity and, thus, to conversion of new recessive alleles into traits expressed in early polyploid generations. Such interactions probably contribute to rapid adaptation and evolution of polyploid plant species.     

A structural-maintenance-of-chromosomes hinge domain-containing protein is required for RNA-directed DNA methylation

RNA-directed DNA methylation (RdDM) is a process in which dicer-generated small RNAs guide de novo cytosine methylation at the homologous DNA region. To identify components of the RdDM machinery important for Arabidopsis thaliana development, we targeted an enhancer active in meristems for methylation, which resulted in silencing of a downstream GFP reporter gene. This silencing system also features secondary siRNAs, which trigger methylation that spreads beyond the targeted enhancer region. A screen for mutants defective in meristem silencing and enhancer methylation retrieved six dms complementation groups, which included the known factors DRD1 (ref. 3; a SNF2-like chromatin-remodeling protein) and Pol IVb subunits. Additionally, we identified a previously unknown gene DMS3 (At3g49250), encoding a protein similar to the hinge-domain region of structural maintenance of chromosomes (SMC) proteins. This finding implicates a putative chromosome architectural protein that can potentially link nucleic acids in facilitating an RNAi-mediated epigenetic modification involving secondary siRNAs and spreading of DNA methylation.     

The common PPARgamma Pro12Ala polymorphism is associated with decreased risk of type 2 diabetes

Genetic association studies are viewed as problematic and plagued by irreproducibility. Many associations have been reported for type 2 diabetes, but none have been confirmed in multiple samples and with comprehensive controls. We evaluated 16 published genetic associations to type 2 diabetes and related sub-phenotypes using a family-based design to control for population stratification, and replication samples to increase power. We were able to confirm only one association, that of the common Pro12Ala polymorphism in peroxisome proliferator-activated receptor-gamma(PPARgamma) with type 2 diabetes. By analysing over 3,000 individuals, we found a modest (1.25-fold) but significant (P=0.002) increase in diabetes risk associated with the more common proline allele (85% frequency). Moreover, our results resolve a controversy about common variation in PPARgamma. An initial study found a threefold effect, but four of five subsequent publications failed to confirm the association. All six studies are consistent with the odds ratio we describe. The data implicate inherited variation in PPARgamma in the pathogenesis of type 2 diabetes. Because the risk allele occurs at such high frequency, its modest effect translates into a large population attributable risk-influencing as much as 25% of type 2 diabetes in the general population.     

Segregation at three loci explains familial and population risk in Hirschsprung disease

Hirschsprung disease (HSCR), the most common hereditary cause of intestinal obstruction, shows considerable variation and complex inheritance. Coding sequence mutations in RET, GDNF, EDNRB, EDN3 and SOX10 lead to long-segment (L-HSCR) and syndromic HSCR but fail to explain the transmission of the much more common short-segment form (S-HSCR). We conducted a genome scan in families with S-HSCR and identified susceptibility loci at 3p21, 10q11 and 19q12 that seem to be necessary and sufficient to explain recurrence risk and population incidence. The gene at 10q11 is probably RET, supporting its crucial role in all forms of HSCR; however, coding sequence mutations are present in only 40% of linked families, suggesting the importance of noncoding variation. Here we show oligogenic inheritance of S-HSCR, the 3p21 and 19q12 loci as RET-dependent modifiers, and a parent-of-origin effect at RET. This study demonstrates by a complete genetic dissection why the inheritance pattern of S-HSCR is nonmendelian.     

Chromosome-scale selective sweeps shape Caenorhabditis elegans genomic diversity

The nematode Caenorhabditis elegans is central to research in molecular, cell and developmental biology, but nearly all of this research has been conducted on a single strain of C. elegans. Little is known about the population genomic and evolutionary history of this species. We characterized C. elegans genetic variation using high-throughput selective sequencing of a worldwide collection of 200 wild strains and identified 41,188 SNPs. Notably, C. elegans genome variation is dominated by a set of commonly shared haplotypes on four of its six chromosomes, each spanning many megabases. Population genetic modeling showed that this pattern was generated by chromosome-scale selective sweeps that have reduced variation worldwide; at least one of these sweeps probably occurred in the last few hundred years. These sweeps, which we hypothesize to be a result of human activity, have drastically reshaped the global C. elegans population in the recent past.     

PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes

DNA microarrays can be used to identify gene expression changes characteristic of human disease. This is challenging, however, when relevant differences are subtle at the level of individual genes. We introduce an analytical strategy, Gene Set Enrichment Analysis, designed to detect modest but coordinate changes in the expression of groups of functionally related genes. Using this approach, we identify a set of genes involved in oxidative phosphorylation whose expression is coordinately decreased in human diabetic muscle. Expression of these genes is high at sites of insulin-mediated glucose disposal, activated by PGC-1alpha and correlated with total-body aerobic capacity. Our results associate this gene set with clinically important variation in human metabolism and illustrate the value of pathway relationships in the analysis of genomic profiling experiments.