The Genotype-Tissue Expression (GTEx) project was established to characterize genetic effects on the transcriptome across human tissues and to link these regulatory mechanisms to trait and disease associations. Here, we present analyses of the version 8 data, examining 15,201 RNA-sequencing samples from 49 tissues of 838 postmortem donors. We comprehensively characterize genetic associations for gene expression and splicing in cis and trans, showing that regulatory associations are found for almost all genes, and describe the underlying molecular mechanisms and their contribution to allelic heterogeneity and pleiotropy of complex traits. Leveraging the large diversity of tissues, we provide insights into the tissue specificity of genetic effects and show that cell type composition is a key factor in understanding gene regulatory mechanisms in human tissues.

Alterations in DNA methylation frequently occur in hepatocellular cancer (HCC). We have previously demonstrated that hypermethylation in candidate genes can be detected in plasma DNA prior to HCC diagnosis. To identify with a genome-wide approach additional genes hypermethylated in HCC that could be used for more accurate analysis of plasma DNA for early diagnosis, we analyzed tumor and adjacent non-tumor tissues from 62 Taiwanese HCC cases using Illumina methylation arrays that screen 26,486 autosomal CpG sites. After Bonferroni adjustment, a total of 2,324 CpG sites significantly differed in methylation level, with 684 CpG sites significantly hypermethylated and 1,640 hypomethylated in tumor compared to non-tumor tissues. Array data were validated with pyrosequencing in a subset of 5 of these genes; correlation coefficients ranged from 0.92 to 0.97. Analysis of plasma DNA from 38 cases demonstrated that 37% to 63% of cases had detectable hypermethylated DNA (≥5% methylation) for these 5 genes individually. At least one of these genes was hypermethylated in 87% of cases, suggesting that measurement of DNA methylation in plasma samples is feasible. The panel of methylated genes indentified in the current study will be further tested in large cohort of prospectively collected samples to determine their utility as early biomarkers of hepatocellular carcinoma.

Arsenic contamination of drinking water is a major public health issue in many countries, increasing risk for a wide array of diseases, including cancer. There is inter-individual variation in arsenic metabolism efficiency and susceptibility to arsenic toxicity; however, the basis of this variation is not well understood. Here, we have performed the first genome-wide association study (GWAS) of arsenic-related metabolism and toxicity phenotypes to improve our understanding of the mechanisms by which arsenic affects health. Using data on urinary arsenic metabolite concentrations and approximately 300,000 genome-wide single nucleotide polymorphisms (SNPs) for 1,313 arsenic-exposed Bangladeshi individuals, we identified genome-wide significant association signals (P<5×10−8) for percentages of both monomethylarsonic acid (MMA) and dimethylarsinic acid (DMA) near the AS3MT gene (arsenite methyltransferase; 10q24.32), with five genetic variants showing independent associations. In a follow-up analysis of 1,085 individuals with arsenic-induced premalignant skin lesions (the classical sign of arsenic toxicity) and 1,794 controls, we show that one of these five variants (rs9527) is also associated with skin lesion risk (P = 0.0005). Using a subset of individuals with prospectively measured arsenic (n = 769), we show that rs9527 interacts with arsenic to influence incident skin lesion risk (P = 0.01). Expression quantitative trait locus (eQTL) analyses of genome-wide expression data from 950 individual's lymphocyte RNA suggest that several of our lead SNPs represent cis-eQTLs for AS3MT (P = 10−12) and neighboring gene C10orf32 (P = 10−44), which are involved in C10orf32-AS3MT read-through transcription. This is the largest and most comprehensive genomic investigation of arsenic metabolism and toxicity to date, the only GWAS of any arsenic-related trait, and the first study to implicate 10q24.32 variants in both arsenic metabolism and arsenical skin lesion risk. The observed patterns of associations suggest that MMA% and DMA% have distinct genetic determinants and support the hypothesis that DMA is the less toxic of these two methylated arsenic species. These results have potential translational implications for the prevention and treatment of arsenic-associated toxicities worldwide.

Telomere shortening is a hallmark of aging. Telomere length (TL) in blood cells has been studied extensively as a biomarker of human aging and disease; however, little is known regarding variability in TL in non-blood, disease-relevant tissue types. Here we characterize variability in TL measurements for 6,391 tissue samples, representing >20 tissue types and 952 individuals from the Genotype-Tissue Expression (GTEx) Project. We describe differences across tissue types, positive correlation among tissue types, and associations with age and ancestry. We show that genetic variation impacts TL in multiple tissue types, and that TL can mediate the effect of age on gene expression. Our results provide the foundational knowledge regarding TL in healthy tissues that is needed to interpret epidemiological studies of TL and human health. ONE SENTENCE SUMMARYTelomere length varies by tissue type but is generally correlated among tissue types (positively) and with age (negatively). MAIN TEXTTelomeres are DNA-protein complexes located at the end of chromosomes that protect chromosome ends from degradation and fusion (1). The length of the DNA component of telomeres shortens as cells divide (2) with short telomeres eventually triggering cellular senescence (3,4). In most human tissues, TL gradually shortens over the life course, and TL shortening is considered a hallmark (and a potential underlying cause) of human aging (5). In human studies, short TL measured in leukocytes is associated with increased risk of aging-related diseases including cardiovascular disease (6) and type II diabetes (7) as well as all-cause mortality (8). However, long TL may increase risk for some types of cancer (9-11). Leukocyte TL is influenced by inherited genetic variation (single nucleotide polymorphisms [SNPs]), some of which reside near genes with roles in telomere maintenance (12)(13)(14)(15). Leukocyte TL is also associated with lifestyle factors (e.g., obesity) and exposures (e.g., cigarette smoking) (16,17).Epidemiologic studies of TL predominantly use blood (occasionally saliva) as a DNA source. Thus, our understanding of variation in TL, its determinants (e.g., demographic, lifestyle, and genetic factors), and its associations with disease phenotypes is based almost entirely on TL measured in leukocytes from whole blood. Few prior studies have compared TL in leukocytes to TL in other human tissue types; these prior studies are relatively small (<100 participants; <5 tissue types) but provide evidence that TL differs across tissue types and that TL measurements from different tissue types are correlated (18,19). However, larger studies of many additional tissue types are needed to gain a comprehensive understanding of variation in TL and its determinants within and across a wide range of human tissues and cell types. In order to address these gaps in our understanding of TL and its role as a biomarker of aging and disease risk, we measured TL in > 6,000 unique tissue samples, representing >20 distinct tissue types and > 950 individual don...

