Hybridization and polyploidy enable genomic plasticity without sex in the most devastating plant-parasitic nematodes.

Root-knot nematodes (genus Meloidogyne) exhibit a diversity of reproductive modes ranging from obligatory sexual to fully asexual reproduction. Intriguingly, the most widespread and devastating species to global agriculture are those that reproduce asexually, without meiosis. To disentangle this surprising parasitic success despite the absence of sex and genetic exchanges, we have sequenced and assembled the genomes of three obligatory ameiotic and asexual Meloidogyne. We have compared them to those of relatives able to perform meiosis and sexual reproduction. We show that the genomes of ameiotic asexual Meloidogyne are large, polyploid and made of duplicated regions with a high within-species average nucleotide divergence of ~8%. Phylogenomic analysis of the genes present in these duplicated regions suggests that they originated from multiple hybridization events and are thus homoeologs. We found that up to 22% of homoeologous gene pairs were under positive selection and these genes covered a wide spectrum of predicted functional categories. To biologically assess functional divergence, we compared expression patterns of homoeologous gene pairs across developmental life stages using an RNAseq approach in the most economically important asexually-reproducing nematode. We showed that >60% of homoeologous gene pairs display diverged expression patterns. These results suggest a substantial functional impact of the genome structure. Contrasting with high within-species nuclear genome divergence, mitochondrial genome divergence between the three ameiotic asexuals was very low, signifying that these putative hybrids share a recent common maternal ancestor. Transposable elements (TE) cover a ~1.7 times higher proportion of the genomes of the ameiotic asexual Meloidogyne compared to the sexual relative and might also participate in their plasticity. The intriguing parasitic success of asexually-reproducing Meloidogyne species could be partly explained by their TE-rich composite genomes, resulting from allopolyploidization events, and promoting plasticity and functional divergence between gene copies in the absence of sex and meiosis.

Reference Transcriptomes and Detection of Duplicated Copies in Hexaploid and Allododecaploid Spartina Species (Poaceae).

In this study, we report the assembly and annotation of five reference transcriptomes for the European hexaploid Spartina species (S. maritima, S. alterniflora and their homoploid hybrids S. x townsendii and S. x neyrautii) and the allododecaploid invasive species S. anglica These transcriptomes were constructed from various leaf and root cDNA libraries that were sequenced using both Roche-454 and Illumina technologies. Considering the high ploidy levels of the Spartina genomes under study, and considering the absence of diploid reference genome and the need of an appropriate analytical strategy, we developed generic bioinformatics tools to (1) detect different haplotypes of each gene within each species and (2) assign a parental origin to haplotypes detected in the hexaploid hybrids and the neo-allopolyploid. The approach described here allows the detection of putative homeologs from sets of short reads. Synonymous substitution rate (KS) comparisons between haplotypes from the hexaploid species revealed the presence of one KS peak (likely resulting from the tetraploid duplication event). The procedure developed in this study can be applied for future differential gene expression or genomics experiments to study the fate of duplicated genes in the invasive allododecaploid S. anglica.

Identification of an Alternative Splicing Product of the Otx2 Gene Expressed in the Neural Retina and Retinal Pigmented Epithelial Cells.

To investigate the complexity of alternative splicing in the retina, we sequenced and analyzed a total of 115,706 clones from normalized cDNA libraries from mouse neural retina (66,217) and rat retinal pigmented epithelium (49,489). Based upon clustering the cDNAs and mapping them with their respective genomes, the estimated numbers of genes were 9,134 for the mouse neural retina and 12,050 for the rat retinal pigmented epithelium libraries. This unique collection of retinal of messenger RNAs is maintained and accessible through a web-base server to the whole community of retinal biologists for further functional characterization. The analysis revealed 3,248 and 3,202 alternative splice events for mouse neural retina and rat retinal pigmented epithelium, respectively. We focused on transcription factors involved in vision. Among the six candidates suitable for functional analysis, we selected Otx2S, a novel variant of the Otx2 gene with a deletion within the homeodomain sequence. Otx2S is expressed in both the neural retina and retinal pigmented epithelium, and encodes a protein that is targeted to the nucleus. OTX2S exerts transdominant activity on the tyrosinase promoter when tested in the physiological environment of primary RPE cells. By overexpressing OTX2S in primary RPE cells using an adeno associated viral vector, we identified 10 genes whose expression is positively regulated by OTX2S. We find that OTX2S is able to bind to the chromatin at the promoter of the retinal dehydrogenase 10 (RDH10) gene.

