A unified model for interpretable latent embedding of multi-sample, multi-condition single-cell data.

Single-cell analysis across multiple samples and conditions requires quantitative modeling of the interplay between the continuum of cell states and the technical and biological sources of sample-to-sample variability. We introduce GEDI, a generative model that identifies latent space variations in multi-sample, multi-condition single-cell datasets and attributes them to sample-level covariates. GEDI enables cross-sample cell state mapping on par with state-of-the-art integration methods, cluster-free differential gene expression analysis along the continuum of cell states, and machine learning-based prediction of sample characteristics from single-cell data. GEDI can also incorporate gene-level prior knowledge to infer pathway and regulatory network activities in single cells. Finally, GEDI extends all these concepts to previously unexplored modalities that require joint consideration of dual measurements, such as the joint analysis of exon inclusion/exclusion reads to model alternative cassette exon splicing, or spliced/unspliced reads to model the mRNA stability landscapes of single cells.

Circulating Metabolite Abundances Associated With Risks of Bipolar Disorder, Schizophrenia, and Depression: A Mendelian Randomization Study.

BACKGROUND: Preventive measures and treatments for psychiatric disorders are limited. Circulating metabolites are potential candidates for biomarker and therapeutic target identification, given their measurability and essential roles in biological processes. METHODS: Leveraging large-scale genome-wide association studies, we conducted Mendelian randomization analyses to assess the associations between circulating metabolite abundances and the risks of bipolar disorder, schizophrenia, and depression. Genetic instruments were selected for 94 metabolites measured in the Canadian Longitudinal Study on Aging cohort (N = 8299). We repeated Mendelian randomization analyses based on the UK Biobank, INTERVAL, and EPIC (European Prospective Investigation into Cancer)-Norfolk studies. RESULTS: After validating Mendelian randomization assumptions and colocalization evidence, we found that a 1 SD increase in genetically predicted circulating abundances of eicosapentaenoate and docosapentaenoate was associated with odds ratios of 0.72 (95% CI, 0.65-0.79) and 0.63 (95% CI, 0.55-0.72), respectively, for bipolar disorder. Genetically increased Ω-3 unsaturated fatty acids abundance and Ω-3-to-total fatty acids ratio, as well as genetically decreased Ω-6-to-Ω-3 ratio, were negatively associated with the risk of bipolar disorder in the UK Biobank. Genetically increased circulating abundances of 3 N-acetyl-amino acids were associated with an increased risk of schizophrenia with a maximum odds ratio of 1.31 (95% CI, 1.18-1.44) per 1 SD increase. Furthermore, a 1 SD increase in genetically predicted circulating abundance of hypotaurine was associated with an odds ratio of 0.85 (95% CI, 0.78-0.93) for depression. CONCLUSIONS: The biological mechanisms that underlie Ω-3 unsaturated fatty acids, NAT8-catalyzed N-acetyl-amino acids, and hypotaurine warrant exploration to identify new biomarkers and potential therapeutic targets.

Vitamin D, Cognition, and Alzheimer's Disease: Observational and Two-Sample Mendelian Randomization Studies.

Background: Observational studies have found that vitamin D supplementation is associated with improved cognition. Further, recent Mendelian randomization (MR) studies have shown that higher vitamin D levels, 25(OH)D, may protect against Alzheimer's disease. Thus, it is possible that 25(OH)D may protect against Alzheimer's disease by improving cognition. Objective: We assessed this hypothesis, by examining the relationship between 25(OH)D levels and seven cognitive measurements. Methods: To mitigate bias from confounding, we performed two-sample MR analyses. We used instruments from three publications: Manousaki et al. (2020), Sutherland et al. (2022), and the Emerging Risk Factors Collaboration/EPIC-CVD/Vitamin D Studies Collaboration (2021). Results: Our observational studies suggested a protective association between 25(OH)D levels and cognitive measures. An increase in the natural log of 25(OH)D by 1 SD was associated with a higher PACC score (BetaPACC score = 0.06, 95% CI = (0.04-0.08); p = 1.8×10-10). However, in the MR analyses, the estimated effect of 25(OH)D on cognitive measures was null. Specifically, per 1 SD increase in genetically estimated natural log of 25(OH)D, the PACC scores remained unchanged in the overall population, (BetaPACC score = -0.01, 95% CI (-0.06 to 0.03); p = 0.53), and amongst individuals aged over 60 (BetaPACC score = 0.03, 95% CI (-0.028 to 0.08); p = 0.35). Conclusions: In conclusion, our MR study found no clear evidence to support a protective role of increased 25(OH)D concentrations on cognitive performance in European ancestry individuals. However, our study cannot entirely dismiss the potential beneficial effect on PACC for individuals over the age of 60.

