Widespread infection, diversification and old host associations of Nosema Microsporidia in European freshwater gammarids (Amphipoda).

The microsporidian genus Nosema is primarily known to infect insects of economic importance stimulating high research interest, while other hosts remain understudied. Nosema granulosis is one of the formally described Nosema species infecting amphipod crustaceans, being known to infect only two host species. Our first aim was to characterize Nosema spp. infections in different amphipod species from various European localities using the small subunit ribosomal DNA (SSU) marker. Second, we aimed to assess the phylogenetic diversity, host specificity and to explore the evolutionary history that may explain the diversity of gammarid-infecting Nosema lineages by performing a phylogenetic reconstruction based on RNA polymerase II subunit B1 (RPB1) gene sequences. For the host species Gammarus balcanicus, we also analyzed whether parasites were in excess in females to test for sex ratio distortion in relation with Nosema infection. We identified Nosema spp. in 316 individuals from nine amphipod species being widespread in Europe. The RPB1-based phylogenetic reconstruction using newly reported sequences and available data from other invertebrates identified 39 haplogroups being associated with amphipods. These haplogroups clustered into five clades (A-E) that did not form a single amphipod-infecting monophyletic group. Closely related sister clades C and D correspond to Nosema granulosis. Clades A, B and E might represent unknown Nosema species infecting amphipods. Host specificity seemed to be variable with some clades being restricted to single hosts, and some that could be found in several host species. We show that Nosema parasite richness in gammarid hosts is much higher than expected, illustrating the advantage of the use of RPB1 marker over SSU. Finally, we found no hint of sex ratio distortion in Nosema clade A infecting G. balcanicus. This study shows that Nosema spp. are abundant, widespread and diverse in European gammarids. Thus, Nosema is as diverse in aquatic as in terrestrial hosts.

Investigating Wolbachia symbiont-mediated host protection against a bacterial pathogen using a natural Wolbachia nuclear insert.

Wolbachia bacterial endosymbionts provide protection against pathogens in various arthropod species but the underlying mechanisms remain misunderstood. By using a natural Wolbachia nuclear insert (f-element) in the isopod Armadillidium vulgare, we explored whether Wolbachia presence is mandatory to observe protection in this species or the presence of its genes is sufficient. We assessed survival of closely related females carrying or lacking the f-element (and lacking Wolbachia) challenged with the bacterial pathogen Salmonella enterica. Despite marginal significant effects, the f-element alone did not appear to confer survival benefits to its host, suggesting that Wolbachia presence in cells is crucial for protection.

Heterogeneous distribution of sex ratio distorters in natural populations of the isopod Armadillidium vulgare.

In the isopod Armadillidium vulgare, many females produce progenies with female-biased sex ratios, owing to two feminizing sex ratio distorters (SRD): Wolbachia endosymbionts and the f element. We investigated the distribution and population dynamics of these SRD and mitochondrial DNA variation in 16 populations from Europe and Japan. Confirming and extending results from the 1990s, we found that the SRD are present at variable frequencies in populations and that the f element is overall more frequent than Wolbachia. The two SRD never co-occur at high frequency in any population, suggesting an apparent mutual exclusion. We also detected Wolbachia or the f element in some males, which probably reflects insufficient titer to induce feminization or presence of masculinizing alleles. Our results are consistent with a single integration event of a Wolbachia genome in the A. vulgare genome at the origin of the f element, which contradicts an earlier hypothesis of frequent losses and gains. We identified strong linkage between Wolbachia strains and mitochondrial haplotypes, but no association between the f element and mitochondrial background. Our results open new perspectives on SRD evolutionary dynamics in A. vulgare, the evolution of genetic conflicts and their impact on the variability of sex determination systems.

Chromosome-level genome assembly reveals homologous chromosomes and recombination in asexual rotifer Adineta vaga.

