Revealing myopathy spectrum: integrating transcriptional and clinical features of human skeletal muscles with varying health conditions

Myopathy refers to a large group of heterogeneous, rare muscle diseases. Bulk RNA-sequencing has been utilized for the diagnosis and research of these diseases for many years. However, the existing valuable sequencing data often lack integration and clinical interpretation. In this study, we integrated bulk RNA-sequencing data from 1221 human skeletal muscles (292 with myopathies, 929 controls) from both databases and our local samples. By applying a method similar to single-cell analysis, we revealed a general spectrum of muscle diseases, ranging from healthy to mild disease, moderate muscle wasting, and severe muscle disease. This spectrum was further partly validated in three specific myopathies (97 muscles) through clinical features including trinucleotide repeat expansion, magnetic resonance imaging fat fraction, pathology, and clinical severity scores. This spectrum helped us identify 234 genuinely healthy muscles as unprecedented controls, providing a new perspective for deciphering the hallmark genes and pathways among different myopathies. The newly identified featured genes of general myopathy, inclusion body myositis, and titinopathy were highly expressed in our local muscles, as validated by quantitative polymerase chain reaction.

Genome-wide polygenic risk scores for colorectal cancer have implications for risk-based screening

BACKGROUND

Hereditary factors, including single genetic variants and family history, can be used for targeting colorectal cancer (CRC) screening, but limited data exist on the impact of polygenic risk scores (PRS) on risk-based CRC screening.

METHODS

Using longitudinal health and genomics data on 453,733 Finnish individuals including 8801 CRC cases, we estimated the impact of a genome-wide CRC PRS on CRC screening initiation age through population-calibrated incidence estimation over the life course in men and women.

RESULTS

Compared to the cumulative incidence of CRC at age 60 in Finland (the current age for starting screening in Finland), a comparable cumulative incidence was reached 5 and 11 years earlier in persons with high PRS (80-99% and >99%, respectively), while those with a low PRS (< 20%) reached comparable incidence 7 years later. The PRS was associated with increased risk of post-colonoscopy CRC after negative colonoscopy (hazard ratio 1.76 per PRS SD, 95% CI 1.54-2.01). Moreover, the PRS predicted colorectal adenoma incidence and improved incident CRC risk prediction over non-genetic risk factors.

CONCLUSIONS

Our findings demonstrate that a CRC PRS can be used for risk stratification of CRC, with further research needed to optimally integrate the PRS into risk-based screening.

Designing patient-oriented combination therapies for acute myeloid leukemia based on efficacy/toxicity integration and bipartite network modeling

Acute myeloid leukemia (AML), a heterogeneous and aggressive blood cancer, does not respond well to single-drug therapy. A combination of drugs is required to effectively treat this disease. Computational models are critical for combination therapy discovery due to the tens of thousands of two-drug combinations, even with approved drugs. While predicting synergistic drugs is the focus of current methods, few consider drug efficacy and potential toxicity, which are crucial for treatment success. To find effective new drug candidates, we constructed a bipartite network using patient-derived tumor samples and drugs. The network is based on drug-response screening and summarizes all treatment response heterogeneity as drug response weights. This bipartite network is then projected onto the drug part, resulting in the drug similarity network. Distinct drug clusters were identified using community detection methods, each targeting different biological processes and pathways as revealed by enrichment and pathway analysis of the drugs' protein targets. Four drugs with the highest efficacy and lowest toxicity from each cluster were selected and tested for drug sensitivity using cell viability assays on various samples. Results show that ruxolitinib-ulixertinib and sapanisertib-LY3009120 are the most effective combinations with the least toxicity and the best synergistic effect on blast cells. These findings lay the foundation for personalized and successful AML therapies, ultimately leading to the development of drug combinations that can be used alongside standard first-line AML treatment.

Neutralizing antibodies after the third COVID-19 vaccination in healthcare workers with or without breakthrough infection

BACKGROUND

