Complex trait associations in rare diseases and impacts on Mendelian variant interpretation.

Emerging evidence implicates common genetic variation - aggregated into polygenic scores (PGS) - in the onset and phenotypic presentation of rare diseases. Here, we comprehensively map individual polygenic liability for 1102 open-source PGS in a cohort of 3059 probands enrolled in the Genomic Answers for Kids (GA4K) rare disease study, revealing widespread associations between rare disease phenotypes and PGSs for common complex diseases and traits, blood protein levels, and brain and other organ morphological measurements. Using this resource, we demonstrate increased polygenic liability in probands with an inherited candidate disease variant (VUS) compared to unaffected carrier parents. Further, we show an enrichment for large-effect rare variants in putative core PGS genes for associated complex traits. Overall, our study supports and expands on previous findings of complex trait associations in rare diseases, implicates polygenic liability as a potential mechanism underlying variable penetrance of candidate causal variants, and provides a framework for identifying novel candidate rare disease genes.

Genomic insights into pediatric intestinal inflammatory and eosinophilic disorders using single-cell RNA-sequencing.

Introduction: Chronic inflammation of the gastrointestinal tissues underlies gastrointestinal inflammatory disorders, leading to tissue damage and a constellation of painful and debilitating symptoms. These disorders include inflammatory bowel diseases (Crohn's disease and ulcerative colitis), and eosinophilic disorders (eosinophilic esophagitis and eosinophilic duodenitis). Gastrointestinal inflammatory disorders can often present with overlapping symptoms necessitating the use of invasive procedures to give an accurate diagnosis. Methods: This study used peripheral blood mononuclear cells from individuals with Crohn's disease, ulcerative colitis, eosinophilic esophagitis, and eosinophilic duodenitis to better understand the alterations to the transcriptome of individuals with these diseases and identify potential markers of active inflammation within the peripheral blood of patients that may be useful in diagnosis. Single-cell RNA-sequencing was performed on peripheral blood mononuclear cells isolated from the blood samples of pediatric patients diagnosed with gastrointestinal disorders, including Crohn's disease, ulcerative colitis, eosinophilic esophagitis, eosinophilic duodenitis, and controls with histologically healthy gastrointestinal tracts. Results: We identified 730 (FDR < 0.05) differentially expressed genes between individuals with gastrointestinal disorders and controls across eight immune cell types. Discussion: There were common patterns among GI disorders, such as the widespread upregulation of MTRNR2L8 across cell types, and many differentially expressed genes showed distinct patterns of dysregulation among the different gastrointestinal diseases compared to controls, including upregulation of XIST across cell types among individuals with ulcerative colitis and upregulation of Th2-associated genes in eosinophilic disorders. These findings indicate both overlapping and distinct alterations to the transcriptome of individuals with gastrointestinal disorders compared to controls, which provide insight as to which genes may be useful as markers for disease in the peripheral blood of patients.

Utilization of geospatial distribution in the measurement of study cohort representativeness.

OBJECTIVE: The ability to apply results from a study to a broader population remains a primary objective in translational science. Distinct from intrinsic elements of scientific rigor, the extrinsic concept of generalization requires there be alignment between a study cohort and population in which results are expected to be applied. Widespread efforts have been made to quantify representativeness of study cohorts. These techniques, however, often consider the study and target cohorts as monolithic collections that can be directly compared. Overlooking known impacts to health from socio-demographic and environmental factors tied to individual's geographical location, and potentially obfuscating misalignment in underrepresented population subgroups. This manuscript introduces several measures to account for geographic information in the assessment of cohort representation. METHODS: Metrics were defined across two themes: First, measures of recruitment, to assess if a study cohort is drawn at an expected rate and in an expected geographical pattern with respect to individuals in a reference cohort. Second, measures of individual characteristics, to assess if the individuals in the study cohort accurately reflect the sociodemographic, clinical, and geographic diversity observed across a reference cohort while accounting for the geospatial proximity of individuals. RESULTS: As an empirical demonstration, methods are applied to an active clinical study examining asthma in Black/African American patients at a US Midwestern pediatric hospital. Results illustrate how areas of over- and under-recruitment can be identified and contextualized in light of study recruitment patterns at an individual-level, highlighting the ability to identify a subset of features for which the study cohort closely resembled the broader population. In addition they provide an opportunity to dive deeper into misalignments, to identify study cohort members that are in some way distinct from the communities for which they are expected to represent. CONCLUSION: Together, these metrics provide a comprehensive spatial assessment of a study cohort with respect to a broader target population. Such an approach offers researchers a toolset by which to target expected generalization of results derived from a given study.

Morphometric Analysis of Spina Bifida after Fetal Repair Shows New Subtypes with Associated Outcomes.

Importance: The binary classification of spina bifida lesions as myelomeningocele (with sac) or myeloschisis (without sac) belies a spectrum of morphologies, which have not been correlated to clinical characteristics and outcomes. Objective: To characterize spina bifida lesion types and correlate them with preoperative presentation and postoperative outcomes. Design: Secondary analysis of images and videos obtained during fetoscopic spina bifida repair surgery from 2020-2023. Setting: Fetal surgery was performed at a quaternary care center. Participants: A prospective cohort of patients referred for fetal spina bifida underwent fetoscopic repair under an FDA-approved protocol. Of 60 lesions repaired, 57 had available images and were included in the analysis. Interventions or Exposures: We evaluated lesion morphology on high-resolution intraoperative images and videos to categorize lesions based on placode exposure and nerve root stretching. Main Outcomes and Measures: The reproducibility of the lesion classification was assessed via Kappa interrater agreement. Preoperative characteristics analyzed include ventricle size, tonsillar herniation level, lower extremities movement, and lesion dimensions. Outcomes included surgical time, need for patch for skin closure, gestational age at delivery, preterm premature rupture of membranes (PPROM), and neonatal cerebrospinal fluid (CSF) diversion. Results: We distinguished five lesion types that differ across a range of sac sizes, nerve root stretching, and placode exposure, with 93% agreement between examiners (p<0.001). Fetal characteristics at preoperative evaluation differed significantly by lesion type, including lesion volume (p<0.001), largest ventricle size (p=0.008), tonsillar herniation (p=0.005), and head circumference (p=0.03). Lesion level, talipes, and lower extremities movement did not differ by type. Surgical and perinatal outcomes differed by lesion type, including need for patch skin closure (p<0.001), gestational age at delivery (p=0.01), and NICU length of stay (p<0.001). PPROM, CSF leakage at birth, and CSF diversion in the NICU did not differ between lesion groups. Linear regression associated severity of ventriculomegaly with lesion type, but not with tonsillar herniation level. Conclusions and Relevance: There is a distinct phenotypic spectrum in open spina bifida with differential baseline presentation and outcomes. Severity of ventriculomegaly is associated with lesion type, rather than tonsillar herniation level. Our findings expand the classification of spina bifida to reveal a spectrum that warrants further study.

Enhanced resolution profiling in twins reveals differential methylation signatures of type 2 diabetes with links to its complications.

BACKGROUND: Type 2 diabetes (T2D) susceptibility is influenced by genetic and environmental factors. Previous findings suggest DNA methylation as a potential mechanism in T2D pathogenesis and progression. METHODS: We profiled DNA methylation in 248 blood samples from participants of European ancestry from 7 twin cohorts using a methylation sequencing platform targeting regulatory genomic regions encompassing 2,048,698 CpG sites. FINDINGS: We find and replicate 3 previously unreported T2D differentially methylated CpG positions (T2D-DMPs) at FDR 5% in RGL3, NGB and OTX2, and 20 signals at FDR 25%, of which 14 replicated. Integrating genetic variation and T2D-discordant monozygotic twin analyses, we identify both genetic-based and genetic-independent T2D-DMPs. The signals annotate to genes with established GWAS and EWAS links to T2D and its complications, including blood pressure (RGL3) and eye disease (OTX2). INTERPRETATION: The results help to improve our understanding of T2D disease pathogenesis and progression and may provide biomarkers for its complications. FUNDING: Funding acknowledgements for each cohort can be found in the Supplementary Note.

Complex trait associations in rare diseases and impacts on Mendelian variant interpretation.

Emerging evidence implicates common genetic variation - aggregated into polygenic scores (PGS) - impacting the onset and phenotypic presentation of rare diseases. In this study, we quantified individual polygenic liability for 1,151 previously published PGS in a cohort of 2,374 probands enrolled in the Genomic Answers for Kids (GA4K) rare disease study, revealing widespread associations between rare disease phenotypes and PGSs for common complex diseases and traits, blood protein levels, and brain and other organ morphological measurements. We observed increased polygenic burden in probands with variants of unknown significance (VUS) compared to unaffected carrier parents. We further observed an enrichment in overlap between diagnostic and candidate rare disease genes and large-effect PGS genes. Overall, our study supports and expands on previous findings of complex trait associations in rare disease phenotypes and provides a framework for identifying novel candidate rare disease genes and in understanding variable penetrance of candidate Mendelian disease variants.

DNA Methylation at ABCG1 and Long-term Changes in Adiposity and Fat Distribution in Response to Dietary Interventions: The POUNDS Lost Trial.

OBJECTIVE: To examine whether participants with different levels of diabetes-related DNA methylation at ABCG1 might respond differently to dietary weight loss interventions with long-term changes in adiposity and body fat distribution. RESEARCH DESIGN AND METHODS: The current study included overweight/obese participants from the POUNDS Lost trial. Blood levels of regional DNA methylation at ABCG1 were profiled by high-resolution methylC-capture sequencing at baseline among 673 participants, of whom 598 were followed up at 6 months and 543 at 2 years. Two-year changes in adiposity and computed tomography-measured body fat distribution were calculated. RESULTS: Regional DNA methylation at ABCG1 showed significantly different associations with long-term changes in body weight and waist circumference at 6 months and 2 years in dietary interventions varying in protein intake (interaction P < 0.05 for all). Among participants assigned to an average-protein (15%) diet, lower baseline regional DNA methylation at ABCG1 was associated with greater reductions in body weight and waist circumference at 6 months and 2 years, whereas opposite associations were found among those assigned to a high-protein (25%) diet. Similar interaction patterns were also observed for body fat distribution, including visceral adipose tissue, subcutaneous adipose tissue, deep subcutaneous adipose tissue, and total adipose tissue at 6 months and 2 years (interaction P < 0.05 for all). CONCLUSIONS: Baseline DNA methylation at ABCG1 interacted with dietary protein intake on long-term decreases in adiposity and body fat distribution. Participants with lower methylation at ABCG1 benefitted more in long-term reductions in body weight, waist circumference, and body fat distribution when consuming an average-protein diet.

Leveraging osteoclast genetic regulatory data to identify genes with a role in osteoarthritis.

There has been a growing interest in the role of the subchondral bone and its resident osteoclasts in the progression of osteoarthritis (OA). A recent genome-wide association study (GWAS) identified 100 independent association signals for OA traits. Most of these signals are led by noncoding variants, suggesting that genetic regulatory effects may drive many of the associations. We have generated a unique human osteoclast-like cell-specific expression quantitative trait locus (eQTL) resource for studying the genetics of bone disease. Considering the potential role of osteoclasts in the pathogenesis of OA, we performed an integrative analysis of this dataset with the recently published OA GWAS results. Summary data-based Mendelian randomization (SMR) and colocalization analyses identified 38 genes with a potential role in OA, including some that have been implicated in Mendelian diseases with joint/skeletal abnormalities, such as BICRA, EIF6, CHST3, and FBN2. Several OA GWAS signals demonstrated colocalization with more than one eQTL peak, including at 19q13.32 (hip OA with BCAM, PRKD2, and BICRA eQTL). We also identified a number of eQTL signals colocalizing with more than one OA trait, including FAM53A, GCAT, HMGN1, MGAT4A, RRP7BP, and TRIOBP. An SMR analysis identified 3 loci with evidence of pleiotropic effects on OA-risk and gene expression: LINC01481, CPNE1, and EIF6. Both CPNE1 and EIF6 are located at 20q11.22, a locus harboring 2 other strong OA candidate genes, GDF5 and UQCC1, suggesting the presence of an OA-risk gene cluster. In summary, we have used our osteoclast-specific eQTL dataset to identify genes potentially involved with the pathogenesis of OA.

Extravillous trophoblast cell lineage development is associated with active remodeling of the chromatin landscape.

The extravillous trophoblast cell lineage is a key feature of placentation and successful pregnancy. Knowledge of transcriptional regulation driving extravillous trophoblast cell development is limited. Here, we map the transcriptome and epigenome landscape as well as chromatin interactions of human trophoblast stem cells and their transition into extravillous trophoblast cells. We show that integrating chromatin accessibility, long-range chromatin interactions, transcriptomic, and transcription factor binding motif enrichment enables identification of transcription factors and regulatory mechanisms critical for extravillous trophoblast cell development. We elucidate functional roles for TFAP2C, SNAI1, and EPAS1 in the regulation of extravillous trophoblast cell development. EPAS1 is identified as an upstream regulator of key extravillous trophoblast cell transcription factors, including ASCL2 and SNAI1 and together with its target genes, is linked to pregnancy loss and birth weight. Collectively, we reveal activation of a dynamic regulatory network and provide a framework for understanding extravillous trophoblast cell specification in trophoblast cell lineage development and human placentation.

Direct haplotype-resolved 5-base HiFi sequencing for genome-wide profiling of hypermethylation outliers in a rare disease cohort.

Long-read HiFi genome sequencing allows for accurate detection and direct phasing of single nucleotide variants, indels, and structural variants. Recent algorithmic development enables simultaneous detection of CpG methylation for analysis of regulatory element activity directly in HiFi reads. We present a comprehensive haplotype resolved 5-base HiFi genome sequencing dataset from a rare disease cohort of 276 samples in 152 families to identify rare (~0.5%) hypermethylation events. We find that 80% of these events are allele-specific and predicted to cause loss of regulatory element activity. We demonstrate heritability of extreme hypermethylation including rare cis variants associated with short (~200 bp) and large hypermethylation events (>1 kb), respectively. We identify repeat expansions in proximal promoters predicting allelic gene silencing via hypermethylation and demonstrate allelic transcriptional events downstream. On average 30-40 rare hypermethylation tiles overlap rare disease genes per patient, providing indications for variation prioritization including a previously undiagnosed pathogenic allele in DIP2B causing global developmental delay. We propose that use of HiFi genome sequencing in unsolved rare disease cases will allow detection of unconventional diseases alleles due to loss of regulatory element activity.

DNA Methylation of Birthweight-Blood Pressure Genes and Changes of Blood Pressure in Response to Weight-Loss Diets in the POUNDS Lost Trial.

BACKGROUND: DNA methylation (DNAm) may play a critical role in bridging prenatal adverse events and cardiometabolic disorders including hypertension in later life. METHODS: We included 672 adult participants with overweight or obesity, who participated in a 2-year randomized weight-loss dietary intervention study. We defined the regional DNAm levels as the average methylation level of 5'-cytosine-phosphate-guanine-3' within 500 bp of LINC00319 (cg01820192), ATP2B1 (cg00508575), and LMNA (cg12593793), respectively. Generalized linear regression models were used to assess the association between the regional DNAm and 2-year blood pressure changes. Trajectory analysis was used to identify subgroups that shared a similar underlying trajectory of 2-year blood pressure changes. RESULTS: The regional DNAm at LINC00319, showed significantly different associations with 2-year changes in systolic blood pressure and diastolic blood pressure among participants assigned to low- or high-fat diets (P for interaction<0.05 for all). In response to the low-fat diet, per SD higher regional DNAm at LINC00319 was associated with greater reductions in both 2-year changes in systolic blood pressure (ß, -1.481; P=0.020) and diastolic blood pressure (ß, -1.096; P=0.009). Three trajectories of changes in systolic blood pressure or diastolic blood pressure were identified, and participants with higher regional DNAm at LINC00319 were more likely to experience and maintain decreased systolic blood pressure and diastolic blood pressure (odds ratio of being in decrease-stable versus stable [95% CI], 1.542 [1.146-2.076] and 1.463 [1.125-1.902]). CONCLUSIONS: Our findings suggest that DNAm could be a metabolic memory bridging early and later life, and an indicator of more benefits from eating a low-fat weight-loss diet.

Rare and Common Variants in GALNT3 May Affect Bone Mass Independently of Phosphate Metabolism.

Anabolic treatment options for osteoporosis remain limited. One approach to discovering novel anabolic drug targets is to identify genetic causes of extreme high bone mass (HBM). We investigated a pedigree with unexplained HBM within the UK HBM study, a national cohort of probands with HBM and their relatives. Whole exome sequencing (WES) in a family with HBM identified a rare heterozygous missense variant (NM_004482.4:c.1657C > T, p.Arg553Trp) in GALNT3, segregating appropriately. Interrogation of data from the UK HBM study and the Anglo-Australasian Osteoporosis Genetics Consortium (AOGC) revealed an unrelated individual with HBM with another rare heterozygous variant (NM_004482.4:c.831 T > A, p.Asp277Glu) within the same gene. In silico protein modeling predicted that p.Arg553Trp would disrupt salt-bridge interactions, causing instability of GALNT3, and that p.Asp277Glu would disrupt manganese binding and consequently GALNT3 catalytic function. Bi-allelic loss-of-function GALNT3 mutations alter FGF23 metabolism, resulting in hyperphosphatemia and causing familial tumoral calcinosis (FTC). However, bone mineral density (BMD) in FTC cases, when reported, has been either normal or low. Common variants in the GALNT3 locus show genome-wide significant associations with lumbar, femoral neck, and total body BMD. However, no significant associations with BMD are observed at loci coding for FGF23, its receptor FGFR1, or coreceptor klotho. Mendelian randomization analysis, using expression quantitative trait loci (eQTL) data from primary human osteoblasts and genome-wide association studies data from UK Biobank, suggested increased expression of GALNT3 reduces total body, lumbar spine, and femoral neck BMD but has no effect on phosphate concentrations. In conclusion, rare heterozygous loss-of-function variants in GALNT3 may cause HBM without altering phosphate concentration. These findings suggest that GALNT3 may affect BMD through pathways other than FGF23 regulation, the identification of which may yield novel anabolic drug targets for osteoporosis. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

DNA Methylation Near CPT1A and Changes in Triglyceride-rich Lipoproteins in Response to Weight-loss Diet Interventions.

CONTEXT: Carnitine palmitoyltransferase-1A, encoded by the CPT1A gene, plays a key role in the oxidation of long-chain fatty acids in the mitochondria and may be important in triglyceride metabolism. Previous work has shown that high fat intake was negatively associated with CPT1A methylation and positively associated with CPT1A expression. OBJECTIVE: We aim to investigate the association of DNA methylation (DNAm) at the CPT1A gene with reductions in triglycerides and triglyceride-rich lipoproteins (TRLs) in response to weight-loss diet interventions. METHODS: The current study included 538 White participants, who were randomly assigned to 1 of 4 diets varying in macronutrient components. We defined the regional DNAm at CPT1A as the average methylation level over CpGs within 500 bp of the 3 triglyceride-related DNAm sites. RESULTS: Dietary fat intake significantly modified the association between baseline DNAm at CPT1A and 2-year changes in total plasma triglycerides, independent of concurrent weight loss. Among participants assigned to a low-fat diet, a higher regional DNAm level at CPT1A was associated with a greater reduction in total plasma triglycerides at 2 years (P = .01), compared with those assigned to a high-fat diet (P = .64) (P interaction = .018). Further investigation on lipids and apolipoproteins in very low-density lipoprotein (VLDL) revealed similar interaction patterns for 2-year changes in VLDL-triglycerides, VLDL-cholesterol, and VLDL-apolipoprotein B (P interaction = .009, .002, and .016, respectively), but not for VLDL-apoC-III (P interaction = .36). CONCLUSION: Participants with a higher regional DNAm level at CPT1A benefit more in long-term improvement in triglycerides, particularly in the TRLs and related apolipoproteins when consuming a low-fat weight-loss diet.

Insurance denials and diagnostic rates in a pediatric genomic research cohort.

PURPOSE: This study aimed to assess the amount and types of clinical genetic testing denied by insurance and the rate of diagnostic and candidate genetic findings identified through research in patients who faced insurance denials. METHODS: Analysis consisted of review of insurance denials in 801 patients enrolled in a pediatric genomic research repository with either no previous genetic testing or previous negative genetic testing result identified through cross-referencing with insurance prior-authorizations in patient medical records. Patients and denials were also categorized by type of insurance coverage. Diagnostic findings and candidate genetic findings in these groups were determined through review of our internal variant database and patient charts. RESULTS: Of the 801 patients analyzed, 147 had insurance prior-authorization denials on record (18.3%). Exome sequencing and microarray were the most frequently denied genetic tests. Private insurance was significantly more likely to deny testing than public insurance (odds ratio = 2.03 [95% CI = 1.38-2.99] P = .0003). Of the 147 patients with insurance denials, 53.7% had at least 1 diagnostic or candidate finding and 10.9% specifically had a clinically diagnostic finding. Fifty percent of patients with clinically diagnostic results had immediate medical management changes (5.4% of all patients experiencing denials). CONCLUSION: Many patients face a major barrier to genetic testing in the form of lack of insurance coverage. A number of these patients have clinically diagnostic findings with medical management implications that would not have been identified without access to research testing. These findings support re-evaluation of insurance carriers' coverage policies.

Neonatal hyperoxia induces activated pulmonary cellular states and sex-dependent transcriptomic changes in a model of experimental bronchopulmonary dysplasia.

Hyperoxia disrupts lung development in mice and causes bronchopulmonary dysplasia (BPD) in neonates. To investigate sex-dependent molecular and cellular programming involved in hyperoxia, we surveyed the mouse lung using single cell RNA sequencing (scRNA-seq), and validated our findings in human neonatal lung cells in vitro. Hyperoxia-induced inflammation in alveolar type (AT) 2 cells gave rise to damage-associated transient progenitors (DATPs). It also induced a new subpopulation of AT1 cells with reduced expression of growth factors normally secreted by AT1 cells, but increased mitochondrial gene expression. Female alveolar epithelial cells had less EMT and pulmonary fibrosis signaling in hyperoxia. In the endothelium, expansion of Car4+ EC (Cap2) was seen in hyperoxia along with an emergent subpopulation of Cap2 with repressed VEGF signaling. This regenerative response was increased in females exposed to hyperoxia. Mesenchymal cells had inflammatory signatures in hyperoxia, with a new distal interstitial fibroblast subcluster characterized by repressed lipid biosynthesis and a transcriptomic signature resembling myofibroblasts. Hyperoxia-induced gene expression signatures in human neonatal fibroblasts and alveolar epithelial cells in vitro resembled mouse scRNA-seq data. These findings suggest that neonatal exposure to hyperoxia programs distinct sex-specific stem cell progenitor and cellular reparative responses that underpin lung remodeling in BPD.

The Role of Dynamic DNA Methylation in Liver Transplant Rejection in Children.

Transcriptional regulation of liver transplant (LT) rejection may reveal novel predictive and therapeutic targets. The purpose of this article is to test the role of differential DNA methylation in children with biopsy-proven acute cellular rejection after LT. Methods: Paired peripheral blood DNA samples were obtained before and after LT from 17 children, including 4 rejectors (Rs) and 13 nonrejectors (NRs), and assayed with MethylC capture sequencing approach covering 5 million CpGs in immune-cell-specific regulatory elements. Differentially methylated CpGs (DMCs) were identified using generalized linear regression models adjusting for sex and age and merged into differentially methylated regions (DMRs) comprising 3 or more DMCs. Results: Contrasting Rs versus NRs, we identified 2238 DMCs in post-LT and 2620 DMCs in pre-LT samples, which clustered in 216 and 282 DMRs, respectively. DMCs associated with R were enriched in enhancers and depleted in promoters. Among DMRs, the proportion of hypomethylated DMRs increased from 61/282 (22%) in pre-LT to 103/216 (48%, P < 0.0001) in post-LT samples. The highest-ranked biological processes enriched in post-LT DMCs were antigen processing and presentation via major histocompatibility complex (MHC) class I, MHC class I complex, and peptide binding (P < 7.92 × 10-17), respectively. Top-ranked DMRs mapped to genes that mediate B-cell receptor signaling (ADAP1) or regulate several immune cells (ARRB2) (P < 3.75 × 10-08). DMRs in MHC class I genes were enriched for single nucleotide polymorphisms (SNPs), which bind transcription factors, affect gene expression and splicing, or alter peptide-binding amino acid sequences. Conclusions: Dynamic methylation in distal regulatory regions reveals known transplant-relevant MHC-dependent rejection pathways and identifies novel loci for future mechanistic evaluations in pediatric transplant subcohorts.

Independent phenotypic plasticity axes define distinct obesity sub-types.

Studies in genetically 'identical' individuals indicate that as much as 50% of complex trait variation cannot be traced to genetics or to the environment. The mechanisms that generate this 'unexplained' phenotypic variation (UPV) remain largely unknown. Here, we identify neuronatin (NNAT) as a conserved factor that buffers against UPV. We find that Nnat deficiency in isogenic mice triggers the emergence of a bi-stable polyphenism, where littermates emerge into adulthood either 'normal' or 'overgrown'. Mechanistically, this is mediated by an insulin-dependent overgrowth that arises from histone deacetylase (HDAC)-dependent ß-cell hyperproliferation. A multi-dimensional analysis of monozygotic twin discordance reveals the existence of two patterns of human UPV, one of which (Type B) phenocopies the NNAT-buffered polyphenism identified in mice. Specifically, Type-B monozygotic co-twins exhibit coordinated increases in fat and lean mass across the body; decreased NNAT expression; increased HDAC-responsive gene signatures; and clinical outcomes linked to insulinemia. Critically, the Type-B UPV signature stratifies both childhood and adult cohorts into four metabolic states, including two phenotypically and molecularly distinct types of obesity.

Large-scale analysis of circulating glutamate and adipose gene expression in relation to abdominal obesity.

Circulating levels of the amino acid glutamate are associated with central fat accumulation, yet the pathophysiology of this relationship remains unknown. We aimed to (i) refine and validate the association between circulating glutamate and abdominal obesity in a large twin cohort, and (ii) investigate whether transcriptomic profiles in adipose tissue could provide insight into the biological mechanisms underlying the association. First, in a cohort of 4665 individuals from the TwinsUK resource, we identified individuals with abdominal obesity and compared prevalence of the latter across circulating glutamate quintiles. Second, we used transcriptomic signatures generated from adipose tissue, both subcutaneous and visceral, to investigate associations with circulating glutamate levels. Individuals in the top circulating glutamate quintile had a sevenfold higher prevalence of abdominal obesity compared to those in the bottom quintile. The adipose tissue transcriptomic analyses identified GLUL, encoding Glutamate-Ammonia Ligase, as being associated with circulating glutamate and abdominal obesity, with pronounced signatures in the visceral depot. In conclusion, circulating glutamate is positively associated with the prevalence of abdominal obesity which relates to dysregulated GLUL expression specifically in visceral adipose tissue.

Adipose methylome integrative-omic analyses reveal genetic and dietary metabolic health drivers and insulin resistance classifiers.

BACKGROUND: There is considerable evidence for the importance of the DNA methylome in metabolic health, for example, a robust methylation signature has been associated with body mass index (BMI). However, visceral fat (VF) mass accumulation is a greater risk factor for metabolic disease than BMI alone. In this study, we dissect the subcutaneous adipose tissue (SAT) methylome signature relevant to metabolic health by focusing on VF as the major risk factor of metabolic disease. We integrate results with genetic, blood methylation, SAT gene expression, blood metabolomic, dietary intake and metabolic phenotype data to assess and quantify genetic and environmental drivers of the identified signals, as well as their potential functional roles. METHODS: Epigenome-wide association analyses were carried out to determine visceral fat mass-associated differentially methylated positions (VF-DMPs) in SAT samples from 538 TwinsUK participants. Validation and replication were performed in 333 individuals from 3 independent cohorts. To assess functional impacts of the VF-DMPs, the association between VF and gene expression was determined at the genes annotated to the VF-DMPs and an association analysis was carried out to determine whether methylation at the VF-DMPs is associated with gene expression. Further epigenetic analyses were carried out to compare methylation levels at the VF-DMPs as the response variables and a range of different metabolic health phenotypes including android:gynoid fat ratio (AGR), lipids, blood metabolomic profiles, insulin resistance, T2D and dietary intake variables. The results from all analyses were integrated to identify signals that exhibit altered SAT function and have strong relevance to metabolic health. RESULTS: We identified 1181 CpG positions in 788 genes to be differentially methylated with VF (VF-DMPs) with significant enrichment in the insulin signalling pathway. Follow-up cross-omic analysis of VF-DMPs integrating genetics, gene expression, metabolomics, diet, and metabolic traits highlighted VF-DMPs located in 9 genes with strong relevance to metabolic disease mechanisms, with replication of signals in FASN, SREBF1, TAGLN2, PC and CFAP410. PC methylation showed evidence for mediating effects of diet on VF. FASN DNA methylation exhibited putative causal effects on VF that were also strongly associated with insulin resistance and methylation levels in FASN better classified insulin resistance (AUC=0.91) than BMI or VF alone. CONCLUSIONS: Our findings help characterise the adiposity-associated methylation signature of SAT, with insights for metabolic disease risk.

Cross-reactive antibodies elicited to conserved epitopes on SARS-CoV-2 spike protein after infection and vaccination.

SARS-CoV-2 is a novel betacoronavirus that caused coronavirus disease 2019 and has resulted in millions of deaths worldwide. Novel coronavirus infections in humans have steadily become more common. Understanding antibody responses to SARS-CoV-2, and identifying conserved, cross-reactive epitopes among coronavirus strains could inform the design of vaccines and therapeutics with broad application. Here, we determined that individuals with previous SARS-CoV-2 infection or vaccinated with the Pfizer-BioNTech BNT162b2 vaccine produced antibody responses that cross-reacted with related betacoronaviruses. Moreover, we designed a peptide-conjugate vaccine with a conserved SARS-CoV-2 S2 spike epitope, immunized mice and determined cross-reactive antibody binding to SARS-CoV-2 and other related coronaviruses. This conserved spike epitope also shared sequence homology to proteins in commensal gut microbiota and could prime immune responses in humans. Thus, SARS-CoV-2 conserved epitopes elicit cross-reactive immune responses to both related coronaviruses and host bacteria that could serve as future targets for broad coronavirus therapeutics and vaccines.