Genome-Wide Association Analyses Reveal the Importance of Alternative Splicing in Diversifying Gene Function and Regulating Phenotypic Variation in Maize

Harnessing neo-domestication of wild pigmented rice for enhanced nutrition and sustainable agriculture

Advances in precision gene editing have enabled the rapid domestication of wild crop relatives, a process known as neo-domestication. During domestication, breeding rice for maximum productivity under optimal growth conditions reduced genetic diversity, eliminating variants for stress tolerance and grain nutrients. Wild rice varieties have rich genetic diversity, including variants for disease resistance, stress tolerance, and grain nutritional quality. For example, the grain of pigmented wild rice has abundant antioxidants (anthocyanins, proanthocyanidins, and flavonoids), but low yield, poor plant architecture, and long life cycle limit its cultivation. In this review, we address the neo-domestication of wild pigmented rice, focusing on recent progress, CRISPR-Cas editing toolboxes, selection of key candidate genes for domestication, identifying species with superior potential via generating genomic and multi-omics resources, efficient crop transformation methods and highlight strategies for the promotion and application pigmented rice. We also address critical outstanding questions and potential solutions to enable efficient neo-domestication of wild pigmented rice and thus enhance food security and nutrition.

Characterization of a longan pericarp browning related peroxidase with a focus on its role in proanthocyanidin and lignin polymerization

The longan pericarp turns brown dramatically after harvesting, but the mechanism is not well understood. In this work, two peroxidases were purified from longan pericarp and found to be identical to the class III peroxidases PRX53-2 and PRX53-3. In vitro, PRX53-2/3 catalyzed the browning of several pericarp abundant proanthocyanidin and lignin monomers, such as (-)-epicatechin (EC), (+)-catechin (CT) and coniferyl alcohol (ConA). PRX53-2 was upregulated and highly-expressed, while PRX53-3 was expressed at low levels after harvesting; thus, PRX53-2 was considered a browning-related gene. The reaction with both proanthocyanidin and lignin presented a greater degree of brown coloration compared to the single substrate reactions. Several procyanidins isomers, EC-ConA and CT-ConA were detected in the double-substrate reaction. These results not only demonstrate that the effects of PRX53-2 on proanthocyanidin and lignin polymerization may be crucial for longan pericarp browning, but also help in developing new strategies or preservatives to delay pericarp browning.

From genes to traits: Trends in RNA-binding proteins and their role in plant trait development: A review

RNA-binding proteins (RBPs) are essential for cellular functions by attaching to RNAs, creating dynamic ribonucleoprotein complexes (RNPs) essential for managing RNA throughout its life cycle. These proteins are critical to all post-transcriptional processes, impacting vital cellular functions during development and adaptation to environmental changes. Notably, in plants, RBPs are critical for adjusting to inconsistent environmental conditions, with recent studies revealing that plants possess, more prominent, and both novel and conserved RBP families compared to other eukaryotes. This comprehensive review delves into the varied RBPs covering their structural attributes, domain base function, and their interactions with RNA in metabolism, spotlighting their role in regulating post-transcription and splicing and their reaction to internal and external stimuli. It highlights the complex regulatory roles of RBPs, focusing on plant trait regulation and the unique functions they facilitate, establishing a foundation for appreciating RBPs' significance in plant growth and environmental response strategies.

TWAS facilitates gene-scale trait genetic dissection through gene expression, structural variations, and alternative splicing in soybean

A genome-wide association study (GWAS) identifies trait-associated loci, but identifying the causal genes can be a bottleneck, due in part to slow decay of linkage disequilibrium (LD). A transcriptome-wide association study (TWAS) addresses this issue by identifying gene expression-phenotype associations or integrating gene expression quantitative trait loci with GWAS results. Here, we used self-pollinated soybean (Glycine max [L.] Merr.) as a model to evaluate the application of TWAS to the genetic dissection of traits in plant species with slow LD decay. We generated RNA sequencing data for a soybean diversity panel and identified the genetic expression regulation of 29 286 soybean genes. Different TWAS solutions were less affected by LD and were robust to the source of expression, identifing known genes related to traits from different tissues and developmental stages. The novel pod-color gene L2 was identified via TWAS and functionally validated by genome editing. By introducing a new exon proportion feature, we significantly improved the detection of expression variations that resulted from structural variations and alternative splicing. As a result, the genes identified through our TWAS approach exhibited a diverse range of causal variations, including SNPs, insertions or deletions, gene fusion, copy number variations, and alternative splicing. Using this approach, we identified genes associated with flowering time, including both previously known genes and novel genes that had not previously been linked to this trait, providing insights complementary to those from GWAS. In summary, this study supports the application of TWAS for candidate gene identification in species with low rates of LD decay.

Population-level exploration of alternative splicing and its unique role in controlling agronomic traits of rice

Alternative splicing (AS) plays crucial roles in regulating various biological processes in plants. However, the genetic mechanisms underlying AS and its role in controlling important agronomic traits in rice (Oryza sativa) remain poorly understood. In this study, we explored AS in rice leaves and panicles using the rice minicore collection. Our analysis revealed a high level of transcript isoform diversity, with approximately one-fifth of the potential isoforms acting as major transcripts in both tissues. Regarding the genetic mechanism of AS, we found that the splicing of 833 genes in the leaf and 1,230 genes in the panicle was affected by cis-genetic variation. Twenty-one percent of these AS events could only be explained by large structural variations. Approximately 77.5% of genes with significant splicing quantitative trait loci (sGenes) exhibited tissue-specific regulation, and AS can cause 26.9% (leaf) and 23.6% (panicle) of sGenes to have altered, lost, or gained functional domains. Additionally, through splicing-phenotype association analysis, we identified phosphate-starvation-induced RING-type E3 ligase (OsPIE1; LOC_Os01g72480), whose splicing ratio was significantly associated with plant height. In summary, this study provides an understanding of AS in rice and its contribution to the regulation of important agronomic traits.

RNA-Binding Protein-Mediated Alternative Splicing Regulates Abiotic Stress Responses in Plants

The alternative splicing of pre-mRNA generates distinct mRNA variants from a pre-mRNA, thereby modulating a gene's function. The splicing of pre-mRNA depends on splice sites and regulatory elements in pre-mRNA, as well as the snRNA and proteins that recognize these sequences. Among these, RNA-binding proteins (RBPs) are the primary regulators of pre-mRNA splicing and play a critical role in the regulation of alternative splicing by recognizing the elements in pre-mRNA. However, little is known about the function of RBPs in stress response in plants. Here, we summarized the RBPs involved in the alternative splicing of pre-mRNA and their recognizing elements in pre-mRNA, and the recent advance in the role of RBP-mediated alternative splicing in response to abiotic stresses in plants. This review proposes that the regulation of pre-mRNA alternative splicing by RBPs is an important way for plants to adapt to abiotic stresses, and the regulation of alternative splicing by RBPs is a promising direction for crop breeding.

