Reconstructing the Human Genetic History of Mainland Southeast Asia: Insights from Genome-Wide Data from Thailand and Laos

Thailand and Laos, located in the center of Mainland Southeast Asia (MSEA), harbor diverse ethnolinguistic groups encompassing all five language families of MSEA: Tai-Kadai (TK), Austroasiatic (AA), Sino-Tibetan (ST), Hmong-Mien (HM), and Austronesian (AN). Previous genetic studies of Thai/Lao populations have focused almost exclusively on uniparental markers and there is a paucity of genome-wide studies. We therefore generated genome-wide SNP data for 33 ethnolinguistic groups, belonging to the five MSEA language families from Thailand and Laos, and analyzed these together with data from modern Asian populations and SEA ancient samples. Overall, we find genetic structure according to language family, albeit with heterogeneity in the AA-, HM-, and ST-speaking groups, and in the hill tribes, that reflects both population interactions and genetic drift. For the TK speaking groups, we find localized genetic structure that is driven by different levels of interaction with other groups in the same geographic region. Several Thai groups exhibit admixture from South Asia, which we date to ∼600–1000 years ago, corresponding to a time of intensive international trade networks that had a major cultural impact on Thailand. An AN group from Southern Thailand shows both South Asian admixture as well as overall affinities with AA-speaking groups in the region, suggesting an impact of cultural diffusion. Overall, we provide the first detailed insights into the genetic profiles of Thai/Lao ethnolinguistic groups, which should be helpful for reconstructing human genetic history in MSEA and selecting populations for participation in ongoing whole genome sequence and biomedical studies.     

Patrilineal perspective on the Austronesian diffusion in Mainland Southeast Asia

The Cham people are the major Austronesian speakers of Mainland Southeast Asia (MSEA) and the reconstruction of the Cham population history can provide insights into their diffusion. In this study, we analyzed non-recombining region of the Y chromosome markers of 177 unrelated males from four populations in MSEA, including 59 Cham, 76 Kinh, 25 Lao, and 17 Thai individuals. Incorporating published data from mitochondrial DNA (mtDNA), our results indicated that, in general, the Chams are an indigenous Southeast Asian population. The origin of the Cham people involves the genetic admixture of the Austronesian immigrants from Island Southeast Asia (ISEA) with the local populations in MSEA. Discordance between the overall patterns of Y chromosome and mtDNA in the Chams is evidenced by the presence of some Y chromosome lineages that prevail in South Asians. Our results suggest that male-mediated dispersals via the spread of religions and business trade might play an important role in shaping the patrilineal gene pool of the Cham people.     

Gene Flow between the Korean Peninsula and Its Neighboring Countries

SNP markers provide the primary data for population structure analysis. In this study, we employed whole-genome autosomal SNPs as a marker set (54,836 SNP markers) and tested their possible effects on genetic ancestry using 320 subjects covering 24 regional groups including Northern ( = 16) and Southern ( = 3) Asians, Amerindians ( = 1), and four HapMap populations (YRI, CEU, JPT, and CHB). Additionally, we evaluated the effectiveness and robustness of 50K autosomal SNPs with various clustering methods, along with their dependencies on recombination hotspots (RH), linkage disequilibrium (LD), missing calls and regional specific markers. The RH- and LD-free multi-dimensional scaling (MDS) method showed a broad picture of human migration from Africa to North-East Asia on our genome map, supporting results from previous haploid DNA studies. Of the Asian groups, the East Asian group showed greater differentiation than the Northern and Southern Asian groups with respect to Fst statistics. By extension, the analysis of monomorphic markers implied that nine out of ten historical regions in South Korea, and Tokyo in Japan, showed signs of genetic drift caused by the later settlement of East Asia (South Korea, Japan and China), while Gyeongju in South East Korea showed signs of the earliest settlement in East Asia. In the genome map, the gene flow to the Korean Peninsula from its neighboring countries indicated that some genetic signals from Northern populations such as the Siberians and Mongolians still remain in the South East and West regions, while few signals remain from the early Southern lineages.     

