A Comprehensive Analysis of Alternative Splicing in Paleopolyploid Maize

Citation: Mei, W. B., Liu, S. Z., Schnable, J. C., Yeh, C. T., Springer, N. M., Schnable, P. S., & Barbazuk, W. B. (2017). A Comprehensive Analysis of Alternative Splicing in Paleopolyploid Maize. Frontiers in Plant Science, 8, 19. https://doi.org/10.3389/fpls.2017.00694Identifying and characterizing alternative splicing (AS) enables our understanding of the biological role of transcript isoform diversity. This study describes the use of publicly available RNA-Seq data to identify and characterize the global diversity of AS isoforms in maize using the inbred lines B73 and Mo17, and a related species, sorghum. Identification and characterization of AS within maize tissues revealed that genes expressed in seed exhibit the largest differential AS relative to other tissues examined. Additionally, differences in AS between the two genotypes B73 and Mo17 are greatest within genes expressed in seed. We demonstrate that changes in the level of alternatively spliced transcripts (intron retention and exon skipping) do not solely reflect differences in total transcript abundance, and we present evidence that intron retention may act to fine-tune gene expression across seed development stages. Furthermore, we have identified temperature sensitive AS in maize and demonstrate that drought-induced changes in AS involve distinct sets of genes in reproductive and vegetative tissues. Examining our identified AS isoforms within B73 X Mo17 recombinant inbred lines (RILs) identified splicing QTL (sQTL). The 43.3% of cis-sQTL regulated junctions are actually identified as alternatively spliced junctions in our analysis, while 10 Mb windows on each side of 48.2% of trans-sQTLs overlap with splicing related genes. Using sorghum as an out-group enabled direct examination of loss or conservation of AS between homeologous genes representing the two subgenomes of maize. We identify several instances where AS isoforms that are conserved between one maize homeolog and its sorghum ortholog are absent from the second maize homeolog, suggesting that these AS isoforms may have been lost after the maize whole genome duplication event. This comprehensive analysis provides new insights into the complexity of AS in maize

Genetic control of maize shoot apical meristem architecture

The shoot apical meristem contains a pool of undifferentiated stem cells and generates all above-ground organs of the plant. During vegetative growth, cells differentiate from the meristem to initiate leaves while the pool of meristematic cells is preserved; this balance is determined in part by genetic regulatory mechanisms. To assess vegetative meristem growth and genetic control in Zea mays, we investigated its morphology at multiple time points and identified three stages of growth. We measured meristem height, width, plastochron internode length, and associated traits from 86 individuals of the intermated B73 · Mo17 recombinant inbred line population. For meristem height-related traits, the parents exhibited markedly different phenotypes, with B73 being very tall, Mo17 short, and the population distributed between. In the outer cell layer, differences appeared to be related to number of cells rather than cell size. In contrast, B73 and Mo17 were similar in meristem width traits and plastochron internode length, with transgressive segregation in the population. Multiple loci (629 for each trait) were mapped, indicating meristem architecture is controlled by many regions; none of these coincided with previously described mutants impacting meristem development. Major loci for height and width explaining 16% and 19% of the variation were identified on chromosomes 5 and 8, respectively. Significant loci for related traits frequently coincided, whereas those for unrelated traits did not overlap. With the use of three near-isogenic lines, a locus explaining 16% of the parental variation in meristem height was validated. Published expression data were leveraged to identify candidate genes in significant regions. © 2014 Thompson et al

Characterization of maize roothairless6 which encodes a D-type cellulose synthase and controls the switch from bulge formation to tip growth

Citation: Li, L., Hey, S., Liu, S. Z., Liu, Q., McNinch, C., Hu, H. C., . . . Hochholdinger, F. (2016). Characterization of maize roothairless6 which encodes a D-type cellulose synthase and controls the switch from bulge formation to tip growth. Scientific Reports, 6, 12. doi:10.1038/srep34395Root hairs are tubular extensions of the epidermis. Root hairs of the monogenic recessive maize mutant roothairless 6 (rth6) are arrested after bulge formation during the transition to tip growth and display a rough cell surface. BSR-Seq in combination with Seq-walking and subsequent analyses of four independently generated mutant alleles established that rth6 encodes CSLD5 a plasma membrane localized 129 kD D-type cellulose synthase with eight transmembrane domains. Cellulose synthases are required for the biosynthesis of cellulose, the most abundant biopolymer of plant cell walls. Phylogenetic analyses revealed that RTH6 is part of a monocot specific clade of D-type cellulose synthases. D-type cellulose synthases are highly conserved in the plant kingdom with five gene family members in maize and homologs even among early land plants such as the moss Physcomitrella patens or the clubmoss Selaginella moellendorffii. Expression profiling demonstrated that rth6 transcripts are highly enriched in root hairs as compared to all other root tissues. Moreover, in addition to the strong knock down of rth6 expression in young primary roots of the mutant rth6, the gene is also significantly down-regulated in rth3 and rth5 mutants, while it is up-regulated in rth2 mutants, suggesting that these genes interact in cell wall biosynthesis

