Genomic Diversity and Evolution of Mycobacterium ulcerans Revealed by Next-Generation Sequencing

Mycobacterium ulcerans is the causative agent of Buruli ulcer, the third most common mycobacterial disease after tuberculosis and leprosy. It is an emerging infectious disease that afflicts mainly children and youths in West Africa. Little is known about the evolution and transmission mode of M. ulcerans, partially due to the lack of known genetic polymorphisms among isolates, limiting the application of genetic epidemiology. To systematically profile single nucleotide polymorphisms (SNPs), we sequenced the genomes of three M. ulcerans strains using 454 and Solexa technologies. Comparison with the reference genome of the Ghanaian classical lineage isolate Agy99 revealed 26,564 SNPs in a Japanese strain representing the ancestral lineage. Only 173 SNPs were found when comparing Agy99 with two other Ghanaian isolates, which belong to the two other types previously distinguished in Ghana by variable number tandem repeat typing. We further analyzed a collection of Ghanaian strains using the SNPs discovered. With 68 SNP loci, we were able to differentiate 54 strains into 13 distinct SNP haplotypes. The average SNP nucleotide diversity was low (average 0.06–0.09 across 68 SNP loci), and 96% of the SNP locus pairs were in complete linkage disequilibrium. We estimated that the divergence of the M. ulcerans Ghanaian clade from the Japanese strain occurred 394 to 529 thousand years ago. The Ghanaian subtypes diverged about 1000 to 3000 years ago, or even much more recently, because we found evidence that they evolved significantly faster than average. Our results offer significant insight into the evolution of M. ulcerans and provide a comprehensive report on genetic diversity within a highly clonal M. ulcerans population from a Buruli ulcer endemic region, which can facilitate further epidemiological studies of this pathogen through the development of high-resolution tools.     

Comparative Genome Structure,Secondary Metabolite,and Effector Coding Capacity across Cochliobolus Pathogens

The genomes of five Cochliobolus heterostrophus strains, two Cochliobolus sativus strains, three additional Cochliobolus species (Cochliobolus victoriae, Cochliobolus carbonum, Cochliobolus miyabeanus), and closely related Setosphaeria turcica were sequenced at the Joint Genome Institute (JGI). The datasets were used to identify SNPs between strains and species, unique genomic regions, core secondary metabolism genes, and small secreted protein (SSP) candidate effector encoding genes with a view towards pinpointing structural elements and gene content associated with specificity of these closely related fungi to different cereal hosts. Whole-genome alignment shows that three to five percent of each genome differs between strains of the same species, while a quarter of each genome differs between species. On average, SNP counts among field isolates of the same C. heterostrophus species are more than 25× higher than those between inbred lines and 50× lower than SNPs between Cochliobolus species. The suites of nonribosomal peptide synthetase (NRPS), polyketide synthase (PKS), and SSP–encoding genes are astoundingly diverse among species but remarkably conserved among isolates of the same species, whether inbred or field strains, except for defining examples that map to unique genomic regions. Functional analysis of several strain-unique PKSs and NRPSs reveal a strong correlation with a role in virulence.     

Development of a SNP genotyping panel for genetic monitoring of the laboratory mouse

We have developed a genotyping system for detecting genetic contamination in the laboratory mouse based on assaying single-nucleotide polymorphism (SNP) markers positioned on all autosomes and the X chromosome. This system provides a fast, reliable, and cost-effective way for genetic monitoring, while maintaining a very high degree of confidence. We describe the allelic distribution of 235 SNPs in 48 mouse strains, thereby creating a database of polymorphisms useful for genotyping purposes. The SNP markers used in this study were chosen from publicly available SNP databases. Four genotyping methods were evaluated, and dynamic two-tube allele-specific PCR assays were developed for each marker and tested on a set of 48 inbred mouse strains. The minimal number of assays sufficient to distinguish groups consisting of different numbers of mouse strains was estimated, and a panel of 28 SNPs sufficient to distinguish virtually all of the inbred strains tested was selected. Amplifluor SNP detection assays were developed for these markers and tested on an extended list of 96 strains. This panel was used as a genetic quality control approach to monitor the genotypes of nearly 300 inbred, wild-derived, congenic, consomic, and recombinant inbred strains maintained at The Jackson Laboratory. We have concluded that this marker panel is sufficient for genetic contamination monitoring in colonies containing a large number of genetically diverse mouse strains and that reduced versions of the panel could be implemented in facilities housing a lower number of strains.     

