Patterns of sequence divergence and evolution of the S orthologous regions between Asian and African cultivated rice species

A strong postzygotic reproductive barrier separates the recently diverged Asian and African cultivated rice species, Oryza sativa and O. glaberrima. Recently a model of genetic incompatibilities between three adjacent loci: S(1)A, S(1) and S(1)B (called together the S(1) regions) interacting epistatically, was postulated to cause the allelic elimination of female gametes in interspecific hybrids. Two candidate factors for the S(1) locus (including a putative F-box gene) were proposed, but candidates for S(1)A and S(1)B remained undetermined. Here, to better understand the basis of the evolution of regions involved in reproductive isolation, we studied the genic and structural changes accumulated in the S(1) regions between orthologous sequences. First, we established an 813 kb genomic sequence in O. glaberrima, covering completely the S(1)A, S(1) and the majority of the S(1)B regions, and compared it with the orthologous regions of O. sativa. An overall strong structural conservation was observed, with the exception of three isolated regions of disturbed collinearity: (1) a local invasion of transposable elements around a putative F-box gene within S(1), (2) the multiple duplication and subsequent divergence of the same F-box gene within S(1)A, (3) an interspecific chromosomal inversion in S(1)B, which restricts recombination in our O. sativa×O. glaberrima crosses. Beside these few structural variations, a uniform conservative pattern of coding sequence divergence was found all along the S(1) regions. Hence, the S(1) regions have undergone no drastic variation in their recent divergence and evolution between O. sativa and O. glaberrima, suggesting that a small accumulation of genic changes, following a Bateson-Dobzhansky-Muller (BDM) model, might be involved in the establishment of the sterility barrier. In this context, genetic incompatibilities involving the duplicated F-box genes as putative candidates, and a possible strengthening step involving the chromosomal inversion might participate to the reproductive barrier between Asian and African rice species.     

中国栽培稻的起源与演化

中国栽培稻起源地主要有4种假说:华南起源说、云南起源说、长江中下游说和长江中游-淮河上游说。对其进行了简要概括,并进行评述:目前没有一种假说能完全符合水稻起源的4个条件,此外有些学者提出水稻的多起源中心;从分子生物学角度综述了栽培稻籼粳分化研究进展,主要有两个起源说:一源论和二次(或多次)起源论,大量研究支持籼稻和粳稻多起源论;从古DNA、水稻落粒性与驯化、植硅石等方面探讨栽培稻起源演化的研究方法。旨在能为水稻的起源演化研究提供有价值的参考。     

Elemental Content in Brown Rice by Inductively Coupled Plasma Atomic Emission Spectroscopy Reveals the Evolution of Asian Cultivated Rice

The phylogenetic relationship for classification traits and eight mineral elements in brown rice （Oryza sativa L.） from Yunnan Province in China was carried out using microwave assisted digestion followed by inductively coupled plasma atomic emission spectroscopy, and the analytical procedures were carefully controlled and validated. In general, the results show that the mean levels of K, Ca, Mg, Fe and Cu in brown rice for 789 accessions of rice landraces was distinctly lower than that of improved cultivars. They further demonstrate that Ca plays an important role in the differentiation of subspecies indica-japonica, especially to enhance adaptation of cold stress, and that five mineral elements in brown rice enhance the eurytopicity from landrace to improved cultivar. Hierarchical cluster analysis, using average linkage from SPSS software based on eight mineral elements in brown rice, showed that Yunnan rice could be grouped into rice landrace and improved cultivar, with the rice landrace being further clustered into five subgroups, and that, interestingly, purple rice does not cluster with either of the groups. Our present data confirm that indica is the closest relative of late rice and white rice, and that they constitute rice landraces together, whereas japonica is the closest relatives of non-nuda, early-mid and glutinous rice. It is further shown that japonica, non-nuda, early-mid, glutinous, white and red rice might be more primitive than indica, nuda, late, non-glutinous and purple rice, respectively.     

Elemental Content in Brown Rice by Inductively Coupled Plasma Atomic Emission Spectroscopy Reveals the Evolution of Asian Cultivated Rice

