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 3 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.

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.After transcription, the majority of eukaryotic premRNA is subjected to a series of posttranscriptional modifications, including the removal of introns to form a mature mRNA (Stamm et al., 2005). Although some introns are removed constitutively, many can be processed in a variety of alternative ways, including exon skipping, intron retention, alternative acceptor, alternative donor, and alternative position (change in both acceptor and donor positions; Lorkovic et al., 2000). These alternative splicing events form a crucial regulatory level and have the ability to alter an mRNA's stability, localization, and protein products. Alternative splicing events are controlled by a variety of cis-elements, including the presence of consensus splice sequences at the intron-exon border and intronic and exonic splicing enhancer sequences (Pertea et al., 2007). Trans-acting factors also exert a strong effect on splicing and are mostly composed of Ser/Arg-rich proteins, which typically promote intron removal, and heterogenous nuclear ribonucleoproteins, which typically inhibit it (Erkelenz et al., 2013). A host of other indirect factors can affect splicing, including transcription rate, methylation status, and any cellular conditions that alter RNA secondary structure (Kornblihtt et al

SummarySorghum is the fifth most widely planted cereal crop in the world and is commonly cultivated in arid and semi‐arid regions such as Africa. Despite its importance as a food source, sorghum genetic improvement through transgenic approaches has been limited because of an inefficient transformation system. Here, we report a ternary vector (also known as cohabitating vector) system using a recently described pVIR accessory plasmid that facilitates efficient Agrobacterium‐mediated transformation of sorghum. We report regeneration frequencies ranging from 6% to 29% in Tx430 using different selectable markers and single copy, backbone free ‘quality events’ ranging from 45% to 66% of the total events produced. Furthermore, we successfully applied this ternary system to develop transformation protocols for popular but recalcitrant African varieties including Macia, Malisor 84‐7 and Tegemeo. In addition, we report the use of this technology to develop the first stable CRISPR/Cas9‐mediated gene knockouts in Tx430.

Background Gene expression is a key determinant of cellular response. Natural variation in gene expression bridges genetic variation to phenotypic alteration. Identification of the regulatory variants controlling the gene expression in response to drought, a major environmental threat of crop production worldwide, is of great value for drought-tolerant gene identification. Results A total of 627 RNA-seq analyses are performed for 224 maize accessions which represent a wide genetic diversity under three water regimes; 73,573 eQTLs are detected for about 30,000 expressing genes with high-density genome-wide single nucleotide polymorphisms, reflecting a comprehensive and dynamic genetic architecture of gene expression in response to drought. The regulatory variants controlling the gene expression constitutively or drought-dynamically are unraveled. Focusing on dynamic regulatory variants resolved to genes encoding transcription factors, a drought-responsive network reflecting a hierarchy of transcription factors and their target genes is built. Moreover, 97 genes are prioritized to associate with drought tolerance due to their expression variations through the Mendelian randomization analysis. One of the candidate genes, Abscisic acid 8′-hydroxylase, is verified to play a negative role in plant drought tolerance. Conclusions This study unravels the effects of genetic variants on gene expression dynamics in drought response which allows us to better understand the role of distal and proximal genetic effects on gene expression and phenotypic plasticity. The prioritized drought-associated genes may serve as direct targets for functional investigation or allelic mining.

