The Arabidopsis-mei2-like genes play a role in meiosis and vegetative growth in Arabidopsis

Uncovering genes involved in pollinator-driven mating system shifts and selfing syndrome evolution in Brassica rapa

Shifts in pollinator occurrence and their pollen transport effectiveness drive the evolution of mating systems in flowering plants. Understanding the genomic basis of these changes is essential for predicting the persistence of a species under environmental changes. We investigated the genomic changes in Brassica rapa over nine generations of pollination by hoverflies associated with rapid morphological evolution toward the selfing syndrome. We combined a genotyping-by-sequencing (GBS) approach with a genome-wide association study (GWAS) to identify candidate genes, and assessed their functional role in the observed morphological changes by studying mutations of orthologous genes in the model plant Arabidopsis thaliana. We found 31 candidate genes involved in a wide range of functions from DNA/RNA binding to transport. Our functional assessment of orthologous genes in A. thaliana revealed that two of the identified genes in B. rapa are involved in regulating the size of floral organs. We found a protein kinase superfamily protein involved in petal width, an important trait in plant attractiveness to pollinators. Moreover, we found a histone lysine methyltransferase (HKMT) associated with stamen length. Altogether, our study shows that hoverfly pollination leads to rapid evolution toward the selfing syndrome mediated by polygenic changes.

A co-located QTL for seven spike architecture-related traits shows promising breeding use potential in common wheat (Triticum aestivum L.)

KEY MESSAGE

A co-located novel QTL for TFS, FPs, FMs, FFS, FFPs, KWS, and KWPs with potential of improving wheat yield was identified and validated. Spike-related traits, including fertile florets per spike (FFS), kernel weight per spike (KWS), total florets per spike (TFS), florets per spikelet (FPs), florets in the middle spikelet (FMs), fertile florets per spikelet (FFPs), and kernel weight per spikelet (KWPs), are key traits in improving wheat yield. In the present study, quantitative trait loci (QTL) for these traits evaluated under various environments were detected in a recombinant inbred line population (msf/Chuannong 16) mainly genotyped using the 16 K SNP array. Ultimately, we identified 60 QTL, but only QFFS.sau-MC-1A for FFS was a major and stably expressed QTL. It was located on chromosome arm 1AS, where loci for TFS, FPs, FMs, FFS, FFPs, KWS, and KWPs were also simultaneously co-mapped. The effect of QFFS.sau-MC-1A was further validated in three independent segregating populations using a Kompetitive Allele-Specific PCR marker. For the co-located QTL, QFFS.sau-MC-1A, the presence of a positive allele from msf was associate with increases for all traits: + 12.29% TFS, + 10.15% FPs, + 13.97% FMs, + 17.12% FFS, + 14.75% FFPs, + 22.17% KWS, and + 19.42% KWPs. Furthermore, pleiotropy analysis showed that the positive allele at QFFS.sau-MC-1A simultaneously increased the spike length, spikelet number per spike, and thousand-kernel weight. QFFS.sau-MC-1A represents a novel QTL for marker-assisted selection with the potential for improving wheat yield. Four genes, TraesCS1A03G0012700, TraesCS1A03G0015700, TraesCS1A03G0016000, and TraesCS1A03G0016300, which may affect spike development, were predicted in the physical interval harboring QFFS.sau-MC-1A. Our results will help in further fine mapping QFFS.sau-MC-1A and be useful for improving wheat yield.

Heat stress during male meiosis impairs cytoskeletal organization, spindle assembly and tapetum degeneration in wheat

The significance of heat stress in agriculture is ever-increasing with the progress of global climate changes. Due to a negative effect on the yield of staple crops, including wheat, the impairment of plant reproductive development triggered by high ambient temperature became a restraint in food production. Although the heat sensitivity of male meiosis and the following gamete development in wheat has long been recognized, a detailed structural characterization combined with a comprehensive gene expression analysis has not been done about this phenomenon. We demonstrate here that heat stress severely alters the cytoskeletal configuration, triggers the failure of meiotic division in wheat. Moreover, it changes the expression of genes related to gamete development in male meiocytes and the tapetum layer in a genotype-dependent manner. 'Ellvis', a heat-tolerant winter wheat cultivar, showed high spikelet fertility rate and only scarce structural aberrations upon exposure to high temperature. In addition, heat shock genes and genes involved in scavenging reactive oxygen species were significantly upregulated in 'Ellvis', and the expression of meiosis-specific and major developmental genes showed high stability in this cultivar. In the heat-sensitive 'Mv 17-09', however, genes participating in cytoskeletal fiber nucleation, the spindle assembly checkpoint genes, and tapetum-specific developmental regulators were downregulated. These alterations may be related to the decreased cytoskeleton content, frequent micronuclei formation, and the erroneous persistence of the tapetum layer observed in the sensitive genotype. Our results suggest that understanding the heat-sensitive regulation of these gene functions would be an essential contribution to the development of new, heat-tolerant cultivars.

SPINDLY mediates O-fucosylation of hundreds of proteins and sugar-dependent growth in Arabidopsis

The recent discovery of SPINDLY (SPY)-catalyzed protein O-fucosylation revealed a novel mechanism for regulating nucleocytoplasmic protein functions in plants. Genetic evidence indicates the important roles of SPY in diverse developmental and physiological processes. However, the upstream signal controlling SPY activity and the downstream substrate proteins O-fucosylated by SPY remain largely unknown. Here, we demonstrated that SPY mediates sugar-dependent growth in Arabidopsis (Arabidopsis thaliana). We further identified hundreds of O-fucosylated proteins using lectin affinity chromatography followed by mass spectrometry. All the O-fucosylation events quantified in our proteomic analyses were undetectable or dramatically decreased in the spy mutants, and thus likely catalyzed by SPY. The O-fucosylome includes mostly nuclear and cytosolic proteins. Many O-fucosylated proteins function in essential cellular processes, phytohormone signaling, and developmental programs, consistent with the genetic functions of SPY. The O-fucosylome also includes many proteins modified by O-linked N-acetylglucosamine (O-GlcNAc) and by phosphorylation downstream of the target of rapamycin (TOR) kinase, revealing the convergence of these nutrient signaling pathways on key regulatory functions such as post-transcriptional/translational regulation and phytohormone responses. Our study identified numerous targets of SPY/O-fucosylation and potential nodes of crosstalk among sugar/nutrient signaling pathways, enabling future dissection of the signaling network that mediates sugar regulation of plant growth and development.

Unraveling Genomic Regions Controlling Root Traits as a Function of Nitrogen Availability in the MAGIC Wheat Population WM-800

An ever-growing world population demands to be fed in the future and environmental protection and climate change need to be taken into account. An important factor here is nitrogen uptake efficiency (NUpE), which is influenced by the root system (the interface between plant and soil). To understand the natural variation of root system architecture (RSA) as a function of nitrogen (N) availability, a subset of the multiparent advanced generation intercross (MAGIC) winter wheat population WM-800 was phenotyped under two contrasting N treatments in a high-throughput phenotyping system at the seedling stage. Fourteen root and shoot traits were measured. Subsequently, these traits were genetically analyzed using 13,060 polymorphic haplotypes and SNPs in a genome-wide association study (GWAS). In total, 64 quantitative trait loci (QTL) were detected; 60 of them were N treatment specific. Candidate genes for the detected QTL included NRT1.1 and genes involved in stress signaling under N-, whereas candidate genes under N+ were more associated with general growth, such as mei2 and TaWOX11b. This finding may indicate (i) a disparity of the genetic control of root development under low and high N supply and, furthermore, (ii) the need for an N specific selection of genes and genotypes in breeding new wheat cultivars with improved NUpE.

The social network of target of rapamycin complex 1 in plants

Target of rapamycin complex 1 (TORC1) is a highly conserved serine-threonine protein kinase crucial for coordinating growth according to nutrient availability in eukaryotes. It works as a central integrator of multiple nutrient inputs such as sugar, nitrogen, and phosphate and promotes growth and biomass accumulation in response to nutrient sufficiency. Studies, especially in the past decade, have identified the central role of TORC1 in regulating growth through interaction with hormones, photoreceptors, and stress signaling machinery in plants. In this review, we comprehensively analyse the interactome and phosphoproteome of the Arabidopsis TORC1 signaling network. Our analysis highlights the role of TORC1 as a central hub kinase communicating with the transcriptional and translational apparatus, ribosomes, chaperones, protein kinases, metabolic enzymes, and autophagy and stress response machinery to orchestrate growth in response to nutrient signals. This analysis also suggests that along with the conserved downstream components shared with other eukaryotic lineages, plant TORC1 signaling underwent several evolutionary innovations and co-opted many lineage-specific components during. Based on the protein-protein interaction and phosphoproteome data, we also discuss several uncharacterized and unexplored components of the TORC1 signaling network, highlighting potential links for future studies.

A grazing-driven positive nutrient feedback loop and active sexual reproduction underpin widespread Noctiluca green tides

