Analysis of incidental findings in Qatar genome participants reveals novel functional variants in LMNA and DSP.

In order to report clinically actionable incidental findings in genetic testing, the American College of Medical Genetics and Genomics (ACMG) recommended the evaluation of variants in 59 genes associated with highly penetrant mutations. However, there is a lack of epidemiological data on medically actionable rare variants in these genes in Arab populations. We used whole genome sequencing data from 6045 participants from the Qatar Genome Programme and integrated it with phenotypic data collected by the Qatar Biobank. We identified novel putative pathogenic variants in the 59 ACMG genes by filtering previously unrecorded variants based on computational prediction of pathogenicity, variant rarity and segregation evidence. We assessed the phenotypic associations of candidate variants in genes linked to cardiovascular diseases. Finally, we used a zebrafish knockdown and synthetic human mRNA co-injection assay to functionally characterize two of these novel variants. We assessed the zebrafish cardiac function in terms of heart rate, rhythm and hemodynamics, as well as the heart structure. We identified 52 492 novel variants, which have not been reported in global and disease-specific databases. A total of 74 novel variants were selected with potentially pathogenic effect. We prioritized two novel cardiovascular variants, DSP c.1841A > G (p.Asp614Gly) and LMNA c.326 T > G (p.Val109Gly) for functional characterization. Our results showed that both variants resulted in abnormal zebrafish heart rate, rhythm and structure. This study highlights medically actionable variants that are specific to the Middle Eastern Qatari population.

BRCA1-specific machine learning model predicts variant pathogenicity with high accuracy.

Identification of novel BRCA1 variants outpaces their clinical annotation which highlights the importance of developing accurate computational methods for risk assessment. Therefore our aim was to develop a BRCA1-specific machine learning model to predict the pathogenicity of all types of BRCA1 variants and to apply this model and our previous BRCA2-specific model to assess BRCA variants of uncertain significance (VUS) among Qatari patients with breast cancer. We developed an XGBoost model that utilizes variant information such as position frequency and consequence as well as prediction scores from numerous in silico tools. We trained and tested the model with BRCA1 variants that were reviewed and classified by the Evidence-Based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) consortium. In addition we tested the model's performance on an independent set of missense variants of uncertain significance with experimentally determined functional scores. The model performed excellently in predicting the pathogenicity of ENIGMA-classified variants (accuracy: 99.9%) and in predicting the functional consequence of the independent set of missense variants (accuracy: 93.4%). Moreover it predicted 2 115 potentially pathogenic variants among the 31 058 unreviewed BRCA1 variants in the BRCA exchange database. Using two BRCA-specific models we did not identify any pathogenic BRCA1 variants among those found in patients in Qatar but predicted four potentially pathogenic BRCA2 variants, which could be prioritized for functional validation.

Genomic and pleiotropic analyses of resting QT interval identifies novel loci and overlap with atrial electrical disorders.

The resting QT interval, an electrocardiographic (ECG) measure of ventricular myocardial repolarization, is a heritable risk marker of cardiovascular mortality, but the mechanisms remain incompletely understood. Previously reported candidate genes have provided insights into the regulatory mechanisms of the QT interval. However, there are still important knowledge gaps. We aimed to gain new insights by (i) providing new candidate genes, (ii) identifying pleiotropic associations with other cardiovascular traits, and (iii) scanning for sexually dimorphic genetic effects. We conducted a genome-wide association analysis for resting QT interval with ~9.8 million variants in 52 107 individuals of European ancestry without known cardiovascular disease from the UK Biobank. We identified 40 loci, 13 of which were novel, including 2 potential sex-specific loci, explaining ~11% of the trait variance. Candidate genes at novel loci were involved in myocardial structure and arrhythmogenic cardiomyopathy. Investigation of pleiotropic effects of QT interval variants using phenome-wide association analyses in 302 000 unrelated individuals from the UK Biobank and pairwise genome-wide comparisons with other ECG and cardiac imaging traits revealed genetic overlap with atrial electrical pathology. These findings provide novel insights into how abnormal myocardial repolarization and increased cardiovascular mortality may be linked.