A plasma protein-based risk score to predict hip fractures.

As there are effective treatments to reduce hip fractures, identification of patients at high risk of hip fracture is important to inform efficient intervention strategies. To obtain a new tool for hip fracture prediction, we developed a protein-based risk score in the Cardiovascular Health Study using an aptamer-based proteomic platform. The proteomic risk score predicted incident hip fractures and improved hip fracture discrimination in two Trøndelag Health Study validation cohorts using the same aptamer-based platform. When transferred to an antibody-based proteomic platform in a UK Biobank validation cohort, the proteomic risk score was strongly associated with hip fractures (hazard ratio per s.d. increase, 1.64; 95% confidence interval 1.53-1.77). The proteomic risk score, but not available polygenic risk scores for fractures or bone mineral density, improved the C-index beyond the fracture risk assessment tool (FRAX), which integrates information from clinical risk factors (C-index, FRAX 0.735 versus FRAX + proteomic risk score 0.776). The developed proteomic risk score constitutes a new tool for stratifying patients according to hip fracture risk; however, its improvement in hip fracture discrimination is modest and its clinical utility beyond FRAX with information on femoral neck bone mineral density remains to be determined.

SharePro: an accurate and efficient genetic colocalization method accounting for multiple causal signals.

MOTIVATION: Colocalization analysis is commonly used to assess whether two or more traits share the same genetic signals identified in genome-wide association studies (GWAS), and is important for prioritizing targets for functional follow-up of GWAS results. Existing colocalization methods can have suboptimal performance when there are multiple causal variants in one genomic locus. RESULTS: We propose SharePro to extend the COLOC framework for colocalization analysis. SharePro integrates linkage disequilibrium (LD) modeling and colocalization assessment by grouping correlated variants into effect groups. With an efficient variational inference algorithm, posterior colocalization probabilities can be accurately estimated. In simulation studies, SharePro demonstrated increased power with a well-controlled false positive rate at a low computational cost. Compared to existing methods, SharePro provided stronger and more consistent colocalization evidence for known lipid-lowering drug target proteins and their corresponding lipid traits. Through an additional challenging case of the colocalization analysis of the circulating abundance of R-spondin 3 GWAS and estimated bone mineral density GWAS, we demonstrated the utility of SharePro in identifying biologically plausible colocalized signals. AVAILABILITY AND IMPLEMENTATION: SharePro for colocalization analysis is written in Python and openly available at https://github.com/zhwm/SharePro_coloc.

Large-scale circulating proteome association study (CPAS) meta-analysis identifies circulating proteins and pathways predicting incident hip fractures.

