Viral population estimation using pyrosequencing

The diversity of virus populations within single infected hosts presents a major difficulty for the natural immune response as well as for vaccine design and antiviral drug therapy. Recently developed pyrophosphate based sequencing technologies (pyrosequencing) can be used for quantifying this diversity by ultra-deep sequencing of virus samples. We present computational methods for the analysis of such sequence data and apply these techniques to pyrosequencing data obtained from HIV populations within patients harboring drug resistant virus strains. Our main result is the estimation of the population structure of the sample from the pyrosequencing reads. This inference is based on a statistical approach to error correction, followed by a combinatorial algorithm for constructing a minimal set of haplotypes that explain the data. Using this set of explaining haplotypes, we apply a statistical model to infer the frequencies of the haplotypes in the population via an EM algorithm. We demonstrate that pyrosequencing reads allow for effective population reconstruction by extensive simulations and by comparison to 165 sequences obtained directly from clonal sequencing of four independent, diverse HIV populations. Thus, pyrosequencing can be used for cost-effective estimation of the structure of virus populations, promising new insights into viral evolutionary dynamics and disease control strategies.Comment: 23 pages, 13 figure

Distribution of genetic diversity in wild European populations of prickly lettuce (Lactuca serriola): implications for plant genetic resources management

Genetic variation in Lactuca serriola, the closest wild relative of cultivated lettuce, was studied across Europe from the Czech Republic to the United Kingdom, using three molecular marker systems, simple sequence repeat (SSR, microsatellites), AFLP and nucleotide-binding site (NBS) profiling. The ‘functional’ marker system NBS profiling, targeting disease resistance genes of the NBS/LRR family, did not show marked differences in genetic diversity parameters to the other systems. The autogamy of the species resulted in low observed heterozygosity and high population differentiation. Intra-population variation ranged from complete homogeneity to nearly complete heterogeneity. The highest genetic diversity was found in central Europe. The SSR results were compared to SSR variation screened earlier in the lettuce collection of the Centre for Genetic Resources, the Netherlands (CGN). In the UK, practically only a single SSR genotype was found. This genotype together with a few other common SSR genotypes comprised a large part of the plants sampled on the continent. Among the ten most frequent SSR genotypes observed, eight were already present in the CGN collection. Overall, the CGN collection appears to already have a fair representation of genetic variation from NW Europe. The results are discussed in relation to sampling strategies for improving genebank collections of crop wild relatives

Impact of sample preservation and manipulation on insect gut microbiome profiling : a test case with fruit flies (Diptera, Tephritidae)

High-throughput sequencing (HTS) techniques are of great value for the investigation of microbial communities, and have been extensively used to study the gut microbiome. While most studies focus on the human gut, many others have investigated insects. However, because of the rapid spread of HTS techniques, a lot of variation exists in the protocols for sample preparation. In the present study, we investigated the impact of two widely adopted sample-processing procedures preceding library preparation, i.e., preservation of insect tissue in 70% ethanol (EtOH) and sample dissection. We used the fruit fly Ceratitis capitata (Diptera: Tephritidae) as a model organism and set up two experiments, one comparing the effects of sample manipulation and preservation across life stages and the other across fruit samples from different sources. The results of this study showed no major effects of dissection on the outcome of HTS. However, EtOH preservation did have effects on the recovered gut microbiome, the main effect being a significant reduction of the dominant genus, Providencia, in EtOH-preserved samples. Less abundant bacterial groups were also affected resulting in altered microbial profiles obtained from samples preserved in 70% EtOH. These results have important implications for the planning of future studies and when comparing studies that used different sample preparation protocols

Limited antigenic diversity of Plasmodium falciparum apical membrane antigen 1 supports the development of effective multi-allele vaccines

