Transcriptome analysis has important applications in many biological fields. However, assembling a transcriptome without a known reference remains a challenging task requiring algorithmic improvements. We present two methods for substantially improving transcriptome de novo assembly. The first method relies on the observation that the use of a single k-mer length by current de novo assemblers is suboptimal to assemble transcriptomes where the sequence coverage of transcripts is highly heterogeneous. We present the Multiple-k method in which various k-mer lengths are used for de novo transcriptome assembly. We demonstrate its good performance by assembling de novo a published next-generation transcriptome sequence data set of Aedes aegypti, using the existing genome to check the accuracy of our method. The second method relies on the use of a reference proteome to improve the de novo assembly. We developed the Scaffolding using Translation Mapping (STM) method that uses mapping against the closest available reference proteome for scaffolding contigs that map onto the same protein. In a controlled experiment using simulated data, we show that the STM method considerably improves the assembly, with few errors. We applied these two methods to assemble the transcriptome of the non-model catfish Loricaria gr. cataphracta. Using the Multiple-k and STM methods, the assembly increases in contiguity and in gene identification, showing that our methods clearly improve quality and can be widely used. The new methods were used to assemble successfully the transcripts of the core set of genes regulating tooth development in vertebrates, while classic de novo assembly failed.
Next-generation sequencing technologies are rapidly generating whole-genome datasets for an increasing number of organisms. However, phylogenetic reconstruction of genomic data remains difficult because de novo assembly for non-model genomes and multi-genome alignment are challenging. To greatly simplify the analysis, we present an Assembly and Alignment-Free (AAF) method ( https://sourceforge.net/projects/aaf-phylogeny ) that constructs phylogenies directly from unassembled genome sequence data, bypassing both genome assembly and alignment. Using mathematical calculations, models of sequence evolution, and simulated sequencing of published genomes, we address both evolutionary and sampling issues caused by direct reconstruction, including homoplasy, sequencing errors, and incomplete sequencing coverage. From these results, we calculate the statistical properties of the pairwise distances between genomes, allowing us to optimize parameter selection and perform bootstrapping. As a test case with real data, we successfully reconstructed the phylogeny of 12 mammals using raw sequencing reads. We also applied AAF to 21 tropical tree genome datasets with low coverage to demonstrate its effectiveness on non-model organisms. Our AAF method opens up phylogenomics for species without an appropriate reference genome or high sequence coverage, and rapidly creates a phylogenetic framework for further analysis of genome structure and diversity among non-model organisms.
Metabarcoding potentially offers a rapid and cheap method of monitoring biodiversity, but real-world applications are few. We investigated its utility in studying patterns of litter arthropod diversity and composition in the tropics. We collected litter arthropods from 35 matched forest-plantation sites across Xishuangbanna, southwestern China. A new primer combination and the MiSeq platform were used to amplify and sequence a wide variety of litter arthropods using simulated and real-world communities. Quality filtered reads were clustered into 3,624 MOTUs at ≥97% similarity and the taxonomy of each MOTU was predicted. We compared diversity and compositional differences between forests and plantations (rubber and tea) for all MOTUs and for eight arthropod groups. We obtained ~100% detection rate after in silico sequencing six mock communities with known arthropod composition. Ordination showed that rubber, tea and forest communities formed distinct clusters. α-diversity declined significantly between forests and adjacent plantations for more arthropod groups in rubber than tea, and diversity of order Orthoptera increased significantly in tea. Turnover was higher in forests than plantations, but patterns differed among groups. Metabarcoding is useful for quantifying diversity patterns of arthropods under different land-uses and the MiSeq platform is effective for arthropod metabarcoding in the tropics.
25 June 2008 in Philosophical Transactions of the Royal Society B: Biological Sciences
The biogeographic patterns in sexually reproducing animals in island archipelagos may be interpreted as reflecting the importance of allopatric speciation. However, as the forms are allopatric, their reproductive isolation is largely untestable. A historical perspective integrating geology and molecular phylogeny reveals specific cases where ancient precursor islands coalesce, which allows the application of population genetics to critically test genetic isolation. The Anolis populations on Martinique in the Lesser Antilles are one such case where species-level populations on ancient precursor islands ( ca 6–8 Myr BP) have met relatively recently. The distribution of the mtDNA lineages is tightly linked to the precursor island, but the population genetic analysis of microsatellite variation in large samples shows no evidence of restricted genetic exchange between these forms in secondary contact. This tests, and rejects, the hypothesis of simple allopatric speciation in these forms. By contrast, Martinique has pronounced environmental zonation, to which anoles are known to adapt. The population genetic analysis shows restricted genetic exchange across the ecotone between xeric coastal habitat and montane rainforest. This does not indicate full ecological speciation in these forms, but it does suggest the relative importance of the role of ecology in speciation in general.
1 March 2001 in Molecular Phylogenetics and Evolution
The evolution of viviparity in squamates has been the focus of much scientific attention in previous years. In particular, the possibility of the transition from viviparity back to oviparity has been the subject of a vigorous debate. Some studies have suggested this reversal is more frequent than previously thought. However, none of them provide conclusive evidence. We investigated this problem by studying the phylogenetic relationships between oviparous and viviparous lineages of the reproductively bimodal lizard species Zootoca vivipara. Our results show that viviparous populations are not monophyletic, and that several evolutionary transitions in parity mode have occurred. The most parsimonious scenario involves a single origin of viviparity followed by a reversal back to oviparity. This is the first study with a strongly supported phylogenetic framework supporting a transition from viviparity to oviparity.