Clustering biological sequences with dynamic sequence similarity threshold.

BACKGROUND: Biological sequence clustering is a complicated data clustering problem owing to the high computation costs incurred for pairwise sequence distance calculations through sequence alignments, as well as difficulties in determining parameters for deriving robust clusters. While current approaches are successful in reducing the number of sequence alignments performed, the generated clusters are based on a single sequence identity threshold applied to every cluster. Poor choices of this identity threshold would thus lead to low quality clusters. There is however little support provided to users in selecting thresholds that are well matched with the input sequences. RESULTS: We present a novel sequence clustering approach called ALFATClust that exploits rapid pairwise alignment-free sequence distance calculations and community detection in graph for clusters generation. Instead of a single threshold applied to every generated cluster, ALFATClust is capable of dynamically determining the cut-off threshold for each individual cluster by considering both cluster separation and intra-cluster sequence similarity. Benchmarking analysis shows that ALFATClust generally outperforms existing approaches by simultaneously maintaining cluster robustness and substantial cluster separation for the benchmark datasets. The software also provides an evaluation report for verifying the quality of the non-singleton clusters obtained. CONCLUSIONS: ALFATClust is able to generate sequence clusters having high intra-cluster sequence similarity and substantial separation between clusters without having users to decide precise similarity cut-off thresholds.

Optimized analysis parameters of variant calling for whole genome-based phylogeny of Mycobacteroides abscessus.

Whole-genome sequence (WGS) analysis provides the best resolution for reconstructing bacterial phylogeny. However, the resulting tree could vary according to parameters used in the WGS pipeline, making it difficult to compare results across multiple studies. This study compares effects on phylogenies when applying different parameter stringencies. We used as the study model to optimize parameters strains of Mycobacteroides abscessus serially isolated at various intervals, isolates known to represent persistent infection (PI) cases or re-infection (RI) cases and isolates from different subspecies. Un-optimized parameters with low stringency provided an excessive number of SNPs (823) compared to the optimized setting (3 SNPs) between paired strains isolated 1 day apart from PI cases, discordant tree topology and misclassification of subspecies and of instances of RI. We demonstrated that using high-quality variants provides more accuracy for recognizing serial isolates of the same clone versus different clones and for phylogenetic analysis of M. abscessus. Our approach might be used as a model for analyses requiring phylogenetic reconstruction of other bacteria.

MIRUReader: MIRU-VNTR typing directly from long sequencing reads.

SUMMARY: Mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) typing is widely used to genotype Mycobacterium tuberculosis complex in epidemiological studies for tracking tuberculosis transmission. Recent long-read sequencing technologies from Pacific Biosciences and Oxford Nanopore Technologies can produce reads that are long enough to cover the entire repeat regions in each MIRU-VNTR locus which was previously not possible using the short reads from Illumina high-throughput sequencing technologies. We thus developed MIRUReader for MIRU-VNTR typing directly from long sequence reads. AVAILABILITY AND IMPLEMENTATION: Source code and documentation for MIRUReader program is freely available at https://github.com/phglab/MIRUReader. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Government trust, perceptions of COVID-19 and behaviour change: cohort surveys, Singapore.

OBJECTIVE: To evaluate how public perceptions and trust in government communications affected the adoption of protective behaviour in Singapore during the coronavirus disease 2019 (COVID-19) pandemic. METHODS: We launched our community-based cohort to assess public perceptions of infectious disease outbreaks in mid-2019. After the first case of COVID-19 was reported in Singapore on 23 January, we launched a series of seven COVID-19 surveys to both existing and regularly enrolled new participants every 2 weeks. As well as sociodemographic properties of the participants, we recorded changing responses to judge awareness of the situation, trust in various information sources and perceived risk. We used multivariable logistic regression models to evaluate associations with perceptions of risk and self-reported adopted frequencies of protective behaviour. FINDINGS: Our cohort of 633 participants provided 2857 unique responses during the seven COVID-19 surveys. Most agreed or strongly agreed that information from official government sources (99.1%; 528/533) and Singapore-based news agencies (97.9%; 522/533) was trustworthy. Trust in government communication was significantly associated with higher perceived threat (odds ratio, OR: 2.2; 95% confidence interval, CI: 1.6-3.0), but inversely associated with perceived risk of infection (OR: 0.6; 95% CI: 0.4-0.8) or risk of death if infected (OR: 0.6; 95% CI: 0.4-0.9). Trust in government communication was also associated with a greater likelihood of adopting protective behaviour. CONCLUSION: Our findings show that trust is a vital commodity when managing an evolving outbreak. Our repeated surveys provided real-time feedback, allowing an improved understanding of the interplay between perceptions, trust and behaviour.