Microsatellite markers reveal population structure and low gene flow among collections of Bactrocera cucurbitae (Diptera: Tephritidae) in Asia

The melon fruit fly, Bactrocera cucurbitae (Coquillett) (Diptera: Tephritidae), is widespread agricultural pest, and it is known to have the potential to establish invasive populations in various tropical and subtropical areas. Despite the economic risk associated with a putative stable presence of this fly, the population genetics of this pest have remained relatively unexplored in Asia, the main area for distribution of this pest. The goals for this study were to employ nuclear markers to examine geographic collections for population genetic structure and quantify the extent of gene flow within these Southeast Asian and Chinese populations. To achieve these goals, we used 12 polymorphic microsatellite markers. A low level of genetic diversity was found among collections from China and higher levels were seen in Southeast Asia collections. Three genetically distinct groups, Southeast Asia, southwest China, and southeast China, were recovered by Bayesian model-based clustering methods, the phylogenetic reconstruction and the principal coordinate analysis. The Mantel test clearly shows geographical distance contributed in the genetic structuring of B. cucurbitae's populations. No recent bottlenecks for any of the populations examined. The results of clustering, migration analyses, and Mantel test, strongly suggest that the regional structure observed may be due to geographical factors such as mountains, rivers, and islands. We found a high rate of migration in some sites from the southwest China region (cluster 1) and the southeast China region (cluster 2), suggesting that China-Guangdong-Guangzhou (GZ) may be the center of melon fruit fly in the southeast China region.     

Evolutionary history of Helicobacter pylori sequences reflect past human migrations in Southeast Asia

The human population history in Southeast Asia was shaped by numerous migrations and population expansions. Their reconstruction based on archaeological, linguistic or human genetic data is often hampered by the limited number of informative polymorphisms in classical human genetic markers, such as the hypervariable regions of the mitochondrial DNA. Here, we analyse housekeeping gene sequences of the human stomach bacterium Helicobacter pylori from various countries in Southeast Asia and we provide evidence that H. pylori accompanied at least three ancient human migrations into this area: i) a migration from India introducing hpEurope bacteria into Thailand, Cambodia and Malaysia; ii) a migration of the ancestors of Austro-Asiatic speaking people into Vietnam and Cambodia carrying hspEAsia bacteria; and iii) a migration of the ancestors of the Thai people from Southern China into Thailand carrying H. pylori of population hpAsia2. Moreover, the H. pylori sequences reflect iv) the migrations of Chinese to Thailand and Malaysia within the last 200 years spreading hspEasia strains, and v) migrations of Indians to Malaysia within the last 200 years distributing both hpAsia2 and hpEurope bacteria. The distribution of the bacterial populations seems to strongly influence the incidence of gastric cancer as countries with predominantly hspEAsia isolates exhibit a high incidence of gastric cancer while the incidence is low in countries with a high proportion of hpAsia2 or hpEurope strains. In the future, the host range expansion of hpEurope strains among Asian populations, combined with human motility, may have a significant impact on gastric cancer incidence in Asia.     

Phylogeography of the Asian rice gall midge <Emphasis Type="Italic">Orseolia oryzae</Emphasis> (Wood Mason) (Diptera: Cecidomyiidae) in Thailand

The Asian rice gall midge (RGM) Orseolia oryzae (Wood Mason) (Diptera: Cecidomyiidae) is a major pest of rice, leading to yield losses in Thailand and many Asian countries. Despite an increasing number of reported midge outbreaks and the presence of many susceptible rice varieties, only a few studies have focused on the genetic variation of the midges. Therefore, we analyzed the phylogeography among Thai RGM populations covering north, northeast and central Thailand. Two mitochondrial DNA genes, cytochrome C oxidase I (COI) and 12S, and a non-coding repeat region (RR) situated just before COI were amplified. Overall, the haplotype diversity for COI and 12S genes of the Thai population was high, but the nucleotide diversity was quite low. Altogether, the phylogenetic tree and pairwise F _st values indicated that Thai RGM populations recently expanded and were homogeneously distributed throughout the country, except for some populations in the north, which most likely became recently isolated from the main population. Two non-coding repeat motifs, that were recently observed in the mitogenome of RGM in India, were absent in Thai populations and replaced by an 89 bp non-coding sequence. Tandem nucleotide repeats of the sequence TA were also observed. The repeat copy number varied from 2 to 11 and was not correlated with geographical repartition of the midge. Finally, COI barcoding divergence between Indian and Thai populations was high (6.3% in average), giving insights into the potential existence of an RGM species complex in Asia.     

