Tracing 100 million years of grass genome evolutionary plasticity.

Grasses derive from a family of monocotyledonous plants that includes crops of major economic importance such as wheat, rice, sorghum and barley, sharing a common ancestor some 100 million years ago. The genomic attributes of plant adaptation remain obscure and the consequences of recurrent whole genome duplications (WGD) or polyploidization events, a major force in plant evolution, remain largely speculative. We conducted a comparative analysis of omics data from ten grass species to unveil structural (inversions, fusions, fissions, duplications, substitutions) and regulatory (expression and methylation) basis of genome plasticity, as possible attributes of plant long lasting evolution and adaptation. The present study demonstrates that diverged polyploid lineages sharing a common WGD event often present the same patterns of structural changes and evolutionary dynamics, but these patterns are difficult to generalize across independent WGD events as a result of non-WGD factors such as selection and domestication of crops. Polyploidy is unequivocally linked to the evolutionary success of grasses during the past 100 million years, although it remains difficult to attribute this success to particular genomic consequences of polyploidization, suggesting that polyploids harness the potential of genome duplication, at least partially, in lineage-specific ways. Overall, the present study clearly demonstrates that post-polyploidization reprogramming is more complex than traditionally reported in investigating single species and calls for a critical and comprehensive comparison across independently polyploidized lineages.

A MITE insertion abolishes the AP3-3 self-maintenance regulatory loop in apetalous flowers of Nigella damascena.

MADS-box transcription factors are important regulators of floral organ identity through their binding to specific motifs, termed CArG, in the promoter of their target genes. Petal initiation and development depend on class A and B genes, but MADS-box genes of the APETALA3 (AP3) clade are key regulators of this process. In the early diverging eudicot Nigella damascena, an apetalous [T] morph is characterized by the lack of expression of the NdAP3-3 gene, with its expression being petal-specific in the wild-type [P] morph. All [T] morph plants are homozygous for an NdAP3-3 allele with a Miniature Inverted-repeat Transposable Element (MITE) insertion in the second intron of the gene. Here, we investigated to which extent the MITE insertion impairs regulation of the NdAP3-3 gene. We found that expression of NdAP3-3 is initiated in the [T] morph, but the MITE insertion prevents its positive self-maintenance by affecting the correct splicing of the mRNA. We also found specific CArG features in the promoter of the NdAP3-3 genes with petal-specific expression. However, they are not sufficient to drive expression only in petals of transgenic Arabidopsis, highlighting the existence of Nigella-specific cis/trans-acting factors in regulating AP3 paralogs.

A genome sequence resource for the genus Passiflora, the genome of the wild diploid species Passiflora organensis.

The genus Passiflora comprises a large group of plants popularly known as passionfruit, much appreciated for their exotic flowers and edible fruits. The species (∼500) are morphologically variable (e.g., growth habit, size, and color of flowers) and are adapted to distinct tropical ecosystems. In this study, we generated the genome of the wild diploid species Passiflora organensis Gardner by adopting a hybrid assembly approach. Passiflora organensis has a small genome of 259 Mbp and a heterozygosity rate of 81%, consistent with its reproductive system. Most of the genome sequences could be integrated into its chromosomes with cytogenomic markers (satellite DNA) as references. The repeated sequences accounted for 58.55% of the total DNA analyzed, and the Tekay lineage was the prevalent retrotransposon. In total, 25,327 coding genes were predicted. Passiflora organensis retains 5,609 singletons and 15,671 gene families. We focused on the genes potentially involved in the locus determining self-incompatibility and the MADS-box gene family, allowing us to infer expansions and contractions within specific subfamilies. Finally, we recovered the organellar DNA. Structural rearrangements and two mitoviruses, besides relics of other mobile elements, were found in the chloroplast and mt-DNA molecules, respectively. This study presents the first draft genome assembly of a wild Passiflora species, providing a valuable sequence resource for genomic and evolutionary studies on the genus, and support for breeding cropped passionfruit species.

A wheat cysteine-rich receptor-like kinase confers broad-spectrum resistance against Septoria tritici blotch.

The poverty of disease resistance gene reservoirs limits the breeding of crops for durable resistance against evolutionary dynamic pathogens. Zymoseptoria tritici which causes Septoria tritici blotch (STB), represents one of the most genetically diverse and devastating wheat pathogens worldwide. No fully virulent Z. tritici isolates against synthetic wheats carrying the major resistant gene Stb16q have been identified. Here, we use comparative genomics, mutagenesis and complementation to identify Stb16q, which confers broad-spectrum resistance against Z. tritici. The Stb16q gene encodes a plasma membrane cysteine-rich receptor-like kinase that was recently introduced into cultivated wheat and which considerably slows penetration and intercellular growth of the pathogen.

