Organellar Transcripts Research Articles

Mitochondrial RNA Helicases: Key Players in the Regulation of Plant Organellar RNA Splicing and Gene Expression.

Mitochondrial genomes of land plants are large and exhibit a complex mode of gene organization and expression, particularly at the post-transcriptional level. The primary organellar transcripts in plants undergo extensive maturation steps, including endo- and/or exo-nucleolytic cleavage, RNA-base modifications (mostly C-to-U deaminations) and both 'cis'- and 'trans'-splicing events. These essential processing steps rely on the activities of a large set of nuclear-encoded factors. RNA helicases serve as key players in RNA metabolism, participating in the regulation of transcription, mRNA processing and translation. They unwind RNA secondary structures and facilitate the formation of ribonucleoprotein complexes crucial for various stages of gene expression. Furthermore, RNA helicases are involved in RNA metabolism by modulating pre-mRNA maturation, transport and degradation processes. These enzymes are, therefore, pivotal in RNA quality-control mechanisms, ensuring the fidelity and efficiency of RNA processing and turnover in plant mitochondria. This review summarizes the significant roles played by helicases in regulating the highly dynamic processes of mitochondrial transcription, RNA processing and translation in plants. We further discuss recent advancements in understanding how dysregulation of mitochondrial RNA helicases affects the splicing of organellar genes, leading to respiratory dysfunctions, and consequently, altered growth, development and physiology of land plants.

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.

Organellar-genome analyses from the lycophyte genus Isoetes L. show one of the highest frequencies of RNA editing in land plants.

RNA editing is a post-transcriptional process that challenges the central dogma of molecular biology by modifying RNA sequences, introducing nucleotide changes at specific sites, and generating functional diversity beyond the genomic code, especially when it concerns organellar transcripts. In plants, this phenomenon is widespread, but its extent varies significantly among species and organellar genomes. Among land plants, the heterosporous lycophytes (i.e., Isoetes and Selaginella) stand out for their exceptionally high numbers of RNA-editing sites, despite their morphological stasis and ancient lineage. In this study, we explore the complete set of organellar protein-coding genes in the aquatic plant group Isoetes, providing a detailed analysis of RNA editing in both the mitochondrial and plastid genomes. Our findings reveal a remarkable abundance of RNA editing, particularly in the mitochondrial genome, with thousands of editing sites identified. Interestingly, the majority of these edits result in non-silent substitutions, suggesting a role in fine-tuning protein structure and function. Furthermore, we observe a consistent trend of increased hydrophobicity in membrane-bound proteins, supporting the notion that RNA editing may confer a selective advantage by preserving gene functionality in Isoetes. The conservation of highly edited RNA sequences over millions of years underscores the evolutionary significance of RNA editing. Additionally, the study sheds light on the dynamic nature of RNA editing, with shared editing sites reflecting common ancestry whereas exclusive edits matching more recent radiation events within the genus. This work advances our understanding of the intricate interplay between RNA editing, adaptation, and evolution in land plants and highlights the unique genomic features of Isoetes, providing a foundation for further investigations into the functional consequences of RNA editing in this enigmatic plant lineage.

Mitochondrial transport of catalytic RNAs and targeting of the organellar transcriptome in human cells.

Mutations in the small genome present in mitochondria often result in severe pathologies. Different genetic strategies have been explored, aiming to rescue such mutations. A number of these strategies were based on the capacity of human mitochondria to import RNAs from the cytosol and designed to repress the replication of the mutated genomes or to provide the organelles with wild-type versions of mutant transcripts. However, the mutant RNAs present in mitochondria turned out to be an obstacle to therapy and little attention has been devoted so far to their elimination. Here, we present the development of a strategy to knockdown mitochondrial RNAs in human cells using the transfer RNA-like structure of Brome mosaic virus or Tobacco mosaic virus as a shuttle to drive trans-cleaving ribozymes into the organelles in human cell lines. We obtained a specific knockdown of the targeted mitochondrial ATP6 mRNA, followed by a deep drop in ATP6 protein and a functional impairment of the oxidative phosphorylation chain. Our strategy provides a powerful approach to eliminate mutant organellar transcripts and to analyse the control and communication of the human organellar genetic system.

Response of the organellar and nuclear (post)transcriptomes of Arabidopsis to drought.