When splicing meets transgene DNA methylation to uncover elite alleles for plant yield

This article is a Commentary on Zhao et al. (2024), 243: 1724–1741.

Alternative splicing as a driver of natural variation in abscisic acid response

Abscisic acid (ABA) is a crucial player in plant responses to the environment. It accumulates under stress, activating downstream signaling to implement molecular responses that restore homeostasis. Natural variance in ABA sensitivity remains barely understood, and the ABA pathway has been mainly studied at the transcriptional level, despite evidence that posttranscriptional regulation, namely, via alternative splicing, contributes to plant stress tolerance. Here, we identified the Arabidopsis accession Kn-0 as less sensitive to ABA than the reference Col-0, as shown by reduced effects of the hormone on seedling establishment, root branching, and stomatal closure, as well as by decreased induction of ABA marker genes. An in-depth comparative transcriptome analysis of the ABA response in the two variants revealed lower expression changes and fewer genes affected for the least ABA-sensitive ecotype. Notably, Kn-0 exhibited reduced levels of the ABA-signaling SnRK2 protein kinases and lower basal expression of ABA-reactivation genes, consistent with our finding that Kn-0 contains less endogenous ABA than Col-0. ABA also markedly affected alternative splicing, primarily intron retention, with Kn-0 being less responsive regarding both the number and magnitude of alternative splicing events, particularly exon skipping. We find that alternative splicing introduces a more ecotype-specific layer of ABA regulation and identify ABA-responsive splicing changes in key ABA pathway regulators that provide a functional and mechanistic link to the differential sensitivity of the two ecotypes. Our results offer new insight into the natural variation of ABA responses and corroborate a key role for alternative splicing in implementing ABA-mediated stress responses.

Genetic insight into a polygenic trait using a novel genome-wide association approach in a wild amphibian population

Body size variation is central in the evolution of life-history traits in amphibians, but the underlying genetic architecture of this complex trait is still largely unknown. Herein, we studied the genetic basis of body size and fecundity of the alternative morphotypes in a wild population of the Greek smooth newt (Lissotriton graecus). By combining a genome-wide association approach with linkage disequilibrium network analysis, we were able to identify clusters of highly correlated loci thus maximizing sequence data for downstream analysis. The putatively associated variants explained 12.8% to 44.5% of the total phenotypic variation in body size and were mapped to genes with functional roles in the regulation of gene expression and cell cycle processes. Our study is the first to provide insights into the genetic basis of complex traits in newts and provides a useful tool to identify loci potentially involved in fitness-related traits in small data sets from natural populations in non-model species.

Divergence in regulatory mechanisms of GR-RBP genes in different plants under abiotic stress

The IVa subfamily of glycine-rich proteins (GRPs) comprises a group of glycine-rich RNA binding proteins referred to as GR-RBPa here. Previous studies have demonstrated functions of GR-RBPa proteins in regulating stress response in plants. However, the mechanisms responsible for the differential regulatory functions of GR-RBPa proteins in different plant species have not been fully elucidated. In this study, we identified and comprehensively studied a total of 34 GR-RBPa proteins from five plant species. Our analysis revealed that GR-RBPa proteins were further classified into two branches, with proteins in branch I being relatively more conserved than those in branch II. When subjected to identical stresses, these genes exhibited intensive and differential expression regulation in different plant species, corresponding to the enrichment of cis-acting regulatory elements involving in environmental and internal signaling in these genes. Unexpectedly, all GR-RBPa genes in branch I underwent intensive alternative splicing (AS) regulation, while almost all genes in branch II were only constitutively spliced, despite having more introns. This study highlights the complex and divergent regulations of a group of conserved RNA binding proteins in different plants when exposed to identical stress conditions. These species-specific regulations may have implications for stress responses and adaptations in different plant species.

Comprehensive deciphering the alternative splicing patterns involved in leaf morphogenesis of Liriodendron chinense

Alternative splicing (AS), a pivotal post-transcriptional regulatory mechanism, profoundly amplifies diversity and complexity of transcriptome and proteome. Liriodendron chinense (Hemsl.) Sarg., an excellent ornamental tree species renowned for its distinctive leaf shape, which resembles the mandarin jacket. Despite the documented potential genes related to leaf development of L. chinense, the underlying post-transcriptional regulatory mechanisms remain veiled. Here, we conducted a comprehensive analysis of the transcriptome to clarify the genome-wide landscape of the AS pattern and the spectrum of spliced isoforms during leaf developmental stages in L. chinense. Our investigation unveiled 50,259 AS events, involving 10,685 genes (32.9%), with intron retention as the most prevalent events. Notably, the initial stage of leaf development witnessed the detection of 804 differentially AS events affiliated with 548 genes. Although both differentially alternative splicing genes (DASGs) and differentially expressed genes (DEGs) were enriched into morphogenetic related pathways during the transition from fishhook (P2) to lobed (P7) leaves, there was only a modest degree of overlap between DASGs and DEGs. Furthermore, we conducted a comprehensively AS analysis on homologous genes involved in leaf morphogenesis, and most of which are subject to post-transcriptional regulation of AS. Among them, the AINTEGUMENTA-LIKE transcript factor LcAIL5 was characterization in detailed, which experiences skipping exon (SE), and two transcripts displayed disparate expression patterns across multiple stages. Overall, these findings yield a comprehensive understanding of leaf development regulation via AS, offering a novel perspective for further deciphering the mechanism of plant leaf morphogenesis.

Synergizing biotechnology and natural farming: pioneering agricultural sustainability through innovative interventions

The world has undergone a remarkable transformation from the era of famines to an age of global food production that caters to an exponentially growing population. This transformation has been made possible by significant agricultural revolutions, marked by the intensification of agriculture through the infusion of mechanical, industrial, and economic inputs. However, this rapid advancement in agriculture has also brought about the proliferation of agricultural inputs such as pesticides, fertilizers, and irrigation, which have given rise to long-term environmental crises. Over the past two decades, we have witnessed a concerning plateau in crop production, the loss of arable land, and dramatic shifts in climatic conditions. These challenges have underscored the urgent need to protect our global commons, particularly the environment, through a participatory approach that involves countries worldwide, regardless of their developmental status. To achieve the goal of sustainability in agriculture, it is imperative to adopt multidisciplinary approaches that integrate fields such as biology, engineering, chemistry, economics, and community development. One noteworthy initiative in this regard is Zero Budget Natural Farming, which highlights the significance of leveraging the synergistic effects of both plant and animal products to enhance crop establishment, build soil fertility, and promote the proliferation of beneficial microorganisms. The ultimate aim is to create self-sustainable agro-ecosystems. This review advocates for the incorporation of biotechnological tools in natural farming to expedite the dynamism of such systems in an eco-friendly manner. By harnessing the power of biotechnology, we can increase the productivity of agro-ecology and generate abundant supplies of food, feed, fiber, and nutraceuticals to meet the needs of our ever-expanding global population.

