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.     

Population genomics of human gene expression

Genetic variation influences gene expression, and this variation in gene expression can be efficiently mapped to specific genomic regions and variants. Here we have used gene expression profiling of Epstein-Barr virus-transformed lymphoblastoid cell lines of all 270 individuals genotyped in the HapMap Consortium to elucidate the detailed features of genetic variation underlying gene expression variation. We find that gene expression is heritable and that differentiation between populations is in agreement with earlier small-scale studies. A detailed association analysis of over 2.2 million common SNPs per population (5% frequency in HapMap) with gene expression identified at least 1,348 genes with association signals in cis and at least 180 in trans. Replication in at least one independent population was achieved for 37% of cis signals and 15% of trans signals, respectively. Our results strongly support an abundance of cis-regulatory variation in the human genome. Detection of trans effects is limited but suggests that regulatory variation may be the key primary effect contributing to phenotypic variation in humans. We also explore several methodologies that improve the current state of analysis of gene expression variation.     

Regulatory changes underlying expression differences within and between Drosophila species

Differences in gene expression are an important source of phenotypic variation, and can be caused by changes in cis and/or trans regulation. cis-regulatory variants alter allele-specific expression, whereas trans-regulatory variants influence expression of both alleles in a diploid cell. Because of this difference, we hypothesize that natural selection may favor one type of change over the other. Here, we investigate contributions of cis- and trans-regulatory changes to variable intra- and interspecific gene expression using four strains of Drosophila melanogaster, three strains of D. simulans and a total of 78 genes. We show that cis-regulatory changes account for a greater proportion of the expression differences observed between rather than within species. These data are inconsistent with a neutral model assuming equal probabilities of fixation for cis- and trans-regulatory polymorphisms, suggesting that natural selection influences the molecular mechanisms underlying divergent gene expression. Specifically, cis-regulatory changes seem to accumulate preferentially over time.     

Screening a large reference sample to identify very low frequency sequence variants: comparisons between two genes

Most human sequence variation is in the form of single-nucleotide polymorphisms (SNPs). It has been proposed that coding-region SNPs (cSNPs) be used for direct association studies to determine the genetic basis of complex traits. The success of such studies depends on the frequency of disease-associated alleles, and their distribution in different ethnic populations. If disease-associated alleles are frequent in most populations, then direct genotyping of candidate variants could show robust associations in manageable study samples. This approach is less feasible if the genetic risk from a given candidate gene is due to many infrequent alleles. Previous studies of several genes demonstrated that most variants are relatively infrequent (<0.05). These surveys genotyped small samples (n<75) and thus had limited ability to identify rare alleles. Here we evaluate the prevalence and distribution of such rare alleles by genotyping an ethnically diverse reference sample that is more than six times larger than those used in previous studies (n=450). We screened for variants in the complete coding sequence and intron-exon junctions of two candidate genes for neuropsychiatric phenotypes: SLC6A4, encoding the serotonin transporter; and SLC18A2, encoding the vesicular monoamine transporter. Both genes have unique roles in neuronal transmission, and variants in either gene might be associated with neurobehavioral phenotypes.     

The genetic basis for differences in leaf form between Arabidopsis thaliana and its wild relative Cardamine hirsuta

A key question in biology is how differences in gene function or regulation produce new morphologies during evolution. Here we investigate the genetic basis for differences in leaf form between two closely related plant species, Arabidopsis thaliana and Cardamine hirsuta. We report that in C. hirsuta, class I KNOTTED1-like homeobox (KNOX) proteins are required in the leaf to delay cellular differentiation and produce a dissected leaf form, in contrast to A. thaliana, in which KNOX exclusion from leaves results in a simple leaf form. These differences in KNOX expression arise through changes in the activity of upstream gene regulatory sequences. The function of ASYMMETRIC LEAVES1/ROUGHSHEATH2/PHANTASTICA (ARP) proteins to repress KNOX expression is conserved between the two species, but in C. hirsuta the ARP-KNOX regulatory module controls new developmental processes in the leaf. Thus, evolutionary tinkering with KNOX regulation, constrained by ARP function, may have produced diverse leaf forms by modulating growth and differentiation patterns in developing leaf primordia.     

