Shared genetic control of expression and methylation in peripheral blood.

BACKGROUND: Expression QTLs and epigenetic marks are often employed to provide an insight into the possible biological mechanisms behind GWAS hits. A substantial proportion of the variation in gene expression and DNA methylation is known to be under genetic control. We address the proportion of genetic control that is shared between these two genomic features. RESULTS: An exhaustive search for pairwise phenotypic correlations between gene expression and DNA methylation in samples from human blood (n = 610) was performed. Of the 5 × 10(9) possible pairwise tests, 0.36 % passed Bonferroni corrected p-value cutoff of 9.9 × 10(-12). We determined that the correlation structure between probe pairs was largely due to blood cell type specificity of the expression and methylation probes. Upon adjustment of the expression and methylation values for observed blood cellular composition (n = 422), the number of probe pairs which survived Bonferroni correction reduced by more than 5400 fold. Of the 614 correlated probe pairs located on the same chromosome, 75 % share at least one methylation and expression QTL at nominal 10(-5) p-value cutoff. Those probe pairs are located within 1Mbp window from each other and have a mean of absolute value of genetic correlation equal to 0.69, further demonstrating the high degree of shared genetic control. CONCLUSIONS: Overall, this study demonstrates notable genetic covariance between DNA methylation and gene expression and reaffirms the importance of correcting for cell-counts in studies on non-homogeneous tissues.

Endometriosis risk alleles at 1p36.12 act through inverse regulation of CDC42 and LINC00339.

Genome-wide association studies (GWAS) have identified markers within the WNT4 region on chromosome 1p36.12 showing consistent and strong association with increasing endometriosis risk. Fine mapping using sequence and imputed genotype data has revealed strong candidates for the causal SNPs within these critical regions; however, the molecular pathogenesis of these SNPs is currently unknown. We used gene expression data collected from whole blood from 862 individuals and endometrial tissue from 136 individuals from independent populations of European descent to examine the mechanism underlying endometriosis susceptibility. Association mapping results from 7,090 individuals (2,594 cases and 4,496 controls) supported rs3820282 as the SNP with the strongest association for endometriosis risk (P = 1.84 × 10−5, OR = 1.244 (1.126-1.375)). SNP rs3820282 is a significant eQTL in whole blood decreasing expression of LINC00339 (also known as HSPC157) and increasing expression of CDC42 (P = 2.0 ×10−54 and 4.5x10−4 respectively). The largest effects were for two LINC00339 probes (P = 2.0 ×10−54; 1.0 × 10−34). The eQTL for LINC00339 was also observed in endometrial tissue (P = 2.4 ×10−8) with the same direction of effect for both whole blood and endometrial tissue. There was no evidence for eQTL effects for WNT4. Chromatin conformation capture provides evidence for risk SNPs interacting with the promoters of both LINC00339 and CDC4 and luciferase reporter assays suggest the risk SNP rs12038474 is located in a transcriptional silencer for CDC42 and the risk allele increases expression of CDC42. However, no effect of rs3820282 was observed in the LINC00339 expression in Ishikawa cells. Taken together, our results suggest that SNPs increasing endometriosis risk in this region act through CDC42, but further functional studies are required to rule out inverse regulation of both LINC00339 and CDC42.

Population genetic differentiation of height and body mass index across Europe.

Across-nation differences in the mean values for complex traits are common, but the reasons for these differences are unknown. Here we find that many independent loci contribute to population genetic differences in height and body mass index (BMI) in 9,416 individuals across 14 European countries. Using discovery data on over 250,000 individuals and unbiased effect size estimates from 17,500 sibling pairs, we estimate that 24% (95% credible interval (CI) = 9%, 41%) and 8% (95% CI = 4%, 16%) of the captured additive genetic variance for height and BMI, respectively, reflect population genetic differences. Population genetic divergence differed significantly from that in a null model (height, P < 3.94 × 10(-8); BMI, P < 5.95 × 10(-4)), and we find an among-population genetic correlation for tall and slender individuals (r = -0.80, 95% CI = -0.95, -0.60), consistent with correlated selection for both phenotypes. Observed differences in height among populations reflected the predicted genetic means (r = 0.51; P < 0.001), but environmental differences across Europe masked genetic differentiation for BMI (P < 0.58).

Seasonal effects on gene expression.

