Daniel Glass

Mapping cis- and trans-regulatory effects across multiple tissues in twins

Sequence-based variation in gene expression is a key driver of disease risk. Common variants regulating expression in cis have been mapped in many expression quantitative trait locus (eQTL) studies, typically in single tissues from unrelated individuals. Here, we present a comprehensive analysis of gene expression across multiple tissues conducted in a large set of mono- and dizygotic twins that allows systematic dissection of genetic (cis and trans) and non-genetic effects on gene expression. Using identity-by-descent estimates, we show that at least 40% of the total heritable cis effect on expression cannot be accounted for by common cis variants, a finding that reveals the contribution of low-frequency and rare regulatory variants with respect to both transcriptional regulation and complex trait susceptibility. We show that a substantial proportion of gene expression heritability is trans to the structural gene, and we identify several replicating trans variants that act predominantly in a tissue-restricted manner and may regulate the transcription of many genes.

Epigenome-Wide Scans Identify Differentially Methylated Regions for Age and Age-Related Phenotypes in a Healthy Ageing Population

Age-related changes in DNA methylation have been implicated in cellular senescence and longevity, yet the causes and functional consequences of these variants remain unclear. To elucidate the role of age-related epigenetic changes in healthy ageing and potential longevity, we tested for association between whole-blood DNA methylation patterns in 172 female twins aged 32 to 80 with age and age-related phenotypes. Twin-based DNA methylation levels at 26,690 CpG-sites showed evidence for mean genome-wide heritability of 18%, which was supported by the identification of 1,537 CpG-sites with methylation QTLs in cis at FDR 5%. We performed genome-wide analyses to discover differentially methylated regions (DMRs) for sixteen age-related phenotypes (ap-DMRs) and chronological age (a-DMRs). Epigenome-wide association scans (EWAS) identified age-related phenotype DMRs (ap-DMRs) associated with LDL (STAT5A), lung function (WT1), and maternal longevity (ARL4A, TBX20). In contrast, EWAS for chronological age identified hundreds of predominantly hyper-methylated age DMRs (490 a-DMRs at FDR 5%), of which only one (TBX20) was also associated with an age-related phenotype. Therefore, the majority of age-related changes in DNA methylation are not associated with phenotypic measures of healthy ageing in later life. We replicated a large proportion of a-DMRs in a sample of 44 younger adult MZ twins aged 20 to 61, suggesting that a-DMRs may initiate at an earlier age. We next explored potential genetic and environmental mechanisms underlying a-DMRs and ap-DMRs. Genome-wide overlap across cis-meQTLs, genotype-phenotype associations, and EWAS ap-DMRs identified CpG-sites that had cis-meQTLs with evidence for genotype–phenotype association, where the CpG-site was also an ap-DMR for the same phenotype. Monozygotic twin methylation difference analyses identified one potential environmentally-mediated ap-DMR associated with total cholesterol and LDL (CSMD1). Our results suggest that in a small set of genes DNA methylation may be a candidate mechanism of mediating not only environmental, but also genetic effects on age-related phenotypes.

The architecture of gene regulatory variation across multiple human tissues: the MuTHER study.

While there have been studies exploring regulatory variation in one or more tissues, the complexity of tissue-specificity in multiple primary tissues is not yet well understood. We explore in depth the role of cis-regulatory variation in three human tissues: lymphoblastoid cell lines (LCL), skin, and fat. The samples (156 LCL, 160 skin, 166 fat) were derived simultaneously from a subset of well-phenotyped healthy female twins of the MuTHER resource. We discover an abundance of cis-eQTLs in each tissue similar to previous estimates (858 or 4.7% of genes). In addition, we apply factor analysis (FA) to remove effects of latent variables, thus more than doubling the number of our discoveries (1,822 eQTL genes). The unique study design (Matched Co-Twin Analysis—MCTA) permits immediate replication of eQTLs using co-twins (93%–98%) and validation of the considerable gain in eQTL discovery after FA correction. We highlight the challenges of comparing eQTLs between tissues. After verifying previous significance threshold-based estimates of tissue-specificity, we show their limitations given their dependency on statistical power. We propose that continuous estimates of the proportion of tissue-shared signals and direct comparison of the magnitude of effect on the fold change in expression are essential properties that jointly provide a biologically realistic view of tissue-specificity. Under this framework we demonstrate that 30% of eQTLs are shared among the three tissues studied, while another 29% appear exclusively tissue-specific. However, even among the shared eQTLs, a substantial proportion (10%–20%) have significant differences in the magnitude of fold change between genotypic classes across tissues. Our results underline the need to account for the complexity of eQTL tissue-specificity in an effort to assess consequences of such variants for complex traits.

Gene expression changes with age in skin, adipose tissue, blood and brain