The mixoplankton green Noctiluca scintillans (gNoctiluca) is known to form extensive green tides in tropical coastal ecosystems prone to eutrophication. In the Arabian Sea, their recent appearance and annual recurrence have upended an ecosystem that was once exclusively dominated by diatoms. Despite evidence of strong links to eutrophication, hypoxia and warming, the mechanisms underlying outbreaks of this mixoplanktonic dinoflagellate remain uncertain. Here we have used eco-physiological measurements and transcriptomic profiling to ascribe gNoctiluca's explosive growth during bloom formation to the form of sexual reproduction that produces numerous gametes. Rapid growth of gNoctiluca coincided with active ammonium and phosphate release from gNoctiluca cells, which exhibited high transcriptional activity of phagocytosis and metabolism generating ammonium. This grazing-driven nutrient flow ostensibly promotes the growth of phytoplankton as prey and offers positive support successively for bloom formation and maintenance. We also provide the first evidence that the host gNoctiluca cell could be manipulating growth of its endosymbiont population in order to exploit their photosynthetic products and meet critical energy needs. These findings illuminate gNoctiluca's little known nutritional and reproductive strategies that facilitate its ability to form intense and expansive gNoctiluca blooms to the detriment of regional water, food and the socio-economic security in several tropical countries.

Sporophytic control of pollen meiotic progression is mediated by tapetum expression of ABORTED MICROSPORES

Pollen development is dependent on the tapetum, a sporophytic anther cell layer surrounding the microspores that functions in pollen wall formation but is also essential for meiosis-associated development. There is clear evidence of crosstalk and co-regulation between the tapetum and microspores, but how this is achieved is currently not characterized. ABORTED MICROSPORES (AMS), a tapetum transcription factor, is important for pollen wall formation, but also has an undefined role in early pollen development. We conducted a detailed investigation of chromosome behaviour, cytokinesis, radial microtubule array (RMA) organization, and callose formation in the ams mutant. Early meiosis initiates normally in ams, shows delayed progression after the pachytene stage, and then fails during late meiosis, with disorganized RMA, defective cytokinesis, abnormal callose formation, and microspore degeneration, alongside abnormal tapetum development. Here, we show that selected meiosis-associated genes are directly repressed by AMS, and that AMS is essential for late meiosis progression. Our findings indicate that AMS has a dual function in tapetum-meiocyte crosstalk by playing an important regulatory role during late meiosis, in addition to its previously characterized role in pollen wall formation. AMS is critical for RMA organization, callose deposition, and therefore cytokinesis, and is involved in the crosstalk between the gametophyte and sporophytic tissues, which enables synchronous development of tapetum and microspores.

Structural insights reveal the specific recognition of meiRNA by the Mei2 protein

In the fission yeast Schizosaccharomyces pombe, Mei2, an RNA-binding protein essential for entry into meiosis, regulates meiosis initiation. Mei2 binds to a specific non-coding RNA species, meiRNA, and accumulates at sme2 gene locus, which encodes meiRNA. Previous research has shown that the Mei2 C-terminal RNA recognition motif (RRM3) physically interacts with meiRNA 5' region in vitro and stimulates meiosis in vivo. However, the underlying mechanism still remains elusive. We first employed an in vitro crosslinking and immunoprecipitation sequencing (CLIP-seq) assay and demonstrated a preference for U-rich motifs of meiRNA by Mei2 RRM3. We then solved the crystal structures of Mei2 RRM3 in the apo form and complex with an 8mer RNA fragment, derived from meiRNA, as detected by in vitro CLIP-seq. These results provide structural insights into Mei2 RRM3-meiRNA complex and reveal that Mei2 RRM3 binds specifically to the UUC(U) sequence. Furthermore, a structure-based Mei2 mutation, Mei2F644A causes defective karyogamy, suggesting an essential role of the RNA-binding ability of Mei2 in regulating meiosis.

Linking genome wide RNA sequencing with physio-biochemical and cytological responses to catalogue key genes and metabolic pathways for alkalinity stress tolerance in lentil (Lens culinaris Medikus)

BACKGROUND

Alkaline soils cause low productivity in crop plants including lentil. Alkalinity adaptation strategies in lentil were revealed when morpho-anatomical and physio-biochemical observations were correlated with transcriptomics analysis in tolerant (PDL-1) and sensitive (L-4076) cultivars at seedling stage.

RESULTS

PDL-1 had lesser salt injury and performed better as compared to L-4076. Latter showed severe wilting symptoms and higher accumulation of Na⁺ and lower K⁺ in roots and shoots. PDL-1 performed better under high alkalinity stress which can be attributed to its higher mitotic index, more accumulation of K⁺ in roots and shoots and less aberrantly dividing cells. Also, antioxidant enzyme activities, osmolytes' accumulation, relative water content, membrane stability index and abscisic acid were higher in this cultivar. Differentially expressed genes (DEGs) related to these parameters were upregulated in tolerant genotypes compared to the sensitive one. Significantly up-regulated DEGs were found to be involved in abscisic acid (ABA) signalling and secondary metabolites synthesis. ABA responsive genes viz. dehydrin 1, 9-cis-epoxycarotenoid dioxygenase, ABA-responsive protein 18 and BEL1-like homeodomain protein 1 had log₂fold change above 4.0. A total of 12,836 simple sequence repeats and 4,438 single nucleotide polymorphisms were identified which can be utilized in molecular studies.

CONCLUSIONS

Phyto-hormones biosynthesis-predominantly through ABA signalling, and secondary metabolism are the most potent pathways for alkalinity stress tolerance in lentil. Cultivar PDL-1 exhibited high tolerance towards alkalinity stress and can be used in breeding programmes for improving lentil production under alkalinity stress conditions.

Histone demethylase IBM1-mediated meiocyte gene expression ensures meiotic chromosome synapsis and recombination

Histone methylation and demethylation play important roles in plant growth and development, but the involvement of histone demethylation during meiosis is poorly understood. Here we show that disruption of Arabidopsis thaliana INCREASE IN BONSAI METHYLATION 1 (IBM1) causes incomplete synapsis, chromosome entanglement and reduction of recombination during meiosis, leading to sterility. Interestingly, these ibm1 meiotic defects are rescued by mutations in either SUVH4/KYP or CMT3. Using transcriptomic analyses we show that mutation of IBM1 down-regulates thousands of genes expressed in meiocytes, and that expression of about 38% of these genes are restored to wild type levels in ibm1 cmt3 double mutants. Changes in the expression of 437 of these, including the ARABIDOPSIS MEI2-LIKE AML3-5 genes, are correlated with a significant reduction of gene body CHG methylation. Consistently, the aml3 aml4 aml5 triple have defects in synapsis and chromosome entanglement similar to ibm1. Genetic analysis shows that aml3 aml4 aml5 ibm1 quadruple mutants resembles the ibm1 single mutant. Strikingly, over expression of AML5 in ibm1 can partially rescue the ibm1 meiotic defects. Taken together, our results demonstrate that histone demethylase IBM1 is required for meiosis likely via coordinated regulation of meiocyte gene expression during meiosis.

The Arabidopsis PeptideAtlas: Harnessing worldwide proteomics data to create a comprehensive community proteomics resource

We developed a resource, the Arabidopsis PeptideAtlas (www.peptideatlas.org/builds/arabidopsis/), to solve central questions about the Arabidopsis thaliana proteome, such as the significance of protein splice forms and post-translational modifications (PTMs), or simply to obtain reliable information about specific proteins. PeptideAtlas is based on published mass spectrometry (MS) data collected through ProteomeXchange and reanalyzed through a uniform processing and metadata annotation pipeline. All matched MS-derived peptide data are linked to spectral, technical, and biological metadata. Nearly 40 million out of ∼143 million MS/MS (tandem MS) spectra were matched to the reference genome Araport11, identifying ∼0.5 million unique peptides and 17,858 uniquely identified proteins (only isoform per gene) at the highest confidence level (false discovery rate 0.0004; 2 non-nested peptides ≥9 amino acid each), assigned canonical proteins, and 3,543 lower-confidence proteins. Physicochemical protein properties were evaluated for targeted identification of unobserved proteins. Additional proteins and isoforms currently not in Araport11 were identified that were generated from pseudogenes, alternative start, stops, and/or splice variants, and small Open Reading Frames; these features should be considered when updating the Arabidopsis genome. Phosphorylation can be inspected through a sophisticated PTM viewer. PeptideAtlas is integrated with community resources including TAIR, tracks in JBrowse, PPDB, and UniProtKB. Subsequent PeptideAtlas builds will incorporate millions more MS/MS data.

Transcriptome and metabolome analyses of cold and darkness-induced pellicle cysts of Scrippsiella trochoidea

BACKGROUND

Dinoflagellates are a group of unicellular organisms that are a major component of aquatic eukaryotes and important contributors to marine primary production. Nevertheless, many dinoflagellates are considered harmful algal bloom (HAB) species due to their detrimental environmental and human health impacts. Cyst formation is widely perceived as an adaptive strategy of cyst-forming dinoflagellates in response to adverse environmental conditions. Dinoflagellate cysts play critical roles in bloom dynamics. However, our insight into the underlying molecular basis of encystment is still limited. To investigate the molecular processes regulating encystment in dinoflagellates, transcriptome and metabolome investigations were performed on cold and darkness-induced pellicle cysts of Scrippsiella trochoidea.

