Genome-wide characterization of circulating metabolic biomarkers.

Genome-wide association analyses using high-throughput metabolomics platforms have led to novel insights into the biology of human metabolism1-7. This detailed knowledge of the genetic determinants of systemic metabolism has been pivotal for uncovering how genetic pathways influence biological mechanisms and complex diseases8-11. Here we present a genome-wide association study for 233 circulating metabolic traits quantified by nuclear magnetic resonance spectroscopy in up to 136,016 participants from 33 cohorts. We identify more than 400 independent loci and assign probable causal genes at two-thirds of these using manual curation of plausible biological candidates. We highlight the importance of sample and participant characteristics that can have significant effects on genetic associations. We use detailed metabolic profiling of lipoprotein- and lipid-associated variants to better characterize how known lipid loci and novel loci affect lipoprotein metabolism at a granular level. We demonstrate the translational utility of comprehensively phenotyped molecular data, characterizing the metabolic associations of intrahepatic cholestasis of pregnancy. Finally, we observe substantial genetic pleiotropy for multiple metabolic pathways and illustrate the importance of careful instrument selection in Mendelian randomization analysis, revealing a putative causal relationship between acetone and hypertension. Our publicly available results provide a foundational resource for the community to examine the role of metabolism across diverse diseases.

A genetic association study of circulating coagulation factor VIII and von Willebrand factor levels.

ABSTRACT: Coagulation factor VIII (FVIII) and its carrier protein von Willebrand factor (VWF) are critical to coagulation and platelet aggregation. We leveraged whole-genome sequence data from the Trans-Omics for Precision Medicine (TOPMed) program along with TOPMed-based imputation of genotypes in additional samples to identify genetic associations with circulating FVIII and VWF levels in a single-variant meta-analysis, including up to 45 289 participants. Gene-based aggregate tests were implemented in TOPMed. We identified 3 candidate causal genes and tested their functional effect on FVIII release from human liver endothelial cells (HLECs) and VWF release from human umbilical vein endothelial cells. Mendelian randomization was also performed to provide evidence for causal associations of FVIII and VWF with thrombotic outcomes. We identified associations (P < 5 × 10-9) at 7 new loci for FVIII (ST3GAL4, CLEC4M, B3GNT2, ASGR1, F12, KNG1, and TREM1/NCR2) and 1 for VWF (B3GNT2). VWF, ABO, and STAB2 were associated with FVIII and VWF in gene-based analyses. Multiphenotype analysis of FVIII and VWF identified another 3 new loci, including PDIA3. Silencing of B3GNT2 and the previously reported CD36 gene decreased release of FVIII by HLECs, whereas silencing of B3GNT2, CD36, and PDIA3 decreased release of VWF by HVECs. Mendelian randomization supports causal association of higher FVIII and VWF with increased risk of thrombotic outcomes. Seven new loci were identified for FVIII and 1 for VWF, with evidence supporting causal associations of FVIII and VWF with thrombotic outcomes. B3GNT2, CD36, and PDIA3 modulate the release of FVIII and/or VWF in vitro.

Normal range CAG repeat size variations in the HTT gene are associated with an adverse lipoprotein profile partially mediated by body mass index.

Tandem cytosine-adenine-guanine (CAG) repeat sizes of 36 or more in the huntingtin gene (HTT) cause Huntington's disease (HD). Apart from neuropsychiatric complications, the disease is also accompanied by metabolic dysregulation and weight loss, which contribute to a progressive functional decline. Recent studies also reported an association between repeats below the pathogenic threshold (<36) for HD and body mass index (BMI), suggesting that HTT repeat sizes in the non-pathogenic range are associated with metabolic dysregulation. In this study, we hypothesized that HTT repeat sizes < 36 are associated with metabolite levels, possibly mediated through reduced BMI. We pooled data from three European cohorts (n = 10 228) with genotyped HTT CAG repeat size and metabolomic measurements. All 145 metabolites were measured on the same targeted platform in all studies. Multilevel mixed-effects analysis using the CAG repeat size in HTT identified 67 repeat size metabolite associations. Overall, the metabolomic profile associated with larger CAG repeat sizes in HTT were unfavorable-similar to those of higher risk of coronary artery disease and type 2 diabetes-and included elevated levels of amino acids, fatty acids, low-density lipoprotein (LDL)-, very low-density lipoprotein- and intermediate density lipoprotein (IDL)-related metabolites while with decreased levels of very large high-density lipoprotein (HDL)-related metabolites. Furthermore, the associations of 50 metabolites, in particular, specific very large HDL-related metabolites, were mediated by lower BMI. However, no mediation effect was found for 17 metabolites related to LDL and IDL. In conclusion, our findings indicate that large non-pathogenic CAG repeat sizes in HTT are associated with an unfavorable metabolomic profile despite their association with a lower BMI.