Minisatellite MS32 Alleles Show Population Specificity Among Thai,Chinese, and Japanese

Lineages of structurally related alleles at minisatellite MS32 in human populations show considerable differentiation at the continental level. However, the regional specificity of these lineages remains unknown. We now describe the comparison of allele structures in Thai, Han Chinese, and Japanese populations with lineages previously established for North Europeans and Africans. The great majority of alignable Asian alleles showed their closest structural relative in Asia, with few instances of preferential alignment of Asian with European alleles and only one isolated incident showing a best match with an African allele. Further, there was a strong tendency, most marked for Japanese, for Asian alleles to align preferentially with other alleles from the same population, indicating strong regional specificity of allele lineages. This rapidly evolving minisatellite can therefore serve as a lineage marker for exploring recent events in human population history and dissecting population structure at the fine-scale level, as well as being an extremely informative DNA marker for personal identification.     

Genetic variation,population structure and identification of yellow catfish, <Emphasis Type="Italic">Mystus nemurus</Emphasis> (C&;V) in Thailand using RAPD,ISSR and SCAR marker

Random amplified polymorphic DNA (RAPD) and inter-simple sequence repeat (ISSR) markers were used to investigate the genetic structure of four subpopulations of Mystus nemurus in Thailand. The 7 RAPD and 7 ISSR primers were selected. Of 83 total RAPD fragments, 80 (96.39%) were polymorphic loci, and of 81 total ISSR fragments, 75 (92.59%) were polymorphic loci. Genetic variation and genetic differentiation obtained from RAPD fragments or ISSR fragments showed similar results. Percentage of polymorphic loci (%P), observed number of alleles, effective number of alleles, Nei’s gene diversity (H) and Shannon’s information index revealed moderate to high level of genetic variations within each M. nemurus subpopulation and overall population. High levels of genetic differentiations were received from pairwise unbiased genetic distance (D) and coefficient of differentiation. Mantel test between D or gene flow and geographical distance showed a low to moderate correlation. Analysis of molecular variance indicated that variations among subpopulations were higher than those within subpopulations. The UPGMA dendrograms, based on RAPD and ISSR, showing the genetic relationship among subpopulations are grouped into three clusters; Songkhla (SK) subpopulation was separated from the other subpopulations. The candidate species-specific and subpopulation-specific RAPD fragments were sequenced and used to design sequence-characterized amplified region primers which distinguished M. nemurus from other species and divided SK subpopulation from the other subpopulations. The markers used in this study should be useful for breeding programs and future aquacultural development of this species in Thailand.     

Further resolution of porcine phylogeny in Southeast Asia by Thai mtDNA haplotypes

The mtDNA control region of 72 Thai native pigs and 11 Thai wild boars indigenous to Northern Thailand was comparatively sequenced. In total, 36 nucleotide variations that accounted for 24 haplotypes have been described (TNH01 to TNH20 and TWH01 to TWH04). These haplotypes and further publicly available mtDNA haplotypes were used to assess phylogenetic relationships. Twenty-three of the 24 haplotypes became integrated into the Asian clade of the phylogenetic tree and eight of them recapitulated another major cluster of haplotypes within this clade (Thai haplogroup, THG). Only haplotype TNH01 fit in with the European clade of the phylogenetic tree. An additional analysis using 510 bp of the mtDNA incorporated the THG haplotypes in to clade MTSEA (mountainous and Southeast Asian distribution) to form haplogroup MTSEA-THG. Recently, MTSEA was renamed in MC3. MC3 contains only signatures of pigs scattered across the Indo-Burma Biodiversity Hotspot (IBBH), a region including Thailand to the Kra Isthmus. Here we propose a putative independent porcine domestication event in South-east Asia (SEA). All haplotypes of haplogroup MTSEA-THG have revealed unique and previously unknown nucleotide signatures at positions 24 (nucleotide A) and at positions 183 (nucleotide C) that differentiate them from all other porcine mtDNA haplotypes.     