The phylogenetic relationship for classification traits and eight mineral elements in brown rice (Oryza sativa L.) from Yunnan Province in China was carried out using microwave assisted digestion followed by inductively coupled plasma atomic emission spectroscopy, and the analytical procedures were carefully controlled and validated. In general, the results show that the mean levels of K, Ca, Mg, Fe and Cu in brown rice for 789 accessions of rice landraces was distinctly lower than that of improved cultivars. They further demonstrate that Ca plays an important role in the differentiation of subspecies indica-japonica, especially to enhance adaptation of cold stress, and that five mineral elements in brown rice enhance the eurytopicity from landrace to improved cultivar. Hierarchical cluster analysis, using average linkage from SPSS software based on eight mineral elements in brown rice, showed that Yunnan rice could be grouped into rice landrace and improved cultivar, with the rice landrace being further clustered into five subgroups, and that, interestingly, purple rice does not cluster with either of the groups. Our present data confirm that indica is the closest relative of late rice and white rice, and that they constitute rice landraces together, whereas japonica is the closest relatives of non-nuda, early-mid and glutinous rice. It is further shown that japonica, non-nuda, early-mid, glutinous, white and red rice might be more primitive than indica, nuda, late, non-glutinous and purple rice, respectively.     

基于线粒体基因片段核苷酸多态性的亚洲栽培稻起源进化研究

亚洲栽培稻的祖先是普通野生稻,已成为世界公认的观点,然而亚洲栽培稻的2个亚种:粳稻和籼稻是一次起源还是二次起源仍存在很大争议,其起源地是国内还是国外依然是国际学者间争论的焦点。本文通过对184份亚洲栽培稻和203份普通野生稻3段基因序列cox3、cox1、orf 224和2段基因间序列ssv-39/178、rps2-trnfM的多样性研究,验证了以下观点:1)粳稻起源于中国,籼稻起源于中国和国外;2)亚洲栽培稻的起源为二次起源,即普通野生稻存在偏籼和偏粳2种类型,亚洲栽培稻的2个亚种籼稻和粳稻在进化过程中分别由偏籼型的普通野生稻和偏粳型的普通野生稻进化而来。     

Evolution of Orthologous Intronless and Intron-Bearing Globin Genes in Two Insect Species

While globin genes ctt-2β and ctt-9.1 in Chironomus thummi thummi each have a single intron, all of the other insect globin genes reported so far are intronless. We analyzed four globin genes linked to the two intron-bearing genes in C. th. thummi. Three have a single intron at the same position as ctt-2β and ctt-9.1; the fourth is intronless and lies between intron bearing genes. Finally, in addition to its intron, one gene (ctt-13RT) was recently interrupted by retrotransposition. Phylogenetic analyses show that the six genes in C. th. thummi share common ancestry with five globin genes in the distantly related species C. tentans, and that a 5-gene ancestral cluster predates the divergence of the two species. One gene in the ancestral cluster gave rise to ctn-ORFB in C. tentans, and duplicated in C. th. thummi to create ctt-11 and ctt-12. From parsimonious calculations of evolutionary distances since speciation, ctt-11, ctt-12, and ctn-ORFB evolved rapidly, while ctn-ORFE in C. tentans evolved slowly compared to other globin genes in the clusters. While these four globins are under selective pressure, we suggest that most chironomid globin genes were not selected for their unique function. Instead, we propose that high gene copy number itself was selected because conditions favored organisms that could synthesize more hemoglobin. High gene copy number selection to produce more of a useful product may be the basis of forming multigene families, all of whose members initially accumulate neutral substitutions while retaining essential function. Maintenance of a large family of globin genes not only ensured high levels of hemoglobin production, but may have facilitated the extensive divergence of chironomids into as many as 5000 species. Received: 31 December 1996 / Accepted: 16 May 1997     

亚洲栽培稻起源、演化中两个重要稻种类型的研究