Bdelloid rotifers are notorious as a speciose ancient clade comprising only asexual lineages. Thanks to their ability to repair highly fragmented DNA, most bdelloid species also withstand complete desiccation and ionizing radiation. Producing a well-assembled reference genome is a critical step to developing an understanding of the effects of long-term asexuality and DNA breakage on genome evolution. To this end, we present the first high-quality chromosome-level genome assemblies for the bdelloid Adineta vaga, composed of six pairs of homologous (diploid) chromosomes with a footprint of paleotetraploidy. The observed large-scale losses of heterozygosity are signatures of recombination between homologous chromosomes, either during mitotic DNA double-strand break repair or when resolving programmed DNA breaks during a modified meiosis. Dynamic subtelomeric regions harbor more structural diversity (e.g., chromosome rearrangements, transposable elements, and haplotypic divergence). Our results trigger the reappraisal of potential meiotic processes in bdelloid rotifers and help unravel the factors underlying their long-term asexual evolutionary success.

Assessing the Impact of a Viral Infection on the Expression of Transposable Elements in the Cabbage Looper Moth (Trichoplusia ni).

Most studies of stress-induced transposable element (TE) expression have so far focused on abiotic sources of stress. Here, we analyzed the impact of an infection by the AcMNPV baculovirus on TE expression in a cell line (Tnms42) and midgut tissues of the cabbage looper moth (Trichoplusia ni). We find that a large fraction of TE families (576/636 in Tnms42 cells and 503/612 in midgut) is lowly expressed or not expressed at all [≤ 4 transcripts per million (TPM)] in the uninfected condition (median TPM of 0.37 in Tnms42 and 0.46 in midgut cells). In the infected condition, a total of 62 and 187 TE families were differentially expressed (DE) in midgut and Tnms42 cells, respectively, with more up- (46) than downregulated (16) TE families in the former and as many up- (91) as downregulated (96) TE families in the latter. Expression log2 fold changes of DE TE families varied from -4.95 to 9.11 in Tnms42 cells and from -4.28 to 7.66 in midgut. Large variations in expression profiles of DE TEs were observed depending on the type of cells and on time after infection. Overall, the impact of AcMNPV on TE expression in T. ni is moderate but potentially sufficient to affect TE activity and genome architecture. Interestingly, one host-derived TE integrated into AcMNPV genomes is highly expressed in infected Tnms42 cells. This result shows that virus-borne TEs can be expressed, further suggesting that they may be able to transpose and that viruses may act as vectors of horizontal transfer of TEs in insects.

Early Response to Dehydration Six-Like Transporter Family: Early Origin in Streptophytes and Evolution in Land Plants.

Carbon management by plants involves the activity of many sugar transporters, which play roles in sugar subcellular partitioning and reallocation at the whole organism scale. Among these transporters, the early response to dehydration six-like (ESL) monosaccharide transporters (MSTs) are still poorly characterized although they represent one of the largest sugar transporter subfamilies. In this study, we used an evolutionary genomic approach to infer the evolutionary history of this multigenic family. No ESL could be identified in the genomes of rhodophytes, chlorophytes, and the brown algae Ectocarpus siliculosus, whereas one ESL was identified in the genome of Klebsormidium nitens providing evidence for the early emergence of these transporters in Streptophytes. A phylogenetic analysis using the 519 putative ESL proteins identified in the genomes of 47 Embryophyta species and being representative of the plant kingdom has revealed that ESL protein sequences can be divided into three major groups. The first and second groups originated in the common ancestor of all spermaphytes [ζ: 340 million years ago (MYA)] and of angiosperms (ε: 170-235 MYA), respectively, and the third group originated before the divergence of rosids and asterids (γ/1R: 117 MYA). In some eudicots (Vitales, Malpighiales, Myrtales, Sapindales, Brassicales, Malvales, and Solanales), the ESL family presents remarkable expansions of gene copies associated with tandem duplications. The analysis of non-synonymous and synonymous substitutions for the dN/dS ratio of the ESL copies of the genus Arabidopsis has revealed that ESL genes are evolved under a purifying selection even though the progressive increase of dN/dS ratios in the three groups suggests subdiversification phenomena. To further explore the possible acquisition of novel functions by ESL MSTs, we identified the gene structure and promoter cis-acting elements for Arabidopsis thaliana ESL genes. The expression profiling of Arabidopsis ESL unraveled some gene copies that are almost constitutively expressed, whereas other gene copies display organ-preferential expression patterns. This study provides an evolving framework to better understand the roles of ESL transporters in plant development and response to environmental constraints.

