Recent spread of a Y-chromosomal lineage in northern China and Mongolia

We have identified a Y-chromosomal lineage that is unusually frequent in northeastern China and Mongolia, in which a haplotype cluster defined by 15 Y short tandem repeats was carried by approximately 3.3% of the males sampled from East Asia. The most recent common ancestor of this lineage lived 590 +/- 340 years ago (mean +/- SD), and it was detected in Mongolians and six Chinese minority populations. We suggest that the lineage was spread by Qing Dynasty (1644-1912) nobility, who were a privileged elite sharing patrilineal descent from Giocangga (died 1582), the grandfather of Manchu leader Nurhaci, and whose documented members formed approximately 0.4% of the minority population by the end of the dynasty.     

Population genetic analysis of the N-acylsphingosine amidohydrolase gene associated with mental activity in humans

To understand the evolution of human mental activity, we performed population genetic analyses of nucleotide sequences ( approximately 11 kb) from a worldwide sample of 60 chromosomes of the N-acylsphingosine amidohydrolase (ASAH1) gene. ASAH1 hydrolyzes ceramides and regulates neuronal development, and its deficiency often results in mental retardation. In the region ( approximately 4.4 kb) encompassing exons 3 and 4 of this gene, two distinct lineages (V and M) have been segregating in the human population for 2.4 +/- 0.4 million years (MY). The persistence of these two lineages is attributed to ancient population structure of humans in Africa. However, all haplotypes belonging to the V lineage exhibit strong linkage disequilibrium, a high frequency (62%), and small nucleotide diversity (pi = 0.05%). These features indicate a signature of positive Darwinian selection for the V lineage. Compared with the orthologs in mammals and birds, it is only Val at amino acid site 72 that is found exclusively in the V lineage in humans, suggesting that this Val is a likely target of positive selection. Computer simulation confirms that demographic models of modern humans except for the ancient population structure cannot explain the presence of two distinct lineages, and neutrality is incompatible with the observed small genetic variation of the V lineage at ASAH1. On the basis of the above observations, it is argued that positive selection is possibly operating on ASAH1 in the modern human population.     

Discordant Phylogeographic Patterns between the Y Chromosome and Mitochondrial DNA in the House Mouse: Selection on the Y Chromosome?

We have compared patterns of geographic variation and molecular divergence of mitochondrial DNA (mtDNA) and Y chromosome over the range of the different subspecies of Mus musculus. MtDNA was typed for 305 nucleotides in the control region, the Y chromosome for 834 base pairs (bp) in Zfy introns and 242 bp in Sry, a Zfy2 18-bp deletion, and two microsatellites. Apparent discrepancies exist between the distributions of the lineages of mtDNA and of the two major Y-chromosome lineages thus defined: some subspecies share the same mtDNA lineage but have different Y-chromosome lineages or vice versa. One microsatellite reveals a geographically clustered variation inside the distribution of each Y-chromosome lineage, showing that new Y-chromosome variants can rapidly spread locally. The two major Y-chromosome lineages have a divergence time only about one fourth of that between mtDNA lineages. Although this recent coalescence of the Y chromosomes between subspecies could partly be due to a lower ancestral polymorphism of the Y chromosome, it suggests that secondary introgression after the radiation of the subspecies might have occurred. There is evidence that the differentiation of the Y-chromosome lineages contributes to partial reproductive isolation between subspecies, and patterns of molecular evolution suggest that selection has played a role in the rapid spread across subspecies.     

A selective sweep driven by pyrimethamine treatment in southeast asian malaria parasites

