RNA Editing Sites Research Articles

Complete Organelle Genome of the Desiccation-Tolerant (DT) Moss Tortula atrovirens and Comparative Analysis of the Pottiaceae Family.

Tortula atrovirens (Sm.) Lindb. is an important component of biological soil crusts and possesses an extraordinary tolerance against desiccation in dryland habitats. However, knowledge of the organelle genome of this desiccation-tolerant (DT) moss is still lacking. Here, we assembled the first reported Tortula organelle genome and conducted a comprehensive analysis within the Pottiaceae family. T. atrovirens exhibited the second largest chloroplast genome (129,646 bp) within the Pottiaceae, whereas its mitogenome (105,877 bp) and those of other mosses were smaller in size compared to other land plants. The chloroplast and mitochondrial genomes of T. atrovirens were characterized by the expansion of IR boundaries and the absence of homologous recombination-mediated by large repeats. A total of 57 RNA editing sites were detected through mapping RNA-seq data. Moreover, the gene content and order were highly conserved among the Pottiaceae organelle genomes. Phylogenetic analysis showed that bryophytes are paraphyletic, with their three lineages (hornworts, mosses, and liverworts) and vascular plants forming successive sister clades. Timmiella anomala is clearly separated from the monophyletic Pottiaceae, and T. atrovirens is closely related to Syntrichia filaris within the Pottioideae. In addition, we detected four hypervariable regions for candidate-molecular markers. Our findings provide valuable insights into the organelle genomes of T. atrovirens and the evolutionary relationships within the Pottiaceae family, facilitating future discovery of DT genetic resources from bryophytes.

Assembly and comparative analysis of the first complete mitochondrial genome of Setaria italica.

In this study, we assembled the first complete mitochondrial genome of Setaria italica and confirmed the multi-branched architecture. The foxtail millet (Setaria italica) holds significant agricultural importance, particularly in arid and semi-arid regions. It plays a pivotal role in diversifying dietary patterns and shaping planting strategies. Although the chloroplast genome of S. italica has been elucidated in recent studies, the complete mitochondrial genome remains largely unexplored. In this study, we employed PacBio HiFi sequencing platforms to sequence and assemble the complete mitochondrial genome. The mitochondrial genome spans a total length of 446,614 base pairs and harbors a comprehensive set of genetic elements, including 33 unique protein-coding genes (PCGs), encompassing 24 unique mitochondrial core genes and 9 variable genes, along with 20 transfer RNA (tRNA) genes and 3 ribosomal RNA (rRNA) genes. Our analysis of mitochondrial PCGs revealed a pronounced codon usage preference. For instance, the termination codon exhibits a marked preference for UAA, while alanine (Ala) exhibits a preference for GCU, and glutamine (Gln) favors CAA. Notably, the maximum Relative Synonymous Codon Usage (RSCU) values for cysteine (Cys) and phenylalanine (Phe) are both below 1.2, indicating a lack of strong codon usage preference for these amino acids. Phylogenetic analyses consistently place S. italica in close evolutionary proximity to Chrysopogon zizanioides, relative to other Panicoideae plants. Collinearity analysis showed that a total of 39 fragments were identified to display homology with both the mitochondrial and chloroplast genomes. A total of 417 potential RNA-editing sites were discovered across the 33 mitochondrial PCGs. Notably, all these editing events involved the conversion of cytosine (C) to uracil (U). Through the employment of PCR validation coupled with Sanger sequencing for the anticipated editing sites of these codons, RNA-editing events were conclusively identified at two specific loci: nad4L-2 and atp6-1030. The results of this study provide a pivotal foundation for advanced genomic breeding research in foxtail millet. Furthermore, they impart essential insights that will be instrumental for forthcoming investigations into the evolutionary and molecular dynamics of Panicoideae species.

Full-length PacBio Amplicon Sequencing to Unveil RNA Editing Sites

