gener: a server‐based analysis of pollen pool structure

The server‐based program gener performs the two‐generation analysis of pollen flow for data consisting of mother/offspring arrays using genetic markers. The gener program decomposes the genetic variance sampled by maternal individuals within and among pollen pool components of genetic variance and is accessible from http://dyerlab.bio.vcu.edu . These estimates are used to construct the test statistic, Φ_ft , whose significance is tested via permutation. The Φ_ft statistic can subsequently be used to quantify genetic effective pollen donor population size (N_ep), effective mating area and dispersal distance. Furthermore, the gener program can calculate Φ_ft values for all pairs of substrata within the data set.     

圈养小熊猫遗传多样性与种群遗传结构分析

小熊猫是亚洲特有的珍稀濒危动物,目前受到栖息地减少、片断化和人类活动干扰等威胁。中国圈养小熊猫已经有60 多年历史,约55 个机构曾经饲养过小熊猫,现今圈养数量有400 多只,评估小熊猫圈养种群的遗传多样性和遗传结构对科学维持圈养种群和保存遗传种质资源意义重大。本研究利用19 个微卫星座位,对中国境内11 个小熊猫圈养种群的116 只个体进行了遗传多样性评估及遗传结构分析。结果显示11 个种群都具有较高的遗传多样性,平均基因丰富度3.505 ± 1.033 (北京)至4.026 ± 1.219 (冕宁),期望杂合度0.631 ± 0.225(黄山)至0.782 ±0.171 (温岭)。其中福州和无锡种群极显著偏离Hardy-Weinberg 平衡。整个圈养群体内各个种群遗传分化系数为0.055,呈显著分化,表明11 个种群遗传分化水平较高。Bayesian 遗传聚类分析将11 个种群聚为三个遗传簇,与野生种群的遗传聚类结果一致。结论:小熊猫圈养种群与野生种群相比,同样具有较高的遗传多样性。因此,圈养小熊猫遗传管理的重点不再是引进野生个体充实圈养种群,应制订科学的繁殖计划,避免近交,从而维持圈养种群的遗传多样性。     

三峡库区特有种疏花水柏枝的保护遗传学研究

采用超薄平板微型聚丙烯酰胺等电聚焦电泳方法对三峡库区特有种疏花水柏枝 (Myricarialaxiflora)的 13个自然居群和 1个人工迁地保护居群的等位酶遗传变异进行了初步研究。检测了 5个酶系统 ,得到 13个等位酶位点 ,遗传多样性及其遗传结构分析结果表明 :疏花水柏枝具有较高水平的遗传多样性 ,平均多态位点比率P =78.7% ,每位点平均等位基因数A =1.8,平均预期杂合度He=0 .317,高于中国植物特有种的平均水平 ,且群体中杂合基因型个体偏多 ;其遗传变异主要发生于居群内 ( 84.86 % ) ,居群间又存在一定的遗传分化 (Gst=0 .15 14) ,居群间平均基因流Nm =1.40 1,居群间遗传距离为 0 .0 0 2～ 0 .176 ;UPGMA聚类分析显示疏花水柏枝在三峡湖北境内的白水河—泄滩一带分为上游和下游 2个居群组 ;武汉植物园迁地保护的混合居群基本保育了其遗传多样性总水平。在分析讨论疏花水柏枝遗传多样性与其繁育系统、生境及其起源进化的关系的基础上 ,探讨了疏花水柏枝濒危的主要原因 ,推断其应为第四纪冰期影响后的古孑遗种。最后 ,在评价迁地保护成果的基础上 ,提出了今后进一步保育的策略。     

群体遗传学研究中的数据处理方法I．RAPD数据的AMOVA分析

近年来,RAPD数据和AMOVA分析广泛地应用于群体遗传学和保护遗传学研究。然而,由于RAPD标记具显性特点。加上目前进行AMOVA分析所依赖的RAPDistance软件不完善,使得对RAPD数据进行AMOVA分析时存在许多不足。本文介绍了AMOVA分析的基本过程,同时引入一个新的程序DCFA用以替代RADistance并详述了将DCFA与WINAMOVA联用,对RAPD数据进行AMOVA分析的具体步骤与注意事项,最后,以产自中国和巴西8个普通野生稻（Oryza furipogon)天然群体为例,演示了对RAPD表型数据进行AMOVA分析的过程,讨论了AMOVA分析结果在群体遗传结构上的意义。通过对AMOVA算法的分析,同时比较4种距离系数所得AMOVA结果,我们认为在进行AMOVA分析时选择NEI－LI距离和欧氏距离平方较为合适,而目前国内使用较多的JACCARD系数不适合AMOVA分析。     

