Evolutionary Relationship of DNA Sequences in Finite Populations

With the aim of analyzing and interpreting data on DNA polymorphism obtained by DNA sequencing or restriction enzyme technique, a mathematical theory on the expected evolutionary relationship among DNA sequences (nucleons) sampled is developed under the assumption that the evolutionary change of nucleons is determined solely by mutation and random genetic drift. The statistical property of the number of nucleotide differences between randomly chosen nucleons and that of heterozygosity or nucleon diversity is investigated using this theory. These studies indicate that the estimates of the average number of nucleotide differences and nucleon diversity have a large variance, and a large part of this variance is due to stochastic factors. Therefore, increasing sample size does not help reduce the variance significantly. The distribution of sample allele (nucleomorph) frequencies is also studied, and it is shown that a small number of samples are sufficient in order to know the distribution pattern.     

A Simple Method for Estimating Average Number of Nucleotide Substitutions within and between Populations from Restriction Data

A simple method is proposed for estimating the average number of nucleotide substitutions per site within and between populations for the case where a large number of individuals are examined for many restriction enzymes. This method gives essentially the same results as those obtained by Nei and Li's method but saves a large amount of computer time. The variances of the quantities estimated can be obtained by the jackknife method, and these variances are very similar to those obtained by Nei and Jin's more sophisticated method. A similar method can also be applied to DNA sequence data.     

Small‐ and large‐scale heterogeneity in genetic variation across the collard flycatcher genome: implications for estimating genetic diversity in nonmodel organisms

Population genetic studies in nonmodel organisms are often hampered by a lack of reference genomes that are essential for whole‐genome resequencing. In the light of this, genotyping methods have been developed to effectively eliminate the need for a reference genome, such as genotyping by sequencing or restriction site‐associated DNA sequencing (RAD‐seq). However, what remains relatively poorly studied is how accurately these methods capture both average and variation in genetic diversity across an organism's genome. In this issue of Molecular Ecology Resources, Dutoit et al. (2016) use whole‐genome resequencing data from the collard flycatcher to assess what factors drive heterogeneity in nucleotide diversity across the genome. Using these data, they then simulate how well different sequencing designs, including RAD sequencing, could capture most of the variation in genetic diversity. They conclude that for evolutionary and conservation‐related studies focused on the estimating genomic diversity, researchers should emphasize the number of loci analysed over the number of individuals sequenced.     

A Cladistic Analysis of Phenotypic Associations with Haplotypes Inferred from Restriction Endonuclease Mapping and DNA Sequence Data. III. Cladogram Estimation

We previously developed a cladistic approach to identify subsets of haplotypes defined by restriction endonuclease mapping or DNA sequencing that are associated with significant phenotypic deviations. Our approach was limited to segments of DNA in which little recombination occurs. In such cases, a cladogram can be constructed from the restriction site or sequence data that represents the evolutionary steps that interrelate the observed haplotypes. The cladogram is used to define a nested statistical design to identify mutational steps associated with significant phenotypic deviations. The central assumption behind this strategy is that any undetected mutation causing a phenotypic effect is embedded within the same evolutionary history that is represented by the cladogram. The power of this approach depends upon the confidence one has in the particular cladogram used to draw inferences. In this paper, we present a strategy for estimating the set of cladograms that are consistent with a particular sample of either restriction site or nucleotide sequence data and that includes the possibility of recombination. We first evaluate the limits of parsimony in constructing cladograms. Once these limits have been determined, we construct the set of parsimonious and nonparsimonious cladograms that is consistent with these limits. Our estimation procedure also identifies haplotypes that are candidates for being products of recombination. If recombination is extensive, our algorithm subdivides the DNA region into two or more subsections, each having little or no internal recombination. We apply this estimation procedure to three data sets to illustrate varying degrees of cladogram ambiguity and recombination.     

Estimation of evolutionary distance between nucleotide sequences

A mathematical formula for estimating the average number of nucleotide substitutions per site (delta) between two homologous DNA sequences is developed by taking into account unequal rates of substitution among different nucleotide pairs. Although this formula is obtained for the equal-input model of nucleotide substitution, computer simulations have shown that it gives a reasonably good estimate for a wide range of nucleotide substitution patterns as long as delta is equal to or smaller than 1. Furthermore, the frequency of cases to which the formula is inapplicable is much lower than that for other similar methods recently proposed. This point is illustrated using insulin genes. A statistical method for estimating the number of nucleotide changes due to deletion and insertion is also developed. Application of this method to globin gene data indicates that the number of nucleotide changes per site increases with evolutionary time but the pattern of the increase is quite irregular.      