Vaccinations against the SARS-CoV-2 are still crucial in combating the ongoing pandemic that has caused more than 700 million infections and claimed almost 7 million lives in the past four years. Omicron (B.1.1.529) variants have incurred mutations that challenge the protection against infection and severe disease by the current vaccines, potentially compromising vaccination efforts.

METHODS

We analyzed serum samples taken up to 9 months post third dose from 432 healthcare workers. Enzyme-linked immunosorbent assays (ELISA) and microneutralization tests (MNT) were used to assess the prevalence of vaccine-induced neutralizing antibodies against various SARS-CoV-2 Omicron variants.

RESULTS

In this serological analysis we show that SARS-CoV-2 vaccine combinations of BNT162b2, mRNA-1273, and ChAdOx1 mount SARS-CoV-2 binding and neutralizing antibodies with similar kinetics, but with differing neutralization capabilities. The most recent Omicron variants, BQ.1.1 and XBB.1.5, show a significant increase in the ability to escape vaccine and infection-induced antibody responses. Breakthrough infections in thrice vaccinated adults were seen in over 50% of the vaccinees, resulting in a stronger antibody response than without infection.

CONCLUSIONS

Different three-dose vaccine combinations seem to induce considerable levels of neutralizing antibodies against most SARS-CoV-2 variants. However, the ability of the newer variants BQ1.1 and XBB 1.5 to escape vaccine-induced neutralizing antibody responses underlines the importance of updating vaccines as new variants emerge.

Spatial analysis of the osteoarthritis microenvironment: techniques, insights, and applications

Osteoarthritis (OA) is a debilitating degenerative disease affecting multiple joint tissues, including cartilage, bone, synovium, and adipose tissues. OA presents diverse clinical phenotypes and distinct molecular endotypes, including inflammatory, metabolic, mechanical, genetic, and synovial variants. Consequently, innovative technologies are needed to support the development of effective diagnostic and precision therapeutic approaches. Traditional analysis of bulk OA tissue extracts has limitations due to technical constraints, causing challenges in the differentiation between various physiological and pathological phenotypes in joint tissues. This issue has led to standardization difficulties and hindered the success of clinical trials. Gaining insights into the spatial variations of the cellular and molecular structures in OA tissues, encompassing DNA, RNA, metabolites, and proteins, as well as their chemical properties, elemental composition, and mechanical attributes, can contribute to a more comprehensive understanding of the disease subtypes. Spatially resolved biology enables biologists to investigate cells within the context of their tissue microenvironment, providing a more holistic view of cellular function. Recent advances in innovative spatial biology techniques now allow intact tissue sections to be examined using various -omics lenses, such as genomics, transcriptomics, proteomics, and metabolomics, with spatial data. This fusion of approaches provides researchers with critical insights into the molecular composition and functions of the cells and tissues at precise spatial coordinates. Furthermore, advanced imaging techniques, including high-resolution microscopy, hyperspectral imaging, and mass spectrometry imaging, enable the visualization and analysis of the spatial distribution of biomolecules, cells, and tissues. Linking these molecular imaging outputs to conventional tissue histology can facilitate a more comprehensive characterization of disease phenotypes. This review summarizes the recent advancements in the molecular imaging modalities and methodologies for in-depth spatial analysis. It explores their applications, challenges, and potential opportunities in the field of OA. Additionally, this review provides a perspective on the potential research directions for these contemporary approaches that can meet the requirements of clinical diagnoses and the establishment of therapeutic targets for OA.

Metabolic changes in response to varying whole-grain wheat and rye intake

