Fat metabolism is associated with telomere length in six population-based studies.

Telomeres are repetitive DNA sequences located at the end of chromosomes, which are associated to biological aging, cardiovascular disease, cancer and mortality. Lipid and fatty acid metabolism have been associated with telomere shortening. We have conducted an in-depth study investigating the association of metabolic biomarkers with telomere length (LTL). We performed an association analysis of 226 metabolic biomarkers with LTL using data from 11 775 individuals from six independent population-based cohorts (BBMRI-NL consortium). Metabolic biomarkers include lipoprotein lipids and subclasses, fatty acids, amino acids, glycolysis measures and ketone bodies. LTL was measured by quantitative polymerase chain reaction or FlowFISH. Linear regression analysis was performed adjusting for age, sex, lipid-lowering medication and cohort-specific covariates (model 1) and additionally for body mass index (BMI) and smoking (model 2), followed by inverse variance-weighted meta-analyses (significance threshold Pmeta = 6.5 × 10-4). We identified four metabolic biomarkers positively associated with LTL, including two cholesterol to lipid ratios in small VLDL (S-VLDL-C % and S-VLDL-CE %) and two omega-6 fatty acid ratios (FAw6/FA and LA/FA). After additionally adjusting for BMI and smoking, these metabolic biomarkers remained associated with LTL with similar effect estimates. In addition, cholesterol esters in very small VLDL (XS-VLDL-CE) became significantly associated with LTL (P = 3.6 × 10-4). We replicated the association of FAw6/FA with LTL in an independent dataset of 7845 individuals (P = 1.9 × 10-4). To conclude, we identified multiple metabolic biomarkers involved in lipid and fatty acid metabolism that may be involved in LTL biology. Longitudinal studies are needed to exclude reversed causation.

Genome-wide association study meta-analysis of dizygotic twinning illuminates genetic regulation of female fecundity.

