Validating Internal Control Genes for the Accurate Normalization of qPCR Expression Analysis of the Novel Model Plant Setaria viridis

GOLDEN2-like1 is sufficient but not necessary for chloroplast biogenesis in mesophyll cells of C4 grasses

Chloroplasts are the site of photosynthesis. In land plants, chloroplast biogenesis is regulated by a family of transcription factors named GOLDEN2-like (GLK). In C4 grasses, it has been hypothesized that genome duplication events led to the sub-functionalization of GLK paralogs (GLK1 and GLK2) to control chloroplast biogenesis in two distinct cell types: mesophyll and bundle sheath cells. Although previous characterization of golden2 (g2) mutants in maize has demonstrated a role for GLK2 paralogs in regulating chloroplast biogenesis in bundle sheath cells, the function of GLK1 has remained elusive. Here we show that, contrary to expectations, GLK1 is not required for chloroplast biogenesis in mesophyll cells of maize. Comparisons between maize and Setaria viridis, which represent two independent C4 origins within the Poales, further show that the role of GLK paralogs in controlling chloroplast biogenesis in mesophyll and bundle sheath cells differs between species. Despite these differences, complementation analysis revealed that GLK1 and GLK2 genes from maize are both sufficient to restore functional chloroplast development in mesophyll and bundle sheath cells of S. viridis mutants. Collectively our results suggest an evolutionary trajectory in C4 grasses whereby both orthologs retained the ability to induce chloroplast biogenesis but GLK2 adopted a more prominent developmental role, particularly in relation to chloroplast activation in bundle sheath cells.

Domestication-associated PHYTOCHROME C is a flowering time repressor and a key factor determining Setaria as a short-day plant

Phytochromes play vital roles in the regulation of flowering time, but little is known in Panicoideae species, especially the C₄ model Setaria. Here, genomic variations of PHYTOCHROME C (PHYC) between wild and cultivated Setaria gene pools were analysed and three SiphyC mutants were identified. The function of SiPHYC was verified by CRISPR-Cas9 approach and transcriptome sequencing. Furthermore, efficiency of indoor cultivation of SiphyC mutants were systematically evaluated. An extreme purified selection of PHYC was detected in wild to cultivated domestication process of Setaria. SiphyC mutants and knockout transgenic plants showed an early heading date and a loss of response to short-day photoperiod. Furthermore, variable expression of SiFTa, SiMADS14 and SiMADS15 might be responsible for promoting flowering of SiphyC mutants. Moreover, SiphyC mutant was four times that of the indoor plot ratio of wild-type and produced over 200 seeds within 45 d per individual. Our results suggest that domestication-associated SiPHYC repressed flowering and determined Setaria as a short-day plant, and SiphyC mutants possess the potential for creating efficient indoor cultivation system suitable for research on Setaria as a model, and either for maize or sorghum as well.

Multi-omics intervention in Setaria to dissect climate-resilient traits: Progress and prospects

Millets constitute a significant proportion of underutilized grasses and are well known for their climate resilience as well as excellent nutritional profiles. Among millets, foxtail millet (Setaria italica) and its wild relative green foxtail (S. viridis) are collectively regarded as models for studying broad-spectrum traits, including abiotic stress tolerance, C₄ photosynthesis, biofuel, and nutritional traits. Since the genome sequence release, the crop has seen an exponential increase in omics studies to dissect agronomic, nutritional, biofuel, and climate-resilience traits. These studies have provided first-hand information on the structure, organization, evolution, and expression of several genes; however, knowledge of the precise roles of such genes and their products remains elusive. Several open-access databases have also been instituted to enable advanced scientific research on these important crops. In this context, the current review enumerates the contemporary trend of research on understanding the climate resilience and other essential traits in Setaria, the knowledge gap, and how the information could be translated for the crop improvement of related millets, biofuel crops, and cereals. Also, the review provides a roadmap for studying other underutilized crop species using Setaria as a model.

Regulatory mechanisms behind the phenotypic plasticity associated with Setaria italica water deficit tolerance