分布子云南及南亚的光稃稻和镰刀谷在研究亚洲栽培稻(O.sativa L.)的起源、演化及稻作育种上有重要意义。它们在分类上的地位往往被误分,进而影响这些类型的开发利用。研究表明,光稃稻是普通野生稻被引上亚洲的热带与亚热带山区在旱地上首先被人类驯化成的原始粳型陆稻,当其进一步被引向高海拔和高纬度时则演化成普通粳稻。镰刀谷迄今多被划为籼稻,我们的初步研究表明,镶刀谷大部分属粳稻,其中也包括一部分籼稻。它的起源及其分类上的地位尚待研究。     

Comparative Sequence Analysis of the Ghd7 Orthologous Regions Revealed Movement of Ghd7 in the Grass Genomes

Ghd7 is an important rice gene that has a major effect on several agronomic traits, including yield. To reveal the origin of Ghd7 and sequence evolution of this locus, we performed a comparative sequence analysis of the Ghd7 orthologous regions from ten diploid Oryza species, Brachypodium distachyon, sorghum and maize. Sequence analysis demonstrated high gene collinearity across the genus Oryza and a disruption of collinearity among non-Oryza species. In particular, Ghd7 was not present in orthologous positions except in Oryza species. The Ghd7 regions were found to have low gene densities and high contents of repetitive elements, and that the sizes of orthologous regions varied tremendously. The large transposable element contents resulted in a high frequency of pseudogenization and gene movement events surrounding the Ghd7 loci. Annotation information and cytological experiments have indicated that Ghd7 is a heterochromatic gene. Ghd7 orthologs were identified in B. distachyon, sorghum and maize by phylogenetic analysis; however, the positions of orthologous genes differed dramatically as a consequence of gene movements in grasses. Rather, we identified sequence remnants of gene movement of Ghd7 mediated by illegitimate recombination in the B. distachyon genome.     

Patterns and Processes of Genome-Wide Divergence Between North American and African Drosophila melanogaster

