Genetic variants for head size share genes and pathways with cancer

The size of the human head is highly heritable, but genetic drivers of its variation within the general population remain unmapped. We perform a genome-wide association study on head size (N = 80,890) and identify 67 genetic loci, of which 50 are novel. Neuroimaging studies show that 17 variants affect specific brain areas, but most have widespread effects. Gene set enrichment is observed for various cancers and the p53, Wnt, and ErbB signaling pathways. Genes harboring lead variants are enriched for macrocephaly syndrome genes (37-fold) and high-fidelity cancer genes (9-fold), which is not seen for human height variants. Head size variants are also near genes preferentially expressed in intermediate progenitor cells, neural cells linked to evolutionary brain expansion. Our results indicate that genes regulating early brain and cranial growth incline to neoplasia later in life, irrespective of height. This warrants investigation of clinical implications of the link between head size and cancer.

Antidepressant Exposure and DNA Methylation: Insights from a Methylome-Wide Association Study

Importance

Understanding antidepressant mechanisms could help design more effective and tolerated treatments.

Objective

Identify DNA methylation (DNAm) changes associated with antidepressant exposure.

Design

Case-control methylome-wide association studies (MWAS) of antidepressant exposure were performed from blood samples collected between 2006-2011 in Generation Scotland (GS). The summary statistics were tested for enrichment in specific tissues, gene ontologies and an independent MWAS in the Netherlands Study of Depression and Anxiety (NESDA). A methylation profile score (MPS) was derived and tested for its association with antidepressant exposure in eight independent cohorts, alongside prospective data from GS.

Setting

Cohorts; GS, NESDA, FTC, SHIP-Trend, FOR2107, LBC1936, MARS-UniDep, ALSPAC, E-Risk, and NTR.

Participants

Participants with DNAm data and self-report/prescription derived antidepressant exposure.

Main Outcomes and Measures

Whole-blood DNAm levels were assayed by the EPIC/450K Illumina array (9 studies, N exposed = 661, N unexposed = 9,575) alongside MBD-Seq in NESDA (N exposed = 398, N unexposed = 414). Antidepressant exposure was measured by self- report and/or antidepressant prescriptions.

Results

The self-report MWAS (N = 16,536, N exposed = 1,508, mean age = 48, 59% female) and the prescription-derived MWAS (N = 7,951, N exposed = 861, mean age = 47, 59% female), found hypermethylation at seven and four DNAm sites (p < 9.42x10 -8 ), respectively. The top locus was cg26277237 ( KANK1, p self-report = 9.3x10 -13 , p prescription = 6.1x10 -3 ). The self-report MWAS found a differentially methylated region, mapping to DGUOK-AS1 ( p adj = 5.0x10 -3 ) alongside significant enrichment for genes expressed in the amygdala, the "synaptic vesicle membrane" gene ontology and the top 1% of CpGs from the NESDA MWAS (OR = 1.39, p < 0.042). The MPS was associated with antidepressant exposure in meta-analysed data from external cohorts (N studies = 9, N = 10,236, N exposed = 661, f3 = 0.196, p < 1x10 -4 ).

Conclusions and Relevance

Antidepressant exposure is associated with changes in DNAm across different cohorts. Further investigation into these changes could inform on new targets for antidepressant treatments.

3 Key Points

Question: Is antidepressant exposure associated with differential whole blood DNA methylation?Findings: In this methylome-wide association study of 16,536 adults across Scotland, antidepressant exposure was significantly associated with hypermethylation at CpGs mapping to KANK1 and DGUOK-AS1. A methylation profile score trained on this sample was significantly associated with antidepressant exposure (pooled f3 [95%CI]=0.196 [0.105, 0.288], p < 1x10 -4 ) in a meta-analysis of external datasets. Meaning: Antidepressant exposure is associated with hypermethylation at KANK1 and DGUOK-AS1 , which have roles in mitochondrial metabolism and neurite outgrowth. If replicated in future studies, targeting these genes could inform the design of more effective and better tolerated treatments for depression.

Genome-wide analysis in over 1 million individuals of European ancestry yields improved polygenic risk scores for blood pressure traits

Hypertension affects more than one billion people worldwide. Here we identify 113 novel loci, reporting a total of 2,103 independent genetic signals (P < 5 × 10-8) from the largest single-stage blood pressure (BP) genome-wide association study to date (n = 1,028,980 European individuals). These associations explain more than 60% of single nucleotide polymorphism-based BP heritability. Comparing top versus bottom deciles of polygenic risk scores (PRSs) reveals clinically meaningful differences in BP (16.9 mmHg systolic BP, 95% CI, 15.5-18.2 mmHg, P = 2.22 × 10-126) and more than a sevenfold higher odds of hypertension risk (odds ratio, 7.33; 95% CI, 5.54-9.70; P = 4.13 × 10-44) in an independent dataset. Adding PRS into hypertension-prediction models increased the area under the receiver operating characteristic curve (AUROC) from 0.791 (95% CI, 0.781-0.801) to 0.826 (95% CI, 0.817-0.836, ∆AUROC, 0.035, P = 1.98 × 10-34). We compare the 2,103 loci results in non-European ancestries and show significant PRS associations in a large African-American sample. Secondary analyses implicate 500 genes previously unreported for BP. Our study highlights the role of increasingly large genomic studies for precision health research.

Association analysis between an epigenetic alcohol risk score and blood pressure

Background Epigenome-wide association studies have revealed multiple DNA methylation sites (CpGs) associated with alcohol consumption, an important lifestyle risk factor for cardiovascular diseases. Results We generated an alcohol consumption epigenetic risk score (ERS) based on previously reported 144 alcohol-associated CpGs and examined the association of the ERS with systolic blood pressure (SBP), diastolic blood pressure (DBP), and hypertension (HTN) in 3,898 Framingham Heart Study (FHS) participants. We found an association of alcohol intake with the ERS in the meta-analysis with 0.09 units higher ERS per drink consumed per day ( p  < 0.0001). Cross-sectional analyses in FHS revealed that a one-unit increment of the ERS was associated with 1.93 mm Hg higher SBP ( p  = 4.64E-07), 0.68 mm Hg higher DBP ( p  = 0.006), and an odds ratio of 1.78 for HTN ( p  < 2E-16). Meta-analysis of the cross-sectional association of the ERS with BP traits in eight independent external cohorts (n = 11,544) showed similar relationships with blood pressure levels, i.e., a one-unit increase in ERS was associated with 0.74 ( p  = 0.002) and 0.50 ( p  = 0.0006) mm Hg higher SBP and DBP, but could not confirm the association with hypertension. Longitudinal analyses in FHS (n = 3,260) and five independent external cohorts (n = 4,021) showed that the baseline ERS was not associated with a change in blood pressure over time or with incident HTN. Conclusions Our findings provide proof-of-concept that utilizing an ERS is a useful approach to capture the recent health consequences of lifestyle behaviors such as alcohol consumption.

Association analysis between an epigenetic alcohol risk score and blood pressure

Background

Epigenome-wide association studies have revealed multiple DNA methylation sites (CpGs) associated with alcohol consumption, an important lifestyle risk factor for cardiovascular diseases.

Results

We generated an alcohol consumption epigenetic risk score (ERS) based on previously reported 144 alcohol-associated CpGs and examined the association of the ERS with systolic blood pressure (SBP), diastolic blood pressure (DBP), and hypertension (HTN) in 3,898 Framingham Heart Study (FHS) participants. We found an association of alcohol intake with the ERS in the meta-analysis with 0.09 units higher ERS per drink consumed per day (p < 0.0001). Cross-sectional analyses in FHS revealed that a one-unit increment of the ERS was associated with 1.93 mm Hg higher SBP (p = 4.64E-07), 0.68 mm Hg higher DBP (p = 0.006), and an odds ratio of 1.78 for HTN (p < 2E-16). Meta-analysis of the cross-sectional association of the ERS with BP traits in eight independent external cohorts (n = 11,544) showed similar relationships with blood pressure levels, i.e., a one-unit increase in ERS was associated with 0.74 (p = 0.002) and 0.50 (p = 0.0006) mm Hg higher SBP and DBP, but could not confirm the association with hypertension. Longitudinal analyses in FHS (n = 3,260) and five independent external cohorts (n = 4,021) showed that the baseline ERS was not associated with a change in blood pressure over time or with incident HTN.

Conclusions

Our findings provide proof-of-concept that utilizing an ERS is a useful approach to capture the recent health consequences of lifestyle behaviors such as alcohol consumption.

The role of thyroid function in borderline personality disorder and schizophrenia: a Mendelian Randomisation study

BACKGROUND