We conducted a case-control study to investigate interindividual variability in susceptibility to health effects of inorganic arsenic due to arsenic metabolism efficiency, genetic factors, and their interaction. A total of 594 cases of arsenicinduced skin lesions and 1,041 controls was selected from baseline participants in a large prospective cohort study in Bangladesh. Adjusted odds ratios (OR) for skin lesions were estimated in relation to the polymorphisms in the glutathione S-transferase w1 and methylenetetrahydrofolate reductase genes, the percentage of monomethylarsonous acid (%MMA) and dimethylarsinic acid (%DMA) in urine, and the ratios of MMA to inorganic arsenic and DMA to MMA. Water arsenic concentration was positively associated with %MMA and inversely associated with %DMA. The doseresponse relationship of risk of skin lesion with %MMA was more apparent than those with other methylation indices; the ORs for skin lesions in relation to increasing %MMA quartiles were 1.00 (reference), 1.33 [95% confidence interval (95% CI), 0.92-1.93], 1.68 (95% CI, 1.17-2.42), and 1.57 (95% CI, 1.10-2.26; P for trend = 0.01). The ORs for skin lesions in relation to the methylenetetrahydrofolate reductase 677TT/ 1298AA and 677CT/1298AA diplotypes (compared with 677CC/1298CC diplotype) were 1.66 (95% CI, 1.00-2.77) and 1.77 (95% CI, 0.61-5.14), respectively. The OR for skin lesions in relation to the glutathione S-transferase w1 diplotype containing all at-risk alleles was 3.91 (95% CI, 1.03-14.79). Analysis of joint effects of genotypes/diplotypes with water arsenic concentration and urinary %MMA suggests additivity of these factors. The findings suggest that arsenic metabolism, particularly the conversion of MMA to DMA, may be saturable and that differences in urinary arsenic metabolites, genetic factors related to arsenic metabolism, and their joint distributions modulate arsenic toxicity. (Cancer Epidemiol Biomarkers Prev 2007;16(6):1270 -8)

Telomere shortening is a hallmark of aging. Telomere length (TL) in blood cells has been studied extensively as a biomarker of human aging and disease; however, little is known regarding variability in TL in nonblood, disease-relevant tissue types. Here, we characterize variability in TLs from 6391 tissue samples, representing >20 tissue types and 952 individuals from the Genotype-Tissue Expression (GTEx) project. We describe differences across tissue types, positive correlation among tissue types, and associations with age and ancestry. We show that genetic variation affects TL in multiple tissue types and that TL may mediate the effect of age on gene expression. Our results provide the foundational knowledge regarding TL in healthy tissues that is needed to interpret epidemiological studies of TL and human health.

No abstract

A large fraction of human genes are regulated by genetic variation near the transcribed sequence (cis-eQTL, expression quantitative trait locus), and many cis-eQTLs have implications for human disease. Less is known regarding the effects of genetic variation on expression of distant genes (trans-eQTLs) and their biological mechanisms. In this work, we use genome-wide data on SNPs and array-based expression measures from mononuclear cells obtained from a population-based cohort of 1,799 Bangladeshi individuals to characterize cis- and trans-eQTLs and determine if observed trans-eQTL associations are mediated by expression of transcripts in cis with the SNPs showing trans-association, using Sobel tests of mediation. We observed 434 independent trans-eQTL associations at a false-discovery rate of 0.05, and 189 of these trans-eQTLs were also cis-eQTLs (enrichment P<0.0001). Among these 189 trans-eQTL associations, 39 were significantly attenuated after adjusting for a cis-mediator based on Sobel P<10-5. We attempted to replicate 21 of these mediation signals in two European cohorts, and while only 7 trans-eQTL associations were present in one or both cohorts, 6 showed evidence of cis-mediation. Analyses of simulated data show that complete mediation will be observed as partial mediation in the presence of mediator measurement error or imperfect LD between measured and causal variants. Our data demonstrates that trans-associations can become significantly stronger or switch directions after adjusting for a potential mediator. Using simulated data, we demonstrate that this phenomenon is expected in the presence of strong cis-trans confounding and when the measured cis-transcript is correlated with the true (unmeasured) mediator. In conclusion, by applying mediation analysis to eQTL data, we show that a substantial fraction of observed trans-eQTL associations can be explained by cis-mediation. Future studies should focus on understanding the mechanisms underlying widespread cis-mediation and their relevance to disease biology, as well as using mediation analysis to improve eQTL discovery.