Identification of superior reference genes for data normalisation of expression studies via quantitative PCR in hybrid roses (Rosa hybrida)

Brevipalpus yothersi Baker (Tenuipalpidae) development in sweet orange plants is influenced by previous mite infestation and the presence of shelters

Citrus leprosis is the most important viral disease affecting citrus. The disease is caused predominantly by CiLV-C and is transmitted by Brevipalpus yothersi Baker mites. This study brings some insight into the colonization of B. yothersi in citrus [(Citrus × sinensis (L.) Osbeck (Rutaceae)] previously infested by viruliferous or non-viruliferous B. yothersi. It also assesses the putative role of shelters on the behavior of B. yothersi. Expression of PR1 and PR4 genes, markers of plant defense mechanisms, were evaluated by RT-qPCR to correlate the role of the plant hormonal changes during the tri-trophic virus-mite-plant interplay. A previous infestation with either non-viruliferous and viruliferous mites positively influenced oviposition and the number of adult individuals in the resulting populations. Mite populations were higher on branches that had received a previous mite infestation than branches that did not. There was an increase in the expression of PR4, a marker gene in the jasmonic acid (JA) pathway, in the treatment with non-viruliferous mites, indicating a response from the plant to their feeding. Conversely, an induced expression of PR1, a marker gene in the salicylic acid (SA) pathway, was observed mainly in the treatment with viruliferous mites, which suggests the activation of a plant response against the pathogen. The earlier mite infestation, as well as the presence of leprosis lesions and a gypsum mixture as artificial shelters, all fostered the growth of the B. yothersi populations after the second infestation, regardless of the presence or absence of CiLV-C. Furthermore, it is suggested that B. yothersi feeding actually induces the JA pathway in plants. At the same time, the CiLV-C represses the JA pathway and induces the SA pathway, which benefits the mite vector.

Indole-3 acetic acid negatively regulates rose black spot disease resistance through antagonizing the salicylic acid signaling pathway via jasmonic acid

MAIN CONCLUSION

IAA cooperates with JA to inhibit SA and negatively regulates rose black spot disease resistance. Black spot disease caused by the fungus Marssonina rosae is the most prevalent and severe ailment in rose cultivation, leading to the appearance of black spots on leaves and eventual leaf fall, significantly impacting the utilization of roses in gardens. Salicylic acid (SA) and jasmonic acid (JA) are pivotal hormones that collaborate with indole-3 acetic acid (IAA) in regulating plant defense responses; however, the detailed mechanisms underlying the induction of black spot disease resistance by IAA, JA, and SA remain unclear. In this study, transcript analysis was conducted on resistant (R13-54) and susceptible (R12-26) lines following M. rosae infection. In addition, the impact of exogenous interference with IAA on SA- and JA-mediated disease resistance was examined. The continuous accumulation of JA, in synergy with IAA, inhibited activation of the SA signaling pathway in the early infection stage, thereby negatively regulating the induction of effective resistance to black spot disease. IAA administration alleviated the inhibition of SA on JA to negatively regulate the resistance of susceptible strains by further enhancing the synthesis and accumulation of JA. However, IAA did not contribute to the negative regulation of black spot resistance when high levels of JA were inhibited. Virus-induced gene silencing of RcTIFY10A, an inhibitor of the JA signaling pathway, further suggested that IAA upregulation led to a decrease in disease resistance, a phenomenon not observed when the JA signal was inhibited. Collectively, these findings indicate that the IAA-mediated negative regulation of black spot disease resistance relies on activation of the JA signaling pathway.

Identification and Validation of Reference Genes for qRT-PCR Analysis of Petal-Color-Related Genes in Rosa praelucens

The flower's color is regarded as one of the most outstanding features of the rose. Rosa praelucens Byhouwer, an endemic and critically endangered decaploid wild rose species, is abundant in phenotypic diversity, especially in flower color variation, from white to different degrees of pink. The mechanism underlying this variation, e.g., the level of petal-color-related genes, is worth probing. Seven candidate reference genes for qRT-PCR analysis, including tubulin α chain (TUBA), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), histone H2B (Histone2A), eukaryotic translation elongation factor 1-α (EEF1A), 60S ribosomal protein (RPL37), eukaryotic translation initiation factor 1-α (EIF1A), and aquaporins (AQP), were detected from the transcriptome datasets of full blooming flowers of white-petaled and pink-petaled individuals, and their expression stabilities were evaluated through qRT-PCR analysis. According to stability rankings analysis, EEF1A showed the highest stability and could be chosen as the most suitable reference gene. Moreover, the reliability of EEF1A was demonstrated via qRT-PCR analysis of six petal-color-related target genes, the expression patterns of which, through EEF1A normalization, were found to be consistent with the findings of transcriptome analysis. The result provides an optimal reference gene for exploring the expression level of petal-color-related genes in R. praelucens, which will accelerate the dissection of petal-color-variation mechanisms in R. praelucens.

Comparative transcriptome and metabolome analysis revealed diversity in the response of resistant and susceptible rose (Rosa hybrida) varieties to Marssonina rosae

Rose black spot disease caused by Marssonina rosae is among the most destructive diseases that affects the outdoor cultivation and production of roses; however, the molecular mechanisms underlying the defensive response of roses to M. rosae have not been clarified. To investigate the diversity of response to M. rosae in resistant and susceptible rose varieties, we performed transcriptome and metabolome analyses of resistant (KT) and susceptible (FG) rose varieties and identified differentially expressed genes (DEGs) and differentially accumulated metabolites (DAMs) in response to M. rosae at different time points. In response to M. rosae, DEGs and DAMs were mainly upregulated compared to the control and transcription factors were concentrated in the WRKY and AP2/ERF families. Gene Ontology analysis showed that the DEGs of FG were mainly enriched in biological processes, such as the abscisic acid-activated signaling pathway, cell wall, and defense response, whereas the DEGs of KT were mainly enriched in Golgi-mediated vesicle transport processes. Kyoto Encyclopedia of Genes and Genomes analysis showed that the DEGs of both varieties were concentrated in plant-pathogen interactions, plant hormone signal transduction, and mitogen-activated protein kinase signaling pathways, with the greatest number of DEGs associated with brassinosteroid (BR) in the plant hormone signal transduction pathway. The reliability of the transcriptome results was verified by qRT-PCR. DAMs of KT were significantly enriched in the butanoate metabolism pathway, whereas DAMs of FG were significantly enriched in BR biosynthesis, glucosinolate biosynthesis, and tryptophan metabolism. Moreover, the DAMs in these pathways were significantly positively correlated with the DEGs. Disease symptoms were aggravated when FG leaves were inoculated with M. rosae after 24-epibrassinolide treatment, indicating that the response of FG to M. rosae involves the BR signaling pathway. Our results provide new insights into the molecular mechanisms underlying rose response to M. rosae and lay a theoretical foundation for formulating rose black spot prevention and control strategies and cultivating resistant varieties.

The rose INFLORESCENCE DEFICIENT IN ABSCISSION-LIKE genes, RbIDL1 and RbIDL4, regulate abscission in an ethylene-responsive manner

KEY MESSAGE

RbIDL1 and RbIDL4 are up-regulated in an ethylene-responsive manner during rose petal abscission and restored the Arabidopsis ida-2 mutant abscission defect suggesting functional conservation of the IDA pathway in rose. Abscission is an ethylene-regulated developmental process wherein plants shed unwanted organs in a controlled manner. The INFLORESCENCE DEFICIENT IN ABSCISSION family has been identified as a key regulator of abscission in Arabidopsis, encoding peptides that interact with receptor-like kinases to activate abscission. Loss of function ida mutants show abscission deficiency in Arabidopsis. Functional conservation of the IDA pathway in other plant abscission processes is a matter of interest given the discovery of these genes in several plants. We have identified four members of the INFLORESCENCE DEFICIENT IN ABSCISSION-LIKE family from the ethylene-sensitive, early-abscising fragrant rose, Rosa bourboniana. All four are conserved in sequence and possess well-defined PIP, mIDa and EPIP motifs. Three of these, RbIDL1, RbIDL2 and RbIDL4 show a three-fourfold increase in transcript levels in petal abscission zones (AZ) during ethylene-induced petal abscission as well as natural abscission. The genes are also expressed in other floral tissues but respond differently to ethylene in these tissues. RbIDL1 and RbIDL4, the more prominently expressed IDL genes in rose, can complement the abscission defect of the Arabidopsis ida-2 mutant; while, promoters of both genes can drive AZ-specific expression in an ethylene-responsive manner even in Arabidopsis silique AZs indicating recognition of AZ-specific and ethylene-responsive cis elements in their promoters by the abscission machinery of rose as well as Arabidopsis.

P Starvation in Roses Leads to Strongly Genotype-Dependent Induction of P-Transporter Genes during Black Spot Leaf Disease

Phosphorous starvation in plants has been reported to have contrasting effects on the interaction with pathogens in different plant pathogen systems and plant species. Both increases and decreases in susceptibility have been observed in numerous reports. Here, we analysed black spot infection and the leaf expression of two plant phosphate transporters and one defence marker gene in roses after phosphorous starvation. We varied three factors: phosphate starvation versus full supply of phosphorous, black spot infection vs. mock inoculation, and different susceptible and resistant progeny of a biparental rose population. Black spot susceptibility or resistance was not significantly changed upon phosphate starvation in either compatible or incompatible interactions. The expression of phosphate transporters was strongly induced upon starvation, but in some genotypes, expression was altered by black spot interaction as well. The marker for pathogenic interactions was exclusively induced by interaction with black spot, but the expression was altered by a combination of phosphate starvation and interaction with the fungus in some genotypes. In summary, phosphate starvation has clear effects on the gene expression of phosphate transporters in rose leaves, and the interaction with a hemibiotrophic leaf pathogen is strongly genotype dependent.

Selection and stability validation of reference gene candidates for transcriptional analysis in Rousettus aegyptiacus

Bats are the only mammals capable of powered flight and their body temperature can reach up to 42 °C during flight. Additionally, bats display robust type I IFN interferon (IFN-I) responses and some species constitutively express IFN-α. Reference genes with stable expression under temperature oscillations and IFN-I release are therefore critical for normalization of quantitative reverse-transcription polymerase chain reaction (qRT-PCR) data in bats. The expression stability of reference genes in Rousettus aegyptiacus remains elusive, although this species is frequently used in the infection research. We selected ACTB, EEF1A1, GAPDH and PGK1 as candidate reference genes and evaluated their expression stability in various tissues and cells from this model bat species upon IFN-I treatment at 35 °C, 37 °C and 40 °C by qRT-PCR. We employed two statistical algorithms, BestKeeper and NormFinder, and found that EEF1A1 exhibited the highest expression stability under all tested conditions. ACTB and GAPDH displayed unstable expression upon temperature change and IFN-I treatment, respectively. By normalizing to EEF1A1, we uncovered that GAPDH expression was significantly induced by IFN-I in R. aegyptiacus. Our study identifies EEF1A1 as the most suitable reference gene for qRT-PCR studies upon temperature changes and IFN-I treatment and unveils the induction of GAPDH expression by IFN-I in R. aegyptiacus. These findings are pertinent to other bat species and may be relevant for non-volant mammals that show physiological fluctuations of core body temperature.

Comparative Transcriptome and Weighted Gene Co-expression Network Analysis Identify Key Transcription Factors of Rosa chinensis 'Old Blush' After Exposure to a Gradual Drought Stress Followed by Recovery

Rose is one of the most fundamental ornamental crops, but its yield and quality are highly limited by drought. The key transcription factors (TFs) and co-expression networks during rose’s response to drought stress and recovery after drought stress are still limited. In this study, the transcriptomes of leaves of 2-year-old cutting seedlings of Rosa chinensis ‘Old Blush’ from three continuous droughted stages (30, 60, 90 days after full watering) and rewatering were analyzed using RNA sequencing. Weighted gene co-expression network analysis (WGCNA) was used to construct a co-expression network, which was associated with the physiological traits of drought response to discovering the hub TFs involved in drought response. More than 45 million high-quality clean reads were generated from the sample and used for comparison with the rose reference genome. A total of 46433 differentially expressed genes (DEGs) were identified. Gene Ontology (GO) term enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that drought stress caused significant changes in signal transduction, plant hormones including ABA, auxin, brassinosteroid (BR), cytokinin, ethylene (ET), jasmonic acid (JA) and salicylic acid (SA), primary and secondary metabolism, and a certain degree of recovery after rewatering. Gene co-expression analysis identified 18 modules, in which four modules showed a high degree of correlation with physiological traits. In addition, 42 TFs including members of NACs, WRKYs, MYBs, AP2/ERFs, ARFs, and bHLHs with high connectivity in navajowhite1 and blue modules were screened. This study provides the transcriptome sequencing report of R. chinensis ‘Old Blush’ during drought stress and rewatering process. The study also identifies the response of candidate TFs to drought stress, providing guidelines for improving the drought tolerance of the rose through molecular breeding in the future.

Prickle morphogenesis in rose is coupled with secondary metabolite accumulation and governed by canonical MBW transcriptional complex

Rose is an economically important flowering plant that holds an essential place in cut flower, medicinal, and aromatic industries. The presence of prickles, epidermal outgrowths resembling trichomes, on rose is highly undesirable as these make harvesting and transportation difficult. Attempts were made for generating rose varieties lacking prickles via breeding and natural selections; however, these approaches obtained only chimeric and genetically unstable prickle-less mutants. The alternative way to get rid of prickles is via genetic manipulations, but the molecular mechanisms of prickle initiation and development in rose are almost unexplored. Therefore, the present study was carried out to understand the morphological, molecular, and correlated metabolic changes underlining prickle morphogenesis in a prickle-bearing Rosa hybrida L. cv. "First Red (FR)". The histological and metabolomic analyses at three distinct stages of the prickle morphogenesis, namely, emerging tiny initiating prickles, partially greenish soft prickles, and brownish hard prickles, demonstrated a gradually increasing deposition of phenolic compounds and lignification with development. Corresponding RNAseq analysis revealed an upregulation of the genes involved in secondary metabolism, especially in the phenylpropanoid biosynthetic pathway. A set of genes encoding a transcriptional network similar to the one regulating epidermal cell differentiation leading to phenylpropanoid accumulation and trichome development, was also upregulated. Differential expression of this transcriptional network in prickle-less R. hybrida L. cv. "Himalayan Wonder" compared to prickly FR plants substantiated its involvement in prickle morphogenesis. The results collectively supported the proposition that prickles are evolved from trichomes and provided molecular clues towards engineering prickle-less roses.

SIGNIFICANCE STATEMENT

Prickles, the vasculature less epidermal outgrowths resembling trichomes, are defense organs protecting plants against herbivory. Despite biological significance, the mechanism of prickle morphogenesis remains obscure. Here, we show that like trichomes, prickles accumulate secondary metabolites, especially lignin and flavonoids, during morphogenesis. Cognate transcriptome analysis demonstrated that upregulation of a hormone-regulated transcriptional activation-inhibition network, known to govern trichome morphogenesis, likely triggers the differentiation of epidermal cells to outgrow into prickle.

Diversity and selection of the continuous-flowering gene, RoKSN, in rose

Blooming seasonality is an important trait in ornamental plants and was selected by humans. Wild roses flower only in spring whereas most cultivated modern roses can flower continuously. This trait is explained by a mutation of a floral repressor gene, RoKSN, a TFL1 homologue. In this work, we studied the origin, the diversity and the selection of the RoKSN gene. We analyzed 270 accessions, including wild and old cultivated Asian and European roses as well as modern roses. By sequencing the RoKSN gene, we proposed that the allele responsible for continuous-flowering, RoKSN^copia, originated from Chinese wild roses (Indicae section), with a recent insertion of the copia element. Old cultivated Asian roses with the RoKSN^copia allele were introduced in Europe, and the RoKSN^copia allele was progressively selected during the 19th and 20th centuries, leading to continuous-flowering modern roses. Furthermore, we detected a new allele, RoKSN^A181, leading to a weak reblooming. This allele encodes a functional floral repressor and is responsible for a moderate accumulation of RoKSN transcripts. A transient selection of this RoKSN^A181 allele was observed during the 19th century. Our work highlights the selection of different alleles at the RoKSN locus for recurrent blooming in rose.

Selection of suitable reference genes for quantitative real time PCR in different Tulasnella isolates and orchid-fungus symbiotic germination system

Under natural conditions, mycorrhizal symbiosis accompanies nearly the entire life cycle of orchids from seed germination through to flowering and fruiting. Tulasnella-like orchid mycorrhizal fungi are the most common mycorrhizal fungi found in association with orchid species. Presently suitable reference genes have not been systematically selected for the quantification of gene expression via Real-Time Quantitative Reverse Transcription PCR (RT-qPCR). We evaluated 12 candidate Tulasnella genes in nine different Tulasnella isolates and in the Dendrobium-fungal symbiotic germination associations followed by statistical analysis using the programs Bestkeeper, geNorm, and Normfinder to analyze the expression stability of the individual genes. The results showed that the EF2, UBC, and PP2A genes had the highest rankings with relatively stable expression levels across the different genotypes and during the symbiotic seed germination process by the three programs, and may be suitable for RT-qPCR normalization. Furthermore, the gene encoding C-5 Sterol desaturase (C5SD) was selected to verify the reliability of EF2, UBC, and PP2A expression during the Tulasnella-Dendrobium symbiotic seed germination process. This study is the first systematic exploration of optimal reference genes for gene expression studies during the colonization of orchid seeds by the mycorrhizal fungus Tulasnella.

An R2R3-MYB Transcription Factor RmMYB108 Responds to Chilling Stress of Rosa multiflora and Conferred Cold Tolerance of Arabidopsis

A sudden cooling in the early spring or late autumn negatively impacts the plant growth and development. Although a number of studies have characterized the role of the transcription factors (TFs) of plant R2R3-myeloblastosis (R2R3-MYB) in response to biotic and abiotic stress, plant growth, and primary and specific metabolisms, much less is known about their role in Rosa multiflora under chilling stress. In the present study, RmMYB108, which encodes a nuclear-localized R2R3-MYB TF with a self-activation activity, was identified based on the earlier published RNA-seq data of R. multiflora plants exposed to short-term low-temperature stress and also on the results of prediction of the gene function referring Arabidopsis. The RmMYB108 gene was induced by stress due to chilling, salt, and drought and was expressed in higher levels in the roots than in the leaves. The heterologous expression of RmMYB108 in Arabidopsis thaliana significantly enhanced the tolerance of transgenic plants to freezing, water deficit, and high salinity, enabling higher survival and growth rates, earlier flowering and silique formation, and better seed quantity and quality compared with the wild-type (WT) plants. When exposed to a continuous low-temperature stress at 4°C, transgenic Arabidopsis lines-overexpressing RmMYB108 showed higher activities of superoxide dismutase and peroxidase, lower relative conductivity, and lower malondialdehyde content than the WT. Moreover, the initial fluorescence (F _o) and maximum photosynthetic efficiency of photosystem II (F _v/F _m) changed more dramatically in the WT than in transgenic plants. Furthermore, the expression levels of cold-related genes involved in the ICE1 (Inducer of CBF expression 1)-CBFs (C-repeat binding factors)-CORs (Cold regulated genes) cascade were higher in the overexpression lines than in the WT. These results suggest that RmMYB108 was positively involved in the tolerance responses when R. multiflora was exposed to challenges against cold, freeze, salt, or drought and improved the cold tolerance of transgenic Arabidopsis by reducing plant damage and promoting plant growth.

Petal abscission in fragrant roses is associated with large scale differential regulation of the abscission zone transcriptome

Flowers of fragrant roses such as Rosa bourboniana are ethylene-sensitive and undergo rapid petal abscission while hybrid roses show reduced ethylene sensitivity and delayed abscission. To understand the molecular mechanism underlying these differences, a comparative transcriptome of petal abscission zones (AZ) of 0 h and 8 h ethylene-treated flowers from R. bourboniana was performed. Differential regulation of 3700 genes (1518 up, 2182 down) representing 8.5% of the AZ transcriptome was observed between 0 and 8 h ethylene-treated R. bourboniana petal AZ. Abscission was associated with large scale up-regulation of the ethylene pathway but prominent suppression of the JA, auxin and light-regulated pathways. Regulatory genes encoding kinases/phosphatases/F-box proteins and transcription factors formed the major group undergoing differential regulation besides genes for transporters, wall modification, defense and phenylpropanoid pathways. Further comparisons with ethylene-treated petals of R. bourboniana and 8 h ethylene-treated AZ (R. hybrida) identified a core set of 255 genes uniquely regulated by ethylene in R. bourboniana AZ. Almost 23% of these encoded regulatory proteins largely conserved with Arabidopsis AZ components. Most of these were up-regulated while an entire set of photosystem genes was prominently down-regulated. The studies provide important information on regulation of petal abscission in roses.

Genome-wide identification of WD40 genes reveals a functional diversification of COP1-like genes in Rosaceae

Genome-wide identification of WD40-like genes reveals a duplication of COP1-like genes, one of the key players involved in regulation of flowering time and photomorphogenesis, with strong functional diversification in Rosaceae. WD40 proteins play crucial roles in a broad spectrum of developmental and physiological processes. Here, we conducted a systematic characterization of this family of genes in Rosa chinensis 'Old Blush' (OB), a founder genotype for modern rose domestication. We identified 187 rose WD40 genes and classified them into 5 clusters and 15 subfamilies with 11 of RcWD40s presumably generated via tandem duplication. We found RcWD40 genes were expressed differentially following stages of vegetative and reproductive development. We detected a duplication of CONSTITUTIVE PHOTOMORPHOGENIC1-like genes in rose (RcCOP1 and RcCOP1L) and other Rosaceae plants. Featuring a distinct expression pattern and a different profile of cis-regulatory-elements in the transcriptional regulatory regions, RcCOP1 seemed being evolutionarily conserved while RcCOP1L did not dimerize with RcHY5 and RcSPA4. Our data thus reveals a functional diversification of COP1-like genes in Rosacaeae plants, and provides a valuable resource to explore the potential function and evolution of WD40-like genes in Rosaceae plants.

Identification of candidate reference genes for quantitative RT-PCR in Miscanthus sinensis subjected to various abiotic stresses

Quantitative real-time PCR (qRT-PCR) has been widely used for studying gene expression at the transcript level. Its accuracy usually relies on the reference genes that are utilized for data normalization. Miscanthus sinensis, a perennial C4 grass with high biomass and strong resistance to adversities, is often utilized as a high value energy crop. However, no reliable reference genes have been investigated for normalizing gene expression for this species. In this study, 12 candidate reference genes were selected to identify their stability under five different abiotic stress treatments (drought, salt, cadmium, chromium and arsenic) by using geNorm, NormFinder, BestKeeper and RefFinder softwares. The results showed that 18S rRNA and Unigene33312 were the best reference genes under drought treatments. Unigene33312 and Unigene33024 were found to be the most stably expressed genes under salt stress and Cd stress. Moreover, Unigene33024 and PP2A were the most suitable reference genes under Cr stress and Unigene33024 and Sb09g019750 were deemed more suitable reference genes under As stress. In total, considering all the samples, Unigene33024 and PP2A were the most stable genes while ACTIN and Unigene26576 were the least stable reference genes for internal control. The expression patterns of two target genes (Cu/Zn SOD and CAT) were used to further verify those selected reference genes under different conditions. The results showed that the most and the least stable reference genes had clearly different expression patterns. This work comprehensively estimated the stability of reference genes in M. sinensis which may give insight to the reference genes selection in other tissues as well as other related varieties. These suggested reference genes would assist in further putative gene expression validation in M. sinensis.

Genome-wide identification and functional analysis of JmjC domain-containing genes in flower development of Rosa chinensis

We genome-wide identified 28 JmjC domain-containing genes, further spatio-temporal expression profiling and genetic analysis defined them as epigenetic regulators in flowering initiation of Rosa chinensis. The JmjC domain-containing histone demethylases play critical roles in maintaining homeostasis of histone methylations, thus are vital for plant growth and development. Genome-wide identification of the JmjC domain-containing genes have been reported in several species, however, no systematic study has been performed in rose plants. In this paper, we identified 28 JmjC domain-containing genes from the newly published genome database of Rosa chinensis. The JmjC domain-containing proteins in R. chinensis were divided into seven groups, KDM3 was the largest group with 13 members, and JmjC domain-only A and KDM5B were the smallest clades both with only one member. Although all the JmjC domain proteins having a conserved JmjC domain, the gene and protein structure experienced differentiation and specification during the evolution, especially in KDM3 clade, one gene (RcJMJ40) was found carrying site deletions for cofactors Fe (II) and α-KG binding which were crucial for demethylase activities, three genes (RcJMJ41, RcJMJ43 and RcJMJ44) had no intron while two of them had tandem JmjC domains. Spatial expression pattern analysis of these JmjC domain-containing genes in different tissues showed most of them were highly expressed in reproductive tissues such as floral meristem and closed flowers than vegetative tissues, demonstrating their important functions in developmental switch from vegetative to reproductive growth of roses. Temporal expression profiling indicated majority of JmjC domain-containing genes from R. chinensis fluctuated along with floral bud differentiation and development, further proving their essential roles in flower organogenesis. VIGS induced silencing of RcJMJ12 led to delayed flowering time, and decreased the expression levels of flowering integrator such as RcFT, RcSOC1, RcFUL, RcLFY and RcAP1, therefore providing the genetic evidence of RcJMJ12 in flowering initiation. Collectively, spatio-temporal expression profiling and genetic analysis defined the JmjC domain-containing genes as important epigenetic regulators in flower development of R. chinensis.

Selection and validation of reference genes for gene expression studies in Pseudomonas brassicacearum GS20 using real-time quantitative reverse transcription PCR

Pseudomonas brassicacearum GS20 is an antagonistic strain of bacteria recently isolated from the rhizosphere of Codonopsis pilosula. No validated reference gene has been indentified from P. brassicacearum to use in the normalization of real-time quantitative reverse transcription-PCR data. Therefore, in this study, nine candidate reference genes (recA, gyrA, rpoD, proC, gmk, rho, 16S, ftsz, and secA) were assessed at different growth phases and under various growth conditions. The expression stability of these candidate genes was evaluated using BestKeeper, NormFinder and GeNorm. In general, the results showed rho, rpoD and gyrA were the most suitable reference genes for P. brassicacearum GS20. The relative expression of iron-regulated gene (fhu) was normalized to verify the reliability of the proposed reference genes under iron-replete and iron-limited conditions. The trend in relative expression was consistent with the change in siderophore production under different iron conditions. This study presents reliable reference genes for transcriptional studies in P. brassicacearum GS20 under the chosen experimental conditions.

Higher Ramie mosaic virus transmission efficiency by females than by males of Bemisia tabaci MED

Begomoviruses can modify their transmission vector, Bemisia tabaci, to benefit their spread, although this may not always be the case. Here, the new begomovirus Ramie mosaic virus (RaMoV) and its vector B. tabaci MED, which is dominant in China and many regions of the world, were used as a model to examine direct and indirect interaction and virus transmission by B. tabaci MED of different sexes. No significant direct or indirect effects of RaMoV were observed in B. tabaci MED females, although RaMoV could shorten the life span of B. tabaci MED females by up to 4 days. A test of RaMoV transmission by different sexes of B. tabaci MED showed that there was higher virus transmission efficiency by females than males. Overall, RaMoV is transmitted by B. tabaci MED in a sex-dependent manner, and further research is needed to uncover the mechanism of the difference in RaMoV transmission by different sexes of B. tabaci.

Selection of Suitable Reference Genes in Pinus massoniana Lamb. Under Different Abiotic Stresses for qPCR Normalization

The normalization of data by choosing suitable reference genes is fundamental for obtaining accurate and reliable results in quantitative real-time polymerase chain reaction (qPCR) analyses. In this study, the expression stability of 12 candidate reference genes of Pinus massoniana under different abiotic stresses was evaluated using four statistical algorithms: geNorm, NormFinder, BestKeeper, and RefFinder. The results indicate that the following genes could be used as reference genes under different treatments: Actin 2 (ACT2) and F-box family gene (F-box) for salinity treatment, cyclophilin (CYP) and alpha-tubulin (TUA) for ABA treatment, actin 7 (ACT7) and CYP for drought treatment, actin 1 (ACT1) and ACT7 for cold treatment, ACT1 and CYP for heat treatment, and TUA and ACT2 for the “Total” group. To validate the suitability of the selected reference genes in this study, the Short-Root protein (SHR), Alpha-pinene synthase (APS), and Pyrabactin resistance-like protein (PYL) gene expression patterns were analyzed. The expression patterns had significant biases when the most unstable reference genes were used for normalization, compared with when the optimum reference gene or gene combinations were used for normalization. These results will be beneficial for further studies on gene transcription in early-stage, unlignified seedlings of P. massoniana.

Evaluation of reference genes for real-time quantitative PCR analysis in southern corn rootworm, Diabrotica undecimpunctata howardi (Barber)

Quantitative reverse transcription PCR (RT-qPCR) is one of the most efficient, reliable and widely used techniques to quantify gene expression. In this study, we evaluated the performance of six southern corn rootworm, Diabrotica undecimpunctata howardi (Barber), housekeeping genes (HKG), β-actin (Actin), β-tubulin (Tubulin), elongation factor 1 alpha (EF1α), glyceraldehyde-3 phosphate dehydrogenase (GAPDH), 40 S ribosomal protein S9 (RpS9) and ubiquitin-conjugating protein (Ubi), under different experimental conditions such as developmental stage, exposure of neonate and adults to dsRNA, exposure of adults to different temperatures, different 3^rd instar larva tissues, and neonate starvation. The HKGs were analyzed with four algorithms, including geNorm, NormFinder, BestKeeper, and delta-CT. Although the six HKGs showed a relatively stable expression pattern among different treatments, some variability was observed. Among the six genes, EF1α exhibited the lowest Ct values for all treatments while Ubi exhibited the highest. Among life stages and across treatments, Ubi exhibited the least stable expression pattern. GAPDH, Actin, and EF1α were among the most stable HKGs in the majority of the treatments. This research provides HKG for accurate normalization of RT-qPCR data in the southern corn rootworm. Furthermore, this information can contribute to future genomic and functional genomic research in Diabrotica species.

Selection of reference genes for flowering pathway analysis in the masting plants, Celmisia lyallii and Chionochloa pallens, under variable environmental conditions

Mast flowering is characterised by mass synchronised flowering at irregular intervals over a wide geographical area. An understanding of the molecular drivers of mast flowering requires expression analysis of key developmentally regulated gene(s). Reverse transcription-quantitative PCR is the gold standard technique used to assess expression of target gene(s) and to validate high-throughput sequencing data. Selection and validation of appropriate reference gene(s), used as normalisation factors in transcript abundance analysis, is an essential step to avoid ambiguous expression results. Eight candidate reference genes were assessed to select the best internal normalisation factors in naturally growing masting plants Chionochloa pallens and Celmisia lyallii. Statistical packages geNorm, Normfinder, BestKeeper, ΔCt and RefFinder were used to determine the expression stability in plants translocated to different altitudes and sampled across the season. GAPDH and PP2a in Celmisia and ExP and THP in Chionochloa were found to be the best pairs of reference genes for normalisation of the gene expression data. Our study revealed environmentally-induced changes in reference gene expression, information that will be utilised as we investigate flowering phenology of masting plants under global climatic change.

Selection and validation of reference genes of Paeonia lactiflora in growth development and light stress

The stem of Paeonia lactiflora will bend when it grows in greenhouse at a low light intensity. It is important to explore causes of morphological changes of peony to improve its quality. Gene expression can be evaluated by quantitative real-time PCR, based on reference gene. However, systematic selection of reference genes under weak lighting for herbaceous peony is lacking. To address this problem, we first selected 10 candidate reference genes based on a coefficient of variation of gene expression from peony stem transcriptome data. Then, geNorm, NormFinder and BestKeeper were applied to assess the stability of the genes, and RankAggreg was used to give a comprehensive ranking. The results show that there are some differences in optimal reference genes among samples from different organs and under the two lighting conditions, and the optimal number of suitable reference genes is distinct. Two selected suitable reference genes were then used to normalize target genes, and the results were compared with transcriptome data. Consistent gene expression trends were obtained, indicating the reliability of the method. To the best of our knowledge, this is the first time reference genes for herbaceous peony were selected in different organs, developmental stages and under two kinds of lighting conditions. The findings can provide a practical method for selecting reference genes for peony under these conditions and demonstrate a useful combination of reference genes.

Differential and reciprocal regulation of ethylene pathway genes regulates petal abscission in fragrant and non-fragrant roses

The fragrant rose, Rosa bourboniana, is highly sensitive to ethylene and shows rapid petal abscission (within 16-18 h) while the non-fragrant hybrid rose, R. hybrida, shows delayed abscission (50-52 h) due to reduced ethylene sensitivity. To understand the molecular basis governing these differences, all components of the ethylene pathway (biosynthesis/ receptor/signalling) were studied for expression during abscission. Transcript accumulation of most ethylene biosynthesis genes (ACS/ACO families) increased rapidly in petal abscission zones of R. bourboniana within 4-8 h of ethylene treatment. The expression of most receptor and signalling genes encoding CTRs, EIN2 and EIN3/EIL homologues also followed similar kinetics. Under natural field conditions where abscission takes longer, there was a temporal delay in transcript accumulation of most ethylene pathway genes while some biosynthesis genes (showing reduced ethylene sensitivity) were more strongly up-regulated by abscission cues. In contrast, in R. hybrida where even ethylene-induced abscission is considerably delayed, transcript accumulation of most ethylene biosynthesis and signalling genes was, surprisingly, reduced by ethylene and showed an opposite regulation compared to R. bourboniana. The results suggest that differential and reciprocal regulation of ethylene pathway is one of the major reasons for differences in petal abscission and vase-life between Rosa bourboniana and R. hybrida.

Interaction of roses with a biotrophic and a hemibiotrophic leaf pathogen leads to differences in defense transcriptome activation

Transcriptomic analysis resulted in the upregulation of the genes related to common defense mechanisms for black spot and the downregulation of the genes related to photosynthesis and cell wall modification for powdery mildew. Plant pathogenic fungi successfully colonize their hosts by manipulating the host defense mechanisms, which is accompanied by major transcriptome changes in the host. To characterize compatible plant pathogen interactions at early stages of infection by the obligate biotrophic fungus Podosphaera pannosa, which causes powdery mildew, and the hemibiotrophic fungus Diplocarpon rosae, which causes black spot, we analyzed changes in the leaf transcriptome after the inoculation of detached rose leaves with each pathogen. In addition, we analyzed differences in the transcriptomic changes inflicted by both pathogens as a first step to characterize specific infection strategies. Transcriptomic changes were analyzed using next-generation sequencing based on the massive analysis of cDNA ends approach, which was validated using high-throughput qPCR. We identified a large number of differentially regulated genes. A common set of the differentially regulated genes comprised of pathogenesis-related (PR) genes, such as of PR10 homologs, chitinases and defense-related transcription factors, such as various WRKY genes, indicating a conserved but insufficient PTI [pathogen associated molecular pattern (PAMP) triggered immunity] reaction. Surprisingly, most of the differentially regulated genes were specific to the interactions with either P. pannosa or D. rosae. Specific regulation in response to D. rosae was detected for genes from the phenylpropanoid and flavonoid pathways and for individual PR genes, such as paralogs of PR1 and PR5, and other factors of the salicylic acid signaling pathway. Differently, inoculation with P. pannosa leads in addition to the general pathogen response to a downregulation of genes related to photosynthesis and cell wall modification.

Selection of suitable reference genes for quantitative real-time PCR gene expression analysis in Mulberry (Morus alba L.) under different abiotic stresses

Mulberry (Morus alba L.) is the sole food source for the mulberry silkworm, Bombyx mori and therefore important for sericulture industry. Different abiotic stress conditions like drought, salt, heat and cold stress adversely affect the productivity and quality of mulberry leaves. Quantitative real time PCR (qPCR) is a reliable and widely used method to identify abiotic stress responsive genes and molecular mechanism in different plant species. Selection of suitable reference genes is important requirement for normalizing the expression of genes through qRT-PCR study. In the present study, we have selected eight candidate reference genes in mulberry for analyzing their expression stability in different abiotic stress treatments including drought, salt, heat and cold stresses. The expression stability of these reference genes was determined using geNorm, NormFinder and RefFinder statistical algorithms. The results showed that Ubiquitin and protein phosphatase 2A regulatory subunit A (PP2A) were the most stable genes across all the treatment samples. However, analysis of individual stresses revealed different expression profiles and stability of reference genes. Actin3 and PP2A were most stable in drought and salt conditions respectively. RPL3 most preferred in heat stress and Ubiquitin was most stable in cold stress. We propose the ubiquitin and PP2A are the preferred reference genes for normalization of gene expression data from abiotic stresses. In addition, Actin3 are preferred for drought stress, PP2A for salt stress, RPL3 for heat stress and Ubiquitin for cold stress studies.

Selection of suitable reference genes for qRT-PCR normalisation under different experimental conditions in Eucommia ulmoides Oliv

Normalisation of data, by choosing the appropriate reference genes, is fundamental for obtaining reliable results in quantitative real-time PCR (qPCR). This study evaluated the expression stability of 11 candidate reference genes with different varieties, developmental periods, tissues, and abiotic stresses by using four statistical algorithms: geNorm, NormFinder, BestKeeper, and RefFinder. The results indicated that ubiquitin-conjugating enzyme S (UBC) and ubiquitin-conjugating enzyme E2 (UBC E2) could be used as reference genes for different E. ulmoides varieties and tissues, UBC and histone H4 (HIS4) for different developmental periods, beta-tubulin (TUB) and UBC for cold treatment, ubiquitin extension protein (UBA80) and HIS4 for drought treatment, and ubiquitin-60S ribosomal protein L40 (UBA52) and UBC E2 for salinity treatment. UBC and UBC E2 for the group “Natural growth” and “Total”, UBA80 and UBC for the group “Abiotic stresses”. To validate the suitability of the selected reference genes in this study, mevalonate kinase (MK), phenylalanine ammonia-lyase (PAL), and 4-coumarate-CoA ligase (4CL) gene expression patterns were analysed. When the most unstable reference genes were used for normalisation, the expression patterns had significant biases compared with the optimum reference gene combinations. These results will be beneficial for more accurate quantification of gene expression levels in E. ulmoides.

Identification of suitable internal control genes for transcriptional studies in Eleusine coracana under different abiotic stress conditions

Finger millet [Eleusine coracana (L.) Gaertn] is an excellent food and forage crop of arid and semiarid areas in Africa and Asia. It is well adapted to drought, heat, high salinity, poor soil fertility and low pH with an efficient C₄ carbon fixation mechanism for high yield potential. To normalize the target gene expression data, the identification of suitable reference genes is essential. Ten candidate reference genes were selected and their expression stability was analyzed in various samples treated with different abiotic stress conditions. Five different statistical algorithms: geNorm, NormFinder, BestKeeper, ΔCt, and RefFinder were used to determine the stability of these genes. Our results revealed GAPDH, EEF1a, ACT and CYC as highly stable reference genes and PP2A and eIF4A as least stable reference genes across all the samples and suggesting that these genes could be used for accurate transcript normalization under abiotic stress. To the best of our knowledge, this is the first report on identification of suitable reference genes for accurate transcript normalization using qRT-PCR in finger millet under abiotic stress.

The TNL gene Rdr1 confers broad-spectrum resistance to Diplocarpon rosae

Black spot disease, which is caused by the ascomycete Diplocarpon rosae, is the most severe disease in field-grown roses in temperate regions and has been distributed worldwide, probably together with commercial cultivars. Here, we present data indicating that muRdr1A is the active Rdr1 gene, a single-dominant TIR-NBS-LRR (Toll/interleukin-1 receptor-nucleotide binding site-leucine rich repeat) (TNL)-type resistance gene against black spot disease, which acts against a broad range of pathogenic isolates independent of the genetic background of the host genotype. Molecular analyses revealed that, compared with the original donor genotype, the multiple integrations that are found in the primary transgenic clone segregate into different integration patterns in its sexual progeny and do not show any sign of overexpression. Rdr1 provides resistance to 13 different single-spore isolates belonging to six different races and broad field mixtures of conidia; thus far, Rdr1 is only overcome by two races. The expression of muRdr1A, the active Rdr1 gene, leads to interaction patterns that are identical in the transgenic clones and the non-transgenic original donor genotype. This finding indicates that the interacting avirulence (Avr) factor on the pathogen side must be widespread among the pathogen populations and may have a central function in the rose-black spot interaction. Therefore, the Rdr1 gene, pyramided with only a few other R genes by sexual crosses, might be useful for breeding roses that are resistant to black spot because the spread of new pathogenic races of the fungus appears to be slow.

Reference gene selection for qRT-PCR analysis of flower development in Lagerstroemia indica and L. speciosa

Quantitative real-time polymerase chain reaction (qRT-PCR) is a prevalent method for gene expression analysis, depending on the stability of the reference genes for data normalization. Lagerstroemia indica and L. speciosa are popular ornamental plants which are famous for the long flowering period. However, no systematic studies on reference genes in Lagerstroemia have yet been conducted. In the present study, we selected nine candidate reference genes (GAPDH, TUA, TUB, 18S, RPII, EF-1α, ATC, EIF5A and CYP) and evaluated their expression stability in different tissues during floral development of L. indica and L. speciosa using four algorithms (geNorm, NormFinder, BestKeeper and, RefFinder). Results showed that RPII and EF-1α were the most stably expressed and suitable reference genes for both of Lagerstroemia species. Moreover, ACT exhibited high expression stability in L. indica and GAPDH was a suitable reference gene for L. speciosa in different flower development stages. TUB was an unsuitable reference gene for gene expression normalization due to significant variations in expression across all samples. Finally, we verified the reliability of the selected candidate reference genes by amplifying an AGAMOUS homolog (LsAG1) of Arabidopsis thaliana. This study provides a list of suitable reference genes, thereby broadening the genetic basis of the gene expression patterns in Lagerstroemia species.

microRNAs as reference genes for quantitative PCR in cotton

Cotton (Gossypium hirsutum) is the most important non-food plant in the world. Studies concerning the fiber quality and plant fitness of cotton at molecular level depend on high sensitive and reproducible gene-expression assays. However, only a few reports have described genes suitable for normalizing gene expression data. In this study, we report for the first time that microRNAs (miRNAs) are reliable reference genes (RGs) for cotton gene expression normalization in quantitative real-time reverse transcription (RT)-PCR. The stability of cotton miRNAs was assayed in root, stem, leaf and flower samples from three different cultivars [FiberMax (FM966), Delta Opal (DO) and Cedro] and under conditions of biotic stress caused by infection with Cotton leafroll dwarf virus (CLRDV). The stability of mRNAs already described as reference genes in cotton was also assessed. The geNorm, NormFinder, BestKeeper and ΔCt algorithms were used to select the best reference genes. In 8 of the 12 sets tested, miRNAs (miR172, 168 and 390) were found to be the best RGs. To validate the best selected RGs, miR159, miR164, miR2118, miR2910, miR3476, GhDCL2 and GhDCL4 expression levels were evaluated under biotic stress conditions, and miR164 and a putative myo-inositol oxigenase gene (GhMIOX) were assessed in leaves and flowers. The RGs selected in this work proved to be excellent reference genes in the two cases studied. Our results support the use of miRNAs as reference genes for miRNA and protein-coding genes.

Reference genes for quantitative real-time PCR analysis in symbiont Entomomyces delphacidicola of Nilaparvata lugens (Stål)

Nilaparvata lugens (Stål) (Hemiptera: Delphacidae) is a major rice pest that harbors an endosymbiont ascomycete fungus, Entomomyces delphacidicola str. NLU (also known as yeast-like symbiont, YLS). Driving by demand of novel population management tactics (e.g. RNAi), the importance of YLS has been studied and revealed, which greatly boosts the interest of molecular level studies related to YLS. The current study focuses on reference genes for RT-qPCR studies related to YLS. Eight previously unreported YLS genes were cloned, and their expressions were evaluated for N. lugens samples of different developmental stages and sexes, and under different nutritional conditions and temperatures. Expression stabilities were analyzed by BestKeeper, geNorm, NormFinder, ΔCt method and RefFinder. Furthermore, the selected reference genes for RT-qPCR of YLS genes were validated using targeted YLS genes that respond to different nutritional conditions (amino acid deprivation) and RNAi. The results suggest that ylsRPS15p/ylsACT are the most suitable reference genes for temporal gene expression profiling, while ylsTUB/ylsACT and ylsRPS15e/ylsGADPH are the most suitable reference gene choices for evaluating nutrition and temperature effects. Validation studies demonstrated the advantage of using endogenous YLS reference genes for YLS studies.

An efficient transient expression system for gene function analysis in rose

BACKGROUND

Roses are widely used as garden ornamental plants and cut flowers. Rosa chinensis cv 'Old Blush' has been used as a model genotype in rose studies due to its contribution to recurrent flowering and tea scent traits of modern roses. The deficiency of efficient genetic transformation systems is a handicap limiting functional genetics studies of roses. Agrobacterium-mediated transient transformation offers a powerful tool for the characterization of gene function in plants.

RESULTS

A convenient and highly efficient Agrobacterium mediated genetic transformation protocol using R. chinensis cv 'Old Blush' seedlings in vitro as an expression system is described in this paper. The most important factor affecting transformation efficiency in this system is seedling age; 3/4-week-old rose shoots with or without roots from sub-culturing are optimal for transformation, requiring no complicated inoculation media, supplements, or carefully tuned plant growth conditions. This transient expression system was successfully applied to analysis of the gene promoter activities, DNA binding capacity of transcription factors, protein-protein interaction in physiological contexts using luciferase as a reporter gene.

CONCLUSIONS

This transient transformation system was validated as a robust and efficient platform, thus providing a new option for gene function and signaling pathway investigation in roses and further extending the utility of R. chinensis cv 'Old Blush' as a model plant to study diverse gene function and signaling pathways in Rosaceae.

Validation of Reference Genes for Quantitative Real-Time PCR during Bicolor Tepal Development in Asiatic Hybrid Lilies (Lilium spp.)

Quantitative real-time PCR (qRT-PCR) is a reliable and high-throughput technique for gene expression studies, but its accuracy depends on the expression stability of reference genes. To date, several reliable reference gene identifications have been reported in Lilium spp., but none has been obtained for lily tepals at different developmental stages. In this study, ten candidate reference genes were selected and evaluated for their expression stability in Lilium 'Tiny Padhye' during the process of bicolor tepal development. The expression stability of these candidates was evaluated by three software programs (geNorm, NormFinder, and BestKeeper) and the comparative ΔCt method, and comprehensive stability rankings were generated by RefFinder. As a result, TIP41-like family gene (TIP41) and actin (ACT) were the best combination of reference genes for tepals at different developmental stages; TIP41 and F-box family gene (F-box) for tepals under shading treatment; ACT, actin11 (ACT11), and elongation factor 1-α (EF1-α) for different tissues; and ACT, TIP41, and ACT11 for all samples. The selected optimal reference genes were further verified by analyzing the expression levels of flavonoid 3'-hydroxylase (LhF3'H) and anthocyanidin 3-O-glucosyltransfersae (LhUFGT) in tepals at different developmental stages. This study provides useful information for gene expression characterization in lilies under different experimental conditions, and can serve as a basis for similar research in other closely related species.

