Selection and Validation of Reference Genes for qRT-PCR Analysis in the Oil-Rich Tuber Crop Tiger Nut (Cyperus esculentus) Based on Transcriptome Data

Genome-wide identification and molecular evolution of Dof transcription factors in Cyperus esculentus

Dof transcription factor family in Cyperus esculentus genome was identified and analyzed using bioinformatics. The analysis results revealed that C.esculentus genome contains 29 Dof genes (CesDof), all of which are located in the nucleus according to subcellular localization prediction. CesDof proteinrs have a range of 124 to 512 amino acids, with most being basic proteins. Their secondary structure was mainly irregular curl. The promoter sequence of CesDof genes contains cis-acting elements that respond to light, drought, hormones, low temperature, and circadian rhythm. Codon preference analysis showed that CesDof genes' codon preference ends in T/A. Collinearity analysis revealed that C.esculentus had three pairs of collinear CesDof genes. Additionally, there were 15 pairs of collinear genes between C.esculentus and Arabidopsis thaliana. The genetic relationship between C.esculentus and Rhynchospora pubera was found to be the closest. Phylogenetic tree analysis revealed that 29 CesDof genes of C.esculentus can be classified into 4 subgroups. Additionally, 144 miRNAs were predicted to target these CesDof genes. Furthermore, protein interaction analysis indicated that 15 Dof proteins in C.esculentus had interactions. The qRT-PCR verification results of drought stress and salt stress treatment experiments showed that most CesDof genes were involved in drought stress and salt stress responses, and the gene expression trends under drought stress and salt stress conditions were consistent. These results lay a theoretical foundation for further studying the molecular functions of Dof gene family in C.esculentus and its molecular mechanisms in regulating the life activities of C.esculentus.

RNA-Seq Reveals That Multiple Pathways Are Involved in Tuber Expansion in Tiger Nuts (Cyperus esculentus L.)

The tiger nut (Cyperus esculentus L.) is a usable tuber and edible oil plant. The size of the tubers is a key trait that determines the yield and the mechanical harvesting of tiger nut tubers. However, little is known about the anatomical and molecular mechanisms of tuber expansion in tiger nut plants. This study conducted anatomical and comprehensive transcriptomics analyses of tiger nut tubers at the following days after sowing: 40 d (S1); 50 d (S2); 60 d (S3); 70 d (S4); 90 d (S5); and 110 d (S6). The results showed that, at the initiation stage of a tiger nut tuber (S1), the primary thickening meristem (PTM) surrounded the periphery of the stele and was initially responsible for the proliferation of parenchyma cells of the cortex (before S1) and then the stele (S2-S3). The increase in cell size of the parenchyma cells occurred mainly from S1 to S3 in the cortex and from S3 to S4 in the stele. A total of 12,472 differentially expressed genes (DEGs) were expressed to a greater extent in the S1-S3 phase than in S4-S6 phase. DEGs related to tuber expansion were involved in cell wall modification, vesicle transport, cell membrane components, cell division, the regulation of plant hormone levels, signal transduction, and metabolism. DEGs involved in the biosynthesis and the signaling of indole-3-acetic acid (IAA) and jasmonic acid (JA) were expressed highly in S1-S3. The endogenous changes in IAA and JAs during tuber development showed that the highest concentrations were found at S1 and S1-S3, respectively. In addition, several DEGs were related to brassinosteroid (BR) signaling and the G-protein, MAPK, and ubiquitin-proteasome pathways, suggesting that these signaling pathways have roles in the tuber expansion of tiger nut. Finally, we come to the conclusion that the cortex development preceding stele development in tiger nut tubers. The auxin signaling pathway promotes the division of cortical cells, while the jasmonic acid pathway, brassinosteroid signaling, G-protein pathway, MAPK pathway, and ubiquitin protein pathway regulate cell division and the expansion of the tuber cortex and stele. This finding will facilitate searches for genes that influence tuber expansion and the regulatory networks in developing tubers.

Selection and Validation of Reference Genes for Gene Expression in Bactericera gobica Loginova under Different Insecticide Stresses