Hip fractures are associated with significant disability, high cost, and mortality. However, the exact biological mechanisms underlying susceptibility to hip fractures remain incompletely understood. In an exploratory search of the underlying biology as reflected through the circulating proteome, we performed a comprehensive Circulating Proteome Association Study (CPAS) meta-analysis for incident hip fractures. Analyses included 6430 subjects from two prospective cohort studies (Cardiovascular Health Study and Trøndelag Health Study) with circulating proteomics data (aptamer-based 5 K SomaScan version 4.0 assay; 4979 aptamers). Associations between circulating protein levels and incident hip fractures were estimated for each cohort using age and sex-adjusted Cox regression models. Participants experienced 643 incident hip fractures. Compared with the individual studies, inverse-variance weighted meta-analyses yielded more statistically significant associations, identifying 23 aptamers associated with incident hip fractures (conservative Bonferroni correction 0.05/4979, P < 1.0 × 10-5). The aptamers most strongly associated with hip fracture risk corresponded to two proteins of the growth hormone/insulin growth factor system (GHR and IGFBP2), as well as GDF15 and EGFR. High levels of several inflammation-related proteins (CD14, CXCL12, MMP12, ITIH3) were also associated with increased hip fracture risk. Ingenuity pathway analysis identified reduced LXR/RXR activation and increased acute phase response signaling to be overrepresented among those proteins associated with increased hip fracture risk. These analyses identified several circulating proteins and pathways consistently associated with incident hip fractures. These findings underscore the usefulness of the meta-analytic approach for comprehensive CPAS in a similar manner as has previously been observed for large-scale human genetic studies. Future studies should investigate the underlying biology of these potential novel drug targets.

Hip fractures are associated with significant disability, high cost, and mortality. However, the exact biological mechanisms underlying susceptibility to hip fractures remain incompletely understood. To increase the understanding of the underlying mechanisms, we performed a meta-analysis of the associations between 4860 circulating proteins and risk of fractures using two large cohorts, including 6430 participants with 643 incident hip fractures. We identified 23 proteins/aptamers associated with incident hip fractures. Two proteins of the growth hormone/insulin growth factor system (GHR and IGFBP2), as well as GDF15 and EGFR were most strongly associated with hip fracture risk. High levels of several inflammation-related proteins were also associated with increased hip fracture risk. Pathway analysis identified reduced LXR/RXR activation and increased acute phase response signaling to be overrepresented among those proteins associated with increased hip fracture risk. Future mechanistic studies should investigate the underlying biology of these novel protein biomarkers which may be potential drug targets.

Conditional Knockout of IL-1R1 in Endothelial Cells Attenuates Seizures and Neurodegeneration via Inhibiting Neuroinflammation Mediated by Nrf2/HO-1/NLRP3 Signaling in Status Epilepticus Model.

Studies on the bench and at bedside have demonstrated that the process of epileptogenesis is involved in neuroinflammatory responses. As the receptor of proinflammatory cytokine IL-1ß, IL-1ß type 1 receptor (IL-1R1) is reported to express abundantly in the endothelial cells in epileptic brains, which is deemed to be implicated in the epileptogenic process. However, whether and how endothelial IL-1R1 modulates neuroinflammatory responses in the pathological process of epileptic seizures and/or status epilepticus (SE) remains obscure. Here, we indicated endothelial IL-1R1 is involved in neuroinflammation, facilitating epilepsy progress via Nrf2/HO-1/NLRP3. In vitro, we observed upregulation of inflammatory cytokines in co-culture model under IL-1ß challenge, as well as in BV2 cells after stimulation with conditional medium (CM) from IL-1ß-stimulated bEnd.3 cells. In vivo, mice with conditional knockout of endothelial IL-1R1 (IL-1R1-CKO) were generated by hybrid IL-1R1flox/flox mice with Tek-Cre mice. IL-1R1-CKO reduced seizure susceptibility in kainic acid (KA)-induced SE model. In addition, IL-1R1-CKO KA mice exhibited lessened hippocampal neuroinflammation, mitigated neuronal damage, and decreased abnormal neurogenesis. In cognitive behavioral tests, IL-1R1-CKO KA mice presented improvement in learning and memory. Furthermore, we also indicated blockage of endothelial IL-1R1 downregulated the expressions of Nrf2/HO-1/NLRP3 pathway-related proteins. Nrf2-siRNA reversed the downregulation of HO-1, NLRP3, caspase-1, and IL-1ß. These results demonstrated CKO of endothelial IL-1R1 reduces seizure susceptibility and attenuates SE-related neurobehavioral damage by suppressing hippocampal neuroinflammation via Nrf2/HO-1/NLRP3.

