Neutrality tests using DNA polymorphism from multiple samples

The polymorphism of a gene or a locus is studied with increasing frequency by multiple laboratories or the same group at different times. Such practice results in polymorphism being revealed by different samples at different regions of the locus. Tests of neutrality have been widely conducted for polymorphism data but commonly used statistical tests cannot be applied directly to such data. This article provides a procedure to conduct a neutrality test and details are given for two commonly used tests. Applying the two new tests to the chemokine-receptor gene (CCR5) in humans, we found that the hypothesis that all mutations are selectively neutral cannot explain the observed pattern of DNA polymorphism.

Molecular evolution of CXCR1, a G protein-coupled receptor involved in signal transduction of neutrophils

Human neutrophils are a type of white blood cell, which forms an early line of defense against bacterial infections. Neutrophils are highly responsive to the chemokine, interleukin-8 (IL-8) due to the abundant distribution of CXCR1, one of the IL-8 receptors on the neutrophil cell surface. As a member of the GPCR family, CXCR1 plays a crucial role in the IL-8 signal transduction pathway in neutrophils. We sequenced the complete coding region of the CXCR1 gene in worldwide human populations and five representative nonhuman primate species. Our results indicate accelerated protein evolution in the human lineage, which was likely caused by Darwinian positive selection. The sliding window analysis and the codon-based neutrality test identified signatures of positive selection at the N-terminal ligand/receptor recognition domain of human CXCR1.

Intra- and interspecific variation of the CCR5 gene in higher primates

We have evaluated the molecular evolution of the chemokine receptor CCR5 in primates. The chemokine receptor CCR5 serves as a major co-receptor for human immunodeficiency virus/simian immunodeficiency virus (HIV/SIV) infection. Knowledge of evolution of the CCR5 molecule and selection on the CCR5 gene may shed light on its functional role. The comparison of differences between intraspecific polymorphisms and interspecific fixed substitutions provides useful information regarding modes of selection during the course of evolution. There is marked polymorphism in the CCR5 gene sequence within different primate species, whereas sequence divergence between different species is small. By using contingency tests, we compared synonymous (SS) and nonsynonymous (NS) CCR5 mutations occurring within and between a broad range of primates. Our results demonstrate that CCR5 evolution did not follow expectations, of strict neutrality at the level of the whole gene. The proportion of NS to SS at the intraspecific level was significantly higher than that observed at the interspecific level. These results suggest that most CCR5 NS polymorphisms are slightly deleterious. However, at domains more closely correlated with its known biological functions, there was no obvious evidence to support deviation from neutrality.

Evidence for Positive Selection on the Osteogenin (BMP3) Gene in Human Populations

Background: Human skeletal system has evolved rapidly since the dispersal of modern humans from Africa, potentially driven by selection and adaptation. Osteogenin (BMP3) plays an important role in skeletal development and bone osteogenesis as an antagonist of the osteogenic bone morphogenetic proteins, and negatively regulates bone mineral density. Methodology/Principal Findings: Here, we resequenced the BMP3 gene from individuals in four geographically separated modern human populations. Features supportive of positive selection in the BMP3 gene were found including the presence of an excess of nonsynonymous mutations in modern humans, and a significantly lower genetic diversity that deviates from neutrality. The prevalent haplotypes of the first exon region in Europeans demonstrated features of long-range haplotype homogeneity. In contrast with findings in European, the derived allele SNP Arg192Gln shows higher extended haplotype homozygosity in East Asian. The worldwide allele frequency distribution of SNP shows not only a high-derived allele frequency in Asians, but also in Americans, which is suggestive of functional adaptation. Conclusions/Significance: In conclusion, we provide evidence for recent positive selection operating upon a crucial gene in skeletal development, which may provide new insight into the evolution of the skeletal system and bone development.

Speciation and Phylogeography of Opsariichthys bidens (Pisces: Cypriniformes: Cyprinidae) in China: Analysis of the Cytochrome b Gene of mtDNA from Diverse Populations

