Haijian Lin

Identification of miRNAs and their target genes in developing maize ears by combined small RNA and degradome sequencing

BackgroundIn plants, microRNAs (miRNAs) are endogenous ~22 nt RNAs that play important regulatory roles in many aspects of plant biology, including metabolism, hormone response, epigenetic control of transposable elements, and stress response. Extensive studies of miRNAs have been performed in model plants such as rice and Arabidopsis thaliana. In maize, most miRNAs and their target genes were analyzed and identified by clearly different treatments, such as response to low nitrate, salt and drought stress. However, little is known about miRNAs involved in maize ear development. The objective of this study is to identify conserved and novel miRNAs and their target genes by combined small RNA and degradome sequencing at four inflorescence developmental stages.ResultsWe used deep-sequencing, miRNA microarray assays and computational methods to identify, profile, and describe conserved and non-conserved miRNAs at four ear developmental stages, which resulted in identification of 22 conserved and 21-maize-specific miRNA families together with their corresponding miRNA*. Comparison of miRNA expression in these developmental stages revealed 18 differentially expressed miRNA families. Finally, a total of 141 genes (251 transcripts) targeted by 102 small RNAs including 98 miRNAs and 4 ta-siRNAs were identified by genomic-scale high-throughput sequencing of miRNA cleaved mRNAs. Moreover, the differentially expressed miRNAs-mediated pathways that regulate the development of ears were discussed.ConclusionsThis study confirmed 22 conserved miRNA families and discovered 26 novel miRNAs in maize. Moreover, we identified 141 target genes of known and new miRNAs and ta-siRNAs. Of these, 72 genes (117 transcripts) targeted by 62 differentially expressed miRNAs may attribute to the development of maize ears. Identification and characterization of these important classes of regulatory genes in maize may improve our understanding of molecular mechanisms controlling ear development.

Genome expression profile analysis of the immature maize embryo during dedifferentiation.

Maize is one of the most important cereal crops worldwide and one of the primary targets of genetic manipulation, which provides an excellent way to promote its production. However, the obvious difference of the dedifferentiation frequency of immature maize embryo among various genotypes indicates that its genetic transformation is dependence on genotype and immature embryo-derived undifferentiated cells. To identify important genes and metabolic pathways involved in forming of embryo-derived embryonic calli, in this study, DGE (differential gene expression) analysis was performed on stages I, II, and III of maize inbred line 18-599R and corresponding control during the process of immature embryo dedifferentiation. A total of ∼21 million cDNA tags were sequenced, and 4,849,453, 5,076,030, 4,931,339, and 5,130,573 clean tags were obtained in the libraries of the samples and the control, respectively. In comparison with the control, 251, 324 and 313 differentially expressed genes (DEGs) were identified in the three stages with more than five folds, respectively. Interestingly, it is revealed that all the DEGs are related to metabolism, cellular process, and signaling and information storage and processing functions. Particularly, the genes involved in amino acid and carbohydrate transport and metabolism, cell wall/membrane/envelope biogenesis and signal transduction mechanism have been significantly changed during the dedifferentiation. To our best knowledge, this study is the first genome-wide effort to investigate the transcriptional changes in dedifferentiation immature maize embryos and the identified DEGs can serve as a basis for further functional characterization.

Genome expression profile analysis reveals important transcripts in maize roots responding to the stress of heavy metal Pb

Lead (Pb) has become one of the most abundant heavy metal pollutants of the environment. With its large biomass, maize could be an important object for studying the phytoremediation of Pb-contaminated soil. In our previous research, we screened 19 inbred lines of maize for Pb concentration, and line 178 was identified to be a hyperaccumulator for Pb in both the roots and aboveground parts. To identify important genes and metabolic pathways related to Pb accumulation and tolerance, line 178 was underwent genome expression profile under Pb stress and a control (CK). A total of approximately 11 million cDNA tags were sequenced and 4 665 539 and 4 936 038 clean tags were obtained from the libraries of the test and CK, respectively. In comparison to CK, 2379 and 1832 genes were identified up- or downregulated, respectively, more than fivefolds under Pb stress. Interestingly, all the genes were related to cellular processes and signaling, information storage and processing or metabolism functions. Particularly, the genes involved in posttranslational modification, protein turnover and chaperones; signal transduction, carbohydrate transport and metabolism; and lipid transport and metabolism significantly changed under the treatment. In addition, seven pathways including ribosome, photosynthesis, and carbon fixation were affected significantly, with 118, 12, 34, 21, 18, 72 and 43 differentially expressed genes involved. The significant upregulation of the ribosome pathway may reveal an important secret for Pb tolerance of line 178. And the sharp increase of laccase transcripts and metal ion transporters were suggested to account in part for Pb hyperaccumulation in the line.

Cytoplasmic male sterility-regulated novel microRNAs from maize