Malaria parasites (Plasmodium falciparum) provide an excellent system in which to study the genomic effects of strong selection in a recombining eukaryote because the rapid spread of resistance to multiple drugs during the last the past 50 years has been well documented, the full genome sequence and a microsatellite map are now available, and haplotype data can be easily generated. We examined microsatellite variation around the dihydrofolate reductase (dhfr) gene on chromosome 4 of P. falciparum. Point mutations in dhfr are known to be responsible for resistance to the antimalarial drug pyrimethamine, and resistance to this drug has spread rapidly in Southeast (SE) Asia after its introduction in 1970s. We genotyped 33 microsatellite markers distributed across chromosome 4 in 61 parasites from a location on the Thailand/Myanmar border. We observed minimal microsatellite length variation in a 12-kb (0.7-cM) region flanking the dhfr gene and diminished variation for approximately 100 kb (6 cM), indicative of a single origin of resistant alleles. Furthermore, we found the same or similar microsatellite haplotypes flanked resistant dhfr alleles sampled from 11 parasite populations in five SE Asian countries indicating recent invasion of a single lineage of resistant dhfr alleles in locations 2000 km apart. Three features of these data are of especially interest. (1). Pyrimethamine resistance is generally assumed to have evolved multiple times because the genetic basis is simple and resistance can be selected easily in the laboratory. Yet our data clearly indicate a single origin of resistant dhfr alleles sampled over a large region of SE Asia. (2). The wide valley ( approximately 6 cM) of reduced variation around dhfr provides "proof-of-principle" that genome-wide association may be an effective way to locate genes under strong recent selection. (3). The width of the selective valley is consistent with predictions based on independent measures of recombination, mutation, and selection intensity, suggesting that we have reasonable estimates of these parameters. We conclude that scanning the malaria parasite genome for evidence of recent selection may prove an extremely effective way to locate genes underlying recently evolved traits such as drug resistance, as well as providing an opportunity to study the dynamics of selective events that have occurred recently or are currently in progress.     

Adaptive evolution of genes underlying schizophrenia

Schizophrenia poses an evolutionary-genetic paradox because it exhibits strongly negative fitness effects and high heritability, yet it persists at a prevalence of approximately 1% across all human cultures. Recent theory has proposed a resolution: that genetic liability to schizophrenia has evolved as a secondary consequence of selection for human cognitive traits. This hypothesis predicts that genes increasing the risk of this disorder have been subject to positive selection in the evolutionary history of humans and other primates. We evaluated this prediction using tests for recent selective sweeps in human populations and maximum-likelihood tests for selection during primate evolution. Significant evidence for positive selection was evident using one or both methods for 28 of 76 genes demonstrated to mediate liability to schizophrenia, including DISC1, DTNBP1 and NRG1, which exhibit especially strong and well-replicated functional and genetic links to this disorder. Strong evidence of non-neutral, accelerated evolution was found for DISC1, particularly for exon 2, the only coding region within the schizophrenia-associated haplotype. Additionally, genes associated with schizophrenia exhibited a statistically significant enrichment in their signals of positive selection in HapMap and PAML analyses of evolution along the human lineage, when compared with a control set of genes involved in neuronal activities. The selective forces underlying adaptive evolution of these genes remain largely unknown, but these findings provide convergent evidence consistent with the hypothesis that schizophrenia represents, in part, a maladaptive by-product of adaptive changes during human evolution.     

Genomic surveillance of SARS-CoV-2 tracks early interstate transmission of P.1 lineage and diversification within P.2 clade in Brazil