Genomic tools and analyses are now being widely used to understand genome-wide patterns and processes associated with speciation and adaptation. In this article, we apply a genomics approach to the model organism Drosophila melanogaster. This species originated in Africa and subsequently spread and adapted to temperate environments of Eurasia and the New World, leading some populations to evolve reproductive isolation, especially between cosmopolitan and Zimbabwean populations. We used tiling arrays to identify highly differentiated regions within and between North America (the United States and Caribbean) and Africa (Cameroon and Zimbabwe) across 63% of the D. melanogaster genome and then sequenced representative fragments to study their genetic divergence. Consistent with previous findings, our results showed that most differentiation was between populations living in Africa vs. outside of Africa (i.e., “out-of-Africa” divergence), with all other geographic differences being less substantial (e.g., between cosmopolitan and Zimbabwean races). The X chromosome was much more strongly differentiated than the autosomes between North American and African populations (i.e., greater X divergence). Overall differentiation was positively associated with recombination rates across chromosomes, with a sharp reduction in regions near centromeres. Fragments surrounding these high F_ST sites showed reduced haplotype diversity and increased frequency of rare and derived alleles in North American populations compared to African populations. Nevertheless, despite sharp deviation from neutrality in North American strains, a small set of bottleneck/expansion demographic models was consistent with patterns of variation at the majority of our high F_ST fragments. Although North American populations were more genetically variable compared to Europe, our simulation results were generally consistent with those previously based on European samples. These findings support the hypothesis that most differentiation between North America and Africa was likely driven by the sorting of African standing genetic variation into the New World via Europe. Finally, a few exceptional loci were identified, highlighting the need to use an appropriate demographic null model to identify possible cases of selective sweeps in species with complex demographic histories.THE study of genetic differentiation between populations and species has recently been empowered by the use of genomic techniques and analysis (e.g., Noor and Feder 2006; Stinchcombe and Hoekstra 2008). In the past decade, genetic studies of adaptation and speciation have taken advantage of emerging molecular techniques to scan the genomes of diverging populations for highly differentiated genetic regions (e.g., Wilding et al. 2001; Emelianov et al. 2003; Beaumont and Balding 2004; Campbell and Bernatchez 2004; Scotti-Saintagne et al. 2004; Achere et al. 2005; Turner et al. 2005; Vasemagi et al. 2005; Bonin et al. 2006, 2007; Murray and Hare 2006; Savolainen et al. 2006; Yatabe et al. 2007; Nosil et al. 2008, 2009; Turner et al. 2008a,b; Kulathinal et al. 2009). As a result, genome scans can identify candidate regions that may be associated with adaptive evolution between diverging populations and, more broadly, are able to describe genome-wide patterns and processes of population differentiation (Begun et al. 2007; Stinchcombe and Hoekstra 2008).Genome scans in well-studied genetic model species such as Drosophila melanogaster gain particular power because differentiated loci are mapped to a well-annotated genome. Moreover, the evolutionary history of D. melanogaster is rich with adaptive and demographic events with many parallels to human evolution. Most notable is the historical out-of-Africa migration and subsequent adaptation to temperate ecological environments of Europe, Asia, North America, and Australia. This has resulted in widespread genetic and phenotypic divergence between African and non-African populations (e.g., David and Capy 1988; Begun and Aquadro 1993; Capy et al. 1994; Colegrave et al. 2000; Rouault et al. 2001; Takahashi et al. 2001; Caracristi and Schlötterer 2003; Baudry et al. 2004; Pool and Aquadro 2006; Schmidt et al. 2008; Yukilevich and True 2008a,b). Further, certain populations in Africa and in the Caribbean vary in their degree of reproductive isolation from populations in more temperate regions (Wu et al. 1995; Hollocher et al. 1997; Yukilevich and True 2008a,b). In particular, the Zimbabwe and nearby populations of southern Africa are strongly sexually isolated from all other populations, designating them as a distinct behavioral race (Wu et al. 1995).D. melanogaster has received a great deal of attention from the population geneticists in studying patterns of sequence variation across African and non-African populations. Many snapshots have been taken of random microsatellite and SNP variants spread across X and autosomes, and these have generated several important conclusions. Polymorphism patterns in European populations are characterized by reduced levels of nucleotide and haplotype diversity, an excess of high frequency-derived polymorphisms, and elevated levels of linkage disequilibrium relative to African populations (e.g., Begun and Aquadro 1993; Andolfatto 2001; Glinka et al. 2003; Haddrill et al. 2005; Ometto et al. 2005; Thornton and Andolfatto 2006; Hutter et al. 2007; Singh et al. 2007). These results have been generally interpreted as compatible with population size reduction/bottlenecks followed by recent population expansions. On the other hand, African populations are generally assumed either to have been relatively constant in size over time or to have experienced population size expansions. They generally show higher levels of nucleotide and haplotype diversity, an excess of rare variants, and a deficit of high frequency-derived alleles (Glinka et al. 2003; Ometto et al. 2005; Pool and Aquadro 2006; Hutter et al. 2007; but see Haddrill et al. 2005 for evidence of bottlenecks in Africa).Previous work also shows that the ratio of X-linked to autosomal polymorphism deviates from neutral expectations in opposite directions in African and European populations with more variation on the X than expected in Africa and less variation on the X than expected in Europe (Andolfatto 2001; Kauer et al. 2002; Hutter et al. 2007; Singh et al. 2007). The deviation from neutrality in the ratio of X-autosome polymorphism may be explained by positive selection being more prevalent on the X in Europe and/or by a combination of bottlenecks and male-biased sex ratios in Europe and female-biased sex ratios in Africa (Charlesworth 2001; Hutter et al. 2007; Singh et al. 2007). The selective explanation stems from the argument that, under the hitchhiking selection model, X-linked loci are likely to be more affected by selective sweeps than autosomal loci (Maynard Smith and Haigh 1974; Charlesworth et al. 1987; Vicoso and Charlesworth 2006, 2009).The relative contribution of selective and demographic processes in shaping patterns of genomic variation and differentiation is highly debated (Wall et al. 2002; Glinka et al. 2003; Haddrill et al. 2005; Ometto et al. 2005; Schöfl and Schlötterer 2004; Thornton and Andolfatto 2006; Hutter et al. 2007; Singh et al. 2007; Shapiro et al. 2007; Stephan and Li 2007; Hahn 2008; Macpherson et al. 2008; Noor and Bennett 2009; Sella et al. 2009). This is especially the case in D. melanogaster because derived non-African populations have likely experienced a complex set of demographic events during their migration out of Africa (e.g., Thornton and Andolfatto 2006; Singh et al. 2007; Stephan and Li 2007), making population genetics signatures of demography and selection difficult to tease apart (e.g., Macpherson et al. 2008). Thus it is still unclear what role selection has played in shaping overall patterns of genomic variation and differentiation relative to demographic processes in this species.While there is a long tradition in studying arbitrarily or opportunistically chosen sequences in D. melanogaster, genomic scans that focus particularly on highly differentiated sites across the genome have received much less attention. Such sites are arguably the best candidates to resolve the debate on which processes have shaped genomic differentiation within species (e.g., Przeworski 2002). Recently, a genome-wide scan of cosmopolitan populations in the United States and in Australia was performed to investigate clinal genomic differentiation on the two continents (Turner et al. 2008a). Many single feature polymorphisms differentiating Northern and Southern Hemisphere populations were identified. Among the most differentiated loci in common between continents, 80% were differentiated in the same orientation relative to the Equator, implicating selection as the likely explanation (Turner et al. 2008a). Larger regions of genomic differentiation within and between African and non-African populations have also been discovered, some of them possibly being driven by divergent selection (e.g., Dopman and Hartl 2007; Emerson et al. 2008; Turner et al. 2008a, Aguade 2009). Despite this recent progress, we still know relatively little about large-scale patterns of genomic differentiation in this species, especially between African and non-African populations, and whether most of this differentiation is consistent with demographic processes alone or if it requires selective explanations.In this work, we explicitly focus on identifying differentiated sites across the genome between U.S., Caribbean, West African, and Zimbabwean populations. This allows us to address several fundamental questions related to genomic evolution in D. melanogaster, such as the following: (1) Do genome-wide patterns of differentiation reflect patterns of reproductive isolation? (2) Is genomic differentiation random across and within chromosomes or are some regions overrepresented? (3) What are the population genetics properties of differentiated sites and their surrounding sequences? (4) Can demographic historical processes alone explain most of the observed differentiation on a genome-wide level or is it necessary to involve selection in their explanation?In general, our findings revealed that most genomic differentiation within D. melanogaster shows an out-of-Africa genetic signature. These results are inconsistent with the notion that most genomic differentiation occurs between cosmopolitan and Zimbabwean reproductively isolated races. Further, we found that the X is more differentiated between North American and African populations and more strongly deviates from pure neutrality in North American populations relative to autosomes. Nevertheless, our article shows that much of this deviation from neutrality is broadly consistent with several demographic null models, with a few notable exceptions. Athough this does not exclude selection as a possible alternative mechanism for the observed patterns, it supports the idea that most differentiation in D. melanogaster was likely driven by the sorting of African standing genetic variation into the New World.     