Roles of Aquaporins in Setaria viridis Stem Development and Sugar Storage

Setaria viridis is a C4 grass used as a model for bioenergy feedstocks. The elongating internodes in developing S. viridis stems grow from an intercalary meristem at the base, and progress acropetally toward fully expanded cells that store sugar. During stem development and maturation, water flow is a driver of cell expansion and sugar delivery. As aquaporin proteins are implicated in regulating water flow, we analyzed elongating and mature internode transcriptomes to identify putative aquaporin encoding genes that had particularly high transcript levels during the distinct stages of internode cell expansion and maturation. We observed that SvPIP2;1 was highly expressed in internode regions undergoing cell expansion, and SvNIP2;2 was highly expressed in mature sugar accumulating regions. Gene co-expression analysis revealed SvNIP2;2 expression was highly correlated with the expression of five putative sugar transporters expressed in the S. viridis internode. To explore the function of the proteins encoded by SvPIP2;1 and SvNIP2;2, we expressed them in Xenopus laevis oocytes and tested their permeability to water. SvPIP2;1 and SvNIP2;2 functioned as water channels in X. laevis oocytes and their permeability was gated by pH. Our results indicate that SvPIP2;1 may function as a water channel in developing stems undergoing cell expansion and SvNIP2;2 is a candidate for retrieving water and possibly a yet to be determined solute from mature internodes. Future research will investigate whether changing the function of these proteins influences stem growth and sugar yield in S. viridis.