Corrigendum: FANCM Limits Meiotic Crossovers in Brassica Crops.

[This corrects the article DOI: 10.3389/fpls.2018.00368.].

The wild grape genome sequence provides insights into the transition from dioecy to hermaphroditism during grape domestication.

BACKGROUND: A key step in domestication of the grapevine was the transition from separate sexes (dioecy) in wild Vitis vinifera ssp. sylvestris (V. sylvestris) to hermaphroditism in cultivated Vitis vinifera ssp. sativa (V. vinifera). It is known that V. sylvestris has an XY system and V. vinifera a modified Y haplotype (Yh) and that the sex locus is small, but it has not previously been precisely characterized. RESULTS: We generate a high-quality de novo reference genome for V. sylvestris, onto which we map whole-genome re-sequencing data of a cross to locate the sex locus. Assembly of the full X, Y, and Yh haplotypes of V. sylvestris and V. vinifera sex locus and examining their gene content and expression profiles during flower development in wild and cultivated accessions show that truncation and deletion of tapetum and pollen development genes on the X haplotype likely causes male sterility, while the upregulation of a Y allele of a cytokinin regulator (APRT3) may cause female sterility. The downregulation of this cytokinin regulator in the Yh haplotype may be sufficient to trigger reversal to hermaphroditism. Molecular dating of X and Y haplotypes is consistent with the sex locus being as old as the Vitis genus, but the mechanism by which recombination was suppressed remains undetermined. CONCLUSIONS: We describe the genomic and evolutionary characterization of the sex locus of cultivated and wild grapevine, providing a coherent model of sex determination in the latter and for transition from dioecy to hermaphroditism during domestication.

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.

A Repertory of Rearrangements and the Loss of an Inverted Repeat Region in Passiflora Chloroplast Genomes.

Chloroplast genomes (cpDNA) in angiosperms are usually highly conserved. Although rearrangements have been observed in some lineages, such as Passiflora, the mechanisms that lead to rearrangements are still poorly elucidated. In the present study, we obtained 20 new chloroplast genomes (18 species from the genus Passiflora, and Dilkea retusa and Mitostemma brevifilis from the family Passifloraceae) in order to investigate cpDNA evolutionary history in this group. Passiflora cpDNAs vary in size considerably, with ∼50 kb between shortest and longest. Large inverted repeat (IR) expansions were identified, and at the extreme opposite, the loss of an IR was detected for the first time in Passiflora, a rare event in angiosperms. The loss of an IR region was detected in Passiflora capsularis and Passiflora costaricensis, a species in which occasional biparental chloroplast inheritance has previously been reported. A repertory of rearrangements such as inversions and gene losses were detected, making Passiflora one of the few groups with complex chloroplast genome evolution. We also performed a phylogenomic study based on all the available cp genomes and our analysis implies that there is a need to reconsider the taxonomic classifications of some species in the group.

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.

Influence of CNV on transcript levels of HvCBF genes at Fr-H2 locus revealed by resequencing in resistant barley cv. 'Nure' and expression analysis.

Resequencing in resistant cultivar 'Nure' and structural comparison with the same region of susceptible 'Morex' was performed in order to gain a better insight into barley Frost-resistance-H2 locus. Accurate annotation showed copy number variation (CNV) in the proximal part of the locus. In 'Nure', two exact copies of the HvCBF4-HvCBF2A region and one of the HvCBF4-HvCBF2B segment were observed, while in 'Morex' the corresponding region harboured a single HvCBF4-HvCBF2A (22 kb) segment. Abundance and diversity of repetitive element classes, gene function gain/losses, regulatory motifs and SNPs in gene sequences were identified. An expression study of key HvCBFs with/without CNV on selected genotypes contrasting for frost resistance and estimated HvCBF4-HvCBF2B copy number (2-10 copies) was also performed. Under light stimulus at warm temperature (23 °C), CNV of HvCBF2A and HvCBF4 correlated with their expression levels and reported frost resistance of genotypes; moreover, expression levels of HvCBF2A and HvCBF14 were strongly correlated (r = 0.908, p < 0.01). On the other hand, frost resistance correlated to HvCBF14 expression (r = 0.871, p < 0.01) only after cold induction (6°C) in the dark. A complex interplay of HvCBFs expression levels under different light/temperature stimuli is discussed in light of CNV and presence/number of regulatory elements that integrate different signal transduction pathways.