Assessing the genetic burden of familial hypercholesterolemia in a large middle eastern biobank.

BACKGROUND: The genetic architecture underlying Familial Hypercholesterolemia (FH) in Middle Eastern Arabs is yet to be fully described, and approaches to assess this from population-wide biobanks are important for public health planning and personalized medicine. METHODS: We evaluate the pilot phase cohort (n = 6,140 adults) of the Qatar Biobank (QBB) for FH using the Dutch Lipid Clinic Network (DLCN) criteria, followed by an in-depth characterization of all genetic alleles in known dominant (LDLR, APOB, and PCSK9) and recessive (LDLRAP1, ABCG5, ABCG8, and LIPA) FH-causing genes derived from whole-genome sequencing (WGS). We also investigate the utility of a globally established 12-SNP polygenic risk score to predict FH individuals in this cohort with Arab ancestry. RESULTS: Using DLCN criteria, we identify eight (0.1%) 'definite', 41 (0.7%) 'probable' and 334 (5.4%) 'possible' FH individuals, estimating a prevalence of 'definite or probable' FH in the Qatari cohort of ~ 1:125. We identify ten previously known pathogenic single-nucleotide variants (SNVs) and 14 putatively novel SNVs, as well as one novel copy number variant in PCSK9. Further, despite the modest sample size, we identify one homozygote for a known pathogenic variant (ABCG8, p. Gly574Arg, global MAF = 4.49E-05) associated with Sitosterolemia 2. Finally, calculation of polygenic risk scores found that individuals with 'definite or probable' FH have a significantly higher LDL-C SNP score than 'unlikely' individuals (p = 0.0003), demonstrating its utility in Arab populations. CONCLUSION: We design and implement a standardized approach to phenotyping a population biobank for FH risk followed by systematically identifying known variants and assessing putative novel variants contributing to FH burden in Qatar. Our results motivate similar studies in population-level biobanks - especially those with globally under-represented ancestries - and highlight the importance of genetic screening programs for early detection and management of individuals with high FH risk in health systems.

Flavin Adenine Dinucleotide (FAD) and Pyridoxal 5'-Phosphate (PLP) Bind to Sox9 and Alter the Expression of Key Pancreatic Progenitor Transcription Factors.

Cofactor flavin adenine dinucleotide (FAD), a compound with flavin moiety and a derivative of riboflavin (vitamin B2), is shown to bind to Sox9 (a key transcription factor in early pancreatic development) and, subsequently, induce a large increase in markers of pancreatic development, including Ngn3 and PTF1a. Pyridoxal 5'-phosphate (PLP), the active form of vitamin B6, also binds to Sox9 and results in a similar increase in pancreatic development markers. Sox9 is known to be specifically important for pancreatic progenitors. Previously, there was no known link between FAD, PLP, or other co-factors and Sox9 for function. Thus, our findings show the mechanism by which FAD and PLP interact with Sox9 and result in the altered expression of pancreatic progenitor transcription factors involved in the pancreas development.

Actionable genomic variants in 6045 participants from the Qatar Genome Program.

In a clinical setting, DNA sequencing can uncover findings unrelated to the purpose of genetic evaluation. The American College of Medical Genetics and Genomics (ACMG) recommends the evaluation and reporting of 59 genes from clinic genomic sequencing. While the prevalence of secondary findings is available from large population studies, these data lack Arab and other Middle Eastern populations. The Qatar Genome Program (QGP) generates whole-genome sequencing (WGS) data and combines it with phenotypic information to create a comprehensive database for studying the Qatari and wider Arab and Middle Eastern populations at the molecular level. This study identified and analyzed medically actionable variants in the 59 ACMG genes using WGS data from 6045 QGP participants. Our results identified a total of 60 pathogenic and likely pathogenic variants in 25 ACMG genes in 141 unique individuals. Overall, 2.3% of the QGP sequenced participants carried a pathogenic or likely pathogenic variant in one of the 59 ACMG genes. We evaluated the QGP phenotype-genotype association of additional nonpathogenic ACMG variants. These variants were found in patients from the Hamad Medical Corporation or reported incidental findings data in Qatar. We found a significant phenotype association for two variants, c.313+3A>C in LDLR, and c.58C>T (p.Gln20*) in the TPM1.

