Global Population Genomics of Two Subspecies of Cryptosporidium hominis during 500 Years of Evolution

Cryptosporidiosis is a major global health problem and a primary cause of diarrhea, particularly in young children in low- and middle-income countries (LMICs). The zoonotic Cryptosporidium parvum and anthroponotic Cryptosporidium hominis cause most human infections. Here, we present a comprehensive whole-genome study of C. hominis, comprising 114 isolates from 16 countries within five continents. We detect two lineages with distinct biology and demography, which diverged circa 500 years ago. We consider these lineages two subspecies and propose the names C. hominis hominis and C. hominis aquapotentis (gp60 subtype IbA10G2). In our study, C. h. hominis is almost exclusively represented by isolates from LMICs in Africa and Asia and appears to have undergone recent population contraction. In contrast, C. h. aquapotentis was found in high-income countries, mainly in Europe, North America, and Oceania, and appears to be expanding. Notably, C. h. aquapotentis is associated with high rates of direct human-to-human transmission, which may explain its success in countries with well-developed environmental sanitation infrastructure. Intriguingly, we detected genomic regions of introgression following secondary contact between the subspecies. This resulted in high diversity and divergence in genomic islands of putative virulence genes, including muc5 (CHUDEA2_430) and a hypothetical protein (CHUDEA6_5270). This diversity is maintained by balancing selection, suggesting a co-evolutionary arms race with the host. Finally, we find that recent gene flow from C. h. aquapotentis to C. h. hominis, likely associated with increased human migration, maybe driving the evolution of more virulent C. hominis variants.     

Interpreting the genomic landscape of speciation: a road map for finding barriers to gene flow

Speciation, the evolution of reproductive isolation among populations, is continuous, complex, and involves multiple, interacting barriers. Until it is complete, the effects of this process vary along the genome and can lead to a heterogeneous genomic landscape with peaks and troughs of differentiation and divergence. When gene flow occurs during speciation, barriers restricting gene flow locally in the genome lead to patterns of heterogeneity. However, genomic heterogeneity can also be produced or modified by variation in factors such as background selection and selective sweeps, recombination and mutation rate variation, and heterogeneous gene density. Extracting the effects of gene flow, divergent selection and reproductive isolation from such modifying factors presents a major challenge to speciation genomics. We argue one of the principal aims of the field is to identify the barrier loci involved in limiting gene flow. We first summarize the expected signatures of selection at barrier loci, at the genomic regions linked to them and across the entire genome. We then discuss the modifying factors that complicate the interpretation of the observed genomic landscape. Finally, we end with a road map for future speciation research: a proposal for how to account for these modifying factors and to progress towards understanding the nature of barrier loci. Despite the difficulties of interpreting empirical data, we argue that the availability of promising technical and analytical methods will shed further light on the important roles that gene flow and divergent selection have in shaping the genomic landscape of speciation.     

Population history of Berthelot's pipit: colonization, gene flow and morphological divergence in Macaronesia

The fauna of oceanic islands provide exceptional models with which to examine patterns of dispersal, isolation and diversification, from incipient speciation to species level radiations. Here, we investigate recent differentiation and microevolutionary change in Berthelot's pipit (Anthus berthelotii), an endemic bird species inhabiting three Atlantic archipelagos. Mitochondrial DNA sequence data and microsatellite markers were used to deduce probable colonization pathway, genetic differentiation, and gene flow among the 12 island populations. Phenotypic differentiation was investigated based on eight biologically important morphological traits. We found little mitochondrial DNA variability, with only one and four haplotypes for the control region and cytochrome b, respectively. However, microsatellite data indicated moderate population differentiation (FST=0.069) between the three archipelagos that were identified as genetically distinct units with limited gene flow. Both results, combined with the estimated time of divergence (2.5 millions years ago) from the Anthus campestris (the sister species), suggest that this species has only recently dispersed throughout these islands. The genetic relationships, patterns of allelic richness and exclusive alleles among populations suggest the species originally colonized the Canary Islands and only later spread from there to the Madeiran archipelago and Selvagen Islands. Differentiation has also occurred within archipelagos, although to a lesser degree. Gene flow was observed more among the eastern and central islands of the Canaries than between these and the western islands or the Madeiran Islands. Morphological differences were also more important between than within archipelagos. Concordance between morphological and genetic differentiation provided ambiguous results suggesting that genetic drift alone was not sufficient to explain phenotypic differentiation. The observed genetic and morphological differences may therefore be the result of differing patterns of selection pressures between populations, with Berthelot's pipit undergoing a process of incipient differentiation.     