Bioinformatic analysis and experimental validation of cuproptosis-related LncRNA as a novel biomarker for prognosis and immunotherapy of oral squamous cell carcinoma.

BACKGROUND: The novel form of regulatory cell death, cuproptosis, is characterized by proteotoxicity, which ultimately leads to cell death. Its targeting has emerged as a promising therapeutic approach for oral squamous cell carcinoma (OSCC). Long noncoding RNAs (lncRNAs) participate in epigenetic regulation and have been linked to the progression, prognosis, and treatment of OSCC. Thus, this study aimed to identify new cuproptosis-related lncRNAs (CRLs), establish predictive models for clinical prognosis, immune response, and drug sensitivity, and provide novel insights into immune escape and tumor drug resistance. METHODS: The present study screened eight CRLs (THAP9-AS1, STARD4-AS1, WDFY3-AS2, LINC00847, CDKN2A-DT, AL132800.1, GCC2-AS1, AC005746.1) using Lasso Cox regression analysis to develop an eight-CRL prognostic model. Patients were categorized into high- and low-risk groups using risk scores. To evaluate the predictive ability of the model, Kaplan-Meier analysis, ROC curves, and nomograms were employed. Furthermore, the study investigated the differences in immune function and anticancer drug sensitivity between the high- and low-risk groups. To validate the expression of CRLs in the model, OSCC cell lines were subjected to quantitative real-time fluorescence PCR (qRT-PCR). RESULTS: The results of the study showed that the high-risk group had a shorter overall survival (OS) time in OSCC patients. Cox regression analysis demonstrated that the high-risk score was an independent risk factor for a poor prognosis. The validity of the model was confirmed using ROC curve analysis, and a nomogram was developed to predict the prognosis of OSCC patients. Furthermore, patients in the high-risk group with high TMB had a poorer prognosis. Patients in the low-risk group responded better to immunotherapy than those in the high-risk group. Additionally, the risk scores were significantly associated with drug sensitivity in OSCC patients. Finally, the findings of qRT-PCR supported the reliability of the proposed risk model. CONCLUSION: The study identified and established the 8-CRL model, which represents a novel pathway of lncRNA regulation of cuproptosis in OSCC. This model provides guidance for the prognosis and treatment of OSCC and offers a new insight into immune escape and tumor drug resistance.

Cross-Ancestry Investigation of Venous Thromboembolism Genomic Predictors.

BACKGROUND: Venous thromboembolism (VTE) is a life-threatening vascular event with environmental and genetic determinants. Recent VTE genome-wide association studies (GWAS) meta-analyses involved nearly 30 000 VTE cases and identified up to 40 genetic loci associated with VTE risk, including loci not previously suspected to play a role in hemostasis. The aim of our research was to expand discovery of new genetic loci associated with VTE by using cross-ancestry genomic resources. METHODS: We present new cross-ancestry meta-analyzed GWAS results involving up to 81 669 VTE cases from 30 studies, with replication of novel loci in independent populations and loci characterization through in silico genomic interrogations. RESULTS: In our genetic discovery effort that included 55 330 participants with VTE (47 822 European, 6320 African, and 1188 Hispanic ancestry), we identified 48 novel associations, of which 34 were replicated after correction for multiple testing. In our combined discovery-replication analysis (81 669 VTE participants) and ancestry-stratified meta-analyses (European, African, and Hispanic), we identified another 44 novel associations, which are new candidate VTE-associated loci requiring replication. In total, across all GWAS meta-analyses, we identified 135 independent genomic loci significantly associated with VTE risk. A genetic risk score of the significantly associated loci in Europeans identified a 6-fold increase in risk for those in the top 1% of scores compared with those with average scores. We also identified 31 novel transcript associations in transcriptome-wide association studies and 8 novel candidate genes with protein quantitative-trait locus Mendelian randomization analyses. In silico interrogations of hemostasis and hematology traits and a large phenome-wide association analysis of the 135 GWAS loci provided insights to biological pathways contributing to VTE, with some loci contributing to VTE through well-characterized coagulation pathways and others providing new data on the role of hematology traits, particularly platelet function. Many of the replicated loci are outside of known or currently hypothesized pathways to thrombosis. CONCLUSIONS: Our cross-ancestry GWAS meta-analyses identified new loci associated with VTE. These findings highlight new pathways to thrombosis and provide novel molecules that may be useful in the development of improved antithrombosis treatments.

