Competitive Metagenomic DNA Hybridization Identifies Host-Specific Microbial Genetic Markers in Cow Fecal Samples

Several PCR methods have recently been developed to identify fecal contamination in surface waters. In all cases, researchers have relied on one gene or one microorganism for selection of host-specific markers. Here we describe the application of a genome fragment enrichment (GFE) method to identify host-specific genetic markers from fecal microbial community DNA. As a proof of concept, bovine fecal DNA was challenged against a porcine fecal DNA background to select for bovine-specific DNA sequences. Bioinformatic analyses of 380 bovine enriched metagenomic sequences indicated a preponderance of Bacteroidales-like regions predicted to encode membrane-associated and secreted proteins. Oligonucleotide primers capable of annealing to select Bacteroidales-like bovine GFE sequences exhibited extremely high specificity (>99%) in PCR assays with total fecal DNAs from 279 different animal sources. These primers also demonstrated a broad distribution of corresponding genetic markers (81% positive) among 148 different bovine sources. These data demonstrate that direct metagenomic DNA analysis by the competitive solution hybridization approach described is an efficient method for identifying potentially useful fecal genetic markers and for characterizing differences between environmental microbial communities.     

二代及三代测序技术在遗传学诊断中的应用进展

遗传病的防治是公共卫生领域的重大课题,而明确病因是遗传病防治的重要环节。高通量测序技术（又称二代测序技术）具有高通量、低成本、高准确度的优点,为遗传诊断及咨询提供了直接证据,已成为遗传学检测不可或缺的有力工具;第三代测序也凭借其长读长的独特优势在临床应用中占据一席之地。二代及三代测序技术各有特点,互为补充,临床中针对不同的检测需求有多种类型的测序方案可供选择。基于此,对二代及三代测序技术的原理、分类及其在遗传学诊断中的应用进展做一综述,以期为临床测序方案的选择提供思路和指导。     

Next Generation Sequencing Uncovers Unexpected Bacterial Pathogens in Ticks in Western Europe

Background and Aims

Ticks are highly susceptible to global environmental and socio-economical changes. Several tick-borne pathogens have been reported in new geographical regions while new species, strains or genetic variants of tick-borne microorganisms are continually being detected. However, tick-borne pathogens are still poorly understood, and it is estimated that half of all human tick-borne disease has an unknown origin. Therefore in order to prevent these diseases, more effort is required to identify unknown or unexpected tick-borne pathogens. Ixodes ricinus is the vector for a broad range of bacterial pathogens and the most prevalent tick in Europe. The aim of the present study was to evaluate the capability of Next Generation Sequencing (NGS) to extend the inventory of pathogenic bacteria carried by this species of tick in France.

Methods

RNA and DNA were extracted from 1450 I. ricinus questing nymphs collected by flagging in Alsace, France. RNA was pooled and used for NGS. Following de novo assembly, bacterial contigs were assigned to the closest known taxonomy. DNA was used for real time PCR to confirm taxonomic species assignment of NGS-derived contigs for the doubtful cases, and for determination of prevalence.

Results

We have generated a global in-depth picture of tick-borne bacteria. We identified RNA from the main pathogenic bacterial species known to be transmitted by I. ricinus. In addition we also identified unanticipated bacterial species for which we have estimated the prevalence within those ticks inhabiting the studied areas.

Conclusions

The data obtained from this study has proven that NGS has an enormous potential to detect the unexpected and provides the means to monitor pathogen occurrence.     

Genome-Wide Characterization of Insertion and Deletion Variation in Chicken Using Next Generation Sequencing

Insertion and deletion (INDEL) is one of the main events contributing to genetic and phenotypic diversity, which receives less attention than SNP and large structural variation. To gain a better knowledge of INDEL variation in chicken genome, we applied next generation sequencing on 12 diverse chicken breeds at an average effective depth of 8.6. Over 1.3 million non-redundant short INDELs (1–49 bp) were obtained, the vast majority (92.48%) of which were novel. Follow-up validation assays confirmed that most (88.00%) of the randomly selected INDELs represent true variations. The majority (95.76%) of INDELs were less than 10 bp. Both the detected number and affected bases were larger for deletions than insertions. In total, INDELs covered 3.8 Mbp, corresponding to 0.36% of the chicken genome. The average genomic INDEL density was estimated as 0.49 per kb. INDELs were ubiquitous and distributed in a non-uniform fashion across chromosomes, with lower INDEL density in micro-chromosomes than in others, and some functional regions like exons and UTRs were prone to less INDELs than introns and intergenic regions. Nearly 620,253 INDELs fell in genic regions, 1,765 (0.28%) of which located in exons, spanning 1,358 (7.56%) unique Ensembl genes. Many of them are associated with economically important traits and some are the homologues of human disease-related genes. We demonstrate that sequencing multiple individuals at a medium depth offers a promising way for reliable identification of INDELs. The coding INDELs are valuable candidates for further elucidation of the association between genotypes and phenotypes. The chicken INDELs revealed by our study can be useful for future studies, including development of INDEL markers, construction of high density linkage map, INDEL arrays design, and hopefully, molecular breeding programs in chicken.     

