Exploring Pandora's box: potential and pitfalls of low coverage genome surveys for evolutionary biology

Characterization of The Complete Mitochondrial Genome of the Bromeliad Crab Metopaulias depressus (Rathbun, 1896) (Crustacea: Decapoda: Brachyura: Sesarmidae)

Metopaulias depressus is a non-marine crab endemic to Jamaica that dwells in rainforest bromeliads and exhibits elaborate active parental care behavior. Current genomic resources on M. depressus are rare, limiting the understanding of its adaptation to terrestrial life in species that evolved from marine ancestors. This study reports the complete mitochondrial genome of M. depressus assembled using Sanger sequencing. The AT-rich mitochondrial genome of M. depressus is 15,765 bp in length and comprises 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22 transfer RNA genes. A single 691 bp-long intergenic space is assumed to be the control region (CR) or D-loop. A set of selective pressure analyses indicate that the entirety of the PCGs experience purifying selection. Cox1, cox2, nad5, cox3, and atp6 experience strong purifying selection, and atp8 experiences weak purifying selection compared to the rest of the PCGs. The secondary structures of most tRNA genes exhibit a standard ‘cloverleaf’ structure, with the exception of trnS1, which lacks the dihydroxyuridine (DHU) arm but not the loop, the trnH gene, which lacks the thymine pseudouracil cytosine (T) loop but not the arm, and trnM, which exhibits an overly developed T loop. A maximum likelihood phylogenetic analysis based on all PCGs indicated that M. depressus is more closely related to the genera Clistocoeloma, Nanosesarma, and Parasesarma than to Chiromantes, Geosesarma, and Orisarma. This study contributes to deciphering the phylogenetic relationships within the family Sesarmidae and represents a new genomic resource for this iconic crab species.

Phylogenomics of the longitarsal Colossendeidae: The evolutionary history of an Antarctic sea spider radiation

Sea spiders (Pycnogonida) constitute a group of marine benthic arthropods that has a particularly high species diversity in the Southern Ocean. The "longitarsal" group of the sea spider family Colossendeidae is especially abundant in this region. However, this group also includes some representatives from other oceans, which raises the question where the group originates from. Therefore, we here investigated the phylogeny of the group with a hybrid enrichment approach that yielded a dataset of 1607 genes and over one million base pairs. We obtained a well-resolved phylogeny of the group, which is mostly consistent with morphological data. The data support an Antarctic origin of the longitarsal Colossendeidae and multiple dispersal events to other regions, which occurred at different timescales. This scenario is consistent with evidence found in other groups of marine invertebrates and highlights the role of the Southern Ocean as a source for non-Antarctic biota, especially of the deep sea. Our results suggest an initially slow rate of diversification followed by a more rapid radiation possibly correlated with the mid-Miocene cooling of Antarctica, similar to what is found in other taxa.

The Draft Genome of the Invasive Walking Stick, Medauroidea extradendata, Reveals Extensive Lineage-Specific Gene Family Expansions of Cell Wall Degrading Enzymes in Phasmatodea

Plant cell wall components are the most abundant macromolecules on Earth. The study of the breakdown of these molecules is thus a central question in biology. Surprisingly, plant cell wall breakdown by herbivores is relatively poorly understood, as nearly all early work focused on the mechanisms used by symbiotic microbes to breakdown plant cell walls in insects such as termites. Recently, however, it has been shown that many organisms make endogenous cellulases. Insects, and other arthropods, in particular have been shown to express a variety of plant cell wall degrading enzymes in many gene families with the ability to break down all the major components of the plant cell wall. Here we report the genome of a walking stick, Medauroidea extradentata, an obligate herbivore that makes uses of endogenously produced plant cell wall degrading enzymes. We present a draft of the 3.3Gbp genome along with an official gene set that contains a diversity of plant cell wall degrading enzymes. We show that at least one of the major families of plant cell wall degrading enzymes, the pectinases, have undergone a striking lineage-specific gene family expansion in the Phasmatodea. This genome will be a useful resource for comparative evolutionary studies with herbivores in many other clades and will help elucidate the mechanisms by which metazoans breakdown plant cell wall components.