Drought is one of the main environmental stresses that negatively impacts vegetative and reproductive yield. Water deficit responses are determined by the duration and intensity of the stress, which, together with plant genotype, will define the chances of plant survival. The metabolic adjustments in response to water deficit are complex and involve gene expression modulation regulated by DNA-binding proteins and epigenetic modifications. This last mechanism may also regulate the activity of transposable elements, which in turn impact the expression of nearby loci. Setaria italica plants submitted to five water deficit regimes were analyzed through a phenotypical approach, including growth, physiological, RNA-seq and sRNA-seq analyses. The results showed a progressive reduction in yield as a function of water deficit intensity associated with signaling pathway modulation and metabolic adjustments. We identified a group of loci that were consistently associated with drought responses, some of which were related to water deficit perception, signaling and regulation. Finally, an analysis of the transcriptome and sRNAome allowed us to identify genes putatively regulated by TE- and sRNA-related mechanisms and an intriguing positive correlation between transcript levels and sRNA accumulation in gene body regions. These findings shed light on the processes that allow S. italica to overcome drought and survive under water restrictive conditions.

A low CO2-responsive mutant of Setaria viridis reveals that reduced carbonic anhydrase limits C4 photosynthesis

In C4 species, β-carbonic anhydrase (CA), localized to the cytosol of the mesophyll cells, accelerates the interconversion of CO2 to HCO3-, the substrate used by phosphoenolpyruvate carboxylase (PEPC) in the first step of C4 photosynthesis. Here we describe the identification and characterization of low CO2-responsive mutant 1 (lcr1) isolated from an N-nitroso-N-methylurea- (NMU) treated Setaria viridis mutant population. Forward genetic investigation revealed that the mutated gene Sevir.5G247800 of lcr1 possessed a single nucleotide transition from cytosine to thymine in a β-CA gene causing an amino acid change from leucine to phenylalanine. This resulted in severe reduction in growth and photosynthesis in the mutant. Both the CO2 compensation point and carbon isotope discrimination values of the mutant were significantly increased. Growth of the mutants was stunted when grown under ambient pCO2 but recovered at elevated pCO2. Further bioinformatics analyses revealed that the mutation has led to functional changes in one of the conserved residues of the protein, situated near the catalytic site. CA transcript accumulation in the mutant was 80% lower, CA protein accumulation 30% lower, and CA activity ~98% lower compared with the wild type. Changes in the abundance of other primary C4 pathway enzymes were observed; accumulation of PEPC protein was significantly increased and accumulation of malate dehydrogenase and malic enzyme decreased. The reduction of CA protein activity and abundance in lcr1 restricts the supply of bicarbonate to PEPC, limiting C4 photosynthesis and growth. This study establishes Sevir.5G247800 as the major CA allele in Setaria for C4 photosynthesis and provides important insights into the function of CA in C4 photosynthesis that would be required to generate a rice plant with a functional C4 biochemical pathway.

Systematic selection and validation of suitable reference genes for quantitative real-time PCR normalization studies of gene expression in Nitraria tangutorum

Suitable reference genes can be used to calibrate the error in quantitative real-time PCR (qPCR) experiments, making the results more credible. However, there are no reference genes suitable for multiple species and under different experimental conditions. Nitraria tangutorum Bobr. is a typical plant native to desert areas. It is drought-resistant, saline-alkali resistant, extreme temperatures-resistant, and has strong adaptability. To date, the importance of this germplasm has not been sufficiently understood; therefore, it is still unclear which genes can be used as reference genes to calibrate qPCR data of N. tangutorum. In this study we analyzed the expression levels of 10 candidate reference genes (ACT, GAPDH, TUA, TUB, CYP, UBC, His, PP2A, HSP, and EF1-α) in N. tangutorum seedlings under a series of experimental conditions, including in different organs (root, stem, and leaf) and under abiotic stresses (salt, drought, heat, and cold) and hormone stimuli (abscisic acid) by qPCR. Three software programs (geNorm, NormFinder, and BestKeeper) were used to evaluate the expression stability of the ten genes. Comprehensive analysis showed that EF1-α and His had the best expression stability, whereas HSP was the least suitable as a reference gene. The expression profile of NtCER7, a gene related to the regulation of cuticular wax biosynthesis in N. tangutorum, verified the accuracy of the experimental results. Based on this study, we recommend EF1-α and His as suitable reference genes for N. tangutorum. This paper provides the first data on stable reference genes in N. tangutorum, which will be beneficial to studying the gene expression of N. tangutorum and other Nitraria species in the future.