Background: RNA editing enriches post-transcriptional sequence changes. Currently detecting RNA editing sites is mostly based on the Sanger sequencing platform and second-generation sequencing. However, detection with Sanger sequencing is limited by the disturbing background peaks using the direct sequencing method and the clone number using the clone sequencing method, while second-generation sequencing detection is constrained by its short read. Objective: We aimed to design a pipeline that can accurately detect RNA editing sites for full-length long-read amplicons to meet the requirement when focusing on a few specific genes of interest. Method: We developed a novel high-throughput RNA editing sites detection pipeline based on the PacBio circular consensus sequences sequencing which is accurate with high-throughput and long-read coverage. We tested the pipeline on cytosolic malate dehydrogenase in the hard-shelled mussel Mytilus coruscus and further validated it using direct Sanger sequencing. Results: Data generated from the PacBio circular consensus sequences (CCS) amplicons in three mussels were first filtered by quality and then selected by open reading frame. After filtering, 225-2047 sequences of the three mussels, respectively, were used to identify RNA editing sites. With corresponding genomic DNA sequences, we extracted 227-799 candidate RNA editing sites excluding heterozygous sites. We further figured out 7-11 final RESs using a new error model specially designed for RNA editing site detection. The resulting RNA editing sites all agree with the validation using the Sanger sequencing. Conclusion: We report a near-zero error rate method in identifying RNA editing sites of long-read amplicons with the use of PacBio CCS sequencing.

De Novo Assembly and Comparative Analysis of Mitochondrial Genomes of Two Pueraria montana Varieties.

Pueraria montana is a species with important medicinal value and a complex genetic background. In this study, we sequenced and assembled the mitochondrial (mt) genomes of two varieties of P. montana. The mt genome lengths of P. montana var. thomsonii and P. montana var. montana were 457,390 bp and 456,731 bp, respectively. Both P. montana mitogenomes showed a multi-branched structure consisting of two circular molecules, with 56 genes annotated, comprising 33 protein-coding genes, 18 tRNA genes (trnC-GCA and trnM-CAU are multi-copy genes), and 3 rRNA genes. Then, 207 pairs of long repeats and 96 simple sequence repeats (SSRs) were detected in the mt genomes of P. montana, and 484 potential RNA-editing sites were found across the 33 mitochondrial protein-coding genes of each variety. Additionally, a syntenic sequence analysis showed a high collinearity between the two mt genomes. This work is the first to analyze the mt genomes of P. montana. It can provide information that can be used to analyze the structure of mt genomes of higher plants and provide a foundation for future comparative genomic studies and evolutionary biology research in related species.

Comprehensive analysis of the mitochondrial genome of Rehmannia glutinosa: insights into repeat-mediated recombinations and RNA editing-induced stop codon acquisition.

Rehmannia glutinosa is an economically significant medicinal plant. Yet, the structure and sequence of its mitochondrial genome has not been published, which plays a crucial role in evolutionary analysis and regulating respiratory-related macromolecule synthesis. In this study, the R. glutinosa mitogenome was sequenced employing a combination of Illumina short reads and Nanopore long reads, with subsequent assembly using a hybrid strategy. We found that the predominant configuration of the R. glutinosa mitogenome comprises two circular chromosomes. The primary structure of the mitogenome encompasses two mitochondrial chromosomes corresponding to the two major configurations, Mac1-1 and Mac1-2. The R. glutinosa mitogenome encoded an angiosperm-typical set of 24 core genes, nine variable genes, three rRNA genes, and 15 tRNA genes. A phylogenetic analysis using the 16 shared protein-coding genes (PCG) yielded a tree consistent with the phylogeny of Lamiales species and two outgroup taxa. Mapping RNA-seq data to the coding sequences (CDS) of the PCGs revealed 507 C-to-U RNA editing sites across 31 PCGs of the R. glutinosa mitogenome. Furthermore, one start codon (nad4L) and two stop codons (rpl10 and atp6) were identified as products of RNA editing events in the R. glutinosa mitogenome.

Mitogenome Assembly Reveals Gene Migration and RNA Editing Events in Plateau Hongliu (Myricaria elegans Royle.)

The Plateau Hongliu (Myricaria elegans Royle.) is a woody shrub halophyte that thrives in arid areas of western Tibet, in the Himalayan Mountains. It is acclaimed as superior in saline stress acclimation and as a critical pharmaceutical resource of the Tibetan traditional herb. Nevertheless, the mitogenome in the genus Myricaria remains unknown. Here, using the Illumina and PacBio sequencing assays, the first complete mitogenome of the M. elegans revealed a multi-branched skeleton with a total length of 416,354 bp and GC content of 44.33%, comprising two circular molecules (M1 and 2). The complete mitogenome annotates 31 unique protein-encoding genes (PEGs), fifteen tRNAs, and three rRNA genes. The UAA exhibits the most prominent codon usage preference as a termination, followed by UUA codons for leucine. The mitogenome contains 99 simple sequence repeats and 353 pairs of dispersed repeats, displaying the most frequent in palindromic repeats. Gene transfer analyses identified 8438 bp of 18 homologous fragments from the plastome, accounting for 2.03% of the total length. Using the PREP suite, 350 C-U RNA editing sites were predicted, of which nad4 and ccmB were on the top frequency. Syntenic and phylogenetic analyses suggested weakly conserved patterns of M. elegans in Caryophyllales owing to the genome rearrangement. In summary, the deciphered unique features and complexities of the mitogenome in M. elegans provide novel insights into understanding the evolution and biological conservation underlying climate resilience in halophytes.

