Validation of reference genes for quantitative real-time PCR during leaf and flower development in Petunia hybrida

Evaluation of reference genes and expression patterns of CONSTANS-LIKE genes in Tetrastigma hemsleyanum under different photoperiods

CONSTANS-LIKE (COL ) genes are a key signalling molecule that regulates plant growth and development during the photoperiod. Our preliminary experiments showed that the photoperiod greatly influence the formation of Tetrastigma hemsleyanum root tubers. In this study, we examined the oscillation patterns and expression characteristics of COL genes in leaves of T. hemsleyanum under different photoperiod conditions. Six genes were selected as candidate reference genes for further analyses: (1) 18S ribosomal RNA (18S rRNA ); (2) α-tubulin (TUBA ); (3) 30S ribosomal RNA (30S rRNA ); (4) TATA binding protein (TBP ); (5) elongation factor 1α (EF-1α ); and (6) RNA polymerase II (RPII ). The geNorm, NormFinder, and BestKeeper software programs were used to evaluate expression stability. Two ThCOL genes were screened in the T. hemsleyanum transcriptome library, and their expression patterns under different photoperiod conditions were analysed using quantitative reverse transcription PCR. The genes EF-1α , TUBA , and 18S rRNA were used to analyse the expression profiles of CONSTANS genes (ThCOL4 and ThCOL5 ) under different photoperiods. The expression peaks of ThCOL4 and ThCOL5 appeared at different times, demonstrating that their oscillation patterns were influenced by the photoperiod. We speculate that these two ThCOL genes may be involved in different biological processes.

Dorsoventrally asymmetric expression of miR319/TCP generates dorsal-specific venation patterning in petunia corolla tube

Vein-associated pigmentation (venation) is a type of floral coloration adopted by plants to attract pollinators. Several petunia (Petunia hybrida) lines generate dorsoventrally asymmetric venation patterning of the corolla tube, in which venation is only present in the dorsal tube. The molecular mechanism underlying this trait is unknown. Here, we demonstrate that miR319 is preferentially expressed in the dorsal corolla tube, leading to dorsoventrally asymmetric expression of its target genes. Transgenic lines overexpressing phy-miR319a generated uniform venation patterning of the corolla tube. Knockout of TCP genes targeted by miR319 promoted venation patterning in the ventral and dorsal tube, while overexpression of the miR319 target gene, PhTCP6, completely inhibited corolla tube venation patterning. In addition, miR319-targeted TCPs negatively regulated venation patterning, probably by repressing the regulator of venation patterning, AN4. Together, our data demonstrate that asymmetric expression of miR319 promotes venation patterning in the petunia dorsal tube alone by repressing the expression of its target TCP genes, which negatively regulate corolla tube venation patterning. These findings provide novel insights into how the dorsoventrally asymmetric distribution of venation patterning is established in zygomorphic flowers.

CRISPR/Cas9-mediated mutagenesis of FT/TFL1 in petunia improves plant architecture and early flowering

Petunias are renowned ornamental species widely cultivated as pot plants for their aesthetic appeal both indoors and outdoors. The preference for pot plants depends on their compact growth habit and abundant flowering. While genome editing has gained significant popularity in many crop plants in addressing growth and development and abiotic and biotic stress factors, relatively less emphasis has been placed on its application in ornamental plant species. Genome editing in ornamental plants opens up possibilities for enhancing their aesthetic qualities, offering innovative opportunities for manipulating plant architecture and visual appeal through precise genetic modifications. In this study, we aimed to optimize the procedure for an efficient genome editing system in petunia plants using the highly efficient multiplexed CRISPR/Cas9 system. Specifically, we targeted a total of six genes in Petunia which are associated with plant architecture traits, two paralogous of FLOWERING LOCUS T (PhFT) and four TERMINAL FLOWER-LIKE1 (PhTFL1) paralogous genes separately in two constructs. We successfully induced homogeneous and heterogeneous indels in the targeted genes through precise genome editing, resulting in significant phenotypic alterations in petunia. Notably, the plants harboring edited PhTFL1 and PhFT exhibited a conspicuously early flowering time in comparison to the wild-type counterparts. Furthermore, mutants with alterations in the PhTFL1 demonstrated shorter internodes than wild-type, likely by downregulating the gibberellic acid pathway genes PhGAI, creating a more compact and aesthetically appealing phenotype. This study represents the first successful endeavor to produce compact petunia plants with increased flower abundance through genome editing. Our approach holds immense promise to improve economically important potting plants like petunia and serve as a potential foundation for further improvements in similar ornamental plant species.

Exploring laccase: a sustainable enzymatic solution for the paper recycling domain

The global advocacy of resource conservation and waste management emphasizes the significance of sustainable practices, particularly in sectors such as paper manufacturing and recycling. Currently, conventional chemical methods are predominant for paper production, necessitating the use of substantial amount of toxic chemicals. This chemical-intensive approach compromises the recycled fiber quality, generates hazardous effluent causing serious ecological threats which triggers regulatory complexities for the mills. To address these challenges modern research suggests adopting sustainable eco-friendly practices such as employing enzymes. This review aims to explore the applicability of 'laccase' enzyme for paper recycling, investigating its properties and contribution to improved recycling practices. By delving into the potential application of laccase integration into the papermaking process, this article sheds light on the limitations inherent in traditional methods surmounted within both research and translational landscapes. Culture and process optimization studies, supporting the technological improvements and the future prospects have been documented.

Reliable reference gene selection for quantitative real-time PCR (qRT-PCR) in floral developmental phases of dioecious species Coccinia grandis

The flowering process is intricate and regulated by a combination of external and internal factors. Delving into gene expression research has the potential to enhance our comprehension of the molecular foundations underlying floral development. Because of its accuracy, specificity, reproducibility, and efficiency, qRT-PCR is now a biological research tool for studying expression pattern of desired genes. The gene expression investigations using qRT-PCR required a reference gene with relatively uniform expression levels in multiple biological samples, including different developmental stages, tissues, and experimental conditions. In this study, experimental sets offloral and floral organ development in the male and female plants of C. grandis, a dioecious Cucurbitaceae species, qRT-PCR profiling was performed using six reference genes as internal control with B-class floral identity gene, PISTILLATA (PI). To analyse the data, algorithms such as geNorm, NormFinder, RefFinder, and BestKeeper were used to pick out the best internal controls from a group of candidates. The optimal reference gene for qRT-PCR studies with floral samples has been recommended as β-actin combined with β-tubulin. This is the first report on the validation of candidate reference genes across flower developmental stages in the dioecious species C. grandis, which will provide basic data for research on the molecular mechanism underlying flower development in this species and lay the groundwork for similar studies in other related species.

Abscisic acid-induced transcription factor PsMYB306 negatively regulates tree peony bud dormancy release

Bud dormancy is a crucial strategy for perennial plants to withstand adverse winter conditions. However, the regulatory mechanism of bud dormancy in tree peony (Paeonia suffruticosa) remains largely unknown. Here, we observed dramatically reduced and increased accumulation of abscisic acid (ABA) and bioactive gibberellins (GAs) GA1 and GA3, respectively, during bud endodormancy release of tree peony under prolonged chilling treatment. An Illumina RNA sequencing study was performed to identify potential genes involved in the bud endodormancy regulation in tree peony. Correlation matrix, principal component, and interaction network analyses identified a downregulated MYB transcription factor gene, PsMYB306, the expression of which positively correlated with 9-CIS-EPOXYCAROTENOID DIOXYGENASE 3 (PsNCED3) expression. Protein modeling analysis revealed 4 residues within the R2R3 domain of PsMYB306 to possess DNA binding capability. Transcription of PsMYB306 was increased by ABA treatment. Overexpression of PsMYB306 in petunia (Petunia hybrida) inhibited seed germination and plant growth, concomitant with elevated ABA and decreased GA contents. Silencing of PsMYB306 accelerated cold-triggered tree peony bud burst and influenced the production of ABA and GAs and the expression of their biosynthetic genes. ABA application reduced bud dormancy release and transcription of ENT-KAURENOIC ACID OXIDASE 1 (PsKAO1), GA20-OXIDASE 1 (PsGA20ox1), and GA3-OXIDASE 1 (PsGA3ox1) associated with GA biosynthesis in PsMYB306-silenced buds. In vivo and in vitro binding assays confirmed that PsMYB306 specifically transactivated the promoter of PsNCED3. Silencing of PsNCED3 also promoted bud break and growth. Altogether, our findings suggest that PsMYB306 negatively modulates cold-induced bud endodormancy release by regulating ABA production.

Characterizing reference genes for high-fidelity gene expression analysis under different abiotic stresses and elicitor treatments in fenugreek leaves

BACKGROUND