The sharp increase of COVID-19 cases in late 2020 has made Brazil the new epicenter of the ongoing SARS-CoV-2 pandemic. The novel viral lineages P.1 (Variant of Concern Gamma) and P.2, respectively identified in the Brazilian states of Amazonas and Rio de Janeiro, have been associated with potentially higher transmission rates and antibody neutralization escape. In this study, we performed the whole-genome sequencing of 185 samples isolated from three out of the five Brazilian regions, including Amazonas (North region), Rio Grande do Norte, Paraíba and Bahia (Northeast region), and Rio de Janeiro (Southeast region) in order to monitor the spread of SARS-CoV-2 lineages in Brazil in the first months of 2021. Here, we showed a widespread dispersal of P.1 and P.2 across Brazilian regions and, except for Amazonas, P.2 was the predominant lineage identified in the sampled states. We estimated the origin of P.2 lineage to have happened in February, 2020 and identified that it has differentiated into new clades. Interstate transmission of P.2 was detected since March, but reached its peak in December, 2020 and January, 2021. Transmission of P.1 was also high in December and its origin was inferred to have happened in August 2020. We also confirmed the presence of lineage P.7, recently described in the southernmost region of Brazil, to have spread across the Northeastern states. P.1, P.2 and P.7 are descended from the ancient B.1.1.28 strain, which co-dominated the first phase of the pandemic in Brazil with the B.1.1.33 strain. We also identified the occurrence of a new lineage descending from B.1.1.33 that convergently carries the E484K mutation, N.9. Indeed, the recurrent report of many novel SARS-CoV-2 genetic variants in Brazil could be due to the absence of effective control measures resulting in high SARS-CoV2 transmission rates. Altogether, our findings provided a landscape of the critical state of SARS-CoV-2 across Brazil and confirm the need to sustain continuous sequencing of the SARS-CoV-2 isolates worldwide in order to identify novel variants of interest and monitor for vaccine effectiveness.     

Evolution of the dentition in prehistoric Ohio Valley Native Americans: II. morphology of the deciduous dentition

In order to evaluate the microevolutionary dynamics of morphological features of the deciduous dentition, I collected data on the variation of 57 features (33 crown and 24 root) from prehistoric Ohio Valley populations. I sampled a total of 370 individuals from 26 populations representing a lineage that inhabited the middle and upper Ohio valley region from approximately 3000 to 350 BP. Evolutionary changes in the frequencies of morphological features of the deciduous teeth in this lineage were limited. Over 80% of the features show no significant differences among the populations. The relatively few features that show consistent differences separate pre- and postmaize agricultural populations. I discuss explanations for this change in terms of selection differences or gene flow. The general pattern of morphological trait expression in the deciduous teeth of this Ohio Valley lineage corresponds to what has been termed the Mongoloid dental complex (sinodonty in the permanent teeth). I suggest additional features that, with further study, may be added to this morphological complex. Am J Phys Anthropol 106:189–205, 1998. © 1998 Wiley-Liss, Inc.     

Emergence of a New Lineage of Dengue Virus Type 2 Identified in Travelers Entering Western Australia from Indonesia, 2010-2012

Dengue virus (DENV) transmission is ubiquitous throughout the tropics. More than 70% of the current global dengue disease burden is borne by people who live in the Asia-Pacific region. We sequenced the E gene of DENV isolated from travellers entering Western Australia between 2010–2012, most of whom visited Indonesia, and identified a diverse array of DENV1-4, including multiple co-circulating viral lineages. Most viruses were closely related to lineages known to have circulated in Indonesia for some time, indicating that this geographic region serves as a major hub for dengue genetic diversity. Most notably, we identified a new lineage of DENV-2 (Cosmopolitan genotype) that emerged in Bali in 2011–2012. The spread of this lineage should clearly be monitored. Surveillance of symptomatic returned travellers provides important and timely information on circulating DENV serotypes and genotypes, and can reveal the herald wave of dengue and other emerging infectious diseases.     

Natural Polymorphisms in Tap2 Influence Negative Selection and CD4∶CD8 Lineage Commitment in the Rat

