Comparative analysis of the male inflorescence transcriptome profiles of an ms22 mutant of maize

In modern agricultural production, maize is the most successful crop utilizing heterosis. 712C-ms22 is an important male sterile material in maize. In this study, we performed transcriptome sequencing analysis of the V10 stage of male inflorescence. Through this analysis, 27.63 million raw reads were obtained, and trimming of the raw data revealed 26.63 million clean reads, with an average match rate of 94.64%. Using Tophat software, we matched these clean reads to the maize reference genome. The abundance of 39,622 genes was measured, and 35,399 genes remained after filtering out the non-expressed genes across all the samples. These genes were classified into 19 categories by clusters of orthologous groups of protein annotation. Transcriptome sequencing analysis of the male sterile and fertile 712C-ms22 maize revealed some key DEGs that may be related to metabolic pathways. qRT-PCR analysis validated the gene expression patterns identified by RNA-seq. This analysis revealed some of the essential genes responsible for pollen development and for pollen tube elongation. Our findings provide useful markers of male sterility and new insights into the global mechanisms mediating male sterility in maize.


1.Introduction
Maize (Zea mays L.)is widely grown across the globe and constitutes the world's largest food crop, exceeding rice, wheat and other cereals in terms of importance since 2012 [1]. Maize is consumed in large quantities by both humans and animals, and has a number of industrial applications such asitsuse as a biofuel [2]. As an exemplar plant species, important research is frequently performed on maize, including studies of domestication, molecular evolution and genome architecture [3].
Over 11 million years ago, the maize genome was duplicated through an ancient polyploidy event, and the modern maize genome is replete with chromosomal duplications and repetitive DNA. In modern agricultural production, maize is the most successful crop that utilizes heterosis [4][5].
With advances in technology and with the completion of the maize B73 self-bred line genome sequence, a number of genomic, transcriptomic, epigenomic and proteomic approaches have increasingly been employed to provide new insights into the molecular a1111111111 a1111111111 a1111111111 a1111111111 a1111111111 mechanisms of plants [6]. Among these techniques, transcriptome analysis is a holistic approach that strives to understand the role and the interaction of individual biological components in phenotypic characteristics [7].
In recent years, advances in transcriptome sequencing, or RNA-seq, have allowed for the high throughput analysis of the transcriptome profiles of model plants [8][9], such as Arabidopsis thaliana, rice and maize, as well as other important crops such as soybeans, wheat, barley, Medicago, Sorghum and Setariaitalica [10][11][12]. By using transcriptome analysis in plants, a complete set of transcripts in a cell at a specific developmental stage or under specific physiological conditions can be obtained [13]; this set can then provide information on differentially expressed gene (DEG) levels and gene regulation for a specific pathway of interest [14]. Transcript analysis technology also allows for the direct comparison and the analysis of thousands of genes within one experiment [15]. By comparing DEGs, global transcriptome profiling analysis provides a mechanism with which to identify gene net works for the discovery of functional connections [16].
Yan and colleagues (2012) conducted genome-wide transcriptome sequencing analysis of the young floral buds of the Brassicanapus cytoplasmic male sterility (CMS)line Nsa and its novel restorer line NR1 and found a group of candidate genes associated with male sterility [17]. Similarly, Liu and colleagues (2013) employed Solexa/Illumina techniques to analyze the DEGs between the male sterile line 121A and its near isogenic line 121C in peppers. The comparison and analysis of these plant materials found a group of key genes and found significant pathways associated with male sterility [18]. Hao(2016)employed transcriptome profile analysis for drought-tolerant RIL70 and drought-sensitive RIL93 maize recombination inbred lines and detected DEGs between different stages of drought stress in these two lines, providing insights into the dynamic mechanisms underlying drought tolerance in maize seedlings [19][20].
New technological advances are required to increase the yield of hybrid seed production in maize [21]. Therefore, a break through in maize heterosis is required, similar to that which occurred for rice breeding in China. CMS and GMS are common among most plants, and male sterility exists widely among crops [22]. Male sterile lines are useful in the study of plant reproductive development, and these materials can be used as effective tools for hybrid seed production through heterosis. Moreover, GMS plays an important role in the production of hybrid maize seeds [23].
Therefore, the study of male sterility is valuable not only in terms of our theoretical understanding but also in terms of practical applications. However, the molecular mechanism of GMS in maize remains unclear. Although RNA-seq has been employed to examine a variety of plant species, few reports exist on the application of this technology in the study of GMS in maize.
In this study, the 712C-ms22mutantline was selected as the research material to examine the molecular genetic mechanism of GMS in maize. We performed comparative transcriptome analysis of 712C-ms22 mutant male inflorescence(V10) using Illumina sequencing technology to identify the DEGs and metabolic pathways that play a role in maize GMS. Our findings advance our understanding of the molecular mechanisms of male sterility in maize and may contribute towards technological advances, for example, in hybrid seed production and in other male sterility lines of maize. and the presence of pollen grains. However, compared with the fertile plant, the staminate inflorescences of sterile 712C-ms22 plants were abnormal in appearance, with tightly folded inflorescences and small, thin florets (Fig 1).
The mutant sterile males were smaller and had fewer branches than the normal fertile males.Male sterile (Clever shell)plants possessed tightly closed staminate flowers, whereas fertile (Clever shell) plants possessed loose staminate flowers that are easily cracked.Normal fertile plants were light green to pale yellow in color, whereas sterile plants were light green.
We also found that the ear and filament of female plants grew faster than those of male sterile plants. The ear growth of the fertile plants was slow, and the filaments were shorter. For cultivable plants,the filaments were not drawn out and were observed to be completely in powder form. The ears of fertile plants were short, whereas the fruit ears of sterile plants were longer with denser filaments. This phenomenon may be related to the nutritional distribution of reproductive growth.
After dissecting the fertile and sterile 712C-ms22 plants, phenotypic differences became evident (Fig 2). The fertile plant has normal anthers. In fertile flowers, the anther was plump,uniform, smooth and yellow in color, whereas in the sterile plant, the anther was shriveled, and the flowers were abnormal in appearance.