The Nuclear and Mitochondrial Genomes of the Facultatively Eusocial Orchid Bee Euglossa dilemma

Bees provide indispensable pollination services to both agricultural crops and wild plant populations, and several species of bees have become important models for the study of learning and memory, plant–insect interactions, and social behavior. Orchid bees (Apidae: Euglossini) are especially important to the fields of pollination ecology, evolution, and species conservation. Here we report the nuclear and mitochondrial genome sequences of the orchid bee Euglossa dilemma Bembé & Eltz. E. dilemma was selected because it is widely distributed, highly abundant, and it was recently naturalized in the southeastern United States. We provide a high-quality assembly of the 3.3 Gb genome, and an official gene set of 15,904 gene annotations. We find high conservation of gene synteny with the honey bee throughout 80 MY of divergence time. This genomic resource represents the first draft genome of the orchid bee genus Euglossa, and the first draft orchid bee mitochondrial genome, thus representing a valuable resource to the research community.

Development of genomic tools in a widespread tropical tree, Symphonia globulifera L.f.: a new low-coverage draft genome, SNP and SSR markers

Population genetic studies in tropical plants are often challenging because of limited information on taxonomy, phylogenetic relationships and distribution ranges, scarce genomic information and logistic challenges in sampling. We describe a strategy to develop robust and widely applicable genetic markers based on a modest development of genomic resources in the ancient tropical tree species Symphonia globulifera L.f. (Clusiaceae), a keystone species in African and Neotropical rainforests. We provide the first low-coverage (11X) fragmented draft genome sequenced on an individual from Cameroon, covering 1.027 Gbp or 67.5% of the estimated genome size. Annotation of 565 scaffolds (7.57 Mbp) resulted in the prediction of 1046 putative genes (231 of them containing a complete open reading frame) and 1523 exact simple sequence repeats (SSRs, microsatellites). Aligning a published transcriptome of a French Guiana population against this draft genome produced 923 high-quality single nucleotide polymorphisms. We also preselected genic SSRs in silico that were conserved and polymorphic across a wide geographical range, thus reducing marker development tests on rare DNA samples. Of 23 SSRs tested, 19 amplified and 18 were successfully genotyped in four S. globulifera populations from South America (Brazil and French Guiana) and Africa (Cameroon and São Tomé island, F_ST  = 0.34). Most loci showed only population-specific deviations from Hardy-Weinberg proportions, pointing to local population effects (e.g. null alleles). The described genomic resources are valuable for evolutionary studies in Symphonia and for comparative studies in plants. The methods are especially interesting for widespread tropical or endangered taxa with limited DNA availability.

Connectivity in the cold: the comparative population genetics of vent-endemic fauna in the Scotia Sea, Southern Ocean

We report the first comparative population genetics study for vent fauna in the Southern Ocean using cytochrome C oxidase I and microsatellite markers. Three species are examined: the kiwaid squat lobster, Kiwa tyleri, the peltospirid gastropod, Gigantopelta chessoia, and a lepetodrilid limpet, Lepetodrilus sp., collected from vent fields 440 km apart on the East Scotia Ridge (ESR) and from the Kemp Caldera on the South Sandwich Island Arc, ~95 km eastwards. We report no differentiation for all species across the ESR, consistent with panmixia or recent range expansions. A lack of differentiation is notable for Kiwa tyleri, which exhibits extremely abbreviated lecithotrophic larval development, suggestive of a very limited dispersal range. Larval lifespans may, however, be extended by low temperature-induced metabolic rate reduction in the Southern Ocean, muting the impact of dispersal strategy on patterns of population structure. COI diversity patterns suggest all species experienced demographic bottlenecks or selective sweeps in the past million years and possibly at different times. ESR and Kemp limpets are divergent, although with evidence of very recent ESR-Kemp immigration. Their divergence, possibility indicative of incipient speciation, along with the absence of the other two species at Kemp, may be the consequence of differing dispersal capabilities across a ~1000 m depth range and/or different selective regimes between the two areas. Estimates of historic and recent limpet gene flow between the ESR and Kemp are consistent with predominantly easterly currents and potentially therefore, cross-axis currents on the ESR, with biogeographic implications for the region.

