Recent changes in the mutational dynamics of the SARS-CoV-2 main protease substantiate the danger of emerging resistance to antiviral drugs

IntroductionThe current coronavirus pandemic is being combated worldwide by nontherapeutic measures and massive vaccination programs. Nevertheless, therapeutic options such as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main-protease (Mpro) inhibitors are essential due to the ongoing evolution toward escape from natural or induced immunity. While antiviral strategies are vulnerable to the effects of viral mutation, the relatively conserved Mpro makes an attractive drug target: Nirmatrelvir, an antiviral targeting its active site, has been authorized for conditional or emergency use in several countries since December 2021, and a number of other inhibitors are under clinical evaluation. We analyzed recent SARS-CoV-2 genomic data, since early detection of potential resistances supports a timely counteraction in drug development and deployment, and discovered accelerated mutational dynamics of Mpro since early December 2021.MethodsWe performed a comparative analysis of 10.5 million SARS-CoV-2 genome sequences available by June 2022 at GISAID to the NCBI reference genome sequence NC_045512.2. Amino-acid exchanges within high-quality regions in 69,878 unique Mpro sequences were identified and time- and in-depth sequence analyses including a structural representation of mutational dynamics were performed using in-house software.ResultsThe analysis showed a significant recent event of mutational dynamics in Mpro. We report a remarkable increase in mutational variability in an eight-residue long consecutive region (R188-G195) near the active site since December 2021.DiscussionThe increased mutational variability in close proximity to an antiviral-drug binding site as described herein may suggest the onset of the development of antiviral resistance. This emerging diversity urgently needs to be further monitored and considered in ongoing drug development and lead optimization

Whole genome sequencing refines stratification and therapy of patients with clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is the most common form of kidney cancer, but a comprehensive description of its genomic landscape is lacking. We report the whole genome sequencing of 778 ccRCC patients enrolled in the 100,000 Genomes Project, providing for a detailed description of the somatic mutational landscape of ccRCC. We identify candidate driver genes, which as well as emphasising the major role of epigenetic regulation in ccRCC highlight additional biological pathways extending opportunities for therapeutic interventions. Genomic characterisation identified patients with divergent clinical outcome; higher number of structural copy number alterations associated with poorer prognosis, whereas VHL mutations were independently associated with a better prognosis. The observations that higher T-cell infiltration is associated with better overall survival and that genetically predicted immune evasion is not common supports the rationale for immunotherapy. These findings should inform personalised surveillance and treatment strategies for ccRCC patients

Neuroimaging-Based Classification of PTSD Using Data-Driven Computational Approaches:A Multisite Big Data Study from the ENIGMA-PGC PTSD Consortium

BACKGROUND: Recent advances in data-driven computational approaches have been helpful in devising tools to objectively diagnose psychiatric disorders. However, current machine learning studies limited to small homogeneous samples, different methodologies, and different imaging collection protocols, limit the ability to directly compare and generalize their results. Here we aimed to classify individuals with PTSD versus controls and assess the generalizability using a large heterogeneous brain datasets from the ENIGMA-PGC PTSD Working group.METHODS: We analyzed brain MRI data from 3,477 structural-MRI; 2,495 resting state-fMRI; and 1,952 diffusion-MRI. First, we identified the brain features that best distinguish individuals with PTSD from controls using traditional machine learning methods. Second, we assessed the utility of the denoising variational autoencoder (DVAE) and evaluated its classification performance. Third, we assessed the generalizability and reproducibility of both models using leave-one-site-out cross-validation procedure for each modality.RESULTS: We found lower performance in classifying PTSD vs. controls with data from over 20 sites (60% test AUC for s-MRI, 59% for rs-fMRI and 56% for d-MRI), as compared to other studies run on single-site data. The performance increased when classifying PTSD from HC without trauma history in each modality (75% AUC). The classification performance remained intact when applying the DVAE framework, which reduced the number of features. Finally, we found that the DVAE framework achieved better generalization to unseen datasets compared with the traditional machine learning frameworks, albeit performance was slightly above chance.CONCLUSION: These results have the potential to provide a baseline classification performance for PTSD when using large scale neuroimaging datasets. Our findings show that the control group used can heavily affect classification performance. The DVAE framework provided better generalizability for the multi-site data. This may be more significant in clinical practice since the neuroimaging-based diagnostic DVAE classification models are much less site-specific, rendering them more generalizable.</p

Regularity Properties and Pathologies of Position-Space Renormalization-Group Transformations