Evolutionary change of restriction cleavage sites and phylogenetic inference for man and apes

A mathematical theory for the evolutionary change of restriction endonuclease cleavage sites is developed, and the probabilities of various types of restriction-site changes are evaluated. A computer simulation is also conducted to study properties of the evolutionary change of restriction sites. These studies indicate that parsimony methods of constructing phylogenetic trees often make erroneous inferences about evolutionary changes of restriction sites unless the number of nucleotide substitutions per site is less than 0.01 for all branches of the tree. This introduces a systematic error in estimating the number of mutational changes for each branch and, consequently, in constructing phylogenetic trees. Therefore, parsimony methods should be used only in cases where nucleotide sequences are closely related. Reexamination of Ferris et al.'s data on restriction-site differences of mitochondrial DNAs does not support Templeton's conclusions regarding the phylogenetic tree for man and apes and the molecular clock hypothesis. Templeton's claim that Nei and Li's method of estimating the number of nucleotide substitutions per site is seriously affected by parallel losses and loss-gains of restriction sites is also unsupported.      

基于线粒体Cyt b基因的龙头鱼群体遗传结构分析

在我国黄海、东海和南海海岸带上, 选取青岛、南通、舟山、三门、宁德、泉州和湛江共7个群体164尾龙头鱼(Harpadon nehereus)为研究对象, 通过PCR扩增共获得长度为1112 bp的线粒体DNA细胞色素b (Cyt b)基因序列, 共检测到32个变异位点, 其中单变异位点27个, 简约信息位点1个。164条序列定义了29个单倍型, 平均单倍型多样性(H_d)和核苷酸多样性(π)分别为0.3026±0.0479和0.000371±0.000379, 其中泉州群体的单倍型多样性和核苷酸多样性均最低。分析比较了不同龙头鱼群体间遗传变异情况, 发现群体间平均遗传距离为0.00035, 遗传分化指数F_ST值均小于0.05, 群体间没有发生明显的遗传分化。AMOVA分析结果显示龙头鱼群体遗传差异主要来源于群体内个体间的变异。中性检验的Tajima’s D和Fu’s F_s统计值均为负值且显著偏离中性, 核苷酸不配对分布图呈现明显的单峰分布, 表明龙头鱼历史上经历了群体扩张事件, 参考已知海洋鱼类Cyt b基因2%每百万年的进化速率, 估算群体扩张发生的时间大约在0.08—0.32百万年前的第四纪更新世中晚期。     

Mitochondrial DNA variation of Salvelinus species found in Finland

Mitochondrial DNA (mtDNA) was purified from the Arctic charr, Salvelinus alpinus , the brook charr, Salvelinus fontinalis , and the lake charr, Salvelinus namaycush , and digested with restriction enzymes Ava II, Hinf I, Eco R V, Pst I and Xba I. Two Arctic charr samples were from natural populations and they represented two different morphotypes of Arctic charr. All other studied populations were hatchery maintained. Eight additional restriction enzymes and double digestions were employed to study morphotypes of Arctic charr. We distinguished two morphotypes with restriction enzyme Nci I. Sequence divergence among mtDNA types was 2.9–3.8% between S. alpinus and S. fontinalis , 3.4–4.6% between S. alpinus and S. namaycush , and 4.7–5.3% between S. fontinalis and S. namaycush . lntraspecific variation was lowest in Arctic charr, the average of nucleon diversity for three populations being 0.179, while for brook charr and for lake charr nucleon diversity was 0.334 and 0.550, respectively. According to the number of mtDNA types, it is obvious that introduction to Finland and hatchery propagation have not greatly affected the mtDNA variation of brook charr or lake charr.     

Phylogeography of the Turkish brown trout Salmo trutta L.: mitochondrial DNA PCR‐RFLP variation