Insecticide resistance has long been a problem in crop pest control. Bactericera gobica is a major pest on the well-known medicinal plants Lycium barbarum L. Investigating insecticide resistance mechanisms of B. gobica will help to identify pesticide reduction strategies to control the pest. Gene expression normalization by RT-qPCR requires the selection and validation of appropriate reference genes (RGs). Here, 15 candidate RGs were selected from transcriptome data of B. gobica. Their expression stability was evaluated with five algorithms (Delta Ct, GeNorm, Normfinder, BestKeeper and RefFinder) for sample types differing in response to five insecticide stresses and in four other experimental conditions. Our results indicated that the RGs RPL10 + RPS15 for Imidacloprid and Abamectin; RPL10 + AK for Thiamethoxam; RPL32 + RPL10 for λ-cyhalothrin; RPL10 + RPL8 for Matrine; and EF2 + RPL32 under different insecticide stresses were the most suitable RGs for RT-qPCR normalization. EF1α + RPL8, EF1α + β-actin, β-actin + EF2 and β-actin + RPS15 were the optimal combination of RGs under odor stimulation, temperature, developmental stages and both sexes, respectively. Overall, EF2 and RPL8 were the two most stable RGs in all conditions, while α-TUB and RPL32 were the least stable RGs. The corresponding suitable RGs and one unstable RG were used to normalize a target cytochrome P450 CYP6a1 gene between adult and nymph stages and under imidacloprid stress. The results of CYP6a1 expression were consistent with transcriptome data. This study is the first research on the most stable RG selection in B. gobica nymphs exposed to different insecticides, which will contribute to further research on insecticide resistance mechanisms in B. gobica.

Selection and Validation of Reference Genes for Gene Expression Studies in Euonymus japonicus Based on RNA Sequencing

Euonymus japonicus is one of the most low-temperature-tolerant evergreen broad-leaved tree species in the world and is widely used in urban greening. However, there are very few molecular biology studies on its low-temperature tolerance mechanism. So far, no researcher has selected and reported on its reference genes. In this study, 21 candidate reference genes (12 traditional housekeeping genes and 9 other genes) were initially selected based on gene expression and coefficient of variation (CV) through RNA-Seq (unpublished data), and qRT-PCR was used to detect the expression levels of candidate reference genes in three different groups of samples (leaves under different temperature stresses, leaves of plants at different growth stages, and different organs). After further evaluating the expression stability of these genes using geNorm, NormFinder, Bestkeeper, and RefFind, the results show that the traditional housekeeping gene eIF5A and the new reference gene RTNLB1 have good stability in the three different groups of samples, so they are reference genes with universality. In addition, we used eIF5A and RTNLB1 as reference genes to calibrate the expression pattern of the target gene EjMAH1, which confirmed this view. This article is the first to select and report on the reference gene of E. japonicus, laying the foundation for its low-temperature tolerance mechanism and other molecular biology research.

Molecular characterization of oleosin genes in Cyperus esculentus, a Cyperaceae plant producing oil in underground tubers

KEY MESSAGE

CeOLE genes exhibit a tuber-predominant expression pattern and their mRNA/protein abundances are positively correlated with oil accumulation during tuber development. Overexpression could significantly increase the oil content of tobacco leaves. Oleosins (OLEs) are abundant structural proteins of lipid droplets (LDs) that function in LD formation and stabilization in seeds of oil crops. However, little information is available on their roles in vegetative tissues. In this study, we present the first genome-wide characterization of the oleosin family in tigernut (Cyperus esculentus L., Cyperaceae), a rare example accumulating high amounts of oil in underground tubers. Six members identified represent three previously defined clades (i.e. U, SL and SH) or six out of seven orthogroups (i.e. U, SL1, SL2, and SH1-3) proposed in this study. Comparative genomics analysis reveals that lineage-specific expansion of Clades SL and SH was contributed by whole-genome duplication and dispersed duplication, respectively. Moreover, presence of SL2 and SH3 in Juncus effuses implies their appearance sometime before Cyperaceae-Juncaceae divergence, whereas SH2 appears to be Cyperaceae specific. Expression analysis showed that CeOLE genes exhibit a tuber-predominant expression pattern and transcript levels are considerably more abundant than homologs in the close relative Cyperus rotundus. Moreover, CeOLE mRNA and protein abundances were shown to positively correlate with oil accumulation during tuber development. Additionally, two dominant isoforms (i.e. CeOLE2 and -5) were shown to locate in LDs as well as the endoplasmic reticulum of tobacco (Nicotiana benthamiana) leaves, and are more likely to function in homo and heteromultimers. Furthermore, overexpression of CeOLE2 and -5 in tobacco leaves could significantly increase the oil content, supporting their roles in oil accumulation. These findings provide insights into lineage-specific family evolution and putative roles of CeOLE genes in oil accumulation of vegetative tissues, which facilitate further genetic improvement for tigernut.

Selection and verification of reliable internal reference genes in stem development of herbaceous peony (Paeonia lactiflora Pall.)

