Integrated multi-omics analysis of zinc-finger proteins uncovers roles in RNA regulation.

RNA interactome studies have revealed that hundreds of zinc-finger proteins (ZFPs) are candidate RNA-binding proteins (RBPs), yet their RNA substrates and functional significance remain largely uncharacterized. Here, we present a systematic multi-omics analysis of the DNA- and RNA-binding targets and regulatory roles of more than 100 ZFPs representing 37 zinc-finger families. We show that multiple ZFPs are previously unknown regulators of RNA splicing, alternative polyadenylation, stability, or translation. The examined ZFPs show widespread sequence-specific RNA binding and preferentially bind proximal to transcription start sites. Additionally, several ZFPs associate with their targets at both the DNA and RNA levels. We highlight ZNF277, a C2H2 ZFP that binds thousands of RNA targets and acts as a multi-functional RBP. We also show that ZNF473 is a DNA/RNA-associated protein that regulates the expression and splicing of cell cycle genes. Our results reveal diverse roles for ZFPs in transcriptional and post-transcriptional gene regulation.

RANK ligand converts the NCoR/HDAC3 co-repressor to a PGC1ß- and RNA-dependent co-activator of osteoclast gene expression.

The nuclear receptor co-repressor (NCoR) complex mediates transcriptional repression dependent on histone deacetylation by histone deacetylase 3 (HDAC3) as a component of the complex. Unexpectedly, we found that signaling by the receptor activator of nuclear factor κB (RANK) converts the NCoR/HDAC3 co-repressor complex to a co-activator of AP-1 and NF-κB target genes that are required for mouse osteoclast differentiation. Accordingly, the dominant function of NCoR/HDAC3 complexes in response to RANK signaling is to activate, rather than repress, gene expression. Mechanistically, RANK signaling promotes RNA-dependent interaction of the transcriptional co-activator PGC1ß with the NCoR/HDAC3 complex, resulting in the activation of PGC1ß and inhibition of HDAC3 activity for acetylated histone H3. Non-coding RNAs Dancr and Rnu12, which are associated with altered human bone homeostasis, promote NCoR/HDAC3 complex assembly and are necessary for RANKL-induced osteoclast differentiation in vitro. These findings may be prototypic for signal-dependent functions of NCoR in other biological contexts.

Exploiting dynamic enhancer landscapes to decode macrophage and microglia phenotypes in health and disease.

The development and functional potential of metazoan cells is dependent on combinatorial roles of transcriptional enhancers and promoters. Macrophages provide exceptionally powerful model systems for investigation of mechanisms underlying the activation of cell-specific enhancers that drive transitions in cell fate and cell state. Here, we review recent advances that have expanded appreciation of the diversity of macrophage phenotypes in health and disease, emphasizing studies of liver, adipose tissue, and brain macrophages as paradigms for other tissue macrophages and cell types. Studies of normal tissue-resident macrophages and macrophages associated with cirrhosis, obese adipose tissue, and neurodegenerative disease illustrate the major roles of tissue environment in remodeling enhancer landscapes to specify the development and functions of distinct macrophage phenotypes. We discuss the utility of quantitative analysis of environment-dependent changes in enhancer activity states as an approach to discovery of regulatory transcription factors and upstream signaling pathways.

A TLR4/TRAF6-dependent signaling pathway mediates NCoR coactivator complex formation for inflammatory gene activation.