In the present study, mitochondrial DNA polymerase chain reaction‐restriction fragment length polymorphism (PCR‐RFLP) assay was used to assess the phylogenetic and phylogeographic relationships among 27 brown trout Salmo trutta populations from Turkey. The complete NADH 5/6 region and a second segment comprising the cytochrome b gene and D‐loop of mtDNA amplified by PCR were digested with six and five restriction enzymes, respectively. A total of 27 haplotypes were observed and divided into three major phylogenetic assemblages, namely Danubian (DA), Adriatic (AD) and a newly proposed Tigris (TI) lineage. The timing of the net nucleotide divergence between the major lineages along with the geological history of Turkey suggested pre‐Pleistocene isolation of the Turkish brown trout and provided evidence that Turkey could be considered as a centre of diversification for these lineages. The average haplotype diversity (0·1397) and the nucleotide diversity (0·000416) within populations were low in comparison to the observed interpopulation nucleotide diversity (0·021266). PCR‐RFLP analysis showed that most of the mtDNA sequence variation found in the Turkish brown trout populations was imputable to differences among lineages. On the other hand, there was also an obvious relationship between geographical distribution of the populations and their clustering. The present study showed that brown trout populations from Turkey are highly divergent and mainly have a unique genetic profile that could be used for conservation and management purposes.     

Geographical variation of mitochondrial DNA between populations of the white-clawed crayfish Austropotamobius pallipes

1. The white-clawed crayfish Austropotamobius pallipes has a widespread distribution in Europe, but since the last century its distribution has been severely affected by a combination of factors, particularly the disease crayfish plague, introduced by foreign crayfish species. At present A. pallipes is considered as vulnerable and endangered and measures are being taken in a number of countries to conserve it.
2. The application of genetics to conservation is of practical value particularly where restocking is being considered. This study assesses levels of genetic differentiation from an analysis of mitochondrial DNA (mtDNA) variation by restriction fragment length polymorphism of four populations sampled on a European scale. Six different composite haplotypes were detected among sixty-four crayfish representing the four populations. mtDNA nucleon diversity values within species ranged from 0.63% to 12.35%. A close genetic similarity was found between British and French populations of A. pallipes .
3. Results show that mtDNA can be as a suitable diagnostic marker for measuring genetic diversity between crayfish populations and that such information is of use in planning crayfish conservation strategies.     

Geographical variation of mitochondrial DNA between populations of the white-clawed crayfish Austropotamobius pallipes

1. The white-clawed crayfish Austropotamobius pallipes has a widespread distribution in Europe, but since the last century its distribution has been severely affected by a combination of factors, particularly the disease crayfish plague, introduced by foreign crayfish species. At present A. pallipes is considered as vulnerable and endangered and measures are being taken in a number of countries to conserve it.
2. The application of genetics to conservation is of practical value particularly where restocking is being considered. This study assesses levels of genetic differentiation from an analysis of mitochondrial DNA (mtDNA) variation by restriction fragment length polymorphism of four populations sampled on a European scale. Six different composite haplotypes were detected among sixty-four crayfish representing the four populations. mtDNA nucleon diversity values within species ranged from 0.63% to 12.35%. A close genetic similarity was found between British and French populations of A. pallipes .
3. Results show that mtDNA can be as a suitable diagnostic marker for measuring genetic diversity between crayfish populations and that such information is of use in planning crayfish conservation strategies.     

分子生态学重要概念——遗传距离及其测度的理论研究概况

综述了遗传距离的概念、背景,有关遗传距离的几种基本的突变模型以及和遗传距离有关的参量和几种常用统计量,指出在处理蛋白质数据、分子数据以及序列数据时,如何选择相应的统计量和可用的软件包,同时还着重指明了各种模型的假设前提,为处理实际的蛋白质或分子数据时选择合适的模型,和对数据的最终解释提供一些帮助。     

Empirical relationship between the number of nucleotide substitutions and interspecific identity of amino acid sequences in some proteins

There are three different methods of estimating the number of nucleotide substitutions between a pair of species from amino acid sequence data, i.e. the Poisson correction method, random evolutionary hit method, and counting the actual but minimum number of nucleotide substitutions. In this paper the relationships among the estimates obtained by these methods are studied empirically. The results obtained indicate that there is a high correlation among these estimates and in practice any of the three methods may be used for constructing evolutionary trees or relating nucleotide substitutions to evolutionary time. The effects of varying rates of nucleotide substition among different sites on the Poisson correction and random evolutionary hit methods are also studied mathematically. It is shown that these two methods are quite insensitive to the variation of the rate of nucleotide substitution.     

A statistical characterization of consistent patterns of human immunodeficiency virus evolution within infected patients