Plants have evolved sophisticated mechanisms to cope with drought, which involve massive changes in nuclear gene expression. However, little is known about the roles of post-transcriptional processing of nuclear or organellar transcripts and how meaningful these changes are. To address these issues, we used RNA-sequencing after ribosomal RNA depletion to monitor (post)transcriptional changes during different times of drought exposure in Arabidopsis Col-0. Concerning the changes detected in the organellar transcriptomes, chloroplast transcript levels were globally reduced, editing efficiency dropped, but splicing was not affected. Mitochondrial transcripts were slightly elevated, while editing and splicing were unchanged. Conversely, alternative splicing (AS) affected nearly 1,500 genes (9% of expressed nuclear genes). Of these, 42% were regulated solely at the level of AS, representing transcripts that would have gone unnoticed in a microarray-based approach. Moreover, we identified 927 isoform switching events. We provide a table of the most interesting candidates, and as proof of principle, increased drought tolerance of the carbonic anhydrase ca1 and ca2 mutants is shown. In addition, altering the relative contributions of the spliced isoforms could increase drought resistance. For example, our data suggest that the accumulation of a nonfunctional FLM (FLOWERING LOCUS M) isoform and not the ratio of FLM-ß and -δ isoforms may be responsible for the phenotype of early flowering under long-day drought conditions. In sum, our data show that AS enhances proteome diversity to counteract drought stress and represent a valuable resource that will facilitate the development of new strategies to improve plant performance under drought.

Emergence of Novel RNA-Editing Sites by Changes in the Binding Affinity of a Conserved PPR Protein.

RNA editing converts cytidines to uridines in plant organellar transcripts. Editing typically restores codons for conserved amino acids. During evolution, specific C-to-U editing sites can be lost from some plant lineages by genomic C-to-T mutations. By contrast, the emergence of novel editing sites is less well documented. Editing sites are recognized by pentatricopeptide repeat (PPR) proteins with high specificity. RNA recognition by PPR proteins is partially predictable, but prediction is often inadequate for PPRs involved in RNA editing. Here we have characterized evolution and recognition of a recently gained editing site. We demonstrate that changes in the RNA recognition motifs that are not explainable with the current PPR code allow an ancient PPR protein, QED1, to uniquely target the ndhB-291 site in Brassicaceae. When expressed in tobacco, the Arabidopsis QED1 edits 33 high-confident off-target sites in chloroplasts and mitochondria causing a spectrum of mutant phenotypes. By manipulating the relative expression levels of QED1 and ndhB-291, we show that the target specificity of the PPR protein depends on the RNA:protein ratio. Finally, our data suggest that the low expression levels of PPR proteins are necessary to ensure the specificity of editing site selection and prevent deleterious off-target editing.

Organellar transcripts dominate the cellular mRNA pool across plants of varying ploidy levels

Mitochondrial and plastid functions depend on coordinated expression of proteins encoded by genomic compartments that have radical differences in copy number of organellar and nuclear genomes. In polyploids, doubling of the nuclear genome may add challenges to maintaining balanced expression of proteins involved in cytonuclear interactions. Here, we use ribo-depleted RNA sequencing (RNA-seq) to analyze transcript abundance for nuclear and organellar genomes in leaf tissue from four different polyploid angiosperms and their close diploid relatives. We find that even though plastid genomes contain <1% of the number of genes in the nuclear genome, they generate the majority (69.9 to 82.3%) of messenger RNA (mRNA) transcripts in the cell. Mitochondrial genes are responsible for a much smaller percentage (1.3 to 3.7%) of the leaf mRNA pool but still produce much higher transcript abundances per gene compared to nuclear genome. Nuclear genes encoding proteins that functionally interact with mitochondrial or plastid gene products exhibit mRNA expression levels that are consistently more than 10-fold lower than their organellar counterparts, indicating an extreme cytonuclear imbalance at the RNA level despite the predominance of equimolar interactions at the protein level. Nevertheless, interacting nuclear and organellar genes show strongly correlated transcript abundances across functional categories, suggesting that the observed mRNA stoichiometric imbalance does not preclude coordination of cytonuclear expression. Finally, we show that nuclear genome doubling does not alter the cytonuclear expression ratios observed in diploid relatives in consistent or systematic ways, indicating that successful polyploid plants are able to compensate for cytonuclear perturbations associated with nuclear genome doubling.