Epidemiological studies have shown associations between whole-grain intake and lowered disease risk. A sufficient level of whole-grain intake to reach the health benefits has not been established, and there is limited knowledge about the impact of whole-grain intake on metabolite levels. In this clinical intervention study, we aimed to identify plasma and urine metabolites associated with two different intake levels of whole-grain wheat and rye and to correlate them with clinical plasma biomarkers. Healthy volunteers (N = 68) were divided into two groups receiving either whole-grain wheat or whole-grain rye in two four-week interventions with 48 and 96 g/d of whole grains consumed. The metabolomics of the plasma samples was performed with UPLC-QTOF-MS. Plasma alkylresorcinols were quantified with GC-MS and plasma and urinary mammalian lignans with HPLC-ECD. The high-dose intervention impacted the metabolite profile, including microbial metabolites, more in the rye-enriched diet compared with wheat. Among the increased metabolites were alkylresorcinol glucuronides, sinapyl alcohol, and pipecolic acid betaine, while the decreased metabolites included acylcarnitines and ether lipids. Plasma alkylresorcinols, urinary enterolactone, and total mammalian lignans reflected the study diets in a dose-dependent manner. Several key metabolites linked with whole-grain consumption and gut microbial metabolism increased in a linear manner between the two interventions. The results reveal that an increase in whole-grain intake, particularly rye, is strongly reflected in the metabolite profile, is correlated with clinical variables, and suggests that a diet rich in whole grains promotes the growth and/or metabolism of microbes producing potentially beneficial microbial metabolites.

Recessive TMOD1 mutation causes childhood cardiomyopathy

Familial cardiomyopathy in pediatric stages is a poorly understood presentation of heart disease in children that is attributed to pathogenic mutations. Through exome sequencing, we report a homozygous variant in tropomodulin 1 (TMOD1; c.565C>T, p.R189W) in three individuals from two unrelated families with childhood-onset dilated and restrictive cardiomyopathy. To decipher the mechanism of pathogenicity of the R189W mutation in TMOD1, we utilized a wide array of methods, including protein analyses, biochemistry and cultured cardiomyocytes. Structural modeling revealed potential defects in the local folding of TMOD1R189W and its affinity for actin. Cardiomyocytes expressing GFP-TMOD1R189W demonstrated longer thin filaments than GFP-TMOD1wt-expressing cells, resulting in compromised filament length regulation. Furthermore, TMOD1R189W showed weakened activity in capping actin filament pointed ends, providing direct evidence for the variant's effect on actin filament length regulation. Our data indicate that the p.R189W variant in TMOD1 has altered biochemical properties and reveals a unique mechanism for childhood-onset cardiomyopathy.

Boosting cytotoxicity of adoptive allogeneic NK cell therapy with an oncolytic adenovirus encoding a human vIL-2 cytokine for the treatment of human ovarian cancer

Despite good results in the treatment of hematological malignancies, Natural killer (NK) cells have shown limited effectiveness in solid tumors, such as ovarian cancer (OvCa). Here, we assessed the potential of an oncolytic adenovirus expressing a variant interleukin-2 (vIL-2) cytokine, Ad5/3-E2F-d24-vIL2 (vIL-2 virus), also known as TILT-452, to enhance NK cell therapy efficacy in human OvCa ex vivo. Human OvCa surgical specimens were processed into single-cell suspensions and NK cells were expanded from healthy blood donors. OvCa sample digests were co-cultured ex vivo with NK cells and vIL-2 virus and cancer cell killing potential assessed in real time through cell impedance measurement. Proposed therapeutic combination was evaluated in vivo with an OvCa patient-derived xenograft (PDX) in mice. Addition of vIL-2 virus significantly enhanced NK cell therapy killing potential in treated OvCa co-cultures. Similarly, vIL-2 virus in combination with NK cell therapy promoted the best in vivo OvCa tumor control. Mechanistically, vIL-2 virus induced higher percentages of granzyme B in NK cells, and CD8+ T cells, while T regulatory cell proportions remained comparable to NK cell monotherapy in vivo. Ad5/3-E2F-d24-vIL2 virus treatment represents a promising strategy to boost adoptive NK cell therapeutic effect in human OvCa.

Inhibition of TGF-β signaling, invasion, and growth of cutaneous squamous cell carcinoma by PLX8394

