Genome-wide SNP genotyping highlights the role of natural selection in Plasmodium falciparum population divergence

Background

The malaria parasite Plasmodium falciparum exhibits abundant genetic diversity, and this diversity is key to its success as a pathogen. Previous efforts to study genetic diversity in P. falciparum have begun to elucidate the demographic history of the species, as well as patterns of population structure and patterns of linkage disequilibrium within its genome. Such studies will be greatly enhanced by new genomic tools and recent large-scale efforts to map genomic variation. To that end, we have developed a high throughput single nucleotide polymorphism (SNP) genotyping platform for P. falciparum.

Results

Using an Affymetrix 3,000 SNP assay array, we found roughly half the assays (1,638) yielded high quality, 100% accurate genotyping calls for both major and minor SNP alleles. Genotype data from 76 global isolates confirm significant genetic differentiation among continental populations and varying levels of SNP diversity and linkage disequilibrium according to geographic location and local epidemiological factors. We further discovered that nonsynonymous and silent (synonymous or noncoding) SNPs differ with respect to within-population diversity, inter-population differentiation, and the degree to which allele frequencies are correlated between populations.

Conclusions

The distinct population profile of nonsynonymous variants indicates that natural selection has a significant influence on genomic diversity in P. falciparum, and that many of these changes may reflect functional variants deserving of follow-up study. Our analysis demonstrates the potential for new high-throughput genotyping technologies to enhance studies of population structure, natural selection, and ultimately enable genome-wide association studies in P. falciparum to find genes underlying key phenotypic traits.     

Genetics of chloroquine-resistant malaria: a haplotypic view

The development and rapid spread of chloroquine resistance (CQR) in Plasmodium falciparum have triggered the identification of several genetic target(s) in the P. falciparum genome. In particular, mutations in the Pfcrt gene, specifically, K76T and mutations in three other amino acids in the region adjoining K76 (residues 72, 74, 75 and 76), are considered to be highly related to CQR. These various mutations form several different haplotypes and Pfcrt gene polymorphisms and the global distribution of the different CQR- Pfcrt haplotypes in endemic and non-endemic regions of P. falciparum malaria have been the subject of extensive study. Despite the fact that the Pfcrt gene is considered to be the primary CQR gene in P. falciparum , several studies have suggested that this may not be the case. Furthermore, there is a poor correlation between the evolutionary implications of the Pfcrt haplotypes and the inferred migration of CQR P. falciparum based on CQR epidemiological surveillance data. The present paper aims to clarify the existing knowledge on the genetic basis of the different CQR- Pfcrt haplotypes that are prevalent in worldwide populations based on the published literature and to analyse the data to generate hypotheses on the genetics and evolution of CQR malaria.     

Genome-wide analysis of the role of GlnR in Streptomyces venezuelae provides new insights into global nitrogen regulation in actinomycetes

Background

Knowledge of the origins, distribution, and inheritance of variation in the malaria parasite (Plasmodium falciparum) genome is crucial for understanding its evolution; however the 81% (A+T) genome poses challenges to high-throughput sequencing technologies. We explore the viability of the Roche 454 Genome Sequencer FLX (GS FLX) high throughput sequencing technology for both whole genome sequencing and fine-resolution characterization of genetic exchange in malaria parasites.

Results

We present a scheme to survey recombination in the haploid stage genomes of two sibling parasite clones, using whole genome pyrosequencing that includes a sliding window approach to predict recombination breakpoints. Whole genome shotgun (WGS) sequencing generated approximately 2 million reads, with an average read length of approximately 300 bp. De novo assembly using a combination of WGS and 3 kb paired end libraries resulted in contigs ≤ 34 kb. More than 8,000 of the 24,599 SNP markers identified between parents were genotyped in the progeny, resulting in a marker density of approximately 1 marker/3.3 kb and allowing for the detection of previously unrecognized crossovers (COs) and many non crossover (NCO) gene conversions throughout the genome.

Conclusions

By sequencing the 23 Mb genomes of two haploid progeny clones derived from a genetic cross at more than 30× coverage, we captured high resolution information on COs, NCOs and genetic variation within the progeny genomes. This study is the first to resequence progeny clones to examine fine structure of COs and NCOs in malaria parasites.     