We reconsider the conceptual foundations of the renormalization-group (RG) formalism, and prove some rigorous theorems on the regularity properties and possible pathologies of the RG map. Regarding regularity, we show that the RG map, defined on a suitable space of interactions (= formal Hamiltonians), is always single-valued and Lipschitz continuous on its domain of definition. This rules out a recently proposed scenario for the RG description of first-order phase transitions. On the pathological side, we make rigorous some arguments of Griffiths, Pearce and Israel, and prove in several cases that the renormalized measure is not a Gibbs measure for any reasonable interaction. This means that the RG map is ill-defined, and that the conventional RG description of first-order phase transitions is not universally valid. For decimation or Kadanoff transformations applied to the Ising model in dimension $d \ge 3$, these pathologies occur in a full neighborhood $\{ \beta > \beta_0 ,\, |h| < \epsilon(\beta) \}$ of the low-temperature part of the first-order phase-transition surface. For block-averaging transformations applied to the Ising model in dimension $d \ge 2$, the pathologies occur at low temperatures for arbitrary magnetic-field strength. Pathologies may also occur in the critical region for Ising models in dimension $d \ge 4$. We discuss in detail the distinction between Gibbsian and non-Gibbsian measures, and give a rather complete catalogue of the known examples. Finally, we discuss the heuristic and numerical evidence on RG pathologies in the light of our rigorous theorems.Comment: 273 pages including 14 figures, Postscript, See also ftp.scri.fsu.edu:hep-lat/papers/9210/9210032.ps.

Human SHBG mRNA Translation Is Modulated by Alternative 5′-Non-Coding Exons 1A and 1B

BACKGROUND: The human sex hormone-binding globulin (SHBG) gene comprises at least 6 different transcription units (TU-1, -1A, -1B, -1C, -1D and -1E), and is regulated by no less than 6 different promoters. The best characterized are TU-1 and TU-1A: TU-1 is responsible for producing plasma SHBG, while TU-1A is transcribed and translated in the testis. Transcription of the recently described TU-1B, -1C, and -1D has been demonstrated in human prostate tissue and prostate cancer cell lines, as well as in other human cell lines such as HeLa, HepG2, HeK 293, CW 9019 and imr 32. However, there are no reported data demonstrating their translation. In the present study, we aimed to determine whether TU-1A and TU-1B are indeed translated in the human prostate and whether 5' UTR exons 1A and 1B differently regulate SHBG translation. RESULTS: Cis-regulatory elements that could potentially regulate translation were identified within the 5'UTRs of SHBG TU-1A and TU-1B. Although full-length SHBG TU-1A and TU-1B mRNAs were present in prostate cancer cell lines, the endogenous SHBG protein was not detected by western blot in any of them. LNCaP prostate cancer cells transfected with several SHBG constructs containing exons 2 to 8 but lacking the 5'UTR sequence did show SHBG translation, whereas inclusion of the 5'UTR sequences of either exon 1A or 1B caused a dramatic decrease in SHBG protein levels. The molecular weight of SHBG did not vary between cells transfected with constructs with or without the 5'UTR sequence, thus confirming that the first in-frame ATG of exon 2 is the translation start site of TU-1A and TU-1B. CONCLUSIONS: The use of alternative SHBG first exons 1A and 1B differentially inhibits translation from the ATG situated in exon 2, which codes for methionine 30 of transcripts that begin with the exon 1 sequence

Functional informed genome-wide interaction analysis of body mass index, diabetes and colorectal cancer risk

Background: Body mass index (BMI) and diabetes are established risk factors for colorectal cancer (CRC), likely through perturbations in metabolic traits (e.g. insulin resistance and glucose homeostasis). Identification of interactions between variation in genes and these metabolic risk factors may identify novel biologic insights into CRC etiology. Methods: To improve statistical power and interpretation for gene-environment interaction (G × E) testing, we tested genetic variants that regulate expression of a gene together for interaction with BMI (kg/m2 ) and diabetes on CRC risk among 26 017 cases and 20 692 controls. Each variant was weighted based on PrediXcan analysis of gene expression data from colon tissue generated in the Genotype-Tissue Expression Project for all genes with heritability ≥1%. We used a mixed-effects model to jointly measure the G × E interaction in a gene by partitioning the interactions into the predicted gene expression levels (fixed effects), and residual G × E effects (random effects). G × BMI analyses were stratified by sex as BMI-CRC associations differ by sex. We used false discovery rates to account for multiple comparisons and reported all results with FDR <0.2. Results: Among 4839 genes tested, genetically predicted expressions of FOXA1 (P = 3.15 × 10-5 ), PSMC5 (P = 4.51 × 10-4 ) and CD33 (P = 2.71 × 10-4 ) modified the association of BMI on CRC risk for men; KIAA0753 (P = 2.29 × 10-5 ) and SCN1B (P = 2.76 × 10-4 ) modified the association of BMI on CRC risk for women; and PTPN2 modified the association between diabetes and CRC risk in both sexes (P = 2.31 × 10-5 ). Conclusions: Aggregating G × E interactions and incorporating functional information, we discovered novel genes that may interact with BMI and diabetes on CRC risk

Mendelian randomisation study of age at menarche and age at menopause and the risk of colorectal cancer.