Quantifying gene expression is a critical aspect of applied genomics research across all organisms, and real-time PCR has emerged as a powerful tool for this purpose. However, selecting appropriate internal control genes for data normalization presents specific challenges. This study aimed to identify suitable reference genes for gene expression analysis under various conditions, encompassing salinity, low and high-temperature stresses, and different elicitor treatments. These treatments included titanium dioxide, cold plasma, 24-epibrassinolide, and melatonin, resulting in a total of 13 unique treatments and 148 treatment combinations applied to fenugreek plants.

RESULTS

As per the analysis performed with the BestKeeper tool, EEF-1α, and GAPDH were recognized as the most stable reference genes under the majority of conditions. Furthermore, the GeNorm and NormFinder tools identified β-tubulin and EEF-1α as the most stable reference genes. The findings of this research demonstrated that, although the stability of three reference genes expression was acceptable in almost all evaluated treatments, fluctuations in their expression were observed under the treatments of cold stress with TiO2 NPs application, cold plasma application with salinity stress, and cold plasma application with high-temperature stress compared to others. Simultaneously, the GeNorm analysis results demonstrated that in the mentioned treatments, relying on only one reference gene is inadequate. To corroborate the results, we examined the expression profile of the SSR gene, a pivotal gene in diosgenin biosynthesis, under all investigated treatments and treatment combinations. The outcomes suggested that employing stable reference genes yielded highly consistent results.

CONCLUSIONS

The varying expression patterns of the target genes emphasize the crucial need for precise optimization of experimental conditions and selecting stable reference genes to achieve accurate results in gene expression studies utilizing real-time PCR. These findings offer valuable insights into the selection of appropriate reference genes for gene expression analysis under diverse conditions using real-time PCR.

Comprehensive Analyses of Four PhNF-YC Genes from Petunia hybrida and Impacts on Flowering Time

Nuclear Factor Y (NF-Y) is a class of heterotrimeric transcription factors composed of three subunits: NF-A, NF-YB, and NF-YC. NF-YC family members play crucial roles in various developmental processes, particularly in the regulation of flowering time. However, their functions in petunia remain poorly understood. In this study, we isolated four PhNF-YC genes from petunia and confirmed their subcellular localization in both the nucleus and cytoplasm. We analyzed the transcript abundance of all four PhNF-YC genes and found that PhNF-YC2 and PhNF-YC4 were highly expressed in apical buds and leaves, with their transcript levels decreasing before flower bud differentiation. Silencing PhNF-YC2 using VIGS resulted in a delayed flowering time and reduced chlorophyll content, while PhNF-YC4-silenced plants only exhibited a delayed flowering time. Furthermore, we detected the transcript abundance of flowering-related genes involved in different signaling pathways and found that PhCO, PhGI, PhFBP21, PhGA20ox4, and PhSPL9b were regulated by both PhNF-YC2 and PhNF-YC4. Additionally, the transcript abundance of PhSPL2, PhSPL3, and PhSPL4 increased only in PhNF-YC2-silenced plants. Overall, these results provide evidence that PhNF-YC2 and PhNF-YC4 negatively regulate flowering time in petunia by modulating a series of flowering-related genes.

Transcriptome-Based Identification of the Optimal Reference Genes for Quantitative Real-Time Polymerase Chain Reaction Analyses of Lingonberry Fruits throughout the Growth Cycle

(1) Background: Vaccinium vitis-idaea is a nutritionally and economically valuable natural wild plant species that produces berries useful for treating various diseases. There is growing interest in lingonberry, but there is limited information regarding lingonberry reference genes suitable for gene expression analyses of different tissues under various abiotic stress conditions. The objective of this study was to identify stable reference genes suitable for different lingonberry tissues in response to abiotic stress. (2) Methods: The delta Ct method and the GeNorm v3.5 and NormFinder v20 programs were used to comprehensively analyze gene expression stability. (3) Results: Actin Unigene23839 was the best reference gene for analyzing different cultivars, whereas Actin CL5740.Contig2 was the most suitable reference gene for analyzing different tissues and alkali stress. In contrast, 18S rRNA CL5051.Contig1 was the most stable reference gene under drought conditions. (4) Conclusions: These suitable reference genes may be used in future qRT-PCR analyses of different lingonberry tissues and the effects of abiotic stresses. Furthermore, the study data may be useful for functional genomics studies and the molecular breeding of lingonberry. In summary, internal reference genes or internal reference gene combinations should be carefully selected according to the experimental conditions to ensure that the generated gene expression data are accurate.

New insights into the genetic manipulation of the R2R3-MYB and CHI gene families on anthocyanin pigmentation in Petunia hybrida

Several R2R3-MYB genes control anthocyanin pigmentation in petunia, and ANTHOCYANIN-2 (AN2) is treated as the main player in petal limbs. However, the actual roles of R2R3-MYBs in the coloration of different floral tissues in the so called "darkly-veined" petunias are still not clear. The genetic background and expression of AN2 paralogs from various petunias with different color patterns were identified. All "darkly-veined" genotypes have the identical mutation in the AN2 gene, but express a different functional paralog - ANTHOCYANIN-4 (AN4) - abundantly in flowers. Constitutive overexpression of PhAN4 in this petunia resulted not only in a fully colored flower but also in a clearly visible pigmentation in the green tissue and roots, which can be rapidly increased by stress conditions. Suppression of AN4 gene resulted in discolored petals and whitish anthers. Interestingly, when a similar white flower phenotype was achieved by knockout of an essential structural gene of anthocyanin biosynthesis - CHALCONE ISOMERASE-A (CHI-A) - the plant responded directly by upregulating of another paralogs - DEEP PURPLE (DPL) and PURPLE HAZE (PHZ). Moreover, we also found that CHI-B can partially substitute for CHI-A in anthers, but not in vegetative tissues. Further, no significant effects on the longevity of white or enhanced colored flowers were observed compared with the wild type. We concluded that endogenous up-regulation of AN4 leads to the restoration of petal color in the "darkly-veined" phenotypes as a result of the breeding process under human selection, and CHI-B is a backup for CHI-A acitvity in some floral tissues.

Selection and validation of reliable reference genes for quantitative real-time PCR in Barnyard millet (Echinochloa spp.) under varied abiotic stress conditions

Quantitative real-time polymerase chain reaction (RT-qPCR) using a stable reference gene is widely used for gene expression research. Barnyard millet (Echinochloa spp.) is an ancient crop in Asia and Africa that is widely cultivated for food and fodder. It thrives well under drought, salinity, cold, and heat environmental conditions, besides adapting to any soil type. To date, there are no gene expression studies performed to identify the potential candidate gene responsible for stress response in barnyard millet, due to lack of reference gene. Here, 10 candidate reference genes, Actin (ACT), α-tubulin (α-TUB), β-tubulin (β-TUB), RNA pol II (RP II), elongation factor-1 alpha (EF-1α), adenine phosphoribosyltransferase (APRT), TATA-binding protein-like factor (TLF), ubiquitin-conjugating enzyme 2 (UBC2), ubiquitin-conjugating enzyme E2L5 (UBC5) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH), were selected from mRNA sequences of E. crus-galli and E. colona var frumentacea. Five statistical algorithms (geNorm, NormFinder, BestKeeper, ΔCt, and RefFinder) were applied to determine the expression stabilities of these genes in barnyard millet grown under four different abiotic stress (drought, salinity, cold and heat) exposed at different time points. The UBC5 and ɑ-TUB in drought, GAPDH in salinity, GAPDH and APRT in cold, and EF-1α and RP II in heat were the most stable reference genes, whereas ß-TUB was the least stable irrespective of stress conditions applied. Further Vn/Vn + 1 analysis revealed two reference genes were sufficient to normalize gene expression across all sample sets. The suitability of identified reference genes was validated with Cu-ZnSOD (SOD1) in the plants exposed to different abiotic stress conditions. The results revealed that the relative quantification of the SOD1 gene varied according to reference genes and the number of reference genes used, thus highlighting the importance of the choice of a reference gene in such experiments. This study provides a foundational framework for standardizing RT-qPCR analyses, enabling accurate gene expression profiling in barnyard millet.

Hormonal Signaling during dPCD: Cytokinin as the Determinant of RNase-Based Self-Incompatibility in Solanaceae

Research into molecular mechanisms of self-incompatibility (SI) in plants can be observed in representatives of various families, including Solanaceae. Earlier studies of the mechanisms of S-RNase-based SI in petunia (Petunia hybrida E. Vilm.) demonstrate that programmed cell death (PCD) is an SI factor. These studies suggest that the phytohormon cytokinin (CK) is putative activator of caspase-like proteases (CLPs). In this work, data confirming this hypothesis were obtained in two model objects-petunia and tomato (six Solanaceae representatives). The exogenous zeatin treatment of tomato and petunia stigmas before a compatible pollination activates CLPs in the pollen tubes in vivo, as shown via the intravital imaging of CLP activities. CK at any concentration slows down the germination and growth of petunia and tomato male gametophytes both in vitro and in vivo; shifts the pH of the cytoplasm (PHc) to the acid region, thereby creating the optimal conditions for CLP to function and inhibiting the F-actin formation and/or destructing the cytoskeleton in pollen tubes to point foci during SI-induced PCD; and accumulates in style tissues during SI response. The activity of the ISOPENTENYLTRANSFERASE 5 (IPT5) gene at this moment exceeds its activity in a cross-compatible pollination, and the levels of expression of the CKX1 and CKX2 genes (CK OXIDASE/DEHYDROGENASE) are significantly lower in self-incompatible pollination. All this suggests that CK plays a decisive role in the mechanism underlying SI-induced PCD.

