Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gene regulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

The regulatory and transcriptional landscape associated with carbon utilization in a filamentous fungus

Filamentous fungi, such as Neurospora crassa , are very efficient in deconstructing plant biomass by the secretion of an arsenal of plant cell wall-degrading enzymes, by remodeling metabolism to accommodate production of secreted enzymes, and by enabling transport and intracellular utilization of plant biomass components. Although a number of enzymes and transcriptional regulators involved in plant biomass utilization have been identified, how filamentous fungi sense and integrate nutritional information encoded in the plant cell wall into a regulatory hierarchy for optimal utilization of complex carbon sources is not understood. Here, we performed transcriptional profiling of N. crassa on 40 different carbon sources, including plant biomass, to provide data on how fungi sense simple to complex carbohydrates. From these data, we identified regulatory factors in N. crassa and characterized one (PDR-2) associated with pectin utilization and one with pectin/hemicellulose utilization (ARA-1). Using in vitro DNA affinity purification sequencing (DAP-seq), we identified direct targets of transcription factors involved in regulating genes encoding plant cell wall-degrading enzymes. In particular, our data clarified the role of the transcription factor VIB-1 in the regulation of genes encoding plant cell wall-degrading enzymes and nutrient scavenging and revealed a major role of the carbon catabolite repressor CRE-1 in regulating the expression of major facilitator transporter genes. These data contribute to a more complete understanding of cross talk between transcription factors and their target genes, which are involved in regulating nutrient sensing and plant biomass utilization on a global level.

59 BASIC BIOLOGICAL SCIENCES↗

Enhanced Resistance Pines for Improved Renewable Biofuel and Chemical Production (Technical Report)

We completed phenotyping constitutive and inducible oleoresin flow across two seasons, constitutive resin canal number and density and wood terpene content in our ADEPT2 and CCLONES populations. We completed genetic association between 19 oleoresin phenotypes and a total of 523,192 SNP markers from ADEPT2 and 13,883 SNP markers in CCLONES using four mixed linear models. A total of 293 significant SNPs (FDR = 0.20) were identified. We used the MENTOR tool to mine mechanistic connections from a multiplex network constructed from poplar multi-omic data to construct a conceptual model for a subset of these significant SNPs. Our model contains 6 transcriptional regulators in addition to 3 monoterpene synthases. To generate more lines of evidence for these significant SNPs, we completed a time course RNAseq experiment after inducing vascular zone cells to differentiate into new resin canals with a methyl jasmonate treatment, a single nuclei RNAseq that identified differentiating resin canal epithelial cells and are completing analysis for a QTL study in a hybrid pine population. The time course identified 4634 significantly down and 1890 significantly up regulated transcripts after treatment with methyl jasmonate, an inducer of new resin canal formation in the vascular cambial meristem. To analyze this large set of differentially regulated genes, we created a predictive expression network and analyzed it with random walk restart using 6 seed genes coding for transcription factors regulating xylem differentiation in poplar. Of the top ranked 200 transcripts, 119 transcripts were significant differentially expressed supporting these transcripts as potential candidates regulating resin canal formation. Analysis of single nuclei sequencing of shoot tips that contain differentiating resin canals, identified 10 clusters. One cluster was highly enriched in transcripts coding for 9 of the enzymes in the MEP pathway 3 prenyl synthetases, and 3 monoterpene synthases strongly suggesting that this cluster represents resin canal epithelial cells. We are mining the additional transcripts to create a trajectory analysis. In summary, we have identified > 10 novel genes that are strongly supported candidates for further analysis in breeding lines and for genetic engineering over- and under- expressing lines to increase wood terpene content to improve resistance to insect and fungal pathogens while simultaneously increasing terpene supplies for renewable chemicals and biofuels.

59 BASIC BIOLOGICAL SCIENCES↗

A Powerful Toolkit for Synthetic Biology: Over 3.8 Billion Years of Evolution

The combination of evolutionary with engineering principles will enhance synthetic biology. Conversely, synthetic biology has the potential to enrich evolutionary biology by explaining why some adaptive space is empty, on Earth or elsewhere. Synthetic biology, the design and construction of artificial biological systems, substitutes bio-engineering for evolution, which is seen as an obstacle. But because evolution has produced the complexity and diversity of life, it provides a proven toolkit of genetic materials and principles available to synthetic biology. Evolution operates on the population level, with the populations composed of unique individuals that are historical entities. The source of genetic novelty includes mutation, gene regulation, sex, symbiosis, and interspecies gene transfer. At a phenotypic level, variation derives from regulatory control, replication and diversification of components, compartmentalization, sexual selection and speciation, among others. Variation is limited by physical constraints such as diffusion, and chemical constraints such as reaction rates and membrane fluidity. While some of these tools of evolution are currently in use in synthetic biology, all ought to be examined for utility. A hybrid approach of synthetic biology coupled with fine-tuning through evolution is suggested

