Huiru Peng

Diverse set of microRNAs are responsive to powdery mildew infection and heat stress in wheat (Triticum aestivum L.)

BackgroundMicroRNAs (miRNAs) are a class of small non-coding regulatory RNAs that regulate gene expression by guiding target mRNA cleavage or translational inhibition. MiRNAs can have large-scale regulatory effects on development and stress response in plants.ResultsTo test whether miRNAs play roles in regulating response to powdery mildew infection and heat stress in wheat, by using Solexa high-throughput sequencing we cloned the small RNA from wheat leaves infected by preponderant physiological strain Erysiphe graminis f. sp. tritici (Egt) or by heat stress treatment. A total of 153 miRNAs were identified, which belong to 51 known and 81 novel miRNA families. We found that 24 and 12 miRNAs were responsive to powdery mildew infection and heat stress, respectively. We further predicted that 149 target genes were potentially regulated by the novel wheat miRNA.ConclusionsOur results indicated that diverse set of wheat miRNAs were responsive to powdery mildew infection and heat stress and could function in wheat responses to both biotic and abiotic stresses.

Identification and characterization of wheat long non-protein coding RNAs responsive to powdery mildew infection and heat stress by using microarray analysis and SBS sequencing

BackgroundBiotic and abiotic stresses, such as powdery mildew infection and high temperature, are important limiting factors for yield and grain quality in wheat production. Emerging evidences suggest that long non-protein coding RNAs (npcRNAs) are developmentally regulated and play roles in development and stress responses of plants. However, identification of long npcRNAs is limited to a few plant species, such as Arabidopsis, rice and maize, no systematic identification of long npcRNAs and their responses to abiotic and biotic stresses is reported in wheat.ResultsIn this study, by using computational analysis and experimental approach we identified 125 putative wheat stress responsive long npcRNAs, which are not conserved among plant species. Among them, some were precursors of small RNAs such as microRNAs and siRNAs, two long npcRNAs were identified as signal recognition particle (SRP) 7S RNA variants, and three were characterized as U3 snoRNAs. We found that wheat long npcRNAs showed tissue dependent expression patterns and were responsive to powdery mildew infection and heat stress.ConclusionOur results indicated that diverse sets of wheat long npcRNAs were responsive to powdery mildew infection and heat stress, and could function in wheat responses to both biotic and abiotic stresses, which provided a starting point to understand their functions and regulatory mechanisms in the future.

TamiR159 Directed Wheat TaGAMYB Cleavage and Its Involvement in Anther Development and Heat Response

In Arabidopsis and rice, miR159-regulated GAMYB-like family transcription factors function in flower development and gibberellin (GA) signaling in cereal aleurone cells. In this study, the involvement of miR159 in the regulation of its putative target TaGAMYB and its relationship to wheat development were investigated. First, we demonstrated that cleavage of TaGAMYB1 and TaGAMYB2 was directed by miR159 using 5′-RACE and a transient expression system. Second, we overexpressed TamiR159, TaGAMYB1 and mTaGAMYB1 (impaired in the miR159 binding site) in transgenic rice, revealing that the accumulation in rice of mature miR159 derived from the precursor of wheat resulted in delayed heading time and male sterility. In addition, the number of tillers and primary branches in rice overexpressing mTaGAMYB1 increased relative to the wild type. Our previous study reported that TamiR159 was downregulated after two hours of heat stress treatment in wheat (Triticum aestivum L.). Most notably, the TamiR159 overexpression rice lines were more sensitive to heat stress relative to the wild type, indicating that the downregulation of TamiR159 in wheat after heat stress might participate in a heat stress-related signaling pathway, in turn contributing to heat stress tolerance.

Whole-genome discovery of miRNAs and their targets in wheat ( Triticum aestivum L.)

BackgroundMicroRNAs (miRNAs) are small, non-coding RNAs playing essential roles in plant growth, development, and stress responses. Sequencing of small RNAs is a starting point for understanding their number, diversity, expression and possible roles in plants.ResultsIn this study, we conducted a genome-wide survey of wheat miRNAs from 11 tissues, characterizing a total of 323 novel miRNAs belonging to 276 families in wheat. A miRNA conservation analysis identified 191 wheat-specific miRNAs, 2 monocot-specific miRNAs, and 30 wheat-specific variants from 9 highly conserved miRNA families. To understand possible roles of wheat miRNAs, we determined 524 potential targets for 124 miRNA families through degradome sequencing, and cleavage of a subset of them was validated via 5′ RACE. Based on the genome-wide identification and characterization of miRNAs and their associated target genes, we further identified 64 miRNAs preferentially expressing in developing or germinating grains, which could play important roles in grain development.ConclusionWe discovered 323 wheat novel miRNAs and 524 target genes for 124 miRNA families in a genome-wide level, and our data will serve as a foundation for future research into the functional roles of miRNAs in wheat.

