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.

Combined genome-wide association study of 136 quantitative ear morphology traits in multiple populations reveal 8 novel loci.

Human ear morphology, a complex anatomical structure represented by a multidimensional set of correlated and heritable phenotypes, has a poorly understood genetic architecture. In this study, we quantitatively assessed 136 ear morphology traits using deep learning analysis of digital face images in 14,921 individuals from five different cohorts in Europe, Asia, and Latin America. Through GWAS meta-analysis and C-GWASs, a recently introduced method to effectively combine GWASs of many traits, we identified 16 genetic loci involved in various ear phenotypes, eight of which have not been previously associated with human ear features. Our findings suggest that ear morphology shares genetic determinants with other surface ectoderm-derived traits such as facial variation, mono eyebrow, and male pattern baldness. Our results enhance the genetic understanding of human ear morphology and shed light on the shared genetic contributors of different surface ectoderm-derived phenotypes. Additionally, gene editing experiments in mice have demonstrated that knocking out the newly ear-associated gene (Intu) and a previously ear-associated gene (Tbx15) causes deviating mouse ear morphology.

Bone mineral density loci specific to the skull portray potential pleiotropic effects on craniosynostosis.

Skull bone mineral density (SK-BMD) provides a suitable trait for the discovery of key genes in bone biology, particularly to intramembranous ossification, not captured at other skeletal sites. We perform a genome-wide association meta-analysis (n ~ 43,800) of SK-BMD, identifying 59 loci, collectively explaining 12.5% of the trait variance. Association signals cluster within gene-sets involved in skeletal development and osteoporosis. Among the four novel loci (ZIC1, PRKAR1A, AZIN1/ATP6V1C1, GLRX3), there are factors implicated in intramembranous ossification and as we show, inherent to craniosynostosis processes. Functional follow-up in zebrafish confirms the importance of ZIC1 on cranial suture patterning. Likewise, we observe abnormal cranial bone initiation that culminates in ectopic sutures and reduced BMD in mosaic atp6v1c1 knockouts. Mosaic prkar1a knockouts present asymmetric bone growth and, conversely, elevated BMD. In light of this evidence linking SK-BMD loci to craniofacial abnormalities, our study provides new insight into the pathophysiology, diagnosis and treatment of skeletal diseases.

Genetic prediction of male pattern baldness based on large independent datasets.

Genetic prediction of male pattern baldness (MPB) is important in science and society. Previous genetic MPB prediction models were limited by sparse marker coverage, small sample size, and/or data dependency in the different analytical steps. Here, we present novel models for genetic prediction of MPB based on a large set of markers and large independent subsample sets drawn among 187,435 European subjects. We selected 117 SNP predictors within 85 distinct loci from a list of 270 previously MPB-associated SNPs in 55,573 males of the UK Biobank Study (UKBB). Based on these 117 SNPs with and without age as additional predictor, we trained, by use of different methods, prediction models in a non-overlapping subset of 104,694 UKBB males and tested them in a non-overlapping subset of 26,177 UKBB males. Estimates of prediction accuracy were similar between methods with AUC ranges of 0.725-0.728 for severe, 0.631-0.635 for moderate, 0.598-0.602 for slight, and 0.708-0.711 for no hair loss with age, and slightly lower without, while prediction of any versus no hair loss gave 0.690-0.711 with age and slightly lower without. External validation in an early-onset enriched MPB dataset from the Bonn Study (N = 991) showed improved prediction accuracy without considering age such as AUC of 0.830 for no vs. any hair loss. Because of the large number of markers and the large independent datasets used for the different analytical steps, the newly presented genetic prediction models are the most reliable ones currently available for MPB or any other human appearance trait.

Codiversification of gut microbiota with humans.

The gut microbiomes of human populations worldwide have many core microbial species in common. However, within a species, some strains can show remarkable population specificity. The question is whether such specificity arises from a shared evolutionary history (codiversification) between humans and their microbes. To test for codiversification of host and microbiota, we analyzed paired gut metagenomes and human genomes for 1225 individuals in Europe, Asia, and Africa, including mothers and their children. Between and within countries, a parallel evolutionary history was evident for humans and their gut microbes. Moreover, species displaying the strongest codiversification independently evolved traits characteristic of host dependency, including reduced genomes and oxygen and temperature sensitivity. These findings all point to the importance of understanding the potential role of population-specific microbial strains in microbiome-mediated disease phenotypes.

