Background/Objectives: Nutritional deficiencies have been proposed as possible etiological causes for autoimmune diseases, among which type 1 diabetes (T1D). Vitamin K (VK) has potentially positive effects on type 2 diabetes, but its role on T1D in humans remains largely unknown. We aimed to examine the presence of a causal association between VK and T1D using a Mendelian randomization (MR) approach. Methods: Genetic variants from a genome-wide association study (GWAS) for VK (N = 2138 Europeans) were used as instruments in our two-sample MR study to investigate whether circulating VK levels are causally associated with the risk of T1D in a large European T1D GWAS cohort (18,942 cases/520,580 controls). Through a multivariable MR (MVMR), the effects of both VK and specific gut microbiota on T1D were investigated given that the gut microbiome synthesizes VK. Results: We found that changes in levels of circulating VK did not affect T1D risk in our univariate two-sample MR, but this study had limited power to detect small effects of VK (OR for T1D of less than 0.8). However, our MVMR indicated a suggestive association of VK with the risk of T1D adjusting for two different gut microbiome populations. Conclusions: In conclusion, VK levels are unlikely to significantly affect the risk of T1D, but small effects cannot be excluded, and the role of gut microbiome in this association should be further investigated.

Cytochromes P450 (CYP450) are hemoproteins generally involved in the detoxification of the body of xenobiotic molecules. They participate in the metabolism of many drugs and genetic polymorphisms in humans have been found to impact drug responses and metabolic functions. In this study, we investigate the genetic diversity of CYP450 genes. We found that two clusters, CYP3A and CYP4F, are notably differentiated across human populations with evidence for selective pressures acting on both clusters: we found signals of recent positive selection in CYP3A and CYP4F genes and signals of balancing selection in CYP4F genes. Furthermore, an extensive amount of unusual linkage disequilibrium is detected in this latter cluster, indicating co-evolution signatures among CYP4F genes. Several of the selective signals uncovered co-localize with expression quantitative trait loci (eQTL), which could suggest epistasis acting on co-regulation in these gene families. In particular, we detected a potential co-regulation event between CYP3A5 and CYP3A43, a gene whose function remains poorly characterized. We further identified a causal relationship between CYP3A5 expression and reticulocyte count through Mendelian randomization analyses, potentially involving a regulatory region displaying a selective signal specific to African populations. Our findings linking natural selection and gene expression in CYP3A and CYP4F subfamilies are of importance in understanding population differences in metabolism of nutrients and drugs.

Abstract Background The role of metabolism in the variation of age at menarche (AAM) and age at natural menopause (ANM) in the female population is not entirely known. We aimed to investigate the causal role of circulating metabolites in AAM and ANM using Mendelian randomization (MR). Methods We combined MR with genetic colocalization to investigate potential causal associations between 658 metabolites and AAM and between 684 metabolites and ANM. We extracted genetic instruments for our exposures from four genome-wide association studies (GWAS) on circulating metabolites and queried the effects of these variants on the outcomes in two large GWAS from the ReproGen consortium. Additionally, we assessed the mediating role of the body mass index (BMI) in these associations, identified metabolic pathways implicated in AAM and ANM, and sought validation for selected metabolites in the Avon Longitudinal Study of Parents and Children (ALSPAC). Results Our analysis identified 10 candidate metabolites for AAM, but none of them colocalized with AAM. For ANM, 76 metabolites were prioritized (FDR-adjusted MR P-value ≤ 0.05), with 17 colocalizing, primarily in the glycerophosphocholines class, including the omega-3 fatty acid and phosphatidylcholine (PC) categories. Pathway analyses and validation in ALSPAC mothers also highlighted the role of omega and polyunsaturated fatty acids levels in delaying age at menopause. Conclusions Our study suggests that metabolites from the glycerophosphocholine and fatty acid families play a causal role in the timing of both menarche and menopause. This underscores the significance of specific metabolic pathways in the biology of female reproductive longevity.

