The UK Adult Twin Registry (TwinsUK Resource)

Insights into non-crossover recombination from long-read sperm sequencing

Meiotic recombination is a fundamental process that generates genetic diversity by creating new combinations of existing alleles. Although human crossovers have been studied at the pedigree, population and single-cell level, the more frequent non-crossover events that lead to gene conversion are harder to study, particularly at the individual level. Here we show that single high-fidelity long sequencing reads from sperm can capture both crossovers and non-crossovers, allowing effectively arbitrary sample sizes for analysis from one male. Using fifteen sperm samples from thirteen donors we demonstrate variation between and within donors for the rates of different types of recombination. Intriguingly, we observe a tendency for non-crossover gene conversions to occur upstream of nearby PRDM9 binding sites, whereas crossover locations have a slight downstream bias. We further provide evidence for two distinct non-crossover processes. One gives rise to the vast majority of non-crossovers with mean conversion tract length under 50bp, which we suggest is an outcome of standard PRDM9-induced meiotic recombination. In contrast ~2% of non-crossovers have much longer mean tract length, and potentially originate from the same process as complex events with more than two haplotype switches, which is not associated with PRDM9 binding sites and is also seen in somatic cells.

Exploring the causal association between serum metabolites and erectile dysfunction: a bidirectional Mendelian randomisation study

Erectile dysfunction is a common sexual disorder in men. Some studies have found a strong association between some serum metabolites and erectile dysfunction. To investigate this association further, we used bidirectional Mendelian randomisation to investigate causality and possible biological mechanisms.Firstly, this study screened the statistics of genome-wide association studies of serum metabolites and erectile dysfunction to obtain instrumental variables. Inverse variance weighting was used as the primary method for causal effect analysis of instrumental variables in forward or reverse Mendelian randomisation, and the results obtained by MR-Egger regression and the weighted median method were used as references. Subsequently, the metabolites causally associated with erectile dysfunction were subjected to replication analyses and meta-analyses, and the results of the meta-analyses were analysed by pathway analyses to find influential pathways. In this process, Mendelian randomisation results need to be assessed for stability and reliability using sensitivity analysis.It was found that a total of six serum metabolites were causally associated with erectile dysfunction in a forward Mendelian randomisation study. 1,3,7-trimethyluraten (0.85 (0.73-0.99), P = 0.0368), ergothioneine (0.65 (0.45-0.94), P = 0.0226) and gamma-glutamylglutamate (0.63 (0.46-0.88), P = 0.0059) were protective against the development of erectile dysfunction, whereas 2-hydroxyhippurate (1.10 (1.02-1.19), P = 0.0152), N2,N2-dimethylguanosine (1.57 (1.02-2.40), P = 0.0395) and octanoylcarnitine (1.38 (1.06-1.82), P = 0.0183) were able to induce the development of erectile dysfunction. In addition, metabolic pathway analysis showed that 1,3,7-trimethylurate was able to influence the development of erectile dysfunction via the caffeine metabolism pathway (P = 0.0454). On the other hand, reverse Mendelian randomisation analysis showed that erectile dysfunction reduced serum homocitrulline levels (0.99 (0.97-1.00), P = 0.0360). Sensitivity analyses, including heterogeneity tests and pleiotropy tests, confirmed the reliability of the results.In conclusion, this study demonstrated a bidirectional causal relationship between serum metabolites and erectile dysfunction using bidirectional Mendelian randomisation analysis and replication meta-analysis. On this basis, this study provides a new direction of thinking and strong evidence for the therapeutic application and adjunctive diagnosis of serum metabolites in erectile dysfunction, and provides a certain reference value for subsequent related studies.

Building national patient registries in Mexico: insights from the MexOMICS Consortium

Objective

To introduce MexOMICS, a Mexican Consortium focused on establishing electronic databases to collect, cross-reference, and share health-related and omics data on the Mexican population.

Methods

Since 2019, the MexOMICS Consortium has established three electronic-based registries: the Mexican Twin Registry (TwinsMX), Mexican Lupus Registry (LupusRGMX), and the Mexican Parkinson's Research Network (MEX-PD), designed and implemented using the Research Electronic Data Capture web-based application. Participants were enrolled through voluntary participation and on-site engagement with medical specialists. We also acquired DNA samples and Magnetic Resonance Imaging scans in subsets of participants.

Results

The registries have successfully enrolled a large number of participants from a variety of regions within Mexico: TwinsMX (n = 2,915), LupusRGMX (n = 1,761) and MEX-PD (n = 750). In addition to sociodemographic, psychosocial, and clinical data, MexOMICS has collected DNA samples to study the genetic biomarkers across the three registries. Cognitive function has been assessed with the Montreal Cognitive Assessment in a subset of 376 MEX-PD participants. Furthermore, a subset of 267 twins have participated in cognitive evaluations with the Creyos platform and in MRI sessions acquiring structural, functional, and spectroscopy brain imaging; comparable evaluations are planned for LupusRGMX and MEX-PD.

Conclusions

The MexOMICS registries offer a valuable repository of information concerning the potential interplay of genetic and environmental factors in health conditions among the Mexican population.

MC1R and age heteroclassification of face phenotypes in the Rio Grande do Sul population

BACKGROUND

Forensic DNA phenotyping (FDP) consists of the use of methodologies for predicting externally visible characteristics (EVCs) from the genetic material of biological samples found in crime scenes and has proven to be a promising tool in aiding human identification in police activities. Currently, methods based on multiplex assays and statistical models of prediction of EVCs related to hair, skin, and iris pigmentation using panels of SNP and INDEL biomarkers have already been developed and validated by the forensic scientific community. As well as traces of pigmentation, an individual's perceived age (PA) can also be considered an EVC and its estimation in unknown individuals can be useful for the progress of investigations. Liu and colleagues (2016) were pioneers in evidencing that, in addition to lifestyle and environmental factors, the presence of SNP and INDEL variants in the MC1R gene - which encodes a transmembrane receptor responsible for regulating melanin production - seems to contribute to an individual's PA. The group highlighted the association between these MC1R gene polymorphisms and the PA in the European population, where carriers of risk haplotypes appeared to be up to 2 years older in comparison to their chronological age (CA).

PURPOSE

Understanding that genotype-phenotype relationships cannot be extrapolated between different population groups, this study aimed to test this hypothesis and verify the applicability of this variant panel in the Rio Grande do Sul admixed population.

METHODS

Based on genomic data from a sample of 261 volunteers representative of gaucho population and using a multiple linear regression (MLR) model, our group was able to verify a significant association among nine intronic variants in loci adjacent to MC1R (e.g., AFG3L1P, TUBB3, FANCA) and facial age appearance, whose PA was defined after age heteroclassification of standard frontal face images through 11 assessors.

RESULTS

Different from that observed in European populations, our results show that the presence of effect alleles (R) of the selected variants in our sample influenced both younger and older face phenotypes. The influence of each variant on PA is expressed as β values.

CONCLUSIONS

There are important molecular mechanisms behind the effects of MC1R locus on PA, and the genomic background of each population seems to be crucial to determine this influence.

Identification of blood metabolites associated with risk of Alzheimer's disease by integrating genomics and metabolomics data

Specific metabolites have been reported to be potentially associated with Alzheimer's disease (AD) risk. However, the comprehensive understanding of roles of metabolite biomarkers in AD etiology remains elusive. We performed a large AD metabolome-wide association study (MWAS) by developing blood metabolite genetic prediction models. We evaluated associations between genetically predicted levels of metabolites and AD risk in 39,106 clinically diagnosed AD cases, 46,828 proxy AD and related dementia (proxy-ADD) cases, and 401,577 controls. We further conducted analyses to determine microbiome features associated with the detected metabolites and characterize associations between predicted microbiome feature levels and AD risk. We identified fourteen metabolites showing an association with AD risk. Five microbiome features were further identified to be potentially related to associations of five of the metabolites. Our study provides new insights into the etiology of AD that involves blood metabolites and gut microbiome, which warrants further investigation.

Metformin and feeding increase levels of the appetite-suppressing metabolite Lac-Phe in humans

Metformin, a widely used first-line treatment for type 2 diabetes (T2D), is known to reduce blood glucose levels and suppress appetite. Here we report a significant elevation of the appetite-suppressing metabolite N-lactoyl phenylalanine (Lac-Phe) in the blood of individuals treated with metformin across seven observational and interventional studies. Furthermore, Lac-Phe levels were found to rise in response to acute metformin administration and post-prandially in patients with T2D or in metabolically healthy volunteers.

Genetic influences on circulating retinol and its relationship to human health

Retinol is a fat-soluble vitamin that plays an essential role in many biological processes throughout the human lifespan. Here, we perform the largest genome-wide association study (GWAS) of retinol to date in up to 22,274 participants. We identify eight common variant loci associated with retinol, as well as a rare-variant signal. An integrative gene prioritisation pipeline supports novel retinol-associated genes outside of the main retinol transport complex (RBP4:TTR) related to lipid biology, energy homoeostasis, and endocrine signalling. Genetic proxies of circulating retinol were then used to estimate causal relationships with almost 20,000 clinical phenotypes via a phenome-wide Mendelian randomisation study (MR-pheWAS). The MR-pheWAS suggests that retinol may exert causal effects on inflammation, adiposity, ocular measures, the microbiome, and MRI-derived brain phenotypes, amongst several others. Conversely, circulating retinol may be causally influenced by factors including lipids and serum creatinine. Finally, we demonstrate how a retinol polygenic score could identify individuals more likely to fall outside of the normative range of circulating retinol for a given age. In summary, this study provides a comprehensive evaluation of the genetics of circulating retinol, as well as revealing traits which should be prioritised for further investigation with respect to retinol related therapies or nutritional intervention.