Evolutionary transition to XY sex chromosomes associated with Y-linked duplication of a male hormone gene in a terrestrial isopod.

Sex chromosomes are highly variable in some taxonomic groups, but the evolutionary mechanisms underlying this diversity are not well understood. In terrestrial isopod crustaceans, evolutionary turnovers in sex chromosomes are frequent, possibly caused by Wolbachia, a vertically-transmitted endosymbiont causing male-to-female sex reversal. Here, we use surgical manipulations and genetic crosses, plus genome sequencing, to examine sex chromosomes in the terrestrial isopod Trachelipus rathkei. Although an earlier cytogenetics study suggested a ZZ/ZW sex chromosome system in this species, we surprisingly find multiple lines of evidence that in our study population, sex is determined by an XX/XY system. Consistent with a recent evolutionary origin for this XX/XY system, the putative male-specific region of the genome is small. The genome shows evidence of Y-linked duplications of the gene encoding the androgenic gland hormone, a major component of male sexual differentiation in isopods. Our analyses also uncover sequences horizontally acquired from past Wolbachia infections, consistent with the hypothesis that Wolbachia may have interfered with the evolution of sex determination in T. rathkei. Overall, these results provide evidence for the co-occurrence of multiple sex chromosome systems within T. rathkei, further highlighting the relevance of terrestrial isopods as models for the study of sex chromosome evolution.

Monitoring Insect Transposable Elements in Large Double-Stranded DNA Viruses Reveals Host-to-Virus and Virus-to-Virus Transposition.

The mechanisms by which transposable elements (TEs) can be horizontally transferred between animals are unknown, but viruses are possible candidate vectors. Here, we surveyed the presence of host-derived TEs in viral genomes in 35 deep sequencing data sets produced from 11 host-virus systems, encompassing nine arthropod host species (five lepidopterans, two dipterans, and two crustaceans) and six different double-stranded (ds) DNA viruses (four baculoviruses and two iridoviruses). We found evidence of viral-borne TEs in 14 data sets, with frequencies of viral genomes carrying a TE ranging from 0.01% to 26.33% for baculoviruses and from 0.45% to 7.36% for iridoviruses. The analysis of viral populations separated by a single replication cycle revealed that viral-borne TEs originating from an initial host species can be retrieved after viral replication in another host species, sometimes at higher frequencies. Furthermore, we detected a strong increase in the number of integrations in a viral population for a TE absent from the hosts' genomes, indicating that this TE has undergone intense transposition within the viral population. Finally, we provide evidence that many TEs found integrated in viral genomes (15/41) have been horizontally transferred in insects. Altogether, our results indicate that multiple large dsDNA viruses have the capacity to shuttle TEs in insects and they underline the potential of viruses to act as vectors of horizontal transfer of TEs. Furthermore, the finding that TEs can transpose between viral genomes of a viral species sets viruses as possible new niches in which TEs can persist and evolve.

Characterization of a Sex-Determining Region and Its Genomic Context via Statistical Estimates of Haplotype Frequencies in Daughters and Sons Sequenced in Pools.

Sex chromosomes are generally derived from a pair of autosomes that have acquired a locus controlling sex. Sex chromosomes may evolve reduced recombination around this locus and undergo a long process of molecular divergence. At that point, the original loci controlling sex may be difficult to pinpoint. This difficulty has affected many model species from mammals to birds to flies, which present highly diverged sex chromosomes. Identifying sex-controlling loci is easier in species with molecularly similar sex chromosomes. Here we aimed at pinpointing the sex-determining region (SDR) of Armadillidium vulgare, a terrestrial isopod with female heterogamety (ZW females and ZZ males) and whose sex chromosomes appear to show low genetic divergence. To locate the SDR, we assessed single-nucleotide polymorphism (SNP) allele frequencies in F1 daughters and sons sequenced in pools (pool-seq) in several families. We developed a Bayesian method that uses the SNP genotypes of individually sequenced parents and pool-seq data from F1 siblings to estimate the genetic distance between a given genomic region (contig) and the SDR. This allowed us to assign more than 43 Mb of contigs to sex chromosomes, and to demonstrate extensive recombination and very low divergence between these chromosomes. By taking advantage of multiple F1 families, we delineated a very short genomic region (∼65 kb) that presented no evidence of recombination with the SDR. In this short genomic region, the comparison of sequencing depths between sexes highlighted female-specific genes that have undergone recent duplication, and which may be involved in sex determination in A. vulgare.