Illness Characteristics of COVID-19 in Children Infected with the SARS-CoV-2 Delta Variant.

BACKGROUND: The Delta (B.1.617.2) SARS-CoV-2 variant was the predominant UK circulating strain between May and November 2021. We investigated whether COVID-19 from Delta infection differed from infection with previous variants in children. METHODS: Through the prospective COVID Symptom Study, 109,626 UK school-aged children were proxy-reported between 28 December 2020 and 8 July 2021. We selected all symptomatic children who tested positive for SARS-CoV-2 and were proxy-reported at least weekly, within two timeframes: 28 December 2020 to 6 May 2021 (Alpha (B.1.1.7), the main UK circulating variant) and 26 May to 8 July 2021 (Delta, the main UK circulating variant), with all children unvaccinated (as per national policy at the time). We assessed illness profiles (symptom prevalence, duration, and burden), hospital presentation, and presence of long (≥28 day) illness, and calculated odds ratios for symptoms presenting within the first 28 days of illness. RESULTS: 694 (276 younger (5-11 years), 418 older (12-17 years)) symptomatic children tested positive for SARS-CoV-2 with Alpha infection and 706 (227 younger and 479 older) children with Delta infection. Median illness duration was short with either variant (overall cohort: 5 days (IQR 2-9.75) with Alpha, 5 days (IQR 2-9) with Delta). The seven most prevalent symptoms were common to both variants. Symptom burden over the first 28 days was slightly greater with Delta compared with Alpha infection (in younger children, 3 (IQR 2-5) symptoms with Alpha, 4 (IQR 2-7) with Delta; in older children, 5 (IQR 3-8) symptoms with Alpha, 6 (IQR 3-9) with Delta infection ). The odds of presenting several symptoms were higher with Delta than Alpha infection, including headache and fever. Few children presented to hospital, and long illness duration was uncommon, with either variant. CONCLUSIONS: COVID-19 in UK school-aged children due to SARS-CoV-2 Delta strain B.1.617.2 resembles illness due to the Alpha variant B.1.1.7., with short duration and similar symptom burden.

Impacts of dietary exposure to pesticides on faecal microbiome metabolism in adult twins.

BACKGROUND: Dietary habits have a profound influence on the metabolic activity of gut microorganisms and their influence on health. Concerns have been raised as to whether the consumption of foodstuffs contaminated with pesticides can contribute to the development of chronic disease by affecting the gut microbiome. We performed the first pesticide biomonitoring survey of the British population, and subsequently used the results to perform the first pesticide association study on gut microbiome composition and function from the TwinsUK registry. METHODS: Dietary exposure of 186 common insecticide, herbicide, or fungicide residues and the faecal microbiome in 65 twin pairs in the UK was investigated. We evaluated if dietary habits, geographic location, or the rural/urban environment, are associated with the excretion of pesticide residues. The composition and metabolic activity of faecal microbiota was evaluated using shotgun metagenomics and metabolomics respectively. We performed a targeted urine metabolomics analysis in order to evaluate whether pesticide urinary excretion was also associated with physiological changes. RESULTS: Pyrethroid and/or organophosphorus insecticide residues were found in all urine samples, while the herbicide glyphosate was found in 53% of individuals. Food frequency questionnaires showed that residues from organophosphates were higher with increased consumption of fruit and vegetables. A total of 34 associations between pesticide residue concentrations and faecal metabolite concentrations were detected. Glyphosate excretion was positively associated with an overall increased bacterial species richness, as well as to fatty acid metabolites and phosphate levels. The insecticide metabolite Br2CA, reflecting deltamethrin exposure, was positively associated with the phytoestrogens enterodiol and enterolactone, and negatively associated with some N-methyl amino acids. Urine metabolomics performed on a subset of samples did not reveal associations with the excretion of pesticide residues. CONCLUSIONS: The consumption of conventionally grown fruit and vegetables leads to higher ingestion of pesticides with unknown long-term health consequences. Our results highlight the need for future dietary intervention studies to understand effects of pesticide exposure on the gut microbiome and possible health consequences.

