The influence of trematode parasite burden on gene expression in a mammalian host.

BACKGROUND: Parasites can profoundly impact their hosts and are responsible for a plethora of debilitating diseases. To identify global changes in host gene expression related to parasite infection, we sequenced, assembled, and annotated the liver transcriptomes of Balb/cj mice infected with the trematode parasite Schistosoma mansoni and compared the results to uninfected mice. We used two different methodologies (i.e. de novo and reference guided) to evaluate the influence of parasite sequences on host transcriptome assembly. RESULTS: Our results demonstrate that the choice of assembly methodology significantly impacted the proportion of parasitic reads detected from the host library, yet the presence of non-target (xenobiotic) sequences did not create significant structural errors in the assembly. After removing parasite sequences from the mouse transcriptomes, we analyzed host gene expression under different parasite infection levels and observed significant differences in the associated immunologic and metabolic responses based on infection level. In particular, genes associated with T-helper type 1 (Th-1) and T-helper type 2 (Th-2) were up-regulated in infected mice whereas genes related to amino acid and carbohydrate metabolism were down-regulated in infected mice. These changes in gene expression scale with infection status and likely impact the evolutionary fitness of hosts. CONCLUSIONS: Overall, our data indicate that a) infected mice reduce the expression of key metabolic genes in direct proportion to their infection level; b) infected mice similarly increase the expression of key immune genes in response to infection; c) patterns of gene expression correspond to the pathological symptoms of schistosomiasis; and d) identifying and filtering out non-target sequences (xenobiotics) improves differential expression prediction. Our findings identify parasite targets for RNAi or other therapies and provide a better understanding of the pathology and host immune repertoire involved in response to S. mansoni infections.

Formalizing the definition of meta-analysis in Molecular Ecology.

Meta-analysis, the statistical synthesis of pertinent literature to develop evidence-based conclusions, is relatively new to the field of molecular ecology, with the first meta-analysis published in the journal Molecular Ecology in 2003 (Slate & Phua 2003). The goal of this article is to formalize the definition of meta-analysis for the authors, editors, reviewers and readers of Molecular Ecology by completing a review of the meta-analyses previously published in this journal. We also provide a brief overview of the many components required for meta-analysis with a more specific discussion of the issues related to the field of molecular ecology, including the use and statistical considerations of Wright's FST and its related analogues as effect sizes in meta-analysis. We performed a literature review to identify articles published as 'meta-analyses' in Molecular Ecology, which were then evaluated by at least two reviewers. We specifically targeted Molecular Ecology publications because as a flagship journal in this field, meta-analyses published in Molecular Ecology have the potential to set the standard for meta-analyses in other journals. We found that while many of these reviewed articles were strong meta-analyses, others failed to follow standard meta-analytical techniques. One of these unsatisfactory meta-analyses was in fact a secondary analysis. Other studies attempted meta-analyses but lacked the fundamental statistics that are considered necessary for an effective and powerful meta-analysis. By drawing attention to the inconsistency of studies labelled as meta-analyses, we emphasize the importance of understanding the components of traditional meta-analyses to fully embrace the strengths of quantitative data synthesis in the field of molecular ecology.

The genomic proliferation of transposable elements in colonizing populations: Schistosoma mansoni in the new world.

Transposable elements (TEs) are mobile genes with an inherent ability to move within and among genomes. Theory predicts that TEs proliferate extensively during physiological stress due to the breakdown of TE repression systems. We tested this hypothesis in Schistosoma mansoni, a widespread trematode parasite that causes the human disease schistosomiasis. According to phylogenetic analysis, S. mansoni invaded the new world during the last 500 years. We hypothesized that new world strains of S. mansoni would have more copies of TEs than old world strains due to the physiological stress associated with invasion of the new world. We quantified the copy number of six TEs (Saci-1, Saci-2 and Saci-3, Perere-1, Merlin-sm1, and SmTRC1) in the genome and the transcriptome of old world and new world strains of S. mansoni, using qPCR relative quantification. As predicted, the genomes of new world parasites contain significantly more copies of class I and class II TEs in both laboratory and field strains. However, such differences are not observed in the transcriptome suggesting that either TE silencing mechanisms have reactivated to control the expression of these elements or the presence of inactive truncated copies of TEs.

A Plasma Sample Preparation for Mass Spectrometry using an Automated Workstation.

Sample preparation for mass spectrometry analysis in proteomics requires enzymatic cleavage of proteins into a peptide mixture. This process involves numerous incubation and liquid transfer steps in order to achieve denaturation, reduction, alkylation, and cleavage. Adapting this workflow onto an automated workstation can increase efficiency and reduce coefficients of variance, thereby providing more reliable data for statistical comparisons between sample types. We previously described an automated proteomic sample preparation workflow1. Here, we report the development of a more efficient and better controlled workflow with the following advantages: 1) The number of liquid transfer steps is reduced from nine to six by combining reagents; 2) Pipetting time is reduced by selective tip pipetting using a 96-position pipetting head with multiple channels; 3) Potential throughput is increased by the availability of up to 45 deck positions; 4) Complete enclosure of the system provides improved temperature and environmental control and reduces the potential for contamination of samples or reagents; and 5) The addition of stable isotope labeled peptides, as well as ß-galactosidase protein, to each sample makes monitoring and quality control possible throughout the entire process. These hardware and process improvements provide good reproducibility and improve intra-assay and inter-assay precision (CV of less than 20%) for LC-MS based protein and peptide quantification. The entire workflow for digesting 96 samples in a 96-well plate can be completed in approximately 5 hours.

