Selection and validation of reference genes for quantitative expression analysis of miRNAs and mRNAs in Poplar

The Mechanism of TRIM21 Inhibiting the Invasion and Migration of ccRCC by Stabilizing ASS1

Clear cell renal cell carcinoma (ccRCC) is characterized by its aggressive invasion and metastasis, presenting significant clinical challenges. Gaining insights into the molecular mechanisms underlying its progression is crucial for the development of effective therapeutic strategies. Addressing a critical knowledge gap in understanding ccRCC tumorigenesis, this study aims to elucidate the expression patterns of TRIM21 in ccRCC, unravel its impact on ccRCC patient prognosis, and investigate the regulatory role of TRIM21 in ASS1 expression and urea cycle dysregulation within the context of ccRCC. The results demonstrate that TRIM21 is downregulated in ccRCC, and low expression of TRIM21 predicts an unfavorable prognosis for ccRCC patients. Furthermore, the upregulation of TRIM21 can inhibit the migration and invasion of ccRCC cells by regulating the ubiquitination modification of ASS1. This not only expands the functional role of TRIM21 in ccRCC tumorigenesis but also demonstrates its ability to reverse urea cycle dysregulation through stabilizing ASS1 expression. Specifically, abnormal downregulation of TRIM21 in ccRCC reduces K63 ubiquitination modification of ASS1, leading to decreased stability of the ASS1 protein, aggravated urea cycle dysregulation, and facilitated migration and invasion of ccRCC cells. Additionally, reduction in ASS1 reverses the depressed migration and invasion caused by overexpression of TRIM21 in ccRCC cells. In summary, our findings contribute to a deeper understanding of the functional role played by TRIM21 in ccRCC progression, pinpoint a unique and novel regulatory mechanism involving ectopic downregulation-mediated ASS1 ubiquitination modification and urea cycle dysfunction during ccRCC progression, and provide fresh insights for further investigation into the pathogenesis and metabolic reprogramming associated with ccRCC.

Development and evaluation of suitable reference genes for qRT-PCR normalization of hybrids derived from Lycium barbarum and Lycium ruthenicum

BACKGROUND

A correct and stably expressing reference gene is prerequisite for successful quantitative real-time PCR (qRT-PCR). Investigating gene expression profiling during flower development could enhance our understanding of the molecular mechanisms of flower formation and fertility in Lycium.

METHODS AND RESULTS

In this study, 11 candidate reference genes in Lycium flower development were selected from transcriptome sequence data and evaluated with five traditional housekeeping genes from previous studies based on qRT-PCR amplification. Comparing the expression stability result of 16 candidate genes using GeNorm, NormFinder, BestKeeper, and Delta Ct algorithms, Lba04g01649 and Lba12g02820 were validated as the optimal reference genes for the flower development of Lycium.

CONCLUSIONS

The reference genes identified in this study would improve the accuracy of qRT-PCR quantification of target gene expression in Lycium flower development and facilitate future functional genomics studies on flower development. This research could lay the foundation for the study of the reproduction and development of the Lycium flower.

Analysis of oral lichen planus severity on micro-RNA linked with malignant transformation risks

OBJECTIVE

The present study evaluated the oral tissue expression of micro-RNA (miRNAs) linked to the potential malignant evolution of oral lichen planus (OLP). Furthermore, the correlation between OLP severity and miRNAs expression was assessed, and possible predictors of miRNAs in OLP patients were identified.

METHODS

The present study enrolled 41 patients with OLP (median age 58 years) and 42 healthy controls (median age 59 years). In each patient, miRNA levels (miR-7a-3p,-7a2-3p,-7a-5p,-21-3p,-21-5p,-100-3p,-100-5p,-125b-2-3p,-125b-5p,-200b-3p,-200b-5p) were assessed and analyzed through reverse transcription polymerase chain reaction. Clinical parameters and the eventual presence of OLP symptoms, signs, and disease severity scores in each patient were reported using an anamnestic questionnaire.

RESULTS

In comparison with healthy controls, OLP patients showed significantly higher miR-7a-3p,-7a-2-3p,-21-3p, miR-21-5p and miR-100-5p levels (p < 0.05) and significantly lower miR-125b-2-3p,-125b-5p,-200b-3p, and -200b-5p levels (p < 0.05). Furthermore, OLP symptoms and signs and disease severity scores were significantly correlated and were also predictors of all analyzed miRNAs (p < 0.05).

CONCLUSIONS

In comparison with healthy subjects, OLP patients exhibited unbalanced oral miRNAs expression linked to the risk of potential malignant evolution of OLP. Furthermore, some miRNAs were correlated with OLP extent and were significant predictors of OLP symptoms, signs, and disease severity scores.

A protocol for identifying universal reference genes within a genus based on RNA-Seq data: a case study of poplar stem gene expression