The nuclear receptor corepressor (NCoR) forms a complex with histone deacetylase 3 (HDAC3) that mediates repressive functions of unliganded nuclear receptors and other transcriptional repressors by deacetylation of histone substrates. Recent studies provide evidence that NCoR/HDAC3 complexes can also exert coactivator functions in brown adipocytes by deacetylating and activating PPARγ coactivator 1α (PGC1α) and that signaling via receptor activator of nuclear factor kappa-B (RANK) promotes the formation of a stable NCoR/HDAC3/PGC1ß complex that coactivates nuclear factor kappa-B (NFκB)- and activator protein 1 (AP-1)-dependent genes required for osteoclast differentiation. Here, we demonstrate that activation of Toll-like receptor (TLR) 4, but not TLR3, the interleukin 4 (IL4) receptor nor the Type I interferon receptor, also promotes assembly of an NCoR/HDAC3/PGC1ß coactivator complex. Receptor-specific utilization of TNF receptor-associated factor 6 (TRAF6) and downstream activation of extracellular signal-regulated kinase 1 (ERK1) and TANK-binding kinase 1 (TBK1) accounts for the common ability of RANK and TLR4 to drive assembly of an NCoR/HDAC3/PGC1ß complex in macrophages. ERK1, the p65 component of NFκB, and the p300 histone acetyltransferase (HAT) are also components of the induced complex and are associated with local histone acetylation and transcriptional activation of TLR4-dependent enhancers and promoters. These observations identify a TLR4/TRAF6-dependent signaling pathway that converts NCoR from a corepressor of nuclear receptors to a coactivator of NFκB and AP-1 that may be relevant to functions of NCoR in other developmental and homeostatic processes.

Robust single-cell discovery of RNA targets of RNA-binding proteins and ribosomes.

RNA-binding proteins (RBPs) are critical regulators of gene expression and RNA processing that are required for gene function. Yet the dynamics of RBP regulation in single cells is unknown. To address this gap in understanding, we developed STAMP (Surveying Targets by APOBEC-Mediated Profiling), which efficiently detects RBP-RNA interactions. STAMP does not rely on ultraviolet cross-linking or immunoprecipitation and, when coupled with single-cell capture, can identify RBP-specific and cell-type-specific RNA-protein interactions for multiple RBPs and cell types in single, pooled experiments. Pairing STAMP with long-read sequencing yields RBP target sites in an isoform-specific manner. Finally, Ribo-STAMP leverages small ribosomal subunits to measure transcriptome-wide ribosome association in single cells. STAMP enables the study of RBP-RNA interactomes and translational landscapes with unprecedented cellular resolution.

Nuclear and cytoplasmic poly(A) binding proteins (PABPs) favor distinct transcripts and isoforms.

The poly(A)-tail appended to the 3'-end of most eukaryotic transcripts plays a key role in their stability, nuclear transport, and translation. These roles are largely mediated by Poly(A) Binding Proteins (PABPs) that coat poly(A)-tails and interact with various proteins involved in the biogenesis and function of RNA. While it is well-established that the nuclear PABP (PABPN) binds newly synthesized poly(A)-tails and is replaced by the cytoplasmic PABP (PABPC) on transcripts exported to the cytoplasm, the distribution of transcripts for different genes or isoforms of the same gene on these PABPs has not been investigated on a genome-wide scale. Here, we analyzed the identity, splicing status, poly(A)-tail size, and translation status of RNAs co-immunoprecipitated with endogenous PABPN or PABPC in human cells. At steady state, many protein-coding and non-coding RNAs exhibit strong bias for association with PABPN or PABPC. While PABPN-enriched transcripts more often were incompletely spliced and harbored longer poly(A)-tails and PABPC-enriched RNAs had longer half-lives and higher translation efficiency, there are curious outliers. Overall, our study reveals the landscape of RNAs bound by PABPN and PABPC, providing new details that support and advance the current understanding of the roles these proteins play in poly(A)-tail synthesis, maintenance, and function.

FLARE: a fast and flexible workflow for identifying RNA editing foci.

BACKGROUND: Fusion of RNA-binding proteins (RBPs) to RNA base-editing enzymes (such as APOBEC1 or ADAR) has emerged as a powerful tool for the discovery of RBP binding sites. However, current methods that analyze sequencing data from RNA-base editing experiments are vulnerable to false positives due to off-target editing, genetic variation and sequencing errors. RESULTS: We present FLagging Areas of RNA-editing Enrichment (FLARE), a Snakemake-based pipeline that builds on the outputs of the SAILOR edit site discovery tool to identify regions statistically enriched for RNA editing. FLARE can be configured to analyze any type of RNA editing, including C to U and A to I. We applied FLARE to C-to-U editing data from a RBFOX2-APOBEC1 STAMP experiment, to show that our approach attains high specificity for detecting RBFOX2 binding sites. We also applied FLARE to detect regions of exogenously introduced as well as endogenous A-to-I editing. CONCLUSIONS: FLARE is a fast and flexible workflow that identifies significantly edited regions from RNA-seq data. The FLARE codebase is available at https://github.com/YeoLab/FLARE .

