The RNA binding protein quaking regulates formation of circRNAs.

Circular RNAs (circRNAs), formed by non-sequential back-splicing of pre-mRNA transcripts, are a widespread form of non-coding RNA in animal cells. However, it is unclear whether the majority of circRNAs represent splicing by-products without function or are produced in a regulated manner to carry out specific cellular functions. We show that hundreds of circRNAs are regulated during human epithelial-mesenchymal transition (EMT) and find that the production of over one-third of abundant circRNAs is dynamically regulated by the alternative splicing factor, Quaking (QKI), which itself is regulated during EMT. Furthermore, by modulating QKI levels, we show the effect on circRNA abundance is dependent on intronic QKI binding motifs. Critically, the addition of QKI motifs is sufficient to induce de novo circRNA formation from transcripts that are normally linearly spliced. These findings demonstrate circRNAs are both purposefully synthesized and regulated by cell-type specific mechanisms, suggesting they play specific biological roles in EMT.

ESRP1 controls biogenesis and function of a large abundant multiexon circRNA.

While the majority of circRNAs are formed from infrequent back-splicing of exons from protein coding genes, some can be produced at quite high level and in a regulated manner. We describe the regulation, biogenesis and function of circDOCK1(2-27), a large, abundant circular RNA that is highly regulated during epithelial-mesenchymal transition (EMT) and whose formation depends on the epithelial splicing regulator ESRP1. CircDOCK1(2-27) synthesis in epithelial cells represses cell motility both by diverting transcripts from DOCK1 mRNA production to circRNA formation and by direct inhibition of migration by the circRNA. HITS-CLIP analysis and CRISPR-mediated deletions indicate ESRP1 controls circDOCK1(2-27) biosynthesis by binding a GGU-containing repeat region in intron 1 and detaining its splicing until Pol II completes its 157 kb journey to exon 27. Proximity-dependent biotinylation (BioID) assay suggests ESRP1 may modify the RNP landscape of intron 1 in a way that disfavours communication of exon 1 with exon 2, rather than physically bridging exon 2 to exon 27. The X-ray crystal structure of RNA-bound ESRP1 qRRM2 domain reveals it binds to GGU motifs, with the guanines embedded in clamp-like aromatic pockets in the protein.

SARSCoV-2 antibody prevalence and titers in persons living with HIV cared for at a large tertiary reference center in Mexico City.

OBJECTIVE: To assess SARS-CoV-2 antibody prevalence and titers in people living with HIV (PLWHIV) on antiretroviral treatment (ART) enrolled at a tertiary reference hospital in Mexico. METHODS: Two plasma aliquots per person, used for HIV viral load follow-up between 01/2020 and 09/2021, were used to assess total anti-N and neutralizing SARS-CoV-2 antibodies. Sociodemographic, clinical, and SARS-CoV-2 exposure risk information were collected. The risk associated with SARS-CoV-2 exposure and associations with antibody titers were analyzed with logistic, Cox, and linear multivariable models. RESULTS: 803 PLWHIV participated; 233 had detectable SARS-CoV-2 antibodies (prevalent cases), and 132 seroconverted (incident cases). Overall, the adjusted prevalence was 46.45%, with an incidence rate of 3.78 cases/100 person-months. Factors associated with prevalent cases included lower age, location (western zone of Mexico City and the neighboring Mexico State), use of public transport, attendance at meetings without social distancing, and higher CD4 + T cell counts (p < 0.05; multivariable logistic model). BNT162b2 vaccination reduced incident cases (Cox adjusted HR = 0.4; p = 0.013). Notably, previously infected and vaccinated individuals showed maximization of neutralizing activity (p < 0.001). No associations between SARS-CoV-2 neutralization and HIV-related variables (CD4 + T cell counts, viral load, number of years in viral suppression, ART regimen) were found in multivariable analysis. CONCLUSIONS: SARS-CoV-2 infection was associated with community risk rather than HIV-associated variables in PLWH on ART and clinical follow-up. Antibody neutralization activity in vaccinated participants was maximized with previous SARS-CoV-2 infection.