Real-time quantitative reverse transcription polymerase chain reaction (RT-qPCR) plays a crucial role in relative gene expression analysis, and accurate normalization relies on suitable reference genes (RGs). In this study, a pipeline for identifying candidate RGs from publicly available stem-related RNA-Seq data of different Populus species under various developmental and abiotic stress conditions is presented. DESeq2's median of ratios yielded the smallest coefficient of variance (CV) values in a total of 292 RNA-Seq samples and was therefore chosen as the method for sample normalization. A total of 541 stably expressed genes were retrieved based on the CV values with a cutoff of 0.3. Universal gene-specific primer pairs were designed based on the consensus sequences of the orthologous genes of each Populus RG candidate. The expression levels of 12 candidate RGs and six reported RGs in stems under different abiotic stress conditions or in different Populus species were assessed by RT-qPCR. The expression stability of selected genes was further evaluated using ΔCt, geNorm, NormFinder, and BestKeeper. All candidate RGs were stably expressed in different experiments and conditions in Populus. A test dataset containing 117 RNA-Seq samples was then used to confirm the expression stability, six candidate RGs and three reported RGs met the requirement of CV ≤ 0.3. In summary, this study was to propose a systematic and optimized protocol for the identification of constitutively and stably expressed genes based on RNA-Seq data, and Potri.001G349400 (CNOT2) was identified as the best candidate RG suitable for gene expression studies in poplar stems.

PagEXPA1 combines with PagCDKB2;1 to regulate plant growth and the elongation of fibers in Populus alba × Populus glandulosa

Expansins are important plant cell wall proteins. They can loosen and soften the cell walls and lead to wall extension and cell expansion. To investigate their role in wood formation and fiber elongation, the PagEXPA1 that highly expressed in cell differentiation and expansion tissues was cloned from 84K poplar (Populus alba × P. glandulosa). The subcellular localization showed that PagEXPA1 located in the cell wall and it was highly expressed in primary stems and young leaves. Compared with non-transgenic 84K poplar, overexpression of PagEXPA1 can promote plant-growth, lignification, and fiber cell elongation, while PagEXPA1 Cas9-editing mutant lines exhibited the opposite phenotype. Transcriptome analysis revealed that DEGs were mainly enriched in some important processes, which are associated with cell wall formation and cellulose synthesis. The protein interaction prediction and expression analysis showed that PagCDKB2:1 and PagEXPA1 might have an interaction relationship. The luciferase complementary assay and bimolecular fluorescence complementary assay validated that PagEXPA1 can combined with PagCDKB2;1. So they promoted the expansion of xylem vascular tissues and the development of poplar though participating in the regulation of cell division and differentiation by programming the cell-cycle. It provides good foundation for molecular breeding of fast-growing and high-quality poplar varieties.

Genome-Wide Identification and Expression Profiling of Velvet Complex Transcription Factors in Populus alba × Populus glandulosa

The discovery of new genes with novel functions is a major driver of adaptive evolutionary innovation in plants. Especially in woody plants, due to genome expansion, new genes evolve to regulate the processes of growth and development. In this study, we characterized the unique VeA transcription factor family in Populus alba × Populus glandulosa, which is associated with secondary metabolism. Twenty VeA genes were characterized systematically on their phylogeny, genomic distribution, gene structure and conserved motif, promoter binding site, and expression profiling. Furthermore, through ChIP-qPCR, Y1H, and effector-reporter assays, it was demonstrated that PagMYB128 directly regulated PagVeA3 to influence the biosynthesis of secondary metabolites. These results provide a basis for further elucidating the function of VeAs gene in poplar and its genetic regulation mechanism.

Growth-regulating factor 15-mediated vascular cambium differentiation positively regulates wood formation in hybrid poplar (Populus alba × P. glandulosa)

Introduction

Hybrid poplars are industrial trees in China. An understanding of the molecular mechanism underlying wood formation in hybrid poplars is necessary for molecular breeding. Although the division and differentiation of vascular cambial cells is important for secondary growth and wood formation, the regulation of this process is largely unclear.

Methods

In this study, mPagGRF15 OE and PagGRF15-SRDX transgenic poplars were generated to investigate the function of PagGRF15. RNA-seq and qRT-PCR were conducted to analyze genome-wide gene expression, while ChIP‒seq and ChIP-PCR were used to identified the downstream genes regulated by PagGRF15.

Results and discussion