Frequency and management of medically actionable incidental findings from genome and exome sequencing data: a systematic review.

The application of whole genome/exome sequencing technologies in clinical genetics and research has resulted in the discovery of incidental findings unrelated to the primary purpose of genetic testing. The American College of Medical Genetics and Genomics published guidelines for reporting pathogenic and likely pathogenic variants that are deemed to be medically actionable, which allowed us to learn about the epidemiology of incidental findings in different populations. However, consensus guidelines for variant reporting and classification are still lacking. We conducted a systematic literature review of incidental findings in whole genome/exome sequencing studies to obtain a comprehensive understanding of variable reporting and classification methods for variants that are deemed to be medically actionable across different populations. The review highlights the elements that demand further consideration or adjustment.

A Comprehensive Overview of the COVID-19 Literature: Machine Learning-Based Bibliometric Analysis.

BACKGROUND: Shortly after the emergence of COVID-19, researchers rapidly mobilized to study numerous aspects of the disease such as its evolution, clinical manifestations, effects, treatments, and vaccinations. This led to a rapid increase in the number of COVID-19-related publications. Identifying trends and areas of interest using traditional review methods (eg, scoping and systematic reviews) for such a large domain area is challenging. OBJECTIVE: We aimed to conduct an extensive bibliometric analysis to provide a comprehensive overview of the COVID-19 literature. METHODS: We used the COVID-19 Open Research Dataset (CORD-19) that consists of a large number of research articles related to all coronaviruses. We used a machine learning-based method to analyze the most relevant COVID-19-related articles and extracted the most prominent topics. Specifically, we used a clustering algorithm to group published articles based on the similarity of their abstracts to identify research hotspots and current research directions. We have made our software accessible to the community via GitHub. RESULTS: Of the 196,630 publications retrieved from the database, we included 28,904 in our analysis. The mean number of weekly publications was 990 (SD 789.3). The country that published the highest number of COVID-19-related articles was China (2950/17,270, 17.08%). The highest number of articles were published in bioRxiv. Lei Liu affiliated with the Southern University of Science and Technology in China published the highest number of articles (n=46). Based on titles and abstracts alone, we were able to identify 1515 surveys, 733 systematic reviews, 512 cohort studies, 480 meta-analyses, and 362 randomized control trials. We identified 19 different topics covered among the publications reviewed. The most dominant topic was public health response, followed by clinical care practices during the COVID-19 pandemic, clinical characteristics and risk factors, and epidemic models for its spread. CONCLUSIONS: We provide an overview of the COVID-19 literature and have identified current hotspots and research directions. Our findings can be useful for the research community to help prioritize research needs and recognize leading COVID-19 researchers, institutes, countries, and publishers. Our study shows that an AI-based bibliometric analysis has the potential to rapidly explore a large corpus of academic publications during a public health crisis. We believe that this work can be used to analyze other eHealth-related literature to help clinicians, administrators, and policy makers to obtain a holistic view of the literature and be able to categorize different topics of the existing research for further analyses. It can be further scaled (for instance, in time) to clinical summary documentation. Publishers should avoid noise in the data by developing a way to trace the evolution of individual publications and unique authors.

The biological impact of blood pressure-associated genetic variants in the natriuretic peptide receptor C gene on human vascular smooth muscle.