Evaluating the molecular diagnostic yield of joint genotyping-based approach for detecting rare germline pathogenic and putative loss-of-function variants.

PURPOSE: Cohort-based germline variant characterization is the standard approach for pathogenic variant discovery in clinical and research samples. However, the impact of cohort size on the molecular diagnostic yield of joint genotyping is largely unknown. METHODS: Head-to-head comparison of the molecular diagnostic yield of joint genotyping in two cohorts of 239 cancer patients in the absence and then in the presence of 100 additional germline exomes. RESULTS: In 239 testicular cancer patients, 4 (7.4%, 95% confidence interval [CI]: 2.1-17.9) of 54 pathogenic variants in the cancer predisposition and American College of Medical Genetics and Genomics (ACMG) genes were missed by one or both computational runs of joint genotyping. Similarly, 8 (12.1%, 95% CI: 5.4-22.5) of 66 pathogenic variants in these genes were undetected by joint genotyping in another independent cohort of 239 breast cancer patients. An exome-wide analysis of putative loss-of-function (pLOF) variants in the testicular cancer cohort showed that 162 (8.2%, 95% CI: 7.1-9.6) pLOF variants were only detected in one analysis run but not the other, while 433 (22.0%, 95% CI: 20.2-23.9%) pLOF variants were filtered out by both analyses despite having sufficient sequencing coverage. CONCLUSION: Our analysis of the standard germline variant detection method highlighted a substantial impact of concurrently analyzing additional genomic data sets on the ability to detect clinically relevant germline pathogenic variants.

Detection of Pathogenic Variants With Germline Genetic Testing Using Deep Learning vs Standard Methods in Patients With Prostate Cancer and Melanoma.

Importance: Less than 10% of patients with cancer have detectable pathogenic germline alterations, which may be partially due to incomplete pathogenic variant detection. Objective: To evaluate if deep learning approaches identify more germline pathogenic variants in patients with cancer. Design, Setting, and Participants: A cross-sectional study of a standard germline detection method and a deep learning method in 2 convenience cohorts with prostate cancer and melanoma enrolled in the US and Europe between 2010 and 2017. The final date of clinical data collection was December 2017. Exposures: Germline variant detection using standard or deep learning methods. Main Outcomes and Measures: The primary outcomes included pathogenic variant detection performance in 118 cancer-predisposition genes estimated as sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV). The secondary outcomes were pathogenic variant detection performance in 59 genes deemed actionable by the American College of Medical Genetics and Genomics (ACMG) and 5197 clinically relevant mendelian genes. True sensitivity and true specificity could not be calculated due to lack of a criterion reference standard, but were estimated as the proportion of true-positive variants and true-negative variants, respectively, identified by each method in a reference variant set that consisted of all variants judged to be valid from either approach. Results: The prostate cancer cohort included 1072 men (mean [SD] age at diagnosis, 63.7 [7.9] years; 857 [79.9%] with European ancestry) and the melanoma cohort included 1295 patients (mean [SD] age at diagnosis, 59.8 [15.6] years; 488 [37.7%] women; 1060 [81.9%] with European ancestry). The deep learning method identified more patients with pathogenic variants in cancer-predisposition genes than the standard method (prostate cancer: 198 vs 182; melanoma: 93 vs 74); sensitivity (prostate cancer: 94.7% vs 87.1% [difference, 7.6%; 95% CI, 2.2% to 13.1%]; melanoma: 74.4% vs 59.2% [difference, 15.2%; 95% CI, 3.7% to 26.7%]), specificity (prostate cancer: 64.0% vs 36.0% [difference, 28.0%; 95% CI, 1.4% to 54.6%]; melanoma: 63.4% vs 36.6% [difference, 26.8%; 95% CI, 17.6% to 35.9%]), PPV (prostate cancer: 95.7% vs 91.9% [difference, 3.8%; 95% CI, -1.0% to 8.4%]; melanoma: 54.4% vs 35.4% [difference, 19.0%; 95% CI, 9.1% to 28.9%]), and NPV (prostate cancer: 59.3% vs 25.0% [difference, 34.3%; 95% CI, 10.9% to 57.6%]; melanoma: 80.8% vs 60.5% [difference, 20.3%; 95% CI, 10.0% to 30.7%]). For the ACMG genes, the sensitivity of the 2 methods was not significantly different in the prostate cancer cohort (94.9% vs 90.6% [difference, 4.3%; 95% CI, -2.3% to 10.9%]), but the deep learning method had a higher sensitivity in the melanoma cohort (71.6% vs 53.7% [difference, 17.9%; 95% CI, 1.82% to 34.0%]). The deep learning method had higher sensitivity in the mendelian genes (prostate cancer: 99.7% vs 95.1% [difference, 4.6%; 95% CI, 3.0% to 6.3%]; melanoma: 91.7% vs 86.2% [difference, 5.5%; 95% CI, 2.2% to 8.8%]). Conclusions and Relevance: Among a convenience sample of 2 independent cohorts of patients with prostate cancer and melanoma, germline genetic testing using deep learning, compared with the current standard genetic testing method, was associated with higher sensitivity and specificity for detection of pathogenic variants. Further research is needed to understand the relevance of these findings with regard to clinical outcomes.

