Integrative identification of non-coding regulatory regions driving metastatic prostate cancer.

Large-scale sequencing efforts of thousands of tumor samples have been undertaken to understand the mutational landscape of the coding genome. However, the vast majority of germline and somatic variants occur within non-coding portions of the genome. These genomic regions do not directly encode for specific proteins, but can play key roles in cancer progression, for example by driving aberrant gene expression control. Here, we designed an integrative computational and experimental framework to identify recurrently mutated non-coding regulatory regions that drive tumor progression. Application of this approach to whole-genome sequencing (WGS) data from a large cohort of metastatic castration-resistant prostate cancer (mCRPC) revealed a large set of recurrently mutated regions. We used (i) in silico prioritization of functional non-coding mutations, (ii) massively parallel reporter assays, and (iii) in vivo CRISPR-interference (CRISPRi) screens in xenografted mice to systematically identify and validate driver regulatory regions that drive mCRPC. We discovered that one of these enhancer regions, GH22I030351, acts on a bidirectional promoter to simultaneously modulate expression of U2-associated splicing factor SF3A1 and chromosomal protein CCDC157. We found that both SF3A1 and CCDC157 are promoters of tumor growth in xenograft models of prostate cancer. We nominated a number of transcription factors, including SOX6, to be responsible for higher expression of SF3A1 and CCDC157. Collectively, we have established and confirmed an integrative computational and experimental approach that enables the systematic detection of non-coding regulatory regions that drive the progression of human cancers.

An mRNA processing pathway suppresses metastasis by governing translational control from the nucleus.

Cancer cells often co-opt post-transcriptional regulatory mechanisms to achieve pathologic expression of gene networks that drive metastasis. Translational control is a major regulatory hub in oncogenesis; however, its effects on cancer progression remain poorly understood. Here, to address this, we used ribosome profiling to compare genome-wide translation efficiencies of poorly and highly metastatic breast cancer cells and patient-derived xenografts. We developed dedicated regression-based methods to analyse ribosome profiling and alternative polyadenylation data, and identified heterogeneous nuclear ribonucleoprotein C (HNRNPC) as a translational controller of a specific mRNA regulon. We found that HNRNPC is downregulated in highly metastatic cells, which causes HNRNPC-bound mRNAs to undergo 3' untranslated region lengthening and, subsequently, translational repression. We showed that modulating HNRNPC expression impacts the metastatic capacity of breast cancer cells in xenograft mouse models. In addition, the reduced expression of HNRNPC and its regulon is associated with the worse prognosis in breast cancer patient cohorts.

A sense-antisense RNA interaction promotes breast cancer metastasis via regulation of NQO1 expression.

Antisense RNAs are ubiquitous in human cells, yet their role is largely unexplored. Here we profiled antisense RNAs in the MDA-MB-231 breast cancer cell line and its highly lung metastatic derivative. We identified one antisense RNA that drives cancer progression by upregulating the redox enzyme NADPH quinone dehydrogenase 1 (NQO1), and named it NQO1-AS. Knockdown of either NQO1 or NQO1-AS reduced lung colonization in a mouse model, and investigation into the role of NQO1 indicated that it is broadly protective against oxidative damage and ferroptosis. Breast cancer cells in the lung are dependent on this pathway, and this dependence can be exploited therapeutically by inducing ferroptosis while inhibiting NQO1. Together, our findings establish a role for NQO1-AS in the progression of breast cancer by regulating its sense mRNA post-transcriptionally. Because breast cancer predominantly affects females, the disease models used in this study are of female origin and the results are primarily applicable to females.

Two isoleucyl tRNAs that decode synonymous codons divergently regulate breast cancer metastatic growth by controlling translation of proliferation-regulating genes.

