A ranking system for reference libraries of DNA barcodes: application to marine fish species from Portugal

Background

The increasing availability of reference libraries of DNA barcodes (RLDB) offers the opportunity to the screen the level of consistency in DNA barcode data among libraries, in order to detect possible disagreements generated from taxonomic uncertainty or operational shortcomings. We propose a ranking system to attribute a confidence level to species identifications associated with DNA barcode records from a RLDB. Here we apply the proposed ranking system to a newly generated RLDB for marine fish of Portugal.

Methodology/Principal Findings

Specimens (n = 659) representing 102 marine fish species were collected along the continental shelf of Portugal, morphologically identified and archived in a museum collection. Samples were sequenced at the barcode region of the cytochrome oxidase subunit I gene (COI-5P). Resultant DNA barcodes had average intra-specific and inter-specific Kimura-2-parameter distances (0.32% and 8.84%, respectively) within the range usually observed for marine fishes. All specimens were ranked in five different levels (A–E), according to the reliability of the match between their species identification and the respective diagnostic DNA barcodes. Grades A to E were attributed upon submission of individual specimen sequences to BOLD-IDS and inspection of the clustering pattern in the NJ tree generated. Overall, our study resulted in 73.5% of unambiguous species IDs (grade A), 7.8% taxonomically congruent barcode clusters within our dataset, but awaiting external confirmation (grade B), and 18.7% of species identifications with lower levels of reliability (grades C/E).

Conclusion/Significance

We highlight the importance of implementing a system to rank barcode records in RLDB, in order to flag taxa in need of taxonomic revision, or reduce ambiguities of discordant data. With increasing DNA barcode records publicly available, this cross-validation system would provide a metric of relative accuracy of barcodes, while enabling the continuous revision and annotation required in taxonomic work.     

Reading the complex skipper butterfly fauna of one tropical place

Background

An intense, 30-year, ongoing biodiversity inventory of Lepidoptera, together with their food plants and parasitoids, is centered on the rearing of wild-caught caterpillars in the 120,000 terrestrial hectares of dry, rain, and cloud forest of Area de Conservacion Guanacaste (ACG) in northwestern Costa Rica. Since 2003, DNA barcoding of all species has aided their identification and discovery. We summarize the process and results for a large set of the species of two speciose subfamilies of ACG skipper butterflies (Hesperiidae) and emphasize the effectiveness of barcoding these species (which are often difficult and time-consuming to identify).

Methodology/Principal Findings

Adults are DNA barcoded by the Biodiversity Institute of Ontario, Guelph, Canada; and they are identified by correlating the resulting COI barcode information with more traditional information such as food plant, facies, genitalia, microlocation within ACG, caterpillar traits, etc. This process has found about 303 morphologically defined species of eudamine and pyrgine Hesperiidae breeding in ACG (about 25% of the ACG butterfly fauna) and another 44 units indicated by distinct barcodes (n = 9,094), which may be additional species and therefore may represent as much as a 13% increase. All but the members of one complex can be identified by their DNA barcodes.

Conclusions/Significance

Addition of DNA barcoding to the methodology greatly improved the inventory, both through faster (hence cheaper) accurate identification of the species that are distinguishable without barcoding, as well as those that require it, and through the revelation of species “hidden” within what have long been viewed as single species. Barcoding increased the recognition of species-level specialization. It would be no more appropriate to ignore barcode data in a species inventory than it would be to ignore adult genitalia variation or caterpillar ecology.     

Multi-Locus Variable Number of Tandem Repeat Analysis for Rapid and Accurate Typing of Virulent Multidrug Resistant Escherichia coli Clones