Paternal gene pool of Malays in Southeast Asia and its applications for the early expansion of Austronesians.

OBJECTIVES: The origin and differentiation of Austronesian populations and their languages have long fascinated linguists, archeologists, and geneticists. However, the founding process of Austronesians and when they separated from their close relatives, such as the Daic and Austro-Asiatic populations in the mainland of Asia, remain unclear. In this study, we explored the paternal origin of Malays in Southeast Asia and the early differentiation of Austronesians. MATERIALS AND METHODS: We generated whole Y-chromosome sequences of 50 Malays and co-analyzed 200 sequences from other Austronesians and related populations. We generated a revised phylogenetic tree with time estimation. RESULTS: We identified six founding paternal lineages among the studied Malays samples. These founding lineages showed a surprisingly coincident expansion age at 5000 to 6000 years ago. We also found numerous mostly close related samples of the founding lineages of Malays among populations from Mainland of Asia. CONCLUSION: Our analyses provided a refined phylogenetic resolution for the dominant paternal lineages of Austronesians found by previous studies. We suggested that the co-expansion of numerous founding paternal lineages corresponds to the initial differentiation of the most recent common ancestor of modern Austronesians. The splitting time and divergence pattern in perspective of paternal Y-chromosome evidence are highly consistent with the previous theories of ethnologists, linguists, and archeologists.

ARGDIT: a validation and integration toolkit for Antimicrobial Resistance Gene Databases.

MOTIVATION: Antimicrobial resistance is currently one of the main challenges in public health due to the excessive use of antimicrobials in medical treatments and agriculture. The advancements in high-throughput next-generation sequencing and development of bioinformatics tools allow simultaneous detection and identification of antimicrobial resistance genes (ARGs) from clinical, food and environment samples, to monitor the prevalence and track the dissemination of these ARGs. Such analyses are however reliant on a comprehensive database of ARGs with accurate sequence content and annotation. Most of the current ARG databases are therefore manually curated, but this is a time-consuming process and the resulting curation errors could be hard to detect. Several secondary ARG databases consolidate contents from different source ARG databases, and hence modifications in the primary databases might not be propagated and updated promptly in the secondary ARG databases. RESULTS: To address these problems, a validation and integration toolkit called ARGDIT was developed to validate ARG database fidelity, and merge multiple primary ARG databases into a single consolidated secondary ARG database with optional automated sequence re-annotation. Experimental results demonstrated the effectiveness of this toolkit in identifying errors such as sequence annotation typos in current ARG databases and generating an integrated non-redundant ARG database with structured annotation. A toolkit-oriented workflow is also proposed to minimize the efforts in validating, curating and merging multiple ARG protein or coding sequence databases. Database developers therefore benefit from faster update cycles and lower costs for database maintenance, while ARG pipeline users can easily evaluate the reference ARG database quality. AVAILABILITY AND IMPLEMENTATION: ARGDIT is available at https://github.com/phglab/ARGDIT. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Analysis of five deep-sequenced trio-genomes of the Peninsular Malaysia Orang Asli and North Borneo populations.