Molecular epidemiology and outcome of Helicobacter pylori infection in Thailand: a cultural cross roads

ABSTRACT Background. Thailand is at the cultural cross roads between East and South Asia. It has been suggested that this is also the region where the predominant Helicobacter pylori (H. pylori) genotype changes from East Asian to South Asian. Methods. We compared the molecular epidemiology and outcome of H. pylori infections among different ethnic groups in Thailand (Thai, Thai-Chinese and Chinese). H. pylori isolates were genotyped by polymerase chain reaction based on cagA, cag right end junction and vacA genotypes. Results. Ninety-eight isolates from 38 ethnic Thai, 20 ethnic Chinese and 40 Thai-Chinese were categorized into East Asian (45%), South/Central Asian (26%), Western (1%) or mixed type (29%). The East Asian genotype was the most common among Chinese (85%) and Thai-Chinese (55%) (p <.01 compared to ethnic Thai). The ethnicity of the mother among mixed Thai-Chinese marriages predicted the genotype of the child's H. pylori (e.g. when the mother was Chinese, 84% had East Asian type vs. 29% when the mother was Thai) (p <.001). Gastric cancer was common among ethnic Chinese with East Asian genotype (e.g. all Chinese with gastric cancer or peptic ulcer disease had East Asian genotype, whereas only 40% of Chinese with gastritis had this genotype). Conclusions. Immigration, intermarriage and the variety of H. pylori genotypes in Thailand suggest that Thailand is an ideal site for epidemiological studies attempting to relate H. pylori genotypes and host factors to outcome. Our data also support the hypothesis that the primary caretaker of the children is most likely the source of the infection.     

Seventeen novel microsatellite markers from an enriched library of the pest species Bactrocera dorsalis sensu stricto

Seventeen novel microsatellite markers were isolated and characterized from an enriched DNA library of the pest species Bactrocera dorsalis sensu stricto. The polymorphism of these loci was tested in individual fruit flies from a Thai mass‐rearing strain and three wild Asian populations: two from Thailand and one from Myanmar. Allele numbers per locus ranged from two to 15. In the natural populations, a high level of polymorphism was detected in all loci suggesting the usefulness of these markers to quantify the genetic variation in this highly invasive insect.     

Genetic structure and gene flow in a metapopulation of an endangered plant species,Silene tatarica

We investigated the distribution of genetic variation within and between seven subpopulations in a riparian population of Silene tatarica in northern Finland by using amplified fragment length polymorphism (AFLP) markers. A Bayesian approach-based clustering program indicated that the marker data contained not only one panmictic population, but consisted of seven clusters, and that each original sample site seems to consist of a distinct subpopulation. A coalescent-based simulation approach shows recurrent gene flow between subpopulations. Relative high FST values indicated a clear subpopulation differentiation. However, amova analysis and UPGMA-dendrogram did not suggest any hierarchical regional structuring among the subpopulations. There was no correlation between geographical and genetic distances among the subpopulations, nor any correlation between the subpopulation census size and amount of genetic variation. Estimates of gene flow suggested a low level of gene flow between the subpopulations, and the assignment tests proposed a few long-distance bidirectional dispersal events between the subpopulations. No apparent difference was found in within-subpopulation genetic diversity among upper, middle and lower regions along the river. Relative high amounts of linkage disequilibrium at subpopulation level indicated recent population bottlenecks or admixture, and at metapopulation levels a high subpopulation turnover rate. The overall pattern of genetic variation within and between subpopulations also suggested a 'classical' metapopulation structure of the species suggested by the ecological surveys.     

Preliminary DNA identification for the tsunami victims in Thailand