Development and characterization of novel microsatellite markers by Next Generation Sequencing for the blue and red shrimp Aristeus antennatus

The blue and red shrimp, Aristeus antennatus, is a commercially important crustacean, in the Mediterranean Sea, which has been listed as a priority species for fishery management. Hypervariable microsatellite markers could be a useful tool to identify genetic stocks among geographically close fishing grounds. Potential microsatellite markers (97) identified from next-generation sequencing of an individual shrimp using a 454 GS Junior Pyrosequencer were tested on a preliminary panel of 15 individuals representing the four worldwide genetic stocks of the species from which 35 polymorphic loci were identified and used to characterize an additional 20 individuals from the Western Mediterranean Sea. In the Western Mediterranean sample, 32 out of 35 were polymorphic loci and the number of alleles per locus ranged from 2 to 14 and expected heterozygosity ranged from 0.050 to 0.968. No linkage disequilibrium was detected, indicating the independence of the loci. These novel microsatellites provide additional tools to address questions relating to genetic diversity, parentage studies and connectivity patterns of A. antennatus populations and help develop effective strategies to ensure long-term sustainability of this resource.

A First Insight into the Genome of the Filter-Feeder Mussel Mytilus galloprovincialis

Mussels belong to the phylum Mollusca, one of the largest and most diverse taxa in the animal kingdom. Despite their importance in aquaculture and in biology in general, genomic resources from mussels are still scarce. To broaden and increase the genomic knowledge in this family, we carried out a whole-genome sequencing study of the cosmopolitan Mediterranean mussel (Mytilus galloprovincialis). We sequenced its genome (32X depth of coverage) on the Illumina platform using three pair-end libraries with different insert sizes. The large number of contigs obtained pointed out a highly complex genome of 1.6 Gb where repeated elements seem to be widespread (~30% of the genome), a feature that is also shared with other marine molluscs. Notwithstanding the limitations of our genome sequencing, we were able to reconstruct two mitochondrial genomes and predict 10,891 putative genes. A comparative analysis with other molluscs revealed a gene enrichment of gene ontology categories related to multixenobiotic resistance, glutamate biosynthetic process, and the maintenance of ciliary structures.

Antarctic krill population genomics: apparent panmixia, but genome complexity and large population size muddy the water

Antarctic krill (Euphausia superba; hereafter krill) are an incredibly abundant pelagic crustacean which has a wide, but patchy, distribution in the Southern Ocean. Several studies have examined the potential for population genetic structuring in krill, but DNA-based analyses have focused on a limited number of markers and have covered only part of their circum-Antarctic range. We used mitochondrial DNA and restriction site-associated DNA sequencing (RAD-seq) to investigate genetic differences between krill from five sites, including two from East Antarctica. Our mtDNA results show no discernible genetic structuring between sites separated by thousands of kilometres, which is consistent with previous studies. Using standard RAD-seq methodology, we obtained over a billion sequences from >140 krill, and thousands of variable nucleotides were identified at hundreds of loci. However, downstream analysis found that markers with sufficient coverage were primarily from multicopy genomic regions. Careful examination of these data highlights the complexity of the RAD-seq approach in organisms with very large genomes. To characterize the multicopy markers, we recorded sequence counts from variable nucleotide sites rather than the derived genotypes; we also examined a small number of manually curated genotypes. Although these analyses effectively fingerprinted individuals, and uncovered a minor laboratory batch effect, no population structuring was observed. Overall, our results are consistent with panmixia of krill throughout their distribution. This result may indicate ongoing gene flow. However, krill's enormous population size creates substantial panmictic inertia, so genetic differentiation may not occur on an ecologically relevant timescale even if demographically separate populations exist.

Characterization of fossilized relatives of the White Spot Syndrome Virus in genomes of decapod crustaceans

Background

The White Spot Syndrome Virus (WSSV) is an important pathogen that infects a variety of decapod species and causes a highly contagious disease in penaeid shrimps. Mass mortalities caused by WSSV have pronounced commercial impact on shrimp aquaculture. Until now WSSV is the only known member of the virus family Nimaviridae, a group with obscure phylogenetic affinities. Its isolated position makes WSSV studies challenging due to large number of genes without homology in other viruses or cellular organisms.

