The Phosphatome of Medicinal and Edible Fungus Wolfiporia cocos.

Wolfiporia cocos is an important medicinal and edible fungus that grows in association with pine trees, and its dried sclerotium has been used as a traditional medicine in China for centuries. However, the commercial production of W. cocos sclerotia is currently limited by shortages in pine wood resources. Since protein phosphatases (PPs) play significant roles in growth, signal transduction, development, metabolism, sexual reproduction, cell cycle, and environmental stress responses in fungi, the phosphatome of W. cocos was analyzed in this study by identifying PP genes, studying transcript profiles and assigning PPs to orthologous groups. Fifty-four putative PP genes were putatively identified in W. cocos genome based on homologous sequences searching using BLASTx program against the Saccharomyces cerevisiae, Fusarium graminearum, and Sclerotinia sclerotiorum databases. Based on known and presumed functions of orthologues of these PP genes found in other fungi, the putative roles of these W. cocos PPs in colonization, hyphal growth, sclerotial formation, secondary metabolism, and stress tolerance to environment were discussed in this study. And the level of transcripts from PP genes in the mycelium and sclerotium stages was also analyzed by qRT-PCR. Our study firstly identified and functional discussed the phosphatome in the medicinal and edible fungus W. cocos. The data from our study contribute to a better understanding of PPs potential roles in various cellar processes of W. cocos, and systematically provide comprehensive and novel insights into W. cocos economically important traits that could be extended to other fungi.

KRN4 Controls Quantitative Variation in Maize Kernel Row Number.

Kernel row number (KRN) is an important component of yield during the domestication and improvement of maize and controlled by quantitative trait loci (QTL). Here, we fine-mapped a major KRN QTL, KRN4, which can enhance grain productivity by increasing KRN per ear. We found that a ~3-Kb intergenic region about 60 Kb downstream from the SBP-box gene Unbranched3 (UB3) was responsible for quantitative variation in KRN by regulating the level of UB3 expression. Within the 3-Kb region, the 1.2-Kb Presence-Absence variant was found to be strongly associated with quantitative variation in KRN in diverse maize inbred lines, and our results suggest that this 1.2-Kb transposon-containing insertion is likely responsible for increased KRN. A previously identified A/G SNP (S35, also known as Ser220Asn) in UB3 was also found to be significantly associated with KRN in our association-mapping panel. Although no visible genetic effect of S35 alone could be detected in our linkage mapping population, it was found to genetically interact with the 1.2-Kb PAV to modulate KRN. The KRN4 was under strong selection during maize domestication and the favorable allele for the 1.2-Kb PAV and S35 has been significantly enriched in modern maize improvement process. The favorable haplotype (Hap1) of 1.2-Kb-PAV-S35 was selected during temperate maize improvement, but is still rare in tropical and subtropical maize germplasm. The dissection of the KRN4 locus improves our understanding of the genetic basis of quantitative variation in complex traits in maize.

De novo characterization of Panax japonicus C. A. Mey transcriptome and genes related to triterpenoid saponin biosynthesis.

Panax japonicus C.A.Mey, the traditional medicinal herb in the Araliaceae family, has been used as a tonic, anti-inflammatory and haemostatic agent in China for more than thousand years. Its clinical effects are mainly due to the presence of triterpenoid saponins. However, little is known at the genetic level about how saponins are biosynthesized in this plant. We have therefore performed the de novo transcriptome assembly and high throughput RNA-seq analysis for P. japonicus. 66,403 unigenes were assembled from 19.6 Gbp raw data, and 34,639 unigenes were annotated. After annotation, 29 unigenes involved in putative backbone biosynthesis of triterpenoid saponin were selected. Additionally, 34 Cytochrome P450 and 18 UDP-glycosyltransferase unigenes were predicted based on the annotation, which were related to the saponin backbone modification. The expression level of related key genes were further verified by qPCR analysis. The results of this study provide the most comprehensive expressed sequence resources for P. japonicus, which will enlarge the available P. japonicus gene pool and facilitate further genome-wide research and analyses in this species.

Multiple loci not only Rf3 involved in the restoration ability of pollen fertility, anther exsertion and pollen shedding to S type cytoplasmic male sterile in maize.