Gao-Yan Li, Xu-Zhen Wang, Ya-Hui Zhao, Jie Zhang, Chun-Guang Zhang, and Shun-Ping He (2009) Speciation and phylogeography of Opsariichthys bidens (Pisces: Cypriniformes: Cyprinidae) in China: analysis of the cytochrome b gene of mtDNA from diverse populations. Zoological Studies 48(4): 569-583. The cyprinid fish Opsariichthys bidens Gunther is distributed in all major river systems of continental East Asia, and represents an attractive model for phylogeographic studies among cyprinid species or within a given species. In this study, we investigated the phylogeographic and demographic history of this species, using partial sequences of the cytochrome (cyt) b gene in mitochondrial (mt)DNA. Fish samples were collected from almost all major river systems where O. bidens is distributed in China. Sequence analysis showed remarkably high polymorphism, with 125 haplotypes in the 234 specimens examined, and with 89.8% of haplotypes occurring in only 1 specimen. A neutrality test indicated that some groups were not at mutation-drift equilibrium, suggesting a past population expansion. These results were supported by a mismatch distribution analysis. Based on our analysis, O. bidens consists of 4 groups belonging to 2 clades. The divergence time of the 2 clades was estimated to be 11.06-8.04 my. This value corresponds to the time of the 2nd uplift of the Qinghai-Tibet Plateau, the emergence of the East Asian monsoon, and the Epoch-6 Event. A two species scheme is proposed. http://zoolstud.sinica.edu.tw/Journals/48.4/569.pdf

香港地区湍蛙遗传多样性研究