In vivo RNA interference demonstrates a role for Nramp1 in modifying susceptibility to type 1 diabetes

Type 1 diabetes is an autoimmune disease influenced by multiple genetic loci. Although more than 20 insulin-dependent diabetes (Idd) loci have been implicated in the nonobese diabetic (NOD) mouse model, few causal gene variants have been identified. Here we show that RNA interference (RNAi) can be used to probe candidate genes in this disease model. Slc11a1 encodes a phagosomal ion transporter, Nramp1, that affects resistance to intracellular pathogens and influences antigen presentation. This gene is the strongest candidate among the 42 genes in the Idd5.2 region; a naturally occurring mutation in the protective Idd5.2 haplotype results in loss of function of the Nramp1 protein. Using lentiviral transgenesis, we generated NOD mice in which Slc11a1 is silenced by RNAi. Silencing reduced the frequency of type 1 diabetes, mimicking the protective Idd5.2 region. Our results demonstrate a role for Slc11a1 in modifying susceptibility to type 1 diabetes and illustrate that RNAi can be used to study causal genes in a mammalian model organism.     

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.     

Characterization of single-nucleotide polymorphisms in coding regions of human genes.

A major goal in human genetics is to understand the role of common genetic variants in susceptibility to common diseases. This will require characterizing the nature of gene variation in human populations, assembling an extensive catalogue of single-nucleotide polymorphisms (SNPs) in candidate genes and performing association studies for particular diseases. At present, our knowledge of human gene variation remains rudimentary. Here we describe a systematic survey of SNPs in the coding regions of human genes. We identified SNPs in 106 genes relevant to cardiovascular disease, endocrinology and neuropsychiatry by screening an average of 114 independent alleles using 2 independent screening methods. To ensure high accuracy, all reported SNPs were confirmed by DNA sequencing. We identified 560 SNPs, including 392 coding-region SNPs (cSNPs) divided roughly equally between those causing synonymous and non-synonymous changes. We observed different rates of polymorphism among classes of sites within genes (non-coding, degenerate and non-degenerate) as well as between genes. The cSNPs most likely to influence disease, those that alter the amino acid sequence of the encoded protein, are found at a lower rate and with lower allele frequencies than silent substitutions. This likely reflects selection acting against deleterious alleles during human evolution. The lower allele frequency of missense cSNPs has implications for the compilation of a comprehensive catalogue, as well as for the subsequent application to disease association.     

Disruption of an imprinted gene cluster by a targeted chromosomal translocation in mice

Genomic imprinting is an epigenetic process in which the activity of a gene is determined by its parent of origin. Mechanisms governing genomic imprinting are just beginning to be understood. However, the tendency of imprinted genes to exist in chromosomal clusters suggests a sharing of regulatory elements. To better understand imprinted gene clustering, we disrupted a cluster of imprinted genes on mouse distal chromosome 7 using the Cre/loxP recombination system. In mice carrying a site-specific translocation separating Cdkn1c and Kcnq1, imprinting of the genes retained on chromosome 7, including Kcnq1, Kcnq1ot1, Ascl2, H19 and Igf2, is unaffected, demonstrating that these genes are not regulated by elements near or telomeric to Cdkn1c. In contrast, expression and imprinting of the translocated Cdkn1c, Slc22a1l and Tssc3 on chromosome 11 are affected, consistent with the hypothesis that elements regulating both expression and imprinting of these genes lie within or proximal to Kcnq1. These data support the proposal that chromosomal abnormalities, including translocations, within KCNQ1 that are associated with the human disease Beckwith-Wiedemann syndrome (BWS) may disrupt CDKN1C expression. These results underscore the importance of gene clustering for the proper regulation of imprinted genes.     

Regulatory evolution across the protein interaction network

Protein-protein interactions may impose constraints on both structural and regulatory evolution. Here we show that protein-protein interactions are negatively associated with evolutionary variation in gene expression. Moreover, interacting proteins have similar levels of variation in expression, and their expression levels are positively correlated across strains. Our results suggest that interacting proteins undergo similar evolutionary dynamics, and that their expression levels are evolutionarily coupled. These patterns hold for organisms as diverse as budding yeast and fruit flies.