Identification of new potential molecular actors related to fiber quality in flax through Omics.

One of the biggest challenges for a more widespread utilization of plant fibers is to better understand the different molecular factors underlying the variability in fineness and mechanical properties of both elementary and scutched fibers. Accordingly, we analyzed genome-wide transcription profiling from bast fiber bearing tissues of seven different flax varieties (4 spring, 2 winter fiber varieties and 1 winter linseed) and identified 1041 differentially expressed genes between varieties, of which 97 were related to cell wall metabolism. KEGG analysis highlighted a number of different enriched pathways. Subsequent statistical analysis using Partial Least-Squares Discriminant Analysis showed that 73% of the total variance was explained by the first 3 X-variates corresponding to 56 differentially expressed genes. Calculation of Pearson correlations identified 5 genes showing a strong correlation between expression and morphometric data. Two-dimensional gel proteomic analysis on the two varieties showing the most discriminant and significant differences in morphometrics revealed 1490 protein spots of which 108 showed significant differential abundance. Mass spectrometry analysis successfully identified 46 proteins representing 32 non-redundant proteins. Statistical clusterization based on the expression level of genes corresponding to the 32 proteins showed clear discrimination into three separate clusters, reflecting the variety type (spring-/winter-fiber/oil). Four of the 32 proteins were also highly correlated with morphometric features. Examination of predicted functions for the 9 (5 + 4) identified genes highlighted lipid metabolism and senescence process. Calculation of Pearson correlation coefficients between expression data and retted fiber mechanical measurements (strength and maximum force) identified 3 significantly correlated genes. The genes were predicted to be connected to cell wall dynamics, either directly (Expansin-like protein), or indirectly (NAD(P)-binding Rossmann-fold superfamily protein). Taken together, our results have allowed the identification of molecular actors potentially associated with the determination of both in-planta fiber morphometrics, as well as ex-planta fiber mechanical properties, both of which are key parameters for elementary fiber and scutched fiber quality in flax.

An improved assembly of the pearl millet reference genome using Oxford Nanopore long reads and optical mapping.

Pearl millet (Pennisetum glaucum (L.)) R. Br. syn. Cenchrus americanus (L.) Morrone) is an important crop in South Asia and sub-Saharan Africa which contributes to ensuring food security. Its genome has an estimated size of 1.76 Gb and displays a high level of repetitiveness above 80%. A first assembly was previously obtained for the Tift 23D2B1-P1-P5 cultivar genotype using short-read sequencing technologies. This assembly is, however, incomplete and fragmented with around 200 Mb unplaced on chromosomes. We report here an improved quality assembly of the pearl millet Tift 23D2B1-P1-P5 cultivar genotype obtained with an approach combining Oxford Nanopore long reads and Bionano Genomics optical maps. This strategy allowed us to add around 200 Mb at the chromosome-level assembly. Moreover, we strongly improved continuity in the order of the contigs and scaffolds within the chromosomes, particularly in the centromeric regions. Notably, we added more than 100 Mb around the centromeric region on chromosome 7. This new assembly also displayed a higher gene completeness with a complete BUSCO score of 98.4% using the Poales database. This more complete and higher quality assembly of the Tift 23D2B1-P1-P5 genotype now available to the community will help in the development of research on the role of structural variants and more broadly in genomics studies and the breeding of pearl millet.

Long-read and chromosome-scale assembly of the hexaploid wheat genome achieves high resolution for research and breeding.

BACKGROUND: The sequencing of the wheat (Triticum aestivum) genome has been a methodological challenge for many years owing to its large size (15.5 Gb), repeat content, and hexaploidy. Many initiatives aiming at obtaining a reference genome of cultivar Chinese Spring have been launched in the past years and it was achieved in 2018 as the result of a huge effort to combine short-read sequencing with many other resources. Reference-quality genome assemblies were then produced for other accessions, but the rapid evolution of sequencing technologies offers opportunities to reach high-quality standards at lower cost. RESULTS: Here, we report on an optimized procedure based on long reads produced on the Oxford Nanopore Technology PromethION device to assemble the genome of the French bread wheat cultivar Renan. CONCLUSIONS: We provide the most contiguous chromosome-scale assembly of a bread wheat genome to date. Coupled with an annotation based on RNA-sequencing data, this resource will be valuable for the crop community and will facilitate the rapid selection of agronomically important traits. We also provide a framework to generate high-quality assemblies of complex genomes using ONT.

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.