ABSTRACT Adolescent idiopathic scoliosis (AIS) is the most common form of pediatric musculoskeletal disorder. Observational studies have pointed to several risk factors for AIS, but almost no evidence exists to support their causal association with AIS. Here, we applied Mendelian randomization (MR), known to limit bias from confounding and reverse causation, to investigate causal associations between body composition and puberty‐related exposures and AIS risk in Europeans and Asians. For our two‐sample MR studies, we used single nucleotide polymorphisms (SNPs) associated with body mass index (BMI), waist‐hip ratio, lean mass, childhood obesity, bone mineral density (BMD), 25‐hydroxyvitamin D (25OHD), age at menarche, and pubertal growth in large European genome‐wide association studies (GWAS), and with adult osteoporosis risk and age of menarche in Biobank Japan. We extracted estimates of the aforementioned SNPs on AIS risk from the European or Asian subsets of the largest multiancestry AIS GWAS ( N = 7956 cases/88,459 controls). The results of our inverse variance‐weighted (IVW) MR estimates suggest no causal association between the aforementioned risk factors and risk of AIS. Pleiotropy‐sensitive MR methods yielded similar results. However, restricting our analysis to European females with AIS, we observed a causal association between estimated BMD and the risk of AIS (IVW odds ratio for AIS = 0.1, 95% confidence interval 0.01 to 0.7, p = 0.02 per SD increase in estimated BMD), but this association was no longer significant after adjusting for BMI, body fat mass, and 25OHD and remained significant after adjusting for age at menarche in multivariable MR. In conclusion, we demonstrated a protective causal effect of BMD on AIS risk in females of European ancestry, but this effect was modified by BMI, body fat mass, and 25OHD levels. Future MR studies using larger AIS GWAS are needed to investigate small effects of the aforementioned exposures on AIS. © 2023 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Genetic variants in drug targets can be used to predict the effect of drugs. Here, we extend this principle to assess how sex and body mass index may modify the effect of a genetically predicted lower CETP levels on biomarkers and cardiovascular outcomes. We found sex and BMI to be modifiers of the association between genetically predicted lower CETP and lipid biomarkers in UK Biobank participants. Female sex and lower BMI were associated with higher HDL-cholesterol and lower LDL-cholesterol for a same genetically predicted reduction in CETP concentration. We found that sex also modulated the effect of genetically lower CETP on cholesterol efflux capacity in samples from the Montreal Heart Institute Biobank. However, these modifying effects did not extend to sex-differences in cardiovascular outcomes in our data. Our results provide insight on the clinical effects of CETP inhibitors in the presence of effect modification based on observational genetic data. The approach can support precision medicine applications and help assess the external validity of clinical trials.

The genome of the Severe Acute Respiratory Syndrome coronavirus 2 (SARS-CoV-2), the pathogen that causes coronavirus disease 2019 (COVID-19), has been sequenced at an unprecedented scale leading to a tremendous amount of viral genome sequencing data. To assist in tracing infection pathways and design preventive strategies, a deep understanding of the viral genetic diversity landscape is needed. We present here a set of genomic surveillance tools from population genetics which can be used to better understand the evolution of this virus in humans. To illustrate the utility of this toolbox, we detail an in depth analysis of the genetic diversity of SARS-CoV-2 in first year of the COVID-19 pandemic. We analyzed 329,854 high-quality consensus sequences published in the GISAID database during the pre-vaccination phase. We demonstrate that, compared to standard phylogenetic approaches, haplotype networks can be computed efficiently on much larger datasets. This approach enables real-time lineage identification, a clear description of the relationship between variants of concern, and efficient detection of recurrent mutations. Furthermore, time series change of Tajima's D by haplotype provides a powerful metric of lineage expansion. Finally, principal component analysis (PCA) highlights key steps in variant emergence and facilitates the visualization of genomic variation in the context of SARS-CoV-2 diversity. The computational framework presented here is simple to implement and insightful for real-time genomic surveillance of SARS-CoV-2 and could be applied to any pathogen that threatens the health of populations of humans and other organisms.