Genetic variation in the major histocompatibility complex (MHC) affects CD4∶CD8 lineage commitment and MHC expression. However, the contribution of specific genes in this gene-dense region has not yet been resolved. Nor has it been established whether the same genes regulate MHC expression and T cell selection. Here, we assessed the impact of natural genetic variation on MHC expression and CD4∶CD8 lineage commitment using two genetic models in the rat. First, we mapped Quantitative Trait Loci (QTLs) associated with variation in MHC class I and II protein expression and the CD4∶CD8 T cell ratio in outbred Heterogeneous Stock rats. We identified 10 QTLs across the genome and found that QTLs for the individual traits colocalized within a region spanning the MHC. To identify the genes underlying these overlapping QTLs, we generated a large panel of MHC-recombinant congenic strains, and refined the QTLs to two adjacent intervals of ∼0.25 Mb in the MHC-I and II regions, respectively. An interaction between these intervals affected MHC class I expression as well as negative selection and lineage commitment of CD8 single-positive (SP) thymocytes. We mapped this effect to the transporter associated with antigen processing 2 (Tap2) in the MHC-II region and the classical MHC class I gene(s) (RT1-A) in the MHC-I region. This interaction was revealed by a recombination between RT1-A and Tap2, which occurred in 0.2% of the rats. Variants of Tap2 have previously been shown to influence the antigenicity of MHC class I molecules by altering the MHC class I ligandome. Our results show that a restricted peptide repertoire on MHC class I molecules leads to reduced negative selection of CD8SP cells. To our knowledge, this is the first study showing how a recombination between natural alleles of genes in the MHC influences lineage commitment of T cells.     

Origin and evolution of Japanese encephalitis virus in southeast Asia

Since it emerged in Japan in the 1870s, Japanese encephalitis has spread across Asia and has become the most important cause of epidemic encephalitis worldwide. Four genotypes of Japanese encephalitis virus (JEV) are presently recognized (representatives of genotypes I to III have been fully sequenced), but its origin is not known. We have determined the complete nucleotide and amino acid sequence of a genotype IV Indonesian isolate (JKT6468) which represents the oldest lineage, compared it with other fully sequenced genomes, and examined the geographical distribution of all known isolates. JKT6468 was the least similar, with nucleotide divergence ranging from 17.4 to 19.6% and amino acid divergence ranging from 4.7 to 6.5%. It included an unusual series of amino acids at the carboxy terminus of the core protein unlike that seen in other JEV strains. Three signature amino acids in the envelope protein (including E327 Leu-->Thr/Ser on the exposed lateral surface of the putative receptor binding domain) distinguished genotype IV strains from more recent genotypes. Analysis of all 290 JEV isolates for which sequence data are available showed that the Indonesia-Malaysia region has all genotypes of JEV circulating, whereas only more recent genotypes circulate in other areas (P < 0.0001). These results suggest that JEV originated from its ancestral virus in the Indonesia-Malaysia region and evolved there into the different genotypes which then spread across Asia. Our data, together with recent evidence on the origins of other emerging viruses, including dengue virus and Nipah virus, imply that tropical southeast Asia may be an important zone for emerging pathogens.     

Processes of diversification and dispersion of rice yellow mottle virus inferred from large-scale and high-resolution phylogeographical studies

Phylogeography of Rice yellow mottle virus (RYMV) was reconstructed from the coat protein gene sequences of a selection of 173 isolates from the 14 countries of mainland Africa where the disease occurred and from the full sequences of 16 representative isolates. Genetic variation was linked to geographical distribution and not to host species as isolates from wild rice always clustered with isolates from cultivated rice of the same region. Genetic variation was not associated to agro-ecology, viral interference and insect vector species. Distinct RYMV lineages occurred in East, Central and West Africa, although the Central African lineage included isolates from Benin, Togo and Niger at the west, adjacent to countries of the West African lineage. Genetic subdivision at finer geographical scales was apparent within lineages of Central and West Africa, although less pronounced than in East Africa. Physical obstacles, but also habitat fragmentation, as exemplified by the small low-lying island of Pemba offshore Tanzania mainland, explained strain localization. Three new highly divergent strains were found in eastern Tanzania. By contrast, intensive surveys in Cote d'Ivoire and Guinea at the west of Africa did not reveal any new variant. Altogether, this supported the view that the Eastern Arc Mountains biodiversity hotspot was the centre of origin of RYMV and that the virus spread subsequently from east to west across Africa. In West Africa, specific strains occurred in the Inner Niger Delta and suggested it was a secondary centre of diversification. Processes for diversification and dispersion of RYMV are proposed.     