Dual RNA-seq reveals Meloidogyne graminicola transcriptome and candidate effectors during the interaction with rice plants.

Root-knot nematodes secrete proteinaceous effectors into plant tissues to facilitate infection by suppressing host defences and reprogramming the host metabolism to their benefit. Meloidogyne graminicola is a major pest of rice (Oryza sativa) in Asia and Latin America, causing important crop losses. The goal of this study was to identify M. graminicola pathogenicity genes expressed during the plant-nematode interaction. Using the dual RNA-sequencing (RNA-seq) strategy, we generated transcriptomic data of M. graminicola samples covering the pre-parasitic J2 stage and five parasitic stages in rice plants, from the parasitic J2 to the adult female. In the absence of a reference genome, a de novo M. graminicola transcriptome of 66 396 contigs was obtained from those reads that were not mapped on the rice genome. Gene expression profiling across the M. graminicola life cycle revealed key genes involved in nematode development and provided insights into the genes putatively associated with parasitism. The development of a 'secreted protein prediction' pipeline revealed a typical set of proteins secreted by nematodes, as well as a large number of cysteine-rich proteins and putative nuclear proteins. Combined with expression data, this pipeline enabled the identification of 15 putative effector genes, including two homologues of well-characterized effectors from cyst nematodes (CLE-like and VAP1) and a metallothionein. The localization of gene expression was assessed by in situ hybridization for a subset of candidates. All of these data represent important molecular resources for the elucidation of M. graminicola biology and for the selection of potential targets for the development of novel control strategies for this nematode species.

Convergent Evolution of Endosymbiont Differentiation in Dalbergioid and Inverted Repeat-Lacking Clade Legumes Mediated by Nodule-Specific Cysteine-Rich Peptides.

Nutritional symbiotic interactions require the housing of large numbers of microbial symbionts, which produce essential compounds for the growth of the host. In the legume-rhizobium nitrogen-fixing symbiosis, thousands of rhizobium microsymbionts, called bacteroids, are confined intracellularly within highly specialized symbiotic host cells. In Inverted Repeat-Lacking Clade (IRLC) legumes such as Medicago spp., the bacteroids are kept under control by an arsenal of nodule-specific cysteine-rich (NCR) peptides, which induce the bacteria in an irreversible, strongly elongated, and polyploid state. Here, we show that in Aeschynomene spp. legumes belonging to the more ancient Dalbergioid lineage, bacteroids are elongated or spherical depending on the Aeschynomene spp. and that these bacteroids are terminally differentiated and polyploid, similar to bacteroids in IRLC legumes. Transcriptome, in situ hybridization, and proteome analyses demonstrated that the symbiotic cells in the Aeschynomene spp. nodules produce a large diversity of NCR-like peptides, which are transported to the bacteroids. Blocking NCR transport by RNA interference-mediated inactivation of the secretory pathway inhibits bacteroid differentiation. Together, our results support the view that bacteroid differentiation in the Dalbergioid clade, which likely evolved independently from the bacteroid differentiation in the IRLC clade, is based on very similar mechanisms used by IRLC legumes.

Transcriptomic response of Manduca sexta immune tissues to parasitization by the bracovirus associated wasp Cotesia congregata.

During oviposition, Cotesia congregata parasitoid wasps inject into their host, Manduca sexta, some biological factors such as venom, ovarian fluid and a symbiotic polydnavirus (PDV) named Cotesia congregata bracovirus (CcBV). During parasitism, complex interactions occur between wasp-derived factors and host targets that lead to important modifications in host physiology. In particular, the immune response leading to wasp egg encapsulation is inhibited allowing wasp survival. To date, the regulation of host genes during the interaction had only been studied for a limited number of genes. In this study, we analysed the global impact of parasitism on host gene regulation 24 h post oviposition by high throughput 454 transcriptomic analyses of two tissues known to be involved in the host immune response (hemocytes and fat body). To identify specific effects of parasitism on host transcription at this time point, transcriptomes were obtained from non-treated and parasitized larvae, and also from larvae injected with heat-killed bacteria and double stimulated larvae that were parasitized prior to bacterial challenge. Results showed that, immune challenge by bacteria leads to induction of certain antimicrobial peptide (AMP) genes in M. sexta larvae whether they were parasitized or not prior to bacterial challenge. These results show that at 24 h post oviposition pathways leading to expression of AMP genes are not all inactivated suggesting wasps are in an antiseptic environment. In contrast, at this time point genes involved in phenoloxidase activation and cellular immune responses were globally down-regulated after parasitism in accordance with the observed inhibition of wasp egg encapsulation.