Punctate vascular expression1 is a novel maize gene required for leaf pattern formation that functions downstream of the trans-acting small interfering RNA pathway

The maize (Zea mays) gene RAGGED SEEDLING2-R (RGD2-R) encodes an ARGONAUTE7-like protein required for the biogenesis of trans-acting small interfering RNA, which regulates the accumulation of AUXIN RESPONSE FACTOR3A transcripts in shoots. Although dorsiventral polarity is established in the narrow and cylindrical leaves of rgd2-R mutant plants, swapping of adaxial/abaxial epidermal identity occurs and suggests a model wherein RGD2 is required to coordinate dorsiventral and mediolateral patterning in maize leaves. Laser microdissection-microarray analyses of the rgd2-R mutant shoot apical meristem identified a novel gene, PUNCTATE VASCULAR EXPRESSION1 (PVE1), that is down-regulated in rgd2-R mutant apices. Transcripts of PVE1 provide an early molecular marker for vascular morphogenesis. Reverse genetic analyses suggest that PVE1 functions during vascular development and in mediolateral and dorsiventral patterning of maize leaves. Molecular genetic analyses of PVE1 and of rgd2-R;pve1-M2 double mutants suggest a model wherein PVE1 functions downstream of RGD2 in a pathway that intersects and interacts with the trans-acting small interfering RNA pathway

Last Frontier of Agricultural Big Data in Rotation?

Over the years new and promising varieties are bred not only for resistance to diseases but also for direct yield by stress tolerance, architectural or other properties. However, these potential yield improvements under optimal circumstances are not met by the growers in their practice, and is even widening. ..

Mendelian and Non-Mendelian Regulation of Gene Expression in Maize

Transcriptome variation plays an important role in affecting the phenotype of an organism. However, an understanding of the underlying mechanisms regulating transcriptome variation in segregating populations is still largely unknown. We sought to assess and map variation in transcript abundance in maize shoot apices in the intermated B73×Mo17 recombinant inbred line population. RNA-based sequencing (RNA-seq) allowed for the detection and quantification of the transcript abundance derived from 28,603 genes. For a majority of these genes, the population mean, coefficient of variation, and segregation patterns could be predicted by the parental expression levels. Expression quantitative trait loci (eQTL) mapping identified 30,774 eQTL including 96 trans-eQTL "hotspots," each of which regulates the expression of a large number of genes. Interestingly, genes regulated by a trans-eQTL hotspot tend to be enriched for a specific function or act in the same genetic pathway. Also, genomic structural variation appeared to contribute to cis-regulation of gene expression. Besides genes showing Mendelian inheritance in the RIL population, we also found genes whose expression level and variation in the progeny could not be predicted based on parental difference, indicating that non-Mendelian factors also contribute to expression variation. Specifically, we found 145 genes that show patterns of expression reminiscent of paramutation such that all the progeny had expression levels similar to one of the two parents. Furthermore, we identified another 210 genes that exhibited unexpected patterns of transcript presence/absence. Many of these genes are likely to be gene fragments resulting from transposition, and the presence/absence of their transcripts could influence expression levels of their ancestral syntenic genes. Overall, our results contribute to the identification of novel expression patterns and broaden the understanding of transcriptional variation in plants. © 2013 Lin et al

Genic and nongenic contributions to natural variation of quantitative traits in maize

