Comparison of tet(X4)-containing contigs assembled from metagenomic sequencing data with plasmid sequences of isolates from a cohort of healthy subjects

Recently discovered tet(X) gene variants have provided new insights into microbial antibiotic resistance mechanisms and their potential consequences for public health. This study focused on detection, analysis, and characterization of Tet(X4)-positive Enterobacterales from the gut microbiota of a healthy cohort of individuals in Singapore using cultivation-dependent and cultivation-independent approaches. Twelve Tet(X4)-positive Enterobacterales strains that were previously obtained from the cohort were fully genome-sequenced and comparatively analyzed. A metagenomic sequencing (MS) data set of the same samples was mined for contigs that harbored the tet(X4) resistance gene. The sequences of tet(X4)-containing contigs and plasmids sequences were compared. The presence of the resistance genes floR and estT (previously annotated as catD) was detected in the same cassette in 10 and 12 out of the 12 tet(X4)-carrying plasmids, respectively. MS detected tet(X4)-containing contigs in 2 out of the 109 subjects, while cultivation-dependent analysis previously reported a prevalence of 10.1%. The tet(X4)-containing sequences assembled from MS data are relatively short (~14 to 33 kb) but show high similarity to the respective plasmid sequences of the isolates. Our findings show that MS can complement efforts in the surveillance of antibiotic resistance genes for clinical samples, while it has a lower sensitivity than a cultivation-based method when the target organism has a low abundance. Further optimization is required if MS is to be utilized in antibiotic resistance surveillance.IMPORTANCEThe global rise in antibiotic resistance makes it necessary to develop and apply new approaches to detect and monitor the prevalence of antibiotic resistance genes in human populations. In this regard, of particular interest are resistances against last-resort antibiotics, such as tigecycline. In this study, we show that metagenomic sequencing can help to detect high abundance of the tigecycline resistance gene tet(X4) in fecal samples from a cohort of healthy human subjects. However, cultivation-based approaches currently remain the most reliable and cost-effective method for detection of antibiotic-resistant bacteria.

Mutual Exclusion of Methanobrevibacter Species in the Human Gut Microbiota Facilitates Directed Cultivation of a Candidatus Methanobrevibacter Intestini Representative

Methanogenic Archaea (methanogens) are a phylogenetically diverse group of microorganisms and are considered to be the most abundant archaeal representatives in the human gut. However, the gut methanogen diversity of human populations in many global regions remains poorly investigated. Here, we report the abundance and diversity of gut methanogenic Archaea in a multi-ethnic cohort of healthy Singaporeans by using a concerted approach of metagenomic sequencing, 16S rRNA gene amplicon sequencing, and quantitative PCR. Our results indicate a mutual exclusion of Methanobrevibacter species, i.e., the highly prevalent Methanobrevibacter smithii and the less prevalent Candidatus Methanobrevibacter intestini in more than 80% of the samples when using an amplicon sequencing-based approach. Leveraging on this finding, we were able to select a fecal sample to isolate a representative strain, TLL-48-HuF1, for Candidatus Methanobrevibacter intestini. The analyzed physiological parameters of M. smithii DSM 861T and strain TLL-48-HuF1 suggest high similarity of the two species. Comparative genome analysis and the mutual exclusion of the Methanobrevibacter species indicate potentially different niche adaptation strategies in the human host, which may support the designation of Candidatus M. intestini as a novel species. IMPORTANCE Methanogens are important hydrogen consumers in the gut and are associated with differing host health. Here, we determine the prevalence and abundance of archaeal species in the guts of a multi-ethnic cohort of healthy Singapore residents. While Methanobrevibacter smithii is the most prevalent and abundant methanogen in the human gut of local subjects, the recently proposed Candidatus Methanobrevibacter intestini is the abundant methanogen in a minority of individuals that harbor them. The observed potential mutual exclusion of M. smithii and Ca. M. intestini provides further support to the proposal that the two physiologically similar strains may belong to different Methanobrevibacter species.