Genome-wide association studies have reported a genetic overlap between borderline personality disorder (BPD) and schizophrenia (SCZ). Epidemiologically, the direction and causality of the association between thyroid function and risk of BPD and SCZ are unclear. We aim to test whether genetically predicted variations in TSH and FT4 levels or hypothyroidism are associated with the risk of BPD and SCZ.

METHODS

We employed Mendelian Randomisation (MR) analyses using genetic instruments associated with TSH and FT4 levels as well as hypothyroidism to examine the effects of genetically predicted thyroid function on BPD and SCZ risk. Bidirectional MR analyses were employed to investigate a potential reverse causal association.

RESULTS

Genetically predicted higher FT4 was not associated with the risk of BPD (OR: 1.18; P = 0.60, IVW) or the risk of SCZ (OR: 0.93; P = 0.19, IVW). Genetically predicted higher TSH was not associated with the risk of BPD (OR: 1.11; P = 0.51, IVW) or SCZ (OR: 0.98, P = 0.55, IVW). Genetically predicted hypothyroidism was not associated with BPD or SCZ. We found no evidence for a reverse causal effect between BPD or SCZ on thyroid function.

CONCLUSIONS

We report evidence for a null association between genetically predicted FT4, TSH or hypothyroidism with BPD or SCZ risk. There was no evidence for reverse causality.

Alpha-1 Antitrypsin Augmentation and the Liver Phenotype of Adults With Alpha-1 Antitrypsin Deficiency (Genotype Pi∗ZZ)

BACKGROUND & AIMS

α1-Antitrypsin (AAT) is a major protease inhibitor produced by hepatocytes. The most relevant AAT mutation giving rise to AAT deficiency (AATD), the 'Pi∗Z' variant, causes harmful AAT protein accumulation in the liver, shortage of AAT in the systemic circulation, and thereby predisposes to liver and lung injury. Although intravenous AAT augmentation constitutes an established treatment of AATD-associated lung disease, its impact on the liver is unknown.

METHODS

Liver-related parameters were assessed in a multinational cohort of 760 adults with severe AATD (Pi∗ZZ genotype) and available liver phenotyping, of whom 344 received augmentation therapy and 416 did not. Liver fibrosis was evaluated noninvasively via the serum test AST-to-platelet ratio index and via transient elastography-based liver stiffness measurement. Histologic parameters were compared in 15 Pi∗ZZ adults with and 35 without augmentation.

RESULTS

Compared with nonaugmented subjects, augmented Pi∗ZZ individuals displayed lower serum liver enzyme levels (AST 71% vs 75% upper limit of normal, P < .001; bilirubin 49% vs 58% upper limit of normal, P = .019) and lower surrogate markers of fibrosis (AST-to-platelet ratio index 0.34 vs 0.38, P < .001; liver stiffness measurement 6.5 vs 7.2 kPa, P = .005). Among biopsied participants, augmented individuals had less pronounced liver fibrosis and less inflammatory foci but no differences in AAT accumulation were noted.

CONCLUSIONS

The first evaluation of AAT augmentation on the Pi∗ZZ-related liver disease indicates liver safety of a widely used treatment for AATD-associated lung disease. Prospective studies are needed to confirm the beneficial effects and to demonstrate the potential efficacy of exogenous AAT in patients with Pi∗ZZ-associated liver disease.

Multi-trait analysis characterizes the genetics of thyroid function and identifies causal associations with clinical implications

To date only a fraction of the genetic footprint of thyroid function has been clarified. We report a genome-wide association study meta-analysis of thyroid function in up to 271,040 individuals of European ancestry, including reference range thyrotropin (TSH), free thyroxine (FT4), free and total triiodothyronine (T3), proxies for metabolism (T3/FT4 ratio) as well as dichotomized high and low TSH levels. We revealed 259 independent significant associations for TSH (61% novel), 85 for FT4 (67% novel), and 62 novel signals for the T3 related traits. The loci explained 14.1%, 6.0%, 9.5% and 1.1% of the total variation in TSH, FT4, total T3 and free T3 concentrations, respectively. Genetic correlations indicate that TSH associated loci reflect the thyroid function determined by free T3, whereas the FT4 associations represent the thyroid hormone metabolism. Polygenic risk score and Mendelian randomization analyses showed the effects of genetically determined variation in thyroid function on various clinical outcomes, including cardiovascular risk factors and diseases, autoimmune diseases, and cancer. In conclusion, our results improve the understanding of thyroid hormone physiology and highlight the pleiotropic effects of thyroid function on various diseases.

Methylation Patterns of the FKBP5 Gene in Association with Childhood Maltreatment and Depressive Disorders

Childhood maltreatment is an important risk factor for adult depression and has been associated with changes in the hypothalamic pituitary adrenal (HPA) axis, including cortisol secretion and methylation of the FKBP5 gene. Furthermore, associations between depression and HPA changes have been reported. This study investigated the associations of whole-blood FKBP5 mRNA levels, serum cortisol levels, childhood maltreatment, and depressive symptoms with the whole-blood methylation status (assessed via target bisulfite sequencing) of 105 CpGs at the FKBP5 locus using data from the general population-based Study of Health in Pomerania (SHIP) (N = 203). Both direct and interaction effects with the rs1360780 single-nucleotide polymorphism were investigated. Nominally significant associations of main effects on methylation of a single CpG site were observed at intron 3, intron 7, and the 3'-end of the gene. Additionally, methylation at two clusters at the 3'-end and intron 7 were nominally associated with childhood maltreatment × rs1360780 and depressive symptoms × rs1360780, respectively. The results add to the understanding of molecular mechanisms underlying the emergence of depression and could aid the development of personalised depression therapy and drug development.

Multi-omics and pathway analyses of genome-wide associations implicate regulation and immunity in verbal declarative memory performance

BACKGROUND

Uncovering the functional relevance underlying verbal declarative memory (VDM) genome-wide association study (GWAS) results may facilitate the development of interventions to reduce age-related memory decline and dementia.

METHODS

We performed multi-omics and pathway enrichment analyses of paragraph (PAR-dr) and word list (WL-dr) delayed recall GWAS from 29,076 older non-demented individuals of European descent. We assessed the relationship between single-variant associations and expression quantitative trait loci (eQTLs) in 44 tissues and methylation quantitative trait loci (meQTLs) in the hippocampus. We determined the relationship between gene associations and transcript levels in 53 tissues, annotation as immune genes, and regulation by transcription factors (TFs) and microRNAs. To identify significant pathways, gene set enrichment was tested in each cohort and meta-analyzed across cohorts. Analyses of differential expression in brain tissues were conducted for pathway component genes.

RESULTS

The single-variant associations of VDM showed significant linkage disequilibrium (LD) with eQTLs across all tissues and meQTLs within the hippocampus. Stronger WL-dr gene associations correlated with reduced expression in four brain tissues, including the hippocampus. More robust PAR-dr and/or WL-dr gene associations were intricately linked with immunity and were influenced by 31 TFs and 2 microRNAs. Six pathways, including type I diabetes, exhibited significant associations with both PAR-dr and WL-dr. These pathways included fifteen MHC genes intricately linked to VDM performance, showing diverse expression patterns based on cognitive status in brain tissues.

CONCLUSIONS

VDM genetic associations influence expression regulation via eQTLs and meQTLs. The involvement of TFs, microRNAs, MHC genes, and immune-related pathways contributes to VDM performance in older individuals.

X-chromosome and kidney function: evidence from a multi-trait genetic analysis of 908,697 individuals reveals sex-specific and sex-differential findings in genes regulated by androgen response elements

X-chromosomal genetic variants are understudied but can yield valuable insights into sexually dimorphic human traits and diseases. We performed a sex-stratified cross-ancestry X-chromosome-wide association meta-analysis of seven kidney-related traits (n = 908,697), identifying 23 loci genome-wide significantly associated with two of the traits: 7 for uric acid and 16 for estimated glomerular filtration rate (eGFR), including four novel eGFR loci containing the functionally plausible prioritized genes ACSL4, CLDN2, TSPAN6 and the female-specific DRP2. Further, we identified five novel sex-interactions, comprising male-specific effects at FAM9B and AR/EDA2R, and three sex-differential findings with larger genetic effect sizes in males at DCAF12L1 and MST4 and larger effect sizes in females at HPRT1. All prioritized genes in loci showing significant sex-interactions were located next to androgen response elements (ARE). Five ARE genes showed sex-differential expressions. This study contributes new insights into sex-dimorphisms of kidney traits along with new prioritized gene targets for further molecular research.

Beyond the Global Brain Differences: Intraindividual Variability Differences in 1q21.1 Distal and 15q11.2 BP1-BP2 Deletion Carriers

BACKGROUND