KEY MESSAGE: Thirty loci for fertility restoration of pollen fertility, anther exsertion and pollen shedding to maize CMS-S were identified by GWAS. S type cytoplasmic male sterile (CMS-S) is the main type of CMS in maize; poor understanding of the genetic architecture of fertility restoration to CMS-S is one of the reasons to impede its utility in hybrid breeding. In this study, genome-wide identification of genetic loci for fertility restoration ability to CMS-S was firstly conducted with a set of testcrossing association mapping panel in three environments. A total of 19, 3 and 8 significant loci (P < 1.8 × 10(-6), α = 1) for pollen fertility, anther exsertion and pollen shedding were identified, respectively, and individual locus explained up to 28.26% of phenotypic variation. Of them, only Rf3, the main restorer-fertility gene of CMS-S, was identified for the three traits simultaneously. In addition, 83 candidate genes within the 100 kb extension regions of these loci were predicted. These results revealed that besides Rf3 multiple genetic loci and mechanisms are involved in the fertility restoration ability to CMS-S. Results in this study would provide important information for understanding the genetic architecture of fertility restoration to CMS-S in maize.

Genetic architecture of maize kernel row number and whole genome prediction.

KEY MESSAGE: Maize kernel row number might be dominated by a set of large additive or partially dominant loci and several small dominant loci and can be accurately predicted by fewer than 300 top KRN-associated SNPs. Kernel row number (KRN) is an important yield component in maize and directly affects grain yield. In this study, we combined linkage and association mapping to uncover the genetic architecture of maize KRN and to evaluate the phenotypic predictability using these detected loci. A genome-wide association study revealed 31 associated single nucleotide polymorphisms (SNPs) representing 17 genomic loci with an effect in at least one of five individual environments and the best linear unbiased prediction (BLUP) over all environments. Linkage mapping in three F2:3 populations identified 33 KRN quantitative trait loci (QTLs) representing 21 QTLs common to several population/environments. The majority of these common QTLs that displayed a large effect were additive or partially dominant. We found 70% KRN-associated genomic loci were mapped in KRN QTLs identified in this study, KRN-associated SNP hotspots detected in NAM population and/or previous identified KRN QTL hotspots. Furthermore, the KRN of inbred lines and hybrids could be predicted by the additive effect of the SNPs, which was estimated using inbred lines as a training set. The prediction accuracy using the top KRN-associated tag SNPs was obviously higher than that of the randomly selected SNPs, and approximately 300 top KRN-associated tag SNPs were sufficient for predicting the KRN of the inbred lines and hybrids. The results suggest that the KRN-associated loci and QTLs that were detected in this study show great potential for improving the KRN with genomic selection in maize breeding.

[Transcriptome profiling and analysis of Panax japonicus var. major].

The rhizome of Panax japonicus var. major have been used as the natural medicinal agent by Chinese traditional doctors for more than thousand years. Most of the therapeutic effects of P. japonicus var. major had been reported due to the presence of tetracyclic or pentacyclic triterpene saponins. In this study, Illumina pair-end RNA-sequencing and de novo splicing were done in order to understand the pathway of triterpenoid saponins in this species. The valid reads data of 15. 6 Gb were obtained. The 62 240 unigenes were finally obtained by de novo splicing. After annotation, we discovered 19 unigenes involved in ginsenoside backbone biosynthesis. Additionally, 69 unigenes and 18 unigenes were predicted to have potential function of cytochrome P450 and UDP-glycosyltransferase based on the annotation results, which may encode enzymes responsible for ginsenoside backbone modification. This study provides global expressed datas for P. japonicus var. major, which will contribute significantly to further genome-wide research and analysis for this species.

ZmGA3ox2, a candidate gene for a major QTL, qPH3.1, for plant height in maize.

Maize plant height is closely associated with biomass, lodging resistance and grain yield. Determining the genetic basis of plant height by characterizing and cloning plant height genes will guide the genetic improvement of crops. In this study, a quantitative trait locus (QTL) for plant height, qPH3.1, was identified on chromosome 3 using populations derived from a cross between Zong3 and its chromosome segment substitution line, SL15. The plant height of the two lines was obviously different, and application of exogenous gibberellin A(3) removed this difference. QTL mapping placed qPH3.1 within a 4.0 cM interval, explaining 32.3% of the phenotypic variance. Furthermore, eight homozygous segmental isolines (SILs) developed from two larger F(2) populations further narrowed down qPH3.1 to within a 12.6 kb interval. ZmGA3ox2, an ortholog of OsGA3ox2, which encodes a GA3 ß-hydroxylase, was positionally cloned. Association mapping identified two polymorphisms in ZmGA3ox2 that were significantly associated with plant height across two experiments. Quantitative RT-PCR showed that SL15 had higher ZmGA3ox2 expression relative to Zong3. The resultant higher GA(1) accumulation led to longer internodes in SL15 because of increased cell lengths. Moreover, a large deletion in the coding region of ZmGA3ox2 is responsible for the dwarf mutant d1-6016. The successfully isolated qPH3.1 enriches our knowledge on the genetic basis of plant height in maize, and provides an opportunity for improvement of plant architecture in maize breeding.