THE CARBON FOOTPRINT OF BIOINFORMATICS

Bioinformatic research relies on large-scale computational infrastructures which have a nonzero carbon footprint but so far, no study has quantified the environmental costs of bioinformatic tools and commonly run analyses. In this work, we estimate the carbon footprint of bioinformatics (in kilograms of CO₂ equivalent units, kgCO₂e) using the freely available Green Algorithms calculator (www.green-algorithms.org, last accessed 2022). We assessed 1) bioinformatic approaches in genome-wide association studies (GWAS), RNA sequencing, genome assembly, metagenomics, phylogenetics, and molecular simulations, as well as 2) computation strategies, such as parallelization, CPU (central processing unit) versus GPU (graphics processing unit), cloud versus local computing infrastructure, and geography. In particular, we found that biobank-scale GWAS emitted substantial kgCO₂e and simple software upgrades could make it greener, for example, upgrading from BOLT-LMM v1 to v2.3 reduced carbon footprint by 73%. Moreover, switching from the average data center to a more efficient one can reduce carbon footprint by approximately 34%. Memory over-allocation can also be a substantial contributor to an algorithm’s greenhouse gas emissions. The use of faster processors or greater parallelization reduces running time but can lead to greater carbon footprint. Finally, we provide guidance on how researchers can reduce power consumption and minimize kgCO₂e. Overall, this work elucidates the carbon footprint of common analyses in bioinformatics and provides solutions which empower a move toward greener research.

Emergence of tigecycline- and eravacycline-resistant Tet(X4)-producing Enterobacteriaceae in the gut microbiota of healthy Singaporeans

OBJECTIVES

The recently discovered tigecycline-inactivating enzyme Tet(X4) can confer high-level tigecycline resistance on its hosts, which makes it a public health concern. This study focused on isolation and screening of Tet(X4)-positive Enterobacteriaceae from the gut microbiota of a cohort of healthy individuals in Singapore.

METHODS

MinION and Illumina sequencing was performed to obtain the complete genome sequences of Escherichia coli 2EC1-1 and 94EC. Subsequently, 109 human faecal samples were screened retrospectively for eravacycline-resistant Enterobacteriaceae strains, which were further tested for tet(X4) by PCR. The taxonomy of the isolated strains was determined by 16S rRNA gene PCR and Sanger sequencing.

RESULTS

Comparative genomic analysis of E. coli 2EC1-1 and 94EC revealed that both carry tet(X4), which is encoded by IncI1-type plasmids p2EC1-1 and p94EC-2, respectively. Retrospective screening of faecal samples collected from 109 healthy individuals showed that the faecal carriage rate of Tet(X4)-producing Enterobacteriaceae is 10.1% (95% CI = 5.1%-17.3%), suggesting that tet(X4) is widely distributed in the gut microbiota of healthy individuals in Singapore.

CONCLUSIONS

To the best of our knowledge, this is the first report on the prevalence of tet(X4) in the gut microbiota of a healthy human cohort, as well as the first description of this resistance mechanism outside of China. Our findings suggest that surveillance of tet(X4) in community settings is vital to monitor the spread of this resistance mechanism.

Genome-wide identification of cis DNA methylation quantitative trait loci in three Southeast Asian Populations

DNA methylation (DNAm) is an epigenetic modification that acts to regulate gene transcription, is essential for cellular processes and plays an important role in complex traits and disease. Variation in DNAm levels is influenced by both genetic and environmental factors. Several studies have examined the extent to which common genetic variation influences DNAm (i.e. mQTLs), however, an improved understanding of mQTLs across diverse human populations is needed to increase their utility in integrative genomic studies in order to further our understanding of complex trait and disease biology. Here, we systematically examine cis-mQTLs in three Southeast Asian populations in the Singapore Integrative Omics (iOmics) Study, comprised of Chinese (n = 93), Indians (n = 83) and Malays (n = 78). A total of 24 851 cis-mQTL probes were associated with at least one SNP in meta- and ethnicity-specific analyses at a stringent significance level. These cis-mQTL probes show significant differences in local SNP heritability between the ethnicities, enrichment in functionally relevant regions using data from the Roadmap Epigenomics Mapping Consortium and are associated with nearby genes and complex traits due to pleiotropy. Importantly, DNAm prediction performance and the replication of cis-mQTLs both within iOmics and between two independent mQTL studies in European and Bangladeshi individuals is best when the genetic distance between the ethnicities is small, with differences in cis-mQTLs likely due to differences in allele frequency and linkage disequilibrium. This study highlights the importance of, and opportunities from, extending investigation of the genetic control of DNAm to Southeast Asian populations.