RESULTS

No significant transcriptional response was observed at 2 h; however, massive transcriptome and metabolome reprogramming occurred at 5 h and in pellicle cysts. The gene-to-metabolite network demonstrated that the initial transformation from vegetative cells into pellicle cysts was highly energy demanding through the activation of catabolism, including glycolysis, β-oxidation, TCA cycle and oxidative phosphorylation, to cope with cold-darkness-induced stress. However, after transformation into pellicle cysts, the metabolism was greatly reduced, and various sugars, polyunsaturated fatty acids and amino acids accumulated to prolong survival. The identification of 56 differentially expressed genes (DEGs) related to signal transduction indicated that S. trochoidea received a cold-darkness signal that activated multiple signal transduction pathways, leading to encystment. The elevated expression of genes encoding enzymes involved in ROS stress suggested that pellicle cysts respond to increased oxidative stress. Several cell cycle-related genes were repressed. Intriguingly, 11 DEGs associated with sexual reproduction suggested that pellicle cysts (or some portion thereof) may be a product of sexual reproduction.

CONCLUSIONS

This study provides the first transcriptome and metabolome analyses conducted during the encystment of S. trochoidea, an event that requires complex regulatory mechanisms and impacts on population dynamics. The results reveal comprehensive molecular regulatory processes underlying life cycle regulation in dinoflagellates involving signal transduction, gene expression and metabolite profile, which will improve our ability to understand and monitor dinoflagellate blooms.

Untangling structural factors driving genome stabilization in nascent Brassica napus allopolyploids

Allopolyploids have globally higher fitness than their diploid progenitors; however, by comparison, most resynthesized allopolyploids have poor fertility and highly unstable genome. Elucidating the evolutionary processes promoting genome stabilization and fertility is thus essential to comprehend allopolyploid success. Using the Brassica model, we mimicked the speciation process of a nascent allopolyploid species by resynthesizing allotetraploid Brassica napus and systematically selecting for euploid individuals over eight generations in four independent allopolyploidization events with contrasted genetic backgrounds, cytoplasmic donors, and polyploid formation type. We evaluated the evolution of meiotic behavior and fertility and identified rearrangements in S₁ to S₉ lineages to explore the positive consequences of euploid selection on B. napus genome stability. Recurrent selection of euploid plants for eight generations drastically reduced the percentage of aneuploid progenies as early as the fourth generation, concomitantly with a decrease in number of newly fixed homoeologous rearrangements. The consequences of homoeologous rearrangements on meiotic behavior and seed number depended strongly on the genetic background and cytoplasm donor. The combined use of both self-fertilization and recurrent euploid selection allowed identification of genomic regions associated with fertility and meiotic behavior, providing complementary evidence to explain B. napus speciation success.

A major quantitative trait locus on chromosome A9, BnaPh1, controls homoeologous recombination in Brassica napus

Ensuring faithful homologous recombination in allopolyploids is essential to maintain optimal fertility of the species. Variation in the ability to control aberrant pairing between homoeologous chromosomes in Brassica napus has been identified. The current study exploited the extremes of such variation to identify genetic factors that differentiate newly resynthesised B. napus, which is inherently unstable, and established B. napus, which has adapted to largely control homoeologous recombination. A segregating B. napus mapping population was analysed utilising both cytogenetic observations and high-throughput genotyping to quantify the levels of homoeologous recombination. Three quantitative trait loci (QTL) were identified that contributed to the control of homoeologous recombination in the important oilseed crop B. napus. One major QTL on BnaA9 contributed between 32 and 58% of the observed variation. This study is the first to assess homoeologous recombination and map associated QTLs resulting from deviations in normal pairing in allotetraploid B. napus. The identified QTL regions suggest candidate meiotic genes that could be manipulated in order to control this important trait and further allow the development of molecular markers to utilise this trait to exploit homoeologous recombination in a crop.

A replication-competent late liver stage-attenuated human malaria parasite

Whole-sporozoite vaccines engender sterilizing immunity against malaria in animal models and importantly, in humans. Gene editing allows for the removal of specific parasite genes, enabling generation of genetically attenuated parasite (GAP) strains for vaccination. Using rodent malaria parasites, we have previously shown that late liver stage-arresting replication-competent (LARC) GAPs confer superior protection when compared with early liver stage-arresting replication-deficient GAPs and radiation-attenuated sporozoites. However, generating a LARC GAP in the human malaria parasite Plasmodium falciparum (P. falciparum) has been challenging. Here, we report the generation and characterization of a likely unprecedented P. falciparum LARC GAP generated by targeted gene deletion of the Mei2 gene: P. falciparum mei2-. Robust exoerythrocytic schizogony with extensive cell growth and DNA replication was observed for P. falciparum mei2- liver stages in human liver-chimeric mice. However, P. falciparum mei2- liver stages failed to complete development and did not form infectious exoerythrocytic merozoites, thereby preventing their transition to asexual blood stage infection. Therefore, P. falciparum mei2- is a replication-competent, attenuated human malaria parasite strain with potentially increased potency, useful for vaccination to protect against P. falciparum malaria infection.

Control of Grain Size and Weight by the GSK2-LARGE1/OML4 Pathway in Rice

Regulation of grain size is crucial for improving crop yield and is also a basic aspect in developmental biology. However, the genetic and molecular mechanisms underlying grain size control in crops remain largely unknown despite their central importance. Here, we report that the MEI2-LIKE PROTEIN4 (OML4) encoded by the LARGE1 gene is phosphorylated by GLYCOGEN SYNTHASE KINASE2 (GSK2) and negatively controls grain size and weight in rice (Oryza sativa). Loss of function of OML4 leads to large and heavy grains, while overexpression of OML4 causes small and light grains. OML4 regulates grain size by restricting cell expansion in the spikelet hull. OML4 is expressed in developing panicles and grains, and the GFP-OML4 fusion protein is localized in the nuclei. Biochemical analyses show that the GSK2 physically interacts with OML4 and phosphorylates it, thereby possibly influencing the stability of OML4. Genetic analyses support that GSK2 and OML4 act, at least in part, in a common pathway to control grain size in rice. These results reveal the genetic and molecular mechanism of a GSK2-OML4 regulatory module in grain size control, suggesting that this pathway is a suitable target for improving seed size and weight in crops.

Differential gene expression profiling of one- and two-dimensional apogamous gametophytes of the fern Dryopteris affinis ssp. affinis

Apomixis was originally defined as the replacement of sexual reproduction by an asexual process that does not involve fertilization but, in angiosperms, it is often used in the more restricted sense of asexual reproduction through seeds. In ferns, apomixis combines the production of unreduced spores (diplospory) and the formation of sporophytes from somatic cells of the prothallium (apogamy). The genes that control the onset of apogamy in ferns are largely unknown. In this study, we describe the gametophyte transcriptome of the apogamous fern Dryopteris affinis ssp. affinis using an RNA-Seq approach to compare the gene expression profiles of one- and two-dimensional gametophytes, the latter containing apogamic centers. After collapsing highly similar de novo transcripts, we obtained 166,191 unigenes, of which 30% could be annotated using public databases. Multiple quality metrics indicate a good quality of the de novo transcriptome with a low level of fragmentation. Our data show a total of 10,679 genes (6% of all genes) to be differentially expressed between gametophytes of filamentous (one-dimensional) and prothallial (two-dimensional) architecture. 6,110 genes were up-regulated in two-dimensional relative to one-dimensional gametophytes, some of which are implicated in the regulation of meristem growth, auxin signaling, reproduction, and sucrose metabolism. 4,570 genes were down-regulated in two-dimensional versus one-dimensional gametophytes, which are enriched in stimulus and defense genes, as well as genes involved in epigenetic gene regulation and ubiquitin degradation. Our results provide insights into free-living gametophyte development, focusing on the filamentous-to-prothallus growth transition, and provide a useful resource for further investigations of asexual reproduction.

Identification and expression analysis of meiosis-related genes in the harmful alga Heterosigma akashiwo (Raphidophyceae)

Sexual reproduction (SR), hallmarked by meiosis, is widespread in eukaryotes. In phytoplankton, SR has been observed in many lineages, but molecular information on SR or meiosis of harmful algal bloom (HAB) species is scarce. The raphidophyte Heterosigma akashiwo is a globally distributed and devastating HAB species, but molecular evidence of its SR or meiosis is lacking. Here, to address the gap of knowledge, the presence of meiotic genes in H. akashiwo were examined. Interestingly, seven meiosis-specific or related genes (SPO11, MND, RAD21, RAD51, MSH2, MSH6 and MEI2) were identified from H. akashiwo transcriptomes. Furthermore, expression patterns of these genes except MSH6 (excluded due to primer failure) were investigated using quantitative reverse-transcription PCR. Results showed that the examined genes exhibited similar diel rhythms, typically, highest in early dark period and then gradually decreasing until mid-day. In addition, the expression of these six genes was not higher in the stationary growth stage than in the exponential stage, as would be expected if meiosis was to form cysts, and their elevated expression in response to colchicine treatment (arresting cells in the G2/M transition) indicated a potential role of these genes in cell division and population growth in H. akashiwo. Consistent with this, we also found a morning to afternoon shift in the expression of these genes during the bloom of H. akashiwo. This study documents a part of the typical SR gene repertoire and its potential role in regulating cell division in H. akashiwo, offering candidates for population growth markers for bloom monitoring although its linkage to meiosis and SR remain to be studied further in the future.

