ZmWAK02 encoding an RD-WAK protein confers maize resistance against gray leaf spot.

Gray leaf spot (GLS) caused by Cercospora zeina or C. zeae-maydis is a major maize disease throughout the world. Although more than 100 QTLs resistant against GLS have been identified, very few of them have been cloned. Here, we identified a major resistance QTL against GLS, qRglsSB, explaining 58.42% phenotypic variation in SB12×SA101 BC1 F1 population. By fine-mapping, it was narrowed down into a 928 kb region. By using transgenic lines, mutants and complementation lines, it was confirmed that the ZmWAK02 gene, encoding an RD wall-associated kinase, is the responsible gene in qRglsSB resistant against GLS. The introgression of the ZmWAK02 gene into hybrid lines significantly improves their grain yield in the presence of GLS pressure and does not reduce their grain yield in the absence of GLS. In summary, we cloned a gene, ZmWAK02, conferring large effect of GLS resistance and confirmed its great value in maize breeding.

A Novel Soybean Diacylglycerol Acyltransferase 1b Variant with Three Amino Acid Substitutions Increases Seed Oil Content.

Improving soybean (Glycine max) seed composition by increasing the protein and oil components will add significant value to the crop and enhance environmental sustainability. Diacylglycerol acyltransferase (DGAT) catalyzes the final rate-limiting step in triacylglycerol (TAG) biosynthesis and has a major impact on seed oil accumulation. We previously identified a soybean DGAT1b variant with 14 amino acid substitutions (GmDGAT1b-MOD) that increases total oil content by 3 percentage points when overexpressed in soybean seeds. In the present study, additional GmDGAT1b variants were generated to further increase oil with a reduced number of substitutions. Variants with one to four amino acid substitutions were screened in the model systems S. cerevisiae and transient N. benthamiana leaf. Promising GmDGAT1b variants resulting in high oil accumulation in the model systems were selected for over-expression in soybeans. One GmDGAT1b variant with three novel amino acid substitutions (GmDGAT1b-3aa) increased total soybean oil to levels near the previously discovered GmDGAT1b-MOD variant. In a multiple location field trial, GmDGAT1b-3aa transgenic events had significantly increased oil and protein by up to 2.3 and 0.6 percentage points, respectively. Modeling of the GmDGAT1b-3aa protein structure provided insights into the potential function of the three substitutions. These findings will guide efforts to improve soybean oil content and overall seed composition by CRISPR editing.

Parallel tuning of semi-dwarfism via differential splicing of Brachytic1 in commercial maize and smallholder sorghum.

In the current genomic era, the search and deployment of new semi-dwarf alleles have continued to develop better plant types in all cereals. We characterized an agronomically optimal semi-dwarf mutation in Zea mays L. and a parallel polymorphism in Sorghum bicolor L. We cloned the maize brachytic1 (br1-Mu) allele by a modified PCR-based Sequence Amplified Insertion Flanking Fragment (SAIFF) approach. Histology and RNA-Seq elucidated the mechanism of semi-dwarfism. GWAS linked a sorghum plant height QTL with the Br1 homolog by resequencing a West African sorghum landraces panel. The semi-dwarf br1-Mu allele encodes an MYB transcription factor78 that positively regulates stalk cell elongation by interacting with the polar auxin pathway. Semi-dwarfism is due to differential splicing and low functional Br1 wild-type transcript expression. The sorghum ortholog, SbBr1, co-segregates with the major plant height QTL qHT7.1 and is alternatively spliced. The high frequency of the Sbbr1 allele in African landraces suggests that African smallholder farmers used the semi-dwarf allele to improve plant height in sorghum long before efforts to introduce Green Revolution-style varieties in the 1960s. Surprisingly, variants for differential splicing of Brachytic1 were found in both commercial maize and smallholder sorghum, suggesting parallel tuning of plant architecture across these systems.

The NLRomes of Zea mays NAM founder lines and Zea luxurians display presence-absence variation, integrated domain diversity, and mobility.