Development and assessment of an RNA editing-based risk model for the prognosis of cervical cancer patients.

RNA editing, as an epigenetic mechanism, exhibits a strong correlation with the occurrence and development of cancers. Nevertheless, few studies have been conducted to investigate the impact of RNA editing on cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC). In order to study the connection between RNA editing and CESC patients' prognoses, we obtained CESC-related information from The Cancer Genome Atlas (TCGA) database and randomly allocated the patients into the training group or testing group. An RNA editing-based risk model for CESC patients was established by Cox regression analysis and least absolute shrinkage and selection operator (LASSO). According to the median score generated by this RNA editing-based risk model, patients were categorized into subgroups with high and low risks. We further constructed the nomogram by risk scores and clinical characteristics and analyzed the impact of RNA editing levels on host gene expression levels and adenosine deaminase acting on RNA. Finally, we also compared the biological functions and pathways of differentially expressed genes (DEGs) between different subgroups by enrichment analysis. In this risk model, we screened out 6 RNA editing sites with significant prognostic value. The constructed nomogram performed well in forecasting patients' prognoses. Furthermore, the level of RNA editing at the prognostic site exhibited a strong correlation with host gene expression. In the high-risk subgroup, we observed multiple biological functions and pathways associated with immune response, cell proliferation, and tumor progression. This study establishes an RNA editing-based risk model that helps forecast patients' prognoses and offers a new understanding of the underlying mechanism of RNA editing in CESC.

Nuclear genomes of dinoflagellates reveal evolutionarily conserved pattern of RNA editing relative to stress response

Dinoflagellates are a group of diverse protists with complex genomes whose gene expression regulation mechanisms remain little known. RNA editing is a post-transcriptional regulatory mechanism of gene expression utilized by diverse species, and has been described primarily in the plastid and mitochondrial genomes of dinoflagellates. Its role in post-transcriptional regulation in the nuclear genomes of dinoflagellates remains largely unexplored. Here, integrating genome and transcriptome data from two dinoflagellate taxa in a comparative analysis, we identified 10,486 and 69,953 putative RNA editing sites in the nuclear genomes of the coral symbiont, Durusdinium trenchii CCMP2556 and the free-living bloom-forming taxon, Prorocentrum cordatum CCMP1329. We recovered all 12 possible types of RNA edits, with more edits representing transitions than transversions. In contrast to other eukaryotes, we found a dominance of A-to-T transversion in non-coding regions, many of which were condition-specific. Overall, the RNA editing sites implicate 7.5% of D. trenchii genes and 13.2% of P. cordatum genes. Some sites (1.5% in D. trenchii and more-substantially 62.3% in P. cordatum) were edited at significantly different frequencies in distinct growth conditions. The distribution of editing types and locations exhibited conserved patterns between the two phylogenetically distant species. Interestingly, A-to-T editing within the untranslated regions appear to be associated with upregulation of the edited genes in response to heat stress. These results lend support to the hypothesis that RNA editing is a key molecular mechanism that underpins regulation of gene expression in dinoflagellates.

De novo assembly of the complete mitochondrial genome of pepino (Solanum muricatum) using PacBio HiFi sequencing: insights into structure, phylogenetic implications, and RNA editing

