The Chalcidoidea bush of life: evolutionary history of a massive radiation of minute wasps.

Chalcidoidea are mostly parasitoid wasps that include as many as 500 000 estimated species. Capturing phylogenetic signal from such a massive radiation can be daunting. Chalcidoidea is an excellent example of a hyperdiverse group that has remained recalcitrant to phylogenetic resolution. We combined 1007 exons obtained with Anchored Hybrid Enrichment with 1048 ultra-conserved elements (UCEs) for 433 taxa including all extant families, >95% of all subfamilies, and 356 genera chosen to represent the vast diversity of the superfamily. Going back and forth between the molecular results and our collective knowledge of morphology and biology, we detected bias in the analyses that was driven by the saturation of nucleotide data. Our final results are based on a concatenated analysis of the least saturated exons and UCE datasets (2054 loci, 284 106 sites). Our analyses support an expected sister relationship with Mymarommatoidea. Seven previously recognized families were not monophyletic, so support for a new classification is discussed. Natural history in some cases would appear to be more informative than morphology, as illustrated by the elucidation of a clade of plant gall associates and a clade of taxa with planidial first-instar larvae. The phylogeny suggests a transition from smaller soft-bodied wasps to larger and more heavily sclerotized wasps, with egg parasitism as potentially ancestral for the entire superfamily. Deep divergences in Chalcidoidea coincide with an increase in insect families in the fossil record, and an early shift to phytophagy corresponds with the beginning of the "Angiosperm Terrestrial Revolution". Our dating analyses suggest a middle Jurassic origin of 174 Ma (167.3-180.5 Ma) and a crown age of 162.2 Ma (153.9-169.8 Ma) for Chalcidoidea. During the Cretaceous, Chalcidoidea may have undergone a rapid radiation in southern Gondwana with subsequent dispersals to the Northern Hemisphere. This scenario is discussed with regard to knowledge about the host taxa of chalcid wasps, their fossil record and Earth's palaeogeographic history.

Sequence, assembly and count datasets of viruses associated to the pine processionary moth Thaumetopoea pityocampa (Denis & Schiffermüller) (Lepidoptera, Notodontidae) identified from transcriptomic high-throughput sequencing.

The pine processionary moth Thaumetopoea pityocampa is a Lepidopteran pest species occurring in the Western Mediterranean. It causes heavy pine defoliations and it is a public and animal health concern because of its urticating caterpillars. Very little is known about the viruses associated to this species, as only two viruses were described so far. We here present a dataset corresponding to 34 viral transcripts, among which 27 could be confidently assigned to 9 RNA and DNA viral families (Iflaviridae, Reoviridae, Partitiviridae, Permutotetraviridae, Flaviviridae, Rhabdoviridae, Parvoviridae, Baculoviridae and PolyDNAviridae). These transcripts were identified from an original transcriptome assembled for the insect host, using both blast search and phylogenetic approaches. The data were acquired from 2 populations in Portugal and 2 populations in Italy. The transcripts were de novo assembled and used to identify viral sequences by homology searches. We also provide information about the populations and life stages in which each virus was identified. The data produced will allow to enrich the virus taxonomy in Lepidopteran hosts, and to develop PCR-based diagnostic tools to screen colonies across the range and determine the distribution and prevalence of the identified viral species.

Ultra-Conserved Elements and morphology reciprocally illuminate conflicting phylogenetic hypotheses in Chalcididae (Hymenoptera, Chalcidoidea).

