Genome-wide Characterization of Shared and Distinct Genetic Components that Influence Blood Lipid Levels in Ethnically Diverse Human Populations

Blood lipid concentrations are heritable risk factors associated with atherosclerosis and cardiovascular diseases. Lipid traits exhibit considerable variation among populations of distinct ancestral origin as well as between individuals within a population. We performed association analyses to identify genetic loci influencing lipid concentrations in African American and Hispanic American women in the Women’s Health Initiative SNP Health Association Resource. We validated one African-specific high-density lipoprotein cholesterol locus at CD36 as well as 14 known lipid loci that have been previously implicated in studies of European populations. Moreover, we demonstrate striking similarities in genetic architecture (loci influencing the trait, direction and magnitude of genetic effects, and proportions of phenotypic variation explained) of lipid traits across populations. In particular, we found that a disproportionate fraction of lipid variation in African Americans and Hispanic Americans can be attributed to genomic loci exhibiting statistical evidence of association in Europeans, even though the precise genes and variants remain unknown. At the same time, we found substantial allelic heterogeneity within shared loci, characterized both by population-specific rare variants and variants shared among multiple populations that occur at disparate frequencies. The allelic heterogeneity emphasizes the importance of including diverse populations in future genetic association studies of complex traits such as lipids; furthermore, the overlap in lipid loci across populations of diverse ancestral origin argues that additional knowledge can be gleaned from multiple populations.     

A stochastic expectation and maximization algorithm for detecting quantitative trait-associated genes

MOTIVATION: Most biological traits may be correlated with the underlying gene expression patterns that are partially determined by DNA sequence variation. The correlations between gene expressions and quantitative traits are essential for understanding the functions of genes and dissecting gene regulatory networks. RESULTS: In the present study, we adopted a novel statistical method, called the stochastic expectation and maximization (SEM) algorithm, to analyze the associations between gene expression levels and quantitative trait values and identify genetic loci controlling the gene expression variations. In the first step, gene expression levels measured from microarray experiments were assigned to two different clusters based on the strengths of their association with the phenotypes of a quantitative trait under investigation. In the second step, genes associated with the trait were mapped to genetic loci of the genome. Because gene expressions are quantitative, the genetic loci controlling the expression traits are called expression quantitative trait loci. We applied the same SEM algorithm to a real dataset collected from a barley genetic experiment with both quantitative traits and gene expression traits. For the first time, we identified genes associated with eight agronomy traits of barley. These genes were then mapped to seven chromosomes of the barley genome. The SEM algorithm and the result of the barley data analysis are useful to scientists in the areas of bioinformatics and plant breeding. Availability and implementation: The R program for the SEM algorithm can be downloaded from our website: http://www.statgen.ucr.edu.     

Methylation of CpG loci in 5'-flanking region alters steady-state expression of adenomatous polyposis coli gene in colon cancer cell lines

The APC genetic locus has been linked to the tumorigenesis and progression of colorectal cancer, although the precise mechanism of its involvement in this disease remains unknown. We used high sensitivity mapping of the methylated cytosine, Northern blot analysis and immunocytochemical staining in six colorectal cancer cell lines (DLD-1, SW480, Colo320, HT29, WiDr, and Colo201) to examine the relationship between the methylation status of the CpG loci in the 5'-flanking region of the APC gene and its expression. APC mRNA expression levels determined by Northern blot analysis correlated well with APC protein levels visualized by immunocytochemistry. In these colorectal cancer cell lines, no major genetic alterations of the APC gene, such as amplification or deletion, were detected. Analysis of the epigenetic control of APC gene expression in these lines revealed that methylation of the CpG loci in the 5'-untranslated region of APC mRNA repressed steady-state expression of the gene. Furthermore, epigenetic alteration of the APC gene was independent of the APC protein truncation and CpG methylation of the hMLH1 promoter. Although less eminent than protein truncation by point mutation within the coding region of the APC gene, epigenetic alteration suppressing APC gene expression may significantly contribute to oncogenesis and the progression of colorectal cancer.     

