Model Of Gene Flow Research Articles

Understanding admixture fractions: theory and estimation of gene-flow.

Estimation of admixture proportions has become one of the most commonly used computational tools in population genomics. However, there is remarkably little population genetic theory on statistical properties of these variables. We develop theoretical results that can accurately predict means and variances of admixture proportions within a population using models with recombination and genetic drift. Based on established theory on measures of multilocus disequilibrium, we show that there is a set of recurrence relations that can be used to derive expectations for higher moments of the admixture proportions distribution. We obtain closed form solutions for some special cases. Using these results, we develop a method for estimating admixture parameters from estimated admixture proportions obtained from programs such as Structure or Admixture. We apply this method to HapMap 3 data and find that the population history of African Americans, as expected, is not best explained by a single admixture event between people of European and African ancestry. The model of constant gene flow starting at 8 generations and ending at 2 generations before present gives the best fit.

Islands of ice: Glacier‐dwelling metazoans form regionally distinct populations despite extensive periods of deglaciation

AbstractAimGlaciers cover considerable portion of land and host diverse life forms from single‐celled organisms to invertebrates. However, the determinants of diversity and community composition of these organisms remain underexplored. This study addresses the biogeography, population connectivity and dispersal of these organisms, especially critical in understanding during the rapid recession of glaciers and increased extinction risk for isolated populations. By reconstructing the Quaternary biogeographic history of Fontourion glacialis, a widespread in Northern Hemisphere glacier obligate species of Tardigrada, we aim to understand how populations of glacier‐dwelling metazoans receive immigrants, respond to disappearing glaciers and to what extent remaining glaciers can serve as refugia.LocationGlaciers across Svalbard, Scandinavia, Greenland and Iceland.MethodsWe analysed mtDNA (COI gene) variability of 263 F. glacialis specimens collected across the distribution range. Phylogeographic and coalescent‐based approaches were used to detect population differentiation patterns, investigate most likely models of gene flow and test the influences of geographical and climatic factors on the distribution of F. glacialis genetic variants.ResultsOur findings indicate that the distribution of F. glacialis genetic variants is primarily influenced by geographical rather than climatic factors. Populations exhibit a dispersal‐limited distribution pattern, influenced by geographical distance and local barriers, even between neighbouring glaciers. Significantly, the genetic structure within Scandinavia suggests the existence of “southern” glacial or low‐temperature refugia, where F. glacialis may have survived a period of extensive deglaciation during the Holocene climatic optimum (8–5 kyr ago).Main ConclusionThe study uncovers complex metapopulation structures in F. glacialis, with impacts of local barriers, population bottlenecks as well as historical ice sheet fluctuations. It suggests that such populations can endure extended periods of deglaciation, highlighting the resilience of glacial refugia. The study highlights the necessity of understanding the diversity and population structure of ice‐dwelling fauna in both spatial and temporal contexts.

Maintenance of a narrow hybrid zone between native and introduced red foxes (Vulpes vulpes) despite conspecificity and high dispersal capabilities.

Human-facilitated introductions of nonnative populations can lead to secondary contact between allopatric lineages, resulting in lineage homogenisation or the formation of stable hybrid zones maintained by reproductive barriers. We investigated patterns of gene flow between the native Sacramento Valley red fox (Vulpes vulpes patwin) and introduced conspecifics of captive-bred origin in California's Central Valley. Considering their recent divergence (20-70 kya), we hypothesised that any observed barriers to gene flow were primarily driven by pre-zygotic (e.g. behavioural differences) rather than post-zygotic (e.g. reduced hybrid fitness) barriers. We also explored whether nonnative genes could confer higher fitness in the human-dominated landscape resulting in selective introgression into the native population. Genetic analysis of red foxes (n = 682) at both mitochondrial (cytochrome b + D-loop) and nuclear (19,051 SNPs) loci revealed narrower cline widths than expected under a simulated model of unrestricted gene flow, consistent with the existence of reproductive barriers. We identified several loci with reduced introgression that were previously linked to behavioural divergence in captive-bred and domestic canids, supporting pre-zygotic, yet possibly hereditary, barriers as a mechanism driving the narrowness and stability of the hybrid zone. Several loci with elevated gene flow from the nonnative into the native population were linked to genes associated with domestication and adaptation to human-dominated landscapes. This study contributes to our understanding of hybridisation dynamics in vertebrates, particularly in the context of species introductions and landscape changes, underscoring the importance of considering how multiple mechanisms may be maintaining lineages at the species and subspecies level.

