Identification of human genetic variants controlling circular RNA expression

An integrative pipeline for circular RNA quantitative trait locus discovery with application in human T cells

MOTIVATION

Molecular quantitative trait locus (QTL) mapping has proven to be a powerful approach for prioritizing genetic regulatory variants and causal genes identified by genome-wide association studies. Recently, this success has been extended to circular RNA (circRNA), a potential group of RNAs that can serve as markers for the diagnosis, prognosis, or therapeutic targets of various human diseases. However, a well-developed computational pipeline for circRNA QTL (circQTL) discovery is still lacking.

RESULTS

We introduce an integrative method for circQTL mapping and implement it as an automated pipeline based on Nextflow, named cscQTL. The proposed method has two main advantages. Firstly, cscQTL improves the specificity by systematically combining outputs of multiple circRNA calling algorithms to obtain highly confident circRNA annotations. Secondly, cscQTL improves the sensitivity by accurately quantifying circRNA expression with the help of pseudo references. Compared to the single method approach, cscQTL effectively identifies circQTLs with an increase of 20%-100% circQTLs detected and recovered all circQTLs that are highly supported by the single method approach. We apply cscQTL to a dataset of human T cells and discover genetic variants that control the expression of 55 circRNAs. By colocalization tests, we further identify circBACH2 and circYY1AP1 as potential candidates for immune disease regulation.

AVAILABILITY AND IMPLEMENTATION

cscQTL is freely available at: https://github.com/datngu/cscQTL and https://doi.org/10.5281/zenodo.7851982.

Circular RNAs arising from synaptic host genes during human neuronal differentiation are modulated by SFPQ RNA-binding protein

BACKGROUND

Circular RNA (circRNA) molecules, generated through non-canonical back-splicing of exon-exon junctions, have recently been implicated in diverse biological functions including transcriptional regulation and modulation of protein interactions. CircRNAs are emerging as a key component of the complex neural transcriptome implicated in brain development. However, the specific expression patterns and functions of circRNAs in human neuronal differentiation have not been explored.

RESULTS

Using total RNA sequencing analysis, we identified expressed circRNAs during the differentiation of human neuroepithelial stem (NES) cells into developing neurons and discovered that many circRNAs originated from host genes associated with synaptic function. Interestingly, when assessing population data, exons giving rise to circRNAs in our dataset had a higher frequency of genetic variants. Additionally, screening for RNA-binding protein sites identified enrichment of Splicing Factor Proline and Glutamine Rich (SFPQ) motifs in increased circRNAs, several of which were reduced by SFPQ knockdown and enriched in SFPQ ribonucleoprotein complexes.

CONCLUSIONS

Our study provides an in-depth characterisation of circRNAs in a human neuronal differentiation model and highlights SFPQ as both a regulator and binding partner of circRNAs elevated during neuronal maturation.

Insight on Non-Coding RNAs from Biofluids in Ovarian Tumors

Ovarian tumors are the most frequent adnexal mass, raising diagnostic and therapeutic issues linked to a large spectrum of tumors, with a continuum from benign to malignant. Thus far, none of the available diagnostic tools have proven efficient in deciding strategy, and no consensus exists on the best strategy between "single test", "dual testing", "sequential testing", "multiple testing options" and "no testing". In addition, there is a need for prognostic tools such as biological markers of recurrence and theragnostic tools to detect women not responding to chemotherapy in order to adapt therapies. Non-coding RNAs are classified as small or long based on their nucleotide count. Non-coding RNAs have multiple biological functions such as a role in tumorigenesis, gene regulation and genome protection. These ncRNAs emerge as new potential tools to differentiate benign from malignant tumors and to evaluate prognostic and theragnostic factors. In the specific setting of ovarian tumors, the goal of the present work is to offer an insight into the contribution of biofluid non-coding RNAs (ncRNA) expression.

Genetic regulation of circular RNA expression in human aortic smooth muscle cells and vascular traits

Circular RNAs (circRNAs) are a class of non-coding RNAs that have cell-type-specific expression and are relevant in cardiovascular disease. Aortic smooth muscle cells (SMCs) play a crucial role in cardiovascular disease. In this study, we employed a systems genetics approach to identify SMC circRNA transcripts and their relevance in cardiovascular traits across the genome. We quantified circRNA expression across 151 quiescent and proliferative human aortic SMCs from donors of various genetic ancestries. We identified 1,589 expressed circRNAs. Between quiescent and proliferative SMCs, we identified 173 differentially expressed circRNAs. To characterize the genetic regulation of circRNA expression, we associated the genotypes of 6.3 million single nucleotide polymorphisms (SNPs) with circRNA abundance and found 96 circRNAs that were associated with genetic loci. Three SNPs were associated with circRNA expression in proliferative SMCs but not quiescent SMCs. We identified six SNPs that had distinct association directions with circRNA isoforms from the same gene. Lastly, to identify the relevance of circRNAs in cardiovascular disease, we overlapped genetic loci associated with circRNA expression with vascular disease-related genome-wide association studies loci. We identified 14 blood pressure, one myocardial infarction, and three coronary artery disease loci, which were associated with a circRNA transcript but not an mRNA transcript. Overall, our results provide insight into the genetic basis of vascular disease traits mediated by circRNA expression.