Population genetic analysis of Plasmodium falciparum parasites using a customized Illumina GoldenGate genotyping assay

The diversity in the Plasmodium falciparum genome can be used to explore parasite population dynamics, with practical applications to malaria control. The ability to identify the geographic origin and trace the migratory patterns of parasites with clinically important phenotypes such as drug resistance is particularly relevant. With increasing single-nucleotide polymorphism (SNP) discovery from ongoing Plasmodium genome sequencing projects, a demand for high SNP and sample throughput genotyping platforms for large-scale population genetic studies is required. Low parasitaemias and multiple clone infections present a number of challenges to genotyping P. falciparum. We addressed some of these issues using a custom 384-SNP Illumina GoldenGate assay on P. falciparum DNA from laboratory clones (long-term cultured adapted parasite clones), short-term cultured parasite isolates and clinical (non-cultured isolates) samples from East and West Africa, Southeast Asia and Oceania. Eighty percent of the SNPs (n = 306) produced reliable genotype calls on samples containing as little as 2 ng of total genomic DNA and on whole genome amplified DNA. Analysis of artificial mixtures of laboratory clones demonstrated high genotype calling specificity and moderate sensitivity to call minor frequency alleles. Clear resolution of geographically distinct populations was demonstrated using Principal Components Analysis (PCA), and global patterns of population genetic diversity were consistent with previous reports. These results validate the utility of the platform in performing population genetic studies of P. falciparum.     

Allelic Diversity and Geographical Distribution of the Gene Encoding Plasmodium falciparum Merozoite Surface Protein-3 in Thailand

Merozoite surface proteins (MSPs) of malaria parasites play critical roles during the erythrocyte invasion and so are potential candidates for malaria vaccine development. However, because MSPs are often under strong immune selection, they can exhibit extensive genetic diversity. The gene encoding the merozoite surface protein-3 (MSP-3) of Plasmodium falciparum displays 2 allelic types, K1 and 3D7. In Thailand, the allelic frequency of the P. falciparum msp-3 gene was evaluated in a single P. falciparum population in Tak at the Thailand and Myanmar border. However, no study has yet looked at the extent of genetic diversity of the msp-3 gene in P. falciparum populations in other localities. Here, we genotyped the msp-3 alleles of 63 P. falciparum samples collected from 5 geographical populations along the borders of Thailand with 3 neighboring countries (Myanmar, Laos, and Cambodia). Our study indicated that the K1 and 3D7 alleles coexisted, but at different proportions in different Thai P. falciparum populations. K1 was more prevalent in populations at the Thailand-Myanmar and Thailand-Cambodia borders, whilst 3D7 was more prevalent at the Thailand-Laos border. Global analysis of the msp-3 allele frequencies revealed that proportions of K1 and 3D7 alleles of msp-3 also varied in different continents, suggesting the divergence of malaria parasite populations. In conclusion, the variation in the msp-3 allelic patterns of P. falciparum in Thailand provides fundamental knowledge for inferring the P. falciparum population structure and for the best design of msp-3 based malaria vaccines.     

Identification of genome-wide single nucleotide polymorphisms in allopolyploid crop Brassica napus

Background

Single nucleotide polymorphisms (SNPs) are the most common type of genetic variation. Identification of large numbers of SNPs is helpful for genetic diversity analysis, map-based cloning, genome-wide association analyses and marker-assisted breeding. Recently, identifying genome-wide SNPs in allopolyploid Brassica napus (rapeseed, canola) by resequencing many accessions has become feasible, due to the availability of reference genomes of Brassica rapa (2n = AA) and Brassica oleracea (2n = CC), which are the progenitor species of B. napus (2n = AACC). Although many SNPs in B. napus have been released, the objective in the present study was to produce a larger, more informative set of SNPs for large-scale and efficient genotypic screening. Hence, short-read genome sequencing was conducted on ten elite B. napus accessions for SNP discovery. A subset of these SNPs was randomly selected for sequence validation and for genotyping efficiency testing using the Illumina GoldenGate assay.

Results

A total of 892,536 bi-allelic SNPs were discovered throughout the B. napus genome. A total of 36,458 putative amino acid variants were located in 13,552 protein-coding genes, which were predicted to have enriched binding and catalytic activity as a result. Using the GoldenGate genotyping platform, 94 of 96 SNPs sampled could effectively distinguish genotypes of 130 lines from two mapping populations, with an average call rate of 92%.