BackgroundSolanum muricatum is an emerging horticultural fruit crop with rich nutritional and antioxidant properties. Although the chromosome-scale genome of this species has been sequenced, its mitochondrial genome sequence has not been reported to date.ResultsPacBio HiFi sequencing was used to assemble the circular mitogenome of S. muricatum, which was 433,466 bp in length. In total, 38 protein-coding, 19 tRNA, and 3 rRNA genes were annotated. The reticulate mitochondrial conformations with multiple junctions were verified by polymerase chain reaction, and codon usage, sequence repeats, and gene migration from chloroplast to mitochondrial genome were determined. A collinearity analysis of eight Solanum mitogenomes revealed high structural variability. Overall, 585 RNA editing sites in protein coding genes were identified based on RNA-seq data. Among them, mttB was the most frequently edited (52 times), followed by ccmB (46 times). A phylogenetic analysis based on the S. muricatum mitogenome and those of 39 other taxa (including 25 Solanaceae species) revealed the evolutionary and taxonomic status of S. muricatum.ConclusionsWe provide the first report of the assembled and annotated S. muricatum mitogenome. This information will help to lay the groundwork for future research on the evolutionary biology of Solanaceae species. Furthermore, the results will assist the development of molecular breeding strategies for S. muricatum based on the most beneficial agronomic traits of this species.

Elucidating the multichromosomal structure within the Brasenia schreberi mitochondrial genome through assembly and analysis

Brasenia schreberi, a plant species traditionally utilized in Chinese medicine and cuisine, represents an early evolutionary stage among flowering plants (angiosperms). While the plastid genome of this species has been published, its mitochondrial genome (mitogenome) has not been extensively explored, with a notable absence of thorough comparative analyses of its organellar genomes. In our study, we had assembled the entire mitogenome of B. schreberi utilizing the sequencing data derived from both Illumina platform and Oxford Nanopore. The B. schreberi mitogenome mostly exists as six circular DNA molecules, with the largest being 628,257 base pairs (bp) and the smallest 110,220 bp, amounting to 1.49 megabases (Mb). Then we annotated the mitogenome of B. schreberi. The mitogenome encompasses a total of 71 genes: 40 of these are coding proteins genes (PCGs), 28 are genes for transfer RNA (tRNA), and the remaining 3 are genes for ribosomal RNA (rRNA). In the analysis of codon usage, we noted a unique codon preference specific to each amino acid. The most commonly used codons exhibited an average RSCU of 1.36, indicating a noticeable bias in codon selection. In the repeat sequence analysis, a total of 553 simple sequence repeats (SSRs) were identified, 1,822 dispersed repeats (comprising 1,015 forward and 807 palindromic repeats), and 608 long terminal repeats (LTRs). Additionally, in the analysis of homologous sequences between organelle genomes, we detected 38 homologous sequences derived from the plastid genome, each exceeding 500 bp, within the B. schreberi mitochondrial genome. Notably, ten tRNA genes (trnC-GCA, trnM-CAU, trnI-CAU, trnQ-UUG, trnN-GUU, trnT-GGU, trnW-CCA, trnA-UGC, trnI-GAU, and trnV-GAC) appear to have been completely transferred from the chloroplast to the mitogenome. Utilizing the Deepred-mt to predict the RNA editing sites in the mitogenome, we have identified 675 high-quality RNA editing sites in the 40 mitochondrial PCGs. In the final stage of our study, we performed an analysis of colinearity and inferred the phylogenetic relationship of B. schreberi with other angiosperms, utilizing the mitochondrial PCGs as a basis. The results showed that the non-coding regions of the B. schreberi mitogenome are characterized by an abundance of repetitive sequences and exogenous sequences, and B. schreberi is more closely related with Euryale ferox.

The first complete chloroplast genome of Thalictrum fargesii: insights into phylogeny and species identification.

Thalictrum fargesii is a medicinal plant belonging to the genus Thalictrum of the Ranunculaceae family and has been used in herbal medicine in the Himalayan regions of China and India. This species is taxonomically challenging because of its morphological similarities to other species within the genus. Thus, herbal drugs from this species are frequently adulterated, substituted, or mixed with other species, thereby endangering consumer safety. The present study aimed to sequence and assemble the entire chloroplast (cp) genome of T. fargesii using the Illumina HiSeq 2500 platform to better understand the genomic architecture, gene composition, and phylogenetic relationships within the Thalictrum. The cp genome was 155,929 bp long and contained large single-copy (85,395 bp) and small single-copy (17,576 bp) regions that were segregated by a pair of inverted repeat regions (26,479 bp) to form a quadripartite structure. The cp genome contains 133 genes, including 88 protein-coding genes (PCGs), 37 tRNA genes, and 8 rRNA genes. Additionally, this genome contains 64 codons that encode 20 amino acids, the most preferred of which are alanine and leucine. We identified 68 SSRs, 27 long repeats, and 242 high-confidence C-to-U RNA-editing sites in the cp genome. Moreover, we discovered seven divergent hotspot regions in the cp genome of T. fargesii, among which ndhD-psaC and rpl16-rps3 may be useful for developing molecular markers for identifying ethnodrug species and their contaminants. A comparative study with eight other species in the genus revealed that pafI and rps19 had highly variable sites in the cp genome of T.fargesii. Additionally, two special features, (i) the shortest length of the ycf1 gene at the IRA-SSC boundary and (ii) the distance between the rps19 fragment and trnH at the IRA-LSC junction, distinguish the cp genome of T. fargesii from those of other species within the genus. Furthermore, phylogenetic analysis revealed that T. fargesii was closely related to T. tenue and T. petaloidium. Considering all these lines of evidence, our findings offer crucial molecular and evolutionary information that could play a significant role in further species identification, evolution, and phylogenetic studies on T. fargesii.