Carriers of the 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants exhibit regional and global brain differences compared with noncarriers. However, interpreting regional differences is challenging if a global difference drives the regional brain differences. Intraindividual variability measures can be used to test for regional differences beyond global differences in brain structure.

METHODS

Magnetic resonance imaging data were used to obtain regional brain values for 1q21.1 distal deletion (n = 30) and duplication (n = 27) and 15q11.2 BP1-BP2 deletion (n = 170) and duplication (n = 243) carriers and matched noncarriers (n = 2350). Regional intra-deviation scores, i.e., the standardized difference between an individual's regional difference and global difference, were used to test for regional differences that diverge from the global difference.

RESULTS

For the 1q21.1 distal deletion carriers, cortical surface area for regions in the medial visual cortex, posterior cingulate, and temporal pole differed less and regions in the prefrontal and superior temporal cortex differed more than the global difference in cortical surface area. For the 15q11.2 BP1-BP2 deletion carriers, cortical thickness in regions in the medial visual cortex, auditory cortex, and temporal pole differed less and the prefrontal and somatosensory cortex differed more than the global difference in cortical thickness.

CONCLUSIONS

We find evidence for regional effects beyond differences in global brain measures in 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants. The results provide new insight into brain profiling of the 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants, with the potential to increase understanding of the mechanisms involved in altered neurodevelopment.

Mendelian randomization indicates causal effects of estradiol levels on kidney function in males

Context

Chronic kidney disease (CKD) is a public health burden worldwide. Epidemiological studies observed an association between sex hormones, including estradiol, and kidney function.

Objective

We conducted a Mendelian randomization (MR) study to assess a possible causal effect of estradiol levels on kidney function in males and females.

Design

We performed a bidirectional two-sample MR using published genetic associations of serum levels of estradiol in men (n = 206,927) and women (n = 229,966), and of kidney traits represented by estimated glomerular filtration rate (eGFR, n = 567,460), urine albumin-to-creatinine ratio (UACR, n = 547,361), and CKD (n = 41,395 cases and n = 439,303 controls) using data obtained from the CKDGen Consortium. Additionally, we conducted a genome-wide association study using UK Biobank cohort study data (n = 11,798 men and n = 6,835 women) to identify novel genetic associations with levels of estradiol, and then used these variants as instruments in a one-sample MR.

Results

The two-sample MR indicated that genetically predicted estradiol levels are significantly associated with eGFR in men (beta = 0.077; p = 5.2E-05). We identified a single locus at chromosome 14 associated with estradiol levels in men being significant in the one-sample MR on eGFR (beta = 0.199; p = 0.017). We revealed significant results with eGFR in postmenopausal women and with UACR in premenopausal women, which did not reach statistical significance in the sensitivity MR analyses. No causal effect of eGFR or UACR on estradiol levels was found.

Conclusions

We conclude that serum estradiol levels may have a causal effect on kidney function. Our MR results provide starting points for studies to develop therapeutic strategies to reduce kidney disease.

Height, Autoimmune Thyroid Disease, and Thyroid Cancer: A Mendelian Randomization Study

Background: Increased height has been associated with increased risk of hypothyroidism or thyroid cancer in epidemiological studies. However, the potential causal association between height and hypothyroidism or thyroid cancer has not been thoroughly explored. Autoimmune thyroid disease (AITD) mainly presents as hypothyroidism, thus we aim to evaluate the causal relationship between height as exposure and its association with AITD or thyroid cancer. Methods: Mendelian randomization (MR) analyses were performed by using genetic instruments associated with height, which were selected from the largest genome-wide association meta-analysis for height in up to 5.4 million individuals. Summary-level data for AITD and thyroid cancer (including 30,234 and 3001 cases, respectively) were collected from the large number of available genome-wide association studies. Bidirectional MR was performed to test for reverse causal association between AITD and adult height. Results: MR analyses showed that increased genetically predicted height was associated with a 4% increased risk of AITD ([CI 1.02 to 1.07], p-value = 1.99E-03) per 1-standard deviation (SD) increase in genetically predicted height. The bidirectional MR did not show any causal association between AITD and adult height. Additionally, increased genetically predicted height was associated with 15% increased risk of thyroid cancer ([CI 1.07 to 1.23], p-value = 2.32E-04) per 1-SD increase in height. Sensitivity analysis confirmed the main results. Conclusions: This MR study showed that 1-SD increase in genetically predicted height was associated with increased risk of AITD and thyroid cancer. In contrast, there was no evidence of a causal association of genetically predicted AITD with height. These results could further aid in investigation of height-related pathways as a means of gaining new mechanistic insights into AITD and thyroid cancer.

Correction: Analysis of epidemiological association patterns of serum thyrotropin by combining random forests and Bayesian networks

[This corrects the article DOI: 10.1371/journal.pone.0271610.].

Genetic Variants in WNT16 and PKD2L1 Locus Affect Heel Ultrasound Bone Stiffness: Analyses from the General Population and Patients Evaluated for Osteoporosis

Osteoporosis, a complex chronic disease with increasing prevalence, is characterised by reduced bone mineral density (BMD) and increased fracture risk. The high heritability of BMD suggests substantial impact of the individual genetic disposition on bone phenotypes and the development of osteoporosis. In the past years, genome-wide association studies (GWAS) identified hundreds of genetic variants associated with BMD or osteoporosis. Here, we analysed 1103 single nucleotide polymorphisms (SNPs), previously identified as associated with estimated BMD (eBMD) in the UK Biobank. We assessed whether these SNPs are related to heel stiffness index obtained by quantitative ultrasound in 5665 adult participants of the Study of Health in Pomerania (SHIP). We confirmed 45 significant associations after correction for multiple testing. Next, we analysed six selected SNPs in 631 patients evaluated for osteoporosis [rs2707518 (CPED1/WNT16), rs3779381 (WNT16), rs115242848 (LOC101927709/EN1), rs10239787 (JAZF1), rs603424 (PKD2L1) and rs6968704 (JAZF1)]. Differences in minor allele frequencies (MAF) of rs2707518 and rs3779381 between SHIP participants (higher MAF) and patients evaluated for osteoporosis (lower MAF) indicated a protective effect of the minor allele on bone integrity. In contrast, differences in MAF of rs603424 indicated a harmful effect. Co-localisation analyses indicated that the rs603424 effect may be mediated via stearoyl-CoA desaturase (SCD) expression, an enzyme highly expressed in adipose tissue with a crucial role in lipogenesis. Taken together, our results support the role of the WNT16 pathway in the regulation of bone properties and indicate a novel causal role of SCD expression in adipose tissue on bone integrity.

High Thyroid-Stimulating Hormone and Low Free Triiodothyronine Levels Are Associated with Chronic Kidney Disease in Three Population-Based Studies from Germany

High serum thyroid-stimulating hormone (TSH) levels have previously been associated with a low estimated glomerular filtration rate (eGFR), but studies associating thyroid hormone levels with albuminuria revealed inconsistent results. We used cross-sectional data from 7933 individuals aged 20 to 93 years of the Berlin Aging Study II and the Study of Health in Pomerania to associate serum TSH, fT3, and fT4 levels with eGFR and albuminuria. In multivariable analyses adjusted for confounding, we found inverse non-linear associations of serum TSH levels with eGFR, while serum fT3 levels showed a positive association with eGFR. High as well as low serum fT4 levels were associated with a lower eGFR. Age but not sex modified the association between thyroid hormone levels and eGFR. The inverse associations between serum TSH levels and eGFR were strongest in the youngest age groups, while the positive associations between serum fT3 levels and eGFR were strongest in older individuals. No significant associations between thyroid hormone levels and albuminuria were found. Our results indicate that hypothyroidism might be associated with a reduced kidney function. Thyroid function might be more tightly related to the eGFR than to albuminuria in the general population.

Novel (sulfated) thyroid hormone transporters in the solute carrier 22 family

Objective

Thyroid hormone (TH) transport represents a critical first step in governing intracellular TH regulation. It is still unknown whether the full repertoire of TH transporters has been identified. Members of the solute carrier (SLC) 22 family have substrates in common with the known TH transporters of the organic anion-transporting peptide family. Therefore, we screened the SLC22 family for TH transporters.

Methods

Uptake of 1 nM of iodothyronines or sulfated iodothyronines in COS1 cells expressing SLC22 proteins was performed.

Results

We first tested 25 mouse (m) SLC22 proteins for TH uptake and found that the majority of the organic anion transporter (OAT) clade were capable of 3,3',5-triiodothyronine and/or thyroxine (T4) transport. Based on phylogenetic tree analysis of the mouse and human (h) SLC22 family, we selected eight hSLC22s that grouped with the newly identified mouse TH transporters. Of these, four tested positive for uptake of one or more substrates, particularly hSLC22A11 showed robust (3-fold over control) uptake of T4. Uptake of sulfated iodothyronines was strongly (up to 17-fold) induced by some SLC22s, most notably SLC22A8, hSLC22A9, mSLC22A27 and mSLC22A29. Finally, the zebrafish orthologues of SLC22A6/8 drOatx and drSlc22a6l also transported almost all (sulfated) iodothyronines tested. The OAT inhibitors lesinurad and probenecid inhibited most SLC22 proteins.

