Average Pairwise Difference Research Articles

Genetic structure of silver pomfret (Pampus argenteus (Euphrasen, 1788)) in the Arabian Sea, Bay of Bengal, and South China Sea as indicated by mitochondrial COI gene sequences

Summary The silver pomfret (Pampus argenteus) is a commercially important species native to the Indo-western Pacific area. In the present study, the genetic structure of five P. argenteus populations from the Arabian Sea, Bay of Bengal, and South China Sea were investigated using sequences of the mitochondrial cytochrome c oxidase subunit I (COI) gene. All populations showed high levels of haplotype diversity (h) and low levels of nucleotide diversity (π). Phylogenetic trees indicated that the haplotypes could be divided into three clusters, indicating the geographic distribution of the five populations. Analysis of molecular variance (amova) indicated the intra-region variation among these populations to be highly significant. The average pairwise differences and Fst values among the three regions were also found to be significant. These results suggest that P. argenteus in these regions has strong geographic genetic structure. The present results provide new information for the genetic assessment, fishery management, and conservation of this species.

Population genetic structure and demographic history of Pampus argenteus in the Indo-West Pacific inferred from mitochondrial cytochrome b sequences

Abstract Population genetic structure and historical demography of Pampus argenteus in three areas (the South China Sea, the Arabian Sea and the Bay of Bengal) of the Indo-West Pacific were evaluated using mitochondrial cytochrome b gene sequences. The net Tumura and Nei genetic distance among the three regions ranged from 0.023 to 0.034. The fixation index analysis showed that significant genetic divergence occurred among the populations from different geographical regions. Analysis of molecular variance (AMOVA) reveals that intra-region variation among the populations was highly significant. The average pairwise differences and φST values (0.885–0.939) among the three regions were significant. The populations have high haplotype diversity (Hd, 0.736–0.947) and low nucleotide diversity (π, 0.141–0.586%). The results suggest that P. argenteus has strong geographic genetic structures. In addition, mismatch distribution analyses and tests of neutral evolution imply that a late Pleistocene population expansion may have happened in all populations.

Influence of finite‐sites mutation, population subdivision and sampling schemes on patterns of nucleotide polymorphism for species with molecular hyperdiversity