丙型肝炎病毒基因变异及候选保守表位分析

HCV分离株主要分为4个基因型(HCV Ⅰ～Ⅳ),各型间的氨基酸及核苷酸组成同源性均小于80％, 氨基酸变异率分别为C 8%,E1 35%,E2/NS1 53%,NS3 27%,NS4 35%,NS5 39%.不同型别的HCV有不同的地区分布特征.根据HCV表达产物多肽的保守性、亲水性、抗原性及空间构型等特性,已在HCV表达产物中鉴定出一些高度保守的候选B细胞表位及T细胞表位, 其中B细胞表位一般为12～40肽, T细胞表位一般为7～9肽, 这些B/T细胞保守表位的鉴定, 将有助于推动HCV的免疫治疗及疫苗研究的发展.     

Sequencing and Genetic Variation of Multidrug Resistance Plasmids in Klebsiella pneumoniae

Background

The development of multidrug resistance is a major problem in the treatment of pathogenic microorganisms by distinct antimicrobial agents. Characterizing the genetic variation among plasmids from different bacterial species or strains is a key step towards understanding the mechanism of virulence and their evolution.

Results

We applied a deep sequencing approach to 206 clinical strains of Klebsiella pneumoniae collected from 2002 to 2008 to understand the genetic variation of multidrug resistance plasmids, and to reveal the dynamic change of drug resistance over time. First, we sequenced three plasmids (70 Kb, 94 Kb, and 147 Kb) from a clonal strain of K. pneumoniae using Sanger sequencing. Using the Illumina sequencing technology, we obtained more than 17 million of short reads from two pooled plasmid samples. We mapped these short reads to the three reference plasmid sequences, and identified a large number of single nucleotide polymorphisms (SNPs) in these pooled plasmids. Many of these SNPs are present in drug-resistance genes. We also found that a significant fraction of short reads could not be mapped to the reference sequences, indicating a high degree of genetic variation among the collection of K. pneumoniae isolates. Moreover, we identified that plasmid conjugative transfer genes and antibiotic resistance genes are more likely to suffer from positive selection, as indicated by the elevated rates of nonsynonymous substitution.

Conclusion

These data represent the first large-scale study of genetic variation in multidrug resistance plasmids and provide insight into the mechanisms of plasmid diversification and the genetic basis of antibiotic resistance.     

A Genome-Wide Survey of Genetic Variation in Gorillas Using Reduced Representation Sequencing

All non-human great apes are endangered in the wild, and it is therefore important to gain an understanding of their demography and genetic diversity. Whole genome assembly projects have provided an invaluable foundation for understanding genetics in all four genera, but to date genetic studies of multiple individuals within great ape species have largely been confined to mitochondrial DNA and a small number of other loci. Here, we present a genome-wide survey of genetic variation in gorillas using a reduced representation sequencing approach, focusing on the two lowland subspecies. We identify 3,006,670 polymorphic sites in 14 individuals: 12 western lowland gorillas (Gorilla gorilla gorilla) and 2 eastern lowland gorillas (Gorilla beringei graueri). We find that the two species are genetically distinct, based on levels of heterozygosity and patterns of allele sharing. Focusing on the western lowland population, we observe evidence for population substructure, and a deficit of rare genetic variants suggesting a recent episode of population contraction. In western lowland gorillas, there is an elevation of variation towards telomeres and centromeres on the chromosomal scale. On a finer scale, we find substantial variation in genetic diversity, including a marked reduction close to the major histocompatibility locus, perhaps indicative of recent strong selection there. These findings suggest that despite their maintaining an overall level of genetic diversity equal to or greater than that of humans, population decline, perhaps associated with disease, has been a significant factor in recent and long-term pressures on wild gorilla populations.     