Comparative analysis of mutational patterns in triple negative breast cancer before and after neoadjuvant chemotherapy in patients with residual disease

Triple-negative breast cancer (TNBC) is a heterogeneous disease with poor survival compared to other subtypes. Patients with residual disease after neoadjuvant chemotherapy (NAC) face an increased risk of relapse and death. We aimed to characterize the mutational landscape of this subset to offer insights into relapse pathogenesis and potential therapeutic targets. We retrospectively analyzed archived paired (pre- and post-NAC) tumor samples from 25 patients with TNBC with residual disease using a targeted 72-gene next-generation sequencing panel. Our findings revealed a stable mutational burden in both pre- and post-NAC samples, with a median count of 12 variants (IQR 7-17.25) per sample. TP53, PMS2, PTEN, ERBB2, and NOTCH1 variants were observed in pre-NAC samples predominantly. Notably, post-NAC samples exhibited a significant increase in AR gene mutations, suggesting potential prognostic and predictive implications. No difference in mutational burden was found between patients who did and did not receive platinum (p = 0.94), or between those with and without recurrence (p = 0.49). We employed K-means clustering to categorize the patients based on their variant profiles, aiding in the prediction of possible patterns associated with recurrence. Our study was limited by its small sample size and retrospective design, suggesting the need for further validation in larger prospective cohorts.

Investigating causal associations among gut microbiota, metabolites, and psoriatic arthritis: a Mendelian randomization study

Background

Currently, there has been observed a significant alteration in the composition of the gut microbiome (GM) and serum metabolites in patients with psoriatic arthritis (PsA) compared to healthy individuals. However, previous observational studies have shown inconsistent results regarding the alteration of gut microbiota/metabolites. In order to shed light on this matter, we utilized Mendelian randomization to determine the causal effect of GM/metabolites on PsA.

Methods

We retrieved summary-level data of GM taxa/metabolites and PsA from publicly available GWAS statistics. Causal relationships between GM/metabolites and PsA were determined using a two-sample MR analysis, with the IVW approach serving as the primary analysis method. To ensure the robustness of our findings, we conducted sensitivity analyses, multivariable MR analysis (MVMR), and additional analysis including replication verification analysis, LDSC regression, and Steiger test analysis. Furthermore, we investigated reverse causality through a reverse MR analysis. Finally, we conducted an analysis of expression quantitative trait loci (eQTLs) involved in the metabolic pathway to explore potential molecular mechanisms of metabolism.

Results

Our findings reveal that eight GM taxa and twenty-three serum metabolites are causally related to PsA (P < 0.05). Notably, a higher relative abundance of Family Rikenellaceae (ORIVW: 0.622, 95% CI: 0.438-0.883, FDR = 0.045) and elevated serum levels of X-11538 (ORIVW: 0.442, 95% CI: 0.250-0.781, FDR = 0.046) maintain significant causal associations with a reduced risk of PsA, even after adjusting for multiple testing correction and conducting MVMR analysis. These findings suggest that Family Rikenellaceae and X-11538 may have protective effects against PsA. Our sensitivity analysis and additional analysis revealed no significant horizontal pleiotropy, reverse causality, or heterogeneity. The functional enrichment analysis revealed that the eQTLs examined were primarily associated with glycerolipid metabolism and the expression of key metabolic factors influenced by bacterial infections (Vibrio cholerae and Helicobacter pylori) as well as the mTOR signaling pathway.

Conclusion

In conclusion, our study demonstrates that Family Rikenellaceae and X-11538 exhibit a strong and negative causal relationship with PsA. These particular GM taxa and metabolites have the potential to serve as innovative biomarkers, offering valuable insights into the treatment and prevention of PsA. Moreover, bacterial infections and mTOR-mediated activation of metabolic factors may play an important role in this process.

Glaucoma Patients Have a Lower Abundance of Butyrate-Producing Taxa in the Gut

Purpose

Glaucoma is an eye disease that is the most common cause of irreversible blindness worldwide. It has been suggested that gut microbiota can produce reactive oxygen species and pro-inflammatory cytokines that may travel from the gastric mucosa to distal sites, for example, the optic nerve head or trabecular meshwork. There is evidence for a gut-eye axis, as microbial dysbiosis has been associated with retinal diseases. We investigated the microbial composition in patients with glaucoma and healthy controls. Moreover, we analyzed the association of the gut microbiome with intraocular pressure (IOP; risk factor of glaucoma) and vertical cup-to-disc ratio (VCDR; quantifying glaucoma severity).

Methods

The discovery analyses included participants of the Rotterdam Study and the Erasmus Glaucoma Cohort. A total of 225 patients with glaucoma and 1247 age- and sex-matched participants without glaucoma were included in our analyses. Stool samples were used to generate 16S rRNA gene profiles. We assessed associations with 233 genera and species. We used data from the TwinsUK and the Study of Health in Pomerania (SHIP) to replicate our findings.

Results

Several butyrate-producing taxa (e.g. Butyrivibrio, Caproiciproducens, Clostridium sensu stricto 1, Coprococcus 1, Ruminococcaceae UCG 007, and Shuttleworthia) were less abundant in people with glaucoma compared to healthy controls. The same taxa were also associated with lower IOP and smaller VCDR. The replication analyses confirmed the findings from the discovery analyses.

Conclusions

Large human studies exploring the link between the gut microbiome and glaucoma are lacking. Our results suggest that microbial dysbiosis plays a role in the pathophysiology of glaucoma.

Trans-ancestry epigenome-wide association meta-analysis of DNA methylation with lifetime cannabis use

Cannabis is widely used worldwide, yet its links to health outcomes are not fully understood. DNA methylation can serve as a mediator to link environmental exposures to health outcomes. We conducted an epigenome-wide association study (EWAS) of peripheral blood-based DNA methylation and lifetime cannabis use (ever vs. never) in a meta-analysis including 9436 participants (7795 European and 1641 African ancestry) from seven cohorts. Accounting for effects of cigarette smoking, our trans-ancestry EWAS meta-analysis revealed four CpG sites significantly associated with lifetime cannabis use at a false discovery rate of 0.05 ( p < 5.85 × 10 - 7 ) : cg22572071 near gene ADGRF1, cg15280358 in ADAM12, cg00813162 in ACTN1, and cg01101459 near LINC01132. Additionally, our EWAS analysis in participants who never smoked cigarettes identified another epigenome-wide significant CpG site, cg14237301 annotated to APOBR. We used a leave-one-out approach to evaluate methylation scores constructed as a weighted sum of the significant CpGs. The best model can explain 3.79% of the variance in lifetime cannabis use. These findings unravel the DNA methylation changes associated with lifetime cannabis use that are independent of cigarette smoking and may serve as a starting point for further research on the mechanisms through which cannabis exposure impacts health outcomes.

Exploring quantitative traits-associated copy number deletions through reanalysis of UK10K consortium whole genome sequencing cohorts

OBJECTIVES

We performed comprehensive association analyses of common high-confidence gnomAD-reported copy number deletions (CNDs) with 60 quantitative traits from UK10K consortium WGS data.

METHODS

The study made use of data generated by the UK10K Consortium. UK10K consortium WGS data consist of TwinsUK (n = 1754, middle-aged females) and ALSPAC (n = 1867, birth to adolescence) cohorts. UK10K consortium called 18,739 CNDs (hg19) with GenomeSTRiP software. After filtering out variants with minor allele frequency < 0.05 or HWE P < 1.0 × 10- 6, 1222 (TwinsUK) and 1211 (ALSPAC) CNDs remained for association analyses with 60 normalized quantitative traits.

RESULTS

We identified 23 genome-wide significant associations at 13 loci, among which 2 associations reached experiment-wide significance. We found that two common deletions in chromosome 4, located between WDR1 and ZNF518B (23.3 kb, dbVar ID:nssv15888957, 4:10211262-10,234,569 and 9.8 kb, dbVar ID:nssv15888975, 4:10392422-10,402,191), were associated with uric acid levels (P = 5.23 × 10- 11 and 2.29 × 10- 8, respectively). We also discovered a novel deletion spanning chromosome 18 (823 bp, dbVar ID: nssv15841628, 8:74347187-74,348,010) associated with low HDL cholesterol levels (P = 4.15 × 10- 7). Additionally, we observed two red blood cell traits-associated loci with genome-wide significance, a 13.2 kb deletion in 7q22.1 (nssv15922542) and a 3.7 kb deletion in 12q24.12 (nssv15813226), both of which were located in regions previously reported to be associated with red blood cell traits. Two deletions in 11q11 (nssv15803200 and nssv15802240), where clusters of multiple olfactory receptor genes exist, and a deletion (nssv15929560) upstream to DOCK5 were associated with childhood obesity. Finally, when defining Trait-Associated copy number Deletions (TADs) as CNDs with phenotype associations at sub-threshold significance (P < 10- 3), we identified 157 (97.5%) out of 161 TADs in non-coding regions, with a mean size of 4 kb (range: 209 - 47,942 bp).

