A high resolution map of the Arabidopsis thaliana developmental transcriptome based on RNA-seq profiling

From tetraploid to diploid, a pangenomic approach to identify genes lost during synthetic diploidization of Eragrostis curvula

INTRODUCTION

In Eragrostis curvula, commonly known as weeping lovegrass, a synthetic diploidization event of the facultative apomictic tetraploid Tanganyika INTA cv. originated from the sexual diploid Victoria cv. Apomixis is an asexual reproduction by seeds in which the progeny is genetically identical to the maternal plant.

METHODS

To assess the genomic changes related to ploidy and to the reproductive mode occurring during diploidization, a mapping approach was followed to obtain the first E. curvula pangenome assembly. In this way, gDNA of Tanganyika INTA was extracted and sequenced in 2x250 Illumina pair-end reads and mapped against the Victoria genome assembly. The unmapped reads were used for variant calling, while the mapped reads were assembled using Masurca software.

RESULTS

The length of the assembly was 28,982,419 bp distributed in 18,032 contigs, and the variable genes annotated in these contigs rendered 3,952 gene models. Functional annotation of the genes showed that the reproductive pathway was differentially enriched. PCR amplification in gDNA and cDNA of Tanganyika INTA and Victoria was conducted to validate the presence/absence variation in five genes related to reproduction and ploidy. The polyploid nature of the Tanganyika INTA genome was also evaluated through the variant calling analysis showing the single nucleotide polymorphism (SNP) coverage and allele frequency distribution with a segmental allotetraploid pairing behavior.

DISCUSSION

The results presented here suggest that the genes were lost in Tanganyika INTA during the diploidization process that was conducted to suppress the apomictic pathway, affecting severely the fertility of Victoria cv.

Identification of structural variations related to drought tolerance in wheat (Triticum aestivum L.)

Structural variations are common in plant genomes, affecting meiotic recombination and distorted segregation in wheat. And presence/absence variations can significantly affect drought tolerance in wheat. Drought is a major abiotic stress limiting wheat production. Common wheat has a complex genome with three sub-genomes, which host large numbers of structural variations (SVs). SVs play critical roles in understanding the genetic contributions of plant domestication and phenotypic plasticity, but little is known about their genomic characteristics and their effects on drought tolerance. In the present study, high-resolution karyotypes of 180 doubled haploids (DHs) were developed. Signal polymorphisms between the parents involved with 8 presence-absence variations (PAVs) of tandem repeats (TR) distributed on the 7 (2A, 4A, 5A, 7A, 3B, 7B, and 2D) of 21 chromosomes. Among them, PAV on chromosome 2D showed distorted segregation, others transmit normal conforming to a 1:1 segregation ration in the population; and a PAVs recombination occurred on chromosome 2A. Association analysis of PAV and phenotypic traits under different water regimes, we found PAVs on chromosomes 4A, 5A, and 7B showed negative effect on grain length (GL) and grain width (GW); PAV.7A had opposite effect on grain thickness (GT) and spike length (SL), with the effect on traits differing under different water regimes. PAVs on linkage group 2A, 4A, 7A, 2D, and 7B associated with the drought tolerance coefficients (DTCs), and significant negative effect on drought resistance values (D values) were detected in PAV.7B. Additionally, quantitative trait loci (QTL) associated with phenotypic traits using the 90 K SNP array showed QTL for DTCs and grain-related traits in chromosomes 4A, and 5A, 3B were co-localized in differential regions of PAVs. These PAVs can cause the differentiation of the target region of SNP and could be used for genetic improvement of agronomic traits under drought stress via marker-assisted selection (MAS) breeding.

Identification of new potential downstream transcriptional targets of the strigolactone pathway including glucosinolate biosynthesis

Strigolactones regulate shoot branching and many aspects of plant growth, development, and allelopathy. Strigolactones are often discussed alongside auxin because they work together to inhibit shoot branching. However, the roles and mechanisms of strigolactones and how they act independently of auxin are still elusive. Additionally, there is still much in general to be discovered about the network of molecular regulators and their interactions in response to strigolactones. Here, we conducted an experiment in Arabidopsis with physiological treatments and strigolactone mutants to determine transcriptional pathways associated with strigolactones. The three physiological treatments included shoot tip removal with and without auxin treatment and treatment of intact plants with the auxin transport inhibitor, N-1-naphthylphthalamic acid (NPA). We identified the glucosinolate biosynthesis pathway as being upregulated across strigolactone mutants indicating strigolactone-glucosinolate crosstalk. Additionally, strigolactone application cannot restore the highly branched phenotype observed in glucosinolate biosynthesis mutants, placing glucosinolate biosynthesis downstream of strigolactone biosynthesis. Oxidative stress genes were enriched across the experiment suggesting that this process is mediated through multiple hormones. Here, we also provide evidence supporting non-auxin-mediated, negative feedback on strigolactone biosynthesis. Increases in strigolactone biosynthesis gene expression seen in strigolactone mutants could not be fully restored by auxin. By contrast, auxin could fully restore auxin-responsive gene expression increases, but not sugar signaling-related gene expression. Our data also point to alternative roles of the strigolactone biosynthesis genes and potential new signaling functions of strigolactone precursors. In this study, we identify a strigolactone-specific regulation of glucosinolate biosynthesis genes indicating that the two are linked and may work together in regulating stress and shoot ranching responses in Arabidopsis. Additionally, we provide evidence for non-auxinmediated feedback on strigolactone biosynthesis and discuss this in the context of sugar signaling.

Charting plant gene functions in the multi-omics and single-cell era

Despite the increased access to high-quality plant genome sequences, the set of genes with a known function remains far from complete. With the advent of novel bulk and single-cell omics profiling methods, we are entering a new era where advanced and highly integrative functional annotation strategies are being developed to elucidate the functions of all plant genes. Here, we review different multi-omics approaches to improve functional and regulatory gene characterization and highlight the power of machine learning and network biology to fully exploit the complementary information embedded in different omics layers. Finally, we discuss the potential of emerging single-cell methods and algorithms to further increase the resolution, allowing generation of functional insights about plant biology.

The Arabidopsis DUF239 gene family encodes Neprosin-like proteins that are widely expressed in seed endosperm

Domain of unknown function 239 (DUF239) is a conserved sequence found in the catalytic site of Neprosins which are specific secreted prolyl endopeptidases found in the Nepenthes genus. Neprosins participate in the nitrogen cycle by digesting preys trapped in the pitcher of these carnivorous plants. Apart from that, DUF239s have been poorly documented in plants. We have identified 50 genes containing DUF239-coding sequences in the Arabidopsis genome that are distributed across six distinct phylogenetic clusters. The chromosomal distribution suggests that several genes are the result of recent duplication events, with up to eight genes found in a strict tandem distribution. In Arabidopsis, most of DUF239-containing sequences are also associated to a Neprosin-activating domain (DUF4409) and an amino-terminal α-helix which corresponds to the typical domain organization of the Neprosins described in the Nepenthes genus. Analysis of Arabidopsis transcriptomic datasets reveals that 39 genes are exclusively expressed in reproductive organs, mainly during seed development and more specifically in the endosperm (23 genes). The peculiar expression pattern of the DUF239 gene family in Arabidopsis suggests new functions of Neprosin-like proteins in plants during seed development.

Promoter activity and transcriptome analyses decipher functions of CgbHLH001 gene (Chenopodium glaucum L.) in response to abiotic stress

BACKGROUND

Our previous studies revealed that CgbHLH001 transcription factor (TF) played an important role in abiotic stress tolerance, suggesting that its promoter was a potential target in response to stress signals. In addition, the regulatory mechanism of CgbHLH001 TF is still limited.

RESULTS

In the present study, a 1512 bp of 5'-flanking sequence of CgbHLH001 gene was identified, and the sequence carried quite a few of cis-acting elements. The gene promoter displayed strong activity and was induced by multiple abiotic stress. A series of 5'-deletions of the promoter sequence resulted in a gradual decrease in its activity, especially, the 5' untranslated region (UTR) was necessary to drive promoter activity. Further, CgbHLH001 promoter drove its own gene overexpression ectopically at the transcriptional and translational levels, which in turn conferred the stress tolerance to transgenic Arabidopsis. Transcriptome analysis showed that salt stress induced a large number of genes involved in multiple biological regulatory processes. Differentially expressed genes (DEGs) that mediate phytohormone signal transduction and mitogen-activated protein kinase (MAPK) signaling pathway were widely induced and mostly upregulated under salt stress, and the transcription levels in P_bHLH::bHLH-overexpressing transgenic lines were higher than that of 35S::bHLH overexpression.

CONCLUSIONS