BackgroundPrevious studies have demonstrated that gene expression levels change with age. These changes are hypothesized to influence the aging rate of an individual. We analyzed gene expression changes with age in abdominal skin, subcutaneous adipose tissue and lymphoblastoid cell lines in 856 female twins in the age range of 39-85 years. Additionally, we investigated genotypic variants involved in genotype-by-age interactions to understand how the genomic regulation of gene expression alters with age.ResultsUsing a linear mixed model, differential expression with age was identified in 1,672 genes in skin and 188 genes in adipose tissue. Only two genes expressed in lymphoblastoid cell lines showed significant changes with age. Genes significantly regulated by age were compared with expression profiles in 10 brain regions from 100 postmortem brains aged 16 to 83 years. We identified only one age-related gene common to the three tissues. There were 12 genes that showed differential expression with age in both skin and brain tissue and three common to adipose and brain tissues.ConclusionsSkin showed the most age-related gene expression changes of all the tissues investigated, with many of the genes being previously implicated in fatty acid metabolism, mitochondrial activity, cancer and splicing. A significant proportion of age-related changes in gene expression appear to be tissue-specific with only a few genes sharing an age effect in expression across tissues. More research is needed to improve our understanding of the genetic influences on aging and the relationship with age-related diseases.

Male-pattern baldness susceptibility locus at 20p11

We conducted a genome-wide association study for androgenic alopecia in 1,125 men and identified a newly associated locus at chromosome 20p11.22, confirmed in three independent cohorts (n = 1,650; OR = 1.60, P = 1.1 × 10−14 for rs1160312). The one man in seven who harbors risk alleles at both 20p11.22 and AR (encoding the androgen receptor) has a sevenfold-increased odds of androgenic alopecia (OR = 7.12, P = 3.7 × 10−15).

Gene-gene and gene-environment interactions detected by transcriptome sequence analysis in twins

Understanding the genetic architecture of gene expression is an intermediate step in understanding the genetic architecture of complex diseases. RNA sequencing technologies have improved the quantification of gene expression and allow measurement of allele-specific expression (ASE). ASE is hypothesized to result from the direct effect of cis regulatory variants, but a proper estimation of the causes of ASE has not been performed thus far. In this study, we take advantage of a sample of twins to measure the relative contributions of genetic and environmental effects to ASE, and we find substantial effects from gene × gene (G×G) and gene × environment (G×E) interactions. We propose a model where ASE requires genetic variability in cis, a difference in the sequence of both alleles, but where the magnitude of the ASE effect depends on trans genetic and environmental factors that interact with the cis genetic variants.

No Correlation Between Childhood Maltreatment and Telomere Length

To the Editor: Telomeres are lengths of repetitive DNA that cap the ends of chromosomes. They protect the ends of the chromosome and shorten with each cell division. Short leukocyte telomere length has been related to a number of age-related diseases (1). In addition, shorter telomere length has been associated with environmental factors such as smoking and lack of exercise (2,3). In a recent issue of Biological Psychiatry, Tyrka et al. (4) published a report suggesting a link between maltreatment in childhood and telomere shortening in 31 subjects. Individuals who had suffered maltreatment had telomere length .70 ± .24 compared with 1.02 ± .52 in individuals who had not been abused. We were interested in this plausible theory and decided to investigate these findings in the larger Twins UK cohort. Twins UK is a registry of ∼11,000 British twins that has been shown to be representative of the general population. Twins have answered questionnaires on two occasions relating to maltreatment, be it sexual, physical, or emotional abuse. One questionnaire asked specifically about maltreatment in childhood (2005) and one referred to maltreatment at any time (2008). The specific questions asked were the following: “As a child, were you ever a victim of physical abuse?” “As a child, were you ever a victim of sexual abuse?” “Have you ever been physically abused?” “Have you ever been emotionally abused?” “Have you ever been sexually abused?” The possible answers were Yes, No, or no answer. We only considered subjects that gave a consistent answer to the physical abuse or the sexual abuse questions in both questionnaires. We felt it further strengthened our phenotype to have corroborative data at two time points. Mean leukocyte telomere length was derived from the mean of the terminal restriction fragment length by using the Southern blot method on DNA extracted from peripheral leukocytes, as described elsewhere (5). Each DNA sample was analyzed twice. If the difference between the duplicates was >5%, a third measurement was performed, and the mean of the two results <5% apart was taken. The coefficient of variation of the terminal restriction fragment length assay in the study was 1.5%. Telomere length was essentially normally distributed. There is no significant difference in telomere length between people who have been abused and those who have not been abused (Figure 1, Table 1). The result holds for parametric and nonparametric tests and including or excluding possible confounding factors (age, gender, smoking, and body mass index). Figure 1 A graph showing no significant difference in mean telomere length in those who had suffered physical or sexual abuse compared with those who had not suffered abuse (t test p value .61 and .10, respectively). LCL, . Table 1 Telomere Length in Individuals with a History of Abuse Compared with Those with No History of Abuse When we looked at individuals who had answered the maltreatment in childhood questionnaire alone, we had a much larger sample size. With 123 cases and 1751 control subjects, there remains no significant difference in the mean telomere length between those subjects who had suffered sexual abuse and those who had not. The mean telomere length in this group is identical in cases and control subjects. This result was unchanged when possible confounders (age, body mass index, smoking status, and gender) were included in the linear model. We have failed to replicate the findings of Tyrka et al. (4), despite a sample size several-fold larger. There are certain methodological differences that may account for this lack of replication. Our telomere measurements were performed with Southern blot as opposed to quantitative polymerase chain reaction used in the previous article. We have, however, previously demonstrated significant changes using these assays with other phenotypes, such as lack of exercise, obesity, smoking, social class, and number of nevi. The assays are significantly correlated and the Southern blot method, although labor intensive, is thought to be the gold standard (6). The questionnaires to assess abuse used by Tyrka et al. (4) and Twins UK are different. We have confidence in our phenotype; by using twins who have reported maltreatment in both questionnaires, we have isolated a robust group of cases for analysis. Many different environmental factors are known to reduce telomeres, such as cigarette smoking, obesity, lack of exercise, and social class, as well as genes (2,3,7). With a small selected sample such as that used by Tyrka et al. (4), rather than the abuse, any of these confounding external factors could actually be responsible for the difference in telomere length they observed or it could have been due to chance.