We report that PagGRF15 from hybrid poplar (Populus alba × P. glandulosa), a growth-regulating factor, plays a critical role in the regulation of vascular cambium activity. PagGRF15 was expressed predominantly in the cambial zone of vascular tissue. Overexpression of mPagGRF15 (the mutated version of GRF15 in the miR396 target sequence) in Populus led to decreased plant height and internode number. Further stem cross sections showed that the mPagGRF15 OE plants exhibited significant changes in vascular pattern with an increase in xylem and a reduction in phloem. In addition, cambium cell files were decreased in the mPagGRF15 OE plants. However, dominant suppression of the downstream genes of PagGRF15 using PagGRF15-SRDX showed an opposite phenotype. Based on the RNA-seq and ChIP-seq results, combining qRT-PCR and ChIP-PCR analysis, candidate genes, such as WOX4b, PXY and GID1.3, were obtained and found to be mainly involved in cambial activity and xylem differentiation. Accordingly, we speculated that PagGRF15 functions as a positive regulator mediating xylem differentiation by repressing the expression of the WOX4a and PXY genes to set the pace of cambial activity. In contrast, PagGRF15 mediated the GA signaling pathway by upregulating GID1.3 expression to stimulate xylem differentiation. This study provides valuable information for further studies on vascular cambium differentiation mechanisms and genetic improvement of the specific gravity of wood in hybrid poplars.

Reference gene selection for quantitative real-time PCR analysis of Hymenopellis radicata under abiotic stress

Hymenopellis radicata (H. radicata) is an edible fungus rich in protein and mineral elements, with high edible and medical value. And reference genes suitable for normalization of qRT-PCR data from this species have not been investigated. In this study, therefore, we selected 11 housekeeping genes common in biology. The expression levels of these housekeeping genes were measured in three different tissues and six different abiotic stress treatments in mycelium. They were evaluated for expression stability using online tools. The results showed that gene ACT could be stable expressed in all samples. The expressions of genes TUB and UBQ10 are the most stable under heat stress, ACT and EF are the most stable genes under salt stress, ACT and TUB are the most stable genes under oxidation stress, RPL6 and EF are the most stable genes under pH condition, ACT and RPB2 are the most stable genes under cadmium stress, and RPB2 and UBC are the most stable genes under drought condition. ACT and PP2A are the most stable genes at different tissue sites. This study is of great help to explore the gene expression pattern of H. radicata, and also provides reference for internal reference gene screening under other conditions.

Screening and validation of optimal miRNA reference genes in different developing stages and tissues of Lilium henryi Baker

Dynamic miRNA detection using the qRT-PCR technique requires appropriate reference genes to ensure data reliability. Previous studies have screened internal reference genes in plants during embryonic development and various stress treatment, involving relatively few tissues and organs. There is no relevant miRNA study in Lilium henryi Baker and limited research on the optimal miRNA reference genes in lilies, such as 5S, 18S, U6 and Actin. Twelve genes were selected as candidate reference genes whose expression stability was analyzed in petals at different developmental stages and other tissues using various algorithms, such as geNorm, NormFinder, BestKeeper, and Delta CT. The results revealed that the optimal combination of reference genes for Lilium henryi Baker petals at different developmental stages was osa-miR166m and osa-miR166a-3p, while that for different tissues of Lilium henryi Baker was osa-miR166g-3p and osa-miR166a-3p.Four important genes related to growth and development regulation, namely, osa-miR156a, osa-miR395b, osa-miR396a-3p, and osa-miR396a-5p, were selected for validation. The findings of the present study could contribute to future investigations onmiRNA expression and the related functions in Lilium henryi Baker while providing important references for the normalization of the miRNA expression in other varieties of lily.

Integrative analysis of negatively regulated miRNA-mRNA axes for esophageal squamous cell carcinoma

BACKGROUND

MicroRNAs regulating mRNA expression by targeting at mRNAs is known constructive in tumor occurrence, immune escape, and metastasis.

OBJECTIVE

This research aims at finding negatively regulatory miRNA-mRNA pairs in esophageal squamous cell carcinoma (ESCC).

METHODS

GENE expression data of The Cancer Genome Atlas (TCGA) and GEO database were employed in differently expressed RNA and miRNA (DE-miRNAs/DE-mRNAs) screening. Function analysis was conducted with DAVID-mirPath. MiRNA-mRNA axes were identified by MiRTarBase and TarBase and verified in esophageal specimen by real-time reverse transcription polymerase chain reaction (RT-qPCR). Receiver operation characteristic (ROC) curve and Decision Curve Analysis (DCA) were applied in miRNA-mRNA pairs predictive value estimation. Interactions between miRNA-mRNA regulatory pairs and immune features were analyzed using CIBERSORT.

RESULTS

Combining TCGA database, 4 miRNA and 10 mRNA GEO datasets, totally 26 DE-miRNAs (13 up and 13 down) and 114 DE-mRNAs (64 up and 50 down) were considered significant. MiRTarBase and TarBase identified 37 reverse regulation miRNA-mRNA pairs, 14 of which had been observed in esophageal tissue or cell line. Through analysis of RT-qPCR outcome, miR-106b-5p/KIAA0232 signature was chosen as characteristic pair of ESCC. ROC and DCA verified the predictive value of model containing miRNA-mRNA axis in ESCC. Via affecting mast cells, miR-106b-5p/KIAA0232 may contribute to tumor microenvironment.

CONCLUSIONS

The diagnostic model of miRNA-mRNA pair in ESCC was established. Their complex role in ESCC pathogenesis especially tumor immunity was partly disclosed.