SplintQuant: a method for accurately quantifying circular RNA transcript abundance without reverse transcription bias.

Reverse transcription of RNA is fallible, introducing biases and confounding the quantification of transcript abundance. We demonstrate that circular RNAs (circRNAs) are more subjective to overestimation of transcript abundance than cognate linear RNAs due to their covalently closed, circular form, producing multiple concatameric products from a single priming of reverse transcriptase. We developed SplintQuant, where custom DNA oligonucleotides are ligated by PBCV-1 DNA ligase only when bound to their target RNA. These circRNA-specific DNA oligonucleotides are terminally tagged with universal primers, allowing SplintQuant to accurately quantify even lowly abundant circRNAs through highly specific quantitative PCR (qPCR) in the absence of reverse transcription. SplintQuant is sensitive, specific, highly reproducible, and applicable to the quantification of canonical and noncanonical RNA transcripts including alternative splice variants, gene fusions, and offers a gold-standard approach for accurately quantifying circRNAs.

Tetramerization of MADS family transcription factors SEPALLATA3 and AGAMOUS is required for floral meristem determinacy in Arabidopsis.

The MADS transcription factors (TF) constitute an ancient family of TF found in all eukaryotes that bind DNA as obligate dimers. Plants have dramatically expanded the functional diversity of the MADS family during evolution by adding protein-protein interaction domains to the core DNA-binding domain, allowing the formation of heterotetrameric complexes. Tetramerization of plant MADS TFs is believed to play a central role in the evolution of higher plants by acting as one of the main determinants of flower formation and floral organ specification. The MADS TF, SEPALLATA3 (SEP3), functions as a central protein-protein interaction hub, driving tetramerization with other MADS TFs. Here, we use a SEP3 splice variant, SEP3Δtet, which has dramatically abrogated tetramerization capacity to decouple SEP3 tetramerization and DNA-binding activities. We unexpectedly demonstrate that SEP3 heterotetramer formation is required for correct termination of the floral meristem, but plays a lesser role in floral organogenesis. The heterotetramer formed by SEP3 and the MADS protein, AGAMOUS, is necessary to activate two target genes, KNUCKLES and CRABSCLAW, which are required for meristem determinacy. These studies reveal unique and highly specific roles of tetramerization in flower development and suggest tetramerization may be required to activate only a subset of target genes in closed chromatin regions.

Structural basis for the oligomerization of the MADS domain transcription factor SEPALLATA3 in Arabidopsis.

In plants, MADS domain transcription factors act as central regulators of diverse developmental pathways. In Arabidopsis thaliana, one of the most central members of this family is SEPALLATA3 (SEP3), which is involved in many aspects of plant reproduction, including floral meristem and floral organ development. SEP3 has been shown to form homo and heterooligomeric complexes with other MADS domain transcription factors through its intervening (I) and keratin-like (K) domains. SEP3 function depends on its ability to form specific protein-protein complexes; however, the atomic level determinants of oligomerization are poorly understood. Here, we report the 2.5-Å crystal structure of a small portion of the intervening and the complete keratin-like domain of SEP3. The domains form two amphipathic alpha helices separated by a rigid kink, which prevents intramolecular association and presents separate dimerization and tetramerization interfaces comprising predominantly hydrophobic patches. Mutations to the tetramerization interface demonstrate the importance of highly conserved hydrophobic residues for tetramer stability. Atomic force microscopy was used to show SEP3-DNA interactions and the role of oligomerization in DNA binding and conformation. Based on these data, the oligomerization patterns of the larger family of MADS domain transcription factors can be predicted and manipulated based on the primary sequence.

Circular RNA in cancer.