The human genome contains 61 codons encoding 20 amino acids. Synonymous codons representing a given amino acid are decoded by a set of transfer RNAs (tRNAs) called isoacceptors. We report the surprising observation that two isoacceptor tRNAs that decode synonymous codons become modulated in opposing directions during breast cancer progression. Specifically, tRNAIleUAU became upregulated, whereas tRNAIleGAU became repressed as breast cancer cells attained enhanced metastatic capacity. Functionally, tRNAIleUAU promoted and tRNAIleGAU suppressed metastatic colonization in mouse xenograft models. These tRNAs mediated opposing effects on codon-dependent translation of growth-promoting genes, consistent with genomic enrichment or depletion of their cognate codons in mitotic genes. Our findings uncover a specific isoacceptor tRNA pair that act in opposition, divergently impacting growth-regulating genes and a disease phenotype. Degeneracy of the genetic code can thus be biologically exploited by human cancer cells via tRNA isoacceptor shifts that causally facilitate the transition toward a growth-promoting state.

Global patterns of antigen receptor repertoire disruption across adaptive immune compartments in COVID-19.

Whereas pathogen-specific T and B cells are a primary focus of interest during infectious disease, we have used COVID-19 to ask whether their emergence comes at a cost of broader B cell and T cell repertoire disruption. We applied a genomic DNA-based approach to concurrently study the immunoglobulin-heavy (IGH) and T cell receptor (TCR) ß and δ chain loci of 95 individuals. Our approach detected anticipated repertoire focusing for the IGH repertoire, including expansions of clusters of related sequences temporally aligned with SARS-CoV-2-specific seroconversion, and enrichment of some shared SARS-CoV-2-associated sequences. No significant age-related or disease severity-related deficiencies were noted for the IGH repertoire. By contrast, whereas focusing occurred at the TCRß and TCRδ loci, including some TCRß sequence-sharing, disruptive repertoire narrowing was almost entirely limited to many patients aged older than 50 y. By temporarily reducing T cell diversity and by risking expansions of nonbeneficial T cells, these traits may constitute an age-related risk factor for COVID-19, including a vulnerability to new variants for which T cells may provide key protection.

Highly multiplexed immune repertoire sequencing links multiple lymphocyte classes with severity of response to COVID-19.

Background: Disease progression of subjects with coronavirus disease 2019 (COVID-19) varies dramatically. Understanding the various types of immune response to SARS-CoV-2 is critical for better clinical management of coronavirus outbreaks and to potentially improve future therapies. Disease dynamics can be characterized by deciphering the adaptive immune response. Methods: In this cross-sectional study we analyzed 117 peripheral blood immune repertoires from healthy controls and subjects with mild to severe COVID-19 disease to elucidate the interplay between B and T cells. We used an immune repertoire Primer Extension Target Enrichment method (immunoPETE) to sequence simultaneously human leukocyte antigen (HLA) restricted T cell receptor beta chain (TRB) and unrestricted T cell receptor delta chain (TRD) and immunoglobulin heavy chain (IgH) immune receptor repertoires. The distribution was analyzed of TRB, TRD and IgH clones between healthy and COVID-19 infected subjects. Using McFadden's Adjusted R2 variables were examined for a predictive model. The aim of this study is to analyze the influence of the adaptive immune repertoire on the severity of the disease (value on the World Health Organization Clinical Progression Scale) in COVID-19. Findings: Combining clinical metadata with clonotypes of three immune receptor heavy chains (TRB, TRD, and IgH), we found significant associations between COVID-19 disease severity groups and immune receptor sequences of B and T cell compartments. Logistic regression showed an increase in shared IgH clonal types and decrease of TRD in subjects with severe COVID-19. The probability of finding shared clones of TRD clonal types was highest in healthy subjects (controls). Some specific TRB clones seems to be present in severe COVID-19 (Figure S7b). The most informative models (McFadden´s Adjusted R2=0.141) linked disease severity with immune repertoire measures across all three cell types, as well as receptor-specific cell counts, highlighting the importance of multiple lymphocyte classes in disease progression. Interpretation: Adaptive immune receptor peripheral blood repertoire measures are associated with COVID-19 disease severity. Funding: The study was funded with grants from the Berlin Institute of Health (BIH).