Identification of Reference Genes for Quantitative Gene Expression Studies in a Non-Model Tree Pistachio (Pistacia vera L.)

The tree species, Pistacia vera (P. vera) is an important commercial product that is salt-tolerant and long-lived, with a possible lifespan of over one thousand years. Gene expression analysis is an efficient method to explore the possible regulatory mechanisms underlying these characteristics. Therefore, having the most suitable set of reference genes is required for transcript level normalization under different conditions in P. vera. In the present study, we selected eight widely used reference genes, ACT, EF1α, α-TUB, β-TUB, GAPDH, CYP2, UBQ10, and 18S rRNA. Using qRT-PCR their expression was assessed in 54 different samples of three cultivars of P. vera. The samples were collected from different organs under various abiotic treatments (cold, drought, and salt) across three time points. Several statistical programs (geNorm, NormFinder, and BestKeeper) were applied to estimate the expression stability of candidate reference genes. Results obtained from the statistical analysis were then exposed to Rank aggregation package to generate a consensus gene rank. Based on our results, EF1α was found to be the superior reference gene in all samples under all abiotic treatments. In addition to EF1α, ACT and β-TUB were the second best reference genes for gene expression analysis in leaf and root. We recommended β-TUB as the second most stable gene for samples under the cold and drought treatments, while ACT holds the same position in samples analyzed under salt treatment. This report will benefit future research on the expression profiling of P. vera and other members of the Anacardiaceae family.

