Geographic population structure of the honeybee microsporidian parasite Vairimorpha (Nosema) ceranae in the South West Indian Ocean.

The microsporidian Vairimorpha (Nosema) ceranae is one of the most common parasites of the honeybee. A single honeybee carries many parasites and therefore multiple alleles of V. ceranae genes that seem to be ubiquitous. As a consequence, nucleotide diversity analyses have not allowed discriminating genetic structure of parasite populations. We performed deep loci-targeted sequencing to monitor the haplotype frequencies of genome markers in isolates from discontinuous territories, namely the tropical islands of the South West Indian Ocean. The haplotype frequency distribution corroborated the suspected tetraploidy of the parasite. Most major haplotypes were ubiquitous in the area but with variable frequency. While oceanic isolates differed from European and Asian outgroups, parasite populations from distinct archipelagoes also differed in their haplotype distribution. Interestingly an original and very divergent Malagasy isolate was detected. The observed population structure allowed formulating hypotheses upon the natural history of V. ceranae in this oceanic area. We also discussed the usefulness of allelic distribution assessment, using multiple informative loci or genome-wide analyses, when parasite population is not clonal within a single host.

The compact genome of the sponge Oopsacas minuta (Hexactinellida) is lacking key metazoan core genes.

BACKGROUND: Explaining the emergence of the hallmarks of bilaterians is a central focus of evolutionary developmental biology-evodevo-and evolutionary genomics. For this purpose, we must both expand and also refine our knowledge of non-bilaterian genomes, especially by studying early branching animals, in particular those in the metazoan phylum Porifera. RESULTS: We present a comprehensive analysis of the first whole genome of a glass sponge, Oopsacas minuta, a member of the Hexactinellida. Studying this class of sponge is evolutionary relevant because it differs from the three other Porifera classes in terms of development, tissue organization, ecology, and physiology. Although O. minuta does not exhibit drastic body simplifications, its genome is among the smallest of animal genomes sequenced so far, and surprisingly lacks several metazoan core genes (including Wnt and several key transcription factors). Our study also provides the complete genome of a symbiotic Archaea dominating the associated microbial community: a new Thaumarchaeota species. CONCLUSIONS: The genome of the glass sponge O. minuta differs from all other available sponge genomes by its compactness and smaller number of encoded proteins. The unexpected loss of numerous genes previously considered ancestral and pivotal for metazoan morphogenetic processes most likely reflects the peculiar syncytial tissue organization in this group. Our work further documents the importance of convergence during animal evolution, with multiple convergent evolution of septate-like junctions, electrical-signaling and multiciliated cells in metazoans.

First Vanilla planifolia High-Density Genetic Linkage Map Provides Quantitative Trait Loci for Resistance to Fusarium oxysporum.

Fusarium oxysporum f. sp. radicis-vanillae (Forv), the causal agent of root and stem rot disease, is the main pathogen affecting vanilla production. Sources of resistance have been reported in Vanilla planifolia G. Jackson ex Andrews, the main cultivated vanilla species. In this study, we developed the first high-density genetic map in this species with 1,804 genotyping-by-sequencing (GBS)-generated single nucleotide polymorphism (SNP) markers using 125 selfed progenies of the CR0040 traditional vanilla cultivar. Sixteen linkage groups (LG) were successfully constructed, with a mean of 113 SNPs and an average length of 207 cM per LG. The map had a high density with an average of 5.45 SNP every 10 cM and an average distance of 1.85 cM between adjacent markers. The first three LG were aligned against the first assembled chromosome of CR0040, and the other 13 LG were correctly associated with the other 13 assembled chromosomes. The population was challenged with the highly pathogenic Forv strain Fo072 using the root-dip inoculation method. Five traits were mapped, and 20 QTLs were associated with resistance to Fo072. Among the genes retrieved in the CR0040 physical regions associated with QTLs, genes potentially involved in biotic resistance mechanisms, coding for kinases, E3 ubiquitin ligases, pentatricopeptide repeat-containing proteins, and one leucine-rich repeat receptor underlying the qFo72_08.1 QTL have been highlighted. This study should provide useful resources for marker-assisted selection in V. planifolia.

A chromosome-level, haplotype-phased Vanilla planifolia genome highlights the challenge of partial endoreplication for accurate whole-genome assembly.