Single-Nucleotide Polymorphism Phylotyping of Escherichia coli

We describe a rapid and easily automated phylogenetic grouping technique based on analysis of bacterial genome single-nucleotide polymorphisms (SNPs). We selected 13 SNPs derived from a complete sequence analysis of 11 essential genes previously used for multilocus sequence typing (MLST) of 30 Escherichia coli strains representing the genetic diversity of the species. The 13 SNPs were localized in five genes, trpA, trpB, putP, icdA, and polB, and were selected to allow recovery of the main phylogenetic groups (groups A, B1, E, D, and B2) and subgroups of the species. In the first step, we validated the SNP approach in silico by extracting SNP data from the complete sequences of the five genes for a panel of 65 pathogenic strains belonging to different E. coli pathovars, which were previously analyzed by MLST. In the second step, we determined these SNPs by dideoxy single-base extension of unlabeled oligonucleotide primers for a collection of 183 commensal and extraintestinal clinical E. coli isolates and compared the SNP phylotyping method to previous well-established typing methods. This SNP phylotyping method proved to be consistent with the other methods for assigning phylogenetic groups to the different E. coli strains. In contrast to the other typing methods, such as multilocus enzyme electrophoresis, ribotyping, or PCR phylotyping using the presence/absence of three genomic DNA fragments, the SNP typing method described here is derived from a solid phylogenetic analysis, and the results obtained by this method are more meaningful. Our results indicate that similar approaches may be used for a wide variety of bacterial species.     

Construction of a single nucleotide polymorphism linkage map for the silkworm, Bombyx mori, based on bacterial artificial chromosome end sequences

We have developed a linkage map for the silkworm Bombyx mori based on single nucleotide polymorphisms (SNPs) between strains p50T and C108T initially found on regions corresponding to the end sequences of bacterial artificial chromosome (BAC) clones. Using 190 segregants from a backcross of a p50T female x an F1 (p50T x C108T) male, we analyzed segregation patterns of 534 SNPs between p50T and C108T, detected among 3840 PCR amplicons, each associated with a p50T BAC end sequence. This enabled us to construct a linkage map composed of 534 SNP markers spanning 1305 cM in total length distributed over the expected 28 linkage groups. Of the 534 BACs whose ends harbored the SNPs used to construct the linkage map, 89 were associated with 107 different ESTs. Since each of the SNP markers is directly linked to a specific genomic BAC clone and to whole-genome sequence data, and some of them are also linked to EST data, the SNP linkage map will be a powerful tool for investigating silkworm genome properties, mutation mapping, and map-based cloning of genes of industrial and agricultural interest.     

Identification of an Interaction between VWF rs7965413 and Platelet Count as a Novel Risk Marker for Metabolic Syndrome: An Extensive Search of Candidate Polymorphisms in a Case-Control Study

Although many single nucleotide polymorphisms (SNPs) have been identified to be associated with metabolic syndrome (MetS), there was only a slight improvement in the ability to predict future MetS by the simply addition of SNPs to clinical risk markers. To improve the ability to predict future MetS, combinational effects, such as SNP—SNP interaction, SNP—environment interaction, and SNP—clinical parameter (SNP × CP) interaction should be also considered. We performed a case-control study to explore novel SNP × CP interactions as risk markers for MetS based on health check-up data of Japanese male employees. We selected 99 SNPs that were previously reported to be associated with MetS and components of MetS; subsequently, we genotyped these SNPs from 360 cases and 1983 control subjects. First, we performed logistic regression analyses to assess the association of each SNP with MetS. Of these SNPs, five SNPs were significantly associated with MetS (P < 0.05): LRP2 rs2544390, rs1800592 between UCP1 and TBC1D9, APOA5 rs662799, VWF rs7965413, and rs1411766 between MYO16 and IRS2. Furthermore, we performed multiple logistic regression analyses, including an SNP term, a CP term, and an SNP × CP interaction term for each CP and SNP that was significantly associated with MetS. We identified a novel SNP × CP interaction between rs7965413 and platelet count that was significantly associated with MetS [SNP term: odds ratio (OR) = 0.78, P = 0.004; SNP × CP interaction term: OR = 1.33, P = 0.001]. This association of the SNP × CP interaction with MetS remained nominally significant in multiple logistic regression analysis after adjustment for either the number of MetS components or MetS components excluding obesity. Our results reveal new insight into platelet count as a risk marker for MetS.     