In higher plants, microRNA (miRNA) is involved in regulation of developmental processes, including sexual organ development. Seven novel miRNA families with one known miRNA were isolated by constructing a small RNA library from a mixture of anther from a cytoplasmic male sterile line and its maintainer. Two miRNAs are conserved in plant species. A total of 18 potential targets were identified for the eight miRNA families, including 15 proteins annotated with function and three unknown proteins. The known proteins include several proteins relevant to cell structure and stress response, transcription factors, and enzymes associated with metabolic and signaling pathways, playing important roles in microspore development. Quantitative real-time PCR assay revealed different expression patterns of the miRNAs between the cytoplasmic male sterile line and its maintainer. Each of the miRNAs tended to be down-regulated after the tetrad stage in a fertile line. However, most of the miRNAs in the cytoplasmic male sterile line were shown to be up-regulated from the tetrad to mononuclear stage, displaying special expression patterns differing from the ones in fertile line. We conclude that additional inactive miRNA pathways are essential during pollen development for a fertile line to ensure male fertility. Contrarily, miRNAs are up-regulated during the period from the tetrad to mononuclear stage, which contributes to pollen abortion for a cytoplasmic male sterile line.

Combined small RNA and degradome sequencing reveals microRNA regulation during immature maize embryo dedifferentiation

Genetic transformation of maize is highly dependent on the development of embryonic calli from the dedifferentiated immature embryo. To better understand the regulatory mechanism of immature embryo dedifferentiation, we generated four small RNA and degradome libraries from samples representing the major stages of dedifferentiation. More than 186 million raw reads of small RNA and degradome sequence data were generated. We detected 102 known miRNAs belonging to 23 miRNA families. In total, we identified 51, 70 and 63 differentially expressed miRNAs (DEMs) in the stage I, II, III samples, respectively, compared to the control. However, only 6 miRNAs were continually up-regulated by more than fivefold throughout the process of dedifferentiation. A total of 87 genes were identified as the targets of 21 DEM families. This group of targets was enriched in members of four significant pathways including plant hormone signal transduction, antigen processing and presentation, ECM-receptor interaction, and alpha-linolenic acid metabolism. The hormone signal transduction pathway appeared to be particularly significant, involving 21 of the targets. While the targets of the most significant DEMs have been proved to play essential roles in cell dedifferentiation. Our results provide important information regarding the regulatory networks that control immature embryo dedifferentiation in maize.

Cloning and characterization of miRNAs from maize seedling roots under low phosphorus stress

MicroRNAs (miRNAs) are a class of small, non-coding regulatory RNAs that regulate gene expression by guiding target mRNA cleavage or translational inhibition in plants and animals. In this study, a small RNA library was constructed to identify conserved miRNAs as well as novel miRNAs in maize seedling roots under low level phosphorus stress. Twelve miRNAs were identified by high throughput sequencing of the library and subsequent analysis, two belong to conserved miRNA families (miRNA399b and miRNA156), and the remaining ten are novel and one of latter is conserved in gramineous species. Based on sequence homology, we predicted 125 potential target genes of these miRNAs and then expression patterns of 7 miRNAs were validated by semi-RT-PCR analysis. MiRNA399b, Zma-miR3, and their target genes (Zmpt1 and Zmpt2) were analyzed by real-time PCR. It is shown that both miRNA399b and Zma-miR3 are induced by low phosphorus stress and regulated by their target genes (Zmpt1 and Zmpt2). Moreover, Zma-miR3, regulated by two maize inorganic phosphate transporters as a newly identified miRNAs, would likely be directly involved in phosphate homeostasis, so was miRNA399b in Arabidopsis and rice. These results indicate that both conserved and maize-specific miRNAs play important roles in stress responses and other physiological processes correlated with phosphate starvation, regulated by their target genes. Identification of these differentially expressed miRNAs will facilitate us to uncover the molecular mechanisms underlying the progression of maize seedling roots development under low level phosphorus stress.

The Dynamics of DNA Methylation in Maize Roots under Pb Stress

Plants adapt to adverse conditions through a series of physiological, cellular, and molecular processes, culminating in stress tolerance. However, little is known about the associated regulatory mechanisms at the epigenetic level in maize under lead (Pb) stress. Therefore, in this study, we aimed to compare DNA methylation profiles during the dynamic development of maize roots following Pb treatment to identify candidate genes involved in the response to Pb stress. Methylated DNA immunoprecipitation-sequencing (MeDIP-seq) was used to investigate the genome-wide DNA methylation patterns in maize roots under normal condition (A1) and 3 mM Pb(NO3)2 stress for 12 h (K2), 24 h (K3) and 48 h (K4). The results showed that the average methylation density was the highest in CpG islands (CGIs), followed by the intergenic regions. Within the gene body, the methylation density of the introns was higher than those of the UTRs and exons. In total, 3857 methylated genes were found in 4 tested samples, including 1805 differentially methylated genes for K2 versus A1, 1508 for K3 versus A1, and 1660 for K4 versus A1. Further analysis showed that 140 genes exhibited altered DNA methylation in all three comparisons, including some well-known stress-responsive transcription factors and proteins, such as MYB, AP2/ERF, bZIP, serine-threonine/tyrosine-proteins, pentatricopeptide repeat proteins, RING zinc finger proteins, F-box proteins, leucine-rich repeat proteins and tetratricopeptide repeat proteins. This study revealed the genome-scale DNA methylation patterns of maize roots in response to Pb exposure and identified candidate genes that potentially regulate root dynamic development under Pb stress at the methylation level.