High-sensitivity in situ capture of endogenous RNA-protein interactions in fixed cells and primary tissues.

RNA-binding proteins (RBPs) have pivotal functions in RNA metabolism, but current methods are limited in retrieving RBP-RNA interactions within endogenous biological contexts. Here, we develop INSCRIBE (IN situ Sensitive Capture of RNA-protein Interactions in Biological Environments), circumventing the challenges through in situ RNA labeling by precisely directing a purified APOBEC1-nanobody fusion to the RBP of interest. This method enables highly specific RNA-binding site identification across a diverse range of fixed biological samples such as HEK293T cells and mouse brain tissue and accurately identifies the canonical binding motifs of RBFOX2 (UGCAUG) and TDP-43 (UGUGUG) in native cellular environments. Applicable to any RBP with available primary antibodies, INSCRIBE enables sensitive capture of RBP-RNA interactions from ultra-low input equivalent to ~5 cells. The robust, versatile, and sensitive INSCRIBE workflow is particularly beneficial for precious tissues such as clinical samples, empowering the exploration of genuine RBP-RNA interactions in RNA-related disease contexts.

Accelerated Aging and Microsatellite Instability in Recessive Dystrophic Epidermolysis Bullosa-Associated Cutaneous Squamous Cell Carcinoma.

Recessive dystrophic epidermolysis bullosa (RDEB) is a severely debilitating disorder caused by pathogenic variants in COL7A1 and is characterized by extreme skin fragility, chronic inflammation, and fibrosis. A majority of patients with RDEB develop squamous cell carcinoma, a highly aggressive skin cancer with limited treatment options currently available. In this study, we utilized an approach leveraging whole-genome sequencing and RNA sequencing across 3 different tissues in a single patient with RDEB to gain insight into possible mechanisms of RDEB-associated squamous cell carcinoma progression and to identify potential therapeutic options. As a result, we identified PLK-1 as a possible candidate for targeted therapy and discovered microsatellite instability and accelerated aging as factors potentially contributing to the aggressive nature and early onset of RDEB squamous cell carcinoma. By integrating multitissue genomic and transcriptomic analyses in a single patient, we demonstrate the promise of bridging the gap between genomic research and clinical applications for developing tailored therapies for patients with rare genetic disorders such as RDEB.

Integrative Analysis of Germline Rare Variants in Clear and Non-clear Cell Renal Cell Carcinoma.