Vanilla planifolia, the species cultivated to produce one of the world's most popular flavors, is highly prone to partial genome endoreplication, which leads to highly unbalanced DNA content in cells. We report here the first molecular evidence of partial endoreplication at the chromosome scale by the assembly and annotation of an accurate haplotype-phased genome of V. planifolia. Cytogenetic data demonstrated that the diploid genome size is 4.09 Gb, with 16 chromosome pairs, although aneuploid cells are frequently observed. Using PacBio HiFi and optical mapping, we assembled and phased a diploid genome of 3.4 Gb with a scaffold N50 of 1.2 Mb and 59 128 predicted protein-coding genes. The atypical k-mer frequencies and the uneven sequencing depth observed agreed with our expectation of unbalanced genome representation. Sixty-seven percent of the genes were scattered over only 30% of the genome, putatively linking gene-rich regions and the endoreplication phenomenon. By contrast, low-coverage regions (non-endoreplicated) were rich in repeated elements but also contained 33% of the annotated genes. Furthermore, this assembly showed distinct haplotype-specific sequencing depth variation patterns, suggesting complex molecular regulation of endoreplication along the chromosomes. This high-quality, anchored assembly represents 83% of the estimated V. planifolia genome. It provides a significant step toward the elucidation of this complex genome. To support post-genomics efforts, we developed the Vanilla Genome Hub, a user-friendly integrated web portal that enables centralized access to high-throughput genomic and other omics data and interoperable use of bioinformatics tools.

Genotyping-by-Sequencing Technology in Plant Taxonomy and Phylogeny.

Genotyping-by-sequencing (GBS) is a method to discover and genotype simultaneous genome-wide high-throughput single nucleotide polymorphisms (SNPs). GBS is based on reducing genome complexity with restriction enzymes. Here we describe a method developed by Elshire et al. for constructing simplified GBS libraries and recent bioinformatic approaches developed to analyze the large volume of polymorphism data generated by this method. GBS approach is suitable for population studies, taxonomic and phylogenic studies, germplasm characterization, and breeding and trait mapping for a wide range of organisms, including plants with complex genomes.

High-throughput sequencing of complete genomes of ipomoviruses associated with an epidemic of cassava brown streak disease in the Comoros Archipelago.

Using high-throughput sequencing of small interfering RNAs (siRNAs), virion-associated nucleic acid (VANA), and double stranded RNAs (dsRNAs), we have determined the complete genome sequences of Comorian isolates of two ipomoviruses, cassava brown streak virus (CBSV) and a divergent isolate of Ugandan cassava brown streak virus (UCBSV-KM) representing a new strain of this virus. While the large ORF of CBSV shares the highest nucleotide sequence identity (95.9%) with a Tanzanian isolate of CBSV, the large UCBSV-KM ORF shares the highest nucleotide sequence identity (77.5%) with a Malawian isolate of UCBSV. This low value is near the species demarcation threshold for the family Potyviridae (<76%). Phylogenetic analysis confirms that UCBSV-KM represents a new lineage that is genetically distinct from the currently described UCBSV strains.

New genomic data and analyses challenge the traditional vision of animal epithelium evolution.

BACKGROUND: The emergence of epithelia was the foundation of metazoan expansion. Epithelial tissues are a hallmark of metazoans deeply rooted in the evolution of their complex developmental morphogenesis processes. However, studies on the epithelial features of non-bilaterians are still sparse and it remains unclear whether the last common metazoan ancestor possessed a fully functional epithelial toolkit or if it was acquired later during metazoan evolution. RESULTS: To investigate the early evolution of animal epithelia, we sequenced the genome and transcriptomes of two new sponge species to characterize epithelial markers such as the E-cadherin complex and the polarity complexes for all classes (Calcarea, Demospongiae, Hexactinellida, Homoscleromorpha) of sponges (phylum Porifera) and compare them with their homologues in Placozoa and in Ctenophora. We found that Placozoa and most sponges possess orthologues of all essential genes encoding proteins characteristic of bilaterian epithelial cells, as well as their conserved interaction domains. In stark contrast, we found that ctenophores lack several major polarity complex components such as the Crumbs complex and Scribble. Furthermore, the E-cadherin ctenophore orthologue exhibits a divergent cytoplasmic domain making it unlikely to interact with its canonical cytoplasmic partners. CONCLUSIONS: These unexpected findings challenge the current evolutionary paradigm on the emergence of epithelia. Altogether, our results raise doubt on the homology of protein complexes and structures involved in cell polarity and adhesive-type junctions between Ctenophora and Bilateria epithelia.

Genotyping by Sequencing Highlights a Polygenic Resistance to Ralstonia pseudosolanacearum in Eggplant (Solanum melongena L.).