Low‐coverage whole‐genome sequencing reveals molecular markers for spawning season and sex identification in Gulf of Maine Atlantic cod (Gadus morhua,Linnaeus 1758)

Atlantic cod (Gadus morhua, Linnaeus 1758) in the western Gulf of Maine are managed as a single stock despite several lines of evidence supporting two spawning groups (spring and winter) that overlap spatially, while exhibiting seasonal spawning isolation. Low‐coverage whole‐genome sequencing was used to evaluate the genomic population structure of Atlantic cod spawning groups in the western Gulf of Maine and Georges Bank using 222 individuals collected over multiple years. Results indicated low total genomic differentiation, while also showing strong differentiation between spring and winter‐spawning groups at specific regions of the genome. Guided regularized random forest and ranked F _ST methods were used to select panels of single nucleotide polymorphisms (SNPs) that could reliably distinguish spring and winter‐spawning Atlantic cod (88.5% assignment rate), as well as males and females (95.0% assignment rate) collected in the western Gulf of Maine. These SNP panels represent a valuable tool for fisheries research and management of Atlantic cod in the western Gulf of Maine that will aid investigations of stock production and support accuracy of future assessments.     

Insights to Genetic Characterization Tools for Epidemiological Tracking of Francisella tularensis in Sweden

Tularaemia, caused by the bacterium Francisella tularensis, is endemic in Sweden and is poorly understood. The aim of this study was to evaluate the effectiveness of three different genetic typing systems to link a genetic type to the source and place of tularemia infection in Sweden. Canonical single nucleotide polymorphisms (canSNPs), MLVA including five variable number of tandem repeat loci and PmeI-PFGE were tested on 127 F. tularensis positive specimens collected from Swedish case-patients. All three typing methods identified two major genetic groups with near-perfect agreement. Higher genetic resolution was obtained with canSNP and MLVA compared to PFGE; F. tularensis samples were first assigned into ten phylogroups based on canSNPs followed by 33 unique MLVA types. Phylogroups were geographically analysed to reveal complex phylogeographic patterns in Sweden. The extensive phylogenetic diversity found within individual counties posed a challenge to linking specific genetic types with specific geographic locations. Despite this, a single phylogroup (B.22), defined by a SNP marker specific to a lone Swedish sequenced strain, did link genetic type with a likely geographic place. This result suggests that SNP markers, highly specific to a particular reference genome, may be found most frequently among samples recovered from the same location where the reference genome originated. This insight compels us to consider whole-genome sequencing (WGS) as the appropriate tool for effectively linking specific genetic type to geography. Comparing the WGS of an unknown sample to WGS databases of archived Swedish strains maximizes the likelihood of revealing those rare geographically informative SNPs.     

Efficient strategies for changing the diapause character of silkworm eggs and for the germline transformation of diapause silkworm strains