Convergent evolution of the UbiA prenyltransferase family underlies the independent acquisition of furanocoumarins in plants.

Furanocoumarins (FCs) are plant-specialized metabolites with potent allelochemical properties. The distribution of FCs is scattered with a chemotaxonomical tendency towards four distant families with highly similar FC pathways. The mechanism by which this pathway emerged and spread in plants has not been elucidated. Furanocoumarin biosynthesis was investigated in Ficus carica (fig, Moraceae), focusing on the first committed reaction catalysed by an umbelliferone dimethylallyltransferase (UDT). Comparative RNA-seq analysis among latexes of different fig organs led to the identification of a UDT. The phylogenetic relationship of this UDT to previously reported Apiaceae UDTs was evaluated. The expression pattern of F. carica prenyltransferase 1 (FcPT1) was related to the FC contents in different latexes. Enzymatic characterization demonstrated that one of the main functions of FcPT1 is UDT activity. Phylogenetic analysis suggested that FcPT1 and Apiaceae UDTs are derived from distinct ancestors, although they both belong to the UbiA superfamily. These findings are supported by significant differences in the related gene structures. This report describes the identification of FcPT1 involved in FC biosynthesis in fig and provides new insights into multiple origins of the FC pathway and, more broadly, into the adaptation of plants to their environments.

A receptor-like kinase enhances sunflower resistance to Orobanche cumana.

Orobanche cumana (sunflower broomrape) is an obligate parasitic plant that infects sunflower roots, causing yield losses. Here, by using a map-based cloning strategy, we identified HaOr7-a gene that confers resistance to O. cumana race F-which was found to encode a leucine-rich repeat receptor-like kinase. The complete HAOR7 protein is present in resistant lines of sunflower and prevents O. cumana from connecting to the vascular system of sunflower roots, whereas susceptible lines encode a truncated protein that lacks transmembrane and kinase domains.

Transposable element discovery and characterization of LTR-retrotransposon evolutionary lineages in the tropical fruit species Passiflora edulis.

A significant proportion of plant genomes is consists of transposable elements (TEs), especially LTR retrotransposons (LTR-RTs) which are known to drive genome evolution. However, not much information is available on the structure and evolutionary role of TEs in the Passifloraceae family (Malpighiales order). Against this backdrop, we identified, characterized, and inferred the potential genomic impact of the TE repertoire found in the available genomic resources for Passiflora edulis, a tropical fruit species. A total of 250 different TE sequences were identified (96% Class I, and 4% Class II), corresponding to ~ 19% of the P. edulis draft genome. TEs were found preferentially in intergenic spaces (70.4%), but also overlapping genes (30.6%). LTR-RTs accounted for 181 single elements corresponding to ~ 13% of the draft genome. A phylogenetic inference of the reverse transcriptase domain of the LTR-RT revealed association of 37 elements with the Copia superfamily (Angela, Ale, Tork, and Sire) and 128 with the Gypsy (Del, Athila, Reina, CRM, and Galadriel) superfamily, and Del elements were the most frequent. Interestingly, according to insertion time analysis, the majority (95.9%) of the LTR-RTs were recently inserted into the P. edulis genome (< 2.0 Mya), and with the exception of the Athila lineage, all LTR-RTs are transcriptionally active. Moreover, functional analyses disclosed that the Angela, Del, CRM and Tork lineages are conserved in wild Passiflora species, supporting the idea of a common expansion of Copia and Gypsy superfamilies. Overall, this is the first study describing the P. edulis TE repertoire, and it also lends weight to the suggestion that LTR-RTs had a recent expansion into the analyzed gene-rich region of the P. edulis genome, possibly along WGD (Whole genome duplication) events, but are under negative selection due to their potential deleterious impact on gene regions.

µLAS technology for DNA isolation coupled to Cas9-assisted targeting for sequencing and assembly of a 30 kb region in plant genome.