The 2004 Southeast Asia Tsunami killed nearly 5,400 people in Southern Thailand, including foreign tourists and local residents. To recover DNA evidence as much as possible from the seriously decomposed bodies, we explored procedures of sample preparation from both bone and tooth samples as well as both mitochondrial and nuclear markers. Despite having failed to recover enough DNA for nuclear marker typing, we succeeded in obtaining fully informative results for mitochondrial markers （HV1 and HV2） from 258 tooth samples with a success rate of 51% （258/507）. Using an organic DNA extraction method coupled with an ultrafiltration step, we obtained 16 STR （including 13 CODIS loci, one sex discrimination locus, and two Identifiler loci） profiles for 834 samples with a success rate of 79% （834/1,062）. In addition, by comparing the allelic frequencies between the typed samples as a group and other index populations, we conclude that the Thai tsunami victims are a combined group of several populations. Our results provide valuable evidence and protocols for the future forensic practice.     

Genetic attributes of midwife toad (Alytes obstetricans) populations do not correlate with degree of species decline

Genetic diversity is crucial for long‐term population persistence. Population loss and subsequent reduction in migration rate among the most important processes that are expected to lead to a reduction in genetic diversity and an increase in genetic differentiation. While the theory behind this is well‐developed, empirical evidence from wild populations is inconsistent. Using microsatellite markers, we compared the genetic structure of populations of an amphibian species, the midwife toad (Alytes obstetricans), in four Swiss regions where the species has suffered variable levels of subpopulation extirpation. We also quantified the effects of several geographic factors on genetic structure and used a model selection approach to ascertain which of the variables were important for explaining genetic variation. Although subpopulation pairwise F_ST‐values were highly significant even over small geographic scales, neither any of the geographic variables nor loss of subpopulations were important factors for predicting spatial genetic structure. The absence of a signature of subpopulation loss on genetic differentiation may suggest that midwife toad subpopulations function as relatively independent units.     

The freshwater Calanoida (Crustacea: Copepoda) of Thailand

14 Species of freshwater Calanoida in Thailand are reported in this study. The most common species is Neodiaptomus botulifer. Most of these belong to the Southeast Asian species, the remainder are either Indian or East Asian copepods. The zoogeographical demarcation of copepods in Thailand is not very distinct due to its interconnecting riverine and reservoir systems as well as frequent flood. The peculiarity of the Thai Calanoida and their distribution is discussed in a Southeast Asia context. A taxonomic key to the Thai Calanoida is formulated.     

High-resolution analysis of Y-chromosomal polymorphisms reveals signatures of population movements from central Asia and West Asia into India

Linguistic evidence suggests that West Asia and Central Asia have been the two major geographical sources of genes in the contemporary Indian gene pool. To test the nature and extent of similarities in the gene pools of these regions we have collected DNA samples from four ethnic populations of northern India, and have screened these samples for a set of 18 Y-chromosome polymorphic markers (12 unique event polymorphisms and six short tandem repeats). These data from Indian populations have been analysed in conjunction with published data from several West Asian and Central Asian populations. Our analyses have revealed traces of population movement from Central Asia and West Asia into India. Two haplogroups, HG-3 and HG-9, which are known to have arisen in the Central Asian region, are found in reasonably high frequencies (41.7% and 14.3% respectively) in the study populations. The ages estimated for these two haplogroups are less in the Indian populations than those estimated from data on Middle Eastern populations. A neighbour-joining tree based on Y-haplogroup frequencies shows that the North Indians are genetically placed between the West Asian and Central Asian populations. This is consistent with gene flow from West Asia and Central Asia into India.     

Genotypic Distribution of Hepatitis C Virus in Thailand and Southeast Asia