Within-patient HIV populations evolve rapidly because of a high mutation rate, short generation time, and strong positive selection pressures. Previous studies have identified "consistent patterns" of viral sequence evolution. Just before HIV infection progresses to AIDS, evolution seems to slow markedly, and the genetic diversity of the viral population drops. This evolutionary slowdown could be caused either by a reduction in the average viral replication rate or because selection pressures weaken with the collapse of the immune system. The former hypothesis (which we denote "cellular exhaustion") predicts a simultaneous reduction in both synonymous and nonsynonymous evolution, whereas the latter hypothesis (denoted "immune relaxation") predicts that only nonsynonymous evolution will slow. In this paper, we present a set of statistical procedures for distinguishing between these alternative hypotheses using DNA sequences sampled over the course of infection. The first component is a new method for estimating evolutionary rates that takes advantage of the temporal information in longitudinal DNA sequence samples. Second, we develop a set of probability models for the analysis of evolutionary rates in HIV populations in vivo. Application of these models to both synonymous and nonsynonymous evolution affords a comparison of the cellular-exhaustion and immune-relaxation hypotheses. We apply the procedures to longitudinal data sets in which sequences of the env gene were sampled over the entire course of infection. Our analyses (1) statistically confirm that an evolutionary slowdown occurs late in infection, (2) strongly support the immune-relaxation hypothesis, and (3) indicate that the cessation of nonsynonymous evolution is associated with disease progression.     

Variability of Some Milk-Associated Genes and Proteins in Several Breeds of Saudi Arabian Camels

To gain knowledge on the molecular basis of diversity of several clans of Saudi camel (Camelus dromedarius) characterization of these animals was conducted at both genetic and protein levels. To this end, blood and milk samples were collected from several camel breeds at different Saudi Arabia locations (northern Jeddah, Riyadh, and Alwagh governorates). Genomic DNA was extracted from blood of four Saudi camel breeds (Majahem, Safra, Wadha, and Hamara), and DNA fragments of the casein and α-lactalbumin genes were amplified. The retrieved DNA sequences were analyzed for genetic variability. The inter-simple sequence repeat technique was used for confirming the relationships among the analyzed camel breeds, and the PCR–RFLP with two restriction enzymes was utilized for exploring their molecular variations. The number of haplotypes, gene diversity, nucleotide diversity, average number of nucleotide differences, and sequence conservation were calculated for all the analyzed DNA sequences. These analyses revealed the presence of several single nucleotide polymorphisms in the analyzed DNA sequences. A group of neighbor joining trees was built for inferring the evolutionary variations among the studied animals. Protein profiling of milk from different camel clans was also conducted, and differences between and within the Saudi camel clans were easily found based on the isoelectric focusing (IEF) profiles using ampholytes with different IEF range. This study revealed that analyzed camel breeds show low levels of genetic differences. This may be a reflection of the evolutionary history of C. dromedarius that was domesticated based on a highly homogeneous ancestor ecotype.     

Estimation of the Amount of DNA Polymorphism When the Neutral Mutation Rate Varies among Sites

Knowing the amount of DNA polymorphism is essential to understand the mechanism of maintaining DNA polymorphism in a natural population. The amount of DNA polymorphism can be measured by the average number of nucleotide differences per site (π), the proportion of segregating (polymorphic) site (s) and the minimum number of mutations per site (s*). Since the latter two quantities depend on the sample size, θ is often used as a measure of the amount of DNA polymorphism, where θ = 4Nμ, N is the effective population size and μ is the neutral mutation rate per site per generation. It is known that θ estimated from π, s and s* under the infinite site model can be biased when the mutation rate varies among sites. We have therefore developed new methods for estimating θ under the finite site model. Using computer simulations, it has been shown that the new methods give almost unbiased estimates even when the mutation rate varies among sites substantially. Furthermore, we have also developed new statistics for testing neutrality by modifying Tajima's D statistic. Computer simulations suggest that the new test statistics can be used even when the mutation rate varies among sites.     

锁阳ITS序列遗传多样性分析