CONCLUSION

We conducted a reanalysis of the UK10K Whole Genome Sequencing cohort, which led to the identification of multiple high confidence copy number deletions associated with quantitative traits. These deletions have standard dbVar IDs and replicate previous findings, as well as reveal novel loci that require further replication studies.

Gene-metabolite annotation with shortest reactional distance enhances metabolite genome-wide association studies results

Metabolite genome-wide association studies (mGWAS) have advanced our understanding of the genetic control of metabolite levels. However, interpreting these associations remains challenging due to a lack of tools to annotate gene-metabolite pairs beyond the use of conservative statistical significance threshold. Here, we introduce the shortest reactional distance (SRD) metric, drawing from the comprehensive KEGG database, to enhance the biological interpretation of mGWAS results. We applied this approach to three independent mGWAS, including a case study on sickle cell disease patients. Our analysis reveals an enrichment of small SRD values in reported mGWAS pairs, with SRD values significantly correlating with mGWAS p values, even beyond the standard conservative thresholds. We demonstrate the utility of SRD annotation in identifying potential false negatives and inaccuracies within current metabolic pathway databases. Our findings highlight the SRD metric as an objective, quantitative and easy-to-compute annotation for gene-metabolite pairs, suitable to integrate statistical evidence to biological networks.

Protective effect of uridine on atrial fibrillation: a Mendelian randomisation study

Uridine, a pyrimidine nucleoside, is crucial in the synthesis of metabolites. According to observational studies, a higher plasma uridine level is associated with a lower risk of atrial fibrillation (AF). However, the casual relationship between uridine and AF is still unknown. In this study, we used the Mendelian randomisation (MR) approach to explore causality. Three genetic variants associated with uridine were identified from the Metabolomics GWAS server (7824 participants); summary-level datasets associated with AF were acquired from a genome-wide association study (GWAS) meta-analysis with 1,030,836 European participants (60,620 AF cases). We duplicated the MR analyses using datasets from AF HRC studies and the FinnGen Consortium, and then conducted a meta-analysis which combined the main results. The risk of AF was significantly associated with the genetically determined plasma uridine level (odds ratio [OR] 0.27; 95% confidence interval [CI] 0.16, 0.47; p = 2.39 × 10-6). The association remained consistent in the meta-analysis of the various datasets (OR 0.27; 95% CI 0.17, 0.42; p = 1.34 × 10-8). In conclusion, the plasma uridine level is inversely associated with the risk of AF. Raising the plasma uridine level may have prophylactic potential against AF.

Speech processing performance of Hungarian-speaking twins and singletons

Studying speech processing in twins versus their singleton peers provides opportunities to study both genetic and environmental effects on how children acquire these aspects of their speech and - by extension - their phonological systems. Our study focused on speech processing in typically developing Hungarian-speaking twins and their singleton peers between 5 and 9 years of age. Participants included 384 monolingual Hungarian-speaking children (192 twins, and 192 singletons). Data from four tasks - repetition of synthesised monosyllables, nonsense words, well-formed noisy sentences, and well-formed phonologically complex sentences - were analysed. There was a main effect for birth status, and singletons outperformed their twin peers on the majority of the speech processing tasks. Age and task also had effects on the performance of the participants, and there was a three-way task by age by twin versus singleton status indicating that the speech processing performance of twins versus singletons is interdependent with the type of task and age. Our results also indicate that monolingual Hungarian-speaking twins may be at higher risk for developmental speech delays relative to their singleton peers.

Reproducible microbiome composition signatures of anxiety and depressive symptoms

The gut microbiome is a significant contributor to mental health, with growing evidence linking its composition to anxiety and depressive disorders. Gut microbiome composition is associated with signs of anxiety and depression both in clinically diagnosed mood disorders and subclinically in the general population and may be influenced by dietary fibre intake and the presence of chronic pain. We provide an update of current evidence on the role of gut microbiome composition in depressive and anxiety disorders or symptoms by reviewing available studies. Analysing data from three independent cohorts (osteoarthritis 1 (OA1); n = 46, osteoarthritis 2 (OA2); n = 58, and healthy controls (CON); n = 67), we identified microbial composition signatures of anxiety and depressive symptoms at genus level and cross-validated our findings performing meta-analyses of our results with results from previously published studies. The genera Bifidobacterium (fixed-effect beta (95% CI) = -0.22 (-0.34, -0.10), p = 3.90e-04) and Lachnospiraceae NK4A136 group (fixed-effect beta (95% CI) = -0.09 (-0.13, -0.05), p = 2.53e-06) were found to be the best predictors of anxiety and depressive symptoms, respectively, across our three cohorts and published literature taking into account demographic and lifestyle covariates, such as fibre intake. The association with anxiety was robust in accounting for heterogeneity between cohorts and supports previous observations of the potential prophylactic effect of Bifidobacterium against anxiety symptoms.

Characterisation of Fasting and Postprandial NMR Metabolites: Insights from the ZOE PREDICT 1 Study

BACKGROUND

Postprandial metabolomic profiles and their inter-individual variability are not well characterised. Here, we describe postprandial metabolite changes, their correlations with fasting values and their inter- and intra-individual variability, following a standardised meal in the ZOE PREDICT 1 cohort.

METHODS

In the ZOE PREDICT 1 study (n = 1002 (NCT03479866)), 250 metabolites, mainly lipids, were measured by a Nightingale NMR panel in fasting and postprandial (4 and 6 h after a 3.7 MJ mixed nutrient meal, with a second 2.2 MJ mixed nutrient meal at 4 h) serum samples. For each metabolite, inter- and intra-individual variability over time was evaluated using linear mixed modelling and intraclass correlation coefficients (ICC) were calculated.

RESULTS