The circular RNA landscape in multiple sclerosis: Disease-specific associated variants and exon methylation shape circular RNA expression profile

BACKGROUND

Circular RNAs (circRNAs) are a class of non-coding RNAs increasingly emerging as crucial actors in the pathogenesis of human diseases, including autoimmune and neurological disorders as multiple sclerosis (MS). Despite several efforts, the mechanisms regulating circRNAs expression are still largely unknown and the circRNA profile and regulation in MS-relevant cell models has not been completely investigated. In this work, we aimed at exploring the global landscape of circRNA expression in MS patients, also evaluating a possible correlation with their genetic and epigenetic background.

METHODS

We performed RNA-seq experiments on circRNA-enriched samples, derived from peripheral blood mononuclear cells (PBMCs) of 10 MS patients and 10 matched controls and performed differential circRNA expression. The genetic background was evaluated using array genotyping, and an expression quantitative trait loci (eQTL) analysis was carried out.

RESULTS

Expression analysis revealed 166 differentially expressed circRNAs in MS patients, 125 of which are downregulated. One of the top dysregulated circRNAs, hsa_circ_0007990, derives from the PGAP3 gene, encoding a protein relevant for the control of autoimmune responses. The downregulation of this circRNA was confirmed in two independent replication cohorts, suggesting its implementation as a possible RNA-based biomarker. The eQTL analysis evidenced a significant association between 89 MS-associated loci and the expression of at least one circRNA, suggesting that MS-associated variants could impact on disease pathogenesis by altering circRNA profiles. Finally, we found a significant correlation between exon methylation and circRNA expression levels, supporting the hypothesis that epigenetic features may play an important role in the definition of the cell circRNA pool.

CONCLUSION

We described the circRNA expression profile of PBMCs in MS patients, suggesting that MS-associated variants may tune the expression levels of circRNAs acting as "circ-QTLs", and proposing a role for exon-based DNA methylation in regulating circRNA expression.

Evidence for correlations between BMI-associated SNPs and circRNAs

Circular RNAs (circRNAs) are regulators of processes like adipogenesis. Their expression can be modulated by SNPs. We analysed links between BMI-associated SNPs and circRNAs. First, we detected an enrichment of BMI-associated SNPs on circRNA genomic loci in comparison to non-significant variants. Analysis of sex-stratified GWAS data revealed that circRNA genomic loci encompassed more genome-wide significant BMI-SNPs in females than in males. To explore whether the enrichment is restricted to BMI, we investigated nine additional GWAS studies. We showed an enrichment of trait-associated SNPs in circRNAs for four analysed phenotypes (body height, chronic kidney disease, anorexia nervosa and autism spectrum disorder). To analyse the influence of BMI-affecting SNPs on circRNA levels in vitro, we examined rs4752856 located on hsa_circ_0022025. The analysis of heterozygous individuals revealed an increased level of circRNA derived from the BMI-increasing SNP allele. We conclude that genetic variation may affect the BMI partly through circRNAs.

Circular RNA circRNF169 functions as a miR-30c-5p sponge to promote cellular senescence

Circular RNAs (circRNAs), characterized as single-stranded closed circular RNA molecules, have been established to exert pivotal functions in various biological or pathological processes. Nonetheless, the effects and underlying mechanisms concerning circRNAs on the aging and aging-related diseases remain elusive. We herein compared the expression patterns of circRNAs in young and senescent mouse embryonic fibroblasts (MEFs), and uncovered that circRNF169 was dramatically up-regulated in senescent MEFs compared with that in young MEFs. Therefore, we further digged into the role and potential mechanisms of circRNF169 in the senescence of MEFs. The results of senescence-associate-β-galactosidase staining and BrdU incorporation assay showed that silencing of circRNF169 significantly delayed MEFs senescence and promoted cell proliferation, while ectopic expression of circRNF169 exhibited the opposite effects. Moreover, the dual-luciferase reporter assay confirmed that circRNF169 acted as an endogenous miR-30c-5p sponge, which accelerated cellular senescence by sequestering and inhibiting miR-30c-5p activity. Taken together, our results suggested that circRNF169 exerted a crucial role in cellular senescence through sponging miR-30c-5p and represented a promising target for aging intervention.