The CgbHLH001 promoter exhibited a positive response to abiotic stress and its 5' UTR sequence enhanced the regulation of gene expression to stress. A few important pathways and putative key genes involved in salt tolerance were identified, which can be used to elucidate the mechanism of salt tolerance and decipher the regulatory mechanism of promoters to develop an adaptation strategy for desert halophytes.

Genome-wide screening of the RNase T2 gene family and functional analyses in jujube (Ziziphus jujuba Mill.)

BACKGROUND

Ribonuclease (RNase T2) plays crucial roles in plant evolution and breeding. However, there have been few studies on the RNase T2 gene family in Ziziphus jujuba Mill., one of important dried fruit tree species. Recently, the released sequences of the reference genome of jujube provide a good chance to perform genome-wide identification and characterization of ZjRNase gene family in the jujube.

RESULTS

In this study, we identified four members of RNase T2 in jujube distributed on three chromosomes and unassembled chromosomes. They all contained two conserved sites (CASI and CASII). Analysis of the phylogenetic relationships revealed that the RNase T2 genes in jujube could be divided into two groups: ZjRNase1 and ZjRNase2 belonged to class I, while ZjRNase3 and ZjRNase4 belonged to class II. Only ZjRNase1 and ZjRNase2 expression were shown by the jujube fruit transcriptome analysis. So ZjRNase1 and ZjRNase2 were selected functional verification by overexpression transformation of Arabidopsis. The overexpression of these two genes led to an approximately 50% reduction in seed number, which deserve further attention. Moreover, the leaves of the ZjRNase1 overexpression transgenic lines were curled and twisted. Overexpression of ZjRNase2 resulted in shortened and crisp siliques and the production of trichomes, and no seeds were produced.

CONCLUSION

In summary, these findings will provide new insights into the molecular mechanisms of low number of hybrid seeds in jujube and a reference for the future molecular breeding of jujube.

Analyses of Cullin1 homologs reveal functional redundancy in S-RNase-based self-incompatibility and evolutionary relationships in eudicots

In Petunia (Solanaceae family), self-incompatibility (SI) is regulated by the polymorphic S-locus, which contains the pistil-specific S-RNase and multiple pollen-specific S-Locus F-box (SLF) genes. SLFs assemble into E3 ubiquitin ligase complexes known as Skp1-Cullin1-F-box complexes (SCFSLF). In pollen tubes, these complexes collectively mediate ubiquitination and degradation of all nonself S-RNases, but not self S-RNase, resulting in cross-compatible, but self-incompatible, pollination. Using Petunia inflata, we show that two pollen-expressed Cullin1 (CUL1) proteins, PiCUL1-P and PiCUL1-B, function redundantly in SI. This redundancy is lost in Petunia hybrida, not because of the inability of PhCUL1-B to interact with SSK1, but due to a reduction in the PhCUL1-B transcript level. This is possibly caused by the presence of a DNA transposon in the PhCUL1-B promoter region, which was inherited from Petunia axillaris, one of the parental species of Pe. hybrida. Phylogenetic and syntenic analyses of Cullin genes in various eudicots show that three Solanaceae-specific CUL1 genes share a common origin, with CUL1-P dedicated to S-RNase-related reproductive processes. However, CUL1-B is a dispersed duplicate of CUL1-P present only in Petunia, and not in the other species of the Solanaceae family examined. We suggest that the CUL1s involved (or potentially involved) in the SI response in eudicots share a common origin.

Shoot maturation strengthens FLS2-mediated resistance to Pseudomonas syringae

A temporal-spatial regulation of immunity components is essential for properly activating plant defense response. Flagellin-sensing 2 (FLS2) is a surface-localized receptor that recognizes bacterial flagellin. The immune function of FLS2 is compromised in early stages of shoot development. However, the underlying mechanism for the age-dependent FLS2 signaling is not clear. Here, we show that the reduced basal immunity of juvenile leaves against Pseudomonas syringae pv. tomato DC3000 is independent of FLS2. The flg22-induced marker gene expression and ROS activation were comparable in juvenile and adult stage, but callose deposition was more evident in the adult stage than that of juvenile stage. We further demonstrated that microRNA156, a master regulator of plant aging, suppressed callose deposition in juvenile leaves in response to flg22 but not the expression of FLS2 and FRK1 (Flg22-induced receptor-like kinase 1) . Altogether, we revealed an intrinsic mechanism that regulates the amplitude of FLS2-mediated resistance during aging.

Gaining Insight into Large Gene Families with the Aid of Bioinformatic Tools

Proteins participating in plant cell morphogenesis are often encoded by large gene families, in some cases comprising paralogs with variable (modular) domain organization, as in the case of the formin (FH2 protein) family of actin nucleators that can have also additional functions. Unravelling the phylogeny of such a complex gene family brings a number of specific challenges but may be crucial for predictions of protein function and for experimental design. Here we present an overview of our "cottage industry" semi-manual bioinformatic approach, based mostly, though not exclusively, on freely available software tools, which we used to obtain insight into the evolutionary history of plant FH2 proteins and some other components of the plant cell morphogenesis apparatus.

Alternative Splicing in the Regulatory Circuit of Plant Temperature Response

As sessile organisms, plants have evolved complex mechanisms to rapidly respond to ever-changing ambient temperatures. Temperature response in plants is modulated by a multilayer regulatory network, including transcriptional and post-transcriptional regulations. Alternative splicing (AS) is an essential post-transcriptional regulatory mechanism. Extensive studies have confirmed its key role in plant temperature response, from adjustment to diurnal and seasonal temperature changes to response to extreme temperatures, which has been well documented by previous reviews. As a key node in the temperature response regulatory network, AS can be modulated by various upstream regulations, such as chromatin modification, transcription rate, RNA binding proteins, RNA structure and RNA modifications. Meanwhile, a number of downstream mechanisms are affected by AS, such as nonsense-mediated mRNA decay (NMD) pathway, translation efficiency and production of different protein variants. In this review, we focus on the links between splicing regulation and other mechanisms in plant temperature response. Recent advances regarding how AS is regulated and the following consequences in gene functional modulation in plant temperature response will be discussed. Substantial evidence suggests that a multilayer regulatory network integrating AS in plant temperature response has been unveiled.

A Revised View of the LSU Gene Family: New Functions in Plant Stress Responses and Phytohormone Signaling

LSUs (RESPONSE TO LOW SULFUR) are plant-specific proteins of unknown function that were initially identified during transcriptomic studies of the sulfur deficiency response in Arabidopsis. Recent functional studies have shown that LSUs are important hubs of protein interaction networks with potential roles in plant stress responses. In particular, LSU proteins have been reported to interact with members of the brassinosteroid, jasmonate signaling, and ethylene biosynthetic pathways, suggesting that LSUs may be involved in response to plant stress through modulation of phytohormones. Furthermore, in silico analysis of the promoter regions of LSU genes in Arabidopsis has revealed the presence of cis-regulatory elements that are potentially responsive to phytohormones such as ABA, auxin, and jasmonic acid, suggesting crosstalk between LSU proteins and phytohormones. In this review, we summarize current knowledge about the LSU gene family in plants and its potential role in phytohormone responses.

A high-throughput yeast approach to characterize aquaporin permeabilities: Profiling the Arabidopsis PIP aquaporin sub-family

INTRODUCTION

Engineering membrane transporters to achieve desired functionality is reliant on availability of experimental data informing structure-function relationships and intelligent design. Plant aquaporin (AQP) isoforms are capable of transporting diverse substrates such as signaling molecules, nutrients, metalloids, and gases, as well as water. AQPs can act as multifunctional channels and their transport function is reliant on many factors, with few studies having assessed transport function of specific isoforms for multiple substrates.

METHODS

High-throughput yeast assays were developed to screen for transport function of plant AQPs, providing a platform for fast data generation and cataloguing of substrate transport profiles. We applied our high-throughput growth-based yeast assays to screen all 13 Arabidopsis PIPs (AtPIPs) for transport of water and several neutral solutes: hydrogen peroxide (H2O2), boric acid (BA), and urea. Sodium (Na+) transport was assessed using elemental analysis techniques.

RESULTS

All AtPIPs facilitated water and H2O2 transport, although their growth phenotypes varied, and none were candidates for urea transport. For BA and Na+ transport, AtPIP2;2 and AtPIP2;7 were the top candidates, with yeast expressing these isoforms having the most pronounced toxicity response to BA exposure and accumulating the highest amounts of Na+. Linking putative AtPIP isoform substrate transport profiles with phylogenetics and gene expression data, enabled us to align possible substrate preferences with known and hypothesized biological roles of AtPIPs.

DISCUSSION