Shared Signature of Recent Positive Selection on the TSBP1-BTNL2-HLA-DRA Genes in Five Native Populations from North Borneo

North Borneo (NB) is home to more than 40 native populations. These natives are believed to have undergone local adaptation in response to environmental challenges such as the mosquito-abundant tropical rainforest. We attempted to trace the footprints of natural selection from the genomic data of NB native populations using a panel of ∼2.2 million genome-wide single nucleotide polymorphisms. As a result, an ∼13-kb haplotype in the Major Histocompatibility Complex Class II region encompassing candidate genes TSBP1–BTNL2–HLA-DRA was identified to be undergoing natural selection. This putative signature of positive selection is shared among the five NB populations and is estimated to have arisen ∼5.5 thousand years (∼220 generations) ago, which coincides with the period of Austronesian expansion. Owing to the long history of endemic malaria in NB, the putative signature of positive selection is postulated to be driven by Plasmodium parasite infection. The findings of this study imply that despite high levels of genetic differentiation, the NB populations might have experienced similar local genetic adaptation resulting from stresses of the shared environment.

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.

Genome-wide association meta-analysis highlights light-induced signaling as a driver for refractive error

Refractive errors, including myopia, are the most frequent eye disorders worldwide and an increasingly common cause of blindness. This genome-wide association meta-analysis in 160,420 participants and replication in 95,505 participants increased the number of established independent signals from 37 to 161 and showed high genetic correlation between Europeans and Asians (>0.78). Expression experiments and comprehensive in silico analyses identified retinal cell physiology and light processing as prominent mechanisms, and also identified functional contributions to refractive-error development in all cell types of the neurosensory retina, retinal pigment epithelium, vascular endothelium and extracellular matrix. Newly identified genes implicate novel mechanisms such as rod-and-cone bipolar synaptic neurotransmission, anterior-segment morphology and angiogenesis. Thirty-one loci resided in or near regions transcribing small RNAs, thus suggesting a role for post-transcriptional regulation. Our results support the notion that refractive errors are caused by a light-dependent retina-to-sclera signaling cascade and delineate potential pathobiological molecular drivers.

Dataset on gene expression in the elderly after Mindfulness Awareness Practice or Health Education Program

It has been reported that relaxation techniques can improve physical health and cognitive function. A number of studies involving different types of relaxation practices showed changes in expression of genes. We investigated the gene expression pattern of a cohort of elderly subjects of Asian descent after weekly (for the first three months) and monthly (for the subsequent six months) intervention. Sixty consenting elderly subjects (aged 60-90 years) with mild cognitive impairment were assigned to either the Mindfulness Awareness Practice (MAP) or Health Education Program (HEP) group in a randomized controlled trial to assess the effectiveness of the programs in preventing further cognitive decline and evaluate the influence on neurological, cellular and biochemical factors. Blood samples were collected before the start of intervention and after nine months for gene expression profiling using Affymetrix Human Genome U133 Plus 2.0 arrays. The dataset is publicly available for further analyses.

The fine-scale genetic structure and evolution of the Japanese population