Functional characterization of three TERMINAL FLOWER 1-like genes from Platanus acerifolia

KEY MESSAGE

TFL1-like genes of the basal eudicot Platanus acerifolia have conserved roles in maintaining vegetative growth and inhibiting flowering, but may act through distinct regulatory mechanism. Three TERMINAL FLOWER 1 (TFL1)-like genes were isolated and characterized from London plane tree (Platanus acerifolia). All genes have conserved genomic organization and characteristic of the phosphatidylethanolamine-binding protein (PEBP) family. Sequence alignment and phylogenetic analysis indicated that two genes belong to the TFL1 clade, designated as PlacTFL1a and PlacTFL1b, while another one was grouped in the BFT clade, named as PlacBFT. qRT-PCR analysis showed that all three genes primarily expressed in vegetative phase, but the expression of PlacTFL1a was much higher and wider than that of PlacTFL1b, with the latter only detected at relatively low expression levels in apical and lateral buds in April. PlacBFT was mainly expressed in young stems of adult trees followed by juvenile tissues. Ectopic expression of any TFL1-like gene in Arabidopsis showed phenotypes of delayed or repressed flowering. Furthermore, overexpression of PlacTFL1a gene in petunia also resulted in extremely delayed flowering. In non-flowering 35:PlacTFL1a transgenic petunia plants, the FT-like gene (PhFT) gene was significantly upregulated and AP1 homologues PFG, FBP26 and FBP29 were significantly down-regulated in leaves. Yeast two-hybrid analysis indicated that only weak interactions were detected between PlacTFL1a and PlacFDL, and PlacTFL1a showed no interaction with PhFDL1/2. These results indicated that the TFL1-like genes of Platanus have conserved roles in repressing flowering, but probably via a distinct regulatory mechanism.

Drought and Oxidative Stress in Flax (Linum usitatissimum L.) Entails Harnessing Non-Canonical Reference Gene for Precise Quantification of qRT-PCR Gene Expression

Flax (Linum usitatissimum L.) is a self-pollinating, annual, diploid crop grown for multi-utility purposes for its quality oil, shining bast fiber, and industrial solvent. Being a cool (Rabi) season crop, it is affected by unprecedented climatic changes such as high temperature, drought, and associated oxidative stress that, globally, impede its growth, production, and productivity. To precisely assess the imperative changes that are inflicted by drought and associated oxidative stress, gene expression profiling of predominant drought-responsive genes (AREB, DREB/CBF, and ARR) was carried out by qRT-PCR. Nevertheless, for normalization/quantification of data obtained from qRT-PCR results, a stable reference gene is mandatory. Here, we evaluated a panel of four reference genes (Actin, EF1a, ETIF5A, and UBQ) and assessed their suitability as stable reference genes for the normalization of gene expression data obtained during drought-induced oxidative stress in flax. Taking together, from the canonical expression of the proposed reference genes in three different genotypes, we report that EF1a as a stand-alone and EF1a and ETIF5A in tandem are suitable reference genes to be used for the real-time visualization of cellular impact of drought and oxidative stress on flax.

Physiological changes besides the enhancement of pigmentation in Petunia hybrida caused by overexpression of PhAN2, an R2R3-MYB transcription factor

Ectopic expression of PhAN2 in vegetative tissue can improve regeneration and adventitious rooting but inhibit axillary bud outgrowth of petunia, while overexpression specifically in flowers could shorten longevity. Anthocyanin 2 has been only treated as a critical positive regulation factor of anthocyanin biosynthesis in petunia flowers. To determine if this gene had other functions in plant growth, we overexpressed this gene in an an2 mutant petunia cultivar driven by promoters with different strengths or tissue specificity. Various physiological processes of transformants in different growth stages and environments were analyzed. Besides the expected pigmentation improvement in different tissues, the results also showed that ectopic expression of AN2 could improve the regeneration skill but inhibit the axillary bud germination of in vitro plants. Moreover, the rooting ability of shoot tips of transformants was significantly improved, while some transgenic lines' flower longevity was shortened. Gene expression analysis showed that the transcripts level of AN2, partner genes anthocyanin 1 (AN1), anthocyanin 11 (AN11), and target gene dihydroflavonol 4-reductase (DFR) was altered in the different transgenic lines. In addition, ethylene biosynthesis-related genes 1-aminocyclopropane-1-carboxylic acid synthase (ACS1) and ACC oxidase (ACO1) were upregulated in rooting and flower senescence processes but at different time points. Overall, our data demonstrate that the critical role of this AN2 gene in plant growth physiology may extend beyond that of a single activator of anthocyanin biosynthesis.

Selection and Validation of Reference Genes for RT-qPCR Analysis of Gene Expression in Nicotiana benthamiana upon Single Infections by 11 Positive-Sense Single-Stranded RNA Viruses from Four Genera

Quantitative real-time PCR (RT-qPCR) is a widely used method for studying alterations in gene expression upon infections caused by diverse pathogens such as viruses. Positive-sense single-stranded (ss(+)) RNA viruses form a major part of all known plant viruses, and some of them are damaging pathogens of agriculturally important crops. Analysis of gene expression following infection by ss(+) RNA viruses is crucial for the identification of potential anti-viral factors. However, viral infections are known to globally affect gene expression and therefore selection and validation of reference genes for RT-qPCR is particularly important. In this study, the expression of commonly used reference genes for RT-qPCR was studied in Nicotiana benthamiana following single infection by 11 ss(+) RNA viruses, including five tobamoviruses, four potyviruses, one potexvirus and one polerovirus. Stability of gene expression was analyzed in parallel by four commonly used algorithms: geNorm, NormFinder, BestKeeper, and Delta CT, and RefFinder was finally used to summarize all the data. The most stably expressed reference genes differed significantly among the viruses, even when those viruses were from the same genus. Our study highlights the importance of the selection and validation of reference genes upon different viral infections.

Identification and validation of reference genes for qRT-PCR analyses under different experimental conditions in Allium wallichii

Himalayan onion (Allium wallichii) is a perennial bulbous herb with high ornamental value and has long been used as traditional medicines in Nepal and China because of the anti-cancer and anti-microbial activities. Wild Allium wallichii features different flower colors, including purple, pink, deep purple and white. However, little is known about the molecular mechanisms of color formation during A. wallichii flower development stages due to the lack of optimal reference genes. Quantitative real-time polymerase chain reaction (qRT-PCR) is a powerful tool for quantifying expression levels of target genes. The accuracy of qRT-PCR analyses is largely dependent on the identification of stable reference genes for data normalization. The stability of reference gene expression may vary with plant species and environmental conditions. The aim of this study was to select stable reference genes for qRT-PCR analyses of target genes at flower development stages, in different flower colors and organs for Allium wallichii. The CDSs of eight potential reference genes (TUB2, ACT1, GAPC, EF1α, UBQ, UBC, SAND and CYP1) were cloned and their stability was evaluated by four programs (Delta Ct, geNorm, NormFinder and BestKeeper), and the results were further integrated into a comprehensive rank by RefFinder. The results showed that TUB2 and GAPC were the most stable two reference genes at different developmental stages of purple- and white-flower genotypes and across all samples. UBC and TUB2 expression was stable at different developmental stages of purple flowers. CYP1 and TUB2 were stably expressed at different developmental stages of white flowers. GAPC and SAND showed the highest rankings in different flower colors. TUB2 and EF1α performed the best in different tissues. ACT1 was the least stable gene in all tested samples. Moreover, DIHYDROFLAVONOL-4-REDUCTASE (DFR) gene that involved in anthocyanin synthesis was used to evaluate the effectiveness of the selected candidates. This study identified the first set of suitable reference genes for qRT-PCR analyses, which will lay the foundation for gene function study in A. wallichii.

Effect of Savirin or Ticagrelor Treatment on the Expression of Commonly Used Reference Genes in Staphylococcus aureus

