The UK Adult Twin Registry (TwinsUK Resource)

Testing the impact of trait prevalence priors in Bayesian-based genetic prediction modeling of human appearance traits

The prediction of appearance traits by use of solely genetic information has become an established approach and a number of statistical prediction models have already been developed for this purpose. However, given limited knowledge on appearance genetics, currently available models are incomplete and do not include all causal genetic variants as predictors. Therefore such prediction models may benefit from the inclusion of additional information that acts as a proxy for this unknown genetic background. Use of priors, possibly informed by trait category prevalence values in biogeographic ancestry groups, in a Bayesian framework may thus improve the prediction accuracy of previously predicted externally visible characteristics, but has not been investigated as of yet. In this study, we assessed the impact of using trait prevalence-informed priors on the prediction performance in Bayesian models for eye, hair and skin color as well as hair structure and freckles in comparison to the respective prior-free models. Those prior-free models were either similarly defined either very close to the already established ones by using a reduced predictive marker set. However, these differences in the number of the predictive markers should not affect significantly our main outcomes. We observed that such priors often had a strong effect on the prediction performance, but to varying degrees between different traits and also different trait categories, with some categories barely showing an effect. While we found potential for improving the prediction accuracy of many of the appearance trait categories tested by using priors, our analyses also showed that misspecification of those prior values often severely diminished the accuracy compared to the respective prior-free approach. This emphasizes the importance of accurate specification of prevalence-informed priors in Bayesian prediction modeling of appearance traits. However, the existing literature knowledge on spatial prevalence is sparse for most appearance traits, including those investigated here. Due to the limitations in appearance trait prevalence knowledge, our results render the use of trait prevalence-informed priors in DNA-based appearance trait prediction currently infeasible.

Reference exome data for a Northern Brazilian population

Exome sequencing is widely used in the diagnosis of rare genetic diseases and provides useful variant data for analysis of complex diseases. There is not always adequate population-specific reference data to assist in assigning a diagnostic variant to a specific clinical condition. Here we provide a catalogue of variants called after sequencing the exomes of 45 babies from Rio Grande do Nord in Brazil. Sequence data were processed using an 'intersect-then-combine' (ITC) approach, using GATK and SAMtools to call variants. A total of 612,761 variants were identified in at least one individual in this Brazilian Cohort, including 559,448 single nucleotide variants (SNVs) and 53,313 insertion/deletions. Of these, 58,111 overlapped with nonsynonymous (nsSNVs) or splice site (ssSNVs) SNVs in dbNSFP. As an aid to clinical diagnosis of rare diseases, we used the American College of Medicine Genetics and Genomics (ACMG) guidelines to assign pathogenic/likely pathogenic status to 185 (0.32%) of the 58,111 nsSNVs and ssSNVs. Our data set provides a useful reference point for diagnosis of rare diseases in Brazil. (169 words).

Genetic variations in estrogen and progesterone pathway genes in preeclampsia patients and controls in Bavaria

PURPOSE

Hypertensive pregnancy disorders and preeclampsia are major causes of maternal and fetal morbidity and mortality worldwide. Many different organs are involved in the diseases' clinical phenotype. The underlying mechanism is still unknown, with a possible genetic component. This case-control study investigated effects on the risk of preeclampsia of genetic variations (single nucleotide polymorphisms, SNPs) in the estrogen and progesterone pathway genes.

METHODS

The study included 167 patients with preeclampsia and 115 healthy controls from the "Franconian Maternal Health Evaluation Studies" (FRAMES). All patients completed an epidemiological questionnaire, data from which were correlated with prospective data on pregnancy and labor. DNA was isolated from blood samples and genotyping was done by PCR. Variants in the aromatase gene CYP19A1 (rs10046, rs4646), progesterone receptor gene (rs1042838, rs10895068), and estrogen receptor-α gene (rs488133) were examined, and the genotype distribution in the two groups was analyzed statistically.

RESULTS

A significant difference in the distribution frequency of genotypes between preeclampsia patients and controls was identified in one of the five SNPs. For rs10895068 in the progesterone receptor gene, genotype G/A was significantly more frequent among cases than controls (P = 0.023). No significant differences between the two cohorts were found in the other SNPs.

CONCLUSIONS

This study showed a significant association between only one SNP in the progesterone receptor and preeclampsia. Other studies have also noted genetic aspects of preeclampsia. The underlying mechanism and causal relationship are not yet known, and further research is needed to explain the extent of genetic variations and the causal relationship in preeclampsia.

FADS genetic and metabolomic analyses identify the ∆5 desaturase (FADS1) step as a critical control point in the formation of biologically important lipids

Humans have undergone intense evolutionary selection to optimize their capacity to generate necessary quantities of long chain (LC-) polyunsaturated fatty acid (PUFA)-containing lipids. To better understand the impact of genetic variation within a locus of three FADS genes (FADS1, FADS2, and FADS3) on a diverse family of lipids, we examined the associations of 247 lipid metabolites (including four major classes of LC-PUFA-containing molecules and signaling molecules) with common and low-frequency genetic variants located within the FADS locus. Genetic variation in the FADS locus was strongly associated (p < 1.2 × 10^–8) with 52 LC-PUFA-containing lipids and signaling molecules, including free fatty acids, phospholipids, lyso-phospholipids, and an endocannabinoid. Notably, the majority (80%) of FADS-associated lipids were not significantly associated with genetic variants outside of this FADS locus. These findings highlight the central role genetic variation at the FADS locus plays in regulating levels of physiologically critical LC-PUFA-containing lipids that participate in innate immunity, energy homeostasis, and brain development/function.

Genetic aetiology of blood pressure relates to aortic stiffness with bi-directional causality: evidence from heritability, blood pressure polymorphisms, and Mendelian randomization

AIMS

Haemodynamic determinants of blood pressure (BP) include cardiac output (CO), systemic vascular resistance (SVR), and arterial stiffness. We investigated the heritability of these phenotypes, their association with BP-related single-nucleotide polymorphisms (SNPs), and the causal association between BP and arterial stiffness.

METHODS AND RESULTS

We assessed BP, central BP components, and haemodynamic properties (during a single visit) including CO, SVR, and pulse wave velocity (PWV, measure of arterial stiffness) in 3531 (1934 monozygotic, 1586 dizygotic) female TwinsUK participants. Heritability was estimated using structural equation modelling. Association with 984 BP-associated SNP was examined using least absolute shrinkage and selection operator (LASSO) and generalized estimating equation regression. One and two-sample Mendelian randomization (MR) was used to estimate the causal direction between BP and arterial stiffness including data on 436 419 UK Biobank participants. We found high heritability for systolic and pulsatile components of BP (>50%) and PWV (65%) with overlapping genes accounting for >50% of their observed correlation. Environmental factors explained most of the variability of CO and SVR (>80%). Regression identified SNPs (n = 5) known to be associated with BP to also be associated with PWV. One-sample MR showed evidence of bi-directional causal association between BP and PWV in TwinsUK participants. Two-sample MR, confirmed a bi-directional causal effect of PWV on BP (inverse variance weighted (IVW) beta = 0.11, P < 0.02) and BP on arterial stiffness (IVW beta = 0.004, P < 0.0001).

CONCLUSION

The genetic basis of BP is mediated not only by genes regulating BP but also by genes that influence arterial stiffness. Mendelian randomization indicates a bi-directional causal association between BP and arterial stiffness.

Impact of host genetics on gut microbiome: Take-home lessons from human and mouse studies

The intestinal microbiome has emerged as an important component involved in various diseases. Therefore, the interest in understanding the factors shaping its composition is growing. The gut microbiome, often defined as a complex trait, contains diverse components and its properties are determined by a combination of external and internal effects. Although much effort has been invested so far, it is still difficult to evaluate the extent to which human genetics shape the composition of the gut microbiota. However, in mouse studies, where the environmental factors are better controlled, the effect of the genetic background was significant. The purpose of this paper is to provide a current assessment of the role of human host genetics in shaping the gut microbiome composition. Despite the inconsistency of the reported results, it can be estimated that the genetic factor affects a portion of the microbiome. However, this effect is currently lower than the initial estimates, and it is difficult to separate the genetic influence from the environmental effect. Additionally, despite the differences between the microbial composition of humans and mice, results from mouse models can strengthen our knowledge of host genetics underlying the human gut microbial variation.