The roles of 6K protein on Getah virus replication and pathogenicity.

Alphavirus is a type of arbovirus that can infect both humans and animals. The amino acid sequence of the 6K protein, being one of the structural proteins of the alphavirus, is not conserved. Deletion of this protein will result in varying effects on different alphaviruses. Our study focuses on the function of the Getah virus (GETV) 6K protein in infected cells and mice. We successfully constructed infectious clone plasmids and created resulting viruses (rGETV and rGETV-Δ6K). Our comprehensive microscopic analysis revealed that the 6 K protein mainly stays in the endoplasmic reticulum. In addition, rGETV-Δ6K has lower thermal stability and sensitivity to temperature than GETV. Although the deletion of the 6K protein does not reduce virion production in ST cells, it affects the release of virions from host cells by inhibiting the process of E2 protein transportation to the plasma membrane. Subsequent in vivo testing demonstrated that neonatal mice infected with rGETV-Δ6K had a lower virus content, less significant pathological changes in tissue slices, and milder disease than those infected with the wild-type virus. Our results indicate that the 6K protein effectively reduces the viral titer by influencing the release of viral particles. Furthermore, the 6K protein play a role in the clinical manifestation of GETV disease.

An atlas of genetic determinants of forearm fracture.

Osteoporotic fracture is among the most common and costly of diseases. While reasonably heritable, its genetic determinants have remained elusive. Forearm fractures are the most common clinically recognized osteoporotic fractures with a relatively high heritability. To establish an atlas of the genetic determinants of forearm fractures, we performed genome-wide association analyses including 100,026 forearm fracture cases. We identified 43 loci, including 26 new fracture loci. Although most fracture loci associated with bone mineral density, we also identified loci that primarily regulate bone quality parameters. Functional studies of one such locus, at TAC4, revealed that Tac4-/- mice have reduced mechanical bone strength. The strongest forearm fracture signal, at WNT16, displayed remarkable bone-site-specificity with no association with hip fractures. Tall stature and low body mass index were identified as new causal risk factors for fractures. The insights from this atlas may improve fracture prediction and enable therapeutic development to prevent fractures.

HLA allele-calling using multi-ancestry whole-exome sequencing from the UK Biobank identifies 129 novel associations in 11 autoimmune diseases.

The human leukocyte antigen (HLA) region on chromosome 6 is strongly associated with many immune-mediated and infection-related diseases. Due to its highly polymorphic nature and complex linkage disequilibrium patterns, traditional genetic association studies of single nucleotide polymorphisms do not perform well in this region. Instead, the field has adopted the assessment of the association of HLA alleles (i.e., entire HLA gene haplotypes) with disease. Often based on genotyping arrays, these association studies impute HLA alleles, decreasing accuracy and thus statistical power for rare alleles and in non-European ancestries. Here, we use whole-exome sequencing (WES) from 454,824 UK Biobank (UKB) participants to directly call HLA alleles using the HLA-HD algorithm. We show this method is more accurate than imputing HLA alleles and harness the improved statistical power to identify 360 associations for 11 auto-immune phenotypes (at least 129 likely novel), leading to better insights into the specific coding polymorphisms that underlie these diseases. We show that HLA alleles with synonymous variants, often overlooked in HLA studies, can significantly influence these phenotypes. Lastly, we show that HLA sequencing may improve polygenic risk scores accuracy across ancestries. These findings allow better characterization of the role of the HLA region in human disease.

Beyond apples and pears: sex-specific genetics of body fat percentage.