The impact of genetically controlled splicing on exon inclusion and protein structure

Common variants affecting mRNA splicing are typically identified though splicing quantitative trait locus (sQTL) mapping and have been shown to be enriched for GWAS signals by a similar degree to eQTLs. However, the specific splicing changes induced by these variants have been difficult to characterize, making it more complicated to analyze the effect size and direction of sQTLs, and to determine downstream splicing effects on protein structure. In this study, we catalogue sQTLs using exon percent spliced in (PSI) scores as a quantitative phenotype. PSI is an interpretable metric for identifying exon skipping events and has some advantages over other methods for quantifying splicing from short read RNA sequencing. In our set of sQTL variants, we find evidence of selective effects based on splicing effect size and effect direction, as well as exon symmetry. Additionally, we utilize AlphaFold2 to predict changes in protein structure associated with sQTLs overlapping GWAS traits, highlighting a potential new use-case for this technology for interpreting genetic effects on traits and disorders.

Natural uORF variation in plants

Taking advantage of natural variation promotes our understanding of phenotypic diversity and trait evolution, ultimately accelerating plant breeding, in which the identification of causal variations is critical. To date, sequence variations in the coding region and transcription level polymorphisms caused by variations in the promoter have been prioritized. An upstream open reading frame (uORF) in the 5' untranslated region (5' UTR) regulates gene expression at the post-transcription or translation level. In recent years, studies have demonstrated that natural uORF variations shape phenotypic diversity. This opinion article highlights recent researches and speculates on future directions for natural uORF variation in plants.

iFLAS: positive-unlabeled learning facilitates full-length transcriptome-based identification and functional exploration of alternatively spliced isoforms in maize

The advent of full-length transcriptome sequencing technologies has accelerated the discovery of novel splicing isoforms. However, existing alternative splicing (AS) tools are either tailored for short-read RNA-Seq data or designed for human and animal studies. The disparities in AS patterns between plants and animals still pose a challenge to the reliable identification and functional exploration of novel isoforms in plants. Here, we developed integrated full-length alternative splicing analysis (iFLAS), a plant-optimized AS toolkit that introduced a semi-supervised machine learning method known as positive-unlabeled (PU) learning to accurately identify novel isoforms. iFLAS also enables the investigation of AS functions from various perspectives, such as differential AS, poly(A) tail length, and allele-specific AS (ASAS) analyses. By applying iFLAS to three full-length transcriptome sequencing datasets, we systematically identified and functionally characterized maize (Zea mays) AS patterns. We found intron retention not only introduces premature termination codons, resulting in lower expression levels of isoforms, but may also regulate the length of 3'UTR and poly(A) tail, thereby affecting the functional differentiation of isoforms. Moreover, we observed distinct ASAS patterns in two genes within heterosis offspring, highlighting their potential value in breeding. These results underscore the broad applicability of iFLAS in plant full-length transcriptome-based AS research.

Genome-Wide Association Study for Maize Hybrid Performance in a Typical Breeder Population

Maize is one of the major crops that has demonstrated success in the utilization of heterosis. Developing high-yield hybrids is a crucial part of plant breeding to secure global food demand. In this study, we conducted a genome-wide association study (GWAS) for 10 agronomic traits using a typical breeder population comprised 442 single-cross hybrids by evaluating additive, dominance, and epistatic effects. A total of 49 significant single nucleotide polymorphisms (SNPs) and 69 significant pairs of epistasis were identified, explaining 26.2% to 64.3% of the phenotypic variation across the 10 traits. The enrichment of favorable genotypes is significantly correlated to the corresponding phenotype. In the confident region of the associated site, 532 protein-coding genes were discovered. Among these genes, the Zm00001d044211 candidate gene was found to negatively regulate starch synthesis and potentially impact yield. This typical breeding population provided a valuable resource for dissecting the genetic architecture of yield-related traits. We proposed a novel mating strategy to increase the GWAS efficiency without utilizing more resources. Finally, we analyzed the enrichment of favorable alleles in the Shaan A and Shaan B groups, as well as in each inbred line. Our breeding practice led to consistent results. Not only does this study demonstrate the feasibility of GWAS in F1 hybrid populations, it also provides a valuable basis for further molecular biology and breeding research.

Alternative Splicing Variation: Accessing and Exploiting in Crop Improvement Programs

Alternative splicing (AS) is a gene regulatory mechanism modulating gene expression in multiple ways. AS is prevalent in all eukaryotes including plants. AS generates two or more mRNAs from the precursor mRNA (pre-mRNA) to regulate transcriptome complexity and proteome diversity. Advances in next-generation sequencing, omics technology, bioinformatics tools, and computational methods provide new opportunities to quantify and visualize AS-based quantitative trait variation associated with plant growth, development, reproduction, and stress tolerance. Domestication, polyploidization, and environmental perturbation may evolve novel splicing variants associated with agronomically beneficial traits. To date, pre-mRNAs from many genes are spliced into multiple transcripts that cause phenotypic variation for complex traits, both in model plant Arabidopsis and field crops. Cataloguing and exploiting such variation may provide new paths to enhance climate resilience, resource-use efficiency, productivity, and nutritional quality of staple food crops. This review provides insights into AS variation alongside a gene expression analysis to select for novel phenotypic diversity for use in breeding programs. AS contributes to heterosis, enhances plant symbiosis (mycorrhiza and rhizobium), and provides a mechanistic link between the core clock genes and diverse environmental clues.

Mapping intron retention events contributing to complex traits using splice quantitative trait locus

BACKGROUND

Alternative splicing (AS) of mRNA plays an important roles in transcriptome diversity, involving regulation of plant growth and stress response. Understanding the variation of AS events underlying GWAS loci in a crop population can provide insight into the molecular mechanisms of complex agronomic traits. To date, genome-wide association studies relating AS events to agronomic traits have rarely been conducted at the population level in crops.

RESULTS

Here, a pipeline was constructed to identify candidate AS events related to complex traits. Firstly, ovule transcriptome data were used to characterize intron retention (IR), the predominant type of AS in plants, on a genome-wide scale. This was done in a natural population consisting of 279 upland cotton lines. Secondly, splice quantitative trait locus (sQTL) analysis was carried out, which yielded a total of 2295 sQTLs involving 1607 genes. Of these, 14.25% (n = 427) were cis-regulatory loci. Integration with expression quantitative trait loci (eQTL) revealed that 53 (21.4%) cis-sGenes were regulated by both cis-sQTLs and cis-eQTLs. Finally, co-localization analysis integrated with GWAS loci in this population showed 32 cis-QTLs to be co-located with genetic regulatory loci related to fiber yield and quality traits, indicating that sQTLs are likely to participate in regulating cotton fiber yield and quality. An in-depth evaluation confirmed that differences in the IR rates of sQTL-regulated candidate genes such as GhLRRK1 and GhGC1 are associated with lint percentage (LP), which has potential in breeding applications.