Identification of genetic variants contributing to cisplatin-induced cytotoxicity by use of a genomewide approach

Cisplatin, a platinating agent commonly used to treat several cancers, is associated with nephrotoxicity, neurotoxicity, and ototoxicity, which has hindered its utility. To gain a better understanding of the genetic variants associated with cisplatin-induced toxicity, we present a stepwise approach integrating genotypes, gene expression, and sensitivity of HapMap cell lines to cisplatin. Cell lines derived from 30 trios of European descent (CEU) and 30 trios of African descent (YRI) were used to develop a preclinical model to identify genetic variants and gene expression that contribute to cisplatin-induced cytotoxicity in two different populations. Cytotoxicity was determined as cell-growth inhibition at increasing concentrations of cisplatin for 48 h. Gene expression in 176 HapMap cell lines (87 CEU and 89 YRI) was determined using the Affymetrix GeneChip Human Exon 1.0 ST Array. We identified six, two, and nine representative SNPs that contribute to cisplatin-induced cytotoxicity through their effects on 8, 2, and 16 gene expressions in the combined, Centre d'Etude du Polymorphisme Humain (CEPH), and Yoruban populations, respectively. These genetic variants contribute to 27%, 29%, and 45% of the overall variation in cell sensitivity to cisplatin in the combined, CEPH, and Yoruban populations, respectively. Our whole-genome approach can be used to elucidate the expression of quantitative trait loci contributing to a wide range of cellular phenotypes.     

Molecular-Marker-Facilitated Investigations of Quantitative-Trait Loci in Maize. I. Numbers, Genomic Distribution and Types of Gene Action

Individual genetic factors which underlie variation in quantitative traits of maize were investigated in each of two F2 populations by examining the mean trait expressions of genotypic classes at each of 17-20 segregating marker loci. It was demonstrated that the trait expression of marker locus classes could be interpreted in terms of genetic behavior at linked quantitative trait loci (QTLs). For each of 82 traits evaluated, QTLs were detected and located to genomic sites. The numbers of detected factors varied according to trait, with the average trait significantly influenced by almost two-thirds of the marked genomic sites. Most of the detected associations between marker loci and quantitative traits were highly significant, and could have been detected with fewer than the 1800-1900 plants evaluated in each population. The cumulative, simple effects of marker-linked regions of the genome explained between 8 and 40% of the phenotypic variation for a subset of 25 traits evaluated. Single marker loci accounted for between 0.3% and 16% of the phenotypic variation of traits. Individual plant heterozygosity, as measured by marker loci, was significantly associated with variation in many traits. The apparent types of gene action at the QTLs varied both among traits and between loci for given traits, although overdominance appeared frequently, especially for yield-related traits. The prevalence of apparent overdominance may reflect the effects of multiple QTLs within individual marker-linked regions, a situation which would tend to result in overestimation of dominance. Digenic epistasis did not appear to be important in determining the expression of the quantitative traits evaluated. Examination of the effects of marked regions on the expression of pairs of traits suggests that genomic regions vary in the direction and magnitudes of their effects on trait correlations, perhaps providing a means of selecting to dissociate some correlated traits. Marker-facilitated investigations appear to provide a powerful means of examining aspects of the genetic control of quantitative traits. Modifications of the methods employed herein will allow examination of the stability of individual gene effects in varying genetic backgrounds and environments.     

Revealing the architecture of gene regulation: the promise of eQTL studies

Expression quantitative trait loci (eQTL) mapping studies have become a widely used tool for identifying genetic variants that affect gene regulation. In these studies, expression levels are viewed as quantitative traits, and gene expression phenotypes are mapped to particular genomic loci by combining studies of variation in gene expression patterns with genome-wide genotyping. Results from recent eQTL mapping studies have revealed substantial heritable variation in gene expression within and between populations. In many cases, genetic factors that influence gene expression levels can be mapped to proximal (putatively cis) eQTLs and, less often, to distal (putatively trans) eQTLs. Beyond providing great insight into the biology of gene regulation, a combination of eQTL studies with results from traditional linkage or association studies of human disease may help predict a specific regulatory role for polymorphic sites previously associated with disease.     