Rothschild, Lynn J.↗

CovS inactivation reduces CovR promoter binding at diverse virulence factor encoding genes in group A Streptococcus

The control of virulence gene regulator (CovR), also called caspsule synthesis regulator (CsrR), is critical to how the major human pathogen group A Streptococcus fine-tunes virulence factor production. CovR phosphorylation (CovR~P) levels are determined by its cognate sensor kinase CovS, and functional abrogating mutations in CovS can occur in invasive GAS isolates leading to hypervirulence. Presently, the mechanism of CovR-DNA binding specificity is unclear, and the impact of CovS inactivation on global CovR binding has not been assessed. Thus, we performed CovR chromatin immunoprecipitation sequencing (ChIP-seq) analysis in the emm1 strain MGAS2221 and its CovS kinase deficient derivative strain 2221-CovS-E281A. We identified that CovR bound in the promoter regions of nearly all virulence factor encoding genes in the CovR regulon. Additionally, direct CovR binding was observed for numerous genes encoding proteins involved in amino acid metabolism, but we found limited direct CovR binding to genes encoding other transcriptional regulators. The consensus sequence AATRANAAAARVABTAAA was present in the promoters of genes directly regulated by CovR, and mutations of highly conserved positions within this motif relieved CovR repression of the hasA and MGAS2221_0187 promoters. Analysis of strain 2221-CovS-E281A revealed that binding of CovR at repressed, but not activated, promoters is highly dependent on CovR~P state. CovR repressed virulence factor encoding genes could be grouped dependent on how CovR~P dependent variation in DNA binding correlated with gene transcript levels. Taken together, the data show that CovR repression of virulence factor encoding genes is primarily direct in nature, involves binding to a newly-identified DNA binding motif, and is relieved by CovS inactivation. These data provide new mechanistic insights into one of the most important bacterial virulence regulators and allow for subsequent focused investigations into how CovR-DNA interaction at directly controlled promoters impacts GAS pathogenesis.

59 BASIC BIOLOGICAL SCIENCES↗

A Foundational Study for Normal F8 -Containing Mouse Models for the miRNA Regulation of Hemophilia A: Identification and Analysis of Mouse miRNAs that Downregulate the Murine F8 Gene

Hemophilia A (HA) is associated with defects in the F8 gene, encoding coagulation factor VIII (FVIII). Our previous studies show that F8-targeting micro RNAs (miRNAs), a group of small RNAs involved in gene regulation, can downregulate F8 expression causing HA in individuals with normal F8-genotypes and increased HA severity in patients with mutations in F8. Understanding the mechanistic underpinnings of human genetic diseases caused or modulated by miRNAs require a small animal model, such as a mouse model. Here, we report a foundational study to develop such a model system. We identified the mouse 3'untranslated region (3'UTR) on murine F8-mRNA (muF8-mRNA) that can bind to murine miRNAs. We then selected three miRNAs for evaluation: miR-208a, miR-351 and miR-125a. We first demonstrate that these three miRNAs directly target the 3'UTR of muF8-mRNA and reduce the expression of a reporter gene (luciferase) mRNA fused to the muF8-3' UTR in mammalian cells. Furthermore, in mouse cells that endogenously express the F8 gene and produce FVIII protein, the ectopic expression of these miRNAs downregulated F8-mRNA and FVIII protein. These results provide proof-of-concept and reagents as a foundation for using a normal F8-containing mouse as a model for the miRNA regulation of normal F8 in causing or aggravating the genetic disease HA.

59 BASIC BIOLOGICAL SCIENCES↗

Regulatory helix plays a key role in genetic ON-OFF switching for the 2’-deoxyguanosine sensing mRNA element