Eggplant cultivation is limited by numerous diseases, including the devastating bacterial wilt (BW) caused by the Ralstonia solanacearum species complex (RSSC). Within the RSSC, Ralstonia pseudosolanacearum (including phylotypes I and III) causes severe damage to all solanaceous crops, including eggplant. Therefore, the creation of cultivars resistant to R. pseudosolanacearum strains is a major goal for breeders. An intraspecific eggplant population, segregating for resistance, was created from the cross between the susceptible MM738 and the resistant EG203 lines. The population of 123 doubled haploid lines was challenged with two strains belonging to phylotypes I (PSS4) and III (R3598), which both bypass the published EBWR9 BW-resistance quantitative trait locus (QTL). Ten and three QTLs of resistance to PSS4 and to R3598, respectively, were detected and mapped. All were strongly influenced by environmental conditions. The most stable QTLs were found on chromosomes 3 and 6. Given their estimated physical position, these newly detected QTLs are putatively syntenic with BW-resistance QTLs in tomato. In particular, the QTLs' position on chromosome 6 overlaps with that of the major broad-spectrum tomato resistance QTL Bwr-6. The present study is a first step towards understanding the complex polygenic system, which underlies the high level of BW resistance of the EG203 line.

Eggplant Resistance to the Ralstonia solanacearum Species Complex Involves Both Broad-Spectrum and Strain-Specific Quantitative Trait Loci.

Bacterial wilt (BW) is a major disease of solanaceous crops caused by the Ralstonia solanacearum species complex (RSSC). Strains are grouped into five phylotypes (I, IIA, IIB, III, and IV). Varietal resistance is the most sustainable strategy for managing BW. Nevertheless, breeding to improve cultivar resistance has been limited by the pathogen's extensive genetic diversity. Identifying the genetic bases of specific and non-specific resistance is a prerequisite to breed improvement. A major gene (ERs1) was previously mapped in eggplant (Solanum melongena L.) using an intraspecific population of recombinant inbred lines derived from the cross of susceptible MM738 (S) × resistant AG91-25 (R). ERs1 was originally found to control three strains from phylotype I, while being totally ineffective against a virulent strain from the same phylotype. We tested this population against four additional RSSC strains, representing phylotypes I, IIA, IIB, and III in order to clarify the action spectrum of ERs1. We recorded wilting symptoms and bacterial stem colonization under controlled artificial inoculation. We constructed a high-density genetic map of the population using single nucleotide polymorphisms (SNPs) developed from genotyping-by-sequencing and added 168 molecular markers [amplified fragment length polymorphisms (AFLPs), simple sequence repeats (SSRs), and sequence-related amplified polymorphisms (SRAPs)] developed previously. The new linkage map based on a total of 1,035 markers was anchored on eggplant, tomato, and potato genomes. Quantitative trait locus (QTL) mapping for resistance against a total of eight RSSC strains resulted in the detection of one major phylotype-specific QTL and two broad-spectrum QTLs. The major QTL, which specifically controls three phylotype I strains, was located at the bottom of chromosome 9 and corresponded to the previously identified major gene ERs1. Five candidate R-genes were underlying this QTL, with different alleles between the parents. The two other QTLs detected on chromosomes 2 and 5 were found to be associated with partial resistance to strains of phylotypes I, IIA, III and strains of phylotypes IIA and III, respectively. Markers closely linked to these three QTLs will be crucial for breeding eggplant with broad-spectrum resistance to BW. Furthermore, our study provides an important contribution to the molecular characterization of ERs1, which was initially considered to be a major resistance gene.

Evolutionary forces affecting synonymous variations in plant genomes.

Base composition is highly variable among and within plant genomes, especially at third codon positions, ranging from GC-poor and homogeneous species to GC-rich and highly heterogeneous ones (particularly Monocots). Consequently, synonymous codon usage is biased in most species, even when base composition is relatively homogeneous. The causes of these variations are still under debate, with three main forces being possibly involved: mutational bias, selection and GC-biased gene conversion (gBGC). So far, both selection and gBGC have been detected in some species but how their relative strength varies among and within species remains unclear. Population genetics approaches allow to jointly estimating the intensity of selection, gBGC and mutational bias. We extended a recently developed method and applied it to a large population genomic dataset based on transcriptome sequencing of 11 angiosperm species spread across the phylogeny. We found that at synonymous positions, base composition is far from mutation-drift equilibrium in most genomes and that gBGC is a widespread and stronger process than selection. gBGC could strongly contribute to base composition variation among plant species, implying that it should be taken into account in plant genome analyses, especially for GC-rich ones.

A large set of 26 new reference transcriptomes dedicated to comparative population genomics in crops and wild relatives.