中国南疆栽培杏群体遗传结构的荧光AFLP分析

以中国新疆维吾尔自治区南部塔里木盆地边缘库车、喀什和和田3个栽培杏群体的85个品种类型为试材,利用荧光标记AFLP对群体遗传结构进行了研究,结果表明:对64对EcoRⅠ/MseⅠ引物(其中MseI引物为FAM荧光标记物)进行了筛选,其中8对荧光标记引物谱带清晰,多态性高;同一引物在不同群体以及同一群体不同引物扩增多态性均存在显著差异,种级水平多态带百分率P>库车>和田>喀什3个群体水平;种级水平Nei基因多样度指数H和Shannon信息指数I>库车>和田>喀什3个群体水平。种级水平大于种下群体,各群体以库车最高;杏遗传多样性主要存在于群体内,群体间的遗传分化系数GST为0.0882,即群体间的遗传变异占总变异的8.82%,群体间的基因流Nm为5.1689;群体的遗传一致度在0.9772～0.9811之间,遗传距离在0.0191～0.0232之间,说明群体间的相似程度较高,遗传距离较小;UPGMA聚类分析结果表明,库车、喀什、和田3个亚群可能是相对独立的孟德尔群体,但同时存在部分基因交流;各项研究指标都表明库车群体的遗传多样性都最高,它可能是野杏向栽培杏过渡的中间群体;南疆栽培杏群体拥有极为丰富的遗传多样性,为进一步的选择育种提供了更多的种质资源,为该地区杏群体生物多样性保护和利用提供理论依据。     

Inference of Gene Flow between Species under Misspecified Models

Genomic sequence data provide a rich source of information about the history of species divergence and interspecific hybridization or introgression. Despite recent advances in genomics and statistical methods, it remains challenging to infer gene flow, and as a result, one may have to estimate introgression rates and times under misspecified models. Here we use mathematical analysis and computer simulation to examine estimation bias and issues of interpretation when the model of gene flow is misspecified in analysis of genomic datasets, for example, if introgression is assigned to the wrong lineages. In the case of two species, we establish a correspondence between the migration rate in the continuous migration model and the introgression probability in the introgression model. When gene flow occurs continuously through time but in the analysis is assumed to occur at a fixed time point, common evolutionary parameters such as species divergence times are surprisingly well estimated. However, the time of introgression tends to be estimated towards the recent end of the period of continuous gene flow. When introgression events are assigned incorrectly to the parental or daughter lineages, introgression times tend to collapse onto species divergence times, with introgression probabilities underestimated. Overall, our analyses suggest that the simple introgression model is useful for extracting information concerning between-specific gene flow and divergence even when the model may be misspecified. However, for reliable inference of gene flow it is important to include multiple samples per species, in particular, from hybridizing species.     

物种形成过程中的分化基因组岛及其形成机制

理解物种形成机制是生态和进化领域的重要任务。得益于测序技术的快速发展, 越来越多研究发现分化种群(亚种、物种)间的基因组常呈现异质性分化景观, 存在分化基因组岛, 这被认为是基因流存在下的歧化选择引起的, 支持基因流存在下的成种假说。然而, 基因渐渗、祖先多态性的差异分选、连锁选择等其他进化过程也可导致分化基因组岛的形成。现有实证研究在解析分化基因组岛的形成机制时, 往往忽略了上述其他进化过程的作用。为此, 本文在辨析分化基因组岛相关概念的基础上, 总结了利用种群基因组数据鉴定分化基因组岛的方法, 对比了不同进化过程形成分化基因组岛的特征, 指出在区分不同机制时联用基因渐渗程度、绝对分化指数(d_XY)、相对节点深度(RND)、重组率等多个指标的必要性, 归纳了物种形成过程中分化基因组岛形成机制解析的研究思路, 并对未来在生殖隔离机制上的深入探索以及实证研究的整合分析等方面进行了展望。     

Genomic differentiation and patterns of gene flow between two long‐tailed tit species (Aegithalos)