Recent technical advances combined with novel computational approaches have promised the acceleration of our understanding of the tree of life. However, when it comes to hyperdiverse and poorly known groups of invertebrates, studies are still scarce. As published phylogenies will be rarely challenged by future taxonomists, careful attention must be paid to potential analytical bias. We present the first molecular phylogenetic hypothesis for the family Chalcididae, a group of parasitoid wasps, with a representative sampling (144 ingroups and seven outgroups) that covers all described subfamilies and tribes, and 82% of the known genera. Analyses of 538 Ultra-Conserved Elements (UCEs) with supermatrix (RAxML and IQTREE) and gene tree reconciliation approaches (ASTRAL, ASTRID) resulted in highly supported topologies in overall agreement with morphology but reveal conflicting topologies for some of the deepest nodes. To resolve these conflicts, we explored the phylogenetic tree space with clustering and gene genealogy interrogation methods, analyzed marker and taxon properties that could bias inferences and performed a thorough morphological analysis (130 characters encoded for 40 taxa representative of the diversity). This joint analysis reveals that UCEs enable attainment of resolution between ancestry and convergent/divergent evolution when morphology is not informative enough, but also shows that a systematic exploration of bias with different analytical methods and a careful analysis of morphological features is required to prevent publication of artifactual results. We highlight a GC content bias for maximum-likelihood approaches, an artifactual mid-point rooting of the ASTRAL tree and a deleterious effect of high percentage of missing data (>85% missing UCEs) on gene tree reconciliation methods. Based on the results we propose a new classification of the family into eight subfamilies and ten tribes that lay the foundation for future studies on the evolutionary history of Chalcididae.

Exploring systematic biases, rooting methods and morphological evidence to unravel the evolutionary history of the genus Ficus (Moraceae).

Despite many attempts in the Sanger sequencing era, the phylogeny of fig trees remains unresolved, which limits our ability to analyze the evolution of key traits that may have contributed to their evolutionary and ecological success. We used restriction-site-associated DNA sequencing (c. 420 kb) and 102 morphological characters to elucidate the relationships between 70 species of Ficus. To increase phylogenetic information for higher-level relationships, we targeted conserved regions and assembled paired reads into long loci to enable the retrieval of homologous loci in outgroup genomes. We compared morphological and molecular results to highlight discrepancies and reveal possible inference bias. For the first time, we recovered a monophyletic subgenus Urostigma (stranglers) and a clade with all gynodioecious Ficus. However, we show, with a new approach based on iterative principal component analysis, that it is not (and will probably never be) possible to homogenize evolutionary rates and GC content for all taxa before phylogenetic inference. Four competing positions for the root of the molecular tree are possible. The placement of section Pharmacosycea as sister to other fig trees is not supported by morphological data and considered a result of a long-branch attraction artefact to the outgroups. Regarding morphological features and indirect evidence from the pollinator tree of life, the topology that divides Ficus into monoecious versus gynodioecious species appears most plausible. It seems most likely that the ancestor of fig trees was a freestanding tree and active pollination is inferred as the ancestral state, contrary to previous hypotheses. However, ambiguity remains on the ancestral breeding system. Despite morphological plasticity, we advocate restoring a central role to morphology in our understanding of the evolution of Ficus, as it can help detect systematic errors that appear more pronounced with larger molecular datasets.

From refugia to contact: Pine processionary moth hybrid zone in a complex biogeographic setting.

Contact zones occur at the crossroad between specific dispersal routes and are facilitated by biogeographic discontinuities. Here, we focused on two Lepidoptera sister species that come in contact near the Turkish Straits System (TSS). We aimed to infer their phylogeographic histories in the Eastern Mediterranean and finely analyze their co-occurrence and hybridization patterns in this biogeographic context. We used molecular mitochondrial and nuclear markers to study 224 individuals from 42 localities. We used discordances between markers and complementary assignment methods to identify and map hybrids and parental individuals. We confirmed the parapatric distribution of Thaumetopoea pityocampa (Lepidoptera: Notodontidae) in the west and Thaumetopoea wilkinsoni in the east and identified a narrow contact zone. We identified several glacial refugia of T. wilkinsoni in southern Turkey with a strong east-west differentiation in this species. Unexpectedly, T. pityocampa crossed the TSS and occur in northern Aegean Turkey and some eastern Greek islands. We found robust evidence of introgression between the two species in a restricted zone in northwestern Turkey, but we did not identify any F1 individuals. The identified hybrid zone was mostly bimodal. The distributions and genetic patterns of the studied species were strongly influenced both by the Quaternary climatic oscillations and the complex geological history of the Aegean region. T. pityocampa and T. wilkinsoni survived the last glacial maximum in disjoint refugia and met in western Turkey at the edge of the recolonization routes. Expanding population of T. wilkinsoni constrained T. pityocampa to the western Turkish shore. Additionally, we found evidence of recurrent introgression by T. wilkinsoni males in several T. pityocampa populations. Our results suggest that some prezygotic isolation mechanisms, such as differences in timing of the adult emergences, might be a driver of the isolation between the sister species.