The making of a rapid plant invader: genetic diversity and differentiation in the native and invaded range of Senecio inaequidens

To become invasive, exotic species have to succeed in the consecutive phases of introduction, naturalization, and invasion. Each of these phases leaves traces in genetic structure, which may affect the species’ success in subsequent phases. We examined this interplay of genetic structure and invasion dynamics in the South African Ragwort (Senecio inaequidens), one of Europe’s fastest plant invaders. We used AFLP and microsatellite markers to analyze 19 native African and 32 invasive European populations. In combination with historic data, we distinguished invasion routes and traced them back to the native source areas. This revealed that different introduction sites had markedly different success in the three invasion phases. Notably, an observed lag‐phase in Northern Germany was evidently not terminated by factors increasing the invasiveness of the resident population but by invasive spread from another introduction centre. The lineage invading Central Europe was introduced to sites in which winters are more benign than in the native source region. Subsequently, this lineage spread into areas in which winter temperatures match the native climate more closely. Genetic diversity clearly increases with population age in Europe and less clearly decreases with spread rate up to population establishment. This indicates that gene flow along well‐connected invasion routes counteracted losses of genetic diversity during rapid spread. In summary, this study suggests that multiple introductions, environmental preadaptation and high gene flow along invasion routes contributed to the success of this rapid invader. More generally, it demonstrates the benefit of combining genetic, historical, and climatic data for understanding biological invasions.     

Evolution in natural and experimental populations of Amphicarpaea bracteata

Genetically divergent lineages (biotypes) often coexist within local populations of the selfing annual Amphicarpaea bracteata (Leguminosae). At one site, a rapid change in lineage frequency was observed, with one biotype decreasing from 41% to 16% of the population over a 2 year period. During 6 additional years of censuses, this lineage showed no tendency to recover to its former abundance. An experimental population was constructed with approximately equal numbers of individuals from the two dominant biotypes in this population, and plants were allowed to grow and reproduce without interference for 3 generations. The same lineage that declined in the natural habitat was also strongly selected against in this experimental population, with a selection coefficient (0.57 per generation) comparable to that inferred from the initial decline in frequency of this lineage in nature (0.47). A greenhouse experiment suggested that differential adaptation to high vs. low light environments could be one factor affecting relative success of these lineages. Differences in reproductive allocation and disease resistance to local pathogens are also key factors likely to affect relative fitness. Overall, these results show that the genetic structure of selfing plant populations can be highly dynamic, with intense selection and sudden shifts in relative abundance of different lineages.     

Multiple cryptic species with divergent substrate affinities in the Serpula himantioides species complex

Serpula himantioides is a widespread saprotrophic morphospecies mainly colonising coniferous wood in nature, but it appears frequently in buildings as well. From an earlier study, it is known that at least three divergent lineages occur within the S. himantioides species complex. In this study, a broader sample of S. himantioides isolates has been analysed by multi-locus sequencing, including new isolates from Asia, North and South America. Altogether five phylogenetical species (PS1-5) were detected, all recognised across independent gene phylogenies. A new southern South American phylogenetic species (PS1) was found, representing an early diverging lineage within the S. himantioides species complex. The two closely related PS2 and PS3 lineages included isolates from North America only, and PS4 was also dominated by North American isolates. Most of the investigated isolates (76%) clustered into PS5, a lineage that has been found on most continents, including North America. Overall, little phylogeographical structure was found in PS5, indicating frequent and recent long-distance dispersal events within this widespread lineage. Our analyses indicate that South and North America are the centres of divergence for the S.?himantioides species complex. Some of the lineages seem adapted to various substrates, but PS5 is able to decay a wide array of angiosperms and gymnosperms, which may have facilitated the spread of this lineage throughout the world.     

Genetic evidence for a high diversity and wide distribution of endemic strains of the pathogenic chytrid fungus Batrachochytrium dendrobatidis in wild Asian amphibians