Association between ABO haplotypes and the risk of venous thrombosis: impact on disease risk estimation.

Genetic risk score (GRS) analysis is a popular approach to derive individual risk prediction models for complex diseases. In venous thrombosis (VT), such type of analysis shall integrate information at the ABO blood group locus, which is one of the major susceptibility loci. However, there is no consensus about which single nucleotide polymorphisms (SNPs) must be investigated when properly assessing association between ABO locus and VT risk. Using comprehensive haplotype analyses of ABO blood group tagging SNPs in 5425 cases and 8445 controls from 6 studies, we demonstrate that using only rs8176719 (tagging O1) to correctly assess the impact of ABO locus on VT risk is suboptimal, because 5% of rs8176719-delG carriers do not have an increased risk of developing VT. Instead, we recommend the use of 4 SNPs, rs2519093 (tagging A1), rs1053878 (A2), rs8176743 (B), and rs41302905 (O2), when assessing the impact of ABO locus on VT risk to avoid any risk misestimation. Compared with the O1 haplotype, the A2 haplotype is associated with a modest increase in VT risk (odds ratio, âˆ¼1.2), the A1 and B haplotypes are associated with an âˆ¼1.8-fold increased risk, whereas the O2 haplotype tends to be slightly protective (odds ratio, âˆ¼0.80). In addition, although the A1 and B blood groups are associated with increased von Willebrand factor and factor VIII plasma levels, only the A1 blood group is associated with ICAM levels, but in an opposite direction, leaving additional avenues to be explored to fully understand the spectrum of biological effects mediated by ABO locus on cardiovascular traits.

Hepatic triglyceride content is intricately associated with numerous metabolites and biochemical pathways.

BACKGROUND AND AIMS: Non-alcoholic fatty liver disease (NAFLD) is characterized by the pathological accumulation of triglycerides in hepatocytes and is associated with insulin resistance, atherogenic dyslipidaemia and cardiometabolic diseases. Thus far, the extent of metabolic dysregulation associated with hepatic triglyceride accumulation has not been fully addressed. In this study, we aimed to identify metabolites associated with hepatic triglyceride content (HTGC) and map these associations using network analysis. METHODS: To gain insight in the spectrum of metabolites associated with hepatic triglyceride accumulation, we performed a comprehensive plasma metabolomics screening of 1363 metabolites in apparently healthy middle aged (age 45-65) individuals (N = 496) in whom HTGC was measured by proton magnetic resonance spectroscopy. An atlas of metabolite-HTGC associations, based on univariate results, was created using correlation-based Gaussian graphical model (GGM) and genome scale metabolic model network analyses. Pathways associated with the clinical prognosis marker fibrosis 4 (FIB-4) index were tested using a closed global test. RESULTS: Our analyses revealed that 118 metabolites were univariately associated with HTGC (p-value <6.59 × 10-5 ), including 106 endogenous, 1 xenobiotic and 11 partially characterized/uncharacterized metabolites. These associations were mapped to several biological pathways including branched amino acids (BCAA), diglycerols, sphingomyelin, glucosyl-ceramide and lactosyl-ceramide. We also identified a novel possible HTGC-related pathway connecting glutamate, metabolonic lactone sulphate and X-15245 using the GGM network. These pathways were confirmed to be associated with the FIB-4 index as well. The full interactive metabolite-HTGC atlas is provided online: https://tofaquih.github.io/AtlasLiver/. CONCLUSIONS: The combined network and pathway analyses indicated extensive associations between BCAA and the lipids pathways with HTGC and the FIB-4 index. Moreover, we report a novel pathway glutamate-metabolonic lactone sulphate-X-15245 with a potential strong association with HTGC. These findings can aid elucidating HTGC metabolomic profiles and provide insight into novel drug targets for fibrosis-related outcomes.

