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

Cloning and functional characterization of the peptide deformylase encoding gene EuPDF1B from Eucommia ulmoides Oliv

Peptide deformylase can catalyse the removal of formyl groups from the N-terminal formyl methionine of the primary polypeptide chain. The peptide deformylase genes of a few herbaceous plants have been studied to some extent, but the peptide deformylase genes of woody plants have not been studied. In this study, we isolated EuPDF1B from Eucommia ulmoides Oliv. The full-length sequence of EuPDF1B is 1176 bp long with a poly-A tail and contains an open reading frame of 831 bp that encodes a protein of 276 amino acids. EuPDF1B was localized to the chloroplast. qRT‒PCR analysis revealed that this gene was expressed in almost all tissues tested but mainly in mature leaves. Moreover, the expression of EuPDF1B was enhanced by ABA, MeJA and GA and inhibited by shading treatment. The expression pattern of EuPDF1B was further confirmed in EuPDF1Bp: GUS transgenic tobacco plants. Among all the transgenic tobacco plants, EuPDF1Bp-3 showed the highest GUS histochemical staining and activity in different tissues. This difference may be related to the presence of enhancer elements in the region from - 891 bp to - 236 bp of the EuPDF1B promoter. In addition, the expression of the chloroplast gene psbA and the net photosynthetic rate, fresh weight and height of tobacco plants overexpressing EuPDF1B were greater than those of the wild-type tobacco plants, suggesting that EuPDF1B may promote the growth of transgenic tobacco plants. This is the first time that PDF and its promoter have been cloned from woody plants, laying a foundation for further analysis of the function of PDF and the regulation of its expression.

Positive regulation of the Eucommia rubber biosynthesis-related gene EuFPS1 by EuWRKY30 in Eucommia ulmoides

Eucommia rubber is a secondary metabolite from Eucommia ulmoides that has attracted much attention because of its unique properties and enormous potential for application. However, the transcriptional mechanism regulating its biosynthesis has not yet been determined. Farnesyl pyrophosphate synthase is a key enzyme in the Eucommia rubber biosynthesis. In this study, the promoter of EuFPS1 was used as bait, EuWRKY30 was screened from the cDNA library of EuFPS1 via a yeast one-hybrid system. EuWRKY30 belongs to the WRKY IIa subfamily and contains a WRKY domain and a C2H2 zinc finger motif, and the expressed protein is located in the nucleus. EuWRKY30 and EuFPS1 exhibited similar tissue expression patterns, and yeast one-hybrid and dual-luciferase experiments confirmed that EuWRKY30 directly binds to the W-box element in the EuFPS1 promoter and activates its expression. Moreover, the overexpression of EuWRKY30 significantly upregulated the expression level of EuFPS1, further increasing the density of the rubber particles and Eucommia rubber content. The results of this study indicated that EuWRKY30 positively regulates EuFPS1, which plays a critical role in the synthesis of Eucommia rubber, provided a basis for further analysis of the underlying transcriptional regulatory mechanisms.

StCoExpNet: a global co-expression network analysis facilitates identifying genes underlying agronomic traits in potatoes

KEY MESSAGE

We constructed a gene expression atlas and co-expression network for potatoes and identified several novel genes associated with various agronomic traits. This resource will accelerate potato genetics and genomics research. Potato (Solanum tuberosum L.) is the world's most crucial non-cereal food crop and ranks third in food production after wheat and rice. Despite the availability of several potato transcriptome datasets at public databases like NCBI SRA, an effort has yet to be put into developing a global transcriptome atlas and a co-expression network for potatoes. The objectives of our study were to construct a global expression atlas for potatoes using publicly available transcriptome datasets, identify housekeeping and tissue-specific genes, construct a global co-expression network and identify co-expression clusters, investigate the transcriptional complexity of genes involved in various essential biological processes related to agronomic traits, and provide a web server (StCoExpNet) to easily access the newly constructed expression atlas and co-expression network to investigate the expression and co-expression of genes of interest. In this study, we used data from 2299 publicly available potato transcriptome samples obtained from 15 different tissues to construct a global transcriptome atlas. We found that roughly 87% of the annotated genes exhibited detectable expression in at least one sample. Among these, we identified 281 genes with consistent and stable expression levels, indicating their role as housekeeping genes. Conversely, 308 genes exhibited marked tissue-specific expression patterns. We exemplarily linked some co-expression clusters to important agronomic traits of potatoes, such as self-incompatibility, anthocyanin biosynthesis, tuberization, and defense responses against multiple pathogens. The dataset compiled here constitutes a new resource (StCoExpNet), which can be accessed at https://stcoexpnet.julius-kuehn.de . This transcriptome atlas and the co-expression network will accelerate potato genetics and genomics research.

Selection and validation of reference genes for normalization of gene expression in Floccularia luteovirens