BACKGROUND: Recent advances in genomic technologies have facilitated genome-wide investigation of human genetic variations. However, most efforts have focused on the major populations, yet trio genomes of indigenous populations from Southeast Asia have been under-investigated. RESULTS: We analyzed the whole-genome deep sequencing data (~ 30×) of five native trios from Peninsular Malaysia and North Borneo, and characterized the genomic variants, including single nucleotide variants (SNVs), small insertions and deletions (indels) and copy number variants (CNVs). We discovered approximately 6.9 million SNVs, 1.2 million indels, and 9000 CNVs in the 15 samples, of which 2.7% SNVs, 2.3% indels and 22% CNVs were novel, implying the insufficient coverage of population diversity in existing databases. We identified a higher proportion of novel variants in the Orang Asli (OA) samples, i.e., the indigenous people from Peninsular Malaysia, than that of the North Bornean (NB) samples, likely due to more complex demographic history and long-time isolation of the OA groups. We used the pedigree information to identify de novo variants and estimated the autosomal mutation rates to be 0.81 × 10- 8 - 1.33 × 10- 8, 1.0 × 10- 9 - 2.9 × 10- 9, and ~ 0.001 per site per generation for SNVs, indels, and CNVs, respectively. The trio-genomes also allowed for haplotype phasing with high accuracy, which serves as references to the future genomic studies of OA and NB populations. In addition, high-frequency inherited CNVs specific to OA or NB were identified. One example is a 50-kb duplication in DEFA1B detected only in the Negrito trios, implying plausible effects on host defense against the exposure of diverse microbial in tropical rainforest environment of these hunter-gatherers. The CNVs shared between OA and NB groups were much fewer than those specific to each group. Nevertheless, we identified a 142-kb duplication in AMY1A in all the 15 samples, and this gene is associated with the high-starch diet. Moreover, novel insertions shared with archaic hominids were identified in our samples. CONCLUSION: Our study presents a full catalogue of the genome variants of the native Malaysian populations, which is a complement of the genome diversity in Southeast Asians. It implies specific population history of the native inhabitants, and demonstrated the necessity of more genome sequencing efforts on the multi-ethnic native groups of Malaysia and Southeast Asia.

Genetic relatedness of indigenous ethnic groups in northern Borneo to neighboring populations from Southeast Asia, as inferred from genome-wide SNP data.

The region of northern Borneo is home to the current state of Sabah, Malaysia. It is located closest to the southern Philippine islands and may have served as a viaduct for ancient human migration onto or off of Borneo Island. In this study, five indigenous ethnic groups from Sabah were subjected to genome-wide SNP genotyping. These individuals represent the "North Borneo"-speaking group of the great Austronesian family. They have traditionally resided in the inland region of Sabah. The dataset was merged with public datasets, and the genetic relatedness of these groups to neighboring populations from the islands of Southeast Asia, mainland Southeast Asia and southern China was inferred. Genetic structure analysis revealed that these groups formed a genetic cluster that was independent of the clusters of neighboring populations. Additionally, these groups exhibited near-absolute proportions of a genetic component that is also common among Austronesians from Taiwan and the Philippines. They showed no genetic admixture with Austro-Melanesian populations. Furthermore, phylogenetic analysis showed that they are closely related to non-Austro-Melansian Filipinos as well as to Taiwan natives but are distantly related to populations from mainland Southeast Asia. Relatively lower heterozygosity and higher pairwise genetic differentiation index (FST ) values than those of nearby populations indicate that these groups might have experienced genetic drift in the past, resulting in their differentiation from other Austronesians. Subsequent formal testing suggested that these populations have received no gene flow from neighboring populations. Taken together, these results imply that the indigenous ethnic groups of northern Borneo shared a common ancestor with Taiwan natives and non-Austro-Melanesian Filipinos and then isolated themselves on the inland of Sabah. This isolation presumably led to no admixture with other populations, and these individuals therefore underwent strong genetic differentiation. This report contributes to addressing the paucity of genetic data on representatives from this strategic region of ancient human migration event(s).

Genomic structure of the native inhabitants of Peninsular Malaysia and North Borneo suggests complex human population history in Southeast Asia.

Southeast Asia (SEA) is enriched with a complex history of peopling. Malaysia, which is located at the crossroads of SEA, has been recognized as one of the hubs for early human migration. To unravel the genomic complexity of the native inhabitants of Malaysia, we sequenced 12 samples from 3 indigenous populations from Peninsular Malaysia and 4 native populations from North Borneo to a high coverage of 28-37×. We showed that the Negritos from Peninsular Malaysia shared a common ancestor with the East Asians, but exhibited some level of gene flow from South Asia, while the North Borneo populations exhibited closer genetic affinity towards East Asians than the Malays. The analysis of time of divergence suggested that ancestors of Negrito were the earliest settlers in the Malay Peninsula, whom first separated from the Papuans ~ 50-33 thousand years ago (kya), followed by East Asian (~ 40-15 kya), while the divergence time frame between North Borneo and East Asia populations predates the Austronesian expansion period implies a possible pre-Neolithic colonization. Substantial Neanderthal ancestry was confirmed in our genomes, as was observed in other East Asians. However, no significant difference was observed, in terms of the proportion of Denisovan gene flow into these native inhabitants from Malaysia. Judging from the similar amount of introgression in the Southeast Asians and East Asians, our findings suggest that the Denisovan gene flow may have occurred before the divergence of these populations and that the shared similarities are likely an ancestral component.