Transcriptional riboswitches, noncoding mRNA elements that operate in cis to regulate gene expression, have a promising potential in medicine, synthetic biology and directed evolution. They bind to cellular metabolites or metal ions with high specificity, leading to conformational rearrangements that facilitate the activation or premature termination of transcription for downstream genes. This elegant mechanism for feedback regulation of metabolic pathways has been identified in prokaryotes and a few in eukaryotes. Our chemical probing of the 2’-deoxyguanosine (2’-dG)-sensing riboswitch demonstrates that the overall conformational state of the full-length riboswitch (dGsw-fl) is unresponsive to the 2’-dG. Although binding proceeds as expected, dGsw-fl exclusively populates an OFF state of transcriptional inhibition. We chemically probed the structure of a known dGsw transcriptional intermediate (dGsw-int) to evaluate the possibility of a cotranscriptional regulatory role. Interestingly, apo dGsw-int adopts an alternative conformation in which a stable anti-terminator helix is formed, leading to an ON state where transcription can proceed. In the presence of 2’-dG, this anti-terminator helix is destabilized to produce a conformation reminiscent of the full-length, OFF-state dGsw. Using a fluorescence quenching assay, we demonstrate that binding 2’-dG to early transcriptional intermediates can inhibit the formation of the anti-terminator helix, locking dGsw in an OFF state. These data suggest that metabolite sensing occurs during a brief window of time between the synthesis of two transcriptional intermediates. Our studies indicate that dGsw does not function as a binary ON−OFF switch, but instead fine-tunes the transcription of downstream genes during RNA synthesis using key intermediates.

59 BASIC BIOLOGICAL SCIENCES↗

The multifaceted role of c-di-AMP signaling in the regulation of Porphyromonas gingivalis lipopolysaccharide structure and function

This study unveils the intricate functional association between cyclic di-3’,5’-adenylic acid (c-di-AMP) signaling, cellular bioenergetics, and the regulation of lipopolysaccharide (LPS) profile in Porphyromonas gingivalis, a Gram-negative obligate anaerobe considered as a keystone pathogen involved in the pathogenesis of chronic periodontitis. Previous research has identified variations in P. gingivalis LPS profile as a major virulence factor, yet the underlying mechanism of its modulation has remained elusive. We employed a comprehensive methodological approach, combining two mutants exhibiting varying levels of c-di-AMP compared to the wild type, alongside an optimized analytical methodology that combines conventional mass spectrometry techniques with a novel approach known as FLAT n . We demonstrate that c-di-AMP acts as a metabolic nexus, connecting bioenergetic status to nuanced shifts in fatty acid and glycosyl profiles within P. gingivalis LPS. Notably, the predicted regulator gene cdaR, serving as a potent regulator of c-di-AMP synthesis, was found essential for producing N-acetylgalactosamine and an unidentified glycolipid class associated with the LPS profile. The multifaceted roles of c-di-AMP in bacterial physiology are underscored, emphasizing its significance in orchestrating adaptive responses to stimuli. Furthermore, our findings illuminate the significance of LPS variations and c-di-AMP signaling in determining the biological activities and immunostimulatory potential of P. gingivalis LPS, promoting a pathoadaptive strategy. The study expands the understanding of c-di-AMP pathways in Gram-negative species, laying a foundation for future investigations into the mechanisms governing variations in LPS structure at the molecular level and their implications for host-pathogen interactions.

59 BASIC BIOLOGICAL SCIENCES↗

Developing multi-gene CRISPRa/i programs to accelerate DBTL cycles in ABF hosts engineered for chemical production

This project developed and implemented a modular CRISPR activation and interference (CRISPRa/i) platform to accelerate strain optimization and pathway development for industrially relevant microbial hosts. By integrating multiplexed transcriptional perturbation tools with data-driven Design–Build–Test–Learn (DBTL) workflows, the team achieved reductions in cycle time and enhanced production of industrial aromatics, particularly 4-aminocinnamic acid (4-ACA), in Pseudomonas putida. Key accomplishments included: ● Development of a robust, tunable CRISPRa/i system in P. putida that enabled efficient multi-target gene regulation via guide RNA (gRNA) programs ● Completion of two full DBTL cycles, guided by machine learning (ML) models trained on transcriptomic and performance data, reducing engineering time by over 30% ● Optimization of multi-gene regulatory programs to balance expression of host and pathway modules, improve 4-ACA titers, and resolve metabolic bottlenecks ● Demonstration of system portability through a limited proof-of-concept extension in Acinetobacter baylyi, underscoring the generalizability of the approach ● Evaluation of strain performance on lignocellulosic biomass-derived substrates, demonstrating the feasibility of converting renewable carbon into aromatic building blocks These results illustrate the feasibility of applying ML-guided CRISPRa/i perturbation strategies to accelerate strain development in complex microbial systems. The resulting tools and datasets contribute to DOE objectives by improving platform predictability, reducing development costs, and enabling broader access to sustainable, economically viable bioproduction technologies.