Towards the Genetic Architecture of Complex Gene Expression Traits: Challenges and Prospects for eQTL Mapping in Humans

The discovery of expression quantitative trait loci (eQTLs) and their target genes (eGenes) has not only compensated for the limitations of genome-wide association studies for complex phenotypes but has also provided a basis for predicting gene expression. Efforts have been made to develop analytical methods in statistical genetics, a key discipline in eQTL analysis. In particular, mixed model– and deep learning–based analytical methods have been extremely beneficial in mapping eQTLs and predicting gene expression. Nevertheless, we still face many challenges associated with eQTL discovery. Here, we discuss two key aspects of these challenges: 1, the complexity of eTraits with various factors such as polygenicity and epistasis and 2, the voluminous work required for various types of eQTL profiles. The properties and prospects of statistical methods, including the mixed model method, Bayesian inference, the deep learning method, and the integration method, are presented as future directions for eQTL discovery. This review will help expedite the design and use of efficient methods for eQTL discovery and eTrait prediction.

circTAF8 Regulates Myoblast Development and Associated Carcass Traits in Chicken

Recent studies have shown that circular RNAs (circRNAs) play important roles in skeletal muscle development. CircRNA biogenesis is dependent on the genetic context. Single-nucleotide polymorphisms in the introns flanking circRNAs may be intermediate-inducible factors between circRNA expression and phenotypic traits. Our previous study showed that circTAF8 is an abundantly and differentially expressed circRNA in leg muscle during chicken embryonic development. Here, we aimed to investigate circTAF8 function in muscle development and the association of the SNPs in the circTAF8 flanking introns with carcass traits. In this study, we observed that overexpression of circTAF8 could promote the proliferation of chicken primary myoblasts and inhibit their differentiation. In addition, the SNPs in the introns flanking the circTAF8 locus and those associated with chicken carcass traits were analyzed in 335 partridge chickens. A total of eight SNPs were found associated with carcass traits such as leg muscle weight, live weight, and half and full-bore weight. The association analysis results of haplotype combinations were consistent with the association analysis of a single SNP. These results suggest that circTAF8 plays a regulatory role in muscle development. These identified SNPs were found correlated with traits to muscle development and carcass muscle weight in chickens.

CircRNAs as potent biomarkers in ovarian cancer: a systematic scoping review

More powerful prognostic and diagnostic tools are urgently needed for identifying and treating ovarian cancer (OC), which is the most fatal malignancy in women in developed countries. Circular RNAs (circRNAs) are conservative and stable looped molecules that can regulate gene expression by competing with other endogenous microRNA sponges. This discovery provided new insight into novel methods for regulating genes that are involved in many disorders and cancers. This review focuses on the dysregulated expression of circRNAs as well as their diagnostic and prognostic values in OC. We found that studies have identified twenty-one downregulated circRNAs and fifty-seven upregulated ones. The results of these studies confirm that circRNAs might be potent biomarkers with diagnostic, prognostic and therapeutic target value for OC. We also consider the connection between circRNAs and OC cell proliferation, apoptosis, metastasis, and chemotherapy resistance and sensitivity.

Inflammation-Related circRNA Polymorphism and Ischemic Stroke Prognosis

CircRNAs belong to a novel class of noncoding RNAs that are generated by exons of genes by alternative mRNA splicing and involved in pathophysiological processes of ischemic stroke by regulating neuro-inflammation. A total of 982 patients were enrolled in our study for stroke recovery analysis. The aim of our study was to first explore the association between the inflammation-related circRNA polymorphism and functional outcome 3 months after ischemic stroke by using multivariate logistic regression model. Next, we further investigated the role of circRNA polymorphism in predicting stroke recurrence by using Cox proportional hazard regression model. Five circRNA polymorphisms were genotyped by using polymerase chain reaction and ligation detection reaction method. We identified circ-STAT3 (signal transducer and activator of transcription) rs2293152 GG genotype to be associated with poorer recovery 90 days after stroke (OR = 1.452, 95% CI: 1.165-4.362, p = 0.016). After adjusting for confound factors, the association for rs2293152 with 3 months outcome after IS was stronger, suggesting a mechanism that rs2293152 is an independent risk factor for stroke recovery (OR = 2.255, 95% CI: 1.034-2.038, p = 0.031). However, no other circRNA polymorphisms (circ-DLGAP4 rs41274714, circ-TRAF2 rs10870141, circ-ITCH rs10485505, rs4911154) were associated with functional outcome 3 months after stroke in any genetic models. Subgroup analysis revealed that the negative effect of rs2293152 GG genotype was greater in female and older patients, subjects with history of hypertension. Additionally, all the circRNA polymorphisms were not correlated with recurrent risk of ischemic stroke. Our results indicated that circ-STAT3 might be a novel biomarker for predicting functional outcome after stroke and an important contributor to the ischemic stroke recovery.