ACE2 expression in adipose tissue is associated with COVID-19 cardio-metabolic risk factors and cell type composition

COVID-19 severity has varied widely, with demographic and cardio-metabolic factors increasing risk of severe reactions to SARS-CoV-2 infection, but the underlying mechanisms for this remain uncertain. We investigated phenotypic and genetic factors associated with subcutaneous adipose tissue expression of Angiotensin I Converting Enzyme 2 ( ACE2 ), which has been shown to act as a receptor for SARS-CoV-2 cellular entry. In a meta-analysis of three independent studies including up to 1,471 participants, lower adipose tissue ACE2 expression was associated with adverse cardio-metabolic health indices including type 2 diabetes (T2D) and obesity status, higher serum fasting insulin and BMI, and lower serum HDL levels (P<5.32x10 -4 ). ACE2 expression levels were also associated with estimated proportions of cell types in adipose tissue; lower ACE2 expression was associated with a lower proportion of microvascular endothelial cells (P=4.25x10 -4 ) and higher macrophage proportion (P=2.74x10 -5 ), suggesting a link to inflammation. Despite an estimated heritability of 32%, we did not identify any proximal or distal genetic variants (eQTLs) associated with adipose tissue ACE2 expression. Our results demonstrate that at-risk individuals have lower background ACE2 levels in this highly relevant tissue. Further studies will be required to establish how this may contribute to increased COVID-19 severity.

Subregional statistical shape modelling identifies lesser trochanter size as a possible risk factor for radiographic hip osteoarthritis, a cross-sectional analysis from the Osteoporotic Fractures in Men Study

OBJECTIVE

Statistical shape modelling (SSM) of hip dual-energy X-ray absorptiometry (DXA) scans has identified relationships between hip shape and radiographic hip OA (rHOA). We aimed to further elucidate shape characteristics related to rHOA by focusing on subregions identified from whole-hip shape models.

METHOD

SSM was applied to hip DXAs obtained in the Osteoporotic Fractures in Men Study. Whole-hip shape modes (HSMs) associated with rHOA were combined to form a composite at-risk-shape. Subsequently, subregional HSMs (cam-type and lesser trochanter modes) were built, and associations with rHOA were examined by logistic regression. Subregional HSMs were further characterised, by examining associations with 3D-HSMs derived from concurrent hip CT scans.

RESULTS

4,098 participants were identified with hip DXAs and radiographs. Composite shapes from whole-hip HSMs revealed that lesser trochanter size and cam-type femoral head are related to rHOA. From sub-regional models, lesser trochanter mode (LTM)1 [OR 0.74; 95%CI 0.63.0.87] and cam-type mode (CTM)3 [OR 1.27; 1.13.1.42] were associated with rHOA, associations being similar to those for whole hip HSMs. 515 MrOS participants had hip DXAs and 3D-HSMs derived from hip CT scans. LTM1 was associated with 3D-HSMs that also represented a larger lesser trochanter [3D-HSM7 (beta (β)-0.23;-0.33,-0.14) and 3D-HSM9 (β0.36; 0.27.0.45)], and CTM3 with 3D-HSMs describing cam morphology [3D-HSM3 (β-0.16;-0.25,-0.07) and 3D-HSM6 (β 0.19; 0.10.0.28)].

CONCLUSION

Subregional SSM of hip DXA scans suggested larger lesser trochanter and cam morphology underlie associations between overall hip shape and rHOA. 3D hip modelling suggests our subregional SSMs represent true anatomical variations in hip shape, warranting further investigation.

MISTIC: A prediction tool to reveal disease-relevant deleterious missense variants

The diffusion of next-generation sequencing technologies has revolutionized research and diagnosis in the field of rare Mendelian disorders, notably via whole-exome sequencing (WES). However, one of the main issues hampering achievement of a diagnosis via WES analyses is the extended list of variants of unknown significance (VUS), mostly composed of missense variants. Hence, improved solutions are needed to address the challenges of identifying potentially deleterious variants and ranking them in a prioritized short list. We present MISTIC (MISsense deleTeriousness predICtor), a new prediction tool based on an original combination of two complementary machine learning algorithms using a soft voting system that integrates 113 missense features, ranging from multi-ethnic minor allele frequencies and evolutionary conservation, to physiochemical and biochemical properties of amino acids. Our approach also uses training sets with a wide spectrum of variant profiles, including both high-confidence positive (deleterious) and negative (benign) variants. Compared to recent state-of-the-art prediction tools in various benchmark tests and independent evaluation scenarios, MISTIC exhibits the best and most consistent performance, notably with the highest AUC value (> 0.95). Importantly, MISTIC maintains its high performance in the specific case of discriminating deleterious variants from benign variants that are rare or population-specific. In a clinical context, MISTIC drastically reduces the list of VUS (<30%) and significantly improves the ranking of "causative" deleterious variants. Pre-computed MISTIC scores for all possible human missense variants are available at http://lbgi.fr/mistic.

Whole-exome sequencing identifies rare variants in STAB2 associated with venous thromboembolic disease

Deep vein thrombosis and pulmonary embolism, collectively defined as venous thromboembolism (VTE), are the third leading cause of cardiovascular death in the United States. Common genetic variants conferring increased varying degrees of VTE risk have been identified by genome-wide association studies (GWAS). Rare mutations in the anticoagulant genes PROC, PROS1 and SERPINC1 result in perinatal lethal thrombosis in homozygotes and markedly increased VTE risk in heterozygotes. However, currently described VTE variants account for an insufficient portion of risk to be routinely used for clinical decision making. To identify new rare VTE risk variants, we performed a whole-exome study of 393 individuals with unprovoked VTE and 6114 controls. This study identified 4 genes harboring an excess number of rare damaging variants in patients with VTE: PROS1, STAB2, PROC, and SERPINC1. At STAB2, 7.8% of VTE cases and 2.4% of controls had a qualifying rare variant. In cell culture, VTE-associated variants of STAB2 had a reduced surface expression compared with reference STAB2. Common variants in STAB2 have been previously associated with plasma von Willebrand factor and coagulation factor VIII levels in GWAS, suggesting that haploinsufficiency of stabilin-2 may increase VTE risk through elevated levels of these procoagulants. In an independent cohort, we found higher von Willebrand factor levels and equivalent propeptide levels in individuals with rare STAB2 variants compared with controls. Taken together, this study demonstrates the utility of gene-based collapsing analyses to identify loci harboring an excess of rare variants with functional connections to a complex thrombotic disease.

Informing retention in longitudinal cohort studies through a social marketing lens: Raine Study Generation 2 participants' perspectives on benefits and barriers to participation

Background

Longitudinal cohort studies have made significant contributions to medical discoveries and provide the impetus for health interventions which reduce the risk of disease. Establishing and maintaining these cohorts is challenging and costly. While some attrition is unavoidable, maintaining a sufficient number of participants ensures that results remain representative and free from bias. Numerous studies have investigated ways to reduce attrition but few studies have sought to understand the experience of participants, and none have examined this through a social marketing framework. This first paper in a two part-series describes participants’ experiences according to: benefits, barriers, motivators and influencers. The second paper uses this understanding to address issues relating to the 4Ps (product, price, place, promotion) of social marketing.

Methods

Participants were recruited from the Raine Study, a pregnancy cohort study that has been running in Western Australia since 1989. Qualitative interviews were conducted with 29 active and inactive participants from the Generation 2 cohort, who were originally enrolled in the Raine Study at birth by their parents (Generation 1). ‘Active’ participants (n = 17) were defined as those who agreed to attend their 27 year follow-up, while ‘inactive’ (n = 12) participants were defined as those who had not attended either of the past two follow-ups (at 22 and 27 years).

Results

There were considerable differences between active and inactive participants, with active participants perceiving far more personal and collective benefits from their participation. Inactive participants described being constrained by structural barriers around work and life, whereas active participants were able to overcome them to attend follow-ups. Inactive participants also described the value of extrinsic incentives which might motivate their attendance, and active participants described the role of their parents as significant influencers in their propensity to remain in the study.

Conclusions

This paper provides rich descriptions of what participation in a long-running study means to participants. Use of a social marketing framework ensured that participants were constructed as ‘human consumers’ who are influenced by individual and broader social systems. Understanding participants in this way means that differentiated strategies can be tailored to enhance retention.

