Gene Flow Events Research Articles

Detectability of Varied Hybridization Scenarios Using Genome-Scale Hybrid Detection Methods

Hybridization events complicate the accurate reconstruction of phylogenies, as they lead to patterns of genetic heritability that are unexpected under traditional, bifurcating models of species trees. This phenomenon has led to the development of methods to infer these varied hybridization events, both methods that reconstruct networks directly, as well as summary methods that predict individual hybridization events from a subset of taxa. However, a lack of empirical comparisons between methods – especially those pertaining to large networks with varied hybridization scenarios – hinders their practical use. Here, we provide a comprehensive review of popular summary methods: TICR, MSCquartets, HyDe, Patterson’s D-Statistic (ABBA-BABA), D3, and Dp. TICR and MSCquartets are based on quartet concordance factors gathered from gene tree topologies and HyDe, Patterson’s D-Statistic, D3, and Dp use site pattern frequencies to identify hybridization events between sets of three taxa. We then use simulated data to address questions of method accuracy and ideal use scenarios by testing methods against complex networks which depict gene flow events that differ in depth (timing), quantity (single vs. multiple, overlapping hybridizations), and rate of gene flow (γ). We find that deeper or multiple hybridization events may introduce noise and weaken the signal of hybridization, leading to higher relative false negative rates across all methods. Despite some forms of hybridization eluding quartet-based detection methods, MSCquartets displays high precision in most scenarios. While HyDe results in high false negative rates when tested on hybridizations involving extinct or unsampled ghost lineages, HyDe is the only method able to identify the direction of hybridization, distinguishing the source parental lineages from recipient hybrid lineages. Lastly, we test the methods on a dataset of ultraconserved elements from the bee subfamily Nomiinae, finding possible hybridization events between clades which correspond to regions of poor support in the species tree estimated in a previous study.

Gelada genomes highlight events of gene flow, hybridisation and local adaptation that track past climatic changes.

Theropithecus gelada, the last surviving species of this genus, occupy a unique and highly specialised ecological niche in the Ethiopian highlands. A subdivision into three geographically defined populations (Northern, Central and Southern) has been tentatively proposed for this species on the basis of genetic analyses, but genomic data have been investigated only for two of these groups (Northern and Central). Here we combined newly generated whole genome sequences of individuals sampled from the population living south of the East Africa Great Rift Valley with available data from the other two gelada populations to reconstruct the evolutionary history of the species. Integrating genomic and paleoclimatic data we found that gene-flow across populations and with Papio species tracked past climate changes. The isolation and climatic conditions experienced by Southern geladas during the Holocene shaped local diversity and generated diet-related genomic signatures.

Genome analysis reveals three distinct lineages of the cosmopolitan white shark

The white shark (Carcharodon carcharias) (Linnaeus, 1758), an iconic apex predator occurring in all oceans,1,2 is classified as Vulnerable globally3-with global abundance having dropped to 63% of 1970s estimates,4-and as Critically Endangered in Europe.5 Identification of evolutionary significant units and their management are crucial for conservation,6 especially as the white shark is facing various but often region-specific anthropogenic threats.7,8,9,10,11 Assessing connectivity in a cosmopolitan marine species requires worldwide sampling and high-resolution genetic markers.12 Both are lacking for the white shark, with studies to date typified by numerous but geographically limited sampling, and analyses relying largely on relatively small numbers of nuclear microsatellites,13,14,15,16,17,18,19 which can be plagued by various genotyping artefacts and thus require cautious interpretation.20 Sequencing and computational advances are finally allowing genomes21,22,23 to be leveraged into population studies,24,25,26,27 with datasets comprising thousands of single-nucleotide polymorphisms (SNPs). Here, combining target gene capture (TGC)28 sequencing (89 individuals, 4,000 SNPs) and whole-genome re-sequencing (17 individuals, 391,000 SNPs) with worldwide sampling across most of the distributional range, we identify three genetically distinct allopatric lineages (North Atlantic, Indo-Pacific, and North Pacific). These diverged 100,000-200,000 years ago during the Penultimate Glaciation, when low sea levels, different ocean currents, and water temperatures produced significant biogeographic barriers. Our results show that without high-resolution genomic analyses of samples representative of a species' range,12 the true extent of diversity, presence of past and contemporary barriers to gene flow, subsequent speciation, and local evolutionary events will remain enigmatic.