RPL13A and EEF1A1 Are Suitable Reference Genes for qPCR during Adipocyte Differentiation of Vascular Stromal Cells from Patients with Different BMI and HOMA-IR

Real-time or quantitative PCR (qPCR) is a useful technique that requires reliable reference genes for data normalization in gene expression analysis. Adipogenesis is among the biological processes suitable for this technique. The selection of adequate reference genes is essential for qPCR gene expression analysis of human Vascular Stromal Cells (hVSCs) during their differentiation into adipocytes. To the best of our knowledge, there are no studies validating reference genes for the analyses of visceral and subcutaneous adipose tissue hVSCs from subjects with different Body Mass Index (BMI) and Homeostatic Model Assessment of Insulin Resistance (HOMA-IR) index. The present study was undertaken to analyze this question. We first analyzed the stability of expression of five potential reference genes: CYC, GAPDH, RPL13A, EEF1A1, and 18S ribosomal RNA, during in vitro adipogenic differentiation, in samples from these types of patients. The expression of RPL13A and EEF1A1 was not affected by differentiation, thus being these genes the most stable candidates, while CYC, GAPDH, and 18S were not suitable for this sort of analysis. This work highlights that RPL13A and EEF1A1 are good candidates as reference genes for qPCR analysis of hVSCs differentiation into adipocytes from subjects with different BMI and HOMA-IR.