Age-dependent regional retinal nerve fibre changes in SIX1/SIX6 polymorphism

SIX1/SIX6 polymorphism has been shown to be associated with glaucoma. Studies have also found that, in older adults, retinal nerve fibre layer (RNFL) thickness is significantly thinned with each copy of the risk allele in SIX1/SIX6. However, it is not known whether these genetic variants exert their effects in younger individuals. Comparing a healthy young adult with an older adult cohort (mean age 20 vs 63 years), both of Northern European descent, we found that there was no significant RNFL thinning in each copy of the risk alleles in SIX1/SIX6 in the eyes of younger individuals. The older cohort showed an unexpectedly thicker RNFL in the nasal sector with each copy of the risk allele for both the SIX1 (rs10483727) and SIX6 (rs33912345) variants. In the temporal sector, thinner RNFL was found with each copy of the risk allele in rs33912345 with a decrease trend observed in rs10483727. Our results suggest that SIX1/SIX6 gene variants exert their influence later in adult life.

Very low-depth whole-genome sequencing in complex trait association studies

Motivation

Very low-depth sequencing has been proposed as a cost-effective approach to capture low-frequency and rare variation in complex trait association studies. However, a full characterization of the genotype quality and association power for very low-depth sequencing designs is still lacking.

Results

We perform cohort-wide whole-genome sequencing (WGS) at low depth in 1239 individuals (990 at 1× depth and 249 at 4× depth) from an isolated population, and establish a robust pipeline for calling and imputing very low-depth WGS genotypes from standard bioinformatics tools. Using genotyping chip, whole-exome sequencing (75× depth) and high-depth (22×) WGS data in the same samples, we examine in detail the sensitivity of this approach, and show that imputed 1× WGS recapitulates 95.2% of variants found by imputed GWAS with an average minor allele concordance of 97% for common and low-frequency variants. In our study, 1× further allowed the discovery of 140 844 true low-frequency variants with 73% genotype concordance when compared to high-depth WGS data. Finally, using association results for 57 quantitative traits, we show that very low-depth WGS is an efficient alternative to imputed GWAS chip designs, allowing the discovery of up to twice as many true association signals than the classical imputed GWAS design.

Availability and implementation

The HELIC genotype and WGS datasets have been deposited to the European Genome-phenome Archive (https://www.ebi.ac.uk/ega/home): EGAD00010000518; EGAD00010000522; EGAD00010000610; EGAD00001001636, EGAD00001001637. The peakplotter software is available at https://github.com/wtsi-team144/peakplotter, the transformPhenotype app can be downloaded at https://github.com/wtsi-team144/transformPhenotype.

Supplementary information

Supplementary data are available at Bioinformatics online.

Human postprandial responses to food and potential for precision nutrition

Metabolic responses to food influence risk of cardiometabolic disease, but large-scale high-resolution studies are lacking. We recruited n = 1,002 twins and unrelated healthy adults in the United Kingdom to the PREDICT 1 study and assessed postprandial metabolic responses in a clinical setting and at home. We observed large inter-individual variability (as measured by the population coefficient of variation (s.d./mean, %)) in postprandial responses of blood triglyceride (103%), glucose (68%) and insulin (59%) following identical meals. Person-specific factors, such as gut microbiome, had a greater influence (7.1% of variance) than did meal macronutrients (3.6%) for postprandial lipemia, but not for postprandial glycemia (6.0% and 15.4%, respectively); genetic variants had a modest impact on predictions (9.5% for glucose, 0.8% for triglyceride, 0.2% for C-peptide). Findings were independently validated in a US cohort (n = 100 people). We developed a machine-learning model that predicted both triglyceride (r = 0.47) and glycemic (r = 0.77) responses to food intake. These findings may be informative for developing personalized diet strategies. The ClinicalTrials.gov registration identifier is NCT03479866.

Segmented Macular Layer Volumes from Spectral Domain Optical Coherence Tomography in 184 Adult Twins: Associations With Age and Heritability

Purpose

To investigate segmented macular layer volumes from a healthy adult twin cohort (TwinsUK), exploring changes with age and heritability.

Methods

Macular spectral domain optical coherence tomography images were acquired from monozygotic (MZ) and dizygotic (DZ) twins in a cross-sectional study. The following layer volumes were derived for circles of 3 and 6 mm diameter around the foveal center, using automated segmentation software: retinal nerve fiber layer (RNFL), ganglion cell–inner plexiform layer (GCIPL), inner nuclear layer (INL), outer plexiform layer (OPL), outer nuclear layer (ONL), photoreceptors (PR), retinal pigment epithelium (RPE), and total retinal volume (TRV). Correlation coefficients (intereye; age; intrapair for MZ and DZ pairs) were quantified; heritability was estimated using structural equation modeling.

Results

Scans from 184 participants were included. Intereye correlation was highest for TRV and GCIPL. Negative correlations with age (for 3- or 6-mm areas, or both) were observed for TRV, RNFL, GCIPL, and INL. Positive correlations were observed for PR, RPE, and OPL. For all layers, intrapair correlation was greater for MZ than DZ pairs. Heritability estimates were highest (>80%) for TRV and GCIPL volume, and lowest for RPE volume.

Conclusions

Although TRV was negatively correlated with age, all layers did not show negative correlation. Some inner layers thinned with age, whereas some outer volumes increased (not the ONL). Reduced RPE phagocytic function with age and remodeling in the OPL could be contributing factors. Heritability estimates were highest for inner retinal layers (particularly GCIPL), and lowest for RPE volume.

Reference exome data for Australian Aboriginal populations to support health-based research

Whole exome sequencing (WES) is a popular and successful technology which is widely used in both research and clinical settings. However, there is a paucity of reference data for Aboriginal Australians to underpin the translation of health-based genomic research. Here we provide a catalogue of variants called after sequencing the exomes of 50 Aboriginal individuals from the Northern Territory (NT) of Australia and compare these to 72 previously published exomes from a Western Australian (WA) population of Martu origin. Sequence data for both NT and WA samples were processed using an ‘intersect-then-combine’ (ITC) approach, using GATK and SAMtools to call variants. A total of 289,829 variants were identified in at least one individual in the NT cohort and 248,374 variants in at least one individual in the WA cohort. Of these, 166,719 variants were present in both cohorts, whilst 123,110 variants were private to the NT cohort and 81,655 were private to the WA cohort. Our data set provides a useful reference point for genomic studies on Aboriginal Australians.

Measurement(s)	Aboriginal Australian • DNA • sequence feature annotation
Technology Type(s)	Whole Exome Sequencing • DNA sequencing • sequence annotation
Factor Type(s)	ancestry • sex • age
Sample Characteristic - Organism	Homo sapiens
Sample Characteristic - Location	Northern Territory

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.12040638

Elucidating the role of maternal environmental exposures on offspring health and disease using two-sample Mendelian randomization

BACKGROUND

There is considerable interest in estimating the causal effect of a range of maternal environmental exposures on offspring health-related outcomes. Previous attempts to do this using Mendelian randomization methodologies have been hampered by the paucity of epidemiological cohorts with large numbers of genotyped mother-offspring pairs.

METHODS

We describe a new statistical model that we have created which can be used to estimate the effect of maternal genotypes on offspring outcomes conditional on offspring genotype, using both individual-level and summary-results data, even when the extent of sample overlap is unknown.

RESULTS

We describe how the estimates obtained from our method can subsequently be used in large-scale two-sample Mendelian randomization studies to investigate the causal effect of maternal environmental exposures on offspring outcomes. This includes studies that aim to assess the causal effect of in utero exposures related to fetal growth restriction on future risk of disease in offspring. We illustrate our framework using examples related to offspring birthweight and cardiometabolic disease, although the general principles we espouse are relevant for many other offspring phenotypes.

CONCLUSIONS

We advocate for the establishment of large-scale international genetics consortia that are focused on the identification of maternal genetic effects and committed to the public sharing of genome-wide summary-results data from such efforts. This information will facilitate the application of powerful two-sample Mendelian randomization studies of maternal exposures and offspring outcomes.

Self-reported hearing loss questions provide a good measure for genetic studies: a polygenic risk score analysis from UK Biobank

Age-related hearing impairment (ARHI) is very common in older adults and has major impact on quality of life. The heritability of ARHI has been estimated to be around 50%. The present study aimed to estimate heritability and environmental contributions to liability of ARHI and the extent to which a polygenic risk score (PRS) derived from a recent genome-wide association study of questionnaire items regarding hearing loss using the UK Biobank is predictive of hearing loss in other samples. We examined (1) a sample from TwinsUK who have had hearing ability measured by pure-tone audiogram and the speech-to-noise ratio test as well as questionnaire measures that are comparable with the UK Biobank questionnaire items and (2) European and non-European samples from the UK Biobank which were not part of the original GWAS. Results indicated that the questionnaire items were over 50% heritable in TwinsUK and comparable with the objective hearing measures. In addition, we found very high genetic correlation (0.30-0.84) between the questionnaire responses and objective hearing measures in the TwinsUK sample. Finally, PRS computed from weighted UK Biobank GWAS results were predictive of both questionnaire and objective measures of hearing loss in the TwinsUK sample, as well as questionnaire-measured hearing loss in Europeans but not non-European subpopulations. These results demonstrate the utility of questionnaire-based methods in genetic association studies of hearing loss in adults and highlight the differences in genetic predisposition to ARHI by ethnic background.

Dysregulated Antibody, Natural Killer Cell and Immune Mediator Profiles in Autoimmune Thyroid Diseases

The pathogenesis of autoimmune thyroid diseases (AITD) is poorly understood and the association between different immune features and the germline variants involved in AITD are yet unclear. We previously observed systemic depletion of IgG core fucosylation and antennary α1,2 fucosylation in peripheral blood mononuclear cells in AITD, correlated with anti-thyroid peroxidase antibody (TPOAb) levels. Fucose depletion is known to potentiate strong antibody-mediated NK cell activation and enhanced target antigen-expressing cell killing. In autoimmunity, this may translate to autoantibody-mediated immune cell recruitment and attack of self-antigen expressing normal tissues. Hence, we investigated the crosstalk between immune cell traits, secreted proteins, genetic variants and the glycosylation patterns of serum IgG, in a multi-omic and cross-sectional study of 622 individuals from the TwinsUK cohort, 172 of whom were diagnosed with AITD. We observed associations between two genetic variants (rs505922 and rs687621), AITD status, the secretion of Desmoglein-2 protein, and the profile of two IgG N-glycan traits in AITD, but further studies need to be performed to better understand their crosstalk in AITD. On the other side, enhanced afucosylated IgG was positively associated with activatory CD335- CD314+ CD158b+ NK cell subsets. Increased levels of the apoptosis and inflammation markers Caspase-2 and Interleukin-1α positively associated with AITD. Two genetic variants associated with AITD, rs1521 and rs3094228, were also associated with altered expression of the thyrocyte-expressed ligands known to recognize the NK cell immunoreceptors CD314 and CD158b. Our analyses reveal a combination of heightened Fc-active IgG antibodies, effector cells, cytokines and apoptotic signals in AITD, and AITD genetic variants associated with altered expression of thyrocyte-expressed ligands to NK cell immunoreceptors. Together, TPOAb responses, dysregulated immune features, germline variants associated with immunoactivity profiles, are consistent with a positive autoreactive antibody-dependent NK cell-mediated immune response likely drawn to the thyroid gland in AITD.

Spector T, S. Small K, Steves CJ. Introducing ExHiBITT – Exploring Host microBIome inTeractions in Twins –, a colon multiomic cohort study

Background: The colon is populated by approximately 1012 microorganisms, but the relationships between this microbiome and the host health status are still not completely understood. Here, our objective is to present the cohort characteristics of ExHiBITT – Exploring Host microBIome inTeractions in Twins – including i) biomedical phenotypes, ii) environmental factors and ii) colonoscopic findings. Methods: Participants from the TwinsUK cohort were recruited to study the interactions between the microbiome and host adaptive immunity. In total, 205 monozygotic twins were recruited from the wider TwinsUK cohort. They completed health questionnaires, and provided saliva, blood, colon biopsies from three different locations, caecal fluid, and two faecal samples. Results: A significant proportion of this apparently normal cohort had colonic polyps (28%), which are of interest as potential precursors of colorectal cancer, and, as expected, the number of polyps found was significantly correlated with BMI and age. Hitherto undiagnosed diverticulosis was also not infrequently found during colonoscopy (26%) and was associated with changes in Hybrid Th1-17 cells in the colon. Twin proband co-occurrence rate for diverticulosis (82%) was much higher than for polyps (42%). Familial factors affecting morphology or tolerance may contribute to the ease of endoscopy, as both the time to reach the caecum and pain perceived were highly concordant (proband concordance: 85% and 56%, respectively). Conclusions: We found the expected positive relationship between BMI and colonoscopic anomalies such as diverticular disease and polyps in the whole population, but within twin pairs this association was reversed. This suggests that familial factors confound these associations. Host and microbial next generation sequencing and metabolomics of the samples collected are planned in this cohort.