Over the past decade, circular RNA (circRNA) research has evolved into a bona fide research field shedding light on the functional consequence of this unique family of RNA molecules in cancer. Although the method of formation and the abundance of circRNAs can differ from their cognate linear mRNA, the spectrum of interacting partners and their resultant cellular functions in oncogenesis are analogous. However, with 10 times more diversity in circRNA variants compared with linear RNA variants, combined with their hyperstability in the cell, circRNAs are equipped to influence every stage of oncogenesis. This is an opportune time to address the breadth of circRNA in cancer focused on their spatiotemporal expression, mutations in biogenesis factors and contemporary functions through each stage of cancer. In this Review, we highlight examples of functional circRNAs in specific cancers, which satisfy critical criteria, including their physical co-association with the target and circRNA abundance at stoichiometrically valid quantities. These considerations are essential to develop strategies for the therapeutic exploitation of circRNAs as biomarkers and targeted anticancer agents.

Native Circular RNA Pulldown Method to Simultaneously Profile RNA and Protein Interactions.

Circular RNAs (circRNAs) are a widespread, cell-, tissue-, and disease-specific class of largely non-coding RNA transcripts. These single-stranded, covalently-closed transcripts arise through non-canonical splicing of pre-mRNA, a process called back-splicing. Back-splicing results in circRNAs which are distinguishable from their cognate mRNA as they possess a unique sequence of nucleic acids called the backsplice junction (BSJ). CircRNAs have been shown to play key functional roles in various cellular contexts and achieve this through their interaction with other macromolecules, particularly other RNA molecules and proteins. To elucidate the molecular mechanisms underlying circRNA function, it is necessary to identify these interacting partners. Herein, we present an optimized strategy for the simultaneous purification of the circRNA interactome within eukaryotic cells, allowing the identification of both circRNA-RNA and circRNA-protein interactions.

AUTO-TUNE: selecting the distance threshold for inferring HIV transmission clusters.

Molecular surveillance of viral pathogens and inference of transmission networks from genomic data play an increasingly important role in public health efforts, especially for HIV-1. For many methods, the genetic distance threshold used to connect sequences in the transmission network is a key parameter informing the properties of inferred networks. Using a distance threshold that is too high can result in a network with many spurious links, making it difficult to interpret. Conversely, a distance threshold that is too low can result in a network with too few links, which may not capture key insights into clusters of public health concern. Published research using the HIV-TRACE software package frequently uses the default threshold of 0.015 substitutions/site for HIV pol gene sequences, but in many cases, investigators heuristically select other threshold parameters to better capture the underlying dynamics of the epidemic they are studying. Here, we present a general heuristic scoring approach for tuning a distance threshold adaptively, which seeks to prevent the formation of giant clusters. We prioritize the ratio of the sizes of the largest and the second largest cluster, maximizing the number of clusters present in the network. We apply our scoring heuristic to outbreaks with different characteristics, such as regional or temporal variability, and demonstrate the utility of using the scoring mechanism's suggested distance threshold to identify clusters exhibiting risk factors that would have otherwise been more difficult to identify. For example, while we found that a 0.015 substitutions/site distance threshold is typical for US-like epidemics, recent outbreaks like the CRF07_BC subtype among men who have sex with men (MSM) in China have been found to have a lower optimal threshold of 0.005 to better capture the transition from injected drug use (IDU) to MSM as the primary risk factor. Alternatively, in communities surrounding Lake Victoria in Uganda, where there has been sustained heterosexual transmission for many years, we found that a larger distance threshold is necessary to capture a more risk factor-diverse population with sparse sampling over a longer period of time. Such identification may allow for more informed intervention action by respective public health officials.

AUTO-TUNE: SELECTING THE DISTANCE THRESHOLD FOR INFERRING HIV TRANSMISSION CLUSTERS.