2.2.Transcriptome profile of maize in 712C-ms22
The cDNA libraries of the staminate inflorescences of 712C-ms22 sterile and fertile plants at the V10stage of development were prepared and sequenced using the IlluminaHiSeq 2500 platform. Three biological replicates for each sample were sequenced to increase the statistical power of the analysis. Reads of low sequencing quality were filtered out, and 53.34 million and 42.63 million 100 bp paired-end reads were obtained for the fertile and sterile materials, respectively. Of these reads, 78.59% mapped to unique positions. In total, 79.33% (209 million) of the total paired-end reads were aligned to the B73 reference genome (ZmB73_RefGen_v3). After trimming the raw data, 27.63million raw reads and 26.42million clean reads were obtained for the male inflorescence sample (Table 1). Uniquely mapped reads were used to estimate the transcript levels. Expression values were expressed in units of reads per kilobase per million reads mapped (RPKM). As expected, in this experiment, 94.65% of reads mapped to protein-coding genes.

Identification of the DEGs in sterile and fertile samples
Genes with similar expression patterns may also share similar functions, such as being involved in the same metabolic processes or pathways in the developmental stage [24]. To identify the function of genes, researchers use similar expression patterns as an analytical strategy to cluster genes, thereby providing clues as to the function of unknown genes [25]. DEG analysis showed that most of the detected genes were exclusively expressed through down-regulation in the male sterile line. The detected DEGs were characterized using biased information, calculated as the differences between male sterile and male fertility in flowers by using the log expression ratios (trimmed mean of M values).
To gain insights into the differential gene expression patterns between the sterile and fertile 712C-ms22 mutants in maize, we performed hierarchical clustering based on gene expression patterns to identify clusters with functional enrichment (Fig 3). DEGs were characterized as those genes with an RPKM value greater than or equal to 1 in at least one of the genotypes and with a false discovery rate (FDR)-adjusted P-value of less than 0.05. Using this significance level,we identified 395 DEGs between the sterile and fertile plants (log2FC>5 and q-value< 0.05)(S1 Table). We found that gene expression levels in fertile plants showed only small differences, whereas large differences were observed in the gene expression levels in sterile plants. Gene expression levels were calculated using the RPKM values. The DEGs between sample groups were defined using the fold-change values of the normalized RPKM expression values.  Most of these genes were expressed as down-regulated in the male sterile line. The relationship of all the significant DEGs was visualized in a heat map as shown in Fig 3.When comparing the sterile and fertile mutant plants of 712C-ms22,384 DEGs (log2FC>5 and P-value< 0.05) were detected, including nine up-regulated genes and 375 down-regulated genes. The DEGs common to both genotypes were mainly involved in signal transduction, metabolism, transport and transcriptional regulation. Most of the DEGs were unique for each comparison group. Both gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were conducted for these DEGs. The functional classification of DEGs was also carried out(S1 Table, Fig 3).