Introduction: Biological sex influences both overall adiposity and fat distribution. Further, testosterone and sex hormone binding globulin (SHBG) influence adiposity and metabolic function, with differential effects of testosterone in men and women. Here, we aimed to perform sex-stratified genome-wide association studies (GWAS) of body fat percentage (BFPAdj) (adjusting for testosterone and sex hormone binding globulin (SHBG)) to increase statistical power. Methods: GWAS were performed in white British individuals from the UK Biobank (157,937 males and 154,337 females). To avoid collider bias, loci associated with SHBG or testosterone were excluded. We investigated association of BFPAdj loci with high density cholesterol (HDL), triglyceride (TG), type 2 diabetes (T2D), coronary artery disease (CAD), and MRI-derived abdominal subcutaneous adipose tissue (ASAT), visceral adipose tissue (VAT) and gluteofemoral adipose tissue (GFAT) using publicly available data from large GWAS. We also performed 2-sample Mendelian Randomization (MR) using identified BFPAdj variants as instruments to investigate causal effect of BFPAdj on HDL, TG, T2D and CAD in males and females separately. Results: We identified 195 and 174 loci explaining 3.35% and 2.60% of the variation in BFPAdj in males and females, respectively at genome-wide significance (GWS, p<5x10-8). Although the direction of effect at these loci was generally concordant in males and females, only 38 loci were common to both sexes at GWS. Seven loci in males and ten loci in females have not been associated with any adiposity/cardiometabolic traits previously. BFPAdj loci generally did not associate with cardiometabolic traits; several had paradoxically beneficial cardiometabolic effects with favourable fat distribution. MR analyses did not find convincing supportive evidence that increased BFPAdj has deleterious cardiometabolic effects in either sex with highly significant heterogeneity. Conclusions: There was limited genetic overlap between BFPAdj in males and females at GWS. BFPAdj loci generally did not have adverse cardiometabolic effects which may reflect the effects of favourable fat distribution and cardiometabolic risk modulation by testosterone and SHBG.

Identifying Rare Genetic Determinants for Improved Polygenic Risk Prediction of Bone Mineral Density and Fracture Risk.

Osteoporosis and fractures severely impact the elderly population. Polygenic risk scores for bone mineral density have demonstrated potential clinical utility. However, the value of rare genetic determinants in risk prediction has not been assessed. With whole-exome sequencing data from 436,824 UK Biobank participants, we assigned White British ancestry individuals into a training data set (n = 317,434) and a test data set (n = 74,825). In the training data set, we developed a common variant-based polygenic risk score for heel ultrasound speed of sound (SOS). Next, we performed burden testing to identify genes harboring rare determinants of bone mineral density, targeting influential rare variants with predicted high deleteriousness. We constructed a genetic risk score, called ggSOS, to incorporate influential rare variants in significant gene burden masks into the common variant-based polygenic risk score. We assessed the predictive performance of ggSOS in the White British test data set, as well as in populations of non-White British European (n = 18,885), African (n = 7165), East Asian (n = 2236), South Asian (n = 9829), and other admixed (n = 1481) ancestries. Twelve genes in pivotal regulatory pathways of bone homeostasis harbored influential rare variants associated with SOS (p < 5.5 × 10-7 ), including AHNAK, BMP5, CYP19A1, FAM20A, FBXW5, KDM5B, KREMEN1, LGR4, LRP5, SMAD6, SOST, and WNT1. Among 4013 (5.4%) individuals in the test data set carrying these variants, a one standard deviation decrease in ggSOS was associated with 1.35-fold (95% confidence interval [CI] 1.16-1.57) increased hazard of major osteoporotic fracture. However, compared with a common variant-based polygenic risk score (C-index = 0.641), ggSOS had only marginally improved prediction accuracy in identifying at-risk individuals (C-index = 0.644), with overlapping confidence intervals. Similarly, ggSOS did not demonstrate substantially improved predictive performance in non-European ancestry populations. In summary, modeling the effects of rare genetic determinants may assist polygenic prediction of fracture risk among carriers of influential rare variants. Nonetheless, improved clinical utility is not guaranteed for population-level risk screening. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Development of risk prediction models for depression combining genetic and early life risk factors.