Conclusions

Despite the polyploid nature of B. napus, nearly 900,000 simple SNPs were identified by whole genome resequencing. These SNPs were predicted to be effective in high-throughput genotyping assays (51% polymorphic SNPs, 92% average call rate using the GoldenGate assay, leading to an estimated >450 000 useful SNPs). Hence, the development of a much larger genotyping array of informative SNPs is feasible. SNPs identified in this study to cause non-synonymous amino acid substitutions can also be utilized to directly identify causal genes in association studies.     

Whole Mitochondrial Genome Sequence of an Indian Plasmodium falciparum Field Isolate

Mitochondrial genome sequence of malaria parasites has served as a potential marker for inferring evolutionary history of the Plasmodium genus. In Plasmodium falciparum, the mitochondrial genome sequences from around the globe have provided important evolutionary understanding, but no Indian sequence has yet been utilized. We have sequenced the whole mitochondrial genome of a single P. falciparum field isolate from India using novel primers and compared with the 3D7 reference sequence and 1 previously reported Indian sequence. While the 2 Indian sequences were highly divergent from each other, the presently sequenced isolate was highly similar to the reference 3D7 strain.     

A Locally Acquired Falciparum Malaria via Nosocomial Transmission in Korea

A 57-year old man who was admitted to an emergency room of a tertiary hospital with hemoptysis developed malarial fever 19 days later and then died from severe falciparum malaria 2 days later. He had not traveled outside of Korea for over 30 years. Through intensive interviews and epidemiological surveys, we found that a foreign patient with a recent history of travel to Africa was transferred to the same hospital with severe falciparum malaria. We confirmed through molecular genotyping of the MSP-1 gene that Plasmodium falciparum genotypes of the 2 patients were identical. It is suggested that a breach of standard infection control precautions resulted in this P. falciparum transmission between 2 patients in a hospital environment. This is the first report of a nosocomial transmission of falciparum malaria in Korea.     

SNP Discovery and Linkage Map Construction in Cultivated Tomato

Few intraspecific genetic linkage maps have been reported for cultivated tomato, mainly because genetic diversity within Solanum lycopersicum is much less than that between tomato species. Single nucleotide polymorphisms (SNPs), the most abundant source of genomic variation, are the most promising source of polymorphisms for the construction of linkage maps for closely related intraspecific lines. In this study, we developed SNP markers based on expressed sequence tags for the construction of intraspecific linkage maps in tomato. Out of the 5607 SNP positions detected through in silico analysis, 1536 were selected for high-throughput genotyping of two mapping populations derived from crosses between ‘Micro-Tom’ and either ‘Ailsa Craig’ or ‘M82’. A total of 1137 markers, including 793 out of the 1338 successfully genotyped SNPs, along with 344 simple sequence repeat and intronic polymorphism markers, were mapped onto two linkage maps, which covered 1467.8 and 1422.7 cM, respectively. The SNP markers developed were then screened against cultivated tomato lines in order to estimate the transferability of these SNPs to other breeding materials. The molecular markers and linkage maps represent a milestone in the genomics and genetics, and are the first step toward molecular breeding of cultivated tomato. Information on the DNA markers, linkage maps, and SNP genotypes for these tomato lines is available at http://www.kazusa.or.jp/tomato/.     

Evolution of Genetic Polymorphisms of Plasmodium falciparum Merozoite Surface Protein (PfMSP) in Thailand

Plasmodium falciparum malaria is a major public health problem in Thailand due to the emergence of multidrug resistance. The understanding of genetic diversity of malaria parasites is essential for developing effective drugs and vaccines. The genetic diversity of the merozoite surface protein-1 (PfMSP-1) and merozoite surface protein-2 (PfMSP-2) genes was investigated in a total of 145 P. falciparum isolates collected from Mae Sot District, Tak Province, Thailand during 3 different periods (1997-1999, 2005-2007, and 2009-2010). Analysis of genetic polymorphisms was performed to track the evolution of genetic change of P. falciparum using PCR. Both individual genes and their combination patterns showed marked genetic diversity during the 3 study periods. The results strongly support that P. falciparum isolates in Thailand are markedly diverse and patterns changed with time. These 2 polymorphic genes could be used as molecular markers to detect multiple clone infections and differentiate recrudescence from reinfection in P. falciparum isolates in Thailand.     