Evaluation of duplicated reference genes for quantitative real-time PCR analysis in genome unknown hexaploid oat (Avena sativa L.)

BACKGROUND

Oat (Avena sativa L.), a hexaploid crop with unknown genome, has valuable nutritional, medicinal and pharmaceutical uses. However, no suitable RGs (reference genes) for qPCR (quantitative real-time PCR) has been documented for oat yet. Single-copy gene is often selected as RG, which is challengeable or impactable in unexplored polyploids.

RESULTS

In this study, eleven candidate RGs, including four duplicated genes, were selected from oat transcriptome. The stability and the optimal combination of these candidate RGs were assessed in 18 oat samples by using four statistical algorithms including the ΔCt method, geNorm, NormFinder and BestKeeper. The most stable RGs for "all samples", "shoots and roots of seedlings", "developing seeds" and "developing endosperms" were EIF4A (Eukaryotic initiation factor 4A-3), UBC21 (Ubiquitin-Conjugating Enzyme 21), EP (Expressed protein) and EIF4A respectively. Among these RGs, UBC21 was a four-copy duplicated gene. The reliability was validated by the expression patterns of four various genes normalized to the most and the least stable RGs in different sample sets.

CONCLUSIONS

Results provide a proof of concept that the duplicated RG is feasible for qPCR in polyploids. To our knowledge, this study is the first systematic research on the optimal RGs for accurate qPCR normalization of gene expression in different organs and tissues of oat.

Bioaccumulation, antioxidative response, and metallothionein expression in Lupinus luteus L. exposed to heavy metals and silver nanoparticles

Yellow-lupin (Lupinus luteus L.) was grown on soils contaminated with heavy metals during two parallel studies. In the first one, the soil was contaminated by industrial activities whereas, in the second one, the soil was artificially contaminated with a single metal including Cd, Pb, Zn, Ni (in nitrate form), and Ag (in nitrate and nanoparticles form). The study was performed to assess a plant's response to contamination including its antioxidative response and molecular mechanisms involved in metal detoxification through the expression level of metallothioneins (MTs). Overall, the study provided insights into identification and validation of housekeeping genes (HKG) in L. luteus under exposure to metal stress and showed the effects of selected heavy metals and silver nanoparticles on the expression of metallothioneins, the activity of guaiacol peroxidase (GPX) and bioaccumulation of metals in leaves of L. luteus. As such, HKG validation using BestKeeper, NormFinder, and geNorm software allowed for the selection of four most stable reference genes in a context metal contamination for the selected plant. Moreover, a significant increase in the expression levels of MT was observed in plants grown under heavy metal stress and none on plants grown on 25 mg kg^-1 of silver nanoparticles. Also, the GPX activity and MT expression showed statistically significant changes between different conditions and doses which means that they can be used as highly sensitive stress markers for planning the phytoremediation process on a large scale.

Reference gene identification for reliable normalisation of quantitative RT-PCR data in Setaria viridis

BACKGROUND

Quantitative real-time polymerase chain reaction (RT-qPCR) is the key platform for the quantitative analysis of gene expression in a wide range of experimental systems and conditions. However, the accuracy and reproducibility of gene expression quantification via RT-qPCR is entirely dependent on the identification of reliable reference genes for data normalisation. Green foxtail (Setaria viridis) has recently been proposed as a potential experimental model for the study of C₄ photosynthesis and is closely related to many economically important crop species of the Panicoideae subfamily of grasses, including Zea mays (maize), Sorghum bicolor (sorghum) and Sacchurum officinarum (sugarcane). Setaria viridis (Accession 10) possesses a number of key traits as an experimental model, namely; (i) a small sized, sequenced and well annotated genome; (ii) short stature and generation time; (iii) prolific seed production, and; (iv) is amendable to Agrobacterium tumefaciens-mediated transformation. There is currently however, a lack of reference gene expression information for Setaria viridis (S. viridis). We therefore aimed to identify a cohort of suitable S. viridis reference genes for accurate and reliable normalisation of S. viridis RT-qPCR expression data.