Reference genes are frequently used for the normalization of quantitative reverse transcriptase PCR (qRTPCR) data in gene expression studies. Staphylococcus aureus is one of the most common causes of biofilm-related infections. Savirin and ticagrelor show in vitro as well as in vivo antibiofilm activity against S. aureus. The main aim of this study was to identify the most stably expressed reference genes to study the effect of these molecules on genes in a strong biofilm producing S. aureus isolate isolated from biofilm-related infection. Quantitative real-time PCR was performed by using relative quantification method. Four different algorithms, delta Ct, normfinder, bestkeeper, and genorm, followed by a comprehensive analysis was used to identify the most stable reference genes from a list of sixteen different candidate reference genes. All four algorithms reported different results, with some comparable findings among some methods. In the comprehensive analysis of the results of all the algorithms used, the most stable reference genes found were spa, rpoD, and pyk for savirin treatment experiment and gapdH, gyrA, and gmk for ticagrelor treatment experiment. The optimal number of reference genes required was two for both the experimental conditions. Despite having some drawbacks, each algorithm can reliably determine an appropriate reference gene independently. However, based on consensus ranking and the required optimal number of reference genes reported, spa and rpoD were the most appropriate reference genes for savirin treatment experiment, and gapdH and gyrA were most appropriate for ticagrelor treatment experiment. This study provides baseline data on reference genes to study the effect of savirin or ticagrelor treatment on the expression of potential reference genes in S. aureus. We recommend prior re-validation of reference genes on a case-by-case basis before they can be used.

Reference genes selection of Paeonia ostii 'Fengdan' under osmotic stresses and hormone treatments by RT-qPCR

BACKGROUND

Tree peony possess significant ornamental, medicinal and oil values. Osmotic stresses including dehydratiuon and salinity limit the expansion of cultivation area of tree peony. Information on reference genes selection under osmotic stress and hormone stimulation of tree peony still limited. This study aimed to determine the stable reference genes suitable for tree peony under osmotic stresses and hormone treatments, and provide a theoretical basis for the molecular biology research.

METHODS AND RESULTS

Twelve candidate reference genes were evaluated in Paeonia ostii 'Fengdan' under osmotic stress and hormone treatments by RT-qPCR. Delta Ct method, geNorm, and NormFinder were used for the comprehensive expression stability ranking comparison. The results revealed that tubulin-α was the preferred internal reference genes for drought and ABA treatment, tubulin-β was identified as the most suitable reference gene under drought and OPDA induction, 18s-rRNA was regarded as the most stable gene for salinity and JA treatment, eIF-5 A was listed as the most stable gene for JA and MeJA treatments. The experiments also displayed that EF1-α were comparatively unstable under ABA and BR hormone treatments.

CONCLUSION

These preferred reference genes could be useful in qPCR studies involving osmotic or hormonal stresses in Paeonia ostii 'Fengdan'. It is anticipated that the results will benefit tree peony functional genomics studies and molecular breeding research in the future.

Precise measurement of gene expression changes in mouse brain areas denervated by injury

Quantitative PCR (qPCR) is a widely used method to study gene expression changes following brain injury. The accuracy of this method depends on the tissue harvested, the time course analyzed and, in particular on the choice of appropriate internal controls, i.e., reference genes (RGs). In the present study we have developed and validated an algorithm for the accurate normalization of qPCR data using laser microdissected tissue from the mouse dentate gyrus after entorhinal denervation at 0, 1, 3, 7, 14 and 28 days postlesion. The expression stabilities of ten candidate RGs were evaluated in the denervated granule cell layer (gcl) and outer molecular layer (oml) of the dentate gyrus. Advanced software algorithms demonstrated differences in stability for single RGs in the two layers at several time points postlesion. In comparison, a normalization index of several stable RGs covered the entire post-lesional time course and showed high stability. Using these RGs, we validated our findings and quantified glial fibrillary acidic protein (Gfap) mRNA and allograft inflammatory factor 1 (Aif1/Iba1) mRNA in the denervated oml. We compared the use of single RGs for normalization with the normalization index and found that single RGs yield variable results. In contrast, the normalization index gave stable results. In sum, our study shows that qPCR can yield precise, reliable, and reproducible datasets even under such complex conditions as brain injury or denervation, provided appropriate RGs for the model are used. The algorithm reported here can easily be adapted and transferred to any other brain injury model.

Single gene mutation in a plant MYB transcription factor causes a major shift in pollinator preference

Understanding the molecular basis of reproductive isolation and speciation is a key goal of evolutionary genetics. In the South American genus Petunia, the R2R3-MYB transcription factor MYB-FL regulates the biosynthesis of UV-absorbing flavonol pigments, a major determinant of pollinator preference. MYB-FL is highly expressed in the hawkmoth-pollinated P. axillaris, but independent losses of its activity in sister taxa P. secreta and P. exserta led to UV-reflective flowers and associated pollinator shifts in each lineage (bees and hummingbirds, respectively). We created a myb-fl CRISPR mutant in P. axillaris and studied the effect of this single gene on innate pollinator preference. The mutation strongly reduced the expression of the two key flavonol-related biosynthetic genes but only affected the expression of few other genes. The mutant flowers were UV reflective as expected but additionally contained low levels of visible anthocyanin pigments. Hawkmoths strongly preferred the wild-type P. axillaris over the myb-fl mutant, whereas both social and solitary bee preference depended on the level of visible color of the mutants. MYB-FL, with its specific expression pattern, small number of target genes, and key position at the nexus of flavonol and anthocyanin biosynthetic pathways, provides a striking example of evolution by single mutations of large phenotypic effect.

Shikimate Kinase Plays Important Roles in Anthocyanin Synthesis in Petunia

In plants, the shikimate pathway is responsible for the production of aromatic amino acids L-tryptophan, L-phenylalanine, and L-tyrosine. L-Phenylalanine is the upstream substrate of flavonoid and anthocyanin synthesis. Shikimate kinase (SK) catalyzes the phosphorylation of the C3 hydroxyl group of shikimate to produce 3-phosphate shikimate (S3P), the fifth step of the shikimate pathway. However, whether SK participates in flavonoid and anthocyanin synthesis is unknown. This study characterized the single-copy PhSK gene in the petunia (Petunia hybrida) genome. PhSK was localized in chloroplasts. PhSK showed a high transcription level in corollas, especially in the coloring stage of flower buds. Suppression of PhSK changed flower color and shape, reduced the content of anthocyanins, and changed the flavonoid metabolome profile in petunia. Surprisingly, PhSK silencing caused a reduction in the shikimate, a substrate of PhSK. Further qPCR analysis showed that PhSK silencing resulted in a reduction in the mRNA level of PhDHQ/SDH, which encodes the protein catalyzing the third and fourth steps of the shikimate pathway, showing a feedback regulation mechanism of gene expression in the shikimate pathway.

The C2H2-type zinc finger protein PhZFP1 regulates cold stress tolerance by modulating galactinol synthesis in Petunia hybrida

The C2H2 zinc finger proteins (ZFPs) play essential roles in regulating cold stress responses. Similarly, raffinose accumulation contributes to freezing stress tolerance. However, the relationship between C2H2 functions and raffinose synthesis in cold tolerance remains uncertain. Here, we report the characterization of the cold-induced C2H2-type zinc finger protein PhZFP1 in Petunia hybrida. PhZFP1 was found to be predominantly localized in the nucleus. Overexpression of PhZFP1 conferred enhanced cold tolerance in transgenic petunia lines. In contrast, RNAi mediated suppression of PhZFP1 led to increased cold susceptibility. PhZFP1 regulated the expression of a range of abiotic stress responsive-genes including genes encoding proteins involved in reactive oxygen species (ROS) scavenging and raffinose metabolism. The accumulation of galactinol and raffinose, and the levels of PhGolS1-1 transcripts, were significantly increased in PhZFP1-overexpressing plants and decreased in PhZFP1-RNAi plants under cold stress. Moreover, the galactinol synthase (GolS)-encoding gene PhGolS1-1 was identified as a direct target of PhZFP1. Taken together, these results demonstrate that PhZFP1 functions in cold stress tolerance by modulation of galactinol synthesis via regulation of PhGolS1-1. This study also provides new insights into the mechanisms underlying C2H2 zinc finger protein-mediated cold stress tolerance, and has identified a candidate gene for improving cold stress tolerance.

Crotonylation versus acetylation in petunia corollas with reduced acetyl-CoA due to PaACL silencing

Protein acetylation and crotonylation are important posttranslational modifications of lysine. In animal cells, the correlation of acetylation and crotonylation has been well characterized and the lysines of some proteins are acetylated or crotonylated depending on the relative concentrations of acetyl-CoA and crotonyl-CoA. However, in plants, the correlation of acetylation and crotonylation and the effects of the relative intracellular concentrations of crotonyl-CoA and acetyl-CoA on protein crotonylation and acetylation are not well known. In our previous study, PaACL silencing changed the content of acetyl-CoA in petunia (Petunia hybrida) corollas, and the effect of PaACL silencing on the global acetylation proteome in petunia was analyzed. In the present study, we found that PaACL silencing did not significantly alter the content of crotonyl-CoA. We performed a global crotonylation proteome analysis of the corollas of PaACL-silenced and control petunia plants; we found that protein crotonylation was closely related to protein acetylation and that proteins with more crotonylation sites often had more acetylation sites. Crotonylated proteins and acetylated proteins were enriched in many common Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. However, PaACL silencing resulted in different KEGG pathway enrichments of proteins with different levels of crotonylation sites and acetylation sites. PaACLB1-B2 silencing did not led to changes in the opposite direction in crotonylation and acetylation levels at the same lysine site in cytoplasmic proteins, which indicated that cytoplasmic lysine acetylation and crotonylation might not depend on the relative concentrations of acetyl-CoA and crotonyl-CoA. Moreover, the global crotonylome and acetylome were weakly positively correlated in the corollas of PaACL-silenced and control plants.