The complex genomes of many economically important crops present tremendous challenges to understand the genetic control of many quantitative traits with great importance in crop production, adaptation, and evolution. Advances in genomic technology need to be integrated with strategic genetic design and novel perspectives to break new ground. Complementary to individual-gene-targeted research, which remains challenging, a global assessment of the genomic distribution of trait-associated SNPs (TASs) discovered from genome scans of quantitative traits can provide insights into the genetic architecture and contribute to the design of future studies. Here we report the first systematic tabulation of the relative contribution of different genomic regions to quantitative trait variation in maize. We found that TASs were enriched in the nongenic regions, particularly within a 5-kb window upstream of genes, which highlights the importance of polymorphisms regulating gene expression in shaping the natural variation. Consistent with these findings, TASs collectively explained 44%-59% of the total phenotypic variation across maize quantitative traits, and on average, 79% of the explained variation could be attributed to TASs located in genes or within 5 kb upstream of genes, which together comprise only 13% of the genome. Our findings suggest that efficient, cost-effective genome-wide association studies (GWAS) in species with complex genomes can focus on genic and promoter regions

Genome-wide discovery and characterization of maize long non-coding RNAs

BACKGROUND: Long non-coding RNAs (lncRNAs) are transcripts that are 200 bp or longer, do not encode proteins, and potentially play important roles in eukaryotic gene regulation. However, the number, characteristics and expression inheritance pattern of lncRNAs in maize are still largely unknown. RESULTS: By exploiting available public EST databases, maize whole genome sequence annotation and RNA-seq datasets from 30 different experiments, we identified 20,163 putative lncRNAs. Of these lncRNAs, more than 90% are predicted to be the precursors of small RNAs, while 1,704 are considered to be high-confidence lncRNAs. High confidence lncRNAs have an average transcript length of 463 bp and genes encoding them contain fewer exons than annotated genes. By analyzing the expression pattern of these lncRNAs in 13 distinct tissues and 105 maize recombinant inbred lines, we show that more than 50% of the high confidence lncRNAs are expressed in a tissue-specific manner, a result that is supported by epigenetic marks. Intriguingly, the inheritance of lncRNA expression patterns in 105 recombinant inbred lines reveals apparent transgressive segregation, and maize lncRNAs are less affected by cis- than by trans- genetic factors. CONCLUSIONS: We integrate all available transcriptomic datasets to identify a comprehensive set of maize lncRNAs, provide a unique annotation resource of the maize genome and a genome-wide characterization of maize lncRNAs, and explore the genetic control of their expression using expression quantitative trait locus mapping

Integrating transposable elements in the 3D genome

Chromosome organisation is increasingly recognised as an essential component of genome regulation, cell fate and cell health. Within the realm of transposable elements (TEs) however, the spatial information of how genomes are folded is still only rarely integrated in experimental studies or accounted for in modelling. Whilst polymer physics is recognised as an important tool to understand the mechanisms of genome folding, in this commentary we discuss its potential applicability to aspects of TE biology. Based on recent works on the relationship between genome organisation and TE integration, we argue that existing polymer models may be extended to create a predictive framework for the study of TE integration patterns. We suggest that these models may offer orthogonal and generic insights into the integration profiles (or "topography") of TEs across organisms. In addition, we provide simple polymer physics arguments and preliminary molecular dynamics simulations of TEs inserting into heterogeneously flexible polymers. By considering this simple model, we show how polymer folding and local flexibility may generically affect TE integration patterns. The preliminary discussion reported in this commentary is aimed to lay the foundations for a large-scale analysis of TE integration dynamics and topography as a function of the three-dimensional host genome

Evolutionary relationships among barley and <i>Arabidopsis</i> core circadian clock and clock-associated genes

The circadian clock regulates a multitude of plant developmental and metabolic processes. In crop species, it contributes significantly to plant performance and productivity and to the adaptation and geographical range over which crops can be grown. To understand the clock in barley and how it relates to the components in the Arabidopsis thaliana clock, we have performed a systematic analysis of core circadian clock and clock-associated genes in barley, Arabidopsis and another eight species including tomato, potato, a range of monocotyledonous species and the moss, Physcomitrella patens. We have identified orthologues and paralogues of Arabidopsis genes which are conserved in all species, monocot/dicot differences, species-specific differences and variation in gene copy number (e.g. gene duplications among the various species). We propose that the common ancestor of barley and Arabidopsis had two-thirds of the key clock components identified in Arabidopsis prior to the separation of the monocot/dicot groups. After this separation, multiple independent gene duplication events took place in both monocot and dicot ancestors. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s00239-015-9665-0) contains supplementary material, which is available to authorized users