The majority of hepatitis C virus (HCV) infection results in chronic infection, which can lead to liver cirrhosis and hepatocellular carcinoma. Global burden of hepatitis C virus (HCV) is estimated at 150 million individuals, or 3% of the world’s population. The distribution of the seven major genotypes of HCV varies with geographical regions. Since Asia has a high incidence of HCV, we assessed the distribution of HCV genotypes in Thailand and Southeast Asia. From 588 HCV-positive samples obtained throughout Thailand, we characterized the HCV 5’ untranslated region, Core, and NS5B regions by nested PCR. Nucleotide sequences obtained from both the Core and NS5B of these isolates were subjected to phylogenetic analysis, and genotypes were assigned using published reference genotypes. Results were compared to the epidemiological data of HCV genotypes identified within Southeast Asian. Among the HCV subtypes characterized in the Thai samples, subtype 3a was the most predominant (36.4%), followed by 1a (19.9%), 1b (12.6%), 3b (9.7%) and 2a (0.5%). While genotype 1 was prevalent throughout Thailand (27–36%), genotype 3 was more common in the south. Genotype 6 (20.9%) constituted subtype 6f (7.8%), 6n (7.7%), 6i (3.4%), 6j and 6m (0.7% each), 6c (0.3%), 6v and 6xa (0.2% each) and its prevalence was significantly lower in southern Thailand compared to the north and northeast (p = 0.027 and p = 0.030, respectively). Within Southeast Asia, high prevalence of genotype 6 occurred in northern countries such as Myanmar, Laos, and Vietnam, while genotype 3 was prevalent in Thailand and Malaysia. Island nations of Singapore, Indonesia and Philippines demonstrated prevalence of genotype 1. This study further provides regional HCV genotype information that may be useful in fostering sound public health policy and tracking future patterns of HCV spread.     

Y-chromosome and mitochondrial DNA studies on the population structure of the Christmas Island community

Christmas Island is a remote Australian territory located close to the main Indonesian island of Java. Y-chromosome and mitochondrial DNA (mtDNA) markers were used to investigate the genetic structure of the population, which comprises communities of mixed ethnic origin. Analysis of 12 Y-chromosome biallelic polymorphisms revealed a high level of gene diversity and haplotype frequencies that were consistent with source populations in southern China and Southeast Asia. mtDNA hypervariable segment I (HVS-I) sequences displayed high levels of haplotype diversity and nucleotide diversity that were comparable to various Asian populations. Genetic distances revealed extremely low mtDNA differentiation among Christmas Islanders and Asian populations. This was supported by the relatively high proportion of sequence types shared among these populations. The most common mtDNA haplogroups were M* and B, followed by D and F, which are prevalent in East/Southeast Asia. Christmas Islanders of European descent were characterized by the Eurasian haplogroup R*, and a limited degree of admixture was observed. In general, analysis of the genetic data indicated population affinities to southern Chinese (in particular from the Yunnan Province) and Southeast Asia (Thailand, Malaysia, and Cambodia), which was consistent with historical records of settlement. The combined use of these different marker systems provides a useful and appropriate model for the study of contemporary populations derived from different ethnic origins.     

Population structure and genetic differentiation associated with breeding history and selection in tomato (Solanum lycopersicum L.)

Tomato (Solanum lycopersicum L.) has undergone intensive selection during and following domestication. We investigated population structure and genetic differentiation within a collection of 70 tomato lines representing contemporary (processing and fresh-market) varieties, vintage varieties and landraces. The model-based Bayesian clustering software, STRUCTURE, was used to detect subpopulations. Six independent analyses were conducted using all marker data (173 markers) and five subsets of markers based on marker type (single-nucleotide polymorphisms, simple sequence repeats and insertion/deletions) and location (exon and intron sequences) within genes. All of these analyses consistently separated four groups predefined by market niche and age into distinct subpopulations. Furthermore, we detected at least two subpopulations within the processing varieties. These subpopulations correspond to historical patterns of breeding conducted for specific production environments. We found no subpopulation within fresh-market varieties, vintage varieties and landraces when using all marker data. High levels of admixture were shown in several varieties representing a transition in the demarcation between processing and fresh-market breeding. The genetic clustering detected by using the STRUCTURE software was confirmed by two statistics, pairwise F(st) (θ) and Nei's standard genetic distance. We also identified a total of 19 loci under positive selection between processing, fresh-market and vintage germplasm by using an F(st)-outlier method based on the deviation from the expected distribution of F(st) and heterozygosity. The markers and genome locations we identified are consistent with known patterns of selection and linkage to traits that differentiate the market classes. These results demonstrate how human selection through breeding has shaped genetic variation within cultivated tomato.