Floccularia luteovirens is one of the rare edible fungi with high nutritional value found on the Qinghai-Tibet Plateau. However, research at the molecular level on this species is currently constrained due to the lack of reliable reference genes for this species. Thirteen potential reference genes (ACT, GAPDH, EF-Tu, SAMDC, UBI, CLN1, β-TUB, γ-TUB, GTP, H3, UBC, UBC-E2, and GTPBP1) were chosen for the present study, and their expression under various abiotic conditions was investigated. Stability of gene expression was tested using GeNorm, NormFinder, BestKeeper, Delta-Ct, and RefFinder. The results showed that the most suitable reference genes for salt treatment were ACT and EF-Tu. Under drought stress, γ-TUB and UBC-E2 would be suitable for normalization. Under oxidative stress, the reference genes H3 and GAPDH worked well. Under heat stress, the reference genes EF-Tu and γ-TUB were suggested. Under extreme pH stress, UBC-E2 and H3 were appropriate reference genes. Under cadmium stress, the reference genes ACT and UBC-E2 functioned well. In different tissues, H3 and GTPBP1 were appropriate reference genes. The optimal internal reference genes when analyzing all samples were H3 and SAMDC. The expression level of HSP90 was studied to further validate the applicability of the genes identified in this study.

Selection and validation of internal control genes for quantitative real-time RT‒qPCR normalization of Phlebopus portentosus gene expression under different conditions

Phlebopus portentosus (Berk. and Broome) Boedijn is an attractive edible mushroom and is considered the only bolete for which artificial cultivation in vitro has been achieved. Gene expression analysis has become widely used in research on edible fungi and is important for elucidating the functions of genes involved in complex biological processes. Selecting appropriate reference genes is crucial to ensuring reliable RT‒qPCR gene expression analysis results. In our study, a total of 12 candidate control genes were selected from 25 traditional housekeeping genes based on their expression stability in 9 transcriptomes of 3 developmental stages. These genes were further evaluated using geNorm, NormFinder, and RefFinder under different conditions and developmental stages. The results revealed that MSF1 domain-containing protein (MSF1), synaptobrevin (SYB), mitogen-activated protein kinase genes (MAPK), TATA-binding protein 1 (TBP1), and SPRY domain protein (SPRY) were the most stable reference genes in all sample treatments, while elongation factor 1-alpha (EF1), actin and ubiquitin-conjugating enzyme (UBCE) were the most unstably expressed. The gene SYB was selected based on the transcriptome results and was identified as a novel reference gene in P. portentosus. This is the first detailed study on the identification of reference genes in this fungus and may provide new insights into selecting genes and quantifying gene expression.

Evaluation of reference genes for qRT-PCR studies in the colchicine producing Gloriosa superba L

The flame lily, Gloriosa superba L., is one of the two primary sources of the anti-inflammatory drug, colchicine. Previous studies have shown that a higher level of colchicine production occurs in the rhizomes than in leaves and roots. Earlier precursor feeding and transcriptome analysis of G. superba have provided a putative pathway and candidate genes involved in colchicine biosynthesis. Comparative analysis of expression levels of candidate pathway genes in different tissues of G. superba using quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR) can reveal highly expressed genes in the rhizome compared to other tissues which could suggest roles of the gene products in colchicine biosynthesis. Normalization is an important step in effectively analyzing differential gene expression by qRT-PCR with broader applications. The current study selected candidate reference genes from the transcriptome datasets and analyzed them to determine the most stable genes for normalization of colchicine biosynthesis-related genes. Using RefFinder, one stable reference gene, UBC22, was selected to normalize gene expression levels of candidate methyltransferase (MT) genes in the leaves, roots, and rhizomes of G. superba. With UBC22 as reference gene, the methyltransferases, GsOMT1, GsOMT3, and GsOMT4 showed significantly higher expression levels in the rhizome of G. superba, while MT31794 was more highly expressed in the roots. In conclusion, the current results showed a viable reference gene expression analysis system that could help elucidate colchicine biosynthesis and its exploitation for increased production of the drug in G. superba.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11816-023-00840-x.

Predatory stink bug, Eocanthecona furcellata (Wolff) responses to oral exposure route of λ-cyhalothrin via sex-specific modulation manner

The toxic effects of insecticides on predatory arthropods have closely related to their exposure routes. However, little is known about the effects of insecticide on reproductive parameters when the route of exposure occurs at a trophic level via prey intake. We therefore conducted current studies assessing whether Eocanthecona furcellata adults would be affected by feeding with λ-cyhalothrin-contaminated prey. Reproductive parameters, i.e. prolonged premating and preoviposition durations, reduced number of egg batches and egg amount, disturbed ovarian development, and suppressed expression of reproductive related genes were observed in E. furcellata females by feeding with treated prey. Moreover, reduced survival rate and altered carbohydrate metabolism parameters were detected in male bugs. Biochemical parameters, including MDA content, the activities of three antioxidant enzymes and three detoxification enzymes exhibited sex-specific responses after oral-exposure to λ-cyhalothrin in E. furcellata. The results indicate that the insecticide affects the fitness and leads to impairing reproductive potential via sex-specific modulation manner in predator insects. Taken together, our results provide a comprehensive assessment about detrimental impacts of λ-cyhalothrin-exposure on predators via prey intake, as well as a solid basis for further research to protect the predators from hazardous impacts of insecticides.