Molecular surveillance of viral pathogens and inference of transmission networks from genomic data play an increasingly important role in public health efforts, especially for HIV-1. For many methods, the genetic distance threshold used to connect sequences in the transmission network is a key parameter informing the properties of inferred networks. Using a distance threshold that is too high can result in a network with many spurious links, making it difficult to interpret. Conversely, a distance threshold that is too low can result in a network with too few links, which may not capture key insights into clusters of public health concern. Published research using the HIV-TRACE software package frequently uses the default threshold of 0.015 substitutions/site for HIV pol gene sequences, but in many cases, investigators heuristically select other threshold parameters to better capture the underlying dynamics of the epidemic they are studying. Here, we present a general heuristic scoring approach for tuning a distance threshold adaptively, which seeks to prevent the formation of giant clusters. We prioritize the ratio of the sizes of the largest and the second largest cluster, maximizing the number of clusters present in the network. We apply our scoring heuristic to outbreaks with different characteristics, such as regional or temporal variability, and demonstrate the utility of using the scoring mechanism's suggested distance threshold to identify clusters exhibiting risk factors that would have otherwise been more difficult to identify. For example, while we found that a 0.015 substitutions/site distance threshold is typical for US-like epidemics, recent outbreaks like the CRF07_BC subtype among men who have sex with men (MSM) in China have been found to have a lower optimal threshold of 0.005 to better capture the transition from injected drug use (IDU) to MSM as the primary risk factor. Alternatively, in communities surrounding Lake Victoria in Uganda, where there has been sustained hetero-sexual transmission for many years, we found that a larger distance threshold is necessary to capture a more risk factor-diverse population with sparse sampling over a longer period of time. Such identification may allow for more informed intervention action by respective public health officials.

The response of the maize nitrate transport system to nitrogen demand and supply across the lifecycle.

An understanding of nitrate (NO3-) uptake throughout the lifecycle of plants, and how this process responds to nitrogen (N) availability, is an important step towards the development of plants with improved nitrogen use efficiency (NUE). NO3- uptake capacity and transcript levels of putative high- and low-affinity NO3- transporters (NRTs) were profiled across the lifecycle of dwarf maize (Zea mays) plants grown at reduced and adequate NO3-. Plants showed major changes in high-affinity NO3- uptake capacity across the lifecycle, which varied with changing relative growth rates of roots and shoots. Transcript abundances of putative high-affinity NRTs (predominantly ZmNRT2.1 and ZmNRT2.2) were correlated with two distinct peaks in high-affinity root NO3- uptake capacity and also N availability. The reduction in NO3- supply during the lifecycle led to a dramatic increase in NO3- uptake capacity, which preceded changes in transcript levels of NRTs, suggesting a model with short-term post-translational regulation and longer term transcriptional regulation of NO3- uptake capacity. These observations offer new insight into the control of NO3- uptake by both plant developmental processes and N availability, and identify key control points that may be targeted by future plant improvement programmes to enhance N uptake relative to availability and/or demand.

Versatile toolkit for highly-efficient and scarless overexpression of circular RNAs.

Circular RNAs (circRNAs) are a class of single-stranded, covalently closed RNA that contain a unique back-splice junction (bsj) sequence created by the ligation of their 5' and 3' ends via spliceosome-catalyzed back-splicing. A key step in illuminating the cellular roles of specific circRNAs is via increasing their expression. This is frequently done by transfecting cells with plasmid DNA containing cloned exons from which the circRNA is transcribed, flanked by sequences that promote back-splicing. We observed that commonly used plasmids lead to the production of circRNAs with molecular scars at the circRNA bsj. Stepwise redesign of the cloning vector corrected this problem, ensuring bona fide circRNAs are produced with their natural bsj at high efficiency. The fidelity of circRNAs produced from this new construct was validated by RNA sequencing and also functionally validated. To increase the utility of this modified resource for expressing circRNA, we developed an expanded set of vectors incorporating this design that (i) enables selection with a variety of antibiotics and fluorescent proteins, (ii) employs a range of promoters varying in promoter strength and (iii) generated a complementary set of lentiviral plasmids for difficult-to-transfect cells. These resources provide a novel and versatile toolkit for high-efficiency and scarless overexpression of circular RNAs that fulfill a critical need for the investigation of circRNA function.

Functional Characterisation of the Circular RNA, circHTT(2-6), in Huntington's Disease.

