Geographic Origin and Vertical Transmission of Leishmania infantum Parasites in Hunting Hounds, United States.

Vertical transmission of leishmaniasis is common but is difficult to study against the background of pervasive vector transmission. We present genomic data from dogs in the United States infected with Leishmania infantum parasites; these infections have persisted in the apparent absence of vector transmission. We demonstrate that these parasites were introduced from the Old World separately and more recently than L. infantum from South America. The parasite population shows unusual genetics consistent with a lack of meiosis: a high level of heterozygous sites shared across all isolates and no decrease in linkage with genomic distance between variants. Our data confirm that this parasite population has been evolving with little or no sexual reproduction. This demonstration of vertical transmission has profound implications for the population genetics of Leishmania parasites. When investigating transmission in complex natural settings, considering vertical transmission alongside vector transmission is vital.

Systematic identification of genes encoding cell surface and secreted proteins that are essential for in vitro growth and infection in Leishmania donovani.

Leishmaniasis is an infectious disease caused by protozoan parasites belonging to the genus Leishmania for which there are no approved human vaccines. Infections localise to different tissues in a species-specific manner with the visceral form of the disease caused by Leishmania donovani and L. infantum being the most deadly in humans. Although Leishmania spp. parasites are predominantly intracellular, the visceral disease can be prevented in dogs by vaccinating with a complex mixture of secreted products from cultures of L. infantum promastigotes. With the logic that extracellular parasite proteins make good subunit vaccine candidates because they are directly accessible to vaccine-elicited host antibodies, here we attempt to discover proteins that are essential for in vitro growth and host infection with the goal of identifying subunit vaccine candidates. Using an in silico analysis of the Leishmania donovani genome, we identified 92 genes encoding proteins that are predicted to be secreted or externally anchored to the parasite membrane by a single transmembrane region or a GPI anchor. By selecting a transgenic L. donovani parasite that expresses both luciferase and the Cas9 nuclease, we systematically attempted to target all 92 genes by CRISPR genome editing and identified four that were required for in vitro growth. For fifty-five genes, we infected cohorts of mice with each mutant parasite and by longitudinally quantifying parasitaemia with bioluminescent imaging, showed that nine genes had evidence of an attenuated infection although all ultimately established an infection. Finally, we expressed two genes as full-length soluble recombinant proteins and tested them as subunit vaccine candidates in a murine preclinical infection model. Both proteins elicited significant levels of protection against the uncontrolled development of a splenic infection warranting further investigation as subunit vaccine candidates against this deadly infectious tropical disease.

Diversity and Within-Host Evolution of Leishmania donovani from Visceral Leishmaniasis Patients with and without HIV Coinfection in Northern Ethiopia.

Visceral leishmaniasis (VL) is a fatal disease and a growing public health problem in East Africa, where Ethiopia has one of the highest VL burdens. The largest focus of VL in Ethiopia is driven by high prevalence in migrant agricultural workers and associated with a high rate of coinfection with HIV. This coinfection makes VL more difficult to treat successfully and is associated with a high rate of relapse, with VL/HIV patients frequently experiencing many relapses of VL before succumbing to this infection. We present genome-wide data on Leishmania donovani isolates from a longitudinal study of cohorts of VL and VL/HIV patients reporting to a single clinic in Ethiopia. Extensive clinical data allow us to investigate the influence of coinfection and relapse on the populations of parasites infecting these patients. We find that the same parasite population is responsible for both VL and VL/HIV infections and that, in most cases, disease relapse is caused by recrudescence of the population of parasites that caused primary VL. Complex, multiclonal infections are present in both primary and relapse cases, but the infrapopulation of parasites within a patient loses genetic diversity between primary disease presentation and subsequent relapses, presumably due to a population bottleneck induced by treatment. These data suggest that VL/HIV relapses are not caused by genetically distinct parasite infections or by reinfection. Treatment of VL does not lead to sterile cure, and in VL/HIV, the infecting parasites are able to reestablish after clinically successful treatment, leading to repeated relapse of VL. IMPORTANCE Visceral leishmaniasis (VL) is the second largest cause of deaths due to parasite infections and a growing problem in East Africa. In Ethiopia, it is particularly associated with migrant workers moving from regions of nonendemicity for seasonal agricultural work and is frequently found as a coinfection with HIV, which leads to frequent VL relapse following treatment. Insight into the process of relapse in these patients is thus key to controlling the VL epidemic in Ethiopia. We show that there is little genetic differentiation between the parasites infecting HIV-positive and HIV-negative VL patients. Moreover, we provide evidence that relapses are caused by the initially infecting parasite population and that treatment induces a loss of genetic diversity in this population. We propose that restoring functioning immunity and improving antiparasitic treatment may be key in breaking the cycle of relapsing VL in VL/HIV patients.