Population genomics of apricots unravels domestication history and adaptive events.

Among crop fruit trees, the apricot (Prunus armeniaca) provides an excellent model to study divergence and adaptation processes. Here, we obtain nearly 600 Armeniaca apricot genomes and four high-quality assemblies anchored on genetic maps. Chinese and European apricots form two differentiated gene pools with high genetic diversity, resulting from independent domestication events from distinct wild Central Asian populations, and with subsequent gene flow. A relatively low proportion of the genome is affected by selection. Different genomic regions show footprints of selection in European and Chinese cultivated apricots, despite convergent phenotypic traits, with predicted functions in both groups involved in the perennial life cycle, fruit quality and disease resistance. Selection footprints appear more abundant in European apricots, with a hotspot on chromosome 4, while admixture is more pervasive in Chinese cultivated apricots. Our study provides clues to the biology of selected traits and targets for fruit tree research and breeding.

Fonio millet genome unlocks African orphan crop diversity for agriculture in a changing climate.

Sustainable food production in the context of climate change necessitates diversification of agriculture and a more efficient utilization of plant genetic resources. Fonio millet (Digitaria exilis) is an orphan African cereal crop with a great potential for dryland agriculture. Here, we establish high-quality genomic resources to facilitate fonio improvement through molecular breeding. These include a chromosome-scale reference assembly and deep re-sequencing of 183 cultivated and wild Digitaria accessions, enabling insights into genetic diversity, population structure, and domestication. Fonio diversity is shaped by climatic, geographic, and ethnolinguistic factors. Two genes associated with seed size and shattering showed signatures of selection. Most known domestication genes from other cereal models however have not experienced strong selection in fonio, providing direct targets to rapidly improve this crop for agriculture in hot and dry environments.

High-quality genome sequence of white lupin provides insight into soil exploration and seed quality.

White lupin (Lupinus albus L.) is an annual crop cultivated for its protein-rich seeds. It is adapted to poor soils due to the production of cluster roots, which are made of dozens of determinate lateral roots that drastically improve soil exploration and nutrient acquisition (mostly phosphate). Using long-read sequencing technologies, we provide a high-quality genome sequence of a cultivated accession of white lupin (2n = 50, 451 Mb), as well as de novo assemblies of a landrace and a wild relative. We describe a modern accession displaying increased soil exploration capacity through early establishment of lateral and cluster roots. We also show how seed quality may have been impacted by domestication in term of protein profiles and alkaloid content. The availability of a high-quality genome assembly together with companion genomic and transcriptomic resources will enable the development of modern breeding strategies to increase and stabilize white lupin yield.

Tube-Wise Diagnostic Microarray for the Multiplex Characterization of the Complex Plant Pathogen Ralstonia solanacearum.

Ralstonia solanacearum is a well-known agricultural and ecological threat worldwide. The complexity of the R. solanacearum species complex (Rssc) represents a challenge for the accurate characterization of epidemiological strains by official services and research laboratories. The majority of protocols only focus on a narrow range of strains; however, this species complex includes strains that represent major constraints and are under strict regulation. The main drawback associated with the current methods of detecting and characterizing Rssc strains is their reliance on combining different protocols to properly characterize the strains at the ecotype level, which require time and money. Therefore, we used microarray technology (ArrayTube) to develop a standard protocol, which characterizes 17 major groups of interest in the Rssc, in a single multiplex reaction. These 17 majors groups are linked with a phylogenetic assignation (phylotypes, sequevars), but also with an ecotype assignation associated with a range of hosts (e.g., brown rot, Moko). Probes were designed with a 50-mer length constraint and thoroughly evaluated for any flaws or secondary structures. The strains are characterized based on a DNA extraction from pure culture. Validation data showed strong intra-repeatability, inter-repeatability, and reproducibility as well as good specificity. A hierarchical analysis of the probe groups is suitable for an accurate characterization. Compared with single marker detection tests, the method described in this paper addresses efficiently the issue of combining several tests by testing a large number of phylogenetic markers in a single reaction assay. This custom microarray (RsscAT) represents a significant improvement in the epidemiological monitoring of Rssc strains worldwide, and it has the potential to provide insights for phylogenetic incongruence of Rssc strains based on the host of isolation and may be used to indicate potentially emergent strains.