Construction of Pseudomolecule Sequences of the aus Rice Cultivar Kasalath for Comparative Genomics of Asian Cultivated Rice

Having a deep genetic structure evolved during its domestication and adaptation, the Asian cultivated rice (Oryza sativa) displays considerable physiological and morphological variations. Here, we describe deep whole-genome sequencing of the aus rice cultivar Kasalath by using the advanced next-generation sequencing (NGS) technologies to gain a better understanding of the sequence and structural changes among highly differentiated cultivars. The de novo assembled Kasalath sequences represented 91.1% (330.55 Mb) of the genome and contained 35 139 expressed loci annotated by RNA-Seq analysis. We detected 2 787 250 single-nucleotide polymorphisms (SNPs) and 7393 large insertion/deletion (indel) sites (>100 bp) between Kasalath and Nipponbare, and 2 216 251 SNPs and 3780 large indels between Kasalath and 93-11. Extensive comparison of the gene contents among these cultivars revealed similar rates of gene gain and loss. We detected at least 7.39 Mb of inserted sequences and 40.75 Mb of unmapped sequences in the Kasalath genome in comparison with the Nipponbare reference genome. Mapping of the publicly available NGS short reads from 50 rice accessions proved the necessity and the value of using the Kasalath whole-genome sequence as an additional reference to capture the sequence polymorphisms that cannot be discovered by using the Nipponbare sequence alone.     