Genome Improvement and Core Gene Set Refinement of Fugacium kawagutii

Cataloging an accurate functional gene set for the Symbiodiniaceae species is crucial for addressing biological questions of dinoflagellate symbiosis with corals and other invertebrates. To improve the gene models of Fugacium kawagutii, we conducted high-throughput chromosome conformation capture (Hi-C) for the genome and Illumina combined with PacBio sequencing for the transcriptome to achieve a new genome assembly and gene prediction. A 0.937-Gbp assembly of F. kawagutii were obtained, with a N50 > 13 Mbp and the longest scaffold of 121 Mbp capped with telomere motif at both ends. Gene annotation produced 45,192 protein-coding genes, among which, 11,984 are new compared to previous versions of the genome. The newly identified genes are mainly enriched in 38 KEGG pathways including N-Glycan biosynthesis, mRNA surveillance pathway, cell cycle, autophagy, mitophagy, and fatty acid synthesis, which are important for symbiosis, nutrition, and reproduction. The newly identified genes also included those encoding O-methyltransferase (O-MT), 3-dehydroquinate synthase, homologous-pairing protein 2-like (HOP2) and meiosis protein 2 (MEI2), which function in mycosporine-like amino acids (MAAs) biosynthesis and sexual reproduction, respectively. The improved version of the gene set (Fugka_Geneset _V3) raised transcriptomic read mapping rate from 33% to 54% and BUSCO match from 29% to 55%. Further differential gene expression analysis yielded a set of stably expressed genes under variable trace metal conditions, of which 115 with annotated functions have recently been found to be stably expressed under three other conditions, thus further developing the "core gene set" of F. kawagutii. This improved genome will prove useful for future Symbiodiniaceae transcriptomic, gene structure, and gene expression studies, and the refined "core gene set" will be a valuable resource from which to develop reference genes for gene expression studies.

Asexual thalli originated from sporophytic thalli via apomeiosis in the green seaweed Ulva

Apomixis is an asexual reproduction system without fertilization, which is an important proliferation strategy for plants and algae. Here, we report on the apomeiosis in the green seaweed Ulva prolifera, which has sexual and obligate asexual populations. Genomic PCR of mating type (MT)-locus genes revealed asexual thalli carrying both MT genomes. Observation of the chromosomes during the formation of each type of reproductive cell revealed that cells in asexual thalli performed apomeiosis without chromosome reduction. Moreover, genotyping revealed that laboratory-cultured sporophytic thalli produced not only each type of gametophyte but also diploid thalli carrying the mt^- and mt⁺ genome (mt^± thallus strains). The mt^± thallus strain released diploid biflagellate zoids, with ultrastructure and behavior similar to mt⁺ gametes. Additionally, a transcriptomic analysis revealed that some meiosis-related genes (Mei2L and RAD1) were highly expressed in the quadriflagellate zoosporoids. Our results strongly suggest that asexual thalli originally evolved via apomeiosis in sporophytic thalli.

Studies of a rice sterile mutant sstl from the TRIM collection

BACKGROUND

Rice (Oryza sativa) is one of the main crops in the world, and more than 3.9 billion people will consume rice by 2025. Sterility significantly affects rice production and leads to yield defects. The undeveloped anthers or abnormal pollen represent serious defects in rice male sterility. Therefore, understanding the mechanism of male sterility is an important task. Here, we investigated a rice sterile mutant according to its developmental morphology and transcriptional profiles.

RESULTS

An untagged T-DNA insertional mutant showed defective pollen and abnormal anthers as compared with its semi-sterile mutant (sstl) progeny segregates. Transcriptomic analysis of sterile sstl-s revealed several biosynthesis pathways, such as downregulated cell wall, lipids, secondary metabolism, and starch synthesis. This downregulation is consistent with the morphological characterization of sstl-s anthers with irregular exine, absence of intine, no starch accumulation in pollen grains and no accumulated flavonoids in anthers. Moreover, defective microsporangia development led to abnormal anther locule and aborted microspores. The downregulated lipids, starch, and cell wall synthesis-related genes resulted in loss of fertility.

CONCLUSIONS

We illustrate the importance of microsporangia in the development of anthers and functional microspores. Abnormal development of pollen grains, pollen wall, anther locule, etc. result in severe yield reduction.

Genetic basis for coordination of meiosis and sexual structure maturation in Cryptococcus neoformans

In the human fungal pathogen Cryptococcus neoformans, sex can benefit its pathogenicity through production of meiospores, which are believed to offer both physical and meiosis-created lineage advantages for its infections. Cryptococcus sporulation occurs following two parallel events, meiosis and differentiation of the basidium, the characteristic sexual structure of the basidiomycetes. However, the circuit integrating these events to ensure subsequent sporulation is unclear. Here, we show the spatiotemporal coordination of meiosis and basidial maturation by visualizing event-specific molecules in developing basidia defined by a quantitative approach. Monitoring of gene induction timing together with genetic analysis reveals co-regulation of the coordinated events by a shared regulatory program. Two RRM family regulators, Csa1 and Csa2, are crucial components that bridge meiosis and basidial maturation, further determining sporulation. We propose that the regulatory coordination of meiosis and basidial development serves as a determinant underlying the production of infectious meiospores in C. neoformans.

Transcriptomic Analyses of Scrippsiella trochoidea Reveals Processes Regulating Encystment and Dormancy in the Life Cycle of a Dinoflagellate, with a Particular Attention to the Role of Abscisic Acid

Due to the vital importance of resting cysts in the biology and ecology of many dinoflagellates, a transcriptomic investigation on Scrippsiella trochoidea was conducted with the aim to reveal the molecular processes and relevant functional genes regulating encystment and dormancy in dinoflagellates. We identified via RNA-seq 3,874 (out of 166,575) differentially expressed genes (DEGs) between resting cysts and vegetative cells; a pause of photosynthesis (confirmed via direct measurement of photosynthetic efficiency); an active catabolism including β-oxidation, glycolysis, glyoxylate pathway, and TCA in resting cysts (tested via measurements of respiration rate); 12 DEGs encoding meiotic recombination proteins and members of MEI2-like family potentially involved in sexual reproduction and encystment; elevated expressions in genes encoding enzymes responding to pathogens (chitin deacetylase) and ROS stress in cysts; and 134 unigenes specifically expressed in cysts. We paid particular attention to genes pertaining to phytohormone signaling and identified 4 key genes regulating abscisic acid (ABA) biosynthesis and catabolism, with further characterization based on their full-length cDNA obtained via RACE-PCR. The qPCR results demonstrated elevated biosynthesis and repressed catabolism of ABA during the courses of encystment and cyst dormancy, which was significantly enhanced by lower temperature (4 ± 1°C) and darkness. Direct measurements of ABA using UHPLC-MS/MS and ELISA in vegetative cells and cysts both fully supported qPCR results. These results collectively suggest a vital role of ABA in regulating encystment and maintenance of dormancy, akin to its function in seed dormancy of higher plants. Our results provided a critical advancement in understanding molecular processes in resting cysts of dinoflagellates.

The TRAF Mediated Gametogenesis Progression (TRAMGaP) Gene Is Required for Megaspore Mother Cell Specification and Gametophyte Development

In plants, the role of TRAF-like proteins with meprin and the TRAF homology (MATH) domain is far from clear. In animals, these proteins serve as adapter molecules to mediate signal transduction from Tumor Necrosis Factor Receptor to downstream effector molecules. A seed-sterile mutant with a disrupted TRAF-like gene (At5g26290) exhibiting aberrant gametogenesis led us to investigate the developmental role of this gene in Arabidopsis (Arabidopsis thaliana). The mutation was semidominant and resulted in pleiotropic phenotypes with such features as short siliques with fewer ovules, pollen and seed sterility, altered Megaspore Mother Cell (MMC) specification, and delayed programmed cell death in megaspores and the tapetum, features that overlapped those in other well-characterized mutants. Seed sterility and reduced transmission frequency of the mutant alleles pointed to a dual role, sporophytic and gametophytic, for the gene on the male side. The mutant also showed altered expression of various genes involved in such cellular and developmental pathways as regulation of transcription, biosynthesis and transport of lipids, hormone-mediated signaling, and gametophyte development. The diverse phenotypes of the mutant and the altered expression of key genes related to gametophyte and seed development could be explained based on the functional similarly between At5g26290 and MATH-BTB domain proteins that modulate gene expression through the ubiquitin-mediated proteasome system. These results show a novel link between a TRAF-like gene and reproductive development in plants.

De novo transcriptome of the cosmopolitan dinoflagellate Amphidinium carterae to identify enzymes with biotechnological potential