Background and objective: Previous germline studies on renal cell carcinoma (RCC) have usually pooled clear and non-clear cell RCCs and have not adequately accounted for population stratification, which might have led to an inaccurate estimation of genetic risk. Here, we aim to analyze the major germline drivers of RCC risk and clinically relevant but underexplored germline variant types. Methods: We first characterized germline pathogenic variants (PVs), cryptic splice variants, and copy number variants (CNVs) in 1436 unselected RCC patients. To evaluate the enrichment of PVs in RCC, we conducted a case-control study of 1356 RCC patients ancestry matched with 16 512 cancer-free controls using approaches accounting for population stratification and histological subtypes, followed by characterization of secondary somatic events. Key findings and limitations: Clear cell RCC patients (n = 976) exhibited a significant burden of PVs in VHL compared with controls (odds ratio [OR]: 39.1, p = 4.95e-05). Non-clear cell RCC patients (n = 380) carried enrichment of PVs in FH (OR: 77.9, p = 1.55e-08) and MET (OR: 1.98e11, p = 2.07e-05). In a CHEK2-focused analysis with European participants, clear cell RCC (n = 906) harbored nominal enrichment of low-penetrance CHEK2 variants-p.Ile157Thr (OR: 1.84, p = 0.049) and p.Ser428Phe (OR: 5.20, p = 0.045), while non-clear cell RCC (n = 295) exhibited nominal enrichment of CHEK2 loss of function PVs (OR: 3.51, p = 0.033). Patients with germline PVs in FH, MET, and VHL exhibited significantly earlier age of cancer onset than patients without germline PVs (mean: 46.0 vs 60.2 yr, p < 0.0001), and more than half had secondary somatic events affecting the same gene (n = 10/15, 66.7%). Conversely, CHEK2 PV carriers exhibited a similar age of onset to patients without germline PVs (mean: 60.1 vs 60.2 yr, p = 0.99), and only 30.4% carried somatic events in CHEK2 (n = 7/23). Finally, pathogenic germline cryptic splice variants were identified in SDHA and TSC1, and pathogenic germline CNVs were found in 18 patients, including CNVs in FH, SDHA, and VHL. Conclusions and clinical implications: This analysis supports the existing link between several RCC risk genes and RCC risk manifesting in earlier age of onset. It calls for caution when assessing the role of CHEK2 due to the burden of founder variants with varying population frequency. It also broadens the definition of the RCC germline landscape of pathogenicity to incorporate previously understudied types of germline variants. Patient summary: In this study, we carefully compared the frequency of rare inherited mutations with a focus on patients' genetic ancestry. We discovered that subtle variations in genetic background may confound a case-control analysis, especially in evaluating the cancer risk associated with specific genes, such as CHEK2. We also identified previously less explored forms of rare inherited mutations, which could potentially increase the risk of kidney cancer.

Expanded palette of RNA base editors for comprehensive RBP-RNA interactome studies.

RNA binding proteins (RBPs) are key regulators of RNA processing and cellular function. Technologies to discover RNA targets of RBPs such as TRIBE (targets of RNA binding proteins identified by editing) and STAMP (surveying targets by APOBEC1 mediated profiling) utilize fusions of RNA base-editors (rBEs) to RBPs to circumvent the limitations of immunoprecipitation (CLIP)-based methods that require enzymatic digestion and large amounts of input material. To broaden the repertoire of rBEs suitable for editing-based RBP-RNA interaction studies, we have devised experimental and computational assays in a framework called PRINTER (protein-RNA interaction-based triaging of enzymes that edit RNA) to assess over thirty A-to-I and C-to-U rBEs, allowing us to identify rBEs that expand the characterization of binding patterns for both sequence-specific and broad-binding RBPs. We also propose specific rBEs suitable for dual-RBP applications. We show that the choice between single or multiple rBEs to fuse with a given RBP or pair of RBPs hinges on the editing biases of the rBEs and the binding preferences of the RBPs themselves. We believe our study streamlines and enhances the selection of rBEs for the next generation of RBP-RNA target discovery.