Genetic analysis and major QTL detection for maize kernel size and weight in multi-environments.

KEY MESSAGE: Twelve major QTL in five optimal clusters and several epistatic QTL are identified for maize kernel size and weight, some with pleiotropic will be promising for fine-mapping and yield improvement. Kernel size and weight are important target traits in maize (Zea mays L.) breeding programs. Here, we report a set of quantitative trait loci (QTL) scattered through the genome and significantly controlled the performance of four kernel traits including length, width, thickness and weight. From the cross V671 (large kernel) × Mc (small kernel), 270 derived F2:3 families were used to identify QTL of maize kernel-size traits and kernel weight in five environments, using composite interval mapping (CIM) for single-environment analysis along with mixed linear model-based CIM for joint analysis. These two mapping strategies identified 55 and 28 QTL, respectively. Among them, 6 of 23 coincident were detected as interacting with environment. Single-environment analysis showed that 8 genetic regions on chromosomes 1, 2, 4, 5 and 9 clustered more than 60 % of the identified QTL. Twelve stable major QTLs accounting for over 10 % of phenotypic variation were included in five optimal clusters on the genetic region of bins 1.02-1.03, 1.04-1.06, 2.05-2.07, 4.07-4.08 and 9.03-9.04; the addition and partial dominance effects of significant QTL play an important role in controlling the development of maize kernel. These putative QTL may have great promising for further fine-mapping with more markers, and genetic improvement of maize kernel size and weight through marker-assisted breeding.

Genome-wide identification, evolution, and expression analysis of RNA-binding glycine-rich protein family in maize.

The RNA-binding glycine-rich protein (RB-GRP) family is characterized by the presence of a glycine-rich domain arranged in (Gly)n-X repeats and an RNA-recognition motif (RRM). RB-GRPs participate in varied physiological and biochemical processes especially in the stress response of plants. In this study, a total of 23 RB-GRPs distributed on 10 chromosomes were identified in maize (Zea mays L.), and they were divided into four subgroups according to their conserved domain architecture. Five pairs of paralogs were identified, while none of them was located on the same chromosomal region, suggesting that segmental duplication is predominant in the duplication events of the RB-GRPs in maize. Comparative analysis of RB-GRPs in maize, Arabidopsis (Arabidopsis thaliana L.), rice (Oryza sativa L.), and wheat (Triticum aestivum) revealed that two exclusive subgroups were only identified in maize. Expression of eight ZmRB-GRPs was significantly regulated by at least two kinds of stresses. In addition, cis-elements predicted in the promoter regions of the ZmRB-GRPs also indicated that these ZmRB-GRPs would be involved in stress response of maize. The preliminary genome-wide analysis of the RB-GRPs in maize would provide useful information for further study on the function of the ZmRB-GRPs.

Digital gene expression analysis of early root infection resistance to Sporisorium reilianum f. sp. zeae in maize.

The maize smut fungus, Sporisorium reilianum f. sp. zeae, which is an important biotrophic pathogen responsible for extensive crop losses, can infect maize by invading the root during the early seedling stage. In order to investigate disease-resistance mechanisms at this early seedling stage, digital gene expression analysis, which applies a dual-enzyme approach, was used to identify the transcriptional changes in the roots of Huangzao4 (susceptible) and Mo17 (resistant) after root inoculation with S. reilianum. During the infection in the roots, the expression pattern of pathogenesis-related genes in Huangzao4 and Mo17 were significantly differentially regulated at different infection stages. The glutathione S-transferase enzyme activity and reactive oxygen species levels also showed changes before and after inoculation. The total lignin contents and the pattern of lignin depositions in the roots differed during root colonization of Huangzao4 and Mo17. These results suggest that the interplay between S. reilianum and maize during the early infection stage involves many important transcriptional and physiological changes, which offer several novel insights to understanding the mechanisms of resistance to the infection of biotrophic fungal pathogens.

Floral transition in maize infected with Sporisorium reilianum disrupts compatibility with this biotrophic fungal pathogen.