The coffee genome provides insight into the convergent evolution of caffeine biosynthesis.

Coffee is a valuable beverage crop due to its characteristic flavor, aroma, and the stimulating effects of caffeine. We generated a high-quality draft genome of the species Coffea canephora, which displays a conserved chromosomal gene order among asterid angiosperms. Although it shows no sign of the whole-genome triplication identified in Solanaceae species such as tomato, the genome includes several species-specific gene family expansions, among them N-methyltransferases (NMTs) involved in caffeine production, defense-related genes, and alkaloid and flavonoid enzymes involved in secondary compound synthesis. Comparative analyses of caffeine NMTs demonstrate that these genes expanded through sequential tandem duplications independently of genes from cacao and tea, suggesting that caffeine in eudicots is of polyphyletic origin.

Plant genetics. Early allopolyploid evolution in the post-Neolithic Brassica napus oilseed genome.

Oilseed rape (Brassica napus L.) was formed ~7500 years ago by hybridization between B. rapa and B. oleracea, followed by chromosome doubling, a process known as allopolyploidy. Together with more ancient polyploidizations, this conferred an aggregate 72× genome multiplication since the origin of angiosperms and high gene content. We examined the B. napus genome and the consequences of its recent duplication. The constituent An and Cn subgenomes are engaged in subtle structural, functional, and epigenetic cross-talk, with abundant homeologous exchanges. Incipient gene loss and expression divergence have begun. Selection in B. napus oilseed types has accelerated the loss of glucosinolate genes, while preserving expansion of oil biosynthesis genes. These processes provide insights into allopolyploid evolution and its relationship with crop domestication and improvement.

Functional annotation of Cotesia congregata bracovirus: identification of viral genes expressed in parasitized host immune tissues.

UNLABELLED: Bracoviruses (BVs) from the Polydnaviridae family are symbiotic viruses used as biological weapons by parasitoid wasps to manipulate lepidopteran host physiology and induce parasitism success. BV particles are produced by wasp ovaries and injected along with the eggs into the caterpillar host body, where viral gene expression is necessary for wasp development. Recent sequencing of the proviral genome of Cotesia congregata BV (CcBV) identified 222 predicted virulence genes present on 35 proviral segments integrated into the wasp genome. To date, the expressions of only a few selected candidate virulence genes have been studied in the caterpillar host, and we lacked a global vision of viral gene expression. In this study, a large-scale transcriptomic analysis by 454 sequencing of two immune tissues (fat body and hemocytes) of parasitized Manduca sexta caterpillar hosts allowed the detection of expression of 88 CcBV genes expressed 24 h after the onset of parasitism. We linked the expression profiles of these genes to several factors, showing that different regulatory mechanisms control viral gene expression in the host. These factors include the presence of signal peptides in encoded proteins, diversification of promoter regions, and, more surprisingly, gene position on the proviral genome. Indeed, most genes for which expression could be detected are localized in particular proviral regions globally producing higher numbers of circles. Moreover, this polydnavirus (PDV) transcriptomic analysis also reveals that a majority of CcBV genes possess at least one intron and an arthropod transcription start site, consistent with an insect origin of these virulence genes. IMPORTANCE: Bracoviruses (BVs) are symbiotic polydnaviruses used by parasitoid wasps to manipulate lepidopteran host physiology, ensuring wasp offspring survival. To date, the expressions of only a few selected candidate BV virulence genes have been studied in caterpillar hosts. We performed a large-scale analysis of BV gene expression in two immune tissues of Manduca sexta caterpillars parasitized by Cotesia congregata wasps. Genes for which expression could be detected corresponded to genes localized in particular regions of the viral genome globally producing higher numbers of circles. Our study thus brings an original global vision of viral gene expression and paves the way to the determination of the regulatory mechanisms enabling the expression of BV genes in targeted organisms, such as major insect pests. In addition, we identify sequence features suggesting that most BV virulence genes were acquired from insect genomes.

Soil metatranscriptomics for mining eukaryotic heavy metal resistance genes.