Elevated blood pressure (BP) is a major global risk factor for cardiovascular disease. Genome-wide association studies have identified several genetic variants at the NPR3 locus associated with BP, but the functional impact of these variants remains to be determined. Here we confirmed, by a genome-wide association study within UK Biobank, the existence of two independent BP-related signals within NPR3 locus. Using human primary vascular smooth muscle cells (VSMCs) and endothelial cells (ECs) from different individuals, we found that the BP-elevating alleles within one linkage disequilibrium block identified by the sentinel variant rs1173771 was associated with lower endogenous NPR3 mRNA and protein levels in VSMCs, together with reduced levels in open chromatin and nuclear protein binding. The BP-elevating alleles also increased VSMC proliferation, angiotensin II-induced calcium flux and cell contraction. However, an analogous genotype-dependent association was not observed in vascular ECs. Our study identifies novel, putative mechanisms for BP-associated variants at the NPR3 locus to elevate BP, further strengthening the case for targeting NPR-C as a therapeutic approach for hypertension and cardiovascular disease prevention.

Discovery of novel heart rate-associated loci using the Exome Chip.

Resting heart rate is a heritable trait, and an increase in heart rate is associated with increased mortality risk. Genome-wide association study analyses have found loci associated with resting heart rate, at the time of our study these loci explained 0.9% of the variation. This study aims to discover new genetic loci associated with heart rate from Exome Chip meta-analyses.Heart rate was measured from either elecrtrocardiograms or pulse recordings. We meta-analysed heart rate association results from 104 452 European-ancestry individuals from 30 cohorts, genotyped using the Exome Chip. Twenty-four variants were selected for follow-up in an independent dataset (UK Biobank, N = 134 251). Conditional and gene-based testing was undertaken, and variants were investigated with bioinformatics methods.We discovered five novel heart rate loci, and one new independent low-frequency non-synonymous variant in an established heart rate locus (KIAA1755). Lead variants in four of the novel loci are non-synonymous variants in the genes C10orf71, DALDR3, TESK2 and SEC31B. The variant at SEC31B is significantly associated with SEC31B expression in heart and tibial nerve tissue. Further candidate genes were detected from long-range regulatory chromatin interactions in heart tissue (SCD, SLF2 and MAPK8). We observed significant enrichment in DNase I hypersensitive sites in fetal heart and lung. Moreover, enrichment was seen for the first time in human neuronal progenitor cells (derived from embryonic stem cells) and fetal muscle samples by including our novel variants.Our findings advance the knowledge of the genetic architecture of heart rate, and indicate new candidate genes for follow-up functional studies.

The pluripotent regulatory circuitry connecting promoters to their long-range interacting elements.

The mammalian genome harbors up to one million regulatory elements often located at great distances from their target genes. Long-range elements control genes through physical contact with promoters and can be recognized by the presence of specific histone modifications and transcription factor binding. Linking regulatory elements to specific promoters genome-wide is currently impeded by the limited resolution of high-throughput chromatin interaction assays. Here we apply a sequence capture approach to enrich Hi-C libraries for >22,000 annotated mouse promoters to identify statistically significant, long-range interactions at restriction fragment resolution, assigning long-range interacting elements to their target genes genome-wide in embryonic stem cells and fetal liver cells. The distal sites contacting active genes are enriched in active histone modifications and transcription factor occupancy, whereas inactive genes contact distal sites with repressive histone marks, demonstrating the regulatory potential of the distal elements identified. Furthermore, we find that coregulated genes cluster nonrandomly in spatial interaction networks correlated with their biological function and expression level. Interestingly, we find the strongest gene clustering in ES cells between transcription factor genes that control key developmental processes in embryogenesis. The results provide the first genome-wide catalog linking gene promoters to their long-range interacting elements and highlight the complex spatial regulatory circuitry controlling mammalian gene expression.

Genome-wide association study of response to methotrexate in early rheumatoid arthritis patients.