STUDY QUESTION: Which genetic factors regulate female propensity for giving birth to spontaneous dizygotic (DZ) twins? SUMMARY ANSWER: We identified four new loci, GNRH1, FSHR, ZFPM1, and IPO8, in addition to previously identified loci, FSHB and SMAD3. WHAT IS KNOWN ALREADY: The propensity to give birth to DZ twins runs in families. Earlier, we reported that FSHB and SMAD3 as associated with DZ twinning and female fertility measures. STUDY DESIGN, SIZE, DURATION: We conducted a genome-wide association meta-analysis (GWAMA) of mothers of spontaneous dizygotic (DZ) twins (8265 cases, 264 567 controls) and of independent DZ twin offspring (26 252 cases, 417 433 controls). PARTICIPANTS/MATERIALS, SETTING, METHODS: Over 700 000 mothers of DZ twins, twin individuals and singletons from large cohorts in Australia/New Zealand, Europe, and the USA were carefully screened to exclude twins born after use of ARTs. Genetic association analyses by cohort were followed by meta-analysis, phenome wide association studies (PheWAS), in silico and in vivo annotations, and Zebrafish functional validation. MAIN RESULTS AND THE ROLE OF CHANCE: This study enlarges the sample size considerably from previous efforts, finding four genome-wide significant loci, including two novel signals and a further two novel genes that are implicated by gene level enrichment analyses. The novel loci, GNRH1 and FSHR, have well-established roles in female reproduction whereas ZFPM1 and IPO8 have not previously been implicated in female fertility. We found significant genetic correlations with multiple aspects of female reproduction and body size as well as evidence for significant selection against DZ twinning during human evolution. The 26 top single nucleotide polymorphisms (SNPs) from our GWAMA in European-origin participants weakly predicted the crude twinning rates in 47 non-European populations (r = 0.23 between risk score and population prevalence, s.e. 0.11, 1-tail P = 0.058) indicating that genome-wide association studies (GWAS) are needed in African and Asian populations to explore the causes of their respectively high and low DZ twinning rates. In vivo functional tests in zebrafish for IPO8 validated its essential role in female, but not male, fertility. In most regions, risk SNPs linked to known expression quantitative trait loci (eQTLs). Top SNPs were associated with in vivo reproductive hormone levels with the top pathways including hormone ligand binding receptors and the ovulation cycle. LARGE SCALE DATA: The full DZT GWAS summary statistics will made available after publication through the GWAS catalog (https://www.ebi.ac.uk/gwas/). LIMITATIONS, REASONS FOR CAUTION: Our study only included European ancestry cohorts. Inclusion of data from Africa (with the highest twining rate) and Asia (with the lowest rate) would illuminate further the biology of twinning and female fertility. WIDER IMPLICATIONS OF THE FINDINGS: About one in 40 babies born in the world is a twin and there is much speculation on why twinning runs in families. We hope our results will inform investigations of ovarian response in new and existing ARTs and the causes of female infertility. STUDY FUNDING/COMPETING INTEREST(S): Support for the Netherlands Twin Register came from the Netherlands Organization for Scientific Research (NWO) and The Netherlands Organization for Health Research and Development (ZonMW) grants, 904-61-193, 480-04-004, 400-05-717, Addiction-31160008, 911-09-032, Biobanking and Biomolecular Resources Research Infrastructure (BBMRI.NL, 184.021.007), Royal Netherlands Academy of Science Professor Award (PAH/6635) to DIB, European Research Council (ERC-230374), Rutgers University Cell and DNA Repository (NIMH U24 MH068457-06), the Avera Institute, Sioux Falls, South Dakota (USA) and the National Institutes of Health (NIH R01 HD042157-01A1) and the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health and Grand Opportunity grants 1RC2 MH089951. The QIMR Berghofer Medical Research Institute (QIMR) study was supported by grants from the National Health and Medical Research Council (NHMRC) of Australia (241944, 339462, 389927, 389875, 389891, 389892, 389938, 443036, 442915, 442981, 496610, 496739, 552485, 552498, 1050208, 1075175). L.Y. is funded by Australian Research Council (Grant number DE200100425). The Minnesota Center for Twin and Family Research (MCTFR) was supported in part by USPHS Grants from the National Institute on Alcohol Abuse and Alcoholism (AA09367 and AA11886) and the National Institute on Drug Abuse (DA05147, DA13240, and DA024417). The Women's Genome Health Study (WGHS) was funded by the National Heart, Lung, and Blood Institute (HL043851 and HL080467) and the National Cancer Institute (CA047988 and UM1CA182913), with support for genotyping provided by Amgen. Data collection in the Finnish Twin Registry has been supported by the Wellcome Trust Sanger Institute, the Broad Institute, ENGAGE-European Network for Genetic and Genomic Epidemiology, FP7-HEALTH-F4-2007, grant agreement number 201413, National Institute of Alcohol Abuse and Alcoholism (grants AA-12502, AA-00145, AA-09203, AA15416, and K02AA018755) and the Academy of Finland (grants 100499, 205585, 118555, 141054, 264146, 308248, 312073 and 336823 to J. Kaprio). TwinsUK is funded by the Wellcome Trust, Medical Research Council, Versus Arthritis, European Union Horizon 2020, Chronic Disease Research Foundation (CDRF), Zoe Ltd and the National Institute for Health Research (NIHR) Clinical Research Network (CRN) and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. For NESDA, funding was obtained from the Netherlands Organization for Scientific Research (Geestkracht program grant 10000-1002), the Center for Medical Systems Biology (CSMB, NVVO Genomics), Biobanking and Biomolecular Resources Research Infrastructure (BBMRI-NL), VU University's Institutes for Health and Care Research (EMGO+) and Neuroscience Campus Amsterdam, University Medical Center Groningen, Leiden University Medical Center, National Institutes of Health (NIH, ROI D0042157-01A, MH081802, Grand Opportunity grants 1 RC2 Ml-1089951 and IRC2 MH089995). Part of the genotyping and analyses were funded by the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health. Computing was supported by BiG Grid, the Dutch e-Science Grid, which is financially supported by NWO. Work in the Del Bene lab was supported by the Programme Investissements d'Avenir IHU FOReSIGHT (ANR-18-IAHU-01). C.R. was supported by an EU Horizon 2020 Marie Sklodowska-Curie Action fellowship (H2020-MSCA-IF-2014 #661527). H.S. and K.S. are employees of deCODE Genetics/Amgen. The other authors declare no competing financial interests. TRIAL REGISTRATION NUMBER: N/A.

Maternal Characteristics in Natural and Medically Assisted Reproduction Dizygotic Twin Pregnancies.

Previous studies have shown that mothers of naturally conceived dizygotic (DZ) twins tend to be taller, older, and smoke more than mothers of naturally conceived monozygotic (MZ) twin and mothers of singletons. Here, we investigate whether mothers of naturally conceived DZ twins differ from mothers who conceived their DZ twins after medically assisted reproduction (MAR) in eight maternal traits related to fertility based on observational survey data. We include data from 33,648 mothers from the Netherlands Twin Register (NTR) and 1660 mothers of twins from the Norwegian Mother, Father and Child Cohort Study (MoBA). We contrast mothers of naturally conceived DZ twins with mothers of MAR DZ twins. Next, we further segment the MAR group into mothers who underwent hormonal induction of ovulation but not in vitro fertilization (IVF) and those who IVF twins, comparing them both to each other and against the mothers of naturally conceived DZ twins. Mothers of naturally conceived DZ twins smoke more often, differ in body composition, have a higher maternal age and have more offspring before the twins than mothers of MZ twins. Compared to MAR DZ twin mothers, mothers of naturally conceived DZ twins have fewer miscarriages, lower maternal age and increased height, more offspring and are more often smokers. BMI before the twin pregnancy is similar in both natural and MAR DZ twin mothers. Mothers who received hormonal induction of ovulation (OI) have a lower maternal age, fewer miscarriages, and a higher number of offspring before their twin pregnancy than twin mothers who received IVF and/or intracytoplasmic sperm injection (ICSI) treatments. Our study shows that twin mothers are a heterogenous group and the differences between twin mothers should be taken into account in epidemiological and genetic research that includes twins.

Genetic meta-analysis of twin birth weight shows high genetic correlation with singleton birth weight.

Birth weight (BW) is an important predictor of newborn survival and health and has associations with many adult health outcomes, including cardiometabolic disorders, autoimmune diseases and mental health. On average, twins have a lower BW than singletons as a result of a different pattern of fetal growth and shorter gestational duration. Therefore, investigations into the genetics of BW often exclude data from twins, leading to a reduction in sample size and remaining ambiguities concerning the genetic contribution to BW in twins. In this study, we carried out a genome-wide association meta-analysis of BW in 42 212 twin individuals and found a positive correlation of beta values (Pearson's r = 0.66, 95% confidence interval [CI]: 0.47-0.77) with 150 previously reported genome-wide significant variants for singleton BW. We identified strong positive genetic correlations between BW in twins and numerous anthropometric traits, most notably with BW in singletons (genetic correlation [rg] = 0.92, 95% CI: 0.66-1.18). Genetic correlations of BW in twins with a series of health-related traits closely resembled those previously observed for BW in singletons. Polygenic scores constructed from a genome-wide association study on BW in the UK Biobank demonstrated strong predictive power in a target sample of Dutch twins and singletons. Together, our results indicate that a similar genetic architecture underlies BW in twins and singletons and that future genome-wide studies might benefit from including data from large twin registers.

Genome-wide Association Analysis in Humans Links Nucleotide Metabolism to Leukocyte Telomere Length.

Leukocyte telomere length (LTL) is a heritable biomarker of genomic aging. In this study, we perform a genome-wide meta-analysis of LTL by pooling densely genotyped and imputed association results across large-scale European-descent studies including up to 78,592 individuals. We identify 49 genomic regions at a false dicovery rate (FDR) < 0.05 threshold and prioritize genes at 31, with five highlighting nucleotide metabolism as an important regulator of LTL. We report six genome-wide significant loci in or near SENP7, MOB1B, CARMIL1, PRRC2A, TERF2, and RFWD3, and our results support recently identified PARP1, POT1, ATM, and MPHOSPH6 loci. Phenome-wide analyses in >350,000 UK Biobank participants suggest that genetically shorter telomere length increases the risk of hypothyroidism and decreases the risk of thyroid cancer, lymphoma, and a range of proliferative conditions. Our results replicate previously reported associations with increased risk of coronary artery disease and lower risk for multiple cancer types. Our findings substantially expand current knowledge on genes that regulate LTL and their impact on human health and disease.

Longitudinal multi-omics study reveals common etiology underlying association between plasma proteome and BMI trajectories in adolescent and young adult twins.

BACKGROUND: The influence of genetics and environment on the association of the plasma proteome with body mass index (BMI) and changes in BMI remains underexplored, and the links to other omics in these associations remain to be investigated. We characterized protein-BMI trajectory associations in adolescents and adults and how these connect to other omics layers. METHODS: Our study included two cohorts of longitudinally followed twins: FinnTwin12 (N = 651) and the Netherlands Twin Register (NTR) (N = 665). Follow-up comprised 4 BMI measurements over approximately 6 (NTR: 23-27 years old) to 10 years (FinnTwin12: 12-22 years old), with omics data collected at the last BMI measurement. BMI changes were calculated in latent growth curve models. Mixed-effects models were used to quantify the associations between the abundance of 439 plasma proteins with BMI at blood sampling and changes in BMI. In FinnTwin12, the sources of genetic and environmental variation underlying the protein abundances were quantified by twin models, as were the associations of proteins with BMI and BMI changes. In NTR, we investigated the association of gene expression of genes encoding proteins identified in FinnTwin12 with BMI and changes in BMI. We linked identified proteins and their coding genes to plasma metabolites and polygenic risk scores (PRS) applying mixed-effects models and correlation networks. RESULTS: We identified 66 and 14 proteins associated with BMI at blood sampling and changes in BMI, respectively. The average heritability of these proteins was 35%. Of the 66 BMI-protein associations, 43 and 12 showed genetic and environmental correlations, respectively, including 8 proteins showing both. Similarly, we observed 7 and 3 genetic and environmental correlations between changes in BMI and protein abundance, respectively. S100A8 gene expression was associated with BMI at blood sampling, and the PRG4 and CFI genes were associated with BMI changes. Proteins showed strong connections with metabolites and PRSs, but we observed no multi-omics connections among gene expression and other omics layers. CONCLUSIONS: Associations between the proteome and BMI trajectories are characterized by shared genetic, environmental, and metabolic etiologies. We observed few gene-protein pairs associated with BMI or changes in BMI at the proteome and transcriptome levels.

Methylome-wide association study of antidepressant use in Generation Scotland and the Netherlands Twin Register implicates the innate immune system.

Antidepressants are an effective treatment for major depressive disorder (MDD), although individual response is unpredictable and highly variable. Whilst the mode of action of antidepressants is incompletely understood, many medications are associated with changes in DNA methylation in genes that are plausibly linked to their mechanisms. Studies of DNA methylation may therefore reveal the biological processes underpinning the efficacy and side effects of antidepressants. We performed a methylome-wide association study (MWAS) of self-reported antidepressant use accounting for lifestyle factors and MDD in Generation Scotland (GS:SFHS, N = 6428, EPIC array) and the Netherlands Twin Register (NTR, N = 2449, 450 K array) and ran a meta-analysis of antidepressant use across these two cohorts. We found ten CpG sites significantly associated with self-reported antidepressant use in GS:SFHS, with the top CpG located within a gene previously associated with mental health disorders, ATP6V1B2 (ß = -0.055, pcorrected = 0.005). Other top loci were annotated to genes including CASP10, TMBIM1, MAPKAPK3, and HEBP2, which have previously been implicated in the innate immune response. Next, using penalised regression, we trained a methylation-based score of self-reported antidepressant use in a subset of 3799 GS:SFHS individuals that predicted antidepressant use in a second subset of GS:SFHS (N = 3360, ß = 0.377, p = 3.12 × 10-11, R2 = 2.12%). In an MWAS analysis of prescribed selective serotonin reuptake inhibitors, we showed convergent findings with those based on self-report. In NTR, we did not find any CpGs significantly associated with antidepressant use. The meta-analysis identified the two CpGs of the ten above that were common to the two arrays used as being significantly associated with antidepressant use, although the effect was in the opposite direction for one of them. Antidepressants were associated with epigenetic alterations in loci previously associated with mental health disorders and the innate immune system. These changes predicted self-reported antidepressant use in a subset of GS:SFHS and identified processes that may be relevant to our mechanistic understanding of clinically relevant antidepressant drug actions and side effects.

Handedness and 23 Early Life Characteristics in 37,495 Dutch Twins.

In studies of singletons, a range of early-life characteristics have been reported to be associated with handedness, but some of these associations have failed to replicate. We examined associations between 23 early life characteristics with handedness in a large sample of 37,495 5-year-old twins. We considered three definitions of handedness: left-handedness (LH), mixed-handedness (MH), and non-right-handedness (NRH). Our main aim was to test whether the associations with sex, birth weight, gestational age, and season of birth - as reported in singletons - replicate in twins, and to examine twin-specific variables, including zygosity, chorionicity, birth order, and intertwin delivery time. Compared to previously published data from adults born as singletons (7.23%), the prevalence of NRH was higher in both twins (16.19%) and their parents (15.09%). In the twins, LH and NRH were associated with parents' LH. Male sex and lower gestational age were associated with NRH, and LH was associated with not being breastfed. MH was related to neurodevelopmental delays and higher externalizing problems later in childhood. Other previously reported associations were not replicated, and no twin-specific characteristics were related to handedness. These results emphasize the importance of considering multiple definitions of handedness and indicate a small number of replicated associations across studies.

More Than Results: The Clinical and Research Relationship in the Evolving Detection and Surveillance of SARS-CoV-2.

INTRODUCTION: During the coronavirus disease 2019 (COVID-19) pandemic, real-time reverse transcription polymerase chain reaction (RT-PCR) became an essential tool for laboratories to provide high-sensitivity qualitative diagnostic testing for patients and real-time data to public health officials. Here we explore the predictive value of quantitative data from RT-PCR cycle threshold (Ct) values in epidemiological measures, symptom presentation, and variant transition. METHODS: To examine the association with hospitalizations and deaths, data from 74,479 patients referred to the Avera Institute for Human Genetics (AIHG) for COVID-19 testing in 2020 were matched by calendar week to epidemiological data reported by the South Dakota Department of Health. We explored the association between symptom data, patient age, and Ct values for 101 patients. We also explored changes in Ct values during variant transition detected by genomic surveillance sequencing of the AIHG testing population during 2021. RESULTS: Measures from AIHG diagnostic testing strongly explain variance in the South Dakota state positivity percentage (R2 = 0.758), a two-week delay in hospitalizations (R2 = 0.856), and a four-week delay in deaths (R2 = 0.854). Based on factor analysis of patient symptoms, three groups could be distinguished which had different presentations of age, Ct value, and time from collection. Additionally, conflicting Ct value results among SARSCoV- 2 variants during variant transition may reflect the community transmission dynamics. CONCLUSIONS: Measures of Ct value in RT-PCR diagnostic assays combined with routine screening have valuable applications in monitoring the dynamics of SARS-CoV-2 within communities.

Healthy Cotwins Share Gut Microbiome Signatures With Their Inflammatory Bowel Disease Twins and Unrelated Patients.

BACKGROUND & AIMS: It is currently unclear whether reported changes in the gut microbiome are cause or consequence of inflammatory bowel disease (IBD). Therefore, we studied the gut microbiome of IBD-discordant and -concordant twin pairs, which offers the unique opportunity to assess individuals at increased risk of developing IBD, namely healthy cotwins from IBD-discordant twin pairs. METHODS: Fecal samples were obtained from 99 twins (belonging to 51 twin pairs), 495 healthy age-, sex-, and body mass index-matched controls, and 99 unrelated patients with IBD. Whole-genome metagenomic shotgun sequencing was performed. Taxonomic and functional (pathways) composition was compared among healthy cotwins, IBD-twins, unrelated patients with IBD, and healthy controls with multivariable (ie, adjusted for potential confounding) generalized linear models. RESULTS: No significant differences were observed in the relative abundance of species and pathways between healthy cotwins and their IBD-twins (false discovery rate <0.10). Compared with healthy controls, 13, 19, and 18 species, and 78, 105, and 153 pathways were found to be differentially abundant in healthy cotwins, IBD-twins, and unrelated patients with IBD, respectively (false discovery rate <0.10). Of these, 8 (42.1%) of 19 and 1 (5.6%) of 18 species, and 37 (35.2%) of 105 and 30 (19.6%) of 153 pathways overlapped between healthy cotwins and IBD-twins, and healthy cotwins and unrelated patients with IBD, respectively. Many of the shared species and pathways have previously been associated with IBD. The shared pathways include potentially inflammation-related pathways, for example, an increase in propionate degradation and L-arginine degradation pathways. CONCLUSIONS: The gut microbiome of healthy cotwins from IBD-discordant twin pairs displays IBD-like signatures. These IBD-like microbiome signatures might precede the onset of IBD. However, longitudinal follow-up studies are needed to infer a causal relationship.

Changing genetic architecture of body mass index from infancy to early adulthood: an individual based pooled analysis of 25 twin cohorts.

BACKGROUND: Body mass index (BMI) shows strong continuity over childhood and adolescence and high childhood BMI is the strongest predictor of adult obesity. Genetic factors strongly contribute to this continuity, but it is still poorly known how their contribution changes over childhood and adolescence. Thus, we used the genetic twin design to estimate the genetic correlations of BMI from infancy to adulthood and compared them to the genetic correlations of height. METHODS: We pooled individual level data from 25 longitudinal twin cohorts including 38,530 complete twin pairs and having 283,766 longitudinal height and weight measures. The data were analyzed using Cholesky decomposition offering genetic and environmental correlations of BMI and height between all age combinations from 1 to 19 years of age. RESULTS: The genetic correlations of BMI and height were stronger than the trait correlations. For BMI, we found that genetic correlations decreased as the age between the assessments increased, a trend that was especially visible from early to middle childhood. In contrast, for height, the genetic correlations were strong between all ages. Age-to-age correlations between environmental factors shared by co-twins were found for BMI in early childhood but disappeared altogether by middle childhood. For height, shared environmental correlations persisted from infancy to adulthood. CONCLUSIONS: Our results suggest that the genes affecting BMI change over childhood and adolescence leading to decreasing age-to-age genetic correlations. This change is especially visible from early to middle childhood indicating that new genetic factors start to affect BMI in middle childhood. Identifying mediating pathways of these genetic factors can open possibilities for interventions, especially for those children with high genetic predisposition to adult obesity.

Gene-by-Crisis Interaction for Optimism and Meaning in Life: The Effects of the COVID-19 Pandemic.

The corona virus disease 2019 (COVID-19) pandemic and the restrictions to reduce the spread of the virus has had a large impact on daily life. We investigated the individual differences in the effect of the COVID-19 pandemic and first lockdown on optimism and meaning in life in a sample from the Netherlands Twin Register. Participants completed surveys before (N = 9964, Mean age: 48.2, SD = 14.4) and during the first months of the pandemic (i.e. April-May 2020, N = 17,464, Mean age: 44.6 SD = 14.8), with a subsample completing both surveys (N = 6461, Mean age T1: 48.8, SD = 14.5). We applied genetic covariance structure models to twin data to investigate changes in the genetic architecture of the outcome traits due to the pandemic and the interaction of genes with the environmental exposure. Although 56% and 35% of the sample was negatively affected by the pandemic in their optimism and meaning in life, many participants were stable (32% and 43%) or even showed increased optimism and meaning in life (11% and 22%). Subgroups, specifically women, higher educated people, and people with poorer health, experienced larger negative effects. During the first months of the pandemic, slightly lower heritability estimates for optimism and meaning in life (respectively 20% and 25%) were obtained compared to pre-pandemic (respectively 26% and 32%), although confidence intervals overlap. The lower than unity genetic correlations across time (.75 and .63) suggest gene-environment interactions, where the expression of genes that influence optimism and meaning in life differs before and during the pandemic. The COVID-19 pandemic is a strong exposure that leads to imbalanced effects on the well-being of individuals. Some people decrease in well-being, while others get more optimistic and consider their lives as more meaningful during the pandemic. These differences are partly explained by individual differences in genetic sensitivity to extreme environmental change. More knowledge on the person-specific response to specific environmental variables underlying these individual differences is urgently needed to prevent further inequality.

Genetic factors explain a significant part of associations between adolescent well-being and the social environment.

Socio-environmental factors play an important role in adolescent well-being, but potential genetic contributions to these associations are rarely assessed. To address this gap in the literature, associations between well-being and family conflict and functioning, number of friends, friendship importance and satisfaction, and leisure time variables were studied in N = ~ 4700 twin pairs from the Netherlands Twin Register, us ing generalized estimating equations and twin-difference scores. When twin-difference scores indicated a role for genetic factors, we used bivariate genetic models to quantify genetic and environmental contributions to these associations. We identify significant associations between well-being and family functioning, family conflict, different leisure time activities, number of friends, and satisfaction with friendships. Additionally, we find evidence for large (73-91%) genetic influence on the associations between well-being and family conflict and functioning, leisure time sport/scouting clubs, and satisfaction with friendships. Finally, findings support the hypothesis of a causal association between well-being and family conflict and functioning. These findings have important implications for research into the social correlates of well-being in adolescence, as not taking genetic factors into account leads to overestimations of the influence of identified correlates and consequently to recommendations of these correlates as intervention targets.

Phenome-wide investigation of health outcomes associated with genetic predisposition to loneliness.

Humans are social animals that experience intense suffering when they perceive a lack of social connection. Modern societies are experiencing an epidemic of loneliness. Although the experience of loneliness is universally human, some people report experiencing greater loneliness than others. Loneliness is more strongly associated with mortality than obesity, emphasizing the need to understand the nature of the relationship between loneliness and health. Although it is intuitive that circumstantial factors such as marital status and age influence loneliness, there is also compelling evidence of a genetic predisposition toward loneliness. To better understand the genetic architecture of loneliness and its relationship with associated outcomes, we extended the genome-wide association study meta-analysis of loneliness to 511 280 subjects, and detect 19 significant genetic variants from 16 loci, including four novel loci, as well as 58 significantly associated genes. We investigated the genetic overlap with a wide range of physical and mental health traits by computing genetic correlations and by building loneliness polygenic scores in an independent sample of 18 498 individuals with EHR data to conduct a PheWAS with. A genetic predisposition toward loneliness was associated with cardiovascular, psychiatric, and metabolic disorders and triglycerides and high-density lipoproteins. Mendelian randomization analyses showed evidence of a causal, increasing, the effect of both BMI and body fat on loneliness. Our results provide a framework for future studies of the genetic basis of loneliness and its relationship to mental and physical health.

Selective maternal seeding and environment shape the human gut microbiome.

Vertical transmission of bacteria from mother to infant at birth is postulated to initiate a life-long host-microbe symbiosis, playing an important role in early infant development. However, only the tracking of strictly defined unique microbial strains can clarify where the intestinal bacteria come from, how long the initial colonizers persist, and whether colonization by other strains from the environment can replace existing ones. Using rare single nucleotide variants in fecal metagenomes of infants and their family members, we show strong evidence of selective and persistent transmission of maternal strain populations to the vaginally born infant and their occasional replacement by strains from the environment, including those from family members, in later childhood. Only strains from the classes Actinobacteria and Bacteroidia, which are essential components of the infant microbiome, are transmitted from the mother and persist for at least 1 yr. In contrast, maternal strains of Clostridia, a dominant class in the mother's gut microbiome, are not observed in the infant. Caesarean-born infants show a striking lack of maternal transmission at birth. After the first year, strain influx from the family environment occurs and continues even in adulthood. Fathers appear to be more frequently donors of novel strains to other family members than receivers. Thus, the infant gut is seeded by selected maternal bacteria, which expand to form a stable community, with a rare but stable continuing strain influx over time.

Investigating the relationships between unfavourable habitual sleep and metabolomic traits: evidence from multi-cohort multivariable regression and Mendelian randomization analyses.

BACKGROUND: Sleep traits are associated with cardiometabolic disease risk, with evidence from Mendelian randomization (MR) suggesting that insomnia symptoms and shorter sleep duration increase coronary artery disease risk. We combined adjusted multivariable regression (AMV) and MR analyses of phenotypes of unfavourable sleep on 113 metabolomic traits to investigate possible biochemical mechanisms linking sleep to cardiovascular disease. METHODS: We used AMV (N = 17,368) combined with two-sample MR (N = 38,618) to examine effects of self-reported insomnia symptoms, total habitual sleep duration, and chronotype on 113 metabolomic traits. The AMV analyses were conducted on data from 10 cohorts of mostly Europeans, adjusted for age, sex, and body mass index. For the MR analyses, we used summary results from published European-ancestry genome-wide association studies of self-reported sleep traits and of nuclear magnetic resonance (NMR) serum metabolites. We used the inverse-variance weighted (IVW) method and complemented this with sensitivity analyses to assess MR assumptions. RESULTS: We found consistent evidence from AMV and MR analyses for associations of usual vs. sometimes/rare/never insomnia symptoms with lower citrate (- 0.08 standard deviation (SD)[95% confidence interval (CI) - 0.12, - 0.03] in AMV and - 0.03SD [- 0.07, - 0.003] in MR), higher glycoprotein acetyls (0.08SD [95% CI 0.03, 0.12] in AMV and 0.06SD [0.03, 0.10) in MR]), lower total very large HDL particles (- 0.04SD [- 0.08, 0.00] in AMV and - 0.05SD [- 0.09, - 0.02] in MR), and lower phospholipids in very large HDL particles (- 0.04SD [- 0.08, 0.002] in AMV and - 0.05SD [- 0.08, - 0.02] in MR). Longer total sleep duration associated with higher creatinine concentrations using both methods (0.02SD per 1 h [0.01, 0.03] in AMV and 0.15SD [0.02, 0.29] in MR) and with isoleucine in MR analyses (0.22SD [0.08, 0.35]). No consistent evidence was observed for effects of chronotype on metabolomic measures. CONCLUSIONS: Whilst our results suggested that unfavourable sleep traits may not cause widespread metabolic disruption, some notable effects were observed. The evidence for possible effects of insomnia symptoms on glycoprotein acetyls and citrate and longer total sleep duration on creatinine and isoleucine might explain some of the effects, found in MR analyses of these sleep traits on coronary heart disease, which warrant further investigation.

Educational attainment of same-sex and opposite-sex dizygotic twins: An individual-level pooled study of 19 twin cohorts.

Comparing twins from same- and opposite-sex pairs can provide information on potential sex differences in a variety of outcomes, including socioeconomic-related outcomes such as educational attainment. It has been suggested that this design can be applied to examine the putative role of intrauterine exposure to testosterone for educational attainment, but the evidence is still disputed. Thus, we established an international database of twin data from 11 countries with 88,290 individual dizygotic twins born over 100 years and tested for differences between twins from same- and opposite-sex dizygotic pairs in educational attainment. Effect sizes with 95% confidence intervals (CI) were estimated by linear regression models after adjusting for birth year and twin study cohort. In contrast to the hypothesis, no difference was found in women (ß = -0.05 educational years, 95% CI -0.11, 0.02). However, men with a same-sex co-twin were slightly more educated than men having an opposite-sex co-twin (ß = 0.14 educational years, 95% CI 0.07, 0.21). No consistent differences in effect sizes were found between individual twin study cohorts representing Europe, the USA, and Australia or over the cohorts born during the 20th century, during which period the sex differences in education reversed favoring women in the latest birth cohorts. Further, no interaction was found with maternal or paternal education. Our results contradict the hypothesis that there would be differences in the intrauterine testosterone levels between same-sex and opposite-sex female twins affecting education. Our findings in men may point to social dynamics within same-sex twin pairs that may benefit men in their educational careers.

Safe Linkage of Cohort and Population-Based Register Data in a Genomewide Association Study on Health Care Expenditure.

There are research questions whose answers require record linkage of multiple databases that may be characterized by limited options for full data sharing. For this purpose, the Open Data Infrastructure for Social Science and Economic Innovations (ODISSEI) consortium has supported the development of the ODISSEI Secure Supercomputer (OSSC) platform that allows researchers to link cohort data to data from Statistics Netherlands and run large-scale analyses in a high-performance computing (HPC) environment. Here, we report a successful record linkage genomewide association (GWA) study on expenditure for total health, mental health, primary and hospital care, and medication. Record linkage for genotype data from 16,726 participants from the Netherlands Twin Register (NTR) with data from Statistics Netherlands was accomplished in the secure OSSC platform, followed by gene-based tests and estimation of total and single nucleotide polymorphism (SNP)-based heritability. The total heritability of expenditure ranged between 29.4% (SE 0.8) and 37.5% (SE 0.8), but GWA analyses did not identify SNPs or genes that were genomewide significantly associated with health care expenditure. SNP-based heritability was between 0.0% (SE 3.5) and 5.4% (SE 4.0) and was different from zero for mental health care and primary care expenditure. We conclude that successfully linking genotype data to administrative health care expenditure data from Statistics Netherlands is feasible and demonstrates a series of analyses on health care expenditure. The OSSC platform offers secure possibilities for analyzing linked data in large scale and realizing sample sizes required for GWA studies, providing invaluable opportunities to answer many new research questions.

DNA methylation age is associated with an altered hemostatic profile in a multiethnic meta-analysis.

Many hemostatic factors are associated with age and age-related diseases; however, much remains unknown about the biological mechanisms linking aging and hemostatic factors. DNA methylation is a novel means by which to assess epigenetic aging, which is a measure of age and the aging processes as determined by altered epigenetic states. We used a meta-analysis approach to examine the association between measures of epigenetic aging and hemostatic factors, as well as a clotting time measure. For fibrinogen, we performed European and African ancestry-specific meta-analyses which were then combined via a random effects meta-analysis. For all other measures we could not estimate ancestry-specific effects and used a single fixed effects meta-analysis. We found that 1-year higher extrinsic epigenetic age as compared with chronological age was associated with higher fibrinogen (0.004 g/L/y; 95% confidence interval, 0.001-0.007; P = .01) and plasminogen activator inhibitor 1 (PAI-1; 0.13 U/mL/y; 95% confidence interval, 0.07-0.20; P = 6.6 × 10-5) concentrations, as well as lower activated partial thromboplastin time, a measure of clotting time. We replicated PAI-1 associations using an independent cohort. To further elucidate potential functional mechanisms, we associated epigenetic aging with expression levels of the PAI-1 protein encoding gene (SERPINE1) and the 3 fibrinogen subunit-encoding genes (FGA, FGG, and FGB) in both peripheral blood and aorta intima-media samples. We observed associations between accelerated epigenetic aging and transcription of FGG in both tissues. Collectively, our results indicate that accelerated epigenetic aging is associated with a procoagulation hemostatic profile, and that epigenetic aging may regulate hemostasis in part via gene transcription.

A Comparison of the ASEBA Adult Self Report (ASR) and the Brief Problem Monitor (BPM/18-59).

The adult self report (ASR) is a well-validated instrument with multiple scales relating to adult psychopathology. Recently, an 18-item version has been introduced, the brief problem monitor (BPM) to measure Internalizing behavior (INT), Externalizing behavior (EXT), and attention problems (ATT). The present study compared the BPM and ASR and investigated how well the BPM can serve as a supplement or an alternative for the ASR for specific clinical and scientific purposes. In a large sample of adult twins (N = 9.835) from the Netherlands Twin Register (NTR), we compared the internal consistency, clinical classification concordance, means, and variances of the ASR and BPM. Using the classical twin design, we investigated the genetic covariance structure. For external validation, the associations between subjective well-being and different subscales of the ASR and BPM were compared. The internal consistency of the BPM scales (around α = 0.75) was somewhat lower than the ASR (α ~ 0.85). The BPM Externalizing scale showed the lowest internal consistency (α = 0.63). ASR and BPM scores showed good clinical classification concordance (0.61-0.80) and high correlations (r > 0.88). A small reversed sex difference in the BPM Externalizing scale appeared (women > men). Genetic (0.34-0.54) and environmental components (0.46-0.66) explained the variance to a similar extent for the ASR and BPM. The phenotypic and genetic associations with well-being were comparable. In situations where sum scores are sufficient, the BPM performs as well as the longer ASR. Depending on the situation and goal, it is worth considering the BPM as an alternative for the ASR to reduce the participant burden.