Glycosylation of immunoglobulin G is regulated by a large network of genes pleiotropic with inflammatory diseases

Effector functions of immunoglobulin G (IgG) are regulated by the composition of a glycan moiety, thus affecting activity of the immune system. Aberrant glycosylation of IgG has been observed in many diseases, but little is understood about the underlying mechanisms. We performed a genome-wide association study of IgG N-glycosylation (N = 8090) and, using a data-driven network approach, suggested how associated loci form a functional network. We confirmed in vitro that knockdown of IKZF1 decreases the expression of fucosyltransferase FUT8, resulting in increased levels of fucosylated glycans, and suggest that RUNX1 and RUNX3, together with SMARCB1, regulate expression of glycosyltransferase MGAT3. We also show that variants affecting the expression of genes involved in the regulation of glycoenzymes colocalize with variants affecting risk for inflammatory diseases. This study provides new evidence that variation in key transcription factors coupled with regulatory variation in glycogenes modifies IgG glycosylation and has influence on inflammatory diseases.

An unsupervised learning approach to identify novel signatures of health and disease from multimodal data

BACKGROUND

Modern medicine is rapidly moving towards a data-driven paradigm based on comprehensive multimodal health assessments. Integrated analysis of data from different modalities has the potential of uncovering novel biomarkers and disease signatures.

METHODS

We collected 1385 data features from diverse modalities, including metabolome, microbiome, genetics, and advanced imaging, from 1253 individuals and from a longitudinal validation cohort of 1083 individuals. We utilized a combination of unsupervised machine learning methods to identify multimodal biomarker signatures of health and disease risk.

RESULTS

Our method identified a set of cardiometabolic biomarkers that goes beyond standard clinical biomarkers. Stratification of individuals based on the signatures of these biomarkers identified distinct subsets of individuals with similar health statuses. Subset membership was a better predictor for diabetes than established clinical biomarkers such as glucose, insulin resistance, and body mass index. The novel biomarkers in the diabetes signature included 1-stearoyl-2-dihomo-linolenoyl-GPC and 1-(1-enyl-palmitoyl)-2-oleoyl-GPC. Another metabolite, cinnamoylglycine, was identified as a potential biomarker for both gut microbiome health and lean mass percentage. We identified potential early signatures for hypertension and a poor metabolic health outcome. Additionally, we found novel associations between a uremic toxin, p-cresol sulfate, and the abundance of the microbiome genera Intestinimonas and an unclassified genus in the Erysipelotrichaceae family.

CONCLUSIONS

Our methodology and results demonstrate the potential of multimodal data integration, from the identification of novel biomarker signatures to a data-driven stratification of individuals into disease subtypes and stages-an essential step towards personalized, preventative health risk assessment.

Identifying causal variants and genes using functional genomics in specialized cell types and contexts

A central goal in human genetics is the identification of variants and genes that influence the risk of polygenic diseases. In the past decade, genome-wide association studies (GWAS) have identified tens of thousands of genetic loci associated with various diseases. Since the majority of such loci lie within non-coding regions and have many candidate variants in linkage disequilibrium, it has been challenging to accurately identify specific causal variants and genes. To aid in their discovery a variety of statistical and experimental approaches have been developed. These approaches often borrow information from functional genomics assays such as ATAC-seq, ChIP-seq and RNA-seq to annotate functional variants and identify regulatory relationships between variants and genes. While such approaches are powerful, given the diversity of cell types and environments, it is paramount to select disease-relevant contexts for follow-up analyses. In this review, we discuss the latest developments, challenges, and best practices for determining the causal mechanisms of polygenic disease risk variants with functional genomics data from specialized cell types.