RNA editing factor SlORRM2 regulates the formation of fruit pointed tips in tomato.

Domestication of tomato (Solanum lycopersicum) has led to large variation in fruit size and morphology. The development of the distal end of the fruit is a critical factor in determining its overall shape. However, the intricate mechanisms underlying distal fruit development require further exploration. This study aimed to investigate the regulatory role of an organelle RNA recognition motif (RRM)-containing protein SlORRM2 in tomato fruit morphology development. Mutant plants lacking SlORRM2 exhibited fruits with pointed tips at the distal end. However, this phenotype could be successfully restored through the implementation of a "functional complementation" strategy. Our findings suggest that the formation of pointed tips in the fruits of the CR-slorrm2 mutants is linked to alterations in the development of the ovary and style. We observed a substantial decrease in the levels of indole-3-acetic acid (IAA) and altered expression of IAA-related response genes in the ovary and style tissues of CR-slorrm2. Moreover, our data demonstrated that SlORRM2 plays a role in regulating mitochondrial RNA editing sites, particularly within genes encoding various respiratory chain subunits. Additionally, the CR-slorrm2 mutants exhibited modified organellar morphology and increased levels of reactive oxygen species. These findings provide valuable insights into the mechanisms underlying the formation of fruit pointed tips in tomato and offer genetic resources for tomato breeding.

Insights into U-to-C RNA editing from the lycophyte Phylloglossum drummondii.

The lycophyte Phylloglossum drummondii is the sole inhabitant of its genus in the Huperzioideae group and one of a small minority of plants which perform uridine to cytidine RNA editing. We assembled the P.drummondii chloroplast and mitochondrial genomes and used RNA sequence data to build a comprehensive profile of organellar RNA editing events. In addition to many C-to-U editing events in both organelles, we found just four U-to-C editing events in the mitochondrial transcripts cob, nad1, nad5 and rpl2. These events are conserved in related lycophytes in the genera Huperzia and Phlegmariurus. De novo transcriptomes for three of these lycophytes were assembled to search for putative U-to-C RNA editing enzymes. Four putative U-to-C editing factors could be matched to the four mitochondrial U-to-C editing sites. Due to the unusually few numbers of U-to-C RNA editing sites, P. drummondii and related lycophytes are useful models for studying this poorly understood mechanism.

New comparative genomic evidence supporting the proteomic diversification role of A-to-I RNA editing in insects.

Adenosine-to-inosine (A-to-I) RNA editing, resembling A-to-G mutation, confers adaptiveness by increasing proteomic diversity in a temporal-spatial manner. This evolutionary theory named "proteomic diversifying hypothesis" has only partially been tested in very few organisms like Drosophila melanogaster, mainly by observing the positive selection on nonsynonymous editing events. To find additional genome-wide evidences supporting this interesting assumption, we retrieved the genomes of four Drosophila species and collected 20 deep-sequenced transcriptomes of different developmental stages and neuron populations of D. melanogaster. We systematically profiled the RNA editomes in these samples and performed meticulous comparative genomic analyses. Further evidences were found to support the diversifying hypothesis. (1) None of the nonsynonymous editing sites in D. melanogaster had ancestral G-alleles, while the silent editing sites had an unignorable fraction of ancestral G-alleles; (2) Only very few nonsynonymous editing sites in D. melanogaster had corresponding G-alleles derived in the genomes of sibling species, and the fraction of such situation was significantly lower than that of silent editing sites; (3) The few nonsynonymous editing with corresponding G-alleles had significantly more variable editing levels (across samples) than other nonsynonymous editing sites in D. melanogaster. The proteomic diversifying nature of RNA editing in Drosophila excludes the restorative role which favors an ancestral G-allele. The few fixed G-alleles in sibling species might facilitate the adaptation to particular environment and the corresponding nonsynonymous editing in D. melanogaster would introduce stronger advantage of flexible proteomic diversification. With multi-Omics data, our study consolidates the nature of evolutionary significance of A-to-I RNA editing sites in model insects.