Tissue-Specific Effects of Genetic and Epigenetic Variation on Gene Regulation and Splicing

Understanding how genetic variation affects distinct cellular phenotypes, such as gene expression levels, alternative splicing and DNA methylation levels, is essential for better understanding of complex diseases and traits. Furthermore, how inter-individual variation of DNA methylation is associated to gene expression is just starting to be studied. In this study, we use the GenCord cohort of 204 newborn Europeans’ lymphoblastoid cell lines, T-cells and fibroblasts derived from umbilical cords. The samples were previously genotyped for 2.5 million SNPs, mRNA-sequenced, and assayed for methylation levels in 482,421 CpG sites. We observe that methylation sites associated to expression levels are enriched in enhancers, gene bodies and CpG island shores. We show that while the correlation between DNA methylation and gene expression can be positive or negative, it is very consistent across cell-types. However, this epigenetic association to gene expression appears more tissue-specific than the genetic effects on gene expression or DNA methylation (observed in both sharing estimations based on P-values and effect size correlations between cell-types). This predominance of genetic effects can also be reflected by the observation that allele specific expression differences between individuals dominate over tissue-specific effects. Additionally, we discover genetic effects on alternative splicing and interestingly, a large amount of DNA methylation correlating to alternative splicing, both in a tissue-specific manner. The locations of the SNPs and methylation sites involved in these associations highlight the participation of promoter proximal and distant regulatory regions on alternative splicing. Overall, our results provide high-resolution analyses showing how genome sequence variation has a broad effect on cellular phenotypes across cell-types, whereas epigenetic factors provide a secondary layer of variation that is more tissue-specific. Furthermore, the details of how this tissue-specificity may vary across inter-relations of molecular traits, and where these are occurring, can yield further insights into gene regulation and cellular biology as a whole.     

Genetic and Epigenetic Regulation of TOX3 Expression in Breast Cancer

Genome wide association studies (GWAS) have identified low penetrance and high frequency single nucleotide polymorphisms (SNPs) that contribute to genetic susceptibility of breast cancer. The SNPs at 16q12, close to the TOX3 and CASC16 genes, represent one of the susceptibility loci identified by GWAS, showing strong evidence for breast cancer association across various populations. To examine molecular mechanisms of TOX3 regulation in breast cancer, we investigated both genetic and epigenetic factors using cell lines and datasets derived from primary breast tumors available through The Cancer Genome Atlas (TCGA). TOX3 expression is highly up-regulated in luminal subtype tumors compared to normal breast tissues or basal-like tumors. Expression quantitative trait loci (eQTL) analyses revealed significant associations of rs3803662 and rs4784227 genotypes with TOX3 expression in breast tumors. Bisulfite sequencing of four CpG islands in the TOX3 promoter showed a clear difference between luminal and basal-like cancer cell lines. 5-Aza-2’-deoxycytidine treatment of a basal-like cancer cell line increased expression of TOX3. TCGA dataset verified significantly lower levels of methylation of the promoter in luminal breast tumors with an inverse correlation between methylation and expression of TOX3. Methylation QTL (mQTL) analyses showed a weak or no correlation of rs3803662 or rs4784227 with TOX3 promoter methylation in breast tumors, indicating an independent relationship between the genetic and epigenetic events. These data suggest a complex system of TOX3 regulation in breast tumors, driven by germline variants and somatic epigenetic modifications in a subtype specific manner.     

Gene-expression variation within and among human populations