Population declines and extinctions of amphibians have been attributed to the chytrid fungus Batrachochytrium dendrobatidis (Bd), especially one globally emerging recombinant lineage (‘Bd‐GPL’). We used PCR assays that target the ribosomal internal transcribed spacer region (ITS) of Bd to determine the prevalence and genetic diversity of Bd in South Korea, where Bd is widely distributed but is not known to cause morbidity or mortality in wild populations. We isolated Korean Bd strains from native amphibians with low infection loads and compared them to known worldwide Bd strains using 19 polymorphic SNP and microsatellite loci. Bd prevalence ranged between 12.5 and 48.0%, in 11 of 17 native Korean species, and 24.7% in the introduced bullfrog Lithobates catesbeianus. Based on ITS sequence variation, 47 of the 50 identified Korean haplotypes formed a group closely associated with a native Brazilian Bd lineage, separated from the Bd‐GPL lineage. However, multilocus genotyping of three Korean Bd isolates revealed strong divergence from both Bd‐GPL and the native Brazilian Bd lineages. Thus, the ITS region resolves genotypes that diverge from Bd‐GPL but otherwise generates ambiguous phylogenies. Our results point to the presence of highly diversified endemic strains of Bd across Asian amphibian species. The rarity of Bd‐GPL‐associated haplotypes suggests that either this lineage was introduced into Korea only recently or Bd‐GPL has been outcompeted by native Bd strains. Our results highlight the need to consider possible complex interactions among native Bd lineages, Bd‐GPL and their associated amphibian hosts when assessing the spread and impact of Bd‐GPL on worldwide amphibian populations.     

Clonal analysis of vertebrate myogenesis. 3. Developmental changes in the muscle-colony-forming cells of the human fetal limb