Upregulation of circ_0000142 promotes multiple myeloma progression by adsorbing miR-610 and upregulating AKT3 expression

Circular RNAs (circRNAs) play an important regulatory role in a variety of malignancies. Nevertheless, the role of circ_0000142 in multiple myeloma (MM) and its regulatory mechanism remains largely unknown. Real-time polymerase chain reaction was employed to detect the expressions of circ_0000142 and miR-610 in MM tissues and cell lines. The expression of AKT3 and apoptosis-related proteins (Bcl-2, Bax) in MM cells was detected by western blot. The correlation between the expression level of circ_0000142 and the clinicopathological parameters of MM patients was analysed. Cell proliferation, apoptosis, migration and invasion were monitored by Cell Counting Kit 8 assay, flow cytometry analysis and Transwell assay, respectively. The dual-luciferase reporter gene assay and RNA immunoprecipitation assay were employed to verify the targeting relationship between circ_0000142 and miR-610. In this study, it was demonstrated that, circ_0000142 was highly expressed in MM patients, and its high expression level was significantly associated with increased International Staging System and Durie-Salmon stage. Overexpression of circ_0000142 enhanced MM cell proliferation, migration, invasion and suppressed cell apoptosis, while knocking down circ_0000142 had the opposite effects. Mechanistically, circ_0000142 functioned as a competitive endogenous RNA, directly targeting miR-610 and positively regulating AKT3 expression. In brief, circ_0000142 enhances the proliferation and metastasis of MM cells by modulating the miR-610/AKT3 axis.

Function of Circular RNAs in Fish and Their Potential Application as Biomarkers

Circular RNAs (circRNAs) are an emerging class of regulatory RNAs with a covalently closed-loop structure formed during pre-mRNA splicing. Recent advances in high-throughput RNA sequencing and circRNA-specific computational tools have driven the development of novel approaches to their identification and functional characterization. CircRNAs are stable, developmentally regulated, and show tissue- and cell-type-specific expression across different taxonomic groups. They play a crucial role in regulating various biological processes at post-transcriptional and translational levels. However, the involvement of circRNAs in fish immunity has only recently been recognized. There is also broad evidence in mammals that the timely expression of circRNAs in muscle plays an essential role in growth regulation but our understanding of their expression and function in teleosts is still very limited. Here, we discuss the available knowledge about circRNAs and their role in growth and immunity in vertebrates from a comparative perspective, with emphasis on cultured teleost fish. We expect that the interest in teleost circRNAs will increase substantially soon, and we propose that they may be used as biomarkers for selective breeding of farmed fish, thus contributing to the sustainability of the aquaculture sector.

Insights Into the Role of CircRNAs: Biogenesis, Characterization, Functional, and Clinical Impact in Human Malignancies

Circular RNAs (circRNAs) are an evolutionarily conserved novel class of non-coding endogenous RNAs (ncRNAs) found in the eukaryotic transcriptome, originally believed to be aberrant RNA splicing by-products with decreased functionality. However, recent advances in high-throughput genomic technology have allowed circRNAs to be characterized in detail and revealed their role in controlling various biological and molecular processes, the most essential being gene regulation. Because of the structural stability, high expression, availability of microRNA (miRNA) binding sites and tissue-specific expression, circRNAs have become hot topic of research in RNA biology. Compared to the linear RNA, circRNAs are produced differentially by backsplicing exons or lariat introns from a pre-messenger RNA (mRNA) forming a covalently closed loop structure missing 3′ poly-(A) tail or 5′ cap, rendering them immune to exonuclease-mediated degradation. Emerging research has identified multifaceted roles of circRNAs as miRNA and RNA binding protein (RBP) sponges and transcription, translation, and splicing event regulators. CircRNAs have been involved in many human illnesses, including cancer and neurodegenerative disorders such as Alzheimer’s and Parkinson’s disease, due to their aberrant expression in different pathological conditions. The functional versatility exhibited by circRNAs enables them to serve as potential diagnostic or predictive biomarkers for various diseases. This review discusses the properties, characterization, profiling, and the diverse molecular mechanisms of circRNAs and their use as potential therapeutic targets in different human malignancies.