MicroRNA miR171b Positively Regulates Resistance to Huanglongbing of Citrus

Huanglongbing (HLB) is one of the most severe citrus diseases in the world, causing huge economic losses. However, efficient methods of protecting citrus from HLB have not yet been developed. microRNA (miRNA)-mediated regulation of gene expression is a useful tool to control plant diseases, but the miRNAs involved in regulating resistance to HLB have not yet been identified. In this study, we found that miR171b positively regulated resistance to HLB in citrus. Upon infection with HLB bacteria, the bacteria were detected in the second month in the control plants. However, in the miR171b-overexpressing transgenic citrus plants, the bacteria could not be detected until the 24th month. RNA-seq data indicated that multiple pathways, such as photosynthesis, plant–pathogen interaction, the MAPK signaling pathway, etc., might be involved in improving the resistance to HLB in miR171b-overexpressing plants compared with the control. Finally, we determined that miR171b could target SCARECROW-like (SCL) genes to downregulate their expression, which then led to promoted resistance to HLB stress. Collectively, our results demonstrate that miR171b plays a positive regulatory role in resistance to citrus HLB, and provides a new insight into the role of miRNAs in the adaptation of citrus to HLB stress.

Selection of reference genes for RT-qPCR analysis of rice with Rhizoctonia solani infection and biocontrol PGPR/KSi application

BACKGROUND

Rhizoctonia solani (AG1 IA) is an important pathogen of rice (Oryza sativa L.) that causes rice sheath blight (RSB). Since control of RSB by breeding and fungicides have had limited success, novel strategies like biocontrol with plant growth-promoting rhizobacteria (PGPR) can be an effective alternative.

METHOD AND RESULTS

Seven commonly used reference genes (RGs), 18SrRNA, ACT1, GAPDH2, UBC5, RPS27, eIF4a and CYP28, were evaluated for their stability in rice-R. solani-PGPR interaction for real-time quantitative PCR (RT-qPCR) analysis. Different algorithms were examined, Delta Ct, geNorm, NormFinder, BestKeeper, and comprehensive ranking by RefFinder, to evaluate RT-qPCR of rice in tissues infected with R. solani and treated with the PGPR strains, Pseudomonas saponiphilia and Pseudomonas protegens, with potassium silicate (KSi) alone or in combination with each PGPR strain. RG stability was affected for each treatment and treatment-specific RG selection was suggested. Validation analysis was done for nonexpressor of PR-1(NPR1) for each treatment.

CONCLUSION

Overall, ACT1 was the most stable RG with R. solani infection alone, GAPDH2 with R. solani infection plus KSi, UBC5 with R. solani infection plus P. saponiphilia, and eIF4a with R. solani infection plus P. protegens. Both ACT1 and RPS27 were the most stable with the combination of KSi and P. saponiphilia, while RPS27 was the most stable with the combination of KSi and P. protegens.

Identification and function of miRNA-mRNA interaction pairs during lateral root development of hemi-parasitic Santalum album L. seedlings

Sandalwood (Santalum album L.) is a hemi-parasitic tree species famous for its santalol and santalene, which are extracted from its heartwood and roots. The ability to understand root functionality within its branched root system would benefit the regulation of sandalwood growth and enhance the commercial value of sandalwood. Phenotypic and anatomical evidence in this study indicated that seed germination stage 4 (SG4) seemed pivotal for lateral root (LR) morphogenesis. Small RNA (sRNA) high-throughput sequencing of root tissues at three sub-stages of SG4 (lateral root primordia initiation (LRPI), lateral root primordia development (LRPD), and lateral root primordia emergence (LRPE)) was performed to identify microRNAs (miRNAs) associated with LR development. A total of 135 miRNAs, including 70 differentially expressed miRNAs (DEMs), were screened. Ten DEMs were selected to investigate transcript abundance in different organs or developmental stages. Among 100 negative DEM-mRNA interaction pairs, four targets (Sa-miR166m_2, 408d, 858a, and novel_Sa-miR8) were selected for studying cleavage sites by 5' RLM-RACE validation. The expression mode of the four miRNA-mRNA pairs was investigated after indole-3-acetic acid (IAA) treatment. IAA enhanced the abundance of homeobox-leucine-zipper protein 32 (HOX32), laccase 12 (LAC12), myeloblastosis86 (MYB86), and pectin methylesterase inhibitor6 (PMEI6) target transcripts by reducing the expression of Sa-miR166m_2, 408d, 858a, and novel_Sa-miR8 in the first 10 min. A schematic model of miRNA-regulated LR development is proposed for this hemi-parasitic species. This novel genetic information for improving sandalwood root growth and development may allow for the cultivation of fast-growing and high-yielding plantations.

Analysis of aberrant miRNA-mRNA interaction networks in prostate cancer to conjecture its molecular mechanisms

BACKGROUND