One hundred E. coli isolates from Norway (n = 37), Sweden (n = 24), UK (n = 20) and Spain (n = 19), producing CTX-M-type - (n = 84), or SHV-12 (n = 4) extended spectrum β-lactamases, or the plasmid mediated AmpC, CMY-2 (n = 12), were typed using multi-locus sequence typing (MLST) and multi-locus variable number of tandem repeat analysis (MLVA). Isolates clustered into 33 Sequence Types (STs) and 14 Sequence Type Complexes (STCs), and 58 MLVA-Types (MTs) and 25 different MLVA-Type Complexes (MTCs). A strong agreement between the MLST profile and MLVA typing results was observed, in which all ST131-isolates (n = 39) and most of the STC-648 (n = 10), STC-38 (n = 9), STC-10 (n = 9), STC-405 (n = 8) and STC-23 (n = 6) isolates were clustered distinctly into MTC-29, -36, -20, -14, -10 and -39, respectively. MLVA is a rapid and accurate tool for genotyping isolates of globally disseminated virulent multidrug resistant E. coli lineages, including ST131.     

Novel loci for adiponectin levels and their influence on type 2 diabetes and metabolic traits: a multi-ethnic meta-analysis of 45,891 individuals

Circulating levels of adiponectin, a hormone produced predominantly by adipocytes, are highly heritable and are inversely associated with type 2 diabetes mellitus (T2D) and other metabolic traits. We conducted a meta-analysis of genome-wide association studies in 39,883 individuals of European ancestry to identify genes associated with metabolic disease. We identified 8 novel loci associated with adiponectin levels and confirmed 2 previously reported loci (P = 4.5×10⁻⁸–1.2×10⁻⁴³). Using a novel method to combine data across ethnicities (N = 4,232 African Americans, N = 1,776 Asians, and N = 29,347 Europeans), we identified two additional novel loci. Expression analyses of 436 human adipocyte samples revealed that mRNA levels of 18 genes at candidate regions were associated with adiponectin concentrations after accounting for multiple testing (p<3×10⁻⁴). We next developed a multi-SNP genotypic risk score to test the association of adiponectin decreasing risk alleles on metabolic traits and diseases using consortia-level meta-analytic data. This risk score was associated with increased risk of T2D (p = 4.3×10⁻³, n = 22,044), increased triglycerides (p = 2.6×10⁻¹⁴, n = 93,440), increased waist-to-hip ratio (p = 1.8×10⁻⁵, n = 77,167), increased glucose two hours post oral glucose tolerance testing (p = 4.4×10⁻³, n = 15,234), increased fasting insulin (p = 0.015, n = 48,238), but with lower in HDL-cholesterol concentrations (p = 4.5×10⁻¹³, n = 96,748) and decreased BMI (p = 1.4×10⁻⁴, n = 121,335). These findings identify novel genetic determinants of adiponectin levels, which, taken together, influence risk of T2D and markers of insulin resistance.     

Characterization of Salmonella occurring at high prevalence in a population of the land iguana Conolophus subcristatus in Galápagos Islands, Ecuador

The aim of the study was to elucidate the association between the zoonotic pathogen Salmonella and a population of land iguana, Colonophus subcristatus, endemic to Galápagos Islands in Ecuador. We assessed the presence of Salmonella subspecies and serovars and estimated the prevalence of the pathogen in that population. Additionally, we investigated the genetic relatedness among isolates and serovars utilising pulsed field gel electrophoresis (PFGE) on XbaI-digested DNA and determined the antimicrobial susceptibility to a panel of antimicrobials. The study was carried out by sampling cloacal swabs from animals (n = 63) in their natural environment on in the island of Santa Cruz. A high prevalence (62/63, 98.4%) was observed with heterogeneity of Salmonella subspecies and serovars, all known to be associated with reptiles and with reptile-associated salomonellosis in humans. Serotyping revealed 14 different serovars among four Salmonella enterica subspecies: S. enterica subsp. enterica (n = 48), S. enterica subsp. salamae (n = 2), S. enterica subsp. diarizonae (n = 1), and S. enterica subsp. houtenae (n = 7). Four serovars were predominant: S. Poona (n = 18), S. Pomona (n = 10), S. Abaetetuba (n = 8), and S.Newport (n = 5). The S. Poona isolates revealed nine unique XbaI PFGE patterns, with 15 isolates showing a similarity of 70%. Nine S. Pomona isolates had a similarity of 84%. One main cluster with seven (88%) indistinguishable isolates of S. Abaetetuba was observed. All the Salmonella isolates were pan-susceptible to antimicrobials representative of the most relevant therapeutic classes. The high prevalence and absence of clinical signs suggest a natural interaction of the different Salmonella serovars with the host species. The interaction may have been established before any possible exposure of the iguanas and the biocenosis to direct or indirect environmental factors influenced by the use of antimicrobials in agriculture, in human medicine or in veterinary medicine.     

