Harbinger transposon insertion in ethylene signaling gene leads to emergence of new sexual forms in cucurbits.

In flowering plants, the predominant sexual morph is hermaphroditism, and the emergence of unisexuality is poorly understood. Using Cucumis melo (melon) as a model system, we explore the mechanisms driving sexual forms. We identify a spontaneous mutant exhibiting a transition from bisexual to unisexual male flower, and identify the causal mutation as a Harbinger transposon impairing the expression of Ethylene Insensitive 2 (CmEIN2) gene. Genetics and transcriptomic analysis reveal a dual role of CmEIN2 in both sex determination and fruit shape formation. Upon expression of CmACS11, EIN2 is recruited to repress the expression of the carpel inhibitor, CmWIP1. Subsequently, EIN2 is recruited to mediate stamina inhibition. Following the sex determination phase, EIN2 promotes fruit shape elongation. Genome-wide analysis reveals that Harbinger transposon mobilization is triggered by environmental cues, and integrates preferentially in active chromatin, particularly within promoter regions. Characterization of a large collection of melon germplasm points to active transpositions in the wild, compared to cultivated accessions. Our study underscores the association between chromatin dynamics and the temporal aspects of mobile genetic element insertions, providing valuable insights into plant adaptation and crop genome evolution.

Ethylene produced in carpel primordia controls CmHB40 expression to inhibit stamen development.

Sex determination evolved to control the development of unisexual flowers. In agriculture, it conditions how plants are cultivated and bred. We investigated how female flowers develop in monoecious cucurbits. We discovered in melon, Cucumis melo, a mechanism in which ethylene produced in the carpel is perceived in the stamen primordia through spatially differentially expressed ethylene receptors. Subsequently, the CmEIN3/CmEIL1 ethylene signalling module, in stamen primordia, activates the expression of CmHB40, a transcription factor that downregulates genes required for stamen development and upregulates genes associated with organ senescence. Investigation of melon genetic biodiversity revealed a haplotype, originating in Africa, altered in EIN3/EIL1 binding to CmHB40 promoter and associated with bisexual flower development. In contrast to other bisexual mutants in cucurbits, CmHB40 mutations do not alter fruit shape. By disentangling fruit shape and sex-determination pathways, our work opens up new avenues in plant breeding.

Spatially expressed WIP genes control Arabidopsis embryonic root development.

Development of plant organs is a highly organized process. In Arabidopsis, proper root development requires that distinct cell types and tissue layers are specified and formed in a restricted manner in space and over time. Despite its importance, genetic controls underlying such regularity remain elusive. Here we found that WIP genes expressed in the embryo and suspensor functionally oppose those expressed in the surrounding maternal tissues to orchestrate cell division orientation and cell fate specification in the embryonic root, thereby promoting regular root formation. The maternal WIPs act non-cell autonomously to repress root cell fate specification through SIMILAR TO RADICAL-INDUCED CELL DEATH ONE (SRO) family members. When losing all WIPs, root cells divide irregularly in the early embryo, but this barely alters their fate specification and the morphology of post-embryonic roots. Our results reveal cross-communication between the embryonic and maternal WIPs in controlling root development.

Cantaloupe melon genome reveals 3D chromatin features and structural relationship with the ancestral cucurbitaceae karyotype.

Cucumis melo displays a large diversity of horticultural groups with cantaloupe melon the most cultivated type. Using a combination of single-molecule sequencing, 10X Genomics link-reads, high-density optical and genetic maps, and chromosome conformation capture (Hi-C), we assembled a chromosome scale C. melo var. cantalupensis Charentais mono genome. Integration of RNA-seq, MeDip-seq, ChIP-seq, and Hi-C data revealed a widespread compartmentalization of the melon genome, segregating constitutive heterochromatin and euchromatin. Genome-wide comparative and evolutionary analysis between melon botanical groups identified Charentais mono genome increasingly more divergent from Harukei-3 (reticulatus), Payzawat (inodorus), and HS (ssp. agrestis) genomes. To assess the paleohistory of the Cucurbitaceae, we reconstructed the ancestral Cucurbitaceae karyotype and compared it to sequenced cucurbit genomes. In contrast to other species that experienced massive chromosome shuffling, melon has retained the ancestral genome structure. We provide comprehensive genomic resources and new insights in the diversity of melon horticultural groups and evolution of cucurbits.

The miR166-SlHB15A regulatory module controls ovule development and parthenocarpic fruit set under adverse temperatures in tomato.