Historical fragmentation and stepping-stone gene flow led to population genetic differentiation in a coastal seabird.

Understanding the forces that shape population genetic structure is fundamental both for understanding evolutionary trajectories and for conservation. Many factors can influence the geographic distribution of genetic variation, and the extent to which local populations differ can be especially difficult to predict in highly mobile organisms. For example, many species of seabirds are essentially panmictic, but some show strong structure. Pigeon Guillemots (Cepphus columba; Charadriiformes: Alcidae) breed in small colonies scattered along the North Pacific coastline and feed in shallow nearshore waters year-round. Given their distribution, gene flow is potentially lower and population genetic structure is stronger than in most other high-latitude Northern Hemisphere seabirds. We screened variation in the mitochondrial control region, four microsatellite loci, and two nuclear introns in 202 Pigeon Guillemots representing three of five subspecies. Mitochondrial sequences and nuclear loci both showed significant population differences, although structure was weaker for the nuclear loci. Genetic differentiation was correlated with geographic distance between sampling locations for both the mitochondrial and nuclear loci. Mitochondrial gene trees and demographic modeling both provided strong evidence for two refugial populations during the Pleistocene glaciations: one in the Aleutian Islands and one farther east and south. We conclude that historical fragmentation combined with a stepping-stone model of gene flow led to the relatively strong population differentiation in Pigeon Guillemots compared to other high-latitude Northern Hemisphere seabird species. Our study adds to growing evidence that Pleistocene glaciation events affected population genetic structure not only in terrestrial species but also in coastal marine animals.

Modelling jaguar gene flow in fragmented landscapes offers insights into functional population connectivity

ContextPreserving functional connectivity is a conservation priority to secure the long-term viability of geographically dispersed subpopulations, such as the jaguar (Panthera onca) populations in Central America. However, managing connectivity in this region is limited due to the scarcity of local assessments of the connectivity between existing populations, some of which exhibit low levels of gene flow and genetic admixture.ObjectivesWe selected the jaguar as a model species to understand how gene flow of large carnivores is shaped in a heavily human-impacted landscape in the Neotropics. We assessed the impact of landscape features and spatial scale on jaguar gene flow across northern Central America, predicted synoptic, landscape-wide functional connectivity across the region; and compared connectivity predictions based on models of gene flow, habitat suitability, and expert knowledge.MethodsWe employed genetic data based on 335 faecal samples collected from 72 individual jaguars. We parameterized gene flow resistance surfaces using linear mixed effects models and the maximum likelihood population-effects method. We evaluated nine landscape variables at six spatial scales, selecting the optimal scale and transformation for each variable according to univariate models and AIC. To predict jaguar gene flow, we developed multivariate models and implemented resistant kernels to forecast functional connectivity between jaguar populations across the study region under three dispersal distance scenarios. Furthermore, we compared the connectivity estimates based on gene flow against those based on habitat suitability and the corridors delineated by expert knowledge in the region.ResultsLow resistance to jaguar gene flow was associated with greater tree cover and vegetation, lower areas of built-up, and intermediate distances from water bodies. Notably, tree cover affected jaguar gene flow on a smaller scale compared to the rest of variables. Higher connectivity, indicated by lower resistance, was found within and around the largest jaguar conservation units (JCUs) such as the Reserva-de-Biosfera-Transfronteriza, Selva-Maya and Maya-Mountains. This contrasted with the smallest JCUs, like Sierra-Santa-Cruz and Sierra-de-las-Minas in eastern Guatemala, and Cordillera-Nombre-de-Dios in Honduras. Across the region, lower connectivity was observed in the Caribbean connection between eastern Guatemala and midwestern Honduras, as well as in the Honduran and Nicaraguan Miskito area. Models based on gene flow and habitat suitability were similar in their predictions of areas of high connectivity; however, the habitat suitability models predicted larger areas of low connectivity than did the gene flow models. Moreover, the expert knowledge corridors were consistent with areas of high and medium connectivity as predicted by the gene flow model.ConclusionsGene flow of jaguars is positively driven by the presence of forest and water sources, while human impact (built-up and non-vegetation areas) has a negative effect. Areas of lowest resistance largely correspond to the location of JCUs, which serve as crucial reservoirs of high-quality jaguar habitat within the region. While the largest JCU in the region displayed high connectivity, the low connectivity between the smallest JCUs (i.e. Sierra-Santa-Cruz, Sierra-de-las-Minas and Cordillera-Nombre-de-Dios) underscores the need for conservation attention in these areas. Conservation and management actions such as habitat loss-prevention/restoration and anthropogenic impact mitigation should be prioritized in the binational region of Guatemala-Honduras, a key connectivity bottleneck between the species’ northern and southern ranges. Similarly, attention is warranted in the Honduran-Nicaraguan Miskito area.