Identification and Validation of Reference Genes for qRT-PCR Studies of Gene Expression in Dioscorea opposita

Quantitative real-time polymerase chain reaction (qRT-PCR) is one of the most common methods for gene expression studies. Data normalization based on reference genes is essential for obtaining reliable results for qRT-PCR assays. This study evaluated potential reference genes of Chinese yam (Dioscorea opposita Thunb.), which is an important tuber crop and medicinal plant in East Asia. The expression of ten candidate reference genes across 20 samples from different organs and development stages was assessed. We identified the most stable genes for qRT-PCR studies using combined samples from different organs. Our results also suggest that different suitable reference genes or combinations of reference genes for normalization should be applied according to different organs and developmental stages. To validate the suitability of the reference genes, we evaluated the relative expression of PE2.1 and PE53, which are two genes that may be associated with microtuber formation. Our results provide the foundation for reference gene(s) selection in D. opposita and will contribute toward more accurate gene analysis studies of the genus Dioscorea.

Selection and Validation of Reference Genes for Quantitative Real-time PCR in Gentiana macrophylla

Real time quantitative PCR (RT-qPCR or qPCR) has been extensively applied for analyzing gene expression because of its accuracy, sensitivity, and high throughput. However, the unsuitable choice of reference gene(s) can lead to a misinterpretation of results. We evaluated the stability of 10 candidates - five traditional housekeeping genes (UBC21, GAPC2, EF-1α4, UBQ10, and UBC10) and five novel genes (SAND1, FBOX, PTB1, ARP, and Expressed1) - using the transcriptome data of Gentiana macrophylla. Common statistical algorithms ΔC t, GeNorm, NormFinder, and BestKeeper were run with samples collected from plants under various experimental conditions. For normalizing expression levels from tissues at different developmental stages, GAPC2 and UBC21 had the highest rankings. Both SAND1 and GAPC2 proved to be the optimal reference genes for roots from plants exposed to abiotic stresses while EF-1α4 and SAND1 were optimal when examining expression data from the leaves of stressed plants. Based on a comprehensive ranking of stability under different experimental conditions, we recommend that SAND1 and EF-1α4 are the most suitable overall. In this study, to find a suitable reference gene and its real-time PCR assay for G. macrophylla DNA content quantification, we evaluated three target genes including WRKY30, G10H, and SLS, through qualitative and absolute quantitative PCR with leaves under elicitors stressed experimental conditions. Arbitrary use of reference genes without previous evaluation can lead to a misinterpretation of the data. Our results will benefit future research on the expression of genes related to secoiridoid biosynthesis in this species under different experimental conditions.