Restriction-site associated DNA markers provide new insights into the evolutionary history of the bark beetle genus Dendroctonus.

The bark beetle genus Dendroctonus contains some of the most economically important pests of conifers worldwide. Despite many attempts, there is no agreement today on the phylogenetic relationships within the genus, which limits our understanding of its evolutionary history. Here, using restriction-site associated DNA (RAD) markers from 70 specimens representing 17 species (85% of the known diversity) we inferred the phylogeny of the genus, its time of origin and biogeographic history, as well as the evolution of key ecological traits (host plants, larval behavior and adults' attack strategies). For all combinations of tested parameters (from 6444 to 23,570 RAD tags analyzed), the same, fully resolved topology was inferred. Our analyses suggest that the most recent common ancestor (mrca) of all extant Dendroctonus species was widely distributed across eastern Palearctic and western Nearctic during the early Miocene, from where species dispersed to other Holarctic regions. A first main inter-continental vicariance event occurred during early Miocene isolating the ancestors of D. armandi in the Palearctic, which was followed by the radiation of the main Dendroctonus lineages in North America. During the Late Miocene, the ancestor of the 'rufipennis' species group colonized north-east Palearctic regions from western North America, which was followed by a second main inter-continental vicariance event isolating Pleistocene populations in Asia (D. micans) and western North America (D. murrayanae and D. punctatus). The present study supports previous hypotheses explaining intercontinental range disjunctions across the Northern Hemisphere by the fragmentation of a continuous distribution due to climatic cooling, host range fragmentation and geological changes during the late Cenozoic. The reconstruction of ancestral ecological traits indicates that the mrca bored individual galleries and mass attacked the boles of pines. The gregarious feeding behavior of the larvae as well as the individual attack of the base of trees have apparently independently evolved twice in North America (in the 'rufipennis' and the 'valens' species groups), which suggests a higher adaptive potential than previously thought and may be of interest for plant protection and biodiversity conservation in a rapidly changing world.

Pushing the limits of whole genome amplification: successful sequencing of RADseq library from a single microhymenopteran (Chalcidoidea, Trichogramma).

A major obstacle to high-throughput genotyping of microhymenoptera is their small size. As species are difficult to discriminate, and because complexes may exist, the sequencing of a pool of specimens is hazardous. Thus, one should be able to sequence pangenomic markers (e.g., RADtags) from a single specimen. To date, whole genome amplification (WGA) prior to library construction is still a necessity as at most 10 ng of DNA can be obtained from single specimens (sometimes less). However, this amount of DNA is not compatible with manufacturer's requirements for commercial kits. Here we test the accuracy of the GenomiPhi kit V2 on Trichogramma wasps by comparing RAD libraries obtained from the WGA of single specimens (F0 and F1 generation, about1 ng input DNA for the WGA (0.17-2.9 ng)) and a biological amplification of genomic material (the pool of the progeny of the F1 generation). Globally, we found that 99% of the examined loci (up to 48,189 for one of the crosses, 109 bp each) were compatible with the mode of reproduction of the studied model (haplodiploidy) and Mendelian inheritance of alleles. The remaining 1% (0.01% of the analysed nucleotides) could represent WGA bias or other experimental/analytical bias. This study shows that the multiple displacement amplification method on which the GenomiPhi kit relies, could also be of great help for the high-throughput genotyping of microhymenoptera used for biological control, or other organisms from which only a very small amount of DNA can be extracted, such as human disease vectors (e.g.,  sandflies, fleas, ticks etc.).

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/.

Isolation, characterization and PCR multiplexing of microsatellite loci for a mite crop pest, Tetranychus urticae (Acari: Tetranychidae).

BACKGROUND: Tetranychus urticae is a highly polyphagous species with a cosmopolitan distribution that has the status of pest in more than 100 economically significant crops all over the world. Despite a number of previous efforts to isolate genetic markers, only a reduced set of microsatellite loci has been published. Taking advantage of the whole genome sequence of T. urticae that recently became available; we isolated and characterized a new set of microsatellite loci and tested the level of polymorphism across populations originating from a wide geographical area. RESULTS: A total of 42 microsatellite sequences widespread in the T. urticae genome were identified, the exact position in the genome recorded, and PCR amplification of microsatellite loci tested with primers defined here. Fourteen loci showed unambiguous genotype patterns and were further characterized. Three multiplex polymerase chain reaction sets were optimized in order to genotype a total of 24 polymorphic loci, including 10 previously published Tetranychus-specific loci. The microsatellite kits successfully amplified 686 individuals from 60 field populations for which we assessed the level of genetic diversity. The number of alleles per locus ranged from 3 to 16 and the expected heterozygosity values ranged from 0.12 to 0.81. Most of the loci displayed a significant excess of homozygous and did not model the Hardy-Weinberg equilibrium. This can be explained by the arrhenotokous mode of reproduction of T. urticae. CONCLUSIONS: These primers represent a valuable resource for robust studies on the genetic structure, dispersal and population biology of T. urticae, that can be used in managing this destructive agricultural pest.

Empirical assessment of RAD sequencing for interspecific phylogeny.

Next-generation sequencing opened up new possibilities in phylogenetics; however, choosing an appropriate method of sample preparation remains challenging. Here, we demonstrate that restriction-site-associated DNA sequencing (RAD-seq) generates useful data for phylogenomics. Analysis of our RAD library using current bioinformatic and phylogenetic tools produced 400× more sites than our Sanger approach (2,262,825 nt/species), fully resolving relationships between 18 species of ground beetles (divergences up to 17 My). This suggests that RAD-seq is promising to infer phylogeny of eukaryotic species, though potential biases need to be evaluated and new methodologies developed to take full advantage of such data.

Developing nuclear DNA phylogenetic markers in the angiosperm genus Leucadendron (Proteaceae): a next-generation sequencing transcriptomic approach.

Despite the recent advances in generating molecular data, reconstructing species-level phylogenies for non-models groups remains a challenge. The use of a number of independent genes is required to resolve phylogenetic relationships, especially for groups displaying low polymorphism. In such cases, low-copy nuclear exons and non-coding regions, such as 3' untranslated regions (3'-UTRs) or introns, constitute a potentially interesting source of nuclear DNA variation. Here, we present a methodology meant to identify new nuclear orthologous markers using both public-nucleotide databases and transcriptomic data generated for the group of interest by using next generation sequencing technology. To identify PCR primers for a non-model group, the genus Leucadendron (Proteaceae), we adopted a framework aimed at minimizing the probability of paralogy and maximizing polymorphism. We anchored when possible the right-hand primer into the 3'-UTR and the left-hand primer into the coding region. Seven new nuclear markers emerged from this search strategy, three of those included 3'-UTRs. We further compared the phylogenetic potential between our new markers and the ribosomal internal transcribed spacer region (ITS). The sequenced 3'-UTRs yielded higher polymorphism rates than the ITS region did. We did not find strong incongruences with the phylogenetic signal contained in the ITS region and the seven new designed markers but they strongly improved the phylogeny of the genus Leucadendron. Overall, this methodology is efficient in isolating orthologous loci and is valid for any non-model group given the availability of transcriptomic data.