A Way Straight-Forward for Leishmania Genetics.

Genetic exchange between Leishmania parasites was demonstrated in sandflies over 10 years ago. Louradour et al. have shown in vitro hybridization of two Leishmania tropica isolates, with the potential to remove a major roadblock to using forward genetics in Leishmania, understanding Leishmania reproductive biology, and analyzing gene flow in natural populations.

Genomic analysis of natural intra-specific hybrids among Ethiopian isolates of Leishmania donovani.

Parasites of the genus Leishmania (Kinetoplastida: Trypanosomatidae) cause widespread and devastating human diseases. Visceral leishmaniasis due to Leishmania donovani is endemic in Ethiopia where it has also been responsible for major epidemics. The presence of hybrid genotypes has been widely reported in surveys of natural populations, genetic variation reported in a number of Leishmania species, and the extant capacity for genetic exchange demonstrated in laboratory experiments. However, patterns of recombination and the evolutionary history of admixture that produced these hybrid populations remain unclear. Here, we use whole-genome sequence data to investigate Ethiopian L. donovani isolates previously characterized as hybrids by microsatellite and multi-locus sequencing. To date there is only one previous study on a natural population of Leishmania hybrids based on whole-genome sequences. We propose that these hybrids originate from recombination between two different lineages of Ethiopian L. donovani occurring in the same region. Patterns of inheritance are more complex than previously reported with multiple, apparently independent, origins from similar parents that include backcrossing with parental types. Analysis indicates that hybrids are representative of at least three different histories. Furthermore, isolates were highly polysomic at the level of chromosomes with differences between parasites recovered from a recrudescent infection from a previously treated individual. The results demonstrate that recombination is a significant feature of natural populations and contributes to the growing body of data that shows how recombination, and gene flow, shape natural populations of Leishmania.

Global genome diversity of the Leishmania donovani complex.

Protozoan parasites of the Leishmania donovani complex - L. donovani and L. infantum - cause the fatal disease visceral leishmaniasis. We present the first comprehensive genome-wide global study, with 151 cultured field isolates representing most of the geographical distribution. L. donovani isolates separated into five groups that largely coincide with geographical origin but vary greatly in diversity. In contrast, the majority of L. infantum samples fell into one globally-distributed group with little diversity. This picture is complicated by several hybrid lineages. Identified genetic groups vary in heterozygosity and levels of linkage, suggesting different recombination histories. We characterise chromosome-specific patterns of aneuploidy and identified extensive structural variation, including known and suspected drug resistance loci. This study reveals greater genetic diversity than suggested by geographically-focused studies, provides a resource of genomic variation for future work and sets the scene for a new understanding of the evolution and genetics of the Leishmania donovani complex.

Reconstruction of Haplotype-Blocks Selected during Experimental Evolution.

The genetic analysis of experimentally evolving populations typically relies on short reads from pooled individuals (Pool-Seq). While this method provides reliable allele frequency estimates, the underlying haplotype structure remains poorly characterized. With small population sizes and adaptive variants that start from low frequencies, the interpretation of selection signatures in most Evolve and Resequencing studies remains challenging. To facilitate the characterization of selection targets, we propose a new approach that reconstructs selected haplotypes from replicated time series, using Pool-Seq data. We identify selected haplotypes through the correlated frequencies of alleles carried by them. Computer simulations indicate that selected haplotype-blocks of several Mb can be reconstructed with high confidence and low error rates, even when allele frequencies change only by 20% across three replicates. Applying this method to real data from D. melanogaster populations adapting to a hot environment, we identify a selected haplotype-block of 6.93 Mb. We confirm the presence of this haplotype-block in evolved populations by experimental haplotyping, demonstrating the power and accuracy of our haplotype reconstruction from Pool-Seq data. We propose that the combination of allele frequency estimates with haplotype information will provide the key to understanding the dynamics of adaptive alleles.