为阐明锁阳(Cynomorium songaricum)的遗传结构与遗传多样性, 以甘肃省河西走廊地区及青海省共18个居群的188个锁阳个体为研究对象, 利用现代分子生物学技术, 采用序列分析方法从核-质基因方面对遗传结构进行了分析。结果显示, 锁阳ITS序列总长度为687 bp, 含有7个变异位点, 定义9个单倍型, 整体单倍型多态性H_d=0.294 20, 核苷酸多样性π=0.000 49。在整个单倍型网络中介图中, 单倍型H1位于中心位置, 并在所有的居群中均有分布, 为核心的古老单倍型。分子方差分析结果显示, 锁阳种群变异主要来源于种群内。根据ITS序列得到的群体间遗传分化系数以及Mantel检验结果, 锁阳种群间的遗传距离与地理距离之间不存在相关性, 表明现存的锁阳居群是相对近期发生生境片段化的产物。中性检验结果表明, 锁阳拒绝中性进化, 群体历经扩张或者基因座位受到负选择作用, 其中性零假说不能被排除。该研究为锁阳的系统分类、资源鉴定以及保护措施的制定提供了分子证据。     

Assessing the nucleotide diversity of three aphid species by RAPD

A method is presented for the estimation of nucleotide diversity and genetic structure of populations from RAPD (random amplified polymorphic DNA) data. It involves a modification of the technique developed by Lynch and Crease (1990) for the case of restriction sites as survey data. As new elements the method incorporates (i) dominance correction, (ii) values of asexual reproduction of the populations sampled, and (iii) an analytical variance of the number of nucleotide substitutions per site. Sampling was carried out at two geographic scales for three aphid species. At a macrogeographic scale, populations of Rhopalosiphum padi did not show statistical genetic differentiation. Aphis gossypii and Myzus persicae, which were sampled at a microgeographic scale, showed a higher genetic differentiation than R. padi, it being statistically significant in M. persicae. The major sources of sampling variance within- and between-populations were found to be nucleotide (i.e., the number of alleles used as a function of the number of primers used) and population (i.e., sample size) sampling. Extremely low estimates of nucleotide diversity were obtained for the species studied here. This result is consistent with previous reports on genetic diversity for the same or other aphid species which were based on allozyme polymorphism, mitochondrial DNA variation and qualitative analyses of RAPDs.     

Genetic Signature of Rapid IHHNV (Infectious Hypodermal and Hematopoietic Necrosis Virus) Expansion in Wild Penaeus Shrimp Populations

Infectious hypodermal and hematopoietic necrosis virus (IHHNV) is a widely distributed single-stranded DNA parvovirus that has been responsible for major losses in wild and farmed penaeid shrimp populations on the northwestern Pacific coast of Mexico since the early 1990''s. IHHNV has been considered a slow-evolving, stable virus because shrimp populations in this region have recovered to pre-epizootic levels, and limited nucleotide variation has been found in a small number of IHHNV isolates studied from this region. To gain insight into IHHNV evolutionary and population dynamics, we analyzed IHHNV capsid protein gene sequences from 89 Penaeus shrimp, along with 14 previously published sequences. Using Bayesian coalescent approaches, we calculated a mean rate of nucleotide substitution for IHHNV that was unexpectedly high (1.39×10⁻⁴ substitutions/site/year) and comparable to that reported for RNA viruses. We found more genetic diversity than previously reported for IHHNV isolates and highly significant subdivision among the viral populations in Mexican waters. Past changes in effective number of infections that we infer from Bayesian skyline plots closely correspond to IHHNV epizootiological historical records. Given the high evolutionary rate and the observed regional isolation of IHHNV in shrimp populations in the Gulf of California, we suggest regular monitoring of wild and farmed shrimp and restriction of shrimp movement as preventative measures for future viral outbreaks.     

Hierarchical Analysis of Nucleotide Diversity in Geographically Structured Populations

Existing methods for analyzing nucleotide diversity require investigators to identify relevant hierarchical levels before beginning the analysis. We describe a method that partitions diversity into hierarchical components while allowing any structure present in the data to emerge naturally. We present an unbiased version of NEI's nucleotide diversity statistics and show that our modification has the same properties as WRIGHT's F(ST). We compare its statistical properties with several other F(ST) estimators, and we describe how to use these statistics to produce a rooted tree of relationships among the sampled populations in which the mean time to coalescence of haplotypes drawn from populations belonging to the same node is smaller than the mean time to coalescence of haplotypes drawn from populations belonging to different nodes. We illustrate the method by applying it to data from a recent survey of restriction site variation in the chloroplast genome of Coreopsis grandiflora.