Plant pathogens cause significant crop loss worldwide, and new resistance genes deployed to combat diseases can be overcome quickly. Understanding the existing resistance gene diversity within the germplasm of major crops, such as maize, is crucial for the development of new disease-resistant varieties. We analysed the nucleotide-binding leucine-rich repeat receptors (NLRs) of 26 recently sequenced diverse founder lines from the maize nested association mapping (NAM) population and compared them to the R gene complement present in a wild relative of maize, Zea luxurians. We found that NLRs in both species contain a large diversity of atypical integrated domains, including many domains that have not previously been found in the NLRs of other species. Additionally, the single Z. luxurians genome was found to have greater integrated atypical domain diversity than all 26 NAM founder lines combined, indicating that this species may represent a rich source of novel resistance genes. NLRs were also found to have very high sequence diversity and presence-absence variation among the NAM founder lines, with a large NLR cluster on Chr10 representing a diversity hotspot. Additionally, NLRs were shown to be mobile within maize genomes, with several putative interchromosomal translocations identified.

The northern corn leaf blight resistance gene Ht1 encodes an nucleotide-binding, leucine-rich repeat immune receptor.

Northern corn leaf blight, caused by the fungal pathogen Exserohilum turcicum, is a major disease of maize. The first major locus conferring resistance to E. turcicum race 0, Ht1, was identified over 50 years ago, but the underlying gene has remained unknown. We employed a map-based cloning strategy to identify the Ht1 causal gene, which was found to be a coiled-coil nucleotide-binding, leucine-rich repeat (NLR) gene, which we named PH4GP-Ht1. Transgenic testing confirmed that introducing the native PH4GP-Ht1 sequence to a susceptible maize variety resulted in resistance to E. turcicum race 0. A survey of the maize nested association mapping genomes revealed that susceptible Ht1 alleles had very low to no expression of the gene. Overexpression of the susceptible B73 allele, however, did not result in resistant plants, indicating that sequence variations may underlie the difference between resistant and susceptible phenotypes. Modelling of the PH4GP-Ht1 protein indicated that it has structural homology to the Arabidopsis NLR resistance gene ZAR1, and probably forms a similar homopentamer structure following activation. RNA sequencing data from an infection time course revealed that 1 week after inoculation there was a threefold reduction in fungal biomass in the PH4GP-Ht1 transgenic plants compared to wild-type plants. Furthermore, PH4GP-Ht1 transgenics had significantly more inoculation-responsive differentially expressed genes than wild-type plants, with enrichment seen in genes associated with both defence and photosynthesis. These results demonstrate that the NLR PH4GP-Ht1 is the causal gene underlying Ht1, which represents a different mode of action compared to the previously reported wall-associated kinase northern corn leaf blight resistance gene Htn1/Ht2/Ht3.

The RppC-AvrRppC NLR-effector interaction mediates the resistance to southern corn rust in maize.

Southern corn rust (SCR), caused by the fungal pathogen Puccinia polysora, is a major threat to maize production worldwide. Efficient breeding and deployment of resistant hybrids are key to achieving durable control of SCR. Here, we report the molecular cloning and characterization of RppC, which encodes an NLR-type immune receptor and is responsible for a major SCR resistance quantitative trait locus. Furthermore, we identified the corresponding avirulence effector, AvrRppC, which is secreted by P. polysora and triggers RppC-mediated resistance. Allelic variation of AvrRppC directly determines the effectiveness of RppC-mediated resistance, indicating that monitoring of AvrRppC variants in the field can guide the rational deployment of RppC-containing hybrids in maize production. Currently, RppC is the most frequently deployed SCR resistance gene in China, and a better understanding of its mode of action is critical for extending its durability.

Developmental and transcriptional responses of maize to drought stress under field conditions.

Drought is a common abiotic stress which significantly limits global crop productivity. Maize is an important staple crop and its yield is determined by successful development of the female inflorescence, the ear. We investigated drought stress responses across several developmental stages of the maize B73 inbred line under field conditions. Drought suppressed plant growth, but had little impact on progression through developmental stages. While ear growth was suppressed by drought, the process of spikelet initiation was not significantly affected. Tassel growth was reduced to a lesser extent compared to the observed reduction in ear growth under stress. Parallel RNA-seq profiling of leaves, ears, and tassels at several developmental stages revealed tissue-specific differences in response to drought stress. High temperature fluctuation was an additional environmental factor that also likely influenced gene expression patterns in the field. Drought induced significant transcriptional changes in leaves and ears but only minor changes in the tassel. Additionally, more genes were drought responsive in ears compared to leaves over the course of drought treatment. Genes that control DNA replication, cell cycle, and cell division were significantly down-regulated in stressed ears, which was consistent with inhibition of ear growth under drought. Inflorescence meristem genes were affected by drought to a lesser degree which was consistent with the minimal impact of drought on spikelet initiation. In contrast, genes that are involved in floret and ovule development were sensitive to stress, which is consistent with the detrimental effect of drought on gynoecium development and kernel set.