Functional divergence in shrimp anti-lipopolysaccharide factors (ALFs): from recognition of cell wall components to antimicrobial activity.

Antilipopolysaccharide factors (ALFs) have been described as highly cationic polypeptides with a broad spectrum of potent antimicrobial activities. In addition, ALFs have been shown to recognize LPS, a major component of the Gram-negative bacteria cell wall, through conserved amino acid residues exposed in the four-stranded ß-sheet of their three dimensional structure. In penaeid shrimp, ALFs form a diverse family of antimicrobial peptides composed by three main variants, classified as ALF Groups A to C. Here, we identified a novel group of ALFs in shrimp (Group D ALFs), which corresponds to anionic polypeptides in which many residues of the LPS binding site are lacking. Both Group B (cationic) and Group D (anionic) shrimp ALFs were produced in a heterologous expression system. Group D ALFs were found to have impaired LPS-binding activities and only limited antimicrobial activity compared to Group B ALFs. Interestingly, all four ALF groups were shown to be simultaneously expressed in an individual shrimp and to follow different patterns of gene expression in response to a microbial infection. Group B was by far the more expressed of the ALF genes. From our results, nucleotide sequence variations in shrimp ALFs result in functional divergence, with significant differences in LPS-binding and antimicrobial activities. To our knowledge, this is the first functional characterization of the sequence diversity found in the ALF family.