Ph5GT silencing alters flower color and flavonoids metabolome profile in petunia

Anthocyanins are important pigments in plants and glycosylation plays an important role in the stability of anthocyanins. Anthocyanin 5-O-glucosyltransferase (5GT) can glycosylate anthocyanin at the 5-O-position. Till now, the enzymatic activity characteristics of 5GT had been studied in vitro in a variety of plants. However, the subcellular localization of 5GT protein still remained unclear, and little genetic evidence on the roles of 5GT in plants has been reported. The full-length Ph5GT gene from petunia (Petunia hybrida) was isolated in this study. Green fluorescent fusion protein assays revealed that Ph5GT protein was localized to the cytoplasm. Ph5GT was found to be highly expressed in flowers, with highest levels of expression occurring during the coloring stage of flower development. Furthermore, Ph5GT silencing led to the change in flower color from purple to light purple and a significant reduction in total anthocyanin content. The metabolome analysis revealed that the content of malvidins and petunidins modified by glycosylation at the 5-O-position was significantly reduced, while the content of their precursor without glycosylation was significantly increased, implying that Ph5GT could glycosylate malvidin and petunidin derivatives and that the substrate types of Ph5GT were expanded in comparison to previous studies.

PhRHMs play important roles in leaf and flower development and anthocyanin synthesis in petunia

Anthocyanins, vital metabolites in plants, are formed by anthocyanidins combined with various monosaccharides, including glucose, rhamnose, and arabinose. Rhamnose contributes greatly to the glycosylation of anthocyanidins. There are two kinds of rhamnose synthase (RS): rhamnose biosynthesis (RHM), and nucleotide-RS/epimerase-reductase (UER1). Nevertheless, no RS isoform was reported to be involved in anthocyanin synthesis. Here, three homologous PhRHM genes, namely PhRHM1, PhRHM2, and PhRHM3, and one PhUER1 gene from petunia were cloned and characterized. Green fluorescent protein fusion protein assays revealed that PhRHMs and PhUER1 are localized in the cytoplasm. We obtained PhRHM1 or/and PhRHM2 or PhUER1 silenced petunia plants and did not attempt to obtain PhRHM3 silenced plants since PhRHM3 mRNA was not detected in petunia organs examined. PhRHM1 and PhRHM2 (PhRHM1-2) silencing induced abnormal plant growth and decreased the contents of l-rhamnose, photosynthetic pigments and total anthocyanins, while PhUER1 silencing did not cause any visible phenotypic changes. Flavonoid metabolome analysis further revealed that PhRHM1-2 silencing reduced the contents of anthocyanins with rhamnose residue. These results revealed that PhRHMs contribute to the biosynthesis of rhamnose and that PhRHMs participate in the anthocyanin rhamnosylation in petunia, while PhUER1 does not.

Systematic screening and validation of reliable reference genes for qRT-PCR analysis in Okra (Abelmoschus esculentus L.)

Quantitative real-time polymerase chain reaction (qRT-PCR) is a sensitive and widely used technique for quantifying gene expression levels, and its accuracy depends on the reference genes used for data normalization. To date, no reference gene has been reported in the nutritious and functional vegetable okra (Abelmoschus esculentus L.). Herein, 11 candidates of reference genes were selected and evaluated for their expression stability in okra in different tissues at different developmental stages by using three software algorithms (geNorm, NormFinder, BestKeeper) and a web-based tool (RefFinder). Among them, eukaryotic initiation factor 4 alpha (eIF4A) and protein phosphatase 2A (PP2A) showed the highest stability, while TUA5 had the lowest stability. The combined usage of these two most stable reference genes was sufficient to normalize gene expression in okra. Then, the above results were further validated by normalizing the expression of the cellulose synthase gene CesA4. This work provides appropriate reference genes for transcript normalization in okra, which will facilitate subsequent functional gene research on this vegetable crop.

Selection and validation of reference genes for quantitative real-time PCR normalization in Psoralea corylifolia (Babchi) under various abiotic stress

Psoralea corylifolia L. is a popular herb and has long been used in traditional Ayurvedic and Chinese medicine owing to its extensive pharmacological activities, especially in the treatment of various shin diseases. To date, the systematic evaluation and selection of the optimum reference genes for gene expression analysis of P. corylifolia were not reported. Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) is a method for gene expression quantification. Selecting appropriate reference genes is the prerequisite for accurate normalization of RT-qPCR results. This study assessed the expression stability of 10 candidate reference genes under different abiotic stresses. First, amplification primers for RT-qPCR were designed and received testing and optimization. Then, expression data from each candidate gene were evaluated based on three statistical algorithms, and their results were further integrated into a comprehensive ranking based on the geometric mean. Additionally, one target gene, i.e., 1-aminocyclopropane-1-carboxylate oxidase (ACO), was used to validate the effectiveness of the selected reference. Our analysis suggested that thioredoxin-like protein YLS8 (YLS8), TIP41-like family protein (TIP41), and cyclophilin 2 (CYP2) genes provided superior expression normalization under different abiotic stresses. Overall, this work constitutes the first effort to select optimal endogenous controls for RT-qPCR studies of gene expression in P. corylifolia. It also provides a reasonable normalization standard and basis for further analysis of the gene expression of bioactive components in P. corylifolia.

Selection and Validation of Reference Genes for RT-qPCR Analysis in Tibetan Medicinal Plant Saussurea Laniceps Callus under Abiotic Stresses and Hormone Treatments

Real-time quantitative PCR (RT-qPCR) is an important technique for studying gene expression analysis, but accurate and reliable results depend on the use of a stable reference gene. This study proposes to test the expression stability of candidate reference genes in the callus of Saussurea laniceps, a unique Tibetan medicinal plant. Based on the S. laniceps callus transcriptome, eleven candidate reference genes, including TUA2, TUB3, TUB8, TIF3B1, TIF3H1, ELF5A, PP2AA2, UEV1D, UBL5, UBC36, and SKIP1), were validated for RT-qPCR normalization in the callus under abiotic stress (salt, cold, and UV) and hormone treatments (abscisic acid, MeJA, and salicylic acid). The stability of the candidate genes was evaluated in all the samples of S. laniceps. Comprehensive analysis of all samples showed that the best reference genes were UBC36 and UBL5. ELF5A and TIF3B1 were ranked as the most stable genes in the sample sets under abiotic stress. For hormone stimulation, UBC36 and TIF3H1 genes had the best stability. This study provides useful guidelines and a starting point for reference gene selection for expression analysis using RT-qPCR techniques in S. laniceps.

CRISPR/Cas9-Mediated Editing of Autophagy Gene 6 in Petunia Decreases Flower Longevity, Seed Yield, and Phosphorus Remobilization by Accelerating Ethylene Production and Senescence-Related Gene Expression

Developmental petal senescence is a type of programmed cell death (PCD), during which the production of ethylene is induced, the expression of PCD-related genes is upregulated, and nutrients are recycled. Autophagy is an intracellular mechanism involved in PCD modulation and nutrient cycling. As a central component of the autophagy pathway, Autophagy Gene 6 (ATG6) was previously shown as a negative regulator of petal senescence. To better understand the role of autophagy in ethylene biosynthesis and nutrient remobilization during petal senescence, we generated and characterized the knockout (KO) mutants of PhATG6 using CRISPR/Cas9 in Petunia × hybrida 'Mitchell Diploid.' PhATG6-KO lines exhibited decreased flower longevity when compared to the flowers of the wild-type or a non-mutated regenerative line (controls), confirming the negative regulatory role of ATG6 in petal senescence. Smaller capsules and fewer seeds per capsule were produced in the KO plants, indicating the crucial function of autophagy in seed production. Ethylene production and ethylene biosynthesis genes were upregulated earlier in the KO lines than the controls, indicating that autophagy affects flower longevity through ethylene. The transcript levels of petal PCD-related genes, including PhATG6, PhATG8d, PhPI3K (Phosphatidylinositol 3-Kinase), and a metacaspase gene PhMC1, were upregulated earlier in the corollas of PhATG6-KO lines, which supported the accelerated PCD in the KO plants. The remobilization of phosphorus was reduced in the KO lines, showing that nutrient recycling was compromised. Our study demonstrated the important role of autophagy in flower lifespan and seed production and supported the interactions between autophagy and various regulatory factors during developmental petal senescence.

ODORANT1 targets multiple metabolic networks in petunia flowers