Selection of Reliable Reference Genes for Gene Expression Analysis under Abiotic Stresses in the Desert Biomass Willow, Salix psammophila

Salix psammophila is a desert shrub willow that has extraordinary adaptation to abiotic stresses and plays an important role in maintaining local ecosystems. Moreover, S. psammophila is regarded as a promising biomass feedstock because of its high biomass yields and short rotation coppice cycle. However, few suitable reference genes (RGs) for quantitative real-time polymerase chain reaction (qRT-PCR) constrain the study on normalization of gene expression in S. psammophila until now. Here, we investigated the expression stabilities of 14 candidate RGs across tissue types and under four abiotic stress treatments, including heat, cold, salt, and drought treatments. After calculation of PCR efficiencies, three different software, NormFinder, geNorm, and BestKeeper were employed to analyze systematically the qRT-PCR data, and the outputs were merged by RankAggreg software. The optimal RGs selected for gene expression analysis were EF1α (Elongation factor-1 alpha) and OTU (OTU-like cysteine protease family protein) for different tissue types, UBC (Ubiquitin-conjugating enzyme E2) and LTA4H (Leukotriene A-4 hydrolase homolog) for heat treatment, HIS (Histone superfamily protein H3) and ARF2 (ADP-ribosylation factor 2) for cold treatment, OTU and ACT7 (Actin 7) for salt treatment, UBC and LTA4H for drought treatment. The expression of UBC, ARF2, and VHAC (V-type proton ATPase subunit C) varied the least across tissue types and under abiotic stresses. Furthermore, the relative genes expression profiles of one tissue-specific gene WOX1a (WUSCHEL-related homeobox 1a), and four stress-inducible genes, including Hsf-A2 (Heat shock transcription factors A2), CBF3 (C-repeat binding factor 3), HKT1 (High-Affinity K⁺ Transporter 1), and GST (Glutathione S-transferase), were conducted to confirm the validity of the RGs in this study. These results provided an important RGs application guideline for gene expression characterization in S. psammophila.