Occurrence of grapevine leafroll-associated virus complex in Napa Valley

Grapevine leafroll disease (GLD) is caused by a complex of several virus species (grapevine leafroll-associated viruses, GLRaV) in the family Closteroviridae. Because of its increasing importance, it is critical to determine which species of GLRaV is predominant in each region where this disease is occurring. A structured sampling design, utilizing a combination of RT-PCR based testing and sequencing methods, was used to survey GLRaVs in Napa Valley (California, USA) vineyards (n = 36). Of the 216 samples tested for GLRaV-1, -2, -3, -4, -5, and -9, 62% (n = 134) were GLRaV positive. Of the positives, 81% (n = 109) were single infections with GLRaV-3, followed by GLRaV-2 (4%, n = 5), while the remaining samples (15%, n = 20) were mixed infections of GLRaV-3 with GLRaV-1, 2, 4, or 9. Additionally, 468 samples were tested for genetic variants of GLRaV-3, and of the 65% (n = 306) of samples positive for GLRaV-3, 22% were infected with multiple GLRaV-3 variants. Phylogenetic analysis utilizing sequence data from the single infection GLRaV-3 samples produced seven well-supported GLRaV-3 variants, of which three represented 71% of all GLRaV-3 positive samples in Napa Valley. Furthermore, two novel variants, which grouped with a divergent isolate from New Zealand (NZ-1), were identified, and these variants comprised 6% of all positive GLRaV-3 samples. Spatial analyses showed that GLRaV-3a, 3b, and 3c were not homogeneously distributed across Napa Valley. Overall, 86% of all blocks (n = 31) were positive for GLRaVs and 90% of positive blocks (n = 28) had two or more GLRaV-3 variants, suggesting complex disease dynamics that might include multiple insect-mediated introduction events.     

The Guinea-Bissau family of Mycobacterium tuberculosis complex revisited

The Guinea-Bissau family of strains is a unique group of the Mycobacterium tuberculosis complex that, although genotypically closely related, phenotypically demonstrates considerable heterogeneity. We have investigated 414 M. tuberculosis complex strains collected in Guinea-Bissau between 1989 and 2008 in order to further characterize the Guinea-Bissau family of strains. To determine the strain lineages present in the study sample, binary outcomes of spoligotyping were compared with spoligotypes existing in the international database SITVIT2. The major circulating M. tuberculosis clades ranked in the following order: AFRI (n = 195, 47.10%), Latin-American-Mediterranean (LAM) (n = 75, 18.12%), ill-defined T clade (n = 53, 12.8%), Haarlem (n = 37, 8.85%), East-African-Indian (EAI) (n = 25, 6.04%), Unknown (n = 12, 2.87%), Beijing (n = 7, 1.68%), X clade (n = 4, 0.96%), Manu (n = 4, 0.97%), CAS (n = 2, 0.48%). Two strains of the LAM clade isolated in 2007 belonged to the Cameroon family (SIT61). All AFRI isolates except one belonged to the Guinea-Bissau family, i.e. they have an AFRI_1 spoligotype pattern, they have a distinct RFLP pattern with low numbers of IS6110 insertions, and they lack the regions of difference RD7, RD8, RD9 and RD10, RD701 and RD702. This profile classifies the Guinea-Bissau family, irrespective of phenotypic biovar, as part of the M. africanum West African 2 lineage, or the AFRI_1 sublineage according to the spoligtyping nomenclature. Guinea-Bissau family strains display a variation of biochemical traits classically used to differentiate M. tuberculosis from M. bovis. Yet, the differential expression of these biochemical traits was not related to any genes so far investigated (narGHJI and pncA). Guinea-Bissau has the highest prevalence of M. africanum recorded in the African continent, and the Guinea-Bissau family shows a high phylogeographical specificity for Western Africa, with Guinea-Bissau being the epicenter. Trends over time however indicate that this family of strains is waning in most parts of Western Africa, including Guinea-Bissau (p = 0.048).     

