Genetics of gene expression in primary immune cells identifies cell type-specific master regulators and roles of HLA alleles

Trans-acting genetic variants have a substantial, albeit poorly characterized, role in the heritable determination of gene expression. Using paired purified primary monocytes and B cells, we identify new predominantly cell type-specific cis and trans expression quantitative trait loci (eQTLs), including multi-locus trans associations to LYZ and KLF4 in monocytes and B cells, respectively. Additionally, we observe a B cell-specific trans association of rs11171739 at 12q13.2, a known autoimmune disease locus, with IP6K2 (P = 5.8 × 10(-15)), PRIC285 (P = 3.0 × 10(-10)) and an upstream region of CDKN1A (P = 2 × 10(-52)), suggesting roles for cell cycle regulation and peroxisome proliferator-activated receptor γ (PPARγ) signaling in autoimmune pathogenesis. We also find that specific human leukocyte antigen (HLA) alleles form trans associations with the expression of AOAH and ARHGAP24 in monocytes but not in B cells. In summary, we show that mapping gene expression in defined primary cell populations identifies new cell type-specific trans-regulated networks and provides insights into the genetic basis of disease susceptibility.     

Integrated transcriptional profiling and linkage analysis for identification of genes underlying disease

Integration of genome-wide expression profiling with linkage analysis is a new approach to identifying genes underlying complex traits. We applied this approach to the regulation of gene expression in the BXH/HXB panel of rat recombinant inbred strains, one of the largest available rodent recombinant inbred panels and a leading resource for genetic analysis of the highly prevalent metabolic syndrome. In two tissues important to the pathogenesis of the metabolic syndrome, we mapped cis- and trans-regulatory control elements for expression of thousands of genes across the genome. Many of the most highly linked expression quantitative trait loci are regulated in cis, are inherited essentially as monogenic traits and are good candidate genes for previously mapped physiological quantitative trait loci in the rat. By comparative mapping we generated a data set of 73 candidate genes for hypertension that merit testing in human populations. Mining of this publicly available data set is expected to lead to new insights into the genes and regulatory pathways underlying the extensive range of metabolic and cardiovascular disease phenotypes that segregate in these recombinant inbred strains.     

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.     

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.     

Fine mapping of regulatory loci for mammalian gene expression using radiation hybrids

We mapped regulatory loci for nearly all protein-coding genes in mammals using comparative genomic hybridization and expression array measurements from a panel of mouse-hamster radiation hybrid cell lines. The large number of breaks in the mouse chromosomes and the dense genotyping of the panel allowed extremely sharp mapping of loci. As the regulatory loci result from extra gene dosage, we call them copy number expression quantitative trait loci, or ceQTLs. The -2log10P support interval for the ceQTLs was <150 kb, containing an average of <2-3 genes. We identified 29,769 trans ceQTLs with -log10P > 4, including 13 hotspots each regulating >100 genes in trans. Further, this work identifies 2,761 trans ceQTLs harboring no known genes, and provides evidence for a mode of gene expression autoregulation specific to the X chromosome.     

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.     

Common variants at CDKAL1 and KLF9 are associated with body mass index in east Asian populations

Obesity is a disorder with a complex genetic etiology, and its epidemic is a worldwide problem. Although multiple genetic loci associated with body mass index, the most common measure of obesity, have been identified in European populations, few studies have focused on Asian populations. Here we report a genome-wide association study and replication studies with 62,245 east Asian subjects, which identified two new body mass index-associated loci in the CDKAL1 locus at 6p22 (rs2206734, P = 1.4 × 10(-11)) and the KLF9 locus at 9q21 (rs11142387, P = 1.3 × 10(-9)), as well as several previously reported loci (the SEC16B, BDNF, FTO, MC4R and GIPR loci, P < 5.0 × 10(-8)). We subsequently performed gene-gene interaction analyses and identified an interaction (P = 2.0 × 10(-8)) between a SNP in the KLF9 locus (rs11142387) and one in the MSTN (also known as GDF8) locus at 2q32 (rs13034723). These findings should provide useful insights into the etiology of obesity.     

Common variants on 1p36 and 1q42 are associated with cutaneous basal cell carcinoma but not with melanoma or pigmentation traits

To search for new sequence variants that confer risk of cutaneous basal cell carcinoma (BCC), we conducted a genome-wide SNP association study of 930 Icelanders with BCC and 33,117 controls. After analyzing 304,083 SNPs, we observed signals from loci at 1p36 and 1q42, and replicated these associations in additional sample sets from Iceland and Eastern Europe. Overall, the most significant signals were from rs7538876 on 1p36 (OR = 1.28, P = 4.4 x 10(-12)) and rs801114 on 1q42 (OR = 1.28, P = 5.9 x 10(-12)). The 1p36 locus contains the candidate genes PADI4, PADI6, RCC2 and ARHGEF10L, and the gene nearest to the 1q42 locus is the ras-homolog RHOU. Neither locus was associated with fair pigmentation traits that are known risk factors for BCC, and no risk was observed for melanoma. Approximately 1.6% of individuals of European ancestry are homozygous for both variants, and their estimated risk of BCC is 2.68 times that of noncarriers.     

A genome-wide association study in Europeans and South Asians identifies five new loci for coronary artery disease