Dinoflagellates are phytoplanktonic organisms found in both freshwater and marine habitats. They are often studied because related to harmful algal blooms but they are also known to produce bioactive compounds for the treatment of human pathologies. The aim of this study was to sequence the full transcriptome of the dinoflagellate Amphidinium carterae in both nitrogen-starved and -replete culturing conditions (1) to evaluate the response to nitrogen starvation at the transcriptional level, (2) to look for possible polyketide synthases (PKSs) in the studied clone (genes that may be involved in the synthesis of bioactive compounds), (3) if present, to evaluate if nutrient starvation can influence PKS expression, (4) to look for other possible enzymes of biotechnological interest and (5) to test strain cytotoxicity on human cell lines. Results showed an increase in nitrogen metabolism and stress response in nitrogen-starved cells and confirmed the presence of a type I β-ketosynthase. In addition, L-asparaginase (used for the treatment of Leukemia and for acrylamide reduction in food industries) and cellulase (useful for biofuel production and other industrial applications) have been identified for the first time in this species, giving new insights into possible biotechnological applications of dinoflagellates.

TORC1-Dependent Phosphorylation Targets in Fission Yeast

Target of rapamycin (TOR) kinase controls cell metabolism and growth in response to environmental cues such as nutrients, growth factors, and stress. TOR kinase is widely conserved across eukaryotes. As in other organisms, the fission yeast Schizosaccharomyces pombe has two types of TOR complex, namely TOR complex 1 (TORC1) and TORC2. It is interesting that the two TOR complexes in S. pombe have opposite roles in sexual differentiation, which is induced by nutrient starvation. TORC1, which contains Tor2 as a catalytic subunit, promotes vegetative growth and represses sexual differentiation in nutrient-rich conditions, while TORC2 is required for the initiation of sexual differentiation. Multiple targets of TORC1 have been identified. Some of these, such as S6 kinase and an autophagy regulator Atg13, are known targets in other organisms. In addition, there is a novel group of TORC1 targets involved in the regulation of sexual differentiation. Here, we review recent findings on phosphorylation targets of TORC1 in S. pombe. Furthermore, we briefly report a novel S. pombe target of TORC1.

Transcriptomic Analysis of Thermally Stressed Symbiodinium Reveals Differential Expression of Stress and Metabolism Genes

Endosymbioses between dinoflagellate algae (Symbiodinium sp.) and scleractinian coral species form the foundation of coral reef ecosystems. The coral symbiosis is highly susceptible to elevated temperatures, resulting in coral bleaching, where the algal symbiont is released from host cells. This experiment aimed to determine the transcriptional changes in cultured Symbiodinium, to better understand the response of cellular mechanisms under future temperature conditions. Cultures were exposed to elevated temperatures (average 31°C) or control conditions (24.5°C) for a period of 28 days. Whole transcriptome sequencing of Symbiodinium cells on days 4, 19, and 28 were used to identify differentially expressed genes under thermal stress. A large number of genes representing 37.01% of the transcriptome (∼23,654 unique genes, FDR < 0.05) with differential expression were detected at no less than one of the time points. Consistent with previous studies of Symbiodinium gene expression, fold changes across the transcriptome were low, with 92.49% differentially expressed genes at ≤2-fold change. The transcriptional response included differential expression of genes encoding stress response components such as the antioxidant network and molecular chaperones, cellular components such as core photosynthesis machinery, integral light-harvesting protein complexes and enzymes such as fatty acid desaturases. Differential expression of genes encoding glyoxylate cycle enzymes were also found, representing the first report of this in Symbiodinium. As photosynthate transfer from Symbiodinium to coral hosts provides up to 90% of a coral's daily energy requirements, the implications of altered metabolic processes from exposure to thermal stress found in this study on coral-Symbiodinium associations are unknown and should be considered when assessing the stability of the symbiotic relationship under future climate conditions.

ANGUSTIFOLIA, a Plant Homolog of CtBP/BARS Localizes to Stress Granules and Regulates Their Formation

The ANGUSTIFOLIA (AN) gene in Arabidopsis is important for a plethora of morphological phenotypes. Recently, AN was also reported to be involved in responses to biotic and abiotic stresses. It encodes a homolog of the animal C-terminal binding proteins (CtBPs). In contrast to animal CtBPs, AN does not appear to function as a transcriptional co-repressor and instead functions outside nucleus where it might be involved in Golgi-associated membrane trafficking. In this study, we report a novel and unexplored role of AN as a component of stress granules (SGs). Interaction studies identified several RNA binding proteins that are associated with AN. AN co-localizes with several messenger ribonucleoprotein granule markers to SGs in a stress dependent manner. an mutants exhibit an altered SG formation. We provide evidence that the NAD(H) binding domain of AN is relevant in this context as proteins carrying mutations in this domain localize to a much higher degree to SGs and strongly reduce AN dimerization and its interaction with one interactor but not the others. Finally, we show that AN is a negative regulator of salt and osmotic stress responses in Arabidopsis suggesting a functional relevance in SGs.

An RRM-containing mei2-like MCT1 plays a negative role in the seed germination and seedling growth of Arabidopsis thaliana in the presence of ABA

Despite an increasing understanding of the essential role of the Mei2 gene encoding an RNA-binding protein (RBP) in premeiotic DNA synthesis and meiosis in yeasts and animals, the functional roles of the mei2-like genes in plant growth and development are largely unknown. Contrary to other mei2-like RBPs that contain three RNA-recognition motifs (RRMs), the mei2 C-terminal RRM only (MCT) is unique in that it harbors only the last C-terminal RRM. Although MCTs have been implicated to play important roles in plants, their functional roles in stress responses as well as plant growth and development are still unknown. Here, we investigated the expression and functional role of MCT1 (At1g37140) in plant response to abscisic acid (ABA). Confocal analysis of MCT1-GFP-expressing plants revealed that MCT1 is localized to the nucleus. The transcript level of MCT1 was markedly increased upon ABA treatment. Analysis of MCT1-overexpressing transgenic Arabidopsis plants and artificial miRNA-mediated mct1 knockdown mutants demonstrated that MCT1 inhibited seed germination and cotyledon greening of Arabidopsis plants under ABA. The transcript levels of ABA signaling-related genes, such as ABI3, ABI4, and ABI5, were markedly increased in the MCT1-overexpressing transgenic plant. Collectively, these results suggest that ABA-upregulated MCT1 plays a negative role in Arabidopsis seed germination and seedling growth under ABA by modulating the expression of ABA signaling-related genes.

Abundant protein phosphorylation potentially regulates Arabidopsis anther development

As the male reproductive organ of flowering plants, the stamen consists of the anther and filament. Previous studies on stamen development mainly focused on single gene functions by genetic methods or gene expression changes using comparative transcriptomic approaches, especially in model plants such as Arabidopsis thaliana However, studies on Arabidopsis anther protein expression and post-translational modifications are still lacking. Here we report proteomic and phosphoproteomic studies on developing Arabidopsis anthers at stages 4-7 and 8-12. We identified 3908 high-confidence phosphorylation sites corresponding to 1637 phosphoproteins. Among the 1637 phosphoproteins, 493 were newly identified, with 952 phosphorylation sites. Phosphopeptide enrichment prior to LC-MS analysis facilitated the identification of low-abundance proteins and regulatory proteins, thereby increasing the coverage of proteomic analysis, and facilitated the analysis of more regulatory proteins. Thirty-nine serine and six threonine phosphorylation motifs were uncovered from the anther phosphoproteome and further analysis supports that phosphorylation of casein kinase II, mitogen-activated protein kinases, and 14-3-3 proteins is a key regulatory mechanism in anther development. Phosphorylated residues were preferentially located in variable protein regions among family members, but they were they were conserved across angiosperms in general. Moreover, phosphorylation might reduce activity of reactive oxygen species scavenging enzymes and hamper brassinosteroid signaling in early anther development. Most of the novel phosphoproteins showed tissue-specific expression in the anther according to previous microarray data. This study provides a community resource with information on the abundance and phosphorylation status of thousands of proteins in developing anthers, contributing to understanding post-translational regulatory mechanisms during anther development.

Regulation of Small RNAs and Corresponding Targets in Nod Factor-Induced Phaseolus vulgaris Root Hair Cells

A genome-wide analysis identified the set of small RNAs (sRNAs) from the agronomical important legume Phaseolus vulgaris (common bean), including novel P. vulgaris-specific microRNAs (miRNAs) potentially important for the regulation of the rhizobia-symbiotic process. Generally, novel miRNAs are difficult to identify and study because they are very lowly expressed in a tissue- or cell-specific manner. In this work, we aimed to analyze sRNAs from common bean root hairs (RH), a single-cell model, induced with pure Rhizobium etli nodulation factors (NF), a unique type of signal molecule. The sequence analysis of samples from NF-induced and control libraries led to the identity of 132 mature miRNAs, including 63 novel miRNAs and 1984 phasiRNAs. From these, six miRNAs were significantly differentially expressed during NF induction, including one novel miRNA: miR-RH82. A parallel degradome analysis of the same samples revealed 29 targets potentially cleaved by novel miRNAs specifically in NF-induced RH samples; however, these novel miRNAs were not differentially accumulated in this tissue. This study reveals Phaseolus vulgaris-specific novel miRNA candidates and their corresponding targets that meet all criteria to be involved in the regulation of the early nodulation events, thus setting the basis for exploring miRNA-mediated improvement of the common bean–rhizobia symbiosis.

