Barbara Feldmeyer

Short read Illumina data for the de novo assembly of a non-model snail species transcriptome (Radix balthica, Basommatophora, Pulmonata), and a comparison of assembler performance

BackgroundUntil recently, read lengths on the Solexa/Illumina system were too short to reliably assemble transcriptomes without a reference sequence, especially for non-model organisms. However, with read lengths up to 100 nucleotides available in the current version, an assembly without reference genome should be possible. For this study we created an EST data set for the common pond snail Radix balthica by Illumina sequencing of a normalized transcriptome. Performance of three different short read assemblers was compared with respect to: the number of contigs, their length, depth of coverage, their quality in various BLAST searches and the alignment to mitochondrial genes.ResultsA single sequencing run of a normalized RNA pool resulted in 16,923,850 paired end reads with median read length of 61 bases. The assemblies generated by VELVET, OASES, and SeqMan NGEN differed in the total number of contigs, contig length, the number and quality of gene hits obtained by BLAST searches against various databases, and contig performance in the mt genome comparison. While VELVET produced the highest overall number of contigs, a large fraction of these were of small size (< 200bp), and gave redundant hits in BLAST searches and the mt genome alignment. The best overall contig performance resulted from the NGEN assembly. It produced the second largest number of contigs, which on average were comparable to the OASES contigs but gave the highest number of gene hits in two out of four BLAST searches against different reference databases. A subsequent meta-assembly of the four contig sets resulted in larger contigs, less redundancy and a higher number of BLAST hits.ConclusionOur results document the first de novo transcriptome assembly of a non-model species using Illumina sequencing data. We show that de novo transcriptome assembly using this approach yields results useful for downstream applications, in particular if a meta-assembly of contig sets is used to increase contig quality. These results highlight the ongoing need for improvements in assembly methodology.

Climate-driven population divergence in sex-determining systems

Sex determination is a fundamental biological process, yet its mechanisms are remarkably diverse. In vertebrates, sex can be determined by inherited genetic factors or by the temperature experienced during embryonic development. However, the evolutionary causes of this diversity remain unknown. Here we show that live-bearing lizards at different climatic extremes of the species’ distribution differ in their sex-determining mechanisms, with temperature-dependent sex determination in lowlands and genotypic sex determination in highlands. A theoretical model parameterized with field data accurately predicts this divergence in sex-determining systems and the consequence thereof for variation in cohort sex ratios among years. Furthermore, we show that divergent natural selection on sex determination across altitudes is caused by climatic effects on lizard life history and variation in the magnitude of between-year temperature fluctuations. Our results establish an adaptive explanation for intra-specific divergence in sex-determining systems driven by phenotypic plasticity and ecological selection, thereby providing a unifying framework for integrating the developmental, ecological and evolutionary basis for variation in vertebrate sex determination.

Gene expression patterns associated with caste and reproductive status in ants: worker-specific genes are more derived than queen-specific ones.

Variation in gene expression leads to phenotypic diversity and plays a central role in caste differentiation of eusocial insect species. In social Hymenoptera, females with the same genetic background can develop into queens or workers, which are characterized by divergent morphologies, behaviours and lifespan. Moreover, many social insects exhibit behaviourally distinct worker castes, such as brood‐tenders and foragers. Researchers have just started to explore which genes are differentially expressed to achieve this remarkable phenotypic plasticity. Although the queen is normally the only reproductive individual in the nest, following her removal, young brood‐tending workers often develop ovaries and start to reproduce. Here, we make use of this ability in the ant Temnothorax longispinosus and compare gene expression patterns in the queens and three worker castes along a reproductive gradient. We found the largest expression differences between the queen and the worker castes (~2500 genes) and the smallest differences between infertile brood‐tenders and foragers (~300 genes). The expression profile of fertile workers is more worker‐like, but to a certain extent intermediate between the queen and the infertile worker castes. In contrast to the queen, a high number of differentially expressed genes in the worker castes are of unknown function, pointing to the derived status of hymenopteran workers within insects.

Unique evolutionary trajectories in repeated adaptation to hydrogen sulphide-toxic habitats of a neotropical fish (Poecilia mexicana).