Expanded palette of RNA base editors for comprehensive RBP-RNA interactome studies.

RNA binding proteins (RBPs) are key regulators of RNA processing and cellular function. Technologies to discover RNA targets of RBPs such as TRIBE (targets of RNA binding proteins identified by editing) and STAMP (surveying targets by APOBEC1 mediated profiling) utilize fusions of RNA base-editors (rBEs) to RBPs to circumvent the limitations of immunoprecipitation (CLIP)-based methods that require enzymatic digestion and large amounts of input material. To broaden the repertoire of rBEs suitable for editing-based RBP-RNA interaction studies, we have devised experimental and computational assays in a framework called PRINTER (protein-RNA interaction-based triaging of enzymes that edit RNA) to assess over thirty A-to-I and C-to-U rBEs, allowing us to identify rBEs that expand the characterization of binding patterns for both sequence-specific and broad-binding RBPs. We also propose specific rBEs suitable for dual-RBP applications. We show that the choice between single or multiple rBEs to fuse with a given RBP or pair of RBPs hinges on the editing biases of the rBEs and the binding preferences of the RBPs themselves. We believe our study streamlines and enhances the selection of rBEs for the next generation of RBP-RNA target discovery.

Integrative Analysis of Germline Rare Variants in Clear and Non-Clear Cell Renal Cell Carcinoma.

IMPORTANCE: RCC encompasses a set of histologically distinct cancers with a high estimated genetic heritability, of which only a portion is currently explained. Previous rare germline variant studies in RCC have usually pooled clear and non-clear cell RCCs and have not adequately accounted for population stratification that may significantly impact the interpretation and discovery of certain candidate risk genes. OBJECTIVE: To evaluate the enrichment of germline PVs in established cancer-predisposing genes (CPGs) in clear cell and non-clear cell RCC patients compared to cancer-free controls using approaches that account for population stratification and to identify unconventional types of germline RCC risk variants that confer an increased risk of developing RCC. DESIGN SETTING AND PARTICIPANTS: In 1,436 unselected RCC patients with sufficient data quality, we systematically identified rare germline PVs, cryptic splice variants, and copy number variants (CNVs). From this unselected cohort, 1,356 patients were ancestry-matched with 16,512 cancer-free controls, and gene-level enrichment of rare germline PVs were assessed in 143 CPGs, followed by an investigation of somatic events in matching tumor samples. MAIN OUTCOMES AND MEASURES: Gene-level burden of rare germline PVs, identification of secondary somatic events accompanying the germline PVs, and characterization of less-explored types of rare germline PVs in RCC patients. RESULTS: In clear cell RCC (n = 976 patients), patients exhibited significantly higher prevalence of PVs in VHL compared to controls (OR: 39.1, 95% CI: 7.01-218.07, p-value:4.95e-05, q-value:0.00584). In non-clear cell RCC (n = 380 patients), patients carried enriched burden of PVs in FH (OR: 77.9, 95% CI: 18.68-324.97, p-value:1.55e-08, q-value: 1.83e-06) and MET (OR: 1.98e11, 95% CI: 0-inf, p-value: 2.07e-05, q-value: 3.50e-07). In a CHEK2-focused analysis with European cases and controls, clear cell RCC patients (n=906 European patients) harbored nominal enrichment of the previously reported low-penetrance CHEK2 variants, p.Ile157Thr (OR:1.84, 95% CI: 1.00-3.36, p-value:0.049) and p.Ser428Phe (OR:5.20, 95% CI: 1.00-26.40, p-value:0.045) while non-clear cell RCC patients (n=295 European patients) exhibited nominal enrichment of CHEK2 LOF germline PVs (OR: 3.51, 95% CI: 1.10-11.10, p-value: 0.033). RCC patients with germline PVs in FH, MET, and VHL exhibited significantly earlier age of cancer onset compared to patients without any germline PVs in CPGs (Mean: 46.0 vs 60.2 years old, Tukey adjusted p-value < 0.0001), and more than half had secondary somatic events affecting the same gene (n=10/15, 66.7%, 95% CI: 38.7-87.0%). Conversely, patients with rare germline PVs in CHEK2 exhibited a similar age of disease onset to patients without any identified germline PVs in CPGs (Mean: 60.1 vs 60.2 years old, Tukey adjusted p-value: 0.99), and only 30.4% of the patients carried secondary somatic events in CHEK2 (n=7/23, 95% CI: 14.1-53.0%). Finally, rare pathogenic germline cryptic splice variants underexplored in RCC were identified in SDHA and TSC1, and rare pathogenic germline CNVs were found in 18 patients, including CNVs in FH, SDHA, and VHL. CONCLUSIONS AND RELEVANCE: This systematic analysis supports the existing link between several RCC risk genes and elevated RCC risk manifesting in earlier age of RCC onset. Our analysis calls for caution when assessing the role of germline PVs in CHEK2 due to the burden of founder variants with varying population frequency in different ancestry groups. It also broadens the definition of the RCC germline landscape of pathogenicity to incorporate previously understudied types of germline variants, such as cryptic splice variants and CNVs.