Postprandially, 85% (of 250 metabolites) significantly changed from fasting at 6 h (47% increased, 53% decreased; Kruskal-Wallis), with 37 measures increasing by >25% and 14 increasing by >50%. The largest changes were observed in very large lipoprotein particles and ketone bodies. Seventy-one percent of circulating metabolites were strongly correlated (Spearman's rho >0.80) between fasting and postprandial timepoints, and 5% were weakly correlated (rho <0.50). The median ICC of the 250 metabolites was 0.91 (range 0.08-0.99). The lowest ICCs (ICC <0.40, 4% of measures) were found for glucose, pyruvate, ketone bodies (β-hydroxybutyrate, acetoacetate, acetate) and lactate.

CONCLUSIONS

In this large-scale postprandial metabolomic study, circulating metabolites were highly variable between individuals following sequential mixed meals. Findings suggest that a meal challenge may yield postprandial responses divergent from fasting measures, specifically for glycolysis, essential amino acid, ketone body and lipoprotein size metabolites.

Genetic control of N-glycosylation of human blood plasma proteins

Glycosylation is an important protein modification, which influences the physical and chemical properties as well as biological function of these proteins. Large-scale population studies have shown that the levels of various plasma protein N-glycans are associated with many multifactorial human diseases. Observed associations between protein glycosylation levels and human diseases have led to the conclusion that N-glycans can be considered a potential source of biomarkers and therapeutic targets. Although biochemical pathways of glycosylation are well studied, the understanding of the mechanisms underlying general and tissue-specific regulation of these biochemical reactions in vivo is limited. This complicates both the interpretation of the observed associations between protein glycosylation levels and human diseases, and the development of glycan-based biomarkers and therapeutics. By the beginning of the 2010s, high-throughput methods of N-glycome profiling had become available, allowing research into the genetic control of N-glycosylation using quantitative genetics methods, including genome-wide association studies (GWAS). Application of these methods has made it possible to find previously unknown regulators of N-glycosylation and expanded the understanding of the role of N-glycans in the control of multifactorial diseases and human complex traits. The present review considers the current knowledge of the genetic control of variability in the levels of N-glycosylation of plasma proteins in human populations. It briefly describes the most popular physical-chemical methods of N-glycome profiling and the databases that contain genes involved in the biosynthesis of N-glycans. It also reviews the results of studies of environmental and genetic factors contributing to the variability of N-glycans as well as the mapping results of the genomic loci of N-glycans by GWAS. The results of functional in vitro and in silico studies are described. The review summarizes the current progress in human glycogenomics and suggests possible directions for further research.

Twins Early Development Study (TEDS): A genetically sensitive investigation of mental health outcomes in the mid-twenties

The Twins Early Development Study (TEDS) is a longitudinal study following a cohort of twins born 1994-1996 in England and Wales. Of the 13,759 families who originally consented to take part, over 10,000 families remain enrolled in the study. The current focus of TEDS is on mental health in the mid-twenties. Making use of over 25 years of genetically sensitive data, TEDS is uniquely placed to explore the longitudinal genetic and environmental influences on common mental health disorders in early adulthood. This paper outlines recent data collection efforts supporting this work, including a cohort-wide mental health assessment at age 26 and a multi-phase Covid-19 study. It will also provide an update on data linkage efforts and the Children of TEDS (CoTEDS) project.

Spectrum and frequency of genetic variants in sporadic amyotrophic lateral sclerosis

Therapy of motoneuron diseases entered a new phase with the use of intrathecal antisense oligonucleotide therapies treating patients with specific gene mutations predominantly in the context of familial amyotrophic lateral sclerosis. With the majority of cases being sporadic, we conducted a cohort study to describe the mutational landscape of sporadic amyotrophic lateral sclerosis. We analysed genetic variants in amyotrophic lateral sclerosis-associated genes to assess and potentially increase the number of patients eligible for gene-specific therapies. We screened 2340 sporadic amyotrophic lateral sclerosis patients from the German Network for motor neuron diseases for variants in 36 amyotrophic lateral sclerosis-associated genes using targeted next-generation sequencing and for the C9orf72 hexanucleotide repeat expansion. The genetic analysis could be completed on 2267 patients. Clinical data included age at onset, disease progression rate and survival. In this study, we found 79 likely pathogenic Class 4 variants and 10 pathogenic Class 5 variants (without the C9orf72 hexanucleotide repeat expansion) according to the American College of Medical Genetics and Genomics guidelines, of which 31 variants are novel. Thus, including C9orf72 hexanucleotide repeat expansion, Class 4, and Class 5 variants, 296 patients, corresponding to ∼13% of our cohort, could be genetically resolved. We detected 437 variants of unknown significance of which 103 are novel. Corroborating the theory of oligogenic causation in amyotrophic lateral sclerosis, we found a co-occurrence of pathogenic variants in 10 patients (0.4%) with 7 being C9orf72 hexanucleotide repeat expansion carriers. In a gene-wise survival analysis, we found a higher hazard ratio of 1.47 (95% confidence interval 1.02-2.1) for death from any cause for patients with the C9orf72 hexanucleotide repeat expansion and a lower hazard ratio of 0.33 (95% confidence interval 0.12-0.9) for patients with pathogenic SOD1 variants than for patients without a causal gene mutation. In summary, the high yield of 296 patients (∼13%) harbouring a pathogenic variant and oncoming gene-specific therapies for SOD1/FUS/C9orf72, which would apply to 227 patients (∼10%) in this cohort, corroborates that genetic testing should be made available to all sporadic amyotrophic lateral sclerosis patients after respective counselling.

The targeted next-generation sequence revealed SMAD4, AKT1, and TP53 mutations from circulating cell-free DNA of breast cancer and its effect on protein structure - A computational approach

Breast cancer biomarkers that detect marginally advanced stages are still challenging. The detection of specific abnormalities, targeted therapy selection, prognosis, and monitoring of treatment effectiveness over time are all made possible by circulating free DNA (cfDNA) analysis. The proposed study will detect specific genetic abnormalities from the plasma cfDNA of a female breast cancer patient by sequencing a cancer-related gene panel (MGM455 - Oncotrack Ultima), including 56 theranostic genes (SNVs and small INDELs). Initially, we determined the pathogenicity of the observed mutations using PredictSNP, iStable, Align-GVGD, and ConSurf servers. As a next step, molecular dynamics (MD) was implemented to determine the functional significance of SMAD4 mutation (V465M). Lastly, the mutant gene relationships were examined using the Cytoscape plug-in GeneMANIA. Using ClueGO, we determined the gene's functional enrichment and integrative analysis. The structural characteristics of SMAD4 V465M protein by MD simulation analysis further demonstrated that the mutation was deleterious. The simulation showed that the native structure was more significantly altered by the SMAD4 (V465M) mutation. Our findings suggest that SMAD4 V465M mutation might be significantly associated with breast cancer, and other patient-found mutations (AKT1-E17K and TP53-R175H) are synergistically involved in the process of SMAD4 translocate to nuclease, which affects the target gene translation. Therefore, this combination of gene mutations could alter the TGF-β signaling pathway in BC. We further proposed that the SMAD4 protein loss may contribute to an aggressive phenotype by inhibiting the TGF-β signaling pathway. Thus, breast cancer's SMAD4 (V465M) mutation might increase their invasive and metastatic capabilities.Communicated by Ramaswamy H. Sarma.

Gene-metabolite annotation with shortest reactional distance enhances metabolite genome-wide association studies results

Studies combining metabolomics and genetics, known as metabolite genome-wide association studies (mGWAS), have provided valuable insights into our understanding of the genetic control of metabolite levels. However, the biological interpretation of these associations remains challenging due to a lack of existing tools to annotate mGWAS gene-metabolite pairs beyond the use of conservative statistical significance threshold. Here, we computed the shortest reactional distance (SRD) based on the curated knowledge of the KEGG database to explore its utility in enhancing the biological interpretation of results from three independent mGWAS, including a case study on sickle cell disease patients. Results show that, in reported mGWAS pairs, there is an excess of small SRD values and that SRD values and p-values significantly correlate, even beyond the standard conservative thresholds. The added-value of SRD annotation is shown for identification of potential false negative hits, exemplified by the finding of gene-metabolite associations with SRD ≤1 that did not reach standard genome-wide significance cut-off. The wider use of this statistic as an mGWAS annotation would prevent the exclusion of biologically relevant associations and can also identify errors or gaps in current metabolic pathway databases. Our findings highlight the SRD metric as an objective, quantitative and easy-to-compute annotation for gene-metabolite pairs that can be used to integrate statistical evidence to biological networks.

Escape from X-inactivation in twins exhibits intra- and inter-individual variability across tissues and is heritable

X-chromosome inactivation (XCI) silences one X in female cells to balance sex-differences in X-dosage. A subset of X-linked genes escape XCI, but the extent to which this phenomenon occurs and how it varies across tissues and in a population is as yet unclear. To characterize incidence and variability of escape across individuals and tissues, we conducted a transcriptomic study of escape in adipose, skin, lymphoblastoid cell lines and immune cells in 248 healthy individuals exhibiting skewed XCI. We quantify XCI escape from a linear model of genes' allelic fold-change and XIST-based degree of XCI skewing. We identify 62 genes, including 19 lncRNAs, with previously unknown patterns of escape. We find a range of tissue-specificity, with 11% of genes escaping XCI constitutively across tissues and 23% demonstrating tissue-restricted escape, including cell type-specific escape across immune cells of the same individual. We also detect substantial inter-individual variability in escape. Monozygotic twins share more similar escape than dizygotic twins, indicating that genetic factors may underlie inter-individual differences in escape. However, discordant escape also occurs within monozygotic co-twins, suggesting environmental factors also influence escape. Altogether, these data indicate that XCI escape is an under-appreciated source of transcriptional differences, and an intricate phenotype impacting variable trait expressivity in females.

Associations between genetically predicted levels of blood metabolites and pancreatic cancer risk

Pancreatic ductal adenocarcinoma (PDAC) is one of the most aggressive solid malignancies, which is featured by systematic metabolism. Thus, a better understanding of metabolic dysregulation in PDAC is important to better characterize its etiology. Here, we performed a large metabolome-wide association study (MWAS) to systematically explore associations between genetically predicted metabolite levels in blood and PDAC risk. Using data from 881 subjects of European descent in the TwinsUK Project, comprehensive genetic models were built to predict serum metabolite levels. These prediction models were applied to the genetic data of 8275 cases and 6723 controls included in the PanScan (I, II and III) and PanC4 consortia. After assessing the metabolite-PDAC risk associations by a slightly modified TWAS/FUSION framework, we identified five metabolites (including two dipeptides) showing significant associations with PDAC risk at false discovery rate (FDR) <0.05. Integrated with gut microbial information, two-sample Mendelian randomization (MR) analyses were further performed to investigate the relationship among serum metabolites, gut microbiome features and PDAC. The flavonoid-degrading bacteria Flavonifractor sp90199495 was found to be associated with metabolite X-21849 and it was also shown to be associated with PDAC risk. Collectively, our study identified novel candidate metabolites for PDAC risk, which could lead to new insights into the etiology of PDAC and improved treatment options.

An integrated investigation of structural and pathway alteration caused by PIK3CA and TP53 mutations identified in cfDNA of metastatic breast cancer

In peripheral blood, cell-free DNA (cfDNA) contains circulating tumor DNA (ctDNA), which indicates molecular abnormalities in metastatic breast tumor tissue. The sequencing of cfDNA of Metastatic Breast Cancer (MBC) patients allows assessment of therapy response and noninvasive treatment. In the proposed study, clinically significant alterations in PIK3CA and TP53 genes associated with MBC resulting in a missense substitution of His1047Arg and Arg282Trp from an next-generation sequencing-based multi-gene panel were reported in a cfDNA of a patient with MBC. To investigate the impact of the reported mutation, we used molecular docking, molecular dynamics simulation, network analysis, and pathway analysis. Molecular Docking analysis determined the distinct binding pattern revealing H1047R-ATP complex has a higher number of Hydrogen bonds (H-bonds) and binding affinity with a slight difference compared to the PIK3CA-ATP complex. Following, molecular dynamics simulation for 200 ns, of which H1047R-ATP complex resulted in the instability of PIK3CA. Similarly, for TP53 mutant R282W, the zinc-free state (apo) and zinc-bounded (holo) complexes were investigated for conformational change between apo and holo complexes, of which the holo complex mutant R282W was unstable. To validate the conformational change of PIK3CA and TP53, 80% mutation of H1047R in the kinase domain of p110α expressed ubiquitously in PIK3CA protein that alters PI3K pathway, while R282W mutation in DNA binding helix (H2) region of P53 protein inhibits the transcription factor in P53 pathway causing MBC. According to our findings, the extrinsic (hypoxia, oxidative stress, and acidosis); intrinsic factors (MYC amplification) in PIK3CA and TP53 mutations will provide potential insights for developing novel therapeutic methods for MBC therapy.

Voluntary risk mitigation behaviour can reduce impact of SARS-CoV-2: a real-time modelling study of the January 2022 Omicron wave in England

BACKGROUND

Predicting the likely size of future SARS-CoV-2 waves is necessary for public health planning. In England, voluntary "plan B" mitigation measures were introduced in December 2021 including increased home working and face coverings in shops but stopped short of restrictions on social contacts. The impact of voluntary risk mitigation behaviours on future SARS-CoV-2 burden is unknown.

METHODS

We developed a rapid online survey of risk mitigation behaviours ahead of the winter 2021 festive period and deployed in two longitudinal cohort studies in the UK (Avon Longitudinal Study of Parents and Children (ALSPAC) and TwinsUK/COVID Symptom Study (CSS) Biobank) in December 2021. Using an individual-based, probabilistic model of COVID-19 transmission between social contacts with SARS-CoV-2 Omicron variant parameters and realistic vaccine coverage in England, we predicted the potential impact of the SARS-CoV-2 Omicron wave in England in terms of the effective reproduction number and cumulative infections, hospital admissions and deaths. Using survey results, we estimated in real-time the impact of voluntary risk mitigation behaviours on the Omicron wave in England, if implemented for the entire epidemic wave.

RESULTS

Over 95% of survey respondents (NALSPAC = 2686 and NTwins = 6155) reported some risk mitigation behaviours, with vaccination and using home testing kits reported most frequently. Less than half of those respondents reported that their behaviour was due to "plan B". We estimate that without risk mitigation behaviours, the Omicron variant is consistent with an effective reproduction number between 2.5 and 3.5. Due to the reduced vaccine effectiveness against infection with the Omicron variant, our modelled estimates suggest that between 55% and 60% of the English population could be infected during the current wave, translating into between 12,000 and 46,000 cumulative deaths, depending on assumptions about severity and vaccine effectiveness. The actual number of deaths was 15,208 (26 November 2021-1 March 2022). We estimate that voluntary risk reduction measures could reduce the effective reproduction number to between 1.8 and 2.2 and reduce the cumulative number of deaths by up to 24%.

CONCLUSIONS

Predicting future infection burden is affected by uncertainty in disease severity and vaccine effectiveness estimates. In addition to biological uncertainty, we show that voluntary measures substantially reduce the projected impact of the SARS-CoV-2 Omicron variant but that voluntary measures alone would be unlikely to completely control transmission.

Challenges and opportunities for omics-based precision medicine in chronic low back pain

PURPOSE

Chronic low back pain (cLBP) is a common health condition worldwide and a leading cause of disability with an estimated lifetime prevalence of 80-90% in industrialized countries. However, we have had limited success in treating cLBP likely due to its non-specific heterogeneous nature that goes beyond detectable anatomical changes. We propose that omics technologies as precision medicine tools are well suited to provide insight into its pathophysiology and provide diagnostic markers and therapeutic targets. Therefore, in this review, we explore the current state of omics technologies in the diagnosis and classification of cLBP. We identify factors that may serve as markers to differentiate between acute and chronic cases of low back pain (LBP). Finally, we also discuss some challenges that must be overcome to successfully apply precision medicine to the diagnosis and treatment of cLBP.

METHODS

A literature search for the current applications of omics technologies to chronic low back pain was performed using the following search terms- "back pain," "low back pain," "proteomics," "transcriptomics", "epigenomics," "genomics," "omics." We reviewed molecular markers identified from 35 studies which hold promise in providing information regarding molecular insights into cLBP.

RESULTS

GWAS studies have found evidence for the role of single nucleotide polymorphisms (SNPs) associated with pain pathways in individuals with cLBP. Epigenomic modifications in patients with cLBP have been found to be enriched among genes involved in immune signaling and inflammation. Transcriptomics profiles of patients with cLBP show multiple lines of evidence for the role of inflammation in cLBP. The glycomics profiles of patients with cLBP are similar to those of patients with inflammatory conditions. Proteomics and microbiomics show promise but have limited studies currently.

CONCLUSION

Omics technologies have identified associations between inflammatory and pain pathways in the pathophysiology of cLBP. However, in order to integrate information across the range of studies, it is important for the field to identify and adopt standardized definitions of cLBP and control patients. Additionally, most papers have applied a single omics method to a sampling of cLBP patients which have yielded limited insight into the pathophysiology of cLBP. Therefore, we recommend a multi-omics approach applied to large global consortia for advancing subphenotyping and better management of cLBP, via improved identification of diagnostic markers and therapeutic targets.

Empowerment Mitigates Gender Differences in Tertius Iungens Brokering

Tertius iungens brokering that brings together people who might not otherwise meet is crucial for organizational effectiveness. But we know little about whether and why women and men differ in their propensity to engage in this brokering. Our paper focuses on the origins and mitigation of gender differences in the propensity to bring people together. In study 1, we showed that the Totterdell et al. [Totterdell P, Holman D, Hukin A (2008) Social networkers: Measuring and examining individual differences in propensity to connect with others. Soc. Networks 30(4):283–296] propensity-to-join-others scale that we used in study 2 and the Obstfeld [Obstfeld D (2005) Social networks, the tertius iungens orientation, and involvement in innovation. Admin. Sci. Quart. 50(1):100–130] tertius iungens scale overlapped not only conceptually, but also empirically, and that these measures of tertius iungens were distinct from mediation- and separation-brokering propensities [Grosser TJ, Obstfeld D, Labianca G, Borgatti SP (2019) Measuring mediation and separation brokerage orientations: A further step toward studying the social network brokerage process. Acad. Management Discoveries 5(2):114–136]. In study 2, we used a natural experiment to examine the tertius iungens brokering propensities of 876 identical and 625 fraternal same-sex twins. We found that brokering propensity was lower for women than for men, although the propensity toward sociability in terms of making friends and acquaintances was lower for men. We also found that for women, relative to men, tertius iungens brokering propensity was largely affected by environmental influences, such as the experience of stereotyping and discrimination, rather than representing an inherited disposition. Moreover, the differences between men and women with respect to brokering were mitigated for empowered samples, such as well-educated or entrepreneurial individuals. Our research asks new questions about how environmental pressures and empowerment affect social networking. Gender differences in brokering may be amenable to mitigation through empowering practices that include education and entrepreneurial experience.

Next-Generation Sequencing Highlights of Diffuse Large B-cell Lymphoma in a Tertiary Care Hospital in North India

Introduction: Next-generation sequencing (NGS) elucidates the diffuse large B-cell lymphoma (DLBCL) genetic characteristics by finding recurrent and novel somatic mutations. This observational study attempted to create an NGS panel with a focus on identifying novel somatic mutations which could have potential clinical and therapeutic implications. This panel was created to look for mutations in 133 genes chosen on basis of a literature review and it was used to sequence the tumor DNA of 20 DLBCL patients after a centralized histopathologic review.

Methods: The study included 20 patients having DLBCL. The quality and quantity of tumor cells were accessed by H&E staining and correlated with histopathology and Immunohistochemistry (IHC) status. Patients were grouped as ABC (activated B-cell), PMBL (primary mediastinal large B-cell lymphoma), and other or unclassified subtypes. The lymphoma panel of 133 was designed on targeted sequencing of multiple genes for the coding regions through NGS. The libraries were prepared and sequenced using the Illumina platform. The alignment of obtained sequences was performed using Burrows-Wheeler Aligner and identification of somatic mutations was done using LoFreq (version 2) variant caller. The mutations were annotated using an annotation pipeline (VariMAT). Previously published literature and databases were used for the annotation of clinically relevant mutations. The common variants were filtered for reporting based on the presence in various population databases (1000G, ExAC, EVS, 1000Japanese, dbSNP, UK10K, MedVarDb). A custom read-depth-based algorithm was used to determine CNV (Copy Number Variants) from targeted sequencing experiments. Rare CNVs were detected using a comparison of the test data read-depths with the matched reference dataset. Reportable mutations were prioritized and prepared based on AMP-ASCO-CAP (Association for Molecular Pathology-American Society of Clinical Oncology-College of American Pathologists), WHO guidelines, and also based on annotation metrics from OncoMD (a knowledge base of genomic alterations).

Results: The informativity of the panel was 95 percent. NOTCH 1 was the most frequently mutated gene in 16.1% of patients followed by 12.9% who had ARID1A mutations. MYD88 and TP53 mutations were detected in 9.6% of the patient while 6.4% of patients had CSF3R mutations. NOTCH 1 and TP 53 are the most frequently reported gene in the middle age group (40-60). Mutation in MYD88 is reported in every age group. MYD88 (51%) is the most common mutation in ABC subtypes of DLBCL, followed by NOTCH 1 (44%) and SOCS 1 (33%) according to our findings. NOTCH 1 mutations are frequent in ABC and PMBL subtypes. Closer investigation reveals missense mutation is the most frequent mutation observed in the total cohort targeting 68.4% followed by frameshift deletion reported in 26.3%. Six novel variants have been discovered in this study.

Conclusions: This study demonstrates the high yield of information in DLBCL using the NGS Lymphoma panel. Results also highlight the molecular heterogeneity of DLBCL subtypes which indicates the need for further studies to make the results of the NGS more clinically relevant.

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

BACKGROUND

COVID-19 severity varies widely. Although some demographic and cardio-metabolic factors, including age and obesity, are associated with increasing risk of severe illness, the underlying mechanism(s) are uncertain.

SUBJECTS/METHODS

In a meta-analysis of three independent studies of 1471 participants in total, we investigated phenotypic and genetic factors associated with subcutaneous adipose tissue expression of Angiotensin I Converting Enzyme 2 (ACE2), measured by RNA-Seq, which acts as a receptor for SARS-CoV-2 cellular entry.

RESULTS

Lower adipose tissue ACE2 expression was associated with multiple adverse cardio-metabolic health indices, including type 2 diabetes (T2D) (P = 9.14 × 10-6), obesity status (P = 4.81 × 10-5), higher serum fasting insulin (P = 5.32 × 10-4), BMI (P = 3.94 × 10-4), and lower serum HDL levels (P = 1.92 × 10-7). ACE2 expression was also associated with estimated proportions of cell types in adipose tissue: lower expression was associated with a lower proportion of microvascular endothelial cells (P = 4.25 × 10-4) and higher proportion of macrophages (P = 2.74 × 10-5). Despite an estimated heritability of 32%, we did not identify any proximal or distal expression quantitative trait loci (eQTLs) associated with adipose tissue ACE2 expression.

CONCLUSIONS

Our results demonstrate that individuals with cardio-metabolic features known to increase risk of severe COVID-19 have lower background ACE2 levels in this highly relevant tissue. Reduced adipose tissue ACE2 expression may contribute to the pathophysiology of cardio-metabolic diseases, as well as the associated increased risk of severe COVID-19.

Variational autoencoders learn transferrable representations of metabolomics data

Dimensionality reduction approaches are commonly used for the deconvolution of high-dimensional metabolomics datasets into underlying core metabolic processes. However, current state-of-the-art methods are widely incapable of detecting nonlinearities in metabolomics data. Variational Autoencoders (VAEs) are a deep learning method designed to learn nonlinear latent representations which generalize to unseen data. Here, we trained a VAE on a large-scale metabolomics population cohort of human blood samples consisting of over 4500 individuals. We analyzed the pathway composition of the latent space using a global feature importance score, which demonstrated that latent dimensions represent distinct cellular processes. To demonstrate model generalizability, we generated latent representations of unseen metabolomics datasets on type 2 diabetes, acute myeloid leukemia, and schizophrenia and found significant correlations with clinical patient groups. Notably, the VAE representations showed stronger effects than latent dimensions derived by linear and non-linear principal component analysis. Taken together, we demonstrate that the VAE is a powerful method that learns biologically meaningful, nonlinear, and transferrable latent representations of metabolomics data.

Explainable Machine Learning for Longitudinal Multi-Omic Microbiome

Over the years, research studies have shown there is a key connection between the microbial community in the gut, genes, and immune system. Understanding this association may help discover the cause of complex chronic idiopathic disorders such as inflammatory bowel disease. Even though important efforts have been put into the field, the functions, dynamics, and causation of dysbiosis state performed by the microbial community remains unclear. Machine learning models can help elucidate important connections and relationships between microbes in the human host. Our study aims to extend the current knowledge of associations between the human microbiome and health and disease through the application of dynamic Bayesian networks to describe the temporal variation of the gut microbiota and dynamic relationships between taxonomic entities and clinical variables. We develop a set of preprocessing steps to clean, filter, select, integrate, and model informative metagenomics, metatranscriptomics, and metabolomics longitudinal data from the Human Microbiome Project. This study accomplishes novel network models with satisfactory predictive performance (accuracy = 0.648) for each inflammatory bowel disease state, validating Bayesian networks as a framework for developing interpretable models to help understand the basic ways the different biological entities (taxa, genes, metabolites) interact with each other in a given environment (human gut) over time. These findings can serve as a starting point to advance the discovery of novel therapeutic approaches and new biomarkers for precision medicine.

Analysis of Tryptophan and Its Main Metabolite Kynurenine and the Risk of Multiple Cancers Based on the Bidirectional Mendelian Randomization Analysis

Background

Tryptophan and its metabolites have been found related to various cancers, but the direction of this relationship is still unclear. The purpose of this study is to explore the causal associations of tryptophan and kynurenine with multiple cancers based on the bidirectional Mendelian randomization analysis.

Methods

The data of a genome-wide association study meta-analysis on 7,824 individuals was used to explore the genetic variants strongly associated with tryptophan and kynurenine. Genetic instruments of four specific cancers were obtained from available summary-level data of 323,590 European participants. Bidirectional Mendelian randomization analysis was conducted to examine possible causality. Sensitivity analysis was performed to test heterogeneity and horizontal pleiotropy. COX regression analysis was conducted to explore associations between dietary tryptophan and cancer mortality in NHANES 1988-1994.

Results

No evidence of any causal association of tryptophan and kynurenine with the risk of four specific cancers was shown, except for weak correlations were suggested between lung or prostate cancer and kynurenine. Multiple sensitivity analyses generated similar results. Our findings from COX regression analysis were consistent with the above results.

Conclusions

Our study did not find any causal relationship between tryptophan and kynurenine and multiple cancers. The associations still need further research.

Genetic and environmental correlational structure among metabolic syndrome endophenotypes

Metabolic syndrome (MetS) is diagnosed by the presence of high scores on three or more metabolic traits, including systolic and diastolic blood pressure (SBP, DBP), glucose and insulin levels, cholesterol and triglyceride (TG) levels, and central obesity. A diagnosis of MetS is associated with increased risk of cardiovascular disease and type 2 diabetes. The components of MetS have long been demonstrated to have substantial genetic components, but their genetic overlap is less well understood. The present paper takes a multi-prong approach to examining the extent of this genetic overlap. This includes the quantitative genetic and additive Bayesian network modeling of the large TwinsUK project and examination of the results of genome-wide association study (GWAS) of UK Biobank data through use of LD score regression and examination of the number of genes and pathways identified in the GWASes which overlap across MetS traits. Results demonstrate a modest genetic overlap, and the genetic correlations obtained from TwinsUK and UK Biobank are nearly identical. However, these correlations imply more genetic dissimilarity than similarity. Furthermore, examination of the extent of overlap in significant GWAS hits, both at the gene and pathway level, again demonstrates only modest but significant genetic overlap. This lends support to the idea that in clinical treatment of MetS, treating each of the components individually may be an important way to address MetS.

Yoghurt consumption is associated with changes in the composition of the human gut microbiome and metabolome

BACKGROUND

Yoghurt contains live bacteria that could contribute via modulation of the gut microbiota to its reported beneficial effects such as reduced body weight gain and lower incidence of type 2 diabetes. To date, the association between yoghurt consumption and the composition of the gut microbiota is underexplored. Here we used clinical variables, metabolomics, 16S rRNA and shotgun metagenomic sequencing data collected on over 1000 predominantly female UK twins to define the link between the gut microbiota and yoghurt-associated health benefits.

RESULTS

According to food frequency questionnaires (FFQ), 73% of subjects consumed yoghurt. Consumers presented a healthier diet pattern (healthy eating index: beta = 2.17 ± 0.34; P = 2.72x10-10) and improved metabolic health characterised by reduced visceral fat (beta = -28.18 ± 11.71 g; P = 0.01). According to 16S rRNA gene analyses and whole shotgun metagenomic sequencing approach consistent taxonomic variations were observed with yoghurt consumption. More specifically, we identified higher abundance of species used as yoghurt starters Streptococcus thermophilus (beta = 0.41 ± 0.051; P = 6.14x10-12) and sometimes added Bifidobacterium animalis subsp. lactis (beta = 0.30 ± 0.052; P = 1.49x10-8) in the gut of yoghurt consumers. Replication in 1103 volunteers from the LifeLines-DEEP cohort confirmed the increase of S. thermophilus among yoghurt consumers. Using food records collected the day prior to faecal sampling we showed than an increase in these two yoghurt bacteria could be transient. Metabolomics analysis revealed that B. animalis subsp. lactis was associated with 13 faecal metabolites including a 3-hydroxyoctanoic acid, known to be involved in the regulation of gut inflammation.

CONCLUSIONS

Yoghurt consumption is associated with reduced visceral fat mass and changes in gut microbiome including transient increase of yoghurt-contained species (i.e. S. thermophilus and B. lactis).

Rare variant association study of veteran twin whole-genomes links severe depression with a nonsynonymous change in the neuronal gene BHLHE22

OBJECTIVES

Major Depressive Disorder (MDD) is a complex neuropsychiatric disease with known genetic associations, but without known links to rare variation in the human genome. Here we aim to identify rare genetic variants associated with MDD using deep whole-genome sequencing data in an independent population.

METHODS

We report the sequencing of 1,688 whole genomes in a large sample of male-male Veteran twins. Depression status was classified based on a structured diagnostic interview according to DSM-III-R diagnostic criteria. Searching only rare variants in genomic regions from recent GWAS on MDD, we used the optimised sequence kernel association test and Fisher's Exact test to fine map loci associated with severe depression.

RESULTS

Our analysis identified one gene associated with severe depression, basic helix loop helix e22 (P_Adjusted = 0.03) via SKAT-O test between unrelated severely depressed cases compared to unrelated non-depressed controls. The same gene BHLHE22 had a non-silent variant rs13279074 (P_Adjusted = 0.032) based on a single variant Fisher's Exact test between unrelated severely depressed cases compared to unrelated non-depressed controls.

CONCLUSION

The gene BHLHE22 shows compelling genetic evidence of directly impacting the severe depression phenotype. Together these results advance understanding of the genetic contribution to major depressive disorder in a new cohort and link a rare variant to severe forms of the disorder.

Multi-Omic Approaches to Identify Genetic Factors in Metabolic Syndrome

Metabolic syndrome (MetS) is a highly heritable disease and a major public health burden worldwide. MetS diagnosis criteria are met by the simultaneous presence of any three of the following: high triglycerides, low HDL/high LDL cholesterol, insulin resistance, hypertension, and central obesity. These diseases act synergistically in people suffering from MetS and dramatically increase risk of morbidity and mortality due to stroke and cardiovascular disease, as well as certain cancers. Each of these component features is itself a complex disease, as is MetS. As a genetically complex disease, genetic risk factors for MetS are numerous, but not very powerful individually, often requiring specific environmental stressors for the disease to manifest. When taken together, all sequence variants that contribute to MetS disease risk explain only a fraction of the heritable variance, suggesting additional, novel loci have yet to be discovered. In this article, we will give a brief overview on the genetic concepts needed to interpret genome-wide association studies (GWAS) and quantitative trait locus (QTL) data, summarize the state of the field of MetS physiological genomics, and to introduce tools and resources that can be used by the physiologist to integrate genomics into their own research on MetS and any of its component features. There is a wealth of phenotypic and molecular data in animal models and humans that can be leveraged as outlined in this article. Integrating these multi-omic QTL data for complex diseases such as MetS provides a means to unravel the pathways and mechanisms leading to complex disease and promise for novel treatments. © 2022 American Physiological Society. Compr Physiol 12:1-40, 2022.

Identical twins carry a persistent epigenetic signature of early genome programming

Monozygotic (MZ) twins and higher-order multiples arise when a zygote splits during pre-implantation stages of development. The mechanisms underpinning this event have remained a mystery. Because MZ twinning rarely runs in families, the leading hypothesis is that it occurs at random. Here, we show that MZ twinning is strongly associated with a stable DNA methylation signature in adult somatic tissues. This signature spans regions near telomeres and centromeres, Polycomb-repressed regions and heterochromatin, genes involved in cell-adhesion, WNT signaling, cell fate, and putative human metastable epialleles. Our study also demonstrates a never-anticipated corollary: because identical twins keep a lifelong molecular signature, we can retrospectively diagnose if a person was conceived as monozygotic twin.

Protein-coding repeat polymorphisms strongly shape diverse human phenotypes

Many human proteins contain domains that vary in size or copy number because of variable numbers of tandem repeats (VNTRs) in protein-coding exons. However, the relationships of VNTRs to most phenotypes are unknown because of difficulties in measuring such repetitive elements. We developed methods to estimate VNTR lengths from whole-exome sequencing data and impute VNTR alleles into single-nucleotide polymorphism haplotypes. Analyzing 118 protein-altering VNTRs in 415,280 UK Biobank participants for association with 786 phenotypes identified some of the strongest associations of common variants with human phenotypes, including height, hair morphology, and biomarkers of health. Accounting for large-effect VNTRs further enabled fine-mapping of associations to many more protein-coding mutations in the same genes. These results point to cryptic effects of highly polymorphic common structural variants that have eluded molecular analyses to date.

An information-based framework for development national twin registry: Scoping review and focus group discussion

BACKGROUND

Registries in various clinical domains have been established in the last decades. The specific genetic structure of twins has enabled researchers to find answers to the role of genetics and the environment in medical sciences. Thus, twin registries were developed across the world to support twin studies. Our main objective was to devise a conceptual model for developing the national twin registry to ensure the success of this registry.

METHODS

In this descriptive and qualitative study, the combination of literature review and focus group discussions was applied to achieve suitable models for developing a national twin registry based on lessons learned from founded registries. The qualitative synthesis and reporting results were conducted based on the COREQ checklist.

RESULTS

According to a systematic literature review, the characteristics and employed strategies employed by established twin registries were recognized. Moreover, based on our objectives, suitable models for registry development were defined. The source of information, the different levels of data, and the information flow were determined based on this model.

CONCLUSION

Suggesting a conceptual framework for twin registry development at the national level based on the experiences of other countries could contribute to a greater understanding of twin registry implementation efficiently.

A genome-wide association study finds genetic variants associated with neck or shoulder pain in UK Biobank

Background

Common types of musculoskeletal conditions include pain in the neck and shoulder areas. This study seeks to identify the genetic variants associated with neck or shoulder pain based on a genome-wide association approach using 203 309 subjects from the UK Biobank cohort and look for replication evidence from the Generation Scotland: Scottish Family Health Study (GS:SFHS) and TwinsUK.

Methods

A genome-wide association study was performed adjusting for age, sex, BMI and nine population principal components. Significant and independent genetic variants were then sent to GS:SFHS and TwinsUK for replication.

Results

We identified three genetic loci that were associated with neck or shoulder pain in the UK Biobank samples. The most significant locus was in an intergenic region in chromosome 17, rs12453010, having P = 1.66 × 10⁻¹¹. The second most significant locus was located in the FOXP2 gene in chromosome 7 with P = 2.38 × 10⁻¹⁰ for rs34291892. The third locus was located in the LINC01572 gene in chromosome 16 with P = 4.50 × 10⁻⁸ for rs62053992. In the replication stage, among four significant and independent genetic variants, rs2049604 in the FOXP2 gene and rs62053992 in the LINC01572 gene were weakly replicated in GS:SFHS (P = 0.0240 and P = 0.0202, respectively).

Conclusions

We have identified three loci associated with neck or shoulder pain in the UK Biobank cohort, two of which were weakly supported in a replication cohort. Further evidence is needed to confirm their roles in neck or shoulder pain.

The genetic architecture of plasma kynurenine includes cardiometabolic disease mechanisms associated with the SH2B3 gene

Inflammation increases the risk of cardiometabolic disease. Delineating specific inflammatory pathways and biomarkers of their activity could identify the mechanistic underpinnings of the increased risk. Plasma levels of kynurenine, a metabolite involved in inflammation, associates with cardiometabolic disease risk. We used genetic approaches to identify inflammatory mechanisms associated with kynurenine variability and their relationship to cardiometabolic disease. We identified single-nucleotide polymorphisms (SNPs) previously associated with plasma kynurenine, including a missense-variant (rs3184504) in the inflammatory gene SH2B3/LNK. We examined the association between rs3184504 and plasma kynurenine in independent human samples, and measured kynurenine levels in SH2B3-knock-out mice and during human LPS-evoked endotoxemia. We conducted phenome scanning to identify clinical phenotypes associated with each kynurenine-related SNP and with a kynurenine polygenic score using the UK-Biobank (n = 456,422), BioVU (n = 62,303), and Electronic Medical Records and Genetics (n = 32,324) databases. The SH2B3 missense variant associated with plasma kynurenine levels and SH2B3-/- mice had significant tissue-specific differences in kynurenine levels.LPS, an acute inflammatory stimulus, increased plasma kynurenine in humans. Mendelian randomization showed increased waist-circumference, a marker of central obesity, associated with increased kynurenine, and increased kynurenine associated with C-reactive protein (CRP). We found 30 diagnoses associated (FDR q < 0.05) with the SH2B3 variant, but not with SNPs mapping to genes known to regulate tryptophan-kynurenine metabolism. Plasma kynurenine may be a biomarker of acute and chronic inflammation involving the SH2B3 pathways. Its regulation lies upstream of CRP, suggesting that kynurenine may be a biomarker of one inflammatory mechanism contributing to increased cardiometabolic disease risk.

Dietary Interventions Reduce Traditional and Novel Cardiovascular Risk Markers by Altering the Gut Microbiome and Their Metabolites

Aims: The current study investigates the role of diet in mediating the gut microbiome-cardiovascular association which has not yet been explored in humans.

Methods and Results: Using a two-arm dietary intervention study in healthy participants (N = 70), we assessed the effects of omega-3 and fibre supplementation on gut microbiome composition and short-chain fatty acid (SCFA) production. We then investigated how changes in gut microbiome composition correlated with changes in traditional cardiovascular risk factors (cholesterol, triglycerides, blood pressure), cytokines, and novel validated markers such as GlycA and ceramides, previously linked to CVD incidence and mortality. Both interventions resulted in significant drops in blood pressure, cholesterol, proinflammatory cytokines, GlycA and ceramides (all P < 0.05). Decreases in the atherogenic low-density lipoprotein triglyceride fraction, in total serum cholesterol were correlated with increases in butyric acid-production [β(SE) = −0.58 (0.06), P < 0.001; −0.53 (0.04), P < 0.001] and nominally associated with increases in some butyrogenic bacteria. Drops in GlycA were linked to increases in Bifidobacterium [β(SE) = −0.32 (0.04), P = 0.02] and other SCFAs including acetic acid [β(SE) = −0.28 (0.04), P = 0.02] and propionic acid [β(SE) = −0.3 (0.04), P = 0.02]. Additionally, we report for the first-time reductions in specific ceramide ratios that have been shown to predict CVD mortality and major adverse cardiovascular events such as d18:1/16:0, d18:0/24:0, and d18:1/24:1 which were associated with the reduction in the abundance in Colinsella and increases in Bifidobacteriuim and Coprococcus 3 and SCFAs (all P < 0.05).

Conclusion: Overall, these findings support the potential of using simple dietary interventions to alter validated biomarkers linked to cardiovascular risk via the gut microbiome composition and its metabolic functions.

Epigenome-wide association meta-analysis of DNA methylation with coffee and tea consumption

Coffee and tea are extensively consumed beverages worldwide which have received considerable attention regarding health. Intake of these beverages is consistently linked to, among others, reduced risk of diabetes and liver diseases; however, the mechanisms of action remain elusive. Epigenetics is suggested as a mechanism mediating the effects of dietary and lifestyle factors on disease onset. Here we report the results from epigenome-wide association studies (EWAS) on coffee and tea consumption in 15,789 participants of European and African-American ancestries from 15 cohorts. EWAS meta-analysis of coffee consumption reveals 11 CpGs surpassing the epigenome-wide significance threshold (P-value <1.1×10-7), which annotated to the AHRR, F2RL3, FLJ43663, HDAC4, GFI1 and PHGDH genes. Among them, cg14476101 is significantly associated with expression of the PHGDH and risk of fatty liver disease. Knockdown of PHGDH expression in liver cells shows a correlation with expression levels of genes associated with circulating lipids, suggesting a role of PHGDH in hepatic-lipid metabolism. EWAS meta-analysis on tea consumption reveals no significant association, only two CpGs annotated to CACNA1A and PRDM16 genes show suggestive association (P-value <5.0×10-6). These findings indicate that coffee-associated changes in DNA methylation levels may explain the mechanism of action of coffee consumption in conferring risk of diseases.

Recent Trends in Patient Registries for Health Services Research

Background  Patient registries are an established methodology in health services research. Since more than 150 years, registries collect information concerning groups of similar patients to answer research questions. Elaborated recommendations about an appropriate development and an efficient operation of registries are available. However, the scene changes rapidly.

Objectives  The aim of the study is to describe current trends in registry research for health services research.

Methods  Registries developed within a German funding scheme for model registries in health services research were analyzed. The observations were compared with recent recommendations of the Agency for Healthcare Research and Quality (AHRQ) on registries in the 21st century.

Results  Analyzing six registries from the funding scheme revealed the following trends: recruiting healthy individuals, representing familial or other interpersonal relationships, recording of patient-reported experiences or outcomes, accepting participants as study sites, active informing of participants, integrating the registry with other data collections, and transferring data from the registry to electronic patient records. This list partly complies with the issues discussed by the AHRQ. The AHRQ structured its ideas in five chapters, increasing focus on the patient, engaging patients as partners, digital health and patient registries, direct-to-patient registry, and registry networks.

Conclusion  For the near future, it can be said that the concept and the design of a registry should place the patient in the center. Registries will be increasingly linked together and interconnected with other data collections. New challenges arise regarding the management of data quality and the interpretation of results from less controlled settings. Here, further research related to the methodology of registries is needed.

Statistical integration of two omics datasets using GO2PLS

BACKGROUND

Nowadays, multiple omics data are measured on the same samples in the belief that these different omics datasets represent various aspects of the underlying biological systems. Integrating these omics datasets will facilitate the understanding of the systems. For this purpose, various methods have been proposed, such as Partial Least Squares (PLS), decomposing two datasets into joint and residual subspaces. Since omics data are heterogeneous, the joint components in PLS will contain variation specific to each dataset. To account for this, Two-way Orthogonal Partial Least Squares (O2PLS) captures the heterogeneity by introducing orthogonal subspaces and better estimates the joint subspaces. However, the latent components spanning the joint subspaces in O2PLS are linear combinations of all variables, while it might be of interest to identify a small subset relevant to the research question. To obtain sparsity, we extend O2PLS to Group Sparse O2PLS (GO2PLS) that utilizes biological information on group structures among variables and performs group selection in the joint subspace.

RESULTS

The simulation study showed that introducing sparsity improved the feature selection performance. Furthermore, incorporating group structures increased robustness of the feature selection procedure. GO2PLS performed optimally in terms of accuracy of joint score estimation, joint loading estimation, and feature selection. We applied GO2PLS to datasets from two studies: TwinsUK (a population study) and CVON-DOSIS (a small case-control study). In the first, we incorporated biological information on the group structures of the methylation CpG sites when integrating the methylation dataset with the IgG glycomics data. The targeted genes of the selected methylation groups turned out to be relevant to the immune system, in which the IgG glycans play important roles. In the second, we selected regulatory regions and transcripts that explained the covariance between regulomics and transcriptomics data. The corresponding genes of the selected features appeared to be relevant to heart muscle disease.

CONCLUSIONS

GO2PLS integrates two omics datasets to help understand the underlying system that involves both omics levels. It incorporates external group information and performs group selection, resulting in a small subset of features that best explain the relationship between two omics datasets for better interpretability.

Welfare Genome Project: A Participatory Korean Personal Genome Project With Free Health Check-Up and Genetic Report Followed by Counseling

The Welfare Genome Project (WGP) provided 1,000 healthy Korean volunteers with detailed genetic and health reports to test the social perception of integrating personal genetic and healthcare data at a large-scale. WGP was launched in 2016 in the Ulsan Metropolitan City as the first large-scale genome project with public participation in Korea. The project produced a set of genetic materials, genotype information, clinical data, and lifestyle survey answers from participants aged 20–96. As compensation, the participants received a free general health check-up on 110 clinical traits, accompanied by a genetic report of their genotypes followed by genetic counseling. In a follow-up survey, 91.0% of the participants indicated that their genetic reports motivated them to improve their health. Overall, WGP expanded not only the general awareness of genomics, DNA sequencing technologies, bioinformatics, and bioethics regulations among all the parties involved, but also the general public’s understanding of how genome projects can indirectly benefit their health and lifestyle management. WGP established a data construction framework for not only scientific research but also the welfare of participants. In the future, the WGP framework can help lay the groundwork for a new personalized healthcare system that is seamlessly integrated with existing public medical infrastructure.

Neural and Genetic Bases for Human Ability Traits

The judgement of human ability is ubiquitous, from school admissions to job performance reviews. The exact make-up of ability traits, however, is often narrowly defined and lacks a comprehensive basis. We attempt to simplify the spectrum of human ability, similar to how five personality traits are widely believed to describe most personalities. Finding such a basis for human ability would be invaluable since neuropsychiatric disease diagnoses and symptom severity are commonly related to such differences in performance. Here, we identified four underlying ability traits within the National Institutes of Health Toolbox normative data (n = 1, 369): (1) Motor-endurance, (2) Emotional processing, (3) Executive and cognitive function, and (4) Social interaction. We used the Human Connectome Project young adult dataset (n = 778) to show that Motor-endurance and Executive and cognitive function were reliably associated with specific brain functional networks (r ² = 0.305 ± 0.021), and the biological nature of these ability traits was also shown by calculating their heritability (31 and 49%, respectively) from twin data.

Estimating pulse wave velocity from the radial pressure wave using machine learning algorithms

One of the European gold standard measurement of vascular ageing, a risk factor for cardiovascular disease, is the carotid-femoral pulse wave velocity (cfPWV), which requires an experienced operator to measure pulse waves at two sites. In this work, two machine learning pipelines were proposed to estimate cfPWV from the peripheral pulse wave measured at a single site, the radial pressure wave measured by applanation tonometry. The study populations were the Twins UK cohort containing 3,082 subjects aged from 18 to 110 years, and a database containing 4,374 virtual subjects aged from 25 to 75 years. The first pipeline uses Gaussian process regression to estimate cfPWV from features extracted from the radial pressure wave using pulse wave analysis. The mean difference and upper and lower limits of agreement (LOA) of the estimation on the 924 hold-out test subjects from the Twins UK cohort were 0.2 m/s, and 3.75 m/s & -3.34 m/s, respectively. The second pipeline uses a recurrent neural network (RNN) to estimate cfPWV from the entire radial pressure wave. The mean difference and upper and lower LOA of the estimation on the 924 hold-out test subjects from the Twins UK cohort were 0.05 m/s, and 3.21 m/s & -3.11m/s, respectively. The percentage error of the RNN estimates on the virtual subjects increased by less than 2% when adding 20% of random noise to the pressure waveform. These results show the possibility of assessing the vascular ageing using a single peripheral pulse wave (e.g. the radial pressure wave), instead of cfPWV. The proposed code for the machine learning pipelines is available from the following online depository (https://github.com/WeiweiJin/Estimate-Cardiovascular-Risk-from-Pulse-Wave-Signal).