Herbaceous peony (Paeonia lactiflora Pall.) has emerged in the cut flower market due to its beautiful appearance. The bending flower stems caused by a lack of mechanical strength is the main problem restricting the development of the cut P. lactiflora industry. So it is of great worth to reveal the basis of stem development changes in P. lactiflora to improve its cut flower quality. Quantitative research on gene expression characteristics can provide clues for understanding their biological functions, and the screening of relatively stable expression genes is a prerequisite for the quantitative study of gene expression characteristics. Thus, it is necessary to find appropriate genes during stem development so as to analyze the qRT‒PCR results. In this study, 10 genes were screened, and these genes expressed stably in stems of different stem strengths at three different developmental stages. Then, their expressions were evaluated by RefFinder, BestKeeper, NormFinder, and GeNorm programs. The results demonstrated that γ-tubulin (γ-TUB) was the most suitable gene, followed by α-tubulin (α-TUB) and β-D-glucosidase (β-GUS), whereas histone H3 (His) was the least suitable gene. Besides, the temporal and spatial expression characteristics of PlCOMT1, the key gene concerned with the synthesis of cell wall fillers in P. lactiflora, were also used to evaluate the suitability of genes. Consequently, γ-TUB and α-TUB are the two best combinations during stem development, and their combination can be used for the stem development of P. lactiflora. These findings will provide a reference for the selection of genes related to stem development and the study of molecular mechanisms related to stem development in P. lactiflora.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01325-5.

Analysis of oil synthesis pathway in Cyperus esculentus tubers and identification of oleosin and caleosin genes

The tubers of the widely distributed Cyperus esculentus are rich in oil, and therefore, the plant is considered to have a high utilization value in the vegetable oil industry. Oleosins and caleosins are lipid-associated proteins found in oil bodies of seeds; however oleosins and caleosins genes have not been identified in C. esculentus. In this study, we performed transcriptome sequencing and lipid metabolome analysis of C. esculentus tubers at four developmental stages to obtain the information on their genetic profile, expression trends, and metabolites in oil accumulation pathways. Overall, 120,881 non-redundant unigenes and 255 lipids were detected; 18 genes belonged to the acetyl-CoA carboxylase (ACC), malonyl-CoA:ACP transacylase (MCAT), β-ketoacyl-ACP synthase (KAS), and fatty acyl-ACP thioesterase (FAT) gene families involved in fatty acid biosynthesis, and 16 genes belonged to the glycerol-3-phosphate acyltransferase (GPAT), diacylglycerol acyltransferase 3 (DGAT3), phospholipid:diacylglycerol acyltransferase (PDAT), FAD2, and lysophosphatidic acid acyltransferase (LPAAT) gene families playing important roles in triacylglycerol synthesis. We also identified 9 oleosin- and 21 caleosin-encoding genes in C. esculentus tubers. These results provide detailed information on the C. esculentus transcriptional and metabolic profiles, which can be used as reference for the development of strategies to increase oil content in C. esculentus tubers.

Chromosome-scale Genome Assembly of the Yellow Nutsedge (Cyperus esculentus)

The yellow nutsedge (Cyperus esculentus L. 1753) is an unconventional oil plant with oil-rich tubers, and a potential alternative for traditional oil crops. Here, we reported the first high-quality and chromosome-level genome assembly of the yellow nutsedge generated by combining PacBio HiFi long reads, Novaseq short reads, and Hi-C data. The final genome size is 225.6 Mb with an N50 of 4.3 Mb. More than 222.9 Mb scaffolds were anchored to 54 pseudochromosomes with a BUSCO score of 96.0%. We identified 76.5 Mb (33.9%) repetitive sequences across the genome. A total of 23,613 protein-coding genes were predicted in this genome, of which 22,847 (96.8%) were functionally annotated. A whole-genome duplication event was found after the divergence of Carex littledalei and Rhynchospora breviuscula, indicating the rich genetic resources of this species for adaptive evolution. Several significantly enriched GO terms were related to invasiveness of the yellow nutsedge, which may explain its plastic adaptability. In addition, several enriched Kyoto Encyclopedia of Genes and Genomes pathways and expanded gene families were closely related with substances in tubers, partially explaining the genomic basis of characteristics of this oil-rich tuber.

Evaluation of Candidate Reference Genes for Gene Expression Analysis in Wild Lamiophlomis rotata