Comparative genotypic and phenotypic characterisation of methicillin-resistant Staphylococcus aureus ST398 isolated from animals and humans

The high prevalence of methicillin-resistant Staphylococcus aureus (MRSA) ST398 among pigs in certain European countries and North America and its occurrence in other animal species raises a question concerning the molecular mechanisms mediating the success of this lineage. In this study a panel of S. aureus strains belonging to sequence type (ST) 5 (n = 4), ST8 (n = 5), ST15 (n = 5), ST22 (n = 8), clonal complex (CC) 30 (n = 8), CC97 (n = 8), CC130 (n = 4), CC151 (n = 4) and ST398 (n = 18) were screened by DNA microarray and PCR for the carriage of virulence and antimicrobial resistance genes. Isolates belonging to the same sequence type/clonal complex (ST/CC) were found to share similar virulence gene profiles. The ST398 lineage displayed the lowest content of virulence genes, which consisted mainly of genes detected among the majority or all of the analysed lineages. All MRSA ST398 isolates lacked accessory virulence genes that were detected in other ST/CC. In contrast to virulence genotype, the antimicrobial resistance genes profiles varied between isolates belonging to the same ST/CC and profile similarities could be observed for isolates from different lineages. MRSA ST398 isolates in particular displayed significant diversity and high content of antimicrobial resistance genes. This was comparable with certain MRSA belonging to other sequence types particularly the equine MRSA ST8. The apparent lack of significant virulence genes among MRSA ST398 strains, demonstrates that the lineage features a unique genetic background but no ST398-specific virulence markers could be identified.     

Genetic association study identifies HSPB7 as a risk gene for idiopathic dilated cardiomyopathy

Dilated cardiomyopathy (DCM) is a structural heart disease with strong genetic background. Monogenic forms of DCM are observed in families with mutations located mostly in genes encoding structural and sarcomeric proteins. However, strong evidence suggests that genetic factors also affect the susceptibility to idiopathic DCM. To identify risk alleles for non-familial forms of DCM, we carried out a case-control association study, genotyping 664 DCM cases and 1,874 population-based healthy controls from Germany using a 50K human cardiovascular disease bead chip covering more than 2,000 genes pre-selected for cardiovascular relevance. After quality control, 30,920 single nucleotide polymorphisms (SNP) were tested for association with the disease by logistic regression adjusted for gender, and results were genomic-control corrected. The analysis revealed a significant association between a SNP in HSPB7 gene (rs1739843, minor allele frequency 39%) and idiopathic DCM (p = 1.06×10⁻⁶, OR = 0.67 [95% CI 0.57–0.79] for the minor allele T). Three more SNPs showed p < 2.21×10⁻⁵. De novo genotyping of these four SNPs was done in three independent case-control studies of idiopathic DCM. Association between SNP rs1739843 and DCM was significant in all replication samples: Germany (n = 564, n = 981 controls, p = 2.07×10⁻³, OR = 0.79 [95% CI 0.67–0.92]), France 1 (n = 433 cases, n = 395 controls, p = 3.73×10⁻³, OR = 0.74 [95% CI 0.60–0.91]), and France 2 (n = 249 cases, n = 380 controls, p = 2.26×10⁻⁴, OR = 0.63 [95% CI 0.50–0.81]). The combined analysis of all four studies including a total of n = 1,910 cases and n = 3,630 controls showed highly significant evidence for association between rs1739843 and idiopathic DCM (p = 5.28×10⁻¹³, OR = 0.72 [95% CI 0.65–0.78]). None of the other three SNPs showed significant results in the replication stage.This finding of the HSPB7 gene from a genetic search for idiopathic DCM using a large SNP panel underscores the influence of common polymorphisms on DCM susceptibility.     