Schmallenberg Virus Circulation in Culicoides in Belgium in 2012: Field Validation of a Real Time RT-PCR Approach to Assess Virus Replication and Dissemination in Midges

Indigenous Culicoides biting midges are suggested to be putative vectors for the recently emerged Schmallenberg virus (SBV) based on SBV RNA detection in field-caught midges. Furthermore, SBV replication and dissemination has been evidenced in C. sonorensis under laboratory conditions. After SBV had been detected in Culicoides biting midges from Belgium in August 2011, it spread all over the country by the end of 2011, as evidenced by very high between-herd seroprevalence rates in sheep and cattle. This study investigated if a renewed SBV circulation in midges occurred in 2012 in the context of high seroprevalence in the animal host population and evaluated if a recently proposed realtime RT-PCR approach that is meant to allow assessing the vector competence of Culicoides for SBV and bluetongue virus under laboratory conditions was applicable to field-caught midges. Therefore midges caught with 12 OVI traps in four different regions in Belgium between May and November 2012, were morphologically identified, age graded, pooled and tested for the presence of SBV RNA by realtime RT-PCR. The results demonstrate that although no SBV could be detected in nulliparous midges caught in May 2012, a renewed but short lived circulation of SBV in parous midges belonging to the subgenus Avaritia occured in August 2012 at all four regions. The infection prevalence reached up to 2.86% in the south of Belgium, the region where a lower seroprevalence was found at the end of 2011 than in the rest of the country. Furthermore, a frequency analysis of the Ct values obtained for 31 SBV-S segment positive pools of Avaritia midges showed a clear bimodal distribution with peaks of Ct values between 21–24 and 33–36. This closely resembles the laboratory results obtained for SBV infection of C. sonorensis and implicates indigenous midges belonging to the subgenus Avaritia as competent vectors for SBV.     

Genetic and Maternal Variation for Heat Resistance in Drosophila from the Field

In Drosophila, field heritability estimates have focused on morphological traits and ignored maternal effects. This study considers heritable variation and maternal effects in a physiological trait, heat resistance. Drosophila were collected from the field in Melbourne, Australia. Resistance was determined using knock-down time at 37°. Drosophila melanogaster was more resistant than Drosophila simulans, and males tended to be more resistant than females. Field heritability and maternal effects were examined in D. simulans using the regression of laboratory-reared F(1) and F(2) onto field-collected parents. Males from the field were crossed to a laboratory stock to obtain progeny. The additive genetic component to variation in heat resistance was large and significant, and heritability was estimated to be around 0.5. A large maternal effect was also evident. Comparisons of regression coefficients suggested that the maternal effect was not associated with cytoplasmic factors. There was no correlation between body size (as measured by wing length) and heat resistance. Unlike in the case of morphological traits, the heritability for heat resistance in nature is not less than that measured in the laboratory.     

High-Throughput Sequencing Identifies Novel and Conserved Cucumber (Cucumis sativus L.) microRNAs in Response to Cucumber Green Mottle Mosaic Virus Infection

Seedlings of Cucumis sativus L. (cv. ''Zhongnong 16'') were artificially inoculated with Cucumber green mottle mosaic virus (CGMMV) at the three-true-leaf stage. Leaf and flower samples were collected at different time points post-inoculation (10, 30 and 50 d), and processed by high throughput sequencing analysis to identify candidate miRNA sequences. Bioinformatic analysis using screening criteria, and secondary structure prediction, indicated that 8 novel and 23 known miRNAs (including 15 miRNAs described for the first time in vivo) were produced by cucumber plants in response to CGMMV infection. Moreover, gene expression profiles (p-value <0.01) validated the expression of 3 of the novel miRNAs and 3 of the putative candidate miRNAs and identified a further 82 conserved miRNAs in CGMMV-infected cucumbers. Gene ontology (GO) analysis revealed that the predicted target genes of these 88 miRNAs, which were screened using the psRNATarget and miRanda algorithms, were involved in three functional categories: 2265 in molecular function, 1362 as cellular components and 276 in biological process. The subsequent Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that the predicted target genes were frequently involved in metabolic processes (166 pathways) and genetic information processes (40 pathways) and to a lesser degree the biosynthesis of secondary metabolites (12 pathways). These results could provide useful clues to help elucidate host-pathogen interactions in CGMMV and cucumber, as well as for the screening of resistance genes.     

Identification of Genetic Alterations,as Causative Genetic Defects in Long QT Syndrome,Using Next Generation Sequencing Technology

Background