High rates of phasing errors in highly polymorphic species with low levels of linkage disequilibrium.

Short read sequencing of diploid individuals does not permit the direct inference of the sequence on each of the two homologous chromosomes. Although various phasing software packages exist, they were primarily tailored for and tested on human data, which differ from other species in factors that influence phasing, such as SNP density, amounts of linkage disequilibrium (LD) and sample sizes. Despite becoming increasingly popular for other species, the reliability of phasing in non-human data has not been evaluated to a sufficient extent. We scrutinized the phasing accuracy for Drosophila melanogaster, a species with high polymorphism levels and reduced LD relative to humans. We phased two D. melanogaster populations and compared the results to the known haplotypes. The performance increased with size of the reference panel and was highest when the reference panel and phased individuals were from the same population. Full genomic SNP data and inclusion of sequence read information also improved phasing. Despite humans and Drosophila having similar switch error rates between polymorphic sites, the distances between switch errors were much shorter in Drosophila with only fragments <300-1500 bp being correctly phased with ≥95% confidence. This suggests that the higher SNP density cannot compensate for the higher recombination rate in D. melanogaster. Furthermore, we show that populations that have gone through demographic events such as bottlenecks can be phased with higher accuracy. Our results highlight that statistically phased data are particularly error prone in species with large population sizes or populations lacking suitable reference panels.

Patterns of linkage disequilibrium and long range hitchhiking in evolving experimental Drosophila melanogaster populations.

Whole-genome resequencing of experimental populations evolving under a specific selection regime has become a popular approach to determine genotype-phenotype maps and understand adaptation to new environments. Despite its conceptual appeal and success in identifying some causative genes, it has become apparent that many studies suffer from an excess of candidate loci. Several explanations have been proposed for this phenomenon, but it is clear that information about the linkage structure during such experiments is needed. Until now only Pool-Seq (whole-genome sequencing of pools of individuals) data were available, which do not provide sufficient information about the correlation between linked sites. We address this problem in two complementary analyses of three replicate Drosophila melanogaster populations evolving to a new hot temperature environment for almost 70 generations. In the first analysis, we sequenced 58 haploid genomes from the founder population and evolved flies at generation 67. We show that during the experiment linkage disequilibrium (LD) increased almost uniformly over much greater distances than typically seen in Drosophila. In the second analysis, Pool-Seq time series data of the three replicates were combined with haplotype information from the founder population to follow blocks of initial haplotypes over time. We identified 17 selected haplotype-blocks that started at low frequencies in the base population and increased in frequency during the experiment. The size of these haplotype-blocks ranged from 0.082 to 4.01 Mb. Moreover, between 42% and 46% of the top candidate single nucleotide polymorphisms from the comparison of founder and evolved populations fell into the genomic region covered by the haplotype-blocks. We conclude that LD in such rising haplotype-blocks results in long range hitchhiking over multiple kilobase-sized regions. LD in such haplotype-blocks is therefore a major factor contributing to an excess of candidate loci. Although modifications of the experimental design may help to reduce the hitchhiking effect and allow for more precise mapping of causative variants, we also note that such haplotype-blocks might be well suited to study the dynamics of selected genomic regions during experimental evolution studies.

Genome-wide transcriptomic responses of the seagrasses Zostera marina and Nanozostera noltii under a simulated heatwave confirm functional types.