The cell lineage of developing human limb muscle has been investigated by means of an in vitro clonal assay. Single cells capable of forming differentiated muscle colonies have been detected within the prospective leg muscle region as early as the 36th day of human development (Streeter's Horizon XVI). At this stage clonable myoblasts account for 14% of the total colony-forming cells. The relative proportion of clonable myoblasts increases rapidly during subsequent fetal development and attains a plateau level of approximately 90% by the 100th day of development. The 90% plateau level persists at least until day 172. By correlating the percent muscle colony differentiation with clonal plating efficiency and with the number of single cells derived from the total limb muscle region of fetuses of different ages, an estimate of the actual number of clonable myoblasts within the developing limb musculature is obtained.Sequential changes within the muscle cell lineage have been further dissected by the temporal analysis of age-dependent medium effects on muscle colony differentiation. The analysis indicates that clonable myoblasts derived from early fetuses are sensitive to medium conditions to which older muscle-colony-forming cells are relatively insensitive. In addition, fusing and nonfusing colonies have been classified into recognizable morphological types whose relative proportions are observed to change during limb development. The results are correlated with human limb morphogenesis and skeletal muscle histogenesis and an operational model of muscle cell lineage is proposed.     

Evolution of the phosphoglycerate mutase processed gene in human and chimpanzee revealing the origin of a new primate gene

Processed genes are created by retroposition from messenger RNA of expressed genes. The estimated amount of processed copies of genes in the human genome is 10,000-14,000. Some of these might be pseudogenes with the expected pattern for nonfunctional sequences, but some others might be an important source of new genes. We have studied the evolution of a Phosphoglycerate mutase processed gene (PGAM3) described in humans and believed to be a pseudogene. We sequenced PGAM3 in chimpanzee and macaque and obtained polymorphism data for human coding region. We found evidence that PGAM3 likely produces a functional protein, as an example of addressing functionality for human processed pseudogenes. First, the open reading frame was intact despite many deletions that occurred in the 3' untranslated region. Second, it appears that the gene is expressed. Finally, interspecies and intraspecies variation for PGAM3 was not consistent with the neutral model proposed for pseudogenes, suggesting that a new functional primate gene has originated. Amino acid divergence was significantly higher than synonymous divergence in PGAM3 lineage, supporting positive selection acting in this gene. This role of selection was further supported by the excess of rare alleles in a population genetic analysis. PGAM3 is located in a region of very low recombination; therefore, it is conceivable that the rapid fixation events in this newly arising gene may have contributed to a selective sweep of variation in the region.     

Intron positions delineate the evolutionary path of a pervasively appended peptide in five human aminoacyl-tRNA synthetases

Recent progress in genome sequencing has revealed a correspondence between the evolution of multicellularity and the appending of new peptides onto age-old enzyme bodies. Indicative of the pervasive nature of these appended peptides, in some cases the same sequences have been appended to a number of different enzymes. By analyzing the positions of introns within one such roaming peptide, an approximately 50-amino acid motif appended to five human aminoacyl-tRNA synthetases, I have delineated its path in eukaryote evolution. The motif was first acquired as an N-terminal extension by histidyl- and glycyl-tRNA synthetases at a very early stage of eukaryote evolution. Later, but not less than 1200 million years ago, the motif spread from histidyl-tRNA synthetase to the C and N terminals of glutamyl- and prolyl-tRNA synthetase, respectively, and then spread further during the evolution of the Chordate lineage to the N terminal of tryptophanyl-tRNA synthetase. In similar fashion, the motif in glycyl-tRNA synthetase spread to the C terminal of methionyl-tRNA synthetase not later than 1000 million years ago.     

Phenotypic traits of the Mediterranean <Emphasis Type="Italic">Phragmites australis</Emphasis> M1 lineage: differences between the native and introduced ranges

The environmental conditions in the new ranges of introduced plant species are often different from the conditions in their native ranges, and invasive plant species have been assumed to adapt to different environmental conditions by rapid ecological evolution in the invasive range after the introduction. Another interpretation of the change in plant traits after their introduction, however, is ecological fitting, which is based on the inherently high phenotypic plasticity of the species rather than on evolution. The Mediterranean haplotype M1 lineage of the wetland grass Phragmites australis was introduced to the coastal wetlands along the Gulf Coast of North America, where it is exposed to a different climate compared to its original range. The climate in the native range is arid or temperate with dry and hot summers, whereas the climate in the introduced range is warmer and has a higher and more uniform precipitation than that in the native range. This warmer and more humid environment is likely to pose different selection pressures to the plants in the introduced range and thus cause rapid evolutionary change and phenotypic differentiation in the introduced range. Here, we compared phenotypic traits of the M1 lineage from the native and introduced ranges in a common garden experiment to study the processes assisting the successful spread in the introduced range. Overall, the native and introduced groups were similar, but we detected a few phenotypic traits that diverged. Ecological fitting could be the fundamental mechanism by which the P. australis M1 lineage survives and spreads in the introduced Gulf Coast region. However, further research is needed to assess how the diverging traits observed in our study in Denmark (lower photosynthetic rates, lower chlorophylls concentration and higher leaf K concentration for the introduced than for the native genotypes) are expressed in the two ranges.     

Detecting substitution-rate heterogeneity among regions of a nucleotide sequence

Likelihood-ratio statistics are proposed to test for heterogeneity in nucleotide substitution rate among regions of a DNA sequence. The tests examine three-sequence phylogenies, and two specific tests are proposed: a test to detect rate heterogeneity among genic regions within a sequence, over all evolutionary lineages; and a test to detect rate heterogeneity among regions in a specific evolutionary lineage. Simulations examine the ability of tests to detect a single region that varies in nucleotide substitution rate relative to the remainder of the sequence. A 50-bp region with a fivefold substitution-rate increase can be detected > or = 90% of the time when it is found in all three lineages of the phylogeny, and a 50-bp region of fivefold rate increase can be detected with approximately 70% power when it is found in only one evolutionary lineage. Simulation also examines the effect of transition- and transversion-rate differences. The tests are applied to published DNA sequences. While the tests are powerful, significant results can be difficult to interpret biologically.