Many health conditions, ranging from psychiatric disorders to cardiovascular disease, display notable seasonal variation in severity and onset. In order to understand the molecular processes underlying this phenomenon, we have examined seasonal variation in the transcriptome of 606 healthy individuals. We show that 74 transcripts associated with a 12-month seasonal cycle were enriched for processes involved in DNA repair and binding. An additional 94 transcripts demonstrated significant seasonal variability that was largely influenced by blood cell count levels. These transcripts were enriched for immune function, protein production, and specific cellular markers for lymphocytes. Accordingly, cell counts for erythrocytes, platelets, neutrophils, monocytes, and CD19 cells demonstrated significant association with a 12-month seasonal cycle. These results demonstrate that seasonal variation is an important environmental regulator of gene expression and blood cell composition. Notable changes in leukocyte counts and genes involved in immune function indicate that immune cell physiology varies throughout the year in healthy individuals.

Interaction of c-Myb with p300 is required for the induction of acute myeloid leukemia (AML) by human AML oncogenes.

The MYB oncogene is widely expressed in acute leukemias and is important for the continued proliferation of leukemia cells, suggesting that MYB may be a therapeutic target in these diseases. However, realization of this potential requires a significant therapeutic window for MYB inhibition, given its essential role in normal hematopoiesis, and an approach for developing an effective therapeutic. We previously showed that the interaction of c-Myb with the coactivator CBP/p300 is essential for its transforming activity. Here, by using cells from Booreana mice which carry a mutant allele of c-Myb, we show that this interaction is essential for in vitro transformation by the myeloid leukemia oncogenes AML1-ETO, AML1-ETO9a, MLL-ENL, and MLL-AF9. We further show that unlike cells from wild-type mice, Booreana cells transduced with AML1-ETO9a or MLL-AF9 retroviruses fail to generate leukemia upon transplantation into irradiated recipients. Finally, we have begun to explore the molecular mechanisms underlying these observations by gene expression profiling. This identified several genes previously implicated in myeloid leukemogenesis and HSC function as being regulated in a c-Myb-p300-dependent manner. These data highlight the importance of the c-Myb-p300 interaction in myeloid leukemogenesis and suggest disruption of this interaction as a potential therapeutic strategy for acute myeloid leukemia.

Detection and replication of epistasis influencing transcription in humans.

Epistasis is the phenomenon whereby one polymorphism's effect on a trait depends on other polymorphisms present in the genome. The extent to which epistasis influences complex traits and contributes to their variation is a fundamental question in evolution and human genetics. Although often demonstrated in artificial gene manipulation studies in model organisms, and some examples have been reported in other species, few examples exist for epistasis among natural polymorphisms in human traits. Its absence from empirical findings may simply be due to low incidence in the genetic control of complex traits, but an alternative view is that it has previously been too technically challenging to detect owing to statistical and computational issues. Here we show, using advanced computation and a gene expression study design, that many instances of epistasis are found between common single nucleotide polymorphisms (SNPs). In a cohort of 846 individuals with 7,339 gene expression levels measured in peripheral blood, we found 501 significant pairwise interactions between common SNPs influencing the expression of 238 genes (P < 2.91 × 10(-16)). Replication of these interactions in two independent data sets showed both concordance of direction of epistatic effects (P = 5.56 × 10(-31)) and enrichment of interaction P values, with 30 being significant at a conservative threshold of P < 9.98 × 10(-5). Forty-four of the genetic interactions are located within 5 megabases of regions of known physical chromosome interactions (P = 1.8 × 10(-10)). Epistatic networks of three SNPs or more influence the expression levels of 129 genes, whereby one cis-acting SNP is modulated by several trans-acting SNPs. For example, MBNL1 is influenced by an additive effect at rs13069559, which itself is masked by trans-SNPs on 14 different chromosomes, with nearly identical genotype-phenotype maps for each cis-trans interaction. This study presents the first evidence, to our knowledge, for many instances of segregating common polymorphisms interacting to influence human traits.

Fine mapping of variants associated with endometriosis in the WNT4 region on chromosome 1p36.

Genome-wide association studies show strong evidence of association with endometriosis for markers on chromosome 1p36 spanning the potential candidate genes WNT4, CDC42 and LINC00339. WNT4 is involved in development of the uterus, and the expression of CDC42 and LINC00339 are altered in women with endometriosis. We conducted fine mapping to examine the role of coding variants in WNT4 and CDC42 and determine the key SNPs with strongest evidence of association in this region. We identified rare coding variants in WNT4 and CDC42 present only in endometriosis cases. The frequencies were low and cannot account for the common signal associated with increased risk of endometriosis. Genotypes for five common SNPs in the region of chromosome 1p36 show stronger association signals when compared with rs7521902 reported in published genome scans. Of these, three SNPs rs12404660, rs3820282, and rs55938609 were located in DNA sequences with potential functional roles including overlap with transcription factor binding sites for FOXA1, FOXA2, ESR1, and ESR2. Functional studies will be required to identify the gene or genes implicated in endometriosis risk.