Conclusions

Our results demonstrated that members of the OAT clade of the SLC22 family constitute a novel, evolutionary conserved group of transporters for (sulfated) iodothyronines. Future studies should reveal the relevance of these transporters in TH homeostasis and physiology.

Distinct genetic liability profiles define clinically relevant patient strata across common diseases

Genome-wide association studies have unearthed a wealth of genetic associations across many complex diseases. However, translating these associations into biological mechanisms contributing to disease etiology and heterogeneity has been challenging. Here, we hypothesize that the effects of disease-associated genetic variants converge onto distinct cell type specific molecular pathways within distinct subgroups of patients. In order to test this hypothesis, we develop the CASTom-iGEx pipeline to operationalize individual level genotype data to interpret personal polygenic risk and identify the genetic basis of clinical heterogeneity. The paradigmatic application of this approach to coronary artery disease and schizophrenia reveals a convergence of disease associated variant effects onto known and novel genes, pathways, and biological processes. The biological process specific genetic liabilities are not equally distributed across patients. Instead, they defined genetically distinct groups of patients, characterized by different profiles across pathways, endophenotypes, and disease severity. These results provide further evidence for a genetic contribution to clinical heterogeneity and point to the existence of partially distinct pathomechanisms across patient subgroups. Thus, the universally applicable approach presented here has the potential to constitute an important component of future personalized medicine concepts.

DNA methylation analysis is used to identify novel genetic loci associated with circulating fibrinogen levels in blood

BACKGROUND

Fibrinogen plays an essential role in blood coagulation and inflammation. Circulating fibrinogen levels may be determined based on interindividual differences in DNA methylation at cytosine-phosphate-guanine (CpG) sites and vice versa.

OBJECTIVES

To perform an EWAS to examine an association between blood DNA methylation levels and circulating fibrinogen levels to better understand its biological and pathophysiological actions.

METHODS

We performed an epigenome-wide association study of circulating fibrinogen levels in 18 037 White, Black, American Indian, and Hispanic participants, representing 14 studies from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium. Circulating leukocyte DNA methylation was measured using the Illumina 450K array in 12 904 participants and using the EPIC array in 5133 participants. In each study, an epigenome-wide association study of fibrinogen was performed using linear mixed models adjusted for potential confounders. Study-specific results were combined using array-specific meta-analysis, followed by cross-replication of epigenome-wide significant associations. We compared models with and without CRP adjustment to examine the role of inflammation.

RESULTS

We identified 208 and 87 significant CpG sites associated with fibrinogen levels from the 450K (p < 1.03 × 10-7) and EPIC arrays (p < 5.78 × 10-8), respectively. There were 78 associations from the 450K array that replicated in the EPIC array and 26 vice versa. After accounting for overlapping sites, there were 83 replicated CpG sites located in 61 loci, of which only 4 have been previously reported for fibrinogen. The examples of genes located near these CpG sites were SOCS3 and AIM2, which are involved in inflammatory pathways. The associations of all 83 replicated CpG sites were attenuated after CRP adjustment, although many remained significant.

CONCLUSION

We identified 83 CpG sites associated with circulating fibrinogen levels. These associations are partially driven by inflammatory pathways shared by both fibrinogen and CRP.

Genome-wide analysis identifies genetic effects on reproductive success and ongoing natural selection at the FADS locus

Identifying genetic determinants of reproductive success may highlight mechanisms underlying fertility and identify alleles under present-day selection. Using data in 785,604 individuals of European ancestry, we identified 43 genomic loci associated with either number of children ever born (NEB) or childlessness. These loci span diverse aspects of reproductive biology, including puberty timing, age at first birth, sex hormone regulation, endometriosis and age at menopause. Missense variants in ARHGAP27 were associated with higher NEB but shorter reproductive lifespan, suggesting a trade-off at this locus between reproductive ageing and intensity. Other genes implicated by coding variants include PIK3IP1, ZFP82 and LRP4, and our results suggest a new role for the melanocortin 1 receptor (MC1R) in reproductive biology. As NEB is one component of evolutionary fitness, our identified associations indicate loci under present-day natural selection. Integration with data from historical selection scans highlighted an allele in the FADS1/2 gene locus that has been under selection for thousands of years and remains so today. Collectively, our findings demonstrate that a broad range of biological mechanisms contribute to reproductive success.

Genomics of perivascular space burden unravels early mechanisms of cerebral small vessel disease

Perivascular space (PVS) burden is an emerging, poorly understood, magnetic resonance imaging marker of cerebral small vessel disease, a leading cause of stroke and dementia. Genome-wide association studies in up to 40,095 participants (18 population-based cohorts, 66.3 ± 8.6 yr, 96.9% European ancestry) revealed 24 genome-wide significant PVS risk loci, mainly in the white matter. These were associated with white matter PVS already in young adults (N = 1,748; 22.1 ± 2.3 yr) and were enriched in early-onset leukodystrophy genes and genes expressed in fetal brain endothelial cells, suggesting early-life mechanisms. In total, 53% of white matter PVS risk loci showed nominally significant associations (27% after multiple-testing correction) in a Japanese population-based cohort (N = 2,862; 68.3 ± 5.3 yr). Mendelian randomization supported causal associations of high blood pressure with basal ganglia and hippocampal PVS, and of basal ganglia PVS and hippocampal PVS with stroke, accounting for blood pressure. Our findings provide insight into the biology of PVS and cerebral small vessel disease, pointing to pathways involving extracellular matrix, membrane transport and developmental processes, and the potential for genetically informed prioritization of drug targets.

OTTERS: a powerful TWAS framework leveraging summary-level reference data

Most existing TWAS tools require individual-level eQTL reference data and thus are not applicable to summary-level reference eQTL datasets. The development of TWAS methods that can harness summary-level reference data is valuable to enable TWAS in broader settings and enhance power due to increased reference sample size. Thus, we develop a TWAS framework called OTTERS (Omnibus Transcriptome Test using Expression Reference Summary data) that adapts multiple polygenic risk score (PRS) methods to estimate eQTL weights from summary-level eQTL reference data and conducts an omnibus TWAS. We show that OTTERS is a practical and powerful TWAS tool by both simulations and application studies.

Epigenome-Wide Association Study Reveals CpG Sites Associated with Thyroid Function and Regulatory Effects on KLF9

Background: Thyroid hormones play a key role in differentiation and metabolism and are known regulators of gene expression through both genomic and epigenetic processes including DNA methylation. The aim of this study was to examine associations between thyroid hormones and DNA methylation. Methods: We carried out a fixed-effect meta-analysis of epigenome-wide association study (EWAS) of blood DNA methylation sites from 8 cohorts from the ThyroidOmics Consortium, incorporating up to 7073 participants of both European and African ancestry, implementing a discovery and replication stage. Statistical analyses were conducted using normalized beta CpG values as dependent and log-transformed thyrotropin (TSH), free thyroxine, and free triiodothyronine levels, respectively, as independent variable in a linear model. The replicated findings were correlated with gene expression levels in whole blood and tested for causal influence of TSH and free thyroxine by two-sample Mendelian randomization (MR). Results: Epigenome-wide significant associations (p-value <1.1E-7) of three CpGs for free thyroxine, five for free triiodothyronine, and two for TSH concentrations were discovered and replicated (combined p-values = 1.5E-9 to 4.3E-28). The associations included CpG sites annotated to KLF9 (cg00049440) and DOT1L (cg04173586) that overlap with all three traits, consistent with hypothalamic-pituitary-thyroid axis physiology. Significant associations were also found for CpGs in FKBP5 for free thyroxine, and at CSNK1D/LINCO1970 and LRRC8D for free triiodothyronine. MR analyses supported a causal effect of thyroid status on DNA methylation of KLF9. DNA methylation of cg00049440 in KLF9 was inversely correlated with KLF9 gene expression in blood. The CpG at CSNK1D/LINC01970 overlapped with thyroid hormone receptor alpha binding peaks in liver cells. The total additive heritability of the methylation levels of the six significant CpG sites was between 25% and 57%. Significant methylation QTLs were identified for CpGs at KLF9, FKBP5, LRRC8D, and CSNK1D/LINC01970. Conclusions: We report novel associations between TSH, thyroid hormones, and blood-based DNA methylation. This study advances our understanding of thyroid hormone action particularly related to KLF9 and serves as a proof-of-concept that integrations of EWAS with other -omics data can provide a valuable tool for unraveling thyroid hormone signaling in humans by complementing and feeding classical in vitro and animal studies.