Understanding patterns of gene-expression variation within and among human populations will provide important insights into the molecular basis of phenotypic diversity and the interpretation of patterns of expression variation in disease. However, little is known about how gene-expression variation is apportioned within and among human populations. Here, we characterize patterns of natural gene-expression variation in 16 individuals of European and African ancestry. We find extensive variation in gene-expression levels and estimate that approximately 83% of genes are differentially expressed among individuals and that approximately 17% of genes are differentially expressed among populations. By decomposing total gene-expression variation into within- versus among-population components, we find that most expression variation is due to variation among individuals rather than among populations, which parallels observations of extant patterns of human genetic variation. Finally, we performed allele-specific quantitative polymerase chain reaction to demonstrate that cis-regulatory variation in the lymphocyte adaptor protein (SH2B adapter protein 3) contributes to differential expression between European and African samples. These results provide the first insight into how human population structure manifests itself in gene-expression levels and will help guide the search for regulatory quantitative trait loci.     

Identification, replication, and fine-mapping of Loci associated with adult height in individuals of african ancestry

Adult height is a classic polygenic trait of high heritability (h ² ∼0.8). More than 180 single nucleotide polymorphisms (SNPs), identified mostly in populations of European descent, are associated with height. These variants convey modest effects and explain ∼10% of the variance in height. Discovery efforts in other populations, while limited, have revealed loci for height not previously implicated in individuals of European ancestry. Here, we performed a meta-analysis of genome-wide association (GWA) results for adult height in 20,427 individuals of African ancestry with replication in up to 16,436 African Americans. We found two novel height loci (Xp22-rs12393627, P = 3.4×10⁻¹² and 2p14-rs4315565, P = 1.2×10⁻⁸). As a group, height associations discovered in European-ancestry samples replicate in individuals of African ancestry (P = 1.7×10⁻⁴ for overall replication). Fine-mapping of the European height loci in African-ancestry individuals showed an enrichment of SNPs that are associated with expression of nearby genes when compared to the index European height SNPs (P<0.01). Our results highlight the utility of genetic studies in non-European populations to understand the etiology of complex human diseases and traits.     

Estimating African American admixture proportions by use of population-specific alleles.

We analyzed the European genetic contribution to 10 populations of African descent in the United States (Maywood, Illinois; Detroit; New York; Philadelphia; Pittsburgh; Baltimore; Charleston, South Carolina; New Orleans; and Houston) and in Jamaica, using nine autosomal DNA markers. These markers either are population-specific or show frequency differences >45% between the parental populations and are thus especially informative for admixture. European genetic ancestry ranged from 6.8% (Jamaica) to 22.5% (New Orleans). The unique utility of these markers is reflected in the low variance associated with these admixture estimates (SEM 1.3%-2.7%). We also estimated the male and female European contribution to African Americans, on the basis of informative mtDNA (haplogroups H and L) and Y Alu polymorphic markers. Results indicate a sex-biased gene flow from Europeans, the male contribution being substantially greater than the female contribution. mtDNA haplogroups analysis shows no evidence of a significant maternal Amerindian contribution to any of the 10 populations. We detected significant nonrandom association between two markers located 22 cM apart (FY-null and AT3), most likely due to admixture linkage disequilibrium created in the interbreeding of the two parental populations. The strength of this association and the substantial genetic distance between FY and AT3 emphasize the importance of admixed populations as a useful resource for mapping traits with different prevalence in two parental populations.     

Change in maternal environment induced by cross-fostering alters genetic and epigenetic effects on complex traits in mice

The interaction between maternally provided environment and offspring genotype is a major determinant of offspring development and fitness in many organisms. Recent research has demonstrated that not only genetic effects, but also epigenetic effects may be subject to modifications by the maternal environment. Genomic imprinting resulting in parent-of-origin-dependent gene expression is among the best studied of epigenetic effects. However, very little is known about the degree to which genomic imprinting effects can be modulated by the maternally provided environment, which has important implications for phenotypic plasticity. In this study, we investigated this unresolved question using a cross-fostering design in which mouse pups were nursed by either their own or an unrelated mother. We scanned the entire genome to search for quantitative trait loci whose effects depend on cross-fostering and detected 10 of such loci. Of the 10 loci, 4 showed imprinting by cross-foster interactions. In most cases, the interaction effect was due to the presence of an effect in either cross-fostered or non-cross-fostered animals. Our results demonstrate that genomic imprinting effects may often be modified by the maternal environment and that such interactions can impact key fitness-related traits suggesting a greater plasticity of genomic imprinting than previously assumed.     