Molecular discrimination of G1 and G3 genotypes of Echinococcus granulosus sensu stricto obtained from human, cattle, and sheep using the mitochondrial NADH dehydrogenase subunit 5 marker

Cystic echinococcosis (CE) is a common zoonotic disease caused by the larval form of Echinococcus granulosus sensu lato. This study determined the genotype and haplotype differences using the NADH dehydrogenase subunit 5 gene in hydatid cyst samples. Human (n = 12), cattle (n = 28), and sheep (n = 31) hydatid cyst isolates were included. Seventy-one genomic DNA samples were successfully extracted, and a 759 bp mitochondrial NADH dehydrogenase subunit 5 gene fragment was amplified by PCR. Following the sequence analysis, E. granulosus sensu stricto isolates were identified as G1 (n = 61) and G3 (n = 10). A total of 23 haplotypes were obtained from the 71 E. granulosus s.s. G1 and G3 samples. The main haplotype was Hap01 (60.56 %), which consisted of the G1 genotype. The second largest haplotype was Hap04, which consisted entirely of the G3 genotype. Hap14 acted as a bridge between the G1 and G3 genotypes. This study identifies G1 as the dominant genotype in humans and farm animals in Turkey. High haplotype and nucleotide diversity in genotypes were observed. Additionally, this is the first report on the phylogeography and gene flow models of the E. granulosus s.s. population in Turkey using the NADH dehydrogenase subunit 5 gene, the best marker distinguishing between G1 and G3 genotypes.

Genetic diversity, demography and connectivity in the marine snail Thaisella chocolata (Duclos, 1832): implications for sustainable fisheries management along the Peruvian coast