Cutaneous squamous cell carcinoma (cSCC) is the most common metastatic skin cancer. The prognosis of patients with metastatic cSCC is poor emphasizing the need for new therapies. We have previously reported that the activation of Ras/MEK/ERK1/2 and transforming growth factor β (TGF-β)/Smad2 signaling in transformed keratinocytes and cSCC cells leads to increased accumulation of laminin-332 and accelerated invasion. Here, we show that the next-generation B-Raf inhibitor PLX8394 blocks TGF-β signaling in ras-transformed metastatic epidermal keratinocytes (RT3 cells) harboring wild-type B-Raf and hyperactive Ras. PLX8394 decreased phosphorylation of TGF-β receptor II and Smad2, as well as p38 activity, MMP-1 and MMP-13 synthesis, and laminin-332 accumulation. PLX8394 significantly inhibited the growth of human cSCC tumors and in vivo collagen degradation in xenograft model. In conclusion, our data indicate that PLX8394 inhibits several serine-threonine kinases in malignantly transformed human keratinocytes and cSCC cells and inhibits cSCC invasion and tumor growth in vitro and in vivo. We identify PLX8394 as a potential therapeutic compound for advanced human cSCC.

Development of mesothelioma-specific oncolytic immunotherapy enabled by immunopeptidomics of murine and human mesothelioma tumors

Malignant pleural mesothelioma (MPM) is an aggressive tumor with a poor prognosis. As the available therapeutic options show a lack of efficacy, novel therapeutic strategies are urgently needed. Given its T-cell infiltration, we hypothesized that MPM is a suitable target for therapeutic cancer vaccination. To date, research on mesothelioma has focused on the identification of molecular signatures to better classify and characterize the disease, and little is known about therapeutic targets that engage cytotoxic (CD8+) T cells. In this study we investigate the immunopeptidomic antigen-presented landscape of MPM in both murine (AB12 cell line) and human cell lines (H28, MSTO-211H, H2452, and JL1), as well as in patients' primary tumors. Applying state-of-the-art immuno-affinity purification methodologies, we identify MHC I-restricted peptides presented on the surface of malignant cells. We characterize in vitro the immunogenicity profile of the eluted peptides using T cells from human healthy donors and cancer patients. Furthermore, we use the most promising peptides to formulate an oncolytic virus-based precision immunotherapy (PeptiCRAd) and test its efficacy in a mouse model of mesothelioma in female mice. Overall, we demonstrate that the use of immunopeptidomic analysis in combination with oncolytic immunotherapy represents a feasible and effective strategy to tackle untreatable tumors.

Improving the cytotoxic response of tumor-infiltrating lymphocytes towards advanced stage ovarian cancer with an oncolytic adenovirus expressing a human vIL-2 cytokine

While the presence of tumor-infiltrating lymphocytes (TILs) associates with improved survival prognosis in ovarian cancer (OvCa) patients, TIL therapy benefit is limited. Here, we evaluated an oncolytic adenovirus coding for a human variant IL-2 (vIL-2) cytokine, Ad5/3-E2F-d24-vIL2 (vIL-2 virus), also known as TILT-452, as an immunotherapeutic strategy to enhance TIL responsiveness towards advanced stage OvCa tumors. Fragments of resected human OvCa tumors were processed into single-cell suspensions, and autologous TILs were expanded from said samples. OvCa tumor specimens were co-cultured with TILs plus vIL-2 virus, and cell killing was assessed in real time through cell impedance measurement. Combination therapy was further evaluated in vivo through a patient-derived xenograft (PDX) ovarian cancer murine model. The combination of vIL-2 virus plus TILs had best cancer cell killing ex vivo compared to TILs monotherapy. These results were supported by an in vivo experiment, where the best OvCa tumor control was obtained when vIL-2 virus was added to TIL therapy. Furthermore, the proposed therapy induced a highly cytotoxic phenotype demonstrated by increased granzyme B intensity in NK cells, CD4+ T, and CD8+ T cells in treated tumors. Our results demonstrate that Ad5/3-E2F-d24-vIL2 therapy consistently improved TILs therapy cytotoxicity in treated human OvCa tumors.