Dek504 Encodes a Mitochondrion-Targeted E+-Type Pentatricopeptide Repeat Protein Essential for RNA Editing and Seed Development in Maize.

In flowering plants, RNA editing is a post-transcriptional process that selectively deaminates cytidines (C) to uridines (U) in organellar transcripts. Pentatricopeptide repeat (PPR) proteins have been identified as site-specific recognition factors for RNA editing. Here, we report the map-based cloning and molecular characterization of the defective kernel mutant dek504 in maize. Loss of Dek504 function leads to delayed embryogenesis and endosperm development, which produce small and collapsed kernels. Dek504 encodes an E+-type PPR protein targeted to the mitochondria, which is required for RNA editing of mitochondrial NADH dehydrogenase 3 at the nad3-317 and nad3-44 sites. Biochemical analysis of mitochondrial protein complexes revealed a significant reduction in the mitochondrial NADH dehydrogenase complex I activity, indicating that the alteration of the amino acid sequence at nad3-44 and nad3-317 through RNA editing is essential for NAD3 function. Moreover, the amino acids are highly conserved in monocots and eudicots, whereas the events of C-to-U editing are not conserved in flowering plants. Thus, our results indicate that Dek504 is essential for RNA editing of nad3, which is critical for NAD3 function, mitochondrial complex I stability, and seed development in maize.

An overview of bioinformatics, genomics, and transcriptomics resources for bryophytes.

Bryophytes are useful models for the study of plant evolution, development, plant-fungal symbiosis, stress responses, and gametogenesis. Additionally, their dominant haploid gametophytic phase makes them great models for functional genomics research, allowing straightforward genome editing and gene knockout via CRISPR or homologous recombination. Until 2016, however, the only bryophyte genome sequence published was that of Physcomitrium patens. Throughout recent years, several other bryophyte genomes and transcriptome datasets became available, enabling better comparative genomics in evolutionary studies. The increase in the number of bryophyte genome and transcriptome resources available has yielded a plethora of annotations, databases, and bioinformatics tools to access the new data, which covers the large diversity of this clade and whose biology comprises features such as association with arbuscular mycorrhiza fungi, sex chromosomes, low gene redundancy, or loss of RNA editing genes for organellar transcripts. Here we provide a guide to resources available for bryophytes with regards to genome and transcriptome databases and bioinformatics tools.

DYW domain structures imply an unusual regulation principle in plant organellar RNA editing catalysis

RNA editosomes selectively deaminate cytidines to uridines in plant organellar transcripts—mostly to restore protein functionality and consequently facilitate mitochondrial and chloroplast function. The RNA editosomal pentatricopeptide repeat proteins serve target RNA recognition, whereas the intensively studied DYW domain elicits catalysis. Here we present structures and functional data of a DYW domain in an inactive ground state and activated. DYW domains harbour a cytidine deaminase fold and a C-terminal DYW motif, with catalytic and structural zinc atoms, respectively. A conserved gating domain within the deaminase fold regulates the active site sterically and mechanistically in a process that we termed gated zinc shutter. Based on the structures, an autoinhibited ground state and its activation are cross-validated by RNA editing assays and differential scanning fluorimetry. We anticipate that, in vivo, the framework of an active plant RNA editosome triggers the release of DYW autoinhibition to ensure a controlled and coordinated cytidine deamination playing a key role in mitochondrial and chloroplast homeostasis.

ZmPPR26, a DYW-type pentatricopeptide repeat protein, is required for C-to-U RNA editing at atpA-1148 in maize chloroplasts.