Parallel Genetic Evolution Within and Between Bacteriophage Species of Varying Degrees of Divergence

Parallel evolution is the acquisition of identical adaptive traits in independently evolving populations. Understanding whether the genetic changes underlying adaptation to a common selective environment are parallel within and between species is interesting because it sheds light on the degree of evolutionary constraints. If parallel evolution is perfect, then the implication is that forces such as functional constraints, epistasis, and pleiotropy play an important role in shaping the outcomes of adaptive evolution. In addition, population genetic theory predicts that the probability of parallel evolution will decline with an increase in the number of adaptive solutions—if a single adaptive solution exists, then parallel evolution will be observed among highly divergent species. For this reason, it is predicted that close relatives—which likely overlap more in the details of their adaptive solutions—will show more parallel evolution. By adapting three related bacteriophage species to a novel environment we find (1) a high rate of parallel genetic evolution at orthologous nucleotide and amino acid residues within species, (2) parallel beneficial mutations do not occur in a common order in which they fix or appear in an evolving population, (3) low rates of parallel evolution and convergent evolution between species, and (4) the probability of parallel and convergent evolution between species is strongly effected by divergence.     

Duplication and Divergence of Two Distinct Pancreatic Ribonuclease Genes in Leaf-Eating African and Asian Colobine Monkeys

doi: 10.1093/molbev/msl025 John     

Molecular Phylogeny,Evolution, and Functional Divergence of the LSD1-Like Gene Family: Inference from the Rice Genome

The identification of LSD1-like genes in parasite, green algae, moss, pine, and monocot and dicot species allowed us to trace the phylogenetic history of this gene family. Computational analysis showed that the diversification of members of this family could be dated back to the early stage of plant evolution. The evolution of plant LSD1-like genes was possibly shaped by two duplication events. These proteins, which contain three copies of the LSD1 zinc finger (zf-LSD1) domain within their entire polypeptides and play crucial roles in modulating disease defense and cell death, resulted from the second duplication. A gain of zf-LSD1 domain model was reasonable for explaining the origination of three-zf-LSD1 domain-containing proteins. The zf-LSD1 domain phylogeny showed that the middle (M) and C-terminal (C) domains originated from a common ancestor; the N-terminal (N) domain might be more ancient than the former two. The divergence of the N, M, and C domains was well before the monocot-dicot split. Coevolution analysis revealed that four intramolecular domain pairs, including the N domain and the interregion between the M and the C domains (INTER2), the M and C domain, the N- and C-terminus, and the M domain and C-terminus, possibly coevolved during the evolution of three-zf-LSD1 domain-containing proteins. The three zf-LSD1 domains are evolutionary conserved. Thus, the differences at the N- and C-terminus would be crucial for functional specificity of LSD1 genes. Strong functional constraints should work on the zf-LSD1 domains, whereas reduced functional constraint was found in the INTER2 region. Functional divergence analysis showed that three-zf-LSD1 domain-containing proteins were significantly functionally divergent from those proteins containing only one zf-LSD1 domain, a result demonstrating that shifted evolutionary rates between the two clusters were significantly different from each other. [Reviewing Editor: Dr. Joshua Plotkin]     

Evolutionary Processes Acting on Candidate cis-Regulatory Regions in Humans Inferred from Patterns of Polymorphism and Divergence