Allometric body shape indices, type 2 diabetes and kidney function: A two-sample Mendelian randomization study

AIM

To examine the association between body mass index (BMI)-independent allometric body shape indices and kidney function.

MATERIALS AND METHODS

We performed a two-sample Mendelian randomization (MR) analysis, using summary statistics from UK Biobank, CKDGen and DIAGRAM. BMI-independent allometric body shape indices were: A Body Shape Index (ABSI), Waist-Hip Index (WHI) and Hip Index (HI). Kidney function outcomes were: urinary albumin-to-creatinine ratio (UACR), estimated glomerular filtration rate and blood urea nitrogen. Furthermore, we investigated type 2 diabetes (T2D) as a potential mediator on the pathway to albuminuria. The main analysis was inverse variance-weighted random-effects MR in participants of European ancestry. We also performed several sensitivity MR analyses.

RESULTS

A 1-standard deviation (SD) increase in genetically predicted ABSI and WHI levels was associated with higher UACR (β = 0.039 [95% confidence interval: 0.016, 0.063] log [UACR], P = 0.001 for ABSI, and β = 0.028 [0.012, 0.044] log [UACR], P = 6 x 10^-4 for WHI) in women, but not in men. Meanwhile, a 1-SD increase in genetically predicted HI was associated with lower UACR in women (β = -0.021 [-0.041, 0.000] log [UACR], P = 0.05) and in men (β = -0.026 [-0.058, 0.005] log [UACR], P = 0.10). Corresponding estimates in individuals with diabetes were substantially augmented. Risk of T2D increased for genetically high ABSI and WHI in women (P < 6 x 10^-19 ) only, but decreased for genetically high HI in both sexes (P < 9 x 10^-3 ). No other associations were observed.

CONCLUSIONS

Genetically high HI was associated with decreased risk of albuminuria, mediated through decreased T2D risk in both sexes. Opposite associations applied to genetically high ABSI and WHI in women only.

Epigenetic and integrative cross-omics analyses of cerebral white matter hyperintensities on MRI

Cerebral white matter hyperintensities on MRI are markers of cerebral small vessel disease, a major risk factor for dementia and stroke. Despite the successful identification of multiple genetic variants associated with this highly heritable condition, its genetic architecture remains incompletely understood. More specifically, the role of DNA methylation has received little attention. We investigated the association between white matter hyperintensity burden and DNA methylation in blood at ∼450 000 cytosine-phosphate-guanine (CpG) sites in 9732 middle-aged to older adults from 14 community-based studies. Single CpG and region-based association analyses were carried out. Functional annotation and integrative cross-omics analyses were performed to identify novel genes underlying the relationship between DNA methylation and white matter hyperintensities. We identified 12 single CpG and 46 region-based DNA methylation associations with white matter hyperintensity burden. Our top discovery single CpG, cg24202936 (P = 7.6 × 10-8), was associated with F2 expression in blood (P = 6.4 × 10-5) and co-localized with FOLH1 expression in brain (posterior probability = 0.75). Our top differentially methylated regions were in PRMT1 and in CCDC144NL-AS1, which were also represented in single CpG associations (cg17417856 and cg06809326, respectively). Through Mendelian randomization analyses cg06809326 was putatively associated with white matter hyperintensity burden (P = 0.03) and expression of CCDC144NL-AS1 possibly mediated this association. Differentially methylated region analysis, joint epigenetic association analysis and multi-omics co-localization analysis consistently identified a role of DNA methylation near SH3PXD2A, a locus previously identified in genome-wide association studies of white matter hyperintensities. Gene set enrichment analyses revealed functions of the identified DNA methylation loci in the blood-brain barrier and in the immune response. Integrative cross-omics analysis identified 19 key regulatory genes in two networks related to extracellular matrix organization, and lipid and lipoprotein metabolism. A drug-repositioning analysis indicated antihyperlipidaemic agents, more specifically peroxisome proliferator-activated receptor-alpha, as possible target drugs for white matter hyperintensities. Our epigenome-wide association study and integrative cross-omics analyses implicate novel genes influencing white matter hyperintensity burden, which converged on pathways related to the immune response and to a compromised blood-brain barrier possibly due to disrupted cell-cell and cell-extracellular matrix interactions. The results also suggest that antihyperlipidaemic therapy may contribute to lowering risk for white matter hyperintensities possibly through protection against blood-brain barrier disruption.

TREML2 Gene Expression and Its Missense Variant rs3747742 Associate with White Matter Hyperintensity Volume and Alzheimer's Disease-Related Brain Atrophy in the General Population

Although the common pathology of Alzheimer's disease (AD) and white matter hyperintensities (WMH) is disputed, the gene TREML2 has been implicated in both conditions: its whole-blood gene expression was associated with WMH volume and its missense variant rs3747742 with AD risk. We re-examined those associations within one comprehensive dataset of the general population, additionally searched for cross-relations and illuminated the role of the apolipoprotein E (APOE) ε4 status in the associations. For our linear regression and linear mixed effect models, we used 1949 participants from the Study of Health in Pomerania (Germany). AD was assessed using a continuous pre-symptomatic MRI-based score evaluating a participant's AD-related brain atrophy. In our study, increased whole-blood TREML2 gene expression was significantly associated with reduced WMH volume but not with the AD score. Conversely, rs3747742-C was significantly associated with a reduced AD score but not with WMH volume. The APOE status did not influence the associations. In sum, TREML2 robustly associated with WMH volume and AD-related brain atrophy on different molecular levels. Our results thus underpin TREML2's role in neurodegeneration, might point to its involvement in AD and WMH via different biological mechanisms, and highlight TREML2 as a worthwhile target for disentangling the two pathologies.

Genome-wide meta-analyses reveal novel loci for verbal short-term memory and learning

Understanding the genomic basis of memory processes may help in combating neurodegenerative disorders. Hence, we examined the associations of common genetic variants with verbal short-term memory and verbal learning in adults without dementia or stroke (N = 53,637). We identified novel loci in the intronic region of CDH18, and at 13q21 and 3p21.1, as well as an expected signal in the APOE/APOC1/TOMM40 region. These results replicated in an independent sample. Functional and bioinformatic analyses supported many of these loci and further implicated POC1. We showed that polygenic score for verbal learning associated with brain activation in right parieto-occipital region during working memory task. Finally, we showed genetic correlations of these memory traits with several neurocognitive and health outcomes. Our findings suggest a role of several genomic loci in verbal memory processes.

A saturated map of common genetic variants associated with human height

Common single-nucleotide polymorphisms (SNPs) are predicted to collectively explain 40-50% of phenotypic variation in human height, but identifying the specific variants and associated regions requires huge sample sizes1. Here, using data from a genome-wide association study of 5.4 million individuals of diverse ancestries, we show that 12,111 independent SNPs that are significantly associated with height account for nearly all of the common SNP-based heritability. These SNPs are clustered within 7,209 non-overlapping genomic segments with a mean size of around 90 kb, covering about 21% of the genome. The density of independent associations varies across the genome and the regions of increased density are enriched for biologically relevant genes. In out-of-sample estimation and prediction, the 12,111 SNPs (or all SNPs in the HapMap 3 panel2) account for 40% (45%) of phenotypic variance in populations of European ancestry but only around 10-20% (14-24%) in populations of other ancestries. Effect sizes, associated regions and gene prioritization are similar across ancestries, indicating that reduced prediction accuracy is likely to be explained by linkage disequilibrium and differences in allele frequency within associated regions. Finally, we show that the relevant biological pathways are detectable with smaller sample sizes than are needed to implicate causal genes and variants. Overall, this study provides a comprehensive map of specific genomic regions that contain the vast majority of common height-associated variants. Although this map is saturated for populations of European ancestry, further research is needed to achieve equivalent saturation in other ancestries.

Genetic loci and prioritization of genes for kidney function decline derived from a meta-analysis of 62 longitudinal genome-wide association studies