Cas9-assisted targeting of DNA fragments in complex genomes is viewed as an essential strategy to obtain high-quality and continuous sequence data. However, the purity of target loci selected by pulsed-field gel electrophoresis (PFGE) has so far been insufficient to assemble the sequence in one contig. Here, we describe the µLAS technology to capture and purify high molecular weight DNA. First, the technology is optimized to perform high sensitivity DNA profiling with a limit of detection of 20 fg/µl for 50 kb fragments and an analytical time of 50 min. Then, µLAS is operated to isolate a 31.5 kb locus cleaved by Cas9 in the genome of the plant Medicago truncatula. Target purification is validated on a Bacterial Artificial Chromosome plasmid, and subsequently carried out in whole genome with µLAS, PFGE or by combining these techniques. PacBio sequencing shows an enrichment factor of the target sequence of 84 with PFGE alone versus 892 by association of PFGE with µLAS. These performances allow us to sequence and assemble one contig of 29 441 bp with 99% sequence identity to the reference sequence.

Genetic basis and timing of a major mating system shift in Capsella.

A crucial step in the transition from outcrossing to self-fertilization is the loss of genetic self-incompatibility (SI). In the Brassicaceae, SI involves the interaction of female and male specificity components, encoded by the genes SRK and SCR at the self-incompatibility locus (S-locus). Theory predicts that S-linked mutations, and especially dominant mutations in SCR, are likely to contribute to loss of SI. However, few studies have investigated the contribution of dominant mutations to loss of SI in wild plant species. Here, we investigate the genetic basis of loss of SI in the self-fertilizing crucifer species Capsella orientalis, by combining genetic mapping, long-read sequencing of complete S-haplotypes, gene expression analyses and controlled crosses. We show that loss of SI in C. orientalis occurred < 2.6 Mya and maps as a dominant trait to the S-locus. We identify a fixed frameshift deletion in the male specificity gene SCR and confirm loss of male SI specificity. We further identify an S-linked small RNA that is predicted to cause dominance of self-compatibility. Our results agree with predictions on the contribution of dominant S-linked mutations to loss of SI, and thus provide new insights into the molecular basis of mating system transitions.

Gene Duplication in the Sugarcane Genome: A Case Study of Allele Interactions and Evolutionary Patterns in Two Genic Regions.

Sugarcane (Saccharum spp.) is highly polyploid and aneuploid. Modern cultivars are derived from hybridization between S. officinarum and S. spontaneum. This combination results in a genome exhibiting variable ploidy among different loci, a huge genome size (~10 Gb) and a high content of repetitive regions. An approach using genomic, transcriptomic, and genetic mapping can improve our knowledge of the behavior of genetics in sugarcane. The hypothetical HP600 and Centromere Protein C (CENP-C) genes from sugarcane were used to elucidate the allelic expression and genomic and genetic behaviors of this complex polyploid. The physically linked side-by-side genes HP600 and CENP-C were found in two different homeologous chromosome groups with ploidies of eight and ten. The first region (Region01) was a Sorghum bicolor ortholog region with all haplotypes of HP600 and CENP-C expressed, but HP600 exhibited an unbalanced haplotype expression. The second region (Region02) was a scrambled sugarcane sequence formed from different noncollinear genes containing partial duplications of HP600 and CENP-C (paralogs). This duplication resulted in a non-expressed HP600 pseudogene and a recombined fusion version of CENP-C and the orthologous gene Sobic.003G299500 with at least two chimeric gene haplotypes expressed. It was also determined that it occurred before Saccharum genus formation and after the separation of sorghum and sugarcane. A linkage map was constructed using markers from nonduplicated Region01 and for the duplication (Region01 and Region02). We compare the physical and linkage maps, demonstrating the possibility of mapping markers located in duplicated regions with markers in nonduplicated region. Our results contribute directly to the improvement of linkage mapping in complex polyploids and improve the integration of physical and genetic data for sugarcane breeding programs. Thus, we describe the complexity involved in sugarcane genetics and genomics and allelic dynamics, which can be useful for understanding complex polyploid genomes.

Development of a Sequence-Based Reference Physical Map of Pea (Pisum sativum L.).

Whole genome profiling (WGP) is a sequence-based physical mapping technology and uses sequence tags generated by next generation sequencing for construction of bacterial artificial chromosome (BAC) contigs of complex genomes. The physical map provides a framework for assembly of genome sequence and information for localization of genes that are difficult to find through positional cloning. To address the challenges of accurate assembly of the pea genome (∼4.2 GB of which approximately 85% is repetitive sequences), we have adopted the WGP technology for assembly of a pea BAC library. Multi-dimensional pooling of 295,680 BAC clones and sequencing the ends of restriction fragments of pooled DNA generated 1,814 million high quality reads, of which 825 million were deconvolutable to 1.11 million unique WGP sequence tags. These WGP tags were used to assemble 220,013 BACs into contigs. Assembly of the BAC clones using the modified Fingerprinted Contigs (FPC) program has resulted in 13,040 contigs, consisting of 213,719 BACs, and 6,294 singleton BACs. The average contig size is 0.33 Mbp and the N50 contig size is 0.62 Mbp. WGPTM technology has proved to provide a robust physical map of the pea genome, which would have been difficult to assemble using traditional restriction digestion based methods. This sequence-based physical map will be useful to assemble the genome sequence of pea. Additionally, the 1.1 million WGP tags will support efficient assignment of sequence scaffolds to the BAC clones, and thus an efficient sequencing of BAC pools with targeted genome regions of interest.

