Long-read RNA sequencing identifies polyadenylation elongation and differential transcript usage of host transcripts during SARS-CoV-2 in vitro infection

Better methods to interrogate host-pathogen interactions during Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infections are imperative to help understand and prevent this disease. Here we implemented RNA-sequencing (RNA-seq) combined with the Oxford Nanopore Technologies (ONT) long-reads to measure differential host gene expression, transcript polyadenylation and isoform usage within various epithelial cell lines permissive and non-permissive for SARS-CoV-2 infection. SARS-CoV-2-infected and mock-infected Vero (African green monkey kidney epithelial cells), Calu-3 (human lung adenocarcinoma epithelial cells), Caco-2 (human colorectal adenocarcinoma epithelial cells) and A549 (human lung carcinoma epithelial cells) were analysed over time (0, 2, 24, 48 hours). Differential polyadenylation was found to occur in both infected Calu-3 and Vero cells during a late time point (48 hpi), with Gene Ontology (GO) terms such as viral transcription and translation shown to be significantly enriched in Calu-3 data. Poly(A) tails showed increased lengths in the majority of the differentially polyadenylated transcripts in Calu-3 and Vero cell lines (up to ~136 nt in mean poly(A) length, padj = 0.029). Of these genes, ribosomal protein genes such as RPS4X and RPS6 also showed downregulation in expression levels, suggesting the importance of ribosomal protein genes during infection. Furthermore, differential transcript usage was identified in Caco-2, Calu-3 and Vero cells, including transcripts of genes such as GSDMB and KPNA2, which have previously been implicated in SARS-CoV-2 infections. Overall, these results highlight the potential role of differential polyadenylation and transcript usage in host immune response or viral manipulation of host mechanisms during infection, and therefore, showcase the value of long-read sequencing in identifying less-explored host responses to disease.

Diamond-Blackfan anemia