Chloroplast RNA-Binding Protein RBD1 Promotes Chilling Tolerance through 23S rRNA Processing in Arabidopsis

Plants have varying abilities to tolerate chilling (low but not freezing temperatures), and it is largely unknown how plants such as Arabidopsis thaliana achieve chilling tolerance. Here, we describe a genome-wide screen for genes important for chilling tolerance by their putative knockout mutants in Arabidopsis thaliana. Out of 11,000 T-DNA insertion mutant lines representing half of the genome, 54 lines associated with disruption of 49 genes had a drastic chilling sensitive phenotype. Sixteen of these genes encode proteins with chloroplast localization, suggesting a critical role of chloroplast function in chilling tolerance. Study of one of these proteins RBD1 with an RNA binding domain further reveals the importance of chloroplast translation in chilling tolerance. RBD1 is expressed in the green tissues and is localized in the chloroplast nucleoid. It binds directly to 23S rRNA and the binding is stronger under chilling than at normal growth temperatures. The rbd1 mutants are defective in generating mature 23S rRNAs and deficient in chloroplast protein synthesis especially under chilling conditions. Together, our study identifies RBD1 as a regulator of 23S rRNA processing and reveals the importance of chloroplast function especially protein translation in chilling tolerance.

Compared to cold acclimation (enhancement of freezing tolerance by a prior exposure to low non-freezing temperature), the tolerance mechanism to non-freezing chilling temperatures is not well understood. Here, we performed a genome-wide mutant screen for chilling sensitive phenotype and identified 49 candidate genes important for chilling tolerance in Arabidopsis. Among the proteins encoded by these 49 genes, 16 are annotated as having chloroplast localization, suggesting a critical role of chloroplast function in chilling tolerance. We further studied RBD1, one of the four RNA-binding proteins localized to chloroplast. RBD1 is only expressed in the green photosynthetic tissues and is localized to nucleoid of chloroplasts. Furthermore, RBD1 is found to be a regulator of 23S rRNA processing likely through direct binding to the precursor of 23S rRNA in a temperature dependent manner. Our study thus reveals the importance of chloroplast function especially protein translation in chilling tolerance at genome-wide scale and suggests an adaptive mechanism involving low temperature enhanced activities from proteins such as RBD1 in chilling tolerance.

TOR signalling in plants

Although the eukaryotic TOR (target of rapamycin) kinase signalling pathway has emerged as a key player for integrating nutrient-, energy- and stress-related cues with growth and metabolic outputs, relatively little is known of how this ancient regulatory mechanism has been adapted in higher plants. Drawing comparisons with the substantial knowledge base around TOR kinase signalling in fungal and animal systems, functional aspects of this pathway in plants are reviewed. Both conserved and divergent elements are discussed in relation to unique aspects associated with an autotrophic mode of nutrition and adaptive strategies for multicellular development exhibited by plants.

Identification of the Arabidopsis RAM/MOR signalling network: adding new regulatory players in plant stem cell maintenance and cell polarization

BACKGROUND AND AIMS

The RAM/MOR signalling network of eukaryotes is a conserved regulatory module involved in co-ordination of stem cell maintenance, cell differentiation and polarity establishment. To date, no such signalling network has been identified in plants.

METHODS

Genes encoding the bona fide core components of the RAM/MOR pathway were identified in Arabidopsis thaliana (arabidopsis) by sequence similarity searches conducted with the known components from other species. The transcriptional network(s) of the arabidopsis RAM/MOR signalling pathway were identified by running in-depth in silico analyses for genes co-regulated with the core components. In situ hybridization was used to confirm tissue-specific expression of selected RAM/MOR genes.

KEY RESULTS

Co-expression data suggested that the arabidopsis RAM/MOR pathway may include genes involved in floral transition, by co-operating with chromatin remodelling and mRNA processing/post-transcriptional gene silencing factors, and genes involved in the regulation of pollen tube polar growth. The RAM/MOR pathway may act upstream of the ROP1 machinery, affecting pollen tube polar growth, based on the co-expression of its components with ROP-GEFs. In silico tissue-specific co-expression data and in situ hybridization experiments suggest that different components of the arabidopsis RAM/MOR are expressed in the shoot apical meristem and inflorescence meristem and may be involved in the fine-tuning of stem cell maintenance and cell differentiation.

CONCLUSIONS

The arabidopsis RAM/MOR pathway may be part of the signalling cascade that converges in pollen tube polarized growth and in fine-tuning stem cell maintenance, differentiation and organ polarity.

Transcriptomic changes during tuber dormancy release process revealed by RNA sequencing in potato

Potato tuber dormancy release is a critical development process that allows potato to produce new plant. The first Illumina RNA sequencing to generate the expressed mRNAs at dormancy tuber (DT), dormancy release tuber (DRT) and sprouting tuber (ST) was performed. We identified 26,639 genes including 5,912 (3,450 up-regulated while 2,462 down-regulated) and 3,885 (2,141 up-regulated while 1,744 down-regulated) genes were differentially expressed from DT vs DRT and DRT vs ST. The RNA-Seq results were further verified using qRT-PCR. We found reserve mobilization events were activated before the bud emergence (DT vs DRT) and highlighted after dormancy release (DRT vs ST). Overexpressed genes related to metabolism of auxin, gibberellic acid, cytokinin and barssinosteriod were dominated in DT vs DRT, whereas overexpressed genes involved in metabolism of ethylene, jasmonate and salicylate were prominent in DRT vs ST. Various histone and cyclin isoforms associated genes involved in cell division/cycle were mainly up-regulated in DT vs DRT. Dormancy release process was also companied by stress response and redox regulation, those genes related to biotic stress, cell wall and second metabolism was preferentially overexpressed in DRT vs ST, which might accelerate dormancy breaking and sprout outgrowth. The metabolic processes activated during tuber dormancy release were also supported by plant seed models. These results represented the first comprehensive picture of a large number of genes involved in tuber dormancy release process.

The root hair "infectome" of Medicago truncatula uncovers changes in cell cycle genes and reveals a requirement for Auxin signaling in rhizobial infection

Nitrogen-fixing rhizobia colonize legume roots via plant-made intracellular infection threads. Genetics has identified some genes involved but has not provided sufficient detail to understand requirements for infection thread development. Therefore, we transcriptionally profiled Medicago truncatula root hairs prior to and during the initial stages of infection. This revealed changes in the responses to plant hormones, most notably auxin, strigolactone, gibberellic acid, and brassinosteroids. Several auxin responsive genes, including the ortholog of Arabidopsis thaliana Auxin Response Factor 16, were induced at infection sites and in nodule primordia, and mutation of ARF16a reduced rhizobial infection. Associated with the induction of auxin signaling genes, there was increased expression of cell cycle genes including an A-type cyclin and a subunit of the anaphase promoting complex. There was also induction of several chalcone O-methyltransferases involved in the synthesis of an inducer of Sinorhizobium meliloti nod genes, as well as a gene associated with Nod factor degradation, suggesting both positive and negative feedback loops that control Nod factor levels during rhizobial infection. We conclude that the onset of infection is associated with reactivation of the cell cycle as well as increased expression of genes required for hormone and flavonoid biosynthesis and that the regulation of auxin signaling is necessary for initiation of rhizobial infection threads.

Transcriptome analysis during berry development provides insights into co-regulated and altered gene expression between a seeded wine grape variety and its seedless somatic variant

BACKGROUND

Seedless grapes are greatly appreciated for fresh and dry fruit consumption. Parthenocarpy and stenospermocarpy have been described as the main phenomena responsible for seedlessness in Vitis vinifera. However, the key genes underpinning molecular and cellular processes that play a significant role in seed development are not well characterized. To identify important regulators and mechanisms that may be altered in the seedless phenotype, we performed a comprehensive transcriptional analysis to compare the transcriptomes of a popular seeded wine cultivar (wild-type) and its seedless somatic variant (mutant) at three key developmental stages.

RESULTS

The transcriptomes revealed by Illumina mRNA-Seq technology had approximately 98% of grapevine annotated transcripts and about 80% of them were commonly expressed in the two lines. Differential gene expression analysis revealed a total of 1075 differentially expressed genes (DE) in the pairwise comparison of developmental stages, which included DE genes specific to the wild-type background, DE genes specific to the mutant background and DE genes commonly shared in both backgrounds. The analysis of differential expression patterns and functional category enrichment of wild-type and mutant DE genes highlighted significant coordination and enrichment of pollen and ovule developmental pathways. The expression of some selected DE genes was further confirmed by real-time RT-PCR analysis.

CONCLUSIONS

This study represents the most comprehensive attempt to characterize the genetic bases of seed formation in grapevine. With a high throughput method, we have shown that a seeded wine grape and its seedless somatic variant are similar in several biological processes. Nevertheless, we could identify an inventory of genes with altered expression in the mutant compared to the wild-type, which may be responsible for the seedless phenotype. The genes located within known genomic regions regulating seed content may be used for the development of molecular tools to assist table grape breeding. Therefore the data reported here have provided a rich genomic resource for practical use and functional characterization of the genes that potentially underpin seedlessness in grapevine.