MicroRNAs (miRNAs) capable of post-transcriptionally regulating mRNA expression are essential to tumor occurrence and progression.

OBJECTIVE

This study aims to find negatively regulatory miRNA-mRNA pairs in prostate adenocarcinoma (PRAD).

METHODS

Combining The Cancer Genome Atlas (TCGA) RNA-Seq data with Gene Expression Omnibus (GEO) mRNA/miRNA expression profiles, differently expressed miRNA/mRNA (DE-miRNAs/DE-mRNAs) were identified. MiRNA-mRNA pairs were screened by miRTarBase and TarBase, databases collecting experimentally confirmed miRNA-mRNA pairs, and verified in 30 paired prostate specimens by real-time reverse transcription polymerase chain reaction (RT-qPCR). The diagnostic values of miRNA-mRNA pairs were measured by receiver operation characteristic (ROC) curve and Decision Curve Analysis (DCA). DAVID-mirPath database and Connectivity Map were employed in GO/KEGG analysis and compounds research. Interactions between miRNA-mRNA pairs and phenotypic features were analyzed with correlation heatmap in hiplot.

RESULTS

Based on TCGA RNA-Seq data, 22 miRNA and 14 mRNA GEO datasets, 67 (20 down and 47 up) miRNAs and 351 (139 up and 212 down) mRNAs were selected. After screening from 2 databases, 8 miRNA (up)-mRNA (down) and 7 miRNA (down)-mRNA (up) pairs were identified with Pearson's correlation in TCGA. By external validation, miR-221-3p (down)/GALNT3 (up) and miR-20a-5p (up)/FRMD6 (down) were chosen. The model combing 4 signatures possessed better diagnostic value. These two miRNA-mRNA pairs were significantly connected with immune cells fraction and tumor immune microenvironment.

CONCLUSIONS

The diagnostic model containing 2 negatively regulatory miRNA-mRNA pairs was established to distinguish PRADs from normal controls.

Selection and validation of appropriate reference genes for RT-qPCR analysis of Nitraria sibirica under various abiotic stresses

BACKGROUND

Nitraria sibirica Pall. is a halophytic shrub with strong environmental adaptability that can survive in extremely saline-alkali and drought-impacted environments. Gene expression analysis aids in the exploration of the molecular mechanisms of plant responses to abiotic stresses. RT-qPCR is the most common technique for studying gene expression. Stable reference genes are a prerequisite for obtaining accurate target gene expression results in RT-qPCR analysis.

RESULTS

In this study, a total of 10 candidate reference genes were selected from the transcriptome of N. sibirica, and their expression stability in leaves and roots under different treatment conditions (salt, alkali, drought, cold, heat and ABA) was evaluated with the geNorm, NormFinder, BestKeeper, comparative ΔCt and RefFinder programs. The results showed that the expression stability of the candidate reference genes was dependent on the tissue and experimental conditions tested. ACT7 combined with R3H, GAPDH, TUB or His were the most stable reference genes in the salt- or alkali-treated leaves, salt-treated roots and drought-treated roots, respectively; R3H and GAPDH were the most suitable combination for drought-treated leaves, heat-treated root samples and ABA-treated leaves; DIM1 and His maintained stable expression in roots under alkali stress; and TUB combined with R3H was stable in ABA-treated roots. TBCB and GAPDH exhibited stable expression in heat-treated leaves; TBCB, R3H, and ERF3A were stable in cold-treated leaves; and the three most stable reference genes for cold-treated roots were TBCB, ACT11 and DIM1. The reliability of the selected reference genes was further confirmed by evaluating the expression patterns of the NsP5CS gene under the six treatment conditions.

CONCLUSION

This study provides a theoretical reference for N. sibirica gene expression standardization and quantification under various abiotic stress conditions and will help to reveal the molecular mechanisms that confer stress tolerance to N. sibirica.

Multiple reaction monitoring mass spectrometry for the discovery of environmentally modulated proteins in an aquatic invertebrate sentinel species, Gammarus fossarum

Multiple reaction monitoring (MRM) mass spectrometry is emerging as a relevant tool for measuring customized molecular markers in freshwater sentinel species. While this technique is typically used for the validation of protein molecular markers preselected from shotgun experiments, recent gains of MRM multiplexing capacity offer new possibilities to conduct large-scale screening of animal proteomes. By combining the strength of active biomonitoring strategies and MRM technologies, this study aims to propose a new strategy for the discovery of candidate proteins that respond to environmental variability. For this purpose, 249 peptides derived from 147 proteins were monitored by MRM in 273 male gammarids caged in 56 environmental sites, representative of the diversity of French water bodies. A methodology is here proposed to identify a set of customized housekeeping peptides (HKPs) used to correct analytical batch effects and allow proper comparison of peptide levels in gammarids. A comparative analysis performed on HKPs-normalized data resulted in the identification of peptides highly modulated in the environment and derived from proteins likely involved in the environmental stress response. Overall, this study proposes a breakthrough approach to screen and identify potential proteins responding to relevant environmental conditions in sentinel species.

