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.     

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.     

Fusion genes and rearranged genes as a linear function of chromosome aberrations in cancer

Cytogenetic aberrations have been reported in 45,000 human neoplasms. Structural balanced rearrangements are associated with distinct tumor subtypes with remarkable specificity and have been essential for identifying genes involved in tumorigenesis. All balanced rearrangements that have been characterized molecularly act by deregulating a gene in one of the breakpoints or by creating a fusion gene. Because most recurrent aberrations and rearranged genes have been found in hematological disorders, whereas numerous genomic imbalances have been identified in solid tumors, it has become generally accepted that there are pathogenetic differences between these neoplasms. We here show that in every tumor type, the numbers of recurrent balanced chromosome abnormalities, fusion genes and genes rearranged as a consequence of balanced aberrations are simply a function of the number of cases with an abnormal karyotype. Hence, there may not be any fundamental tissue-specific differences in the genetic mechanisms by which neoplasia is initiated.     

Transposition of maize Ac/Ds transposable elements in the yeast Saccharomyces cerevisiae

Excision by transposons is associated with chromosome breaks; generally, host-cell proteins repair this damage, often introducing mutations. Many transposons also use host proteins in the transposition mechanism or in regulation. Transposition in systems lacking host factors that influence the behaviour of these transpositions is useful in determining what those factors are and how they work. In addition, features of transposition and regulation intrinsic to the element itself can be determined. Maize Activator/Dissociation (Ac/Ds) elements transpose in a wide variety of heterologous plants, but their characteristics in these other systems differ from those in maize, including their response to increasing genetic dosage and the types of repair products recovered following excision. Two Arabidopsis thaliana mutants (iae1 and iae2) show increased Ac transposition frequencies. These mutants, and the differences mentioned above, suggest the involvement of host proteins in Ac/Ds activity and potential differences between these proteins among plant species. Here we report that Ac/Ds elements, members of the hAT (hobo, Ac, Tam3) superfamily, transpose in the yeast Saccharomyces cerevisiae, an organism lacking class II ('cut and paste') transposons. This demonstrates that plant-specific proteins are not essential for Ac/Ds transposition. The yeast system is valuable for dissecting the Ac/Ds transposition mechanism and identifying host factors that can influence transposition and the repair of DNA damage induced by Ac/Ds. Mutations caused by Ds excision in yeast suggest formation of a DNA-hairpin intermediate, and reinsertions occur throughout the genome with a frequency similar to that in plants. The high proportion of Ac/Ds reinsertions also makes this system an in vivo mutagenesis and reverse genetics tool in yeast and, presumably, other eukaryotic systems.     

Extensive allelic variation and ultrashort telomeres in senescent human cells

By imposing a limit on the proliferative lifespan of most somatic cells, telomere erosion represents an innate mechanism for tumor suppression and may contribute to age-related disease. A detailed understanding of the pathways that link shortened telomeres to replicative senescence has been severely hindered by the inability of current methods to analyze telomere dynamics in detail. Here we describe single telomere length analysis (STELA), a PCR-based approach that accurately measures the full spectrum of telomere lengths from individual chromosomes. STELA analysis of human XpYp telomeres in fibroblasts identifies several features of telomere biology. We observe bimodal distributions of telomeres in normal fibroblasts; these distributions result from inter-allelic differences of up to 6.5 kb, indicating that unexpectedly large-scale differences in zygotic telomere length are maintained throughout development. Most telomeres shorten in a gradual fashion consistent with simple losses through end replication, and the rates of erosion are independent of allele size. Superimposed on this are occasional, more substantial changes in length, which may be the consequence of additional mutational mechanisms. Notably, some alleles show almost complete loss of TTAGGG repeats at senescence.     

Common variation in three genes, including a noncoding variant in CFH, strongly influences risk of age-related macular degeneration