The contemporary Japanese populations largely consist of three genetically distinct groups—Hondo, Ryukyu and Ainu. By principal-component analysis, while the three groups can be clearly separated, the Hondo people, comprising 99% of the Japanese, form one almost indistinguishable cluster. To understand fine-scale genetic structure, we applied powerful haplotype-based statistical methods to genome-wide single nucleotide polymorphism data from 1600 Japanese individuals, sampled from eight distinct regions in Japan. We then combined the Japanese data with 26 other Asian populations data to analyze the shared ancestry and genetic differentiation. We found that the Japanese could be separated into nine genetic clusters in our dataset, showing a marked concordance with geography; and that major components of ancestry profile of Japanese were from the Korean and Han Chinese clusters. We also detected and dated admixture in the Japanese. While genetic differentiation between Ryukyu and Hondo was suggested to be caused in part by positive selection, genetic differentiation among the Hondo clusters appeared to result principally from genetic drift. Notably, in Asians, we found the possibility that positive selection accentuated genetic differentiation among distant populations but attenuated genetic differentiation among close populations. These findings are significant for studies of human evolution and medical genetics.

Establishing multiple omics baselines for three Southeast Asian populations in the Singapore Integrative Omics Study

The Singapore Integrative Omics Study provides valuable insights on establishing population reference measurement in 364 Chinese, Malay, and Indian individuals. These measurements include > 2.5 millions genetic variants, 21,649 transcripts expression, 282 lipid species quantification, and 284 clinical, lifestyle, and dietary variables. This concept paper introduces the depth of the data resource, and investigates the extent of ethnic variation at these omics and non-omics biomarkers. It is evident that there are specific biomarkers in each of these platforms to differentiate between the ethnicities, and intra-population analyses suggest that Chinese and Indians are the most biologically homogeneous and heterogeneous, respectively, of the three groups. Consistent patterns of correlations between lipid species also suggest the possibility of lipid tagging to simplify future lipidomics assays. The Singapore Integrative Omics Study is expected to allow the characterization of intra-omic and inter-omic correlations within and across all three ethnic groups through a systems biology approach.The Singapore Genome Variation projects characterized the genetics of Singapore's Chinese, Malay, and Indian populations. The Singapore Integrative Omics Study introduced here goes further in providing multi-omic measurements in individuals from these populations, including genetic, transcriptome, lipidome, and lifestyle data, and will facilitate the study of common diseases in Asian communities.

Phylogeography of Y-chromosome haplogroup O3a2b2-N6 reveals patrilineal traces of Austronesian populations on the eastern coastal regions of Asia

Austronesian diffusion is considered one of the greatest dispersals in human history; it led to the peopling of an extremely vast region, ranging from Madagascar in the Indian Ocean to Easter Island in Remote Oceania. The Y-chromosome haplogroup O3a2b*-P164(xM134), a predominant paternal lineage of Austronesian populations, is found at high frequencies in Polynesian populations. However, the internal phylogeny of this haplogroup remains poorly investigated. In this study, we analyzed -seventeen Y-chromosome sequences of haplogroup O3a2b*-P164(xM134) and generated a revised phylogenetic tree of this lineage based on 310 non-private Y-chromosome polymorphisms. We discovered that all available O3a2b*-P164(xM134) samples belong to the newly defined haplogroup O3a2b2-N6 and samples from Austronesian populations belong to the sublineage O3a2b2a2-F706. Additionally, we genotyped a series of Y-chromosome polymorphisms in a large collection of samples from China. We confirmed that the sublineage O3a2b2a2b-B451 is unique to Austronesian populations. We found that O3a2b2-N6 samples are widely distributed on the eastern coastal regions of Asia, from Korea to Vietnam. Furthermore, we propose- that the O3a2b2a2b-B451 lineage represents a genetic connection between ancestors of Austronesian populations and ancient populations in North China, where foxtail millet was domesticated about 11,000 years ago. The large number of newly defined Y-chromosome polymorphisms and the revised phylogenetic tree of O3a2b2-N6 will be helpful to explore the origin of proto-Austronesians and the early diffusion process of Austronesian populations.

Characterising private and shared signatures of positive selection in 37 Asian populations