Diamond-Blackfan anemia (DBA) is an inherited bone marrow failure syndrome, characterized as a rare congenital bone marrow erythroid hypoplasia (OMIM#105650). Erythroid defect in DBA results in erythroblastopenia in bone marrow as a consequence of maturation blockade between the burst forming unit–erythroid and colony forming unit–erythroid developmental stages, leading to moderate to severe usually macrocytic aregenerative (<20 × 109/L of reticulocytes) anemia. Congenital malformations localized mostly in the cephalic area and in the extremities (thumbs), as well as short stature and cardiac and urogenital tract abnormalities, are a feature of 50% of the DBA-affected patients. A significant increased risk for malignancy has been reported. DBA is due to a defect in the ribosomal RNA (rRNA) maturation as a consequence of a heterozygous mutation in 1 of the 20 ribosomal protein genes. Besides classical DBA, some DBA-like diseases have been identified. The relation between the defect in rRNA maturation and the erythroid defect in DBA has yet to be fully defined. However, recent studies have identified a role for GATA1 either due to a specific defect in its translation or due to its defective regulation by its chaperone HSP70. In addition, excess free heme-induced reactive oxygen species and apoptosis have been implicated in the DBA erythroid phenotype. Current treatment options are either regular transfusions with appropriate iron chelation or treatment with corticosteroids starting at 1 year of age. The only curative treatment for the anemia of DBA to date is bone marrow transplantation. Use of gene therapy as a therapeutic strategy is currently being explored.

Recruitment of MLL1 Complex Is Essential for SETBP1 to Induce Myeloid Transformation

Abnormal activation of SETBP1 due to overexpression or missense mutations occurs frequently in various myeloid neoplasms and associates with poor prognosis. Direct activation of Hoxa9/Hoxa10/Myb transcription by SETBP1 and its missense mutants is essential for their transforming capability; however, the underlying mechanisms for such activation remain elusive. We found that knockdown of Mll1 in mouse myeloid progenitors immortalized by SETBP1 or its missense mutant SETBP1(D/N) caused significant reduction in the mRNA levels of Hoxa9/Hoxa10/Myb, suggesting that Mll1 is critical for their transcriptional activation induced by SETBP1 and its missense mutants. Physical association of MLL1 with SETBP1/SETBP1(D/N) was readily detected by co-immunoprecipitation in nuclear extracts of these cells, further suggesting that they may form a complex in myeloid cells to activate transcription. This complex formation is likely mediated by direct interactions between SETBP1/SETBP1(D/N) and MLL1 as both SETBP1 and SETBP1(D/N) are capable of interacting with multiple regions of MLL1 in binding assays using proteins synthesized by in vitro transcription and translation. To better understand the extent of SETBP1/SETBP1(D/N)-MLL1 interaction in regulating gene transcription, we carried out both ChIP-seq and RNA-seq analysis in mouse Lin -Sca-1 +c-Kit + (LSK) cells transduced by pMYs retrovirus expressing SETBP1 or SETBP1(D/N) or empty pMYs virus. These analyses revealed extensive overlap in genomic occupancy for MLL1 and SETBP1/SETBP1(D/N) and their cooperation in activating many oncogenic transcription factor genes in addition to Hoxa9/Hoxa10/Myb, including additional HoxA genes (Hoxa1, Hoxa3, Hoxa5, Hoxa6, and Hoxa7), Myc, Eya1, Mef2c, Meis1, Sox4, Mecom, and Lmo2. A large group of ribosomal protein genes were also found to be directly activated by MLL1 and SETBP1/SETBP1(D/N), identifying ribosomal biogenesis as another significant pathway induced by their cooperation. To further assess the requirement for MLL1 in SETBP1-induced transformation using a genetic approach, we also generated SETBP1/SETBP1(D/N)-induced immortalized myeloid progenitors and AMLs using LSK cells from Mll1 conditional knockout mice. Mll1 deletion in immortalized progenitors significantly decreased SETBP1/SETBP1(D/N)-induced transcriptional activation and their colony-forming potential. More importantly, Mll1 deletion significantly extended the survival of mice transplanted with SETBP1/SETBP1(D/N)-induced AMLs, indicating that Mll1 is essential for the maintenance of such leukemias in vivo. We further found that pharmacological inhibition of MLL1 complex using a WDR5 inhibitor OICR-9429 efficiently abrogated SETBP1/SETBP1(D/N)-induced transcriptional activation and transformation. Thus, MLL1 complex plays a critical role in Setbp1-induced transcriptional activation and transformation and represents a promising target for treating myeloid neoplasms with SETBP1 activation. Disclosures Maciejewski: Novartis: Consultancy; Regeneron: Consultancy; Alexion: Consultancy; Bristol Myers Squibb/Celgene: Consultancy.

ORPER: A Workflow for Constrained SSU rRNA Phylogenies

The continuous increase in sequenced genomes in public repositories makes the choice of interesting bacterial strains for future sequencing projects ever more complicated, as it is difficult to estimate the redundancy between these strains and the already available genomes. Therefore, we developed the Nextflow workflow “ORPER”, for “ORganism PlacER”, containerized in Singularity, which allows the determination the phylogenetic position of a collection of organisms in the genomic landscape. ORPER constrains the phylogenetic placement of SSU (16S) rRNA sequences in a multilocus reference tree based on ribosomal protein genes extracted from public genomes. We demonstrate the utility of ORPER on the Cyanobacteria phylum, by placing 152 strains of the BCCM/ULC collection.

ORPER: a workflow for constrained SSU rRNA phylogenies

The continuous increase of sequenced genomes in public repositories makes the choice of interesting bacterial strains for future sequencing projects evermore complicated, as it is difficult to estimate the redundancy between these strains and the already available genomes. Therefore, we developed the Nextflow workflow ORPER (containerized in Singularity), which allows determining the phylogenetic position of a collection of organisms in the genomic landscape. ORPER constrains the phylogenetic placement of SSU (16S) rRNA sequences in a multilocus reference tree based on ribosomal protein genes extracted from public genomes. We demonstrate the utility of ORPER on the Cyanobacteria phylum, by placing 152 strains of the BCCM/ULC collection.

Ribosome Proteins Represented by RPL27A Mark the Development and Metastasis of Triple-Negative Breast Cancer in Mouse and Human

Triple-negative breast cancer (TNBC) is known to have a poor prognosis and limited treatment options. The lack of targeted therapies and poor prognosis of patients with TNBC have made it urgent to discover novel critical diagnosis and therapeutic targets in the TNBC field. Here, in the current study, we integrated the single-cell RNA-sequencing (scRNA-seq) data from four normal mouse mammary tissues and four mouse breast tumors. Comparative analysis was conducted to identify the gene profiles of normal epithelial cells and cancer cells at different models. Surprisingly, two ribosomal protein genes, Rpl27a and Rpl15, were significantly upregulated in the cancer cells in all the TNBC models. Next, we accessed the scRNA-seq data from human primary and metastatic TNBC tissues, and comparative analysis revealed gene profiles of human primary and metastatic TNBC cancer cells. Ribosomal protein genes, represented by RPL27A and RPL15, showed significantly upregulated expression in metastatic TNBC cancer cells. Pathway analysis on the upregulated genes of the metastatic TNBC cancer cells identified the key regulators and signaling pathways that were driving the metastasis of the TNBC cancer cells. Specifically, EIF2 signaling was significantly activated, and major member genes of this signaling pathway were upregulated. In vitro study revealed that targeting RPL27A or EIF2 signaling in a TNBC cell line, MDA-MB-231, significantly reduced cell migration and invasion. Altogether, these data suggested that the RPL27A gene is conducting critical functions in TNBC cancer development and metastasis and is a potential therapeutic target for TNBC.

Gds1 interacts with NuA4 to promote H4 acetylation at ribosomal protein genes

In our previously published studies, RNA polymerase II transcription initiation complexes were assembled from yeast nuclear extracts onto immobilized transcription templates and analyzed by quantitative mass spectrometry. In addition to the expected basal factors and coactivators, we discovered that the uncharacterized protein Gds1 showed activator- stimulated association with promoter DNA. Gds1 co-precipitated with the histone H4 acetyltransferase NuA4, and its levels often tracked with NuA4 in immobilized template experiments. GDS1 deletion led to reduction in H4 acetylation in vivo and other phenotypes consistent with partial loss of NuA4 activity. Genome-wide chromatin immunoprecipitation revealed that the reduction in H4 acetylation was strongest at ribosomal protein gene promoters and other genes with high NuA4 occupancy. Therefore, while Gds1 is not a stoichiometric subunit of NuA4, we propose that it interacts with and modulates NuA4 in specific promoter contexts. Gds1 has no obvious metazoan homolog, but structural predictions suggest it may be distantly related to the DEK protein.

Insights into molecular structure, genome evolution and phylogenetic implication through mitochondrial genome sequence of Gleditsia sinensis

Gleditsia sinensis is an endemic species widely distributed in China with high economic and medicinal value. To explore the genomic evolution and phylogenetic relationships of G. sinensis, the complete mitochondrial (mt) genome of G. sinensis was sequenced and assembled, which was firstly reported in Gleditsia. The mt genome was circular and 594,121 bp in length, including 37 protein-coding genes (PCGs), 19 transfer RNA (tRNA) genes and 3 ribosomal RNA (rRNA) genes. The overall base composition of the G. sinensis mt genome was 27.4% for A, 27.4% for T, 22.6% for G, 22.7% for C. The comparative analysis of PCGs in Fabaceae species showed that most of the ribosomal protein genes and succinate dehydrogenase genes were lost. In addition, we found that the rps4 gene was only lost in G. sinensis, whereas it was retained in other Fabaceae species. The phylogenetic analysis based on shared PCGs of 24 species (22 Fabaceae and 2 Solanaceae) showed that G. sinensis is evolutionarily closer to Senna species. In general, this research will provide valuable information for the evolution of G. sinensis and provide insight into the phylogenetic relationships within the family Fabaceae.

Transcriptional control of ribosome biogenesis in yeast: links to growth and stress signals

Ribosome biogenesis requires prodigious transcriptional output in rapidly growing yeast cells and is highly regulated in response to both growth and stress signals. This minireview focuses on recent developments in our understanding of this regulatory process, with an emphasis on the 138 ribosomal protein genes (RPGs) themselves and a group of >200 ribosome biogenesis (RiBi) genes whose products contribute to assembly but are not part of the ribosome. Expression of most RPGs depends upon Rap1, a pioneer transcription factor (TF) required for the binding of a pair of RPG-specific TFs called Fhl1 and Ifh1. RPG expression is correlated with Ifh1 promoter binding, whereas Rap1 and Fhl1 remain promoter-associated upon stress-induced down regulation. A TF called Sfp1 has also been implicated in RPG regulation, though recent work reveals that its primary function is in activation of RiBi and other growth-related genes. Sfp1 plays an important regulatory role at a small number of RPGs where Rap1–Fhl1–Ifh1 action is subsidiary or non-existent. In addition, nearly half of all RPGs are bound by Hmo1, which either stabilizes or re-configures Fhl1–Ifh1 binding. Recent studies identified the proline rotamase Fpr1, known primarily for its role in rapamycin-mediated inhibition of the TORC1 kinase, as an additional TF at RPG promoters. Fpr1 also affects Fhl1–Ifh1 binding, either independently or in cooperation with Hmo1. Finally, a major recent development was the discovery of a protein homeostasis mechanism driven by unassembled ribosomal proteins, referred to as the Ribosome Assembly Stress Response (RASTR), that controls RPG transcription through the reversible condensation of Ifh1.

The roles of ribosomal proteins in nasopharyngeal cancer: culprits, sentinels or both

Ribosomal protein genes encode products that are essential for cellular protein biosynthesis and are major components of ribosomes. Canonically, they are involved in the complex system of ribosome biogenesis pivotal to the catalysis of protein translation. Amid this tightly organised process, some ribosomal proteins have unique spatial and temporal physiological activity giving rise to their extra-ribosomal functions. Many of these extra-ribosomal roles pertain to cellular growth and differentiation, thus implicating the involvement of some ribosomal proteins in organogenesis. Consequently, dysregulated functions of these ribosomal proteins could be linked to oncogenesis or neoplastic transformation of human cells. Their suspected roles in carcinogenesis have been reported but not specifically explained for malignancy of the nasopharynx. This is despite the fact that literature since one and half decade ago have documented the association of ribosomal proteins to nasopharyngeal cancer. In this review, we explain the association and contribution of dysregulated expression among a subset of ribosomal proteins to nasopharyngeal oncogenesis. The relationship of these ribosomal proteins with the cancer are explained. We provide information to indicate that the dysfunctional extra-ribosomal activities of specific ribosomal proteins are tightly involved with the molecular pathogenesis of nasopharyngeal cancer albeit mechanisms yet to be precisely defined. The complete knowledge of this will impact future applications in the effective management of nasopharyngeal cancer.