Transcriptome comparison of susceptible and resistant wheat in response to powdery mildew infection.

Powdery mildew (Pm) caused by the infection of Blumeria graminis f. sp. tritici (Bgt) is a worldwide crop disease resulting in significant loss of wheat yield. To profile the genes and pathways responding to the Bgt infection, here, using Affymetrix wheat microarrays, we compared the leaf transcriptomes before and after Bgt inoculation in two wheat genotypes, a Pm-susceptible cultivar Jingdong 8 (S) and its near-isogenic line (R) carrying a single Pm resistant gene Pm30. Our analysis showed that the original gene expression status in the S and R genotypes of wheat was almost identical before Bgt inoculation, since only 60 genes exhibited differential expression by P = 0.01 cutoff. However, 12 h after Bgt inoculation, 3014 and 2800 genes in the S and R genotype, respectively, responded to infection. A wide range of pathways were involved, including cell wall fortification, flavonoid biosynthesis and metabolic processes. Furthermore, for the first time, we show that sense-antisense pair genes might be participants in wheat-powdery mildew interaction. In addition, the results of qRT-PCR analysis on several candidate genes were consistent with the microarray data in their expression patterns. In summary, this study reveals leaf transcriptome changes before and after powdery mildew infection in wheat near-isogenic lines, suggesting that powdery mildew resistance is a highly complex systematic response involving a large amount of gene regulation.

Comparative proteomic analysis of embryos between a maize hybrid and its parental lines during early stages of seed germination.

In spite of commercial use of heterosis in agriculture, the molecular basis of heterosis is poorly understood. It was observed that maize hybrid Zong3/87-1 exhibited an earlier onset or heterosis in radicle emergence. To get insights into the underlying mechanism of heterosis in radicle emergence, differential proteomic analysis between hybrid and its parental lines was performed. In total, the number of differentially expressed protein spots between hybrid and its parental lines in dry and 24 h imbibed seed embryos were 134 and 191, respectively, among which 47.01% (63/134) and 34.55% (66/191) protein spots displayed nonadditively expressed pattern. Remarkably, 54.55% of nonadditively accumulated proteins in 24 h imbibed seed embryos displayed above or equal to the level of the higher parent patterns. Moreover, 155 differentially expressed protein spots were identified, which were grouped into eight functional classes, including transcription & translation, energy & metabolism, signal transduction, disease & defense, storage protein, transposable element, cell growth & division and unclassified proteins. In addition, one of the upregulated proteins in F1 hybrids was ZmACT2, a homolog of Arabidopsis thaliana ACT7 (AtACT7). Expressing ZmACT2 driven by the AtACT7 promoter partially complemented the low germination phenotype in the Atact7 mutant. These results indicated that hybridization between two parental lines can cause changes in the expression of a variety of proteins, and it is concluded that the altered pattern of gene expression at translational level in the hybrid may be responsible for the observed heterosis.

High-Density Genetic Linkage Map Construction and QTL Mapping of Grain Shape and Size in the Wheat Population Yanda1817 × Beinong6

High-density genetic linkage maps are necessary for precisely mapping quantitative trait loci (QTLs) controlling grain shape and size in wheat. By applying the Infinium iSelect 9K SNP assay, we have constructed a high-density genetic linkage map with 269 F 8 recombinant inbred lines (RILs) developed between a Chinese cornerstone wheat breeding parental line Yanda1817 and a high-yielding line Beinong6. The map contains 2431 SNPs and 128 SSR & EST-SSR markers in a total coverage of 3213.2 cM with an average interval of 1.26 cM per marker. Eighty-eight QTLs for thousand-grain weight (TGW), grain length (GL), grain width (GW) and grain thickness (GT) were detected in nine ecological environments (Beijing, Shijiazhuang and Kaifeng) during five years between 2010–2014 by inclusive composite interval mapping (ICIM) (LOD≥2.5). Among which, 17 QTLs for TGW were mapped on chromosomes 1A, 1B, 2A, 2B, 3A, 3B, 3D, 4A, 4D, 5A, 5B and 6B with phenotypic variations ranging from 2.62% to 12.08%. Four stable QTLs for TGW could be detected in five and seven environments, respectively. Thirty-two QTLs for GL were mapped on chromosomes 1B, 1D, 2A, 2B, 2D, 3B, 3D, 4A, 4B, 4D, 5A, 5B, 6B, 7A and 7B, with phenotypic variations ranging from 2.62% to 44.39%. QGl.cau-2A.2 can be detected in all the environments with the largest phenotypic variations, indicating that it is a major and stable QTL. For GW, 12 QTLs were identified with phenotypic variations range from 3.69% to 12.30%. We found 27 QTLs for GT with phenotypic variations ranged from 2.55% to 36.42%. In particular, QTL QGt.cau-5A.1 with phenotypic variations of 6.82–23.59% was detected in all the nine environments. Moreover, pleiotropic effects were detected for several QTL loci responsible for grain shape and size that could serve as target regions for fine mapping and marker assisted selection in wheat breeding programs.