Abstract To overcome the disadvantages of current silkworm Bombyx mori transgenic technology, such as costly and time‐consuming to maintain non‐diapause transgenic silkworms, we report here on the development of treatments for the germline transformation of diapause silkworm strains. Our results showed that HCl treatment within 3 h of oviposition was able to prevent the diapause of eggs from Japanese lineage diapause silkworm strains and was also suitable for germline transformation of the same strains. By incubating developing mother eggs from Chinese lineage diapause silkworm strains at 15°C (15°C‐IME), we were able to prevent the diapause of their daughter eggs; a similar strategy (15°C‐IMES) for the germline transformation of the same strains was that the mother eggs were incubated at 15°C, and the daughter eggs were then microinjected according to the conventional microinjection methods used for non‐diapause eggs. By combining temperature and light controls, the improved 15°C‐IMES strategy prevented diapause in daughter eggs, and also enabled the germline transformation of both Japanese and Chinese lineage diapause silkworm strains. Although each of the strategies developed here has advantages and disadvantages, we suggest that the 15°C‐IMES strategy is a good reference for the establishment of germline transformation technologies of other egg diapause insects. These new strategies for the efficient germline transformation of diapause silkworm strains are likely to improve the practical use of silkworm transgenic lines in sericulture and also highlight silkworm functional genomics research and its modeling.     

香菇三个功能基因部分序列的单核苷酸多态性分析

以15个香菇栽培品种为材料,对尿嘧啶核苷酸-胞嘧啶核苷酸激酶基因(UMP-CMP kinase gene,uck1)、分裂原活化蛋白激酶基因(mitogen-activated protein kinase gene,mapk)和外切β-1,3葡聚糖酶基因(exo-β-1,3-glucanase-encoding gene,exg1)进行了部分序列的单核苷酸多态性(single nucleotide polymorphism,SNP)分析。结果表明,测序中出现的双峰,是菌丝双核体细胞中两细胞核之间的差异造成的。在采用uck1、mapk和exg1的3,126bp中,共发现48处多态性位点,发生频率为1/65bp,其中36个属于转换、12个为颠换。从群体发生频率上,38个属于超过10%的常见SNP,10个属于罕见SNP。不同基因的SNP发生频率不同,uck1、mapk和exg1的SNP发生频率分别为1.40%、0.82%和2.41%。外显子区SNP数量高于内含子,3个基因在外显子区域分布28个,内含子分布20个。外显子的28个SNP位点中,11个为错义突变,17个为同义突变。错义突变引起了编码氨基酸的改变。对SNP位点的聚类分析表明,15个栽培品种间存在的多态性位点在1–36之间,15个品种的SNP类型不同。uck1,mapk,exg1的SNP可用于香菇栽培品种的鉴别。     

CNVrd,a Read-Depth Algorithm for Assigning Copy-Number at the FCGR Locus: Population-Specific Tagging of Copy Number Variation at FCGR3B

The extent of contribution from common gene copy number (CN) variants in human disease is currently unresolved. Part of the reason for this is the technical difficulty in directly measuring CN variation (CNV) using molecular methods, and the lack of single nucleotide polymorphisms (SNPs) that can tag complex CNV that has arisen multiple times on different SNP haplotypes. One CNV locus implicated in human disease is FCGR. Here we aimed to use next-generation sequencing (NGS) data from the 1000 Genomes Project to assign CN at FCGR3A and FCGR3B and to comprehensively assess the ability of SNPs to tag specific CN variants. A read-depth algorithm was developed (CNVrd) and validated on a subset of HapMap samples using CN assignments that had previously been determined using molecular and microarray methods. At 7 out of 9 other complex loci there was >90% concordance with microarray data. However, given that some prior knowledge of CN is required, the generalizability of CNVrd is limited and should be applied to other complex CNV loci with caution. Subsequently, CN was assigned et FCGR3B using CNVrd in a total of 952 samples from the 1000 Genomes Project, using three classes and SNPs that correlated with duplication were identified. The best tag SNP was observed in the Mexican-American sample set for duplication at FCGR3B. This SNP (rs117435514, r² = 0.79) also tagged similar duplication in Chinese and Japanese (r² = 0.35–0.60), but not in Caucasian or African. No tag SNP for duplication at FCGR3A or deletion at FCGR3B was identified in any population. We conclude that it is possible to tag CNV at the FCGR locus, but CN and SNPs have to be characterized and correlated on a population-specific basis.     

Finding Markers That Make a Difference: DNA Pooling and SNP-Arrays Identify Population Informative Markers for Genetic Stock Identification