Draft Genome Sequences of Thelohania contejeani and Cucumispora dikerogammari, Pathogenic Microsporidia of Freshwater Crustaceans.

We announce the draft genome sequences of two pathogenic microsporidia of European freshwater crustaceans, Thelohania contejeani (the causative agent of porcelain disease) and Cucumispora dikerogammari Both species are implicated in mass mortalities in natural populations of their crayfish and amphipod hosts, respectively.

Comparative Genomics of Strictly Vertically Transmitted, Feminizing Microsporidia Endosymbionts of Amphipod Crustaceans.

Microsporidia are obligate intracellular eukaryotic parasites of vertebrates and invertebrates. Microsporidia are usually pathogenic and undergo horizontal transmission or a mix of horizontal and vertical transmission. However, cases of nonpathogenic microsporidia, strictly vertically transmitted from mother to offspring, have been reported in amphipod crustaceans. Some of them further evolved the ability to feminize their nontransmitting male hosts into transmitting females. However, our understanding of the evolution of feminization in microsporidia is hindered by a lack of genomic resources. We report the sequencing and analysis of three strictly vertically transmitted microsporidia species for which feminization induction has been demonstrated (Nosema granulosis) or is strongly suspected (Dictyocoela muelleri and Dictyocoela roeselum), along with a draft genome assembly of their host Gammarus roeselii. Contrary to horizontally transmitted microsporidia that form environmental spores that can be purified, feminizing microsporidia cannot be easily isolated from their host cells. Therefore, we cosequenced symbiont and host genomic DNA and devised a computational strategy to obtain genome assemblies for the different partners. Genomic comparison with feminizing Wolbachia bacterial endosymbionts of isopod crustaceans indicated independent evolution of feminization in microsporidia and Wolbachia at the molecular genetic level. Feminization thus represents a remarkable evolutionary convergence of eukaryotic and prokaryotic microorganisms. Furthermore, a comparative genomics analysis of microsporidia allowed us to identify several candidate genes for feminization, involving functions such as DNA binding and membrane fusion. The genomic resources we generated contribute to establish Gammarus roeselii and its microsporidia symbionts as a new model to study the evolution of symbiont-mediated feminization.

Transposable Elements and the Evolution of Insects.

Insects are major contributors to our understanding of the interaction between transposable elements (TEs) and their hosts, owing to seminal discoveries, as well as to the growing number of sequenced insect genomes and population genomics and functional studies. Insect TE landscapes are highly variable both within and across insect orders, although phylogenetic relatedness appears to correlate with similarity in insect TE content. This correlation is unlikely to be solely due to inheritance of TEs from shared ancestors and may partly reflect preferential horizontal transfer of TEs between closely related species. The influence of insect traits on TE landscapes, however, remains unclear. Recent findings indicate that, in addition to being involved in insect adaptations and aging, TEs are seemingly at the cornerstone of insect antiviral immunity. Thus, TEs are emerging as essential insect symbionts that may have deleterious or beneficial consequences on their hosts, depending on context.

Characterization of a new case of XMLV (Bxv1) contamination in the human cell line Hep2 (clone 2B).

The use of misidentified cell lines contaminated by other cell lines and/or microorganisms has generated much confusion in the scientific literature. Detailed characterization of such contaminations is therefore crucial to avoid misinterpretation and ensure robustness and reproducibility of research. Here we use DNA-seq data produced in our lab to first confirm that the Hep2 (clone 2B) cell line (Sigma-Aldrich catalog number: 85011412-1VL) is indistinguishable from the HeLa cell line by mapping integrations of the human papillomavirus 18 (HPV18) at their expected loci on chromosome 8. We then show that the cell line is also contaminated by a xenotropic murine leukemia virus (XMLV) that is nearly identical to the mouse Bxv1 provirus and we characterize one Bxv1 provirus, located in the second intron of the pseudouridylate synthase 1 (PUS1) gene. Using an RNA-seq dataset, we confirm the high expression of the E6 and E7 HPV18 oncogenes, show that the entire Bxv1 genome is moderately expressed, and retrieve a Bxv1 splicing event favouring expression of the env gene. Hep2 (clone 2B) is the fourth human cell line so far known to be contaminated by the Bxv1 XMLV. This contamination has to be taken into account when using the cell line in future experiments.