CONCLUSION

This study provides a clue that AS of mRNA has an impact on crop yield, along with functional sQTLs are new genetic resources for cotton precision breeding.

Design, execution, and interpretation of plant RNA-seq analyses

Genomics has transformed our understanding of the genetic architecture of traits and the genetic variation present in plants. Here, we present a review of how RNA-seq can be performed to tackle research challenges addressed by plant sciences. We discuss the importance of experimental design in RNA-seq, including considerations for sampling and replication, to avoid pitfalls and wasted resources. Approaches for processing RNA-seq data include quality control and counting features, and we describe common approaches and variations. Though differential gene expression analysis is the most common analysis of RNA-seq data, we review multiple methods for assessing gene expression, including detecting allele-specific gene expression and building co-expression networks. With the production of more RNA-seq data, strategies for integrating these data into genetic mapping pipelines is of increased interest. Finally, special considerations for RNA-seq analysis and interpretation in plants are needed, due to the high genome complexity common across plants. By incorporating informed decisions throughout an RNA-seq experiment, we can increase the knowledge gained.

Plant serine/arginine-rich proteins: versatile players in RNA processing

MAIN CONCLUSION

Serine/arginine-rich (SR) proteins participate in RNA processing by interacting with precursor mRNAs or other splicing factors to maintain plant growth and stress responses. Alternative splicing is an important mechanism involved in mRNA processing and regulation of gene expression at the posttranscriptional level, which is the main reason for the diversity of genes and proteins. The process of alternative splicing requires the participation of many specific splicing factors. The SR protein family is a splicing factor in eukaryotes. The vast majority of SR proteins' existence is an essential survival factor. Through its RS domain and other unique domains, SR proteins can interact with specific sequences of precursor mRNA or other splicing factors and cooperate to complete the correct selection of splicing sites or promote the formation of spliceosomes. They play essential roles in the composition and alternative splicing of precursor mRNAs, providing pivotal functions to maintain growth and stress responses in animals and plants. Although SR proteins have been identified in plants for three decades, their evolutionary trajectory, molecular function, and regulatory network remain largely unknown compared to their animal counterparts. This article reviews the current understanding of this gene family in eukaryotes and proposes potential key research priorities for future functional studies.

A role for heritable transcriptomic variation in maize adaptation to temperate environments

Background

Transcription bridges genetic information and phenotypes. Here, we evaluated how changes in transcriptional regulation enable maize (Zea mays), a crop originally domesticated in the tropics, to adapt to temperate environments.

Result

We generated 572 unique RNA-seq datasets from the roots of 340 maize genotypes. Genes involved in core processes such as cell division, chromosome organization and cytoskeleton organization showed lower heritability of gene expression, while genes involved in anti-oxidation activity exhibited higher expression heritability. An expression genome-wide association study (eGWAS) identified 19,602 expression quantitative trait loci (eQTLs) associated with the expression of 11,444 genes. A GWAS for alternative splicing identified 49,897 splicing QTLs (sQTLs) for 7614 genes. Genes harboring both cis-eQTLs and cis-sQTLs in linkage disequilibrium were disproportionately likely to encode transcription factors or were annotated as responding to one or more stresses. Independent component analysis of gene expression data identified loci regulating co-expression modules involved in oxidation reduction, response to water deprivation, plastid biogenesis, protein biogenesis, and plant-pathogen interaction. Several genes involved in cell proliferation, flower development, DNA replication, and gene silencing showed lower gene expression variation explained by genetic factors between temperate and tropical maize lines. A GWAS of 27 previously published phenotypes identified several candidate genes overlapping with genomic intervals showing signatures of selection during adaptation to temperate environments.

Conclusion

Our results illustrate how maize transcriptional regulatory networks enable changes in transcriptional regulation to adapt to temperate regions.

Supplementary information

The online version contains supplementary material available at 10.1186/s13059-023-02891-3.

De Novo Domestication in the Multi-Omics Era

Most cereal crops were domesticated within the last 12,000 years and subsequently spread around the world. These crops have been nourishing the world by supplying a primary energy and nutrient source, thereby playing a critical role in determining the status of human health and sustaining the global population. Here, we review the major challenges of future agriculture and emphasize the utilization of wild germplasm. De novo domestication is one of the most straightforward strategies to manipulate domestication-related and/or other genes with known function, and thereby introduce desired traits into wild plants. We also summarize known causal variations and their corresponding pathways in order to better understand the genetic basis of crop evolution, and how this knowledge could facilitate de novo domestication. Indeed knowledge-driven de novo domestication has great potential for the development of new sustainable crops that have climate-resilient high yield with low resource input and meet individual nutrient needs. Finally, we discuss current opportunities for and barriers to knowledge-driven de novo domestication.

Population-level transcriptomes reveal gene expression and splicing underlying cadmium accumulation in barley

Genetic variation is an important determinant of gene transcription, which in turn contributes to functional and phenotypic diversity. Identification of the genetic variants controlling gene expression and alternative splicing in crops responding to cadmium (Cd), an important issue for food safety and human health, is of great value to improve our understanding of Cd accumulation-related genes. Here we report an in-depth survey of population-level transcriptome variation of barley (Hordeum vulgare) core accessions under Cd exposure. We reveal marked transcriptomic changes in response to Cd exposure, and these are largely independent of tissues. A genome-wide association study (GWAS) revealed 59 498 expression quantitative trait loci (eQTLs) and 23 854 splicing quantitative trait loci (sQTLs), leading to a complex network that covers 66.6% of the expressed genes, including 68 metal transporter genes. On average, 41.0% of sQTLs overlapped with eQTLs across different tissues, indicating that these two dimensions of transcript variation are largely independent. Moreover, we found that 34.5% of GWAS QTLs that underlie 10 Cd accumulation traits in barley are co-localized with eQTLs and sQTLs, which could imply a mechanistic role of different genetic variants affecting gene expression and alternative splicing in these traits. This study highlights the role of distal and proximal genetic effects on gene expression, splicing, and phenotypic plasticity. We anticipate that our results on the genetic control of expression and splicing underlying Cd accumulation provide a bridge to better understand genetic variation and phenotypic diversity to elucidate the mechanisms underlying Cd accumulation in plants.

Differential alternative splicing genes and isoform co-expression networks of Brassica napus under multiple abiotic stresses