Genome-Wide Analysis of Alternative Splicing during Development and Drought Stress in Maize.

Alternative splicing plays a crucial role in plant development as well as stress responses. Although alternative splicing has been studied during development and in response to stress, the interplay between these two factors remains an open question. To assess the effects of drought stress on developmentally regulated splicing in maize (Zea mays), 94 RNA-seq libraries from ear, tassel, and leaf of the B73 public inbred line were constructed at four developmental stages under both well-watered and drought conditions. This analysis was supplemented with a publicly available series of 53 libraries from developing seed, embryo, and endosperm. More than 48,000 novel isoforms, often with stage- or condition-specific expression, were uncovered, suggesting that developmentally regulated alternative splicing occurs in thousands of genes. Drought induced large developmental splicing changes in leaf and ear but relatively few in tassel. Most developmental stage-specific splicing changes affected by drought were tissue dependent, whereas stage-independent changes frequently overlapped between leaf and ear. A linear relationship was found between gene expression changes in splicing factors and alternative spicing of other genes during development. Collectively, these results demonstrate that alternative splicing is strongly associated with tissue type, developmental stage, and stress condition.

Transcriptome analysis of near-isogenic lines provides molecular insights into starch biosynthesis in maize kernel.

Starch is the major component in maize kernels, providing a stable carbohydrate source for humans and livestock as well as raw material for the biofuel industry. Increasing maize kernel starch content will help meet industry demands and has the potential to increase overall yields. We developed a pair of maize near-isogenic lines (NILs) with different alleles for a starch quantitative trait locus on chromosome 3 (qHS3), resulting in different kernel starch content. To investigate the candidate genes for qHS3 and elucidate their effects on starch metabolism, RNA-Seq was performed for the developing kernels of the NILs at 14 and 21 d after pollination (DAP). Analysis of genomic and transcriptomic data identified 76 genes with nonsynonymous single nucleotide polymorphisms and 384 differentially expressed genes (DEGs) in the introgressed fragment, including a hexokinase gene, ZmHXK3a, which catalyzes the conversion of glucose to glucose-6-phosphate and may play a key role in starch metabolism. The expression pattern of all DEGs in starch metabolism shows that altered expression of the candidate genes for qHS3 promoted starch synthesis, with positive consequences for kernel starch content. These results expand the current understanding of starch biosynthesis and accumulation in maize kernels and provide potential candidate genes to increase starch content.

Differential expression of miRNAs and their target genes in senescing leaves and siliques: insights from deep sequencing of small RNAs and cleaved target RNAs.

MicroRNAs (miRNAs) are a class of small RNAs, which typically function by guiding cleavage of target mRNAs. They are known to play roles in a variety of plant processes including development, responses to environmental stresses and senescence. To identify senescence regulation of miRNAs in Arabidopsis thaliana, eight small RNA libraries were constructed and sequenced at four different stages of development and senescence from both leaves and siliques, resulting in more than 200 million genome-matched sequences. Parallel analysis of RNA ends libraries, which enable the large-scale examination of miRNA-guided cleavage products, were constructed and sequenced, resulting in over 750 million genome-matched sequences. These large datasets led to the identification a new senescence-inducible small RNA locus, as well as new regulation of known miRNAs and their target genes during senescence, many of which have established roles in nutrient responsiveness and cell structural integrity. In keeping with remobilization of nutrients thought to occur during senescence, many miRNAs and targets had opposite expression pattern changes between leaf and silique tissues during the progression of senescence. Taken together, these findings highlight the integral role that miRNAs may play in the remobilization of resources and alteration of cellular structure that is known to occur in senescence.

Genome-wide analysis of alternative splicing in Zea mays: landscape and genetic regulation.

Alternative splicing enhances transcriptome diversity in all eukaryotes and plays a role in plant tissue identity and stress adaptation. To catalog new maize (Zea mays) transcripts and identify genomic loci that regulate alternative splicing, we analyzed over 90 RNA-seq libraries from maize inbred lines B73 and Mo17, as well as Syn10 doubled haploid lines (progenies from B73 × Mo17). Transcript discovery was augmented with publicly available data from 14 maize tissues, expanding the maize transcriptome by more than 30,000 and increasing the percentage of intron-containing genes that undergo alternative splicing to 40%. These newly identified transcripts greatly increase the diversity of the maize proteome, sometimes coding for entirely different proteins compared with their most similar annotated isoform. In addition to increasing proteome diversity, many genes encoding novel transcripts gained an additional layer of regulation by microRNAs, often in a tissue-specific manner. We also demonstrate that the majority of genotype-specific alternative splicing can be genetically mapped, with cis-acting quantitative trait loci (QTLs) predominating. A large number of trans-acting QTLs were also apparent, with nearly half located in regions not shown to contain genes associated with splicing. Taken together, these results highlight the currently underappreciated role that alternative splicing plays in tissue identity and genotypic variation in maize.