Analysis of polymorphism and divergence in the non-coding portion of the human genome yields crucial information about factors driving the evolution of gene regulation. Candidate cis-regulatory regions spanning more than 15,000 genes in 15 African Americans and 20 European Americans were re-sequenced and aligned to the chimpanzee genome in order to identify potentially functional polymorphism and to characterize and quantify departures from neutral evolution. Distortions of the site frequency spectra suggest a general pattern of selective constraint on conserved non-coding sites in the flanking regions of genes (CNCs). Moreover, there is an excess of fixed differences that cannot be explained by a Gamma model of deleterious fitness effects, suggesting the presence of positive selection on CNCs. Extensions of the McDonald-Kreitman test identified candidate cis-regulatory regions with high probabilities of positive and negative selection near many known human genes, the biological characteristics of which exhibit genome-wide trends that differ from patterns observed in protein-coding regions. Notably, there is a higher probability of positive selection in candidate cis-regulatory regions near genes expressed in the fetal brain, suggesting that a larger portion of adaptive regulatory changes has occurred in genes expressed during brain development. Overall we find that natural selection has played an important role in the evolution of candidate cis-regulatory regions throughout hominid evolution.     

Estimating the Rate of Adaptive Molecular Evolution When the Evolutionary Divergence Between Species is Small

We investigate the extent by which the estimates of the rate of adaptive molecular evolution obtained by extending the McDonald-Kreitman test are biased if the species, subjected to analysis, diverged recently. We show that estimates can be biased if the nucleotide divergence between the species is low relative to within species variation, and that the magnitude of the bias depends on the rate of adaptive evolution and the distribution of fitness effects of new mutations. Bias appears to be because of three factors: (1) misattribution of polymorphism to divergence; (2) the contribution of ancestral polymorphism to divergence; and (3) different rates of fixation of neutral and advantageous mutations. If there is little adaptive molecular evolution, then slightly deleterious mutations inflate estimates of the rate of adaptive evolution, because these contribute proportionately more to polymorphism than to nucleotide divergence than neutral mutations. However, if there is substantial adaptive evolution, polymorphism contributing to apparent divergence may downwardly bias estimates. We propose a simple method for correcting the different contributions of slightly deleterious and neutral mutations to polymorphism and divergence, and apply it to datasets from several species. We find that estimates of the rate of adaptive molecular evolution from closely related species may be underestimates by ~10% or more. However, after the contribution of polymorphism to divergence is removed, the rate of adaptive evolution may still be overestimated as a consequence of ancestral polymorphism and time for fixation effects. This bias may be substantial if branch lengths are less than 10N (e) generations.     

Affinities between Asian non-human Schistosoma species,the S. indicum group,and the African human schistosomes

Schistosoma species have traditionally been arranged in groups based on egg morphology, geographical origins, and the genus or family of snail intermediate host. One of these groups is the 'S. indicum group' comprising species from Asia that use pulmonate snails as intermediate hosts. DNA sequences were obtained from the four members of this group (S. indicum, S. spindale, S. nasale and S. incognitum) to provide information concerning their phylogenetic relationships with other Asian and African species and species groups. The sequences came from the second internal transcribed spacer (ITS2) of the ribosomal gene repeat, part of the 28S ribosomal RNA gene (28S), and part of the mitochondrial cytochrome c oxidase subunit 1 (CO1) gene. Tree analyses using both distance and parsimony methods showed the S. indicum group not to be monophyletic. Schistosoma indicum, S. spindale and S. nasale were clustered among African schistosomes, while S. incognitum was placed as sister to the African species (using ITS2 and 28S nucleotide sequences and CO1 amino acid sequences), or as sister to all other species of Schistosoma (CO1 nucleotide sequences). Based on the present molecular data, a scenario for the evolution of the S. indicum group is discussed.     

The Nature of Nucleotide Sequence Divergence between Barley and Maize Chloroplast DNA

