Western Southeast Asia is hosting one of the world's most diverse faunas, and one of the reasons for this huge diversity is the complex geologic past of the area, increasing the frequency of isolation and expansion events over evolutionary time scale. As an example case, the present study reveals the phylogeny and biogeographic history of the Paracanthocobitis zonalternans species complex, small benthic freshwater fish (Teleostei: Nemacheilidae) that are commonly occurring across western Southeast Asia (from central Myanmar through western and southern Thailand to northern Malaysia). The group is particularly interesting since it occurs in three biogeographic subdivisions (Indian, Indochinese, Malay/Sundaic) and across all of the major biogeographic barriers in the region. Basing on mitochondrial and nuclear sequence data of 93 samples from about 50 localities we found six major clades, most with exclusive geographic distribution. Divergence time dated the origin of the P. zonalternans species complex to early Miocene (17.8 MYA) and a biogeographic analysis identified the Tenasserim region as the ancestral region. From this region the fish spread during periods of lowered global sea level, particularly during late Miocene (11-8 MYA) northwards into all Burmese river basins and southwards into south Thailand and northern Malaysia. Besides lowered global sea level periods, local stream capture events allowed the complex to expand, e.g. into the Mae Klong basin. Strong fragmentations during periods with elevated sea level during the Pliocene and Pleistocene repeatedly restricted populations to refuges and shaped the observed major lineages. Our results document a higher diversity within the P. zonalternans species complex than formerly believed and a strong impact of global sea level on its evolutionary history. Low sea levels promoted dispersal and elevated sea levels fragmentation events. A very similar impact of sea level changes can be expected in all stationary fauna (freshwater and terrestrial) in all non-mountainous coastal regions worldwide.
Frullania subgenus Microfrullania is a clade of ca. 15 liverwort species occurring in Australasia, Malesia, and southern South America. We used combined nuclear and chloroplast sequence data from 265 ingroup accessions to test species circumscriptions and estimate the biogeographic history of the subgenus. With dense infra-specific sampling, we document an important role of long-distance dispersal in establishing phylogeographic patterns of extant species. At deeper time scales, a combination of phylogenetic analyses, divergence time estimation and ancestral range estimation were used to reject vicariance and to document the role of long-distance dispersal in explaining the evolution and biogeography of the clade across the southern Hemisphere. A backbone phylogeny for the subgenus is proposed, providing insight into evolution of morphological patterns and establishing the basis for an improved sectional classification of species within Microfrullania. Several species complexes are identified, the presence of two undescribed but genetically and morphologically distinct species is noted, and previously neglected names are discussed.
Platanthera is one of the largest genera of temperate orchids in the Holarctic and exemplifies a lineage that has adaptively radiated into diverse habitats within North America, Asia, Europe, North Africa, Borneo, and Sarawak. Major centers of diversity in this genus are North America and eastern Asia. Despite its diversity, a thorough phylogenetic hypothesis for the genus is lacking because no studies have yet sampled taxa exhaustively or developed a robust molecular toolkit. While there is strong evidence that suggests monophyly of subgenus Limnorchis, most taxa in this group have not been included in a phylogenetic analysis. In this study, we developed a new toolkit for Platanthera consisting of genomic information from 617 low-copy nuclear loci. Using a targeted enrichment approach, we collected high-throughput sequence data in 23 accessions of nine of the 12 diploid species of subgenus Limnorchis and outgroup species across Platanthera. A maximum likelihood analysis resolved a strongly supported monophyletic clade for subgenus Limnorchis. Ancestral biogeographic reconstruction indicated that subgenus Limnorchis originated in western North America ca. 3-4.5 Mya from an ancestor that was widespread in western North America and eastern Asia and subsequently diversified in western North America, followed by dispersal of some species to eastern North America. Our results indicate complex biogeographic connections between Asia and North America, and therefore it suggests that Platanthera is a suitable system to test biogeographic hypotheses over time and space in the Holarctic. Our results are also expected to facilitate further study of diversification and biogeographic spread across Platanthera and lay the groundwork for understanding independent origins, biogeography, and morphological diversification of polyploid species within subgenus Limnorchis.
The genus Rana, notably diversified in Oriental regions from China to Southeast Asia, includes a group of cascade frogs assigned to subgenera Odorrana and Eburana. Among them, R. ishikawae and the R. narina complex represent the northernmost members occurring from Taiwan to the Ryukyu Archipelago of Japan. Relationships of these frogs with the continental members, as well as the history of their invasions to islands, have been unclear. The taxonomic status of Odorrana and related genera varies among authors and no phylogenetic reassessment has been done. Using partial sequences of mitochondrial 12S and 16S rRNA genes, we estimated phylogenetic relationships among 17 species of the section Hylarana including Odorrana and Eburana, and related species from the Ryukyus, Taiwan, China, Thailand, Malaysia, and Indonesia. We estimate that (1) Odorrana is monophyletic and encompasses species of Eburana and R. hosii, which is now placed in Chalcorana, (2) the ancestor of R. ishikawae separated from other Rana in the middle to late Miocene prior to its entry to the Ryukyu Archipelago, (3) the ancestor of the R. narina complex later diversified in continental Asia, and invaded the Ryukyu Archipelago through Taiwan, (4) the R. narina complex attained its current distribution within the Ryukyus through niche segregations, and (5) vicariance of R. hosii between Malay Peninsula and Borneo occurred much later than the divergence events in the R. narina complex. Current subgeneric classification of Rana, at least of Southeast Asian members, requires full reassessment in the light of phylogenetic relationships.
Inferring interfamilial relationships within the eudicot order Ericales has remained one of the more recalcitrant problems in angiosperm phylogenetics, likely due to a rapid, ancient radiation. As a result, no comprehensive time-calibrated tree or biogeographical analysis of the order has been published. Here, we elucidate phylogenetic relationships within the order and then conduct time-dependent biogeographical and diversification analyses by using a taxon and locus-rich supermatrix approach on one-third of the extant species diversity calibrated with 23 macrofossils and two secondary calibration points. Our results corroborate previous studies and also suggest several new but poorly supported relationships. Newly suggested relationships are: (1) holoparasitic Mitrastemonaceae is sister to Lecythidaceae, (2) the clade formed by Mitrastemonaceae + Lecythidaceae is sister to Ericales excluding balsaminoids, (3) Theaceae is sister to the styracoids + sarracenioids + ericoids, and (4) subfamilial relationships with Ericaceae suggest that Arbutoideae is sister to Monotropoideae and Pyroloideae is sister to all subfamilies excluding Arbutoideae, Enkianthoideae, and Monotropoideae. Our results indicate Ericales began to diversify 110 Mya, within Indo-Malaysia and the Neotropics, with exchange between the two areas and expansion out of Indo-Malaysia becoming an important area in shaping the extant diversity of many families. Rapid cladogenesis occurred along the backbone of the order between 104 and 106 Mya. Jump dispersal is important within the order in the last 30 My, but vicariance is the most important cladogenetic driver of disjunctions at deeper levels of the phylogeny. We detect between 69 and 81 shifts in speciation rate throughout the order, the vast majority of which occurred within the last 30 My. We propose that range shifting may be responsible for older shifts in speciation rate, but more recent shifts may be better explained by morphological innovation.
The mountains of Borneo are well known for their high endemicity and historical role in preserving Southeast Asian rainforest biodiversity, but the diversification of populations inhabiting these mountains is poorly studied. Here we examine the genetic structure of 12 Bornean montane passerines by comparing complete mtDNA ND2 gene sequences of populations spanning the island. Maximum likelihood and Bayesian phylogenetic trees and haplotype networks are examined for common patterns that might signal important historical events or boundaries to dispersal. Morphological and ecological characteristics of each species are also examined using phylogenetic generalized least-squares (PGLS) for correlation with population structure. Populations in only four of the 12 species are subdivided into distinct clades or haplotype groups. Although this subdivision occurred at about the same time in each species (ca. 0.6-0.7Ma), the spatial positioning of the genetic break differs among the species. In two species, northeastern populations are genetically divergent from populations elsewhere on the island. In the other two species, populations in the main Bornean mountain chain, including the northeast, are distinct from those on two isolated peaks in northwestern Borneo. We suggest different historical forces played a role in shaping these two distributions, despite commonality in timing. PGLS analysis showed that only a single characteristic-hand-wing index-is correlated with population structure. Birds with longer wings, and hence potentially more dispersal power, have less population structure. To understand historical forces influencing montane population structure on Borneo, future studies must compare populations across the entirety of Sundaland.
Mitochondrial genomes provided the first widely used sequences that were sufficiently informative to resolve relationships among animals across a wide taxonomic domain, from within species to between phyla. However, mitogenome studies supported several anomalous relationships and fell partly out of favour as sequencing multiple, independent nuclear loci proved to be highly effective. A tendency to blame mitochondrial DNA (mtDNA) has overshadowed efforts to understand and ameliorate underlying model misspecification. Here we find that influential assessments of the infidelity of mitogenome phylogenies have often been overstated, but nevertheless, substitution saturation and compositional non-stationarity substantially mislead reconstruction. We show that RY coding the mtDNA, excluding protein-coding 3rd codon sites, partitioning models based on amino acid hydrophobicity and enhanced taxon sampling improve the accuracy of mitogenomic phylogeny reconstruction for placental mammals, almost to the level of multi-gene nuclear datasets. Indeed, combined analysis of mtDNA with 3-fold longer nuclear sequence data either maintained or improved upon the nuclear support for all generally accepted clades, even those that mtDNA alone did not favour, thus indicating "hidden support". Confident mtDNA phylogeny reconstruction is especially important for understanding the evolutionary dynamics of mitochondria themselves, and for merging extinct taxa into the tree of life, with ancient DNA often only accessible as mtDNA. Our ancient mtDNA analyses lend confidence to the relationships of three extinct megafaunal taxa: glyptodonts are nested within armadillos, the South American ungulate, Macrauchenia is sister to horses and rhinoceroses, and sabre-toothed and scimitar cats are the monophyletic sister-group of modern cats.
As opposed to angiosperms, moss species richness is similar among tropical regions of the world, in line with the hypothesis that tropical bryophytes are extremely good dispersers. Here, we reconstructed the phylogeny of the pantropical moss genus Pelekium to test the hypothesis that high migration rates erase any difference in species richness among tropical regions. In contrast with this hypothesis, several species considered to have a pantropical range were resolved as a complex of species with a strong geographic structure. Consequently, a significant phylogeographical signal was found in the data, evidencing that cladogenetic diversification within regions takes place at a faster rate than intercontinental migration. The shape of the Pelekium phylogeny, along with the selection of a constant-rate model of diversification among species in the genus, suggests, however, that the cladogenetic speciation patterns observed in Pelekium are not comparable to some of the spectacular examples of tropical radiations reported in angiosperms. Rather, the results presented here point to the constant accumulation of diversity through time in Pelekium. This, combined with evidence for long-distance dispersal limitations in the genus, suggests that the similar patterns of species richness among tropical areas are better explained in terms of comparable rates of diversification across tropical regions than by the homogenization of species richness by recurrent migrations.
This study provides a first description of the phylogeographic patterns and evolutionary history of two species of the mudskipper genus Periophthalmus. These amphibious gobies are distributed throughout the whole Indo-Pacific region and Atlantic coast of Africa, in peritidal habitats of soft-bottom coastal ecosystems. Three sequence datasets of two widely distributed species, Periophthalmus argentilineatus and P. kalolo, were obtained by amplifying and sequencing two mtDNA markers (D-loop and 16S rDNA) and the nDNA rag1 region. The three datasets were then used to perform phylogeographic, demographic and population genetic analyses. Our results indicate that tectonic events and past climatic oscillations strongly contributed to shape present genetic differentiation, phylogeographic and demographic patterns. We found support for the monophyly of P. kalolo, and only shallow genetic differentiation between East-African and Indo-Malayan populations of this species. However, our collections of the morphospecies P. argentilineatus include three molecularly distinct lineages, one of them more closely related to P. kalolo. The presence of Miocenic timings for the most recent common ancestors of some of these morphologically similar clades, suggests the presence of strong stabilising selection in mudskippers' habitats. At population level, demographic analyses and palaeoecological records of mangrove ecosystems suggest that Pleistocene bottlenecks and expansion plus secondary contact events of the studied species were associated with recurrent sea transgressions during interglacials, and sea regressions or stable regimes during glacials, respectively.
The genus Myotis includes the largest number of species in the family Vespertilionidae (Chiroptera), and its members are distributed throughout most of the world. To re-evaluate the phylogenetic position of East Asian Myotis species with respect to Myotis species worldwide, we analyzed mitochondrial gene sequences of NADH dehydrogenase subunit 1 and cytochrome b from 24 East Asian individuals as well as 42 vespertilionid bats determined previously. The results suggest that: (1) some individuals having the same species name in Europe and Japan do not form a monophyletic clade, indicating that some bat species exhibit morphological convergence, (2) Japanese Myotis mystacinus forms a sister relationship with Myotis brandtii (Palaearctic), and both species are included in the American clade implying that an ancestor of these species originated in North America, and (3) the Black whiskered bat, Myotis pruinosus, is endemic to Japan and forms sister relationships with Myotis yanbarensis and Myotis montivagus collected from Okinawa (Japan) and Selangor (Malaysia), respectively, implying that M. pruinosus originated from the south. The systematics of Japanese and East Asian Myotis bats were revisited by considering their phylogenetic relationships. Our study provides the first extensive phylogenetic hypothesis of the genus Myotis that includes East Asian and Japanese species.
We examined the phylogeny and biogeography of the glassperch family Ambassidae (Teleostei), which is widely distributed in the freshwater, brackish and marine coastal habitats across the Indo-West Pacific region. We first built a comprehensive time-calibrated phylogeny of Ambassidae using five genes. We then used this tree to reconstruct the evolution of the salinity preference and ancestral areas. Our results indicate that the two largest genera of Ambassidae, Ambassis and Parambassis, are each not monophyletic. The most recent common ancestor of Ambassidae was freshwater adapted and lived in Australia about 56 million years ago. Three independent freshwater-to-marine transitions are inferred, but no marine-to-freshwater ones. To explain the distribution of ambassids, we hypothesise two long-distance marine dispersal events from Australia. A first event was towards Southeast Asia during the early Cenozoic, followed by a second one towards Africa during mid-Cenozoic. The phylogenetic signal associated with the salinity adaptation of these events was not detected, possibly because of the selective extinction of intermediate marine lineages. The Ambassidae shares two characteristics with other freshwater fish groups distributed in continental regions surrounding the Indian Ocean: They are too young to support the hypothesis that their distribution is the result of the fragmentation of Gondwana, but they did not retain the phylogenetic signal of their marine dispersal.
To elucidate the evolution of one of the most species-rich ant-plant symbiotic systems, the association between Crematogaster (Myrmicinae) and Macaranga (Euphorbiaceae) in South-East Asia, we conducted a phylogenetic analysis of the ant partners. For the phylogenetic analysis partial mitochondrial cytochrome oxidase I and II were sequenced and Maximum Parsimony analysis was performed. The analyzed Crematogaster of the subgenus Decacrema fell into three distinct clades which are also characterized by specific morphological and ecological traits (queen morphology, host-plants, and colony structure). Our results supported the validity of our currently used morphospecies concept for Peninsula Malaysia. However, on a wider geographic range (including North and North-East Borneo) some morphospecies turned out to be species complexes with genetically quite distinct taxa. Our phylogenetic analysis and host association studies do not indicate strict cocladogenesis between the subgenus Decacrema and their Macaranga host-plants because multiple ant taxa occur on quite distinct host-plants belonging to different clades within in the genus Macaranga. These results support the view that host-shifting or host-expansion is common in the ants colonizing Macaranga. Additionally, the considerable geographic substructuring found in the phylogenetic trees of the ants suggests that allopatric speciation has also played a role in the diversification and the current distribution of the Decacrema ants.
Myrteae (c. 2500 species; 51 genera) is the largest tribe of Myrtaceae and an ecologically important groups of angiosperms in the Neotropics. Systematic relationships in Myrteae are complex, hindering conservation initiatives and jeopardizing evolutionary modelling. A well-supported and robust phylogenetic hypothesis was here targeted towards a comprehensive understanding of the relationships within the tribe. The resultant topology was used as a base for key evolutionary analyses such as age estimation, historical biogeography and diversification rate patterns. One nuclear (ITS) and seven chloroplast (psbA-trnH, matK, ndhF, trnl-trnF, trnQ-rps16, rpl16 and rpl32-trnL) DNA regions for 115 taxa representing 46 out of the 51 genera in the tribe were accessed and analysed using maximum likelihood and Bayesian inference tools for phylogenetic reconstruction. Dates of diversification events were estimated and contrasted using two distinct fossil sets (macro and pollen) in BEAST. The subsequent dated phylogenies were compared and analysed for biogeographical patterns using BioGeoBEARS and diversification rates using BAMM. Myrteae phylogeny presents strong statistical support for three major clades within the tribe: Australasian group, Myrtus group and Main Neotropical Lineage. Dating results from calibration using macrofossil are an average of 20 million years older and show an early Paleocene origin of Myrteae, against a mid-Eocene one from the pollen fossil calibration. Biogeographic analysis shows the origin of Myrteae in Zealandia in both calibration approaches, followed by a widespread distribution throughout the still-linked Gondwana continents and diversification of Neotropical endemic lineages by later vicariance. Best configuration shift indicates three points of acceleration in diversification rates, all of them occurring in the Main Neotropical Lineage. Based on the reconstructed topology, several new taxonomic placements were recovered, including: the relative position of Myrtus communis, the placement of the Blepharocalyx group, the absence of generic endemism in the Caribbean, and the paraphyletism of the former Pimenta group. Distinct calibration approaches affect biogeography interpretation, increasing the number of necessary long distance dispersal events in the topology with older nodes. It is hypothesised that biological intrinsic factors such as modifications of embryo type and polyploidy might have played a role in accelerating shifts of diversification rates in Neotropical lineages. Future perspectives include formal subtribal classification, standardization of fossil calibration approaches and better links between diversification shifts and trait evolution.
Several anuran groups of Laurasian origin are each co-distributed in four isolated regions of the Northern Hemisphere: central/southern Europe and adjacent areas, Korean Peninsula and adjacent areas, Indo-Malaya, and southern North America. Similar distribution patterns have been observed in diverse animal and plant groups. Savage [Savage, J.M., 1973. The geographic distribution of frogs: patterns and predictions. In: Vial, J.L. (Ed.), Evolutionary Biology of the Anurans. University of Missouri Press, Columbia, pp. 351-445] hypothesized that the Miocene global cooling and increasing aridities in interiors of Eurasia and North America caused a southward displacement and range contraction of Laurasian frogs (and other groups). We use the frog genus Bombina to test Savage's biogeographical hypothesis. A phylogeny of Bombina is reconstructed based on three mitochondrial and two nuclear gene fragments. The genus is divided into three major clades: an Indo-Malaya clade includes B. fortinuptialis, B. lichuanensis, B. maxima, and B. microdeladigitora; a European clade includes B. bombina, B. pachypus, and B. variegata; and a Korean clade contains B. orientalis. The European and Korean clades form sister-group relationship. Molecular dating of the phylogenetic tree using the penalized likelihood and Bayesian analyses suggests that the divergence between the Indo-Malaya clade and other Bombina species occurred 5.9-28.6 million years ago. The split time between the European clade and the Korean clade is estimated at 5.1-20.9 million years ago. The divergence times of these clades are not significantly later than the timing of Miocene cooling and drying, and therefore can not reject Savage's hypothesis. Some other aspects of biogeography of Bombina also are discussed. The Korean Peninsula and the Shandong Peninsula might have supplied distinct southern refugia for B. orientalis during the Pleistocene glacial maxima. In the Indo-Malaya clade, the uplift of the Tibetan Plateau might have promoted the split between B. maxima and the other species.
Taxonomic relationships within the Old World fruit bat genus, Cynopterus, have been equivocal for the better part of a century. While nomenclature has been revised multiple times on the basis of phenotypic characters, evolutionary relationships among taxa representing the entire geographic range of the genus have not been determined. We used mitochondrial DNA sequence data to infer phylogenetic relationships among the three most broadly distributed members of the genus: C. brachyotis, C. horsfieldi, and C. sphinx, and to assess whether C. brachyotis represents a single widespread species, or a complex of distinct lineages. Results clearly indicate that C. brachyotis is a complex of lineages. C. sphinx and C. horsfieldi haplotypes formed monophyletic groups nested within the C. brachyotis species complex. We identified six divergent mitochondrial lineages that are currently referred to C. brachyotis. Lineages from India, Myanmar, Sulawesi, and the Philippines are geographically well-defined, while in Malaysia two lineages, designated Sunda and Forest, are broadly sympatric and may be ecologically distinct. Demographic analyses of the Sunda and Forest lineages suggest strikingly different population histories, including a recent and rapid range expansion in the Sunda lineage, possibly associated with changes in sea levels during the Pleistocene. The resolution of the taxonomic issues raised in this study awaits combined analysis of morphometric characters and molecular data. However, since both the Indian and Malaysian Forest C. brachyotis lineages are apparently ecologically restricted to increasingly fragmented forest habitat, we suggest that reevaluation of the conservation status of populations in these regions should be an immediate goal.
Coptotermes gestroi, the Asian subterranean termite (AST), is an economically important structural and agricultural pest that has become established in many areas of the world. For the first time, phylogeography was used to illuminate the origins of new found C. gestroi in the US Commonwealth of Puerto Rico; Ohio, USA; Florida, USA; and Brisbane, Australia. Phylogenetic relationships of C. gestroi collected in indigenous locations within Malaysia, Thailand, and Singapore as well as from the four areas of introduction were investigated using three genes (16S rRNA, COII, and ITS) under three optimality criteria encompassing phenetic and cladistic assumptions (maximum parsimony, maximum likelihood, and neighbor-joining). All three genes showed consistent support for a close genetic relationship between C. gestroi samples from Singapore and Ohio, whereas termite samples from Australia, Puerto Rico, and Key West, FL were more closely related to those from Malaysia. Shipping records further substantiated that Singapore and Malaysia were the likely origin of the Ohio and Australia C. gestroi, respectively. These data provide support for using phylogeography to understand the dispersal history of exotic termites. Serendipitously, we also gained insights into concerted evolution in an ITS cluster from rhinotermitid species in two genera.
By investigating genealogical relationships, we estimated the phylogenetic history and biogeography in the megophryid genus Leptobrachium (sensu lato, including Vibrissaphora) from southern China, Indochina, Thailand and the Sundaland. The genealogical relationships among the 30 named and unnamed taxa were estimated using 2009 bp of sequences from the mitochondrial DNA genes 12S rRNA, tRNA(val), and 16S rRNA using maximum parsimony, maximum likelihood, and Bayesian inference methods. The genus Leptobrachium was a well-supported monophyletic group that contained two major clades. One clade had three subclades primarily from disjunct regions including Borneo, Peninsular Malaysia and Java, and Thailand. The Bornean subclade included one species each from the Philippines and Sumatra. The other major clade consisted of two subclades, one from Indochina and the other from southern China (Vibrissaphora). Divergence times estimated an old evolutionary history of each subclade, one that could not be explained by the geohistory of Southeast Asian major landmasses.
Wolbachia pipientis is a diverse, ubiquitous and most prevalent intracellular bacterial group of alpha-Proteobacteria that is concerned with many biological processes in arthropods. The coconut hispine beetle (CHB), Brontispa longissima (Gestro) is an economically important pest of palm cultivation worldwide. In the present study, we comprehensively surveyed the Wolbachia-infection prevalence and mitochondrial DNA (mtDNA) polymorphism in CHB from five different geographical locations, including China's Mainland and Taiwan, Vietnam, Thailand, Malaysia and Indonesia. A total of 540 sequences were screened in this study through three different genes, i.e., cytochrome oxidase subunit I (COI), Wolbachia outer surface protein (wsp) and multilocus sequencing type (MLST) genes. The COI genetic divergence ranges from 0.08% to 0.67%, and likewise, a significant genetic diversity (π = 0.00082; P = 0.049) was noted within and between all analyzed samples. In the meantime, ten different haplotypes (H) were characterized (haplotype diversity = 0.4379) from 21 different locations, and among them, H6 (46 individuals) have shown a maximum number of population clusters than others. Subsequently, Wolbachia-prevalence results indicated that all tested specimens of CHB were found positive (100%), which suggested that CHB was naturally infected with Wolbachia. Wolbachia sequence results (wsp gene) revealed a high level of nucleotide diversity (π = 0.00047) under Tajima's D test (P = 0.049). Meanwhile, the same trend of nucleotide diversity (π = 0.00041) was observed in Wolbachia concatenated MLST locus. Furthermore, phylogenetic analysis (wsp and concatenated MLST genes) revealed that all collected samples of CHB attributed to same Wolbachia B-supergroup. Our results strongly suggest that Wolbachia bacteria and mtDNA were highly concordant with each other and Wolbachia can affect the genetic structure and diversity within the CHB populations.