Age-related macular degeneration (AMD) is a common, late-onset disease with seemingly typical complexity: recurrence ratios for siblings of an affected individual are three- to sixfold higher than in the general population, and family-based analysis has resulted in only modestly significant evidence for linkage. In a case-control study drawn from a US-based population of European descent, we have identified a previously unrecognized common, noncoding variant in CFH, the gene encoding complement factor H, that substantially increases the influence of this locus on AMD, and we have strongly replicated the associations of four other previously reported common alleles in three genes (P values ranging from 10(-6) to 10(-70)). Despite excellent power to detect epistasis, we observed purely additive accumulation of risk from alleles at these genes. We found no differences in association of these loci with major phenotypic categories of advanced AMD. Genotypes at these five common SNPs define a broad spectrum of interindividual disease risk and explain about half of the classical sibling risk of AMD in our study population.     

Allergie et environnement

Allergic symptoms are the result of a predisposed background (very often genetic background) and an environmental exposure. The combinations between these two parameters are very variable, and these variations explain the differences observed by the clinician. For example, for an allergen, one exposed patient can become allergic within a few weeks though another will in a few years. This kind of difference is necessarily explained by the genetic background. To become allergic, we can postulate that a lightly < programmed ⪢ patient would need to be exposed for a longer period that the strongly < programmed ⪢ one.On the other hand, the modifications in our environment probably play a great role in the increasing process of the allergic illness and genetic parameters cannot be involved.The environment will be described, including indoor and outdoor allergens, but also foods and stinging insects.     

Recurrent DNA copy number variation in the laboratory mouse

Different species, populations and individuals vary considerably in the copy number of discrete segments of their genomes. The manner and frequency with which these genetic differences arise over generational time is not well understood. Taking advantage of divergence among lineages sharing a recent common ancestry, we have conducted a genome-wide analysis of spontaneous copy number variation (CNV) in the laboratory mouse. We used high-resolution microarrays to identify 38 CNVs among 14 colonies of the C57BL/6 strain spanning approximately 967 generations of inbreeding, and we examined these loci in 12 additional strains. It is clear from our results that many CNVs arise through a highly nonrandom process: 18 of 38 were the product of recurrent mutation, and rates of change varied roughly four orders of magnitude across different loci. Recurrent CNVs are found throughout the genome, affect 43 genes and fluctuate in copy number over mere hundreds of generations, observations that raise questions about their contribution to natural variation.     

A genome-wide approach accounting for body mass index identifies genetic variants influencing fasting glycemic traits and insulin resistance

Recent genome-wide association studies have described many loci implicated in type 2 diabetes (T2D) pathophysiology and β-cell dysfunction but have contributed little to the understanding of the genetic basis of insulin resistance. We hypothesized that genes implicated in insulin resistance pathways might be uncovered by accounting for differences in body mass index (BMI) and potential interactions between BMI and genetic variants. We applied a joint meta-analysis approach to test associations with fasting insulin and glucose on a genome-wide scale. We present six previously unknown loci associated with fasting insulin at P < 5 × 10(-8) in combined discovery and follow-up analyses of 52 studies comprising up to 96,496 non-diabetic individuals. Risk variants were associated with higher triglyceride and lower high-density lipoprotein (HDL) cholesterol levels, suggesting a role for these loci in insulin resistance pathways. The discovery of these loci will aid further characterization of the role of insulin resistance in T2D pathophysiology.     

Complex haplotypes, copy number polymorphisms and coding variation in two recently divergent mouse strains

Inbred mouse strains provide the foundation for mouse genetics. By selecting for phenotypic features of interest, inbreeding drives genomic evolution and eliminates individual variation, while fixing certain sets of alleles that are responsible for the trait characteristics of the strain. Mouse strains 129Sv (129S5) and C57BL/6J, two of the most widely used inbred lines, diverged from common ancestors within the last century, yet very little is known about the genomic differences between them. By comparative genomic hybridization and sequence analysis of 129S5 short insert libraries, we identified substantial structural variation, a complex fine-scale haplotype pattern with a continuous distribution of diversity blocks, and extensive nucleotide variation, including nonsynonymous coding SNPs and stop codons. Collectively, these genomic changes denote the level and direction of allele fixation that has occurred during inbreeding and provide a basis for defining what makes these mouse strains unique.     