Genome-wide association studies have identified 11 common variants convincingly associated with coronary artery disease (CAD)1??, a modest number considering the apparent heritability of CAD?. All of these variants have been discovered in European populations. We report a meta-analysis of four large genome-wide association studies of CAD, with ～575,000 genotyped SNPs in a discovery dataset comprising 15,420 individuals with CAD (cases) (8,424 Europeans and 6,996 South Asians) and 15,062 controls. There was little evidence for ancestry-specific associations, supporting the use of combined analyses. Replication in an independent sample of 21,408 cases and 19,185 controls identified five loci newly associated with CAD (P < 5 × 10?? in the combined discovery and replication analysis): LIPA on 10q23, PDGFD on 11q22, ADAMTS7-MORF4L1 on 15q25, a gene rich locus on 7q22 and KIAA1462 on 10p11. The CAD-associated SNP in the PDGFD locus showed tissue-specific cis expression quantitative trait locus effects. These findings implicate new pathways for CAD susceptibility.     

Dense genotyping identifies and localizes multiple common and rare variant association signals in celiac disease

Using variants from the 1000 Genomes Project pilot European CEU dataset and data from additional resequencing studies, we densely genotyped 183 non-HLA risk loci previously associated with immune-mediated diseases in 12,041 individuals with celiac disease (cases) and 12,228 controls. We identified 13 new celiac disease risk loci reaching genome-wide significance, bringing the number of known loci (including the HLA locus) to 40. We found multiple independent association signals at over one-third of these loci, a finding that is attributable to a combination of common, low-frequency and rare genetic variants. Compared to previously available data such as those from HapMap3, our dense genotyping in a large sample collection provided a higher resolution of the pattern of linkage disequilibrium and suggested localization of many signals to finer scale regions. In particular, 29 of the 54 fine-mapped signals seemed to be localized to single genes and, in some instances, to gene regulatory elements. Altogether, we define the complex genetic architecture of the risk regions of and refine the risk signals for celiac disease, providing the next step toward uncovering the causal mechanisms of the disease.     

Functional screening of an asthma QTL in YAC transgenic mice.

Many quantitative trait loci (QTLs) contributing to genetically complex conditions have been discovered, but few causative genes have been identified. This is mainly due to the large size of QTLs and the subtle connection between genotype and quantitative phenotype associated with these conditions. Transgenic mice have been successfully used to analyse well-characterized genes suspected of contributing to quantitative traits. Although this approach is powerful for examining one gene at a time, it can be impractical for surveying the large genomic intervals containing many genes that are typically associated with QTLs. To screen for genes contributing to an asthma QTL mapped to human chromosome 5q3 (refs 6,7), we characterized a panel of large-insert 5q31 transgenics based on studies demonstrating that altering gene dosage frequently affects quantitative phenotypes normally influenced by that gene. This panel of human YAC transgenics, propagating a 1-Mb interval of chromosome 5q31 containing 6 cytokine genes and 17 partially characterized genes, was screened for quantitative changes in several asthma-associated phenotypes. Multiple independent transgenic lines with altered IgE response to antigen treatment shared a 180-kb region containing 5 genes, including those encoding human interleukin 4 (IL4) and interleukin 13 (IL13 ), which induce IgE class switching in B cells. Further analysis of these mice and mice transgenic for mouse Il4 and Il13 demonstrated that moderate changes in Il4 and Il13 expression affect asthma-associated phenotypes in vivo. This functional screen of large-insert transgenics enabled us to identify genes that influence the QTL phenotype in vivo.     

Efficient mapping of mendelian traits in dogs through genome-wide association

With several hundred genetic diseases and an advantageous genome structure, dogs are ideal for mapping genes that cause disease. Here we report the development of a genotyping array with approximately 27,000 SNPs and show that genome-wide association mapping of mendelian traits in dog breeds can be achieved with only approximately 20 dogs. Specifically, we map two traits with mendelian inheritance: the major white spotting (S) locus and the hair ridge in Rhodesian ridgebacks. For both traits, we map the loci to discrete regions of <1 Mb. Fine-mapping of the S locus in two breeds refines the localization to a region of approximately 100 kb contained within the pigmentation-related gene MITF. Complete sequencing of the white and solid haplotypes identifies candidate regulatory mutations in the melanocyte-specific promoter of MITF. Our results show that genome-wide association mapping within dog breeds, followed by fine-mapping across multiple breeds, will be highly efficient and generally applicable to trait mapping, providing insights into canine and human health.     

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.     

Network modeling links breast cancer susceptibility and centrosome dysfunction

Many cancer-associated genes remain to be identified to clarify the underlying molecular mechanisms of cancer susceptibility and progression. Better understanding is also required of how mutations in cancer genes affect their products in the context of complex cellular networks. Here we have used a network modeling strategy to identify genes potentially associated with higher risk of breast cancer. Starting with four known genes encoding tumor suppressors of breast cancer, we combined gene expression profiling with functional genomic and proteomic (or 'omic') data from various species to generate a network containing 118 genes linked by 866 potential functional associations. This network shows higher connectivity than expected by chance, suggesting that its components function in biologically related pathways. One of the components of the network is HMMR, encoding a centrosome subunit, for which we demonstrate previously unknown functional associations with the breast cancer-associated gene BRCA1. Two case-control studies of incident breast cancer indicate that the HMMR locus is associated with higher risk of breast cancer in humans. Our network modeling strategy should be useful for the discovery of additional cancer-associated genes.