Methotrexate (MTX) monotherapy is a common first treatment for rheumatoid arthritis (RA), but many patients do not respond adequately. In order to identify genetic predictors of response, we have combined data from two consortia to carry out a genome-wide study of response to MTX in 1424 early RA patients of European ancestry. Clinical endpoints were change from baseline to 6 months after starting treatment in swollen 28-joint count, tender 28-joint count, C-reactive protein and the overall 3-component disease activity score (DAS28). No single nucleotide polymorphism (SNP) reached genome-wide statistical significance for any outcome measure. The strongest evidence for association was with rs168201 in NRG3 (p = 10-7 for change in DAS28). Some support was also seen for association with ZMIZ1, previously highlighted in a study of response to MTX in juvenile idiopathic arthritis. Follow-up in two smaller cohorts of 429 and 177 RA patients did not support these findings, although these cohorts were more heterogeneous.

Capturing genomic relationships that matter.

There is a strong interrelationship within the cell nucleus between form and function of the genome. This connection is exhibited across multiple hierarchies, ranging from grand-scale positioning of chromosomes and their intersection with specific nuclear functional activities, the segregation of chromosome structure into distinct domains and long-range regulatory contacts that drive spatial and temporal expression patterns of genes. Fifteen years ago, the development of the chromosome conformation capture method placed the nature of specific, long-range regulatory interactions under scrutiny. However, its development and integration with next-generation sequencing technologies has greatly expanded the breadth and scope of what is detected. The sheer scale of data offered by these important advances has come with new and challenging bottlenecks that are both experimental and bioinformatical. Here, we discuss the recent and prospective development and implementation of new methodologies and analytical tools that are allowing an in-depth, yet focussed characterisation of genomic contacts that are associated with functional activities in the nucleus.

Analysis of 14,392 whole genomes reveals 3.5% of Qataris carry medically actionable variants.

Arabic populations are underrepresented in large genome projects; therefore, the frequency of clinically actionable variants among Arabs is largely unknown. Here, we investigated genetic variation in 14,392 whole genomes from the Qatar Genome Program (QGP) across the list of 78 actionable genes (v3.1) determined by the American College of Medical Genetics and Genomics (ACMG). Variants were categorized into one of the following groups: (1) Pathogenic (P), (2) Likely pathogenic (LP), and (3) Rare variants of uncertain significance with evidence of pathogenicity. For the classification, we used variant databases, effect predictors, and the disease-relevant phenotypes available for the cohort. Data on cardiovascular disease, cancer, and hypercholesterolemia allowed us to assess the disease-relevant phenotype association of rare missense variants. We identified 248 distinct variants in 50 ACMG genes that fulfilled our criteria to be included in one of the three groups affecting 1036 genotype-positive participants of the QGP cohort. The most frequent variants were in TTN, followed by RYR1 and ATP7B. The prevalence of reportable secondary findings was 3.5%. A further 46 heterozygous variants in six genes with an autosomal recessive mode of inheritance were detected in 200 individuals, accounting for an additional 1.4%. Altogether, they affect 5% of the population. Due to the high consanguinity rate in the QGP cohort (28% in spouses and 60% in parents), P and LP variants both in genes with dominant and recessive inheritance are important for developing better treatment options and preventive strategies in Qatar and the Arabic population of the Middle East.

Loss of the yeast transporter Agp2 upregulates the pleiotropic drug-resistant pump Pdr5 and confers resistance to the protein synthesis inhibitor cycloheximide.