The NAC transcription factor family in Eucommia ulmoides: Genome-wide identification, characterization, and network analysis in relation to the rubber biosynthetic genes

The NAC transcription factor family is a large plant gene family, participating in plant growth and development, secondary metabolite synthesis, biotic and abiotic stresses responses, and hormone signaling. Eucommia ulmoides is a widely planted economic tree species in China that can produce trans-polyisoprene: Eucommia rubber (Eu-rubber). However, genome-wide identification of the NAC gene family has not been reported in E. ulmoides. In this study, 71 NAC proteins were identified based on genomic database of E. ulmoides. Phylogenetic analysis showed that the EuNAC proteins were distributed in 17 subgroups based on homology with NAC proteins in Arabidopsis, including the E. ulmoides-specific subgroup Eu_NAC. Gene structure analysis suggested that the number of exons varied from 1 to 7, and multitudinous EuNAC genes contained two or three exons. Chromosomal location analysis revealed that the EuNAC genes were unevenly distributed on 16 chromosomes. Three pairs of genes of tandem duplicates genes and 12 segmental duplications were detected, which indicated that segmental duplications may provide the primary driving force of expansion of EuNAC. Prediction of cis-regulatory elements indicated that the EuNAC genes were involved in development, light response, stress response and hormone response. For the gene expression analysis, the expression levels of EuNAC genes in various tissues were quite different. To explore the effect of EuNAC genes on Eu-rubber biosynthesis, a co-expression regulatory network between Eu-rubber biosynthesis genes and EuNAC genes was constructed, which indicated that six EuNAC genes may play an important role in the regulation of Eu-rubber biosynthesis. In addition, this six EuNAC genes expression profiles in E. ulmoides different tissues were consistent with the trend in Eu-rubber content. Quantitative real-time PCR analysis showed that EuNAC genes were responsive to different hormone treatment. These results will provide a useful reference for further studies addressing the functional characteristics of the NAC genes and its potential role in Eu-rubber biosynthesis.

Selection and validation of reference genes for RT-qPCR analysis of different organs at various development stages in Caragana intermedia

Reverse transcription quantitative PCR (RT-qPCR) is a technique widely used to investigate the expression of genes. An appropriate reference gene (RG) is essential for RT-qPCR analysis to obtain accurate and reliable results. Caragana intermedia plays an important role in afforestation as a bush. However, due to the lack of appropriate RGs, the research on development-related genes is limited. In this study, the selection for suitable RGs of different organs at various development stages to normalize the results of RT-qPCR about development-related genes was performed. To test the expression stability across all samples, we used the software algorithms such as geNorm, NormFinder, BestKeeper, and RefFinder to evaluate all the candidate RGs. Our results showed that CiEF1α was the most stable RG with little fluctuation among all samples. In addition, CiGAPDH in roots, CiSKIP1 in stems and leaves, and CiEF1α in different organs were selected as the most stable RGs. To confirm the applicability of the most stable RGs, the relative expression of CiWRKY17 was normalized using different candidate RGs. Taken together, our research laid a foundation for the study of development-related genes in C. intermedia.

Selection and Validation of Reference Genes for qRT-PCR Gene Expression Analysis in Kengyilia melanthera

Kengyilia is a newly established genus. Most species in this genus survive in hash environment, which might be an indicator of an acquirement of stress resistance genes and the potential for molecular breeding in Triticeae species. Quantitative real-time PCR (qRT-PCR) is a widely used technique with varied sensitivity heavily dependent on the optimal level of the reference genes. K. melanthera is a typical psammophyte species which has high drought resistance. The reference genes of K. melanthera are not yet reported. This study aims to evaluate the expression stability of 14 candidate reference genes (EF1A, GAPDH, ACT1, UBI, TUBB3, TIPRL, CACS, PPP2R1B, TUBA1A, EIF4A1, CYPA3, TCTP, ABCG11L, and FBXO6L) under five treatments (drought, heat, cold, salt, and ABA) and find the most stable and suitable one even upon stressed conditions. The software NormFinder, GeNorm, BestKeeper, and RefFinder were used for data analysis. In general, the genes CACS and PPP2R1B are concluded to have the best overall performance under the various treatments. With the ABA treatment, TCTP and TIPRL show the best stability. CACS and TCTP, as well as TIPRL and CYPA3, were most stable under the treatments of cold and salt, respectively. CACS and FBXO6L were ranked the highest with the heat treatment and drought treatment, respectively. Finally, the Catalase-1 (CAT1) gene was used to verify the reliability of the above reference genes. Accordingly, CAT1’s expression pattern remained unchanged after normalization with stable reference genes. This study provides beneficial information about the stability and reliability of potential reference genes for qRT-PCR in K. melanthera.

Prognostic value of SOX9 in cervical cancer: Bioinformatics and experimental approaches