Induction of a peptide with activity against a broad spectrum of pathogens in the Aedes aegypti salivary gland, following Infection with Dengue Virus.

The ultimate stage of the transmission of Dengue Virus (DENV) to man is strongly dependent on crosstalk between the virus and the immune system of its vector Aedes aegypti (Ae. aegypti). Infection of the mosquito's salivary glands by DENV is the final step prior to viral transmission. Therefore, in the present study, we have determined the modulatory effects of DENV infection on the immune response in this organ by carrying out a functional genomic analysis of uninfected salivary glands and salivary glands of female Ae. aegypti mosquitoes infected with DENV. We have shown that DENV infection of salivary glands strongly up-regulates the expression of genes that encode proteins involved in the vector's innate immune response, including the immune deficiency (IMD) and Toll signalling pathways, and that it induces the expression of the gene encoding a putative anti-bacterial, cecropin-like, peptide (AAEL000598). Both the chemically synthesized non-cleaved, signal peptide-containing gene product of AAEL000598, and the cleaved, mature form, were found to exert, in addition to antibacterial activity, anti-DENV and anti-Chikungunya viral activity. However, in contrast to the mature form, the immature cecropin peptide was far more effective against Chikungunya virus (CHIKV) and, furthermore, had strong anti-parasite activity as shown by its ability to kill Leishmania spp. Results from circular dichroism analysis showed that the immature form more readily adopts a helical conformation which would help it to cause membrane permeabilization, thus permitting its transfer across hydrophobic cell surfaces, which may explain the difference in the anti-pathogenic activity between the two forms. The present study underscores not only the importance of DENV-induced cecropin in the innate immune response of Ae. aegypti, but also emphasizes the broad-spectrum anti-pathogenic activity of the immature, signal peptide-containing form of this peptide.

Evidence of a bactericidal permeability increasing protein in an invertebrate, the Crassostrea gigas Cg-BPI.

A cDNA sequence with homologies to members of the LPS-binding protein and bactericidal/permeability-increasing protein (BPI) family was identified in the oyster Crassostrea gigas. The recombinant protein was found to bind LPS, to display bactericidal activity against Escherichia coli, and to increase the permeability of the bacterial cytoplasmic membrane. This indicated that it is a BPI rather than an LPS-binding protein. By in situ hybridization, the expression of the C. gigas BPI (Cg-bpi) was found to be induced in hemocytes after oyster bacterial challenge and to be constitutive in various epithelia of unchallenged oysters. Thus, Cg-bpi transcripts were detected in the epithelial cells of tissues/organs in contact with the external environment (mantle, gills, digestive tract, digestive gland diverticula, and gonad follicles). Therefore, Cg-BPI, whose expression profile and biological properties are reminiscent of mammalian BPIs, may provide a first line of defense against potential bacterial invasion. To our knowledge, this is the first characterization of a BPI in an invertebrate.