Molecular Epidemiology of Bacterial Wilt in the Madagascar Highlands Caused by Andean (Phylotype IIB-1) and African (Phylotype III) Brown Rot Strains of the Ralstonia solanacearum Species Complex.

The Ralstonia solanacearum species complex (RSSC) is a highly diverse cluster of bacterial strains found worldwide, many of which are destructive and cause bacterial wilt (BW) in a wide range of host plants. In 2009, potato production in Madagascar was dramatically affected by several BW epidemics. Controlling this disease is critical for Malagasy potato producers. The first important step toward control is the characterization of strains and their putative origins. The genetic diversity and population structure of the RSSC were investigated in the major potato production areas of the Highlands. A large collection of strains (n = 1224) was assigned to RSSC phylotypes based on multiplex polymerase chain reaction (PCR). Phylotypes I and III have been present in Madagascar for a long time but rarely associated with major potato BW outbreaks. The marked increase of BW prevalence was found associated with phylotype IIB sequevar 1 (IIB-1) strains (n = 879). This is the first report of phylotype IIB-1 strains in Madagascar. In addition to reference strains, epidemic IIB-1 strains (n = 255) were genotyped using the existing MultiLocus Variable-Number Tandem Repeat Analysis (MLVA) scheme RS2-MLVA9, producing 31 haplotypes separated into two related clonal complexes (CCs). One major CC included most of the worldwide haplotypes distributed across wide areas. A regional-scale investigation suggested that phylotype IIB-1 strains were introduced and massively spread via latently infected potato seed tubers. Additionally, the genetic structure of phylotype IIB-1 likely resulted from a bottleneck/founder effect. The population structure of phylotype III, described here for the first time in Madagascar, exhibited a different pattern. Phylotype III strains (n = 217) were genotyped using the highly discriminatory MLVA scheme RS3-MLVA16. High genetic diversity was uncovered, with 117 haplotypes grouped into 11 CCs. Malagasy phylotype III strains were highly differentiated from continental African strains, suggesting no recent migration from the continent. Overall, population structure of phylotype III involves individual small CCs that correlate to restricted geographic areas in Madagascar. The evidence suggests, if at all, that African phylotype III strains are not efficiently transmitted through latently infected potato seed tubers.

Ectopic lignification in the flax lignified bast fiber1 mutant stem is associated with tissue-specific modifications in gene expression and cell wall composition.

Histochemical screening of a flax ethyl methanesulfonate population led to the identification of 93 independent M2 mutant families showing ectopic lignification in the secondary cell wall of stem bast fibers. We named this core collection the Linum usitatissimum (flax) lbf mutants for lignified bast fibers and believe that this population represents a novel biological resource for investigating how bast fiber plants regulate lignin biosynthesis. As a proof of concept, we characterized the lbf1 mutant and showed that the lignin content increased by 350% in outer stem tissues containing bast fibers but was unchanged in inner stem tissues containing xylem. Chemical and NMR analyses indicated that bast fiber ectopic lignin was highly condensed and rich in G-units. Liquid chromatography-mass spectrometry profiling showed large modifications in the oligolignol pool of lbf1 inner- and outer-stem tissues that could be related to ectopic lignification. Immunological and chemical analyses revealed that lbf1 mutants also showed changes to other cell wall polymers. Whole-genome transcriptomics suggested that ectopic lignification of flax bast fibers could be caused by increased transcript accumulation of (1) the cinnamoyl-CoA reductase, cinnamyl alcohol dehydrogenase, and caffeic acid O-methyltransferase monolignol biosynthesis genes, (2) several lignin-associated peroxidase genes, and (3) genes coding for respiratory burst oxidase homolog NADPH-oxidases necessary to increase H2O2 supply.

Bridging physiological and evolutionary time-scales in a gene regulatory network.

Gene regulatory networks (GRNs) govern phenotypic adaptations and reflect the trade-offs between physiological responses and evolutionary adaptation that act at different time-scales. To identify patterns of molecular function and genetic diversity in GRNs, we studied the drought response of the common sunflower, Helianthus annuus, and how the underlying GRN is related to its evolution. We examined the responses of 32,423 expressed sequences to drought and to abscisic acid (ABA) and selected 145 co-expressed transcripts. We characterized their regulatory relationships in nine kinetic studies based on different hormones. From this, we inferred a GRN by meta-analyses of a Gaussian graphical model and a random forest algorithm and studied the genetic differentiation among populations (FST ) at nodes. We identified two main hubs in the network that transport nitrate in guard cells. This suggests that nitrate transport is a critical aspect of the sunflower physiological response to drought. We observed that differentiation of the network genes in elite sunflower cultivars is correlated with their position and connectivity. This systems biology approach combined molecular data at different time-scales and identified important physiological processes. At the evolutionary level, we propose that network topology could influence responses to human selection and possibly adaptation to dry environments.