Replicated ecological gradients are prime systems to study processes of molecular evolution underlying ecological divergence. Here, we investigated the repeated adaptation of the neotropical fish Poecilia mexicana to habitats containing toxic hydrogen sulphide (H2S) and compared two population pairs of sulphide‐adapted and ancestral fish by sequencing population pools of >200 individuals (Pool‐Seq). We inferred the evolutionary processes shaping divergence and tested the hypothesis of increase of parallelism from SNPs to molecular pathways. Coalescence analyses showed that the divergence occurred in the face of substantial bidirectional gene flow. Population divergence involved many short, widely dispersed regions across the genome. Analyses of allele frequency spectra suggest that differentiation at most loci was driven by divergent selection, followed by a selection‐mediated reduction of gene flow. Reconstructing allelic state changes suggested that selection acted mainly upon de novo mutations in the sulphide‐adapted populations. Using a corrected Jaccard index to quantify parallel evolution, we found a negligible proportion of statistically significant parallel evolution of Jcorr = 0.0032 at the level of SNPs, divergent genome regions (Jcorr = 0.0061) and genes therein (Jcorr = 0.0091). At the level of metabolic pathways, the overlap was Jcorr = 0.2545, indicating increasing parallelism with increasing level of biological integration. The majority of pathways contained positively selected genes in both sulphide populations. Hence, adaptation to sulphidic habitats necessitated adjustments throughout the genome. The largely unique evolutionary trajectories may be explained by a high proportion of de novo mutations driving the divergence. Our findings favour Goulds view that evolution is often the unrepeatable result of stochastic events with highly contingent effects.

Characterizing a hybrid zone between a cryptic species pair of freshwater snails

Characterizing hybrid zones and their dynamics is a central goal in evolutionary biology, but this is particularly challenging for morphologically cryptic species. The lack of conspicuous divergence between parental types means intermediate hybrid forms often go undetected. We aimed to detect and characterize a suspected hybrid zone between a pair of morphologically cryptic lineages of the freshwater snail, Radix. We sampled Radix from across a contact zone between two mitochondrial lineages (Radix balthica and an undescribed lineage termed ‘MOTU3’) and detected admixture between two nuclear genotype clusters, which were significantly but not categorically associated with the mitochondrial lineages. Using a model selection approach, we show that the admixture cline is best explained by an interaction between precipitation and temperature gradients over the area, rather than geographic distance. We thus hypothesize that the correlation with climatic gradients suggests environmental selection has played a role in maintaining the hybrid zone. In a 2050 climate change scenario, we furthermore predict an expansion of one of the nuclear clusters and a widening of the hybrid zone as the climate warms and dries.

Sponge communities of the Antarctic Peninsula: influence of environmental variables on species composition and richness

Sponge communities on the Antarctic continental shelf currently represent one of the most extensive sponge grounds in the world, and all sponge classes are known to occur in the Southern Ocean. Main objectives of this study conducted at the tip of the Antarctic Peninsula were (1) to identify all sampled sponges and (2) to investigate whether the species composition and species richness of Southern Ocean sponge communities in the area of the Antarctic Peninsula are significantly influenced by environmental variables. The studied material originated from 25 AGT catches and was sampled during the expedition ANT-XXIX/3 of RV Polarstern. Samples were collected in three large-scale areas in the vicinity of the Antarctic Peninsula: Bransfield Strait, Drake Passage and Weddell Sea. The following six environmental variables were measured from bottom water samples (except for sea-ice cover): depth (m), light transmission (%), oxygen (µmol/kg), salinity, sea-ice cover (%) and temperature (°C). Two hundred and sixty-three sponge samples were analyzed, and 81 species of 33 genera from all Porifera classes (Calcarea, Demospongiae, Hexactinellida and Homoscleromorpha) were identified. Total numbers of sponge species per sample station ranged from 1 to 29. A detrended correspondence analysis and a backward-stepwise model selection were performed to check whether species composition and richness were significantly influenced by environmental variables. The analyses revealed that none of the measured environmental variables significantly influenced species composition but that species richness was significantly influenced by (1) temperature and (2) the combination of temperature and depth. Results of this study are of crucial importance for development, performance and assessment of future protection strategies in case of ongoing climatic changes at the Antarctic Peninsula.

An Annotated Draft Genome for Radix auricularia (Gastropoda, Mollusca)