We produced a unique large data set of reference transcriptomes to obtain new knowledge about the evolution of plant genomes and crop domestication. For this purpose, we validated a RNA-Seq data assembly protocol to perform comparative population genomics. For the validation, we assessed and compared the quality of de novo Illumina short-read assemblies using data from two crops for which an annotated reference genome was available, namely grapevine and sorghum. We used the same protocol for the release of 26 new transcriptomes of crop plants and wild relatives, including still understudied crops such as yam, pearl millet and fonio. The species list has a wide taxonomic representation with the inclusion of 15 monocots and 11 eudicots. All contigs were annotated using BLAST, prot4EST and Blast2GO. A strong originality of the data set is that each crop is associated with close relative species, which will permit whole-genome comparative evolutionary studies between crops and their wild-related species. This large resource will thus serve research communities working on both crops and model organisms. All the data are available at http://arcad-bioinformatics.southgreen.fr/.

Lineage-Specific Evolutionary Histories and Regulation of Major Starch Metabolism Genes during Banana Ripening.

Starch is the most widespread and abundant storage carbohydrate in plants. It is also a major feature of cultivated bananas as it accumulates to large amounts during banana fruit development before almost complete conversion to soluble sugars during ripening. Little is known about the structure of major gene families involved in banana starch metabolism and their evolution compared to other species. To identify genes involved in banana starch metabolism and investigate their evolutionary history, we analyzed six gene families playing a crucial role in plant starch biosynthesis and degradation: the ADP-glucose pyrophosphorylases (AGPases), starch synthases (SS), starch branching enzymes (SBE), debranching enzymes (DBE), α-amylases (AMY) and ß-amylases (BAM). Using comparative genomics and phylogenetic approaches, these genes were classified into families and sub-families and orthology relationships with functional genes in Eudicots and in grasses were identified. In addition to known ancestral duplications shaping starch metabolism gene families, independent evolution in banana and grasses also occurred through lineage-specific whole genome duplications for specific sub-families of AGPase, SS, SBE, and BAM genes; and through gene-scale duplications for AMY genes. In particular, banana lineage duplications yielded a set of AGPase, SBE and BAM genes that were highly or specifically expressed in banana fruits. Gene expression analysis highlighted a complex transcriptional reprogramming of starch metabolism genes during ripening of banana fruits. A differential regulation of expression between banana gene duplicates was identified for SBE and BAM genes, suggesting that part of starch metabolism regulation in the fruit evolved in the banana lineage.

Draft Genome Sequence of an Alphaproteobacterium Associated with the Mediterranean Sponge Oscarella lobularis.

While sequencing DNA purified from the homoscleromorph sponge Oscarella lobularis, we detected a large number of reads with strong similarity to available alphaproteobacteria gene sequences of family Rhodobacteraceae. Here, we present the genome sequence of this putative sponge symbiont that we propose to designate as "Candidatus Rhodobacter lobularis."

Mitochondrial Genome Sequence of the Glass Sponge Oopsacas minuta.

We report the complete mitochondrial genome sequence of the Mediterranean glass sponge Oopsacas minuta. This 19-kb mitochondrial genome has 24 noncoding genes (22 tRNAs and 2 rRNAs) and 14 protein-encoding genes coding for 11 subunits of respiratory chain complexes and 3 ATP synthase subunits.

Expansion of banana (Musa acuminata) gene families involved in ethylene biosynthesis and signalling after lineage-specific whole-genome duplications.

Whole-genome duplications (WGDs) are widespread in plants, and three lineage-specific WGDs occurred in the banana (Musa acuminata) genome. Here, we analysed the impact of WGDs on the evolution of banana gene families involved in ethylene biosynthesis and signalling, a key pathway for banana fruit ripening. Banana ethylene pathway genes were identified using comparative genomics approaches and their duplication modes and expression profiles were analysed. Seven out of 10 banana ethylene gene families evolved through WGD and four of them (1-aminocyclopropane-1-carboxylate synthase (ACS), ethylene-insensitive 3-like (EIL), ethylene-insensitive 3-binding F-box (EBF) and ethylene response factor (ERF)) were preferentially retained. Banana orthologues of AtEIN3 and AtEIL1, two major genes for ethylene signalling in Arabidopsis, were particularly expanded. This expansion was paralleled by that of EBF genes which are responsible for control of EIL protein levels. Gene expression profiles in banana fruits suggested functional redundancy for several MaEBF and MaEIL genes derived from WGD and subfunctionalization for some of them. We propose that EIL and EBF genes were co-retained after WGD in banana to maintain balanced control of EIL protein levels and thus avoid detrimental effects of constitutive ethylene signalling. In the course of evolution, subfunctionalization was favoured to promote finer control of ethylene signalling.