Background: Both genetic and early life risk factors play important roles in the pathogenesis and progression of adult depression. However, the interplay between these risk factors and their added value to risk prediction models have not been fully elucidated. Methods: Leveraging a meta-analysis of major depressive disorder genome-wide association studies (N = 45,591 cases and 97,674 controls), we developed and optimized a polygenic risk score for depression using LDpred in a model selection dataset from the UK Biobank (N = 130,092 European ancestry individuals). In a UK Biobank test dataset (N = 278,730 European ancestry individuals), we tested whether the polygenic risk score and early life risk factors were associated with each other and compared their associations with depression phenotypes. Finally, we conducted joint predictive modeling to combine this polygenic risk score with early life risk factors by stepwise regression, and assessed the model performance in identifying individuals at high risk of depression. Results: In the UK Biobank test dataset, the polygenic risk score for depression was moderately associated with multiple early life risk factors. For instance, a one standard deviation increase in the polygenic risk score was associated with 1.16-fold increased odds of frequent domestic violence (95% CI: 1.14-1.19) and 1.09-fold increased odds of not having access to medical care as a child (95% CI: 1.05-1.14). However, the polygenic risk score was more strongly associated with depression phenotypes than most early life risk factors. A joint predictive model integrating the polygenic risk score, early life risk factors, age and sex achieved an AUROC of 0.6766 for predicting strictly defined major depressive disorder, while a model without the polygenic risk score and a model without any early life risk factors had an AUROC of 0.6593 and 0.6318, respectively. Conclusion: We have developed a polygenic risk score to partly capture the genetic liability to depression. Although genetic and early life risk factors can be correlated, joint predictive models improved risk stratification despite limited improvement in magnitude, and may be explored as tools to better identify individuals at high risk of depression.

Colocalization of expression transcripts with COVID-19 outcomes is rare across cell states, cell types and organs.

Identifying causal genes at GWAS loci can help pinpoint targets for therapeutic interventions. Expression studies can disentangle such loci but signals from expression quantitative trait loci (eQTLs) often fail to colocalize-which means that the genetic control of measured expression is not shared with the genetic control of disease risk. This may be because gene expression is measured in the wrong cell type, physiological state, or organ. We tested whether Mendelian randomization (MR) could identify genes at loci influencing COVID-19 outcomes and whether the colocalization of genetic control of expression and COVID-19 outcomes was influenced by cell type, cell stimulation, and organ. We conducted MR of cis-eQTLs from single cell (scRNA-seq) and bulk RNA sequencing. We then tested variables that could influence colocalization, including cell type, cell stimulation, RNA sequencing modality, organ, symptoms of COVID-19, and SARS-CoV-2 status among individuals with symptoms of COVID-19. The outcomes used to test colocalization were COVID-19 severity and susceptibility as assessed in the Host Genetics Initiative release 7. Most transcripts identified using MR did not colocalize when tested across cell types, cell state and in different organs. Most that did colocalize likely represented false positives due to linkage disequilibrium. In general, colocalization was highly variable and at times inconsistent for the same transcript across cell type, cell stimulation and organ. While we identified factors that influenced colocalization for select transcripts, identifying 33 that mediate COVID-19 outcomes, our study suggests that colocalization of expression with COVID-19 outcomes is partially due to noisy signals even after following quality control and sensitivity testing. These findings illustrate the present difficulty of linking expression transcripts to disease outcomes and the need for skepticism when observing eQTL MR results, even accounting for cell types, stimulation state and different organs.

eXclusionarY: 10 years later, where are the sex chromosomes in GWASs?