Analysis of clonal type-specific antibody reactions in Toxoplasma gondii seropositive humans from Germany by peptide-microarray

Background

Different clonal types of Toxoplasma gondii are thought to be associated with distinct clinical manifestations of infections. Serotyping is a novel technique which may allow to determine the clonal type of T. gondii humans are infected with and to extend typing studies to larger populations which include infected but non-diseased individuals.

Methodology

A peptide-microarray test for T. gondii serotyping was established with 54 previously published synthetic peptides, which mimic clonal type-specific epitopes. The test was applied to human sera (n = 174) collected from individuals with an acute T. gondii infection (n = 21), a latent T. gondii infection (n = 53) and from T. gondii-seropositive forest workers (n = 100).

Findings

The majority (n = 124; 71%) of all T. gondii seropositive human sera showed reactions against synthetic peptides with sequences specific for clonal type II (type II peptides). Type I and type III peptides were recognized by 42% (n = 73) or 16% (n = 28) of the human sera, respectively, while type II–III, type I–III or type I–II peptides were recognized by 49% (n = 85), 36% (n = 62) or 14% (n = 25) of the sera, respectively. Highest reaction intensities were observed with synthetic peptides mimicking type II-specific epitopes. A proportion of the sera (n = 22; 13%) showed no reaction with type-specific peptides. Individuals with acute toxoplasmosis reacted with a statistically significantly higher number of peptides as compared to individuals with latent T. gondii infection or seropositive forest workers.

Conclusions

Type II-specific reactions were overrepresented and higher in intensity in the study population, which was in accord with genotyping studies on T. gondii oocysts previously conducted in the same area. There were also individuals with type I- or type III-specific reactions. Well-characterized reference sera and further specific peptide markers are needed to establish and to perform future serotyping approaches with higher resolution.     

The spectrum of central nervous system infections in an adult referral hospital in hanoi, Vietnam

Objectives

To determine prospectively the causative pathogens of central nervous system (CNS) infections in patients admitted to a tertiary referral hospital in Hanoi, Vietnam.

Methods

From May 2007 to December 2008, cerebrospinal fluid (CSF) samples from 352 adults with suspected meningitis or encephalitis underwent routine testing, staining (Gram, Ziehl-Nielsen, India ink), bacterial culture and polymerase chain reaction targeting Neisseria meningitidis, Streptococcus pneumoniae, S. suis, Haemophilus influenzae type b, Herpes simplex virus (HSV), Varicella Zoster virus (VZV), enterovirus, and 16S ribosomal RNA. Blood cultures and clinically indicated radiology were also performed. Patients were classified as having confirmed or suspected bacterial (BM), tuberculous (TBM), cryptococcal (CRM), eosinophilic (EOM) meningitis, aseptic encephalitis/meningitis (AEM), neurocysticercosis and others.

Results

352 (male: 66%) patients were recruited: median age 34 years (range 13–85). 95/352 (27.3%) diagnoses were laboratory confirmed and one by cranial radiology: BM (n = 62), TBM (n = 9), AEM (n = 19), CRM (n = 5), and neurocysticercosis (n = 1, cranial radiology). S. suis predominated as the cause of BM [48/62 (77.4%)]; Listeria monocytogenese (n = 1), S. pasteurianus (n = 1) and N. meningitidis (n = 2) were infrequent. AEM viruses were: HSV (n = 12), VZV (n = 5) and enterovirus (n = 2). 5 patients had EOM. Of 262/352 (74.4%) patients with full clinical data, 209 (79.8%) were hospital referrals and 186 (71%) had been on antimicrobials. 21 (8%) patients died: TBM (15.2%), AEM (10%), and BM (2.8%).