Identification of Putative Candidate Genes from Galphimia spp. Encoding Enzymes of the Galphimines Triterpenoids Synthesis Pathway with Anxiolytic and Sedative Effects

Galphimia spp. is popularly used in Mexican traditional medicine. Some populations of Galphimia exert anxiolytic and sedative effects due to the presence of the modified triterpenoids galphimines. However, the galphimine synthesis pathway has not yet been elucidated. Hence, in this study, a comparative transcriptome analysis between two contrasting populations of Galphimia spp., a galphimine-producer, and a non-galphimine-producer, is performed using RNA-Seq in the Illumina Next Seq 550 platform to identify putative candidates genes that encode enzymes of this metabolic pathway. Transcriptome functional annotation was performed using the Blast2GO in levels of gene ontology. For differential expression analysis, edgeR, pheatmap, and Genie3 library were used. To validate transcriptome data, qPCR was conducted. In producer and non-producer plants of both populations of Galphimia spp., most of the transcripts were grouped in the Molecular Function level of gene ontology. A total of 680 differentially expressed transcripts between producer and non-producer plants were detected. In galphimine-producer plants, a larger number of highly expressed transcripts related to acyclic and polycyclic terpene synthesis were identified. As putative candidate genes involved in the galphimine synthesis pathway, P450 family members and enzymes with kinase activity were identified.

Identification of suitable reference genes for quantitative reverse transcription PCR in Luffa (Luffa cylindrica)

Reverse transcription real-time quantitative PCR is widely used to quantify gene expression. Reference genes are usually used as internal controls to measure the target gene expression level. To date, there is no consensus on the use of systematically validated reference genes in different tissues of Luffa. This study evaluated the expression stability of 11 candidate reference genes in different tissues using five algorithms (BestKeeper, comparative delta-Ct method, GeNorm, NormFinder, and RefFinder). Protein phosphatase 2A was the most stable gene, while alpha Tubulin was the least stable. The relative expression of ethylene-related genes in different tissues was also analyzed to reveal their role in sex determination. This study provides the basis for using suitable reference genes to evaluate targeted gene expression.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01182-8.

Identification and validation of reference genes for qRT-PCR based studies in horse gram (Macrotyloma uniflorum)

The quantitative real-time polymerase chain reaction (qRT-PCR) is the most sensitive and commonly used technique for gene expression studies in biological systems. However, the reliability of qRT-PCR results depends on the selection of reference gene(s) for data normalization. Horse gram (Macrotyloma uniflorum) is an important legume crop on which several molecular studies have been reported. However, the stability of reference genes has not been evaluated. In the present study, nine candidate reference genes were identified from horse gram RNA-seq data and evaluated in two horse gram genotypes, HPK4 and HPKM317 under six abiotic stresses viz. cold, drought, salinity, heat, abscisic acid and methyl viologen-induced oxidative stress. The results were evaluated using geNorm, Bestkeeper, Normfinder and delta-delta Ct methods and comprehensive ranking was assigned using RefFinder and RankAggreg software. The overall result showed that TCTP was one of the most stable genes in all samples and in genotype HPK4, while in HPKM317 profilin was most stably expressed. However, PSMA5 was identified as least stable in all the experimental conditions. Expression of target genes dehydrin and early response to dehydration 6 under drought stress was also validated using TCTP and profilin for data normalization, either alone or in combination, which confirmed their suitability for qRT-PCR data normalization. Thus, TCTP and profilin genes may be used for qRT-PCR data normalization for molecular and genomic studies in horse gram.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01104-0.

Efficient knockout of the phytoene desaturase gene in a hybrid poplar (Populus alba × Populus glandulosa) using the CRISPR/Cas9 system with a single gRNA

The CRISPR/Cas9 system has been used for genome editing in several plant species; however, there are few reports on its use in trees. Here, CRISPR/Cas9 was used to mutate a target gene in Populus alba × Populus glandulosa hybrid poplars. The hybrid poplar is routinely used in molecular biological studies due to the well-established Agrobacterium-mediated transformation method. A single guide RNA (sgRNA) with reported high mutation efficiency in other popular species was designed with a protospacer adjacent motif sequence for the phytoene desaturase 1 (PagPDS1) gene. The pHSE/Cas9-PagPDS1 sgRNA vector was delivered into hybrid poplar cells using Agrobacterium-mediated transformation. The transgenic plants were propagated and classified them into three groups according to their phenotypes. Among a total of 110 lines of transgenic hybrid poplars, 82 lines showed either an albino or a pale green phenotype, indicating around 74.5% phenotypic mutation efficiency of the PagPDS1 gene. The albino phenotypes were observed when the CRISPR/Cas9-mediated mutations in both PagPDS1 alleles in the transgenic plants. There was no off-target modification of the PagPDS2 gene, which has a potential sgRNA target sequence with two mismatches. The results confirmed that the sgRNA can specifically edit PagPDS1 rather than PagPDS2, indicating that CRISPR/Cas9-mediated genome editing can effectively induce target mutations in the hybrid poplar. This technique will be useful to improve tree quality in hybrid poplars (P. alba × P. glandulosa); for example, by enhancing biomass or stress tolerance.