依据线粒体上ND2和CO1两个变异较大的基因序列分析了香港地区香港湍蛙7种群、华南湍蛙1种群，以及大陆其他地区华南湍蛙7种群，戴云湍蛙1种群，武夷湍蛙1种群的系统发育关系，进而探讨香港湍蛙的遗传多样性、香港湍蛙特有性、如何确定香港湍蛙最佳保护单元以及这四种湍蛙的物种分类地位。 1． 香港湍蛙保护遗传学研究 香港湍蛙核苷酸传多样性较低，从其遗传多样性信息、单倍型网络分析、中性检验值以及岐点分布结果一致显示香港湍蛙很可能经历了瓶颈后的扩张，种群正在由一个较小的有效种群大小迅速增长, 有足够的时间通过变异用于积累单倍型的多态性, 而对于提高核苷酸多样化而言, 时间尚短（Nei M et al，1975，Avise J C，2000；李明等，2003）。 分子变异分析结果显示香港湍蛙种群间存在较多的基因交流，且系统发育树上各种群间交叉在一起，没有形成与地理单元相关的分支，而从其单倍型网络看，他们源于共同的祖先，是一个单系群，与地理单元间没有形成显著的遗传分化。因此应作为一个进化显著单元（ESU）。结合其与其他湍蛙发育关系及遗传距离以及野外采集信息认为香港湍蛙只在香港地区有分布，属于香港特有种。该物种内遗传多样性较低，又属于世界自然保护联盟红皮书中的近危种，同时也是《野生动物保护条例》中的受保护野生动物，且由于香港城市建设等使得其栖息环境受到威胁，因此在香港特别行政区应该受到重点保护。 从单倍型分布和核苷酸多样性可以看出大榄涌种群和城门种群具有较高的单倍型多样性和核苷酸多样性，应该作为保护的重点区域。 2. 华南湍蛙东、南沿海种群间系统关系 华南湍蛙分布广，各种群存在着丰富的遗传多样性信息且中部种群广西龙胜和湖南张家界种群核苷酸多样性明显高于其他边缘种群华南湍蛙。种群间几乎没有基因交流，且各种群间无共享单倍型，可见已形成了显著的遗传分化。各种群间遗传距离都较远，其中广东南昆山种群以及福建三港种群与其他种群距离最远，因此可以推测其他种群（广东深圳、香港大屿山、广西龙胜和防城以及湖南张家界种群）可能为独立进化的种群。但是否是一新种或一隐存种，还需要结合形态学进行更深入的研究。 本研究中无论从系统关系看还是从遗传距离看，大屿山种群与深圳种群最近，支持陈坚峰等将其定为华南湍蛙，即华南湍蛙新增一个分布点：香港大屿山。 系统树上广西防城种群（支B）与龙胜和湖南种群（支A）形成姐妹群。香港大屿山种群与深圳种群先形成姐妹群（支C），但却没有与其距离很近的广东南岭及南昆山种群（支D）形成姐妹群，可能粤北和粤中的环境及气候较复杂因此与粤南其他种群形成了明显的隔离。同时可以看出华南湍蛙种群遗传分化与地理距离没有显著的相关性。 3. 四种湍蛙间的系统关系 根据线粒体CO1基因建立四种湍蛙间的系统关系及其遗传距离，很清楚地看到，香港湍蛙与戴云湍蛙关系很近，而华南湍蛙则与武夷湍蛙较近。然而，戴云湍蛙同一个种群内部共有两个单倍型DY1和DY2，且两个单倍型间遗传距离大于DY1与香港湍蛙间遗传距离，更远远大于香港湍蛙种群内部的距离，即戴云湍蛙内部两个单倍型间遗传距离达到了种级水平，同样在系统发育树上这两个单倍型与香港湍蛙形成并系。但是，戴云湍蛙种内在形态上差异不显著。因此，其是否属于萌芽物种分化形成（budding speciation）或已经完全分化为两个不同的种值得进一步研究？ 与戴云湍蛙香港湍蛙关系类似，从系统树上看华南湍蛙不形成单系，而是分成两个大支，与武夷湍蛙形成并系，且福建和南昆山的华南湍蛙与武夷湍蛙遗传距离远大于武夷湍蛙种内福建种群与浙江种群的遗传距离，达到了种级分化水平。由此，可以推断武夷湍蛙是有效种。系统树上广东深圳、香港大屿山、广西防城和龙胜以及湖南张家界种群与华南湍蛙福建及南昆山各种群间遗传距离已超出了种内各种群间的遗传距离，但是至于这一支是否应为另外一个种，有必要扩大采样，并结合核基因及形态信息进行进一步研究。 MtDNA of ND2 and CO1 gene were used to investigate genetic diversity of Amolops in Hongkong .We collected seven populations of A. hongkongensis,，one population of A.ricketti from Hong Kong and other seven populations of A.ricketti from East and South of Chinese mainland. As well as one population of A. daiyunensis and one population of A.wuyiensis Phylogenetic relationship were analyzed of four species. Discussed whether A.hongkongensis is an endemic species and how can we make the conservation and management decisions. 1. Conservation Genetics of A. hongkongensis A. hongkongensis has a low nucleotide diversity, the results of genetic diversity, haplotype network, neutrality test and the mismatch distributions indicate that A. hongkongensis experienced a recent expansion after a bottle neck. They had enough time to accumulated haplotype diversity, but it’s too short to have a high nucleotide diversity（Nei M et al，1975，Avise J C，2000；Li et al，2003）. The result of AMOVA reveals that it has much gene exchange among the populations of A. hongkongensis. The clades of the phylogenetic tree were mixed together, no significant genetic differentiation among 8 populations and they share the same ancestor from the network analysis, these indicate that they are monophyly and should be protected as one ESU. Combined with the information of relationships of interspecies, genetic distance and distribution investigate, We conclude that A. hongkongensis is an endemic species of Hong Kong. Considering on the status of low genetic diversity in A.hongkongensis, and this species was listed in the IUCN red list as near threatened, as well as listed in the . Furthermore, it’s habitat loss and degradation more rapidly as the human activity got higher and higher. So it’s urgent to protect them in Hong Kong. Our results suggest that Tai Lam Wu and TAI MO Shan -Shing Mun populations have the higher priority to be protected because their higher genetic diversity. 2．Phylogenetic relationships among populations of Amolops ricketti from the Southern and eastern China A. ricketti has the considerable genetic diversity of mitochondrial haplotypes within and among populations, and Mitochondrial DNA diversity was higher in populations at the central area of the present distribution range of the frog，i. e. the Longsheng population and Zhangjiajie population, than at the edges of their distribution range. They have no share haplotype among populations, and have a significant genetic differentiation. Genetic distance is high among the populations, especially the distance of Nankun and Sangang group with others, which suggested that they evolved independently. May be there is a cryptic species or a new species, a further study is needed. The results of gene tree and the genetic distance clearly demonstrate that the population from LanTau island is A. ricketti, so we agree with Chen et al(2005) . That means A.ricketti have a new distribution site: LanTau island, HongKong. Phylogenetic relationships were analyzed through NJ and Mrbayes methods and got a consistent topological structure, the structure indicated that the ingroup were comprised four groups. Populations Longsheng and Zhangjiajie were first clustered as clade A; Populations Fangcheng was clustered together (clade B) as a sister group to clade A；Populations Shenzhen and Lantau island were sister groups and clustered as clade C；Then the clade D included populations Nankunshan and Nanling in Guangdong province and Sangang in Fujian province. 3. Phylogenetic Relationships among these four specises Phylogenetic relationships based on 1503bp CO1 gene and the genetic distance show that A. hongkongensis close to A. daiyunensis whereas A.ricketti near to A.wuyiensis. Nevertheless, there are two haplotypes in A.daiyunensis and the genetic distance between them higher than the distance between DY1 with A. hongkongensis. A. hongkongensis is nested in the paraphyletic ancestral species A. daiyunensis. Without significant difference in the morphological characters, So, we considered both A.daiyunensis and A.hongkongensis are valid species, may be this represents a case of ‘budding speciation’ like Batrachuperus pinchonii(Fu and Zeng,2008) in the population of A. daiyunensis. Just like two species above A. wuyiensis and A. ricketti are not monophyly, instead, A.wuyiensis is nested in the paraphyletic ancestral species A.ricketti. We need do more research to make sure whether they are new species.

Competition of lamellar orientation in thin films of a symmetric poly(styrene)-b-poly(L-lactide)diblock copolymer in melt state