Altered high-density lipoprotein composition is associated with risk for complications in type 2 diabetes mellitus in South Asian descendants: A cross-sectional, case-control study on lipoprotein subclass profiling.

BACKGROUND: Composition of high-density lipoproteins (HDL) is emerging as an important determinant in the development of microvascular complications in type 2 diabetes mellitus (T2DM). Dutch South Asian (DSA) individuals with T2DM display an increased risk of microvascular complications compared with Dutch white Caucasian (DwC) individuals with T2DM. In this study, we aimed to investigate whether changes in HDL composition associate with increased microvascular risk in this ethnic group and lead to new lipoprotein biomarkers. MATERIALS AND METHODS: Using 1 H nuclear magnetic resonance spectroscopy and Bruker IVDr Lipoprotein Subclass Analysis (B.I.LISA) software, plasma lipoprotein changes were determined in 51 healthy individuals (30 DwC, 21 DSA) and 92 individuals with T2DM (45 DwC, 47 DSA) in a cross-sectional, case-control study. Differential HDL subfractions were investigated using multinomial logistic regression analyses, adjusting for possible confounders including BMI and diabetes duration. RESULTS: We identified HDL compositional differences between healthy and diabetic individuals in both ethnic groups. Specifically, levels of apolipoprotein A2 and HDL-4 subfractions were lower in DSA compared with DwC with T2DM. Apolipoprotein A2 and HDL-4 subfractions also negatively correlated with waist circumference, waist-to-hip ratio, haemoglobin A1c, glucose levels and disease duration in DSA with T2DM, and associated with increased incidence of microvascular complications. CONCLUSION: While HDL composition differed between controls and T2DM in both ethnic groups, the lower levels of lipid content in the smallest HDL subclass (HDL-4) in DSA with T2DM appeared to be more clinically relevant, with higher odds of having diabetes-related pan-microvascular complications such as retinopathy and neuropathy. These typical differences in HDL could be used as ethnicity-specific T2DM biomarkers.

Network pharmacological analysis and experimental study of cucurbitacin B in oral squamous cell carcinoma.

Oral squamous cell carcinoma (OSCC) is a malignant tumor with a high incidence and poor prognosis. Cucurbitacin B (CuB) is a tetracyclic triterpenoid small-molecule compound extracted from plants, such as Cucurbitaceae and Brassicaceae, which has powerful anticancer effects. However, the effect and mechanism of CuB on OSCC remain unclear. Within the framework of the current study, network pharmacology was used to analyze the relationship between CuB and OSCC. The network pharmacology analysis showed that CuB and OSCC share 134 common targets; among them, PIK3R1, SRC, STAT3, AKT1, and MAPK1 are the key targets. The molecular docking analysis showed that CuB binds five target proteins. The results of the enrichment analysis showed that CuB exerted effects on OSCC through various pathways; of these pathways, PI3K-AKT was the most important pathway. The results of the in vitro cell experiments showed that CuB could inhibit the proliferation and migration of SCC25 and CAL27 cells, block the cell cycle in the G2 phase, induce cell apoptosis, and regulate the protein expression of the PI3K-AKT signaling pathway. The results of the in vivo animal experiments showed that CuB could inhibit 4NQO-induced oral cancer in mice. Therefore, network pharmacology, molecular docking, cell experiments, and animal experiments showed that CuB could play a role in OSCC by regulating multiple targets and pathways.