Comprehensive investigation of microRNAs enhanced by analysis of sequence variants, expression patterns, ARGONAUTE loading, and target cleavage.

MicroRNAs (miRNAs) are a class of small RNAs that typically function by guiding the cleavage of target messenger RNAs. They have been shown to play major roles in a variety of plant processes, including development, and responses to pathogens and environmental stresses. To identify new miRNAs and regulation in Arabidopsis (Arabidopsis thaliana), 27 small RNA libraries were constructed and sequenced from various tissues, stresses, and small RNA biogenesis mutants, resulting in 95 million genome-matched sequences. The use of rdr2 to enrich the miRNA population greatly enhanced this analysis and led to the discovery of new miRNAs arising from both known and new precursors, increasing the total number of Arabidopsis miRNAs by about 10%. Parallel Analysis of RNA Ends data provide evidence that the majority guide target cleavage. Many libraries represented novel stress/tissue conditions, such as submergence-stressed flowers, which enabled the identification of new stress regulation of both miRNAs and their targets, all of which were validated in wild-type plants. By combining small RNA expression analysis with ARGONAUTE immunoprecipitation data and global target cleavage data from Parallel Analysis of RNA Ends, a much more complete picture of Arabidopsis miRNAs was obtained. In particular, the discovery of ARGONAUTE loading and target cleavage biases gave important insights into tissue-specific expression patterns, pathogen responses, and the role of sequence variation among closely related miRNA family members that would not be evident without this combinatorial approach.

Abiotic stress-associated miRNAs: detection and functional analysis.

MicroRNAs (miRNAs) are small regulatory noncoding RNAs varying in length between 20 and 24 nucleotides. They play a key role during plant development by negatively regulating gene expression at the posttranscriptional level. Moreover, recent studies reported several miRNAs associated with abiotic stress responses. Small RNA cloning and high-throughput deep sequencing methods provide expression profiles of not only known miRNAs, but also novel miRNAs. In this chapter, we describe the methods used to identify and characterize abiotic stress-associated miRNAs and their target genes.

Filtering of deep sequencing data reveals the existence of abundant Dicer-dependent small RNAs derived from tRNAs.

Deep sequencing technologies such as Illumina, SOLiD, and 454 platforms have become very powerful tools in discovering and quantifying small RNAs in diverse organisms. Sequencing small RNA fractions always identifies RNAs derived from abundant RNA species such as rRNAs, tRNAs, snRNA, and snoRNA, and they are widely considered to be random degradation products. We carried out bioinformatic analysis of deep sequenced HeLa RNA and after quality filtering, identified highly abundant small RNA fragments, derived from mature tRNAs that are likely produced by specific processing rather than from random degradation. Moreover, we showed that the processing of small RNAs derived from tRNA(Gln) is dependent on Dicer in vivo and that Dicer cleaves the tRNA in vitro.

Genome-wide analysis for discovery of rice microRNAs reveals natural antisense microRNAs (nat-miRNAs).

Small RNAs (21-24 nt) are involved in gene regulation through translation inhibition, mRNA cleavage, or directing chromatin modifications. In rice, currently approximately 240 microRNAs (miRNAs) have been annotated. We sequenced more than four million small RNAs from rice and identified another 24 miRNA genes. Among these, we found a unique class of miRNAs that derive from natural cis-antisense transcript pairs. This configuration generates miRNAs that can perfectly match their targets. We provide evidence that the miRNAs function by inducing mRNA cleavage in the middle of their complementary site. Their production requires Dicer-like 1 (DCL1) activity, which is essential for canonical miRNA biogenesis. All of the natural antisense miRNAs (nat-miRNAs) identified in this study have large introns in their precursors that appear critical for nat-miRNA evolution and for the formation of functional miRNA loci. These findings suggest that other natural cis-antisense loci with similar exon-intron arrangements could be another source of miRNA genes.