GreeNC 2.0: a comprehensive database of plant long non-coding RNAs.

The Green Non-Coding Database (GreeNC) is one of the reference databases for the study of plant long non-coding RNAs (lncRNAs). Here we present our most recent update where 16 species have been updated, while 78 species have been added, resulting in the annotation of more than 495 000 lncRNAs. Moreover, sequence clustering was applied providing information about sequence conservation and gene families. The current version of the database is available at: http://greenc.sequentiabiotech.com/wiki2/Main_Page.

PRGdb 4.0: an updated database dedicated to genes involved in plant disease resistance process.

The Plant Resistance Genes database (PRGdb; http://prgdb.org/prgdb4/) has been greatly expanded, keeping pace with the increasing amount of available knowledge and data (sequenced proteomes, cloned genes, public analysis data, etc.). The easy-to-use style of the database website has been maintained, while an updated prediction tool, more data and a new section have been added. This new section will contain plant resistance transcriptomic experiments, providing additional easy-to-access experimental information. DRAGO3, the tool for automatic annotation and prediction of plant resistance genes behind PRGdb, has been improved in both accuracy and sensitivity, leading to more reliable predictions. PRGdb offers 199 reference resistance genes and 586.652 putative resistance genes from 182 sequenced proteomes. Compared to the previous release, PRGdb 4.0 has increased the number of reference resistance genes from 153 to 199, the number of putative resistance genes from 177K from 76 proteomes to 586K from 182 sequenced proteomes. A new section has been created that collects plant-pathogen transcriptomic data for five species of agricultural interest. Thereby, with these improvements and data expansions, PRGdb 4.0 aims to serve as a reference to the plant scientific community and breeders worldwide, helping to further study plant resistance mechanisms that contribute to fighting pathogens.

Transcriptomic profiles of the ruminal wall in Italian Mediterranean dairy buffaloes fed green forage.

BACKGROUND: Green feed diet in ruminants exerts a beneficial effect on rumen metabolism and enhances the content of milk nutraceutical quality. At present, a comprehensive analysis focused on the identification of genes, and therefore, biological processes modulated by the green feed in buffalo rumen has never been reported. We performed RNA-sequencing in the rumen of buffaloes fed a total mixed ration (TMR) + the inclusion of 30% of ryegrass green feed (treated) or TMR (control), and identified differentially expressed genes (DEGs) using EdgeR and NOISeq tools. RESULTS: We found 155 DEGs using EdgeR (p-values < 0.05) and 61 DEGs using NOISeq (prob ≥0.8), 30 of which are shared. The rt-qPCR validation suggested a higher reliability of EdgeR results as compared with NOISeq data, in our biological context. Gene Ontology analysis of DEGs identified using EdgeR revealed that green feed modulates biological processes relevant for the rumen physiology and, then, health and well-being of buffaloes, such as lipid metabolism, response to the oxidative stress, immune response, and muscle structure and function. Accordingly, we found: (i) up-regulation of HSD17B13, LOC102410803 (or PSAT1) and HYKK, and down-regulation of CDO1, SELENBP1 and PEMT, encoding factors involved in energy, lipid and amino acid metabolism; (ii) enhanced expression of SIM2 and TRIM14, whose products are implicated in the immune response and defense against infections, and reduced expression of LOC112585166 (or SAAL1), ROR2, SMOC2, and S100A11, encoding pro-inflammatory factors; (iii) up-regulation of NUDT18, DNAJA4 and HSF4, whose products counteract stressful conditions, and down-regulation of LOC102396388 (or UGT1A9) and LOC102413340 (or MRP4/ABCC4), encoding detoxifying factors; (iv) increased expression of KCNK10, CACNG4, and ATP2B4, encoding proteins modulating Ca2+ homeostasis, and reduced expression of the cytoskeleton-related MYH11 and DES. CONCLUSION: Although statistically unpowered, this study suggests that green feed modulates the expression of genes involved in biological processes relevant for rumen functionality and physiology, and thus, for welfare and quality production in Italian Mediterranean dairy buffaloes. These findings, that need to be further confirmed through the validation of additional DEGs, allow to speculate a role of green feed in the production of nutraceutical molecules, whose levels might be enhanced also in milk.

Comparison of Tomato Transcriptomic Profiles Reveals Overlapping Patterns in Abiotic and Biotic Stress Responses.

Until a few years ago, many studies focused on the transcriptomic response to single stresses. However, tomato cultivations are often constrained by a wide range of biotic and abiotic stress that can occur singularly or in combination, and several genes can be involved in the defensive mechanism response. Therefore, we analyzed and compared the transcriptomic responses of resistant and susceptible genotypes to seven biotic stresses (Cladosporium fulvum, Phytophthora infestans, Pseudomonas syringae, Ralstonia solanacearum, Sclerotinia sclerotiorum, Tomato spotted wilt virus (TSWV) and Tuta absoluta) and five abiotic stresses (drought, salinity, low temperatures, and oxidative stress) to identify genes involved in response to multiple stressors. With this approach, we found genes encoding for TFs, phytohormones, or participating in signaling and cell wall metabolic processes, participating in defense against various biotic and abiotic stress. Moreover, a total of 1474 DEGs were commonly found between biotic and abiotic stress. Among these, 67 DEGs were involved in response to at least four different stresses. In particular, we found RLKs, MAPKs, Fasciclin-like arabinogalactans (FLAs), glycosyltransferases, genes involved in the auxin, ET, and JA pathways, MYBs, bZIPs, WRKYs and ERFs genes. Detected genes responsive to multiple stress might be further investigated with biotechnological approaches to effectively improve plant tolerance in the field.

Exploitation of epigenetic variation of crop wild relatives for crop improvement and agrobiodiversity preservation.

Crop wild relatives (CWRs) are recognized as the best potential source of traits for crop improvement. However, successful crop improvement using CWR relies on identifying variation in genes controlling desired traits in plant germplasms and subsequently incorporating them into cultivars. Epigenetic diversity may provide an additional layer of variation within CWR and can contribute novel epialleles for key traits for crop improvement. There is emerging evidence that epigenetic variants of functional and/or agronomic importance exist in CWR gene pools. This provides a rationale for the conservation of epigenotypes of interest, thus contributing to agrobiodiversity preservation through conservation and (epi)genetic monitoring. Concepts and techniques of classical and modern breeding should consider integrating recent progress in epigenetics, initially by identifying their association with phenotypic variations and then by assessing their heritability and stability in subsequent generations. New tools available for epigenomic analysis offer the opportunity to capture epigenetic variation and integrate it into advanced (epi)breeding programmes. Advances in -omics have provided new insights into the sources and inheritance of epigenetic variation and enabled the efficient introduction of epi-traits from CWR into crops using epigenetic molecular markers, such as epiQTLs.

isoCNV: in silico optimization of copy number variant detection from targeted or exome sequencing data.

BACKGROUND: Accurate copy number variant (CNV) detection is especially challenging for both targeted sequencing (TS) and whole-exome sequencing (WES) data. To maximize the performance, the parameters of the CNV calling algorithms should be optimized for each specific dataset. This requires obtaining validated CNV information using either multiplex ligation-dependent probe amplification (MLPA) or array comparative genomic hybridization (aCGH). They are gold standard but time-consuming and costly approaches. RESULTS: We present isoCNV which optimizes the parameters of DECoN algorithm using only NGS data. The parameter optimization process is performed using an in silico CNV validated dataset obtained from the overlapping calls of three algorithms: CNVkit, panelcn.MOPS and DECoN. We evaluated the performance of our tool and showed that increases the sensitivity in both TS and WES real datasets. CONCLUSIONS: isoCNV provides an easy-to-use pipeline to optimize DECoN that allows the detection of analysis-ready CNV from a set of DNA alignments obtained under the same conditions. It increases the sensitivity of DECoN without the need for orthogonal methods. isoCNV is available at https://gitlab.com/sequentiateampublic/isocnv .

Recombination suppression in heterozygotes for a pericentric inversion induces the interchromosomal effect on crossovers in Arabidopsis.

During meiosis, recombination ensures allelic exchanges through crossovers (COs) between the homologous chromosomes. Advances in our understanding of the rules of COs have come from studies of mutations including structural chromosomal rearrangements that, when heterozygous, are known to impair COs in various organisms. In this work, we investigate the effect of a large heterozygous pericentric inversion on male and female recombination in Arabidopsis. The inversion was discovered in the Atmcc1 mutant background and was characterized through genetic and next-generation sequencing analysis. Reciprocal backcross populations, each consisting of over 400 individuals, obtained from the mutant and the wild type, both crossed with Landsberg erecta, were analyzed genome-wide by 143 single-nucleotide polymorphisms. The negative impact of inversion became evident in terms of CO loss in the rearranged chromosome in both male and female meiosis. No single-CO event was detected within the inversion, consistent with a post-meiotic selection operating against unbalanced gametes. Cytological analysis of chiasmata in F1 plants confirmed that COs were reduced in male meiosis in the chromosome with inversion. Crossover suppression on the rearranged chromosome is associated with a significant increase of COs in the other chromosomes, thereby maintaining unchanged the number of COs per cell. The CO pattern observed in our study is consistent with the interchromosomal (IC) effect as first described in Drosophila. In contrast to male meiosis, in female meiosis no IC effect is visible. This may be related to the greater strength of interference that constrains the CO number in excess of the minimum value imposed by CO assurance in Arabidopsis female meiosis.

The zoospores of the thraustochytrid Aurantiochytrium limacinum: Transcriptional reprogramming and lipid metabolism associated to their specific functions.

Aurantiochytrium limacinum (Thraustochytriaceae, class Labyrinthulomycetes) is a marine Stramenopile and a pioneering mangrove decomposer. Its life cycle involves a non-motile stage and zoospore production. We observed that the composition of the medium, the presence of amino acids in particular, affects the release of zoospores. Two opposite conditions were defined, one with a cell population mainly composed of zoospores and another one with almost only non-motile cells. In silico allelic frequency analysis and flow cytometry suggest that zoospores and non-motile cells share the same ploidy level and are diploid. Through an RNA-seq approach, the transcriptional reprogramming accompanying the formation of zoospores was investigated, with a particular focus on their lipid metabolism. Based on a differential expression analysis, zoospores are characterized by high motility, very active signal transduction, an arrest of the cell division, a low amino acid metabolism and low glycolysis. Focusing on lipid metabolism, genes involved in lipase activities and peroxisomal ß-oxidation are upregulated. qRT-PCR of selected lipid genes and lipid analyses during the life span of zoospores confirmed these observations. These results highlight the importance of the lipid dynamics in zoospores and show the metabolic processes required to use these energy-dense molecules as fuel for zoospore survival during their quest of new territories.

PRGdb 3.0: a comprehensive platform for prediction and analysis of plant disease resistance genes.

The Plant Resistance Genes database (PRGdb; http://prgdb.org) has been redesigned with a new user interface, new sections, new tools and new data for genetic improvement, allowing easy access not only to the plant science research community but also to breeders who want to improve plant disease resistance. The home page offers an overview of easy-to-read search boxes that streamline data queries and directly show plant species for which data from candidate or cloned genes have been collected. Bulk data files and curated resistance gene annotations are made available for each plant species hosted. The new Gene Model view offers detailed information on each cloned resistance gene structure to highlight shared attributes with other genes. PRGdb 3.0 offers 153 reference resistance genes and 177 072 annotated candidate Pathogen Receptor Genes (PRGs). Compared to the previous release, the number of putative genes has been increased from 106 to 177 K from 76 sequenced Viridiplantae and algae genomes. The DRAGO 2 tool, which automatically annotates and predicts (PRGs) from DNA and amino acid with high accuracy and sensitivity, has been added. BLAST search has been implemented to offer users the opportunity to annotate and compare their own sequences. The improved section on plant diseases displays useful information linked to genes and genomes to connect complementary data and better address specific needs. Through, a revised and enlarged collection of data, the development of new tools and a renewed portal, PRGdb 3.0 engages the plant science community in developing a consensus plan to improve knowledge and strategies to fight diseases that afflict main crops and other plants.

Whole-Genome Sequencing of an Uncommon Cryptococcus neoformans MLST43 Genotype Isolated in Nigeria.

Cryptococcosis is a human infection caused mainly by two species of the Cryptococcus genus, Cryptococcus neoformans and Cryptococcus gattii, whose populations contain several phylogenetically related haploid (VN/VG-types) and hybrid genotypes. Here, we report the whole-genome sequencing of a Nigerian C. neoformans VNII, Mat-α, strain with a rare multilocus-sequence-type (MLST) genotype (ST43).

GREENC: a Wiki-based database of plant lncRNAs.

Long non-coding RNAs (lncRNAs) are functional non-translated molecules greater than 200 nt. Their roles are diverse and they are usually involved in transcriptional regulation. LncRNAs still remain largely uninvestigated in plants with few exceptions. Experimentally validated plant lncRNAs have been shown to regulate important agronomic traits such as phosphate starvation response, flowering time and interaction with symbiotic organisms, making them of great interest in plant biology and in breeding. There is still a lack of lncRNAs in most sequenced plant species, and in those where they have been annotated, different methods have been used, so making the lncRNAs less useful in comparisons within and between species. We developed a pipeline to annotate lncRNAs and applied it to 37 plant species and six algae, resulting in the annotation of more than 120 000 lncRNAs. To facilitate the study of lncRNAs for the plant research community, the information gathered is organised in the Green Non-Coding Database (GreeNC, http://greenc.sciencedesigners.com/).

Maize RNA PolIV affects the expression of genes with nearby TE insertions and has a genome-wide repressive impact on transcription.

BACKGROUND: RNA-directed DNA methylation (RdDM) is a plant-specific epigenetic process that relies on the RNA polymerase IV (Pol IV) for the production of 24 nucleotide small interfering RNAs (siRNA) that guide the cytosine methylation and silencing of genes and transposons. Zea mays RPD1/RMR6 gene encodes the largest subunit of Pol IV and is required for normal plant development, paramutation, transcriptional repression of certain transposable elements (TEs) and transcriptional regulation of specific alleles. RESULTS: In this study we applied a total RNA-Seq approach to compare the B73 and rpd1/rmr6 leaf transcriptomes. Although previous studies indicated that loss of siRNAs production in RdDM mutants provokes a strong loss of CHH DNA methylation but not massive gene or TEs transcriptional activation in both Arabidopsis and maize, our total RNA-Seq analysis of rpd1/rmr6 transcriptome reveals that loss of Pol IV activity causes a global increase in the transcribed fraction of the maize genome. Our results point to the genes with nearby TE insertions as being the most strongly affected by Pol IV-mediated gene silencing. TEs modulation of nearby gene expression is linked to alternative methylation profiles on gene flanking regions, and these profiles are strictly dependent on specific characteristics of the TE member inserted. Although Pol IV is essential for the biogenesis of siRNAs, the genes with associated siRNA loci are less affected by the pol IV mutation. CONCLUSIONS: This deep and integrated analysis of gene expression, TEs distribution, smallRNA targeting and DNA methylation levels, reveals that loss of Pol IV activity globally affects genome regulation, pointing at TEs as modulator of nearby gene expression and indicating the existence of multiple level epigenetic silencing mechanisms. Our results also suggest a predominant role of the Pol IV-mediated RdDM pathway in genome dominance regulation, and subgenome stability and evolution in maize.

Effects of Mecp2 loss of function in embryonic cortical neurons: a bioinformatics strategy to sort out non-neuronal cells variability from transcriptome profiling.

BACKGROUND: Mecp2 null mice model Rett syndrome (RTT) a human neurological disorder affecting females after apparent normal pre- and peri-natal developmental periods. Neuroanatomical studies in cerebral cortex of RTT mouse models revealed delayed maturation of neuronal morphology and autonomous as well as non-cell autonomous reduction in dendritic complexity of postnatal cortical neurons. However, both morphometric parameters and high-resolution expression profile of cortical neurons at embryonic developmental stage have not yet been studied. Here we address these topics by using embryonic neuronal primary cultures from Mecp2 loss of function mouse model. RESULTS: We show that embryonic primary cortical neurons of Mecp2 null mice display reduced neurite complexity possibly reflecting transcriptional changes. We used RNA-sequencing coupled with a bioinformatics comparative approach to identify and remove the contribution of variable and hard to quantify non-neuronal brain cells present in our in vitro cell cultures. CONCLUSIONS: Our results support the need to investigate both Mecp2 morphological as well as molecular effect in neurons since prenatal developmental stage, long time before onset of Rett symptoms.

Spermatozoa from infertile patients exhibit differences of DNA methylation associated with spermatogenesis-related processes: an array-based analysis.

The influence of aberrant sperm DNA methylation on the reproductive capacity of couples has been postulated as a cause of infertility. This study compared the DNA methylation of spermatozoa of 19 fertile donors and 42 infertile patients using the Illumina 450K array. Clustering analysis of methylation data arranged fertile and infertile patients into two groups. Bivariate clustering analysis identified a differential distribution of samples according to the characteristics of seminogram and age, suggesting a possible link between these parameters and specific methylation profiles. The study identified 696 differentially methylated cytosine-guanine dinucleotides (CpG) associated with 501 genes between fertile donors and infertile patients. Ontological enrichment analysis revealed 13 processes related to spermatogenesis. Data filtering identified a set of 17 differentially methylated genes, some of which had functions relating to spermatogenesis. A significant association was identified between RPS6KA2 hypermethylation and advanced age (P = 0.016); APCS hypermethylation and oligozoospermia (P = 0.041); JAM3/NCAPD3 hypermethylation and numerical chromosome sperm anomalies (P = 0.048); and ANK2 hypermethylation and lower pregnancy rate (P = 0.040). This description of a set of differentially methylated genes provides a framework for further investigation into the influence of such variation in male fertility in larger patient cohorts.

LOXL2-mediated chromatin compaction is required to maintain the oncogenic properties of triple-negative breast cancer cells.

Oxidation of histone H3 at lysine 4 (H3K4ox) is catalyzed by lysyl oxidase homolog 2 (LOXL2). This histone modification is enriched in heterochromatin in triple-negative breast cancer (TNBC) cells and has been linked to the maintenance of compacted chromatin. However, the molecular mechanism underlying this maintenance is still unknown. Here, we show that LOXL2 interacts with RuvB-Like 1 (RUVBL1), RuvB-Like 2 (RUVBL2), Actin-like protein 6A (ACTL6A), and DNA methyltransferase 1associated protein 1 (DMAP1), a complex involved in the incorporation of the histone variant H2A.Z. Our experiments indicate that this interaction and the active form of RUVBL2 are required to maintain LOXL2-dependent chromatin compaction. Genome-wide experiments showed that H2A.Z, RUVBL2, and H3K4ox colocalize in heterochromatin regions. In the absence of LOXL2 or RUVBL2, global levels of the heterochromatin histone mark H3K9me3 were strongly reduced, and the ATAC-seq signal in the H3K9me3 regions was increased. Finally, we observed that the interplay between these series of events is required to maintain H3K4ox-enriched heterochromatin regions, which in turn is key for maintaining the oncogenic properties of the TNBC cell line tested (MDA-MB-231).

Genome-wide analysis of histone modifiers in tomato: gaining an insight into their developmental roles.

BACKGROUND: Histone post-translational modifications (HPTMs) including acetylation and methylation have been recognized as playing a crucial role in epigenetic regulation of plant growth and development. Although Solanum lycopersicum is a dicot model plant as well as an important crop, systematic analysis and expression profiling of histone modifier genes (HMs) in tomato are sketchy. RESULTS: Based on recently released tomato whole-genome sequences, we identified in silico 32 histone acetyltransferases (HATs), 15 histone deacetylases (HDACs), 52 histone methytransferases (HMTs) and 26 histone demethylases (HDMs), and compared them with those detected in Arabidopsis (Arabidopsis thaliana), maize (Zea mays) and rice (Oryza sativa) orthologs. Comprehensive analysis of the protein domain architecture and phylogeny revealed the presence of non-canonical motifs and new domain combinations, thereby suggesting for HATs the existence of a new family in plants. Due to species-specific diversification during evolutionary history tomato has fewer HMs than Arabidopsis. The transcription profiles of HMs within tomato organs revealed a broad functional role for some HMs and a more specific activity for others, suggesting key HM regulators in tomato development. Finally, we explored S. pennellii introgression lines (ILs) and integrated the map position of HMs, their expression profiles and the phenotype of ILs. We thereby proved that the strategy was useful to identify HM candidates involved in carotenoid biosynthesis in tomato fruits. CONCLUSIONS: In this study, we reveal the structure, phylogeny and spatial expression of members belonging to the classical families of HMs in tomato. We provide a framework for gene discovery and functional investigation of HMs in other Solanaceae species.

Phylogenetic and Structural Analyses of VPS13 Proteins in Archaeplastida Reveal Their Complex Evolutionary History in Viridiplantae.

VPS13 is a lipid transfer protein family conserved among Eukaryotes and playing roles in fundamental processes involving vesicular transport and membrane expansion including autophagy and organelle biogenesis. VPS13 folds into a long hydrophobic tunnel, allowing lipid transport, decorated by distinct domains involved in protein localization and regulation. Whereas VPS13 organization and function have been extensively studied in yeast and mammals, information in organisms originating from primary endosymbiosis is scarce. In the higher plant Arabidopsis thaliana, four paralogs, AtVPS13S, X, M1, and M2, were identified, AtVPS13S playing a role in the regulation of root growth, cell patterning, and reproduction. In this work, we performed phylogenetic, as well as domain and structural modeling of VPS13 proteins in Archaeplastida in order to understand their general organization and evolutionary history. We confirmed the presence of human VPS13B orthologues in some phyla and described two new VPS13 families presenting a particular domain arrangement: VPS13R in Rhodophytes and VPS13Y in Chlorophytes and Streptophytes. By focusing on Viridiplantae, we were able to draw the evolutionary history of these proteins made by multiple gene gains and duplications as well as domain rearrangements. We showed that some Chlorophytes have only three (AtVPS13M, S, Y) whereas some Charophytes have up to six VPS13 paralogs (AtVPS13M1, M2, S, Y, X, B). We also highlighted specific structural features of VPS13M and X paralogs. This study reveals the complex evolution of VPS13 family and opens important perspectives for their functional characterization in photosynthetic organisms.

A transcriptomic insight into the human sperm microbiome through next-generation sequencing.

The purpose of this study is to provide novel information through Next Generation Sequencing (NGS) for the characterization of viral and bacterial RNA cargo of human sperm cells from healthy fertile donors. For this, RNA-seq raw data of poly(A) RNA from 12 sperm samples from fertile donors were aligned to microbiome databases using the GAIA software. Species of viruses and bacteria were quantified in Operational Taxonomic Units (OTU) and filtered by minimal expression level (>1% OTU in at least one sample). Mean expression values (and their standard deviation) of each species were estimated. A Hierarchical Cluster Analysis (HCA) and a Principal Component Analysis (PCA) were performed to detect common microbiome patterns among samples. Sixteen microbiome species, families, domains, and orders surpassed the established expression threshold. Of the 16 categories, nine corresponded to viruses (23.07% OTU) and seven to bacteria (2.77% OTU), among which the Herperviriales order and Escherichia coli were the most abundant, respectively. HCA and PCA displayed four clusters of samples with a differentiated microbiome fingerprint. This work represents a pilot study into the viruses and bacteria that make up the human sperm microbiome. Despite the high variability observed, some patterns of similarity among individuals were identified. Further NGS studies under standardized methodological procedures are necessary to achieve a deep knowledge of the semen microbiome and its implications in male fertility.

A Window of Vulnerability: Chronic Environmental Stress Does Not Impair Reproduction in the Swordfish Xiphias gladius.

Xiphias gladius is an important fishing resource. The Mediterranean stock is affected by overfishing and is declining. In this light, the aim of this study was to evaluate the cross-talk among metabolism, stress response, immune system and reproduction in immature and mature females, coupling histological and transcriptomic approaches. The transcriptome of livers from 3 immature and 3 mature females was analyzed using the Artificial Intelligence RNA-Seq. For the histological analysis, ovary and liver samples were collected from 50 specimens caught during the reproductive season in the Mediterranean Sea. A total of 750 genes were differentially expressed between the livers. The gene ontologtabey analysis showed 91 upregulated and 161 downregulated biological process GO terms. Instead, the KEGG enrichment analysis revealed 15 enriched pathways. Furthermore, the binding occurring between estrogen receptors and aryl hydrocarbon receptor nuclear translocator, upregulated in mature females, could be liable for the inhibition of detoxification pathway. Indeed, at the histological level, mature females showed a higher density and number of melanomacrophage centers, biomarkers of stress. The present findings reveal the cross-talk among response to environmental stressors, metabolism and reproduction, highlighting that mature females invest a lot of energy in reproduction instead of immune response and detoxification.

The Chloroplast Genome of Endive (Cichorium endivia L.): Cultivar Structural Variants and Transcriptome Responses to Stress Due to Rain Extreme Events.

The chloroplast (cp) genome diversity has been used in phylogeny studies, breeding, and variety protection, and its expression has been shown to play a role in stress response. Smooth- and curly-leafed endives (Cichorium endivia var. latifolium and var. crispum) are of nutritional and economic importance and are the target of ever-changing breeding programmes. A reference cp genome sequence was assembled and annotated (cultivar 'Confiance'), which was 152,809 base pairs long, organized into the angiosperm-typical quadripartite structure, harboring two inverted repeats separated by the large- and short- single copy regions. The annotation included 136 genes, 90 protein-coding genes, 38 transfer, and 8 ribosomal RNAs and the sequence generated a distinct phyletic group within Asteraceae with the well-separated C. endivia and intybus species. SSR variants within the reference genome were mostly of tri-nucleotide type, and the cytosine to uracil (C/U) RNA editing recurred. The cp genome was nearly fully transcribed, hence sequence polymorphism was investigated by RNA-Seq of seven cultivars, and the SNP number was higher in smooth- than curly-leafed ones. All cultivars maintained C/U changes in identical positions, suggesting that RNA editing patterns were conserved; most cultivars shared SNPs of moderate impact on protein changes in the ndhD, ndhA, and psbF genes, suggesting that their variability may have a potential role in adaptive response. The cp transcriptome expression was investigated in leaves of plants affected by pre-harvest rainfall and rainfall excess plus waterlogging events characterized by production loss, compared to those of a cycle not affected by extreme rainfall. Overall, the analyses evidenced stress- and cultivar-specific responses, and further revealed that genes of the Cytochrome b6/f, and PSI-PSII systems were commonly affected and likely to be among major targets of extreme rain-related stress.