Lamiophlomis rotata (Benth.) Kudo is a perennial and unique medicinal plant of the Qinghai–Tibet Plateau. It has the effects of diminishing inflammation, activating blood circulation, removing blood stasis, reducing swelling, and relieving pain. However, thus far, reliable reference gene identifications have not been reported in wild L. rotata. In this study, we identified suitable reference genes for the analysis of gene expression related to the medicinal compound synthesis in wild L. rotata subjected to five different-altitude habitats. Based on the RNA-Seq data of wild L. rotata from five different regions, the stability of 15 candidate internal reference genes was analyzed using geNorm, NormFinder, BestKeeper, and RefFinder. TFIIS, EF-1α, and CYP22 were the most suitable internal reference genes in the leaves of L. rotata from different regions, while OBP, TFIIS, and CYP22 were the optimal reference genes in the roots of L. rotata. The reference genes identified here would be very useful for gene expression studies with different tissues in L. rotata from different habitats.

Analysis of Lhc family genes reveals development regulation and diurnal fluctuation expression patterns in Cyperus esculentus, a Cyperaceae plant

Sixteen Lhc genes representing 13 phylogenetic groups were identified from the full-length transcriptome of tigernut, exhibiting development regulation and diurnal fluctuation expression patterns in leaves. Nuclear encoded light-harvesting chlorophyll a/b-binding (Lhc) proteins play indispensable roles in oxygenic photosynthesis. In this study, we present the first transcriptome-based characterization of Lhc family genes in tigernut (Cyperus esculentus L.), a Cyperaceae C₄ plant producing oil in underground tubers. A number of 16 Lhc genes representing 13 phylogenetic groups identified from the full-length tigernut transcriptome are equal to that found in both Carex littledalei (another Cyperaceae plant) and papaya, slightly more than 15 members present in both rice and jatropha, but relatively less than 18, 20, and 21 members present in sorghum, cassava, and Arabidopsis, respectively. Nevertheless, nearly one-vs-one orthologous relationship was observed in most groups, though some of them are no longer located in syntenic blocks and species-specific expansion was frequently found in Lhcb1. Comparative genomics analysis revealed that the loss of two groups (i.e., Lhca2 and Lhca5) in C. littledalei is species-specific, sometime after the split with tigernut, and the expansion of Lhcb1 was mainly contributed by tandem duplication as observed in most species. Interestingly, a transposed duplication, which appears to be shared by monocots, was also identified in Lhcb1. Further transcriptome profiling revealed a predominant expression pattern of most CeLhc family genes in photosynthetic tissues and enhanced transcription during leaf maturation, reflecting their key roles in light absorption. Moreover, qRT-PCR analysis revealed an apparent diurnal fluctuation expression pattern of 11 dominant CeLhc genes. These findings not only highlight species-specific evolution of Lhc genes in the Cyperaceae family as well as the monocot lineage, but also provide valuable information for further functional analysis and genetic improvement in tigernut.

Transcriptome-Based Selection and Validation of Reference Genes for Gene Expression in Goji Fruit Fly (Neoceratitis asiatica Becker) under Developmental Stages and Five Abiotic Stresses

Goji fruit fly, Neoceratitis asiatica, is a major pest on the well-known medicinal plant Lycium barbarum. Dissecting molecular mechanisms of infestation and host selection of N. asiatica will contribute to the determination of best management practices for pest fly control. Gene expression normalization by Real-time quantitative PCR (qPCR) requires the selection and validation of appropriate reference genes (RGs). Hence, 15 candidate RGs were selected from transcriptome data of N. asiatica. Their expression stability was evaluated with five algorithms (∆Ct, Normfinder, GeNorm, BestKeeper, and RefFinder) for sample types differing in the developmental stage, sex, tissue type, and in response to five different abiotic stresses. Our results indicated that the RGs β-Actin + GST for sex, RPL32 + EF1α for tissue type, RPS13+ EF1α for developmental stages along with odor stimulation, color induction, and starvation-refeeding stresses, EF1α + GAPDH under insecticide stress, RPS13 + RPS18 under temperature stress, respectively, were selected as the most suitable RGs for qPCR normalization. Overall, RPS18 and EF1α were the two most stable RGs in all conditions, while RPS15 and EF1β were the least stable RGs. The corresponding suitable RGs and one unstable RG were used to normalize a target odorant-binding protein OBP56a gene in male and female antennae, different tissues, and under odor stimulation. The results of OBP56a expression were consistent with transcriptome data. Our study is the first research on the most stable RGs selection in N. asiatica, which will facilitate further studies on the mechanisms of host selection and insecticide resistance in N. asiatica.