Pentatricopeptide repeat (PPR) proteins are involved in the C-to-U RNA editing of organellar transcripts. The maize genome contains over 600 PPR proteins and few have been found to function in the C-to-U RNA editing in chloroplasts. Here, we report the function of ZmPPR26 in the C-to-U RNA editing and chloroplast biogenesis in maize. ZmPPR26 encodes a DYW-type PPR protein targeted to chloroplasts. The zmppr26 mutant exhibits albino seedling-lethal phenotype. Loss of function of ZmPPR26 abolishes the editing at atpA-1148 site, and decreases the editing at ndhF-62, rpl20-308, rpl2-2, rpoC2-2774, petB-668, rps8-182, and ndhA-50 sites. Overexpression of ZmPPR26 in zmppr26 restores the editing efficiency and rescues the albino seedling-lethal phenotype. Abolished editing at atpA-1148 causes a Leu to Ser change at AtpA-383 that leads to a reduction in the abundance of chloroplast ATP synthase in zmppr26. The accumulation of photosynthetic complexes are also markedly reduced in zmppr26, providing an explanation for the albino seedling-lethal phenotype. These results indicate that ZmPPR26 is required for the editing at atpA-1148 and is important for editing at the other seven sites in maize chloroplasts. The editing at atpA-1148 is critical for AtpA function, assembly of ATP synthase complex, and chloroplast biogenesis in maize.

Differentially Expressed Genes Shared by Two Distinct Cytoplasmic Male Sterility (CMS) Types of Silene vulgaris Suggest the Importance of Oxidative Stress in Pollen Abortion.

Cytoplasmic male sterility (CMS), encoded by the interacting mitochondrial and nuclear genes, causes pollen abortion or non-viability. CMS is widely used in agriculture and extensively studied in crops. Much less is known about CMS in wild species. We performed a comparative transcriptomic analysis of male sterile and fertile individuals of Silene vulgaris, a model plant for the study of gynodioecy, to reveal the genes responsible for pollen abortion in this species. We used RNA-seq datasets previously employed for the analysis of mitochondrial and plastid transcriptomes of female and hermaphrodite flower buds, making it possible to compare the transcriptomes derived from three genomes in the same RNA specimen. We assembled de novo transcriptomes for two haplotypes of S. vulgaris and identified differentially expressed genes between the females and hermaphrodites, associated with stress response or pollen development. The gene for alternative oxidase was downregulated in females. The genetic pathways controlling CMS in S. vulgaris are similar to those in crops. The high number of the differentially expressed nuclear genes contrasts with the uniformity of organellar transcriptomes across genders, which suggests these pathways are evolutionarily conserved and that selective mechanisms may shield organellar transcription against changes in the cytoplasmic transcriptome.

Plant organellar RNA editing: what 30years of research has revealed.

The central dogma in biology defines the flow of genetic information from DNA to RNA to protein. Accordingly, RNA molecules generally accurately follow the sequences of the genes from which they are transcribed. This rule is transgressed by RNA editing, which creates RNA products that differ from their DNA templates. Analyses of the RNA landscapes of terrestrial plants have indicated that RNA editing (in the form of C-U base transitions) is highly prevalent within organelles (that is, mitochondria and chloroplasts). Numerous C→U conversions (and in some plants also U→C) alter the coding sequences of many of the organellar transcripts and can also produce translatable mRNAs by creating AUG start sites or eliminating premature stop codons, or affect the RNA structure, influence splicing and alter the stability of RNAs. RNA-binding proteins are at the heart of post-transcriptional RNA expression. The C-to-U RNA editing process in plant mitochondria involves numerous nuclear-encoded factors, many of which have been identified as pentatricopeptide repeat (PPR) proteins that target editing sites in a sequence-specific manner. In this review we report on major discoveries on RNA editing in plant organelles, since it was first documented 30years ago.

SMALL KERNEL4is required for mitochondrial cox1 transcript editing and seed development in maize.

In land plants, cytidine-to-uridine (C-to-U) editing of organellar transcripts is an important post-transcriptional process, which is considered to remediate DNA genetic mutations to restore the coding of functional proteins. Pentatricopeptide repeat (PPR) proteins have key roles in C-to-U editing. Owing to its large number, however, the biological functions of many PPR proteins remain to be identified. Through characterizing a small kernel4 (smk4) mutant, here we report the function of Smk4 and its role in maize growth and development. Null mutation of Smk4 slows plant growth and development, causing small plants, delayed flowering time, and small kernels. Cloning revealed that Smk4 encodes a new E-subclass PPR protein, and localization indicated that SMK4 is exclusively localized in mitochondria. Loss of Smk4 function abolishes C-to-U editing at position 1489 of the cytochrome c oxidase1 (cox1) transcript, causing an amino acid change from serine to proline at 497 in Cox1. Cox1 is a core component of mitochondrial complex IV. Indeed, complex IV activity is reduced in the smk4, along with drastically elevated expression of alternative oxidases (AOX). These results indicate that SMK4 functions in the C-to-U editing of cox1-1489, and this editing is crucial for mitochondrial complex IV activity, plant growth, and kernel development in maize.