Distinct clinical symptom patterns in patients hospitalised with COVID-19 in an analysis of 59,011 patients in the ISARIC-4C study.

COVID-19 is clinically characterised by fever, cough, and dyspnoea. Symptoms affecting other organ systems have been reported. However, it is the clinical associations of different patterns of symptoms which influence diagnostic and therapeutic decision-making. In this study, we applied clustering techniques to a large prospective cohort of hospitalised patients with COVID-19 to identify clinically meaningful sub-phenotypes. We obtained structured clinical data on 59,011 patients in the UK (the ISARIC Coronavirus Clinical Characterisation Consortium, 4C) and used a principled, unsupervised clustering approach to partition the first 25,477 cases according to symptoms reported at recruitment. We validated our findings in a second group of 33,534 cases recruited to ISARIC-4C, and in 4,445 cases recruited to a separate study of community cases. Unsupervised clustering identified distinct sub-phenotypes. First, a core symptom set of fever, cough, and dyspnoea, which co-occurred with additional symptoms in three further patterns: fatigue and confusion, diarrhoea and vomiting, or productive cough. Presentations with a single reported symptom of dyspnoea or confusion were also identified, alongside a sub-phenotype of patients reporting few or no symptoms. Patients presenting with gastrointestinal symptoms were more commonly female, had a longer duration of symptoms before presentation, and had lower 30-day mortality. Patients presenting with confusion, with or without core symptoms, were older and had a higher unadjusted mortality. Symptom sub-phenotypes were highly consistent in replication analysis within the ISARIC-4C study. Similar patterns were externally verified in patients from a study of self-reported symptoms of mild disease. The large scale of the ISARIC-4C study enabled robust, granular discovery and replication. Clinical interpretation is necessary to determine which of these observations have practical utility. We propose that four sub-phenotypes are usefully distinct from the core symptom group: gastro-intestinal disease, productive cough, confusion, and pauci-symptomatic presentations. Importantly, each is associated with an in-hospital mortality which differs from that of patients with core symptoms.

Impact of insufficient sleep on dysregulated blood glucose control under standardised meal conditions.

AIMS/HYPOTHESIS: Sleep, diet and exercise are fundamental to metabolic homeostasis. In this secondary analysis of a repeated measures, nutritional intervention study, we tested whether an individual's sleep quality, duration and timing impact glycaemic response to a breakfast meal the following morning. METHODS: Healthy adults' data (N = 953 [41% twins]) were analysed from the PREDICT dietary intervention trial. Participants consumed isoenergetic standardised meals over 2 weeks in the clinic and at home. Actigraphy was used to assess sleep variables (duration, efficiency, timing) and continuous glucose monitors were used to measure glycaemic variation (>8000 meals). RESULTS: Sleep variables were significantly associated with postprandial glycaemic control (2 h incremental AUC), at both between- and within-person levels. Sleep period time interacted with meal type, with a smaller effect of poor sleep on postprandial blood glucose levels when high-carbohydrate (low fat/protein) (pinteraction = 0.02) and high-fat (pinteraction = 0.03) breakfasts were consumed compared with a reference 75 g OGTT. Within-person sleep period time had a similar interaction (high carbohydrate: pinteraction = 0.001, high fat: pinteraction = 0.02). Within- and between-person sleep efficiency were significantly associated with lower postprandial blood glucose levels irrespective of meal type (both p < 0.03). Later sleep midpoint (time deviation from midnight) was found to be significantly associated with higher postprandial glucose, in both between-person and within-person comparisons (p = 0.035 and p = 0.051, respectively). CONCLUSIONS/INTERPRETATION: Poor sleep efficiency and later bedtime routines are associated with more pronounced postprandial glycaemic responses to breakfast the following morning. A person's deviation from their usual sleep pattern was also associated with poorer postprandial glycaemic control. These findings underscore sleep as a modifiable, non-pharmacological therapeutic target for the optimal regulation of human metabolic health. Trial registration ClinicalTrials.gov NCT03479866.