Molecular hyperdiversity has been documented in viruses, prokaryotes and eukaryotes. Such organisms undermine the assumptions of the infinite-sites mutational model, because multiple mutational events at a site comprise a non-negligible portion of polymorphisms. Moreover, different sampling schemes of individuals from species with subdivided populations can profoundly influence resulting patterns and interpretations of molecular variation. Inspired by molecular hyperdiversity in the nematode Caenorhabditis sp. 5, which exhibits average pairwise differences among synonymous sites of >5% as well as modest population structure, we investigated via coalescent simulation the joint effects of a finite-sites mutation (FSM) process and population subdivision on the variant frequency spectrum. From many demes interconnected through a stepping-stone migration model, we constructed local samples from a single deme, pooled samples from several demes and scattered samples of a single individual from numerous demes. Compared with a single panmictic population at equilibrium, we find that high population mutation rates induce a deficit of rare variants (positive Tajima's D) under a FSM model. Population structure also induces such a skew for local samples when migration is high and for pooled samples when migration is low. Contrasts of sampling schemes for C. sp. 5 imply high mutational input coupled with high migration. We propose that joint analysis of local, pooled and scattered samples for species with subdivided populations provides a means of improving inference of demographic history, by virtue of the partially distinct patterns of polymorphism that manifest when sequences are analyzed according to differing sampling schemes.

Genomic impacts of chromosomal inversions in parapatricDrosophilaspecies

Chromosomal inversions impact genetic variation and facilitate speciation in part by reducing recombination in heterokaryotypes. We generated multiple whole-genome shotgun sequences of the parapatric species pair Drosophila pseudoobscura and Drosophila persimilis and their sympatric outgroup (Drosophila miranda) and compared the average pairwise differences for neutral sites within, just outside and far outside of the three large inversions. Divergence between D. pseudoobscura and D. persimilis is high inside the inversions and in the suppressed recombination regions extending 2.5 Mb outside of inversions, but significantly lower in collinear regions further from the inversions. We observe little evidence of decreased divergence predicted to exist in the centre of inversions, suggesting that gene flow through double crossovers or gene conversion is limited within the inversion, or selection is acting within the inversion to maintain divergence in the face of gene flow. In combination with past studies, we provide evidence that inversions in this system maintain areas of high divergence in the face of hybridization, and have done so for a substantial period of time. The left arm of the X chromosome and chromosome 2 inversions appear to have arisen in the lineage leading to D. persimilis approximately 2 Ma, near the time of the split of D. persimilis-D. pseudoobscura-D. miranda, but likely fixed within D. persimilis much more recently, as diversity within D. persimilis is substantially reduced inside and near these two inversions. We also hypothesize that the inversions in D. persimilis may provide an empirical example of the 'mixed geographical mode' theory of inversion origin and fixation, whereby allopatry and secondary contact both play a role.

Three-Dimensional Solution Structure of Bottromycin A2: A Potent Antibiotic Active against Methicillin-Resistant Staphylococcus aureus and Vancomycin-Resistant Enterococci

The three-dimensional (3D) structure of bottromycin A(2), a natural anti-methicillin-resistant Staphylococcus aureus (MRSA) and anti-vancomycin-resistant Enterococci (VRE) agent consisting of seven amino acids, has been investigated through NMR spectroscopy. On the basis of 57 experimental constraints, a total of 34 converged structures were obtained. The average pairwise atomic root mean square difference is 0.74±0.59 Å for all heavy atoms. The resulting structure indicates an interesting feature in that the three C-terminal residues of bottromycin A(2) fold back on the 12-membered cyclic skeleton made by the four N-terminal residues. Thus, MePro(2) and Thia-β-Ala-OMe(7), modification of which significantly affects the antibacterial activities of bottromycin A(2), are located on one side of its 3D structure. These distinct structural features might be important for the binding of bottromycin A(2) with the bacterial ribosome.

Genetic diversity of simian immunodeficiency virus encoding HIV-1 reverse transcriptase persists in macaques despite antiretroviral therapy.

The impact of antiretroviral therapy (ART) on the genetics of simian immunodeficiency virus (SIV) or human immunodeficiency virus (HIV) populations has been incompletely characterized. We analyzed SIV genetic variation before, during, and after ART in a macaque model. Six pigtail macaques were infected with an SIV/HIV chimeric virus, RT-SHIV(mne), in which SIV reverse transcriptase (RT) was replaced by HIV-1 RT. Three animals received a short course of efavirenz (EFV) monotherapy before combination ART was started. All macaques received 20 weeks of tenofovir, emtricitabine, and EFV. Plasma virus populations were analyzed by single-genome sequencing. Population diversity was measured by average pairwise difference, and changes in viral genetics were assessed by phylogenetic and panmixia analyses. After 20 weeks of ART, viral diversity was not different from pretherapy viral diversity despite more than 10,000-fold declines in viremia, indicating that, within this range, there is no relationship between diversity and plasma viremia. In two animals with consistent SIV RNA suppression to <15 copies/ml during ART, there was no evidence of viral evolution. In contrast, in the four macaques with viremias >15 copies/ml during therapy, there was divergence between pre- and during-ART virus populations. Drug resistance mutations emerged in two of these four animals, resulting in virologic failure in the animal with the highest level of pretherapy viremia. Taken together, these findings indicate that viral diversity does not decrease with suppressive ART, that ongoing replication occurs with viremias >15 copies/ml, and that in this macaque model of ART drug resistance likely emerges as a result of incomplete suppression and preexisting drug resistance mutations.

Developing equine mtDNA profiling for forensic application

Horse mtDNA profiling can be useful in forensic work investigating degraded samples, hair shafts or highly dilute samples. Degraded DNA often does not allow sequencing of fragments longer than 200 nucleotides. In this study we therefore search for the most discriminatory sections within the hypervariable horse mtDNA control region. Among a random sample of 39 horses, 32 different sequences were identified in a stretch of 921 nucleotides. The sequences were assigned to the published mtDNA types A-G, and to a newly labelled minor type H. The random match probability within the analysed samples is 3.61%, and the average pairwise sequence difference is 15 nucleotides. In a "sliding window" analysis of 200-nucleotide sections of the mtDNA control region, we find that the known repetitive central motif divides the mtDNA control region into a highly diverse segment and a markedly less discriminatory segment.

Spatial Patterns of Genetic Differentiation in Brachionus calyciflorus Species Complex

To understand spatial patterns of genetic differentiation in the Brachionus calyciflorus species complex, the rDNA ITS regions of 124 rotifer clones collected from eight geographic regions in East China were sequenced and analyzed. A total of 73 haplotypes were defined, and were grouped into three clades by the phylogenetic trees. The divergences of ITS sequence among the three clades ranged from 4.2% to 25.3%, indicating the occurrence of three cryptic species (cryptic species I, cryptic species II and cryptic species III). The sympatric co-occurrence of these cryptic species was observed at Guangzhou, Danzhou and Wuhu. The nucleotide diversity in B. calyciflorus complex was higher, and so was haplotype diversity. Within cryptic species III, the higher genetic differentiation might be attributed to fragmentation and restricted gene flow. In addition, the rapid increase of colonizers after an ambiguous glacial bottleneck might interfere with effective gene flow and accelerate genetic differentiation of geographic populations. The nested clade analysis suggested to some extent the geographic structure, and the non-significant correlation existed between the geographic distance and the average pairwise difference between populations (Da) as well as Fst. The Younger Dryas(YD) Event after the last full glacial period might has considerable effect on the patterns of geographic structure and distribution of cryptic B. calyciflorus species. During the YD event, the three cryptic species probably retreated and co-occurred in multiple relict refugia to the south of Nanling area. Long-distance colonization possibly coupled with subsequent fragmentation might be the main reason of the current patterns of geographic structure. The geographic distribution of B. calyciflorus complex on a global scale might be attributed to the motion of continental plates.

Population genetic structure and historical demography of grey mullet, Mugil cephalus, along the coast of China, inferred by analysis of the mitochondrial control region

Population genetic structure and historical demography of grey mullet, Mugil cephalus, along the coast of China were examined with a 389-bp segment of mtDNA control region. In total, 117 samples were collected from seven locations and 77 haplotypes were obtained. No haplotype was concurrently presented in all the samples, and samples in Southern and Northern regions had distinctly different haplotypes. Sorting of the data cladistically into phylogenetic trees indicated that these haplotypes can be divided into two main clades. Analysis of molecular variance (AMOVA) indicated significant genetic differentiation among subpopulations in different regions ( φ CT = 0.9505, p < 0.01). The average pairwise differences and φ ST values (0.8943–0.9464) between regional populations were significant. The net Tumura and Nei genetic distance (0.5577 ± 0.1421) between the two groups was large. The deep and unique divergence between the two groups suggested that each group is likely to represent a distinct species. In contrast to the great genetic divergence on a broad scale, small but significant population subdivisions at a less spatial scale were also detected both in Southern and Northern regions, which support the conclusion that M. cephalus has limited dispersal potential. High haplotype diversity ( h = 0.8716 ± 0.1462) and low nucleotide diversity ( π = 0.0046 ± 0.0030) were found among the three Northern samples. Among the four Southern samples, however, extraordinarily high levels of both haplotype diversity ( h = 0.9886 ± 0.0057) and nucleotide diversity ( π = 0.0425 ± 0.0212) were detected. The different polymorphisms suggest an apparent distinguishing demographic history between regionally Southern and Northern populations, though mismatch distribution analyses and neutrality tests implied a late Pleistocene population expansion of both Northern and Southern populations.

Rate heterogeneity in the evolution of Helicobacter pylori and the behavior of homoplastic sites

Helicobacter pylori are bacteria with substantial inter-strain variability and phylogenetic reconstructions of sequence data from the organism have common homoplastic sites. Although frequent recombination events have been proposed to contribute to the variation, the effects of nucleotide substitution rate heterogeneities on the reconstruction of H. pylori genealogies have not been studied. We analyzed the substitution pattern of a housekeeping gene, a homologue of the ribonuclease reductase gene ( rnr), to characterize rate heterogeneities between 11 H. pylori isolates. Evidence of limited recombination was demonstrated by the Sawyer's runs test, but the homoplasy test and site-by-site compatibility tests indicated frequent recombination events. Within the 1935 nucleotide gene, 292 sites were polymorphic with an average pair-wise difference of 5.01%. Xia's distances for amino acids at non-synonymous codon substitution sites were smaller at homoplastic sites than at sites that were not homoplastic. Transitions were significantly more common among homoplastic than among non-homoplastic nucleotide substitutions. Simulations of evolution with or without recombination indicated the transition/transversion ratio is expected to be higher in homoplastic sites with no recombination. Despite evidence of recombination, analyses of the rnr genealogy does not show a random tree but rather base substitution behaviors characteristic of both recombination and substitution saturation at some sites. Analyses of sequences in the H. pylori multilocus sequence-typing database provided similar evidence for substitution saturation in multiple housekeeping genes.

DNA based cladograms augment the discovery of a new Ips species from China (Coleoptera: Curculionidae: Scolytinae).

The implementation of DNA in taxonomic study is in its infancy because the association of the amount and type of nucleotide change with species boundaries has not been fully examined for most taxa. Mitochondrial cytochrome c oxidase I (COI) nucleotide data is currently the most popular molecular marker for delimiting species boundaries and a standard pair-wise nucleotide divergence between groups of individuals has been suggested for the recognition of new species. It is unlikely that such a standard would be applicable across animal species, but the association of the amount and type of nucleotide change with species boundaries could help with the establishment of a taxon-specific DNA taxonomy. This study utilizes DNA data from nuclear and mitochondrial genes to improve the taxonomy of an important forest beetle pest, Ips. Amount and type of nucleotide difference are associated with monophyletic species based on a cladistic analysis of these data. As a result, a new species from China is described for a clade of beetles whose nucleotide differences exceeded the amount of evolutionary change observed within currently recognized species. The COI data are analyzed independently with an expanded taxon data set, including pair-wise nucleotide differences between recognized sister species. The wide range of average intraspecific pair-wise nucleotide difference (0-10.0%) suggests limitations to the application of a standard percent nucleotide difference as a means to identify species boundaries. At most, average COI nucleotide intraspecific difference provides an informal guide to identify potential clades that may warrant further systematic investigation.

Evaluating Neanderthal Genetics and Phylogeny

The retrieval of Neanderthal (Homo neanderthalsensis) mitochondrial DNA is thought to be among the most significant ancient DNA contributions to date, allowing conflicting hypotheses on modern human (Homo sapiens) evolution to be tested directly. Recently, however, both the authenticity of the Neanderthal sequences and their phylogenetic position outside contemporary human diversity have been questioned. Using Bayesian inference and the largest dataset to date, we find strong support for a monophyletic Neanderthal clade outside the diversity of contemporary humans, in agreement with the expectations of the Out-of-Africa replacement model of modern human origin. From average pairwise sequence differences, we obtain support for claims that the first published Neanderthal sequence may include errors due to postmortem damage in the template molecules for PCR. In contrast, we find that recent results implying that the Neanderthal sequences are products of PCR artifacts are not well supported, suffering from inadequate experimental design and a presumably high percentage (>68%) of chimeric sequences due to "jumping PCR" events.

Composition of archaeal, bacterial, and eukaryal RuBisCO genotypes in three Western Pacific arc hydrothermal vent systems

We studied the diversity of all forms of the RuBisCO large subunit-encoding gene cbbL in three RuBisCO uncharacterized hydrothermal vent communities. This diversity included the archaeal cbbL and the forms IC and ID, which have not previously been studied in the deep-sea environment, in addition to the forms IA, IB and II. Vent plume sites were Fryer and Pika in the Mariana arc and the Suiyo Seamount, Izu-Bonin, Japan. The cbbL forms were PCR amplified from plume bulk microbial DNA and then cloned and sequenced. Archaeal cbbL was detected in the Mariana samples only. Both forms IA and II were amplified from all samples, while the form IC was amplified only from the Pika and Suiyo samples. Only the Suiyo sample showed amplification of the form ID. The form IB was not recorded in any sample. Based on rarefaction analysis, nucleotide diversity and average pairwise difference, the archaeal cbbL was the most diverse form in Mariana samples, while the bacterial form IA was the most diverse form in the Suiyo sample. Also, the Pika sample harbored the highest diversity of cbbL phylogenetic lineages. Based on pairwise reciprocal library comparisons, the Fryer and Pika archaeal cbbL libraries showed the most significant difference, while Pika and Suiyo showed the highest similarity for forms IA and II libraries. This suggested that the Fryer supported the most divergent sequences. All archaeal cbbL sequences formed unique phylogenetic lineages within the branches of anaerobic thermophilic archaea of the genera Pyrococcus, Archaeoglobus, and Methanococcus. The other cbbL forms formed novel phylogenetic clusters distinct from any recorded previously in other deep-sea habitats. This is the first evidence for the diversity of archaeal cbbL in environmental samples.

Listeria monocytogenes internalins are highly diverse and evolved by recombination and positive selection

To probe the evolution of internalins with confirmed or suspected roles in Listeria monocytogenes virulence we sequenced the full inlB, inlC2, inlC, inlD, inlE, inlF, inlG, and inlH ORFs from 40 L. monocytogenes isolated from human ( n = 10) and animal ( n = 10) clinical cases, foods ( n = 10), and the natural environment ( n = 10). inlB and inlE were present in all isolates, representing 26 and 20 alleles, respectively. inlC was found in all lineage I and II isolates and represented 21 alleles. inlC2 and inlD represented 22 and 24 alleles, respectively, and were found in all L. monocytogenes isolates, with the exception of three lineage II isolates, which carried inlH, an apparent fusion of the 5′ end of inlC2 with the 3′ end of inlD. inlF and inlG were absent from lineage I isolates and represented 16 and 11 alleles, respectively. Average pairwise nucleotide differences per site ( π) ranged from 0.00849 ( inlF) to 0.07020 ( inlE). Phylogenetic trees generally showed clustering of internalin genes into two major evolutionary lineages consistent with lineages I and II previously assigned by ribotyping. In addition to detection of recombination events within each internalin gene, inlB, inlC, inlC2, and inlF showed significant evidence for positive selection (i.e., selection for an advantageous mutant allele). Overall, our data indicated that (i) internalin genes are highly diverse, (ii) internalin gene sequences cluster consistent with the phylogenetic lineages of L. monocytogenes, (iii) both intragenic recombination and positive selection have contributed to the evolution of L. monocytogenes internalins, and (iv) L. monocytogenes internalins show distinct evolutionary histories.

TheompAGene inChlamydia trachomatisDiffers in Phylogeny and Rate of Evolution from Other Regions of the Genome

Strains of Chlamydia trachomatis are classified into serovars based on nucleotide sequence differences in ompA, the gene that encodes the major outer membrane protein. Phylogenetic characterization of strains based on ompA, however, results in serovar groupings that are inconsistent with the distinguishing features of C. trachomatis pathobiology, e.g., tissue tropisms and disease presentation. We have compared nucleotide sequences at multiple sites distributed around the chlamydial genome from 18 strains representing 16 serovars; sampled regions included genes encoding housekeeping enzymes (totaling 2,073 bp), intergenic noncoding segments (1,612 bp), and a gene encoding a second outer membrane protein (porB; 1,023 bp), with the ompA sequence (1,194 bp) used for reference. These comparative analyses revealed substantial variation in nucleotide substitution patterns among the sampled regions, with average pairwise sequence differences ranging from 0.15% for the housekeeping genes to 12.1% for ompA. Phylogenetic characterization of the sampled genomic sequences yielded a strongly supported tree that divides the strains into groupings consistent with C. trachomatis biology and which has a topology quite distinct from the ompA tree. This phylogenetic incongruity can be accounted for by recombination of the ompA gene between different genomic backgrounds. We found, however, no evidence of recombination within or between any of the sampled regions around the C. trachomatis genome apart from ompA. Parallel analysis of published sequence data on four members of the pmp gene family are consistent with the phylogenetic analyses reported here.

Molecular diversity of local Norwegian meadow fescue (Festuca pratensis Huds.) populations and Nordic cultivars—consequences for management and utilisation

Genetic diversity and relatedness were studied in 30 Norwegian local populations of meadow fescue (Festuca pratensis Huds.) using amplified fragment length polymorphism (AFLP) markers. The populations were also compared with 13 Nordic meadow fescue cultivars in order to analyse the distribution of variation in local populations and cultivars and to elucidate relationships between local populations and cultivars. Analysis of molecular variance (AMOVA: ) analysis showed that most of the variation was present within populations and that little variation was found between local populations and cultivars. Separate AMOVA: analyses of local populations and cultivars revealed a higher level of variation within registered cultivars than within local populations. A cluster analysis based on corrected average pairwise differences between populations showed that the populations could be divided into three clusters, of which one also contained the cultivars. These results were supported by principal coordinates analysis. The results indicate that the Norwegian meadow fescue has a narrow genetic basis and that the local populations in Norway can be divided into three groups following the most probable routes of introduction of the species into Norway. The inland populations are closely related to the cultivars and have most probably been established as a result of migration from sown meadows. The western and southern populations probably originate from human activity-for example, trade-to the coastal western and northern parts of the country and to the central parts of southern Norway.

Molecular analysis of carbon monoxide-oxidizing bacteria associated with recent Hawaiian volcanic deposits.

Genomic DNA extracts from four sites at Kilauea Volcano were used as templates for PCR amplification of the large subunit (coxL) of aerobic carbon monoxide dehydrogenase. The sites included a 42-year-old tephra deposit, a 108-year-old lava flow, a 212-year-old partially vegetated ash-and-tephra deposit, and an approximately 300-year-old forest. PCR primers amplified coxL sequences from the OMP clade of CO oxidizers, which includes isolates such as Oligotropha carboxidovorans, Mycobacterium tuberculosis, and Pseudomonas thermocarboxydovorans. PCR products were used to create clone libraries that provide the first insights into the diversity and phylogenetic affiliations of CO oxidizers in situ. On the basis of phylogenetic and statistical analyses, clone libraries for each site were distinct. Although some clone sequences were similar to coxL sequences from known organisms, many sequences appeared to represent phylogenetic lineages not previously known to harbor CO oxidizers. On the basis of average nucleotide diversity and average pairwise difference, a forested site supported the most diverse CO-oxidizing populations, while an 1894 lava flow supported the least diverse populations. Neither parameter correlated with previous estimates of atmospheric CO uptake rates, but both parameters correlated positively with estimates of microbial biomass and respiration. Collectively, the results indicate that the CO oxidizer functional group associated with recent volcanic deposits of the remote Hawaiian Islands contains substantial and previously unsuspected diversity.

Genomic diversity of Erwinia carotovora subsp. carotovora and its correlation with virulence.

We used genetic and biochemical methods to examine the genomic diversity of the enterobacterial plant pathogen Erwinia carotovora subsp. carotovora. The results obtained with each method showed that E. carotovora subsp. carotovora strains isolated from one ecological niche, potato plants, are surprisingly diverse compared to related pathogens. A comparison of 23 partial mdh sequences revealed a maximum pairwise difference of 10.49% and an average pairwise difference of 2.13%, values which are much greater than the maximum variation (1.81%) and average variation (0.75%) previously reported for Escherichia coli. Pulsed-field gel electrophoresis analysis of I-CeuI-digested genomic DNA revealed seven rrn operons in all E. carotovora subsp. carotovora strains examined except strain WPP17, which had only six copies. We identified 26 I-CeuI restriction fragment length polymorphism patterns and observed significant polymorphism in fragment sizes ranging from 100 to 450 kb for all strains. We detected large plasmids in two strains, including the model strain E. carotovora subsp. carotovora 71. The two least virulent strains had an unusual chromosomal structure, suggesting that a particular pulsotype is correlated with virulence. To compare chromosomal organization of multiple enterobacterial genomes, several genes were mapped onto I-CeuI fragments. We identified portions of the genome that appear to be conserved across enterobacteria and portions that have undergone genome rearrangements. We found that the least virulent strain, WPP17, failed to oxidize cellobiose and was missing several hrp and hrc genes. The unexpected variability among isolates obtained from clonal hosts in one region and in one season suggests that factors other than the host plant, potato, drive the evolution of this common environmental bacterium and key plant pathogen.

Analysis of facultative lithotroph distribution and diversity on volcanic deposits by use of the large subunit of ribulose 1,5-bisphosphate carboxylase/oxygenase.

A 492- to 495-bp fragment of the gene coding for the large subunit of the form I ribulose 1,5-bisphosphate carboxylase/oxygenase (RubisCO) (rbcL) was amplified by PCR from facultatively lithotrophic aerobic CO-oxidizing bacteria, colorless and purple sulfide-oxidizing microbial mats, and genomic DNA extracts from tephra and ash deposits from Kilauea volcano, for which atmospheric CO and hydrogen have been previously documented as important substrates. PCR products from the mats and volcanic sites were used to construct rbcL clone libraries. Phylogenetic analyses showed that the rbcL sequences from all isolates clustered with form IC rbcL sequences derived from facultative lithotrophs. In contrast, the microbial mat clone sequences clustered with sequences from obligate lithotrophs representative of form IA rbcL. Clone sequences from volcanic sites fell within the form IC clade, suggesting that these sites were dominated by facultative lithotrophs, an observation consistent with biogeochemical patterns at the sites. Based on phylogenetic and statistical analyses, clone libraries differed significantly among volcanic sites, indicating that they support distinct lithotrophic assemblages. Although some of the clone sequences were similar to known rbcL sequences, most were novel. Based on nucleotide diversity and average pairwise difference, a forested site and an 1894 lava flow were found to support the most diverse and least diverse lithotrophic populations, respectively. These indices of diversity were not correlated with rates of atmospheric CO and hydrogen uptake but were correlated with estimates of respiration and microbial biomass.

Paternal Population History of East Asia: Sources, Patterns, and Microevolutionary Processes

Asia has served as a focal point for human migration during much of the Late Pleistocene and Holocene. Clarification of East Asia's role as a source and/or transit point for human dispersals requires that this region's own settlement history be understood. To this end, we examined variation at 52 polymorphic sites on the nonrecombining portion of the Y chromosome (NRY) in 1,383 unrelated males, representing 25 populations from southern East Asia (SEAS), northern East Asia (NEAS), and central Asia (CAS). The polymorphisms defined 45 global haplogroups, 28 of which were present in these three regions. Although heterozygosity levels were similar in all three regions, the average pairwise difference among haplogroups was noticeably smaller in SEAS. Multidimensional scaling analysis indicated a general separation of SEAS versus NEAS and CAS populations, and analysis of molecular variance produced very different values of Phi(ST) in NEAS and SEAS populations. In spatial autocorrelation analyses, the overall correlogram exhibited a clinal pattern; however, the NEAS populations showed evidence of both isolation by distance and ancient clines, whereas there was no evidence of structure in SEAS populations. Nested cladistic analysis demonstrated that population history events and ongoing demographic processes both contributed to the contrasting patterns of NRY variation in NEAS and SEAS. We conclude that the peopling of East Asia was more complex than earlier models had proposed-that is, a multilayered, multidirectional, and multidisciplinary framework is necessary. For instance, in addition to the previously recognized genetic and dental dispersal signals from SEAS to NEAS populations, CAS has made a significant contribution to the contemporary gene pool of NEAS, and the Sino-Tibetan expansion has left traces of a genetic trail from northern to southern China.