Novel approaches to the analysis of family data in genetic epidemiology

Association between genes regulating neural pathways for quantitative traits of speech and language disorders

Speech sound disorders (SSD) manifest as difficulties in phonological memory and awareness, oral motor function, language, vocabulary, reading, and spelling. Families enriched for SSD are rare, and typically display a cluster of deficits. We conducted a genome-wide association study (GWAS) in 435 children from 148 families in the Cleveland Family Speech and Reading study (CFSRS), examining 16 variables representing 6 domains. Replication was conducted using the Avon Longitudinal Study of Parents and Children (ALSPAC). We identified 18 significant loci (combined p < 10-8) that we pursued bioinformatically. We prioritized 5 novel gene regions with likely functional repercussions on neural pathways, including those which colocalized with differentially methylated regions in our sample. Polygenic risk scores for receptive language, expressive vocabulary, phonological awareness, phonological memory, spelling, and reading decoding associated with increasing clinical severity. In summary, neural-genetic influence on SSD is primarily multigenic and acts on genomic regulatory elements, similar to other neurodevelopmental disorders.

Genome-wide association analysis of metabolic syndrome quantitative traits in the GENNID multiethnic family study

BACKGROUND

To identify genetic associations of quantitative metabolic syndrome (MetS) traits and characterize heterogeneity across ethnic groups.

METHODS

Data was collected from GENetics of Noninsulin dependent Diabetes Mellitus (GENNID), a multiethnic resource of Type 2 diabetic families and included 1520 subjects in 259 African-American, European-American, Japanese-Americans, and Mexican-American families. We focused on eight MetS traits: weight, waist circumference, systolic and diastolic blood pressure, high-density lipoprotein, triglycerides, fasting glucose, and insulin. Using genotyped and imputed data from Illumina's Multiethnic array, we conducted genome-wide association analyses with linear mixed models for all ethnicities, except for the smaller Japanese-American group, where we used additive genetic models with gene-dropping.

RESULTS

Findings included ethnic-specific genetic associations and heterogeneity across ethnicities. Most significant associations were outside our candidate linkage regions and were coincident within a gene or intergenic region, with two exceptions in European-American families: (a) within previously identified linkage region on chromosome 2, two significant GLI2-TFCP2L1 associations with weight, and (b) one chromosome 11 variant near CADM1-LINC00900 with pleiotropic blood pressure effects.

CONCLUSIONS

This multiethnic family study found genetic heterogeneity and coincident associations (with one case of pleiotropy), highlighting the importance of including diverse populations in genetic research and illustrating the complex genetic architecture underlying MetS.

DNA methylation profiling to explore colorectal tumor differences according to menopausal hormone therapy use in women

Aim: Use of menopausal hormone therapy (MHT) has been associated with a reduced risk for colorectal cancer, but mechanisms underlying this relationship are not well understood. In the colon, MHT appears to act through estrogen receptor β (ERβ) which may influence DNA methylation by binding to DNA. Using genome-wide methylation profiling data, we aimed to identify genes that may be differentially methylated according to MHT use. Materials & methods: DNA methylation was measured using Illumina HumanMethylation450k arrays in two independent tumor sample sets of colorectal cancer patients. Differential methylation was determined using R/limma. Results: In the discovery analysis, two CpG sites showed differential DNA methylation according to MHT use, both were not replicated. In stratified analyses, 342 CpG sites were associated with current MHT use only in ERβ-positive tumors. Conclusion: The suggestive findings of differential methylation according to current MHT use in ERβ-positive tumors warrant further investigation.

Genome-wide DNA methylation differences according to oestrogen receptor beta status in colorectal cancer

Involvement of sex hormones in colorectal cancer (CRC) development has been linked to oestrogen receptor β (ERβ). Expression of ERβ is found reduced in tumour tissue and inversely related to mortality. However, mechanisms are not well understood. Our study aimed to detect differentially methylated genes associated with ERβ expression, which could point to mechanisms by which ERβ could influence risk and prognosis of CRC. Epigenome-wide DNA methylation profiling was performed using Illumina HumanMethylation450k BeadChip arrays in two independent tumour sample sets of CRC patients recruited in 2003-2010 by the German DACHS study (discovery cohort n = 917, replication cohort n = 907). ERβ expression was measured using immunohistochemistry and scored as negative, moderate and high. Differentially methylated CpG sites and genomic regions were determined using limma in the R-package RnBeads. For the comparison of tumours with moderate/high ERβ versus negative expression, differentially methylated CpG sites were identified but not confirmed by replication. Comparing tumours of high with tumours of negative ERβ expression revealed 2,904 differentially methylated CpG sites of which 403 were replicated (FDR adjusted p-value<0.05). Replicated CpGs were annotated to genes such as CD36, HK1 or LRP5. A survival analysis indicates that 30 of the replicated CpGs are also associated with overall survival (FDR-adjusted p-value<0.05). The regional analysis identified 60 differentially methylated promotor regions. The epigenome-wide analysis identified both novel genes as well as genes already implicated in CRC. Follow-up mechanistic studies to better understand the regulatory role of ERβ could inform potential targets for improving treatment or prevention of CRC.

Nightshift work, chronotype, and genome-wide DNA methylation in blood

Molecular mechanisms underlying the negative health effects of shift work are poorly understood, which remains a barrier to developing intervention strategies to protect the long-term health of shift workers. We evaluated genome-wide differences in DNA methylation (measured in blood) between 111 actively employed female nightshift and 86 actively employed female dayshift workers from the Seattle metropolitan area. We also explored the effect of chronotype (i.e., measure of preference for activity earlier or later in the day) on DNA methylation among 110 of the female nightshift workers and an additional group of 131 male nightshift workers. Methylation data were generated using the Illumina Infinium HumanMethylation450 BeadChip (450K) Array. After applying the latest methylation data processing methods, we compared methylation levels at 361,210 CpG loci between the groups using linear regression models adjusted for potential confounders and applied the false-discovery rate (FDR) ≤ 0.05 to account for multiple comparisons. No statistically significant associations at the genome-wide level were observed with shift work or chronotype, though based on raw P values and absolute effect sizes, there were suggestive associations in genes that have been previously linked with cancer (e.g., BACH2, JRK, RPS6KA2) and type-2 diabetes (e.g., KCNQ1). Given that our study was underpowered to detect moderate effects, examining these suggestive results in well-powered independent studies or in pooled data sets may improve our understanding of the pathways underlying the negative health effects of shift work and the influence of personal factors such as chronotype. Such an approach may help identify potential interventions that can be used to protect the long-term health of shift workers.