Molluscs are the second most species-rich phylum in the animal kingdom, yet only 11 genomes of this group have been published so far. Here, we present the draft genome sequence of the pulmonate freshwater snail Radix auricularia. Six whole genome shotgun libraries with different layouts were sequenced. The resulting assembly comprises 4,823 scaffolds with a cumulative length of 910 Mb and an overall read coverage of 72×. The assembly contains 94.6% of a metazoan core gene collection, indicating an almost complete coverage of the coding fraction. The discrepancy of ∼690 Mb compared with the estimated genome size of R. auricularia (1.6 Gb) results from a high repeat content of 70% mainly comprising DNA transposons. The annotation of 17,338 protein coding genes was supported by the use of publicly available transcriptome data. This draft will serve as starting point for further genomic and population genetic research in this scientifically important phylum.

Gene expression patterns underlying parasite-induced alterations in host behaviour and life history.

Many parasites manipulate their hosts’ phenotype. In particular, parasites with complex life cycles take control of their intermediate hosts’ behaviour and life history to increase transmission to their definitive host. The proximate mechanisms underlying these parasite‐induced alterations are poorly understood. The cestode Anomotaenia brevis affects the behaviour, life history and morphology of parasitized Temnothorax nylanderi ants and indirectly of their unparasitized nestmates. To gain insights on how parasites alter host phenotypes, we contrast brain gene expression patterns of T. nylanderi workers parasitized with the cestode, their unparasitized nestmates and unparasitized workers from unparasitized colonies. Over 400 differentially expressed genes between the three groups were identified, with most uniquely expressed genes detected in parasitized workers. Among these are genes that can be linked to the increased lifespan of parasitized workers. Furthermore, many muscle (functionality) genes are downregulated in these workers, potentially causing the observed muscular deformations and their inactive behaviour. Alterations in lifespan and activity could be adaptive for the parasite by increasing the likelihood that infected workers residing in acorns are eaten by their definitive host, a woodpecker. Our transcriptome analysis reveals numerous gene expression changes in parasitized workers and their uninfected nestmates and indicates possible routes of parasite manipulation. Although causality still needs to be established, parasite‐induced alterations in lifespan and host behaviour appear to be partly explained by morphological muscle atrophy instead of central nervous system interference, which is often the core of behavioural regulation. Results of this study will shed light upon the molecular basis of antagonistic species interactions.

Life history evolution in social insects: A female perspective

Social insects are known for their unusual life histories with fecund, long-lived queens and sterile, short-lived workers. We review ultimate factors underlying variation in life history strategies in female social insects, whose social life reshapes common trade-offs, such as the one between fecundity and longevity. Interspecific life history variation is associated with colony size, mediated by changes in division of labour and extrinsic mortality. In addition to the ratio of juvenile to adult mortality, social factors such as queen number influence life history trajectories. We discuss two hypotheses explaining why queen fecundity and lifespan is higher in single-queen societies and suggest further research directions on the evolution of life history variation in social insects.

Positive selection in development and growth rate regulation genes involved in species divergence of the genus Radix

BackgroundLife history traits like developmental time, age and size at maturity are directly related to fitness in all organisms and play a major role in adaptive evolution and speciation processes. Comparative genomic or transcriptomic approaches to identify positively selected genes involved in species divergence can help to generate hypotheses on the driving forces behind speciation. Here we use a bottom-up approach to investigate this hypothesis by comparative analysis of orthologous transcripts of four closely related European Radix species.ResultsSnails of the genus Radix occupy species specific distribution ranges with distinct climatic niches, indicating a potential for natural selection driven speciation based on ecological niche differentiation. We then inferred phylogenetic relationships among the four Radix species based on whole mt-genomes plus 23 nuclear loci. Three different tests to infer selection and changes in amino acid properties yielded a total of 134 genes with signatures of positive selection. The majority of these genes belonged to the functional gene ontology categories “reproduction” and “genitalia” with an overrepresentation of the functions “development” and “growth rate”.ConclusionsWe show here that Radix species divergence may be primarily enforced by selection on life history traits such as (larval-) development and growth rate. We thus hypothesise that life history differences may confer advantages under the according climate regimes, e.g., species occupying warmer and dryer habitats might have a fitness advantage with fast developing susceptible life stages, which are more tolerant to habitat desiccation.