Common deletions and SNPs are in linkage disequilibrium in the human genome

Humans show great variation in phenotypic traits such as height, eye color and susceptibility to disease. Genomic DNA sequence differences among individuals are responsible for the inherited components of these complex traits. Reports suggest that intermediate and large-scale DNA copy number and structural variations are prevalent enough to be an important source of genetic variation between individuals. Because association studies to identify genomic loci associated with particular phenotypic traits have focused primarily on genotyping SNPs, it is important to determine whether common structural polymorphisms are in linkage disequilibrium with common SNPs, and thus can be assessed indirectly in SNP-based studies. Here we examine 100 deletion polymorphisms ranging from 70 bp to 7 kb. We show that common deletions and SNPs ascertained with similar criteria have essentially the same distribution of linkage disequilibrium with surrounding SNPs, indicating that these polymorphisms may share evolutionary history and that most deletion polymorphisms are effectively assayed by proxy in SNP-based association studies.     

Natural variation in light sensitivity of Arabidopsis.

Because plants depend on light for growth, their development and physiology must suit the particular light environment. Plants native to different environments show heritable, apparently adaptive, changes in their response to light. As a first step in unraveling the genetic and molecular basis of these naturally occurring differences, we have characterized intraspecific variation in a light-dependent developmental process-seedling emergence. We examined 141 Arabidopsis thaliana accessions for their response to four light conditions, two hormone conditions and darkness. There was significant variation in all conditions, confirming that Arabidopsis is a rich source of natural genetic diversity. Hierarchical clustering revealed that some accessions had response patterns similar to known photoreceptor mutants, suggesting changes in specific signaling pathways. We found that the unusual far-red response of the Lm-2 accession is due to a single amino-acid change in the phytochrome A (PHYA) protein. This change stabilizes the light-labile PHYA protein in light and causes a 100-fold shift in the threshold for far-red light sensitivity. Purified recombinant Lm-2 PHYA also shows subtle photochemical differences and has a reduced capacity for autophosphorylation. These biochemical changes contrast with previously characterized natural alleles in loci controlling plant development, which result in altered gene expression or loss of gene function.     

Insertion of Inhbb into the Inhba locus rescues the Inhba-null phenotype and reveals new activin functions

The activins (dimers of betaA or betaB subunits, encoded by the genes Inhba and Inhbb, respectively) are TGF-beta superfamily members that have roles in reproduction and development. Whereas mice homozygous for the Inhba-null allele demonstrate disruption of whisker, palate and tooth development, leading to neonatal lethality, homozygous Inhbb-null mice are viable, fertile and have eye defects. To determine if these phenotypes were due to spatiotemporal expression differences of the ligands or disruption of specific ligand-receptor interactions, we replaced the region of Inhba encoding the mature protein with Inhbb, creating the allele Inhbatm2Zuk (hereafter designated InhbaBK). Although the craniofacial phenotypes of the Inhba-null mutation were rescued by the InhbaBK allele, somatic, testicular, genital and hair growth were grossly affected and influenced by the dosage and bioactivity of the allele. Thus, functional compensation within the TGF-beta superfamily can occur if the replacement gene is expressed appropriately. The novel phenotypes in these mice further illustrate the usefulness of insertion strategies for defining protein function.     

Common variants near TARDBP and EGR2 are associated with susceptibility to Ewing sarcoma