Unveiling a potential threat to forest ecosystems: molecular diagnosis of Alliaria petiolata, a newly introduced alien plant in Korea.

Identifying stages of a species invasion in a new habitat (i.e., colonization, establishment, and landscape spread) and their primary determinants in biological invasion warrants attention, as it provides vital insights for preventing non-native species from becoming pervasive invaders. However, delineating invasion stages and their associated factors can pose significant challenges due to the ambiguous distinctions between these stages. Alliaria petiolata, one of the most noxious weeds in woodland habitats, has recently been introduced to Korea and observed in a few distant locations. Although the plant's spread has been relatively slow thus far, rapid spread is highly likely in the future, given the high invasive potential reported elsewhere. We indirectly diagnose the current status of A. petiolata invasion in Korea through the assessment of genetic diversity and phylogenetic inferences using genome-wide molecular markers and cytological data. We analyzed 86 individual samples collected from two native and six introduced populations, employing 1,172 SNPs. Our analysis estimated within- and among-population genetic diversity and included two clustering analyses. Furthermore, we investigated potential gene flow and reticulation events among the sampled populations. Our data unraveled that Korean garlic mustard exhibits a hexaploid ploidy level with two distinct chromosome numbers, 2n = 36 and 42. The extent of genetic diversity measured in Korean populations was comparable to that of native populations. Using genome-wide SNP data, we identified three distinct clusters with minor gene flow, while failing to detect indications of reticulation among Korean populations. Based on the multifaceted analyses, our study provides valuable insights into the colonization process and stressed the importance of closely monitoring A. petiolata populations in Korea.

The relevance of Evolutionary Significant Units for the conservation of island-restricted reptiles: Tarentola boettgeri bischoffi as a case study

Abstract Within vertebrates, reptiles are good island colonisers, often leading to considerable levels of intraspecific diversity among populations inhabiting different islands/archipelagos. This study explores the mitochondrial phylogeographic structure of Tarentola boettgeri, a gecko species endemic to the Macaronesian archipelagos of Selvagens and the Canary Islands. Our research introduces a novel monophyletic group, comprising the populations from the islands of Selvagem Pequena and Ilhéu de Fora. Furthermore, we confirm the previously identified genetic clusters associated with Selvagem Grande, Gran Canaria and El Hierro. We estimate that the origin of T. boettgeri dates to the upper Miocene (ca. 6.4 Mya), and that the separation of T. boettgeri bischoffi on Selvagem Grande, Selvagem Pequena, and Ilhéu de Fora, occurred ca. 0.5 Mya. The absence of genetic differences between Selvagem Pequena and Ilhéu de Fora suggests recent gene flow or founder events, possibly facilitated by land connections during major glaciations. Conversely, the geographic barriers between Selvagem Grande and Selvagem Pequena likely persisted, preventing genetic admixing. The significant genetic distances observed among all populations underscore the necessity of an integrative taxonomic revision for T. boettgeri. In light of our findings, and with particular consideration of the small population sizes of T. boettgeri bischoffi on Selvagem Pequena and Ilhéu de Fora, we propose that the identified monophyletic groups should be managed as Evolutionarily Significant Units (ESUs). Accordingly, our study highlights the importance of recognizing ESUs in island-restricted reptile populations for targeted conservation efforts, especially given their unique intraspecific diversity and the vulnerability of their habitats.

Recurrent gene flow events occurred during the diversification of clownfishes of the skunk complex.

Clownfish (subfamily Amphiprioninae) are an iconic group of coral reef fish that evolved a mutualistic interaction with sea anemones, which triggered the adaptive radiation of the clade. Within clownfishes, the "skunk complex" is particularly interesting. Besides ecological speciation, interspecific gene flow and hybrid speciation are thought to have shaped the evolution of the group. We investigated the mechanisms characterizing the diversification of this complex. By taking advantage of their disjunct geographical distribution, we obtained whole-genome data of sympatric and allopatric populations of the three main species of the complex (Amphiprion akallopisos, A. perideraion and A. sandaracinos). We examined population structure, genomic divergence and introgression signals and performed demographic modelling to identify the most realistic diversification scenario. We excluded scenarios of strict isolation or hybrid origin of A. sandaracinos. We discovered moderate gene flow from A. perideraion to the ancestor of A. akallopisos + A. sandaracinos and weak gene flow between the species in the Indo-Australian Archipelago throughout the diversification of the group. We identified introgressed regions in A. sandaracinos and detected in A. perideraion two large regions of high divergence from the two other species. While we found that gene flow has occurred throughout the species' diversification, we also observed that recent admixture was less pervasive than initially thought, suggesting a role of host repartition or behavioural barriers in maintaining the genetic identity of the species in sympatry.