Development and evaluation of a 9K SNP array for peach by internationally coordinated SNP detection and validation in breeding germplasm

Although a large number of single nucleotide polymorphism (SNP) markers covering the entire genome are needed to enable molecular breeding efforts such as genome wide association studies, fine mapping, genomic selection and marker-assisted selection in peach [Prunus persica (L.) Batsch] and related Prunus species, only a limited number of genetic markers, including simple sequence repeats (SSRs), have been available to date. To address this need, an international consortium (The International Peach SNP Consortium; IPSC) has pursued a coordinated effort to perform genome-scale SNP discovery in peach using next generation sequencing platforms to develop and characterize a high-throughput Illumina Infinium® SNP genotyping array platform. We performed whole genome re-sequencing of 56 peach breeding accessions using the Illumina and Roche/454 sequencing technologies. Polymorphism detection algorithms identified a total of 1,022,354 SNPs. Validation with the Illumina GoldenGate® assay was performed on a subset of the predicted SNPs, verifying ∼75% of genic (exonic and intronic) SNPs, whereas only about a third of intergenic SNPs were verified. Conservative filtering was applied to arrive at a set of 8,144 SNPs that were included on the IPSC peach SNP array v1, distributed over all eight peach chromosomes with an average spacing of 26.7 kb between SNPs. Use of this platform to screen a total of 709 accessions of peach in two separate evaluation panels identified a total of 6,869 (84.3%) polymorphic SNPs.The almost 7,000 SNPs verified as polymorphic through extensive empirical evaluation represent an excellent source of markers for future studies in genetic relatedness, genetic mapping, and dissecting the genetic architecture of complex agricultural traits. The IPSC peach SNP array v1 is commercially available and we expect that it will be used worldwide for genetic studies in peach and related stone fruit and nut species.     

Population-level genome-wide STR discovery and validation for population structure and genetic diversity assessment of Plasmodium species