Functional classification of genes using GO and KEGG analyses
GO is an internationally standardized gene function classification system used to describe the properties of genes and their products in any organism. This classification system contains three ontologies:biological process, cellular component and molecular function [26].The application of the GO assignment system to DEGs can assist in understanding the distribution of gene functions at a macro level [27][28]. In this study, plant GO slim annotation was conducted using Blast2GO software (http://www.blast2go.com/b2ghome). The 1206 DEGs (P-value< 0.05)were enriched in 37 functional terms.The five most highly enriched GO terms were found to be catalytic activity(GO:0043086), metabolic process(GO:0008152), cellular process (GO:0009987), binding(GO:0005488) and single-organism process(GO:0044699).All unigenes were then assigned to the three main GO categories (biological process, cellular component and molecular function)(S2 Table, Fig 4).
Based on sequence homology, the 1206 DEGs were annotated into 36 functional categories, including 17 biological processes, 11 cellular components and eight molecular functions( Fig  4). Among the cellular component categories, "cell part and membrane" was the main functional group. Among the molecular function categories, "catalytic activity" was the main functional group, followed by "molecular function regulator" and "carbohydrate binding".
A GO term was considered significantly enriched if the corrected P-value was below 0.05. In total, 242 DEGs(P-value<0.05) were annotated by gene ontology to belong to 19 functional categories. These genes were similarly classified into 19 categories by clusters of orthologous groups of proteins annotation. Fig 4 shows the data from the GO analysis of the regulated genes in the 712C-ms22plants(S1 and S2 Tables, Fig 4).
Other significantly enriched GO terms were identified when comparing the sterile and fertile plants. Some DEGs were involved in carbohydrate metabolism, and these may be transcription factors that are able to regulate processes such as pollen development, the elimination of reactive oxygen species, cellular signal transduction and programmed cell death in plants.These included "pollen tube growth" (GO: 0009860), "pollen tube development" (GO: 0048868),"sexual reproduction" (GO:0019953),"cell differentiation" (GO: 0030154), "cell tip growth"(GO: 0009932),"developmental growth" (GO: 0048589), "peroxidase activity" (GO: 0004601) and"antioxidant activity"(GO: 0016209) in the male inflorescence mutants of 712C-ms22(S4 Table).
The DEGs common to both genotypes were mainly involved in signal transduction, metabolism, transport and transcriptional regulation. Most of the DEGs were unique for each comparison group. Both GO and KEGG analyses were conducted for these DEGs. The functional classification of DEGs was also carried out. KEGG pathway analysis provides valuable classifications for studying the complex biological functions of genes.

2.5.Analysis of DEGs potentially related to male sterility
Carbohydrate and energy level metabolism is an important pathway during the process of plant growth, asit provides the required energy and carbon sources for the plant's main  Comparative analysis of the male inflorescence transcriptome profiles of an ms22 mutant of maize physiological functions [29]. Defects in metabolic pathways may be a cause of male sterility in plants [30][31][32][33]. In this study, many DEGs were found to be involved in carbohydrate and energy metabolism as well as hydrolase activity and enzyme regulator activity. For example, 28 DEGs participated in carbohydrate and energy metabolism,26 in enzyme regulator activity and 35 in hydrolase activityin 712C-ms22 when comparing sterile and fertile plants.
To further investigate these results, line plots of the top 30 genes were drawnand confirmed that thesegenes had similar patterns that were down-regulatedin expression in the male sterile line(S4 Table, Fig 5, Fig 6).  Table, S3 Table).
We also found other DEGs potentially related to male sterility in 712C-ms22. In total, 21 DEGs were related to the elimination of reactive oxygen species, among which 18 were downregulated and three were up-regulated in 712C-ms22 when comparing the sterile and fertile plants. In addition, 10 DEGs associated with actin binding were detected, and 12 DEGs associated with calcium ion binding were found to be down-regulated (Table 3).