10 years ago, a detailed analysis showed that only 33% of genome-wide association study (GWAS) results included the X chromosome. Multiple recommendations were made to combat such exclusion. Here, we re-surveyed the research landscape to determine whether these earlier recommendations had been translated. Unfortunately, among the genome-wide summary statistics reported in 2021 in the NHGRI-EBI GWAS Catalog, only 25% provided results for the X chromosome and 3% for the Y chromosome, suggesting that the exclusion phenomenon not only persists but has also expanded into an exclusionary problem. Normalizing by physical length of the chromosome, the average number of studies published through November 2022 with genome-wide-significant findings on the X chromosome is ∼1 study/Mb. By contrast, it ranges from ∼6 to ∼16 studies/Mb for chromosomes 4 and 19, respectively. Compared with the autosomal growth rate of ∼0.086 studies/Mb/year over the last decade, studies of the X chromosome grew at less than one-seventh that rate, only ∼0.012 studies/Mb/year. Among the studies that reported significant associations on the X chromosome, we noted extreme heterogeneities in data analysis and reporting of results, suggesting the need for clear guidelines. Unsurprisingly, among the 430 scores sampled from the PolyGenic Score Catalog, 0% contained weights for sex chromosomal SNPs. To overcome the dearth of sex chromosome analyses, we provide five sets of recommendations and future directions. Finally, until the sex chromosomes are included in a whole-genome study, instead of GWASs, we propose such studies would more properly be referred to as "AWASs," meaning "autosome-wide scans."

Dose-dependent Association of Alcohol Consumption With Obesity and Type 2 Diabetes: Mendelian Randomization Analyses.

CONTEXT: Effects of modest alcohol consumption remain controversial. Mendelian randomization (MR) can help to mitigate biases due to confounding and reverse causation in observational studies, and evaluate the potential causal role of alcohol consumption. OBJECTIVE: This work aimed to evaluate dose-dependent effect of alcohol consumption on obesity and type 2 diabetes. METHODS: Assessing 408 540 participants of European ancestry in the UK Biobank, we first tested the association between self-reported alcohol intake frequency and 10 anthropometric measurements, obesity, and type 2 diabetes. We then conducted MR analyses both in the overall population and in subpopulations stratified by alcohol intake frequency. RESULTS: Among individuals having more than 14 drinks per week, a 1-drink-per-week increase in genetically predicted alcohol intake frequency was associated with a 0.36-kg increase in fat mass (SD = 0.03 kg), a 1.08-fold increased odds of obesity (95% CI, 1.06-1.10), and a 1.10-fold increased odds of type 2 diabetes (95% CI, 1.06-1.13). These associations were stronger in women than in men. Furthermore, no evidence was found supporting the association between genetically increased alcohol intake frequency and improved health outcomes among individuals having 7 or fewer drinks per week, as MR estimates largely overlapped with the null. These results withstood multiple sensitivity analyses assessing the validity of MR assumptions. CONCLUSION: As opposed to observational associations, MR results suggest there may not be protective effects of modest alcohol consumption on obesity traits and type 2 diabetes. Heavy alcohol consumption could lead to increased measures of obesity as well as increased risk of type 2 diabetes.

Novel insights into systemic sclerosis using a sensitive computational method to analyze whole-genome bisulfite sequencing data.

BACKGROUND: Abnormal DNA methylation is thought to contribute to the onset and progression of systemic sclerosis. Currently, the most comprehensive assay for profiling DNA methylation is whole-genome bisulfite sequencing (WGBS), but its precision depends on read depth and it may be subject to sequencing errors. SOMNiBUS, a method for regional analysis, attempts to overcome some of these limitations. Using SOMNiBUS, we re-analyzed WGBS data previously analyzed using bumphunter, an approach that initially fits single CpG associations, to contrast DNA methylation estimates by both methods. METHODS: Purified CD4+ T lymphocytes of 9 SSc and 4 control females were sequenced using WGBS. We separated the resulting sequencing data into regions with dense CpG data, and differentially methylated regions (DMRs) were inferred with the SOMNiBUS region-level test, adjusted for age. Pathway enrichment analysis was performed with ingenuity pathway analysis (IPA). We compared the results obtained by SOMNiBUS and bumphunter. RESULTS: Of 8268 CpG regions of ≥ 60 CpGs eligible for analysis with SOMNiBUS, we identified 131 DMRs and 125 differentially methylated genes (DMGs; p-values less than Bonferroni-corrected threshold of 6.05-06 controlling family-wise error rate at 0.05; 1.6% of the regions). In comparison, bumphunter identified 821,929 CpG regions, 599 DMRs (of which none had ≥ 60 CpGs) and 340 DMGs (q-value of 0.05; 0.04% of all regions). The top ranked gene identified by SOMNiBUS was FLT4, a lymphangiogenic orchestrator, and the top ranked gene on chromosome X was CHST7, known to catalyze the sulfation of glycosaminoglycans in the extracellular matrix. The top networks identified by IPA included connective tissue disorders. CONCLUSIONS: SOMNiBUS is a complementary method of analyzing WGBS data that enhances biological insights into SSc and provides novel avenues of investigation into its pathogenesis.