Estimated glomerular filtration rate (eGFR) reflects kidney function. Progressive eGFR-decline can lead to kidney failure, necessitating dialysis or transplantation. Hundreds of loci from genome-wide association studies (GWAS) for eGFR help explain population cross section variability. Since the contribution of these or other loci to eGFR-decline remains largely unknown, we derived GWAS for annual eGFR-decline and meta-analyzed 62 longitudinal studies with eGFR assessed twice over time in all 343,339 individuals and in high-risk groups. We also explored different covariate adjustment. Twelve genome-wide significant independent variants for eGFR-decline unadjusted or adjusted for eGFR-baseline (11 novel, one known for this phenotype), including nine variants robustly associated across models were identified. All loci for eGFR-decline were known for cross-sectional eGFR and thus distinguished a subgroup of eGFR loci. Seven of the nine variants showed variant-by-age interaction on eGFR cross section (further about 350,000 individuals), which linked genetic associations for eGFR-decline with age-dependency of genetic cross-section associations. Clinically important were two to four-fold greater genetic effects on eGFR-decline in high-risk subgroups. Five variants associated also with chronic kidney disease progression mapped to genes with functional in-silico evidence (UMOD, SPATA7, GALNTL5, TPPP). An unfavorable versus favorable nine-variant genetic profile showed increased risk odds ratios of 1.35 for kidney failure (95% confidence intervals 1.03-1.77) and 1.27 for acute kidney injury (95% confidence intervals 1.08-1.50) in over 2000 cases each, with matched controls). Thus, we provide a large data resource, genetic loci, and prioritized genes for kidney function decline, which help inform drug development pipelines revealing important insights into the age-dependency of kidney function genetics.

Genome-wide association analyses of physical activity and sedentary behavior provide insights into underlying mechanisms and roles in disease prevention

Although physical activity and sedentary behavior are moderately heritable, little is known about the mechanisms that influence these traits. Combining data for up to 703,901 individuals from 51 studies in a multi-ancestry meta-analysis of genome-wide association studies yields 99 loci that associate with self-reported moderate-to-vigorous intensity physical activity during leisure time (MVPA), leisure screen time (LST) and/or sedentary behavior at work. Loci associated with LST are enriched for genes whose expression in skeletal muscle is altered by resistance training. A missense variant in ACTN3 makes the alpha-actinin-3 filaments more flexible, resulting in lower maximal force in isolated type IIA muscle fibers, and possibly protection from exercise-induced muscle damage. Finally, Mendelian randomization analyses show that beneficial effects of lower LST and higher MVPA on several risk factors and diseases are mediated or confounded by body mass index (BMI). Our results provide insights into physical activity mechanisms and its role in disease prevention.

Virtual Ontogeny of Cortical Growth Preceding Mental Illness

BACKGROUND

Morphology of the human cerebral cortex differs across psychiatric disorders, with neurobiology and developmental origins mostly undetermined. Deviations in the tangential growth of the cerebral cortex during pre/perinatal periods may be reflected in individual variations in cortical surface area later in life.

METHODS

Interregional profiles of group differences in surface area between cases and controls were generated using T1-weighted magnetic resonance imaging from 27,359 individuals including those with attention-deficit/hyperactivity disorder, autism spectrum disorder, bipolar disorder, major depressive disorder, schizophrenia, and high general psychopathology (through the Child Behavior Checklist). Similarity of interregional profiles of group differences in surface area and prenatal cell-specific gene expression was assessed.

RESULTS

Across the 11 cortical regions, group differences in cortical area for attention-deficit/hyperactivity disorder, schizophrenia, and Child Behavior Checklist were dominant in multimodal association cortices. The same interregional profiles were also associated with interregional profiles of (prenatal) gene expression specific to proliferative cells, namely radial glia and intermediate progenitor cells (greater expression, larger difference), as well as differentiated cells, namely excitatory neurons and endothelial and mural cells (greater expression, smaller difference). Finally, these cell types were implicated in known pre/perinatal risk factors for psychosis. Genes coexpressed with radial glia were enriched with genes implicated in congenital abnormalities, birth weight, hypoxia, and starvation. Genes coexpressed with endothelial and mural genes were enriched with genes associated with maternal hypertension and preterm birth.

CONCLUSIONS

Our findings support a neurodevelopmental model of vulnerability to mental illness whereby prenatal risk factors acting through cell-specific processes lead to deviations from typical brain development during pregnancy.

Limited evidence for blood eQTLs in human sexual dimorphism

BACKGROUND

The genetic underpinning of sexual dimorphism is very poorly understood. The prevalence of many diseases differs between men and women, which could be in part caused by sex-specific genetic effects. Nevertheless, only a few published genome-wide association studies (GWAS) were performed separately in each sex. The reported enrichment of expression quantitative trait loci (eQTLs) among GWAS-associated SNPs suggests a potential role of sex-specific eQTLs in the sex-specific genetic mechanism underlying complex traits.

METHODS

To explore this scenario, we combined sex-specific whole blood RNA-seq eQTL data from 3447 European individuals included in BIOS Consortium and GWAS data from UK Biobank. Next, to test the presence of sex-biased causal effect of gene expression on complex traits, we performed sex-specific transcriptome-wide Mendelian randomization (TWMR) analyses on the two most sexually dimorphic traits, waist-to-hip ratio (WHR) and testosterone levels. Finally, we performed power analysis to calculate the GWAS sample size needed to observe sex-specific trait associations driven by sex-biased eQTLs.

RESULTS

Among 9 million SNP-gene pairs showing sex-combined associations, we found 18 genes with significant sex-biased cis-eQTLs (FDR 5%). Our phenome-wide association study of the 18 top sex-biased eQTLs on >700 traits unraveled that these eQTLs do not systematically translate into detectable sex-biased trait-associations. In addition, we observed that sex-specific causal effects of gene expression on complex traits are not driven by sex-specific eQTLs. Power analyses using real eQTL- and causal-effect sizes showed that millions of samples would be necessary to observe sex-biased trait associations that are fully driven by sex-biased cis-eQTLs. Compensatory effects may further hamper their detection.

CONCLUSIONS

Our results suggest that sex-specific eQTLs in whole blood do not translate to detectable sex-specific trait associations of complex diseases, and vice versa that the observed sex-specific trait associations cannot be explained by sex-specific eQTLs.

Analysis of epidemiological association patterns of serum thyrotropin by combining random forests and Bayesian networks

BACKGROUND

Approaching epidemiological data with flexible machine learning algorithms is of great value for understanding disease-specific association patterns. However, it can be difficult to correctly extract and understand those patterns due to the lack of model interpretability.

METHOD

We here propose a machine learning workflow that combines random forests with Bayesian network surrogate models to allow for a deeper level of interpretation of complex association patterns. We first evaluate the proposed workflow on synthetic data. We then apply it to data from the large population-based Study of Health in Pomerania (SHIP). Based on this combination, we discover and interpret broad patterns of individual serum TSH concentrations, an important marker of thyroid functionality.

RESULTS

Evaluations using simulated data show that feature associations can be correctly recovered by combining random forests and Bayesian networks. The presented model achieves predictive accuracy that is similar to state-of-the-art models (root mean square error of 0.66, mean absolute error of 0.55, coefficient of determination of R2 = 0.15). We identify 62 relevant features from the final random forest model, ranging from general health variables over dietary and genetic factors to physiological, hematological and hemostasis parameters. The Bayesian network model is used to put these features into context and make the black-box random forest model more understandable.

CONCLUSION

We demonstrate that the combination of random forest and Bayesian network analysis is helpful to reveal and interpret broad association patterns of individual TSH concentrations. The discovered patterns are in line with state-of-the-art literature. They may be useful for future thyroid research and improved dosing of therapeutics.

Gene-mapping study of extremes of cerebral small vessel disease reveals TRIM47 as a strong candidate

Cerebral small vessel disease is a leading cause of stroke and a major contributor to cognitive decline and dementia, but our understanding of specific genes underlying the cause of sporadic cerebral small vessel disease is limited. We report a genome-wide association study and a whole-exome association study on a composite extreme phenotype of cerebral small vessel disease derived from its most common MRI features: white matter hyperintensities and lacunes. Seventeen population-based cohorts of older persons with MRI measurements and genome-wide genotyping (n = 41 326), whole-exome sequencing (n = 15 965), or exome chip (n = 5249) data contributed 13 776 and 7079 extreme small vessel disease samples for the genome-wide association study and whole-exome association study, respectively. The genome-wide association study identified significant association of common variants in 11 loci with extreme small vessel disease, of which the chr12q24.11 locus was not previously reported to be associated with any MRI marker of cerebral small vessel disease. The whole-exome association study identified significant associations of extreme small vessel disease with common variants in the 5' UTR region of EFEMP1 (chr2p16.1) and one probably damaging common missense variant in TRIM47 (chr17q25.1). Mendelian randomization supports the causal association of extensive small vessel disease severity with increased risk of stroke and Alzheimer's disease. Combined evidence from summary-based Mendelian randomization studies and profiling of human loss-of-function allele carriers showed an inverse relation between TRIM47 expression in the brain and blood vessels and extensive small vessel disease severity. We observed significant enrichment of Trim47 in isolated brain vessel preparations compared to total brain fraction in mice, in line with the literature showing Trim47 enrichment in brain endothelial cells at single cell level. Functional evaluation of TRIM47 by small interfering RNAs-mediated knockdown in human brain endothelial cells showed increased endothelial permeability, an important hallmark of cerebral small vessel disease pathology. Overall, our comprehensive gene-mapping study and preliminary functional evaluation suggests a putative role of TRIM47 in the pathophysiology of cerebral small vessel disease, making it an important candidate for extensive in vivo explorations and future translational work.