09 BIOMASS FUELS↗

Integrative genomics reveals paths to sex dimorphism in Salix purpurea L

Sex dimorphism and gene expression were studied in developing catkins in 159 F 2 individuals from the bioenergy crop Salix purpurea , and potential mechanisms and pathways for regulating sex development were explored. Differential expression, eQTL, bisulfite sequencing, and network analysis were used to characterize sex dimorphism, detect candidate master regulator genes, and identify pathways through which the sex determination region (SDR) may mediate sex dimorphism. Eleven genes are presented as candidates for master regulators of sex, supported by gene expression and network analyses. These include genes putatively involved in hormone signaling, epigenetic modification, and regulation of transcription. eQTL analysis revealed a suite of transcription factors and genes involved in secondary metabolism and floral development that were predicted to be under direct control of the sex determination region. Furthermore, data from bisulfite sequencing and small RNA sequencing revealed strong differences in expression between males and females that would implicate both of these processes in sex dimorphism pathways. These data indicate that the mechanism of sex determination in Salix purpurea is likely different from that observed in the related genus Populus . This further demonstrates the dynamic nature of SDRs in plants, which involves a multitude of mechanisms of sex determination and a high rate of turnover.

59 BASIC BIOLOGICAL SCIENCES↗

Life: Origin and evolution on Earth--How can we escape?

Exploitation of gene regulation rather than the creation of new genes has been predominantly responsible for the evolutionary advances in animals and plants that are widely recognized today. Until very recently it was not possible to examine life in the absence of gravity. We can now imagine forms of life in the universe adapting to circumstances different from those found on Earth. Our own life forms would surely become different in time if they were transferred to other planets with different conditions, including much lower or higher gravity.

Markert, Clement L.↗

Transcriptomic and Proteomic Insights into Host Immune Responses in Pediatric Severe Malarial Anemia: Dysregulation in HSP60-70-TLR2/4 Signaling and Altered Glutamine Metabolism

Severe malarial anemia (SMA, Hb < 6.0 g/dL) is a leading cause of childhood morbidity and mortality in holoendemic Plasmodium falciparum transmission zones. This study explored the entire expressed human transcriptome in whole blood from 66 Kenyan children with non-SMA (Hb ≥ 6.0 g/dL, n = 41) and SMA (n = 25), focusing on host immune response networks. RNA-seq analysis revealed 6862 differentially expressed genes, with equally distributed up-and down-regulated genes, indicating a complex host immune response. Deconvolution analyses uncovered leukocytic immune profiles indicative of a diminished antigenic response, reduced immune priming, and polarization toward cellular repair in SMA. Weighted gene co-expression network analysis revealed that immune-regulated processes are central molecular distinctions between non-SMA and SMA. A top dysregulated immune response signaling network in SMA was the HSP60-HSP70-TLR2/4 signaling pathway, indicating altered pathogen recognition, innate immune activation, stress responses, and antigen recognition. Validation with high-throughput gene expression from a separate cohort of Kenyan children (n = 50) with varying severities of malarial anemia (n = 38 non-SMA and n = 12 SMA) confirmed the RNA-seq findings. Proteomic analyses in 35 children with matched transcript and protein abundance (n = 19 non-SMA and n = 16 SMA) confirmed dysregulation in the HSP60-HSP70-TLR2/4 signaling pathway. Additionally, glutamine transporter and glutamine synthetase genes were differentially expressed, indicating altered glutamine metabolism in SMA. This comprehensive analysis underscores complex immune dysregulation and novel pathogenic features in SMA.

Microbiology↗

A CRISPRi-dCas9 System for Archaea and Its Use To Examine Gene Function during Nitrogen Fixation by Methanosarcina acetivorans