Results

Here we report the discovery of an unusually large amount of sequences with high similarity to WSSV in a genomic library from the Jamaican bromeliad crab Metopaulias depressus. De novo assembly of these sequences allowed for the partial reconstruction of the genome of this endogenized virus with total length of 200 kbp encompassed in three scaffolds. The genome includes at least 68 putative open reading frames with homology in WSSV, most of which are intact. Among these, twelve orthologs of WSSV genes coding for non-structural proteins and nine genes known to code for the major components of the WSSV virion were discovered. Together with reanalysis of two similar cases of WSSV-like sequences in penaeid shrimp genomic libraries, our data allowed comparison of gene composition and gene order between different lineages related to WSSV. Furthermore, screening of published sequence databases revealed sequences with highest similarity to WSSV and the newly described virus in genomic libraries of at least three further decapod species. Analysis of the viral sequences detected in decapods suggests that they are less a result of contemporary WSSV infection, but rather originate from ancestral infection events. Phylogenetic analyses suggest that genes were acquired repeatedly by divergent viruses or viral strains of the Nimaviridae.

Conclusions

Our results shed new light on the evolution of the Nimaviridae and point to a long association of this viral group with decapod crustaceans.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0380-7) contains supplementary material, which is available to authorized users.

Regional differentiation and extensive hybridization between mitochondrial clades of the Southern Ocean giant sea spider Colossendeis megalonyx

Assessing the enormous diversity of Southern Ocean benthic species and their evolutionary histories is a central task in the era of global climate change. Based on mitochondrial markers, it was recently suggested that the circumpolar giant sea spider Colossendeis megalonyx comprises a complex of at least six cryptic species with mostly small and non-overlapping distribution ranges. Here, we expand the sampling to include over 500 mitochondrial COI sequences of specimens from around the Antarctic. Using multiple species delimitation approaches, the number of distinct mitochondrial OTUs increased from six to 15-20 with our larger dataset. In contrast to earlier studies, many of these clades show almost circumpolar distributions. Additionally, analysis of the nuclear internal transcribed spacer region for a subset of these specimens showed incongruence between nuclear and mitochondrial results. These mito-nuclear discordances suggest that several of the divergent mitochondrial lineages can hybridize and should not be interpreted as cryptic species. Our results suggest survival of C. megalonyx during Pleistocene glaciations in multiple refugia, some of them probably located on the Antarctic shelf, and emphasize the importance of multi-gene datasets to detect the presence of cryptic species, rather than their inference based on mitochondrial data alone.

The mitochondrial genomes of the caddisflies Sericostoma personatum and Thremma gallicum (Insecta: Trichoptera)

The mitochondrial genomes of the caddisfly species Sericostoma personatum and Thremma gallicum were sequenced on a 454 FLX and Illumina MiSeq platform, respectively. Reads were assembled de novo and remaining gaps in the S. personatum mitogenome closed by Sanger sequencing. The lengths of the assembled mitogenomes were 15,260 bp and 15,343 bp for S. personatum and T. gallicum, respectively. Both mitogenomes contained all 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, and the control region. The mitochondrial gene order of both caddisflies is identical with the typical insect gene order. These are the third and fourth published mitogenomes of the order Trichoptera of two formerly unexplored families and thus will be useful in future phylogenetic analysis.

QDD version 3.1: a user-friendly computer program for microsatellite selection and primer design revisited: experimental validation of variables determining genotyping success rate

Microsatellite marker development has been greatly simplified by the use of high-throughput sequencing followed by in silico microsatellite detection and primer design. However, the selection of markers designed by the existing pipelines depends either on arbitrary criteria, or older studies on PCR success. Based on wet laboratory experiments, we have identified the following factors that are most likely to influence genotyping success rate: alignment score between the primers and the amplicon; the distance between primers and microsatellites; the length of the PCR product; target region complexity and the number of reads underlying the sequence. The QDD pipeline has been modified to include these most pertinent factors in the output to help the selection of markers. Furthermore, new features are also included in the present version: (i) not only raw sequencing reads are accepted as input, but also contigs, allowing the analysis of assembled high-coverage data; (ii) input data can be both in fasta and fastq format to facilitate the use of Illumina and IonTorrent reads; (iii) A comparison to known transposable elements allows their detection; (iv) A contamination check can be carried out by BLASTing potential markers against the nucleotide (nt) database of NCBI; (v) QDD3 is now also available imbedded into a virtual machine making installation easier and operating system independent. It can be used both on command-line version as well as integrated into a Galaxy server, providing a user-friendly interface, as well as the possibility to utilize a large variety of NGS tools.