Functional Characterization of a (E)-ß-Ocimene Synthase Gene Contributing to the Defense against Spodoptera litura.

Soybean is a worldwide crop that offers valuable proteins, fatty acids, and phytonutrients to humans but is always damaged by insect pests or pathogens. Plants have captured sophisticated defense mechanisms in resisting the attack of insects and pathogens. How to protect soybean in an environment- or human-friendly way or how to develop plant-based pest control is a hotpot. Herbivore-induced plant volatiles that are released by multiple plant species have been assessed in multi-systems against various insects, of which (E)-ß-ocimene has been reported to show anti-insect function in a variety of plants, including soybean. However, the responsible gene in soybean is unknown, and its mechanism of synthesis and anti-insect properties lacks comprehensive assessment. In this study, (E)-ß-ocimene was confirmed to be induced by Spodoptera litura treatment. A plastidic localized monoterpene synthase gene, designated as GmOCS, was identified to be responsible for the biosynthesis of (E)-ß-ocimene through genome-wide gene family screening and in vitro and in vivo assays. Results from transgenic soybean and tobacco confirmed that (E)-ß-ocimene catalyzed by GmOCS had pivotal roles in repelling a S. litura attack. This study advances the understanding of (E)-ß-ocimene synthesis and its function in crops, as well as provides a good candidate for further anti-insect soybean improvement.

Associations of Mitochondrial and Nuclear Mitochondrial Variants and Genes with Seven Metabolic Traits.