Calculating statistical power in Mendelian randomization studies.

In Mendelian randomization (MR) studies, where genetic variants are used as proxy measures for an exposure trait of interest, obtaining adequate statistical power is frequently a concern due to the small amount of variation in a phenotypic trait that is typically explained by genetic variants. A range of power estimates based on simulations and specific parameters for two-stage least squares (2SLS) MR analyses based on continuous variables has previously been published. However there are presently no specific equations or software tools one can implement for calculating power of a given MR study. Using asymptotic theory, we show that in the case of continuous variables and a single instrument, for example a single-nucleotide polymorphism (SNP) or multiple SNP predictor, statistical power for a fixed sample size is a function of two parameters: the proportion of variation in the exposure variable explained by the genetic predictor and the true causal association between the exposure and outcome variable. We demonstrate that power for 2SLS MR can be derived using the non-centrality parameter (NCP) of the statistical test that is employed to test whether the 2SLS regression coefficient is zero. We show that the previously published power estimates from simulations can be represented theoretically using this NCP-based approach, with similar estimates observed when the simulation-based estimates are compared with our NCP-based approach. General equations for calculating statistical power for 2SLS MR using the NCP are provided in this note, and we implement the calculations in a web-based application.

GWAS of 126,559 individuals identifies genetic variants associated with educational attainment.

A genome-wide association study (GWAS) of educational attainment was conducted in a discovery sample of 101,069 individuals and a replication sample of 25,490. Three independent single-nucleotide polymorphisms (SNPs) are genome-wide significant (rs9320913, rs11584700, rs4851266), and all three replicate. Estimated effects sizes are small (coefficient of determination R(2) ≈ 0.02%), approximately 1 month of schooling per allele. A linear polygenic score from all measured SNPs accounts for ≈2% of the variance in both educational attainment and cognitive function. Genes in the region of the loci have previously been associated with health, cognitive, and central nervous system phenotypes, and bioinformatics analyses suggest the involvement of the anterior caudate nucleus. These findings provide promising candidate SNPs for follow-up work, and our effect size estimates can anchor power analyses in social-science genetics.

Genetically encoded fluorescent indicator for intracellular hydrogen peroxide.

We developed a genetically encoded, highly specific fluorescent probe for detecting hydrogen peroxide (H(2)O(2)) inside living cells. This probe, named HyPer, consists of circularly permuted yellow fluorescent protein (cpYFP) inserted into the regulatory domain of the prokaryotic H(2)O(2)-sensing protein, OxyR. Using HyPer we monitored H(2)O(2) production at the single-cell level in the cytoplasm and mitochondria of HeLa cells treated with Apo2L/TRAIL. We found that an increase in H(2)O(2) occurs in the cytoplasm in parallel with a drop in the mitochondrial transmembrane potential (DeltaPsi) and a change in cell shape. We also observed local bursts in mitochondrial H(2)O(2) production during DeltaPsi oscillations in apoptotic HeLa cells. Moreover, sensitivity of the probe was sufficient to observe H(2)O(2) increase upon physiological stimulation. Using HyPer we detected temporal increase in H(2)O(2) in the cytoplasm of PC-12 cells stimulated with nerve growth factor.

Human RFP2 gene promoter: unique structure and unusual strength.

Human gene RFP2 is a candidate tumor suppressor located at 13q14.3 and deleted in multiple tumor types. To explore regulation of RFP2, we determined structure of the 5'-untranslated region of RFP2 gene and its promoter. RFP2 promoter area is TATA-less, highly enriched in G and C nucleotides, and contains multiple quadruplex forming GGGGA-repeats. Deletion analysis of 5'-flanking sequences demonstrated that repeat containing fragment possesses activity seven times exceeding that of the combined SV40 promoter/enhancer. Other unusual features of the RFP2 promoter include anomalously high electrostatic fields induced by sequence-dependent dipoles and very low nucleosome forming potential. A "minimized" version of the RFP2 promoter could be used for overexpression of the various transgenes in the mammalian cells.