RESULTS

Eleven putative candidate reference genes were identified and examined across thirteen different S. viridis tissues. Of these, the geNorm and NormFinder analysis software identified SERINE/THERONINE-PROTEIN PHOSPHATASE 2A (PP2A), 5'-ADENYLYLSULFATE REDUCTASE 6 (ASPR6) and DUAL SPECIFICITY PHOSPHATASE (DUSP) as the most suitable combination of reference genes for the accurate and reliable normalisation of S. viridis RT-qPCR expression data. To demonstrate the suitability of the three selected reference genes, PP2A, ASPR6 and DUSP, were used to normalise the expression of CINNAMYL ALCOHOL DEHYDROGENASE (CAD) genes across the same tissues.

CONCLUSIONS

This approach readily demonstrated the suitably of the three selected reference genes for the accurate and reliable normalisation of S. viridis RT-qPCR expression data. Further, the work reported here forms a highly useful platform for future gene expression quantification in S. viridis and can also be potentially directly translatable to other closely related and agronomically important C₄ crop species.

Sparse panicle1 is required for inflorescence development in Setaria viridis and maize

Setaria viridis is a rapid-life-cycle model panicoid grass. To identify genes that may contribute to inflorescence architecture and thus have the potential to influence grain yield in related crops such as maize, we conducted an N-nitroso-N-methylurea (NMU) mutagenesis of S. viridis and screened for visible inflorescence mutant phenotypes. Of the approximately 2,700 M₂ families screened, we identified four recessive sparse panicle mutants (spp1-spp4) characterized by reduced and uneven branching of the inflorescence. To identify the gene underlying the sparse panicle1 (spp1) phenotype, we performed bulked segregant analysis and deep sequencing to fine map it to an approximately 1 Mb interval. Within this interval, we identified disruptive mutations in two genes. Complementation tests between spp1 and spp3 revealed they were allelic, and deep sequencing of spp3 identified an independent disruptive mutation in SvAUX1 (AUXIN1), one of the two genes in the ∼1 Mb interval and the only gene disruption shared between spp1 and spp3. SvAUX1 was found to affect both inflorescence development and root gravitropism in S. viridis. A search for orthologous mutant alleles in maize confirmed a very similar role of ZmAUX1 in maize, which highlights the utility of S. viridis in accelerating functional genomic studies in maize.

Leaf transcriptome of two highly divergent genotypes of Urochloa humidicola (Poaceae), a tropical polyploid forage grass adapted to acidic soils and temporary flooding areas

BACKGROUND

Urochloa humidicola (Koronivia grass) is a polyploid (6x to 9x) species that is used as forage in the tropics. Facultative apospory apomixis is present in most of the genotypes of this species, although one individual has been described as sexual. Molecular studies have been restricted to molecular marker approaches for genetic diversity estimations and linkage map construction. The objectives of the present study were to describe and compare the leaf transcriptome of two important genotypes that are highly divergent in terms of their phenotypes and reproduction modes: the sexual BH031 and the aposporous apomictic cultivar BRS Tupi.

RESULTS