Genome-wide transcription analysis between related species occurring in overlapping ranges can provide insights into the molecular basis underlying different ecological niches. The co-occurring seagrass species, Zostera marina and Nanozostera noltii, are found in marine coastal environments throughout the northern hemisphere. Z. marina is often dominant in subtidal environments and subjected to fewer temperature extremes compared to the predominately intertidal and more stress-tolerant N. noltii. We exposed plants of both species to a realistic heat wave scenario in a common-stress-garden experiment. Using RNA-seq (~7million reads/library), four Z. marina and four N. noltii libraries were compared representing northern (Denmark) and southern (Italy) locations within the co-occurring range of the species' European distribution. A total of 8977 expressed genes were identified, of which 78 were directly related to heat stress. As predicted, both species were negatively affected by the heat wave, but showed markedly different molecular responses. In Z. marina the heat response was similar across locations in response to the heatwave at 26°C, with a complex response in functions related to protein folding, synthesis of ribosomal chloroplast proteins, proteins involved in cell wall modification and heat shock proteins (HSPs). In N. noltii the heat response markedly differed between locations, while HSP genes were not induced in either population. Our results suggest that as coastal seawater temperatures increase, Z. marina will disappear along its southern most ranges, whereas N. noltii will continue to move north. As a consequence, sub- and intertidal habitat partitioning may weaken in more northern regions because the higher thermal tolerance of N. noltii provides a competitive advantage in both habitats. Although previous studies have focused on HSPs, the present study clearly demonstrates that a broader examination of stress related genes is necessary.

Massive habitat-specific genomic response in D. melanogaster populations during experimental evolution in hot and cold environments.

Experimental evolution in combination with whole-genome sequencing (evolve and resequence [E&R]) is a promising approach to define the genotype-phenotype map and to understand adaptation in evolving populations. Many previous studies have identified a large number of putative selected sites (i.e., candidate loci), but it remains unclear to what extent these loci are genuine targets of selection or experimental noise. To address this question, we exposed the same founder population to two different selection regimes-a hot environment and a cold environment-and quantified the genomic response in each. We detected large numbers of putative selected loci in both environments, albeit with little overlap between the two sets of candidates, indicating that most resulted from habitat-specific selection. By quantifying changes across multiple independent biological replicates, we demonstrate that most of the candidate SNPs were false positives that were linked to selected sites over distances much larger than the typical linkage disequilibrium range of Drosophila melanogaster. We show that many of these mid- to long-range associations were attributable to large segregating inversions and confirm by computer simulations that such patterns could be readily replicated when strong selection acts on rare haplotypes. In light of our findings, we outline recommendations to improve the performance of future Drosophila E&R studies which include using species with negligible inversion loads, such as D. mauritiana and D. simulans, instead of D. melanogaster.

Allelic imbalance metre (Allim), a new tool for measuring allele-specific gene expression with RNA-seq data.

Estimating differences in gene expression among alleles is of high interest for many areas in biology and medicine. Here, we present a user-friendly software tool, Allim, to estimate allele-specific gene expression. Because mapping bias is a major problem for reliable estimates of allele-specific gene expression using RNA-seq, Allim combines two different strategies to account for the mapping biases. In order to reduce the mapping bias, Allim first generates a polymorphism-aware reference genome that accounts for the sequence variation between the alleles. Then, a sequence-specific simulation tool estimates the residual mapping bias. Statistical tests for allelic imbalance are provided that can be used with the bias corrected RNA-seq data.

Intra-specific regulatory variation in Drosophila pseudoobscura.

It is generally accepted that gene regulation serves an important role in determining the phenotype. To shed light on the evolutionary forces operating on gene regulation, previous studies mainly focused on the expression differences between species and their inter-specific hybrids. Here, we use RNA-Seq to study the intra-specific distribution of cis- and trans-regulatory variation in Drosophila pseudoobscura. Consistent with previous results, we find almost twice as many genes (26%) with significant trans-effects than genes with significant cis-effects (18%). While this result supports the previous suggestion of a larger mutational target of trans-effects, we also show that trans-effects may be subjected to purifying selection. Our results underline the importance of intra-specific analyses for the understanding of the evolution of gene expression.

Identifying core features of adaptive metabolic mechanisms for chronic heat stress attenuation contributing to systems robustness.