Clinical and Molecular Epidemiology of Carbapenem-Resistant Enterobacteriaceae Among Adult Inpatients in Singapore.

BACKGROUND: Since 2010, the incidence of carbapenem-resistant Enterobacteriaceae (CRE) has been increasing in Singapore. We analyzed the clinical and molecular epidemiology of CRE among adult inpatients in Singapore. METHODS: Quarterly incidence of unique subjects (per 100000 patient-days) with positive clinical and surveillance cultures for CRE were estimated based on mandatory data submitted to the National Public Health Laboratory by public hospitals between 2010 and 2015. CRE-positive adult inpatients were prospectively recruited from 6 public sector hospitals between December 2013 and April 2015. Subjects answered a standardized epidemiologic questionnaire and provided samples for this study. Further clinical information was extracted from subjects' electronic medical records. Whole-genome sequencing was performed on study isolates to determine transmission clusters. RESULTS: Incidence of CRE clinical cultures among adult inpatients plateaued from 2013 (range: 7.73 to 10.32 per 100000 patient-days) following an initial increase between 2010 and end-2012. We prospectively recruited 249 subjects. Their median age was 65 years, 108 (43%) were female, and 161 (64.7%) had carbapenemase-producing Enterobacteriaceae (CPE). On multivariate analysis, prior carbapenem exposure (OR: 3.23; 95% CI: 1.67-6.25) and hematological malignancies (OR: 2.85; 95% CI: 1.10-7.41) were associated with non-carbapenemase-producing CRE (NCPE) (n = 88) compared with CPE (n = 161) subjects. Among 430 CRE isolates from the 249 subjects, 307(71.3%) were CPE, of which 154(50.2%) were blaKPC-positive, 97(31.6%) blaNDM-positive, and 42 (13.7%) blaOXA-positive. Klebsiella pneumoniae (n = 180, 41.9%), Escherichia coli (n = 129, 30.0%) and Enterobacter cloacae (n = 62, 14.4%) were the main Enterobacteriaceae species. WGS (n = 206) revealed diverse bacterial strain type (STs). The predominant blaKPC-positive plasmid was pHS102707 (n = 62, 55.4%) and the predominant blaNDM-positive plasmid was pNDM-ECS01 (n = 46, 48.9%). Five transmission clusters involving 13 subjects were detected. CONCLUSIONS: Clinical CRE trend among adult inpatients showed stabilization following a rapid rise since introduction in 2010 potentially due to infection prevention measures and antimicrobial stewardship. More work is needed on understanding CPE transmission dynamics.

Insights into the genetic structure and diversity of 38 South Asian Indians from deep whole-genome sequencing.

South Asia possesses a significant amount of genetic diversity due to considerable intergroup differences in culture and language. There have been numerous reports on the genetic structure of Asian Indians, although these have mostly relied on genotyping microarrays or targeted sequencing of the mitochondria and Y chromosomes. Asian Indians in Singapore are primarily descendants of immigrants from Dravidian-language-speaking states in south India, and 38 individuals from the general population underwent deep whole-genome sequencing with a target coverage of 30X as part of the Singapore Sequencing Indian Project (SSIP). The genetic structure and diversity of these samples were compared against samples from the Singapore Sequencing Malay Project and populations in Phase 1 of the 1,000 Genomes Project (1 KGP). SSIP samples exhibited greater intra-population genetic diversity and possessed higher heterozygous-to-homozygous genotype ratio than other Asian populations. When compared against a panel of well-defined Asian Indians, the genetic makeup of the SSIP samples was closely related to South Indians. However, even though the SSIP samples clustered distinctly from the Europeans in the global population structure analysis with autosomal SNPs, eight samples were assigned to mitochondrial haplogroups that were predominantly present in Europeans and possessed higher European admixture than the remaining samples. An analysis of the relative relatedness between SSIP with two archaic hominins (Denisovan, Neanderthal) identified higher ancient admixture in East Asian populations than in SSIP. The data resource for these samples is publicly available and is expected to serve as a valuable complement to the South Asian samples in Phase 3 of 1 KGP.

In silico region of difference (RD) analysis of Mycobacterium tuberculosis complex from sequence reads using RD-Analyzer.