The Asian Diversity Project (ADP) assembled 37 cosmopolitan and ethnic minority populations in Asia that have been densely genotyped across over half a million markers to study patterns of genetic diversity and positive natural selection. We performed population structure analyses of the ADP populations and divided these populations into four major groups based on their genographic information. By applying a highly sensitive algorithm haploPS to locate genomic signatures of positive selection, 140 distinct genomic regions exhibiting evidence of positive selection in at least one population were identified. We examined the extent of signal sharing for regions that were selected in multiple populations and observed that populations clustered in a similar fashion to that of how the ancestry clades were phylogenetically defined. In particular, populations predominantly located in South Asia underwent considerably different adaptation as compared with populations from the other geographical regions. Signatures of positive selection present in multiple geographical regions were predicted to be older and have emerged prior to the separation of the populations in the different regions. In contrast, selection signals present in a single population group tended to be of lower frequencies and thus can be attributed to recent evolutionary events.

Mapping the genetic diversity of HLA haplotypes in the Japanese populations

Japan has often been viewed as an Asian country that possesses a genetically homogenous community. The basis for partitioning the country into prefectures has largely been geographical, although cultural and linguistic differences still exist between some of the districts/prefectures, especially between Okinawa and the mainland prefectures. The Major Histocompatibility Complex (MHC) region has consistently emerged as the most polymorphic region in the human genome, harbouring numerous biologically important variants; nevertheless the presence of population-specific long haplotypes hinders the imputation of SNPs and classical HLA alleles. Here, we examined the extent of genetic variation at the MHC between eight Japanese populations sampled from Okinawa, and six other prefectures located in or close to the mainland of Japan, specifically focusing at the haplotypes observed within each population, and what the impact of any variation has on imputation. Our results indicated that Okinawa was genetically farther to the mainland Japanese than were Gujarati Indians from Tamil Indians, while the mainland Japanese from six prefectures were more homogeneous than between northern and southern Han Chinese. The distribution of haplotypes across Japan was similar, although imputation was most accurate for Okinawa and several mainland prefectures when population-specific panels were used as reference.

Evaluating the coverage and potential of imputing the exome microarray with next-generation imputation using the 1000 Genomes Project

Next-generation genotyping microarrays have been designed with insights from large-scale sequencing of exomes and whole genomes. The exome genotyping arrays promise to query the functional regions of the human genome at a fraction of the sequencing cost, thus allowing large number of samples to be genotyped. However, two pertinent questions exist: firstly, how representative is the content of the exome chip for populations not involved in the design of the chip; secondly, can the content of the exome chip be imputed with the reference data from the 1000 Genomes Project (1KGP). By deep whole-genome sequencing two Asian populations that are not part of the 1KGP, comprising 96 Southeast Asian Malays and 36 South Asian Indians for which the same samples have also been genotyped on both the Illumina 2.5 M and exome microarrays, we discovered the exome chip is a poor representation of exonic content in our two populations. However, up to 94.1% of the variants on the exome chip that are polymorphic in our populations can be confidently imputed with existing non-exome-centric microarrays using the 1KGP panel. The coverage further increases if there exists population-specific reference data from whole-genome sequencing. There is thus limited gain in using the exome chip for populations not involved in the microarray design. Instead, for the same cost of genotyping 2,000 samples on the exome chip, performing whole-genome sequencing of at least 35 samples in that population to complement the 1KGP may yield a higher coverage of the exonic content from imputation instead.

Risk for ACPA-positive rheumatoid arthritis is driven by shared HLA amino acid polymorphisms in Asian and European populations

