Genomic architecture of endogenous ichnoviruses reveals distinct evolutionary pathways leading to virus domestication in parasitic wasps.

BACKGROUND: Polydnaviruses (PDVs) are mutualistic endogenous viruses inoculated by some lineages of parasitoid wasps into their hosts, where they facilitate successful wasp development. PDVs include the ichnoviruses and bracoviruses that originate from independent viral acquisitions in ichneumonid and braconid wasps respectively. PDV genomes are fully incorporated into the wasp genomes and consist of (1) genes involved in viral particle production, which derive from the viral ancestor and are not encapsidated, and (2) proviral segments harboring virulence genes, which are packaged into the viral particle. To help elucidating the mechanisms that have facilitated viral domestication in ichneumonid wasps, we analyzed the structure of the viral insertions by sequencing the whole genome of two ichnovirus-carrying wasp species, Hyposoter didymator and Campoletis sonorensis. RESULTS: Assemblies with long scaffold sizes allowed us to unravel the organization of the endogenous ichnovirus and revealed considerable dispersion of the viral loci within the wasp genomes. Proviral segments contained species-specific sets of genes and occupied distinct genomic locations in the two ichneumonid wasps. In contrast, viral machinery genes were organized in clusters showing highly conserved gene content and order, with some loci located in collinear wasp genomic regions. This genomic architecture clearly differs from the organization of PDVs in braconid wasps, in which proviral segments are clustered and viral machinery elements are more dispersed. CONCLUSIONS: The contrasting structures of the two types of ichnovirus genomic elements are consistent with their different functions: proviral segments are vehicles for virulence proteins expected to adapt according to different host defense systems, whereas the genes involved in virus particle production in the wasp are likely more stable and may reflect ancestral viral architecture. The distinct genomic architectures seen in ichnoviruses versus bracoviruses reveal different evolutionary trajectories that have led to virus domestication in the two wasp lineages.

Cotesia congregata Bracovirus Circles Encoding PTP and Ankyrin Genes Integrate into the DNA of Parasitized Manduca sexta Hemocytes.

Polydnaviruses (PDVs) are essential for the parasitism success of tens of thousands of species of parasitoid wasps. PDVs are present in wasp genomes as proviruses, which serve as the template for the production of double-stranded circular viral DNA carrying virulence genes that are injected into lepidopteran hosts. PDV circles do not contain genes coding for particle production, thereby impeding viral replication in caterpillar hosts during parasitism. Here, we investigated the fate of PDV circles of Cotesia congregata bracovirus during parasitism of the tobacco hornworm, Manduca sexta, by the wasp Cotesia congregata Sequences sharing similarities with host integration motifs (HIMs) of Microplitis demolitor bracovirus (MdBV) circles involved in integration into DNA could be identified in 12 CcBV circles, which encode PTP and VANK gene families involved in host immune disruption. A PCR approach performed on a subset of these circles indicated that they persisted in parasitized M. sexta hemocytes as linear forms, possibly integrated in host DNA. Furthermore, by using a primer extension capture method based on these HIMs and high-throughput sequencing, we could show that 8 out of 9 circles tested were integrated in M. sexta hemocyte genomic DNA and that integration had occurred specifically using the HIM, indicating that an HIM-mediated specific mechanism was involved in their integration. Investigation of BV circle insertion sites at the genome scale revealed that certain genomic regions appeared to be enriched in BV insertions, but no specific M. sexta target site could be identified.IMPORTANCE The identification of a specific and efficient integration mechanism shared by several bracovirus species opens the question of its role in braconid parasitoid wasp parasitism success. Indeed, results obtained here show massive integration of bracovirus DNA in somatic immune cells at each parasitism event of a caterpillar host. Given that bracoviruses do not replicate in infected cells, integration of viral sequences in host DNA might allow the production of PTP and VANK virulence proteins within newly dividing cells of caterpillar hosts that continue to develop during parasitism. Furthermore, this integration process could serve as a basis to understand how PDVs mediate the recently identified gene flux between parasitoid wasps and Lepidoptera and the frequency of these horizontal transfer events in nature.