BACKGROUND: Whole-genome sequencing is increasingly used in clinical diagnosis of tuberculosis and study of Mycobacterium tuberculosis complex (MTC). MTC consists of several genetically homogenous mycobacteria species which can cause tuberculosis in humans and animals. Regions of difference (RDs) are commonly regarded as gold standard genetic markers for MTC classification. RESULTS: We develop RD-Analyzer, a tool that can accurately infer the species and lineage of MTC isolates from sequence reads based on the presence and absence of a set of 31 RDs. Applied on a publicly available diverse set of 377 sequenced MTC isolates from known major species and lineages, RD-Analyzer achieved an accuracy of 98.14 % (370/377) in species prediction and a concordance of 98.47 % (257/261) in Mycobacterium tuberculosis lineage prediction compared to predictions based on single nucleotide polymorphism markers. By comparing respective sequencing read depths on each genomic position between isolates of different sublineages, we were able to identify the known RD markers in different sublineages of Lineage 4 and provide support for six potential delineating markers having high sensitivities and specificities for sublineage prediction. An extended version of RD-Analyzer was thus developed to allow user-defined RDs for lineage prediction. CONCLUSIONS: RD-Analyzer is a useful and accurate tool for species, lineage and sublineage prediction using known RDs of MTC from sequence reads and is extendable to accepting user-defined RDs for analysis. RD-Analyzer is written in Python and is freely available at https://github.com/xiaeryu/RD-Analyzer .

Local transmission and global dissemination of New Delhi Metallo-Beta-Lactamase (NDM): a whole genome analysis.

BACKGROUND: New Delhi metallo-ß-lactamase (bla NDM), a plasmid-borne carbapenemase gene associated with significant mortality and severely limited treatment options, is of global public health concern as it is found in extremely diverse Gram-negative bacterial strains. This study thus aims to genetically characterize local and global spread of bla NDM. METHODS: To investigate local transmission patterns in the context of a single hospital, whole genome sequencing data of the first 11 bla NDM-positive bacteria isolated in a local hospital were analyzed to: (1) identify and compare bla NDM-positive plasmids; and (2) study the phylogenetic relationship of the bacteria chromosomes. The global analysis was conducted by analyzing 2749 complete plasmid sequences (including 39 bla NDM-positive plasmids) in the NCBI database, where: (1) the plasmids were clustered based on their gene composition similarity; (2) phylogenetic study was conducted for each bla NDM-positive plasmid cluster to infer the phylogenetic relationship within each cluster; (3) gene transposition events introducing bla NDM into different plasmid backbones were identified; and (4) clustering pattern was correlated with the plasmids' incompatibility group and geographical distribution. RESULTS: Analysis of the first 11 bla NDM-positive isolates from a single hospital revealed very low bla NDM-positive plasmid diversity. Local transmission was characterized by clonal spread of a predominant plasmid with 2 sporadic instances of plasmid introduction. In contrast to the low diversity locally, global bla NDM spread involved marked plasmid diversity with no predominant bacterial clone. Thirty-nine (1.4 %) out of the 2749 complete plasmid sequences were bla NDM-positive, and could be resolved into 7 clusters, which were associated with plasmid incompatibility group and geographical distribution. The bla NDM gene module was witnessed to mobilize between different plasmid backbones on at least 6 independent occasions. CONCLUSIONS: Our analysis revealed the complex genetic pathways of bla NDM spread, with global dissemination characterized mainly by transposition of the bla NDM gene cassette into varied plasmids. Early local transmission following plasmid introduction is characterized by plasmid conjugation and bacterial spread. Our findings emphasize the importance of plasmid molecular epidemiology in understanding bla NDM spread.

Predicting HLA alleles from high-resolution SNP data in three Southeast Asian populations.