Mitochondria (MT), the major site of cellular energy production, are under dual genetic control by 37 mitochondrial DNA (mtDNA) genes and numerous nuclear genes (MT-nDNA). In the CHARGEmtDNA+ Consortium, we studied genetic associations of mtDNA and MT-nDNA associations with body mass index (BMI), waist-hip-ratio (WHR), glucose, insulin, HOMA-B, HOMA-IR, and HbA1c. This 45-cohort collaboration comprised 70,775 (insulin) to 170,202 (BMI) pan-ancestry individuals. Validation and imputation of mtDNA variants was followed by single-variant and gene-based association testing. We report two significant common variants, one in MT-ATP6 associated (p ≤ 5E-04) with WHR and one in the D-loop with glucose. Five rare variants in MT-ATP6, MT-ND5, and MT-ND6 associated with BMI, WHR, or insulin. Gene-based meta-analysis identified MT-ND3 associated with BMI (p ≤ 1E-03). We considered 2,282 MT-nDNA candidate gene associations compiled from online summary results for our traits (20 unique studies with 31 dataset consortia's genome-wide associations [GWASs]). Of these, 109 genes associated (p ≤ 1E-06) with at least 1 of our 7 traits. We assessed regulatory features of variants in the 109 genes, cis- and trans-gene expression regulation, and performed enrichment and protein-protein interactions analyses. Of the identified mtDNA and MT-nDNA genes, 79 associated with adipose measures, 49 with glucose/insulin, 13 with risk for type 2 diabetes, and 18 with cardiovascular disease, indicating for pleiotropic effects with health implications. Additionally, 21 genes related to cholesterol, suggesting additional important roles for the genes identified. Our results suggest that mtDNA and MT-nDNA genes and variants reported make important contributions to glucose and insulin metabolism, adipocyte regulation, diabetes, and cardiovascular disease.

Exome-Derived Adiponectin-Associated Variants Implicate Obesity and Lipid Biology.

Circulating levels of adiponectin, an adipocyte-secreted protein associated with cardiovascular and metabolic risk, are highly heritable. To gain insights into the biology that regulates adiponectin levels, we performed an exome array meta-analysis of 265,780 genetic variants in 67,739 individuals of European, Hispanic, African American, and East Asian ancestry. We identified 20 loci associated with adiponectin, including 11 that had been reported previously (p < 2 × 10-7). Comparison of exome array variants to regional linkage disequilibrium (LD) patterns and prior genome-wide association study (GWAS) results detected candidate variants (r2 > .60) spanning as much as 900 kb. To identify potential genes and mechanisms through which the previously unreported association signals act to affect adiponectin levels, we assessed cross-trait associations, expression quantitative trait loci in subcutaneous adipose, and biological pathways of nearby genes. Eight of the nine loci were also associated (p < 1 × 10-4) with at least one obesity or lipid trait. Candidate genes include PRKAR2A, PTH1R, and HDAC9, which have been suggested to play roles in adipocyte differentiation or bone marrow adipose tissue. Taken together, these findings provide further insights into the processes that influence circulating adiponectin levels.

Gene-educational attainment interactions in a multi-ancestry genome-wide meta-analysis identify novel blood pressure loci.

Educational attainment is widely used as a surrogate for socioeconomic status (SES). Low SES is a risk factor for hypertension and high blood pressure (BP). To identify novel BP loci, we performed multi-ancestry meta-analyses accounting for gene-educational attainment interactions using two variables, "Some College" (yes/no) and "Graduated College" (yes/no). Interactions were evaluated using both a 1 degree of freedom (DF) interaction term and a 2DF joint test of genetic and interaction effects. Analyses were performed for systolic BP, diastolic BP, mean arterial pressure, and pulse pressure. We pursued genome-wide interrogation in Stage 1 studies (N = 117 438) and follow-up on promising variants in Stage 2 studies (N = 293 787) in five ancestry groups. Through combined meta-analyses of Stages 1 and 2, we identified 84 known and 18 novel BP loci at genome-wide significance level (P < 5 × 10-8). Two novel loci were identified based on the 1DF test of interaction with educational attainment, while the remaining 16 loci were identified through the 2DF joint test of genetic and interaction effects. Ten novel loci were identified in individuals of African ancestry. Several novel loci show strong biological plausibility since they involve physiologic systems implicated in BP regulation. They include genes involved in the central nervous system-adrenal signaling axis (ZDHHC17, CADPS, PIK3C2G), vascular structure and function (GNB3, CDON), and renal function (HAS2 and HAS2-AS1, SLIT3). Collectively, these findings suggest a role of educational attainment or SES in further dissection of the genetic architecture of BP.

Genome-wide associations for birth weight and correlations with adult disease.

Birth weight (BW) has been shown to be influenced by both fetal and maternal factors and in observational studies is reproducibly associated with future risk of adult metabolic diseases including type 2 diabetes (T2D) and cardiovascular disease. These life-course associations have often been attributed to the impact of an adverse early life environment. Here, we performed a multi-ancestry genome-wide association study (GWAS) meta-analysis of BW in 153,781 individuals, identifying 60 loci where fetal genotype was associated with BW (P < 5 × 10-8). Overall, approximately 15% of variance in BW was captured by assays of fetal genetic variation. Using genetic association alone, we found strong inverse genetic correlations between BW and systolic blood pressure (Rg = -0.22, P = 5.5 × 10-13), T2D (Rg = -0.27, P = 1.1 × 10-6) and coronary artery disease (Rg = -0.30, P = 6.5 × 10-9). In addition, using large -cohort datasets, we demonstrated that genetic factors were the major contributor to the negative covariance between BW and future cardiometabolic risk. Pathway analyses indicated that the protein products of genes within BW-associated regions were enriched for diverse processes including insulin signalling, glucose homeostasis, glycogen biosynthesis and chromatin remodelling. There was also enrichment of associations with BW in known imprinted regions (P = 1.9 × 10-4). We demonstrate that life-course associations between early growth phenotypes and adult cardiometabolic disease are in part the result of shared genetic effects and identify some of the pathways through which these causal genetic effects are mediated.

Vitamin D/VDR signaling inhibits colitis by suppressing HIF-1α activation in colonic epithelial cells.

Vitamin D/vitamin D receptor (VDR) signaling is reported to have a protective effect on the onset or progression of inflammatory bowel diseases (IBD), and hypoxia-inducible factor 1α (HIF-1α) activation is demonstrated to be closely associated with chemical-induced colitis. However, the association between vitamin D/VDR signaling and HIF-1α on IBD development remains a mystery. Here, we showed that HIF-1α expression was largely increased in the colonic epithelial cells of diseased tissues from patients with ulcerative colitis (UC). Consistently, HIF-1α activation was also improved in colonic epithelial cells upon TNFα treatment in a NF-κB pathway-dependent manner. HIF-1α inhibitors treatments ameliorated 2,4,6-trinitrobenzenesulfonic acid (TNBS)- or dextran sulfate sodium (DSS)-induced colitis in animal models. In cell or colitis animal models, vitamin D/VDR signaling suppressed HIF-1α overexpression in colonic epithelial cells via regulating NF-κB pathway, resulting in the inhibition of IFNγ and IL-1ß overproductions in these cells. Collectively, these data suggest that vitamin D/VDR signaling relieves colitis development in animal models, at least in part, by suppressing HIF-1α expression in colonic epithelial cells.NEW & NOTEWORTHY This study demonstrates vitamin D/VDR signaling inhibits colitis by suppressing HIF-1α activation in colonic epithelial cells. Since the effect of vitamin D/VDR signaling is only apparent on patients who seem to be vitamin D deficient, the benefits of vitamin D supplementation in patients who are not vitamin D deficient need to be proven.

Life-Course Genome-wide Association Study Meta-analysis of Total Body BMD and Assessment of Age-Specific Effects.

Bone mineral density (BMD) assessed by DXA is used to evaluate bone health. In children, total body (TB) measurements are commonly used; in older individuals, BMD at the lumbar spine (LS) and femoral neck (FN) is used to diagnose osteoporosis. To date, genetic variants in more than 60 loci have been identified as associated with BMD. To investigate the genetic determinants of TB-BMD variation along the life course and test for age-specific effects, we performed a meta-analysis of 30 genome-wide association studies (GWASs) of TB-BMD including 66,628 individuals overall and divided across five age strata, each spanning 15 years. We identified variants associated with TB-BMD at 80 loci, of which 36 have not been previously identified; overall, they explain approximately 10% of the TB-BMD variance when combining all age groups and influence the risk of fracture. Pathway and enrichment analysis of the association signals showed clustering within gene sets implicated in the regulation of cell growth and SMAD proteins, overexpressed in the musculoskeletal system, and enriched in enhancer and promoter regions. These findings reveal TB-BMD as a relevant trait for genetic studies of osteoporosis, enabling the identification of variants and pathways influencing different bone compartments. Only variants in ESR1 and close proximity to RANKL showed a clear effect dependency on age. This most likely indicates that the majority of genetic variants identified influence BMD early in life and that their effect can be captured throughout the life course.

Higher thyrotropin leads to unfavorable lipid profile and somewhat higher cardiovascular disease risk: evidence from multi-cohort Mendelian randomization and metabolomic profiling.

BACKGROUND: Observational studies suggest interconnections between thyroid status, metabolism, and risk of coronary artery disease (CAD), but causality remains to be proven. The present study aimed to investigate the potential causal relationship between thyroid status and cardiovascular disease and to characterize the metabolomic profile associated with thyroid status. METHODS: Multi-cohort two-sample Mendelian randomization (MR) was performed utilizing genome-wide significant variants as instruments for standardized thyrotropin (TSH) and free thyroxine (fT4) within the reference range. Associations between TSH and fT4 and metabolic profile were investigated in a two-stage manner: associations between TSH and fT4 and the full panel of 161 metabolomic markers were first assessed hypothesis-free, then directional consistency was assessed through Mendelian randomization, another metabolic profile platform, and in individuals with biochemically defined thyroid dysfunction. RESULTS: Circulating TSH was associated with 52/161 metabolomic markers, and fT4 levels were associated with 21/161 metabolomic markers among 9432 euthyroid individuals (median age varied from 23.0 to 75.4 years, 54.5% women). Positive associations between circulating TSH levels and concentrations of very low-density lipoprotein subclasses and components, triglycerides, and triglyceride content of lipoproteins were directionally consistent across the multivariable regression, MR, metabolomic platforms, and for individuals with hypo- and hyperthyroidism. Associations with fT4 levels inversely reflected those observed with TSH. Among 91,810 CAD cases and 656,091 controls of European ancestry, per 1-SD increase of genetically determined TSH concentration risk of CAD increased slightly, but not significantly, with an OR of 1.03 (95% CI 0.99-1.07; p value 0.16), whereas higher genetically determined fT4 levels were not associated with CAD risk (OR 1.00 per SD increase of fT4; 95% CI 0.96-1.04; p value 0.59). CONCLUSIONS: Lower thyroid status leads to an unfavorable lipid profile and a somewhat increased cardiovascular disease risk.

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

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

miR-3065-3p promotes stemness and metastasis by targeting CRLF1 in colorectal cancer.

BACKGROUND: Colorectal cancer is one of the most common malignancy in the world. It has been reported that cancer stem cells (CSCs) serve as the primary drivers of tumorigenesis and tumor progression. There is an urgent need to explore novel molecules that regulate CSCs or their signatures. Increasing evidence has shown that miRNAs are involved in tumorigenesis and progression. Here, we aim to explore the regulatory effect and mechanism of miR-3065-3p on the stemness of colorectal cancer. METHODS: The expression of miR-3065-3p in colorectal cancer and the association of miR-3065-3p expression with prognosis of patients with colorectal cancer were analyzed using TCGA dataset or clinical cases. Gain or loss of function in different models, including colorectal cancer cell lines and orthotopic xenograft or liver metastatic mouse model, were used to investigate the effects of miR-3065-3p on colorectal cancer stemness and metastasis in vitro and in vivo. Cancer stemness was analyzed by detecting the ability of migration and invasion, NANOG, OCT4, and SOX2 expression, ALDH activity and sphere formation. In addition, the interaction of miR-3065-3p and cytokine receptor-like factor 1 (CRLF1) was analyzed theoretically and identified by the luciferase reporter assay. Moreover, the correlation between CRLF1 expression and miR-3065-3p was analyzed in colorectal cancer tissues. Finally, the effect of CRLF1 on the stemness and metastasis of colorectal cancer in vitro and in vivo was assessed. RESULTS: In this report, we found that miR-3065-3p was overexpressed in colorectal cancer and that its high expression was associated with poor prognosis of patients with colorectal cancer. miR-3065-3p promotes the stemness and metastasis of colorectal cancer. Furthermore, CRLF1 was the downstream target of miR-3065-3p and inhibited the stemness of colorectal cancer. In addition, CRLF1 expression was negatively correlated with miR-3065-3p in colorectal cancer tissues. And, CRLF1 mediated the effects of miR-3065-3p on promoting stemness of colorectal cancer cells. CONCLUSION: Our data suggest that miR-3065-3p promoted the stemness and metastasis of colorectal cancer by targeting CRLF1. miR-3065-3p might serve as a promising prognostic marker as well as a therapeutic target for colorectal cancer.

MYB transcription factors GmMYBA2 and GmMYBR function in a feedback loop to control pigmentation of seed coat in soybean.

Soybean has undergone extensive selection pressures for seed nutrient composition and seed color during domestication, but the major genetic loci controlling seed coat color have not been completely understood, and the transcriptional regulation relationship among the loci remains elusive. Here, two major regulators, GmMYBA2 and GmMYBR, were functionally characterized as an anthocyanin activator and repressor, respectively. Ectopic expression of GmMYBA2 in soybean hairy roots conferred the enhanced accumulation of delphinidin and cyanidin types of anthocyanins in W1t and w1T backgrounds, respectively, through activating anthocyanin biosynthetic genes in the reported loci. The seed coat pigmentation of GmMYBA2-overexpressing transgenic plants in the W1 background mimicked the imperfect black phenotype (W1/w1, i, R, t), suggesting that GmMYBA2 was responsible for the R locus. Molecular and biochemical analysis showed that GmMYBA2 interacted with GmTT8a to directly activate anthocyanin biosynthetic genes. GmMYBA2 and GmMYBR might form a feedback loop to fine-tune seed coat coloration, which was confirmed in transgenic soybeans. Both GmTT8a and GmMYBR that were activated by GmMYBA2 in turn enhanced and obstructed the formation of the GmMYBA2-GmTT8a module, respectively. The results revealed the sophisticated regulatory network underlying the soybean seed coat pigmentation loci and shed light on the understanding of the seed coat coloration and other seed inclusions.