Identification of reliable reference genes for quantitative real-time PCR analysis of the Rhus chinensis Mill. leaf response to temperature changes

Rhus chinensis Mill. (RCM) is the host plant of Galla chinensis, which is valued in traditional medicine. Environmental temperature directly determines the probability of gallnut formation and RCM growth. At present, there is no experiment to systematically analyse the stability of internal reference gene (RG) expression in RCM. In this experiment, leaves that did not form gallnuts were used as the control group, while leaves that formed gallnuts were used as the experimental group. First, we conducted transcriptome experiments on RCM leaves to obtain 45 103 differential genes and functional enrichment annotations between the two groups. On this basis, this experiment established a transcriptional gene change model of leaves in the process of gallnut formation after being bitten by aphids, and RCM reference candidate genes were screened from RNA sequencing (RNA‐seq) data. This study is based on RCM transcriptome data and evaluates the stability of 11 potential reference genes under cold stress (4 °C) and heat stress (34 °C), using three statistical algorithms (geNorm, NormFinder, and BestKeeper). The results show that GAPDH1 + PP2A2/UBQ are stable reference genes under heat stress, while GAPDH1 + ACT are the most stable under cold stress. This study is the first to screen candidate reference genes in RCM and could help guide future molecular studies in this genus.

Knockdown of miR393 Promotes the Growth and Biomass Production in Poplar

Short tandem target mimic (STTM), which is composed of two short sequences mimicking small RNA target sites, separated by a linker of optimal size, can block the functions of all members in a miRNA family. microRNA393 (miR393), which is one of the conserved miRNA families in plants, can regulate plant root growth, leaf development, plant architecture, and stress resistance. In order to verify the role of miR393 in the secondary growth of trees, we created its STTM transgenic poplar lines (STTM393). The expression of miR393 in STTM393 lines was reduced by over 10 times compared with the control plants. STTM393 lines showed promoted growth with about 20% higher, 15% thicker, and 2-4 more internodes than the control plants after 3 months of growth. The cross-section of the stems showed that STTM393 lines had wider phloem, xylem, and more cambium cell layers than control plants, and the lignin content in STTM393 lines was also higher as revealed by staining and chemical determination. Based on the transcriptome analysis, the genes related to the auxin signaling pathway, cell cyclin, cell expansion, and lignin synthesis had higher expression in STTM393 lines than that in control plants. The higher expression levels of FBL family members suggested that the auxin signaling pathway was strengthened in STTM393 lines to promote plant growth. Therefore, the knockdown of miR393 using the STTM approach provides a way to improve poplar growth and biomass production.

Selection and Optimization of Reference Genes for MicroRNA Expression Normalization by qRT-PCR in Chinese Cedar (Cryptomeria fortunei) under Multiple Stresses

MicroRNA (miRNA) expression analysis is very important for investigating its functions. To date, no research on reference genes (RGs) for miRNAs in gymnosperms, including Cryptomeria fortunei, has been reported. Here, ten miRNAs (i.e., pab-miR159a, cln-miR162, cas-miR166d, pab-miR395b, ppt-miR894, cln-miR6725, novel1, novel6, novel14 and novel16) and three common RGs (U6, 5S and 18S) were selected as candidate RGs. qRT-PCR was used to analyse their expressions in C. fortunei under various experimental conditions, including multiple stresses (cold, heat, drought, salt, abscisic acid and gibberellin) and in various tissues (roots, stems, tender needles, cones and seeds). Four algorithms (delta Ct, geNorm, NormFinder and BestKeeper) were employed to assess the stability of candidate RG expression; the geometric mean and RefFinder program were used to comprehensively evaluate RG stability. According to the results, novel16, cln-miR6725, novel1 and U6 were the most stable RGs for studying C. fortunei miRNA expression. In addition, the expression of three target miRNAs (aly-miR164c-5p, aly-miR168a-5p and smo-miR396) was examined to verify that the selected RGs are suitable for miRNA expression normalisation. This study may aid further investigations of miRNA expression/function in the response of C. fortunei to abiotic stress and provides an important basis for the standardisation of miRNA expression in other gymnosperm species.

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

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

Identification and validation of miRNA reference genes in poplar under pathogen stress