The major histocompatibility complex (MHC) containing the classical human leukocyte antigen (HLA) Class I and Class II genes is among the most polymorphic and diverse regions in the human genome. Despite the clinical importance of identifying the HLA types, very few databases jointly characterize densely genotyped single nucleotide polymorphisms (SNPs) and HLA alleles in the same samples. To date, the HapMap presents the only public resource that provides a SNP reference panel for predicting HLA alleles, constructed with four collections of individuals of north-western European, northern Han Chinese, cosmopolitan Japanese and Yoruba Nigerian ancestry. Owing to complex patterns of linkage disequilibrium in this region, it is unclear whether the HapMap reference panels can be appropriately utilized for other populations. Here, we describe a public resource for the Singapore Genome Variation Project with: (i) dense genotyping across ∼ 9000 SNPs in the MHC; (ii) four-digit HLA typing for eight Class I and Class II loci, in 96 southern Han Chinese, 89 Southeast Asian Malays and 83 Tamil Indians. This resource provides population estimates of the frequencies of HLA alleles at these eight loci in the three population groups, particularly for HLA-DPA1 and HLA-DPB1 that were not assayed in HapMap. Comparing between population-specific reference panels and a cosmopolitan panel created from all four HapMap populations, we demonstrate that more accurate imputation is obtained with population-specific panels than with the cosmopolitan panel, especially for the Malays and Indians but even when imputing between northern and southern Han Chinese. As with SNP imputation, common HLA alleles were imputed with greater accuracy than low-frequency variants.

Detecting and characterizing genomic signatures of positive selection in global populations.

Natural selection is a significant force that shapes the architecture of the human genome and introduces diversity across global populations. The question of whether advantageous mutations have arisen in the human genome as a result of single or multiple mutation events remains unanswered except for the fact that there exist a handful of genes such as those that confer lactase persistence, affect skin pigmentation, or cause sickle cell anemia. We have developed a long-range-haplotype method for identifying genomic signatures of positive selection to complement existing methods, such as the integrated haplotype score (iHS) or cross-population extended haplotype homozygosity (XP-EHH), for locating signals across the entire allele frequency spectrum. Our method also locates the founder haplotypes that carry the advantageous variants and infers their corresponding population frequencies. This presents an opportunity to systematically interrogate the whole human genome whether a selection signal shared across different populations is the consequence of a single mutation process followed subsequently by gene flow between populations or of convergent evolution due to the occurrence of multiple independent mutation events either at the same variant or within the same gene. The application of our method to data from 14 populations across the world revealed that positive-selection events tend to cluster in populations of the same ancestry. Comparing the founder haplotypes for events that are present across different populations revealed that convergent evolution is a rare occurrence and that the majority of shared signals stem from the same evolutionary event.

Deep whole-genome sequencing of 100 southeast Asian Malays.

Whole-genome sequencing across multiple samples in a population provides an unprecedented opportunity for comprehensively characterizing the polymorphic variants in the population. Although the 1000 Genomes Project (1KGP) has offered brief insights into the value of population-level sequencing, the low coverage has compromised the ability to confidently detect rare and low-frequency variants. In addition, the composition of populations in the 1KGP is not complete, despite the fact that the study design has been extended to more than 2,500 samples from more than 20 population groups. The Malays are one of the Austronesian groups predominantly present in Southeast Asia and Oceania, and the Singapore Sequencing Malay Project (SSMP) aims to perform deep whole-genome sequencing of 100 healthy Malays. By sequencing at a minimum of 30× coverage, we have illustrated the higher sensitivity at detecting low-frequency and rare variants and the ability to investigate the presence of hotspots of functional mutations. Compared to the low-pass sequencing in the 1KGP, the deeper coverage allows more functional variants to be identified for each person. A comparison of the fidelity of genotype imputation of Malays indicated that a population-specific reference panel, such as the SSMP, outperforms a cosmopolitan panel with larger number of individuals for common SNPs. For lower-frequency (<5%) markers, a larger number of individuals might have to be whole-genome sequenced so that the accuracy currently afforded by the 1KGP can be achieved. The SSMP data are expected to be the benchmark for evaluating the value of deep population-level sequencing versus low-pass sequencing, especially in populations that are poorly represented in population-genetics studies.

Tracking inter-institutional spread of NDM and identification of a novel NDM-positive plasmid, pSg1-NDM, using next-generation sequencing approaches.