Two evolutionarily distinct classes of paleopolyploidy.

Whole genome duplications (WGDs) occurred in the distant evolutionary history of many lineages and are particularly frequent in the flowering plant lineages. Following paleopolyploidization in plants, most duplicated genes are deleted by intrachromosomal recombination, a process referred to as fractionation. In the examples studied so far, genes are disproportionately lost from one of the parental subgenomes (biased fractionation) and the subgenome having lost the lowest number of genes is more expressed (genome dominance). In the present study, we analyzed the pattern of gene deletion and gene expression following the most recent WGD in banana (alpha event) and extended our analyses to seven other sequenced plant genomes: poplar, soybean, medicago, arabidopsis, sorghum, brassica, and maize. We propose a new class of ancient WGD, with Musa (alpha), poplar, and soybean as members, where genes are both deleted and expressed to an equal extent (unbiased fractionation and genome equivalence). We suggest that WGDs with genome dominance and biased fractionation (Class I) may result from ancient allotetraploidies, while WGDs without genome dominance or biased fractionation (Class II) may result from ancient autotetraploidies.

The banana (Musa acuminata) genome and the evolution of monocotyledonous plants.

Bananas (Musa spp.), including dessert and cooking types, are giant perennial monocotyledonous herbs of the order Zingiberales, a sister group to the well-studied Poales, which include cereals. Bananas are vital for food security in many tropical and subtropical countries and the most popular fruit in industrialized countries. The Musa domestication process started some 7,000 years ago in Southeast Asia. It involved hybridizations between diverse species and subspecies, fostered by human migrations, and selection of diploid and triploid seedless, parthenocarpic hybrids thereafter widely dispersed by vegetative propagation. Half of the current production relies on somaclones derived from a single triploid genotype (Cavendish). Pests and diseases have gradually become adapted, representing an imminent danger for global banana production. Here we describe the draft sequence of the 523-megabase genome of a Musa acuminata doubled-haploid genotype, providing a crucial stepping-stone for genetic improvement of banana. We detected three rounds of whole-genome duplications in the Musa lineage, independently of those previously described in the Poales lineage and the one we detected in the Arecales lineage. This first monocotyledon high-continuity whole-genome sequence reported outside Poales represents an essential bridge for comparative genome analysis in plants. As such, it clarifies commelinid-monocotyledon phylogenetic relationships, reveals Poaceae-specific features and has led to the discovery of conserved non-coding sequences predating monocotyledon-eudicotyledon divergence.

Characterization of a cinnamoyl-CoA reductase 1 (CCR1) mutant in maize: effects on lignification, fibre development, and global gene expression.

Cinnamoyl-CoA reductase (CCR), which catalyses the first committed step of the lignin-specific branch of monolignol biosynthesis, has been extensively characterized in dicot species, but few data are available in monocots. By screening a Mu insertional mutant collection in maize, a mutant in the CCR1 gene was isolated named Zmccr1(-). In this mutant, CCR1 gene expression is reduced to 31% of the residual wild-type level. Zmccr1(-) exhibited enhanced digestibility without compromising plant growth and development. Lignin analysis revealed a slight decrease in lignin content and significant changes in lignin structure. p-Hydroxyphenyl units were strongly decreased and the syringyl/guaiacyl ratio was slightly increased. At the cellular level, alterations in lignin deposition were mainly observed in the walls of the sclerenchymatic fibre cells surrounding the vascular bundles. These cell walls showed little to no staining with phloroglucinol. These histochemical changes were accompanied by an increase in sclerenchyma surface area and an alteration in cell shape. In keeping with this cell type-specific phenotype, transcriptomics performed at an early stage of plant development revealed the down-regulation of genes specifically associated with fibre wall formation. To the present authors' knowledge, this is the first functional characterization of CCR1 in a grass species.

Evolution and expression of class III peroxidases.

Class III peroxidases are members of a large multigenic family, only detected in the plant kingdom and absent from green algae sensu stricto (chlorophyte algae or Chlorophyta). Their evolution is thought to be related to the emergence of the land plants. However class III peroxidases are present in a lower copy number in some basal Streptophytes (Charapyceae), which predate land colonization. Gene structures are variable among organisms and within species with respect to the number of introns, but their positions are highly conserved. Their high copy number, as well as their conservation could be related to plant complexity and adaptation to increasing stresses. No specific function has been assigned to respective isoforms, but in large multigenic families, particular structure-function relations can be expected. Plant peroxidase sequences contain highly conserved residues and motifs, variable domains surrounded by conserved residues and present a low identity level among their promoter regions, further suggesting the existence of sub-functionalization of the different isoforms.