ABSTRACT Thaisella chocolata is a marine gastropod (family Muricidae) occurring from northern Peru to central Chile. It is of economic importance in these countries, where hundreds of tons are caught every year. In combination with El Niño–Southern Oscillation events, fishing activity has led to the collapse of Chilean populations, raising concerns about the state of Peruvian populations. The aim of this study is to assess the demography, genetic diversity and connectivity of T. chocolata along the Peruvian coast. This would provide important insights into the resilience of this species to anthropogenic disturbance. In this study, four locations along the coast of Peru were investigated for the first time. Mitochondrial cytochrome c oxidase subunit I (COI) sequences were obtained from 116 individuals and four microsatellites from 95 individuals. Both COI (combined with 156 publicly available sequences) and microsatellite data revealed very low genetic differentiation, despite a few significant φST and FST values between population pairs. The most likely gene flow model (MIGRATE-N) is panmixia between all locations, with gene flow occurring equally northwards and southwards. This pattern may be related to the surface currents along the Peruvian coastline combined with the high pelagic larval duration of T. chocolata. Hyperdiversity was found for the COI marker (πsyn > 5%). Microsatellite loci were also polymorphic, with a mean diversity of 12.69 alleles per locus per population. No recent population bottleneck was observed (microsatellite data), while at least one past demographic expansion was suggested (COI data). Given the high genetic diversity and connectivity, T. chocolata is expected to show high evolutionary resilience and can be considered as one management unit in terms of fisheries management. However, further studies, on, for example, the characterization of effective population size (using numerous nuclear loci), are required to confirm this hypothesis. In addition, precautionary fisheries measures should be maintained, given the recent collapse of T. chocolata stocks in Chile.

Highland islands: assessing gene flow among populations of an aquatic highland endemic

Abstract A variety of hierarchical gene flow models have been proposed to explain the distribution of genetic variation in aquatic environments. These models, including the Stream Hierarchy (SH), Death Valley (DV), Headwater (HW), and Widespread Gene flow (WG) models, provide testable hypotheses that focus on the degree of within-basin or within-stream network connectivity and an organism’s dispersal abilities. We tested these models using Tlaloc hildebrandi (Profundulidae, Cyprinodontiformes), a freshwater fish endemic to the highlands of the Grijalva and Usumacinta River basins in southern Mexico. Data from ultraconserved elements (UCEs) showed that although T. hildebrandi was recovered as monophyletic, the sub-basins were not. We generated single nucleotide polymorphisms from the UCEs to analyse the data in a population genetic framework. These results differed between analyses such that two (STRUCTURE analysis) or three (DAPC analysis) genetic groupings were recovered. Overall, the results of this study provide support for the SH model. Some individuals from the Jataté sub-basin, however, conformed to the HW model, due to historical connections among headwater streams and rivers of the Jataté and Amarillo–Chenalhó clades/clusters. The greatest degree of gene flow has occurred from the Grijlava to the Usumacinta basins, two geographically proximate basins that have been hypothesized to have had previous connections.

Drones Do Not Drift between Nests in a Wild Population of Apis cerana.

The modes through which individuals disperse prior to reproduction has important consequences for gene flow in populations. In honey bees (Apis sp.), drones (males) reproduce within a short flight range of their natal nest, leaving and returning each afternoon within a narrow mating window. Drones are assumed to return to their natal nests as they depend on workers to feed them. However, in apiaries, drones are reported to regularly make navigation errors and return to a non-natal nest, where they are accepted and fed by unrelated workers. If such a "drone drift" occurred in wild populations, it could facilitate some further degree of dispersal for males, particularly if drones drift into host nests some distance away from their natal nest. Here, we investigated whether drone drift occurs in an invasive population of the Asian honey bee (Apis cerana). Based on the genotypes of 1462 drones from 19 colonies, we found only a single drone that could be considered a candidate drifter (~0.07%). In three other colonies, drones whose genotypes differed from the inferred queen were best explained by recent queen turnover or worker-laying. We concluded that drone drift in this population is low at best, and A. cerana drones either rarely make navigation errors in wild populations or are not accepted into foreign nests when they do so. We therefore confirm that drone dispersal distance is limited to the distance of daily drone flights from natal nests, a key assumption of both colony density estimates based on sampling of drone congregation areas and population genetic models of gene flow in honey bees.

Population structure of mud flounder Paralichthys orbignyanus from the south-western Atlantic Ocean.