This testing framework enables efficient cataloguing of putative transport functionality of diverse AQPs at a scale that can help accelerate our understanding of AQP biology through big data approaches (e.g. association studies). The principles of the individual assays could be further adapted to test additional substrates. Data generated from this framework could inform future testing of AQP physiological roles, and address knowledge gaps in structure-function relationships to improve engineering efforts.

The FANCC-FANCE-FANCF complex is evolutionarily conserved and regulates meiotic recombination

At meiosis, programmed meiotic DNA double-strand breaks are repaired via homologous recombination, resulting in crossovers (COs). From a large excess of DNA double-strand breaks that are formed, only a small proportion gets converted into COs because of active mechanisms that restrict CO formation. The Fanconi anemia (FA) complex proteins AtFANCM, MHF1 and MHF2 were previously identified in a genetic screen as anti-CO factors that function during meiosis in Arabidopsis thaliana. Here, pursuing the same screen, we identify FANCC as a new anti-CO gene. FANCC was previously only identified in mammals because of low primary sequence conservation. We show that FANCC, and its physical interaction with FANCE-FANCF, is conserved from vertebrates to plants. Further, we show that FANCC, together with its subcomplex partners FANCE and FANCF, regulates meiotic recombination. Mutations of any of these three genes partially rescues CO-defective mutants, which is particularly marked in female meiosis. Functional loss of FANCC, FANCE, or FANCF results in synthetic meiotic catastrophe with the pro-CO factor MUS81. This work reveals that FANCC is conserved outside mammals and has an anti-CO role during meiosis together with FANCE and FANCF.

Tissue-preferential recruitment of electron transfer chains for cytochrome P450-catalyzed phenolic biosynthesis

Cytochrome P450 system consists of P450 monooxygenase and redox pattern(s). While the importance of monooxygenases in plant metabolism is well documented, the metabolic roles of the related redox components have been largely overlooked. Here, we show that distinct electron transfer chains are recruited in phenylpropanoid-monolignol P450 systems to support the synthesis and distribution of different classes of phenolics in different plant tissues. While Arabidopsis cinnamate 4-hydroxylase adopts conventional NADPH-cytochrome P450 oxidoreductase (CPR) electron transfer chain for its para-hydroxylation reaction, ferulate 5-hydroxylase uses both NADPH-CPR-cytochrome b₅ (CB5) and NADH-cytochrome b₅ reductase-CB5 chains to support benzene ring 5-hydroxylation, in which the former route is primarily recruited in the stem for syringyl lignin synthesis, while the latter dominates in the syntheses of 5-hydroxylated phenolics in seeds and seed coat suberin. Our study unveils an additional layer of complexity and versatility of P450 system that the plants evolved for diversifying phenolic repertoires.

Pollen Coat Proteomes of Arabidopsis thaliana, Arabidopsis lyrata, and Brassica oleracea Reveal Remarkable Diversity of Small Cysteine-Rich Proteins at the Pollen-Stigma Interface

The pollen coat is the outermost domain of the pollen grain and is largely derived from the anther tapetum, which is a secretory tissue that degenerates late in pollen development. By being localised at the interface of the pollen-stigma interaction, the pollen coat plays a central role in mediating early pollination events, including molecular recognition. Amongst species of the Brassicaceae, a growing body of data has revealed that the pollen coat carries a range of proteins, with a number of small cysteine-rich proteins (CRPs) being identified as important regulators of the pollen-stigma interaction. By utilising a state-of-the-art liquid chromatography/tandem mass spectrometry (LC-MS/MS) approach, rich pollen coat proteomic profiles were obtained for Arabidopsis thaliana, Arabidopsis lyrata, and Brassica oleracea, which greatly extended previous datasets. All three proteomes revealed a strikingly large number of small CRPs that were not previously reported as pollen coat components. The profiling also uncovered a wide range of other protein families, many of which were enriched in the pollen coat proteomes and had functions associated with signal transduction, cell walls, lipid metabolism and defence. These proteomes provide an excellent source of molecular targets for future investigations into the pollen-stigma interaction and its potential evolutionary links to plant-pathogen interactions.

Spliceostatin C, a component of a microbial bioherbicide, is a potent phytotoxin that inhibits the spliceosome

Spliceostatin C (SPC) is a component of a bioherbicide isolated from the soil bacterium Burkholderia rinojensis. The chemical structure of SPC closely resembles spliceostatin A (SPA) which was characterized as an anticancer agent and splicing inhibitor. SPC inhibited the growth of Arabidopsis thaliana seedlings with an IC50 value of 2.2 µM. The seedlings exposed to SPC displayed a significant response with decreased root length and number and inhibition of gravitropism. Reverse transcriptase semi-quantitative PCR (RT-sqPCR) analyses of 19 selected genes demonstrated the active impact of SPC on the quality and quantity of transcripts that underwent intron rearrangements as well as up or down expression upon exposure to SPC. Qualitative and quantitative proteomic profiles identified 66 proteins that were significantly affected by SPC treatment. Further proteomics data analysis revealed that spliceostatin C induces hormone-related responses in Arabidopsis seedlings. In silico binding studies showed that SPC binds to a pocket between the SF3B3 and PF5A of the spliceosome.

Silencing of a Pectin Acetylesterase (PAE) Gene Highly Expressed in Tobacco Pistils Negatively Affects Pollen Tube Growth

Successful plant reproduction and fruit formation depend on adequate pollen and pistil development, and pollen-pistil interactions. In Nicotiana tabacum, pollen tubes grow through the intercellular spaces of pistil-specialized tissues, stigmatic secretory zone, and stylar transmitting tissue (STT). These intercellular spaces are supposed to be formed by the modulation of cell wall pectin esterification. Previously we have identified a gene preferentially expressed in pistils encoding a putative pectin acetylesterase (PAE), named NtPAE1. Here, we characterized the NtPAE1 gene and performed genome-wide and phylogenetic analyses of PAEs. We identified 30 PAE sequences in the N. tabacum genome, distributed in four clades. The expression of NtPAE1 was assessed by RT-qPCR and in situ hybridization. We confirmed NtPAE1 preferential expression in stigmas/styles and ovaries and demonstrated its high expression in the STT. Structural predictions and comparisons between NtPAE1 and functional enzymes validated its identity as a PAE. Transgenic plants were produced, overexpressing and silencing the NtPAE1 gene. Overexpressed plants displayed smaller flowers while silencing plants exhibited collapsed pollen grains, which hardly germinate. NtPAE1 silencing plants do not produce fruits, due to impaired pollen tube growth in their STTs. Thus, NtPAE1 is an essential enzyme regulating pectin modifications in flowers and, ultimately, in plant reproduction.

Cytosolic phosphoglucose isomerase is essential for microsporogenesis and embryogenesis in Arabidopsis

Phosphoglucose isomerase (PGI) catalyzes the interconversion of fructose-6-phosphate and glucose-6-phosphate, which impacts cell carbon metabolic flow. Arabidopsis (Arabidopsis thaliana) contains two nuclear PGI genes respectively encoding plastidial PGI1 and cytosolic PGI (cPGI). The loss of PGI1 impairs the conversion of F6P of the Calvin-Benson cycle to G6P for the synthesis of transitory starch in leaf chloroplasts. Since cpgi knockout mutants have not yet been obtained, they are thought to be lethal. The cpgi lethality can be rescued by expressing CaMV 35S promoter (p35S)-driven cPGI; however, the complemented line is completely sterile due to pollen degeneration. Here, we generated a cpgi mutant expressing p35S::cPGI-YFP in which YFP fluorescence in developing anthers was undetectable specifically in the tapetum and in pollen, which could be associated with male sterility. We also generated RNAi-cPGI knockdown lines with strong cPGI repression in floral buds that exhibited reduced male fertility due to the degeneration of most pollen. Histological analyses indicated that the synthesis of intersporal callose walls was impaired, causing microsporocytes to fail to separate haploid daughter nuclei to form tetrads, which might be responsible for subsequent pollen degeneration. We successfully isolated cpgi knockout mutants in the progeny of a heterozygous cpgi mutant floral-dipped with sugar solutions. The rescued cpgi mutants exhibited diminished young vegetative growth, reduced female fertility, and impaired intersporal callose wall formation in a meiocyte, and, thus, male sterility. Collectively, our data suggest that cPGI plays a vital role in carbohydrate partitioning, which is indispensable for microsporogenesis and early embryogenesis.

Isolation of Nuclei Tagged in Specific Cell Types (INTACT) in Arabidopsis