Heavy metals are pollutants which affect all organisms. Since a small number of eukaryotes have been investigated with respect to metal resistance, we hypothesize that many genes that control this phenomenon remain to be identified. This was tested by screening soil eukaryotic metatranscriptomes which encompass RNA from organisms belonging to the main eukaryotic phyla. Soil-extracted polyadenylated mRNAs were converted into cDNAs and 35 of them were selected for their ability to rescue the metal (Cd or Zn) sensitive phenotype of yeast mutants. Few of the genes belonged to families known to confer metal resistance when overexpressed in yeast. Several of them were homologous to genes that had not been studied in the context of metal resistance. For instance, the BOLA ones, which conferred cross metal (Zn, Co, Cd, Mn) resistance may act by interfering with Fe homeostasis. Other genes, such as those encoding 110- to 130-amino-acid-long, cysteine-rich polypeptides, had no homologues in databases. This study confirms that functional metatranscriptomics represents a powerful approach to address basic biological processes in eukaryotes. The selected genes can be used to probe new pathways involved in metal homeostasis and to manipulate the resistance level of selected organisms.

Genome structure and metabolic features in the red seaweed Chondrus crispus shed light on evolution of the Archaeplastida.

Red seaweeds are key components of coastal ecosystems and are economically important as food and as a source of gelling agents, but their genes and genomes have received little attention. Here we report the sequencing of the 105-Mbp genome of the florideophyte Chondrus crispus (Irish moss) and the annotation of the 9,606 genes. The genome features an unusual structure characterized by gene-dense regions surrounded by repeat-rich regions dominated by transposable elements. Despite its fairly large size, this genome shows features typical of compact genomes, e.g., on average only 0.3 introns per gene, short introns, low median distance between genes, small gene families, and no indication of large-scale genome duplication. The genome also gives insights into the metabolism of marine red algae and adaptations to the marine environment, including genes related to halogen metabolism, oxylipins, and multicellularity (microRNA processing and transcription factors). Particularly interesting are features related to carbohydrate metabolism, which include a minimalistic gene set for starch biosynthesis, the presence of cellulose synthases acquired before the primary endosymbiosis showing the polyphyly of cellulose synthesis in Archaeplastida, and cellulases absent in terrestrial plants as well as the occurrence of a mannosylglycerate synthase potentially originating from a marine bacterium. To explain the observations on genome structure and gene content, we propose an evolutionary scenario involving an ancestral red alga that was driven by early ecological forces to lose genes, introns, and intergenetic DNA; this loss was followed by an expansion of genome size as a consequence of activity of transposable elements.

454-pyrosequencing of Coffea arabica leaves infected by the rust fungus Hemileia vastatrix reveals in planta-expressed pathogen-secreted proteins and plant functions in a late compatible plant-rust interaction.

Coffee (Coffea arabica L.), one of the key export and cash crops in tropical and subtropical countries, suffers severe losses from the rust fungus Hemileia vastatrix. The transcriptome of H. vastatrix was analysed during a compatible interaction with coffee to obtain an exhaustive repertoire of the genes expressed during infection and to identify potential effector genes. Large-scale sequencing (454-GS-FLEX Titanium) of mixed coffee and rust cDNAs obtained from 21-day rust-infected leaves generated 352 146 sequences which assembled into 22 774 contigs. In the absence of any reference genomic sequences for Coffea or Hemileia, specific trinucleotide frequencies within expressed sequence tags (ESTs) and blast homology against a set of dicots and basidiomycete genomes were used to distinguish pathogen from plant sequences. About 30% (6763) of the contigs were assigned to H. vastatrix and 61% (13 951) to C. arabica. The majority (60%) of the rust sequences did not show homology to any genomic database, indicating that they were potential novel fungal genes. In silico analyses of the 6763 H. vastatrix contigs predicted 382 secreted proteins and identified homologues of the flax rust haustorially expressed secreted proteins (HESPs) and bean rust transferred protein 1 (RTP1). These rust candidate effectors showed conserved amino-acid domains and conserved patterns of cysteine positions suggestive of conserved functions during infection of host plants. Quantitative reverse transcription-polymerase chain reaction profiling of selected rust genes revealed dynamic expression patterns during the time course of infection of coffee leaves. This study provides the first valuable genomic resource for the agriculturally important plant pathogen H. vastatrix and the first comprehensive C. arabica EST dataset.

Stage-specific gene expression during spermatogenesis in the dogfish (Scyliorhinus canicula).

