Identification of Accurate Reference Genes for qRT-PCR Analysis of Gene Expression in Eremochloa ophiuroides under Multiple Stresses of Phosphorus Deficiency and/or Aluminum Toxicity

Centipedegrass (Eremochloa ophiuroides (Munro.) Hack.) is a species originating in China and is an excellent warm-season turfgrass. As a native species in southern China, it is naturally distributed in the phosphorus-deficient and aluminum-toxic acid soil areas. It is important to research the molecular mechanism of centipedegrass responses to phosphorus-deficiency and/or aluminum-toxicity stress. Quantitative Real-Time PCR (qRT-PCR) is a common method for gene expression analysis, and the accuracy of qRT-PCR results depends heavily on the stability of internal reference genes. However, there are still no reported stable and effective reference genes for qRT-PCR analysis of target genes under the acid-soil-related stresses in different organs of centipedegrass. For scientific rigor, the gene used as a reference for any plant species and/or any stress conditions should be first systematically screened and evaluated. This study is the first to provide a group of reliable reference genes to quantify the expression levels of functional genes of Eremochloa ophiuroides under multiple stresses of P deficiency and/or aluminum toxicity. In this study, centipedegrass seedlings of the acid-soil-resistant strain ‘E041’ and acid-soil-sensitive strain ‘E089’ were used for qRT-PCR analysis. A total of 11 candidate reference genes (ACT, TUB, GAPDH, TIP41, CACS, HNR, EP, EF1α, EIF4α, PP2A and actin) were detected by qRT-PCR technology, and the stability of candidate genes was evaluated with the combination of four internal stability analysis software programs. The candidate reference genes exhibited differential stability of expression in roots, stems and leaves under phosphorus-deficiency and/or aluminum-toxicity stress. On the whole, the results showed that GAPDH, TIP41 and HNR were the most stable in the total of samples. In addition, for different tissues under various stresses, the selected reference genes were also different. CACS and PP2A were identified as two stable reference genes in roots through all three stress treatments (phosphate deficiency, aluminum toxicity, and the multiple stress treatment of aluminum toxicity and phosphate deficiency). Moreover, CACS was also stable as a reference gene in roots under each treatment (phosphate deficiency, aluminum toxicity, or multiple stresses of aluminum toxicity and phosphate deficiency). In stems under all three stress treatments, GAPDH and EIF4α were the most stable reference genes; for leaves, PP2A and TIP41 showed the two highest rankings in all three stress treatments. Finally, qRT-PCR analysis of the expression patterns of the target gene ALMT1 was performed to verify the selected reference genes. The application of the reference genes identified as internal controls for qRT-PCR analysis will enable accurate analysis of the target gene expression levels and expression patterns in centipedegrass under acid-soil-related stresses.