Genome-wide Identification of Chile Pepper Pentatricopeptide Repeat Domains Provides Insight into Fertility Restoration

Pentatricopeptide repeat (PPR) proteins constitute one of the largest protein families in plants. They are typically targeted to the mitochondria or chloroplast and bind to one or more organellar transcripts, influencing expression. Genes responsible for inhibiting the mitochondrial genes that cause cytoplasmic male sterility (CMS), restorer-of-fertility (Rf) genes, often encode PPRs. Using in silico techniques, 552 PPR domains were identified throughout the chile pepper (Capsicum annuum) genome. The domains were mapped across 12 chromosomes and were found to be largely distally or proximally located. About 28% of the chile pepper PPR domains identified in this study have high structural similarity to previously reported PPRs in arabidopsis (Arabidopsis thaliana). In addition, 11 candidate Rf genes clustered on chromosome 6, and 1 on chromosome 1 were identified that were characterized in 16 A- (S rfrf), B- (N rfrf), and C-line (N RfRf) backgrounds. These findings support a multigene model for fertility restoration and broaden our understanding of the restoration of fertility. This may be an explanation for the lack of widely applicable molecular markers for this important trait. With this new information, specific Rf markers may be developed and will facilitate the implementation of hybrid breeding programs in chile pepper. In addition, this work provides a basis for future research in PPRs, an increasingly important gene family.

Towards a comprehensive picture of C-to-U RNA editing sites in angiosperm mitochondria.

Our understanding of the dynamic and evolution of RNA editing in angiosperms is in part limited by the few editing sites identified to date. This study identified 10,217 editing sites from 17 diverse angiosperms. Our analyses confirmed the universality of certain features of RNA editing, and offer new evidence behind the loss of editing sites in angiosperms. RNA editing is a post-transcriptional process that substitutes cytidines (C) for uridines (U) in organellar transcripts of angiosperms. These substitutions mostly take place in mitochondrial messenger RNAs at specific positions called editing sites. By means of publicly available RNA-seq data, this study identified 10,217 editing sites in mitochondrial protein-coding genes of 17 diverse angiosperms. Even though other types of mismatches were also identified, we did not find evidence of non-canonical editing processes. The results showed an uneven distribution of editing sites among species, genes, and codon positions. The analyses revealed that editing sites were conserved across angiosperms but there were some species-specific sites. Non-synonymous editing sites were particularly highly conserved (~ 80%) across the plant species and were efficiently edited (80% editing extent). In contrast, editing sites at third codon positions were poorly conserved (~ 30%) and only partially edited (~ 40% editing extent). We found that the loss of editing sites along angiosperm evolution is mainly occurring by replacing editing sites with thymidines, instead of a degradation of the editing recognition motif around editing sites. Consecutive and highly conserved editing sites had been replaced by thymidines as result of retroprocessing, by which edited transcripts are reverse transcribed to cDNA and then integrated into the genome by homologous recombination. This phenomenon was more pronounced in eudicots, and in the gene cox1. These results suggest that retroprocessing is a widespread driving force underlying the loss of editing sites in angiosperm mitochondria.

High-throughput sequencing of the chloroplast and mitochondrion of Chlamydomonas reinhardtii to generate improved de novo assemblies, analyze expression patterns and transcript speciation, and evaluate diversity among laboratory strains and wild isolates.