We have studied the lamellar orientation in thin films of a model diblock copolymer, symmetric poly(styrene)-b-poly(L-lactide) (PS-PLLA), in the melt state on supported silicon wafer surface. In this system, while the PLLA block prefers to wet the polymer/substrate interface, the polymer/air as well as polymer/polymer interface is neutral for both blocks due to the similar surface energies of PS and PLLA in melt state. Our results demonstrate that the interplay of the interfaces during phase separation results in a series of structures before approaching the equilibrium state. Lamellar orientation of thin films with different initial film thicknesses at different annealing stages has been investigated using atomic force microscopy (AFM), transmission electronic microscopy (TEM), and X-ray photoelectron spectroscopy (XPS). It is found that in the early stage (annealing time t < 10 min), the polymer/substrate interface dominates the structure evolution, leading to a parallel lamellar structure with holes or islands formed depending on the initial film thickness. Later on, the neutral air interface becomes important and leads to a transition of lamellar orientation from parallel to perpendicular. It is interesting to see that for films with thickness h > 2L, where L is the bulk lamellar period, the lamellar orientation transition can occur independently in different parallel lamellar domains due to the neutrality of polymer/polymer interface.

Study of cupric hexacyanoferrate-modified platinum electrodes using probe beam deflection and electrochemical quartz crystal microbalance techniques

The ion exchange mechanism accompanying the oxidation/reduction processes of cupric hexacyanoferrate-modified platinum electrodes in different aqueous electrolyte solutions has been studied by means of in situ probe beam deflection and the electrochemical quartz crystal microbalance technique. The results demonstrate that the charge neutrality of the film during the reoxidation/reduction process is accomplished predominantly by the movement of cations, but anions and/or solvent are also participator(s). Moreover, in KHC8H4O4 (potassium biphthalate) solution, the EQCM data obtained from chronoamperometry experiment are more complicated than those in KCl and K2SO4 solutions. (C) 1997 Elsevier Science Ltd.

Population genetic structure of crimson snapper Lutjanus erythropterus in East Asia, revealed by analysis of the mitochondrial control region

The population genetic structure of the crimson snapper Lutjanus erythropterus in East Asia was examined with a 427-bp hypervariable portion of the mtDNA control region. A total of 262 samples were collected and 75 haplotypes were obtained. Neutrality tests (Tajima's and Fu's) suggested that Lutjanus erythropterus in East Asia had experienced a bottleneck followed by population expansion since the late Pleistocene. Despite the low phylogeographic structures in mtDNA haplotypes, a hierarchical examination of populations in 11 localities from four geographical regions using analysis of molecular variance (AMOVA) indicated significant genetic differentiation among regions (Phi(CT) = 0.08564, p < 0.01). Limited gene flow between the eastern region (including a locality in the western Pacific Ocean and two localities in the East Sea) and three geographic regions of the South China Sea largely contributed to the genetic subdivision. However, comparisons among three geographic regions of the South China Sea showed little to no genetic difference. Populations of Lutjanus erythropterus in East Asia are inferred to be divided into two major groups: an eastern group, including populations of the western Pacific Ocean and the East Sea, and a South China Sea group, consisting of populations from northern Malaysia to South China. The results suggest that fishery management should reflect the genetic differentiation and diversity in East Asia. (c) 2006 International Council for the Exploration of the Sea. Published by Elsevier Ltd. All rights reserved.

Comparative molecular population genetics of phycoerythrin locus in Prochlorococcus

As the only remainder type of phycobiliproteins in Prochlorococcus, the actual role of phycoerythrin still remains unknown. Previous studies revealed that two different forms of phycoerythrin gene were found in two ecotypes of Prochlorococcus that are specifically adapted to either high light (HL) or low light (LL) conditions. Here we analyze patterns of phycoerythrin nucleotide variation in the HL- and LL-Prochlorococcus populations. Our analyses reveal a significantly greater number of non-synonymous fixed substitutions in peB and peA than expected based on interspecific comparisons. This pattern of excess non-synonymous fixed substitutions is not seen in other five phycoerythrin-related genes (peZ/V/Y/T/S). Several neutrality statistical tests indicate an excess of rare frequency polymorphisms in the LL-Prochlorococcus data, but an excess of intermediate frequency polymorphisms in the HL-Prochlorococcus data. Distributions of the positively selected sites identified using the likelihood ratio test, when mapped onto the phycoerythrin tertiary structure, reveal that HL- and LL-phycoerythrin should be under different selective patterns. These findings may provide insights into the likely role of selection at the phycoerythrin locus and motivate further research to unveil the function of phycoerythrin in Prochlorococcus.