The importance of including imperfect detection models in eDNA experimental design.

Environmental DNA (eDNA) is DNA that has been isolated from field samples, and it is increasingly used to infer the presence or absence of particular species in an ecosystem. However, the combination of sampling procedures and subsequent molecular amplification of eDNA can lead to spurious results. As such, it is imperative that eDNA studies include a statistical framework for interpreting eDNA presence/absence data. We reviewed published literature for studies that utilized eDNA where the species density was known and compared the probability of detecting the focal species to the sampling and analysis protocols. Although biomass of the target species and the volume per sample did not impact detectability, the number of field replicates and number of samples from each replicate were positively related to detection. Additionally, increased number of PCR replicates and increased primer specificity significantly increased detectability. Accordingly, we advocate for increased use of occupancy modelling as a method to incorporate effects of sampling effort and PCR sensitivity in eDNA study design. Based on simulation results and the hierarchical nature of occupancy models, we suggest that field replicates, as opposed to molecular replicates, result in better detection probabilities of target species.

Horizontal gene transfer in schistosomes: A critical assessment.

Horizontal gene transfer (HGT), the movement of genetic material between distinct evolutionary lineages, has long been known as a principal force of diversification and adaptation of prokaryotes. More recently, genomic and transcriptomic datasets have suggested gene transfers among various eukaryotic taxa (e.g., Porifera, Cnidaria, Nematoda, Arthropoda, Rotifera, Craniata, and Plantae). Although the exact mechanism of eukaryotic HGT is often unknown, host-parasite interactions may provide ample opportunities for HGT. Schistosomes are trematode blood parasites with complex life cycles that have been repeatedly implicated in HGT. We employed molecular, bioinformatic and phylogenetic approaches to critically analyze 13 published reports of direct HGTs between schistosomes and their hosts to better understand host-parasite co-evolution. Our research suggests that reported cases of schistosome-associated HGT may be due to technical artifacts as opposed to biological reality as we were unable to substantiate them. HGT clearly occurs in eukaryotic organisms, but the burden of proof is high and we emphasize the importance of multiple lines of evidence to conclusively document HGT.

The genome sequence of a widespread apex predator, the golden eagle (Aquila chrysaetos).

Biologists routinely use molecular markers to identify conservation units, to quantify genetic connectivity, to estimate population sizes, and to identify targets of selection. Many imperiled eagle populations require such efforts and would benefit from enhanced genomic resources. We sequenced, assembled, and annotated the first eagle genome using DNA from a male golden eagle (Aquila chrysaetos) captured in western North America. We constructed genomic libraries that were sequenced using Illumina technology and assembled the high-quality data to a depth of ∼40x coverage. The genome assembly includes 2,552 scaffolds >10 Kb and 415 scaffolds >1.2 Mb. We annotated 16,571 genes that are involved in myriad biological processes, including such disparate traits as beak formation and color vision. We also identified repetitive regions spanning 92 Mb (∼6% of the assembly), including LINES, SINES, LTR-RTs and DNA transposons. The mitochondrial genome encompasses 17,332 bp and is ∼91% identical to the Mountain Hawk-Eagle (Nisaetus nipalensis). Finally, the data reveal that several anonymous microsatellites commonly used for population studies are embedded within protein-coding genes and thus may not have evolved in a neutral fashion. Because the genome sequence includes ∼800,000 novel polymorphisms, markers can now be chosen based on their proximity to functional genes involved in migration, carnivory, and other biological processes.

Hosts, parasites, and horizontal gene transfer.

Mendelian inheritance transfers genes vertically within lineages, whereas horizontal gene transfer (HGT) moves genetic material between or among lineages. Herein, we explore possible mechanisms of HGT between parasites and their hosts, as their intimate contact affords substantial opportunities for HGT. We review studies of host-parasite HGT, discussing their merits, their shortcomings, and the multiple lines of evidence needed to conclusively document HGT while avoiding false positives. We focus mainly on schistosomes and other parasites with complex life cycles as they provide numerous opportunities for HGT among the parasite, intermediate, and definitive host genomes. We also highlight future research directions that could prove illuminating with regard to the occurrence, prevalence, and overall importance of HGT in host-parasite coevolution.

Of contigs and quagmires: next-generation sequencing pitfalls associated with transcriptomic studies.

Molecular ecologists have good reasons to be excited about the newest DNA/RNA sequencing technologies. However, this exuberance should be tempered with a hefty dose of reality: new sequencing technologies come with significant new challenges. Herein, we offer a brief overview of some practical problems encountered during transcriptomics studies conducted in our laboratory, and of nontrivial issues that prospective practitioners should consider. These include template contamination (e.g. from xenobiotics) and the cutting-room floor problem, whereby most of the data are often unassembled, unannotated and unused. We also highlight computational requirements, including hardware, personnel time and associated skill sets. We are very optimistic about the future of molecular ecology, but we hope this cautionary overview will help neophytes better recognize some key challenges associated with new technologies.