2.6.Analysis of DEGs by Qrt-PCR
According to the results of our functional and metabolic pathway analysis combined with the previously reported literature, 12 DEGs were randomly selected for Qrt-PCR analysis using the same sample as that used for RNA-seq. The expression patterns of the 12 DEGs were found to be consistent between Qrt-PCR and RNA-seq (Table 3; Fig 7), with a coincidence rate of 87.51% between the datasets, confirming that the RNA-seq results were reliable.

3.Discussion
GMS is common among most plants, and GMS plant material plays an important role in the production of hybrid crop seeds. To increase the yield as well as to overcome hybridization barriers, studies of maize gamete development, the pollen tube journey and fertilization mechanisms were initiated more than a century ago[33-35].
Previous studies have used transcriptional profiling or proteomic analysis to detect DEGs from sterile and fertile male lines [36][37][38][39]. During the process of plant growth and development, aberrant gene expression in any of the developmental stages relating to fertility (i.e., stamen development, male gametophyte development, pollen formation, pollen nutrient support and pollen tube elongation) can induce male sterility [40][41]. In the present study, RNA-seq was performed, and the transcripts of the staminate inflorescences between sterile and fertile males in ms22 were compared using Illumina sequencing technology. We identified 395 DEGs between the sterile and fertile plants (log2FC>5 and P-value<0.05), among which 386 DEGs were down-regulated and 9 DEGs were up-regulated in 712C-ms22(S1). To validate the RNAseq data, qRT-PCR was performed on 12 of these genes and consistent results were obtained on both platforms.
The spatial and temporal regulation of the expression of genes relating to programmed cell death metabolic pathways, which mediate tapetum degeneration and anther dehiscence, has been proven in many studies [42][43][44]. In our analysis, we found that cell wall loosening was a prominent feature in ms22,involving pectatelyase, polygalacturonase, expansins, cellulose and glucan endo-1 3-beta-glucosidase5. Polygalacturonaseis known to be involved in anther dehiscence, pollen mother cell wall degradation and pollen tube growth [45]. Pectatelyase plays a role in pollen wall development and pollen tube growth [46].
In this study, arabinogalactan protein (AGP)-encoding genes were detected among the DEGs that play an important role in the development of pollen, such as GPI-anchored proteins, GPI transamidase, COBRA-like proteins encoding GPI-anchored proteins and fasciclinlike AGP6,AGP7,AGP16,AGP20 and AGP17.AGPsparticipate in cell expansion, division, seed germination and pollen tube growth and guidance [52][53][54][55]. Calcium gradient-mediated pollen tube growth is one of the best-characterized metabolic systems [56]. Several studies have demonstrated that Ca2+, a messenger in cellular signal transduction, functions asapivotal regulator of the cell life cycle,including cell division, differentiation and apoptosis [57]. The germination of the male gametophyte (the pollen grain) and proper elongation of the growing tube are essential processes in the sexual reproduction of plants, which involves many signaling pathways that guide male sperm cells to their target, the haploid egg cell (the female gametophyte) [58]. In this study,the expression of genes relating to calcium gradient-mediated pollen tube growth was detected in ms22, such as calcium-dependent protein kinase 2, calcium-binding protein 39, calcium-dependent protein kinase ZmCPK11, calcium-binding EF-hand protein, putative calcium-dependent lipid-binding (CaLB domain) family protein and calcium ion-binding protein. The calcium-binding proteins act as calcium sensors and relay calcium signals [59]. Major calcium sensors have been reported, such as calcium-dependent protein kinase, CaM, CBL proteins, CDPK protein and CML protein.Pollen development depended on multiple signaling pathways, in which calmodulin was a key element. In Arabidopsis, CIPK19 is expressed specifically in pollen grains and controls pollen polarity [60][61].
The MADS-box transcription factor controls the specification of stamen primordia.MADSbox transcription factor 1 was detected by our transcriptome analysis. Zinc-finger family proteins play an important role in flower development [62]. The C2H2 zinc-finger family protein and Cys2/His2 zinc-finger transcription factor were found to be significantly activated during anther development in plants [63].
The final developmental stage of pollen maturation includes a programmed desiccation process for enhancing pollen-geminating efficiency. In this study, genes related to this process were detected, including pollen-specific proteins (C13, SF3), late embryogenesis abundant proteins (EMB564, D-34, Lea14-A) and dehydrin proteins (DHN1, DHN13). Additionally, the NADP reductase, ribonuclease, flap endonuclease, exosome complex exonuclease and sucrose cleavage protein were all detected by our transcriptome analysis. The study of the molecular mechanisms mediating male sterility in maize by comprehensive and systematic transcriptome profiling will provide valuable insights into the biological and developmental processes that occur throughout the life cycle of a plant66. We report DEGs that might be associated with the male sterility of ms22 in maize, but their specific functions require further studies.
Male sterility has been widely studied in response to both its biological significance and commercial use in hybrid seed production. Genic male sterility occurs when lesions in nuclear-encoded genes disrupt normal male gametogenesis. Anthers develop through a series of highly coordinated steps that are well-defined at the morphological level [64].
Recent research in maizehas revealed that some mutants with anther identity defects,such as ms-si-355 and ems71990, lack anthers within spikelets at the time when immature anthers normally exist. Furthermore, the mutants ms25, ms26 and ems-63131 exhibit highly irregular anther differentiation.It was found that microsporangia and all cell layers typical of the wildtype anther wall fail to differentiate.The ms8, ms23, ms25, ms26, ms32 and msca1 mutants each showed abnormalities in anther wall development. Mutant ms8 showed meiocyte and tapetum degradation, ms32showed anther wall layer defects and undifferentiated cell layers, and msca1-ms6064showed anther identity defects [65] MS23 regulates the specification and development of the tapetum in maize. Encoding an anther-specific predicted basic helix-loop-helix (bHLH) transcription factor required for tapetal differentiation.Microarray transcript profiling demonstrates a more severe developmental disruption in ms23-ref than in ms23 anthers, which possess a different bHLH defect [66].
In the ms8 mutant, tapetal cells become vacuolated and subsequently degrade. Interestingly, the ms8 mutant also exhibits several mild defects, including an excess number of epidermal cells that are shorter than normal, fewer tapetal cells that are larger than normal, and excess callose accumulation during meiosis [67].
In this study, 39 DEGs relating to actin cytoskeleton organization were expressed in ms22 including actin-depolymerizing factor 1 (ADF1), actin-depolymerizing factor 2 (ADF2) and actin-depolymerizing factor 3 (ADF3). ADF regulates the actin cytoskeleton in Arabidopsis and is involved in polarized tip growth [68]. As an actin-binding protein, it affects the cytoskeleton architecture dynamics and plays a role in the regulation of F-actin filament assembly [69].
The development of the male reproductive organ includes sporophytic cell division through mitosis and male gametogenesis through meiosis. Dehydrins or late embryogenesis abundant proteins protect the pollen from desiccation. Mitotic and meiotic-specific gene destruction during the cell division process can lead to male sterility [70][71]. In this study, we identified 17 DEGs relating to reproduction in the ms22mutant and 25 DEGs relating to sexual reproduction in this mutant,including cell division cycle proteins(cdt2, cdt48, cdt50),elongation factor Ts and cell division protein AAA of the ATPase family. Interestingly, the expression of this DEG was down-regulated in the male sterile maize in our study.