Genetic footprints of adaptive divergence in the bracovirus of Cotesia sesamiae identified by targeted resequencing.

The African parasitoid wasp Cotesia sesamiae is a generalist species structured in locally adapted populations showing differences in host range. The recent discovery of Cotesia typhae, a specialist, sister species to C. sesamiae, provides a good framework to study the genetic determinants of parasitoid host range. To investigate the genomic bases of divergence between these populations and species, we used a targeted sequencing approach on 24 samples. We targeted the bracovirus genomic region encoding virulence genes involved in the interaction with the lepidopteran hosts of the wasps. High sequencing coverage was obtained for all samples, allowing the study of genetic variation between wasp populations and species. By combining population genetic estimations, such as nucleotide diversity (π), relative differentiation (FST ) and absolute divergence (dxy ), with branch-site dN/dS measures, we identified six of 98 bracovirus genes showing significant divergence and evidence of positive selection. These genes, belonging to different gene families, are potentially involved in host adaptation and in the specialization process. Fine-scale analyses of genetic variation also revealed mutations and large deletions in certain genes inducing pseudogenization and loss of function. The image emerging from these results is that adaptation mediated by bracovirus genes happens through selection of particularly adaptive alleles and loss of nonadaptive genes. These results highlight the central role of the bracovirus in the molecular interactions between the wasps and their hosts and in the evolutionary processes of specialization.

Deciphering the Theobroma cacao self-incompatibility system: from genomics to diagnostic markers for self-compatibility.

Cocoa self-compatibility is an important yield factor and has been described as being controlled by a late gameto-sporophytic system expressed only at the level of the embryo sac. It results in gametic non-fusion and involves several loci. In this work, we identified two loci, located on chromosomes 1 and 4 (CH1 and CH4), involved in cocoa self-incompatibility by two different processes. Both loci are responsible for gametic selection, but only one (the CH4 locus) is involved in the main fruit drop. The CH1 locus acts prior to the gamete fusion step and independently of the CH4 locus. Using fine-mapping and genome-wide association studies, we focused analyses on restricted regions and identified candidate genes. Some of them showed a differential expression between incompatible and compatible reactions. Immunolocalization experiments provided evidence of CH1 candidate genes expressed in ovule and style tissues. Highly polymorphic simple sequence repeat (SSR) diagnostic markers were designed in the CH4 region that had been identified by fine-mapping. They are characterized by a strong linkage disequilibrium with incompatibility alleles, thus allowing the development of efficient diagnostic markers predicting self-compatibility and fruit setting according to the presence of specific alleles or genotypes. SSR alleles specific to self-compatible Amelonado and Criollo varieties were also identified, thus allowing screening for self-compatible plants in cocoa populations.

Mutations in ACTRT1 and its enhancer RNA elements lead to aberrant activation of Hedgehog signaling in inherited and sporadic basal cell carcinomas.

Basal cell carcinoma (BCC), the most common human cancer, results from aberrant activation of the Hedgehog signaling pathway. Although most cases of BCC are sporadic, some forms are inherited, such as Bazex-Dupré-Christol syndrome (BDCS)-a cancer-prone genodermatosis with an X-linked, dominant inheritance pattern. We have identified mutations in the ACTRT1 gene, which encodes actin-related protein T1 (ARP-T1), in two of the six families with BDCS that were examined in this study. High-throughput sequencing in the four remaining families identified germline mutations in noncoding sequences surrounding ACTRT1. These mutations were located in transcribed sequences encoding enhancer RNAs (eRNAs) and were shown to impair enhancer activity and ACTRT1 expression. ARP-T1 was found to directly bind to the GLI1 promoter, thus inhibiting GLI1 expression, and loss of ARP-T1 led to activation of the Hedgehog pathway in individuals with BDCS. Moreover, exogenous expression of ACTRT1 reduced the in vitro and in vivo proliferation rates of cell lines with aberrant activation of the Hedgehog signaling pathway. In summary, our study identifies a disease mechanism in BCC involving mutations in regulatory noncoding elements and uncovers the tumor-suppressor properties of ACTRT1.