The mud flounder Paralichthys orbignyanus (Pleuronectiformes, Paralichthyidae) inhabits shallow waters of low salinities and mud bottoms in the temperate marine coastal regions of the Bonaerensean Ecoregion of the Argentinean Biogeographic Province in the south-western Atlantic Ocean. Specimens of P. orbignyanus were collected from Lagoa dos Patos (LDP) (southern Brazil), Mar Chiquita (MCH) and Marisol (MAR) both located in Buenos Aires (Argentina), and San Antonio Oeste (SAO) in the San Matías Gulf, Rio Negro (Argentina). A fragment of the mitochondrial DNA of the Control Region and seven microsatellite loci were characterized. In the Control Region, P. orbignyanus showed high variability, low nucleotide diversity, mild population expansion and a coalescence time of 35,000 years before the present. Flounders provided evidence of a genetic structure between the sampling sites LDP, MCH, MAR vs. SAO. On the other hand, P. orbignyanus displayed a lower to moderate contemporary genetic structure among all samples except between LDP and MCH. With no evidence of isolation by distance, this analysis supports a model of limited gene flow that is likely to be associated with a consistent larvae retention in all sampling sites. In addition, the present connectivity is ascribed to a lower migration process from SAO in the San Matías Gulf congruent with the prevailing littoral drift.

Ice Age effects on genetic divergence of the American crocodile (Crocodylus acutus) in Panama: reconstructing limits of gene flow and environmental ranges: a reply to O'Dea et al.

O'Dea et al. (2022) (Pleistocene sea level changes and crocodile population histories on the isthmus of panama: a comment on Avila-Cervantes et al. (2020). Evolution, 76(11), 2778-2783. https://doi.org/10.1111/evo.14610) question our hypothesis that sea-level changes during the past glaciation played a role in restricting gene flow between Pacific and Caribbean Crocodylus acutus in Panama. They argue that an error in sea-level high-stand reconstruction during the last interglacial period (118-130 ka) does not support our hypothesis. Although they are correct in our high-stand reconstruction error, overlooked the point in that we were presenting a model of restricted gene flow across the Panamanian Isthmus during low sea levels. We review the assumptions of gene demographic methods, emphasizing that we were focusing on times of genetic divergence. We expand here why gene flow between these coastal populations could have been restricted during the last glacial maximum (19-26.5 ka) and the 50,000 years preceding it when sea levels were lower than today. O'Dea et al. suggest local climates may have played larger roles than sea levels. We demonstrate that paleoclimate estimates for the past 3.3 Ma in Panama are within the bounds of extant C. acutus. The importance of Ice Age Sea-level dynamics on Neotropical species was likely profound and should be incorporated into evolutionary studies of these taxa.

Complex cycles of divergence and migration shape lineage structure in the common kingsnake species complex.

The Nearctic is a complex patchwork of habitats and geologic features that form barriers to gene flow resulting in phylogeographic structure and speciation in many lineages. Habitats are rarely stable over geologic time, and the Nearctic has undergone major climatic changes in the past few million years. We use the common kingsnake species complex to study how climate, geography, and history influence lineage formation over a large, complex landscape. Nearctic/North America. Common kingsnake, Lampropeltis getula, species complex. We analyzed genome-wide sequence data from 51 snakes spanning the majority of the species complex's range. We used population clustering, generalized dissimilarity modeling, and coalescent methods to identify the number of genetic clusters within the L. getula complex, infer the environmental correlates of genetic differentiation, and estimate models of divergence and gene flow among lineages. We identified three major lineages within the L. getula complex and further continuous spatial structure within lineages. The most important ecological correlates of genetic distance in the complex are related to aridity and precipitation, consistent with lineage breaks at the Great Plains/Desert ecotone and the Cochise Filter Barrier. Lineages are estimated to have undergone multiple rounds of isolation and secondary contact, with highly asymmetric migration occurring at present. Changing climates combined with a large and geologically complex landscape have resulted in a mosaic of discrete and spatially continuous genetic structure. Multiple rounds of isolation and secondary contact as climate fluctuated over the past ~4.4 My have likely driven the evolution of discrete lineages that maintain high levels of gene flow. Continuous structure is strongly shaped by aridity and precipitation, suggesting roles for major precipitation gradients in helping to maintain lineage identity in the face of gene flow when lineages are in geographic contact.

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.