Patterns of heterogeneous genomic differentiation have been well documented between closely related species, with some highly differentiated genomic regions (“genomic differentiation islands”) spread throughout the genome. Differential levels of gene flow are proposed to account for this pattern, as genomic differentiation islands are suggested to be resistant to gene flow. Recent studies have also suggested that genomic differentiation islands could be explained by linked selection acting on genomic regions with low recombination rates. Here, we investigate genomic differentiation and gene‐flow patterns for autosomes using RAD‐seq data between two closely related species of long‐tailed tits (Aegithalos bonvaloti and A. fuliginosus) in both allopatric and contact zone populations. The results confirm recent or ongoing gene flow between these two species. However, there is little evidence that the genomic regions that were found to be highly differentiated between the contact zone populations are resistant to gene flow, suggesting that differential levels of gene flow is not the cause of the heterogeneous genomic differentiation. Linked selection may be the cause of genomic differentiation islands between the allopatric populations with no or very limited gene flow, but this could not account for the heterogeneous genomic differentiation between the contact zone populations, which show evidence of recent or ongoing gene flow.     

Gene flow and selection interact to promote adaptive divergence in regions of low recombination

Adaptation to new environments often occurs in the face of gene flow. Under these conditions, gene flow and recombination can impede adaptation by breaking down linkage disequilibrium between locally adapted alleles. Theory predicts that this decay can be halted or slowed if adaptive alleles are tightly linked in regions of low recombination, potentially favouring divergence and adaptive evolution in these regions over others. Here, we compiled a global genomic data set of over 1,300 individual threespine stickleback from 52 populations and compared the tendency for adaptive alleles to occur in regions of low recombination between populations that diverged with or without gene flow. In support of theory, we found that putatively adaptive alleles (F_ST and d_XY outliers) tend to occur more often in regions of low recombination in populations where divergent selection and gene flow have jointly occurred. This result remained significant when we employed different genomic window sizes, controlled for the effects of mutation rate and gene density, controlled for overall genetic differentiation, varied the genetic map used to estimate recombination and used a continuous (rather than discrete) measure of geographic distance as proxy for gene flow/shared ancestry. We argue that our study provides the first statistical evidence that the interaction of gene flow and selection biases divergence toward regions of low recombination.     

Population Genetic Structure of the Ark Shell <Emphasis Type="Italic">Scapharca broughtonii</Emphasis> Schrenck from Korea,China, and Russia Based on <Emphasis Type="Italic">COI</Emphasis> Gene Sequences

Haplotype distribution, gene flow, and population genetic structure of the ark shell (Scapharca broughtonii) were studied using a partial sequence of a mitochondrial COI gene. The sequence analysis of 100 specimens obtained from a total of seven localities-five in Korea, one in China, and one in Russia- revealed 29 haplotypes, ranging in sequence divergence from 0.1% to 2.1%. Among these, the most frequent haplotype, SB16, was extensively distributed over study areas, especially in all Korean localities. This extensive distribution consequently resulted in the near absence of statistically significant genetic distance. Also, a high rate of gene flow was characteristic among localities in Korea. A test of genetic population structure showed that the ark shell in Korea formed a large genetic group. Moreover, an AMOVA test to determine the allocation of the genetic variance showed that most of the variance was distributed between localities, instead of within localities. However, a significant population differentiation was found between geographic populations [i.e., Jinhae (locality 6) in Korea and Sangdong (locality 5) in China and Vladivostok (locality 7) in Russia] based on geographic distance and population structure. These distinct groups may be associated with geographic characteristics and barriers. The results suggest that most of the ark shell populations in Korea caused considerable distribution to form a genetically homogeneous and intermixing structure, whereas some of the Korean and Chinese and Russian populations had a significantly different genetic structure.     

Genetic structure of Polytrichum formosum in relation to the breeding system as revealed by microsatellites

Microsatellite variation was determined for three Danish and three Dutch populations of the haploid moss species Polytrichum formosum to gain insight into the relative importance of sexual vs. asexual reproduction for the amount and structure of genetic variation. In general, low levels of microsatellite variation were observed within this species. Even when estimated for polymorphic loci only, the levels of microsatellite variability (P=90.6, A=4.3 and H_S=0.468) within populations were on average lower than those reported for most other plant species. In contrast, genotypic diversity was high within each of the examined populations, indicating that sexual reproduction is a very important determinant of the genetic structure of P. formosum within populations. In agreement with previous findings for allozyme data, no significant genetic differentiation (F_ST=0.028, R_ST=0.015) was observed neither between populations nor between regions approximately 450 km apart (Denmark vs. the Netherlands). These low levels of population differentiation observed for both types of genetic markers are probably best explained by a high level of effective spore dispersal (gene flow) between populations. Therefore, also on a large geographical scale sexual reproduction is the most important determinant of the genetic structure of P. formosum, despite the high potential to reproduce clonally.