The GenomeAsia 100K Project enables genetic discoveries across Asia

The underrepresentation of non-Europeans in human genetic studies so far has limited the diversity of individuals in genomic datasets and led to reduced medical relevance for a large proportion of the world’s population. Population-specific reference genome datasets as well as genome-wide association studies in diverse populations are needed to address this issue. Here we describe the pilot phase of the GenomeAsia 100K Project. This includes a whole-genome sequencing reference dataset from 1,739 individuals of 219 population groups and 64 countries across Asia. We catalogue genetic variation, population structure, disease associations and founder effects. We also explore the use of this dataset in imputation, to facilitate genetic studies in populations across Asia and worldwide.

Using whole-genome sequencing data from 1,739 individuals, the GenomeAsia 100K Project catalogues genetic variation, population structure and disease associations to facilitate genetic studies in Asian populations and increase representation in genetics studies worldwide.

Heritability of skewed X-inactivation in female twins is tissue-specific and associated with age

Female somatic X-chromosome inactivation (XCI) balances the X-linked transcriptional dosages between the sexes. Skewed XCI toward one parental X has been observed in several complex human traits, but the extent to which genetics and environment influence skewed XCI is largely unexplored. To address this, we quantify XCI-skew in multiple tissues and immune cell types in a twin cohort. Within an individual, XCI-skew differs between blood, fat and skin tissue, but is shared across immune cell types. XCI skew increases with age in blood, but not other tissues, and is associated with smoking. XCI-skew is increased in twins with Rheumatoid Arthritis compared to unaffected identical co-twins. XCI-skew is heritable in blood of females >55 years old (h² = 0.34), but not in younger individuals or other tissues. This results in a Gene x Age interaction that shifts the functional dosage of all X-linked heterozygous loci in a tissue-restricted manner.

Skewing of X chromosome inactivation (XCI) occurs when the silencing of one parental X chromosome is non-random. Here, Zito et al. report XCI patterns in lymphoblastoid cell lines, blood, subcutaneous adipose tissue samples and skin samples of monozygotic and dizygotic twins and find XCI skew to associate with tissue and age.

Systems genetics applications in metabolism research

The common forms of metabolic diseases are highly complex, involving hundreds of genes, environmental and lifestyle factors, age-related changes, sex differences and gut-microbiome interactions. Systems genetics is a population-based approach to address this complexity. In contrast to commonly used 'reductionist' approaches, such as gain or loss of function, that examine one element at a time, systems genetics uses high-throughput 'omics' technologies to quantitatively assess the many molecular differences among individuals in a population and then to relate these to physiologic functions or disease states. Unlike genome-wide association studies, systems genetics seeks to go beyond the identification of disease-causing genes to understand higher-order interactions at the molecular level. The purpose of this review is to introduce the systems genetics applications in the areas of metabolic and cardiovascular disease. Here, we explain how large clinical and omics-level data and databases from both human and animal populations are available to help researchers place genes in the context of pathways and networks and formulate hypotheses that can then be experimentally examined. We provide lists of such databases and examples of the integration of reductionist and systems genetics data. Among the important applications emerging is the development of improved nutritional and pharmacological strategies to address the rise of metabolic diseases.

GWAS Identifies 44 Independent Associated Genomic Loci for Self-Reported Adult Hearing Difficulty in UK Biobank

Age-related hearing impairment (ARHI) is the most common sensory impairment in the aging population; a third of individuals are affected by disabling hearing loss by the age of 65. It causes social isolation and depression and has recently been identified as a risk factor for dementia. The genetic risk factors and underlying pathology of ARHI are largely unknown, meaning that targets for new therapies remain elusive, yet heritability estimates range between 35% and 55%. We performed genome-wide association studies (GWASs) for two self-reported hearing phenotypes, using more than 250,000 UK Biobank (UKBB) volunteers aged between 40 and 69 years. Forty-four independent genome-wide significant loci (p < 5E-08) were identified, considerably increasing the number of established trait loci. Thirty-four loci are novel associations with hearing loss of any form, and only one of the ten known hearing loci has a previously reported association with an ARHI-related trait. Gene sets from these loci are enriched in auditory processes such as synaptic activities, nervous system processes, inner ear morphology, and cognition, while genetic correlation analysis revealed strong positive correlations with multiple personality and psychological traits for the first time. Immunohistochemistry for protein localization in adult mouse cochlea implicate metabolic, sensory, and neuronal functions for NID2, CLRN2, and ARHGEF28. These results provide insight into the genetic landscape underlying ARHI, opening up novel therapeutic targets for further investigation. In a wider context, our study also highlights the viability of using self-report phenotypes for genetic discovery in very large samples when deep phenotyping is unavailable.

TwinsUK: The UK Adult Twin Registry Update