Speciation Hypotheses from Phylogeographic Delimitation Yield an Integrative Taxonomy for Seal Salamanders (Desmognathus monticola).

Significant advances have been made in species delimitation and numerous methods can test precisely defined models of speciation, though the synthesis of phylogeography and taxonomy is still sometimes incomplete. Emerging consensus treats distinct genealogical clusters in genome-scale data as strong initial evidence of speciation in most cases, a hypothesis that must therefore be falsified under an explicit evolutionary model. We can now test speciation hypotheses linking trait differentiation to specific mechanisms of divergence with increasingly large data sets. Integrative taxonomy can, therefore, reflect an understanding of how each axis of variation relates to underlying speciation processes, with nomenclature for distinct evolutionary lineages. We illustrate this approach here with Seal Salamanders (Desmognathus monticola) and introduce a new unsupervised machine-learning approach for species delimitation. Plethodontid salamanders are renowned for their morphological conservatism despite extensive phylogeographic divergence. We discover 2 geographic genetic clusters, for which demographic and spatial models of ecology and gene flow provide robust support for ecogeographic speciation despite limited phenotypic divergence. These data are integrated under evolutionary mechanisms (e.g., spatially localized gene flow with reduced migration) and reflected in emergent properties expected under models of reinforcement (e.g., ethological isolation and selection against hybrids). Their genetic divergence is prima facie evidence for species-level distinctiveness, supported by speciation models and divergence along axes such as behavior, geography, and climate that suggest an ecological basis with subsequent reinforcement through prezygotic isolation. As data sets grow more comprehensive, species-delimitation models can be tested, rejected, or corroborated as explicit speciation hypotheses, providing for reciprocal illumination of evolutionary processes and integrative taxonomies. [Desmognathus; integrative taxonomy; machine learning; species delimitation.].

Constrained control of gene-flow models

In ecology and population dynamics, gene flow refers to the transfer of a trait (e.g. genetic material) from one population to another. This phenomenon is of great relevance in studying the spread of diseases or the evolution of social features, such as languages. From the mathematical point of view, gene flow is modeled using bistable reaction–diffusion equations. The unknown is the proportion p of the population that possesses a certain trait, within an overall population N . In such models, gene flow is taken into account by assuming that the population density N depends either on p (if the trait corresponds to fitter individuals) or on the location x (if some zones in the domain can carry more individuals). Recent applications stemming from mosquito-borne-disease control problems or from the study of bilingualism have called for the investigation of the controllability properties of these models. At the mathematical level, this corresponds to boundary control problems and, since we are working with proportions, the control u has to satisfy the constraint 0\leq u \leq 1 . In this article, we provide a thorough analysis of the influence of the gene-flow effect on boundary controllability properties. We prove that, when the population density N only depends on the trait proportion p , the geometry of the domain is the only criterion that has to be considered. We then tackle the case of population densities N varying in x . We first prove that, when N varies slowly in x and when the domain is narrow enough, controllability always holds. This result is proved using a robust domain perturbation method. We then consider the case of sharp fluctuations in N : we first give examples that prove that controllability may fail. Conversely, we give examples of heterogeneities N such that controllability will always be guaranteed: in other words the controllability properties of the equation are very strongly influenced by the variations of N . All negative controllability results are proved by showing the existence of nontrivial stationary states, which act as barriers. The existence of such solutions and the methods of proof are of independent interest. Our article is completed by several numerical experiments that confirm our analysis.

Population genomic structure in Goodman's mouse lemur reveals long-standing separation of Madagascar's Central Highlands and eastern rainforests.