OBJECTIVES: Owing to gene transposition and plasmid conjugation, New Delhi metallo-ß-lactamase (NDM) is typically identified among varied Enterobacteriaceae species and STs. We used WGS to characterize the chromosomal and plasmid molecular epidemiology of NDM transmission involving four institutions in Singapore. METHODS: Thirty-three Enterobacteriaceae isolates (collection years 2010-14) were sequenced using short-read sequencing-by-synthesis and analysed. Long-read single molecule, real-time sequencing (SMRTS) was used to characterize genetically a novel plasmid pSg1-NDM carried on Klebsiella pneumoniae ST147. RESULTS: In 20 (61%) isolates, blaNDM was located on the pNDM-ECS01 plasmid in the background of multiple bacterial STs, including eight K. pneumoniae STs and five Escherichia coli STs. In six (18%) isolates, a novel blaNDM-positive plasmid, pSg1-NDM, was found only in K. pneumoniae ST147. The pSg1-NDM-K. pneumoniae ST147 clone (Sg1-NDM) was fully sequenced using SMRTS. pSg1-NDM, a 90 103 bp IncR plasmid, carried genes responsible for resistance to six classes of antimicrobials. A large portion of pSg1-NDM had no significant homology to any known plasmids in GenBank. pSg1-NDM had no conjugative transfer region. Combined chromosomal-plasmid phylogenetic analysis revealed five clusters of clonal bacterial NDM-positive plasmid transmission, of which two were inter-institution clusters. The largest inter-institution cluster involved six K. pneumoniae ST147-pSg1-NDM isolates. Fifteen patients were involved in transmission clusters, of which four had ward contact, six had hospital contact and five had an unknown transmission link. CONCLUSIONS: A combined sequencing-by-synthesis and SMRTS approach can determine effectively the transmission clusters of blaNDM and genetically characterize novel plasmids. Plasmid molecular epidemiology is important to understanding NDM spread as blaNDM-positive plasmids can conjugate extensively across species and STs.

Comparing methods for performing trans-ethnic meta-analysis of genome-wide association studies.

Genome-wide association studies (GWASs) have discovered thousands of variants that are associated with human health and disease. Whilst early GWASs have primarily focused on genetically homogeneous populations of European, East Asian and South Asian ancestries, the next-generation genome-wide surveys are starting to pool studies from ethnically diverse populations within a single meta-analysis. However, classical epidemiological strategies for meta-analyses that assume fixed- or random-effects may not be the most suitable approaches to combine GWAS findings as these either confer low statistical power or identify mostly loci where the variants carry homogeneous effect sizes that are present in most of the studies. In a trans-ethnic meta-analysis, it is likely that some genetic loci will exhibit heterogeneous effect sizes across the populations. This may be due to differences in study designs, differences arising from the interactions with other genetic variants, or genuine biological differences attributed to environmental, dietary or lifestyle factors that modulate the influence of the genes. Here we compare different strategies for meta-analyzing GWAS across genetically diverse populations, where we intentionally vary the effect sizes present across the different populations. We subsequently applied the methods that yielded the highest statistical power to a trans-ethnic meta-analysis of seven GWAS in type 2 diabetes, and showed that these methods identified bona fide associations that would otherwise have been missed by the classical strategies.

Multiple Genetic Mutations Associated with Polymyxin Resistance in Acinetobacter baumannii.

We studied polymyxin B resistance in 10 pairs of clinical Acinetobacter baumannii isolates, two of which had developed polymyxin B resistance in vivo. All polymyxin B-resistant isolates had lower growth rates than and substitution mutations in the lpx or pmrB gene compared to their parent isolates. There were significant differences in terms of antibiotic susceptibility and genetic determinants of resistance in A. baumannii isolates that had developed polymyxin B resistance in vivo compared to isolates that had developed polymyxin B resistance in vitro.

Differential positive selection of malaria resistance genes in three indigenous populations of Peninsular Malaysia.

The indigenous populations from Peninsular Malaysia, locally known as Orang Asli, continue to adopt an agro-subsistence nomadic lifestyle, residing primarily within natural jungle habitats. Leading a hunter-gatherer lifestyle in a tropical jungle environment, the Orang Asli are routinely exposed to malaria. Here we surveyed the genetic architecture of individuals from four Orang Asli tribes with high-density genotyping across more than 2.5 million polymorphisms. These tribes reside in different geographical locations in Peninsular Malaysia and belong to three main ethno-linguistic groups, where there is minimal interaction between the tribes. We first dissect the genetic diversity and admixture between the tribes and with neighboring urban populations. Later, by implementing five metrics, we investigated the genome-wide signatures for positive natural selection of these Orang Asli, respectively. Finally, we searched for evidence of genomic adaptation to the pressure of malaria infection. We observed that different evolutionary responses might have emerged in the different Orang Asli communities to mitigate malaria infection.