We sequenced the leaf transcriptome of Koronivia grass using an Illumina GAIIx system, which produced 13.09 Gb of data that consisted of 163,575,526 paired-end reads between the two libraries. We de novo-assembled 76,196 transcripts with an average length of 1,152 bp and filtered 35,093 non-redundant unigenes. A similarity search against the non-redundant National Center of Biotechnology Information (NCBI) protein database returned 65 % hits. We annotated 24,133 unigenes in the Phytozome database and 14,082 unigenes in the UniProtKB/Swiss-Prot database, assigned 108,334 gene ontology terms to 17,255 unigenes and identified 5,324 unigenes in 327 known metabolic pathways. Comparisons with other grasses via a reciprocal BLAST search revealed a larger number of orthologous genes for the Panicum species. The unigenes were involved in C4 photosynthesis, lignocellulose biosynthesis and flooding stress responses. A search for functional molecular markers revealed 4,489 microsatellites and 560,298 single nucleotide polymorphisms (SNPs). A quantitative real-time PCR analysis validated the RNA-seq expression analysis and allowed for the identification of transcriptomic differences between the two evaluated genotypes. Moreover, 192 unannotated sequences were classified as containing complete open reading frames, suggesting that the new, potentially exclusive genes should be further investigated.

CONCLUSION

The present study represents the first whole-transcriptome sequencing of U. humidicola leaves, providing an important public information source of transcripts and functional molecular markers. The qPCR analysis indicated that the expression of certain transcripts confirmed the differential expression observed in silico, which demonstrated that RNA-seq is useful for identifying differentially expressed and unique genes. These results corroborate the findings from previous studies and suggest a hybrid origin for BH031.

Screening and Validation of Housekeeping Genes of the Root and Cotyledon of Cunninghamia lanceolata under Abiotic Stresses by Using Quantitative Real-Time PCR

Cunninghamia lanceolata (Chinese fir) is a fast-growing and commercially important conifer of the Cupressaceae family. Due to the unavailability of complete genome sequences and relatively poor genetic background information of the Chinese fir, it is necessary to identify and analyze the expression levels of suitable housekeeping genes (HKGs) as internal reference for precise analysis. Based on the results of database analysis and transcriptome sequencing, we have chosen five candidate HKGs (Actin, GAPDH, EF1a, 18S rRNA, and UBQ) with conservative sequences in the Chinese fir and related species for quantitative analysis. The expression levels of these HKGs in roots and cotyledons under five different abiotic stresses in different time intervals were measured by qRT-PCR. The data were statistically analyzed using the following algorithms: NormFinder, BestKeeper, and geNorm. Finally, RankAggreg was applied to merge the sequences generated from three programs and rank these according to consensus sequences. The expression levels of these HKGs showed variable stabilities under different abiotic stresses. Among these, Actin was the most stable internal control in root, and GAPDH was the most stable housekeeping gene in cotyledon. We have also described an experimental procedure for selecting HKGs based on the de novo sequencing database of other non-model plants.

Selection of reliable reference genes for RT-qPCR analysis during developmental stages and abiotic stress in Setaria viridis

Real-time PCR (RT-qPCR) expression analysis is a powerful analytical technique, but reliable results depend on the use of stable reference genes for proper normalization. This study proposed to test the expression stability of 13 candidate reference genes in Setaria viridis, a monocot species recently proposed as a new C4 model plant. Gene expression stability of these genes was assayed across different tissues and developmental stages of Setaria and under drought or aluminum stress. In general, our results showed Protein Kinase, RNA Binding Protein and SDH as the most stable genes. Moreover, pairwise analysis showed that two reference genes were sufficient to normalize the gene expression data under each condition. By contrast, GAPDH and ACT were the least stably expressed genes tested. Validation of suitable reference genes was carried out to profile the expression of P5CS and GolS during abiotic stress. In addition, normalization of gene expression of SuSy, involved in sugar metabolism, was assayed in the developmental dataset. This study provides a list of reliable reference genes for transcript normalization in S. viridis in different tissues and stages of development and under abiotic stresses, which will facilitate genetic studies in this monocot model plant.

Validation of Reference Genes for Relative Quantitative Gene Expression Studies in Cassava (Manihot esculenta Crantz) by Using Quantitative Real-Time PCR