CRISPR-based systems are emerging as the premier method to manipulate many cellular processes. In this study, a simple and efficient CRISPR interference (CRISPRi) system for targeted gene repression in archaea was developed. The Methanosarcina acetivorans CRISPR-Cas9 system was repurposed by replacing Cas9 with the catalytically dead Cas9 (dCas9) to generate a CRISPRi-dCas9 system for targeted gene repression. To test the utility of the system, genes involved in nitrogen (N2) fixation were targeted for dCas9-mediated repression. First, the nif operon (nifHI 1 I 2 DKEN) that encodes molybdenum nitrogenase was targeted by separate guide RNAs (gRNAs), one targeting the promoter and the other targeting nifD. Remarkably, growth of M. acetivorans with N 2 was abolished by dCas9-mediated repression of the nif operon with each gRNA. The abundance of nif transcripts was >90% reduced in both strains expressing the gRNAs, and NifD was not detected in cell lysate. Next, we targeted NifB, which is required for nitrogenase cofactor biogenesis. Expression of a gRNA targeting the coding sequence of NifB decreased nifB transcript abundance >85% and impaired but did not abolish growth of M. acetivorans with N2. Finally, to ascertain the ability to study gene regulation using CRISPRi-dCas9, nrpR1, encoding a subunit of the repressor of the nif operon, was targeted. The nrpR1 repression strain grew normally with N2 but had increased nif operon transcript abundance, consistent with NrpR1 acting as a repressor. Furthermore, these results highlight the utility of the system, whereby a single gRNA when expressed with dCas9 can block transcription of targeted genes and operons in M. acetivorans.

59 BASIC BIOLOGICAL SCIENCES↗

Changes in DNA Methylation in Response to 6-Benzylaminopurine Affect Allele-Specific Gene Expression in Populus Tomentosa

Cytokinins play important roles in the growth and development of plants. Physiological and photosynthetic characteristics are common indicators to measure the growth and development in plants. However, few reports have described the molecular mechanisms of physiological and photosynthetic changes in response to cytokinin, particularly in woody plants. DNA methylation is an essential epigenetic modification that dynamically regulates gene expression in response to the external environment. In this study, we examined genome-wide DNA methylation variation and transcriptional variation in poplar (Populus tomentosa) after short-term treatment with the synthetic cytokinin 6-benzylaminopurine (6-BA). We identified 460 significantly differentially methylated regions (DMRs) in response to 6-BA treatment. Transcriptome analysis showed that 339 protein-coding genes, 262 long non-coding RNAs (lncRNAs), and 15,793 24-nt small interfering RNAs (siRNAs) were differentially expressed under 6-BA treatment. Among these, 79% were differentially expressed between alleles in P. tomentosa, and 102,819 allele-specific expression (ASE) loci in 19,200 genes were detected showing differences in ASE levels after 6-BA treatment. Combined DNA methylation and gene expression analysis demonstrated that DNA methylation plays an important role in regulating allele-specific gene expression. To further investigate the relationship between these 6-BA-responsive genes and phenotypic variation, we performed SNP analysis of 460 6-BA-responsive DMRs via re-sequencing using a natural population of P. tomentosa, and we identified 206 SNPs that were significantly associated with growth and wood properties. Association analysis indicated that 53% of loci with allele-specific expression had primarily dominant effects on poplar traits. Our comprehensive analyses of P. tomentosa DNA methylation and the regulation of allele-specific gene expression suggest that DNA methylation is an important regulator of imbalanced expression between allelic loci.

54 ENVIRONMENTAL SCIENCES↗

Adaptive laboratory evolution and metabolic engineering of Cupriavidus necator for improved catabolism of volatile fatty acids

Bioconversion of high-volume waste streams into value-added products will be an integral component of the growing bioeconomy. Volatile fatty acids (VFAs) (e.g., butyrate, valerate, and hexanoate) are an emerging and promising waste-derived feedstock for microbial carbon upcycling. Cupriavidus necator H16 is a favorable host for conversion of VFAs into various bioproducts due to its diverse carbon metabolism, ease of metabolic engineering, and use at industrial scales. Here, in this study, we report that a common strategy to improve product titers in C. necator, deletion of the polyhydroxybutyrate (PHB) biosynthetic operon, results in a significant growth defect on VFA substrates. Using adaptive laboratory evolution, we identify mutations to the regulator gene phaR, the two-component response regulator-histidine kinase pair encoded by H16_A1372/H16_A1373, and the tripartite transporter assembly encoded by H16_A2296-A2298 as causative for improved growth on VFA substrates. Deletion of phaR and H16_A1373 led to significantly reduced NADH abundance accompanied by large changes to expression of genes involved in carbon metabolism, balance of electron carriers, and oxidative stress tolerance that may be responsible for improved growth of these engineered strains. These results provide insight into the role of PHB biosynthesis in carbon and energy metabolism and highlight a key role for the regulator PhaR in global regulatory networks. By combining mutations, we generated platform strains with significant growth improvements on VFAs, which can enable improved conversion of waste-derived VFA substrates to target bioproducts.