The transmembrane protein Agp2, initially shown as a transporter of L-carnitine, mediates the high-affinity transport of polyamines and the anticancer drug bleomycin-A5. Cells lacking Agp2 are hyper-resistant to polyamine and bleomycin-A5. In these earlier studies, we showed that the protein synthesis inhibitor cycloheximide blocked the uptake of bleomycin-A5 into the cells suggesting that the drug uptake system may require de novo synthesis. However, our recent findings demonstrated that cycloheximide, instead, induced rapid degradation of Agp2, and in the absence of Agp2 cells are resistant to cycloheximide. These observations raised the possibility that the degradation of Agp2 may allow the cell to alter its drug resistance network to combat the toxic effects of cycloheximide. In this study, we show that membrane extracts from agp2Δ mutants accentuated several proteins that were differentially expressed in comparison to the parent. Mass spectrometry analysis of the membrane extracts uncovered the pleiotropic drug efflux pump, Pdr5, involved in the efflux of cycloheximide, as a key protein upregulated in the agp2Δ mutant. Moreover, a global gene expression analysis revealed that 322 genes were differentially affected in the agp2Δ mutant versus the parent, including the prominent PDR5 gene and genes required for mitochondrial function. We further show that Agp2 is associated with the upstream region of the PDR5 gene, leading to the hypothesis that cycloheximide resistance displayed by the agp2Δ mutant is due to the derepression of the PDR5 gene.

Multi-omics and machine learning reveal context-specific gene regulatory activities of PML::RARA in acute promyelocytic leukemia.

The PML::RARA fusion protein is the hallmark driver of Acute Promyelocytic Leukemia (APL) and disrupts retinoic acid signaling, leading to wide-scale gene expression changes and uncontrolled proliferation of myeloid precursor cells. While known to be recruited to binding sites across the genome, its impact on gene regulation and expression is under-explored. Using integrated multi-omics datasets, we characterize the influence of PML::RARA binding on gene expression and regulation in an inducible PML::RARA cell line model and APL patient ex vivo samples. We find that genes whose regulatory elements recruit PML::RARA are not uniformly transcriptionally repressed, as commonly suggested, but also may be upregulated or remain unchanged. We develop a computational machine learning implementation called Regulatory Element Behavior Extraction Learning to deconvolute the complex, local transcription factor binding site environment at PML::RARA bound positions to reveal distinct signatures that modulate how PML::RARA directs the transcriptional response.

Gene-specific machine learning model to predict the pathogenicity of BRCA2 variants.

Background: Existing BRCA2-specific variant pathogenicity prediction algorithms focus on the prediction of the functional impact of a subtype of variants alone. General variant effect predictors are applicable to all subtypes, but are trained on putative benign and pathogenic variants and do not account for gene-specific information, such as hotspots of pathogenic variants. Local, gene-specific information have been shown to aid variant pathogenicity prediction; therefore, our aim was to develop a BRCA2-specific machine learning model to predict pathogenicity of all types of BRCA2 variants. Methods: We developed an XGBoost-based machine learning model to predict pathogenicity of BRCA2 variants. The model utilizes general variant information such as position, frequency, and consequence for the canonical BRCA2 transcript, as well as deleteriousness prediction scores from several tools. We trained the model on 80% of the expert reviewed variants by the Evidence-Based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) consortium and tested its performance on the remaining 20%, as well as on an independent set of variants of uncertain significance with experimentally determined functional scores. Results: The novel gene-specific model predicted the pathogenicity of ENIGMA BRCA2 variants with an accuracy of 99.9%. The model also performed excellently on predicting the functional consequence of the independent set of variants (accuracy was up to 91.3%). Conclusion: This new, gene-specific model is an accurate method for interpreting the pathogenicity of variants in the BRCA2 gene. It is a valuable addition for variant classification and can prioritize unreviewed variants for functional analysis or expert review.

Comparison of Capture Hi-C Analytical Pipelines.