The evolution of genetic architectures underlying quantitative traits

In the classic view introduced by R. A. Fisher, a quantitative trait is encoded by many loci with small, additive effects. Recent advances in quantitative trait loci mapping have begun to elucidate the genetic architectures underlying vast numbers of phenotypes across diverse taxa, producing observations that sometimes contrast with Fisher''s blueprint. Despite these considerable empirical efforts to map the genetic determinants of traits, it remains poorly understood how the genetic architecture of a trait should evolve, or how it depends on the selection pressures on the trait. Here, we develop a simple, population-genetic model for the evolution of genetic architectures. Our model predicts that traits under moderate selection should be encoded by many loci with highly variable effects, whereas traits under either weak or strong selection should be encoded by relatively few loci. We compare these theoretical predictions with qualitative trends in the genetics of human traits, and with systematic data on the genetics of gene expression levels in yeast. Our analysis provides an evolutionary explanation for broad empirical patterns in the genetic basis for traits, and it introduces a single framework that unifies the diversity of observed genetic architectures, ranging from Mendelian to Fisherian.     

Parental Analysis of Introgressive Hybridization between African and European Honeybees Using Nuclear DNA Rflps

African honeybees, introduced into Brazil 33 years ago, have spread through most of South and Central America and have largely replaced the extant European bees. Due to a paucity of genetic markers, genetic interactions between European and African bees are not well understood. Three restriction fragment length polymorphisms (RFLPs), detected with random, nuclear DNA probes, are described. The polymorphisms are specific to bees of European descent, possibly specific to certain European races. Each European marker was found present at a high frequency in U.S. colonies but absent in South African bees. Previous mitochondrial DNA studies of neotropical bees have revealed negligible maternal gene flow from managed European apiaries into feral African populations. The findings reported here with nuclear DNA show paternal gene flow between the two but suggest asymmetries in levels of introgressive hybridization. Managed colonies in southern Mexico, derived from European maternal lines, showed diminished levels of the European nuclear markers, reflecting significant hybridization with African drones. The European alleles were present only at low frequencies in feral swarms from the same area. The swarms were of African maternal descent. In Venezuelan colonies, also derived from African maternal lines, the European markers were almost totally absent. The results point to limited paternal introgression from European colonies into the African honeybee populations. These findings dispute other views regarding modes of Africanization.     

Multiethnic Genetic Association Studies Improve Power for Locus Discovery

To date, genome-wide association studies have focused almost exclusively on populations of European ancestry. These studies continue with the advent of next-generation sequencing, designed to systematically catalog and test low-frequency variation for a role in disease. A complementary approach would be to focus further efforts on cohorts of multiple ethnicities. This leverages the idea that population genetic drift may have elevated some variants to higher allele frequency in different populations, boosting statistical power to detect an association. Based on empirical allele frequency distributions from eleven populations represented in HapMap Phase 3 and the 1000 Genomes Project, we simulate a range of genetic models to quantify the power of association studies in multiple ethnicities relative to studies that exclusively focus on samples of European ancestry. In each of these simulations, a first phase of GWAS in exclusively European samples is followed by a second GWAS phase in any of the other populations (including a multiethnic design). We find that nontrivial power gains can be achieved by conducting future whole-genome studies in worldwide populations, where, in particular, African populations contribute the largest relative power gains for low-frequency alleles (<5%) of moderate effect that suffer from low power in samples of European descent. Our results emphasize the importance of broadening genetic studies to worldwide populations to ensure efficient discovery of genetic loci contributing to phenotypic trait variability, especially for those traits for which large numbers of samples of European ancestry have already been collected and tested.