09 BIOMASS FUELS↗

Spatiotemporal DNA methylome dynamics of the developing mouse fetus

Cytosine DNA methylation is essential for mammalian development but understanding of its spatiotemporal distribution in the developing embryo remains limited1,2. Here, as part of the mouse Encyclopedia of DNA Elements (ENCODE) project, we profiled 168 methylomes from 12 mouse tissues or organs at 9 developmental stages from embryogenesis to adulthood. We identified 1,808,810 genomic regions that showed variations in CG methylation by comparing the methylomes of different tissues or organs from different developmental stages. These DNA elements predominantly lose CG methylation during fetal development, whereas the trend is reversed after birth. During late stages of fetal development, non-CG methylation accumulated within the bodies of key developmental transcription factor genes, coinciding with their transcriptional repression. Integration of genome-wide DNA methylation, histone modification and chromatin accessibility data enabled us to predict 461,141 putative developmental tissue-specific enhancers, the human orthologues of which were enriched for disease-associated genetic variants. These spatiotemporal epigenome maps provide a resource for studies of gene regulation during tissue or organ progression, and a starting point for investigating regulatory elements that are involved in human developmental disorders.

59 BASIC BIOLOGICAL SCIENCES↗

miRNAs: Critical mediators of breast cancer metastatic programming

Highlights: • This review highlights the crucial role played by miRNAs in the process of metastatic evolution of breast cancer. MicroRNA mediated aberrant gene regulation has been implicated in several diseases including cancer. Recent research has highlighted the role of epigenetic modulation of the complex process of breast cancer metastasis by miRNAs. miRNAs play a crucial role in the process of metastatic evolution by facilitating alterations in the phenotype of tumor cells and the tumor microenvironment that promote this process. They act as critical determinants of the multi-step progression starting from carcinogenesis all the way to organotropism. In this review, we focus on the current understanding of the compelling role of miRNAs in breast cancer metastasis.

60 APPLIED LIFE SCIENCES↗

Noncoding variants alter GATA2 expression in rhombomere 4 motor neurons and cause dominant hereditary congenital facial paresis

Hereditary congenital facial paresis type 1 (HCFP1) is an autosomal dominant disorder of absent or limited facial movement that maps to chromosome 3q21-q22 and is hypothesized to result from facial branchial motor neuron (FBMN) maldevelopment. In the present study, we report that HCFP1 results from heterozygous duplications within a neuron-specific GATA2 regulatory region that includes two enhancers and one silencer, and from noncoding single-nucleotide variants (SNVs) within the silencer. Some SNVs impair binding of NR2F1 to the silencer in vitro and in vivo and attenuate in vivo enhancer reporter expression in FBMNs. Gata2 and its effector Gata3 are essential for inner-ear efferent neuron (IEE) but not FBMN development. A humanized HCFP1 mouse model extends Gata2 expression, favors the formation of IEEs over FBMNs and is rescued by conditional loss of Gata3. These findings highlight the importance of temporal gene regulation in development and of noncoding variation in rare mendelian disease.

59 BASIC BIOLOGICAL SCIENCES↗

Long- and short-read sequencing methods discover distinct circular RNA pools in Lotus japonicus

Circular RNAs (circRNAs) are covalently closed single-stranded RNAs, generated through a back-splicing process that links a downstream 5′ site to an upstream 3′ end. The only distinction in the sequence between circRNA and their linear cognate RNA is the back splice junction. Their low abundance and sequence similarity with their linear origin RNA have made the discovery and identification of circRNA challenging. We have identified almost 6000 novel circRNAs from Lotus japonicus leaf tissue using different enrichment, amplification, and sequencing methods as well as alternative bioinformatics pipelines. The different methodologies identified different pools of circRNA with little overlap. We validated circRNA identified by the different methods using reverse transcription polymerase chain reaction and characterized sequence variations using nanopore sequencing. We compared validated circRNA identified in L. japonicus to other plant species and showed conservation of high-confidence circRNA-expressing genes. This is the first identification of L. japonicus circRNA and provides a resource for further characterization of their function in gene regulation. CircRNAs identified in this study originated from genes involved in all biological functions of eukaryotic cells. The comparison of methodologies and technologies to sequence, identify, analyze, and validate circRNA from plant tissues will enable further research to characterize the function and biogenesis of circRNA in L. japonicus.

59 BASIC BIOLOGICAL SCIENCES↗