Genome-Wide Identification and Validation of Reference Genes in Infected Tomato Leaves for Quantitative RT-PCR Analyses

The Gram-negative bacterium Xanthomonas campestris pv. vesicatoria (Xcv) causes bacterial spot disease of pepper and tomato by direct translocation of type III effector proteins into the plant cell cytosol. Once in the plant cell the effectors interfere with host cell processes and manipulate the plant transcriptome. Quantitative RT-PCR (qRT-PCR) is usually the method of choice to analyze transcriptional changes of selected plant genes. Reliable results depend, however, on measuring stably expressed reference genes that serve as internal normalization controls. We identified the most stably expressed tomato genes based on microarray analyses of Xcv-infected tomato leaves and evaluated the reliability of 11 genes for qRT-PCR studies in comparison to four traditionally employed reference genes. Three different statistical algorithms, geNorm, NormFinder and BestKeeper, concordantly determined the superiority of the newly identified reference genes. The most suitable reference genes encode proteins with homology to PHD finger family proteins and the U6 snRNA-associated protein LSm7. In addition, we identified pepper orthologs and validated several genes as reliable normalization controls for qRT-PCR analysis of Xcv-infected pepper plants. The newly identified reference genes will be beneficial for future qRT-PCR studies of the Xcv-tomato and Xcv-pepper pathosystems, as well as for the identification of suitable normalization controls for qRT-PCR studies of other plant-pathogen interactions, especially, if related plant species are used in combination with bacterial pathogens.