Microbial Degradation of a Recalcitrant Pesticide: Chlordecone.

Chlordecone (Kepone®) is a synthetic organochlorine insecticide (C10Cl10O) used worldwide mostly during the 1970 and 1980s. Its intensive application in the French West Indies to control the banana black weevil Cosmopolites sordidus led to a massive environmental pollution. Persistence of chlordecone in soils and water for numerous decades even centuries causes global public health and socio-economic concerns. In order to investigate the biodegradability of chlordecone, microbial enrichment cultures from soils contaminated by chlordecone or other organochlorines and from sludge of a wastewater treatment plant have been conducted. Different experimental procedures including original microcosms were carried out anaerobically over long periods of time. GC-MS monitoring resulted in the detection of chlorinated derivatives in several cultures, consistent with chlordecone biotransformation. More interestingly, disappearance of chlordecone (50 µg/mL) in two bacterial consortia was concomitant with the accumulation of a major metabolite of formula C9Cl5H3 (named B1) as well as two minor metabolites C10Cl9HO (named A1) and C9Cl4H4 (named B3). Finally, we report the isolation and the complete genomic sequences of two new Citrobacter isolates, closely related to Citrobacter amalonaticus, and that were capable of reproducing chlordecone transformation. Further characterization of these Citrobacter strains should yield deeper insights into the mechanisms involved in this transformation process.

Extensive transcription analysis of the Hyposoter didymator Ichnovirus genome in permissive and non-permissive lepidopteran host species.

Ichnoviruses are large dsDNA viruses that belong to the Polydnaviridae family. They are specifically associated with endoparasitic wasps of the family Ichneumonidae and essential for host parasitization by these wasps. We sequenced the Hyposoter didymator Ichnovirus (HdIV) encapsidated genome for further analysis of the transcription pattern of the entire set of HdIV genes following the parasitization of four different lepidopteran host species. The HdIV genome was found to consist of at least 50 circular dsDNA molecules, carrying 135 genes, 98 of which formed 18 gene families. The HdIV genome had general features typical of Ichnovirus (IV) genomes and closely resembled that of the IV carried by Hyposoter fugitivus. Subsequent transcriptomic analysis with Illumina technology during the course of Spodoptera frugiperda parasitization led to the identification of a small subset of less than 30 genes with high RPKM values in permissive hosts, consisting with these genes encoding crucial virulence proteins. Comparisons of HdIV expression profiles between host species revealed differences in transcript levels for given HdIV genes between two permissive hosts, S. frugiperda and Pseudoplusia includens. However, we found no evident intrafamily gene-specific transcription pattern consistent with the presence of multigenic families within IV genomes reflecting an ability of the wasps concerned to exploit different host species. Interestingly, in two non-permissive hosts, Mamestra brassiccae and Anticarsia gemmatalis (most of the parasitoid eggs were eliminated by the host cellular immune response), HdIV genes were generally less strongly transcribed than in permissive hosts. This suggests that successful parasitism is dependent on the expression of given HdIV genes exceeding a particular threshold value. These results raise questions about the mecanisms involved in regulating IV gene expression according to the nature of the lepidopteran host species encountered.

Functional endogenous viral elements in the genome of the parasitoid wasp Cotesia congregata: insights into the evolutionary dynamics of bracoviruses.