图论在中国三尖杉属植物区系地理成分划分中的应用

本文根据图论中树与生成树的理论,选用中国三尖杉属(Cephalotaxus)8种植物的地理分布资料,采用图论分类法——Prim法,对国产三尖杉属植物地理分布进行定量研究。结果表明,本文应用图论分类法划分植物区系地理成分基本上是成功的,其结果与传统分析基本相符,尤其是采用不对称信息系数构成相异矩阵的方法;此外,本方法还有以下特点:(1)结果直观、简单;(2)能处理非常复杂的植物地理分布问题;(3)能避免因人为主观意识所引起的片面性甚至错误;(4)能准确地反映各类区系地理成分之间及其在地质时期植物地理联系的“亲缘关系”。     

A program to compare genetic differentiation statistics across loci using resampling of individuals and loci

Comparisons of genetic differentiation across populations based on different loci can provide insight into the evolutionary patterns acting on various regions of genomes. Here, we develop a program to statistically compare population genetic differentiation statistics (F(ST) or G'(ST) ) calculated from different loci. The program employs a routine that resamples either or both of individuals and loci and calculates a bootstrap confidence interval in the statistics. Resampling individuals is important when fewer than 25 individuals are sampled per population and when confidence intervals are required for individual loci. Resampling loci provides confidence intervals for sets of loci, such as a set presumed to be neutral, but can be anticonservative if fewer than 20 loci are analysed. We demonstrate the program using previously published data on the genetic differentiation at a major histocompatibility complex locus and at microsatellite loci across 10 populations of the guppy (Poecilia reticulata).     

The computer program structure for assigning individuals to populations: easy to use but easier to misuse

The computer program Structure implements a Bayesian method, based on a population genetics model, to assign individuals to their source populations using genetic marker data. It is widely applied in the fields of ecology, evolutionary biology, human genetics and conservation biology for detecting hidden genetic structures, inferring the most likely number of populations (K), assigning individuals to source populations and estimating admixture and migration rates. Recently, several simulation studies repeatedly concluded that the program yields erroneous inferences when samples from different populations are highly unbalanced in size. Analysing both simulated and empirical data sets, this study confirms that Structure indeed yields poor individual assignments to source populations and gives frequently incorrect estimates of K when sampling is unbalanced. However, this poor performance is mainly caused by the adoption of the default ancestry prior, which assumes all source populations contribute equally to the pooled sample of individuals. When the alternative ancestry prior, which allows for unequal representations of the source populations by the sample, is adopted, accurate individual assignments could be obtained even if sampling is highly unbalanced. The alternative prior also improves the inference of K by two estimators, albeit the improvement is not as much as that in individual assignments to populations. For the difficult case of many populations and unbalanced sampling, a rarely used parameter combination of the alternative ancestry prior, an initial ALPHA value much smaller than the default and the uncorrelated allele frequency model is required for Structure to yield accurate inferences. I conclude that Structure is easy to use but is easier to misuse because of its complicated genetic model and many parameter (prior) options which may not be obvious to choose, and suggest using multiple plausible models (parameters) and K estimators in conducting comparative and exploratory Structure analysis.     

F2breed,a New Program for Construction of Genetic Linkage Maps

This paper presents a new program F2breed for Windows and Linux operating systems, which was designed for estimation of recombination frequency between the genetic loci and construction of genetic maps based on the analysis of inheritance in a F2 population. The program implements the approach of the orientation of n points corresponding to n analyzed loci in the (n–1)-dimensional space for relative arrangement of the loci on the genetic map by consideration of various projections of this structure. The program is characterized by a convenient, intuitive interface, is easy to operate, and is suitable for working with small groups of loci (up to 100).     

湖北河岸带植物中华蚊母树遗传多样性的SRAP分析

利用SRAP分子标记技术,对湖北省河岸带植物中华蚊母树的4个自然居群和1个迁地居群的遗传多样性和遗传结构进行了分析。结果表明,中华蚊母树物种具有较高水平的遗传多样性,7对SRAP引物进行PCR扩增的多态性位点百分率(PPF)为80.43%,每个位点的等位基因数(A)为2,有效等位基因数(Ae)为1.34,总遗传多样性Nei’s基因多样性指数(Hp)为0.215 9,Shannon信息多样性指数(I)为0.350 9。在居群水平上,5个居群总的遗传变异Ht为0.218 8,居群内的遗传变异Hs为0.193 4,居群间的遗传分化系数Gst为0.116 1,表明在总的遗传变异中有88.39%变异存在于居群内,仅11.61%存在于居群间,居群间的基因流Nm为3.807 2,表明居群间有较大程度的基因交流。UPGMA聚类分析和主成分分析显示中华蚊母树主要分为两个居群组,在长江三峡沿岸香溪和乐天溪由于遗传距离比较近聚为一小类再与高家堰聚为一大类,而沿渡河和三峡植物园居群聚为另一大类,表明迁地居群三峡植物园的中华蚊母树与来自巴东沿渡河居群的样本亲缘关系最近,且三峡植物园迁地保护居群基本保育了其遗传多样性总水平。同时在分析讨论了中华蚊母树遗传多样性与其繁育系统、生境及其起源进化的关系的基础上,评价了中华蚊母树的保护策略,并在评价保护成果的基础上,提出了今后进一步保育的策略。结果还表明SRAP标记是分析中华蚊母树遗传多样性和遗传结构非常可靠的一种标记,而且这是使用SRAP标记研究中华蚊母树的首次报道。     

Genetic variability,management, and conservation implications of the critically endangered Brazilian pitviper Bothrops insularis

Information on demographic, genetic, and environmental parameters of wild and captive animal populations has proven to be crucial to conservation programs and strategies. Genetic approaches in conservation programs of Brazilian snakes remain scarce despite their importance for critically endangered species, such as Bothrops insularis, the golden lancehead, which is endemic to Ilha da Queimada Grande, coast of São Paulo State, Brazil. This study aims to (a) characterize the genetic diversity of ex situ and in situ populations of B. insularis using heterologous microsatellites; (b) investigate genetic structure among and within these populations; and (c) provide data for the conservation program of the species. Twelve informative microsatellites obtained from three species of the B. neuwiedi group were used to access genetic diversity indexes of ex situ and in situ populations. Low‐to‐medium genetic diversity parameters were found. Both populations showed low—albeit significant—values of system of mating inbreeding coefficient, whereas only the in situ population showed a significant value of pedigree inbreeding coefficient. Significant values of genetic differentiation indexes suggest a small differentiation between the two populations. Discriminant analysis of principal components (DAPC) recovered five clusters. No geographic relationship was found in the island, suggesting the occurrence of gene flow. Also, our data allowed the establishment of six preferential breeding couples, aiming to minimize inbreeding and elucidate uncertain parental relationships in the captive population. In a conservation perspective, continuous monitoring of both populations is demanded: it involves the incorporation of new individuals from the island into the captive population to avoid inbreeding and to achieve the recommended allelic similarity between the two populations. At last, we recommend that the genetic data support researches as a base to maintain a viable and healthy captive population, highly genetically similar to the in situ one, which is crucial for considering a reintroduction process into the island.     

Combining multiple analytical approaches for the identification of population structure and genetic delineation of two subspecies of the endemic Arabian burnet moth Reissita simonyi (Zygaenidae; Lepidoptera)

Habitat fragmentation and landscape topology may influence the genetic structure and connectivity between natural populations. Six microsatellite loci were used to infer the population structure of 35 populations (N = 788) of the alpine Arabian burnet moth Reissita simonyi (Lepidoptera, Zygaenidae) in Yemen and Oman. Due to the patchy distribution of larval food plants, R. simonyi is not continuously distributed throughout the studied area and the two recognized subspecies of this endemic species (Reissita s. simonyi/R. s. yemenicola) are apparently discretely distributed. All microsatellites showed prevalence of null alleles and therefore a thorough investigation of the impact of null alleles on different population genetic parameters (F _ST, inbreeding coefficients, and Population Graph topologies) is given. In general, null alleles reduced genetic covariance and independence of allele frequencies resulting in a more connected genetic topology in Population Graphs and an overestimation of pairwise F _ST values and inbreeding coefficients. Despite the presence of null alleles, Population Graphs also showed a much higher genetic connectivity within subspecies (and lower genetic differentiation (via F _ST)) than between; supporting existing taxonomic distinction. Partial Mantel tests showed that both geographical distance and altitude were highly correlated with the observed distribution of genetic structure within R. simonyi. In conclusion, we identified geographical and altitudinal distances in R. simonyi as well as an intervening desert area to be the main factors for spatial genetic structure in this species and show that the taxonomic division into two subspecies is confirmed by genetic analysis.     

Graph theoretic topology of the Great but small Barrier Reef world

The transport of larvae between coral reefs is critical to the functioning of Australia’s Great Barrier Reef (GBR) because it determines recruitment rates and genetic exchange. One way of modelling the transport of larvae from one reef to another is to use information about currents. However the connectivity relationships of the entire system have not been fully examined. Graph theory provides a framework for the representation and analysis of connections via larval transport. In the past, the geometric arrangement (topology) of biological systems, such as food webs and neural networks, has revealed a common set of characteristics known as the ‘small world’ property. We use graph theory to examine and describe the topology and connectivity of a species living in 321 reefs in the central section of the GBR over 32 years. This section of the GBR can be described by a directional weighted graph, and we discovered that it exhibits scale-free small-world characteristics. The conclusion that the GBR is a small-world network for biological organisms is robust to variation in both the life history of the species modelled and yearly variation in hydrodynamics. The GBR is the first reported mesoscale biological small-world network.     

综合分离分析的微机程序设计及其应用

基于单基因遗传模型,我们在IBM PC/XT微型计算机上,用BASIC语言设计了综合分离分析程序。并应用该计算机程序分析了133个用家族法和家族史法调查的寻常性鱼鳞癣核心家系,98个仅用家族法调查的寻常性鱼鳞癣核心家系和101个原因不明的父母表型正常的感音神经性耳聋核心家系。结果表明,该计算机程序可用于遗传方式分析和遗传咨询。     

Genetic analyses of the federally endangered Echinacea laevigata using amplified fragment length polymorphisms (AFLP)—Inferences in population genetic structure and mating system

Echinacea laevigata (Boynton and Beadle) Blake is a federally endangered flowering plant species restricted to four states in the southeastern United States. To determine the population structure and outcrossing rate across the range of the species, we conducted AFLP analysis using four primer combinations for 22 populations. The genetic diversity of this species was high based on the level of polymorphic loci (200 of 210 loci; 95.24%) and Nei’s gene diversity (ranging from 0.1398 to 0.2606; overall 0.2611). There was significant population genetic differentiation (G_ST = 0.294; Ө^II = 0.218 from the Bayesian f = 0 model). Results from the AMOVA analysis suggest that a majority of the genetic variance is attributed to variation within populations (70.26%), which is also evident from the PCoA. However, 82% of individuals were assigned back to the original population based on the results of the assignment test. An isolation by distance analysis indicated that genetic differentiation among populations was a function of geographic distance, although long-distance gene dispersal between some populations was suggested from an analysis of relatedness between populations using the neighbor-joining method. An estimate of the outcrossing rate based on genotypes of progenies from six of the 22 populations using the multilocus method from the program MLTR ranged from 0.780 to 0.912, suggesting that the species is predominantly outcrossing. These results are encouraging for conservation, signifying that populations may persist due to continued genetic exchange sustained by the outcrossing mating system of the species.     

Do rivers function as genetic barriers for the plateau wood frog at high elevations?

Using the plateau wood frog Rana kukunoris from the Hengduan Mountains as a model system, we tested whether rivers form significant genetic barriers (the riverine barrier hypothesis) to high elevation amphibians. Samples were collected from eight sites across three major river drainages, the Min, the Dadu and the Yalong Rivers, and the population genetic structure of these frogs was evaluated with data from eight microsatellite DNA loci. A large amount of genetic structure was found, and the pairwise F _ST ranged from 0.022 to 0.508 and a global F _ST was 0.215. Both analysis of molecular variance and isolation by distance analysis suggested that rivers, mountain ridges and geographic distances all contributed significantly to the population structure. However, no single landscape has prominent barrier effect to the plateau wood frog populations. An assignment analysis using the computer program Structure grouped the eight populations into four population clusters, and no single type of landscape can sufficiently explain the clustering. In conclusion, rivers do not appear to be the leading genetic barriers for the plateau wood frog. The strong population genetic structure is likely the consequence of attributes of the species, as opposed to environmental fragmentation, and the barrier effect of the landscapes is largely swamped by the large amount of intrinsic population structure.     

hybridlab (version 1.0): a program for generating simulated hybrids from population samples

We present the computer program hybridlab 1.0 for simulating intraspecific hybrids from population samples of nuclear genetic markers such as microsatellites, allozymes or SNPs (single nucleotide polymorphisms). The program generates a user‐specified number of multilocus F1 hybrid genotypes between any pair of potentially hybridizing populations included in a standard input‐file of multilocus genotypes for population genetic analysis. This simple, user‐friendly program has a wide range of applications for studying natural and artificial hybridization; in particular, for evaluating the statistical power for individual assignment of parental and hybrid individuals. An example of application for Atlantic cod populations is given.     

Identifying the number of population clusters with structure: problems and solutions

The program structure has been used extensively to understand and visualize population genetic structure. It is one of the most commonly used clustering algorithms, cited over 11 500 times in Web of Science since its introduction in 2000. The method estimates ancestry proportions to assign individuals to clusters, and post hoc analyses of results may indicate the most likely number of clusters, or populations, on the landscape. However, as has been shown in this issue of Molecular Ecology Resources by Puechmaille ( 2016 ), when sampling is uneven across populations or across hierarchical levels of population structure, these post hoc analyses can be inaccurate and identify an incorrect number of population clusters. To solve this problem, Puechmaille ( 2016 ) presents strategies for subsampling and new analysis methods that are robust to uneven sampling to improve inferences of the number of population clusters.     

Mitochondrial genetic variation and population structure of the striped snakehead,Channa striata in Malaysia and Sumatra,Indonesia

We explored the genetic diversity and structure of the striped snakehead (Channa striata) across Malaysia and Sumatra (Indonesia) using the partial mtDNA CO1 gene. Twenty five populations (n = 345) were assayed and subdivided into six regions, following the physiogeographical barriers. Populations Sega (SG), Tanjung Tambutan (TR), Kajang (KJ) and Linggi (LG) are highly diversified (Hd: 0.484–0.762, π: 0.0033–0.0059) which could serve as candidates for a selective breeding programme. The only population that contributed to the total allelic richness is Takengon (CS) as it is highly differentiated from other populations and genetically variable within population. We detected two major phylogenies: 1) northwest Peninsular Malaysia and 2) all other regions of Malaysia and Sumatra, Indonesia. They are products of the physical restriction to gene flow between the two lineages by the Bintang Mountain Range. A total of 92.4% of the population pairwise comparison F_ST showed significant structuring, yet several geographically distant populations showed a close genetic relationship. The discrepancy is due to ancient population dispersal and human-mediated translocation. These major findings provide an important base study for initiating a selective breeding program. The high population genetic diversity requires independent conservation as they contain most of the total diversity in this area.     

StructureSelector: A web‐based software to select and visualize the optimal number of clusters using multiple methods

Inferences of population genetic structure are of great importance to the fields of ecology and evolutionary biology. The program structure has been widely used to infer population genetic structure. However, previous studies demonstrated that uneven sampling often leads to wrong inferences on hierarchical structure. The most widely used ΔK method tends to identify the uppermost hierarchy of population structure. Recently, four alternative statistics (medmedk , medmeak , maxmedk and maxmeak ) were proposed, which appear to be more accurate than the previously used methods for both even and uneven sampling data. However, the lack of easy‐to‐use software limits the use of these appealing new estimators. Here, we developed a web‐based user‐friendly software structureselector to calculate the four appealing alternative statistics together with the commonly used Ln Pr(X|K) and ΔK statistics. structureselector accepts the result files of structure , admixture or faststructure as input files. It reports the “best” K for each estimator, and the results are available as HTML or tab separated tables. The program can also generate graphical representations for specific K, which can be easily downloaded from the server. The software is freely available at http://lmme.qdio.ac.cn/StructureSelector/ .