Trinucleotide repeat disorders comprise ~20 severe, inherited, human neuromuscular and neurodegenerative disorders, which result from an abnormal expansion of repetitive sequences in the DNA. The most common of these, Huntington's disease (HD), results from expansion of the CAG repeat region in exon 1 of the HTT gene via an unknown mechanism. Since non-coding RNAs have been implicated in the initiation and progression of many diseases, herein we focused on a circular RNA (circRNA) molecule arising from non-canonical splicing (backsplicing) of HTT pre-mRNA. The most abundant circRNA from HTT, circHTT(2-6), was found to be more highly expressed in the frontal cortex of HD patients, compared with healthy controls, and positively correlated with CAG repeat tract length. Furthermore, the mouse orthologue (mmu_circHTT(2-6)) was found to be enriched within the brain and specifically the striatum, a region enriched for medium spiny neurons that are preferentially lost in HD. Transgenic overexpression of circHTT(2-6) in two human cell lines-SH-SY5Y and HEK293-reduced cell proliferation and nuclear size without affecting cell cycle progression or cellular size, or altering the CAG repeat region length within HTT. CircHTT(2-6) overexpression did not alter total HTT protein levels, but reduced its nuclear localisation. As these phenotypic and genotypic changes resemble those observed in HD patients, our results suggest that circHTT(2-6) may play a functional role in the pathophysiology of this disease.

Circular RNAs drive oncogenic chromosomal translocations within the MLL recombinome in leukemia.

The first step of oncogenesis is the acquisition of a repertoire of genetic mutations to initiate and sustain the malignancy. An important example of this initiation phase in acute leukemias is the formation of a potent oncogene by chromosomal translocations between the mixed lineage leukemia (MLL) gene and one of 100 translocation partners, known as the MLL recombinome. Here, we show that circular RNAs (circRNAs)-a family of covalently closed, alternatively spliced RNA molecules-are enriched within the MLL recombinome and can bind DNA, forming circRNA:DNA hybrids (circR loops) at their cognate loci. These circR loops promote transcriptional pausing, proteasome inhibition, chromatin re-organization, and DNA breakage. Importantly, overexpressing circRNAs in mouse leukemia xenograft models results in co-localization of genomic loci, de novo generation of clinically relevant chromosomal translocations mimicking the MLL recombinome, and hastening of disease onset. Our findings provide fundamental insight into the acquisition of chromosomal translocations by endogenous RNA carcinogens in leukemia.

Seroepidemiology of SARS-CoV-2 in healthcare personnel working at the largest tertiary COVID-19 referral hospitals in Mexico City.

INTRODUCTION: We performed a longitudinal SARS-CoV-2 seroepidemiological study in healthcare personnel of the two largest tertiary COVID-19 referral hospitals in Mexico City. METHODS: All healthcare personnel, including staff physicians, physicians in training, nurses, laboratory technicians, researchers, students, housekeeping, maintenance, security, and administrative staff were invited to voluntarily participate, after written informed consent. Participants answered a computer-assisted self-administered interview and donated blood samples for antibody testing every three weeks from October 2020 to June 2021. RESULTS: A total of 883 participants (out of 3639 registered employees) contributed with at least one blood sample. The median age was 36 years (interquartile range: 28-46) and 70% were women. The most common occupations were nurse (28%), physician (24%), and administrative staff (22%). Two hundred and ninety participants (32.8%) had a positive-test result in any of the visits, yielding an overall adjusted prevalence of 33.5% for the whole study-period. Two hundred and thirty-five positive tests were identified at the baseline visit (prevalent cases), the remaining 55 positive tests were incident cases. Prevalent cases showed associations with both occupational (institution 2 vs. 1: adjusted odds ratio [aOR] = 2.24, 95% confidence interval [CI]: 1.54-3.25; laboratory technician vs. physician: aOR = 4.38, 95% CI: 1.75-10.93) and community (municipality of residence Xochimilco vs. Tlalpan: aOR = 2.03, 95% CI: 1.09-3.79) risk-factors. The incidence rate was 3.0 cases per 100 person-months. Incident cases were associated with community-acquired risk, due to contact with suspect/confirmed COVID-19 cases (HR = 2.45, 95% CI: 1.21-5.00). CONCLUSIONS: We observed that between October 2020 and June 2021, healthcare workers of the two largest tertiary COVID-19 referral centers in Mexico City had similar level of exposure to SARS-CoV-2 than the general population. Most variables associated with exposure in this setting pointed toward community rather than occupational risk. Our observations are consistent with successful occupational medicine programs for SARS-CoV-2 infection control in the participating institutions but suggest the need to strengthen mitigation strategies in the community.