Bracoviruses represent the most complex endogenous viral elements (EVEs) described to date. Nudiviral genes have been hosted within parasitoid wasp genomes since approximately 100 Ma. They play a crucial role in the wasp life cycle as they produce bracovirus particles, which are injected into parasitized lepidopteran hosts during wasp oviposition. Bracovirus particles encapsidate multiple dsDNA circles encoding virulence genes. Their expression in parasitized caterpillars is essential for wasp parasitism success. Here, we report on the genomic organization of the proviral segments (i.e. master sequences used to produce the encapsidated dsDNA circles) present in the Cotesia congregata parasitoid wasp genome. The provirus is composed of a macrolocus, comprising two-thirds of the proviral segments and of seven dispersed loci, each containing one to three segments. Comparative genomic analyses with closely related species gave insights into the evolutionary dynamics of bracovirus genomes. Conserved synteny in the different wasp genomes showed the orthology of the proviral macrolocus across different species. The nudiviral gene odv-e66-like1 is conserved within the macrolocus, suggesting an ancient co-localization of the nudiviral genome and bracovirus proviral segments. By contrast, the evolution of proviral segments within the macrolocus has involved a series of lineage-specific duplications.

Germline BAP1 mutations predispose to renal cell carcinomas.

The genetic cause of some familial nonsyndromic renal cell carcinomas (RCC) defined by at least two affected first-degree relatives is unknown. By combining whole-exome sequencing and tumor profiling in a family prone to cases of RCC, we identified a germline BAP1 mutation c.277A>G (p.Thr93Ala) as the probable genetic basis of RCC predisposition. This mutation segregated with all four RCC-affected relatives. Furthermore, BAP1 was found to be inactivated in RCC-affected individuals from this family. No BAP1 mutations were identified in 32 familial cases presenting with only RCC. We then screened for germline BAP1 deleterious mutations in familial aggregations of cancers within the spectrum of the recently described BAP1-associated tumor predisposition syndrome, including uveal melanoma, malignant pleural mesothelioma, and cutaneous melanoma. Among the 11 families that included individuals identified as carrying germline deleterious BAP1 mutations, 6 families presented with 9 RCC-affected individuals, demonstrating a significantly increased risk for RCC. This strongly argues that RCC belongs to the BAP1 syndrome and that BAP1 is a RCC-predisposition gene.

Alteration of fatty-acid-metabolizing enzymes affects mitochondrial form and function in hereditary spastic paraplegia.

Hereditary spastic paraplegia (HSP) is considered one of the most heterogeneous groups of neurological disorders, both clinically and genetically. The disease comprises pure and complex forms that clinically include slowly progressive lower-limb spasticity resulting from degeneration of the corticospinal tract. At least 48 loci accounting for these diseases have been mapped to date, and mutations have been identified in 22 genes, most of which play a role in intracellular trafficking. Here, we identified mutations in two functionally related genes (DDHD1 and CYP2U1) in individuals with autosomal-recessive forms of HSP by using either the classical positional cloning or a combination of whole-genome linkage mapping and next-generation sequencing. Interestingly, three subjects with CYP2U1 mutations presented with a thin corpus callosum, white-matter abnormalities, and/or calcification of the basal ganglia. These genes code for two enzymes involved in fatty-acid metabolism, and we have demonstrated in human cells that the HSP pathophysiology includes alteration of mitochondrial architecture and bioenergetics with increased oxidative stress. Our combined results focus attention on lipid metabolism as a critical HSP pathway with a deleterious impact on mitochondrial bioenergetic function.

KIF1A missense mutations in SPG30, an autosomal recessive spastic paraplegia: distinct phenotypes according to the nature of the mutations.

The hereditary spastic paraplegias (HSPs) are a clinically and genetically heterogeneous group of neurodegenerative diseases characterised by progressive spasticity in the lower limbs. The nosology of autosomal recessive forms is complex as most mapped loci have been identified in only one or a few families and account for only a small percentage of patients. We used next-generation sequencing focused on the SPG30 chromosomal region on chromosome 2q37.3 in two patients from the original linked family. In addition, wide genome scan and candidate gene analysis were performed in a second family of Palestinian origin. We identified a single homozygous mutation, p.R350G, that was found to cosegregate with the disease in the SPG30 kindred and was absent in 970 control chromosomes while affecting a strongly conserved amino acid at the end of the motor domain of KIF1A. Homozygosity and linkage mapping followed by mutation screening of KIF1A allowed us to identify a second mutation, p.A255V, in the second family. Comparison of the clinical features with the nature of the mutations of all reported KIF1A families, including those reported recently with hereditary sensory and autonomic neuropathy, suggests phenotype-genotype correlations that may help to understand the mechanisms involved in motor neuron degeneration. We have shown that mutations in the KIF1A gene are responsible for SPG30 in two autosomal recessive HSP families. In published families, the nature of the KIF1A mutations seems to be of good predictor of the underlying phenotype and vice versa.