Estrogen and COVID-19 symptoms: Associations in women from the COVID Symptom Study.

It has been widely observed that adult men of all ages are at higher risk of developing serious complications from COVID-19 when compared with women. This study aimed to investigate the association of COVID-19 positivity and severity with estrogen exposure in women, in a population based matched cohort study of female users of the COVID Symptom Study application in the UK. Analyses included 152,637 women for menopausal status, 295,689 women for exogenous estrogen intake in the form of the combined oral contraceptive pill (COCP), and 151,193 menopausal women for hormone replacement therapy (HRT). Data were collected using the COVID Symptom Study in May-June 2020. Analyses investigated associations between predicted or tested COVID-19 status and menopausal status, COCP use, and HRT use, adjusting for age, smoking and BMI, with follow-up age sensitivity analysis, and validation in a subset of participants from the TwinsUK cohort. Menopausal women had higher rates of predicted COVID-19 (P = 0.003). COCP-users had lower rates of predicted COVID-19 (P = 8.03E-05), with reduction in hospital attendance (P = 0.023). Menopausal women using HRT or hormonal therapies did not exhibit consistent associations, including increased rates of predicted COVID-19 (P = 2.22E-05) for HRT users alone. The findings support a protective effect of estrogen exposure on COVID-19, based on positive association between predicted COVID-19 with menopausal status, and negative association with COCP use. HRT use was positively associated with COVID-19, but the results should be considered with caution due to lack of data on HRT type, route of administration, duration of treatment, and potential unaccounted for confounders and comorbidities.

Diet quality and risk and severity of COVID-19: a prospective cohort study.

OBJECTIVE: Poor metabolic health and unhealthy lifestyle factors have been associated with risk and severity of COVID-19, but data for diet are lacking. We aimed to investigate the association of diet quality with risk and severity of COVID-19 and its interaction with socioeconomic deprivation. DESIGN: We used data from 592 571 participants of the smartphone-based COVID-19 Symptom Study. Diet information was collected for the prepandemic period using a short food frequency questionnaire, and diet quality was assessed using a healthful Plant-Based Diet Score, which emphasises healthy plant foods such as fruits or vegetables. Multivariable Cox models were fitted to calculate HRs and 95% CIs for COVID-19 risk and severity defined using a validated symptom-based algorithm or hospitalisation with oxygen support, respectively. RESULTS: Over 3 886 274 person-months of follow-up, 31 815 COVID-19 cases were documented. Compared with individuals in the lowest quartile of the diet score, high diet quality was associated with lower risk of COVID-19 (HR 0.91; 95% CI 0.88 to 0.94) and severe COVID-19 (HR 0.59; 95% CI 0.47 to 0.74). The joint association of low diet quality and increased deprivation on COVID-19 risk was higher than the sum of the risk associated with each factor alone (Pinteraction=0.005). The corresponding absolute excess rate per 10 000 person/months for lowest vs highest quartile of diet score was 22.5 (95% CI 18.8 to 26.3) among persons living in areas with low deprivation and 40.8 (95% CI 31.7 to 49.8) among persons living in areas with high deprivation. CONCLUSIONS: A diet characterised by healthy plant-based foods was associated with lower risk and severity of COVID-19. This association may be particularly evident among individuals living in areas with higher socioeconomic deprivation.

Genome-wide association study in almost 195,000 individuals identifies 50 previously unidentified genetic loci for eye color.

Human eye color is highly heritable, but its genetic architecture is not yet fully understood. We report the results of the largest genome-wide association study for eye color to date, involving up to 192,986 European participants from 10 populations. We identify 124 independent associations arising from 61 discrete genomic regions, including 50 previously unidentified. We find evidence for genes involved in melanin pigmentation, but we also find associations with genes involved in iris morphology and structure. Further analyses in 1636 Asian participants from two populations suggest that iris pigmentation variation in Asians is genetically similar to Europeans, albeit with smaller effect sizes. Our findings collectively explain 53.2% (95% confidence interval, 45.4 to 61.0%) of eye color variation using common single-nucleotide polymorphisms. Overall, our study outcomes demonstrate that the genetic complexity of human eye color considerably exceeds previous knowledge and expectations, highlighting eye color as a genetically highly complex human trait.