Reverse transcription quantitative real-time polymerase chain reaction (real-time PCR, also referred to as quantitative RT-PCR or RT-qPCR) is a highly sensitive and high-throughput method used to study gene expression. Despite the numerous advantages of RT-qPCR, its accuracy is strongly influenced by the stability of internal reference genes used for normalizations. To date, few studies on the identification of reference genes have been performed on cassava (Manihot esculenta Crantz). Therefore, we selected 26 candidate reference genes mainly via the three following channels: reference genes used in previous studies on cassava, the orthologs of the most stable Arabidopsis genes, and the sequences obtained from 32 cassava transcriptome sequence data. Then, we employed ABI 7900 HT and SYBR Green PCR mix to assess the expression of these genes in 21 materials obtained from various cassava samples under different developmental and environmental conditions. The stability of gene expression was analyzed using two statistical algorithms, namely geNorm and NormFinder. geNorm software suggests the combination of cassava4.1_017977 and cassava4.1_006391 as sufficient reference genes for major cassava samples, the union of cassava4.1_014335 and cassava4.1_006884 as best choice for drought stressed samples, and the association of cassava4.1_012496 and cassava4.1_006391 as optimal choice for normally grown samples. NormFinder software recommends cassava4.1_006884 or cassava4.1_006776 as superior reference for qPCR analysis of different materials and organs of drought stressed or normally grown cassava, respectively. Results provide an important resource for cassava reference genes under specific conditions. The limitations of these findings were also discussed. Furthermore, we suggested some strategies that may be used to select candidate reference genes.

Validation of Reference Genes for Relative Quantitative Gene Expression Studies in Cassava (Manihot esculenta Crantz) by Using Quantitative Real-Time PCR

Reverse transcription quantitative real-time polymerase chain reaction (real-time PCR, also referred to as quantitative RT-PCR or RT-qPCR) is a highly sensitive and high-throughput method used to study gene expression. Despite the numerous advantages of RT-qPCR, its accuracy is strongly influenced by the stability of internal reference genes used for normalizations. To date, few studies on the identification of reference genes have been performed on cassava (Manihot esculenta Crantz). Therefore, we selected 26 candidate reference genes mainly via the three following channels: reference genes used in previous studies on cassava, the orthologs of the most stable Arabidopsis genes, and the sequences obtained from 32 cassava transcriptome sequence data. Then, we employed ABI 7900 HT and SYBR Green PCR mix to assess the expression of these genes in 21 materials obtained from various cassava samples under different developmental and environmental conditions. The stability of gene expression was analyzed using two statistical algorithms, namely geNorm and NormFinder. geNorm software suggests the combination of cassava4.1_017977 and cassava4.1_006391 as sufficient reference genes for major cassava samples, the union of cassava4.1_014335 and cassava4.1_006884 as best choice for drought stressed samples, and the association of cassava4.1_012496 and cassava4.1_006391 as optimal choice for normally grown samples. NormFinder software recommends cassava4.1_006884 or cassava4.1_006776 as superior reference for qPCR analysis of different materials and organs of drought stressed or normally grown cassava, respectively. Results provide an important resource for cassava reference genes under specific conditions. The limitations of these findings were also discussed. Furthermore, we suggested some strategies that may be used to select candidate reference genes.

Selection of Reference Genes for Quantitative Real-Time RT-PCR Studies in Tomato Fruit of the Genotype MT-Rg1

Quantitative real-time RT-PCR (qRT-PCR) has become one of the most widely used methods for accurate quantification of gene expression. Since there are no universal reference genes for normalization, the optimal strategy to normalize raw qRT-PCR data is to perform an initial comparison of a set of independent reference genes to assess the most stable ones in each biological model. Normalization of a qRT-PCR experiment helps to ensure that the results are both statistically significant and biologically meaningful. Tomato is the model of choice to study fleshy fruit development. The miniature tomato (Solanum lycopersicum L.) cultivar Micro-Tom (MT) is considered a model system for tomato genetics and functional genomics. A new genotype, containing the Rg1 allele, improves tomato in vitro regeneration. In this work, we evaluated the expression stability of four tomato reference genes, namely CAC, SAND, Expressed, and ACTIN2. We showed that the genes CAC and Exp are the best reference genes of the four we tested during fruit development in the MT-Rg1 genotype. Furthermore, we validated the reference genes by showing that the expression profiles of the transcription factors FRUITFULL1 and APETALA2c during fruit development are comparable to previous reports using other tomato cultivars.