Long QT Syndrome is an inherited channelopathy leading to sudden cardiac death due to ventricular arrhythmias. Despite that several genes have been associated with the disease, nearly 20% of cases remain without an identified genetic cause. Other genetic alterations such as copy number variations have been recently related to Long QT Syndrome. Our aim was to take advantage of current genetic technologies in a family affected by Long QT Syndrome in order to identify the cause of the disease.

Methods

Complete clinical evaluation was performed in all family members. In the index case, a Next Generation Sequencing custom-built panel, including 55 sudden cardiac death-related genes, was used both for detection of sequence and copy number variants. Next Generation Sequencing variants were confirmed by Sanger method. Copy number variations variants were confirmed by Multiplex Ligation dependent Probe Amplification method and at the mRNA level. Confirmed variants and copy number variations identified in the index case were also analyzed in relatives.

Results

In the index case, Next Generation Sequencing revealed a novel variant in TTN and a large deletion in KCNQ1, involving exons 7 and 8. Both variants were confirmed by alternative techniques. The mother and the brother of the index case were also affected by Long QT Syndrome, and family cosegregation was observed for the KCNQ1 deletion, but not for the TTN variant.

Conclusions

Next Generation Sequencing technology allows a comprehensive genetic analysis of arrhythmogenic diseases. We report a copy number variation identified using Next Generation Sequencing analysis in Long QT Syndrome. Clinical and familiar correlation is crucial to elucidate the role of genetic variants identified to distinguish the pathogenic ones from genetic noise.     

Genetic Variation of Human Papillomavirus Type 16 in Individual Clinical Specimens Revealed by Deep Sequencing

Viral genetic diversity within infected cells or tissues, called viral quasispecies, has been mostly studied for RNA viruses, but has also been described among DNA viruses, including human papillomavirus type 16 (HPV16) present in cervical precancerous lesions. However, the extent of HPV genetic variation in cervical specimens, and its involvement in HPV-induced carcinogenesis, remains unclear. Here, we employ deep sequencing to comprehensively analyze genetic variation in the HPV16 genome isolated from individual clinical specimens. Through overlapping full-circle PCR, approximately 8-kb DNA fragments covering the whole HPV16 genome were amplified from HPV16-positive cervical exfoliated cells collected from patients with either low-grade squamous intraepithelial lesion (LSIL) or invasive cervical cancer (ICC). Deep sequencing of the amplified HPV16 DNA enabled de novo assembly of the full-length HPV16 genome sequence for each of 7 specimens (5 LSIL and 2 ICC samples). Subsequent alignment of read sequences to the assembled HPV16 sequence revealed that 2 LSILs and 1 ICC contained nucleotide variations within E6, E1 and the non-coding region between E5 and L2 with mutation frequencies of 0.60% to 5.42%. In transient replication assays, a novel E1 mutant found in ICC, E1 Q381E, showed reduced ability to support HPV16 origin-dependent replication. In addition, partially deleted E2 genes were detected in 1 LSIL sample in a mixed state with the intact E2 gene. Thus, the methods used in this study provide a fundamental framework for investigating the influence of HPV somatic genetic variation on cervical carcinogenesis.     

Integration of Experiments across Diverse Environments Identifies the Genetic Determinants of Variation in Sorghum bicolor Seed Element Composition