Mixed effects of a national protected area network on terrestrial and freshwater biodiversity

Protected areas are considered fundamental to counter biodiversity loss. However, evidence for their effectiveness in averting local extinctions remains scarce and taxonomically biased. We employ a robust counterfactual multi-taxon approach to compare occupancy patterns of 638 species, including birds (150), mammals (23), plants (39) and phytoplankton (426) between protected and unprotected sites across four decades in Finland. We find mixed impacts of protected areas, with only a small proportion of species explicitly benefiting from protection-mainly through slower rates of decline inside protected areas. The benefits of protection are enhanced for larger protected areas and are traceable to when the sites were protected, but are mostly unrelated to species conservation status or traits (size, climatic niche and threat status). Our results suggest that the current protected area network can partly contribute to slow down declines in occupancy rates, but alone will not suffice to halt the biodiversity crisis. Efforts aimed at improving coverage, connectivity and management will be key to enhance the effectiveness of protected areas towards bending the curve of biodiversity loss.

Transposable elements as tissue-specific enhancers in cancers of endodermal lineage

Transposable elements (TE) are repetitive genomic elements that harbor binding sites for human transcription factors (TF). A regulatory role for TEs has been suggested in embryonal development and diseases such as cancer but systematic investigation of their functions has been limited by their widespread silencing in the genome. Here, we utilize unbiased massively parallel reporter assay data using a whole human genome library to identify TEs with functional enhancer activity in two human cancer types of endodermal lineage, colorectal and liver cancers. We show that the identified TE enhancers are characterized by genomic features associated with active enhancers, such as epigenetic marks and TF binding. Importantly, we identify distinct TE subfamilies that function as tissue-specific enhancers, namely MER11- and LTR12-elements in colon and liver cancers, respectively. These elements are bound by distinct TFs in each cell type, and they have predicted associations to differentially expressed genes. In conclusion, these data demonstrate how different cancer types can utilize distinct TEs as tissue-specific enhancers, paving the way for comprehensive understanding of the role of TEs as bona fide enhancers in the cancer genomes.

Signal peptide mimicry primes Sec61 for client-selective inhibition

Preventing the biogenesis of disease-relevant proteins is an attractive therapeutic strategy, but attempts to target essential protein biogenesis factors have been hampered by excessive toxicity. Here we describe KZR-8445, a cyclic depsipeptide that targets the Sec61 translocon and selectively disrupts secretory and membrane protein biogenesis in a signal peptide-dependent manner. KZR-8445 potently inhibits the secretion of pro-inflammatory cytokines in primary immune cells and is highly efficacious in a mouse model of rheumatoid arthritis. A cryogenic electron microscopy structure reveals that KZR-8445 occupies the fully opened Se61 lateral gate and blocks access to the lumenal plug domain. KZR-8445 binding stabilizes the lateral gate helices in a manner that traps select signal peptides in the Sec61 channel and prevents their movement into the lipid bilayer. Our results establish a framework for the structure-guided discovery of novel therapeutics that selectively modulate Sec61-mediated protein biogenesis.

Structural mechanism for inhibition of PP2A-B56α and oncogenicity by CIP2A

The protein phosphatase 2A (PP2A) heterotrimer PP2A-B56α is a human tumour suppressor. However, the molecular mechanisms inhibiting PP2A-B56α in cancer are poorly understood. Here, we report molecular level details and structural mechanisms of PP2A-B56α inhibition by an oncoprotein CIP2A. Upon direct binding to PP2A-B56α trimer, CIP2A displaces the PP2A-A subunit and thereby hijacks both the B56α, and the catalytic PP2Ac subunit to form a CIP2A-B56α-PP2Ac pseudotrimer. Further, CIP2A competes with B56α substrate binding by blocking the LxxIxE-motif substrate binding pocket on B56α. Relevant to oncogenic activity of CIP2A across human cancers, the N-terminal head domain-mediated interaction with B56α stabilizes CIP2A protein. Functionally, CRISPR/Cas9-mediated single amino acid mutagenesis of the head domain blunted MYC expression and MEK phosphorylation, and abrogated triple-negative breast cancer in vivo tumour growth. Collectively, we discover a unique multi-step hijack and mute protein complex regulation mechanism resulting in tumour suppressor PP2A-B56α inhibition. Further, the results unfold a structural determinant for the oncogenic activity of CIP2A, potentially facilitating therapeutic modulation of CIP2A in cancer and other diseases.