Comparative transcriptomics of the model mushroom Coprinopsis cinerea reveals tissue-specific armories and a conserved circuitry for sexual development

Background

It is well known that mushrooms produce defense proteins and secondary metabolites against predators and competitors; however, less is known about the correlation between the tissue-specific expression and the target organism (antagonist) specificity of these molecules. In addition, conserved transcriptional circuitries involved in developing sexual organs in fungi are not characterized, despite the growing number of gene expression datasets available from reproductive and vegetative tissue. The aims of this study were: first, to evaluate the tissue specificity of defense gene expression in the model mushroom Coprinopsis cinerea and, second, to assess the degree of conservation in transcriptional regulation during sexual development in basidiomycetes.

Results

In order to characterize the regulation in the expression of defense loci and the transcriptional circuitries controlling sexual reproduction in basidiomycetes, we sequenced the poly (A)-positive transcriptome of stage 1 primordia and vegetative mycelium of C. cinerea A43mutB43mut. Our data show that many genes encoding predicted and already characterized defense proteins are differentially expressed in these tissues. The predicted specificity of these proteins with regard to target organisms suggests that their expression pattern correlates with the type of antagonists these tissues are confronted with. Accordingly, we show that the stage 1 primordium-specific protein CC1G_11805 is toxic to insects and nematodes. Comparison of our data to analogous data from Laccaria bicolor and Schizophyllum commune revealed that the transcriptional regulation of nearly 70 loci is conserved and probably subjected to stabilizing selection. A Velvet domain-containing protein was found to be up-regulated in all three fungi, providing preliminary evidence of a possible role of the Velvet protein family in sexual development of basidiomycetes. The PBS-soluble proteome of C. cinerea primordia and mycelium was analyzed by shotgun LC-MS. This proteome data confirmed the presence of intracellular defense proteins in primordia.

Conclusions

This study shows that the exposure of different tissues in fungi to different types of antagonists shapes the expression pattern of defense loci in a tissue-specific manner. Furthermore, we identify a transcriptional circuitry conserved among basidiomycetes during fruiting body formation that involves, amongst other transcription factors, the up-regulation of a Velvet domain-containing protein.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-492) contains supplementary material, which is available to authorized users.

Recombination-independent mechanisms and pairing of homologous chromosomes during meiosis in plants

Meiosis is the specialized eukaryotic cell division that permits the halving of ploidy necessary for gametogenesis in sexually reproducing organisms. This involves a single round of DNA replication followed by two successive divisions. To ensure balanced segregation, homologous chromosome pairs must migrate to opposite poles at the first meiotic division and this means that they must recognize and pair with each other beforehand. Although understanding of the mechanisms by which meiotic chromosomes find and pair with their homologs has greatly advanced, it remains far from being fully understood. With some notable exceptions such as male Drosophila, the recognition and physical linkage of homologs at the first meiotic division involves homologous recombination. However, in addition to this, it is clear that many organisms, including plants, have also evolved a series of recombination-independent mechanisms to facilitate homolog recognition and pairing. These implicate chromosome structure and dynamics, telomeres, centromeres, and, most recently, small RNAs. With a particular focus on plants, we present here an overview of understanding of these early, recombination-independent events that act in the pairing of homologous chromosomes during the first meiotic division.

Functional analysis of sporophytic transcripts repressed by the female gametophyte in the ovule of Arabidopsis thaliana

To investigate the genetic and molecular regulation that the female gametophyte could exert over neighboring sporophytic regions of the ovule, we performed a quantitative comparison of global expression in wild-type and nozzle/sporocyteless (spl) ovules of Arabidopsis thaliana (Arabidopsis), using Massively Parallel Signature Sequencing (MPSS). This comparison resulted in 1517 genes showing at least 3-fold increased expression in ovules lacking a female gametophyte, including those encoding 89 transcription factors, 50 kinases, 25 proteins containing a RNA-recognition motif (RRM), and 20 WD40 repeat proteins. We confirmed that eleven of these genes are either preferentially expressed or exclusive of spl ovules lacking a female gametophyte as compared to wild-type, and showed that six are also upregulated in determinant infertile1 (dif1), a meiotic mutant affected in a REC8-like cohesin that is also devoided of female gametophytes. The sporophytic misexpression of IOREMPTE, a WD40/transducin repeat gene that is preferentially expressed in the L1 layer of spl ovules, caused the arrest of female gametogenesis after differentiation of a functional megaspore. Our results show that in Arabidopsis, the sporophytic-gametophytic cross talk includes a negative regulation of the female gametophyte over specific genes that are detrimental for its growth and development, demonstrating its potential to exert a repressive control over neighboring regions in the ovule.

Transcriptome database resource and gene expression atlas for the rose

BACKGROUND

For centuries roses have been selected based on a number of traits. Little information exists on the genetic and molecular basis that contributes to these traits, mainly because information on expressed genes for this economically important ornamental plant is scarce.

RESULTS

Here, we used a combination of Illumina and 454 sequencing technologies to generate information on Rosa sp. transcripts using RNA from various tissues and in response to biotic and abiotic stresses. A total of 80714 transcript clusters were identified and 76611 peptides have been predicted among which 20997 have been clustered into 13900 protein families. BLASTp hits in closely related Rosaceae species revealed that about half of the predicted peptides in the strawberry and peach genomes have orthologs in Rosa dataset. Digital expression was obtained using RNA samples from organs at different development stages and under different stress conditions. qPCR validated the digital expression data for a selection of 23 genes with high or low expression levels. Comparative gene expression analyses between the different tissues and organs allowed the identification of clusters that are highly enriched in given tissues or under particular conditions, demonstrating the usefulness of the digital gene expression analysis. A web interface ROSAseq was created that allows data interrogation by BLAST, subsequent analysis of DNA clusters and access to thorough transcript annotation including best BLAST matches on Fragaria vesca, Prunus persica and Arabidopsis. The rose peptides dataset was used to create the ROSAcyc resource pathway database that allows access to the putative genes and enzymatic pathways.

CONCLUSIONS

The study provides useful information on Rosa expressed genes, with thorough annotation and an overview of expression patterns for transcripts with good accuracy.

Major transitions in the evolution of early land plants: a bryological perspective

Background Molecular phylogeny has resolved the liverworts as the earliest-divergent clade of land plants and mosses as the sister group to hornworts plus tracheophytes, with alternative topologies resolving the hornworts as sister to mosses plus tracheophytes less well supported. The tracheophytes plus fossil plants putatively lacking lignified vascular tissue form the polysporangiophyte clade. Scope This paper reviews phylogenetic, developmental, anatomical, genetic and paleontological data with the aim of reconstructing the succession of events that shaped major land plant lineages. Conclusions Fundamental land plant characters primarily evolved in the bryophyte grade, and hence the key to a better understanding of the early evolution of land plants is in bryophytes. The last common ancestor of land plants was probably a leafless axial gametophyte bearing simple unisporangiate sporophytes. Water-conducting tissue, if present, was restricted to the gametophyte and presumably consisted of perforate cells similar to those in the early-divergent bryophytes Haplomitrium and Takakia. Stomata were a sporophyte innovation with the possible ancestral functions of producing a transpiration-driven flow of water and solutes from the parental gametophyte and facilitating spore separation before release. Stomata in mosses, hornworts and polysporangiophytes are viewed as homologous, and hence these three lineages are collectively referred to as the 'stomatophytes'. An indeterminate sporophyte body (the sporophyte shoot) developing from an apical meristem was the key innovation in polysporangiophytes. Poikilohydry is the ancestral condition in land plants; homoiohydry evolved in the sporophyte of polysporangiophytes. Fungal symbiotic associations ancestral to modern arbuscular mycorrhizas evolved in the gametophytic generation before the separation of major present-living lineages. Hydroids are imperforate water-conducting cells specific to advanced mosses. Xylem vascular cells in polysporangiophytes arose either from perforate cells or de novo. Food-conducting cells were a very early innovation in land plant evolution. The inferences presented here await testing by molecular genetics.

The function of the RNA-binding protein TEL1 in moss reveals ancient regulatory mechanisms of shoot development

The shoot represents the basic body plan in land plants. It consists of a repeated structure composed of stems and leaves. Whereas vascular plants generate a shoot in their diploid phase, non-vascular plants such as mosses form a shoot (called the gametophore) in their haploid generation. The evolution of regulatory mechanisms or genetic networks used in the development of these two kinds of shoots is unclear. TERMINAL EAR1-like genes have been involved in diploid shoot development in vascular plants. Here, we show that disruption of PpTEL1 from the moss Physcomitrella patens, causes reduced protonema growth and gametophore initiation, as well as defects in gametophore development. Leafy shoots formed on ΔTEL1 mutants exhibit shorter stems with more leaves per shoot, suggesting an accelerated leaf initiation (shortened plastochron), a phenotype shared with the Poaceae vascular plants TE1 and PLA2/LHD2 mutants. Moreover, the positive correlation between plastochron length and leaf size observed in ΔTEL1 mutants suggests a conserved compensatory mechanism correlating leaf growth and leaf initiation rate that would minimize overall changes in plant biomass. The RNA-binding protein encoded by PpTEL1 contains two N-terminus RNA-recognition motifs, and a third C-terminus non-canonical RRM, specific to TEL proteins. Removal of the PpTEL1 C-terminus (including this third RRM) or only 16-18 amino acids within it seriously impairs PpTEL1 function, suggesting a critical role for this third RRM. These results show a conserved function of the RNA-binding PpTEL1 protein in the regulation of shoot development, from early ancestors to vascular plants, that depends on the third TEL-specific RRM.