circVAR database: genome-wide archive of genetic variants for human circular RNAs

BACKGROUND

Circular RNAs (circRNAs) play important roles in regulating gene expression through binding miRNAs and RNA binding proteins. Genetic variation of circRNAs may affect complex traits/diseases by changing their binding efficiency to target miRNAs and proteins. There is a growing demand for investigations of the functions of genetic changes using large-scale experimental evidence. However, there is no online genetic resource for circRNA genes.

RESULTS

We performed extensive genetic annotation of 295,526 circRNAs integrated from circBase, circNet and circRNAdb. All pre-computed genetic variants were presented at our online resource, circVAR, with data browsing and search functionality. We explored the chromosome-based distribution of circRNAs and their associated variants. We found that, based on mapping to the 1000 Genomes and ClinVAR databases, chromosome 17 has a relatively large number of circRNAs and associated common and health-related genetic variants. Following the annotation of genome wide association studies (GWAS)-based circRNA variants, we found many non-coding variants within circRNAs, suggesting novel mechanisms for common diseases reported from GWAS studies. For cancer-based somatic variants, we found that chromosome 7 has many highly complex mutations that have been overlooked in previous research.

CONCLUSION

We used the circVAR database to collect SNPs and small insertions and deletions (INDELs) in putative circRNA regions and to identify their potential phenotypic information. To provide a reusable resource for the circRNA research community, we have published all the pre-computed genetic data concerning circRNAs and associated genes together with data query and browsing functions at http://soft.bioinfo-minzhao.org/circvar .

CircFOXO3 rs12196996, a polymorphism at the gene flanking intron, is associated with circFOXO3 levels and the risk of coronary artery disease

CircFOXO3 plays an important role in the pathogenesis of coronary artery disease (CAD). Single nucleotide polymorphisms (SNPs) at circRNA flanking introns may change its back-splicing and influence circRNA formation. Here, we aimed to investigate the influence of the polymorphisms at the circFOXO3 flanking introns on individual susceptibility to CAD. A total of 1185 individuals were included in the case-control study. In a multivariate logistic regression analysis, we determined that the rs12196996 G variant was significantly associated with increased CAD risk (OR = 1.36, P = 0.014). A similar trend of the association was observed in the recessive model (OR = 2.57, P = 0.003). Stratified analysis revealed a more significant association with CAD risk among younger subjects and non-smokers. Consistent with these results, the haplotype rs12196996G-rs9398171C containing rs12196996G allele was also associated with increased CAD risk (OR = 1.31, P = 0.013). Further investigation revealed that the rs12196996 GG genotype was associated with decreased circFOXO3 expression, but not linear FOXO3 levels. Taken together, our data provide the first evidence that the rs12196996 polymorphism at the circFOXO3 gene flanking intron is associated with CAD risk in the Chinese Han population, which is probably due to influence circFOXO3 levels.

Role of non-coding RNA networks in leukemia progression, metastasis and drug resistance

Early-stage detection of leukemia is a critical determinant for successful treatment of the disease and can increase the survival rate of leukemia patients. The factors limiting the current screening approaches to leukemia include low sensitivity and specificity, high costs, and a low participation rate. An approach based on novel and innovative biomarkers with high accuracy from peripheral blood offers a comfortable and appealing alternative to patients, potentially leading to a higher participation rate.

Recently, non-coding RNAs due to their involvement in vital oncogenic processes such as differentiation, proliferation, migration, angiogenesis and apoptosis have attracted much attention as potential diagnostic and prognostic biomarkers in leukemia. Emerging lines of evidence have shown that the mutational spectrum and dysregulated expression of non-coding RNA genes are closely associated with the development and progression of various cancers, including leukemia. In this review, we highlight the expression and functional roles of different types of non-coding RNAs in leukemia and discuss their potential clinical applications as diagnostic or prognostic biomarkers and therapeutic targets.

A Multi-Omics Perspective of Quantitative Trait Loci in Precision Medicine

Quantitative trait loci (QTL) analysis is an important approach to investigate the effects of genetic variants identified through an increasing number of large-scale, multidimensional 'omics data sets. In this 'big data' era, the research community has identified a significant number of molecular QTLs (molQTLs) and increased our understanding of their effects. Herein, we review multiple categories of molQTLs, including those associated with transcriptome, post-transcriptional regulation, epigenetics, proteomics, metabolomics, and the microbiome. We summarize approaches to identify molQTLs and to infer their causal effects. We further discuss the integrative analysis of molQTLs through a multi-omics perspective. Our review highlights future opportunities to better understand the functional significance of genetic variants and to utilize the discovery of molQTLs in precision medicine.