Quantitative real time polymerase chain reaction (qRT-PCR) is a common method to analyze gene expression. Due to differences in RNA quantity, quality, and reverse transcription efficiency between qRT-PCR samples, reference genes are used as internal standards to normalize gene expression. However, few universal genes, especially miRNAs, have been identified as reference so far. Therefore, it is essential to identify reference genes that can be used across various experimental conditions, stress treatments, or tissues. In this study, 14 microRNAs (miRNAs) and 5.8S rRNA were assessed for expression stability in poplar trees infected with canker pathogen. Using geNorm, NormFinder and Bestkeeper reference gene analysis programs, we found that miR156g and miR156a exhibited stable expression throughout the infection process. miR156g, miR156a and 5.8S rRNA were then tested as internal standards to measure the expression of miR1447 and miR171c, and the results were compared to small RNA sequencing (RNA-seq) data. We found that when miR156a and 5.8S rRNA were used as the reference gene, the expression of miR1447 and miR171c were consistent with the small RNA-seq expression profiles. Therefore, miR156a was the most stable miRNAs examined in this study, and could be used as a reference gene in poplar under canker pathogen stress, which should enable comprehensive comparisons of miRNAs expression and avoid the bias caused by different length between detected miRNAs and traditional reference genes. The present study has expanded the miRNA reference genes available for gene expression studies in trees under biotic stress.

Selection of reference genes for gene expression analysis in Liriodendron hybrids' somatic embryogenesis and germinative tissues

The differential expression of genes is crucial for plant somatic embryogenesis (SE), and the accurate quantification of gene expression levels relies on choosing appropriate reference genes. To select the most suitable reference genes for SE studies, 10 commonly used reference genes were examined in synchronized somatic embryogenic and subsequent germinative cultures of Liriodendron hybrids by using quantitative real-time reverse transcription PCR. Four popular normalization algorithms: geNorm, NormFinder, Bestkeeper and Delta-Ct were used to select and validate the suitable reference genes. The results showed that elongation factor 1-gamma, histone H1 linker protein, glyceraldehyde-3-phosphate dehydrogenase and α-tubulin were suitable for SE tissues, while elongation factor 1-gamma and actin were best for the germinative organ tissues. Our work will benefit future studies of gene expression and functional analyses of SE in Liriodendron hybrids. It is also serves as a guide of reference gene selection in early embryonic gene expression analyses for other woody plant species.

Optimization of reference genes for qRT-PCR analysis of microRNA expression under abiotic stress conditions in sweetpotato

Sweetpotato (Ipomoea batatas. L) is an important food crop, harvested for its nutrient-rich tuberous roots. Drought and salt stresses are two major factors limiting the sweetpotato production. Since microRNAs (miRNAs) are well known to play crucial roles in regulation of plant stress responses, quantitative profiling of miRNA expression under stress conditions will facilitate identification and genetic manipulation of novel miRNAs to improve stress tolerance. Real-time quantitative reverse transcription PCR (qRT-PCR) is a commonly used tool for this purpose, but not without challenges. Although stem-loop and poly(A)-tail modified qRT-PCR methods were developed for characterizing miRNA expression, accurate profiling of miRNAs is still difficult in many plant species because of a lack of reliable reference genes for normalizing miRNA transcripts. To identify reference genes that are suitable for normalizing miRNA expression in sweetpotato, the expression stability of eight candidate miRNAs and two commonly used reference genes were tested in 96 samples involving four tissues and two cultivars under drought and salt stress treatments. Data analysis using the geNorm, NormFinder and Bestkeeper algorithms demonstrated that miRn60, miR482, and their combination were reliable references. We further validated the reference genes by expression analysis of the well-characterized miR319 and miR156 that regulate drought and salt stress responses, respectively. The reference genes identified in this study will facilitate future miRNA analysis under abiotic stress conditions in sweetpotato.

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.

miR-26a-5p is a Stable Reference Gene for miRNA Studies in Chondrocytes from Developing Human Cartilage

miRNAs are emerging as key regulators of complex biological systems in several developmental processes. qRT-PCR is a powerful tool to quantitatively assess the profiles and modulation of miRNA expression. In the emerging field of cartilage maturation studies, from precursor to hypertrophic chondrocytes, few data about miRNA regulation are available, and no consensus on the best reference gene (RG) has been reached. This is a crucial pitfall since reliable outcomes depend on proper data normalization. The aim of this work was to identify reliable and stable miRNA RGs, basing the analysis on available high throughput qRT-PCR miRNA data (from the NCBI Gene Expression Omnibus database, GSE49152) obtained from human embryonic cartilage tissues enriched in the precursor, differentiated, and hypertrophic chondrocytes. Four normalization approaches were used, and the stability was quantified by combining BestKeeper, delta-Ct, geNorm, and NormFinder statistical tools. An integrated approach allowed to identify miR-26a-5p as the most stable RG and miR-212-3p as the worst one. RNU44, used in original dataset analysis, performed as second best RG. Applications of different normalization strategies significantly impacted the profiles and modulation of miRNA expression. Herein presented results point out the crucial need of a consensus on data normalization studies aimed at dissecting miRNA role in human cartilage development, to avoid the postulation of unreliable biological conclusions.