Leucyl-tRNA synthetase is a tumour suppressor in breast cancer and regulates codon-dependent translation dynamics.

Tumourigenesis and cancer progression require enhanced global protein translation1-3. Such enhanced translation is caused by oncogenic and tumour-suppressive events that drive the synthesis and activity of translational machinery4,5. Here we report the surprising observation that leucyl-tRNA synthetase (LARS) becomes repressed during mammary cell transformation and in human breast cancer. Monoallelic genetic deletion of LARS in mouse mammary glands enhanced breast cancer tumour formation and proliferation. LARS repression reduced the abundance of select leucine tRNA isoacceptors, leading to impaired leucine codon-dependent translation of growth suppressive genes, including epithelial membrane protein 3 (EMP3) and gamma-glutamyltransferase 5 (GGT5). Our findings uncover a tumour-suppressive tRNA synthetase and reveal that dynamic repression of a specific tRNA synthetase-along with its downstream cognate tRNAs-elicits a downstream codon-biased translational gene network response that enhances breast tumour formation and growth.

A stress-induced tyrosine-tRNA depletion response mediates codon-based translational repression and growth suppression.

Eukaryotic transfer RNAs can become selectively fragmented upon various stresses, generating tRNA-derived small RNA fragments. Such fragmentation has been reported to impact a small fraction of the tRNA pool and thus presumed to not directly impact translation. We report that oxidative stress can rapidly generate tyrosine-tRNAGUA fragments in human cells-causing significant depletion of the precursor tRNA. Tyrosine-tRNAGUA depletion impaired translation of growth and metabolic genes enriched in cognate tyrosine codons. Depletion of tyrosine tRNAGUA or its translationally regulated targets USP3 and SCD repressed proliferation-revealing a dedicated tRNA-regulated growth-suppressive pathway for oxidative stress response. Tyrosine fragments are generated in a DIS3L2 exoribonuclease-dependent manner and inhibit hnRNPA1-mediated transcript destabilization. Moreover, tyrosine fragmentation is conserved in C. elegans. Thus, tRNA fragmentation can coordinately generate trans-acting small RNAs and functionally deplete a tRNA. Our findings reveal the existence of an underlying adaptive codon-based regulatory response inherent to the genetic code.

Androgen Signaling Regulates SARS-CoV-2 Receptor Levels and Is Associated with Severe COVID-19 Symptoms in Men.

SARS-CoV-2 infection has led to a global health crisis, and yet our understanding of the disease and potential treatment options remains limited. The infection occurs through binding of the virus with angiotensin converting enzyme 2 (ACE2) on the cell membrane. Here, we established a screening strategy to identify drugs that reduce ACE2 levels in human embryonic stem cell (hESC)-derived cardiac cells and lung organoids. Target analysis of hit compounds revealed androgen signaling as a key modulator of ACE2 levels. Treatment with antiandrogenic drugs reduced ACE2 expression and protected hESC-derived lung organoids against SARS-CoV-2 infection. Finally, clinical data on COVID-19 patients demonstrated that prostate diseases, which are linked to elevated androgen, are significant risk factors and that genetic variants that increase androgen levels are associated with higher disease severity. These findings offer insights on the mechanism of disproportionate disease susceptibility in men and identify antiandrogenic drugs as candidate therapeutics for COVID-19.

Androgen Regulates SARS-CoV-2 Receptor Levels and Is Associated with Severe COVID-19 Symptoms in Men.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection has led to a global health crisis, and yet our understanding of the disease pathophysiology and potential treatment options remains limited. SARS-CoV-2 infection occurs through binding and internalization of the viral spike protein to angiotensin converting enzyme 2 (ACE2) on the host cell membrane. Lethal complications are caused by damage and failure of vital organs that express high levels of ACE2, including the lungs, the heart and the kidneys. Here, we established a high-throughput drug screening strategy to identify therapeutic candidates that reduce ACE2 levels in human embryonic stem cell (hESC) derived cardiac cells. Drug target analysis of validated hit compounds, including 5 alpha reductase inhibitors, revealed androgen signaling as a key modulator of ACE2 levels. Treatment with the 5 alpha reductase inhibitor dutasteride reduced ACE2 levels and internalization of recombinant spike receptor binding domain (Spike-RBD) in hESC-derived cardiac cells and human alveolar epithelial cells. Finally, clinical data on coronavirus disease 2019 (COVID-19) patients demonstrated that abnormal androgen states are significantly associated with severe disease complications and cardiac injury as measured by blood troponin T levels. These findings provide important insights on the mechanism of increased disease susceptibility in male COVID-19 patients and identify androgen receptor inhibition as a potential therapeutic strategy.