Conclusions

Most infections lacked microbiological confirmation. S. suis was the most common cause of BM in this setting. Improved diagnostics are needed for meningoencephalitic syndromes to inform treatment and prevention strategies.     

Hepatitis C virus in Vietnam: high prevalence of infection in dialysis and multi-transfused patients involving diverse and novel virus variants

Hepatitis C virus (HCV) is a genetically diverse pathogen infecting approximately 2–3% of the world''s population. Herein, we describe results of a large, multicentre serological and molecular epidemiological study cataloguing the prevalence and genetic diversity of HCV in five regions of Vietnam; Ha Noi, Hai Phong, Da Nang, Khanh Hoa and Can Tho. Individuals (n = 8654) with varying risk factors for infection were analysed for the presence of HCV Ab/Ag and, in a subset of positive specimens, for HCV RNA levels (n = 475) and genotype (n = 282). In lower risk individuals, including voluntary blood donors, military recruits and pregnant women, the prevalence of infection was 0.5% (n = 26/5250). Prevalence rates were significantly higher (p<0.001) in intravenous drug users (IDUs; 55.6%, n = 556/1000), dialysis patients (26.6%, n = 153/575) commercial sex workers (CSWs; 8.7%, n = 87/1000), and recipients of multiple blood transfusions (6.0%, n = 32/529). The prevalence of HCV in dialysis patients varied but remained high in all regions (11–43%) and was associated with the receipt of blood transfusions [OR: 2.08 (1.85–2.34), p = 0.001], time from first transfusion [OR: 1.07 (1.01–1.13), p = 0.023], duration of dialysis [OR: 1.31 (1.19–1.43), p<0.001] and male gender [OR: 1.60 (1.06–2.41), p = 0.026]. Phylogenetic analysis revealed high genetic diversity, particularly amongst dialysis and multi-transfused patients, identifying subtypes 1a (33%), 1b (27%), 2a (0.4%), 3a (0.7%), 3b (1.1%), 6a (18.8%), 6e (6.0%), 6h (4.6%), 6l (6.4%) and 2 clusters of novel genotype 6 variants (2.1%). HCV genotype 1 predominated in Vietnam (60%, n = 169/282) but the proportion of infections attributable to genotype 1 varied between regions and risk groups and, in the Southern part of Vietnam, genotype 6 viruses dominated in dialysis and multi-transfused patients (73.9%). This study confirms a high prevalence of HCV infection in Vietnamese IDUs and, notably, reveals high levels of HCV infection associated with dialysis and blood transfusion.     

Population structure of Hispanics in the United States: the multi-ethnic study of atherosclerosis

Using ∼60,000 SNPs selected for minimal linkage disequilibrium, we perform population structure analysis of 1,374 unrelated Hispanic individuals from the Multi-Ethnic Study of Atherosclerosis (MESA), with self-identification corresponding to Central America (n = 93), Cuba (n = 50), the Dominican Republic (n = 203), Mexico (n = 708), Puerto Rico (n = 192), and South America (n = 111). By projection of principal components (PCs) of ancestry to samples from the HapMap phase III and the Human Genome Diversity Panel (HGDP), we show the first two PCs quantify the Caucasian, African, and Native American origins, while the third and fourth PCs bring out an axis that aligns with known South-to-North geographic location of HGDP Native American samples and further separates MESA Mexican versus Central/South American samples along the same axis. Using k-means clustering computed from the first four PCs, we define four subgroups of the MESA Hispanic cohort that show close agreement with self-identification, labeling the clusters as primarily Dominican/Cuban, Mexican, Central/South American, and Puerto Rican. To demonstrate our recommendations for genetic analysis in the MESA Hispanic cohort, we present pooled and stratified association analysis of triglycerides for selected SNPs in the LPL and TRIB1 gene regions, previously reported in GWAS of triglycerides in Caucasians but as yet unconfirmed in Hispanic populations. We report statistically significant evidence for genetic association in both genes, and we further demonstrate the importance of considering population substructure and genetic heterogeneity in genetic association studies performed in the United States Hispanic population.     