Fruit set is inhibited by adverse temperatures, with consequences on yield. We isolated a tomato mutant producing fruits under non-permissive hot temperatures and identified the causal gene as SlHB15A, belonging to class III homeodomain leucine-zipper transcription factors. SlHB15A loss-of-function mutants display aberrant ovule development that mimics transcriptional changes occurring in fertilized ovules and leads to parthenocarpic fruit set under optimal and non-permissive temperatures, in field and greenhouse conditions. Under cold growing conditions, SlHB15A is subjected to conditional haploinsufficiency and recessive dosage sensitivity controlled by microRNA 166 (miR166). Knockdown of SlHB15A alleles by miR166 leads to a continuum of aberrant ovules correlating with parthenocarpic fruit set. Consistent with this, plants harboring an Slhb15a-miRNA166-resistant allele developed normal ovules and were unable to set parthenocarpic fruit under cold conditions. DNA affinity purification sequencing and RNA-sequencing analyses revealed that SlHB15A is a bifunctional transcription factor expressed in the ovule integument. SlHB15A binds to the promoters of auxin-related genes to repress auxin signaling and to the promoters of ethylene-related genes to activate their expression. A survey of tomato genetic biodiversity identified pat and pat-1, two historical parthenocarpic mutants, as alleles of SlHB15A. Taken together, our findings demonstrate the role of SlHB15A as a sentinel to prevent fruit set in the absence of fertilization and provide a mean to enhance fruiting under extreme temperatures.

A Venomics Approach Coupled to High-Throughput Toxin Production Strategies Identifies the First Venom-Derived Melanocortin Receptor Agonists.

Animal venoms are rich in hundreds of toxins with extraordinary biological activities. Their exploitation is difficult due to their complexity and the small quantities of venom available from most venomous species. We developed a Venomics approach combining transcriptomic and proteomic characterization of 191 species and identified 20,206 venom toxin sequences. Two complementary production strategies based on solid-phase synthesis and recombinant expression in Escherichia coli generated a physical bank of 3597 toxins. Screened on hMC4R, this bank gave an incredible hit rate of 8%. Here, we focus on two novel toxins: N-TRTX-Preg1a, exhibiting an inhibitory cystine knot (ICK) motif, and N-BUTX-Ptr1a, a short scorpion-CSαß structure. Neither N-TRTX-Preg1a nor N-BUTX-Ptr1a affects ion channels, the known targets of their toxin scaffolds, but binds to four melanocortin receptors with low micromolar affinities and activates the hMC1R/Gs pathway. Phylogenetically, these two toxins form new groups within their respective families and represent novel hMC1R agonists, structurally unrelated to the natural agonists.

Integrative genome-wide analysis reveals the role of WIP proteins in inhibition of growth and development.

In cucurbits, CmWIP1 is a master gene controlling sex determination. To bring new insight in the function of CmWIP1, we investigated two Arabidopsis WIP transcription factors, AtWIP1/TT1 and AtWIP2/NTT. Using an inducible system we showed that WIPs are powerful inhibitor of growth and inducer of cell death. Using ChIP-seq and RNA-seq we revealed that most of the up-regulated genes bound by WIPs display a W-box motif, associated with stress signaling. In contrast, the down-regulated genes contain a GAGA motif, a known target of polycomb repressive complex. To validate the role of WIP proteins in inhibition of growth, we expressed AtWIP1/TT1 in carpel primordia and obtained male flowers, mimicking CmWIP1 function in melon. Using other promoters, we further demonstrated that WIPs can trigger growth arrest of both vegetative and reproductive organs. Our data supports an evolutionary conserved role of WIPs in recruiting gene networks controlling growth and adaptation to stress.

LeGOO: An Expertized Knowledge Database for the Model Legume Medicago truncatula.