TwinsUK is the largest cohort of community-dwelling adult twins in the UK. The registry comprises over 14,000 volunteer twins (14,838 including mixed, single and triplets); it is predominantly female (82%) and middle-aged (mean age 59). In addition, over 1800 parents and siblings of twins are registered volunteers. During the last 27 years, TwinsUK has collected numerous questionnaire responses, physical/cognitive measures and biological measures on over 8500 subjects. Data were collected alongside four comprehensive phenotyping clinical visits to the Department of Twin Research and Genetic Epidemiology, King’s College London. Such collection methods have resulted in very detailed longitudinal clinical, biochemical, behavioral, dietary and socioeconomic cohort characterization; it provides a multidisciplinary platform for the study of complex disease during the adult life course, including the process of healthy aging. The major strength of TwinsUK is the availability of several ‘omic’ technologies for a range of sample types from participants, which includes genomewide scans of single-nucleotide variants, next-generation sequencing, metabolomic profiles, microbiomics, exome sequencing, epigenetic markers, gene expression arrays, RNA sequencing and telomere length measures. TwinsUK facilitates and actively encourages sharing the ‘TwinsUK’ resource with the scientific community — interested researchers may request data via the TwinsUK website (http://twinsuk.ac.uk/resources-for-researchers/access-our-data/) for their own use or future collaboration with the study team. In addition, further cohort data collection is planned via the Wellcome Open Research gateway (https://wellcomeopenresearch.org/gateways). The current article presents an up-to-date report on the application of technological advances, new study procedures in the cohort and future direction of TwinsUK.

Twin Registries Moving Forward and Meeting the Future: A Review

Twin registries have developed as a valuable resource for the study of many aspects of disease and society over the years in many different countries. A number of these registries include large numbers of twins with data collected at varying information levels for twin cohorts over the past several decades. More recent expansion of twin datasets has allowed for the collection of genetic data, together with many other levels of 'omic' information along with multiple demographic, physiological, health outcomes and other measures typically used in epidemiologic research. Other twin data sources outside these registries reflect research interests in particular aspects of disease or specific phenotypic assessment. Twin registries have the potential to play a key role in many aspects of the artificial intelligence/machine learning-driven projects of the future and will continue to keep adapting to the changing research landscape.

Dissecting the role of the gut microbiota and diet on visceral fat mass accumulation

Both gut microbiota and diet have been shown to impact visceral fat mass (VFM), a major risk factor for cardiometabolic disease, but their relative contribution has not been well characterised. We aimed to estimate and separate the effect of gut microbiota composition from that of nutrient intake on VFM in 1760 older female twins. Through pairwise association analyses, we identified 93 operational taxonomic units (OTUs) and 10 nutrients independently linked to VFM (FDR < 5%). Conditional analyses revealed that the majority (87%) of the 93 VFM-associated OTUs remained significantly associated with VFM irrespective of nutrient intake correction. In contrast, we observed that the effect of fibre, magnesium, biotin and vitamin E on VFM was partially mediated by OTUs. Moreover, we estimated that OTUs were more accurate predictors of VFM than nutrients and accounted for a larger percentage of its variance. Our results suggest that while the role of certain nutrients on VFM appears to depend on gut microbiota composition, specific gut microbes may affect host adiposity regardless of dietary intake. The findings imply that the gut microbiota may have a greater contribution towards shaping host VFM than diet alone. Thus, microbial-based therapy should be prioritised for VFM reduction in overweight and obese subjects.

Fine-Mapping of the Human Blood Plasma N-Glycome onto Its Proteome

Most human proteins are glycosylated. Attachment of complex oligosaccharides to the polypeptide part of these proteins is an integral part of their structure and function and plays a central role in many complex disorders. One approach towards deciphering this human glycan code is to study natural variation in experimentally well characterized samples and cohorts. High-throughput capable large-scale methods that allow for the comprehensive determination of blood circulating proteins and their glycans have been recently developed, but so far, no study has investigated the link between both traits. Here we map for the first time the blood plasma proteome to its matching N-glycome by correlating the levels of 1116 blood circulating proteins with 113 N-glycan traits, determined in 344 samples from individuals of Arab, South-Asian, and Filipino descent, and then replicate our findings in 46 subjects of European ancestry. We report protein-specific N-glycosylation patterns, including a correlation of core fucosylated structures with immunoglobulin G (IgG) levels, and of trisialylated, trigalactosylated, and triantennary structures with heparin cofactor 2 (SERPIND2). Our study reveals a detailed picture of protein N-glycosylation and suggests new avenues for the investigation of its role and function in the associated complex disorders.

Defining the genetic control of human blood plasma N-glycome using genome-wide association study

Glycosylation is a common post-translational modification of proteins. Glycosylation is associated with a number of human diseases. Defining genetic factors altering glycosylation may provide a basis for novel approaches to diagnostic and pharmaceutical applications. Here we report a genome-wide association study of the human blood plasma N-glycome composition in up to 3811 people measured by Ultra Performance Liquid Chromatography (UPLC) technology. Starting with the 36 original traits measured by UPLC, we computed an additional 77 derived traits leading to a total of 113 glycan traits. We studied associations between these traits and genetic polymorphisms located on human autosomes. We discovered and replicated 12 loci. This allowed us to demonstrate an overlap in genetic control between total plasma protein and IgG glycosylation. The majority of revealed loci contained genes that encode enzymes directly involved in glycosylation (FUT3/FUT6, FUT8, B3GAT1, ST6GAL1, B4GALT1, ST3GAL4, MGAT3 and MGAT5) and a known regulator of plasma protein fucosylation (HNF1A). However, we also found loci that could possibly reflect other more complex aspects of glycosylation process. Functional genomic annotation suggested the role of several genes including DERL3, CHCHD10, TMEM121, IGH and IKZF1. The hypotheses we generated may serve as a starting point for further functional studies in this research area.

Rare variant association testing for multicategory phenotype

Genetic association studies have provided new insights into the genetic variability of human complex traits with a focus mainly on continuous or binary traits. Methods have been proposed to take into account disease heterogeneity between subgroups of patients when studying common variants but none was specifically designed for rare variants. Because rare variants are expected to have stronger effects and to be more heterogeneously distributed among cases than common ones, subgroup analyses might be particularly attractive in this context. To address this issue, we propose an extension of burden tests by using a multinomial regression model, which enables association tests between rare variants and multicategory phenotypes. We evaluated the type I error and the power of two burden tests, CAST and WSS, by simulating data under different scenarios. In the case of genetic heterogeneity between case subgroups, we showed an advantage of multinomial regression over logistic regression, which considers all the cases against the controls. We replicated these results on real data from Moyamoya disease where the burden tests performed better when cases were stratified according to age-of-onset. We implemented the functions for association tests in the R package "Ravages" available on Github.

Self-reported Anxiety Sensitivity Index in the TwinsUK cohort

The anxiety sensitivity (AS) construct has received considerable attention in anxiety research and is considered to be a cognitive vulnerability factor for the study of anxiety related disorders. The Anxiety Sensitivity Index (ASI) is the most widely used instrument for the study of AS. The present Data Note provides an overview of all the 16-item ASI questionnaires filled and returned by the twins in the TwinsUK registry. This work does not provide any multidimensional or factor structure analysis of the responses provided. TwinsUK registry encompasses a wide range of clinical and self-reported data that can be used as confounding factors in the study of cognitive and mental health.

Genome-wide meta-analysis identifies genetic locus on chromosome 9 associated with Modic changes

BACKGROUND

Low back pain (LBP) is a common disabling condition. Lumbar disc degeneration (LDD) may be a contributing factor for LBP. Modic change (MC), a distinct phenotype of LDD, is presented as a pathological bone marrow signal change adjacent to vertebral endplate on MRI. It is strongly associated with LBP and has heritability around 30%. Our objective was to identify genetic loci associated with MC using a genome-wide meta-analysis.

METHODS

Presence of MC was evaluated in lumbar MRI in the Northern Finland Birth Cohort 1966 (n=1182) and TwinsUK (n=647). Genome-wide association analyses were carried out using linear regression model. Inverse-variance weighting approach was used in the meta-analysis.

RESULTS

A locus associated with MC (p<5e-8) was found on chromosome 9 with the lead SNP rs1934268 in an intron of the PTPRD gene. It is located in the binding region of BCL11A, SPI1 and PBX3 transcription factors. The SNP was nominally associated with LBP in TwinsUK (p=0.001) but not associated in the UK Biobank (p=0.914). Suggestive signals (p<1e-5) were identified near XKR4, SCIN, MGMT, DLG2, ZNF184 and OPRK1.

CONCLUSION

PTPRD is a novel candidate gene for MC that may act via the development of cartilage or nervous system; further work is needed to define the mechanisms underlying the pathways leading to development of MC. This is the first genome-wide meta-analysis of MC, and the results pave the way for further studies on the genetic factors underlying the various features of spine degeneration and LBP.

The analysis of causal relationships between blood lipid levels and BMD

Bone mineral density (BMD) and lipid levels are two of the most extensively studied risk factors for common diseases of aging, such as cardiovascular disease (CVD) and osteoporosis (OP). These two risk factors are also correlated with each other, but little is known about the molecular mechanisms behind this correlation. Recent studies revealed that circulating levels of several metabolites involved in the biosynthesis of androsterone correlate significantly with BMD and have the capacity to affect cholesterol and lipids levels. A main aim of the present study was to investigate the hypothesis that androsterone-related metabolites could provide a link between CVD and OP, as a common cause of lipid levels and BMD. The present study employed data from the NIHR BRC TwinsUK BioResource, comprising 1909 and 1994 monozygotic and dizygotic twin pairs, respectively, to address the causal relationships among BMD and lipids, and their associated metabolites, using reciprocal causation twin modelling, as well as Mendelian randomization (MR) using large publicly-available GWAS datasets on lipids and BMD, in conjunction with TwinsUK metabolite data. While results involving the twin modelling and MR analyses with metabolites were unable to establish a causal link between metabolite levels and either lipids or BMD, MR analyses of BMD and lipids suggest that lipid levels have a causal impact on BMD, which is consistent with findings from clinical trials of lipid-lowering drugs, which have also increased BMD.

Profound Perturbation of the Metabolome in Obesity Is Associated with Health Risk

Obesity is a heterogeneous phenotype that is crudely measured by body mass index (BMI). There is a need for a more precise yet portable method of phenotyping and categorizing risk in large numbers of people with obesity to advance clinical care and drug development. Here, we used non-targeted metabolomics and whole-genome sequencing to identify metabolic and genetic signatures of obesity. We find that obesity results in profound perturbation of the metabolome; nearly a third of the assayed metabolites associated with changes in BMI. A metabolome signature identifies the healthy obese and lean individuals with abnormal metabolomes-these groups differ in health outcomes and underlying genetic risk. Specifically, an abnormal metabolome associated with a 2- to 5-fold increase in cardiovascular events when comparing individuals who were matched for BMI but had opposing metabolome signatures. Because metabolome profiling identifies clinically meaningful heterogeneity in obesity, this approach could help select patients for clinical trials.

An overview of the TwinsUK cohort’s anxiety and depression assessment, using the self-reported Hospital Anxiety and Depression Scale

Anxiety and depression are common mental health conditions that frequently coexist. The Hospital Anxiety and Depression (HADS) questionnaire, is a widely used questionnaire, both in the clinic and in research. HADS provides an assessment of anxiety and depression based on non-physical symptoms, hence it can be used in patients with other physical pain and symptoms. This resource provides an overview of all the HADS questionnaires completed by the twins in the TwinsUK registry. This includes both monozygotic and dizygotic twins, over a wide range of age and socioeconomic backgrounds. TwinsUK holds a wide range of phenotype and genotyping data that can be used as confounding factors in the study of mental health.

Genome-wide association study of offspring birth weight in 86 577 women identifies five novel loci and highlights maternal genetic effects that are independent of fetal genetics

Genome-wide association studies of birth weight have focused on fetal genetics, whereas relatively little is known about the role of maternal genetic variation. We aimed to identify maternal genetic variants associated with birth weight that could highlight potentially relevant maternal determinants of fetal growth. We meta-analysed data on up to 8.7 million SNPs in up to 86 577 women of European descent from the Early Growth Genetics (EGG) Consortium and the UK Biobank. We used structural equation modelling (SEM) and analyses of mother–child pairs to quantify the separate maternal and fetal genetic effects. Maternal SNPs at 10 loci (MTNR1B, HMGA2, SH2B3, KCNAB1, L3MBTL3, GCK, EBF1, TCF7L2, ACTL9, CYP3A7) were associated with offspring birth weight at P < 5 × 10⁻⁸. In SEM analyses, at least 7 of the 10 associations were consistent with effects of the maternal genotype acting via the intrauterine environment, rather than via effects of shared alleles with the fetus. Variants, or correlated proxies, at many of the loci had been previously associated with adult traits, including fasting glucose (MTNR1B, GCK and TCF7L2) and sex hormone levels (CYP3A7), and one (EBF1) with gestational duration. The identified associations indicate that genetic effects on maternal glucose, cytochrome P450 activity and gestational duration, and potentially on maternal blood pressure and immune function, are relevant for fetal growth. Further characterization of these associations in mechanistic and causal analyses will enhance understanding of the potentially modifiable maternal determinants of fetal growth, with the goal of reducing the morbidity and mortality associated with low and high birth weights.

Definitive Zygosity Scores in the Peas in the Pod Questionnaire is a Sensitive and Accurate Assessment of the Zygosity of Adult Twins

Twin researchers face the challenge of accurately determining the zygosity of twins for research. As part of the annual questionnaire between 1999 and 2006, 8,307 twins from the TwinsUK registry were asked to complete five questions (independently from their co-twin) to ascertain their self-perceived zygosity during childhood on up to five separate occasions. This questionnaire is known as the 'peas in the pod' questionnaire (PPQ), but there is little evidence of its validation. Answers were scored and classified as monozygotic (MZ), dizygotic (DZ), or unknown zygosity (UZ) and were compared with 4,484 twins with genotyping data who had not been selected for zygosity. Of these, 3,859 individuals (46.5% of those who had a zygosity from PPQ) had zygosity classified by both the PPQ and genotyping. Of the 708 individual twins whose answers meant that they were consistently classed as MZ in the PPQ, 683 (96.5%) were MZ within the genotype data. Of the 945 individual twins consistently classed as DZ within questionnaire, 936 (99.0%) were DZ in the genotype data. Where both twins scored MZ consistently across multiple questionnaires, 99.6% were MZ on genotyping, 99.7% were DZ on genotyping if both twins consistently scored DZ. However, for the initial questionnaire, 88.6% of those scoring as MZ were genotypically MZ and 98.7% DZ. For twin pairs where both scored UZ, 94.7% were DZ. Using the PPQ on a single occasion provided a definitive classification of whether the twin was MZ or DZ with an overall accuracy of 86.9%, increasing to 97.9% when there was a consistent classification of zygosity across multiple questionnaires. This study has shown that the PPQ questionnaire is an excellent proxy indicator of zygosity in the absence of genotyping information.

An Investigation Into Physical Frailty as a Link Between the Gut Microbiome and Cognitive Health

The preservation of cognitive abilities with aging is a priority both for individuals and nations given the aging populations of many countries. Recently the gut microbiome has been identified as a new territory to explore in relation to cognition. Experiments using rodents have identified a link between the gut microbiome and cognitive function, particularly that low microbial diversity leads to poor cognition function. Similar studies in humans could identify novel targets to encourage healthy cognition in an aging population. Here, we investigate the association of gut microbiota and cognitive function in a human cohort considering the influence of physical frailty. We analyzed 16S rRNA gene sequence data, derived from fecal samples obtained from 1,551 individuals over the age of 40. Cognitive data was collected using four cognitive tests: verbal fluency (n = 1,368), Deary-Liewald Reaction Time Test (DLRT; n = 873), Mini Mental State Examination (recall; n = 1,374) and Paired Associates Learning from the Cambridge Neuropsychological Test Automated Battery (CANTAB-PAL; n = 405). We use mixed effects models to identify associations with alpha diversity, operational taxonomic units (OTUs) and taxa and performed further analyses adjusting for physical frailty. We then repeated the analyses in a subset of individuals with dietary data, also excluding those using medications shown to influence gut microbiome composition. DLRT and verbal fluency were negatively associated with alpha diversity of the gut microbiota (False-Discovery Rate, FDR, p < 0.05). However, when considering frailty as a covariate, only associations between the DLRT and diversity measures remained. Repeating analyses excluding Proton pump inhibitor (PPI) and antibiotic users and accounting for diet, we similarly observe significant negative associations between the DLRT and alpha diversity measures and a further negative association between DLRT and the abundance of the order Burkholderiales that remains significant after adjusting for host frailty. This highlights the importance of considering concurrent differences in physical health in studies of cognitive performance and suggests that physical health has a relatively larger association with the gut microbiome. However, the frailty independent cognitive-gut microbiota associations that were observed might represent important targets for further research, with potential for use in diagnostic surveillance in cognitive aging and interventions to improve vitality.

Endplate Defect Is Heritable, Associated With Low Back Pain and Triggers Intervertebral Disc Degeneration: A Longitudinal Study From TwinsUK

STUDY DESIGN

Longitudinal study of spine magnetic resonance imaging (MRI) in a large-scale population-based study.

OBJECTIVE

To determine the order of appearance of degenerative change in vertebral bodies and intervertebral discs. We also sought to define the influence of endplate defect on low back pain (LBP) and to determine whether there is a genetic influence on endplate defect.

SUMMARY OF BACKGROUND DATA

Endplate defect is a magnetic resonance imaging trait, found to be associated with intervertebral disc degeneration. There is a lack of understanding regarding the mechanism underlying lumbar disc degeneration (LDD). Recent attention has shifted to vertebral endplate defects and their role in spine degeneration pathology.

METHODS

Individuals from the TwinsUK spine study having longitudinal T2-weighted lumbar MR scans at baseline (n = 996) and a decade later (n = 438) were included. LDD, vertebral endplate defect by calculating a total endplate score, and Modic change (MC) were assessed using standard techniques. Mixed-effects models were used to determine the association between the features of spine pathology, adjusted for covariates. Endplate defect heritability was estimated using variance component analysis.

RESULTS

Significant association was found between endplate defect, LDD, MRI features of LDD and MC was observed. Endplate defect was associated with severe disabling LBP (P ≤ 0.013) in multivariate analysis. An association between disc degeneration (DD) at baseline and MC at follow-up was shown at upper lumbar levels. Total endplate score was heritable with estimated additive genetic component A = 55.3% (95% CI 43.0-65.4).

CONCLUSION

Endplate defect, LDD, and MC are all independent risk factors for episodes of severe and disabling LBP. Longitudinal analysis showed DD is followed by MC. Endplate defect has significant heritability of 55%. However, whether endplate defect triggers DD or these pathological changes occur concurrently could not be conclusively determined.

LEVEL OF EVIDENCE

2.

Social and population health science approaches to understand the human microbiome

The microbiome is now considered our 'second genome' with potentially comparable importance to the genome in determining human health. There is, however, a relatively limited understanding of the broader environmental factors, particularly social conditions, that shape variation in human microbial communities. Fulfilling the promise of microbiome research - particularly the microbiome's potential for modification - will require collaboration between biologists and social and population scientists. For life scientists, the plasticity and adaptiveness of the microbiome calls for an agenda to understand the sensitivity of the microbiome to broader social environments already known to be powerful predictors of morbidity and mortality. For social and population scientists, attention to the microbiome may help answer nagging questions about the underlying biological mechanisms that link social conditions to health. We outline key substantive and methodological advances that can be made if collaborations between social and population health scientists and life scientists are strategically pursued.

A GWAS meta-analysis from 5 population-based cohorts implicates ion channel genes in the pathogenesis of irritable bowel syndrome

BACKGROUND

Irritable bowel syndrome (IBS) shows genetic predisposition, however, large-scale, powered gene mapping studies are lacking. We sought to exploit existing genetic (genotype) and epidemiological (questionnaire) data from a series of population-based cohorts for IBS genome-wide association studies (GWAS) and their meta-analysis.

METHODS

Based on questionnaire data compatible with Rome III Criteria, we identified a total of 1335 IBS cases and 9768 asymptomatic individuals from 5 independent European genotyped cohorts. Individual GWAS were carried out with sex-adjusted logistic regression under an additive model, followed by meta-analysis using the inverse variance method. Functional annotation of significant results was obtained via a computational pipeline exploiting ontology and interaction networks, and tissue-specific and gene set enrichment analyses.

KEY RESULTS

Suggestive GWAS signals (P ≤ 5.0 × 10^-6 ) were detected for 7 genomic regions, harboring 64 gene candidates to affect IBS risk via functional or expression changes. Functional annotation of this gene set convincingly (best FDR-corrected P = 3.1 × 10^-10 ) highlighted regulation of ion channel activity as the most plausible pathway affecting IBS risk.

CONCLUSION & INFERENCES

Our results confirm the feasibility of population-based studies for gene-discovery efforts in IBS, identify risk genes and loci to be prioritized in independent follow-ups, and pinpoint ion channels as important players and potential therapeutic targets warranting further investigation.

Genetically low vitamin D concentrations and myopic refractive error: a Mendelian randomization study

Background

Myopia prevalence has increased in the past 20 years, with many studies linking the increase to reduced time spent outdoors. A number of recent observational studies have shown an inverse association between vitamin D [25(OH)D] serum levels and myopia. However, in such studies it is difficult to separate the effects of time outdoors and vitamin D levels. In this work we use Mendelian randomization (MR) to assess if genetically determined 25(OH)D levels contribute to the degree of myopia.

Methods

We performed MR using results from a meta-analysis of refractive error (RE) genome-wide association study (GWAS) that included 37 382 and 8 376 adult participants of European and Asian ancestry, respectively, published by the Consortium for Refractive Error And Myopia (CREAM). We used single nucleotide polymorphisms (SNPs) in the DHCR7, CYP2R1, GC and CYP24A1 genes with known effects on 25(OH)D concentration as instrumental variables (IV). We estimated the effect of 25(OH)D on myopia level using a Wald-type ratio estimator based on the effect estimates from the CREAM GWAS.

Results

Using the combined effect attributed to the four SNPs, the estimate for the effect of 25(OH)D on refractive error was -0.02 [95% confidence interval (CI) -0.09, 0.04] dioptres (D) per 10 nmol/l increase in 25(OH)D concentration in Caucasians and 0.01 (95% CI -0.17, 0.19) D per 10 nmol/l increase in Asians.

Conclusions

The tight confidence intervals on our estimates suggest the true contribution of vitamin D levels to degree of myopia is very small and indistinguishable from zero. Previous findings from observational studies linking vitamin D levels to myopia were likely attributable to the effects of confounding by time spent outdoors.

Effect of varying skin surface electrode position on electroretinogram responses recorded using a handheld stimulating and recording system

Purpose

A handheld device (the RETeval system, LKC Technologies) aims to increase the ease of electroretinogram (ERG) recording by using specially designed skin electrodes, rather than corneal electrodes. We explored effects of electrode position on response parameters recorded using this device.

Methods

Healthy adult twins were recruited from the TwinsUK cohort and underwent recording of light-adapted flicker ERGs (corresponding to international standard stimuli). In Group 1, skin electrodes were placed in a “comfortable” position, which was up to 20 mm below the lid margin. For subsequent participants (Group 2), the electrode was positioned 2 mm from the lid margin as recommended by the manufacturer. Amplitudes and peak times (averaged from both eyes) were compared between groups after age-matching and inclusion of only one twin per pair. Light-adapted flicker and flash ERGs were recorded for an additional 10 healthy subjects in two consecutive recording sessions: in the test eye, electrode position was varied from 2 to 10–20 mm below the lid margin between sessions; in the fellow (control) eye, the electrode was 2 mm below the lid margin throughout. Amplitudes and peak times (test eye normalised to control eye) were compared for the two sessions.

Results

Including one twin per pair, and age-matching yielded 28 individuals per group. Flicker ERG amplitudes were significantly lower for Group 1 than Group 2 participants (p = 0.0024). However, mean peak times did not differ between groups (p = 0.54). For the subjects in whom electrode position was changed between recording sessions, flash and flicker amplitudes were significantly lower when positioned further from the lid margin (p < 0.005), but peak times were similar (p > 0.5).

Conclusions

Moving the skin electrodes further from the lid margin significantly reduces response amplitudes, highlighting the importance of consistent electrode positioning. However, this does not significantly affect peak times. Thus, it may be feasible to adopt a more comfortable position in participants who cannot tolerate the recommended position if analysis is restricted to peak time parameters.

Gut microbiota associations with common diseases and prescription medications in a population-based cohort

The human gut microbiome has been associated with many health factors but variability between studies limits exploration of effects between them. Gut microbiota profiles are available for >2700 members of the deeply phenotyped TwinsUK cohort, providing a uniform platform for such comparisons. Here, we present gut microbiota association analyses for 38 common diseases and 51 medications within the cohort. We describe several novel associations, highlight associations common across multiple diseases, and determine which diseases and medications have the greatest association with the gut microbiota. These results provide a reference for future studies of the gut microbiome and its role in human health.

The human gut microbiome has been associated with many health factors, but variability between studies limits exploration of these effects. Here, Jackson et al. analyse gut microbiota associations for 38 common diseases and 51 medications within >2700 members of the TwinsUK cohort.

Sex-specific glioma genome-wide association study identifies new risk locus at 3p21.31 in females, and finds sex-differences in risk at 8q24.21

Incidence of glioma is approximately 50% higher in males. Previous analyses have examined exposures related to sex hormones in women as potential protective factors for these tumors, with inconsistent results. Previous glioma genome-wide association studies (GWAS) have not stratified by sex. Potential sex-specific genetic effects were assessed in autosomal SNPs and sex chromosome variants for all glioma, GBM and non-GBM patients using data from four previous glioma GWAS. Datasets were analyzed using sex-stratified logistic regression models and combined using meta-analysis. There were 4,831 male cases, 5,216 male controls, 3,206 female cases and 5,470 female controls. A significant association was detected at rs11979158 (7p11.2) in males only. Association at rs55705857 (8q24.21) was stronger in females than in males. A large region on 3p21.31 was identified with significant association in females only. The identified differences in effect of risk variants do not fully explain the observed incidence difference in glioma by sex.

Association of Raynaud's phenomenon with a polymorphism in the NOS1 gene

Background

Raynaud’s phenomenon (RP) describes the phenomenon of recurrent vasospasm of digital arteries, associated with skin colour changes: pallor, cyanosis and erythema. Twin studies have indicated a genetic predisposition for RP; however, the precise aetiology of RP remains unknown. It is thought that genetic variation in temperature-responsive or vasospastic genes might underlie RP so performed a candidate gene study in a large, population based sample. We assessed the association between RP and single nucleotide polymorphisms (SNPs) in the TRPA1, TRPM8, CALCA, CALCB and NOS1 genes.

Methods

Analysis included a total of 4276 individuals from the TwinsUK database. RP status had been determined using validated, self-administered questionnaires and was diagnosed in 640 individuals (17.6%). 66 tag SNPs across the candidate genes were tested for association with RP status using a linear regression model, accounting for covariates. Adjustment was made for multiple testing. RegulomeDB and GTEx databases were used to assess possible functional effects of the polymorphisms.

Results

Nominally significant associations between RP and four SNPs in NOS1 and one in CALCB were identified. After permutation testing, rs527590 SNP in NOS1 passed the significance threshold. RegulomeDB scores indicated an unlikely functional effect of this variant, while the survey of the GTEx database found the SNP and several variants in linkage disequilibrium to be cis-eQTLs in skin.

Conclusion

Results indicate that RP is associated with variation in gene NOS1. This finding may be related to the observation that the significant SNP in NOS1 is known to exhibit functional influence on the gene expression.