Genetic stock identification (GSI) using molecular markers is an important tool for management of migratory species. Here, we tested a cost-effective alternative to individual genotyping, known as allelotyping, for identification of highly informative SNPs for accurate genetic stock identification. We estimated allele frequencies of 2880 SNPs from DNA pools of 23 Atlantic salmon populations using Illumina SNP-chip. We evaluated the performance of four common strategies (global F _ST, pairwise F _ST, Delta and outlier approach) for selection of the most informative set of SNPs and tested their effectiveness for GSI compared to random sets of SNP and microsatellite markers. For the majority of cases, SNPs selected using the outlier approach performed best followed by pairwise F _ST and Delta methods. Overall, the selection procedure reduced the number of SNPs required for accurate GSI by up to 53% compared with randomly chosen SNPs. However, GSI accuracy was more affected by populations in the ascertainment group rather than the ranking method itself. We demonstrated for the first time the compatibility of different large-scale SNP datasets by compiling the largest population genetic dataset for Atlantic salmon to date. Finally, we showed an excellent performance of our top SNPs on an independent set of populations covering the main European distribution range of Atlantic salmon. Taken together, we demonstrate how combination of DNA pooling and SNP arrays can be applied for conservation and management of salmonids as well as other species.     

Proteomic analysis of nucleopolyhedrovirus infection resistance in the silkworm, Bombyx mori (Lepidoptera: Bombycidae)

Silkworm hemolymph is an important defense tissue to resist bacteria and virus infections. To study the response of silkworm hemolymph in the resistance of Bombyx mori L. nucleopolyhedrovirus (BmNPV), we constructed a near-isogenic silkworm line with BmNPV resistance using highly resistant and highly susceptible parental strains. In this paper, two-dimensional gel electrophoresis (2-DE) and Matrix-Assisted Laser Desorption/Ionization (MALDI)-mass spectrometry were employed to investigate the differences of protein patterns in the hemolymph of the highly resistant, highly susceptible and near-isogenic silkworm strains after BmNPV was administrated to the larvae. A comparison between the proteomes of these three silkworm strains led us to identify two differentially expressed proteins, beta-N-acetylglucosaminidase 2 and aminoacylase. The expression levels of these proteins were higher in the BmNPV resistant strains.     

Haplotypic analysis of the TNF locus by association efficiency and entropy

Background

To understand the causal basis of TNF associations with disease, it is necessary to understand the haplotypic structure of this locus. We genotyped 12 single-nucleotide polymorphisms (SNPs) distributed over 4.3 kilobases in 296 healthy, unrelated Gambian and Malawian adults. We generated 592 high-quality haplotypes by integrating family- and population-based reconstruction methods.

Results

We found 32 different haplotypes, of which 13 were shared between the two populations. Both populations were haplotypically diverse (gene diversity = 0.80, Gambia; 0.85, Malawi) and significantly differentiated (p < 10^-5 by exact test). More than a quarter of marker pairs showed evidence of intragenic recombination (29% Gambia; 27% Malawi). We applied two new methods of analyzing haplotypic data: association efficiency analysis (AEA), which describes the ability of each SNP to detect every other SNP in a case-control scenario; and the entropy maximization method (EMM), which selects the subset of SNPs that most effectively dissects the underlying haplotypic structure. AEA revealed that many SNPs in TNF are poor markers of each other. The EMM showed that 8 of 12 SNPs (Gambia) and 7 of 12 SNPs (Malawi) are required to describe 95% of the haplotypic diversity.

Conclusions

The TNF locus in the Gambian and Malawi sample is haplotypically diverse and has a rich history of intragenic recombination. As a consequence, a large proportion of TNF SNPs must be typed to detect a disease-modifying SNP at this locus. The most informative subset of SNPs to genotype differs between the two populations.     

Molecular Analysis of the Relatedness of Five Domesticated Turkey Strains