Among gynecological cancers, cervical cancer is a common malignancy and remains the leading cause of cancer-related death for women. However, the exact molecular pathogenesis of cervical cancer is not known. Hence, understanding the molecular mechanisms underlying cervical cancer pathogenesis will aid in the development of effective treatment modalities. In this research, we attempted to discern candidate biomarkers for cervical cancer by using multiple bioinformatics approaches. First, we performed differential expression analysis based on cervical squamous cell carcinoma and endocervical adenocarcinoma data from The Cancer Genome Atlas database, then used differentially expressed genes for weighted gene co-expression network construction to find the most relevant gene module for cervical cancer. Next, the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses were performed on the module genes, followed by using protein–protein interaction network analysis and Cytoscape to find the key gene. Finally, we validated the key gene by using multiple online sites and experimental methods. Through weighted gene co-expression network analysis, we found the turquoise module was the highest correlated module with cervical cancer diagnosis. The biological process of the module genes focused on cell proliferation, cell adhesion, and protein binding processes, while the Kyoto Encyclopedia of Genes and Genomes pathway of the module significantly enriched pathways related to cancer and cell circle. Among the module genes, SOX9 was identified as the hub gene, and its expression was associated with cervical cancer prognosis. We found the expression of SOX9 correlates with cancer-associated fibroblast immune infiltration in immune cells by Timer2.0. Furthermore, cancer-associated fibroblast infiltration is linked to cervical cancer patients’ prognosis. Compared to those in normal adjacent, immunohistochemical and real-time quantitative polymerase chain reaction (qPCR) showed that the protein and mRNA expression of SOX9 in cervical cancer were higher. Therefore, the SOX9 gene acts as an oncogene in cervical cancer, interactive with immune infiltration of cancer-associated fibroblasts, thereby affecting the prognosis of patients with cervical cancer.

Selection and validation of reference genes for quantitative real-time PCR of Quercus mongolica Fisch. ex Ledeb under abiotic stresses

Quercus mongolica Fisch. ex Ledeb is the main species of coniferous and broadleaved mixed forests in northeast and north China, which has high ornamental, economic, and ecological value. The appropriate reference genes must be selected for quantitative real-time PCR to reveal the molecular mechanisms of stress responses and their contribution to breeding of Q. mongolica. In the present study, we chose 11 candidate reference genes (TUA, CYP18, HIS4, RPS13, ACT97, TUB1, UBQ10, UBC5, SAND, PP2A, and SAMDC) and used four programs (GeNorm, NormFinder, BestKeeper, and RefFinder) to assess the expression stability of the above genes in roots, stems, and leaves under five abiotic stress factors (cold, salt, drought, weak light, and heavy metal). The findings revealed that under various experimental environments, the most stable genes were different; CYP18, ACT97, and RPS13 ranked the highest under most experimental environments. Moreover, two genes induced by stress, CMO and P5CS2, were chosen to demonstrate the reliability of the selected reference genes in various tissues under various stress conditions. Our research provides a significant basis for subsequent gene function studies of Q. mongolica.

Absolute quantification of E. coli virulence and housekeeping genes to determine pathogen loads in enumerated environmental samples

Quantifying pathogenic genes with q-PCR in complex samples to determine the pathogen loads is influenced by a wide range of factors, including choice of extraction method, standard curve, and the decision to use relative versus absolute quantification of the genes. The aim was to investigate the standardisation of q-PCR methods to determine enumerated E. coli gene ratios grown with the IDEXX Colilert^® Quanti-Trays^® using enteropathogenic E. coli as the model pathogen. q-PCR targeting the eaeA and gadAB genes was used to calculate the eaeA: gadAB ratios for clinical strains collected between [2005–2006 (n = 55)] and [2008–2009 (n = 19)] using the LinRegPCR software and Corbett Research Thermal cycler software. Both programs grouped the isolates into two distinct groups based on the gene ratios although the Corbett Research Thermal cycler software gave results one log higher than the LinRegPCR program. Although the eaeA: gadAB ratio range was determined using extracted E. coli DNA, the impact of free DNA and other bacteria present in the sample needed to be understood. Standard curve variations using serially diluted extracted E. coli DNA, serially diluted pure E. coli culture followed by DNA extraction from each dilution with or without other bacteria was tested using the eaeA q-PCR to quantify the genes. Comparison of the standard curves showed no significant difference between standard curves prepared with diluted DNA or with cells diluted before the DNA is extracted (P = 0.435). Significant differences were observed when background DNA was included in the diluent or Coliform cells added to the diluent to dilute cells before the DNA is extracted (P < 0.001). The “carrier” DNA and Coliform cells enhanced the DNA extraction results resulting in better PCR efficiency. This will have an influence on the quantification of gene ratios and pathogen load in samples containing lower numbers of E. coli.

Fulvic acid enhances drought resistance in tea plants by regulating the starch and sucrose metabolism and certain secondary metabolism