Identification of Suitable Reference Genes for Gene Expression Normalization in the Quantitative Real-Time PCR Analysis of Sweet Osmanthus (Osmanthus fragrans Lour.)

Quantitative real-time PCR (RT-qPCR), a sensitive technique for quantifying gene expression, depends on the stability of the reference gene(s) used for data normalization. Several studies examining the selection of reference genes have been performed in ornamental plants but none in sweet osmanthus (Osmanthus fragrans Lour.). Based on transcriptomic sequencing data from O. fragrans buds at four developmental stages, six reference genes (OfACT, OfEF1α, OfIDH, OfRAN1, OfTUB, and OfUBC2) with stable expression (0.5 to 2 fold change in expression levels between any two developmental stages), as well as the commonly used reference gene Of18S, were selected as candidates for gene expression normalization in the RT-qPCR analysis of O. fragrans. For the normalization of RT-qPCR with two dyes, SYBR Green and EvaGreen, the expressional stability of seven candidate reference genes in 43 O. fragrans samples was analyzed using geNorm, NormFinder and BestKeeper. For RT-qPCR using SYBR Green, OfRAN1 and OfUBC2 were the optimal reference genes for all samples and different cultivars, OfACT and OfEF1α were suitable for different floral developmental stages, and OfACT was the optimal reference gene for different temperature treatments. The geometric mean values of the optimal reference gene pairs for the normalization of RT-qPCR are recommended to be used for all samples, different cultivars and different floral developmental stages in O. fragrans. For RT-qPCR using EvaGreen, OfUBC2 was the optimal reference gene for all samples and different cultivars, and OfACT was the optimal reference gene for different floral developmental stages and different temperature treatments. As the worst reference gene, Of18S should not be used as a reference gene in O. fragrans in the future. Our results provide a reference gene application guideline for O. fragrans gene expression characterization using RT-qPCR.

Sucrose is an early modulator of the key hormonal mechanisms controlling bud outgrowth in Rosa hybrida

Sugar has only recently been identified as a key player in triggering bud outgrowth, while hormonal control of bud outgrowth is already well established. To get a better understanding of sugar control, the present study investigated how sugar availability modulates the hormonal network during bud outgrowth in Rosa hybrida. Other plant models, for which mutants are available, were used when necessary. Buds were grown in vitro to manipulate available sugars. The temporal patterns of the hormonal regulatory network were assessed in parallel with bud outgrowth dynamics. Sucrose determined bud entrance into sustained growth in a concentration-dependent manner. Sustained growth was accompanied by sustained auxin production in buds, and sustained auxin export in a DR5::GUS-expressing pea line. Several events occurred ahead of sucrose-stimulated bud outgrowth. Sucrose upregulated early auxin synthesis genes (RhTAR1, RhYUC1) and the auxin efflux carrier gene RhPIN1, and promoted PIN1 abundance at the plasma membrane in a pPIN1::PIN1-GFP-expressing tomato line. Sucrose downregulated both RwMAX2, involved in the strigolactone-transduction pathway, and RhBRC1, a repressor of branching, at an early stage. The presence of sucrose also increased stem cytokinin content, but sucrose-promoted bud outgrowth was not related to that pathway. In these processes, several non-metabolizable sucrose analogues induced sustained bud outgrowth in R. hybrida, Pisum sativum, and Arabidopsis thaliana, suggesting that sucrose was involved in a signalling pathway. In conclusion, we identified potential hormonal candidates for bud outgrowth control by sugar. They are central to future investigations aimed at disentangling the processes that underlie regulation of bud outgrowth by sugar.

Heat shock factors in carrot: genome-wide identification, classification, and expression profiles response to abiotic stress

Heat shock factors (HSFs) play key roles in the response to abiotic stress in eukaryotes. In this study, 35 DcHSFs were identified from carrot (Daucus carota L.) based on the carrot genome database. All 35 DcHSFs were divided into three classes (A, B, and C) according to the structure and phylogenetic relationships of four different plants, namely, Arabidopsis thaliana, Vitis vinifera, Brassica rapa, and Oryza sativa. Comparative analysis of algae, gymnosperms, and angiosperms indicated that the numbers of HSF transcription factors were related to the plant's evolution. The expression profiles of five DcHsf genes (DcHsf 01, DcHsf 02, DcHsf 09, DcHsf 10, and DcHsf 16), which selected from each subfamily (A, B, and C), were detected by quantitative real-time PCR under abiotic stresses (cold, heat, high salinity, and drought) in two carrot cultivars, D. carota L. cvs. Kurodagosun and Junchuanhong. The expression levels of DcHsfs were markedly increased by heat stress, except that of DcHsf 10, which was down regulated. The expression profiles of different DcHsfs in the same class also differed under various stress treatments. The expression profiles of these DcHsfs were also different in tissues of two carrot cultivars. This study is the first to identify and characterize the DcHSF family transcription factors in plants of Apiaceae using whole-genome analysis. The results of this study provide an in-depth understanding of the DcHSF family transcription factors' structure, function, and evolution in carrot.

Validation of reference housekeeping genes for gene expression studies in western corn rootworm (Diabrotica virgifera virgifera)

Quantitative Real-time PCR (qRT-PCR) is a powerful technique to investigate comparative gene expression. In general, normalization of results using a highly stable housekeeping gene (HKG) as an internal control is recommended and necessary. However, there are several reports suggesting that regulation of some HKGs is affected by different conditions. The western corn rootworm (WCR), Diabrotica virgifera virgifera LeConte (Coleoptera: Chrysomelidae), is a serious pest of corn in the United States and Europe. The expression profile of target genes related to insecticide exposure, resistance, and RNA interference has become an important experimental technique for study of western corn rootworms; however, lack of information on reliable HKGs under different conditions makes the interpretation of qRT-PCR results difficult. In this study, four distinct algorithms (Genorm, NormFinder, BestKeeper and delta-CT) and five candidate HKGs to genes of reference (β-actin; GAPDH, glyceraldehyde-3-phosphate dehydrogenase; β-tubulin; RPS9, ribosomal protein S9; EF1a, elongation factor-1α) were evaluated to determine the most reliable HKG under different experimental conditions including exposure to dsRNA and Bt toxins and among different tissues and developmental stages. Although all the HKGs tested exhibited relatively stable expression among the different treatments, some differences were noted. Among the five candidate reference genes evaluated, β-actin exhibited highly stable expression among different life stages. RPS9 exhibited the most similar pattern of expression among dsRNA treatments, and both experiments indicated that EF1a was the second most stable gene. EF1a was also the most stable for Bt exposure and among different tissues. These results will enable researchers to use more accurate and reliable normalization of qRT-PCR data in WCR experiments.

Reference gene selection for quantitative real-time PCR normalization in Reaumuria soongorica

Despite its superiority for evaluating gene expression, real-time quantitative polymerase chain reaction (qPCR) results can be significantly biased by the use of inappropriate reference genes under different experimental conditions. Reaumuria soongorica is a dominant species of desert ecosystems in arid central Asia. Given the increasing interest in ecological engineering and potential genetic resources for arid agronomy, it is important to analyze gene function. However, systematic evaluation of stable reference genes should be performed prior to such analyses. In this study, the stabilities of 10 candidate reference genes were analyzed under 4 kinds of abiotic stresses (drought, salt, dark, and heat) within 4 accessions (HG010, HG020, XGG030, and XGG040) from 2 different habitats using 3 algorithms (geNorm, NormFinder, and BestKeeper). After validation of the ribulose-1,5-bisphosphate carboxylase/oxygenase large unite (rbcL) expression pattern, our data suggested that histone H2A (H2A) and eukaryotic initiation factor 4A-2 (EIF4A2) were the most stable reference genes, cyclophilin (CYCL) was moderate, and elongation factor 1α (EF1α) was the worst choice. This first systematic analysis for stably expressed genes will facilitate future functional analyses and deep mining of genetic resources in R. soongorica and other species of the Reaumuria genus.