Alternative splicing (AS) is an important regulatory process that affects plant development and stress responses by greatly increasing the complexity of transcriptome and proteome. To understand how the AS landscape of B. napus changes in response to abiotic stresses, we investigated 26 RNA-seq libraries, including control and treatments with cold, dehydration, salt, and abscisic acid (ABA) at two different time points, to perform comparative alternative splicing analysis. Apparently, AS events increased under all stresses except dehydration for 1 h, and intron retention was the most common AS mode. In addition, a total of 357 differential alternative splicing (DAS) genes were identified under four abiotic stresses, among which 81 DAS genes existed in at least two stresses, and 276 DAS genes were presented under only one stress. A weighted gene co-expression network analysis (WGCNA) based on the splicing isoforms, rather than the genes, pinpointed out 23 co-expression modules associated with different abiotic stresses. Among them, a number of significant hub genes were also found to be DAS genes, which encode key isoforms involved in responses to single stress or multiple stresses, including RNA-binding proteins, transcription factors, and other important genes, such as RBP45C, LHY, MYB59, SCL30A, RS40, MAJ23.10, and DWF4. The splicing isoforms of candidate genes identified in this study could be a valuable resource for improving tolerance of B. napus against multiple abiotic stresses.

Regulation of alternative splicing of PaFT and PaFDL1, the FT and FD homologs in Platanus acerifolia

Alternative splicing (AS) selects different alternative splice sites and produces a variety of transcripts with different exon/intron combinations, which may result in multiple protein isoforms. The splicing signals include cis-elements and RNA structures; however, the mechanisms of AS regulation in plants have yet to be elucidated. Previous studies have shown that in Platanus acerifolia, the FLOWERING LOCUS T (FT) homolog PaFT has a unique and complex AS pattern, in which most of the splice forms of PaFT involve the first and/or second intron, and the FD homolog PaFDL1 produces two transcripts via AS, whereas the other FT homolog PaFTL is not regulated by AS. In this study, the regulatory mechanism of the AS of PaFT was demonstrated to be conserved in different plant species. To define the distribution of the AS regulatory signals, the intron-swap, site-directed mutagenesis of alternative splice sites, and deletion experiment were performed. For the PaFT gene, all the signals that regulate the AS of the first intron were located within this intron, while the usage of the first alternative splice site in the second intron was determined by the first intron. Meanwhile, the AS of PaFDL1 might be co-regulated by exons and the first intron. Additionally, the first alternative splice site and adjacent region in PaFT intron 1 might contain cis-elements and/or RNA structures that affect the use of the other sites. This study had provided a deeper insight into the distribution of AS signals in plants, namely the AS signals of different splice sites might exist in the intron where the sites were present, and might also be distributed in exons or other introns.

Genome-Wide Association Analysis for Candidate Genes Contributing to Kernel-Related Traits in Maize

Maize grain size is the main factor determining grain yield. Dissecting the genetic basis of maize grain size may help reveal the regulatory mechanism of maize seed development and yield formation. In this study, two associated populations were used for genome-wide association analysis of kernel length, kernel width, kernel thickness, and hundred-kernel weight from multiple locations in AM122 and AM180, respectively. Then, genome-wide association mapping was performed based on the maize 6H90K SNP chip. A total of 139 loci were identified as associated with the four traits with p < 1 × 10^-4 using two models (FarmCPU and MLM). The transcriptome data showed that 15 of them were expressed differentially in two maize-inbred lines KB182 (small kernel) and KB020 (big kernel) during kernel development. These candidate genes were enriched in regulating peroxidase activity, oxidoreductase, and leaf senescence. The molecular function was major in binding and catalytic activity. This study provided important reference information for exploring maize kernel development mechanisms and applying molecular markers in high-yield breeding.

Emerging Function of Ecotype-Specific Splicing in the Recruitment of Commensal Microbiome

In recent years, host-microbiome interactions in both animals and plants has emerged as a novel research area for studying the relationship between host organisms and their commensal microbial communities. The fitness advantages of this mutualistic interaction can be found in both plant hosts and their associated microbiome, however, the driving forces mediating this beneficial interaction are poorly understood. Alternative splicing (AS), a pivotal post-transcriptional mechanism, has been demonstrated to play a crucial role in plant development and stress responses among diverse plant ecotypes. This natural variation of plants also has an impact on their commensal microbiome. In this article, we review the current progress of plant natural variation on their microbiome community, and discuss knowledge gaps between AS regulation of plants in response to their intimately related microbiota. Through the impact of this article, an avenue could be established to study the biological mechanism of naturally varied splicing isoforms on plant-associated microbiome assembly.

The Importance of a Genome-Wide Association Analysis in the Study of Alternative Splicing Mutations in Plants with a Special Focus on Maize

Alternative splicing is an important mechanism for regulating gene expressions at the post-transcriptional level. In eukaryotes, the genes are transcribed in the nucleus to produce pre-mRNAs and alternative splicing can splice a pre-mRNA to eventually form multiple different mature mRNAs, greatly increasing the number of genes and protein diversity. Alternative splicing is involved in the regulation of various plant life activities, especially the response of plants to abiotic stresses and is also an important process of plant growth and development. This review aims to clarify the usefulness of a genome-wide association analysis in the study of alternatively spliced variants by summarizing the application of alternative splicing, genome-wide association analyses and genome-wide association analyses in alternative splicing, as well as summarizing the related research progress.

The Intron Retention Variant CsClpP3m Is Involved in Leaf Chlorosis in Some Tea Cultivars

Tea products made from chlorotic or albino leaves are very popular for their unique flavor. Probing into the molecular mechanisms underlying the chlorotic leaf phenotype is required to better understand the formation of these tea cultivars and aid in future practical breeding. In this study, transcriptional alterations of multiple subunit genes of the caseinolytic protease complex (Clp) in the chlorotic tea cultivar 'Yu-Jin-Xiang' (YJX) were found. Cultivar YJX possessed the intron retention variant of ClpP3, named as CsClpP3m, in addition to the non-mutated ClpP3. The mutated variant results in a truncated protein containing only 166 amino acid residues and lacks the catalytic triad S182-H206-D255. Quantitative analysis of two CsClpP3 variants in different leaves with varying degrees of chlorosis in YJX and analyses of different chlorotic tea cultivars revealed that the transcript ratios of CsClpP3m over CsClpP3 were negatively correlated with leaf chlorophyll contents. The chlorotic young leaf phenotype was also generated in the transgenic tobacco by suppressing ClpP3 using the RNAi method; complementation with non-mutated CsClpP3 rescued the wild-type phenotype, whereas CsClpP3m failed to complement. Taken together, CsClpP3m is involved in leaf chlorosis in YJX and some other tea cultivars in a dose-dependent manner, likely resulting from the failure of Clp complex assembly due to the truncated sequence of CsClpP3m. Our data shed light on the mechanisms controlling leaf chlorosis in tea plants.