Whole-genome landscape of Medicago truncatula symbiotic genes.

Advances in deciphering the functional architecture of eukaryotic genomes have been facilitated by recent breakthroughs in sequencing technologies, enabling a more comprehensive representation of genes and repeat elements in genome sequence assemblies, as well as more sensitive and tissue-specific analyses of gene expression. Here we show that PacBio sequencing has led to a substantially improved genome assembly of Medicago truncatula A17, a legume model species notable for endosymbiosis studies1, and has enabled the identification of genome rearrangements between genotypes at a near-base-pair resolution. Annotation of the new M. truncatula genome sequence has allowed for a thorough analysis of transposable elements and their dynamics, as well as the identification of new players involved in symbiotic nodule development, in particular 1,037 upregulated long non-coding RNAs (lncRNAs). We have also discovered that a substantial proportion (~35% and 38%, respectively) of the genes upregulated in nodules or expressed in the nodule differentiation zone colocalize in genomic clusters (270 and 211, respectively), here termed symbiotic islands. These islands contain numerous expressed lncRNA genes and display differentially both DNA methylation and histone marks. Epigenetic regulations and lncRNAs are therefore attractive candidate elements for the orchestration of symbiotic gene expression in the M. truncatula genome.

A gene-rich fraction analysis of the Passiflora edulis genome reveals highly conserved microsyntenic regions with two related Malpighiales species.

Passiflora edulis is the most widely cultivated species of passionflowers, cropped mainly for industrialized juice production and fresh fruit consumption. Despite its commercial importance, little is known about the genome structure of P. edulis. To fill in this gap in our knowledge, a genomic library was built, and now completely sequenced over 100 large-inserts. Sequencing data were assembled from long sequence reads, and structural sequence annotation resulted in the prediction of about 1,900 genes, providing data for subsequent functional analysis. The richness of repetitive elements was also evaluated. Microsyntenic regions of P. edulis common to Populus trichocarpa and Manihot esculenta, two related Malpighiales species with available fully sequenced genomes were examined. Overall, gene order was well conserved, with some disruptions of collinearity identified as rearrangements, such as inversion and translocation events. The microsynteny level observed between the P. edulis sequences and the compared genomes is surprising, given the long divergence time that separates them from the common ancestor. P. edulis gene-rich segments are more compact than those of the other two species, even though its genome is much larger. This study provides a first accurate gene set for P. edulis, opening the way for new studies on the evolutionary issues in Malpighiales genomes.

The Genomic Basis of Color Pattern Polymorphism in the Harlequin Ladybird.

Many animal species comprise discrete phenotypic forms. A common example in natural populations of insects is the occurrence of different color patterns, which has motivated a rich body of ecological and genetic research [1-6]. The occurrence of dark, i.e., melanic, forms displaying discrete color patterns is found across multiple taxa, but the underlying genomic basis remains poorly characterized. In numerous ladybird species (Coccinellidae), the spatial arrangement of black and red patches on adult elytra varies wildly within species, forming strikingly different complex color patterns [7, 8]. In the harlequin ladybird, Harmonia axyridis, more than 200 distinct color forms have been described, which classic genetic studies suggest result from allelic variation at a single, unknown, locus [9, 10]. Here, we combined whole-genome sequencing, population-based genome-wide association studies, gene expression, and functional analyses to establish that the transcription factor Pannier controls melanic pattern polymorphism in H. axyridis. We show that pannier is necessary for the formation of melanic elements on the elytra. Allelic variation in pannier leads to protein expression in distinct domains on the elytra and thus determines the distinct color patterns in H. axyridis. Recombination between pannier alleles may be reduced by a highly divergent sequence of ∼170 kb in the cis-regulatory regions of pannier, with a 50 kb inversion between color forms. This most likely helps maintain the distinct alleles found in natural populations. Thus, we propose that highly variable discrete color forms can arise in natural populations through cis-regulatory allelic variation of a single gene.