Pharmacogenomic studies have revealed associations between rs1967309 in the adenylyl cyclase type 9 (ADCY9) gene and clinical responses to the cholesteryl ester transfer protein (CETP) modulator dalcetrapib, however, the mechanism behind this interaction is still unknown. Here, we characterized selective signals at the locus associated with the pharmacogenomic response in human populations and we show that rs1967309 region exhibits signatures of positive selection in several human populations. Furthermore, we identified a variant in CETP, rs158477, which is in long-range linkage disequilibrium with rs1967309 in the Peruvian population. The signal is mainly seen in males, a sex-specific result that is replicated in the LIMAA cohort of over 3400 Peruvians. Analyses of RNA-seq data further suggest an epistatic interaction on CETP expression levels between the two SNPs in multiple tissues, which also differs between males and females. We also detected interaction effects of the two SNPs with sex on cardiovascular phenotypes in the UK Biobank, in line with the sex-specific genotype associations found in Peruvians at these loci. We propose that ADCY9 and CETP coevolved during recent human evolution due to sex-specific selection, which points toward a biological link between dalcetrapib's pharmacogene ADCY9 and its therapeutic target CETP.

The genome of the Severe Acute Respiratory Syndrome coronavirus 2 (SARS-CoV-2), the pathogen that causes coronavirus disease 2019 (COVID-19), has been sequenced at an unprecedented scale, leading to a tremendous amount of viral genome sequencing data. To understand the evolution of this virus in humans, and to assist in tracing infection pathways and designing preventive strategies, we present a set of computational tools that span phylogenomics, population genetics and machine learning approaches. To illustrate the utility of this toolbox, we detail an in depth analysis of the genetic diversity of SARS-CoV-2 in first year of the COVID-19 pandemic, using 329,854 high-quality consensus sequences published in the GISAID database during the pre-vaccination phase. We demonstrate that, compared to standard phylogenetic approaches, haplotype networks can be computed efficiently on much larger datasets, enabling real-time analyses. Furthermore, time series change of Tajima’s D provides a powerful metric of population expansion. Unsupervised learning techniques further highlight key steps in variant detection and facilitate the study of the role of this genomic variation in the context of SARS-CoV-2 infection, with Multiscale PHATE methodology identifying fine-scale structure in the SARS-CoV-2 genetic data that underlies the emergence of key lineages. The computational framework presented here is useful for real-time genomic surveillance of SARS-CoV-2 and could be applied to any pathogen that threatens the health of worldwide populations of humans and other organisms.

Rothmund-Thomson syndrome (RTS) is an autosomal-recessive disorder characterized by poikiloderma, sparse hair, short stature, and skeletal anomalies. Type 2 RTS, which is defined by the presence of bi-allelic mutations in RECQL4, is characterized by increased cancer susceptibility and skeletal anomalies, whereas the genetic basis of RTS type 1, which is associated with juvenile cataracts, is unknown. We studied ten individuals, from seven families, who had RTS type 1 and identified a deep intronic splicing mutation of the ANAPC1 gene, a component of the anaphase-promoting complex/cyclosome (APC/C), in all affected individuals, either in the homozygous state or in trans with another mutation. Fibroblast studies showed that the intronic mutation causes the activation of a 95 bp pseudoexon, leading to mRNAs with premature termination codons and nonsense-mediated decay, decreased ANAPC1 protein levels, and prolongation of interphase. Interestingly, mice that were heterozygous for a knockout mutation have an increased incidence of cataracts. Our results demonstrate that deficiency in the APC/C is a cause of RTS type 1 and suggest a possible link between the APC/C and RECQL4 helicase because both proteins are involved in DNA repair and replication.