Molecular characterization of events is an integral part of the advancement process during genetically modified (GM) crop product development. Assessment of these events is traditionally accomplished by polymerase chain reaction (PCR) and Southern blot analyses. Southern blot analysis can be time-consuming and comparatively expensive and does not provide sequence-level detail. We have developed a sequence-based application, Southernby-Sequencing (SbS), utilizing sequence capture coupled with nextgeneration sequencing (NGS) technology to replace Southern blot analysis for event selection in a high-throughput molecular characterization environment. SbS is accomplished by hybridizing indexed and pooled whole-genome DNA libraries from GM plants to biotinylated probes designed to target the sequence of transformation plasmids used to generate events within the pool. This sequence capture process enriches the sequence data obtained for targeted regions of interest (transformation plasmid DNA). Taking advantage of the DNA adjacent to the targeted bases (referred to as nextto-target sequence) that accompanies the targeted transformation plasmid sequence, the data analysis detects plasmid-to-genome and plasmid-to-plasmid junctions introduced during insertion into the plant genome. Analysis of these junction sequences provides sequence-level information as to the following: the number of insertion loci including detection of unlinked, independently segregating, small DNA fragments; copy number; rearrangements, truncations, or deletions of the intended insertion DNA; and the presence of transformation plasmid backbone sequences. This molecular evidence from SbS analysis is used to characterize and select GM plants meeting optimal molecular characterization criteria. SbS technology has proven to be a robust event screening tool for use in a highthroughput molecular characterization environment.

CRISPR-Cas9 enabled genome engineering has great potential for improving agriculture productivity, but the possibility of unintended off-target edits has evoked some concerns. Here we employ a three-step strategy to investigate Cas9 nuclease specificity in a complex plant genome. Our approach pairs computational prediction with genome-wide biochemical off-target detection followed by validation in maize plants. Our results reveal high frequency (up to 90%) on-target editing with no evidence of off-target cleavage activity when guide RNAs were bioinformatically predicted to be specific. Predictable off-target edits were observed but only with a promiscuous guide RNA intentionally designed to validate our approach. Off-target editing can be minimized by designing guide RNAs that are different from other genomic locations by at least three mismatches in combination with at least one mismatch occurring in the PAM proximal region. With well-designed guides, genetic variation from Cas9 off-target cleavage in plants is negligible, and much less than inherent variation.

Class 2 CRISPR systems are exceptionally diverse, nevertheless, all share a single effector protein that contains a conserved RuvC-like nuclease domain. Interestingly, the size of these CRISPR-associated (Cas) nucleases ranges from >1000 amino acids (aa) for Cas9/Cas12a to as small as 400-600 aa for Cas12f. For in vivo genome editing applications, compact RNA-guided nucleases are desirable and would streamline cellular delivery approaches. Although miniature Cas12f effectors have been shown to cleave double-stranded DNA, targeted DNA modification in eukaryotic cells has yet to be demonstrated. Here, we biochemically characterize two miniature type V-F Cas nucleases, SpCas12f1 (497 aa) and AsCas12f1 (422 aa), and show that SpCas12f1 functions in both plant and human cells to produce targeted modifications with outcomes in plants being enhanced with short heat pulses. Our findings pave the way for the development of miniature Cas12f1-based genome editing tools.

tassel-less1 (tls1) is a classical maize (Zea mays) inflorescence mutant. Homozygous mutant plants have no tassels or very small tassels, and ear development is also impaired. Using a positional cloning approach, ZmNIP3;1 (a NOD26-like intrinsic protein) was identified as the candidate gene for tls1. The ZmNIP3;1 gene is completely deleted in the tls1 mutant genome. Two Mutator-insertional TUSC alleles of ZmNIP3;1 exhibited tls1-like phenotypes, and allelism tests confirmed that the tls1 gene encodes ZmNIP3;1. Transgenic plants with an RNA interference (RNAi) construct to down-regulate ZmNIP3;1 also showed tls1-like phenotypes, further demonstrating that TLS1 is ZmNIP3;1. Sequence analysis suggests that ZmNIP3;1 is a boron channel protein. Foliar application of boron could rescue the tls1 phenotypes and restore the normal tassel and ear development. Gene expression analysis indicated that in comparison with that of the wild type or tls1 plants treated with boron, the transition from the vegetative to reproductive phase or the development of the floral meristem is impaired in the shoot apical meristem of the tls1 mutant plants. It is concluded that the tls1 mutant phenotypes are caused by impaired boron transport, and boron is essential for inflorescence development in maize.