Magnesium transporters, MGT2/MRS2-1 and MGT3/MRS2-5, are important for magnesium partitioning within Arabidopsis thaliana mesophyll vacuoles.

• Magnesium accumulates at high concentrations in dicotyledonous leaves but it is not known in which leaf cell types it accumulates, by what mechanism this occurs and the role it plays when stored in the vacuoles of these cell types. • Cell-specific vacuolar elemental profiles from Arabidopsis thaliana (Arabidopsis) leaves were analysed by X-ray microanalysis under standard and serpentine hydroponic growth conditions and correlated with the cell-specific complement of magnesium transporters identified through microarray analysis and quantitative polymerase chain reaction (qPCR). • Mesophyll cells accumulate the highest vacuolar concentration of magnesium in Arabidopsis leaves and are enriched for members of the MGT/MRS2 family of magnesium transporters. Specifically, AtMGT2/AtMRS2-1 and AtMGT3/AtMRS2-5 were shown to be targeted to the tonoplast and corresponding T-DNA insertion lines had perturbed mesophyll-specific vacuolar magnesium accumulation under serpentine conditions. Furthermore, transcript abundance of these genes was correlated with the accumulation of magnesium under serpentine conditions, in a low calcium-accumulating mutant and across 23 Arabidopsis ecotypes varying in their leaf magnesium concentrations. • We implicate magnesium as a key osmoticum required to maintain growth in low calcium concentrations in Arabidopsis. Furthermore, two tonoplast-targeted members of the MGT/MRS2 family are shown to contribute to this mechanism under serpentine conditions.

Characteristics and growth of the genetic HIV transmission network of Mexico City during 2020.

INTRODUCTION: Molecular surveillance systems could provide public health benefits to focus strategies to improve the HIV care continuum. Here, we infer the HIV genetic network of Mexico City in 2020, and identify actively growing clusters that could represent relevant targets for intervention. METHODS: All new diagnoses, referrals from other institutions, as well as persons returning to care, enrolling at the largest HIV clinic in Mexico City were invited to participate in the study. The network was inferred from HIV pol sequences, using pairwise genetic distance methods, with a locally hosted, secure version of the HIV-TRACE tool: Seguro HIV-TRACE. Socio-demographic, clinical and behavioural metadata were overlaid across the network to design focused prevention interventions. RESULTS: A total of 3168 HIV sequences from unique individuals were included. One thousand and one-hundred and fifty (36%) sequences formed 1361 links within 386 transmission clusters in the network. Cluster size varied from 2 to 14 (63% were dyads). After adjustment for covariates, lower age (adjusted odds ratio [aOR]: 0.37, p<0.001; >34 vs. <24 years), being a man who has sex with men (MSM) (aOR: 2.47, p = 0.004; MSM vs. cisgender women), having higher viral load (aOR: 1.28, p<0.001) and higher CD4+ T cell count (aOR: 1.80, p<0.001; ≥500 vs. <200 cells/mm3 ) remained associated with higher odds of clustering. Compared to MSM, cisgender women and heterosexual men had significantly lower education (none or any elementary: 59.1% and 54.2% vs. 16.6%, p<0.001) and socio-economic status (low income: 36.4% and 29.0% vs. 18.6%, p = 0.03) than MSM. We identified 10 (2.6%) clusters with constant growth, for prioritized intervention, that included intersecting sexual risk groups, highly connected nodes and bridge nodes between possible sub-clusters with high growth potential. CONCLUSIONS: HIV transmission in Mexico City is strongly driven by young MSM with higher education level and recent infection. Nevertheless, leveraging network inference, we identified actively growing clusters that could be prioritized for focused intervention with demographic and risk characteristics that do not necessarily reflect the ones observed in the overall clustering population. Further studies evaluating different models to predict growing clusters are warranted. Focused interventions will have to consider structural and risk disparities between the MSM and the heterosexual populations.