RBMS1 Suppresses Colon Cancer Metastasis through Targeted Stabilization of Its mRNA Regulon.

Identifying master regulators that drive pathologic gene expression is a key challenge in precision oncology. Here, we have developed an analytic framework, named PRADA, that identifies oncogenic RNA-binding proteins through the systematic detection of coordinated changes in their target regulons. Application of this approach to data collected from clinical samples, patient-derived xenografts, and cell line models of colon cancer metastasis revealed the RNA-binding protein RBMS1 as a suppressor of colon cancer progression. We observed that silencing RBMS1 results in increased metastatic capacity in xenograft mouse models, and that restoring its expression blunts metastatic liver colonization. We have found that RBMS1 functions as a posttranscriptional regulator of RNA stability by directly binding its target mRNAs. Together, our findings establish a role for RBMS1 as a previously unknown regulator of RNA stability and as a suppressor of colon cancer metastasis with clinical utility for risk stratification of patients. SIGNIFICANCE: By applying a new analytic approach to transcriptomic data from clinical samples and models of colon cancer progression, we have identified RBMS1 as a suppressor of metastasis and as a post-transcriptional regulator of RNA stability. Notably, RBMS1 silencing and downregulation of its targets are negatively associated with patient survival.See related commentary by Carter, p. 1261.This article is highlighted in the In This Issue feature, p. 1241.

Toward high-resolution population genomics using archaeological samples.

The term 'ancient DNA' (aDNA) is coming of age, with over 1,200 hits in the PubMed database, beginning in the early 1980s with the studies of 'molecular paleontology'. Rooted in cloning and limited sequencing of DNA from ancient remains during the pre-PCR era, the field has made incredible progress since the introduction of PCR and next-generation sequencing. Over the last decade, aDNA analysis ushered in a new era in genomics and became the method of choice for reconstructing the history of organisms, their biogeography, and migration routes, with applications in evolutionary biology, population genetics, archaeogenetics, paleo-epidemiology, and many other areas. This change was brought by development of new strategies for coping with the challenges in studying aDNA due to damage and fragmentation, scarce samples, significant historical gaps, and limited applicability of population genetics methods. In this review, we describe the state-of-the-art achievements in aDNA studies, with particular focus on human evolution and demographic history. We present the current experimental and theoretical procedures for handling and analysing highly degraded aDNA. We also review the challenges in the rapidly growing field of ancient epigenomics. Advancement of aDNA tools and methods signifies a new era in population genetics and evolutionary medicine research.

Deficiency of the adaptor protein SLy1 results in a natural killer cell ribosomopathy affecting tumor clearance.

Individuals with robust natural killer (NK) cell function incur lower rates of malignancies. To expand our understanding of genetic factors contributing to this phenomenon, we analyzed NK cells from cancer resistant and susceptible strains of mice. We identified a correlation between NK levels of the X-chromosome-located adaptor protein SLy1 and immunologic susceptibility to cancer. Unlike the case for T or B lymphocytes, where SLy1 shuttles between the cytoplasm and nucleus to facilitate signal transduction, in NK cells SLy1 functions as a ribosomal protein and is located solely in the cytoplasm. In its absence, ribosomal instability results in p53-mediated NK cell senescence and decreased clearance of malignancies. NK defects are reversible under inflammatory conditions and viral clearance is not impacted by SLy1 deficiency. Our work defines a previously unappreciated X-linked ribosomopathy that results in a specific and subtle NK cell dysfunction leading to immunologic susceptibility to cancer.