In the dogfish testis, the cystic arrangement and polarization of germ cell stages make it possible to observe all stages of spermatogenesis in a single transverse section. By taking advantage of the zonation of this organ, we have used suppressive subtractive libraries construction, real-time PCR, and in situ hybridization to identify 32 dogfish genes showing differential expressions during spermatogenesis. These include homologs of genes already known to be expressed in the vertebrate testis, but found here to be specifically expressed either in pre-meiotic and/or meiotic zones (ribosomal protein S8, high-mobility group box 3, ubiquitin carboxyl-terminal esterase L3, 20beta-hydroxysteroid dehydrogenase, or cyclophilin B) or in post-meiotic zone (speriolin, Soggy, zinc finger protein 474, calreticulin, or phospholipase c-zeta). We also report, for the first time, testis-specific expression patterns for dogfish genes coding for A-kinase anchor protein 5, ring finger protein 152, or F-box only protein 7. Finally, the study highlights the differential expression of new sequences whose identity remains to be assessed. This study provides the first molecular characterization of spermatogenesis in a chondrichthyan, a key species to gain insight into the evolution of this process in gnathostomes.

Convergent origins and rapid evolution of spliced leader trans-splicing in metazoa: insights from the ctenophora and hydrozoa.

Replacement of mRNA 5' UTR sequences by short sequences trans-spliced from specialized, noncoding, spliced leader (SL) RNAs is an enigmatic phenomenon, occurring in a set of distantly related animal groups including urochordates, nematodes, flatworms, and hydra, as well as in Euglenozoa and dinoflagellates. Whether SL trans-splicing has a common evolutionary origin and biological function among different organisms remains unclear. We have undertaken a systematic identification of SL exons in cDNA sequence data sets from non-bilaterian metazoan species and their closest unicellular relatives. SL exons were identified in ctenophores and in hydrozoan cnidarians, but not in other cnidarians, placozoans, or sponges, or in animal unicellular relatives. Mapping of SL absence/presence obtained from this and previous studies onto current phylogenetic trees favors an evolutionary scenario involving multiple origins for SLs during eumetazoan evolution rather than loss from a common ancestor. In both ctenophore and hydrozoan species, multiple SL sequences were identified, showing high sequence diversity. Detailed analysis of a large data set generated for the hydrozoan Clytia hemisphaerica revealed trans-splicing of given mRNAs by multiple alternative SLs. No evidence was found for a common identity of trans-spliced mRNAs between different hydrozoans. One feature found specifically to characterize SL-spliced mRNAs in hydrozoans, however, was a marked adenosine enrichment immediately 3' of the SL acceptor splice site. Our findings of high sequence divergence and apparently indiscriminate use of SLs in hydrozoans, along with recent findings in other taxa, indicate that SL genes have evolved rapidly in parallel in diverse animal groups, with constraint on SL exon sequence evolution being apparently rare.

Are Hox genes ancestrally involved in axial patterning? Evidence from the hydrozoan Clytia hemisphaerica (Cnidaria).

BACKGROUND: The early evolution and diversification of Hox-related genes in eumetazoans has been the subject of conflicting hypotheses concerning the evolutionary conservation of their role in axial patterning and the pre-bilaterian origin of the Hox and ParaHox clusters. The diversification of Hox/ParaHox genes clearly predates the origin of bilaterians. However, the existence of a "Hox code" predating the cnidarian-bilaterian ancestor and supporting the deep homology of axes is more controversial. This assumption was mainly based on the interpretation of Hox expression data from the sea anemone, but growing evidence from other cnidarian taxa puts into question this hypothesis. METHODOLOGY/PRINCIPAL FINDINGS: Hox, ParaHox and Hox-related genes have been investigated here by phylogenetic analysis and in situ hybridisation in Clytia hemisphaerica, an hydrozoan species with medusa and polyp stages alternating in the life cycle. Our phylogenetic analyses do not support an origin of ParaHox and Hox genes by duplication of an ancestral ProtoHox cluster, and reveal a diversification of the cnidarian HOX9-14 genes into three groups called A, B, C. Among the 7 examined genes, only those belonging to the HOX9-14 and the CDX groups exhibit a restricted expression along the oral-aboral axis during development and in the planula larva, while the others are expressed in very specialised areas at the medusa stage. CONCLUSIONS/SIGNIFICANCE: Cross species comparison reveals a strong variability of gene expression along the oral-aboral axis and during the life cycle among cnidarian lineages. The most parsimonious interpretation is that the Hox code, collinearity and conservative role along the antero-posterior axis are bilaterian innovations.