Sporisorium reilianum f. sp. zeae is an important biotrophic pathogen that causes head smut disease in maize. Head smut is not obvious until the tassels and ears emerge. S. reilianum has a very long life cycle that spans almost the entire developmental program of maize after the pathogen successfully invades the root. The aim of this study was to understand at a molecular level how this pathogen interacts with the host during its long life cycle, and how this interaction differs between susceptible and resistant varieties of maize after hyphal invasion. We investigated transcriptional changes in the resistant maize line Mo17 at four developmental stages using a maize 70mer-oligonucleotide microarray. We found that there was a lengthy compatible relationship between the pathogen and host until the early eighth-leaf stage. The resistance in Mo17 relied on the assignment of auxin and regulation of flavonoids in the early floral primordium during the early floral transition stage. We propose a model describing the putative mechanism of head smut resistance in Mo17 during floral transition. In the model, the synergistic regulations among auxin, flavonoids, and hyphal growth play a key role in maintaining compatibility with S. reilianum in the resistant maize line.

Uncovering small RNA-mediated responses to cold stress in a wheat thermosensitive genic male-sterile line by deep sequencing.

The male sterility of thermosensitive genic male sterile (TGMS) lines of wheat (Triticum aestivum) is strictly controlled by temperature. The early phase of anther development is especially susceptible to cold stress. MicroRNAs (miRNAs) play an important role in plant development and in responses to environmental stress. In this study, deep sequencing of small RNA (smRNA) libraries obtained from spike tissues of the TGMS line under cold and control conditions identified a total of 78 unique miRNA sequences from 30 families and trans-acting small interfering RNAs (tasiRNAs) derived from two TAS3 genes. To identify smRNA targets in the wheat TGMS line, we applied the degradome sequencing method, which globally and directly identifies the remnants of smRNA-directed target cleavage. We identified 26 targets of 16 miRNA families and three targets of tasiRNAs. Comparing smRNA sequencing data sets and TaqMan quantitative polymerase chain reaction results, we identified six miRNAs and one tasiRNA (tasiRNA-ARF [for Auxin-Responsive Factor]) as cold stress-responsive smRNAs in spike tissues of the TGMS line. We also determined the expression profiles of target genes that encode transcription factors in response to cold stress. Interestingly, the expression of cold stress-responsive smRNAs integrated in the auxin-signaling pathway and their target genes was largely noncorrelated. We investigated the tissue-specific expression of smRNAs using a tissue microarray approach. Our data indicated that miR167 and tasiRNA-ARF play roles in regulating the auxin-signaling pathway and possibly in the developmental response to cold stress. These data provide evidence that smRNA regulatory pathways are linked with male sterility in the TGMS line during cold stress.

Identification of combining ability loci for five yield-related traits in maize using a set of testcrosses with introgression lines.

Combining ability is essential for hybrid breeding in crops. However, the genetic basis of combining ability remains unclear and has been seldom investigated. Identifying molecular markers associated with this complex trait would help to understand its genetic basis and provide useful information for hybrid breeding in maize. In this study, we identified genetic loci of general combining ability (GCA) and specific combining ability (SCA) for five yield-related traits under three environments using a set of testcrosses with introgression lines (ILs). GCA or SCA of the five yield-related traits of the ILs was estimated by the performance of testcrosses with four testers from different heterotic groups. Genetic correlations between GCA of the traits and the corresponding traits per se were not significant or not strong, suggesting that the genetic basis between them is different. A total of 56 significant loci for GCA and 21 loci for SCA were commonly identified in at least two environments, and only 5 loci were simultaneously controlling GCA and SCA, indicating that the genetic basis of GCA and SCA is different. For all of the traits investigated, positive and significant correlations between the number of GCA loci in the ILs and the performance of the corresponding GCA of the ILs were detected, implying that pyramiding GCA loci would have positive effect on the performance of GCA. Results in this study would be useful for maize hybrid breeding.

Genome-wide association analysis for nine agronomic traits in maize under well-watered and water-stressed conditions.

Drought can cause severe reduction in maize production, and strongly threatens crop yields. To dissect this complex trait and identify superior alleles, 350 tropical and subtropical maize inbred lines were genotyped using a 1536-SNP array developed from drought-related genes and an array of 56,110 random SNPs. The inbred lines were crossed with a common tester, CML312, and the testcrosses were phenotyped for nine traits under well-watered and water-stressed conditions in seven environments. Using genome-wide association mapping with correction for population structure, 42 associated SNPs (P ≤ 2.25 × 10(-6) 0.1/N) were identified, located in 33 genes for 126 trait × environment × treatment combinations. Of these genes, three were co-localized to drought-related QTL regions. Gene GRMZM2G125777 was strongly associated with ear relative position, hundred kernel weight and timing of male and female flowering, and encodes NAC domain-containing protein 2, a transcription factor expressed in different tissues. These results provide some good information for understanding the genetic basis for drought tolerance and further studies on identified candidate genes should illuminate mechanisms of drought tolerance and provide tools for designing drought-tolerant maize cultivars tailored to different environmental scenarios.