Processing body (P-body) and its mediators in cancer

In recent years, processing bodies (P-bodies) formed by liquid-liquid phase separation, have attracted growing scientific attention due to their involvement in numerous cellular activities, including the regulation of mRNAs decay or storage. These cytoplasmic dynamic membraneless granules contain mRNA storage and decay components such as deadenylase and decapping factors. In addition, different mRNA metabolic regulators, including m⁶A readers and gene-mediated miRNA-silencing, are also associated with such P-bodies. Cancerous cells may profit from these mRNA decay shredders by up-regulating the expression level of oncogenes and down-regulating tumor suppressor genes. The main challenges of cancer treatment are drug resistance, metastasis, and cancer relapse likely associated with cancer stem cells, heterogeneity, and plasticity features of different tumors. The mRNA metabolic regulators based on P-bodies play a great role in cancer development and progression. The dysregulation of P-bodies mediators affects mRNA metabolism. However, less is known about the relationship between P-bodies mediators and cancerous behavior. The current review summarizes the recent studies on P-bodies mediators, their contribution to tumor development, and their potential in the clinical setting, particularly highlighting the P-bodies as potential drug-carriers such as exosomes to anticancer in the future.

Alternative splicing dynamics and evolutionary divergence during embryogenesis in wheat species

Among polyploid species with complex genomic architecture, variations in the regulation of alternative splicing (AS) provide opportunities for transcriptional and proteomic plasticity and the potential for generating trait diversities. However, the evolution of AS and its influence on grain development in diploid grass and valuable polyploid wheat crops are poorly understood. To address this knowledge gap, we developed a pipeline for the analysis of alternatively spliced transcript isoforms, which takes the high sequence similarity among polyploid wheat subgenomes into account. Through analysis of synteny and detection of collinearity of homoeologous subgenomes, conserved and specific AS events across five wheat and grass species were identified. A global analysis of the regulation of AS in diploid grass and polyploid wheat grains revealed diversity in AS events not only between the endosperm, pericarp and embryo overdevelopment, but also between subgenomes. Analysis of AS in homoeologous triads of polyploid wheats revealed evolutionary divergence between gene-level and transcript-level regulation of embryogenesis. Evolutionary age analysis indicated that the generation of novel transcript isoforms has occurred in young genes at a more rapid rate than in ancient genes. These findings, together with the development of comprehensive AS resources for wheat and grass species, advance understanding of the evolution of regulatory features of AS during embryogenesis and grain development in wheat.

Genome-wide discovery of natural variation in pre-mRNA splicing and prioritising causal alternative splicing to salt stress response in rice

Pre-mRNA splicing is an essential step for the regulation of gene expression. In order to specifically capture splicing variants in plants for genome-wide association studies (GWAS), we developed a software tool to quantify and visualise Variations of Splicing in Population (VaSP). VaSP can quantify splicing variants from short-read RNA-seq datasets and discover genotype-specific splicing (GSS) events, which can be used to prioritise causal pre-mRNA splicing events in GWAS. We applied our method to an RNA-seq dataset with 328 samples from 82 genotypes from a rice diversity panel exposed to optimal and saline growing conditions. In total, 764 significant GSS events were identified in salt stress conditions. GSS events were used as markers for a GWAS with the shoot Na⁺ accumulation, which identified six GSS events in five genes significantly associated with the shoot Na⁺ content. Two of these genes, OsNUC1 and OsRAD23 emerged as top candidate genes with splice variants that exhibited significant divergence between the variants for shoot growth under salt stress conditions. VaSP is a versatile tool for alternative splicing analysis in plants and a powerful tool for prioritising candidate causal pre-mRNA splicing and corresponding genomic variations in GWAS.

Transcriptome and metabolome profiling provide insights into molecular mechanism of pseudostem elongation in banana

BACKGROUND

Banana plant height is an important trait for horticultural practices and semi-dwarf cultivars show better resistance to damages by wind and rain. However, the molecular mechanisms controlling the pseudostem height remain poorly understood. Herein, we studied the molecular changes in the pseudostem of a semi-dwarf banana mutant Aifen No. 1 (Musa spp. Pisang Awak sub-group ABB) as compared to its wild-type dwarf cultivar using a combined transcriptome and metabolome approach.

RESULTS

A total of 127 differentially expressed genes and 48 differentially accumulated metabolites were detected between the mutant and its wild type. Metabolites belonging to amino acid and its derivatives, flavonoids, lignans, coumarins, organic acids, and phenolic acids were up-regulated in the mutant. The transcriptome analysis showed the differential regulation of genes related to the gibberellin pathway, auxin transport, cell elongation, and cell wall modification. Based on the regulation of gibberellin and associated pathway-related genes, we discussed the involvement of gibberellins in pseudostem elongation in the mutant banana. Genes and metabolites associated with cell wall were explored and their involvement in cell extension is discussed.

CONCLUSIONS

The results suggest that gibberellins and associated pathways are possibly developing the observed semi-dwarf pseudostem phenotype together with cell elongation and cell wall modification. The findings increase the understanding of the mechanisms underlying banana stem height and provide new clues for further dissection of specific gene functions.

CRISPR-Cas technology in corn: a new key to unlock genetic knowledge and create novel products

Since its inception in 2012, CRISPR-Cas technologies have taken the life science community by storm. Maize genetics research is no exception. Investigators around the world have adapted CRISPR tools to advance maize genetics research in many ways. The principle application has been targeted mutagenesis to confirm candidate genes identified using map-based methods. Researchers are also developing tools to more effectively apply CRISPR-Cas technologies to maize because successful application of CRISPR-Cas relies on target gene identification, guide RNA development, vector design and construction, CRISPR-Cas reagent delivery to maize tissues, and plant characterization, each contributing unique challenges to CRISPR-Cas efficacy. Recent advances continue to chip away at major barriers that prevent more widespread use of CRISPR-Cas technologies in maize, including germplasm-independent delivery of CRISPR-Cas reagents and production of high-resolution genomic data in relevant germplasm to facilitate CRISPR-Cas experimental design. This has led to the development of novel breeding tools to advance maize genetics and demonstrations of how CRISPR-Cas technologies might be used to enhance maize germplasm.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01200-9.

Genome-wide identification, evolution, expression, and alternative splicing profiles of peroxiredoxin genes in cotton