Assembly and comparative analysis of the complete multichromosomal mitochondrial genome of Cymbidium ensifolium, an orchid of high economic and ornamental value

BackgroundOrchidaceae is one of the largest groups of angiosperms, and most species have high economic value and scientific research value due to their ornamental and medicinal properties. In China, Chinese Cymbidium is a popular ornamental orchid with high economic value and a long history. However, to date, no detailed information on the mitochondrial genome of any species of Chinese Cymbidium has been published.ResultsHere, we present the complete assembly and annotation of the mitochondrial genome of Cymbidium ensifolium (L.) Sw. The mitogenome of C. ensifolium was 560,647 bp in length and consisted of 19 circular subgenomes ranging in size from 21,995 bp to 48,212 bp. The genome encoded 35 protein-coding genes, 36 tRNAs, 3 rRNAs, and 3405 ORFs. Repeat sequence analysis and prediction of RNA editing sites revealed a total of 915 dispersed repeats, 162 simple repeats, 45 tandem repeats, and 530 RNA editing sites. Analysis of codon usage showed a preference for codons ending in A/T. Interorganellar DNA transfer was identified in 13 of the 19 chromosomes, with plastid-derived DNA fragments representing 6.81% of the C. ensifolium mitochondrial genome. The homologous fragments of the mitochondrial genome and nuclear genome were also analysed. Comparative analysis showed that the GC content was conserved, but the size, structure, and gene content of the mitogenomes varied greatly among plants with multichromosomal mitogenome structure. Phylogenetic analysis based on the mitogenomes reflected the evolutionary and taxonomic statuses of C. ensifolium. Interestingly, compared with the mitogenomes of Cymbidium lancifolium Hook. and Cymbidium macrorhizon Lindl., the mitogenome of C. ensifolium lost 8 ribosomal protein-coding genes.ConclusionIn this study, we assembled and annotated the mitogenome of C. ensifolium and compared it with the mitogenomes of other Liliidae and plants with multichromosomal mitogenome structures. Our findings enrich the mitochondrial genome database of orchid plants and reveal the rapid structural evolution of Cymbidium mitochondrial genomes, highlighting the potential for mitochondrial genes to help decipher plant evolutionary history.

Is RNA editing truly absent in the complex thalloid liverworts (Marchantiopsida)? Evidence of extensive RNA editing from Cyathodium cavernarum.

RNA editing is a crucial modification in plants' organellar transcripts that converts cytidine to uridine (C-to-U; and sometimes uridine to cytidine) in RNA molecules. This post-transcriptional process is controlled by the PLS-class protein with a DYW domain, which belongs to the pentatricopeptide repeat (PPR) protein family. RNA editing is widespread in land plants; however, complex thalloid liverworts (Marchantiopsida) are the only group reported to lack both RNA editing and DYW-PPR protein. The liverwort Cyathodium cavernarum (Marchantiopsida, Cyathodiaceae), typically found in cave habitats, was newly found to have 129 C-to-U RNA editing sites in its chloroplast and 172 sites in its mitochondria. The Cyathodium genus, specifically C. cavernarum, has a large number of PPR editing factor genes, including 251 DYW-type PPR proteins. These DYW-type PPR proteins may be responsible for C-to-U RNA editing in C. cavernarum. Cyathodium cavernarum possesses both PPR DYW proteins and RNA editing. Our analysis suggests that the remarkable RNA editing capability of C. cavernarum may have been acquired alongside the emergence of DYW-type PPR editing factors. These findings provide insight into the evolutionary pattern of RNA editing in land plants.

Abstract LB360: Establishment of predictive biomarkers for carcinogenesis in ulcerative colitis using RNA editing analyses