Telomere-independent ageing in the longest-lived non-colonial animal, Arctica islandica

The shortening of telomeres as a causative factor in ageing is a widely discussed hypothesis in ageing research. The study of telomere length and its regenerating enzyme telomerase in the longest-lived non-colonial animal on earth, Arctica islandica, should inform whether the maintenance of telomere length plays a role in reaching the extreme maximum lifespan (MLSP) of >500years in this species. Since longitudinal measurements on living animals cannot be achieved, a cross-sectional analysis of a short-lived (MLSP 40years from the Baltic Sea) and a long-lived population (MLSP 226years Northeast of Iceland) and in different tissues of young and old animals from the Irish Sea was performed. A high heterogeneity of telomere length was observed in investigated A. islandica over a wide age range (10-36years for the Baltic Sea, 11-194years for Irish Sea, 6-226years for Iceland). Constant telomerase activity and telomere lengths were detected at any age and in different tissues; neither correlated with age or population habitat. Stable telomere maintenance might contribute to the long lifespan of A. islandica. Telomere dynamics are no explanation for the distinct MLSPs of the examined populations and thus the cause of it remains to be investigated.

Random DNA libraries from three species of the stick insect genus Bacillus (Insecta: Phasmida): repetitive DNA characterization and first observation of polyneopteran MITEs

The repetitive DNA content of the stick insect species Bacillus rossius (facultative parthenogenetic), Bacillus grandii (gonochoric), and Bacillus atticus (obligate parthenogenetic) was analyzed through the survey of random genomic libraries roughly corresponding to 0.006% of the genome. By repeat masking, 19 families of transposable elements were identified (two LTR and six non-LTR retrotransposons; 11 DNA transposons). Moreover, a de novo analysis revealed, among the three libraries, the first MITE family observed in polyneopteran genomes. On the whole, transposable element abundance represented 23.3% of the genome in B. rossius, 22.9% in B. atticus, and 18% in B. grandii. Tandem repeat content in the three libraries is much lower: 1.32%, 0.64%, and 1.86% in B. rossius, B. grandii, and B. atticus, respectively. Microsatellites are the most abundant in all species. Minisatellites were only found in B. rossius and B. atticus, and five monomers belonging to the Bag320 satellite family were detected in B. atticus. Assuming the survey provides adequate representation of the relative genome, the obligate parthenogenetic species (B. atticus), compared with the other two species analyzed, does not show a lower transposable element content, as expected from some theoretical and empirical studies.

Isolation and characterization of nine polymorphic microsatellite markers for the deep-sea shrimp Nematocarcinus lanceopes (Crustacea: Decapoda: Caridea)

Background

The shrimp Nematocarcinus lanceopes Bate, 1888 is found in the deep sea around Antarctica and sub-Antarctic islands. Previous studies on mitochondrial data and species distribution models provided evidence for a homogenous circum-Antarctic population of N. lanceopes. However, to analyze the fine-scale population genetic structure and to examine influences of abiotic environmental conditions on population composition and genetic diversity, a set of fast evolving nuclear microsatellite markers is required.

Findings

We report the isolation and characterization of nine polymorphic microsatellite markers from the Antarctic deep-sea shrimp species Nematocarcinus lanceopes (Crustacea: Decapoda: Caridea). Microsatellite markers were screened in 55 individuals from different locations around the Antarctic continent. All markers were polymorphic with 9 to 25 alleles per locus. The observed heterozygosity ranged from 0.545 to 0.927 and the expected heterozygosity from 0.549 to 0.934.

Conclusions

The reported markers provide a novel tool to study genetic structure and diversity in Nematocarcinus lanceopes populations in the Southern Ocean and monitor effects of ongoing climate change in the region on the populations inhabiting these.