Predicting ExWAS findings from GWAS data: a shorter path to causal genes.

GWAS has identified thousands of loci associated with disease, yet the causal genes within these loci remain largely unknown. Identifying these causal genes would enable deeper understanding of the disease and assist in genetics-based drug development. Exome-wide association studies (ExWAS) are more expensive but can pinpoint causal genes offering high-yield drug targets, yet suffer from a high false-negative rate. Several algorithms have been developed to prioritize genes at GWAS loci, such as the Effector Index (Ei), Locus-2-Gene (L2G), Polygenic Prioritization score (PoPs), and Activity-by-Contact score (ABC) and it is not known if these algorithms can predict ExWAS findings from GWAS data. However, if this were the case, thousands of associated GWAS loci could potentially be resolved to causal genes. Here, we quantified the performance of these algorithms by evaluating their ability to identify ExWAS significant genes for nine traits. We found that Ei, L2G, and PoPs can identify ExWAS significant genes with high areas under the precision recall curve (Ei: 0.52, L2G: 0.37, PoPs: 0.18, ABC: 0.14). Furthermore, we found that for every unit increase in the normalized scores, there was an associated 1.3-4.6-fold increase in the odds of a gene reaching exome-wide significance (Ei: 4.6, L2G: 2.5, PoPs: 2.1, ABC: 1.3). Overall, we found that Ei, L2G, and PoPs can anticipate ExWAS findings from widely available GWAS results. These techniques are therefore promising when well-powered ExWAS data are not readily available and can be used to anticipate ExWAS findings, allowing for prioritization of genes at GWAS loci.

Transcriptional reprogramming of skeletal muscle stem cells by the niche environment.

Adult stem cells are indispensable for tissue regeneration, but their function declines with age. The niche environment in which the stem cells reside plays a critical role in their function. However, quantification of the niche effect on stem cell function is lacking. Using muscle stem cells (MuSC) as a model, we show that aging leads to a significant transcriptomic shift in their subpopulations accompanied by locus-specific gain and loss of chromatin accessibility and DNA methylation. By combining in vivo MuSC transplantation and computational methods, we show that the expression of approximately half of all age-altered genes in MuSCs from aged male mice can be restored by exposure to a young niche environment. While there is a correlation between gene reversibility and epigenetic alterations, restoration of gene expression occurs primarily at the level of transcription. The stem cell niche environment therefore represents an important therapeutic target to enhance tissue regeneration in aging.

Proteome-wide Mendelian randomization implicates nephronectin as an actionable mediator of the effect of obesity on COVID-19 severity.

Obesity is a major risk factor for Coronavirus disease (COVID-19) severity; however, the mechanisms underlying this relationship are not fully understood. As obesity influences the plasma proteome, we sought to identify circulating proteins mediating the effects of obesity on COVID-19 severity in humans. Here, we screened 4,907 plasma proteins to identify proteins influenced by body mass index using Mendelian randomization. This yielded 1,216 proteins, whose effect on COVID-19 severity was assessed, again using Mendelian randomization. We found that an s.d. increase in nephronectin (NPNT) was associated with increased odds of critically ill COVID-19 (OR = 1.71, P = 1.63 × 10-10). The effect was driven by an NPNT splice isoform. Mediation analyses supported NPNT as a mediator. In single-cell RNA-sequencing, NPNT was expressed in alveolar cells and fibroblasts of the lung in individuals who died of COVID-19. Finally, decreasing body fat mass and increasing fat-free mass were found to lower NPNT levels. These findings provide actionable insights into how obesity influences COVID-19 severity.