Previous studies have emphasized ethnically heterogeneous human leukocyte antigen (HLA) classical allele associations to rheumatoid arthritis (RA) risk. We fine-mapped RA risk alleles within the major histocompatibility complex (MHC) in 2782 seropositive RA cases and 4315 controls of Asian descent. We applied imputation to determine genotypes for eight class I and II HLA genes to Asian populations for the first time using a newly constructed pan-Asian reference panel. First, we empirically measured high imputation accuracy in Asian samples. Then we observed the most significant association in HLA-DRβ1 at amino acid position 13, located outside the classical shared epitope (Pomnibus = 6.9 × 10(-135)). The individual residues at position 13 have relative effects that are consistent with published effects in European populations (His > Phe > Arg > Tyr ≅ Gly > Ser)--but the observed effects in Asians are generally smaller. Applying stepwise conditional analysis, we identified additional independent associations at positions 57 (conditional Pomnibus = 2.2 × 10(-33)) and 74 (conditional Pomnibus = 1.1 × 10(-8)). Outside of HLA-DRβ1, we observed independent effects for amino acid polymorphisms within HLA-B (Asp9, conditional P = 3.8 × 10(-6)) and HLA-DPβ1 (Phe9, conditional P = 3.0 × 10(-5)) concordant with European populations. Our trans-ethnic HLA fine-mapping study reveals that (i) a common set of amino acid residues confer shared effects in European and Asian populations and (ii) these same effects can explain ethnically heterogeneous classical allelic associations (e.g. HLA-DRB1*09:01) due to allele frequency differences between populations. Our study illustrates the value of high-resolution imputation for fine-mapping causal variants in the MHC.

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.

Indians of South Asia has long been a population of interest to a wide audience, due to its unique diversity. We have deep-sequenced 38 individuals of Indian descent residing in Singapore (SSIP) in an effort to illustrate their diversity from a whole-genome standpoint. Indeed, among Asians in our population panel, SSIP was most diverse, followed by the Malays in Singapore (SSMP). Their diversity is further observed in the population's chromosome Y haplogroup and mitochondria haplogroup profiles; individuals with European-dominant haplogroups had greater proportion of European admixture. Among variants (single nucleotide polymorphism and small insertions/deletions) discovered in SSIP, 21.69% were novel with respect to previous sequencing projects. In addition, some 14 loss-of-function variants (LOFs) were associated to cancer, Type II diabetes, and cholesterol levels. Finally, D statistic test with ancient hominids concurred that there was gene flow to East Asians compared to South Asians.

Evaluating the possibility of detecting evidence of positive selection across Asia with sparse genotype data from the HUGO Pan-Asian SNP Consortium

Background

The HUGO Pan-Asian SNP Consortium (PASNP) has generated a genetic resource of almost 55,000 autosomal single nucleotide polymorphisms (SNPs) across more than 1,800 individuals from 73 urban and indigenous populations in Asia. This has offered valuable insights into the correlation between the genetic ancestry of these populations with major linguistic systems and geography. Here, we attempt to understand whether adaptation to local climate, diet and environment partly explains the genetic variation present in these populations by investigating the genomic signatures of positive selection.

Results

To evaluate the impact to the selection analyses due to the considerably lower SNP density as compared to other population genetics resources such as the International HapMap Project (HapMap) or the Singapore Genome Variation Project, we evaluated the extent of haplotype phasing switch errors and the consistency of selection signals from three haplotype-based approaches (iHS, XP-EHH, haploPS) when the HapMap data is thinned to a similar density as PASNP. We subsequently applied haploPS to detect and characterize positive selection in the PASNP populations, identifying 59 genomics regions that were selected in at least one PASNP populations. A cluster analysis on the basis of these 59 signals showed that indigenous populations such as the Negrito from Malaysia and Philippines, the China Hmong, and the Taiwan Ami and Atayal shared more of these signals. We also reported evidence of a positive selection signal encompassing the beta globin gene in the Taiwan Ami and Atayal that was distinct from the signal in the HapMap Africans, suggesting the possibility of convergent evolution at this locus due to malarial selection.

Conclusions

We established that the lower SNP content of the PASNP data conferred weaker ability to detect signatures of positive selection, but the availability of the new approach haploPS retained modest power. Out of all the populations in PASNP, we identified only 59 signals, suggesting a strong need for high-density population-level genotyping data or sequencing data in order to achieve a comprehensive survey of positive selection in Asian populations.

Electronic supplementary material

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