Analysis of a 2175-base pair (bp) SmaI-HindIII fragment of barley chloroplast DNA revealed that rbcL (the gene for the large subunit of ribulose 1,5-bisphosphate carboxylase) and atpB (the gene for the beta subunit of ATPase) are transcribed divergently and are separated by an untranscribed region of 155-166 bp. The rbcL mRNA has a 320-residue untranslated leader region, whereas the atpB mRNA has a 296- to 309-residue leader region. The sequence of these regions, together with the initial 113 bp of the atpB-coding region and the initial 1279 bp of the rbcL-coding region, is compared with the analogous maize chloroplast DNA sequences. Two classes of nucleotide differences are present, substitutions and insertions/deletions. Nucleotide substitutions show a 1.9-fold bias toward transitions in the rbcL-coding region and a 1.5-fold bias toward transitions in the noncoding region. The level of nucleotide substitutions between the barley and maize sequences is about 0.065/bp. Seventy-one percent of the substitutions in the rbcL-coding region are at the third codon position, and 95% of these are synonymous changes. Insertion/deletion events, which are confined to the noncoding regions, are not randomly distributed in these regions and are often associated with short repeated sequences. The extent of change for the noncoding regions (about 0.093 events/bp) is less than the extent of change at the third codon positions in the rbcL-coding region (about 0.135 events/bp), including insertion/delection events. Limited sequence analysis of the analogous DNA from a wild line ( Hordeum spontaneum) and a primitive Iranian barley (H. vulgare) suggested a low rate of chloroplast DNA evolution. Compared to spinach chloroplast DNA, the barley rbcL-atpB untranslated region is extremely diverged, with only the putative rbcL promoters and ribosome-binding site being extensively conserved.     

Sequence Divergence of Microsatellites and Phylogeny Analysis in Tetraploid Cotton Species and Their Putative Diploid Ancestors

To determine the level of microsatellite sequence differences and to use the information to construct a phylogenetic relationship for cultivated tetraploid cotton （Gossypium spp.） species and their putative diploid ancestors, 10 genome-derived microsatellite primer pairs were used to amplify eight species, including two tetraploid and six diploid species, in Gossypium. A total of 92 unique amplicons were resolved using polyacrylamide gel electrophoresis. Each amplicon was cloned, sequenced, and analyzed using standard phylogenetic software. Allelic diversities were caused mostly by changes in the number of simple sequence repeat （SSR） motif repeats and only a small proportion resulted from interruption of the SSR motif within the locus for the same genome. The frequency of base substitutions was 0.5%-1.0% in different genomes, with only few indels found. Based on the combined 10 SSR flanking sequence data, the homology of A-genome diploid species averaged 98.9%, even though most of the amplicons were of the same size, and the sequence homology between G. gossypioides （Ulbr.） Standl. and three other D-genome species （G. raimondii Ulbr., G. davidsonii Kell., and G. thurberi Tod.） was 98.5%, 98.6%, and 98.5%, respectively. Phylogenetic trees of the two allotetraploid species and their putative diploid progenitors showed that homoelogous sequences from the A- and D-subgenome were still present in the polyploid subgenomes and they evolved independently. Meanwhile, homoelogous sequence interaction that duplicated loci in the polyploid subgenomes became phylogenetic sisters was also found in the evolutionary history of tetraploid cotton species. The results of the present study suggest that evaluation of SSR variation at the sequence level can be effective in exploring the evolutionary relationships among Gossypuim species.     

Kin selection and the Evolution of Mutualisms between Species

Hamilton's theory of kin selection has revolutionized and inspired fifty years of additional theories and experiments on social evolution. Whereas Hamilton's broader intent was to explain the evolutionary stability of cooperation, his focus on shared genetic history appears to have limited the application of his theory to populations within a single species rather than across interacting species. The evolutionary mechanisms for cooperation between species require both spatial and temporal correlations among interacting partners for the benefits to be not only predictable but of sufficient duration to be reliably delivered. As a consequence when the benefits returned by mutualistic partners are redirected to individuals other than the original donor, cooperation usually becomes unstable and parasitism may evolve. However, theoretically, such redirection of mutualistic benefits may actually reinforce, rather than undermine, mutualisms between species when the recipients of these redirected benefits are genetically related to the original donor. Here, I review the few mathematical models that have used Hamilton's theory of kin selection to predict the evolution of mutualisms between species. I go on to examine the applicability of these models to the most well‐studied case of mutualism, pollinating seed predators, where the role of kin selection may have been previously overlooked. Future detailed studies of the direct, and indirect, benefits of mutualism are likely to reveal additional possibilities for applying Hamilton's theory of kin selection to mutualisms between species.