A competitive precision CRISPR method to identify the fitness effects of transcription factor binding sites

Here we describe a competitive genome editing method that measures the effect of mutations on molecular functions, based on precision CRISPR editing using template libraries with either the original or altered sequence, and a sequence tag, enabling direct comparison between original and mutated cells. Using the example of the MYC oncogene, we identify important transcriptional targets and show that E-box mutations at MYC target gene promoters reduce cellular fitness.

Comparative analysis of COVID-19 vaccine responses and third booster dose-induced neutralizing antibodies against Delta and Omicron variants

Two COVID-19 mRNA (of BNT162b2, mRNA-1273) and two adenovirus vector vaccines (ChAdOx1 and Janssen) are licensed in Europe, but optimization of regime and dosing is still ongoing. Here we show in health care workers (n = 328) that two doses of BNT162b2, mRNA-1273, or a combination of ChAdOx1 adenovirus vector and mRNA vaccines administrated with a long 12-week dose interval induce equally high levels of anti-SARS-CoV-2 spike antibodies and neutralizing antibodies against D614 and Delta variant. By contrast, two doses of BNT162b2 with a short 3-week interval induce 2-3-fold lower titers of neutralizing antibodies than those from the 12-week interval, yet a third BNT162b2 or mRNA-1273 booster dose increases the antibody levels 4-fold compared to the levels after the second dose, as well as induces neutralizing antibody against Omicron BA.1 variant. Our data thus indicates that a third COVID-19 mRNA vaccine may induce cross-protective neutralizing antibodies against multiple variants.

Human cell transformation by combined lineage conversion and oncogene expression

Cancer is the most complex genetic disease known, with mutations implicated in more than 250 genes. However, it is still elusive which specific mutations found in human patients lead to tumorigenesis. Here we show that a combination of oncogenes that is characteristic of liver cancer (CTNNB1, TERT, MYC) induces senescence in human fibroblasts and primary hepatocytes. However, reprogramming fibroblasts to a liver progenitor fate, induced hepatocytes (iHeps), makes them sensitive to transformation by the same oncogenes. The transformed iHeps are highly proliferative, tumorigenic in nude mice, and bear gene expression signatures of liver cancer. These results show that tumorigenesis is triggered by a combination of three elements: the set of driver mutations, the cellular lineage, and the state of differentiation of the cells along the lineage. Our results provide direct support for the role of cell identity as a key determinant in transformation and establish a paradigm for studying the dynamic role of oncogenic drivers in human tumorigenesis.

Gaussian Light Model in Brightfield Optical Projection Tomography

This study focuses on improving the reconstruction process of the brightfield optical projection tomography (OPT). OPT is often described as the optical equivalent of X-ray computed tomography, but based on visible light. The detection optics used to collect light in OPT focus on a certain distance and induce blurring in those features out of focus. However, the conventionally used inverse Radon transform assumes an absolute focus throughout the propagation axis. In this study, we model the focusing properties of the detection by coupling Gaussian beam model (GBM) with the Radon transform. The GBM enables the construction of a projection operator that includes modeling of the blurring caused by the light beam. We also introduce the concept of a stretched GBM (SGBM) in which the Gaussian beam is scaled in order to avoid the modeling errors related to the determination of the focal plane. Furthermore, a thresholding approach is used to compress memory usage. We tested the GBM and SGBM approaches using simulated and experimental data in mono- and multifocal modes. When compared with the traditionally used filtered backprojection algorithm, the iteratively computed reconstructions, including the Gaussian models GBM and SGBM, provided smoother images with higher contrast.