The role of historical biogeography in shaping colour morph diversity in the common wall lizard.

The maintenance of polymorphisms often depends on multiple selective forces, but less is known on the role of stochastic or historical processes in maintaining variation. The common wall lizard (Podarcis muralis) is a colour polymorphic species in which local colour morph frequencies are thought to be modulated by natural and sexual selection. Here, we used genome-wide single-nucleotide polymorphism data to investigate the relationships between morph composition and population biogeography at a regional scale, by comparing morph composition with patterns of genetic variation of 54 populations sampled across the Pyrenees. We found that genetic divergence was explained by geographic distance but not by environmental features. Differences in morph composition were associated with genetic and environmental differentiation, as well as differences in sex ratio. Thus, variation in colour morph frequencies could have arisen via historical events and/or differences in the permeability to gene flow, possibly shaped by the complex topography and environment. In agreement with this hypothesis, colour morph diversity was positively correlated with genetic diversity and rates of gene flow and inversely correlated with the likelihood of the occurrence of bottlenecks. Concurrently, we did not find conclusive evidence for selection in the two colour loci. As an illustration of these effects, we observed that populations with higher proportions of the rarer yellow and yellow-orange morphs had higher genetic diversity. Our results suggest that processes involving a decay in overall genetic diversity, such as reduced gene flow and/or bottleneck events have an important role in shaping population-specific morph composition via non-selective processes.

Exploring the Distribution of Phylogenetic Networks Generated Under a Birth-Death-Hybridization Process

Gene-flow processes such as hybridization and introgression play important roles in shaping diversity across the tree of life. Recent studies extending birth-death models have made it possible to investigate patterns of reticulation in a macroevolutionary context. These models allow for different macroevolutionary patterns of gene flow events that can either add, maintain, or remove lineages—with the gene flow itself possibly being dependent on the relatedness between species—thus creating complex diversification scenarios. Further, many reticulate phylogenetic inference methods assume specific reticulation structures or phylogenies belonging to certain network classes. However, the distributions of phylogenetic networks under reticulate birth-death processes are poorly characterized, and it is unknown whether they violate common methodological assumptions. We use simulation techniques to explore phylogenetic network space under a birth-death-hybridization process where the hybridization rate can have a linear dependence on genetic distance. Specifically, we measured the number of lineages through time and role of hybridization in diversification along with the proportion of phylogenetic networks that belong to commonly used network classes (e.g., tree-child, tree-based, or level-1 networks). We find that the growth of phylogenetic networks and class membership are largely affected by assumptions about macroevolutionary patterns of gene flow. In accordance with previous studies, a lower proportion of networks belonged to these classes based on type and density of reticulate events. However, under a birth-death-hybridization process, these factors form an antagonistic relationship; the type of reticulation events that cause high membership proportions also lead to the highest reticulation density, consequently lowering the overall proportion of phylogenies in some classes. Further, we observed that genetic distance–dependent gene flow and incomplete sampling increase the proportion of class membership, primarily due to having fewer reticulate events. Our results can inform studies if their biological expectations of gene flow are associated with evolutionary histories that satisfy the assumptions of current methodology and aid in finding phylogenetic classes that are relevant for methods development.

Impact of persistent barrier to gene flow and catastrophic events on red algae evolutionary history along the Chilean coast.