Peroxiredoxin (PRX) is a ubiquitous thioredoxin-dependent peroxidase that can eliminate excessive free radicals produced by stress and protect cells from oxidative damage. PRXs are also involved in reactive oxygen species (ROS)- and redox-dependent signaling by performing redox interactions with other proteins and modify their redox status. At present, PRX family identification, evolution and regulation research has been conducted in some plants; however, systematic research about this family is lacking in cotton. In this study, a total of 44 PRXs were identified in the cotton genome. Phylogenetic and conserved active site analyses showed that the PRXs were divided into six subfamilies according to the conserved site (PxxxTxxC…S…W/F) and conserved cysteinyl residues positions. Segmental duplication and polyploid events were the main methods for PRX family expansion, and the PRXs of diploid G. arboreum were the donors of PRXs in the D subgenomes of allotetraploid G. hirsutum and G. barbadense during the evolution of the PRX family. qRT-PCR analysis confirmed that cis-acting elements play important roles in regulating the expression of PRXs. Alternative splicing events occurred in GhPRX14-D that can increased the complexity of transcripts in G. hirsutum. Subcellular localization showed that most PRX members were located in chloroplasts, the cytoplasmic membrane and the nucleus. Our results provide systematic support for a better understanding of PRXs in cotton and a starting point for further studies of the specific functions of PRXs in cotton.

Viral Perturbation of Alternative Splicing of a Host Transcript Benefits Infection

Pathogens disturb alternative splicing patterns of infected eukaryotic hosts. However, in plants it is unknown if this is incidental to infection or represents a pathogen-induced remodeling of host gene expression needed to support infection. Here, we compared changes in transcription and protein accumulation with changes in transcript splicing patterns in maize (Zea mays) infected with the globally important pathogen sugarcane mosaic virus (SCMV). Our results suggested that changes in alternative splicing play a major role in determining virus-induced proteomic changes. Focusing on maize phytoene synthase1 (ZmPSY1), which encodes the key regulatory enzyme in carotenoid biosynthesis, we found that although SCMV infection decreases total ZmPSY1 transcript accumulation, the proportion of splice variant T001 increases by later infection stages so that ZmPSY1 protein levels are maintained. We determined that ZmPSY1 has two leaf-specific transcripts, T001 and T003, distinguished by differences between the respective 3'-untranslated regions (UTRs). The shorter 3'-UTR of T001 makes it the more efficient mRNA. Nonsense ZmPSY1 mutants or virus-induced silencing of ZmPSY1 expression suppressed SCMV accumulation, attenuated symptoms, and decreased chloroplast damage. Thus, ZmPSY1 acts as a proviral host factor that is required for virus accumulation and pathogenesis. Taken together, our findings reveal that SCMV infection-modulated alternative splicing ensures that ZmPSY1 synthesis is sustained during infection, which supports efficient virus infection.

Natural variation at FLM splicing has pleiotropic effects modulating ecological strategies in Arabidopsis thaliana

Investigating the evolution of complex phenotypes and the underlying molecular bases of their variation is critical to understand how organisms adapt to their environment. Applying classical quantitative genetics on a segregating population derived from a Can-0xCol-0 cross, we identify the MADS-box transcription factor FLOWERING LOCUS M (FLM) as a player of the phenotypic variation in plant growth and color. We show that allelic variation at FLM modulates plant growth strategy along the leaf economics spectrum, a trade-off between resource acquisition and resource conservation, observable across thousands of plant species. Functional differences at FLM rely on a single intronic substitution, disturbing transcript splicing and leading to the accumulation of non-functional FLM transcripts. Associations between this substitution and phenotypic and climatic data across Arabidopsis natural populations, show how noncoding genetic variation at a single gene might be adaptive through pleiotropic effects.

FLOWERING LOCUS M (FLM) is known as a repressor of Arabidopsis flowering. Here, the authors show that a single intronic substitution of FLM modulates leaf color and plant growth strategy along the leaf economics spectrum, as well as plays a role in plant adaptation.

Identification of Ear Morphology Genes in Maize (Zea mays L.) Using Selective Sweeps and Association Mapping

The performance of maize hybrids largely depend on two parental inbred lines. Improving inbred lines using artificial selection is a key task in breeding programs. However, it is important to elucidate the effects of this selection on inbred lines. Altogether, 208 inbred lines from two maize heterosis groups, named Shaan A and Shaan B, were sequenced by the genotype-by-sequencing to detect genomic changes under selection pressures. In addition, we completed genome-wide association analysis in 121 inbred lines to identify candidate genes for ear morphology related traits. In a genome-wide selection scan, the inbred lines from Shaan A and Shaan B groups showed obvious population divergences and different selective signals distributed in 337 regions harboring 772 genes. Meanwhile, functional enrichment analysis showed those selected genes are mainly involved in regulating cell development. Interestingly, some ear morphology related traits showed significant differentiation between the inbred lines from the two heterosis groups. The genome-wide association analysis of ear morphology related traits showed that four associated genes were co-localized in the selected regions with high linkage disequilibrium. Our spatiotemporal pattern and gene interaction network results for the four genes further contribute to our understanding of the mechanisms behind ear and fruit length development. This study provides a novel insight into digging a candidate gene for complex traits using breeding materials. Our findings in relation to ear morphology will help accelerate future maize improvement.

Exploring transcriptional switches from pairwise, temporal and population RNA-Seq data using deepTS

Transcriptional switch (TS) is a widely observed phenomenon caused by changes in the relative expression of transcripts from the same gene, in spatial, temporal or other dimensions. TS has been associated with human diseases, plant development and stress responses. Its investigation is often hampered by a lack of suitable tools allowing comprehensive and flexible TS analysis for high-throughput RNA sequencing (RNA-Seq) data. Here, we present deepTS, a user-friendly web-based implementation that enables a fully interactive, multifunctional identification, visualization and analysis of TS events for large-scale RNA-Seq datasets from pairwise, temporal and population experiments. deepTS offers rich functionality to streamline RNA-Seq-based TS analysis for both model and non-model organisms and for those with or without reference transcriptome. The presented case studies highlight the capabilities of deepTS and demonstrate its potential for the transcriptome-wide TS analysis of pairwise, temporal and population RNA-Seq data. We believe deepTS will help research groups, regardless of their informatics expertise, perform accessible, reproducible and collaborative TS analyses of large-scale RNA-Seq data.

Molecular characterization and functional expression of voltage-gated sodium channel variants in Apolygus lucorum (Meyer-Dür)

BACKGROUND

Apolygus lucorum (Meyer-Dür) is a serious worldwide agricultural pest, especially for Bt cotton in China. Pyrethroids, neonicotinoids and organophosphates are the most effective insecticides to control piercing and sucking insects, including A. lucorum. The voltage-gated sodium channel (Na_v ) is major target site of pyrethroids. Extensive alternative splicing and RNA editing, two major post-transcriptional mechanisms, contribute to generate different functional sodium channel variants. In our research, we characterized the sodium channel variants of A. lucorum.

RESULTS