Discovery and functional interrogation of SARS-CoV-2 protein-RNA interactions.

The COVID-19 pandemic is caused by severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2). The betacoronvirus has a positive sense RNA genome which encodes for several RNA binding proteins. Here, we use enhanced crosslinking and immunoprecipitation to investigate SARS-CoV-2 protein interactions with viral and host RNAs in authentic virus-infected cells. SARS-CoV-2 proteins, NSP8, NSP12, and nucleocapsid display distinct preferences to specific regions in the RNA viral genome, providing evidence for their shared and separate roles in replication, transcription, and viral packaging. SARS-CoV-2 proteins expressed in human lung epithelial cells bind to 4773 unique host coding RNAs. Nine SARS-CoV-2 proteins upregulate target gene expression, including NSP12 and ORF9c, whose RNA substrates are associated with pathways in protein N-linked glycosylation ER processing and mitochondrial processes. Furthermore, siRNA knockdown of host genes targeted by viral proteins in human lung organoid cells identify potential antiviral host targets across different SARS-CoV-2 variants. Conversely, NSP9 inhibits host gene expression by blocking mRNA export and dampens cytokine productions, including interleukin-1α/ß. Our viral protein-RNA interactome provides a catalog of potential therapeutic targets and offers insight into the etiology of COVID-19 as a safeguard against future pandemics.

Discovery and functional interrogation of SARS-CoV-2 protein-RNA interactions.

The COVID-19 pandemic is caused by severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2). The betacoronvirus has a positive sense RNA genome which encodes for several RNA binding proteins. Here, we use enhanced crosslinking and immunoprecipitation to investigate SARS-CoV-2 protein interactions with viral and host RNAs in authentic virus-infected cells. SARS-CoV-2 proteins, NSP8, NSP12, and nucleocapsid display distinct preferences to specific regions in the RNA viral genome, providing evidence for their shared and separate roles in replication, transcription, and viral packaging. SARS-CoV-2 proteins expressed in human lung epithelial cells bind to 4773 unique host coding RNAs. Nine SARS-CoV-2 proteins upregulate target gene expression, including NSP12 and ORF9c, whose RNA substrates are associated with pathways in protein N-linked glycosylation ER processing and mitochondrial processes. Furthermore, siRNA knockdown of host genes targeted by viral proteins in human lung organoid cells identify potential antiviral host targets across different SARS-CoV-2 variants. Conversely, NSP9 inhibits host gene expression by blocking mRNA export and dampens cytokine productions, including interleukin-1α/ß. Our viral protein-RNA interactome provides a catalog of potential therapeutic targets and offers insight into the etiology of COVID-19 as a safeguard against future pandemics.

Integrating molecular profiles into clinical frameworks through the Molecular Oncology Almanac to prospectively guide precision oncology.