TTC21B contributes both causal and modifying alleles across the ciliopathy spectrum.

Ciliary dysfunction leads to a broad range of overlapping phenotypes, collectively termed ciliopathies. This grouping is underscored by genetic overlap, where causal genes can also contribute modifier alleles to clinically distinct disorders. Here we show that mutations in TTC21B, which encodes the retrograde intraflagellar transport protein IFT139, cause both isolated nephronophthisis and syndromic Jeune asphyxiating thoracic dystrophy. Moreover, although resequencing of TTC21B in a large, clinically diverse ciliopathy cohort and matched controls showed a similar frequency of rare changes, in vivo and in vitro evaluations showed a significant enrichment of pathogenic alleles in cases (P < 0.003), suggesting that TTC21B contributes pathogenic alleles to ∼5% of ciliopathy cases. Our data illustrate how genetic lesions can be both causally associated with diverse ciliopathies and interact in trans with other disease-causing genes and highlight how saturated resequencing followed by functional analysis of all variants informs the genetic architecture of inherited disorders.

High level of structural polymorphism driven by mobile elements in the Hox genomic region of the Chaetognath Spadella cephaloptera.

Little is known about the relationships between genome polymorphism, mobile element dynamics, and population size among animal populations. The chaetognath species Spadella cephaloptera offers a unique perspective to examine this issue because they display a high level of genetic polymorphism at the population level. Here, we have investigated in detail the extent of nucleotide and structural polymorphism in a region harboring Hox1 and several coding genes and presumptive functional elements. Sequencing of several bacterial artificial chromosome inserts representative of this nuclear region uncovered a high level of structural heterogeneity, which is mainly caused by the polymorphic insertion of a diversity of genetic mobile elements. By anchoring this variation through individual genotyping, we demonstrated that sequence diversity could be attributed to the allelic pool of a single population, which was confirmed by detection of extensive recombination within the genomic region studied. The high average level of nucleotide heterozygosity provides clues of selection in both coding and noncoding domains. This pattern stresses how selective processes remarkably cope with intense sequence turnover due to substitutions, mobile element insertions, and recombination to preserve the integrity of functional landscape. These findings suggest that genome polymorphism could provide pivotal information for future functional annotation of genomes.

High-throughput sequencing of a 4.1 Mb linkage interval reveals FLVCR2 deletions and mutations in lethal cerebral vasculopathy.

Rare lethal disease gene identification remains a challenging issue, but it is amenable to new techniques in high-throughput sequencing (HTS). Cerebral proliferative glomeruloid vasculopathy (PGV), or Fowler syndrome, is a severe autosomal recessive disorder of brain angiogenesis, resulting in abnormally thickened and aberrant perforating vessels leading to hydranencephaly. In three multiplex consanguineous families, genome-wide SNP analysis identified a locus of 14 Mb on chromosome 14. In addition, 280 consecutive SNPs were identical in two Turkish families unknown to be related, suggesting a founder mutation reducing the interval to 4.1 Mb. To identify the causative gene, we then specifically enriched for this region with sequence capture and performed HTS in a proband of seven families. Due to technical constraints related to the disease, the average coverage was only 7×. Nonetheless, iterative bioinformatic analyses of the sequence data identified mutations and a large deletion in the FLVCR2 gene, encoding a 12 transmembrane domain-containing putative transporter. A striking absence of alpha-smooth muscle actin immunostaining in abnormal vessels in fetal PGV brains, suggests a deficit in pericytes, cells essential for capillary stabilization and remodeling during brain angiogenesis. This is the first lethal disease-causing gene to be identified by comprehensive HTS of an entire linkage interval.