BackgroundPolymorphism in antigens is a common mechanism for immune evasion used by many important pathogens, and presents major challenges in vaccine development. In malaria, many key immune targets and vaccine candidates show substantial polymorphism. However, knowledge on antigenic diversity of key antigens, the impact of polymorphism on potential vaccine escape, and how sequence polymorphism relates to antigenic differences is very limited, yet crucial for vaccine development. Plasmodium falciparum apical membrane antigen 1 (AMA1) is an important target of naturally-acquired antibodies in malaria immunity and a leading vaccine candidate. However, AMA1 has extensive allelic diversity with more than 60 polymorphic amino acid residues and more than 200 haplotypes in a single population. Therefore, AMA1 serves as an excellent model to assess antigenic diversity in malaria vaccine antigens and the feasibility of multi-allele vaccine approaches. While most previous research has focused on sequence diversity and antibody responses in laboratory animals, little has been done on the cross-reactivity of human antibodies.MethodsWe aimed to determine the extent of antigenic diversity of AMA1, defined by reactivity with human antibodies, and to aid the identification of specific alleles for potential inclusion in a multi-allele vaccine. We developed an approach using a multiple-antigen-competition enzyme-linked immunosorbent assay (ELISA) to examine cross-reactivity of naturally-acquired antibodies in Papua New Guinea and Kenya, and related this to differences in AMA1 sequence.ResultsWe found that adults had greater cross-reactivity of antibodies than children, although the patterns of cross-reactivity to alleles were the same. Patterns of antibody cross-reactivity were very similar between populations (Papua New Guinea and Kenya), and over time. Further, our results show that antigenic diversity of AMA1 alleles is surprisingly restricted, despite extensive sequence polymorphism. Our findings suggest that a combination of three different alleles, if selected appropriately, may be sufficient to cover the majority of antigenic diversity in polymorphic AMA1 antigens. Antigenic properties were not strongly related to existing haplotype groupings based on sequence analysis.ConclusionsAntigenic diversity of AMA1 is limited and a vaccine including a small number of alleles might be sufficient for coverage against naturally-circulating strains, supporting a multi-allele approach for developing polymorphic antigens as malaria vaccines

Natural variation in abiotic stress responsive gene expression and local adaptation to climate in Arabidopsis thaliana.

Gene expression varies widely in natural populations, yet the proximate and ultimate causes of this variation are poorly known. Understanding how variation in gene expression affects abiotic stress tolerance, fitness, and adaptation is central to the field of evolutionary genetics. We tested the hypothesis that genes with natural genetic variation in their expression responses to abiotic stress are likely to be involved in local adaptation to climate in Arabidopsis thaliana. Specifically, we compared genes with consistent expression responses to environmental stress (expression stress responsive, "eSR") to genes with genetically variable responses to abiotic stress (expression genotype-by-environment interaction, "eGEI"). We found that on average genes that exhibited eGEI in response to drought or cold had greater polymorphism in promoter regions and stronger associations with climate than those of eSR genes or genomic controls. We also found that transcription factor binding sites known to respond to environmental stressors, especially abscisic acid responsive elements, showed significantly higher polymorphism in drought eGEI genes in comparison to eSR genes. By contrast, eSR genes tended to exhibit relatively greater pairwise haplotype sharing, lower promoter diversity, and fewer nonsynonymous polymorphisms, suggesting purifying selection or selective sweeps. Our results indicate that cis-regulatory evolution and genetic variation in stress responsive gene expression may be important mechanisms of local adaptation to climatic selective gradients

Local adaptation drives the diversification of effectors in the fungal wheat pathogen Parastagonospora nodorum in the United States