Scent bouquets produced by the flowers of Petunia spp. (petunia) are composed of a complex mixture of floral volatile benzenoid and phenylpropanoid compounds (FVBPs), which are specialized metabolites derived from phenylalanine (Phe) through an interconnected network of enzymes. The biosynthesis and emission of high levels of these volatiles requires coordinated transcriptional activation of both primary and specialized metabolic networks. The petunia R2R3-MYB transcription factor ODORANT 1 (ODO1) was identified as a master regulator of FVBP production and emission; however, our knowledge of the direct regulatory targets of ODO1 has remained limited. Using chromatin immunoprecipitation followed by sequencing (ChIP-seq) in petunia flowers, we identify genome-wide ODO1-bound genes that are enriched not only in genes involved in the biosynthesis of the Phe precursor, as previously reported, but also genes associated with the specialized metabolic pathways involved in generating phenylpropanoid intermediates for FVBPs. ODO1-bound genes are also involved in methionine and S-adenosylmethionine metabolism, which could modulate methyl group supplies for certain FVBPs. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) and RNA-seq analysis in an ODO1 RNAi knockdown line revealed that ODO1-bound targets are expressed at lower levels when ODO1 is suppressed. A cis-regulatory motif, CACCAACCCC, was identified as a potential binding site for ODO1 in the promoters of genes that are both bound and activated by ODO1, which was validated by in planta promoter reporter assays with wild-type and mutated promoters. Overall, our work presents a mechanistic model for ODO1 controlling an extensive gene regulatory network that contributes to FVBP production to give rise to floral scent.

Evaluation of Reference Genes for Quantitative PCR in Eustoma grandiflorum under Different Experimental Conditions

Eustoma grandiflorum, commonly known as prairie gentian or Texas bluebells, is among the most popular agriculturally propagated species of cut flowers. Due to its widespread appeal, there is increasing interest in understanding the molecular genetic factors underlying floral development and resistance to abiotic stresses. We analyzed 18 potential reference genes in different organs, at different floral developmental stages and under drought- and salt-stress treatments, for use in RT-qPCR analysis. A total of four analytical tool packages, including geNorm, NormFinder, BestKeeper, and RefFinder were employed to determine the most appropriate reference genes under each treatment condition. The results demonstrate that different reference genes should be used for normalization under different experimental treatments. EgPP and EgPP2A2 were the most stable internal control genes across different organ types, EgPP and Eg18S were the most stable under salt-stress, EgPP and EgACT1 were the most stable across different floral development stages, and EgEF1A and EgTUA were the most stable reference genes under drought-stress. Additional gene expression analyses of EgMIXTA1, EgTOE1, and EgP5CS1 further confirmed the applicability of these reference genes. The results represent a significant contribution to future studies of reference gene selection for the normalization of gene expression in Eustoma grandiflorum.

Transcriptome Profiling of Cu Stressed Petunia Petals Reveals Candidate Genes Involved in Fe and Cu Crosstalk

Copper (Cu) is an essential element for most living plants, but it is toxic for plants when present in excess. To better understand the response mechanism under excess Cu in plants, especially in flowers, transcriptome sequencing on petunia buds and opened flowers under excess Cu was performed. Interestingly, the transcript level of FIT-independent Fe deficiency response genes was significantly affected in Cu stressed petals, probably regulated by basic-helix-loop-helix 121 (bHLH121), while no difference was found in Fe content. Notably, the expression level of bHLH121 was significantly down-regulated in petals under excess Cu. In addition, the expression level of genes related to photosystem II (PSII), photosystem I (PSI), cytochrome b₆/f complex, the light-harvesting chlorophyll II complex and electron carriers showed disordered expression profiles in petals under excess Cu, thus photosynthesis parameters, including the maximum PSII efficiency (F_V/F_M), nonphotochemical quenching (NPQ), quantum yield of the PSII (ΦPS(II)) and photochemical quenching coefficient (qP), were reduced in Cu stressed petals. Moreover, the chlorophyll a content was significantly reduced, while the chlorophyll b content was not affected, probably caused by the increased expression of chlorophyllide a oxygenase (CAO). Together, we provide new insight into excess Cu response and the Cu–Fe crosstalk in flowers.

Evaluation of Angelica decursiva reference genes under various stimuli for RT-qPCR data normalization

Angelica decursiva is one of the lending traditional Chinese medicinal plants producing coumarins. Notably, several studies have focused on the biosynthesis and not the RT-qPCR (quantitative real-time reverse transcription polymerase chain reaction) study of coumarins. This RT-qPCR technique has been extensively used to investigate gene expression levels in plants and the selection of reference genes which plays a crucial role in standardizing the data form the RT-qPCR analysis. In our study, 11 candidate reference genes were selected from the existing transcriptome data of Angelica decursiva. Here, four different types of statistical algorithms (geNorm, NormFinder, BestKeeper, and Delta Ct) were used to calculate and evaluate the stability of gene expression under different external treatments. Subsequently, RefFinder analysis was used to determine the geometric average of each candidate gene ranking, and to perform comprehensive index ranking. The obtained results showed that among all the 11 candidate reference genes, SAND family protein (SAND), protein phosphatase 2A gene (PP2A), and polypyrimidine tract-binding protein (PTBP) were the most stable reference genes, where Nuclear cap binding protein 2 (NCBP2), TIP41-like protein (TIP41), and Beta-6-tubulin (TUBA) were the least stable genes. To the best of our knowledge, this work is the first to evaluate the stability of reference genes in the Angelica decursiva which has provided an important foundation on the use of RT-qPCR for an accurate and far-reaching gene expression analysis in this medicinal plant.

Complex evolution of novel red floral color in Petunia

Red flower color has arisen multiple times and is generally associated with hummingbird pollination. The majority of evolutionary transitions to red color proceeded from purple lineages and tend to be genetically simple, almost always involving a few loss-of-function mutations of major phenotypic effect. Here we report on the complex evolution of a novel red floral color in the hummingbird-pollinated Petunia exserta (Solanaceae) from a colorless ancestor. The presence of a red color is remarkable because the genus cannot synthesize red anthocyanins and P. exserta retains a nonfunctional copy of the key MYB transcription factor AN2. We show that moderate upregulation and a shift in tissue specificity of an AN2 paralog, DEEP PURPLE, restores anthocyanin biosynthesis in P. exserta. An essential shift in anthocyanin hydroxylation occurred through rebalancing the expression of three hydroxylating genes. Furthermore, the downregulation of an acyltransferase promotes reddish hues in typically purple pigments by preventing acyl group decoration of anthocyanins. This study presents a rare case of a genetically complex evolutionary transition toward the gain of a novel red color.

Expression stability of internal reference gene in response to Trichoderma polysporum infection in Avena fatua L

Trichoderma polysporum was a pathogenic fungi which showed strong pathogenicity to Avena fatua L. in recent study. The stress response of A. fatua to T. polysporum is mediated by the regulation of gene expression. Quantification of the expression of genes requires normalizing RT-qPCR data using reference genes with stable expression in the system studied as internal standards. To construct a RT-qPCR system suitable for response of A. fatua to T. polysporum, and screen stable internal reference genes, GeNorm, NormFinder, BestKeeper and RefFinde were used to perform SYBR Green-based RT-qPCR analysis on eight candidate internal reference genes (18S, 28S, TUA, UBC, ACT, GAPDH, TBP and EF-1α) in A. fatua samples after inoculation of T. polysporum Strain HZ-31. The results showed that TBP, 18S and UBC were the most stable internal reference genes, TBP and TUA, TBP and GAPDH, 18S and TBP, UBC and 18S were the most suitable combination of the two internal reference genes, which could be used as internal reference genes for functional gene expression analysis during the interaction between T. polysporum and A. fatua. This is the first study describing a set of reference genes with a stable expression under fungi stress in A. fatua. These genes are also candidate reference genes of choice for studies seeking to identify stress-responsive genes in A. fatua.

Validation of reference genes for whole blood gene expression analysis in cord blood of preterm and full-term neonates and peripheral blood of healthy adults

Background

Preterm birth is the leading cause of neonatal morbidity and mortality, but research efforts in neonatology are complicated due to the unavailability of large volume blood samples. Whole blood assays can be used to overcome this problem by performing both functional and gene expression studies using small amounts of blood. Gene expression studies using RT-qPCR estimate mRNA-levels of target genes normalized to reference genes. The goal of this study was to identify and validate stable reference genes applicable to cord blood samples obtained from developing neonates of different gestational age groups as well as to adult peripheral blood samples. Eight reference gene candidates (ACTB, B2M, GAPDH, GUSB, HPRT, PPIB, RPLP0, RPL13) were analyzed using the three published software algorithms Bestkeeper, GeNorm and NormFinder.

Results

A normalization factor consisting of ACTB and PPIB allows for comparative expression analyses of neonatal samples from different gestational age groups. Normalization factors consisting of GAPDH and PPIB or ACTB and GAPDH are suitable when samples from preterm and full-term neonates and adults are compared. However, all candidate reference genes except RPLP0 exhibited significant intergroup gene expression variance and a higher gene expression towards an older age which resulted in a small but statistically significant systematic bias. Systematic analysis of RNA-seq data revealed new reference gene candidates with potentially superior stability.