Analysis of virion structural components reveals vestiges of the ancestral ichnovirus genome.

Many thousands of endoparasitic wasp species are known to inject polydnavirus (PDV) particles into their caterpillar host during oviposition, causing immune and developmental dysfunctions that benefit the wasp larva. PDVs associated with braconid and ichneumonid wasps, bracoviruses and ichnoviruses respectively, both deliver multiple circular dsDNA molecules to the caterpillar. These molecules contain virulence genes but lack core genes typically involved in particle production. This is not completely unexpected given that no PDV replication takes place in the caterpillar. Particle production is confined to the wasp ovary where viral DNAs are generated from proviral copies maintained within the wasp genome. We recently showed that the genes involved in bracovirus particle production reside within the wasp genome and are related to nudiviruses. In the present work we characterized genes involved in ichnovirus particle production by analyzing the components of purified Hyposoter didymator Ichnovirus particles by LC-MS/MS and studying their organization in the wasp genome. Their products are conserved among ichnovirus-associated wasps and constitute a specific set of proteins in the virosphere. Strikingly, these genes are clustered in specialized regions of the wasp genome which are amplified along with proviral DNA during virus particle replication, but are not packaged in the particles. Clearly our results show that ichnoviruses and bracoviruses particles originated from different viral entities, thus providing an example of convergent evolution where two groups of wasps have independently domesticated viruses to deliver genes into their hosts.

Identification of bracovirus particle proteins and analysis of their transcript levels at the stage of virion formation.

Polydnaviruses (PDVs) are unique symbiotic viruses associated with parasitic wasps; they replicate only in the calyx cells of a wasp's ovaries and are transferred at oviposition along with the parasitoid egg into the lepidopteran host. The DNA packaged in the viral particles encodes factors that manipulate the host's immune defences and development to benefit the parasitoid. PDVs are found in two subfamilies of ichneumonids (ichnoviruses) and in braconids of the microgastroid complex (bracoviruses). We recently showed that the latter derive from an ancestral nudivirus, as 24 nudivirus-related genes were identified in ovaries of two distantly related braconids at the stage of virion formation. Here, we present a comprehensive analysis of the viral particle proteins of the Chelonus inanitus bracovirus (CiBV). Proteins of purified CiBV particles were analysed by mass spectrometry and amino acid sequences matched to the existing ovarian-cDNA database. In addition, transcript quantities of identified genes were measured by quantitative real-time PCR in female pupae at the onset and peak of virion formation and at corresponding stages in male pupae. This combined approach allowed the identification of 44 CiBV particle proteins: 16 were nudivirus-related, three had similarity to ovarian proteins of another braconid, 11 had similarity to cellular proteins and 14 had no similarity to known proteins. The transcripts of all of them increased in female, but not male, pupae. These data confirm the important contribution of nudivirus genes but also indicate the presence of many lineage- or species-specific proteins possibly involved in the parasitoid-host interaction.

Linkage analysis of the C4A/C4B copy number variation and polymorphisms of the adjacent steroid 21-hydroxylase gene in a healthy population.

Genes encoding the steroid 21-hydroxylase (CYP21A2) and the complement component C4 proteins (C4A and C4B) are located in the MHC region in a strongly linked structure named RCCX module. Previous studies found that carriers of C4B gene deficiency (C4B*Q0) have higher risk for cardiovascular diseases. A potential explanation is that lacking the C4B gene may result in altered function of the neighboring CYP21A2 gene. Therefore we sequenced the CYP21A2 gene in 96 healthy individuals to identify polymorphisms and to characterize their linkage pattern. Fifty-three variations were detected including a new one which alters the TATA-box of the gene. Only three known mutations (V281L, Q318X and R479L) associated with congenital adrenal hyperplasia, were found in 7, 2 and 1 subjects, respectively. Linkage analysis revealed that some variations exhibit strong correlation with the C4 copy number polymorphism and constituents of the MHC III region. Rare alleles of three polymorphisms were identified as components of the 8.1 ancestral haplotype. Haplotyping and family study confirmed that the variant alleles of two intronic SNPs were constituents of haplotype blocks lacking the C4B gene. These results suggest that variations of CYP21A2 gene can be involved in disease associations of the 8.1 haplotype and the C4B*Q0 genotype.