Historical vicariance events, linked to the existence of stable physical barriers to gene flow, generate concordant genetic breaks in co-distributed species while stochastic processes (e.g., costal uplift) could cause species-specific genetic breaks as a result of local strong demographic bottlenecks or extinction. In Chile, previous studies show that the area of the 30°S-33°S could correspond to a stable barrier to gene flow that have affected the genetic structure of various algae and marine invertebrates. Here we sequenced two organellar genes (COI and rbcL) in four taxonomically accepted co-distributed red seaweeds species characterized by a low dispersal potential: Mazzaella laminarioides, M. membranacea, Asterfilopsis disciplinalis, and Ahnfeltiopsis vermicularis. Our results revealed the existence of ten strongly differentiated linages in the taxa studied. Strong genetic breaks, concordant in both space and time (divergence estimated to have occurred some 2.9-12.4 million years ago), were observed between taxa distributed across the 33°S. Conversely, in the Central/South part of the Chilean coast, the localization of the genetic breaks/sub-structure observed varied widely (36°S, 38°S, 39°S, and 40°S). These results suggest that a major historical vicariance event has modeled the genetic structure of several Chilean marine organisms in the north of the Chilean coast during the mid-Miocene, while more recent stochastic events and genetic drift could be the driving forces of genetic divergence/structuration in the central-southern part of the coast.

The contribution of gene flow, selection, and genetic drift to five thousand years of human allele frequency change

Genomic time series from experimental evolution studies and ancient DNA datasets offer us a chance to directly observe the interplay of various evolutionary forces. We show how the genome-wide variance in allele frequency change between two time points can be decomposed into the contributions of gene flow, genetic drift, and linked selection. In closed populations, the contribution of linked selection is identifiable because it creates covariances between time intervals, and genetic drift does not. However, repeated gene flow between populations can also produce directionality in allele frequency change, creating covariances. We show how to accurately separate the fraction of variance in allele frequency change due to admixture and linked selection in a population receiving gene flow. We use two human ancient DNA datasets, spanning around 5,000 y, as time transects to quantify the contributions to the genome-wide variance in allele frequency change. We find that a large fraction of genome-wide change is due to gene flow. In both cases, after correcting for known major gene flow events, we do not observe a signal of genome-wide linked selection. Thus despite the known role of selection in shaping long-term polymorphism levels, and an increasing number of examples of strong selection on single loci and polygenic scores from ancient DNA, it appears to be gene flow and drift, and not selection, that are the main determinants of recent genome-wide allele frequency change. Our approach should be applicable to the growing number of contemporary and ancient temporal population genomics datasets.

Comparative genomic analyses provide new insights into evolutionary history and conservation genomics of gorillas

Genome sequencing is a powerful tool to understand species evolutionary history, uncover genes under selection, which could be informative of local adaptation, and infer measures of genetic diversity, inbreeding and mutational load that could be used to inform conservation efforts. Gorillas, critically endangered primates, have received considerable attention and with the recently sequenced Bwindi mountain gorilla population, genomic data is now available from all gorilla subspecies and both mountain gorilla populations. Here, we reanalysed this rich dataset with a focus on evolutionary history, local adaptation and genomic parameters relevant for conservation. We estimate a recent split between western and eastern gorillas of 150,000–180,000 years ago, with gene flow around 20,000 years ago, primarily between the Cross River and Grauer’s gorilla subspecies. This gene flow event likely obscures evolutionary relationships within eastern gorillas: after excluding putatively introgressed genomic regions, we uncover a sister relationship between Virunga mountain gorillas and Grauer’s gorillas to the exclusion of Bwindi mountain gorillas. This makes mountain gorillas paraphyletic. Eastern gorillas are less genetically diverse and more inbred than western gorillas, yet we detected lower genetic load in the eastern species. Analyses of indels fit remarkably well with differences in genetic diversity across gorilla taxa as recovered with nucleotide diversity measures. We also identified genes under selection and unique gene variants specific for each gorilla subspecies, encoding, among others, traits involved in immunity, diet, muscular development, hair morphology and behavior. The presence of this functional variation suggests that the subspecies may be locally adapted. In conclusion, using extensive genomic resources we provide a comprehensive overview of gorilla genomic diversity, including a so-far understudied Bwindi mountain gorilla population, identify putative genes involved in local adaptation, and detect population-specific gene flow across gorilla species.

Phylogenetic Conflict Between Species Tree and Maternally Inherited Gene Trees in a Clade of Emberiza Buntings (Aves: Emberizidae).