Genome-wide identification and analysis of microRNA responding to long-term waterlogging in crown roots of maize seedlings.

MicroRNAs (miRNAs) are critical post-transcriptional modulators of gene expression involving in plant responses to abiotic stress. However, the regulation of miRNA in the morphological response to waterlogging is poorly understood in maize. In this study, we detected miRNAs and their targets that expressed in waterlogged crown roots of maize seedlings in two inbred lines (Hz32 and Mo17) by RNA sequencing. A total of 61 mature miRNAs were found including 36 known maize (zma) miRNAs and 25 potential novel miRNA candidates. Comparison of miRNA expression in both waterlogged and control crown roots revealed 32 waterlogging-responsive miRNAs, most were consistently downregulated under waterlogging in the two inbred lines. We identified the miRNA targets through degradome sequencing. Many known miRNA targets involving in transcription regulation and reactive oxygen species elimination were found in the degradome libraries, and 17 targets of 10 newly detected miRNAs were identified as well. Moreover, the miRNA-mediated pathways that respond to waterlogging and regulate the induction of crown roots were discussed. This study is a comprehensive survey of responsive miRNAs in waterlogged maize crown roots. The results will help to understand the miRNA expression in response to waterlogging and miRNA-mediated regulation of morphological adaptation to waterlogging in maize.

Cold stress contributes to aberrant cytokinesis during male meiosis I in a wheat thermosensitive genic male sterile line.

The male sterility of a wheat thermosensitive genic male sterile (TGMS) line is strictly controlled by temperature. When the TGMS line BS366 was exposed to 10 °C from the pollen mother cell stage to the meiosis stage, a few pollen grains were formed and devoid of starch. We report here a large-scale transcriptomic study using the Affymetrix wheat GeneChip to follow gene expression in BS366 line anthers in response to cold stress. Notably, many cytoskeletal signaling components were gradually induced in response to cold stress in BS366 line anthers. However, the cytoskeleton-associated genes that play key roles in the dynamic organization of the cytoskeleton were dramatically repressed. Histological studies revealed that the separation of dyads occurred abnormally during male meiosis I, indicating defective male meiotic cytokinesis. Fluorescence labelling and subcellular histological observations revealed that the phragmoplast was defectively formed and the cell plate was abnormally assembled during meiosis I under cold stress. Based on the transcriptomic analysis and observations of characterized histological changes, our results suggest that cold stress repressed transcription of cytoskeleton dynamic factors and subsequently caused the defective cytokinesis during meiosis I. The results may explain the male sterility caused by low temperature in wheat TGMS lines.

Prolyl 4-hydroxylase genes are subjected to alternative splicing in roots of maize seedlings under waterlogging.

BACKGROUND: In animals, prolyl 4-hydroxylases (P4Hs) are regarded as oxygen sensors under hypoxia stress, but little is known about their role in the response to waterlogging in maize. METHODS: A comprehensive genome-wide analysis of P4H genes of maize (zmP4H genes) was carried out, including gene structures, phylogeny, protein motifs, chromosomal locations and expression patterns under waterlogging. KEY RESULTS: Nine zmP4H genes were identified in maize, of which five were alternatively spliced into at least 19 transcripts. Different alternative splicing (AS) events were revealed in different inbred lines, even for the same gene, possibly because of organ and developmental specificities or different stresses. The signal strength of splice sites was strongly correlated with selection of donor and receptor sites, and ambiguous junction sites due to small direct repeats at the exon/intron junction frequently resulted in the selection of unconventional splicing sites. Eleven out of 14 transcripts resulting from AS harboured a premature termination codon, rendering them potential candidates for nonsense-mediated RNA degradation. Reverse transcription-PCR (RT-PCR) indicated that zmP4H genes displayed different expression patterns under waterlogging. The diverse transcripts generated from AS were expressed at different levels, suggesting that zmP4H genes were under specific control by post-transcriptional regulation under waterlogging stress in the line HZ32. CONCLUSIONS: Our results provide a framework for future dissection of the function of the emerging zmP4H family and suggest that AS might have an important role in the regulation of the expression profile of this gene family under waterlogging stress.