Conclusions

The current study identified suitable normalization factors and proposed the use of the additional single gene RPLP0 to avoid systematic bias. This combination will enable comparative analyses not only between neonates of different gestational ages, but also between neonates and adults, as it facilitates more detailed investigations of developmental gene expression changes. The use of software algorithms did not prevent unintended systematic bias. This generally highlights the need for careful validation of such results to prevent false interpretation of potential age-dependent changes in gene expression. To identify the most stable reference genes in the future, RNA-seq based global approaches are recommended.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07801-0.

CRISPR/Cas9-mediated targeted mutation reveals a role for AN4 rather than DPL in regulating venation formation in the corolla tube of Petunia hybrida

Venation is a common anthocyanin pattern displayed in flowers that confers important ornamental traits to plants. An anthocyanin-related R2R3-MYB transcription factor, DPL, has been proposed to regulate corolla tube venation in petunia plants. Here, however, we provide evidence redefining the role of DPL in petunia. A CRISPR/Cas9-mediated mutation of DPL resulted in the absence of the vein-associated anthocyanin pattern above the abaxial surface of the flower bud, but not corolla tube venation, thus indicating that DPL did not regulate the formation of corolla tube venation. Alternately, quantitative real-time PCR analysis demonstrated that the spatiotemporal expression pattern of another R2R3-MYB gene, AN4, coincided with the formation of corolla tube venation in petunia. Furthermore, overexpression of AN4 promoted anthocyanin accumulation by increasing the expression of anthocyanin biosynthesis genes. CRISPR/Cas9-mediated mutation of AN4 led to an absence of corolla tube venation, suggesting that this gene in fact determines this key plant trait. Taken together, the results presented here redefine the prime regulator of corolla tube venation, paving the way for further studies on the molecular mechanisms underlying the various venation patterns in petunia.

Selection and validation of appropriate reference genes for RT-qPCR analysis of flowering stages and different genotypes of Iris germanica L

Iris germanica L. is a perennial herbaceous plant that has been widely cultivated worldwide and is popular for its elegant and vibrantly colorful flowers. Selection of appropriate reference genes is the prerequisite for accurate normalization of target gene expression by quantitative real-time PCR. However, to date, the most suitable reference genes for flowering stages have not been elucidated in I. germanica. In this study, eight candidate reference genes were examined for the normalization of RT-qPCR in three I. germanica cultivars, and their stability were evaluated by four different algorithms (GeNorm, NormFinder, BestKeeper, and Ref-finder). The results revealed that IgUBC and IgGAPDH were the most stable reference genes in '00246' and 'Elizabeth', and IgTUB and IgUBC showed stable expression in '2010200'. IgUBC and IgGAPDH were the most stable in all samples, while IgUBQ showed the least stability. Finally, to validate the reliability of the selected reference genes, the expression patterns of IgFT (Flowering Locus T gene) was analyzed and emphasized the importance of appropriate reference gene selection. This work presented the first systematic study of reference genes selection during flower bud development and provided guidance to research of the molecular mechanisms of flowering stages in I. germanica.

Time-course transcriptomic analysis of Petunia ×hybrida leaves under water deficit stress using RNA sequencing

Water deficit limits plant growth and development, resulting in quality loss of horticultural crops. However, there is limited information on gene regulation and signaling pathways related to water deficit stress response at multiple time points. The objective of this research was to investigate global gene expression patterns under water deficit stress to provide an insight into how petunia (Petunia ×hybrida 'Mitchell Diploid') responded in the process of stress. Nine-week-old petunias were irrigated daily or placed under water stress by withholding water. Stressed plants reduced stomatal conductance after five days of water deficit, indicating they perceived stress and initiated stress response mechanisms. To analyze transcriptomic changes at the early stage of water deficit, leaf tissue samples were collected 1, 3, and 5 days after water was withheld for RNA sequencing. Under water deficit stress, 154, 3611, and 980 genes were upregulated and 41, 2806, and 253 genes were downregulated on day 1, 3, and 5, respectively. Gene Ontology analysis revealed that redox homeostasis processes through sulfur and glutathione metabolism pathways, and hormone signal transduction, especially abscisic acid and ethylene, were enriched under water deficit stress. Thirty-four transcription factor families were identified, including members of AP2/ERF, NAC, MYB-related, C2H2, and bZIP families, and TFs in AP2/ERF family was the most abundant in petunia. Interestingly, only one member of GRFs was upregulated on day 1, while most of TFs were differentially expressed on day 3 and/or 5. The transcriptome data from this research will provide valuable molecular resources for understanding the early stages of water stress-responsive networks as well as engineering petunia with enhanced water stress tolerance.

Mitochondrial citrate synthase plays important roles in anthocyanin synthesis in petunia

Anthocyanins are important flavonoid pigments in plants. Malonyl CoA is an important intermediate in anthocyanin synthesis, and citrate, formed by citrate synthase (CS) catalysing oxaloacetate, is the precursor for the formation of malonyl-CoA. CS is composed of two isoforms, mitochondrial citrate synthase (mCS), a key enzyme of the tricarboxylic acid (TCA) cycle, and citrate synthase (CSY) localizated in microbodies in plants. However, no CS isoform involvement in anthocyanin synthesis has been reported. In this study, we identified the entire CS family in petunia (Petunia hybrida): PhmCS, PhCSY1 and PhCSY2. We obtained petunia plants silenced for the three genes. PhmCS silencing resulted in abnormal development of leaves and flowers. The contents of citrate and anthocyanins were significantly reduced in flowers in PhmCS-silenced plants. However, silencing of PhCSY1 and/or PhCSY2 did not cause a visible phenotype change in petunia. These results showed that PhmCS is involved in anthocyanin synthesis and the development of leaves and flowers, and that the citrate involved in anthocyanin synthesis mainly derived from mitochondria rather than microbodies in petunia.

Two plastidic glycolate/glycerate translocator 1 isoforms function together to transport photorespiratory glycolate and glycerate in rice chloroplasts

The photorespiratory pathway is highly compartmentalized. As such, metabolite shuttles between organelles are critical to ensure efficient photorespiratory carbon flux. Arabidopsis plastidic glycolate/glycerate translocator 1 (PLGG1) has been reported as a key chloroplastic glycolate/glycerate transporter. Two homologous genes, OsPLGG1a and OsPLGG1b, have been identified in the rice genome, although their distinct functions and relationships remain unknown. Herein, our analysis of exogenous expression in oocytes and yeast shows that both OsPLGG1a and OsPLGG1b have the ability to transport glycolate and glycerate. Furthermore, we demonstrate in planta that the perturbation of OsPLGG1a or OsPLGG1b expression leads to extensive accumulation of photorespiratory metabolites, especially glycolate and glycerate. Under ambient CO2 conditions, loss-of-function osplgg1a or osplgg1b mutant plants exhibited significant decreases in photosynthesis efficiency, starch accumulation, plant height, and crop productivity. These morphological defects were almost entirely recovered when the mutant plants were grown under elevated CO2 conditions. In contrast to osplgg1a, osplgg1b mutant alleles produced a mild photorespiratory phenotype and had reduced accumulation of photorespiratory metabolites. Subcellular localization analysis showed that OsPLGG1a and OsPLGG1b are located in the inner and outer membranes of the chloroplast envelope, respectively. In vitro and in vivo experiments revealed that OsPLGG1a and OsPLGG1b have a direct interaction. Our results indicate that both OsPLGG1a and OsPLGG1b are chloroplastic glycolate/glycerate transporters required for photorespiratory metabolism and plant growth, and that they may function as a singular complex.

Validation of suitable genes for normalization of diurnal gene expression studies in Chenopodium quinoa

Quinoa depicts high nutritional quality and abiotic stress resistance, attracting strong interest in the last years. To unravel the function of candidate genes for agronomically relevant traits, studying their transcriptional activities by RT-qPCR is an important experimental approach. The accuracy of such experiments strongly depends on precise data normalization. To date, validation of potential candidate genes for normalization of diurnal expression studies has not been performed in C. quinoa. We selected eight candidate genes based on transcriptome data and literature survey, including conventionally used reference genes. We used three statistical algorithms (BestKeeper, geNorm and NormFinder) to test their stability and added further validation by a simulation-based strategy. We demonstrated that using different reference genes, including those top ranked by stability, causes significant differences among the resulting diurnal expression patterns. Our results show that isocitrate dehydrogenase enzyme (IDH-A) and polypyrimidine tract-binding protein (PTB) are suitable genes to normalize diurnal expression data of two different quinoa accessions. Moreover, we validated our reference genes by normalizing two known diurnally regulated genes, BTC1 and BBX19. The validated reference genes obtained in this study will improve the accuracy of RT-qPCR data normalization and facilitate gene expression studies in quinoa.