------Widespread conservation and lineage-specific diversification of genome-wide DNA methylation patterns across arthropods.

Cytosine methylation is an ancient epigenetic modification yet its function and extent within genomes is highly variable across eukaryotes. In mammals, methylation controls transposable elements and regulates the promoters of genes. In insects, DNA methylation is generally restricted to a small subset of transcribed genes, with both intergenic regions and transposable elements (TEs) depleted of methylation. The evolutionary origin and the function of these methylation patterns are poorly understood. Here we characterise the evolution of DNA methylation across the arthropod phylum. While the common ancestor of the arthropods had low levels of TE methylation and did not methylate promoters, both of these functions have evolved independently in centipedes and mealybugs. In contrast, methylation of the exons of a subset of transcribed genes is ancestral and widely conserved across the phylum, but has been lost in specific lineages. A similar set of genes is methylated in all species that retained exon-enriched methylation. We show that these genes have characteristic patterns of expression correlating to broad transcription initiation sites and well-positioned nucleosomes, providing new insights into potential mechanisms driving methylation patterns over hundreds of millions of years.

Impact of transposable elements on genome size variation between two closely related crustacean species.

Identifying and quantifying genome size variation among species and understanding the underlying causes is a long-standing objective in evolutionary biology. Here, we investigated the basis of genome size variation between two closely related species of terrestrial isopods: Armadillidium vulgare and Armadillidium nasatum. The two species diverged 25 million years ago and the A. vulgare genome is ~500 megabases larger than the A. nasatum genome (1.7 vs. 1.2 gigabases, respectively). Our analyses indicated that genome size difference is essentially attributed to transposable elements (TEs). We found that the deletion rate may be slightly higher in A. nasatum than in A. vulgare, but it is unlikely to explain the observed genome size difference. As the two genomes largely share the same TE families, differential transpositional activity also contributes to the observed variation. Analyses of TE expression suggested that the cumulative expression level of all expressed TEs was higher in A. nasatum than in A. vulgare. Assuming TE expression level is a good proxy for TE transpositional activity, our results suggest that the two species may have recently been experiencing different TE transposition dynamics. Overall, our results illustrate the important impact TEs can have on genome structure and evolution between closely related species.

Dataset for sequencing and de novo assembly of the European endangered white-clawed crayfish (Austropotamobius pallipes) abdominal muscle transcriptome.

The white-clawed crayfish (Austropotamobius pallipes) is an endangered species in Europe with limited genomic information. Despite its conservation status there is no transcriptomic data available for A. pallipes in public databases. The data here represents the first transcriptome profile of the white-clawed crayfish generated using Illumina stranded RNA sequencing. Pair-end reads were assembled de novo with three separate transcriptome assemblers (Trinity, RNABloom, and RNASpades) followed by transcript assembly reduction and gene reconstruction using the EvidentialGene pipeline. The transcriptome was functionally annotated using InterProScan and genes coding for carbohydrate-active enzymes were identified through the dbCAN2 server. Raw fastq reads and the final version of the transcriptome assembly have been deposited in the NCBI-SRA (SRR10549898) and NCBI-TSA (GICG01) databases.

Wide spectrum and high frequency of genomic structural variation, including transposable elements, in large double-stranded DNA viruses.