Tumor molecular profiling of single gene-variant ('first-order') genomic alterations informs potential therapeutic approaches. Interactions between such first-order events and global molecular features (for example, mutational signatures) are increasingly associated with clinical outcomes, but these 'second-order' alterations are not yet accounted for in clinical interpretation algorithms and knowledge bases. We introduce the Molecular Oncology Almanac (MOAlmanac), a paired clinical interpretation algorithm and knowledge base to enable integrative interpretation of multimodal genomic data for point-of-care decision making and translational-hypothesis generation. We benchmarked MOAlmanac to a first-order interpretation method across multiple retrospective cohorts and observed an increased number of clinical hypotheses from evaluation of molecular features and profile-to-cell line matchmaking. When applied to a prospective precision oncology trial cohort, MOAlmanac nominated a median of two therapies per patient and identified therapeutic strategies administered in 47% of patients. Overall, we present an open-source computational method for integrative clinical interpretation of individualized molecular profiles.

Genomics of response to immune checkpoint therapies for cancer: implications for precision medicine.

Immune checkpoint blockade (ICB) therapies, which potentiate the body's natural immune response against tumor cells, have shown immense promise in the treatment of various cancers. Currently, tumor mutational burden (TMB) and programmed death ligand 1 (PD-L1) expression are the primary biomarkers evaluated for clinical management of cancer patients across histologies. However, the wide range of responses has demonstrated that the specific molecular and genetic characteristics of each patient's tumor and immune system must be considered to maximize treatment efficacy. Here, we review the various biological pathways and emerging biomarkers implicated in response to PD-(L)1 and cytotoxic T lymphocyte-associated antigen 4 (CTLA-4) therapies, including oncogenic signaling pathways, human leukocyte antigen (HLA) variability, mutation and neoantigen burden, microbiome composition, endogenous retroviruses (ERV), and deficiencies in chromatin remodeling and DNA damage repair (DDR) machinery. We also discuss several mechanisms that have been observed to confer resistance to ICB, such as loss of phosphatase and tensin homolog (PTEN), loss of major histocompatibility complex (MHC) I/II expression, and activation of the indoleamine 2,3-dioxygenase 1 (IDO1) and transforming growth factor beta (TGFß) pathways. Clinical trials testing the combination of PD-(L)1 or CTLA-4 blockade with molecular mediators of these pathways are becoming more common and may hold promise for improving treatment efficacy and response. Ultimately, some of the genes and molecular mechanisms highlighted in this review may serve as novel biological targets or therapeutic vulnerabilities to improve clinical outcomes in patients.

Integrated genomic characterization of oral carcinomas in post-hematopoietic stem cell transplantation survivors.

OBJECTIVES: Secondary oral squamous cell carcinoma (OSCC) is a late complication in allogeneic hematopoietic stem cell transplantation (HSCT) patients, but little is known about long-term outcomes and prognostication. Additionally, molecular alterations and immunologic insights unique to this disease remain largely unexplored. METHODS: We present a cohort of 31 patients with post-HSCT OSCC and reported on clinicopathologic predictors of survival. Whole-exome sequencing was performed on 6 (19%) matched pairs of peripheral blood (post-conditioning, pre-HSCT) and tumor samples. The entire cohort had archival tumor available for immunoprofiling with PD-1/L1 immunohistochemistry. RESULTS: Five-year overall survival (OS) was 57% (95% CI: 46.1-69.8) with a median disease-free survival (DFS) of 13.3 months. Advanced initial staging, a buccal or oral tongue subsite, chronic oral graft-versus-host disease (GVHD) and smoking all negatively impacted survival. High tumor mutational burden (TMB) (median 11.3 vs. 5.0) and unique mutational signatures were noted between unrelated and related donor groups - with a strong correlation between infiltrating PD-1+ lymphocytes and TMB (R = 0.98, p < 0.01). Some differences were observed when comparing commonly mutated genes among our cohort and TCGA, with a predominance of TP53 events. CONCLUSION: Survival outcomes appear similar in HSCT survivors with OSCC compared with non-HSCT OSCC populations. We identified somatic alterations in genes with therapeutic potential unique to this subpopulation of oral cancers.