Seedling establishment and seed nutritional quality require the sequestration of sufficient element nutrients. The identification of genes and alleles that modify element content in the grains of cereals, including sorghum (Sorghum bicolor), is fundamental to developing breeding and selection methods aimed at increasing bioavailable element content and improving crop growth. We have developed a high-throughput work flow for the simultaneous measurement of multiple elements in sorghum seeds. We measured seed element levels in the genotyped Sorghum Association Panel, representing all major cultivated sorghum races from diverse geographic and climatic regions, and mapped alleles contributing to seed element variation across three environments by genome-wide association. We observed significant phenotypic and genetic correlation between several elements across multiple years and diverse environments. The power of combining high-precision measurements with genome-wide association was demonstrated by implementing rank transformation and a multilocus mixed model to map alleles controlling 20 element traits, identifying 255 loci affecting the sorghum seed ionome. Sequence similarity to genes characterized in previous studies identified likely causative genes for the accumulation of zinc, manganese, nickel, calcium, and cadmium in sorghum seeds. In addition to strong candidates for these five elements, we provide a list of candidate loci for several other elements. Our approach enabled the identification of single-nucleotide polymorphisms in strong linkage disequilibrium with causative polymorphisms that can be evaluated in targeted selection strategies for plant breeding and improvement.Sorghum (Sorghum bicolor) is a globally cultivated source of food, feed, and fiber. Contrasting needs for elemental nutrient accumulation limit crop yield and quality for sorghum marketed to different sectors. The seed-bearing reproductive organs, or panicles, in sorghum represent up to 30% of the total dry matter yield (Amaducci et al., 2004). Plant-based diets, in which grains compose the major food source, require the accumulation of bioavailable essential elements in the plant seeds. Currently, iron (Fe) and zinc (Zn) deficiencies negatively affect the health of over two billion people worldwide (World Health Organization, 2002). Increased bioavailable elemental nutrient content in the edible portions of sorghum for human and animal nutrition could ameliorate this nutritional crisis (Graham et al., 1999; World Health Organization, 2002). Additional global health benefits could be achieved by increasing magnesium (Mg), selenium (Se), calcium (Ca), and copper (Cu; White and Broadley, 2005) while reducing the concentration of toxic elements, including arsenic (As) and cadmium (Cd; Ma et al., 2008).Seed element accumulation results from interconnected biological processes, including element uptake by the roots, translocation and remobilization within the plant, and ultimately import, deposition, and assimilation/storage in the seeds. Element availability is further affected by the accumulation of metabolites in seeds (Vreugdenhil et al., 2004). High-throughput ionomic analysis, or concurrent measurement of multiple elements, allows for the quantitative and simultaneous measurement of an organism’s elemental composition, providing a snapshot of the functional state of an organism under different experimental conditions (Salt et al., 2008). Most studies of the plant ionome utilize inductively coupled plasma mass spectroscopy (ICP-MS). Briefly, inductively coupled plasma (ICP) functions to ionize the analyte into atoms, which are then detected by mass spectroscopy. Reference standards are used to identify and quantitate each element of interest in the sample. ICP-MS analysis can be accomplished in as little as 1 min per sample, which allows for high-throughput processing of thousands of samples (Salt et al., 2008). Previous studies have demonstrated that several elements, including Fe, manganese (Mn), Zn, cobalt (Co), and Cd, share mechanisms of accumulation (Yi and Guerinot, 1996; Vert et al., 2002; Connolly et al., 2003). Ionomic signatures derived from multiple elements also have been shown to better predict plant physiological status for some elements than the measure of the element’s concentration, including essential nutrients like Fe (Baxter et al., 2008). Holistically examining the ionome provides significant insights into the networks underlying ion homeostasis beyond single-element studies (Baxter and Dilkes, 2012).There are over 45,000 catalogued lines of sorghum at the U.S. Department of Agriculture Germplasm Resource Information Network. This diverse collection of sorghum germplasm contains genetic variation with undiscovered impact on seed element composition (Das et al., 1997). Mapping quantitative trait loci for seed element concentration has been successful in a number of species, including Arabidopsis (Arabidopsis thaliana; Vreugdenhil et al., 2004; Waters and Grusak, 2008; Buescher et al., 2010), rice (Oryza sativa; Norton et al., 2010; Zhang et al., 2014), wheat (Triticum aestivum; Shi et al., 2008; Peleg et al., 2009), and maize (Zea mays; Simić et al., 2012; Baxter et al., 2013, 2014). Genome-wide association (GWA) mapping is well suited for uncovering the genetic basis for complex traits, including seed element accumulation. One of the key strengths of association mapping is that a priori knowledge is not necessary to identify new loci associated with the trait of interest. Furthermore, a GWA mapping population is composed of lines that have undergone numerous recombination events, allowing for a narrower mapping interval. Previous GWA studies in maize (Tian et al., 2011), rice (Huang et al., 2010), and sorghum (Morris et al., 2013) have been successful in identifying the genetic basis for various agronomic traits. Here, we analyzed the seed ionome from a community-generated association panel to identify potential loci underlying seed element accumulation in sorghum.     

Variation in Field Isolates of Measles Virus during an 8-Year Period in Japan