Half of the European fruit fly species barcoded (Diptera,Tephritidae); a feasibility test for molecular?identification

A feasibility test of molecular identification of European fruit flies (Diptera: Tephritidae) based on COI barcode sequences has been executed. A dataset containing 555 sequences of 135 ingroup species from three subfamilies and 42 genera and one single outgroup species has been analysed. 73.3% of all included species could be identified based on their COI barcode gene, based on similarity and distances. The low success rate is caused by singletons as well as some problematic groups: several species groups within the genus Terellia and especially the genus Urophora. With slightly more than 100 sequences – almost 20% of the total – this genus alone constitutes the larger part of the failure for molecular identification for this dataset. Deleting the singletons and Urophora results in a success-rate of 87.1% of all queries and 93.23% of the not discarded queries as correctly identified. Urophora is of special interest due to its economic importance as beneficial species for weed control, therefore it is desirable to have alternative markers for molecular identification.We demonstrate that the success of DNA barcoding for identification purposes strongly depends on the contents of the database used to BLAST against. Especially the necessity of including multiple specimens per species of geographically distinct populations and different ecologies for the understanding of the intra- versus interspecific variation is demonstrated. Furthermore thresholds and the distinction between true and false positives and negatives should not only be used to increase the reliability of the success of molecular identification but also to point out problematic groups, which should then be flagged in the reference database suggesting alternative methods for identification.     

DNA Barcode Identification of Podocarpaceae—The Second Largest Conifer Family

We have generated matK, rbcL, and nrITS2 DNA barcodes for 320 specimens representing all 18 extant genera of the conifer family Podocarpaceae. The sample includes 145 of the 198 recognized species. Comparative analyses of sequence quality and species discrimination were conducted on the 159 individuals from which all three markers were recovered (representing 15 genera and 97 species). The vast majority of sequences were of high quality (B ₃₀ = 0.596–0.989). Even the lowest quality sequences exceeded the minimum requirements of the BARCODE data standard. In the few instances that low quality sequences were generated, the responsible mechanism could not be discerned. There were no statistically significant differences in the discriminatory power of markers or marker combinations (p = 0.05). The discriminatory power of the barcode markers individually and in combination is low (56.7% of species at maximum). In some instances, species discrimination failed in spite of ostensibly useful variation being present (genotypes were shared among species), but in many cases there was simply an absence of sequence variation. Barcode gaps (maximum intraspecific p–distance > minimum interspecific p–distance) were observed in 50.5% of species when all three markers were considered simultaneously. The presence of a barcode gap was not predictive of discrimination success (p = 0.02) and there was no statistically significant difference in the frequency of barcode gaps among markers (p = 0.05). In addition, there was no correlation between number of individuals sampled per species and the presence of a barcode gap (p = 0.27).     

Decrease in pneumococcal co-colonization following vaccination with the seven-valent pneumococcal conjugate vaccine

Understanding the epidemiology of pneumococcal co-colonization is important for monitoring vaccine effectiveness and the occurrence of horizontal gene transfer between pneumococcal strains. In this study we aimed to evaluate the impact of the seven-valent pneumococcal conjugate vaccine (PCV7) on pneumococcal co-colonization among Portuguese children. Nasopharyngeal samples from children up to 6 years old yielding a pneumococcal culture were clustered into three groups: pre-vaccine era (n = 173), unvaccinated children of the vaccine era (n = 169), and fully vaccinated children (4 doses; n = 150). Co-colonization, serotype identification, and relative serotype abundance were detected by analysis of DNA of the total bacterial growth of the primary culture plate using the plyNCR-RFLP method and a molecular serotyping microarray-based strategy. The plyNCR-RFLP method detected an overall co-colonization rate of 20.1%. Microarray analysis confirmed the plyNCR-RFLP results. Vaccination status was the only factor found to be significantly associated with co-colonization: co-colonization rates were significantly lower (p = 0.004; Fisher''s exact test) among fully vaccinated children (8.0%) than among children from the pre-PCV7 era (17.3%) or unvaccinated children of the PCV7 era (18.3%). In the PCV7 era there were significantly less non-vaccine type (NVT) co-colonization events than would be expected based on the NVT distribution observed in the pre-PCV7 era (p = 0.024). In conclusion, vaccination with PCV7 resulted in a lower co-colonization rate due to an asymmetric distribution between NVTs found in single and co-colonized samples. We propose that some NVTs prevalent in the PCV7 era are more competitive than others, hampering their co-existence in the same niche. This result may have important implications since a decrease in co-colonization events is expected to translate in decreased opportunities for horizontal gene transfer, hindering pneumococcal evolution events such as acquisition of antibiotic resistance determinants or capsular switch. This might represent a novel potential benefit of conjugate vaccines.     