Madagascar's Central Highlands are largely composed of grasslands, interspersed with patches of forest. The historical perspective was that Madagascar's grasslands had anthropogenic origins, but emerging evidence suggests that grasslands were a component of the pre-human Central Highlands vegetation. Consequently, there is now vigorous debate regarding the extent to which these grasslands have expanded due to anthropogenic pressures. Here, we shed light on the temporal dynamics of Madagascar's vegetative composition by conducting a population genomic investigation of Goodman's mouse lemur (Microcebus lehilahytsara; Cheirogaleidae). These small-bodied primates occur both in Madagascar's eastern rainforests and in the Central Highlands, making them a valuable indicator species. Population divergences among forest-dwelling mammals will reflect changes to their habitat, including fragmentation, whereas patterns of post-divergence gene flow can reveal formerly wooded migration corridors. To explore these patterns, we used RADseq data to infer population genetic structure, demographic models of post-divergence gene flow, and population size change through time. The results offer evidence that open habitats are an ancient component of the Central Highlands, and that widespread forest fragmentation occurred naturally during a period of decreased precipitation near the last glacial maximum. Models of gene flow suggest that migration across the Central Highlands has been possible from the Pleistocene through the recent Holocene via riparian corridors. Though our findings support the hypothesis that Central Highland grasslands predate human arrival, we also find evidence for human-mediated population declines. This highlights the extent to which species imminently threatened by human-mediated deforestation may already be vulnerable from paleoclimatic conditions.

Combining Species Delimitation, Species Trees, and Tests for Gene Flow Clarifies Complex Speciation in Scrub-Jays.

Complex speciation, involving rapid divergence and multiple bouts of post-divergence gene flow, can obfuscate phylogenetic relationships and species limits. In North America, cases of complex speciation are common, due at least in part to the cyclical Pleistocene glacial history of the continent. Scrub-Jays in the genus Aphelocoma provide a useful case study in complex speciation because their range throughout North America is structured by phylogeographic barriers with multiple cases of secondary contact between divergent lineages. Here, we show that a comprehensive approach to genomic reconstruction of evolutionary history, i.e., synthesizing results from species delimitation, species tree reconstruction, demographic model testing, and tests for gene flow, is capable of clarifying evolutionary history despite complex speciation. We find concordant evidence across all statistical approaches for the distinctiveness of an endemic southern Mexico lineage (A. w. sumichrasti), culminating in support for the species status of this lineage under any commonly applied species concept. We also find novel genomic evidence for the species status of a Texas endemic lineage A. w. texana, for which equivocal species delimitation results were clarified by demographic modeling and spatially explicit models of gene flow. Finally, we find that complex signatures of both ancient and modern gene flow between the non-sister California Scrub-Jay (Aphelocoma californica) and Woodhouse's Scrub-Jay (Aphelocoma woodhouseii) result in discordant gene trees throughout the species' genomes despite clear support for their overall isolation and species status. In sum, we find that a multi-faceted approach to genomic analysis can increase our understanding of complex speciation histories, even in well-studied groups. Given the emerging recognition that complex speciation is relatively commonplace, the comprehensive framework that we demonstrate for interrogation of species limits and evolutionary history using genomic data can provide a necessary roadmap for disentangling the impacts of gene flow and incomplete lineage sorting (ILS) to better understand the systematics of other groups with similarly complex evolutionary histories. [Gene flow; genomics; introgression; phylogenomics; population genetics; species delimitation; species trees; systematics.].

Long Distance Runners in the Marine Realm: New Insights Into Genetic Diversity, Kin Relationships and Social Fidelity of Indian Ocean Male Sperm Whales