Identification of transcriptome induced in roots of maize seedlings at the late stage of waterlogging.

BACKGROUND: Plants respond to low oxygen stress, particularly that caused by waterlogging, by altering transcription and translation. Previous studies have mostly focused on revealing the mechanism of the response at the early stage, and there is limited information about the transcriptional profile of genes in maize roots at the late stage of waterlogging. The genetic basis of waterlogging tolerance is largely unknown. In this study, the transcriptome at the late stage of waterlogging was assayed in root cells of the tolerant inbred line HZ32, using suppression subtractive hybridization (SSH). A forward SSH library using RNA populations from four time points (12 h, 16 h, 20 h and 24 h) after waterlogging treatment was constructed to reveal up-regulated genes, and transcriptional and linkage data was integrated to identify candidate genes for waterlogging tolerance. RESULTS: Reverse Northern analysis of a set of 768 cDNA clones from the SSH library revealed a large number of genes were up-regulated by waterlogging. A total of 465 ESTs were assembled into 296 unigenes. Bioinformatic analysis revealed that the genes were involved in complex pathways, such as signal transduction, protein degradation, ion transport, carbon and amino acid metabolism, and transcriptional and translational regulation, and might play important roles at the late stage of the response to waterlogging. A significant number of unigenes were of unknown function. Approximately 67% of the unigenes could be aligned on the maize genome and 63 of them were co-located within reported QTLs. CONCLUSION: The late response to waterlogging in maize roots involves a broad spectrum of genes, which are mainly associated with two response processes: defense at the early stage and adaption at the late stage. Signal transduction plays a key role in activating genes related to the tolerance mechanism for survival during prolonged waterlogging. The crosstalk between carbon and amino acid metabolism reveals that amino acid metabolism performs two main roles at the late stage: the regulation of cytoplasmic pH and energy supply through breakdown of the carbon skeleton.

Association mapping of agronomic traits on chromosome 2A of wheat.

Association mapping is a method to test the association between molecular markers and quantitative trait loci (QTL) based on linkage disequilibrium (LD). In this study, the collection of 108 wheat germplasm accessions form China were evaluated for their plant heights, spike length, spikelets per spike, grains per spike, thousand kernel weight and spikelets density in 3 years at three locations. And they were genotyped with 85 SSR markers and 40 EST-SSR markers. The population structure was inferred on the basis of unlinked 48 SSR markers and 40 EST-SSR markers. The extent of LD on chromosome 2A was 2.3 cM. Association of 37 SSR loci on chromosomes 2A with six agronomic traits was analysed with a mixed linear model. A total of 14 SSR loci were significantly associated with agronomic traits. Some of the associated markers were located in the QTL region detected in previous linkage mapping analysis. Our results demonstrated that association mapping can enhance QTL information and achieves higher resolution with short LD extent.

Differential expression of miRNAs in response to salt stress in maize roots.

BACKGROUND AND AIMS: Corn (Zea mays) responds to salt stress via changes in gene expression, metabolism and physiology. This adaptation is achieved through the regulation of gene expression at the transcriptional and post-transcriptional levels. MicroRNAs (miRNAs) have been found to act as key regulating factors of post-transcriptional gene expression. However, little is known about the role of miRNAs in plants' responses to abiotic stresses. METHODS: A custom microparaflo microfluidic array containing release version 10.1 plant miRNA probes (http://microrna.sanger.ac.uk/) was used to discover salt stress-responsive miRNAs using the differences in miRNA expression between the salt-tolerant maize inbred line 'NC286' and the salt-sensitive maize line 'Huangzao4'. Key Results miRNA microarray hybridization revealed that a total of 98 miRNAs, from 27 plant miRNA families, had significantly altered expression after salt treatment. These miRNAs displayed different activities in the salt response, and miRNAs belonging to the same miRNA family showed the same behaviour. Interestingly, 18 miRNAs were found which were only expressed in the salt-tolerant maize line, and 25 miRNAs that showed a delayed regulation pattern in the salt-sensitive line. A gene model was proposed that showed how miRNAs could regulate the abiotic stress-associated process and the gene networks coping with the stress. CONCLUSIONS: Salt-responsive miRNAs are involved in the regulation of metabolic, morphological and physiological adaptations of maize seedlings at the post-transcriptional level. The miRNA genotype-specific expression model might explain the distinct salt sensitivities between maize lines.