Differential and shared genetic effects on kidney function between diabetic and non-diabetic individuals

Reduced glomerular filtration rate (GFR) can progress to kidney failure. Risk factors include genetics and diabetes mellitus (DM), but little is known about their interaction. We conducted genome-wide association meta-analyses for estimated GFR based on serum creatinine (eGFR), separately for individuals with or without DM (nDM = 178,691, nnoDM = 1,296,113). Our genome-wide searches identified (i) seven eGFR loci with significant DM/noDM-difference, (ii) four additional novel loci with suggestive difference and (iii) 28 further novel loci (including CUBN) by allowing for potential difference. GWAS on eGFR among DM individuals identified 2 known and 27 potentially responsible loci for diabetic kidney disease. Gene prioritization highlighted 18 genes that may inform reno-protective drug development. We highlight the existence of DM-only and noDM-only effects, which can inform about the target group, if respective genes are advanced as drug targets. Largely shared effects suggest that most drug interventions to alter eGFR should be effective in DM and noDM.

The Effects of Common Genetic Variation in 96 Genes Involved in Thyroid Hormone Regulation on TSH and FT4 Concentrations

OBJECTIVE

While most of the variation in thyroid function is determined by genetic factors, single nucleotide polymorphisms (SNPs) identified via genome-wide association analyses have only explained ~5% to 9% of this variance so far. Most SNPs were in or nearby genes with no known role in thyroid hormone (TH) regulation. Therefore, we performed a large-scale candidate gene study investigating the effect of common genetic variation in established TH regulating genes on serum thyrotropin [thyroid-stimulating hormone (TSH)] and thyroxine (FT4) concentrations.

METHODS

SNPs in or within 10 kb of 96 TH regulating genes were included (30 031 TSH SNPs, and 29 962 FT4 SNPs). Associations were studied in 54 288 individuals from the ThyroidOmics Consortium. Linkage disequilibrium-based clumping was used to identify independently associated SNPs. SNP-based explained variances were calculated using SumHer software.

RESULTS

We identified 23 novel TSH-associated SNPs in predominantly hypothalamic-pituitary-thyroid axis genes and 25 novel FT4-associated SNPs in mainly peripheral metabolism and transport genes. Genome-wide SNP variation explained ~21% (SD 1.7) of the total variation in both TSH and FT4 concentrations, whereas SNPs in the 96 TH regulating genes explained 1.9% to 2.6% (SD 0.4).

CONCLUSION

Here we report the largest candidate gene analysis on thyroid function, resulting in a substantial increase in the number of genetic variants determining TSH and FT4 concentrations. Interestingly, these candidate gene SNPs explain only a minor part of the variation in TSH and FT4 concentrations, which substantiates the need for large genetic studies including common and rare variants to unravel novel, yet unknown, pathways in TH regulation.

DNA methylation signature of chronic low-grade inflammation and its role in cardio-respiratory diseases

We performed a multi-ethnic Epigenome Wide Association study on 22,774 individuals to describe the DNA methylation signature of chronic low-grade inflammation as measured by C-Reactive protein (CRP). We find 1,511 independent differentially methylated loci associated with CRP. These CpG sites show correlation structures across chromosomes, and are primarily situated in euchromatin, depleted in CpG islands. These genomic loci are predominantly situated in transcription factor binding sites and genomic enhancer regions. Mendelian randomization analysis suggests altered CpG methylation is a consequence of increased blood CRP levels. Mediation analysis reveals obesity and smoking as important underlying driving factors for changed CpG methylation. Finally, we find that an activated CpG signature significantly increases the risk for cardiometabolic diseases and COPD.

Genetic variants associated with longitudinal changes in brain structure across the lifespan

Human brain structure changes throughout the lifespan. Altered brain growth or rates of decline are implicated in a vast range of psychiatric, developmental and neurodegenerative diseases. In this study, we identified common genetic variants that affect rates of brain growth or atrophy in what is, to our knowledge, the first genome-wide association meta-analysis of changes in brain morphology across the lifespan. Longitudinal magnetic resonance imaging data from 15,640 individuals were used to compute rates of change for 15 brain structures. The most robustly identified genes GPR139, DACH1 and APOE are associated with metabolic processes. We demonstrate global genetic overlap with depression, schizophrenia, cognitive functioning, insomnia, height, body mass index and smoking. Gene set findings implicate both early brain development and neurodegenerative processes in the rates of brain changes. Identifying variants involved in structural brain changes may help to determine biological pathways underlying optimal and dysfunctional brain development and aging.

Identifying the Common Genetic Basis of Antidepressant Response

Background

Antidepressants are a first-line treatment for depression. However, only a third of individuals experience remission after the first treatment. Common genetic variation, in part, likely regulates antidepressant response, yet the success of previous genome-wide association studies has been limited by sample size. This study performs the largest genetic analysis of prospectively assessed antidepressant response in major depressive disorder to gain insight into the underlying biology and enable out-of-sample prediction.

Methods

Genome-wide analysis of remission (n remit = 1852, n nonremit = 3299) and percentage improvement (n = 5218) was performed. Single nucleotide polymorphism-based heritability was estimated using genome-wide complex trait analysis. Genetic covariance with eight mental health phenotypes was estimated using polygenic scores/AVENGEME. Out-of-sample prediction of antidepressant response polygenic scores was assessed. Gene-level association analysis was performed using MAGMA and transcriptome-wide association study. Tissue, pathway, and drug binding enrichment were estimated using MAGMA.

Results

Neither genome-wide association study identified genome-wide significant associations. Single nucleotide polymorphism-based heritability was significantly different from zero for remission (h 2 = 0.132, SE = 0.056) but not for percentage improvement (h 2 = -0.018, SE = 0.032). Better antidepressant response was negatively associated with genetic risk for schizophrenia and positively associated with genetic propensity for educational attainment. Leave-one-out validation of antidepressant response polygenic scores demonstrated significant evidence of out-of-sample prediction, though results varied in external cohorts. Gene-based analyses identified ETV4 and DHX8 as significantly associated with antidepressant response.

Conclusions

This study demonstrates that antidepressant response is influenced by common genetic variation, has a genetic overlap schizophrenia and educational attainment, and provides a useful resource for future research. Larger sample sizes are required to attain the potential of genetics for understanding and predicting antidepressant response.

Effects of Apolipoprotein E polymorphism on carotid intima-media thickness, incident myocardial infarction and incident stroke

The Apolipoprotein E (APOE) gene polymorphism (rs429358 and rs7412) shows a well-established association with lipid profiles, but its effect on cardiovascular disease is still conflicting. Therefore, we examined the association of different APOE alleles with common carotid artery intima-media thickness (CCA-IMT), carotid plaques, incident myocardial infarction (MI) and stroke. We analyzed data from 3327 participants aged 20-79 years of the population-based Study of Health in Pomerania (SHIP) from Northeast Germany with a median follow-up time of 14.5 years. Linear, logistic, and Cox-regression models were used to assess the associations of the APOE polymorphism with CCA-IMT, carotid plaques, incident MI and stroke, respectively. In our study, the APOE E2 allele was associated with lower CCA-IMT at baseline compared to E3 homozygotes (β: - 0.02 [95% CI - 0.04, - 0.004]). Over the follow-up, 244 MI events and 218 stroke events were observed. APOE E2 and E4 allele were not associated with incident MI (E2 HR: 1.06 [95% CI 0.68, 1.66]; E4 HR: 1.03 [95% CI 0.73, 1.45]) and incident stroke (E2 HR: 0.79 [95% CI 0.48, 1.30]; E4 HR: 0.96 [95% CI 0.66, 1.38]) in any of the models adjusting for potential confounders. However, the positive association between CCA-IMT and incident MI was more pronounced in E2 carriers than E3 homozygotes. Thus, our study suggests that while APOE E2 allele may predispose individuals to lower CCA-IMT, E2 carriers may be more prone to MI than E3 homozygotes as the CCA-IMT increases. APOE E4 allele had no effect on CCA-IMT, plaques, MI or stroke.

Obesity and Kidney Function: A Two-Sample Mendelian Randomization Study

BACKGROUND

Obesity and type 2 diabetes (T2D) are correlated risk factors for chronic kidney disease (CKD).

METHODS