Is gene transcription involved in seed dry after-ripening?

Orthodox seeds are living organisms that survive anhydrobiosis and may display dormancy, an inability to germinate at harvest. Seed germination potential can be acquired during a prolonged period of dry storage called after-ripening. The aim of this work was to determine if gene transcription is an underlying regulatory mechanism for dormancy alleviation during after-ripening. To identify changes in gene transcription strictly associated with the acquisition of germination potential but not with storage, we used seed storage at low relative humidity that maintains dormancy as control. Transcriptome profiling was performed using DNA microarray to compare change in gene transcript abundance between dormant (D), after-ripened non-dormant (ND) and after-ripened dormant seeds (control, C). Quantitative real-time polymerase chain reaction (qPCR) was used to confirm gene expression. Comparison between D and ND showed the differential expression of 115 probesets at cut-off values of two-fold change (p<0.05). Comparisons between both D and C with ND in transcript abundance showed that only 13 transcripts, among 115, could be specific to dormancy alleviation. qPCR confirms the expression pattern of these transcripts but without significant variation between conditions. Here we show that sunflower seed dormancy alleviation in the dry state is not related to regulated changes in gene expression.

A gene-phenotype network based on genetic variability for drought responses reveals key physiological processes in controlled and natural environments.

Identifying the connections between molecular and physiological processes underlying the diversity of drought stress responses in plants is key for basic and applied science. Drought stress response involves a large number of molecular pathways and subsequent physiological processes. Therefore, it constitutes an archetypical systems biology model. We first inferred a gene-phenotype network exploiting differences in drought responses of eight sunflower (Helianthus annuus) genotypes to two drought stress scenarios. Large transcriptomic data were obtained with the sunflower Affymetrix microarray, comprising 32423 probesets, and were associated to nine morpho-physiological traits (integrated transpired water, leaf transpiration rate, osmotic potential, relative water content, leaf mass per area, carbon isotope discrimination, plant height, number of leaves and collar diameter) using sPLS regression. Overall, we could associate the expression patterns of 1263 probesets to six phenotypic traits and identify if correlations were due to treatment, genotype and/or their interaction. We also identified genes whose expression is affected at moderate and/or intense drought stress together with genes whose expression variation could explain phenotypic and drought tolerance variability among our genetic material. We then used the network model to study phenotypic changes in less tractable agronomical conditions, i.e. sunflower hybrids subjected to different watering regimes in field trials. Mapping this new dataset in the gene-phenotype network allowed us to identify genes whose expression was robustly affected by water deprivation in both controlled and field conditions. The enrichment in genes correlated to relative water content and osmotic potential provides evidence of the importance of these traits in agronomical conditions.

Targeted mRNA oxidation regulates sunflower seed dormancy alleviation during dry after-ripening.

After-ripening is the mechanism by which dormant seeds become nondormant during their dry storage after harvest. The absence of free water in mature seeds does not allow detectable metabolism; thus, the processes associated with dormancy release under these conditions are largely unknown. We show here that sunflower (Helianthus annuus) seed alleviation of dormancy during after-ripening is associated with mRNA oxidation and that this oxidation is prevented when seeds are maintained dormant. In vitro approaches demonstrate that mRNA oxidation results in artifacts in cDNA-amplified fragment length polymorphim analysis and alters protein translation. The oxidation of transcripts is not random but selective, and, using microarrays, we identified 24 stored mRNAs that became highly oxidized during after-ripening. Oxidized transcripts mainly correspond to genes involved in responses to stress and in cell signaling. Among them, protein phosphatase 2C PPH1, mitogen-activated protein kinase phosphatase 1, and phenyl ammonia lyase 1 were identified. We propose that targeted mRNA oxidation during dry after-ripening of dormant seeds could be a process that governs cell signaling toward germination in the early steps of seed imbibition.