Field isolates of measles virus (MV) during an 8-year period in four areas of Japan, i.e., Osaka, Nagoya, Tokyo and Akita, were classified into three types in regard to the electrophoretic mobility of the hemagglutinin (HA) proteins: S type with small (78 K) HA, M type with intermediate (80 K) HA and L type with large (82K) HA. The type of field isolates was closely related with the geographical location and the year of virus isolation. The S type strain was isolated only in an outbreak from 1983 to 1984, whereas the M and L type strains were isolated between 1983 and 1990. The HA genes of the M and L type strains of MV were found to have a nucleotide substitution which introduces a new potential glycosylation site. In addition, the matrix proteins of all field strains isolated after 1977 showed slower electrophoretic mobility of 42 K than 39 K of the Edmonston and Toyoshima strains. These results indicate that MV strains of different HA types existed concomitantly and that major populations of MV currently circulating in Japan are changing from those prevalent in 1983-1984.     

Novel Virus Discovery and Genome Reconstruction from Field RNA Samples Reveals Highly Divergent Viruses in Dipteran Hosts

We investigated whether small RNA (sRNA) sequenced from field-collected mosquitoes and chironomids (Diptera) can be used as a proxy signature of viral prevalence within a range of species and viral groups, using sRNAs sequenced from wild-caught specimens, to inform total RNA deep sequencing of samples of particular interest. Using this strategy, we sequenced from adult Anopheles maculipennis s.l. mosquitoes the apparently nearly complete genome of one previously undescribed virus related to chronic bee paralysis virus, and, from a pool of Ochlerotatus caspius and Oc. detritus mosquitoes, a nearly complete entomobirnavirus genome. We also reconstructed long sequences (1503-6557 nt) related to at least nine other viruses. Crucially, several of the sequences detected were reconstructed from host organisms highly divergent from those in which related viruses have been previously isolated or discovered. It is clear that viral transmission and maintenance cycles in nature are likely to be significantly more complex and taxonomically diverse than previously expected.     

Whole Genome Sequencing of Field Isolates Provides Robust Characterization of Genetic Diversity in Plasmodium vivax

Background

An estimated 2.85 billion people live at risk of Plasmodium vivax transmission. In endemic countries vivax malaria causes significant morbidity and its mortality is becoming more widely appreciated, drug-resistant strains are increasing in prevalence, and an increasing number of reports indicate that P. vivax is capable of breaking through the Duffy-negative barrier long considered to confer resistance to blood stage infection. Absence of robust in vitro propagation limits our understanding of fundamental aspects of the parasite''s biology, including the determinants of its dormant hypnozoite phase, its virulence and drug susceptibility, and the molecular mechanisms underlying red blood cell invasion.

Methodology/Principal Findings

Here, we report results from whole genome sequencing of five P. vivax isolates obtained from Malagasy and Cambodian patients, and of the monkey-adapted Belem strain. We obtained an average 70–400 X coverage of each genome, resulting in more than 93% of the Sal I reference sequence covered by 20 reads or more. Our study identifies more than 80,000 SNPs distributed throughout the genome which will allow designing association studies and population surveys. Analysis of the genome-wide genetic diversity in P. vivax also reveals considerable allele sharing among isolates from different continents. This observation could be consistent with a high level of gene flow among parasite strains distributed throughout the world.

Conclusions

Our study shows that it is feasible to perform whole genome sequencing of P. vivax field isolates and rigorously characterize the genetic diversity of this parasite. The catalogue of polymorphisms generated here will enable large-scale genotyping studies and contribute to a better understanding of P. vivax traits such as drug resistance or erythrocyte invasion, partially circumventing the lack of laboratory culture that has hampered vivax research for years.     

Whole Genome Sequencing of Field Isolates Reveals Extensive Genetic Diversity in Plasmodium vivax from Colombia

Plasmodium vivax is the most prevalent malarial species in South America and exerts a substantial burden on the populations it affects. The control and eventual elimination of P. vivax are global health priorities. Genomic research contributes to this objective by improving our understanding of the biology of P. vivax and through the development of new genetic markers that can be used to monitor efforts to reduce malaria transmission. Here we analyze whole-genome data from eight field samples from a region in Cordóba, Colombia where malaria is endemic. We find considerable genetic diversity within this population, a result that contrasts with earlier studies suggesting that P. vivax had limited diversity in the Americas. We also identify a selective sweep around a substitution known to confer resistance to sulphadoxine-pyrimethamine (SP). This is the first observation of a selective sweep for SP resistance in this species. These results indicate that P. vivax has been exposed to SP pressure even when the drug is not in use as a first line treatment for patients afflicted by this parasite. We identify multiple non-synonymous substitutions in three other genes known to be involved with drug resistance in Plasmodium species. Finally, we found extensive microsatellite polymorphisms. Using this information we developed 18 polymorphic and easy to score microsatellite loci that can be used in epidemiological investigations in South America.