Pigmentation and vitamin D metabolism in Caucasians: low vitamin D serum levels in fair skin types in the UK.

Background Vitamin D may play a protective role in many diseases. Public health messages are advocating sun avoidance to reduce skin cancer risk but the potential deleterious effects of these recommendations for vitamin D metabolism have been poorly investigated. Methodology/Principal Findings We investigated the association between 25-hydroxy-vitamin D (25(OH)D), skin type and ultraviolet exposure in 1414 Caucasian females in the UK. Mean age of the cohort was 47 years (18–79) and mean 25(OH)D levels were 77 nmol/L (6–289). 25(OH)D levels were strongly associated with season of sampling with higher levels in the spring and summer months (p<0.0001). Light skin types (skin type 1 and 2) have lower levels of 25(OH)D (mean 71 nmol/L) compared to darker skin types (skin type 3 and 4) (mean 82 nmol/L) after adjusting for multiple confounders (p<0.0001). The trend for increasing risk of low vitamin D with fairer skin types was highly significant despite adjustment for all confounders (p = 0.001). Conclusions/Significance Contrary to previous studies across different ethnic backgrounds, this study within Caucasian UK females shows that fair skin types have lower levels of 25(OH)D compared to darker skin types with potential detrimental health effects. Public health campaigns advocating sun avoidance in fair skinned individuals may need to be revised in view of their risk of vitamin D deficiency.

Eight common genetic variants associated with serum DHEAS levels suggest a key role in ageing mechanisms.

Dehydroepiandrosterone sulphate (DHEAS) is the most abundant circulating steroid secreted by adrenal glands—yet its function is unknown. Its serum concentration declines significantly with increasing age, which has led to speculation that a relative DHEAS deficiency may contribute to the development of common age-related diseases or diminished longevity. We conducted a meta-analysis of genome-wide association data with 14,846 individuals and identified eight independent common SNPs associated with serum DHEAS concentrations. Genes at or near the identified loci include ZKSCAN5 (rs11761528; p = 3.15×10−36), SULT2A1 (rs2637125; p = 2.61×10−19), ARPC1A (rs740160; p = 1.56×10−16), TRIM4 (rs17277546; p = 4.50×10−11), BMF (rs7181230; p = 5.44×10−11), HHEX (rs2497306; p = 4.64×10−9), BCL2L11 (rs6738028; p = 1.72×10−8), and CYP2C9 (rs2185570; p = 2.29×10−8). These genes are associated with type 2 diabetes, lymphoma, actin filament assembly, drug and xenobiotic metabolism, and zinc finger proteins. Several SNPs were associated with changes in gene expression levels, and the related genes are connected to biological pathways linking DHEAS with ageing. This study provides much needed insight into the function of DHEAS.

Genetics of skin color variation in Europeans: genome‑wide association studies with functional follow‑up

AbstractIn the International Visible Trait Genetics (VisiGen) Consortium, we investigated the genetics of human skin color by combining a series of genome-wide association studies (GWAS) in a total of 17,262 Europeans with functional follow-up of discovered loci. Our GWAS provide the first genome-wide significant evidence for chromosome 20q11.22 harboring the ASIP gene being explicitly associated with skin color in Europeans. In addition, genomic loci at 5p13.2 (SLC45A2), 6p25.3 (IRF4), 15q13.1 (HERC2/OCA2), and 16q24.3 (MC1R) were confirmed to be involved in skin coloration in Europeans. In follow-up gene expression and regulation studies of 22 genes in 20q11.22, we highlighted two novel genes EIF2S2 and GSS, serving as competing functional candidates in this region and providing future research lines. A genetically inferred skin color score obtained from the 9 top-associated SNPs from 9 genes in 940 worldwide samples (HGDP-CEPH) showed a clear gradual pattern in Western Eurasians similar to the distribution of physical skin color, suggesting the used 9 SNPs as suitable markers for DNA prediction of skin color in Europeans and neighboring populations, relevant in future forensic and anthropological investigations.