Epigenetic modification contributes to the expression divergence of three TaEXPA1 homoeologs in hexaploid wheat (Triticum aestivum)

Common wheat is a hexaploid species with most of the genes present as triplicate homoeologs. Expression divergences of homoeologs are frequently observed in wheat, as well as in other polyploid plants. However, the mechanisms underlying this phenomenon are poorly understood. Expansin genes play important roles in the regulation of cell size, as well as organ size. We found that all three TaEXPA1 homoeologs were silenced in seedling roots. In seedling leaves, TaEXPA1-A and TaEXPA1-D were expressed, but TaEXPA1-B was silenced. Further analysis revealed that silencing of TaEXPA1-B in leaves occurred after the formation of the hexaploid. Chromatin immunoprecipitation assays revealed that the transcriptional silencing of three TaEXPA1 homoeologs in roots was correlated with an increased level of H3K9 dimethylation and decreased levels of H3K4 trimethylation and H3K9 acetylation. Reactivation of TaEXPA1-A and TaEXPA1-D expression in leaves was correlated with increased levels of H3K4 trimethylation and H3K9 acetylation, and decreased levels of H3K9 dimethylation in their promoters, respectively. Moreover, a higher level of cytosine methylation was detected in the promoter region of TaEXPA1-B, which may contribute to its silencing in leaves. We demonstrated that epigenetic modifications contribute to the expression divergence of three TaEXPA1 homoeologs during wheat development.

Genome-wide identification of gibberellins metabolic enzyme genes and expression profiling analysis during seed germination in maize.

Gibberellin (GA) is an essential phytohormone that controls many aspects of plant development. To enhance our understanding of GA metabolism in maize, we intensively screened and identified 27 candidate genes encoding the seven GA metabolic enzymes including ent-copalyl diphosphate synthase (CPS), ent-kaurene synthase (KS), ent-kaurene oxidase (KO), ent-kaurenoic acid oxidase (KAO), GA 20-oxidase (GA20ox), GA 3-oxidase (GA3ox), and GA 2-oxidase (GA2ox), using all available public maize databases. The results indicate that maize genome contains three CPS, four KS, two KO and one KAO genes, and most of them are arranged separately on the maize genome, which differs from that in rice. In addition, the enzymes catalyzing the later steps (ZmGA20ox, ZmGA3ox and ZmGA2ox) are also encoded by gene families in maize, but GA3ox enzyme is likely to be encoded by single gene. Expression profiling analysis exhibited that transcripts of 15 GA metabolic genes could be detected during maize seed germination, which provides further evidence for the notion that increased synthesis of active GA in the embryo is required for triggering germination events. Moreover, a variety of temporal genes expression patterns of GA metabolic genes were detected, which revealed the complexity of underlying mechanism for GA regulated seed germination.

Comparative high-resolution mapping of the wax inhibitors Iw1 and Iw2 in hexaploid wheat.

The wax (glaucousness) on wheat leaves and stems is mainly controlled by two sets of genes: glaucousness loci (W1 and W2) and non-glaucousness loci (Iw1 and Iw2). The non-glaucousness (Iw) loci act as inhibitors of the glaucousness loci (W). High-resolution comparative genetic linkage maps of the wax inhibitors Iw1 originating from Triticum dicoccoides, and Iw2 from Aegilops tauschii were developed by comparative genomics analyses of Brachypodium, sorghum and rice genomic sequences corresponding to the syntenic regions of the Iw loci in wheat. Eleven Iw1 and eight Iw2 linked EST markers were developed and mapped to linkage maps on the distal regions of chromosomes 2BS and 2DS, respectively. The Iw1 locus mapped within a 0.96 cM interval flanked by the BE498358 and CA499581 EST markers that are collinear with 122 kb, 202 kb, and 466 kb genomic regions in the Brachypodium 5S chromosome, the sorghum 6S chromosome and the rice 4S chromosome, respectively. The Iw2 locus was located in a 4.1 to 5.4-cM interval in chromosome 2DS that is flanked by the CJ886319 and CJ519831 EST markers, and this region is collinear with a 2.3 cM region spanning the Iw1 locus on chromosome 2BS. Both Iw1 and Iw2 co-segregated with the BF474014 and CJ876545 EST markers, indicating they are most likely orthologs on 2BS and 2DS. These high-resolution maps can serve as a framework for chromosome landing, physical mapping and map-based cloning of the wax inhibitors in wheat.