It is now evident that DNA forms an organized nuclear architecture, which is essential to maintain the structural and functional integrity of the genome. Chromatin organization can be systematically studied due to the recent boom in chromosome conformation capture technologies (e.g., 3C and its successors 4C, 5C and Hi-C), which is accompanied by the development of computational pipelines to identify biologically meaningful chromatin contacts in such data. However, not all tools are applicable to all experimental designs and all structural features. Capture Hi-C (CHi-C) is a method that uses an intermediate hybridization step to target and select predefined regions of interest in a Hi-C library, thereby increasing effective sequencing depth for those regions. It allows researchers to investigate fine chromatin structures at high resolution, for instance promoter-enhancer loops, but it introduces additional biases with the capture step, and therefore requires specialized pipelines. Here, we compare multiple analytical pipelines for CHi-C data analysis. We consider the effect of retaining multi-mapping reads and compare the efficiency of different statistical approaches in both identifying reproducible interactions and determining biologically significant interactions. At restriction fragment level resolution, the number of multi-mapping reads that could be rescued was negligible. The number of identified interactions varied widely, depending on the analytical method, indicating large differences in type I and type II error rates. The optimal pipeline depends on the project-specific tolerance level of false positive and false negative chromatin contacts.

Structure-Based Virtual Screening and Functional Validation of Potential Hit Molecules Targeting the SARS-CoV-2 Main Protease.

The recent global health emergency caused by the coronavirus disease 2019 (COVID-19) pandemic has taken a heavy toll, both in terms of lives and economies. Vaccines against the disease have been developed, but the efficiency of vaccination campaigns worldwide has been variable due to challenges regarding production, logistics, distribution and vaccine hesitancy. Furthermore, vaccines are less effective against new variants of the SARS-CoV-2 virus and vaccination-induced immunity fades over time. These challenges and the vaccines' ineffectiveness for the infected population necessitate improved treatment options, including the inhibition of the SARS-CoV-2 main protease (Mpro). Drug repurposing to achieve inhibition could provide an immediate solution for disease management. Here, we used structure-based virtual screening (SBVS) to identify natural products (from NP-lib) and FDA-approved drugs (from e-Drug3D-lib and Drugs-lib) which bind to the Mpro active site with high-affinity and therefore could be designated as potential inhibitors. We prioritized nine candidate inhibitors (e-Drug3D-lib: Ciclesonide, Losartan and Telmisartan; Drugs-lib: Flezelastine, Hesperidin and Niceverine; NP-lib: three natural products) and predicted their half maximum inhibitory concentration using DeepPurpose, a deep learning tool for drug-target interactions. Finally, we experimentally validated Losartan and two of the natural products as in vitro Mpro inhibitors, using a bioluminescence resonance energy transfer (BRET)-based Mpro sensor. Our study suggests that existing drugs and natural products could be explored for the treatment of COVID-19.

The immune landscape of solid pediatric tumors.

BACKGROUND: Large immunogenomic analyses have demonstrated the prognostic role of the functional orientation of the tumor microenvironment in adult solid tumors, this variable has been poorly explored in the pediatric counterpart. METHODS: We performed a systematic analysis of public RNAseq data (TARGET) for five pediatric tumor types (408 patients): Wilms tumor (WLM), neuroblastoma (NBL), osteosarcoma (OS), clear cell sarcoma of the kidney (CCSK) and rhabdoid tumor of the kidney (RT). We assessed the performance of the Immunologic Constant of Rejection (ICR), which captures an active Th1/cytotoxic response. We also performed gene set enrichment analysis (ssGSEA) and clustered more than 100 well characterized immune traits to define immune subtypes and compared their outcome. RESULTS: A higher ICR score was associated with better survival in OS and high risk NBL without MYCN amplification but with poorer survival in WLM. Clustering of immune traits revealed the same five principal modules previously described in adult tumors (TCGA). These modules divided pediatric patients into six immune subtypes (S1-S6) with distinct survival outcomes. The S2 cluster showed the best overall survival, characterized by low enrichment of the wound healing signature, high Th1, and low Th2 infiltration, while the reverse was observed in S4. Upregulation of the WNT/Beta-catenin pathway was associated with unfavorable outcomes and decreased T-cell infiltration in OS. CONCLUSIONS: We demonstrated that extracranial pediatric tumors could be classified according to their immune disposition, unveiling similarities with adults' tumors. Immunological parameters might be explored to refine diagnostic and prognostic biomarkers and to identify potential immune-responsive tumors.