HIV Pretreatment Drug Resistance Trends in Mexico City, 2017-2020.

In response to increasing pretreatment drug resistance (PDR), Mexico changed its national antiretroviral treatment (ART) policy, recommending and procuring second-generation integrase strand-transfer inhibitor (INSTI)-based regimens as preferred first-line options since 2019. We present a four-year observational study describing PDR trends across 2017-2020 at the largest HIV diagnosis and primary care center in Mexico City. A total of 6688 baseline protease-reverse transcriptase and 6709 integrase sequences were included. PDR to any drug class was 14.4% (95% CI, 13.6-15.3%). A significant increasing trend for efavirenz/nevirapine PDR was observed (10.3 to 13.6%, p = 0.02). No increase in PDR to second-generation INSTI was observed, remaining under 0.3% across the study period. PDR was strongly associated with prior exposure to ART (aOR: 2.9, 95% CI: 1.9-4.6, p < 0.0001). MSM had higher odds of PDR to efavirenz/nevirapine (aOR: 2.0, 95% CI: 1.0-3.7, p = 0.04), reflecting ongoing transmission of mutations such as K103NS and E138A. ART restarters showed higher representation of cisgender women and injectable drug users, higher age, and lower education level. PDR to dolutegravir/bictegravir remained low in Mexico City, although further surveillance is warranted given the short time of ART optimization. Our study identifies demographic characteristics of groups with higher risk of PDR and lost to follow-up, which may be useful to design differentiated interventions locally.

SRRM4 Expands the Repertoire of Circular RNAs by Regulating Microexon Inclusion.

High-throughput RNA sequencing (RNA-seq) and dedicated bioinformatics pipelines have synergized to identify an expansive repertoire of unique circular RNAs (circRNAs), exceeding 100,000 variants. While the vast majority of these circRNAs comprise canonical exonic and intronic sequences, microexons (MEs)-which occur in 30% of functional mRNA transcripts-have been entirely overlooked. CircRNAs which contain these known MEs (ME-circRNAs) could be identified with commonly utilized circRNA prediction pipelines, CIRCexplorer2 and CIRI2, but were not previously recognized as ME-circRNAs. In addition, when employing a bespoke bioinformatics pipeline for identifying RNA chimeras, called Hyb, we could also identify over 2000 ME-circRNAs which contain novel MEs at their backsplice junctions, that are uncalled by either CIRCexplorer2 or CIRI2. Analysis of circRNA-seq datasets from gliomas of varying clinical grades compared with matched control tissue has shown circRNAs have potential as prognostic markers for stratifying tumor from healthy tissue. Furthermore, the abundance of microexon-containing circRNAs (ME-circRNAs) between tumor and normal tissues is correlated with the expression of a splicing associated factor, Serine/arginine repetitive matrix 4 (SRRM4). Overexpressing SRRM4, known for regulating ME inclusion in mRNAs critical for neural differentiation, in human HEK293 cells resulted in the biogenesis of over 2000 novel ME-circRNAs, including ME-circEIF4G3, and changes in the abundance of many canonical circRNAs, including circSETDB2 and circLBRA. This shows SRRM4, in which its expression is correlated with poor prognosis in gliomas, acts as a bona fide circRNA biogenesis factor. Given the known roles of MEs and circRNAs in oncogenesis, the identification of these previously unrecognized ME-circRNAs further increases the complexity and functional purview of this non-coding RNA family.