The aim of this work was to gain insight into the molecular mechanisms underlying the effect of fulvic acid on drought-exposed tea plants. We performed proteomic analysis of fulvic acid-treated tea leaves from the target plants using tandem mass tag quantitative labeling technology and compared the results with those of a previous transcriptomic analysis. We identified 48 and 611 differentially abundant proteins in the leaves of tea plants treated with fulvic acid compared with the control under mild and severe drought, respectively. Comparative analysis showed that, under severe drought, 55 genes had similar expression patterns at the transcriptome and proteome levels, such as PAL, GBE, GBSS and bAS. Bioinformatic analysis revealed that those genes were mainly related to the starch and sucrose metabolism, phenylpropanoid biosynthesis and triterpenoid biosynthesis. SIGNIFICANCE: This study broadens the understanding of the molecular mechanisms underlying the improved drought resistance seen in tea plants in the presence of fulvic acid and provides a basis for further research on the genomics of drought tolerance in these plants. In addition, these findings could be used to develop new guidance strategies for improved drought management systems in tea plantation.

Multinuclei Occurred Under Cryopreservation and Enhanced the Pathogenicity of Melampsora larici-populina

Melampsora larici-populina is a macrocyclic rust, and the haploid stage with two nuclei and the diploid of mononuclear sequentially occur annually. During the preservation of dry urediniospores at −80°C, we found that one isolate, ΔTs₀₆, was different from the usual wild-type isolate Ts₀₆ at −20°C because it has mixed polykaryotic urediniospores. However, the other spores, including the 0, I, III, and IV stages of a life cycle, were the same as Ts₀₆. After five generations of successive inoculation and harvest of urediniospores from the compatible host Populus purdomii, the isolate ΔTs₀₆ steadily maintained more than 20% multiple nucleus spores. To test the pathogenesis variation of ΔTs₀₆, an assay of host poplars was applied to evaluate the differences between ΔTs₀₆ and Ts₀₆. After ΔTs₀₆ and Ts₀₆ inoculation, leaves of P. purdomii were used to detect the expression of small secreted proteins (SSPs) and fungal biomasses using quantitative real-time PCR (qRT-PCR) and trypan blue staining. ΔTs₀₆ displayed stronger expression of five SSPs and had a shorter latent period, a higher density of uredinia, and higher DNA mass. A transcriptomic comparison between ΔTs₀₆ and Ts₀₆ revealed that 3,224 were differentially expressed genes (DEGs), 55 of which were related to reactive oxygen species metabolism, the Mitogen-activated protein kinase (MAPK) signaling pathway, and the meiosis pathway. Ten genes in the mitotic and meiotic pathways and another two genes associated with the “response to DNA damage stimulus” all had an upward expression, which were detected by qRT-PCR in ΔTs₀₆ during cryopreservation. Gas chromatography–mass spectrometry (GC-MS) confirmed that the amounts of hexadecanoic acid and octadecadienoic acid were much more in ΔTs₀₆ than in Ts₀₆. In addition, using spectrophotometry, hydrogen peroxide (H₂O₂) was also present in greater quantities in ΔTs₀₆ compared with those found in Ts₀₆. Increased fatty acids metabolism could prevent damage to urediniospores in super-low temperatures, but oxidant species that involved H₂O₂ may destroy tube proteins of mitosis and meiosis, which could cause abnormal nuclear division and lead to multinucleation, which has a different genotype. Therefore, the multinuclear isolate is different from the wild-type isolate in terms of phenotype and genotype; this multinucleation phenomenon in urediniospores improves the pathogenesis and environmental fitness of M. larici-populina.

Selection and Verification of Appropriate Reference Genes for Expression Normalization in Cryptomeria fortunei under Abiotic Stress and Hormone Treatments