4.Conclusions
Using RNA-seq technology, we performed global transcriptome sequencing analysis of the ms22 mutant of male inflorescenceto compare sterility and fertility at the V10 developmental stage in maize. Our results demonstrated that additive gene expression patterns were prevalent in the staminate inflorescence development stages. This complementary effect may optimize gene expression levels and may play a fundamental role in the male inflorescence of maize.
Our analysis revealed that the male sterility of the ms22mutant in maize may be related to the regulated expression of some key DEGs, such as those genes involved in carbohydrate and energy metabolism, transcription factors, ribonuclease, pollen-specific kinase, reactive oxygen species elimination, cellular signal transduction, programmed cell death and reproductive development regulation. Our findings offer new insights into the molecular mechanisms of male sterility in maize. Future studies will focus on cloning and on determining the transgenic function of candidate genes associated with maize GMS. This analysis revealed some essential genes responsible for pollen development and pollen tube elongation. Our findings provide useful reference data for male sterility in maize and new insights into the global mechanisms responsible for this phenotype.

Preparation of plant materials
The 712C-ms22mutant strain (also known as 7712C-ms22 and U140G) was provided by the Maize Genetics Cooperation Stock Center. Three biological replicates of each of the ms22 mutant materials, 712C-ms22andB73, were planted at the Jiangpu Experimental Station, Nanjing Agricultural University, Nanjing, Jiangsu, China. At developmental stage V10, the male ear was relatively tender, and the anther in the glume could be detected by microscopic examination of the pollen grains, which could be distinguished between the sterile and fertile male ears. Young sterile and fertile floral material was clipped from 712C-ms22plants and then immediately frozen in liquid nitrogen and stored at −80˚C prior to use in RNA-seq and subsequent experiments.