Many functionally distinct plant tissues have relatively low numbers of cells that are embedded within complex tissues. For example, the shoot apical meristem (SAM) consists of a small population of pluripotent stem cells surrounded by developing leaves and/or flowers at the growing tip of the plant. It is technically challenging to collect enough high-quality SAM samples for molecular analyses. Isolation of Nuclei Tagged in specific Cell Types (INTACT) is an easily reproducible method that allows the enrichment of biotin-tagged cell-type-specific nuclei from the total nuclei pool using biotin-streptavidin affinity purification. Here, we provide a detailed INTACT protocol for isolating nuclei from the Arabidopsis SAM. One can also adapt this protocol to isolate nuclei from other tissues and cell types for investigating tissue/cell-type-specific transcriptome and epigenome and their changes during developmental programs at a high spatiotemporal resolution. Furthermore, due to its low cost and simple procedures, INTACT can be conducted in any standard molecular laboratory.

Canalization of genome-wide transcriptional activity in Arabidopsis thaliana accessions by MET1-dependent CG methylation

BACKGROUND

Despite its conserved role on gene expression and transposable element (TE) silencing, genome-wide CG methylation differs substantially between wild Arabidopsis thaliana accessions.

RESULTS

To test our hypothesis that global reduction of CG methylation would reduce epigenomic, transcriptomic, and phenotypic diversity in A. thaliana accessions, we knock out MET1, which is required for CG methylation, in 18 early-flowering accessions. Homozygous met1 mutants in all accessions suffer from common developmental defects such as dwarfism and delayed flowering, in addition to accession-specific abnormalities in rosette leaf architecture, silique morphology, and fertility. Integrated analysis of genome-wide methylation, chromatin accessibility, and transcriptomes confirms that MET1 inactivation greatly reduces CG methylation and alters chromatin accessibility at thousands of loci. While the effects on TE activation are similarly drastic in all accessions, the quantitative effects on non-TE genes vary greatly. The global expression profiles of accessions become considerably more divergent from each other after genome-wide removal of CG methylation, although a few genes with diverse expression profiles across wild-type accessions tend to become more similar in mutants. Most differentially expressed genes do not exhibit altered chromatin accessibility or CG methylation in cis, suggesting that absence of MET1 can have profound indirect effects on gene expression and that these effects vary substantially between accessions.

CONCLUSIONS

Systematic analysis of MET1 requirement in different A. thaliana accessions reveals a dual role for CG methylation: for many genes, CG methylation appears to canalize expression levels, with methylation masking regulatory divergence. However, for a smaller subset of genes, CG methylation increases expression diversity beyond genetically encoded differences.

The PROSCOOP10 Gene Encodes Two Extracellular Hydroxylated Peptides and Impacts Flowering Time in Arabidopsis