Blue poo: impact of gut transit time on the gut microbiome using a novel marker.

BACKGROUND AND AIMS: Gut transit time is a key modulator of host-microbiome interactions, yet this is often overlooked, partly because reliable methods are typically expensive or burdensome. The aim of this single-arm, single-blinded intervention study is to assess (1) the relationship between gut transit time and the human gut microbiome, and (2) the utility of the 'blue dye' method as an inexpensive and scalable technique to measure transit time. METHODS: We assessed interactions between the taxonomic and functional potential profiles of the gut microbiome (profiled via shotgun metagenomic sequencing), gut transit time (measured via the blue dye method), cardiometabolic health and diet in 863 healthy individuals from the PREDICT 1 study. RESULTS: We found that gut microbiome taxonomic composition can accurately discriminate between gut transit time classes (0.82 area under the receiver operating characteristic curve) and longer gut transit time is linked with specific microbial species such as Akkermansia muciniphila, Bacteroides spp and Alistipes spp (false discovery rate-adjusted p values <0.01). The blue dye measure of gut transit time had the strongest association with the gut microbiome over typical transit time proxies such as stool consistency and frequency. CONCLUSIONS: Gut transit time, measured via the blue dye method, is a more informative marker of gut microbiome function than traditional measures of stool consistency and frequency. The blue dye method can be applied in large-scale epidemiological studies to advance diet-microbiome-health research. Clinical trial registry website https://clinicaltrials.gov/ct2/show/NCT03479866 and trial number NCT03479866.

A GWAS in Latin Americans identifies novel face shape loci, implicating VPS13B and a Denisovan introgressed region in facial variation.

To characterize the genetic basis of facial features in Latin Americans, we performed a genome-wide association study (GWAS) of more than 6000 individuals using 59 landmark-based measurements from two-dimensional profile photographs and ~9,000,000 genotyped or imputed single-nucleotide polymorphisms. We detected significant association of 32 traits with at least 1 (and up to 6) of 32 different genomic regions, more than doubling the number of robustly associated face morphology loci reported until now (from 11 to 23). These GWAS hits are strongly enriched in regulatory sequences active specifically during craniofacial development. The associated region in 1p12 includes a tract of archaic adaptive introgression, with a Denisovan haplotype common in Native Americans affecting particularly lip thickness. Among the nine previously unidentified face morphology loci we identified is the VPS13B gene region, and we show that variants in this region also affect midfacial morphology in mice.

Strain-Level Analysis of Bifidobacterium spp. from Gut Microbiomes of Adults with Differing Lactase Persistence Genotypes.

One of the strongest associations between human genetics and the gut microbiome is a greater relative abundance of Bifidobacterium in adults with lactase gene (LCT) single nucleotide polymorphisms (SNPs) associated with lactase nonpersistence (GG genotypes), versus lactase persistence (AA/AG genotypes). To gain a finer-grained phylogenetic resolution of this association, we interrogated 1,680 16S rRNA libraries and 245 metagenomes from gut microbiomes of adults with various lactase persistence genotypes. We further employed a novel genome-capture-based enrichment of Bifidobacterium DNA from a subset of these metagenomes, including monozygotic (MZ) twin pairs, each sampled 2 or 3 times. B. adolescentis and B. longum were the most abundant Bifidobacterium species regardless of host LCT genotype. LCT genotypes could not be discriminated based on relative abundances of Bifidobacterium species or Bifidobacterium community structure. Three distinct metagenomic analysis methods of Bifidobacterium-enriched DNA revealed intraindividual temporal stability of B. longum, B. adolescentis, and B. bifidum strains against the background of a changeable microbiome. Two of our three methods also observed greater strain sharing within MZ twin pairs than within unrelated individuals for B. adolescentis, while no method revealed an effect of host LCT genotype on Bifidobacterium strain composition. Our results support a "rising tide lifts all boats" model for the dominant bifidobacteria in the adult gut: their higher abundance in lactase-nonpersistent than in lactase-persistent individuals results from an expansion at the genus level. Bifidobacterium species are known to be transmitted from mother to child and stable within individuals in infancy and childhood: our results extend this stability into adulthood.IMPORTANCE When humans domesticated animals, some adapted genetically to digest milk into adulthood (lactase persistence). The gut microbiomes of people with lactase-persistent genotypes (AA or AG) differ from those with lactase-nonpersistent genotypes (GG) by containing fewer bacteria belonging to the bifidobacteria, a group which contains beneficial species. Here, we asked if the gut microbiomes of adults with GG and AA/AG genotypes differ in the species of bifidobacteria present. In particular, we used a novel technique which allowed us to compare bifidobacteria in adults at the strain level, without the traditional need for culturing. Our results show that the GG genotype enhances the abundance of bifidobacteria regardless of species. We also noted that a person's specific strains are recoverable several years later, and twins can share the same ones. Given that bifidobacteria are inherited from mother to child, strain stability over time in adulthood suggests long-term, multigenerational inheritance.