Medicago truncatula was proposed, about three decades ago, as a model legume to study the Rhizobium-legume symbiosis. It has now been adopted to study a wide range of biological questions, including various developmental processes (in particular root, symbiotic nodule and seed development), symbiotic (nitrogen-fixing and arbuscular mycorrhizal endosymbioses) and pathogenic interactions, as well as responses to abiotic stress. With a number of tools and resources set up in M. truncatula for omics, genetics and reverse genetics approaches, massive amounts of data have been produced, as well as four genome sequence releases. Many of these data were generated with heterogeneous tools, notably for transcriptomics studies, and are consequently difficult to integrate. This issue is addressed by the LeGOO (for Legume Graph-Oriented Organizer) knowledge base (https://www.legoo.org), which finds the correspondence between the multiple identifiers of the same gene. Furthermore, an important goal of LeGOO is to collect and represent biological information from peer-reviewed publications, whatever the technical approaches used to obtain this information. The information is modeled in a graph-oriented database, which enables flexible representation, with currently over 200,000 relations retrieved from 298 publications. LeGOO also provides the user with mining tools, including links to the Mt5.0 genome browser and associated information (on gene functional annotation, expression, methylome, natural diversity and available insertion mutants), as well as tools to navigate through different model species. LeGOO is, therefore, an innovative database that will be useful to the Medicago and legume community to better exploit the wealth of data produced on this model species.

A miR172 target-deficient AP2-like gene correlates with the double flower phenotype in roses.

One of the well-known floral abnormalities in flowering plants is the double-flower phenotype, which corresponds to flowers that develop extra petals, sometimes even containing entire flowers within flowers. Because of their highly priced ornamental value, spontaneous double-flower variants have been found and selected for in a wide range of ornamental species. Previously, double flower formation in roses was associated with a restriction of AGAMOUS expression domain toward the centre of the meristem, leading to extra petals. Here, we characterized the genomic region containing the mutation associated with the switch from simple to double flowers in the rose. An APETALA2-like gene (RcAP2L), a member of the Target Of EAT-type (TOE-type) subfamily, lies within this interval. In the double flower rose, two alleles of RcAP2L are present, one of which harbours a transposable element inserted into intron 8. This insertion leads to the creation of a miR172 resistant RcAP2L variant. Analyses of the presence of this variant in a set of simple and double flower roses demonstrate a correlation between the presence of this allele and the double flower phenotype. These data suggest a role of this miR172 resistant RcAP2L variant in regulating RcAGAMOUS expression and double flower formation in Rosa sp.

The Rosa genome provides new insights into the domestication of modern roses.

Roses have high cultural and economic importance as ornamental plants and in the perfume industry. We report the rose whole-genome sequencing and assembly and resequencing of major genotypes that contributed to rose domestication. We generated a homozygous genotype from a heterozygous diploid modern rose progenitor, Rosa chinensis 'Old Blush'. Using single-molecule real-time sequencing and a meta-assembly approach, we obtained one of the most comprehensive plant genomes to date. Diversity analyses highlighted the mosaic origin of 'La France', one of the first hybrids combining the growth vigor of European species and the recurrent blooming of Chinese species. Genomic segments of Chinese ancestry identified new candidate genes for recurrent blooming. Reconstructing regulatory and secondary metabolism pathways allowed us to propose a model of interconnected regulation of scent and flower color. This genome provides a foundation for understanding the mechanisms governing rose traits and should accelerate improvement in roses, Rosaceae and ornamentals.

Diversity in sequences, post-translational modifications and expected pharmacological activities of toxins from four Conus species revealed by the combination of cutting-edge proteomics, transcriptomics and bioinformatics.

Venomous animals have developed a huge arsenal of reticulated peptides for defense and predation. Based on various scaffolds, they represent a colossal pharmacological diversity, making them top candidates for the development of innovative drugs. Instead of relying on the classical, low-throughput bioassay-guided approach to identify innovative bioactive peptides, this work exploits a recent paradigm to access to venom diversity. This strategy bypasses the classical approach by combining high-throughput transcriptomics, proteomics and bioinformatics cutting-edge technologies to generate reliable peptide sequences. The strategy employed to generate hundreds of reliable sequences from Conus venoms is deeply described. The study led to the discovery of (i) conotoxins that belong to known pharmacological families targeting various GPCRs or ion-gated channels, and (ii) new families of conotoxins, never described to date. It also focusses on the diversity of genes, sequences, folds, and PTM's provided by such species.

De Novo sequencing and transcriptome analysis for Tetramorium bicarinatum: a comprehensive venom gland transcriptome analysis from an ant species.

BACKGROUND: Arthropod venoms are invaluable sources of bioactive substances with biotechnological application. The limited availability of some venoms, such as those from ants, has restricted the knowledge about the composition and the potential that these biomolecules could represent. In order to provide a global insight on the transcripts expressed in the venom gland of the Brazilian ant species Tetramorium bicarinatum and to unveil the potential of its products, high-throughput approach using Illumina technology has been applied to analyze the genes expressed in active venom glands of this ant species. RESULTS: A total of 212,371,758 pairs of quality-filtered, 100-base-pair Illumina reads were obtained. The de novo assemblies yielded 36,042 contigs for which 27,873 have at least one predicted ORF among which 59.77% produce significant hits in the available databases. The investigation of the reads mapping toxin class revealed a high diversification with the major part consistent with the classical hymenopteran venom protein signature represented by venom allergen (33.3%), followed by a diverse toxin-expression profile including several distinct isoforms of phospholipase A1 and A2, venom serine protease, hyaluronidase, protease inhibitor and secapin. Moreover, our results revealed for the first time the presence of toxin-like peptides that have been previously identified from unrelated venomous animals such as waprin-like (snakes) and agatoxins (spiders and conus).The non-toxin transcripts were mainly represented by contigs involved in protein folding and translation, consistent with the protein-secretory function of the venom gland tissue. Finally, about 40% of the generated contigs have no hits in the databases with 25% of the predicted peptides bearing signal peptide emphasizing the potential of the investigation of these sequences as source of new molecules. Among these contigs, six putative novel peptides that show homologies with previously identified antimicrobial peptides were identified. CONCLUSIONS: To the best of our knowledge, this work reports the first large-scale analysis of genes transcribed by the venomous gland of the ant species T. bicarinatum and helps with the identification of Hymenoptera toxin arsenal. In addition, results from this study demonstrate that de novo transcriptome assembly allows useful venom gene expression analysis in a species lacking a genome sequence database.

Transcriptomics and venomics: implications for medicinal chemistry.

Over the last three decades, transcriptomic studies of venom gland cells have continuously evolved, opening up new possibilities for exploring the molecular diversity of animal venoms, a prerequisite for the discovery of new drug candidates and molecular phylogenetics. The molecular complexity of animal venoms is much greater than initially thought. In this review, we describe the different technologies available for transcriptomic studies of venom, from the original individual cloning approaches to the more recent global Next Generation Sequencing strategies. Our understanding of animal venoms is evolving, with the discovery of complex and diverse bio-optimized cocktails of compounds, including mostly peptides and proteins, which are now beginning to be studied by academic and industrial researchers.

High-throughput production of two disulphide-bridge toxins.

A quick and efficient production method compatible with high-throughput screening was developed using 36 toxins belonging to four different families of two disulphide-bridge toxins. Final toxins were characterized using HPLC co-elution, CD and pharmacological studies.

Atractaspis aterrima toxins: the first insight into the molecular evolution of venom in side-stabbers.

Although snake venoms have been the subject of intense research, primarily because of their tremendous potential as a bioresource for design and development of therapeutic compounds, some specific groups of snakes, such as the genus Atractaspis, have been completely neglected. To date only limited number of toxins, such as sarafotoxins have been well characterized from this lineage. In order to investigate the molecular diversity of venom from Atractaspis aterrima-the slender burrowing asp, we utilized a high-throughput transcriptomic approach completed with an original bioinformatics analysis pipeline. Surprisingly, we found that Sarafotoxins do not constitute the major ingredient of the transcriptomic cocktail; rather a large number of previously well-characterized snake venom-components were identified. Notably, we recovered a large diversity of three-finger toxins (3FTxs), which were found to have evolved under the significant influence of positive selection. From the normalized and non-normalized transcriptome libraries, we were able to evaluate the relative abundance of the different toxin groups, uncover rare transcripts, and gain new insight into the transcriptomic machinery. In addition to previously characterized toxin families, we were able to detect numerous highly-transcribed compounds that possess all the key features of venom-components and may constitute new classes of toxins.

Profiling the venom gland transcriptome of Tetramorium bicarinatum (Hymenoptera: Formicidae): the first transcriptome analysis of an ant species.

Animal venoms are complex mixtures containing a range of bioactive elements with potential pharmacological and therapeutic use. Even though ants account among the most diverse zoological group, little information is available regarding their venom composition. To initiate the characterization of the transcriptomic venom gland expression of the ant species Tetramorium bicarinatum, 400 randomly selected clones from cDNA library were sequenced and a total of 364 high quality expressed sequence tags (ESTs) were generated. Based on the results of BLAST searches, these sequences were clustered and assembled into 83 contigs (22 multiple sequences) and 61 singletons. About 74% (267) of the contigs matched BLASTx hits with an interesting diversity together with an unusual abundance of cellular transcripts related to gene expression regulation (29% of the total library) reflecting the specialization of this tissue. About eighteen per cent of the ESTs were categorized as Hymenoptera venom compounds, the major part represented by allergens (62% of the total venom compounds). In addition, a high number of sequences (26%) had no similarity to any known sequences. This study provides a first insight of the gene expression scenario of the venom gland of T. bicarinatum which might contribute to acquiring a more comprehensive view on the origin and functional diversity of venom proteins among ants and more broadly among Hymenopteran insects.

MtbHLH1, a bHLH transcription factor involved in Medicago truncatula nodule vascular patterning and nodule to plant metabolic exchanges.

This study aimed at defining the role of a basic helix-loop-helix (bHLH) transcription factor gene from Medicago truncatula, MtbHLH1, whose expression is upregulated during the development of root nodules produced upon infection by rhizobia bacteria. We used MtbHLH1 promoter::GUS fusions and quantitative reverse-transcription polymerase chain reaction analyses to finely characterize the MtbHLH1 expression pattern. We altered MtbHLH1 function by expressing a dominantly repressed construct (CRES-T approach) and looked for possible MtbHLH1 target genes by transcriptomics. We found that MtbHLH1 is expressed in nodule primordia cells derived from pericycle divisions, in nodule vascular bundles (VBs) and in uninfected cells of the nitrogen (N) fixation zone. MtbHLH1 is also expressed in root tips, lateral root primordia cells and root VBs, and induced upon auxin treatment. Altering MtbHLH1 function led to an unusual phenotype, with a modified patterning of nodule VB development and a reduced growth of aerial parts of the plant, even though the nodules were able to fix atmospheric N. Several putative MtbHLH1 regulated genes were identified, including an asparagine synthase and a LOB (lateral organ boundary) transcription factor. Our results suggest that the MtbHLH1 gene is involved in the control of nodule vasculature patterning and nutrient exchanges between nodules and roots.

Transcription reprogramming during root nodule development in Medicago truncatula.

Many genes which are associated with root nodule development and activity in the model legume Medicago truncatula have been described. However information on precise stages of activation of these genes and their corresponding transcriptional regulators is often lacking. Whether these regulators are shared with other plant developmental programs also remains an open question. Here detailed microarray analyses have been used to study the transcriptome of root nodules induced by either wild type or mutant strains of Sinorhizobium meliloti. In this way we have defined eight major activation patterns in nodules and identified associated potential regulatory genes. We have shown that transcription reprogramming during consecutive stages of nodule differentiation occurs in four major phases, respectively associated with (i) early signalling events and/or bacterial infection; plant cell differentiation that is either (ii) independent or (iii) dependent on bacteroid differentiation; (iv) nitrogen fixation. Differential expression of several genes involved in cytokinin biosynthesis was observed in early symbiotic nodule zones, suggesting that cytokinin levels are actively controlled in this region. Taking advantage of databases recently developed for M. truncatula, we identified a small subset of gene expression regulators that were exclusively or predominantly expressed in nodules, whereas most other regulators were also activated under other conditions, and notably in response to abiotic or biotic stresses. We found evidence suggesting the activation of the jasmonate pathway in both wild type and mutant nodules, thus raising questions about the role of jasmonate during nodule development. Finally, quantitative RT-PCR was used to analyse the expression of a series of nodule regulator and marker genes at early symbiotic stages in roots and allowed us to distinguish several early stages of gene expression activation or repression.

Environmental regulation of lateral root emergence in Medicago truncatula requires the HD-Zip I transcription factor HB1.

The adaptation of root architecture to environmental constraints is a major agricultural trait, notably in legumes, the third main crop worldwide. This root developmental plasticity depends on the formation of lateral roots (LRs) emerging from primary roots. In the model legume Medicago truncatula, the HD-Zip I transcription factor HB1 is expressed in primary and lateral root meristems and induced by salt stress. Constitutive expression of HB1 in M. truncatula roots alters their architecture, whereas hb1 TILLING mutants showed increased lateral root emergence. Electrophoretic mobility shift assay, promoter mutagenesis, and chromatin immunoprecipitation-PCR assays revealed that HB1 directly recognizes a CAATAATTG cis-element present in the promoter of a LOB-like (for Lateral Organ Boundaries) gene, LBD1, transcriptionally regulated by auxin. Expression of these genes in response to abscisic acid and auxin and their behavior in hb1 mutants revealed an HB1-mediated repression of LBD1 acting during LR emergence. M. truncatula HB1 regulates an adaptive developmental response to minimize the root surface exposed to adverse environmental stresses.