Cryptomeria fortunei has become one of the main timber afforestation species in subtropical high-altitude areas of China due to its fast growth, good material quality, and strong adaptability, showing broad application prospects. Quantitative real-time PCR (qRT-PCR) is the most accurate and widely used gene expression evaluation technique, and selecting appropriate reference genes (RGs) is essential for normalizing qRT-PCR results. However, suitable RGs for gene expression normalization in C. fortunei have not been reported. Here, we tested the expression stability for 12 RGs in C. fortunei under various experimental conditions (simulated abiotic stresses (cold, heat, drought, and salinity) and hormone treatments (methyl jasmonate, abscisic acid, salicylic acid, and gibberellin) and in different tissues (stems, tender needles, needles, cones, and seeds) using four algorithms (delta Ct, geNorm, NormFinder, and BestKeeper). Then, geometric mean rankings from these algorithms and the RefFinder program were used to comprehensively evaluate RG stability. The results indicated CYP, actin, UBC, and 18S as good choices for studying C. fortunei gene expression. qRT-PCR analysis of the expression patterns of three target genes (CAT and MAPK1/6) further verified that the selected RGs were suitable for gene expression normalization. This study provides an important basis for C. fortunei gene expression standardization and quantification.

Identification, evolution and expression analysis of WRKY gene family in Eucommia ulmoides

The WRKY transcription factors is one of the largest families of transcription factors (TFs) in plants and involved in multiple biological processes. However, the role of the WRKY family had not been reported in Eucommia ulmoides. In this study, 45 WRKY genes (EuWRKY1-45) with conserved WRKY domain were identified in E. ulmoides and classified into three groups. The group II was further divided into five subgroups based on phylogenetic analysis, and each clade was well supported by the conserved motifs. All the genes were located on 34 different scaffolds respectively. A number of development-, light-, hormone-, and stress-related elements were randomly distributed in the promoter sequences of EuWRKYs. Expression profiles indicated that EuWRKY genes were involved in leaf development, and majority of EuWRKYs genes were highly expressed in leaf buds. Co-expression analysis of WRKYs suggested an intricate interplay of growth-related responses. EuWRKY4 was involved in a complex proteins interaction network. Collectively, our results provide extensive insights into the WRKY gene family, thereby contributing to the screening of additional candidate genes in E. ulmoides.

Reference gene selection for qRT-PCR analyses of luffa (Luffa cylindrica) plants under abiotic stress conditions

Selecting suitable internal reference genes is an important prerequisite for the application of quantitative real-time PCR (qRT-PCR). However, no systematic studies have been conducted on reference genes in luffa. In this study, seven reference genes were selected, and their expression levels in luffa plants exposed to various simulated abiotic stresses [i.e., cold, drought, heat, salt, H₂O₂, and abscisic acid (ABA) treatments] were analyzed by qRT-PCR. The stability of the reference gene expression levels was validated using the geNorm, NormFinder, BestKeeper, and RefFinder algorithms. The results indicated that EF-1α was the most stably expressed and suitable reference gene overall and for the heat, cold, and ABA treatments. Additionally, UBQ expression was stable following the salt treatment, whereas TUB was identified as a suitable reference gene for H₂O₂ and drought treatments. The reliability of the selected reference genes was verified by analyzing the expression of copper/zinc superoxide dismutase (Cu/Zn-SOD) gene in luffa. When the most unstable reference genes were used for data normalizations, the resulting expression patterns had obvious biases when compared with the expression patterns for the most ideal reference genes used alone or combined. These results will be conducive to more accurate quantification of gene expression levels in luffa.

Selection and validation of reliable reference genes for gene expression studies from Monilinia vaccinii-corymbosi infected wild blueberry phenotypes

Monilinia blight disease caused by Monilinia vaccinii-corymbosi (Reade) Honey (M.vc) causes severe damage and economic losses in wild blueberry growing regions. Molecular mechanisms regulating defence responses of wild blueberry phenotypes towards this causal fungus are not yet fully known. A reliable quantification of gene expression using quantitative real time PCR (qPCR) is fundamental for measuring changes in target gene expression. A crucial aspect of accurate normalisation is the choice of appropriate reference genes. This study evaluated the expression stability of seven candidate reference genes (GAPDH, UBC9, UBC28, TIP41, CaCSa, PPR and RH8) in floral tissues of diploid and tetraploid wild blueberry phenotypes challenged with M.vc. The expression stability was calculated using five algorithms: geNorm, NormFinder, BestKeeper, deltaCt and RefFinder. The results indicated that UBC9 and GAPDH were the most stable reference genes, while RH8 and PPR were the least stable ones. To further validate the suitability of the analyzed reference genes, the expression level of a pathogenesis related protein gene (i.e., PR3) was analysed for both phenotypes at four time points of infection. Our results may be beneficial for future studies involving the quantification of relative gene expression levels in wild blueberry species.

Transcriptome analysis of terpenoid biosynthetic genes and simple sequence repeat marker screening in Eucommia ulmoides

Trans-polyisoprene rubber is produced in the tissues of leaves, bark, and fruit of Eucommia ulmoides and is considered an important energy source. Transcript profiles of two tissues from E. ulmoides cv. Qinzhong No. 3, leaf and fruit, were analysed using the Illumina HiSeq 2000 system. In total, 104 million clean reads were obtained and assembled into 58,863 unigenes. Through gene functional classification, 28,091 unigenes (47.72%) were annotated and 65 unigenes have been hypothesized to encode proteins involved in terpenoid biosynthesis. In addition, 10,041 unigenes were detected as differentially expressed unigenes, and 29 of them were putatively related to terpenoid biosynthesis. The synthesis of trans-polyisoprene rubbers in E. ulmoides was hypothesised to be dominated by the mevalonate pathway. Farnesyl diphosphate synthase 2 (FPPS2) was considered a key component in the biosynthesis of trans-polyprenyl diphosphate. Rubber elongation factor 3 (REF3) might be involved in stabilising the membrane of rubber particles in E. ulmoides. To date, 351 simple sequence repeats (SSRs) were validated as polymorphisms from eight E. ulmoides plants (two parent plants and six F₁ individuals), and these could act as molecular markers for genetic map density increase and breeding improvement of E. ulmoides.

Identification and evaluation of suitable reference genes for gene expression analysis in rubber tree leaf

Gene expression profiles are increasingly applied to investigate molecular mechanism for which, normalization with suitable reference genes is critical. Previously we have reported several suitable reference genes for laticifer samples from rubber tree, however, little is known in leaf. The main objective of this current study was to identify some stable expression reference genes at various developmental stages of leaf, as well as during abiotic (high and low temperature extremes) and biotic stresses (pathogen stress). Gene expression profilings identified the ubiquitin-proteasome system as excellent potential as reference genes for rubber tree leaf. Among a total of 30 tested genes investigated, 24 new candidate (including 11 genes involved in the ubiquitin-proteasome system), 4 previously identified and 2 specific genes, were further evaluated using quantitative real-time PCR. Our results indicated that the new candidate genes had better expression stability comparing with others. For instance, an ubiquitin conjugating enzyme (RG0099) and three ubiquitin-protein ligases (RG0928, RG2190 and RG0118) expressed stably in all samples, and were confirmed to be suitable reference genes for rubber tree leaf under four different conditions. Finally, we suggest that using more than one reference gene may be appropriate in gene expression studies when employing different software to normalize gene expression data. Our findings have significant implications for the reliability of data obtained from genomics studies in rubber tree and perhaps in other species.

Identification and evaluation of reference genes for quantitative real-time PCR analysis in Polygonum cuspidatum based on transcriptome data

BACKGROUND

Polygonum cuspidatum of the Polygonaceae family is a traditional medicinal plant with many bioactive compounds that play important roles in human health and stress responses. Research has attempted to identify biosynthesis genes and metabolic pathways in this species, and quantitative real-time PCR (RT-qPCR) has commonly been used to detect gene expression because of its speed, sensitivity, and specificity. However, no P. cuspidatum reference genes have been identified, which hinders gene expression studies. Here, we aimed to identify suitable reference genes for accurate and reliable normalization of P. cuspidatum RT-qPCR data.

RESULTS

Twelve candidate reference genes, including nine common (ACT, TUA, TUB, GAPDH, EF-1γ, UBQ, UBC, 60SrRNA, and eIF6A) and three novel (SKD1, YLS8, and NDUFA13), were analyzed in different tissues (root, stem, and leaf) without treatment and in leaves under abiotic stresses (salt, ultraviolet [UV], cold, heat, and drought) and hormone stimuli (abscisic acid [ABA], ethylene [ETH], gibberellin [GA₃], methyl jasmonate [MeJA], and salicylic acid [SA]). Expression stability in 65 samples was calculated using the △CT method, geNorm, NormFinder, BestKeeper, and RefFinder. Two reference genes (NDUFA13 and EF-1γ) were sufficient to normalize gene expression across all sample sets. They were also the two most stable genes for abiotic stresses and different tissues, whereas NDUFA13 and SKD1 were the top two choices for hormone stimuli. Considering individual experimental sets, GAPDH was the top-ranked gene under ABA, ETH, and GA₃ treatments, while 60SrRNA showed good stability under MeJA and cold treatments. ACT, UBC, and TUB were suitable genes for drought, UV, and ABA treatments, respectively. TUA was not suitable because of its considerable variation in expression under different conditions. The expression patterns of PcPAL, PcSTS, and PcMYB4 under UV and SA treatments and in different tissues normalized by stable and unstable reference genes demonstrated the suitability of the optimal reference genes.

CONCLUSIONS

We propose NDUFA13 and EF-1γ as reference genes to normalize P. cuspidatum expression data. To our knowledge, this is the first systematic study of reference genes in P. cuspidatum which could help advance molecular biology research in P. cuspidatum and allied species.

Selection of Reference Genes for Normalization of Real-Time PCR Data in Calliptamus italicus (Orthoptera: Acrididae) Under Different Temperature Conditions

Global warming has dominated worldwide climate change trends, and adaptability to high temperatures is the main factor underlying the spread of the pest Calliptamus italicus in Xinjiang Province, China. However, knowledge about the molecular mechanisms responsible for this adaptability and other related biological properties of C. italicus remain relatively unclear. Real-time quantitative polymerase chain reaction (RT-qPCR) is a key tool for gene expression analysis associated with various biological processes. Reference genes are necessary for normalizing gene expression levels across samples taken from specific experimental conditions. In this study, transcript level of five genes (GAPDH, 18S, TUB, ACT, and EF1α), commonly used as reference genes, were evaluated under nine different temperatures (27, 30, 33, 36, 39, 42, 45, 48, and 51°C) to assess their expression stability and further select the most suitable to be used on normalization of target gene expression data. Gene expression profiles were analyzed using geNorm, NormFinder, and BestKeeper software packages. The combined results demonstrated that the best-ranked reference genes for C. italicus are EF1α, GAPDH, and ACT under different thermal stress conditions. This is the first study that assesses gene expression analysis across a range of temperatures to select the most appropriate reference genes for RT-qPCR data normalization in C. italicus. These results should assist target gene expression analysis associated with heat stress in C. italicus.

Integration of Transcriptomes, Small RNAs, and Degradome Sequencing to Identify Putative miRNAs and Their Targets Related to Eu-Rubber Biosynthesis in Eucommia ulmoides

Eucommia ulmoides has attracted much attention as a valuable natural rubber (Eu-rubber) production tree. As a strategic material, Eu-rubber plays a vital role in general and defence industries. However, the study of Eu-rubber biosynthesis at a molecular level is scarce, and the regulatory network between microRNAs (miRNAs) and messenger RNAs (mRNAs) in Eu-rubber biosynthesis has not been assessed. In this study, we comprehensively analyzed the transcriptomes, small RNAs (sRNAs) and degradome to reveal the regulatory network of Eu-rubber biosynthesis in E. ulmoides. A total of 82,065 unigenes and 221 miRNAs were identified using high-throughput sequencing; 20,815 targets were predicted using psRNATarget software. Of these targets, 779 miRNA-target pairs were identified via degradome sequencing. Thirty-one miRNAs were differentially expressed; 22 targets of 34 miRNAs were annotated in the terpenoid backbone biosynthesis pathway (ko00900) based on the Kyoto Encyclopedia of Genes and Genomes (KEGG). These miRNAs were putatively related to Eu-rubber biosynthesis. A regulatory network was constructed according to the expression profiles of miRNAs and their targets. These results provide a comprehensive analysis of transcriptomics, sRNAs and degradome to reveal the Eu-rubber accumulation, and provide new insights into genetic engineering techniques which may improve the content of Eu-rubber in E. ulmoides.

Dynamic Changes in Metabolite Accumulation and the Transcriptome during Leaf Growth and Development in Eucommia ulmoides

Eucommia ulmoides Oliver is widely distributed in China. This species has been used mainly in medicine due to the high concentration of chlorogenic acid (CGA), flavonoids, lignans, and other compounds in the leaves and barks. However, the categories of metabolites, dynamic changes in metabolite accumulation and overall molecular mechanisms involved in metabolite biosynthesis during E. ulmoides leaf growth and development remain unknown. Here, a total of 515 analytes, including 127 flavonoids, 46 organic acids, 44 amino acid derivatives, 9 phenolamides, and 16 vitamins, were identified from four E. ulmoides samples using ultraperformance liquid chromatography-mass spectrometry (UPLC-MS) (for widely targeted metabolites). The accumulation of most flavonoids peaked in growing leaves, followed by old leaves. UPLC-MS analysis indicated that CGA accumulation increased steadily to a high concentration during leaf growth and development, and rutin showed a high accumulation level in leaf buds and growing leaves. Based on single-molecule long-read sequencing technology, 69,020 transcripts and 2880 novel loci were identified in E. ulmoides. Expression analysis indicated that isoforms in the flavonoid biosynthetic pathway and flavonoid metabolic pathway were highly expressed in growing leaves and old leaves. Co-expression network analysis suggested a potential direct link between the flavonoid and phenylpropanoid biosynthetic pathways via the regulation of transcription factors, including MYB (v-myb avian myeloblastosis viral oncogene homolog) and bHLH (basic/helix-loop-helix). Our study predicts dynamic metabolic models during leaf growth and development and will support further molecular biological studies of metabolite biosynthesis in E. ulmoides. In addition, our results significantly improve the annotation of the E. ulmoides genome.

Transcriptome analysis reveals novel insights into the continuous cropping induced response in Codonopsis tangshen, a medicinal herb

Codonopsis tangshen Oliv. (C. tangshen Oliv.), a famous medicinal herb in China, is seriously affected by continuous cropping (C-cro). The physiological and biochemical results indicated that C-cro significantly affected the malonaldehyde (MDA) and chlorophyll content, as well as activities of catalase (CAT) and superoxide dismutase (SOD) when compared with the non-continuous cropping (NC-cro) group. Transcriptome profiling found 762 differentially expressed genes, including 430 up-regulated and 332 down-regulated genes by C-cro. In addition, pathway enrichment analysis revealed that genes related to 'Tyrosine degradation I', 'Glycogen synthesis' and 'Phenylalanine and tyrosine catabolism' were up-regulated, and genes associated with 'Signal transduction', 'Immune system', etc. were down-regulated by C-cro. The expression of target genes was further validated by Q-PCR. In this study, we demonstrated the effects of C-cro on C. tangshen at the transcriptome level, and found possible C-cro responsive candidate genes. These findings could be further beneficial for improving the continuous cropping tolerance.

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.

Abnormally expressed long noncoding RNA B3GALT5-AS1 may serve as a biomarker for the diagnostic and prognostic of gastric cancer

Early diagnosis of gastric cancer (GC) is an effective method to improve prognosis. Increasing number of long noncoding RNAs (lncRNAs) have been reported as biomarkers for several cancers. We aim to detect the level of lncRNA B3GALT5-AS1 and its association with clinical parameters and to further explore its application value in GC. We measured serum B3GALT5-AS1 expression in 107 patients with GC, 40 polyp patients, and 87 normal controls to explore the significance of serum B3GALT5-AS1 in GC using the quantitative real-time polymerase chain reaction method. The result demonstrated that B3GALT5-AS1 level was markedly richer in GC patients than that in normal people (P < .001). B3GALT5-AS1 may be served as a diagnostic marker for distinguishing GC patients from healthy people, and the proportion under the receiver operating characteristics curve is 0.816 (95% confidence interval, 0.758-0.874; P = .03). Further exploration validated that high serum B3GALT5-AS1 level was related to TNM stage (P = .024), and lymph node metastasis (P = .023). Our study suggested that serum B3GALT5-AS1 may be employed as an ideal biomarker for early screening of GC.