In this study, we isolated numerous sodium channel variants that cover the entire coding region of the VGSC gene from A. lucorum. All clones could be grouped into 47 splice types based on the presence of nine alternative exons (exons j, n, o, a, p, b, s, q and t). Exons j, b and t were located independently, while exons n, o, a and p were located adjacently, as were exons s and q. We also found 35 nucleotide changes in different positions in individual variants, of which 18 nucleotide changes were A-to-I RNA editing, 11 nucleotide changes were likely due to U-to-C or C-to-U editing, and the others were likely natural sequence polymorphisms in the population. Furthermore, we expressed all of the variants in Xenopus oocytes. Eighteen of them were expressed in oocytes and sensitive to tetrodotoxin.

CONCLUSION

Our results provide a functional basis for understanding how A. lucorum sodium channel variants work in regulating channel expression, pharmacology and gating properties for agricultural insects. Apolygus lucorum is widely distributed in cotton production. Our results suggest how AlNa_v (the sodium channel of A.lucorum) variants work in regulating channel expression, pharmacology and sodium channel gating for agricultural insects in the future. © 2020 Society of Chemical Industry.

TeoNAM: A Nested Association Mapping Population for Domestication and Agronomic Trait Analysis in Maize

Recombinant inbred lines (RILs) are an important resource for mapping genes controlling complex traits in many species. While RIL populations have been developed for maize, a maize RIL population with multiple teosinte inbred lines as parents has been lacking. Here, we report a teosinte nested association mapping (TeoNAM) population, derived from crossing five teosinte inbreds to the maize inbred line W22. The resulting 1257 BC1S4 RILs were genotyped with 51,544 SNPs, providing a high-density genetic map with a length of 1540 cM. On average, each RIL is 15% homozygous teosinte and 8% heterozygous. We performed joint linkage mapping (JLM) and a genome-wide association study (GWAS) for 22 domestication and agronomic traits. A total of 255 QTL from JLM were identified, with many of these mapping near known genes or novel candidate genes. TeoNAM is a useful resource for QTL mapping for the discovery of novel allelic variation from teosinte. TeoNAM provides the first report that PROSTRATE GROWTH1, a rice domestication gene, is also a QTL associated with tillering in teosinte and maize. We detected multiple QTL for flowering time and other traits for which the teosinte allele contributes to a more maize-like phenotype. Such QTL could be valuable in maize improvement.

Transcriptome-wide alternative splicing modulation during plant-pathogen interactions in wheat

Alternative splicing (AS) enhances the diversities of both transcripts and proteins in eukaryotes, which contribute to stress adaptation. To catalog wheat (Triticum aestivum L.) AS genes, we characterized 45 RNA-seq libraries from wheat seedlings infected by powdery mildew, Blumeria graminis f. sp. tritici (Bgt) or stripe rust fungus, Puccinia striiformis f. sp. tritici (Pst). We discovered that 11.2% and 10.4% of the multiexon genes had AS transcripts during Bgt and Pst infections, respectively. In response to fungal infection, wheat modulated AS not only in disease resistance proteins, but also in splicing related factors. Apart from the stress induced or activated splicing variants by pathogen, the differential expression profiles were fold increased through changing the ratio of full spliced transcripts versus intron retention (IR) transcripts. Comparing AS transcripts produced by the same gene in Bgt with Pst stress, the spliced terminal exons and the stranded introns are independent and different. This demonstrated that differential induction of specific splice variants were activated against two fungal pathogens. The specific induced AS genes in the Pst-resistant plants were enriched in improving the membrane permeability and protein modification ability, whereas gene expression involved in protein translation and transport were strengthened in Pst-susceptible plants.

Genome-Wide Identification of Splicing Quantitative Trait Loci (sQTLs) in Diverse Ecotypes of Arabidopsis thaliana

Alternative splicing (AS) of pre-mRNAs contributes to transcriptome diversity and enables plants to generate different protein isoforms from a single gene and/or fine-tune gene expression during different development stages and environmental changes. Although AS is pervasive, the genetic basis for differential isoform usage in plants is still emerging. In this study, we performed genome-wide analysis in 666 geographically distributed diverse ecotypes of Arabidopsis thaliana to identify genomic regions [splicing quantitative trait loci (sQTLs)] that may regulate differential AS. These ecotypes belong to different microclimatic conditions and are part of the relict and non-relict populations. Although sQTLs were spread across the genome, we observed enrichment for trans-sQTL (trans-sQTLs hotspots) on chromosome one. Furthermore, we identified several sQTL (911) that co-localized with trait-linked single nucleotide polymorphisms (SNP) identified in the Arabidopsis genome-wide association studies (AraGWAS). Many sQTLs were enriched among circadian clock, flowering, and stress-responsive genes, suggesting a role for differential isoform usage in regulating these important processes in diverse ecotypes of Arabidopsis. In conclusion, the current study provides a deep insight into SNPs affecting isoform ratios/genes and facilitates a better mechanistic understanding of trait-associated SNPs in GWAS studies. To the best of our knowledge, this is the first report of sQTL analysis in a large set of Arabidopsis ecotypes and can be used as a reference to perform sQTL analysis in the Brassicaceae family. Since whole genome and transcriptome datasets are available for these diverse ecotypes, it could serve as a powerful resource for the biological interpretation of trait-associated loci, splice isoform ratios, and their phenotypic consequences to help produce more resilient and high yield crop varieties.

De Novo Domestication: An Alternative Route toward New Crops for the Future

Current global agricultural production must feed over 7 billion people. However, productivity varies greatly across the globe and is under threat from both increased competitions for land and climate change and associated environmental deterioration. Moreover, the increase in human population size and dietary changes are putting an ever greater burden on agriculture. The majority of this burden is met by the cultivation of a very small number of species, largely in locations that differ from their origin of domestication. Recent technological advances have raised the possibility of de novo domestication of wild plants as a viable solution for designing ideal crops while maintaining food security and a more sustainable low-input agriculture. Here we discuss how the discovery of multiple key domestication genes alongside the development of technologies for accurate manipulation of several target genes simultaneously renders de novo domestication a route toward crops for the future.

Crop genome-wide association study: a harvest of biological relevance

With the advent of rapid genotyping and next-generation sequencing technologies, genome-wide association study (GWAS) has become a routine strategy for decoding genotype-phenotype associations in many species. More than 1000 such studies over the last decade have revealed substantial genotype-phenotype associations in crops and provided unparalleled opportunities to probe functional genomics. Beyond the many 'hits' obtained, this review summarizes recent efforts to increase our understanding of the genetic architecture of complex traits by focusing on non-main effects including epistasis, pleiotropy, and phenotypic plasticity. We also discuss how these achievements and the remaining gaps in our knowledge will guide future studies. Synthetic association is highlighted as leading to false causality, which is prevalent but largely underestimated. Furthermore, validation evidence is appealing for future GWAS, especially in the context of emerging genome-editing technologies.