Adult male sperm whales (Physeter macrocephalus) are long distance runners of the marine realm, feeding in high latitudes and mating in tropical and subtropical waters where stable social groups of females and immatures live. Several areas of uncertainty still limit our understanding of their social and breeding behavior, in particular concerning the potential existence of geographical and/or social fidelities. In this study, using underwater observation and sloughed-skin sampling, we looked for male social fidelity to a specific matrilineal sperm whale group near Mauritius. In addition, we captured a wider picture of kin relationships and genetic diversity of male sperm whales in the Indian Ocean thanks to biopsies of eight individuals taken in a feeding ground near the Kerguelen and Crozet Archipelagos (Southern Indian Ocean). Twenty-six adult male sperm whales were identified when socializing with adult females and immatures off Mauritius. Sloughed-skin samples were taken from thirteen of them for genetic analysis. Long-term underwater observation recorded several noteworthy social interactions between adult males and adult females and/or immatures. We identified seven possible male recaptures over different years (three by direct observation, and four at the gametic level), which supports a certain level of male social fidelity. Two probable first- and thirty second-degree kin relationships were highlighted between members of the social unit and adult males, confirming that some of the adult males observed in Mauritian waters are reproductive. Male social philopatry to their natal group can be excluded, as none of the males sampled shared the haplotype characteristic of the matrilineal social group. Mitochondrial DNA control region haplotype and nucleotide diversities calculated over the 21 total male sperm whales sampled were similar to values found by others in the Indian Ocean. Our study strongly supports the existence of some levels of male sperm whale social fidelity, not directed to their social group of birth, in the Indian Ocean. Males sampled in breeding and feeding grounds are linked by kin relationships. Our results support a model of male mediated gene flow occurring at the level of the whole Indian Ocean, likely interconnected with large-scale geographical fidelity to ocean basin, and a small-scale social fidelity to matrilineal social groups.

Application of the maximum threshold distances to reduce gene flow frequency in the coexistence between genetically modified (GM) and non-GM maize.

On the coexistence of genetically modified (GM) and non‐GM maize, the isolation distance plays an important role in controlling the transgenic flow. In this study, maize gene flow model was used to quantify the MTD0.1% and MTD1% in the main maize‐planting regions of China; those were the maximum threshold distance for the gene flow frequency equal to or lower than 1% and 0.1%. The model showed that the extreme MTD1% and MTD0.1% were 187 and 548 m, respectively. The regions of northern China and the coastal plain, including Hainan crop winter‐season multiplication base, showed a significantly high risk for maize gene flow, while the west‐south of China was the largest low‐risk areas. Except for a few sites, the isolation distance of 500 m could yield a seed purity of better than 0.1% and meet the production needs of breeder seeds. The parameters of genetic competitiveness (cp) were introduced to assess the effects of hybrid compatibility between the donor and recipient. The results showed that hybrid incompatibility could minimize the risk. When cp = 0.05, MTD1% and MTD0.1% could be greatly reduced within 19 m and 75 m. These data were helpful to provide scientific data to set the isolation distance between GM and non‐GM maize and select the right place to produce the hybrid maize seeds.

Gene Flow Patterns of the Aedes aegypti (Diptera: Culicidae) Mosquito in Colombia: a Continental Comparison Suggests Multiple Invasion Routes and Gene Exchange

In Colombia Aedes aegypti is present in 80% of the country up to 2,300 m; however, little is known of its genetic relations within a country context and, hence, within a broader context, for example, America. The aforementioned, herein, analyzed the gene flow within a context of the Americas, its directionality and genetic diversity of the mitochondrial lineages in the A. aegypti populations for Colombia. This called for the use of the sequences for A. aegypti available of the mitochondrial ND4 gene in the GenBank for Colombia and the American continent. No presence was found of nuclear mitochondrial pseudogenes (NUMTs) for Colombia. It is estimated that in Colombia the gene flow of the A. aegypti populations is occurring from the southeast and northeast toward the center of the country. In comparison with the mitochondrial sequences for America, the vector’s haplotypes in Colombia suggest connections between the populations of mosquitoes from the south with those from the north of the continent. The gene flow model at continental scale suggests bidirectional connections between the populations from the north of the continent with those from the south, while at South American scale it proposes the gene flow in all the directions with respect to the Colombian. The Colombian A. aegypti vector monitoring and control strategies must pay special attention to the vector’s points of entry into Colombia related with Peru, Venezuela, Brazil, Mexico, and North America to avoid the entry of populations with characteristics like resistance to insecticides or vector competition.