Evolutionary genomics of Culex pipiens: global and local adaptations associated with climate, life-history traits and anthropogenic factors.

We present the first genome-wide study of recent evolution in Culex pipiens species complex focusing on the genomic extent, functional targets and likely causes of global and local adaptations. We resequenced pooled samples of six populations of C. pipiens and two populations of the outgroup Culex torrentium. We used principal component analysis to systematically study differential natural selection across populations and developed a phylogenetic scanning method to analyse admixture without haplotype data. We found evidence for the prominent role of geographical distribution in shaping population structure and specifying patterns of genomic selection. Multiple adaptive events, involving genes implicated with autogeny, diapause and insecticide resistance were limited to specific populations. We estimate that about 5-20% of the genes (including several histone genes) and almost half of the annotated pathways were undergoing selective sweeps in each population. The high occurrence of sweeps in non-genic regions and in chromatin remodelling genes indicated the adaptive importance of gene expression changes. We hypothesize that global adaptive processes in the C. pipiens complex are potentially associated with South to North range expansion, requiring adjustments in chromatin conformation. Strong local signature of adaptation and emergence of hybrid bridge vectors necessitate genomic assessment of populations before specifying control agents.

Wolbachia is not all about sex: male-feminizing Wolbachia alters the leafhopper Zyginidia pullula transcriptome in a mainly sex-independent manner.

Wolbachia causes the feminization of chromosomally male embryos in several species of crustaceans and insects, including the leafhopper Zyginidia pullula. In contrast to the relatively well-established ecological aspects of male feminization (e.g., sex ratio distortion and its consequences), the underlying molecular mechanisms remain understudied and unclear. We embarked on an exploratory study to investigate the extent and nature of Wolbachia's effect on gene expression pattern in Z. pullula. We sequenced whole transcriptomes from Wolbachia-infected and uninfected adults. 18147 loci were assembled de novo, including homologs of several Drosophila sex determination genes. A number of transcripts were flagged as candidate Wolbachia sequences. Despite the resemblance of Wolbachia-infected chromosomal males to uninfected and infected chromosomal females in terms of sexual morphology and behavior, principal component analysis revealed that gene expression patterns did not follow these sexual phenotype categories. The principal components generated by differentially expressed genes specified a strong sex-independent Wolbachia effect, followed by a weaker Wolbachia-sexual karyotype interaction effect. Approaches to further examine the molecular mechanism of Wolbachia-host interactions have been suggested based on the presented findings.

Cytochrome c oxidase subunit 1 barcode data of fish of the Nayband National Park in the Persian Gulf and analysis using meta-data flag several cryptic species.

We provide cytochrome c oxidase subunit 1 (COI) barcode sequences of fishes of the Nayband National Park, Persian Gulf, Iran. Industrial activities, ecological considerations and goals of The Fish Barcode of Life campaign make it crucial that fish species residing in the park be identified. To the best of our knowledge, this is the first report of barcoding data on fishes of the Persian Gulf. We examined 187 individuals representing 76 species, 56 genera and 32 families. The data flagged potentially cryptic species of Gerres filamentosus and Plectorhinchus schotaf. 16S rDNA data on these species are provided. Exclusion of these two potential cryptic species resulted in a mean COI intraspecific distance of 0.18%, and a mean inter- to intraspecific divergence ratio of 66.7. There was no overlap between maximum Kimura 2-parameter distances among conspecifics (1.66%) and minimum distance among congeneric species (6.19%). Barcodes shared among species were not observed. Neighbour-joining analysis showed that most species formed cohesive sequence units with little variation. Finally, the comparison of 16 selected species from this study with meta-data of conspecifics from Australia, India, China and South Africa revealed high interregion divergences and potential existence of six cryptic species. Pairwise interregional comparisons were more informative than global divergence assessments with regard to detection of cryptic variation. Our analysis exemplifies optimal use of the expanding barcode data now becoming available.