The contribution of metabolism to heat stress may play a significant role in defining robustness and recovery of systems; either by providing the energy and metabolites required for cellular homeostasis, or through the generation of protective osmolytes. However, the mechanisms by which heat stress attenuation could be adapted through metabolic processes as a stabilizing strategy against thermal stress are still largely unclear. We address this issue through metabolomic and transcriptomic profiles for populations along a thermal cline where two seagrass species, Zostera marina and Zostera noltii, were found in close proximity. Significant changes captured by these profile comparisons could be detected, with a larger response magnitude observed in northern populations to heat stress. Sucrose, fructose, and myo-inositol were identified to be the most responsive of the 29 analyzed organic metabolites. Many key enzymes in the Calvin cycle, glycolysis and pentose phosphate pathways also showed significant differential expression. The reported comparison suggests that adaptive mechanisms are involved through metabolic pathways to dampen the impacts of heat stress, and interactions between the metabolome and proteome should be further investigated in systems biology to understand robust design features against abiotic stress.

Transcriptomic resilience to global warming in the seagrass Zostera marina, a marine foundation species.

Large-scale transcription profiling via direct cDNA sequencing provides important insights as to how foundation species cope with increasing climatic extremes predicted under global warming. Species distributed along a thermal cline, such as the ecologically important seagrass Zostera marina, provide an opportunity to assess temperature effects on gene expression as a function of their long-term adaptation to heat stress. We exposed a southern and northern European population of Zostera marina from contrasting thermal environments to a realistic heat wave in a common-stress garden. In a fully crossed experiment, eight cDNA libraries, each comprising ~125 000 reads, were obtained during and after a simulated heat wave, along with nonstressed control treatments. Although gene-expression patterns during stress were similar in both populations and were dominated by classical heat-shock proteins, transcription profiles diverged after the heat wave. Gene-expression patterns in southern genotypes returned to control values immediately, but genotypes from the northern site failed to recover and revealed the induction of genes involved in protein degradation, indicating failed metabolic compensation to high sea-surface temperature. We conclude that the return of gene-expression patterns during recovery provides critical information on thermal adaptation in aquatic habitats under climatic stress. As a unifying concept for ecological genomics, we propose transcriptomic resilience, analogous to ecological resilience, as an important measure to predict the tolerance of individuals and hence the fate of local populations in the face of global warming.

Comprehensive transcriptome analysis of the highly complex Pisum sativum genome using next generation sequencing.

BACKGROUND: The garden pea, Pisum sativum, is among the best-investigated legume plants and of significant agro-commercial relevance. Pisum sativum has a large and complex genome and accordingly few comprehensive genomic resources exist. RESULTS: We analyzed the pea transcriptome at the highest possible amount of accuracy by current technology. We used next generation sequencing with the Roche/454 platform and evaluated and compared a variety of approaches, including diverse tissue libraries, normalization, alternative sequencing technologies, saturation estimation and diverse assembly strategies. We generated libraries from flowers, leaves, cotyledons, epi- and hypocotyl, and etiolated and light treated etiolated seedlings, comprising a total of 450 megabases. Libraries were assembled into 324,428 unigenes in a first pass assembly.A second pass assembly reduced the amount to 81,449 unigenes but caused a significant number of chimeras. Analyses of the assemblies identified the assembly step as a major possibility for improvement. By recording frequencies of Arabidopsis orthologs hit by randomly drawn reads and fitting parameters of the saturation curve we concluded that sequencing was exhaustive. For leaf libraries we found normalization allows partial recovery of expression strength aside the desired effect of increased coverage. Based on theoretical and biological considerations we concluded that the sequence reads in the database tagged the vast majority of transcripts in the aerial tissues. A pathway representation analysis showed the merits of sampling multiple aerial tissues to increase the number of tagged genes. All results have been made available as a fully annotated database in fasta format. CONCLUSIONS: We conclude that the approach taken resulted in a high quality - dataset which serves well as a first comprehensive reference set for the model legume pea. We suggest future deep sequencing transcriptome projects of species lacking a genomics backbone will need to concentrate mainly on resolving the issues of redundancy and paralogy during transcriptome assembly.