Identification and analysis of phosphorylation status of proteins in dormant terminal buds of poplar

BACKGROUND

Although there has been considerable progress made towards understanding the molecular mechanisms of bud dormancy, the roles of protein phosphorylation in the process of dormancy regulation in woody plants remain unclear.

RESULTS

We used mass spectrometry combined with TiO₂ phosphopeptide-enrichment strategies to investigate the phosphoproteome of dormant terminal buds (DTBs) in poplar (Populus simonii × P. nigra). There were 161 unique phosphorylated sites in 161 phosphopeptides from 151 proteins; 141 proteins have orthologs in Arabidopsis, and 10 proteins are unique to poplar. Only 34 sites in proteins in poplar did not match well with the equivalent phosphorylation sites of their orthologs in Arabidopsis, indicating that regulatory mechanisms are well conserved between poplar and Arabidopsis. Further functional classifications showed that most of these phosphoproteins were involved in binding and catalytic activity. Extraction of the phosphorylation motif using Motif-X indicated that proline-directed kinases are a major kinase group involved in protein phosphorylation in dormant poplar tissues.

CONCLUSIONS

This study provides evidence about the significance of protein phosphorylation during dormancy, and will be useful for similar studies on other woody plants.

Comprehensive network analysis of anther-expressed genes in rice by the combination of 33 laser microdissection and 143 spatiotemporal microarrays

Co-expression networks systematically constructed from large-scale transcriptome data reflect the interactions and functions of genes with similar expression patterns and are a powerful tool for the comprehensive understanding of biological events and mining of novel genes. In Arabidopsis (a model dicot plant), high-resolution co-expression networks have been constructed from very large microarray datasets and these are publicly available as online information resources. However, the available transcriptome data of rice (a model monocot plant) have been limited so far, making it difficult for rice researchers to achieve reliable co-expression analysis. In this study, we performed co-expression network analysis by using combined 44 K agilent microarray datasets of rice, which consisted of 33 laser microdissection (LM)-microarray datasets of anthers, and 143 spatiotemporal transcriptome datasets deposited in RicexPro. The entire data of the rice co-expression network, which was generated from the 176 microarray datasets by the Pearson correlation coefficient (PCC) method with the mutual rank (MR)-based cut-off, contained 24,258 genes and 60,441 genes pairs. Using these datasets, we constructed high-resolution co-expression subnetworks of two specific biological events in the anther, “meiosis” and “pollen wall synthesis”. The meiosis network contained many known or putative meiotic genes, including genes related to meiosis initiation and recombination. In the pollen wall synthesis network, several candidate genes involved in the sporopollenin biosynthesis pathway were efficiently identified. Hence, these two subnetworks are important demonstrations of the efficiency of co-expression network analysis in rice. Our co-expression analysis included the separated transcriptomes of pollen and tapetum cells in the anther, which are able to provide precise information on transcriptional regulation during male gametophyte development in rice. The co-expression network data presented here is a useful resource for rice researchers to elucidate important and complex biological events.

A novel RNA-recognition-motif protein is required for premeiotic G1/S-phase transition in rice (Oryza sativa L.)

The molecular mechanism for meiotic entry remains largely elusive in flowering plants. Only Arabidopsis SWI1/DYAD and maize AM1, both of which are the coiled-coil protein, are known to be required for the initiation of plant meiosis. The mechanism underlying the synchrony of male meiosis, characteristic to flowering plants, has also been unclear in the plant kingdom. In other eukaryotes, RNA-recognition-motif (RRM) proteins are known to play essential roles in germ-cell development and meiosis progression. Rice MEL2 protein discovered in this study shows partial similarity with human proline-rich RRM protein, deleted in Azoospermia-Associated Protein1 (DAZAP1), though MEL2 also possesses ankyrin repeats and a RING finger motif. Expression analyses of several cell-cycle markers revealed that, in mel2 mutant anthers, most germ cells failed to enter premeiotic S-phase and meiosis, and a part escaped from the defect and underwent meiosis with a significant delay or continued mitotic cycles. Immunofluorescent detection revealed that T7 peptide-tagged MEL2 localized at cytoplasmic perinuclear region of germ cells during premeiotic interphase in transgenic rice plants. This study is the first report of the plant RRM protein, which is required for regulating the premeiotic G1/S-phase transition of male and female germ cells and also establishing synchrony of male meiosis. This study will contribute to elucidation of similarities and diversities in reproduction system between plants and other species.

Meiosis is a pivotal event to produce haploid spores and gametes in all sexually reproducing species and is a fundamentally different type of cell cycle from mitosis. Thus, the molecular mechanisms to switch the cell cycle from mitosis to meiosis have been studied by many researchers. In yeast and metazoans, RNA-binding proteins are known to play important roles in the post-transcriptional regulation of genes implicated in the meiotic entry and meiosis. In contrast, in the plant kingdom, the mechanisms to control the meiotic entry have largely remained elusive. In this study, we discover a novel RNA-recognition-motif (RRM) protein in rice (Oryza sativa L.), designated MEL2, and demonstrate that MEL2 is required for the faithful transition of germ cells from mitosis to meiotic cell cycle. Rice MEL2 shows partial similarity with human DAZAP1, which is an RRM protein and relates to Azoospermia syndrome in human, while there are critical structural differences between germline-specific RRM proteins of mammals and plants. Our findings will lead the molecular-biological studies of plant meiotic entry to the next steps and will enable a comparison of the systems of meiotic entry between animals and plants.

Long Nonprotein-Coding RNAs in Plants

In recent years, nonprotein-coding RNAs (or npcRNAs) have emerged as a major part of the eukaryotic transcriptome. Many new regulatory npcRNAs or riboregulators riboregulators have been discovered and characterized due to the advent of new genomic approaches. This growing number suggests that npcRNAs could play a more important role than previously believed and significantly contribute to the generation of evolutionary complexity in multicellular organisms. Regulatory npcRNAs range from small RNAs (si/miRNAs) to very large transcripts (or long npcRNAs) and play diverse functions in development and/or environmental stress responses. Small RNAs include an expanding number of 20-40 nt RNAs that function in the regulation of gene expression by affecting mRNA decay and translational inhibition or lead to DNA methylation and gene silencing. They generally involve double-stranded RNA or stem loops and imply transcriptional or posttranscriptional gene silencing (PTGS). RNA silencing besides small interfering RNA and microRNA, gene silencing in plants is also mediated by tasiRNAs (trans-acting siRNAs) and nat-siRNAs (natural antisense mediated siRNAs). In contrast to small RNAs, much less is known about the large and diverse population of long npcRNAs, and only a few have been implicated in diverse functions such as abiotic stress responses, nodulation and flower development, and sex chromosome-specific expression. Moreover, many long npcRNAs act as antisense transcripts or are substrates of the small RNA pathways, thus interfering with a variety of RNA-related metabolisms. An emerging hypothesis is that long npcRNAs, as shown for small si/miRNAs, integrate into ribonucleoprotein particles (RNPs) to modulate their function, localization, or stability to act on target mRNAs. As plants show a remarkable developmental plasticity to adapt their growth to changing environmental conditions, understanding how npcRNAs work may reveal novel mechanisms involved in growth control and differentiation and help to design new tools for biotechnological applications.

Phylogenetic and expression analysis of RNA-binding proteins with triple RNA recognition motifs in plants

The superfamily of RNA binding proteins (RBPs) is vastly expanded in plants compared to other eukaryotes. A subfamily of RBPs that contain three RNA recognition motifs (RRMs) from the Arabidopsis (24), rice (19) and poplar (37) genomes was analyzed in this study. Phylogenetic analysis with full-length protein sequences of 80 RBPs identified nine clades. The largest clade, comprising 23 members, showed high homology to human RBPs involved in oxidative signaling. Digital northern analysis revealed that Arabidopsis RBPs are transcriptionally responsive to biotic, abiotic and hormonal treatments. Northern blot analysis of eight Arabidopsis RBPs belonging to the tobacco RBP45/47 family showed that these genes respond to ozone stress. AtRBP45b, which shows closest homology to the yeast oxidative stress regulatory protein, CSX1, was expressed in multiple tissues. Two novel splice variant forms of AtRBP45b were identified by 3'RACE analysis. Based on RT-PCR, splice variant AtRBP45b-SV1 was observed only in response to mechanical wounding caused by pathogen or chemical infiltrations and was not detectable in response to salt or temperature stress. Electrophoretic mobility shift assay demonstrated that recombinant full-length and splice variant forms of AtRBP45b bound synthetic RNA. Identifying in vivo RNA targets of AtRBP45b will aid in determining the precise functional role of these proteins during oxidative signaling.