Our knowledge of the genetic relatedness among the eight existing domesticated turkey strains is limited. To begin to address this paucity, genetic relatedness among five turkey strains (Blue Slate, Bourbon Red, Narragansett, Royal Palm, and Spanish Black) was investigated using three molecular marker systems: randomly amplified polymorphic DNA (RAPD), microsatellite, and SNPs derived from a sequence tagged site and a cloned RAPD fragment. The RAPD analyses were based on five primers that revealed a total of 14 informative DNA fragments in all five populations. The microsatellite analyses involved two informative alleles from three primer-pairs. A total of nine SNPs were detected, one of which appeared to be strain specific. This SNP formed the basis of a PCR-RFLP genotyping procedure developed to distinguish one of the strains from the other four. Evidence from these analyses including the SNP-based RFLP-PCR suggests that Royal Palm is distinct from the other four strains, though more closely related to Narragansett. These data provide, for the first time, molecular evidence of the potential relationships among noncommercial domesticated turkey strains.     

Effect of a glucose impulse on the CcpA regulon in Staphylococcus aureus

Background

A low genetic diversity in Francisella tularensis has been documented. Current DNA based genotyping methods for typing F. tularensis offer a limited and varying degree of subspecies, clade and strain level discrimination power. Whole genome sequencing is the most accurate and reliable method to identify, type and determine phylogenetic relationships among strains of a species. However, lower cost typing schemes are necessary in order to enable typing of hundreds or even thousands of isolates.

Results

We have generated a high-resolution phylogenetic tree from 40 Francisella isolates, including 13 F. tularensis subspecies holarctica (type B) strains, 26 F. tularensis subsp. tularensis (type A) strains and a single F. novicida strain. The tree was generated from global multi-strain single nucleotide polymorphism (SNP) data collected using a set of six Affymetrix GeneChip^® resequencing arrays with the non-repetitive portion of LVS (type B) as the reference sequence complemented with unique sequences of SCHU S4 (type A). Global SNP based phylogenetic clustering was able to resolve all non-related strains. The phylogenetic tree was used to guide the selection of informative SNPs specific to major nodes in the tree for development of a genotyping assay for identification of F. tularensis subspecies and clades. We designed and validated an assay that uses these SNPs to accurately genotype 39 additional F. tularensis strains as type A (A1, A2, A1a or A1b) or type B (B1 or B2).

Conclusion

Whole-genome SNP based clustering was shown to accurately identify SNPs for differentiation of F. tularensis subspecies and clades, emphasizing the potential power and utility of this methodology for selecting SNPs for typing of F. tularensis to the strain level. Additionally, whole genome sequence based SNP information gained from a representative population of strains may be used to perform evolutionary or phylogenetic comparisons of strains, or selection of unique strains for whole-genome sequencing projects.     

Pseudomonas syringae pv. actinidiae from Recent Outbreaks of Kiwifruit Bacterial Canker Belong to Different Clones That Originated in China

A recently emerged plant disease, bacterial canker of kiwifruit (Actinidia deliciosa and A. chinensis), is caused by Pseudomonas syringae pv. actinidiae (PSA). The disease was first reported in China and Japan in the 1980s. A severe outbreak of PSA began in Italy in 2008 and has spread to other European countries. PSA was found in both New Zealand and Chile in 2010. To study the evolution of the pathogen and analyse the transmission of PSA between countries, genomes of strains from China and Japan (where the genus Actinidia is endemic), Italy, New Zealand and Chile were sequenced. The genomes of PSA strains are very similar. However, all strains from New Zealand share several single nucleotide polymorphisms (SNPs) that distinguish them from all other PSA strains. Similarly, all the PSA strains from the 2008 Italian outbreak form a distinct clonal group and those from Chile form a third group. In addition to the rare SNPs present in the core genomes, there is abundant genetic diversity in a genomic island that is part of the accessory genome. The island from several Chinese strains is almost identical to the island present in the New Zealand strains. The island from a different Chinese strain is identical to the island present in the strains from the recent Italian outbreak. The Chilean strains of PSA carry a third variant of this island. These genomic islands are integrative conjugative elements (ICEs). Sequencing of these ICEs provides evidence of three recent horizontal transmissions of ICE from other strains of Pseudomonas syringae to PSA. The analyses of the core genome SNPs and the ICEs, combined with disease history, all support the hypothesis of an independent Chinese origin for both the Italian and the New Zealand outbreaks and suggest the Chilean strains also originate from China.