Bacillus pumilus improved drought tolerance in Glycyrrhiza uralensis G5 seedlings through enhancing primary and secondary metabolisms

It has been reported that drought stress adversely affects the growth and yield of Glycyrrhiza uralensis, Chinese liquorice, in agricultural production. Bacillus pumilus, an important plant growth-promoting bacterium, play a significant role in improving plant tolerance to abiotic stress. However, the role of Bacillus pumilus G5 in resisting drought stress is largely unknown. In the present study, we found that drought stress significantly inhibited the growth and reduced the biomass of G. uralensis seedlings by restraining C- and N-metabolism, while this could be effectively reversed by B. pumilus G5 inoculation. Specifically, B. pumilus G5 significantly increased the content of primary metabolites such as soluble sugar, soluble protein, and free amino acids by regulating the C and N metabolic processes in G. uralensis seedlings. Moreover, B. pumilus G5 increased the content of glycyrrhizic acid, one of the important secondary metabolites, likely mediated through the increased content of primary metabolites and by recovering the expression of three key enzymes, HMGR, SQS, and β-AS, in the biosynthesis of glycyrrhizic acid. Interestingly, the regulating effect of B. pumilus G5 inoculation on promoting the accumulation of glycyrrhizic acid and increasing the expression of synthesis-related genes is spatially selective. In summary, our findings suggest that B. pumilus G5 could alleviate adverse effects induced by drought stress on the growth of G. uralensis seedlings by regulating C- and N-metabolisms that further triggered the accumulation of secondary metabolites, and this finally improved the drought tolerance of cultivated G. uralensis seedlings.

Selection and Validation of Reference Genes for RT-qPCR Analysis in Spinacia oleracea under Abiotic Stress

Reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) is an accurate and convenient method for mRNA quantification. Selection of optimal reference gene(s) is an important step in RT-qPCR experiments. However, the stability of housekeeping genes in spinach (Spinacia oleracea) under various abiotic stresses is unclear. Evaluating the stability of candidate genes and determining the optimal gene(s) for normalization of gene expression in spinach are necessary to investigate the gene expression patterns during development and stress response. In this study, ten housekeeping genes, 18S ribosomal RNA (18S rRNA), actin, ADP ribosylation factor (ARF), cytochrome c oxidase subunit 5C (COX), cyclophilin (CYP), elongation factor 1-alpha (EF1α), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), histone H3 (H3), 50S ribosomal protein L2 (RPL2), and tubulin alpha chain (TUBα) from spinach, were selected as candidates in roots, stems, leaves, flowers, and seedlings in response to high temperature, CdCl₂, NaCl, NaHCO₃, and Na₂CO₃ stresses. The expression of these genes was quantified by RT-qPCR and evaluated by NormFinder, BestKeeper, and geNorm. 18S rRNA, actin, ARF, COX, CYP, EF1α, GAPDH, H3, and RPL2 were detected as optimal reference genes for gene expression analysis of different organs and stress responses. The results were further confirmed by the expression pattern normalized with different reference genes of two heat-responsive genes. Here, we optimized the detection method of the gene expression pattern in spinach. Our results provide the optimal candidate reference genes which were crucial for RT-qPCR analysis.

Silencing ATG6 and PI3K accelerates petal senescence and reduces flower number and shoot biomass in petunia

Petal senescence is a form of developmental programmed cell death (PCD) that is regulated by internal and environmental signals. Autophagy, a metabolic pathway that regulates intercellular nutrient recycling, is thought to play an important role in the regulation of petal senescence-associated PCD. To characterize the function of two central autophagy genes in petal senescence, we down-regulated Autophagy Gene 6 (PhATG6) and Phosphoinositide 3-Kinase (PhPI3K) using Virus-Induced Gene Silencing (VIGS) in Petunia × hybrida. The silencing of PhATG6 and PhPI3K accelerated petal senescence, thereby reducing flower longevity. Both PhATG6- and PhPI3K-silenced petunias had reduced flower numbers, flower biomass, and vegetative shoot biomass. These phenotypes were intensified when plants were grown under low nutrient conditions. Additionally, two important regulators of senescence, an ethylene biosynthesis gene (PhACS) and a type I metacaspase gene (PhMC1), were suppressed in senescing petals of PhATG6- and PhPI3K-silenced plants. In conclusion, our study identified PhATG6 and PhPI3K as negative regulators of flower senescence and demonstrated the influence of nutrient limitation on the function of autophagy during petal senescence. Our study also found that autophagy genes potentially influence the transcriptional regulation of metacaspases and ethylene biosynthetic genes during petal senescence. The results of this project will be fundamental for future studies of petal senescence and will provide genetic information for future crop improvement.

The Regulation of Floral Colour Change in Pleroma raddianum (DC.) Gardner

Floral colour change is a widespread phenomenon in angiosperms, but poorly understood from the genetic and chemical point of view. This article investigates this phenomenon in Pleroma raddianum, a Brazilian endemic species whose flowers change from white to purple. To this end, flavonoid compounds and their biosynthetic gene expression were profiled. By using accurate techniques (Ultra Performance Liquid Chromatography-High-Resolution Mass Spectrometry (UPLC-HRMS)), thirty phenolic compounds were quantified. Five key genes of the flavonoid biosynthetic pathway were partially cloned, sequenced, and the mRNA levels were analysed (RT-qPCR) during flower development. Primary metabolism was also investigated by gas chromatography coupled to mass spectrometry (GC-EIMS), where carbohydrates and organic acids were identified. Collectively, the obtained results suggest that the flower colour change in P. raddianum is determined by petunidin and malvidin whose accumulation coincides with the transcriptional upregulation of early and late biosynthetic genes of the flavonoid pathway, mainly CHS and ANS, respectively. An alteration in sugars, organic acids and phenolic co-pigments is observed together with the colour change. Additionally, an increment in the content of Fe³⁺ ions in the petals, from the pink to purple stage, seemed to influence the saturation of the colour.

Selection and validation of reference genes for measuring gene expression in Toona ciliata under different experimental conditions by quantitative real-time PCR analysis

BACKGROUND

Before studying gene expression of different organisms, it is important to determine the best reference gene. At present, the most accurate method of detecting gene expression is quantitative real-time PCR (RT-qPCR). With this method, reference genes that are stable in different biological systems and under different conditions can be obtained. Toona ciliata Roem (T. ciliata). is a valuable and fast-growing timber specie. In this study, 20 reference genes were identified using RT-qPCR, as a primary prerequisite for future gene expression analysis. Four different methods, geNorm, NormFinder, BestKeeper, and RankAggreg were used to evaluate the expression stability of the 20 candidate reference genes in various tissues under different conditions.

RESULTS

The experimental results showed that TUB-α was the most stably expressed reference gene across all samples and UBC17 was the most stable in leaves and young stems under Hypsipyla robusta (H. robusta) and methyl jasmonate (MeJA) treatments. In addition, PP2C59 and UBC5B were the best-performing genes in leaves under H. robusta treatment, while HIS1 and ACT7 were the best reference genes in young stems. The two best reference genes were 60S-18 and TUB-α after treatment at 4 °C. The expression of HIS6 and MUB1 was the most stable under PEG6000 treatment. The accuracy of the selected reference genes was verified using the transcription factor MYB3 (TcMYB3) gene.

CONCLUSIONS

This is the first report to verify the best reference genes for normalizing gene expression in T. ciliata under different conditions, which will facilitate future elucidation of gene regulations in this species.

PaACL silencing accelerates flower senescence and changes the proteome to maintain metabolic homeostasis in Petunia hybrida

Cytosolic acetyl-CoA is an intermediate of the synthesis of most secondary metabolites and the source of acetyl for protein acetylation. The formation of cytosolic acetyl-CoA from citrate is catalysed by ATP-citrate lyase (ACL). However, the function of ACL in global metabolite synthesis and global protein acetylation is not well known. Here, four genes, PaACLA1, PaACLA2, PaACLB1, and PaACLB2, which encode the ACLA and ACLB subunits of ACL in Petunia axillaris, were identified as the same sequences in Petunia hybrida 'Ultra'. Silencing of PaACLA1-A2 and PaACLB1-B2 led to abnormal leaf and flower development, reduced total anthocyanin content, and accelerated flower senescence in petunia 'Ultra'. Metabolome and acetylome analysis revealed that PaACLB1-B2 silencing increased the content of many downstream metabolites of acetyl-CoA metabolism and the levels of acetylation of many proteins in petunia corollas. Mechanistically, the metabolic stress induced by reduction of acetyl-CoA in PaACL-silenced petunia corollas caused global and specific changes in the transcriptome, the proteome, and the acetylome, with the effect of maintaining metabolic homeostasis. In addition, the global proteome and acetylome were negatively correlated under acetyl-CoA deficiency. Together, our results suggest that ACL acts as an important metabolic regulator that maintains metabolic homeostasis by promoting changes in the transcriptome, proteome. and acetylome.