BACKGROUND: Substantial evidence supports an association between use of menopausal hormone therapy and decreased colorectal cancer (CRC) risk, indicating a role of exogenous sex hormones in CRC development. However, findings on endogenous oestrogen exposure and CRC are inconsistent. METHODS: We used a Mendelian randomisation approach to test for a causal effect of age at menarche and age at menopause as surrogates for endogenous oestrogen exposure on CRC risk. Weighted genetic risk scores based on 358 single-nucleotide polymorphisms associated with age at menarche and 51 single-nucleotide polymorphisms associated with age at menopause were used to estimate the association with CRC risk using logistic regression in 12,944 women diagnosed with CRC and 10,741 women without CRC from three consortia. Sensitivity analyses were conducted to address pleiotropy and possible confounding by body mass index. RESULTS: Genetic risk scores for age at menarche (odds ratio per year 0.98, 95% confidence interval: 0.95-1.02) and age at menopause (odds ratio 0.98, 95% confidence interval: 0.94-1.01) were not significantly associated with CRC risk. The sensitivity analyses yielded similar results. CONCLUSIONS: Our study does not support a causal relationship between genetic risk scores for age at menarche and age at menopause and CRC risk

Mendelian randomization analysis of C-reactive protein on colorectal cancer risk

Background: Chronic inflammation is a risk factor for colorectal cancer (CRC). Circulating C-reactive protein (CRP) is also moderately associated with CRC risk. However, observational studies are susceptible to unmeasured confounding or reverse causality. Using genetic risk variants as instrumental variables, we investigated the causal relationship between genetically elevated CRP concentration and CRC risk, using a Mendelian randomization approach. Methods: Individual-level data from 30 480 CRC cases and 22 844 controls from 33 participating studies in three international consortia were used: the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO), the Colorectal Transdisciplinary Study (CORECT) and the Colon Cancer Family Registry (CCFR). As instrumental variables, we included 19 single nucleotide polymorphisms (SNPs) previously associated with CRP concentration. The SNP-CRC associations were estimated using a logistic regression model adjusted for age, sex, principal components and genotyping phases. An inverse-variance weighted method was applied to estimate the causal effect of CRP on CRC risk. Results: Among the 19 CRP-associated SNPs, rs1260326 and rs6734238 were significantly associated with CRC risk (P = 7.5 × 10-4, and P = 0.003, respectively). A genetically predicted one-unit increase in the log-transformed CRP concentrations (mg/l) was not associated with increased risk of CRC [odds ratio (OR) = 1.04; 95% confidence interval (CI): 0.97, 1.12; P = 0.256). No evidence of association was observed in subgroup analyses stratified by other risk factors. Conclusions: In spite of adequate statistical power to detect moderate association, we found genetically elevated CRP concentration was not associated with increased risk of CRC among individuals of European ancestry. Our findings suggested that circulating CRP is unlikely to be a causal factor in CRC development

Mendelian randomization of circulating polyunsaturated fatty acids and colorectal cancer risk

Background: Results from epidemiologic studies examining polyunsaturated fatty acids (PUFA) and colorectal cancer risk are inconsistent. Mendelian randomization may strengthen causal inference from observational studies. Given their shared metabolic pathway, examining the combined effects of aspirin/NSAID use with PUFAs could help elucidate an association between PUFAs and colorectal cancer risk. Methods: Information was leveraged from genome-wide association studies (GWAS) regarding PUFA-associated SNPs to create weighted genetic scores (wGS) representing genetically predicted circulating blood PUFAs for 11,016 non-Hispanic white colorectal cancer cases and 13,732 controls in the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO). Associations per SD increase in the wGS were estimated using unconditional logistic regression. Interactions between PUFA wGSs and aspirin/NSAID use on colorectal cancer risk were also examined. Results: Modest colorectal cancer risk reductions were observed per SD increase in circulating linoleic acid [ORLA = 0.96; 95% confidence interval (CI) = 0.93-0.98; P = 5.2 × 10-4] and α-linolenic acid (ORALA = 0.95; 95% CI = 0.92-0.97; P = 5.4 × 10-5), whereas modest increased risks were observed for arachidonic (ORAA = 1.06; 95% CI = 1.03-1.08; P = 3.3 × 10-5), eicosapentaenoic (OREPA = 1.04; 95% CI = 1.01-1.07; P = 2.5 × 10-3), and docosapentaenoic acids (ORDPA = 1.03; 95% CI = 1.01-1.06; P = 1.2 × 10-2). Each of these effects was stronger among aspirin/NSAID nonusers in the stratified analyses. Conclusions: Our study suggests that higher circulating shorter-chain PUFAs (i.e., LA and ALA) were associated with reduced colorectal cancer risk, whereas longer-chain PUFAs (i.e., AA, EPA, and DPA) were associated with an increased colorectal cancer risk. Impact: The interaction of PUFAs with aspirin/NSAID use indicates a shared colorectal cancer inflammatory pathway. Future research should continue to improve PUFA genetic instruments to elucidate the independent effects of PUFAs on colorectal cancer