The Arabidopsis PROSCOOP genes belong to a family predicted to encode secreted pro-peptides, which undergo maturation steps to produce peptides named SCOOP. Some of them are involved in defence signalling through their perception by a receptor complex including MIK2, BAK1 and BKK1. Here, we focused on the PROSCOOP10 gene, which is highly and constitutively expressed in aerial organs. The MS/MS analyses of leaf apoplastic fluids allowed the identification of two distinct peptides (named SCOOP10#1 and SCOOP10#2) covering two different regions of PROSCOOP10. They both possess the canonical S-X-S family motif and have hydroxylated prolines. This identification in apoplastic fluids confirms the biological reality of SCOOP peptides for the first time. NMR and molecular dynamics studies showed that the SCOOP10 peptides, although largely unstructured in solution, tend to assume a hairpin-like fold, exposing the two serine residues previously identified as essential for the peptide activity. Furthermore, PROSCOOP10 mutations led to an early-flowering phenotype and increased expression of the floral integrators SOC1 and LEAFY, consistent with the de-regulated transcription of PROSCOOP10 in several other mutants displaying early- or late-flowering phenotypes. These results suggest a role for PROSCOOP10 in flowering time, highlighting the functional diversity within the PROSCOOP family.

Analysis of BRCT5 domain-containing proteins reveals a new component of DNA damage repair in Arabidopsis

The integrity of plant genetic information is constantly challenged by various internal and external factors. Therefore, plants use a sophisticated molecular network to identify, signal and repair damaged DNA. Here, we report on the identification and analysis of four uncharacterized Arabidopsis BRCT5 DOMAIN CONTAINING PROTEINs (BCPs). Proteins with the BRCT5 domain are frequently involved in the maintenance of genome stability across eukaryotes. The screening for sensitivity to induced DNA damage identified BCP1 as the most interesting candidate. We show that BCP1 loss of function mutants are hypersensitive to various types of DNA damage and accumulate an increased number of dead cells in root apical meristems upon DNA damage. Analysis of publicly available sog1 transcriptomic and SOG1 genome-wide DNA binding data revealed that BCP1 is inducible by gamma radiation and is a direct target of this key DNA damage signaling transcription factor. Importantly, bcp1 plants showed a reduced frequency of somatic homologous recombination in response to both endogenous and induced DNA damage. Altogether, we identified a novel plant-specific DNA repair factor that acts downstream of SOG1 in homology-based repair.

New insights on the expression patterns of specific Arabinogalactan proteins in reproductive tissues of Arabidopsis thaliana

Arabinogalactan proteins (AGPs) are hydroxyproline-rich glycoproteins containing a high proportion of carbohydrates, widely distributed in the plant kingdom and ubiquitously present in land plants. AGPs have long been suggested to play important roles in plant reproduction and there is already evidence that specific glycoproteins are essential for male and female gametophyte development, pollen tube growth and guidance, and successful fertilization. However, the functions of many of these proteins have yet to be uncovered, mainly due to the difficulty to study individual AGPs. In this work, we generated molecular tools to analyze the expression patterns of a subgroup of individual AGPs in different Arabidopsis tissues, focusing on reproductive processes. This study focused on six AGPs: four classical AGPs (AGP7, AGP25, AGP26, AGP27), one AG peptide (AGP24) and one chimeric AGP (AGP31). These AGPs were first selected based on their predicted expression patterns along the reproductive tissues from available RNA-seq data. Promoter analysis using β-glucuronidase fusions and qPCR in different Arabidopsis tissues allowed to confirm these predictions. AGP7 was mainly expressed in female reproductive tissues, more precisely in the style, funiculus, and integuments near the micropyle region. AGP25 was found to be expressed in the style, septum and ovules with higher expression in the chalaza and funiculus tissues. AGP26 was present in the ovules and pistil valves. AGP27 was expressed in the transmitting tissue, septum and funiculus during seed development. AGP24 was expressed in pollen grains, in mature embryo sacs, with highest expression at the chalazal pole and in the micropyle. AGP31 was expressed in the mature embryo sac with highest expression at the chalaza and, occasionally, in the micropyle. For all these AGPs a co-expression analysis was performed providing new hints on its possible functions. This work confirmed the detection in Arabidopsis male and female tissues of six AGPs never studied before regarding the reproductive process. These results provide novel evidence on the possible involvement of specific AGPs in plant reproduction, as strong candidates to participate in pollen-pistil interactions in an active way, which is significant for this field of study.

Revisiting AGAMOUS-LIKE15, a Key Somatic Embryogenesis Regulator, Using Next Generation Sequencing Analysis in Arabidopsis

AGAMOUS-like 15 (AGL15) is a member of the MADS-domain transcription factor (TF) family. MADS proteins are named for a conserved domain that was originally from an acronym derived from genes expressed in a variety of eukaryotes (MCM1-AGAMOUS-DEFICIENS-SERUM RESPONSE FACTOR). In plants, this family has expanded greatly, with more than one-hundred members generally found in dicots, and the proteins encoded by these genes have often been associated with developmental identity. AGL15 transcript and protein accumulate primarily in embryos and has been found to promote an important process called plant regeneration via somatic embryogenesis (SE). To understand how this TF performs this function, we have previously used microarray technologies to assess direct and indirect responsive targets of this TF. We have now revisited this question using next generation sequencing (NGS) to both characterize in vivo binding sites for AGL15 as well as response to the accumulation of AGL15. We compared these data to the prior microarray results to evaluate the different platforms. The new NGS data brought to light an interaction with brassinosteroid (BR) hormone signaling that was "missed" in prior Gene Ontology analysis from the microarray studies.

Active DNA demethylation in plants: 20 years of discovery and beyond

Maintaining proper DNA methylation levels in the genome requires active demethylation of DNA. However, removing the methyl group from a modified cytosine is chemically difficult and therefore, the underlying mechanism of demethylation had remained unclear for many years. The discovery of the first eukaryotic DNA demethylase, Arabidopsis thaliana REPRESSOR OF SILENCING 1 (ROS1), led to elucidation of the 5-methylcytosine base excision repair mechanism of active DNA demethylation. In the 20 years since ROS1 was discovered, our understanding of this active DNA demethylation pathway, as well as its regulation and biological functions in plants, has greatly expanded. These exciting developments have laid the groundwork for further dissecting the regulatory mechanisms of active DNA demethylation, with potential applications in epigenome editing to facilitate crop breeding and gene therapy.

Biochemical characterization of Pectin Methylesterase Inhibitor 3 from Arabidopsis thaliana

The de-methylesterification of the pectic polysaccharide homogalacturonan (HG) by pectin methylesterases (PMEs) is a critical step in the control of plant cell expansion and morphogenesis. Plants have large gene families encoding PMEs but also PME inhibitors (PMEIs) with differ in their biochemical properties. The Arabidopsis thaliana PECTIN METHYLESTERASE INHIBITOR 3 (PMEI3) gene is frequently used as a tool to manipulate pectin methylesterase activity in studies assessing its role in the control of morphogenesis. One limitation of these studies is that the exact biochemical activity of this protein has not yet been determined. In this manuscript we produced the protein in Pichia pastoris and characterized its activity in vitro. Like other PMEIs, PMEI3 inhibits PME activity at acidic pH in a variety of cell wall extracts and in purified PME preparations, but does not affect the much stronger PME activity at neutral pH. The protein is remarkable heat stable and shows higher activity against PME3 than against PME2, illustrating how different members of the large PMEI family can differ in their specificities towards PME targets. Finally, growing Arabidopsis thaliana seedlings in the presence of purified PMEI3 caused a dose-dependent inhibition of root growth associated with the overall inhibition of HG de-methylesterification of the root surface. This suggests an essential in vivo role for PME activity at acidic pH in HG de-methylesterification and growth control. These results show that purified recombinant PMEI3 is a powerful tool to study the connection between pectin de-methylesterification and cell expansion.

Polymerization of the backbone of the pectic polysaccharide rhamnogalacturonan I

Rhamnogalacturonan I (RG-I) is a major plant cell wall pectic polysaccharide defined by its repeating disaccharide backbone structure of [4)-α-D-GalA-(1,2)-α-L-Rha-(1,]. A family of RG-I:Rhamnosyltransferases (RRT) has previously been identified, but synthesis of the RG-I backbone has not been demonstrated in vitro because the identity of Rhamnogalacturonan I:Galaturonosyltransferase (RG-I:GalAT) was unknown. Here a putative glycosyltransferase, At1g28240/MUCI70, is shown to be an RG-I:GalAT. The name RGGAT1 is proposed to reflect the catalytic activity of this enzyme. When incubated together with the rhamnosyltransferase RRT4, the combined activities of RGGAT1 and RRT4 result in elongation of RG-I acceptors in vitro into a polymeric product. RGGAT1 is a member of a new GT family categorized as GT116, which does not group into existing GT-A clades and is phylogenetically distinct from the GALACTURONOSYLTRANSFERASE (GAUT) family of GalA transferases that synthesize the backbone of the pectin homogalacturonan. RGGAT1 has a predicted GT-A fold structure but employs a metal-independent catalytic mechanism that is rare among glycosyltransferases with this fold type. The identification of RGGAT1 and the 8-member Arabidopsis GT116 family provides a new avenue for studying the mechanism of RG-I synthesis and the function of RG-I in plants.

Evolution of isoform-level gene expression patterns across tissues during lotus species divergence

Both gene duplication and alternative splicing (AS) drive the functional diversity of gene products in plants, yet the relative contributions of the two key mechanisms to the evolution of gene function are largely unclear. Here, we studied AS in two closely related lotus plants, Nelumbo lutea and Nelumbo nucifera, and the outgroup Arabidopsis thaliana, for both single-copy and duplicated genes. We show that most splicing events evolved rapidly between orthologs and that the origin of lineage-specific splice variants or isoforms contributed to gene functional changes during species divergence within Nelumbo. Single-copy genes contain more isoforms, have more AS events conserved across species, and show more complex tissue-dependent expression patterns than their duplicated counterparts. This suggests that expression divergence through isoforms is a mechanism to extend the expression breadth of genes with low copy numbers. As compared to isoforms of local, small-scale duplicates, isoforms of whole-genome duplicates are less conserved and display a less conserved tissue bias, pointing towards their contribution to subfunctionalization. Through comparative analysis of isoform expression networks, we identified orthologous genes of which the expression of at least some of their isoforms displays a conserved tissue bias across species, indicating a strong selection pressure for maintaining a stable expression pattern of these isoforms. Overall, our study shows that both AS and gene duplication contributed to the diversity of gene function during the evolution of lotus.

Finding new friends and revisiting old ones - how plant lipid droplets connect with other subcellular structures

The number of described contact sites between different subcellular compartments and structures in eukaryotic cells has increased dramatically in recent years and, as such, has substantially reinforced the well-known premise that these kinds of connections are essential for overall cellular organization and the proper functioning of cellular metabolic and signaling pathways. Here, we discuss contact sites involving plant lipid droplets (LDs), including LD-endoplasmic reticulum (ER) connections that mediate the biogenesis of new LDs at the ER, LD-peroxisome connections, that facilitate the degradation of LD-stored triacylglycerols (TAGs), and the more recently discovered LD-plasma membrane connections, which involve at least three novel proteins, but have a yet unknown physiological function(s).

Enhanced growth rate under elevated CO2 conditions was observed for transgenic lines of genes identified by intraspecific variation analyses in Arabidopsis thaliana

Using the whole genome and growth data of Arabidopsis thaliana ecotypes, we identified two genes associated with enhancement of the growth rate in response to elevated CO₂ conditions. Improving plant growth under elevated CO₂ conditions may contribute to enhanced agricultural yield under future global climate change. In this study, we examined the genes implicated in the enhancement of growth rates under elevated CO₂ conditions by analyzing the growth rates of Arabidopsis thaliana ecotypes originating from various latitudes and altitudes throughout the world. We also performed a genome-wide association study and a transcriptome study to identify single nucleic polymorphisms that were correlated with the relative growth rate (RGR) under elevated CO₂ conditions or with CO₂ response of RGR. We then selected 43 candidate genes and generated their overexpression and/or RNA interference (RNAi) transgenic mutants for screening. After screening, we have found that RNAi lines of AT3G4000 and AT5G50900 showed significantly higher growth rates under the elevated CO₂ condition. As per our findings, we conclude that natural variation includes genetic variation associated with the enhancement of plant productivity under elevated CO₂ conditions.

High-gradient magnetic fields and starch metabolism: results from a space experiment

Directing plant growth in weightlessness requires understanding the processes that establish plant orientation and how to manipulate them. Both gravi- and phototropism determine directional growth and previous experiments showed that high gradient magnetic fields (HGMF) can induce curvature in roots and shoots. Experiments with Brassica rapa verified that that gravitropism-like induction of curvature is possible in space and that the HGMF-responsive organelles are amyloplasts. We assessed the effect of space and HGMF based on 16 genes and compared their transcription with static growth and clinorotation. Amyloplasts size in root tips increased under weightlessness but decreased under clinorotation but not in response to magnetic fields. Amyloplast size changes were correlated with reduced amylase transcription in space samples and enhanced transcription after clinorotation. Mechanostimulation and weightlessness have opposite effects on the size of amyloplasts. The data show that plants perceive weightlessness, and that their metabolism adjusts to microgravity and mechanostimulation. Thus, clinorotation as surrogate for space research may lead to incorrect interpretations.

MCU proteins dominate in vivo mitochondrial Ca2+ uptake in Arabidopsis roots

Ca2+ signaling is central to plant development and acclimation. While Ca2+-responsive proteins have been investigated intensely in plants, only a few Ca2+-permeable channels have been identified, and our understanding of how intracellular Ca2+ fluxes is facilitated remains limited. Arabidopsis thaliana homologs of the mammalian channel-forming mitochondrial calcium uniporter (MCU) protein showed Ca2+ transport activity in vitro. Yet, the evolutionary complexity of MCU proteins, as well as reports about alternative systems and unperturbed mitochondrial Ca2+ uptake in knockout lines of MCU genes, leave critical questions about the in vivo functions of the MCU protein family in plants unanswered. Here, we demonstrate that MCU proteins mediate mitochondrial Ca2+ transport in planta and that this mechanism is the major route for fast Ca2+ uptake. Guided by the subcellular localization, expression, and conservation of MCU proteins, we generated an mcu triple knockout line. Using Ca2+ imaging in living root tips and the stimulation of Ca2+ transients of different amplitudes, we demonstrated that mitochondrial Ca2+ uptake became limiting in the triple mutant. The drastic cell physiological phenotype of impaired subcellular Ca2+ transport coincided with deregulated jasmonic acid-related signaling and thigmomorphogenesis. Our findings establish MCUs as a major mitochondrial Ca2+ entry route in planta and link mitochondrial Ca2+ transport with phytohormone signaling.

An LCO-responsive homolog of NODULE INCEPTION positively regulates lateral root formation in Populus sp

The transcription factor NODULE INCEPTION (NIN) has been studied extensively for its multiple roles in root nodule symbiosis within plants of the nitrogen-fixing clade (NFC) that associate with soil bacteria, such as rhizobia and Frankia. However, NIN homologs are present in plants outside the NFC, suggesting a role in other developmental processes. Here, we show that the biofuel crop Populus sp., which is not part of the NFC, contains eight copies of NIN with diversified protein sequence and expression patterns. Lipo-chitooligosaccharides (LCOs) are produced by rhizobia and a wide range of fungi, including mycorrhizal ones, and act as symbiotic signals that promote lateral root formation. RNAseq analysis of Populus sp. treated with purified LCO showed induction of the PtNIN2 subfamily. Moreover, the expression of PtNIN2b correlated with the formation of lateral roots and was suppressed by cytokinin treatment. Constitutive expression of PtNIN2b overcame the inhibition of lateral root development by cytokinin under high nitrate conditions. Lateral root induction in response to LCOs likely represents an ancestral function of NIN retained and repurposed in nodulating plants, as we demonstrate that the role of NIN in LCO-induced root branching is conserved in both Populus sp. and legumes. We further established a visual marker of LCO perception in Populus sp. roots, the putative sulfotransferase PtSS1 that can be used to study symbiotic interactions with the bacterial and fungal symbionts of Populus sp.

Glutamate 1-semialdehyde aminotransferase is connected to GluTR by GluTR-binding protein and contributes to the rate-limiting step of 5-aminolevulinic acid synthesis

Tetrapyrroles play fundamental roles in crucial processes including photosynthesis, respiration, and catalysis. In plants, 5-aminolevulinic acid (ALA) is the common precursor of tetrapyrroles. ALA is synthesized from activated glutamate by the enzymes glutamyl-tRNA reductase (GluTR) and glutamate-1-semialdehyde aminotransferase (GSAAT). ALA synthesis is recognized as the rate-limiting step in this pathway. We aimed to explore the contribution of GSAAT to the control of ALA synthesis and the formation of a protein complex with GluTR. In Arabidopsis thaliana, two genes encode GSAAT isoforms: GSA1 and GSA2. A comparison of two GSA knockout mutants with the wild-type revealed the correlation of reduced GSAAT activity and ALA-synthesizing capacity in leaves with lower chlorophyll content. Growth and green pigmentation were more severely impaired in gsa2 than in gsa1, indicating the predominant role of GSAAT2 in ALA synthesis. Interestingly, GluTR accumulated to higher levels in gsa2 than in the wild-type and was mainly associated with the plastid membrane. We propose that the GSAAT content modulates the amount of soluble GluTR available for ALA synthesis. Several different biochemical approaches revealed the GSAAT-GluTR interaction through the assistance of GluTR-binding protein (GBP). A modeled structure of the tripartite protein complex indicated that GBP mediates the stable association of GluTR and GSAAT for adequate ALA synthesis.

Prioritized candidate causal haplotype blocks in plant genome-wide association studies

Genome wide association studies (GWAS) can play an essential role in understanding genetic basis of complex traits in plants and animals. Conventional SNP-based linear mixed models (LMM) that marginally test single nucleotide polymorphisms (SNPs) have successfully identified many loci with major and minor effects in many GWAS. In plant, the relatively small population size in GWAS and the high genetic diversity found in many plant species can impede mapping efforts on complex traits. Here we present a novel haplotype-based trait fine-mapping framework, HapFM, to supplement current GWAS methods. HapFM uses genotype data to partition the genome into haplotype blocks, identifies haplotype clusters within each block, and then performs genome-wide haplotype fine-mapping to prioritize the candidate causal haplotype blocks of trait. We benchmarked HapFM, GEMMA, BSLMM, GMMAT, and BLINK in both simulated and real plant GWAS datasets. HapFM consistently resulted in higher mapping power than the other GWAS methods in high polygenicity simulation setting. Moreover, it resulted in smaller mapping intervals, especially in regions of high LD, achieved by prioritizing small candidate causal blocks in the larger haplotype blocks. In the Arabidopsis flowering time (FT10) datasets, HapFM identified four novel loci compared to GEMMA’s results, and the average mapping interval of HapFM was 9.6 times smaller than that of GEMMA. In conclusion, HapFM is tailored for plant GWAS to result in high mapping power on complex traits and improved on mapping resolution to facilitate crop improvement.

Genome-wide association studies (GWAS) are commonly used in human and plant studies to identify genetic variants responsible for the phenotype of interest and provide foundations for studying disease mechanisms and crop improvement. Most GWAS models are developed and optimized using human datasets. However, the difference between human and plant datasets essentially limits their applications in plant studies, especially when mapping complex traits such as drought resistance and yield. In this study, we present a novel GWAS method, HapFM, tailored for plant datasets to overcome the difficulties of many conventional GWAS methods. HapFM resulted in higher statistical power than conventional GWAS methods for mapping complex traits in our simulation and real dataset analyses. In addition, HapFM reduced the mapping interval by prioritizing candidate causal regions in the genome, which benefits the downstream experimental studies. Last but not least, HapFM can incorporate biological annotations to increase statistical power further. Overall, HapFM balances statistical power, result interpretability, and downstream experimental verifiability.

Modulation of receptor-like transmembrane kinase 1 nuclear localization by DA1 peptidases in Arabidopsis

Signals are often perceived by proteins in one cellular location and transduced to other locations such as the nucleus. Signaling proteins can be cleaved by peptidases to facilitate this movement, but the peptidases involved in this are poorly understood despite their widespread role. We describe a role for the ubiquitin-activated peptidase DA1 in cleaving the membrane-localized receptor-like kinase transmembrane kinase 1 (TMK1) in Arabidopsis. TMK1 is phosphorylated in response to auxin and mediates several auxin responses including growth induction by cell expansion. DA1-mediated cleavage of TMK1 facilitates nuclear localization of its intracellular kinase domain to repress auxin-mediated gene expression, facilitating differential cell expansion during growth. These analyses establish a wider role for DA1 family activities in cell growth.

The cleavage of intracellular domains of receptor-like kinases (RLKs) has an important functional role in the transduction of signals from the cell surface to the nucleus in many organisms. However, the peptidases that catalyze protein cleavage during signal transduction remain poorly understood despite their crucial roles in diverse signaling processes. Here, we report in the flowering plant Arabidopsis thaliana that members of the DA1 family of ubiquitin-regulated Zn metallopeptidases cleave the cytoplasmic kinase domain of transmembrane kinase 1 (TMK1), releasing it for nuclear localization where it represses auxin-responsive cell growth during apical hook formation by phosphorylation and stabilization of the transcriptional repressors IAA32 and IAA34. Mutations in DA1 family members exhibited reduced apical hook formation, and DA1 family-mediated cleavage of TMK1 was promoted by auxin treatment. Expression of the DA1 family-generated intracellular kinase domain of TMK1 by an auxin-responsive promoter fully restored apical hook formation in a tmk1 mutant, establishing the function of DA1 family peptidase activities in TMK1-mediated differential cell growth and apical hook formation. DA1 family peptidase activity therefore modulates TMK1 kinase activity between a membrane location where it stimulates acid cell growth and initiates an auxin-dependent kinase cascade controlling cell proliferation in lateral roots and a nuclear localization where it represses auxin-mediated gene expression and growth.

A seed-like proteome in oil-rich tubers

There are numerous examples of plant organs or developmental stages that are desiccation-tolerant and can withstand extended periods of severe water loss. One prime example are seeds and pollen of many spermatophytes. However, in some plants, also vegetative organs can be desiccation-tolerant. One example are the tubers of yellow nutsedge (Cyperus esculentus), which also store large amounts of lipids similar to seeds. Interestingly, the closest known relative, purple nutsedge (Cyperus rotundus), generates tubers that do not accumulate oil and are not desiccation-tolerant. We generated nanoLC-MS/MS-based proteomes of yellow nutsedge in five replicates of four stages of tuber development and compared them to the proteomes of roots and leaves, yielding 2257 distinct protein groups. Our data reveal a striking upregulation of hallmark proteins of seeds in the tubers. A deeper comparison to the tuber proteome of the close relative purple nutsedge (C. rotundus) and a previously published proteome of Arabidopsis seeds and seedlings indicates that indeed a seed-like proteome was found in yellow but not purple nutsedge. This was further supported by an analysis of the proteome of a lipid droplet-enriched fraction of yellow nutsedge, which also displayed seed-like characteristics. One reason for the differences between the two nutsedge species might be the expression of certain transcription factors homologous to ABSCISIC ACID INSENSITIVE3, WRINKLED1, and LEAFY COTYLEDON1 that drive gene expression in Arabidopsis seed embryos.

Shedding light on AT1G29480 of Arabidopsis thaliana-An enigmatic locus restricted to Brassicacean genomes

A key feature of C₄ Kranz anatomy is the presence of an enlarged, photosynthetically highly active bundle sheath whose cells contain large numbers of chloroplasts. With the aim to identify novel candidate regulators of C₄ bundle sheath development, we performed an activation tagging screen with Arabidopsis thaliana. The reporter gene used encoded a chloroplast-targeted GFP protein preferentially expressed in the bundle sheath, and the promoter of the C₄ phosphoenolpyruvate carboxylase gene from Flaveria trinervia served as activation tag because of its activity in all chlorenchymatous tissues of A. thaliana. Primary mutants were selected based on their GFP signal intensity, and one stable mutant named kb-1 with a significant increase in GFP fluorescence intensity was obtained. Despite the increased GFP signal, kb-1 showed no alterations to bundle sheath anatomy. The causal locus, AT1G29480, is specific to the Brassicaceae with its second exon being conserved. Overexpression and reconstitution studies confirmed that AT1G29480, and specifically its second exon, were sufficient for the enhanced GFP phenotype, which was not dependent on translation of the locus or its parts into protein. We conclude, therefore, that the AT1G29480 locus enhances the GFP reporter gene activity via an RNA-based mechanism.

Seed dormancy 4 like1 of Arabidopsis is a key regulator of phase transition from embryo to vegetative development

Seed dormancy is an adaptive trait that enables plants to survive adverse conditions and restart growth in a season and location suitable for vegetative and reproductive growth. Control of seed dormancy is also important for crop production and food quality because it can help induce uniform germination and prevent preharvest sprouting. Rice preharvest sprouting quantitative trait locus analysis has identified Seed dormancy 4 (Sdr4) as a positive regulator of dormancy development. Here, we analyzed the loss-of-function mutant of the Arabidopsis ortholog, Sdr4 Like1 (SFL1), and found that the sfl1-1 seeds showed precocious germination at the mid- to late-maturation stage similar to rice sdr4 mutant, but converted to become more dormant than the wild type during maturation drying. Coordinated with the dormancy levels, expression levels of the seed maturation and dormancy master regulator genes, ABI3, FUS3, and DOG1 in sfl1-1 seeds were lower than in wild type at early- and mid-maturation stages, but higher at the late-maturation stage. In addition to the seed dormancy phenotype, sfl1-1 seedlings showed a growth arrest phenotype and heterochronic expression of LAFL (LEC1, ABI3, FUS3, LEC2) and DOG1 in the seedlings. These data suggest that SFL1 is a positive regulator of initiation and termination of the seed dormancy program. We also found genetic interaction between SFL1 and the SFL2, SFL3, and SFL4 paralogs of SFL1, which impacts on the timing of the phase transition from embryo maturation to seedling growth.

A combined transcriptome - miRNAome approach revealed that a kinesin gene is differentially targeted by a novel miRNA in an apomictic genotype of Eragrostis curvula

Weeping lovegrass (Eragrostis curvula [Shrad.] Nees) is a perennial grass typically established in semi-arid regions, with good adaptability to dry conditions and sandy soils. This polymorphic complex includes both sexual and apomictic cytotypes, with different ploidy levels (2x-8x). Diploids are known to be sexual, while most polyploids are facultative apomicts, and full apomicts have also been reported. Plant breeding studies throughout the years have focused on achieving the introgression of apomixis into species of agricultural relevance, but, given the complexity of the trait, a deeper understanding of the molecular basis of regulatory mechanisms of apomixis is still required. Apomixis is thought to be associated with silencing or disruption of the sexual pathway, and studies have shown it is influenced by epigenetic mechanisms. In a previous study, we explored the role of miRNA-mRNA interactions using two contrasting E. curvula phenotypes. Here, the sexual OTA-S, the facultative Don Walter and the obligate apomictic Tanganyika cDNA and sRNA libraries were inquired, searching for miRNA discovery and miRNA expression regulation of genes related to the reproductive mode. This allowed for the characterization of seven miRNAs and the validation of their miRNA-target interactions. Interestingly, a kinesin gene was found to be repressed in the apomictic cultivar Tanganyika, targeted by a novel miRNA that was found to be overexpressed in this genotype, suggestive of an involvement in the reproductive mode expression. Our work provided additional evidence of the contribution of the epigenetic regulation of the apomictic pathway.

Circadian autonomy and rhythmic precision of the Arabidopsis female reproductive organ

The plant circadian clock regulates essential biological processes including flowering time or petal movement. However, little is known about how the clock functions in flowers. Here, we identified the circadian components and transcriptional networks contributing to the generation of rhythms in pistils, the female reproductive organ. When detached from the rest of the flower, pistils sustain highly precise rhythms, indicating organ-specific circadian autonomy. Analyses of clock mutants and chromatin immunoprecipitation assays showed distinct expression patterns and specific regulatory functions for clock activators and repressors in pistils. Genetic interaction studies also suggested a hierarchy of the repressing activities that provide robustness and precision to the pistil clock. Globally, the circadian function in pistils primarily governs responses to environmental stimuli and photosynthesis and controls pistil growth and seed weight and production. Understanding the circadian intricacies in reproductive organs may prove useful for optimizing plant reproduction and productivity.

Metabolism and Signaling of Plant Mitochondria in Adaptation to Environmental Stresses

The interaction of mitochondria with cellular components evolved differently in plants and mammals; in plants, the organelle contains proteins such as ALTERNATIVE OXIDASES (AOXs), which, in conjunction with internal and external ALTERNATIVE NAD(P)H DEHYDROGENASES, allow canonical oxidative phosphorylation (OXPHOS) to be bypassed. Plant mitochondria also contain UNCOUPLING PROTEINS (UCPs) that bypass OXPHOS. Recent work revealed that OXPHOS bypass performed by AOXs and UCPs is linked with new mechanisms of mitochondrial retrograde signaling. AOX is functionally associated with the NO APICAL MERISTEM transcription factors, which mediate mitochondrial retrograde signaling, while UCP1 can regulate the plant oxygen-sensing mechanism via the PRT6 N-Degron. Here, we discuss the crosstalk or the independent action of AOXs and UCPs on mitochondrial retrograde signaling associated with abiotic stress responses. We also discuss how mitochondrial function and retrograde signaling mechanisms affect chloroplast function. Additionally, we discuss how mitochondrial inner membrane transporters can mediate mitochondrial communication with other organelles. Lastly, we review how mitochondrial metabolism can be used to improve crop resilience to environmental stresses. In this respect, we particularly focus on the contribution of Brazilian research groups to advances in the topic of mitochondrial metabolism and signaling.

Altered collective mitochondrial dynamics in the Arabidopsis msh1 mutant compromising organelle DNA maintenance

Mitochondria form highly dynamic populations in the cells of plants (and almost all eukaryotes). The characteristics and benefits of this collective behaviour, and how it is influenced by nuclear features, remain to be fully elucidated. Here, we use a recently developed quantitative approach to reveal and analyse the physical and collective 'social' dynamics of mitochondria in an Arabidopsis msh1 mutant where the organelle DNA maintenance machinery is compromised. We use a newly created line combining the msh1 mutant with mitochondrially targeted green fluorescent protein (GFP), and characterize mitochondrial dynamics with a combination of single-cell time-lapse microscopy, computational tracking, and network analysis. The collective physical behaviour of msh1 mitochondria is altered from that of the wild type in several ways: mitochondria become less evenly spread, and networks of inter-mitochondrial encounters become more connected, with greater potential efficiency for inter-organelle exchange-reflecting a potential compensatory mechanism for the genetic challenge to the mitochondrial DNA population, supporting more inter-organelle exchange. We find that these changes are similar to those observed in friendly, where mitochondrial dynamics are altered by a physical perturbation, suggesting that this shift to higher connectivity may reflect a general response to mitochondrial challenges, where physical dynamics of mitochondria may be altered to control the genetic structure of the mtDNA population.

Subfunctionalization of D27 Isomerase Genes in Saffron

Chromoplasts and chloroplasts contain carotenoid pigments as all-trans- and cis-isomers, which function as accessory light-harvesting pigments, antioxidant and photoprotective agents, and precursors of signaling molecules and plant hormones. The carotenoid pathway involves the participation of different carotenoid isomerases. Among them, D27 is a β-carotene isomerase showing high specificity for the C9-C10 double bond catalyzing the interconversion of all-trans- into 9-cis-β-carotene, the precursor of strigolactones. We have identified one D27 (CsD27-1) and two D27-like (CsD27-2 and CsD27-3) genes in saffron, with CsD27-1 and CsD27-3, clearly differing in their expression patterns; specifically, CsD27-1 was mainly expressed in the undeveloped stigma and roots, where it is induced by Rhizobium colonization. On the contrary, CsD27-2 and CsD27-3 were mainly expressed in leaves, with a preferential expression of CsD27-3 in this tissue. In vivo assays show that CsD27-1 catalyzes the isomerization of all-trans- to 9-cis-β-carotene, and could be involved in the isomerization of zeaxanthin, while CsD27-3 catalyzes the isomerization of all-trans- to cis-ζ-carotene and all-trans- to cis-neurosporene. Our data show that CsD27-1 and CsD27-3 enzymes are both involved in carotenoid isomerization, with CsD27-1 being specific to chromoplast/amyloplast-containing tissue, and CsD27-3 more specific to chloroplast-containing tissues. Additionally, we show that CsD27-1 is co-expressed with CCD7 and CCD8 mycorrhized roots, whereas CsD27-3 is expressed at higher levels than CRTISO and Z-ISO and showed circadian regulation in leaves. Overall, our data extend the knowledge about carotenoid isomerization and their implications in several physiological and ecological processes.

WAKL8 Regulates Arabidopsis Stem Secondary Wall Development

Wall-associated kinases/kinase-likes (WAKs/WAKLs) are plant cell surface sensors. A variety of studies have revealed the important functions of WAKs/WAKLs in regulating cell expansion and defense in cells with primary cell walls. Less is known about their roles during the development of the secondary cell walls (SCWs) that are present in xylem vessel (XV) and interfascicular fiber (IF) cells. In this study, we used RNA-seq data to screen Arabidopsis thaliana WAKs/WAKLs members that may be involved in SCW development and identified WAKL8 as a candidate. We obtained T-DNA insertion mutants wakl8-1 (inserted at the promoter region) and wakl8-2 (inserted at the first exon) and compared the phenotypes to wild-type (WT) plants. Decreased WAKL8 transcript levels in stems were found in the wakl8-2 mutant plants, and the phenotypes observed included reduced stem length and thinner walls in XV and IFs compared with those in the WT plants. Cell wall analysis showed no significant changes in the crystalline cellulose or lignin content in mutant stems compared with those in the WT. We found that WAKL8 had alternative spliced versions predicted to have only extracellular regions, which may interfere with the function of the full-length version of WAKL8. Our results suggest WAKL8 can regulate SCW thickening in Arabidopsis stems.

Unsupervised and semi-supervised learning: the next frontier in machine learning for plant systems biology

Advances in high-throughput omics technologies are leading plant biology research into the era of big data. Machine learning (ML) performs an important role in plant systems biology because of its excellent performance and wide application in the analysis of big data. However, to achieve ideal performance, supervised ML algorithms require large numbers of labeled samples as training data. In some cases, it is impossible or prohibitively expensive to obtain enough labeled training data; here, the paradigms of unsupervised learning (UL) and semi-supervised learning (SSL) play an indispensable role. In this review, we first introduce the basic concepts of ML techniques, as well as some representative UL and SSL algorithms, including clustering, dimensionality reduction, self-supervised learning (self-SL), positive-unlabeled (PU) learning and transfer learning. We then review recent advances and applications of UL and SSL paradigms in both plant systems biology and plant phenotyping research. Finally, we discuss the limitations and highlight the significance and challenges of UL and SSL strategies in plant systems biology.

Rapid evolution of pollen and pistil traits as a response to sexual selection in the post-pollination phase of mating

Sexual selection is the basis of some of the most striking phenotypic variation in nature.^1,2 In animals, sexual selection in males can act on traits that improve access to mates prior to copulation,^3-8 but also on sperm traits filtered by sperm competition,^9-14 or female choice expressed simply by the morphology and physiology of genital tracts.^14-16 Although long overlooked as a mode of selection on plant traits, sexual selection should act on land plants too because they are anisogamous: males produce more, and smaller, gametes than females.^17-19 Numerical asymmetry in gamete production is thought to play a central role in selection on traits that affect pollen transfer to mates,^20,21 but very little is known about how pollen competition or cryptic female choice might affect the evolution of traits expressed after pollination.^22,23 Here, we report the divergence of pollen and pistil traits of the dioecious wind-pollinated annual herb Mercurialis annua during evolution over three generations between populations at low versus high plant density, corresponding to low versus higher levels of polyandry;²⁴ we expected selection under higher polyandry to strengthen competition among pollen donors for fertilizing ovules. We found that populations at high density evolved faster-growing pollen tubes (an equivalent of greater sperm velocity), greater expression of pollen proteins involved in pollen growth, and larger stigmas (a trait likely enhancing the number of pollen donors and thus competition for ovules). Our results identify the post-pollination phase of plant mating as an important arena for the action of sexual selection.

Identification and functional annotation of long intergenic non-coding RNAs in Brassicaceae

Long intergenic noncoding RNAs (lincRNAs) are a large yet enigmatic class of eukaryotic transcripts that can have critical biological functions. The wealth of RNA-sequencing (RNA-seq) data available for plants provides the opportunity to implement a harmonized identification and annotation effort for lincRNAs that enables cross-species functional and genomic comparisons as well as prioritization of functional candidates. In this study, we processed >24 Tera base pairs of RNA-seq data from >16,000 experiments to identify ∼130,000 lincRNAs in four Brassicaceae: Arabidopsis thaliana, Camelina sativa, Brassica rapa, and Eutrema salsugineum. We used nanopore RNA-seq, transcriptome-wide structural information, peptide data, and epigenomic data to characterize these lincRNAs and identify conserved motifs. We then used comparative genomic and transcriptomic approaches to highlight lincRNAs in our data set with sequence or transcriptional conservation. Finally, we used guilt-by-association analyses to assign putative functions to lincRNAs within our data set. We tested this approach on a subset of lincRNAs associated with germination and seed development, observing germination defects for Arabidopsis lines harboring T-DNA insertions at these loci. LincRNAs with Brassicaceae-conserved putative miRNA binding motifs, small open reading frames, or abiotic-stress modulated expression are a few of the annotations that will guide functional analyses into this cryptic portion of the transcriptome.

Transcriptome analysis revealed the expression levels of genes related to abscisic acid and auxin biosynthesis in grapevine (Vitis vinifera L.) under root restriction

The root system is essential for the stable growth of plants. Roots help anchor plants in the soil and play a crucial role in water uptake, mineral nutrient absorption and endogenous phytohormone formation. Root-restriction (RR) cultivation, a powerful technique, confines plant roots to a specific soil space. In the present study, roots of one-year-old "Muscat Hamburg" grapevine under RR and control (nR) treatments harvested at 70 and 125 days after planting were used for transcriptome sequencing, and in total, 2031 (nR7 vs. nR12), 1445 (RR7 vs. RR12), 1532 (nR7 vs. RR7), and 2799 (nR12 vs. RR12) differentially expressed genes (DEGs) were identified. Gene Ontology (GO) enrichment analysis demonstrated that there were several genes involved in the response to different phytohormones, including abscisic acid (ABA), auxin (IAA), ethylene (ETH), gibberellins (GAs), and cytokinins (CTKs). Among them, multiple genes, such as PIN2 and ERF113, are involved in regulating vital plant movements by various phytohormone pathways. Moreover, following RR cultivation, DEGs were enriched in the biological processes of plant-type secondary cell wall biosynthesis, the defense response, programmed cell death involved in cell development, and the oxalate metabolic process. Furthermore, through a combined analysis of the transcriptome and previously published microRNA (miRNA) sequencing results, we found that multiple differentially expressed miRNAs (DEMs) and DEG combinations in different comparison groups exhibited opposite trends, indicating that the expression levels of miRNAs and their target genes were negatively correlated. Furthermore, RR treatment indeed significantly increased the ABA content at 125 days after planting and significantly decreased the IAA content at 70 days after planting. Under RR cultivation, most ABA biosynthesis-related genes were upregulated, while most IAA biosynthesis-related genes were downregulated. These findings lay a solid foundation for further establishing the network through which miRNAs regulate grapevine root development through target genes and for further exploring the molecular mechanism through which endogenous ABA and IAA regulate root architecture development in grapevine.