The association between low birth weight, childhood recollections of parental response to illness, and irritable bowel syndrome: a twin study.

BACKGROUND: The aetiology of irritable bowel syndrome (IBS) is multifactorial, including genetic and environmental factors. Previous studies have suggested that low birth weight and family environment during childhood are associated with developing IBS. METHODS: A survey was sent to all individuals in a UK twin registry. Questions included IBS diagnosed by the Rome IV criteria and if a doctor had previously diagnosed them with IBS. Subjects were categorized as having IBS by Rome IV criteria, a medical diagnosis of IBS or no IBS. Further questions included subjects' recollections of their parents' responses to illness in both the respondent as a child and in the parents themselves. Information regarding birth weight and gestational age have been collected previously. KEY RESULTS: 4258 subjects responded to the questionnaire (51.7%), mean age of 52 (SD 14) years, of whom 98.5% were white and 89.6% female. The mean birth weight was 2.4  (0.6) kg. 5.1% satisfied the Rome IV IBS criteria, the same prevalence as the UK population. However, 14.1% had a previous medical diagnosis of IBS. There was no association found between birth weight and IBS or a medical diagnosis of IBS. On multivariable regression analysis, including parental responses to illness, subjects recalling a parent responding to the parent's bowel symptoms by excusing themselves from household chores were associated with a Rome IV diagnosis of IBS (OR 2.19 (95% CI 1.17-4.10), P = .013). CONCLUSIONS AND INFERENCES: There was no association between birth weight and IBS. However, observing their parents excuse themselves from household chores when they had bowel symptoms was associated with IBS in later life.

Genome-wide scan identifies novel genetic loci regulating salivary metabolite levels.

Saliva, as a biofluid, is inexpensive and non-invasive to obtain, and provides a vital tool to investigate oral health and its interaction with systemic health conditions. There is growing interest in salivary biomarkers for systemic diseases, notably cardiovascular disease. Whereas hundreds of genetic loci have been shown to be involved in the regulation of blood metabolites, leading to significant insights into the pathogenesis of complex human diseases, little is known about the impact of host genetics on salivary metabolites. Here we report the first genome-wide association study exploring 476 salivary metabolites in 1419 subjects from the TwinsUK cohort (discovery phase), followed by replication in the Study of Health in Pomerania (SHIP-2) cohort. A total of 14 distinct locus-metabolite associations were identified in the discovery phase, most of which were replicated in SHIP-2. While only a limited number of the loci that are known to regulate blood metabolites were also associated with salivary metabolites in our study, we identified several novel saliva-specific locus-metabolite associations, including associations for the AGMAT (with the metabolites 4-guanidinobutanoate and beta-guanidinopropanoate), ATP13A5 (with the metabolite creatinine) and DPYS (with the metabolites 3-ureidopropionate and 3-ureidoisobutyrate) loci. Our study suggests that there may be regulatory pathways of particular relevance to the salivary metabolome. In addition, some of our findings may have clinical significance, such as the utility of the pyrimidine (uracil) degradation metabolites in predicting 5-fluorouracil toxicity and the role of the agmatine pathway metabolites as biomarkers of oral health.

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.