Different genomic regions may reflect conflicting phylogenetic topologies primarily due to incomplete lineage sorting and/or gene flow. Genomic data are necessary to reconstruct the true species tree and explore potential causes of phylogenetic conflict. Here, we investigate the phylogenetic relationships of 4 Emberiza species (Aves: Emberizidae) and discuss the potential causes of the observed mitochondrial non-monophyly of Emberiza godlewskii (Godlewski's bunting) using phylogenomic analyses based on whole genome resequencing data from 41 birds. Analyses based on both the whole mitochondrial genome and ~39 kilobases from the non-recombining W chromosome reveal sister relationships between each the northern and southern populations of E. godlewskii with E. cioides and E. cia, respectively. In contrast, the monophyly of E. godlewskii is reflected by the phylogenetic signal of autosomal and Z chromosomal sequence data as well as demographic inference analyses, which-in combination-support the following tree topology: ([{E. godlewskii, E. cia}, E. cioides], E. jankowskii). Using D-statistics, we detected multiple gene flow events among different lineages, indicating pervasive introgressive hybridization within this clade. Introgression from an unsampled lineage that is sister to E. cioides or introgression from an unsampled mitochondrial + W chromosomal lineage of E. cioides into northern E. godlewskii may explain the phylogenetic conflict between the species tree estimated from genome-wide data versus mtDNA/W tree topologies. These results underscore the importance of using genomic data for phylogenetic reconstruction and species delimitation.

Slendr: a framework for spatio-temporal population genomic simulations on geographic landscapes.

One of the goals of population genetics is to understand how evolutionary forces shape patterns of genetic variation over time. However, because populations evolve across both time and space, most evolutionary processes also have an important spatial component, acting through phenomena such as isolation by distance, local mate choice, or uneven distribution of resources. This spatial dimension is often neglected, partly due to the lack of tools specifically designed for building and evaluating complex spatio-temporal population genetic models. To address this methodological gap, we present a new framework for simulating spatially-explicit genomic data, implemented in a new R package called slendr (www.slendr.net), which leverages a SLiM simulation back-end script bundled with the package. With this framework, the users can programmatically and visually encode spatial population ranges and their temporal dynamics (i.e., population displacements, expansions, and contractions) either on real Earth landscapes or on abstract custom maps, and schedule splits and gene-flow events between populations using a straightforward declarative language. Additionally, slendr can simulate data from traditional, non-spatial models, either with SLiM or using an alternative built-in coalescent msprime back end. Together with its R-idiomatic interface to the tskit library for tree-sequence processing and analysis, slendr opens up the possibility of performing efficient, reproducible simulations of spatio-temporal genomic data entirely within the R environment, leveraging its wealth of libraries for geospatial data analysis, statistics, and visualization. Here, we present the design of the slendr R package and demonstrate its features on several practical example workflows.

Phylogenomics and divergence pattern of Polygonatum (Asparagaceae: Polygonateae) in the north temperate region

Polygonatum is the largest genus of tribe Polygonateae (Asparagaceae) and is widely distributed in the temperate Northern Hemisphere, especially well diversified in southwestern China to northeastern Asia. Phylogenetic relationships of many species are still controversial. Hence it is necessary to clarify their phylogenetic relationships and infer possible reticulate relationships for the genus. In this study, genome-wide data of 43 species from Polygonatum and its closely related taxa were obtained by Hyb-Seq sequencing. The phylogenetic trees constructed from genome-wide nuclear and chloroplast sequences strongly supported the monophyly of Polygonatum with division into three major clades. A high level of incongruence was detected between nuclear and chloroplast trees as well as among gene trees within the genus, but all occurred within each major clade. However, introgression tests and reticulate evolution analyses revealed low level of gene flow and weak introgression events in the genus, suggesting hybridization and introgression were not dominant during the evolutionary diversification of Polygonatum in the Northern Hemisphere. This study provides important insights into reconstructing evolutionary relationships and speciation pattern of taxa from the north temperate flora.

Investigation of phylogenetic relationships within Saxifraga diversifolia complex (Saxifragaceae) based on restriction-site associated DNA sequence markers.

Subsect. Hirculoideae Engl. & Irmsch., belonging to Saxifraga sect. Ciliatae Haw., has high species richness. It can be divided into S. diversifolia, S. pseudohirculus, and S. sinomontana complexes based on morphological characteristics. The species with prominent leaf veins on the posterior leaf edge were placed in the S. diversifolia complex, which is mainly distributed on the eastern and southern margins of the Qinghai-Tibetan Plateau. In this study, 53 samples, representing 15 of the 33 described species in the S. diversifolia complex, were sequenced using the Restriction-site Associated DNA Sequence (RAD-seq) technique. A total of 111,938 high-quality SNP loci were screened to investigate the phylogenetic relationships within the S. diversifolia complex. The result of the neighbor-joining (NJ) tree shows that the S. diversifolia complex is a paraphyletic group. Despite of some inconsistencies as revealed by genetic structural analysis, clustering results of representative species reconstructed by both NJ and principal component analysis analyses support previous biogeographic and morphological evidences. In addition, long-distance gene flow events for 11 taxa were detected in the S. diversifolia complex, respectively from S. implicans 1 to S. implicans 2, S. diversifolia and S. maxionggouensis, and from S. maxionggouensis to S. nigroglandulifera. These findings may improve our comprehension of the phylogeny, classification, and evolution of the S. diversifolia complex.

First evidence of gonadal hybrid dysgenesis in Chagas disease vectors (Hemiptera, Triatominae): gonad atrophy prevents events of interspecific gene flow and introgression

BackgroundHybridization events between Triatoma spp. have been observed under both natural and laboratory conditions. The ability to produce hybrids can influence different aspects of the parent species, and may even result in events of introgression, speciation and extinction. Hybrid sterility is caused by unviable gametes (due to errors in chromosomal pairing [meiosis]) or by gonadal dysgenesis (GD). All of the triatomine hybrids analyzed so far have not presented GD. We describe here for the first time GD events in triatomine hybrids and highlight these taxonomic and evolutionary implications of these events.MethodsReciprocal experimental crosses were performed between Triatoma longipennis and Triatoma mopan. Intercrosses were also performed between the hybrids, and backcrosses were performed between the hybrids and the parent species. In addition, morphological and cytological analyzes were performed on the atrophied gonads of the hybrids.ResultsHybrids were obtained only for the crosses T. mopan♀ × T. longipennis♂. Intercrosses and backcrosses did not result in offspring. Morphological analyses of the male gonads of the hybrids confirmed that the phenomenon that resulted in sterility of the hybrid was bilateral GD (the gonads of the hybrids were completely atrophied). Cytological analyses of the testes of the hybrids also confirmed GD, with no germ cells observed (only somatic cells, which make up the peritoneal sheath).ConclusionsThe observations made during this study allowed us to characterize, for the first time, GD in triatomines and demonstrated that gametogenesis does not occur in atrophied gonads. The characterization of GD in male hybrids resulting from the crossing of T. mopan♀ × T. longipennis♂ highlights the importance of evaluating both the morphology and the cytology of the gonads to confirm which event resulted in the sterility of the hybrid: GD (which results in no gamete production) or meiotic errors (which results in non-viable gametes).Graphical

Episodic population fragmentation and gene flow reveal a trade-off between heterozygosity and allelic richness.

In episodic environments like deserts, populations of some animal species exhibit irregular fluctuations such that populations are alternately large and connected or small and isolated. Such dynamics are typically driven by periodic resource pulses due, for example, to large but infrequent rainfall events. The repeated population bottlenecks resulting from fragmentation should lower genetic diversity over time, yet species undergoing these fluctuations appear to maintain high levels of genetic diversity. To resolve this apparent paradox, we simulated a metapopulation of constant size undergoing repeat episodes of fragmentation and change in gene flow to mimic outcomes experienced by mammals in an Australian desert. We show that episodic fragmentation and gene flow have contrasting effects on two measures of genetic diversity: heterozygosity and allelic richness. Specifically, fragmentation into many, small subpopulations, coupled with periods of infrequent gene flow, preserves allelic richness at the expense of heterozygosity. In contrast, fragmentation into a few, large subpopulations maintains heterozygosity at the expense of allelic richness. The strength of the trade-off between heterozygosity and allelic richness depends on the amount of gene flow and the frequency of gene flow events. Our results imply that the type of genetic diversity maintained among species living in strongly fluctuating environments will depend on the way populations fragment, with our results highlighting different mechanisms for maintaining allelic richness and heterozygosity in small, fragmented populations.

Comprehensive evaluation of germplasm resources of nine Macrobrachium rosenbergii strains in China

Macrobrachium rosenbergii, as an introduced species, is extensively cultured in China. In this study, 176 individuals from nine domestic strains in China, i.e., Florida (FD), Guangxi (GX), Hefu (HF), Myanmar (MD), “Nantaihu No. 2″ (NT), Thailand (TG), Taiwan (TW), “Shufeng No. 1″ (SF), and Zhengda (ZD), were collectedfor genome resequencing. We evaluated the genetic diversity, population structure and population history of nine cultured strains using SNP markers. We found that genetic diversity of the M. rosenbergii strains was at a low level, among which FD and TG strains showed relatively higher genetic diversity, inferring that the genetic diversity of domestic populations decreased after artificial selection. Population genetic structure analysis showed that the individuals GX, HF, SF and NT strains were irregularly clustered together, and the ZD and TW strains exhibited the specific ancestors compared with other strains. The phylogenetic tree revealed that the individuals from TW, MD, ZD strains grouped into population-specific clusters. The principal components analysis (PCA) also showed that the nine strains were mixed together except the individuals from ZD and TW strains which were separated into specific subgroups. The Fst analysis showed that the Fst values detected from nine strains were less than 0.05, with the lowest Fst observed between SF and NT strains. In addition, two positively selected genes involved in metal transport and stress response were identified in the nine domestic strains by using selective sweep analyses with composite likelihood ratio (CLR), Fst, Pi, and Tajima’ D methods. The TreeMix phylogeny inferred four gene flow events among the nine strains, which occurred in the following directions: (1) from the common ancestor of ZD and TG to SF, (2) from SF to MD, (3) from SF to NT, and (4) from HF to GX. This study contributes to the better understanding genetic background and germplasm status of domestic M. rosenbergii strains in China.

Genomic signatures reveal selection in Chinese and European domesticated geese.

The Swan goose and Greylag goose are species of geese native to East Asia and Europe, respectively, and are widely believed to be the ancestors of Chinese and European domesticated geese. The Yili goose (YL) and European domestic geese originated from the Greylag goose, but the history of domestication is unclear. In this study, we sequenced and analyzed the genome of the YL goose and the Hortobagy goose to combine with other previously sequenced goose populations for in-depth analysis. The population genetic variations in Stone geese, East Zhejiang White Geese, Taihu geese and Zi geese were also identified and compared. The results showed that admixture gene flow existed in the YL geese population, which was introgressed by Chinese geese, suggesting that gene flow events were frequent and widespread among domesticated geese. Further selected sweep analysis identified candidate genes and metabolic pathways that may be related to the differences in morphology. Several genes such as TGFBR3L, CMYA5, FOXD1, ARHGEF28 and SUCLG2 are associated with growth, reproduction and fertility traits. The results of this study will help to understand the genetic characteristics of domestic geese and the genes affecting important traits and provide a basis for the improved breed of domestic geese.

Phylogeographic relationships and morphological evolution between cave and surface Astyanax mexicanus populations (De Filippi 1853) (Actinopterygii, Characidae).

The Astyanax mexicanus complex includes two different morphs, a surface- and a cave-adapted ecotype, found at three mountain ranges in Northeastern Mexico: Sierra de El Abra, Sierra de Guatemala and Sierra de la Colmena (Micos). Since their discovery, multiple studies have attempted to characterize the timing and the number of events that gave rise to the evolution of these cave-adapted ecotypes. Here, using RADseq and genome-wide sequencing, we assessed the phylogenetic relationships, genetic structure and gene flow events between the cave and surface Astyanax mexicanus populations, to estimate the tempo and mode of evolution of the cave-adapted ecotypes. We also evaluated the body shape evolution across different cave lineages using geometric morphometrics to examine the role of phylogenetic signal versus environmental pressures. We found strong evidence of parallel evolution of cave-adapted ecotypes derived from two separate lineages of surface fish and hypothesize that there may be up to four independent invasions of caves from surface fish. Moreover, a strong congruence between the genetic structure and geographic distribution was observed across the cave populations, with the Sierra de Guatemala the region exhibiting most genetic drift among the cave populations analysed. Interestingly, we found no evidence of phylogenetic signal in body shape evolution, but we found support for parallel evolution in body shape across independent cave lineages, with cavefish from the Sierra de El Abra reflecting the most divergent morphology relative to surface and other cavefish populations.