Monitoring the systemic human memory B cell compartment of melanoma patients for anti-tumor IgG antibodies

Melanoma, a potentially lethal skin cancer, is widely thought to be immunogenic in nature. While there has been much focus on T cell-mediated immune responses, limited knowledge exists on the role of mature B cells. We describe an approach, including a cell-based ELISA, to evaluate mature IgG antibody responses to melanoma from human peripheral blood B cells. We observed a significant increase in antibody responses from melanoma patients (n = 10) to primary and metastatic melanoma cells compared to healthy volunteers (n = 10) (P<0.0001). Interestingly, we detected a significant reduction in antibody responses to melanoma with advancing disease stage in our patient cohort (n = 21) (P<0.0001). Overall, 28% of melanoma patient-derived B cell cultures (n = 1,800) compared to 2% of cultures from healthy controls (n = 600) produced antibodies that recognized melanoma cells. Lastly, a patient-derived melanoma-specific monoclonal antibody was selected for further study. This antibody effectively killed melanoma cells in vitro via antibody-mediated cellular cytotoxicity. These data demonstrate the presence of a mature systemic B cell response in melanoma patients, which is reduced with disease progression, adding to previous reports of tumor-reactive antibodies in patient sera, and suggesting the merit of future work to elucidate the clinical relevance of activating humoral immune responses to cancer.     

Methylation of the CpG sites only on the sense strand of the APC gene is specific for hepatocellular carcinoma

Hypermethylation of the promoter of the tumor suppressor gene, adenomatous polyposis coli (APC), occurs in various malignancies, including hepatocellular carcinoma (HCC). However, reports on the specificity of the methylation of the APC gene for HCC have varied. To gain insight into how these variations occur, bisulfite PCR sequencing was performed to analyze the methylation status of both sense and antisense strands of the APC gene in samples of HCC tissue, matched adjacent non-HCC liver tissue, hepatitis, cirrhosis, and normal liver tissues. DNA derived from fetal liver and 12 nonhepatic normal tissue was also examined. These experiments revealed liver-specific, antisense strand-biased CpG methylation of the APC gene and suggested that, although methylation of the antisense strand of the APC gene exists in normal liver and other non-HCC disease liver tissue, methylation of the sense strand of the APC gene occurs predominantly in HCC. To determine the effect of the DNA strand on the specificity of the methylated APC gene as a biomarker for HCC detection, quantitative methylation-specific PCR assays for sense and antisense strand DNA were developed and performed on DNA isolated from HCC (n = 58), matched adjacent non-HCC (n = 58), cirrhosis (n = 41), and hepatitis (n = 39). Receiver operating characteristic curves were constructed. With the cutoff value set at the limit of detection, the specificity of sense and antisense strand methylation was 84% and 43%, respectively, and sensitivity was 67.2% and 72.4%, respectively. This result demonstrated that the identity of the methylated DNA strand impacted the specificity of APC for HCC detection. Interestingly, methylation of the sense strand of APC occurred in 40% of HCCs from patients with serum AFP levels less than 20 ng/mL, suggesting a potential role for APC as a biomarker to complement AFP in HCC screening.