Ewing sarcoma, a pediatric tumor characterized by EWSR1-ETS fusions, is predominantly observed in populations of European ancestry. We performed a genome-wide association study (GWAS) of 401 French individuals with Ewing sarcoma, 684 unaffected French individuals and 3,668 unaffected individuals of European descent and living in the United States. We identified candidate risk loci at 1p36.22, 10q21 and 15q15. We replicated these loci in two independent sets of cases and controls. Joint analysis identified associations with rs9430161 (P = 1.4 × 10(-20); odds ratio (OR) = 2.2) located 25 kb upstream of TARDBP, rs224278 (P = 4.0 × 10(-17); OR = 1.7) located 5 kb upstream of EGR2 and, to a lesser extent, rs4924410 at 15q15 (P = 6.6 × 10(-9); OR = 1.5). The major risk haplotypes were less prevalent in Africans, suggesting that these loci could contribute to geographical differences in Ewing sarcoma incidence. TARDBP shares structural similarities with EWSR1 and FUS, which encode RNA binding proteins, and EGR2 is a target gene of EWSR1-ETS. Variants at these loci were associated with expression levels of TARDBP, ADO (encoding cysteamine dioxygenase) and EGR2.     

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.     

Treatment-specific changes in gene expression discriminate in vivo drug response in human leukemia cells

To elucidate the genomics of cellular responses to cancer treatment, we analyzed the expression of over 9,600 human genes in acute lymphoblastic leukemia cells before and after in vivo treatment with methotrexate and mercaptopurine given alone or in combination. Based on changes in gene expression, we identified 124 genes that accurately discriminated among the four treatments. Discriminating genes included those involved in apoptosis, mismatch repair, cell cycle control and stress response. Only 14% of genes that changed when these medications were given as single agents also changed when they were given together. These data indicate that lymphoid leukemia cells of different molecular subtypes share common pathways of genomic response to the same treatment, that changes in gene expression are treatment-specific and that gene expression can illuminate differences in cellular response to drug combinations versus single agents.     

Evaluating coverage of genome-wide association studies

Genome-wide association studies involving hundreds of thousands of SNPs in thousands of cases and controls are now underway. The first of many analytical challenges in these studies involves the choice of SNPs to genotype. It is not practical to construct a different panel of tag SNPs for each study, so the first generation of genome-wide scans will use predefined, commercially available marker panels, which will in part dictate their success or failure. We compare different approaches in use today, and show that although many of them provide substantial coverage of common variation in non-African populations, the precise extent is strongly dependent on the frequencies of alleles of interest and on specific considerations of study design. Overall, despite substantial differences in genotyping technologies, marker selection strategies and number of markers assayed, the first-generation high-throughput platforms all offer similar levels of genome coverage.     

The codon 72 polymorphic variants of p53 have markedly different apoptotic potential

The gene TP53, encoding p53, has a common sequence polymorphism that results in either proline or arginine at amino-acid position 72. This polymorphism occurs in the proline-rich domain of p53, which is necessary for the protein to fully induce apoptosis. We found that in cell lines containing inducible versions of alleles encoding the Pro72 and Arg72 variants, and in cells with endogenous p53, the Arg72 variant induces apoptosis markedly better than does the Pro72 variant. Our data indicate that at least one source of this enhanced apoptotic potential is the greater ability of the Arg72 variant to localize to the mitochondria; this localization is accompanied by release of cytochrome c into the cytosol. These data indicate that the two polymorphic variants of p53 are functionally distinct, and these differences may influence cancer risk or treatment.     

Differential effects of oncogenic K-Ras and N-Ras on proliferation, differentiation and tumor progression in the colon

Kras is commonly mutated in colon cancers, but mutations in Nras are rare. We have used genetically engineered mice to determine whether and how these related oncogenes regulate homeostasis and tumorigenesis in the colon. Expression of K-Ras(G12D) in the colonic epithelium stimulated hyperproliferation in a Mek-dependent manner. N-Ras(G12D) did not alter the growth properties of the epithelium, but was able to confer resistance to apoptosis. In the context of an Apc-mutant colonic tumor, activation of K-Ras led to defects in terminal differentiation and expansion of putative stem cells within the tumor epithelium. This K-Ras tumor phenotype was associated with attenuated signaling through the MAPK pathway, and human colon cancer cells expressing mutant K-Ras were hypersensitive to inhibition of Raf, but not Mek. These studies demonstrate clear phenotypic differences between mutant Kras and Nras, and suggest that the oncogenic phenotype of mutant K-Ras might be mediated by noncanonical signaling through Ras effector pathways.