Chlamydomonas reinhardtii is a unicellular chlorophyte alga that is widely studied as a reference organism for understanding photosynthesis, sensory and motile cilia, and for development of an algal-based platform for producing biofuels and bio-products. Its highly repetitive, ~205-kbp circular chloroplast genome and ~15.8-kbp linear mitochondrial genome were sequenced prior to the advent of high-throughput sequencing technologies. Here, high coverage shotgun sequencing was used to assemble both organellar genomes de novo. These new genomes correct dozens of errors in the prior genome sequences and annotations. Genome sequencing coverage indicates that each cell contains on average 83 copies of the chloroplast genome and 130 copies of the mitochondrial genome. Using protocols and analyses optimized for organellar transcripts, RNA-Seq was used to quantify their relative abundances across 12 different growth conditions. Forty-six percent of total cellular mRNA is attributable to high expression from a few dozen chloroplast genes. RNA-Seq data were used to guide gene annotation, to demonstrate polycistronic gene expression, and to quantify splicing of psaA and psbA introns. In contrast to a conclusion from a recent study, we found that chloroplast transcripts are not edited. Unexpectedly, cytosine-rich polynucleotide tails were observed at the 3'-end of all mitochondrial transcripts. A comparative genomics analysis of eight laboratory strains and 11 wild isolates of C.reinhardtii identified 2658 variants in the organellar genomes, which is 1/10th as much genetic diversity as is found in the nucleus.

Pervasive, Genome-Wide Transcription in the Organelle Genomes of Diverse Plastid-Bearing Protists.

Organelle genomes are among the most sequenced kinds of chromosome. This is largely because they are small and widely used in molecular studies, but also because next-generation sequencing technologies made sequencing easier, faster, and cheaper. However, studies of organelle RNA have not kept pace with those of DNA, despite huge amounts of freely available eukaryotic RNA-sequencing (RNA-seq) data. Little is known about organelle transcription in nonmodel species, and most of the available eukaryotic RNA-seq data have not been mined for organelle transcripts. Here, we use publicly available RNA-seq experiments to investigate organelle transcription in 30 diverse plastid-bearing protists with varying organelle genomic architectures. Mapping RNA-seq data to organelle genomes revealed pervasive, genome-wide transcription, regardless of the taxonomic grouping, gene organization, or noncoding content. For every species analyzed, transcripts covered ≥85% of the mitochondrial and/or plastid genomes (all of which were ≤105 kb), indicating that most of the organelle DNA—coding and noncoding—is transcriptionally active. These results follow earlier studies of model species showing that organellar transcription is coupled and ubiquitous across the genome, requiring significant downstream processing of polycistronic transcripts. Our findings suggest that noncoding organelle DNA can be transcriptionally active, raising questions about the underlying function of these transcripts and underscoring the utility of publicly available RNA-seq data for recovering complete genome sequences. If pervasive transcription is also found in bigger organelle genomes (>105 kb) and across a broader range of eukaryotes, this could indicate that noncoding organelle RNAs are regulating fundamental processes within eukaryotic cells.

Protein Complexes Implicated in RNA Editing in Plant Organelles

Even after 30 years of active research, RNA editing in plant organelles remains an enigma. Why do most plants post-transcriptionally alter hundreds or thousands of nucleotides in organellar transcripts? How do they achieve such an effective and targeted process? What is known is that in almost all plants, RNA editing is essential for correct synthesis of the photosynthetic and respiratory machinery in mitochondria and chloroplasts. The sites to be edited are defined by a large class of sequence-specific RNA binding proteins: pentatricopeptide repeat (PPR) proteins (Barkan and Small, 2014).

Evaluation of reference genes for reverse transcription quantitative real-time PCR (RT-qPCR) studies in Silene vulgaris considering the method of cDNA preparation.

Accurate gene expression measurements are essential in studies of both crop and wild plants. Reverse transcription quantitative real-time PCR (RT-qPCR) has become a preferred tool for gene expression estimation. A selection of suitable reference genes for the normalization of transcript levels is an essential prerequisite of accurate RT-qPCR results. We evaluated the expression stability of eight candidate reference genes across roots, leaves, flower buds and pollen of Silene vulgaris (bladder campion), a model plant for the study of gynodioecy. As random priming of cDNA is recommended for the study of organellar transcripts and poly(A) selection is indicated for nuclear transcripts, we estimated gene expression with both random-primed and oligo(dT)-primed cDNA. Accordingly, we determined reference genes that perform well with oligo(dT)- and random-primed cDNA, making it possible to estimate levels of nucleus-derived transcripts in the same cDNA samples as used for organellar transcripts, a key benefit in studies of cyto-nuclear interactions. Gene expression variance was estimated by RefFinder, which integrates four different analytical tools. The SvACT and SvGAPDH genes were the most stable candidates across various organs of S. vulgaris, regardless of whether pollen was included or not.