Total RNA extraction and Illumina deep sequencing
Total RNA was extracted using RNAiso Plus Total RNA extraction reagent (TaKaRa, Dalian, China) following the manufacturer's instructions, and the RIN number was measured as an indicator of RNA integrity using an Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA). Qualified total RNA was further purified by anRNAClean XP Kit (Cat A63987, Beckman Coulter,Inc.Kraemer Boulevard, Brea, CA,USA)and an RNase-Free DNase Set (Qiagen, GmBH, Hilden, Germany). RNA samples with an RNA integrity number!8, as assigned by the BioAnalyzer, were used for microarray analysis [27]. After purification, total RNA was sequenced and then subjected to mRNA isolation, frag mentation, first strand and second strand cDNA synthesis, terminal repair, 3'terminal plus A, connection joint and enrichment.Finally, a complete sequence sample library was constructed.

Data analysis of RNA-seq
In this study,TopHat (version2.0.9) was applied for splice mapping of the pre-processed reads to the genome.Two mismatches were allowed, and each read was allowed multiple hits 2, based on the length of the gene and the number of RPKM calculated for each gene [15]. For the three biological replicates, the uniquely mapped reads in each gene were counted using HTSeqversion5.4 software and the genes with RPKM values consistently >0 or 0.25 were used for further analysis. The mapped reads were then assembled with Cufflinks (http://cufflinks. cbcb.umd.edu/) [8][9].

Identification of DEGs in ms22
The number of fragments for each gene was counted by Stringtie (version1.3.0),after comparison with TopHat and were then normalized by the TMM (trimmed mean of the M values) method. Edge R was applied to analyze gene differences between samples. P-value multiple hypothesis tests were corrected to control for the p-value threshold by controlling the false discovery rate (FDR). To identify the DEGs, the DESeq R package (1.10.1) was employed using a model based on a negative binomial distribution with default settings [27]. The multiplier of differential expression was calculated according to the FPKM value.The genes with RPKM values consistently >0 or 0.25 were included for further analysis. An FDR value 0.05 and a log2-fold-change!1 were used as the thresholds to judge the significance of gene expression differences [12].

GO annotation, COG annotation and KEGG enrichment pathway analysis
A hypergeometric test was applied to select GO entries that were significantly enriched in DEGs compared with the whole genome background. Using the GO database, genes can be classified according to the biological processes they participate in, the components of the cells, and the function of the molecules. GO and functional enrichment(http://www.geneontology. org/) analyses were conducted on all of the DEGs using the Goatools software(https://github. com/tanghaibao/goatools), with a q-value 0.05 defining significant enrichment of a GO term in DEGs [26].

qRT-PCR analysis
qRT-PCR analysis was used to verify the RNA-seq gene expression pattern. Reactions comprised 500 ng of total RNA, 1 μl of 10 mMdNTPs (Bioline, London, UK) and 1 μl of oligo-dT (20) primers (Invitrogen, Carlsbad, CA, USA), up to a final volume of 11 μl by the addition of RNase-free H 2 O. Reactions were performed at 65˚C for 5 min, followed by 50˚C for 60 min and were then inactivated at 70˚C for 15 min.
The qRT-PCR primers were designed based on reference unigene sequences with the Primer Premier5.0 software. The Fast Sybr R Green Master Mix (Applied Biosystems,Foster City, CA, USA) was employed, according to the manufacturer's instructions, in a reaction volume of 10 μl. Real-time PCR was conducted on an ABI 7500/7500 Fast Real-Time PCR system(Applied Biosystems).PCR conditions included initial denaturation for 2 min at 95˚C, followed by 40 cycles of denaturation at 95˚C for 30 s,hybridizationat 60˚Cfor 40 s, and elongation at 68˚C for 10 s. The actin2 gene was used as an internal control. The 2 -ΔΔct method was used to calculate the relative level of gene expression, and the B73 sample served as a control. A relative level of gene expression greater than 1 was considered to indicate up-regulation, and less than 1 indicateddown-regulation. All qRT-PCR reactions were performed with the three biological replicates.
Supporting information S1