Filamentous fungi rapidly evolve in response to environmental selection pressures in part due to their genomic plasticity. Parastagonospora nodorum, a fungal pathogen of wheat and causal agent of septoria nodorum blotch, responds to selection pressure exerted by its host, influencing the gain, loss, or functional diversification of virulence determinants, known as effector genes. Whole genome resequencing of 197 P. nodorum isolates collected from spring, durum, and winter wheat production regions of the United States enabled the examination of effector diversity and genomic regions under selection specific to geographically discrete populations. 1,026,859 SNPs/InDels were used to identify novel loci, as well as SnToxA and SnTox3 as factors in disease. Genes displaying presence/absence variation, predicted effector genes, and genes localized on an accessory chromosome had significantly higher pN/pS ratios, indicating a higher rate of sequence evolution. Population structure analyses indicated two P. nodorum populations corresponding to the Upper Midwest (Population 1) and Southern/Eastern United States (Population 2). Prevalence of SnToxA varied greatly between the two populations which correlated with presence of the host sensitivity gene Tsn1 in the most prevalent cultivars in the corresponding regions. Additionally, 12 and 5 candidate effector genes were observed to be under diversifying selection among isolates from Population 1 and 2, respectively, but under purifying selection or neutrally evolving in the opposite population. Selective sweep analysis revealed 10 and 19 regions that had recently undergone positive selection in Population 1 and 2, respectively, involving 92 genes in total. When comparing genes with and without presence/absence variation, those genes exhibiting this variation were significantly closer to transposable elements. Taken together, these results indicate that P. nodorum is rapidly adapting to distinct selection pressures unique to spring and winter wheat production regions by rapid adaptive evolution and various routes of genomic diversification, potentially facilitated through transposable element activity

Outstanding intraindividual genetic diversity in fissiparous planarians (Dugesia, Platyhelminthes) with facultative sex.

Predicted genetic consequences of asexuality include high intraindividual genetic diversity (i.e., the Meselson effect) and accumulation of deleterious mutations (i.e., Muller’s Ratchet), among others. These consequences have been largely studied in parthenogenetic organisms, but studies on fissiparous species are scarce. Differing from parthenogens, fissiparous organisms inherit part of the soma of the progenitor, including somatic mutations. Thus, in the long term, fissiparous reproduction may also result in genetic mosaicism, besides the presence of the Meselson effect and Muller’s Ratchet. Dugesiidae planarians show outstanding regeneration capabilities, allowing them to naturally reproduce by fission, either strictly or combined with sex (facultative). Therefore, they are an ideal model to analyze the genetic footprint of fissiparous reproduction, both when it is alternated with sex and when it is the only mode of reproduction

Ensemble Learning for Free with Evolutionary Algorithms ?

Evolutionary Learning proceeds by evolving a population of classifiers, from which it generally returns (with some notable exceptions) the single best-of-run classifier as final result. In the meanwhile, Ensemble Learning, one of the most efficient approaches in supervised Machine Learning for the last decade, proceeds by building a population of diverse classifiers. Ensemble Learning with Evolutionary Computation thus receives increasing attention. The Evolutionary Ensemble Learning (EEL) approach presented in this paper features two contributions. First, a new fitness function, inspired by co-evolution and enforcing the classifier diversity, is presented. Further, a new selection criterion based on the classification margin is proposed. This criterion is used to extract the classifier ensemble from the final population only (Off-line) or incrementally along evolution (On-line). Experiments on a set of benchmark problems show that Off-line outperforms single-hypothesis evolutionary learning and state-of-art Boosting and generates smaller classifier ensembles

The development of a novel SNP genotyping assay to differentiate cacao clones

In this study, a double-mismatch allele-specific (DMAS) qPCR SNP genotyping method has been designed, tested and validated specifically for cacao, using 65 well annotated international cacao reference accessions retrieved from the Center for Forestry Research and Technology Transfer (CEFORTT) and the International Cocoa Quarantine Centre (ICQC). In total, 42 DMAS-qPCR SNP genotyping assays have been validated, with a 98.05% overall efficiency in calling the correct genotype. In addition, the test allowed for the identification of 15.38% off-types and two duplicates, highlighting the problem of mislabeling in cacao collections and the need for conclusive genotyping assays. The developed method showed on average a high genetic diversity (He = 0.416) and information index (I = 0.601), making it applicable to assess intra-population variation. Furthermore, only the 13 most informative markers were needed to achieve maximum differentiation. This simple, effective method provides robust and accurate genotypic data which allows for more efficient resource management (e.g. tackling mislabeling, conserving valuable genetic material, parentage analysis, genetic diversity studies), thus contributing to an increased knowledge on the genetic background of cacao worldwide. Notably, the described method can easily be integrated in other laboratories for a wide range of objectives and organisms