Using summary data from GIANT (Genetic Investigation of Anthropometric Traits), DIAGRAM (DIAbetes Genetics Replication And Meta-analysis), and CKDGen (CKD Genetics), we examined causality and directionality of the association between obesity and kidney function. Bidirectional 2-sample Mendelian randomization (MR) estimated the total causal effects of body mass index (BMI) and waist-to-hip ratio (WHR) on kidney function, and vice versa. Effects of adverse obesity and T2D were examined by stratifying BMI variants by their association with WHR and T2D. Multivariable MR estimated the direct causal effects of BMI and WHR on kidney function. The inverse variance weighted random-effects MR for Europeans was the main analysis, accompanied by several sensitivity MR analyses.

RESULTS

One standard deviation (SD ≈ 4.8 kg/m2) genetically higher BMI was associated with decreased estimated glomerular filtration rate (eGFR) [β=-0.032 (95% confidence intervals: -0.036, -0.027) log[eGFR], P = 1 × 10-43], increased blood urea nitrogen (BUN) [β = 0.010 (0.005, 0.015) log[BUN], P = 3 × 10-6], increased urinary albumin-to-creatinine ratio [β = 0.199 (0.067, 0.332) log[urinary albumin-to-creatinine ratio (UACR)], P = 0.003] in individuals with diabetes, and increased risk of microalbuminuria [odds ratios (OR) = 1.15 [1.04-1.28], P = 0.009] and CKD [1.13 (1.07-1.19), P = 3 × 10-6]. Corresponding estimates for WHR and for trans-ethnic populations were overall similar. The associations were driven by adverse obesity, and for microalbuminuria additionally by T2D. While genetically high BMI, unlike WHR, was directly associated with eGFR, BUN, and CKD, the pathway to albuminuria was likely through T2D. Genetically predicted kidney function was not associated with BMI or WHR.

CONCLUSIONS

Genetically high BMI is associated with impaired kidney function, driven by adverse obesity, and for albuminuria additionally by T2D.

Association of Telomere Length With Risk of Disease and Mortality

IMPORTANCE

Telomeres protect DNA from damage. Because they shorten with each mitotic cycle, leukocyte telomere length (LTL) serves as a mitotic clock. Reduced LTL has been associated with multiple human disorders.

OBJECTIVE

To determine the association between LTL and overall as well as disease-specific mortality and morbidity.

DESIGN, SETTING, AND PARTICIPANTS

This multicenter, community-based cohort study conducted from March 2006 to December 2010 included longitudinal follow-up (mean [SD], 12 [2] years) for 472 432 English participants from the United Kingdom Biobank (UK Biobank) and analyzed morbidity and mortality. The data were analyzed in 2021.

MAIN OUTCOMES AND MEASURES

Hazard ratios (HRs) and odds ratios for mortality and morbidity associated with a standard deviation change in LTL, adjusted for age, sex, body mass index (calculated as weight in kilograms divided by height in meters squared), and ethnicity.

RESULTS

This study included a total of 472 432 English participants, of whom 54% were women (mean age, 57 years). Reduced LTL was associated with increased overall (HR, 1.08; 95% CI, 1.07-1.09), cardiovascular (HR, 1.09; 95% CI, 1.06-1.12), respiratory (HR, 1.40; 95% CI, 1.34-1.45), digestive (HR, 1.26; 95% CI, 1.19-1.33), musculoskeletal (HR, 1.51; 95% CI, 1.35-1.92), and COVID-19 (HR, 1.15; 95% CI, 1.07-1.23) mortality, but not cancer-related mortality. A total of 214 disorders were significantly overrepresented and 37 underrepresented in participants with shorter LTL. Respiratory (11%), digestive/liver-related (14%), circulatory (18%), and musculoskeletal conditions (6%), together with infections (5%), accounted for most positive associations, whereas (benign) neoplasms and endocrinologic/metabolic disorders were the most underrepresented entities. Malignant tumors, esophageal cancer, and lymphoid and myeloid leukemia were significantly more common in participants with shorter LTL, whereas brain cancer and melanoma were less prevalent. While smoking and alcohol consumption were associated with shorter LTL, additional adjustment for both factors, as well as cognitive function/major comorbid conditions, did not significantly alter the results.

CONCLUSIONS AND RELEVANCE

This cohort study found that shorter LTL was associated with a small risk increase of overall mortality, but a higher risk of mortality was associated with specific organs and diseases.

Longitudinal association of Apolipoprotein E polymorphism with lipid profile, type 2 diabetes and metabolic syndrome: Results from a 15 year follow-up study

AIMS

To examine the association of different APOE alleles with type 2 diabetes mellitus (T2DM) and metabolic syndrome (MetS) as well as the influence of high-sensitive C-reactive protein (hs-CRP) on these associations.

METHODS

We analyzed data from 3917 participants aged 20-81 years of the population-based Study of Health in Pomerania (SHIP) from Northeast Germany with a median follow-up time of 10.8 years. Linear and logistic mixed models were performed to test the association of APOE alleles with T2DM and MetS.

RESULTS

We observed 393 T2DM and 1411 MetS events at baseline, and 576 T2DM and 1342 MetS events over the follow-up. The E4 carriers had a lower odds of developing T2DM (OR: 0.47 [0.24, 0.94]) than E3 homozygotes even after adjustment for potential confounders. The E2 carriers showed no associations. The inverse association between E4 alleles and T2DM moderately attenuated after adjustment for hs-CRP levels. The lower odds of developing T2DM in E4 carriers was more pronounced in participants without obesity, hypertension or MetS. However, both E2 and E4 carriers had higher odds of developing MetS (E2 OR: 1.45 [1.03, 2.03]; E4 OR: 1.56 [1.17, 2.09]) than E3 homozygotes.

CONCLUSIONS

While the presence of APOE E4 allele might increase the chance of MetS through its major action on lipids, E4 allele might offer a protection towards T2DM through its influence on inflammation.

Genome-wide interaction study with major depression identifies novel variants associated with cognitive function

Major Depressive Disorder (MDD) often is associated with significant cognitive dysfunction. We conducted a meta-analysis of genome-wide interaction of MDD and cognitive function using data from four large European cohorts in a total of 3510 MDD cases and 6057 controls. In addition, we conducted analyses using polygenic risk scores (PRS) based on data from the Psychiatric Genomics Consortium (PGC) on the traits of MDD, Bipolar disorder (BD), Schizophrenia (SCZ), and mood instability (MIN). Functional exploration contained gene expression analyses and Ingenuity Pathway Analysis (IPA®). We identified a set of significantly interacting single nucleotide polymorphisms (SNPs) between MDD and the genome-wide association study (GWAS) of cognitive domains of executive function, processing speed, and global cognition. Several of these SNPs are located in genes expressed in brain, with important roles such as neuronal development (REST), oligodendrocyte maturation (TNFRSF21), and myelination (ARFGEF1). IPA® identified a set of core genes from our dataset that mapped to a wide range of canonical pathways and biological functions (MPO, FOXO1, PDE3A, TSLP, NLRP9, ADAMTS5, ROBO1, REST). Furthermore, IPA® identified upstream regulator molecules and causal networks impacting on the expression of dataset genes, providing a genetic basis for further clinical exploration (vitamin D receptor, beta-estradiol, tadalafil). PRS of MIN and meta-PRS of MDD, MIN and SCZ were significantly associated with all cognitive domains. Our results suggest several genes involved in physiological processes for the development and maintenance of cognition in MDD, as well as potential novel therapeutic agents that could be explored in patients with MDD associated cognitive dysfunction.

Discovery of novel eGFR-associated multiple independent signals using a quasi-adaptive method

A decreased estimated glomerular filtration rate (eGFR) leading to chronic kidney disease is a significant public health problem. Kidney function is a heritable trait, and recent application of genome-wide association studies (GWAS) successfully identified multiple eGFR-associated genetic loci. To increase statistical power for detecting independent associations in GWAS loci, we improved our recently developed quasi-adaptive method estimating SNP-specific alpha levels for the conditional analysis, and applied it to the GWAS meta-analysis results of eGFR among 783,978 European-ancestry individuals. Among known eGFR loci, we revealed 19 new independent association signals that were subsequently replicated in the United Kingdom Biobank (n = 408,608). These associations have remained undetected by conditional analysis using the established conservative genome-wide significance level of 5 × 10^-8. Functional characterization of known index SNPs and novel independent signals using colocalization of conditional eGFR association results and gene expression in cis across 51 human tissues identified two potentially causal genes across kidney tissues: TSPAN33 and TFDP2, and three candidate genes across other tissues: SLC22A2, LRP2, and CDKN1C. These colocalizations were not identified in the original GWAS. By applying our improved quasi-adaptive method, we successfully identified additional genetic variants associated with eGFR. Considering these signals in colocalization analyses can increase the precision of revealing potentially functional genes of GWAS loci.