Transcriptomic changes during maize roots development responsive to Cadmium (Cd) pollution using comparative RNAseq-based approach

The heavy metal cadmium (Cd), acts as a widespread environmental contaminant, which has shown to adversely affect human health, food safety and ecosystem safety in recent years. However, research on how plant respond to various kinds of heavy metal stress is scarcely reported, especially for understanding of complex molecular regulatory mechanisms and elucidating the gene networks of plant respond to Cd stress. Here, transcriptomic changes during Mo17 and B73 seedlings development responsive to Cd pollution were investigated and comparative RNAseq-based approach in both genotypes were performed. 115 differential expression genes (DEGs) with significant alteration in expression were found co-modulated in both genotypes during the maize seedling development; of those, most of DGEs were found comprised of stress and defense responses proteins, transporters, as well as transcription factors, such as thaumatin-like protein, ZmOPR2 and ZmOPR5. More interestingly, genotype-specific transcriptional factors changes induced by Cd stress were found contributed to the regulatory mechanism of Cd sensitivity in both different genotypes. Moreover, 12 co-expression modules associated with specific biological processes or pathways (M1 to M12) were identified by consensus co-expression network. These results will expand our understanding of complex molecular mechanism of response and defense to Cd exposure in maize seedling roots.

Heterosis in Early Maize Ear Inflorescence Development: A Genome-Wide Transcription Analysis for Two Maize Inbred Lines and Their Hybrid

Heterosis, or hybrid vigor, contributes to superior agronomic performance of hybrids compared to their inbred parents. Despite its importance, little is known about the genetic and molecular basis of heterosis. Early maize ear inflorescences formation affects grain yield, and are thus an excellent model for molecular mechanisms involved in heterosis. To determine the parental contributions and their regulation during maize ear-development-genesis, we analyzed genome-wide digital gene expression profiles in two maize elite inbred lines (B73 and Mo17) and their F1 hybrid using deep sequencing technology. Our analysis revealed 17,128 genes expressed in these three genotypes and 22,789 genes expressed collectively in the present study. Approximately 38% of the genes were differentially expressed in early maize ear inflorescences from heterotic cross, including many transcription factor genes and some presence/absence variations (PAVs) genes, and exhibited multiple modes of gene action. These different genes showing differential expression patterns were mainly enriched in five cellular component categories (organelle, cell, cell part, organelle part and macromolecular complex), five molecular function categories (structural molecule activity, binding, transporter activity, nucleic acid binding transcription factor activity and catalytic activity), and eight biological process categories (cellular process, metabolic process, biological regulation, regulation of biological process, establishment of localization, cellular component organization or biogenesis, response to stimulus and localization). Additionally, a significant number of genes were expressed in only one inbred line or absent in both inbred lines. Comparison of the differences of modes of gene action between previous studies and the present study revealed only a small number of different genes had the same modes of gene action in both maize seedlings and ear inflorescences. This might be an indication that in different tissues or developmental stages, different global expression patterns prevail, which might nevertheless be related to heterosis. Our results support the hypotheses that multiple molecular mechanisms (dominance and overdominance modes) contribute to heterosis.

Genome-Wide Identification and Analysis of Drought-Responsive Genes and MicroRNAs in Tobacco

Drought stress response is a complex trait regulated at transcriptional and post-transcriptional levels in tobacco. Since the 1990s, many studies have shown that miRNAs act in many ways to regulate target expression in plant growth, development and stress response. The recent draft genome sequence of Nicotiana benthamiana has provided a framework for Digital Gene Expression (DGE) and small RNA sequencing to understand patterns of transcription in the context of plant response to environmental stress. We sequenced and analyzed three Digital Gene Expression (DGE) libraries from roots of normal and drought-stressed tobacco plants, and four small RNA populations from roots, stems and leaves of control or drought-treated tobacco plants, respectively. We identified 276 candidate drought responsive genes (DRGs) with sequence similarities to 64 known DRGs from other model plant crops, 82 were transcription factors (TFs) including WRKY, NAC, ERF and bZIP families. Of these tobacco DRGs, 54 differentially expressed DRGs included 21 TFs, which belonged to 4 TF families such as NAC (6), MYB (4), ERF (10), and bZIP (1). Additionally, we confirmed expression of 39 known miRNA families (122 members) and five conserved miRNA families, which showed differential regulation under drought stress. Targets of miRNAs were further surveyed based on a recently published study, of which ten targets were DRGs. An integrated gene regulatory network is proposed for the molecular mechanisms of tobacco root response to drought stress using differentially expressed DRGs, the changed expression profiles of miRNAs and their target transcripts. This network analysis serves as a reference for future studies on tobacco response stresses such as drought, cold and heavy metals.