Polydnaviruses of braconid wasps derive from an ancestral nudivirus.

Many species of parasitoid wasps inject polydnavirus particles in order to manipulate host defenses and development. Because the DNA packaged in these particles encodes almost no viral structural proteins, their relation to viruses has been debated. Characterization of complementary DNAs derived from braconid wasp ovaries identified genes encoding subunits of a viral RNA polymerase and structural components of polydnavirus particles related most closely to those of nudiviruses--a sister group of baculoviruses. The conservation of this viral machinery in different braconid wasp lineages sharing polydnaviruses suggests that parasitoid wasps incorporated a nudivirus-related genome into their own genetic material. We found that the nudiviral genes themselves are no longer packaged but are actively transcribed and produce particles used to deliver genes essential for successful parasitism in lepidopteran hosts.

Chaetognath transcriptome reveals ancestral and unique features among bilaterians.

BACKGROUND: The chaetognaths (arrow worms) have puzzled zoologists for years because of their astonishing morphological and developmental characteristics. Despite their deuterostome-like development, phylogenomic studies recently positioned the chaetognath phylum in protostomes, most likely in an early branching. This key phylogenetic position and the peculiar characteristics of chaetognaths prompted further investigation of their genomic features. RESULTS: Transcriptomic and genomic data were collected from the chaetognath Spadella cephaloptera through the sequencing of expressed sequence tags and genomic bacterial artificial chromosome clones. Transcript comparisons at various taxonomic scales emphasized the conservation of a core gene set and phylogenomic analysis confirmed the basal position of chaetognaths among protostomes. A detailed survey of transcript diversity and individual genotyping revealed a past genome duplication event in the chaetognath lineage, which was, surprisingly, followed by a high retention rate of duplicated genes. Moreover, striking genetic heterogeneity was detected within the sampled population at the nuclear and mitochondrial levels but cannot be explained by cryptic speciation. Finally, we found evidence for trans-splicing maturation of transcripts through splice-leader addition in the chaetognath phylum and we further report that this processing is associated with operonic transcription. CONCLUSION: These findings reveal both shared ancestral and unique derived characteristics of the chaetognath genome, which suggests that this genome is likely the product of a very original evolutionary history. These features promote chaetognaths as a pivotal model for comparative genomics, which could provide new clues for the investigation of the evolution of animal genomes.

Discovery and characterization of a new bacterial candidate division by an anaerobic sludge digester metagenomic approach.

We have constructed a large fosmid library from a mesophilic anaerobic digester and explored its 16S rDNA diversity using a high-density filter DNA-DNA hybridization procedure. We identified a group of 16S rDNA sequences forming a new bacterial lineage named WWE3 (Waste Water of Evry 3). Only one sequence from the public databases shares a sequence identity above 80% with the WWE3 group which hence cannot be affiliated to any known or candidate prokaryotic division. Despite representing a non-negligible fraction (5% of the 16S rDNA sequences) of the bacterial population of this digester, the WWE3 bacteria could not have been retrieved using the conventional 16S rDNA amplification procedure due to their unusual 16S rDNA gene sequence. WWE3 bacteria were detected by polymerase chain reaction (PCR) in various environments (anaerobic digesters, swine lagoon slurries and freshwater biofilms) using newly designed specific PCR primer sets. Fluorescence in situ hybridization (FISH) analysis of sludge samples showed that WWE3 microorganisms are oval-shaped and located deep inside sludge flocs. Detailed phylogenetic analysis showed that WWE3 bacteria form a distinct monophyletic group deeply branching apart from all known bacterial divisions. A new bacterial candidate division status is proposed for this group.