Short tandem repeats (STRs) are highly informative genetic markers that have been used extensively in population genetics analysis. They are an important source of genetic diversity and can also have functional impact. Despite the availability of bioinformatic methods that permit large-scale genome-wide genotyping of STRs from whole genome sequencing data, they have not previously been applied to sequencing data from large collections of malaria parasite field samples. Here, we have genotyped STRs using HipSTR in more than 3,000 Plasmodium falciparum and 174 Plasmodium vivax published whole-genome sequence data from samples collected across the globe. High levels of noise and variability in the resultant callset necessitated the development of a novel method for quality control of STR genotype calls. A set of high-quality STR loci (6,768 from P. falciparum and 3,496 from P. vivax) were used to study Plasmodium genetic diversity, population structures and genomic signatures of selection and these were compared to genome-wide single nucleotide polymorphism (SNP) genotyping data. In addition, the genome-wide information about genetic variation and other characteristics of STRs in P. falciparum and P. vivax have been available in an interactive web-based R Shiny application PlasmoSTR (https://github.com/bahlolab/PlasmoSTR).     

Characterization of 18 SNP markers for sperm whale (Physeter macrocephalus)

We report the characterization of 18 new single nucleotide polymorphism (SNP) markers for an endangered species, the sperm whale (Physeter macrocephalus), developed using a targeted gene approach. SNP markers were derived from autosomal regions of the genome using primers originally characterized for genome mapping in other mammals. These SNP markers are the first to be designed for genotyping sperm whale populations and will provide a necessary addition to the genetic tools employed for understanding population structure on a global scale and for developing a conservation management strategy for this endangered species.     

Evaluation of the Accuracy of the EasyTest? Malaria Pf/Pan Ag,a Rapid Diagnostic Test,in Uganda

In recent years, rapid diagnostic tests (RDTs) have been widely used for malaria detection, primarily because of their simple operation, fast results, and straightforward interpretation. The Asan EasyTest™ Malaria Pf/Pan Ag is one of the most commonly used malaria RDTs in several countries, including Korea and India. In this study, we tested the diagnostic performance of this RDT in Uganda to evaluate its usefulness for field diagnosis of malaria in this country. Microscopic and PCR analyses, and the Asan EasyTest™ Malaria Pf/Pan Ag rapid diagnostic test, were performed on blood samples from 185 individuals with suspected malaria in several villages in Uganda. Compared to the microscopic analysis, the sensitivity of the RDT to detect malaria infection was 95.8% and 83.3% for Plasmodium falciparum and non-P. falciparum, respectively. Although the diagnostic sensitivity of the RDT decreased when parasitemia was ≤500 parasites/µl, it showed 96.8% sensitivity (98.4% for P. falciparum and 93.8% for non-P. falciparum) in blood samples with parasitemia ≥100 parasites/µl. The specificity of the RDT was 97.3% for P. falciparum and 97.3% for non-P. falciparum. These results collectively suggest that the accuracy of the Asan EasyTest™ Malaria Pf/Pan Ag makes it an effective point-of-care diagnostic tool for malaria in Uganda.     

Single nucleotide polymorphisms generated by genotyping by sequencing to characterize genome-wide diversity,linkage disequilibrium,and selective sweeps in cultivated watermelon

Background

A large single nucleotide polymorphism (SNP) dataset was used to analyze genome-wide diversity in a diverse collection of watermelon cultivars representing globally cultivated, watermelon genetic diversity. The marker density required for conducting successful association mapping depends on the extent of linkage disequilibrium (LD) within a population. Use of genotyping by sequencing reveals large numbers of SNPs that in turn generate opportunities in genome-wide association mapping and marker-assisted selection, even in crops such as watermelon for which few genomic resources are available. In this paper, we used genome-wide genetic diversity to study LD, selective sweeps, and pairwise F_ST distributions among worldwide cultivated watermelons to track signals of domestication.

Results

We examined 183 Citrullus lanatus var. lanatus accessions representing domesticated watermelon and generated a set of 11,485 SNP markers using genotyping by sequencing. With a diverse panel of worldwide cultivated watermelons, we identified a set of 5,254 SNPs with a minor allele frequency of ≥ 0.05, distributed across the genome. All ancestries were traced to Africa and an admixture of various ancestries constituted secondary gene pools across various continents. A sliding window analysis using pairwise F_ST values was used to resolve selective sweeps. We identified strong selection on chromosomes 3 and 9 that might have contributed to the domestication process. Pairwise analysis of adjacent SNPs within a chromosome as well as within a haplotype allowed us to estimate genome-wide LD decay. LD was also detected within individual genes on various chromosomes. Principal component and ancestry analyses were used to account for population structure in a genome-wide association study. We further mapped important genes for soluble solid content using a mixed linear model.

Conclusions

Information concerning the SNP resources, population structure, and LD developed in this study will help in identifying agronomically important candidate genes from the genomic regions underlying selection and for mapping quantitative trait loci using a genome-wide association study in sweet watermelon.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-767) contains supplementary material, which is available to authorized users.     

Complete telomere-to-telomere de novo assembly of the Plasmodium falciparum genome through long-read (>11?kb), single molecule,real-time sequencing

The application of next-generation sequencing to estimate genetic diversity of Plasmodium falciparum, the most lethal malaria parasite, has proved challenging due to the skewed AT-richness [∼80.6% (A + T)] of its genome and the lack of technology to assemble highly polymorphic subtelomeric regions that contain clonally variant, multigene virulence families (Ex: var and rifin). To address this, we performed amplification-free, single molecule, real-time sequencing of P. falciparum genomic DNA and generated reads of average length 12 kb, with 50% of the reads between 15.5 and 50 kb in length. Next, using the Hierarchical Genome Assembly Process, we assembled the P. falciparum genome de novo and successfully compiled all 14 nuclear chromosomes telomere-to-telomere. We also accurately resolved centromeres [∼90–99% (A + T)] and subtelomeric regions and identified large insertions and duplications that add extra var and rifin genes to the genome, along with smaller structural variants such as homopolymer tract expansions. Overall, we show that amplification-free, long-read sequencing combined with de novo assembly overcomes major challenges inherent to studying the P. falciparum genome. Indeed, this technology may not only identify the polymorphic and repetitive subtelomeric sequences of parasite populations from endemic areas but may also evaluate structural variation linked to virulence, drug resistance and disease transmission.