Our knowledge of the diversity and frequency of genomic structural variation segregating in populations of large double-stranded (ds) DNA viruses is limited. Here, we sequenced the genome of a baculovirus (Autographa californica multiple nucleopolyhedrovirus [AcMNPV]) purified from beet armyworm (Spodoptera exigua) larvae at depths >195,000× using both short- (Illumina) and long-read (PacBio) technologies. Using a pipeline relying on hierarchical clustering of structural variants (SVs) detected in individual short- and long-reads by six variant callers, we identified a total of 1,141 SVs in AcMNPV, including 464 deletions, 443 inversions, 160 duplications, and 74 insertions. These variants are considered robust and unlikely to result from technical artifacts because they were independently detected in at least three long reads as well as at least three short reads. SVs are distributed along the entire AcMNPV genome and may involve large genomic regions (30,496 bp on average). We show that no less than 39.9 per cent of genomes carry at least one SV in AcMNPV populations, that the vast majority of SVs (75%) segregate at very low frequency (<0.01%) and that very few SVs persist after ten replication cycles, consistent with a negative impact of most SVs on AcMNPV fitness. Using short-read sequencing datasets, we then show that populations of two iridoviruses and one herpesvirus are also full of SVs, as they contain between 426 and 1,102 SVs carried by 52.4-80.1 per cent of genomes. Finally, AcMNPV long reads allowed us to identify 1,757 transposable elements (TEs) insertions, 895 of which are truncated and occur at one extremity of the reads. This further supports the role of baculoviruses as possible vectors of horizontal transfer of TEs. Altogether, we found that SVs, which evolve mostly under rapid dynamics of gain and loss in viral populations, represent an important feature in the biology of large dsDNA viruses.

Shedding Light on the Antimicrobial Peptide Arsenal of Terrestrial Isopods: Focus on Armadillidins, a New Crustacean AMP Family.

In crustaceans, antimicrobial peptides (AMPs) are clustered into four major groups according to their amino acid composition and structure: (1) single-domain peptides containing cysteine residues such as anti-lipopolysaccharide-factor (ALF), (2) multi-domain or chimeric AMPs such as crustins, (3) non-conventional AMPs, and (4) linear single-domain AMPs. The majority of AMPs has been described in commercially exploited crustaceans, particularly decapods living in aquatic environments (crab, shrimp, lobster, and crayfish). Here, we aimed at establishing the AMPs repertoire of terrestrial isopods (Oniscidea), an original suborder of crustaceans adapted to life outside of the aquatic environment. Using transcriptomic data from 21 species, we identified 110 ALF and 73 crustin sequences. We also characterized the full-length sequence of armadillidins from 17 species, similar to the AMP previously described in the terrestrial isopod Armadillidium vulgare. Furthermore, we tested the antimicrobial activity of three armadillidin peptides characterized from three distantly related species. This analysis revealed similar activity spectra against pathogens, despite extensive structural variation among the tested peptides. In addition to conventional crustacean AMPs, our work highlights armadillidins as a new and independent family of AMPs specific to the Oniscidea, thus opening new perspectives concerning the study of the immune system of terrestrial isopods.

Sex chromosomes control vertical transmission of feminizing Wolbachia symbionts in an isopod.

Microbial endosymbiosis is widespread in animals, with major ecological and evolutionary implications. Successful symbiosis relies on efficient vertical transmission through host generations. However, when symbionts negatively affect host fitness, hosts are expected to evolve suppression of symbiont effects or transmission. Here, we show that sex chromosomes control vertical transmission of feminizing Wolbachia endosymbionts in the isopod Armadillidium nasatum. Theory predicts that the invasion of an XY/XX species by cytoplasmic sex ratio distorters is unlikely because it leads to fixation of the unusual (and often lethal or infertile) YY genotype. We demonstrate that A. nasatum X and Y sex chromosomes are genetically highly similar and that YY individuals are viable and fertile, thereby enabling Wolbachia spread in this XY-XX species. Nevertheless, we show that Wolbachia cannot drive fixation of YY individuals, because infected YY females do not transmit Wolbachia to their offspring, unlike XX and XY females. The genetic basis fits the model of a Y-linked recessive allele (associated with an X-linked dominant allele), in which the homozygous state suppresses Wolbachia transmission. Moreover, production of all-male progenies by infected YY females restores a balanced sex ratio at the host population level. This suggests that blocking of Wolbachia transmission by YY females may have evolved to suppress feminization, thereby offering a whole new perspective on the evolutionary interplay between microbial symbionts and host sex chromosomes.