Abstract Background: Ulcerative colitis (UC) is a chronic inflammatory disease continuous from the rectum. We previously reported that chronic inflammation may induce the RNA-editing enzyme ADAR1 and promote carcinogenesis (Journal of Crohn's and Colitis, 2023). However, the relationship between RNA editing and carcinogenesis in UC remains unresolved because we could not identify specific RNA editing sites. Therefore, in this study, we performed a new comprehensive RNA analysis using RNA sequencing to clarify the role of RNA editing in the carcinogenesis process in UC and to construct a prediction model of carcinogenic risk. Methods: We performed immunostaining of ADAR1, an RNA-editing enzyme, in 139 cases of ulcerative colitis. In addition, normal colorectal tissue, UC inflammatory phase, normal rectal mucosa of UC carcinoma cases, and UC carcinoma sites were submitted to RNA-Seq testing in 6 cases each (24 cases total). Bioinformatics analysis was performed to compare RNA editing sites in each group. Results: RNA-seq analysis showed a significant increase in global RNA editing sites in the inflammatory phase compared to normal colon tissue (p = 0.02). This phenomenon was due to increased ADAR1 expression associated with type 1 interferon. We hypothesized that global RNA editing was decreased in normal colonic mucosa during the remission phase but that some sites of prolonged RNA editing might promote carcinogenesis. We attempted to identify carcinogenesis-specific RNA editing sites. We identified 72 cancer-specific RNA editing sites, mainly in long noncoding RNAs such as NEAT1 and in the 3'UTR region of mRNAs such as CTSS. Assuming that the accumulation of these RNA editing leads to carcinogenesis, pathway analysis revealed that pathways promoting cancer invasion and stemness are activated. When ROC curves were generated using these editing sites as markers, the AUC was 0.935, suggesting that analysis of RNA editing in the colonic mucosa of UC patients may be able to predict the risk of carcinogenesis. Conclusions: We identified potential editing sites for carcinogenesis in UC, allowing the identification of carcinogenic mechanisms in UC and the establishment of predictive markers for carcinogenic risk. Citation Format: Toshiaki Takahashi, Kunitoshi Shigeyasu, Yoshitaka Kondo, Hibiki Umeda, Kazuya Moriwake, Masashi Kayano, Kaori Nitta, Fuminori Teraishi, Yuzo Umeda, Hiroshi Tazawa, Toshiyoshi Fujiwara. Establishment of predictive biomarkers for carcinogenesis in ulcerative colitis using RNA editing analyses [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 2 (Late-Breaking, Clinical Trial, and Invited Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(7_Suppl):Abstract nr LB360.

Assembly and comparative analysis of the first complete mitochondrial genome of the invasive water hyacinth, Eichhornia crassipes

Eichhornia crassipes is an aquatic plant in tropical and subtropical regions, renowned for its notorious invasive tendencies. In this study, we assembled the complete mitogenome of E. crassipes into a single circle molecule of 397,361 bp. The mitogenome has 58 unique genes, including 37 protein-coding genes (PCGs), 18 tRNA genes, three rRNA genes, and 47 % GC content. Sixteen (6.93 %) homologous fragments, ranging from 31 bp to 8548 bp, were identified, indicating the transfer of genetic material from chloroplasts to mitochondria. In addition, we detected positive selection in six PCGs (ccmB, ccmC, ccmFC, nad3, nad4 and sdh4), along with the identification of 782 RNA editing sites across 37 mt-PCGs. These findings suggest a potential contribution to the robust adaptation of this invasive plant to the stressful environment. Lastly, we inferred that phylogenetic conflicts of E. crassipes between the plastome and mitogenome may be attributed to the difference in nucleotide substitution rates between the two organelle genomes. In conclusion, our study provided vital genomic resources for further understanding the invasive mechanism of this species and exploring the dynamic evolution of mitogenomes within the monocot clade.

A Comprehensive Analysis of the Effect of A>I(G) RNA-Editing Sites on Genotoxic Drug Response and Progression in Breast Cancer.

Dysregulated A>I(G) RNA editing, which is mainly catalyzed by ADAR1 and is a type of post-transcriptional modification, has been linked to cancer. A low response to therapy in breast cancer (BC) is a significant contributor to mortality. However, it remains unclear if there is an association between A>I(G) RNA-edited sites and sensitivity to genotoxic drugs. To address this issue, we employed a stringent bioinformatics approach to identify differentially RNA-edited sites (DESs) associated with low or high sensitivity (FDR 0.1, log2 fold change 2.5) according to the IC50 of PARP inhibitors, anthracyclines, and alkylating agents using WGS/RNA-seq data in BC cell lines. We then validated these findings in patients with basal subtype BC. These DESs are mainly located in non-coding regions, but a lesser proportion in coding regions showed predicted deleterious consequences. Notably, some of these DESs are previously reported as oncogenic variants, and in genes related to DNA damage repair, drug metabolism, gene regulation, the cell cycle, and immune response. In patients with BC, we uncovered DESs predominantly in immune response genes, and a subset with a significant association (log-rank test p < 0.05) between RNA editing level in LSR, SMPDL3B, HTRA4, and LL22NC03-80A10.6 genes, and progression-free survival. Our findings provide a landscape of RNA-edited sites that may be involved in drug response mechanisms, highlighting the value of A>I(G) RNA editing in clinical outcomes for BC.

Abstract 263: RNA Editing ADAR1 Deaminase Activation Responsive (READAR) platforms for dynamic detection of clonal myeloproliferative neoplasm stem cell fitness

Abstract Introduction Myeloproliferative neoplasms (MPNs) are clonal hematopoietic stem cell (HSC) derived disorders that can progress to acute myeloid leukemia (AML) at variable rates depending on macroenvironmental stressors and microenvironmental inflammatory cytokines. The total symptom score (TSS) and variant allele frequency (VAF) of mutations have proven valuable predictors of therapeutic response. Recently, we discovered inflammatory-cytokine-responsive ADAR1 RNA editing enzyme overexpression during MPN progression. However, dynamic prediction of clonal MPN stem cell fitness and progression to AML remains challenging. Methods We performed longitudinal 150-gene next-generation sequencing (NGS) analyses for 129 MPN patients. We developed single-cell RNA Editing ADAR1 Deaminase Activation Responsive (READAR) platforms to quantify ADAR1 activity with a novel ADAR1 nanoluc-GFP reporter assay, whole transcriptome RNA editome analyses, and MPN stem cell replating assays. To identify novel RNA editing sites, we performed whole genome sequencing on CD34+ cells from 43 MPN patients and whole transcriptome sequencing on FACS-purified HSC and progenitor cells from 32 MPN patients and 24 non-MPN controls. Results Sequential NGS analysis was performed with a median follow-up of 958 days (range 0-4214; interquartile range (IQR) 314-1350). On average, 6 NGS analyses (range 1-19; IQR 3-9) were performed for each patient. In our cohort, 58.9% of patients were JAK2 V617F+, while mutations in CALR and ASXL1 were observed in 11.6% and 14.0% of cases, respectively. In total, 44.2% (N=57) of patients in our cohort received Inrebic, with 56% reaching the maximum dose of 400mg. Notably, most patients on Inrebic were JAK2 V617F+ (N=38, 67.7%). In JAK2 V617F+ patients on Inrebic, 50% (n=19) experienced a reduction in the JAK2 V617F VAF from their peak value with a median drop of 16% (range 4-57). Our analyses using READAR platforms revealed increased RNA editing levels in JAK2 V617F+ patients from whole transcriptome sequencing analyses compared with JAK2 wild-type patients. Transduction of our ADAR1 nanoluc-GFP reporter into CD34+ cells from MPN patients (N=7), enabled us to quantify ADAR1 activity in an ex vivo setting. Moreover, CD34+ cells from MPN patient samples (n=7) treated in stromal co-cultures with ADAR1 inhibitor, Rebecsinib, demonstrated inhibition of self-renewal in replating assays. A significant reduction in ADAR1-GFP reporter activity was observed in MPN CD34+ cells following Rebecsinib treatment of stromal co-cultures. Also, humanized ADAR1-GFP-luciferase reporter AML mouse models showed a significant reduction in ADAR1 activity and leukemia stem cell self-renewal following Rebecsinib treatment. Conclusion These READAR platform results suggest that ADAR1 activity can be utilized to predict clonal MPN stem cell fitness and progression. Citation Format: Inge van der Werf, Larisa Balaian, Jessica Pham, Wenxue Ma, Athena Mohebbi, Emma Klacking, Antonio Ruiz, Karla Mack, John Mascarenhas, Thomas Whisenant, Ludmil Alexandrov, Sheldon Morris, Catriona Jamieson. RNA Editing ADAR1 Deaminase Activation Responsive (READAR) platforms for dynamic detection of clonal myeloproliferative neoplasm stem cell fitness [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 263.