Gene Flow Research Articles

Chromosome-level baobab genome illuminates its evolutionary trajectory and environmental adaptation

Baobab (Adansonia digitata) is a long-lived tree endemic to Africa with economic, ecological, and cultural importance, yet its genomic features are underexplored. Here, we report a chromosome-level reference genome anchored to 42 chromosomes for A. digitata, alongside draft assemblies for a sibling tree, two trees from distinct locations in Africa, and A. za from Madagascar. The baobab genome is uniquely rich in DNA transposons, which make up 33%, while LTR retrotransposons account for 10%. A. digitata experienced whole genome multiplication (WGM) around 30 million years ago (MYA), followed by a second WGM event 3–11 MYA, likely linked to autotetraploidy. Resequencing of 25 trees identify three subpopulations, with gene flow across West Africa distinct from East Africa. Gene enrichment and fixation index (Fst) analyses show baobab retained multiple circadian, flowering, and light-responsive genes, which likely support longevity through the UV RESISTANCE LOCUS 8 (UVR8) pathway. In sum, we provide genomic resources and insights for baobab breeding and conservation.

Genetic Diversity and Population Structure of Endangered Orchid Cypripedium flavum in Fragmented Habitat Using Fluorescent AFLP Markers

Genetic diversity is crucial for determining the evolutionary potential of a species and is essential for developing optimal conservation strategies. The impact of habitat fragmentation on the genetic diversity of food-deceptive orchids seems to be unpredictable because of their specialized seed and pollen dispersal mechanisms. The habitat of deceptive Cypripedium flavum was severely fragmented during the past half century. This study investigated the genetic diversity and structure of seven fragmented Cypripedium flavum populations in Shangrila County using AFLP markers. A total of 376 alleles were identified, with a range of 70 to 81 alleles per locus. The species exhibited considerable genetic diversity, as evidenced by an average Nei’s gene diversity (H) of 0.339 and a Shannon’s information index (I) of 0.505, with all loci being polymorphic. Based on Molecular Variance (AMOVA), 8.75% of the genetic differentiation was found among populations, while the remaining 91.25% of genetic variation occurred within populations. Population structure analysis revealed that the C. flavum germplasm can be categorized into five distinct groups, among which there was significant gene flow. Despite habitat fragmentation, C. flavum still retained a high level of genetic diversity, and the substantial gene flow (5.0826) is a key factor in maintaining the genetic diversity. These findings offer valuable insights for the conservation and potential use of C. flavum genetic resources.

Population connectivity and size reductions in the Anthropocene: the consequence of landscapes and historical bottlenecks in white forsythia fragmented habitats

BackgroundWhite forsythia (Abeliophyllum distichum) is an endangered Korean Peninsula endemic that has been subjected to recent population genomics studies using SNPs via RAD sequencing. Here, we primarily employed the often underutilized haplotype information from RAD loci to further describe the species’ previously uninvestigated haplotype-based genomic variation and structure, and genetic-geographic characteristics and gene flow patterns among its five earlier identified genetic groups. We also inferred the time of past events that may have impacted the effective population size of these groups, as well as the species’ potential future distribution amidst the warming climate and anthropogenic threats.ResultsOur findings emphasized the recognition of the species’ regional patterns of genetic structure, and the role of topography and its associated gene flow patterns as some of the possible factors that may have influenced the species’ present-day fragmented population distribution. The inferred bottleneck events during the Anthropocene, some of which aligned with the time of historical catastrophic events on the Peninsula (e.g., the Korean War), were revealed to have contributed to the generally low effective population size of its five lineages, particularly those with marginal distributional range. Future distribution under both optimistic and pessimistic climatic scenarios suggests unlikely suitable habitats for these populations to expand from their current range limits, at least in the next 80 years.ConclusionsThe small effective population size and landscape-driven limited gene flow among white forsythia populations will remain a big challenge for the conservation management of the species’ already fragmented population distribution. To help mitigate these impacts, the merging of various research approaches and the use of genomic data to their full potential is recommended to provide the optimized knowledge-based tools for the conservation of this endangered species, and other similar plants under pressure.

Implications of newly developed SNPs for conservation programmes for the threatened Nubian ibex (Capra nubiana) in Oman

AbstractConservation of elusive species affected by habitat degradation, population fragmentation and poaching is challenging. The remaining wild populations of a desert-adapted ungulate, Nubian ibex (Capra nubiana), within Oman are small and fragmented. The appropriateness of captive insurance populations for reinforcing existing, or establishing new, wild populations remains uncertain for Oman due to ambiguity regarding their genetic provenance. For effective management of this threatened species, it is essential to assess the genetic relationships between the wild and captive populations, and to investigate hybridisation with domestic goats (Capra hircus). We identified 5,775 high-quality SNPs using double digest restriction-site associated DNA (ddRAD), to assess genetic structure, gene flow and divergence between Oman’s wild populations of Nubian ibex and in captivity, which are likely of North African provenance. We detected hybridisation with goats in captivity and recommend that genetic assessments of captive individuals are routinely used to evaluate their suitability for conservation programs. Building on previous mitochondrial evidence, substantial nuclear divergence (FST = 0.540) was found between wild Oman and captive populations, providing further evidence that Nubian ibex may be composed of multiple species and urgently needs a taxonomic review. Additionally, an appropriate insurance population should be established for Oman’s threatened wild population. The data provided here will be invaluable for developing marker systems to assess wild populations using low-quality DNA from non-invasive sampling. Consequently, it will support further research into Nubian ibex throughout their range and highlights the need to integrate genetic information for effective conservation management of Nubian ibex.

Genetic characterisation of Annona crassiflora as basis for conservation and breeding

Abstract Annona crassiflora is a fruit-bearing tree species native to the Cerrado that has ecological and economic potential, mainly due to the production of attractive and useful fruits for a number of species, including humans. To provide a basis for its conservation and breeding, the objective was to assess the diversity and genetic structure of natural populations of the species using inter-simple sequence repeat (ISSR) molecular markers. Eight populations were analysed in the state of Minas Gerais, Brazil, with 24 individuals randomly sampled in each population, reaching a total of 192 trees. Ten ISSR primers were used, resulting in high genetic diversity for the combined data (H* = 0.35 and I* = 0.52). However, a wide range of representative values was observed in the individual assessment of populations, with JAN, GM and MC standing out for their low genetic diversity, resulting in H* of 0.19, 0.23 and 0.24 and I* of 0.28, 0.34 and 0.35, respectively. An analysis of molecular variance showed greater variation within populations, indicating gene flow (Nm), but genetic differentiation between populations was moderate. Bayesian analysis, although resulting in four genetic groups, revealed the presence of a majority group for the GM and JAN populations. We propose measures to maintain these populations, such as raising awareness of local extractivism and planting genetically divergent seedlings. Furthermore, we recommend including of all populations in conservation and breeding programmes, aiming to cover the maximum genetic variation for the species.

Genetic Diversity Analysis and Fingerprint Construction for 87 Passionfruit (Passiflora spp.) Germplasm Accessions on the Basis of SSR Fluorescence Markers

A comprehensive genetic diversity analysis of 87 Passiflora germplasm accessions domesticated and cultivated for several years in the karst region of Guizhou, China, was conducted utilizing simple sequence repeat (SSR) fluorescent markers. These Passiflora species, renowned for their culinary and medicinal value, could bring significant economic and ecological benefits to the region. This study aimed to assess the genetic resources of these species and facilitate the selection of superior cultivars adapted to the karst environment. Our analysis revealed an abundance of SSR loci within the Passiflora transcriptome, with single-base repeats being the most prevalent type. Through rigorous primer screening and amplification, we successfully identified 27 SSR primer pairs exhibiting robust polymorphisms. Further interrogation at eight microsatellite loci revealed 68 alleles, underscoring the high level of genetic diversity present in the cultivated accessions. The average expected heterozygosity was 0.202, with the ssr18 locus exhibiting the highest value of 0.768, indicating significant genetic variation. The mean polymorphic information content (PIC) of 0.657 indicates the informativeness of these SSR markers. Comparative analyses of the cultivated and potential wild progenitors revealed distinct genetic variations among the different Passiflora types. Genetic structure and clustering analyses of the 87 accessions revealed seven distinct groups, suggesting gene flow and similarities among the resources. Notably, a DNA fingerprinting system was established using eight SSR primer pairs, effectively distinguishing the selected cultivars that had adapted to the karst mountainous region. This study not only deepens our understanding of Passiflora genetic resources in the karst environment but also provides a valuable reference for conservation, genetic improvement, and cultivar selection. The rich genetic diversity of the Passiflora germplasm underscores their potential for sustainable utilization in breeding programs aimed at enhancing the economic and ecological viability of these valuable plant species.

Mapping Ecological Security Patterns Based on Ecosystem Service Valuation in the Qinling-Daba Mountain Area, China: A Multi-Scenario Study for Development and Conservation Tradeoffs

When focusing on biodiversity maintenance, ecological security pattern (ESP) planning gradually becomes a multi-objective planning strategy for sustainable development; wildlife conservation and ecosystem health maintenance should be balanced with local economic development and people’s livelihood enhancement goals. This study focuses on ESP mapping in the Qinling-Daba Mountain area, which is an ecologically significant and socioeconomically underdeveloped area. The tradeoff between conservation and development is made by varying the area of ecological sources and incorporating ecosystem service tradeoffs into ecological source identification through multi-scenario designation. ESPs under six scenarios were generated based on the minimum resistance model, and the important ecological corridors and strategic points in each scenario were identified and compared. The results show the following: (1) The scenario that sets around 30 percent of the study area as ecological sources maintains the integrity of natural ecosystems and leaves space for food and material supply to residents. (2) In this scenario, the ecological sources are connected by 60 corridors that cross 137 townships with high population densities (>100 people/km2) and intersect with major traffic lines at 71 points. Engineering, management, or education strategies must be taken in these townships or intersections to avoid human–wildlife conflicts. (3) The study area needs to construct both short (north–south) corridors linking proximate ecological sources for species’ daily movement and long (west–east) corridors connecting large and distant sources for species’ seasonal migration and gene flow. (4) The multi-scenario approach turns out to be an effective strategy for ESP planning with considerations for development–conservation tradeoffs.

Extensive Population Structure Highlights an Apparent Paradox of Stasis in the Impala (Aepyceros melampus).

Impalas are unusual among bovids because they have remained morphologically similar over millions of years-a phenomenon referred to as evolutionary stasis. Here, we sequenced 119 whole genomes from the two extant subspecies of impala, the common (Aepyceros melampus melampus) and black-faced (A. m. petersi) impala. We investigated the evolutionary forces working within the species to explore how they might be associated with its evolutionary stasis as a taxon. Despite being one of the most abundant bovid species, we found low genetic diversity overall, and a phylogeographic signal of spatial expansion from southern to eastern Africa. Contrary to expectations under a scenario of evolutionary stasis, we found pronounced genetic structure between and within the two subspecies with indications of ancient, but not recent, gene flow. Black-faced impala and eastern African common impala populations had more runs of homozygosity than common impala in southern Africa, and, using a proxy for genetic load, we found that natural selection is working less efficiently in these populations compared to the southern African populations. Together with the fossil record, our results are consistent with a fixed-optimum model of evolutionary stasis, in which impalas in the southern African core of the range are able to stay near their evolutionary fitness optimum as a generalist ecotone species, whereas eastern African impalas may struggle to do so due to the effects of genetic drift and reduced adaptation to the local habitat, leading to recurrent local extinction in eastern Africa and re-colonisation from the South.

Small-scale geographic differences in multiple-driver environmental variability can modulate contrasting phenotypic plasticity despite high levels of gene flow

Climate change is altering not only the mean conditions of marine environments, but also their temporal variability and predictability. As these alterations are not uniform across seascapes, their biological effects are expected to accentuate intra-specific differences in the adaptive capacity (e.g., plasticity and evolutionary potential) of natural populations. To test this theoretical framework, we assessed the phenotypic and genetic profiles of mussel from three study sites across a multi-driver heterogeneous environmental mosaic in Chilean Patagonia. Our study reveals that temporal variability, predictability, and exposure to extreme events (low pH/low salinity), collectively, can modulate the plasticity and optimal conditions of mussels. Despite these phenotypic differences, we observed low genetic differentiation, likely resulting from significant gene flow induced by aquaculture, ultimately diminishing variation among individuals from different geographic areas. Our findings underscore how variability and predictability are essential factors shaping phenotypic diversity, even at small spatial scales. Balancing these factors could enhance species resilience and ecological success, crucial for biodiversity conservation amidst climate change.

Phylogeography of the Red‐Headed Manakin Supports the River‐Refuge Hypothesis

ABSTRACTAimThe Amazon rainforest is one of the most biodiverse regions on earth, but our understanding of the processes that have shaped its patterns of diversity remains incomplete. One hypothesis for Amazonian speciation is the river‐refuge hypothesis, which suggests that retraction of forests away from the periphery of Amazonia, where rivers are narrow, facilitated divergence of populations separated by wide rivers. Later re‐expansion of forests would have allowed secondary contact between these populations, and co‐occurring hybrid zones may reveal the location where expanding forests first reconnected. Here, we test whether a widespread Amazonian songbird species shows evidence of population contact zones in the eastern headwaters of the Tapajós river, hypothesised to be an area of secondary contact under the river‐refuge hypothesis.LocationAmazon and Atlantic forests of South America.TaxonCeratopipra rubrocapilla (Pipridae, Passeriformes).MethodsWe sampled 147 C. rubrocapilla (Red‐headed Manakin) across its vast range, with 70 samples sequenced using reduced‐representation genomic sequencing and 139 sequenced at the mitochondrial gene cytochrome b. We use population genetic and phylogenetic analyses to investigate patterns of gene flow and population structure across the range, with a particular focus on samples from the headwaters of the Tapajós River to evaluate the potential role of the river‐refuge hypothesis.ResultsWe observe shallow population structure and evidence for widespread gene flow across the range of C. rubrocapilla. Patterns in the Tapajós headwaters match predictions of the river‐refuge hypothesis, with eastern headwater samples most similar to samples farther west on the opposite bank of the Tapajós River.Main ConclusionsThe close co‐occurrence between population contact zones within C. rubrocapilla and other previously studied hybrid zones reinforces the hypothesis that the eastern Tapajós headwaters is where rainforest populations reconnected in the past, following predictions of the river‐refuge hypothesis of Amazonian speciation.

Intergroup variability of the khakas based on non-metric cranial traits data

Introduction. The formation of the gene pool of modern Khakas people occurred through the interaction of various groups of ancient inhabitants of the Minusinsk Basin. This article is focused on the inter-group analysis of craniological series representing different Khakassian subethnic groups. Materials and Methods. The cranial series of the Kachins, Koibals, Sagais, and Beltirs was studied using the battery of 36 cranial non-metric characteristics. The obtained data was compared with data on other populations of Southern and Western Siberia: Shors, Tuvans, Telengits, Selkups, Khants, and Mansis. The analysis of the biological diversity of populations was carried out using Smith's mean measure of divergence (MMD) followed by multidimensional scaling and cluster analysis. Results. The Kachins are the most distant from the other Khakas subethnic groups, forming a cluster with the Telengits and Tuvans. The Koibals and Beltirs were positioned quite close to the Shors. The Sagais occupied a separate position in the Khakas cluster, presumably due to a greater proportion of Europoid admixture compared to other Khakas groups. The pooled Khakas sample shares similar cranial non-metric characteristics with Turkic-speaking ethnic groups of Southern Siberia: Shors, Telengits, and Tuvans. Discussion. The nature of phenetic differentiation of the Khakas sub-ethnic groups presumably reveals their complex population history. The position of the Kachins outside the Khakas cluster based on non-metric traits resulted from specificity of their phenofund, which is also supported by molecular genetic data. The similarity in the phenofund of the Khakas groups to those of the Tuvans and Telengits may result from gene flow between Khakas populations and the peoples of Southern Siberia or/and from their common episodes in their ethnogenesis. The similarity between pooled Khakassian sample and the Shors presumably suggests involvement of common genetic components in the gene pools of these peoples, which is consistent with molecular genetic data. Conclusion. The correspondence of the obtained results with genetic data suggests the possibility of using cranial non-metric traits to identify genetic relationships between ancient populations in the absence of direct genetic information.

Differentiation of the Anammox core microbiome: Unraveling the evolutionary impetus of scalable gene flow

Anaerobic ammonium oxidation bacteria (AAOB), distinguished by their unique autotrophic nitrogen metabolism, hold pivotal positions in the global nitrogen cycle and environmental biotechnologies. However, the ecophysiology and evolution of AAOB remain poorly understood, attributed to the absence of monocultures. Hence, a comprehensive elucidation of the AAOB-dominated core microbiome, anammox core, is imperative to further completing the theory of engineered nitrogen removal and ecological roles of anammox. Performing taxonomic and phylogenetic analyses on collected genome repertoires, we show here that Candidatus Brocadia and Candidatus Kuenenia possesses a more compact core than Candidatus Jettenia, which partly explains why the latter has a less common ecological presence. Evidence of gene flow is particularly striking in functions related to biosynthesis and oxygen detoxification, underscoring the evolutionary forces driving lineage and core differentiation. Furthermore, CRISPR spacer traceback of the AAOB metagenome-assembled genomes (MAGs) reveals a series of genetic traces for the concealed phages. By reconceptualizing the functional divergence of AAOB with the historical role of phages, we ultimately propose a coevolutionary framework to understand the evolutionary trajectory of anammox microecology. The discoveries provided in this study offer new insights into understanding the evolution of AAOB and the ecology of anammox.

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.

Breeding season temporal and spatial trends in continental-scale migration of the monarch butterfly.

The monarch butterfly (Danaus plexippus) is a vagile species that undertakes an annual, multi-generational migration across North America. The abundance of this species at both eastern and western overwintering sites in Central Mexico and California indicates a population decline. Success of continental-scale conservation programs for a migratory species depends on providing, maintaining, and protecting habitats at appropriate temporal and spatial scales. Here, dynamics of monarch continental-scale migration and gene flow were obtained by combined stable isotope, morphological, and genetic analyses. These analyses were applied to temporal monarch samples collected from May to September during 2016-2021 at locations in Iowa, USA and spatial collections from Pennsylvania, Delaware, Iowa, Ohio, Nevada, Idaho, Hawaii, 3 Australian locations during July and August 2016, and Texas in April 2021. Evidence of seasonal multi-generational migration was obtained through δ2H analyses of spatial collections, which was corroborated by decreased wing hue (a morphological marker for non-migratory individuals). In Iowa, 10-15% of monarchs represented migrants from southern areas throughout the breeding season and 6% were migrants from the North in midsummer. Limited sequence variation detected across the mitochondrial genome impacted the capability to detect significant population genetic variation in our North American samples. However, 2 novel substitutions were identified and predicted to be fixed among Australia samples, contributing to intercontinental differentiation from counterparts in North America. Our assessment of temporal and spatial population dynamics across the North American monarch breeding range provides insight into continental-scale migration and previously undetected mitochondrial DNA variation among globally distributed populations.

Analysis of genetic diversity and environmental associations of wild citron (Citrus medica L.) in northeast India

Citron (Citrus medica L.) is a major citrus fruit crop known for its medicinal, nutritional, and aromatic uses. Native to northeast India, Myanmar, and south China, citron is cultivated mostly in the Mediterranean and South Asian countries. The truly wild populations of citron are found in the hills and valleys of northeast India. The extent of genetic diversity in wild citron in relation to environmental variables has not been documented. In the present study, we analyzed genetic diversity and environmental association of wild populations of C. medica in northeast India using 15 polymorphic nuclear SSR (simple sequence repeat) markers and 21 environmental variables. SSR profiling revealed a moderate genetic diversity (Ho = 0.40; He = 0.41; I = 0.70). Neighbor-joining (NJ) and principal coordinate analysis (PCoA) grouped 16 natural populations into two genetic clusters. Bayesian structure results showed two genetic populations with a high admixture of individuals. Analysis of molecular variance (AMOVA) revealed low genetic differentiation (Fst = 0.14; P < 0.001), indicative of high gene flow (Nm = 1.53) in citron. The isolation by distance (IBD) analysis showed a weak correlation, indicating that geographic distance was not a determining factor for genetic distance in C. medica populations. The higher metrics of outcrossing rate (tm > 0.9) revealed outbreeding nature of the species. The genetic diversity of citron showed a positive correlation with precipitation of driest month and precipitation of coldest quarter, but negative with mean diurnal range and isothermality. Overall, the results of the present study could benefit the conservation and management of citron resources and their natural habitat in northeast India. Genetically highly diverse populations identified in the present study can be used for conservation and utilization of citron genetic resources.

Myoglobin inhibits breast cancer cell fatty acid oxidation and migration via heme-dependent oxidant production and not fatty acid binding

The monomeric heme protein myoglobin (Mb) is aberrantly expressed in approximately 40 % of breast tumors. Mb expression is associated with better patient prognosis, yet the molecular mechanisms underlying this effect are unclear. In muscle, Mb's heme moiety confers oxygen storage and delivery. However, prior studies demonstrate that low levels of Mb in cancer cells preclude this function. Several studies propose a fatty acid binding function for Mb via lysine residue K46. Because cancer cells can upregulate fatty acid oxidation (FAO) to fuel cell migration, we tested whether Mb-mediated fatty acid binding modulates FAO and migration. We demonstrate that stable expression of human Mb in MDA-MB-231 breast cancer cells decreases cell migration and FAO. Site-directed mutagenesis of Mb K46 disrupted fatty acid binding but did not improve FAO or migration. Conversely, cells expressing Apo-Mb (with disrupted heme binding) did not show impaired FAO or migration rates, suggesting Mb attenuates FAO and migration via a heme-dependent mechanism rather than through fatty acid binding. Mb's heme-dependent oxidant generation dysregulates migratory gene expression, which is reversed by catalase treatment. Collectively, these data demonstrate that Mb's heme-dependent oxidant production decreases breast cancer cell migration, prompting therapeutic strategies to modulate oxidant production and Mb in tumors.

Agro-morphological genetic diversity assessment of Amaranthus genotypes from Ethiopia based on qualitative traits

Amaranths are a promising plant in the family of Amaranthaceae because of their nutritional and functional properties, such as their high antioxidant content and dietary fiber content. However, it's being disregarded for several reasons, such as ignorance, a lack of in-depth research, and the plant's long-term genetic growth in Ethiopia, among other factors. In the current work, we described the genetic diversity of 120 amaranth genotypes using qualitative criteria. The experiment, which was configured with an alpha lattice design and duplicated twice, was run for two seasons in 2020 and 2021. Twenty qualitative descriptors were looked for in the gathered data. Among the 20 qualitative variables that were assessed, the chi-square test result indicated the presence of prevailing phenotypic variation. The results of the agro-morphological characterization also revealed a significant amount of variance. The overall mean of the Shannon diversity indices (H') was 0.61. The indices for germination rate, leaf margin, prominence of leaf veins, and the existence of auxiliary inflorescence varied from 0.12 to 0.99. The estimated diversity indices showed more intra-regional diversity (0.66) than inter-regional diversity (0.34), demonstrating the existence of gene flow between growing regions. Shannon–Weaver Diversity Index, ranged from 0.00 for auxiliary inflorescence to 1.94 for leaf coloration, with an overall mean of 19 characters (95%) that were found to have high diversity (> 0.76) while auxiliary inflorescence was invariant. Except for auxiliary inflorescence, all qualitative features showed a lot of variation. Additionally, Amaranthus hybridus L. subsp. cruentus (L.) Theil recorded the greatest Shannon diversity index (0.47) while Amaranthus spinosus L. recorded the lowest (0.00). The hierarchical clustering grouped all the genotypes into three clusters. The first cluster included the most genotypes (58), followed by the second (47), and the third cluster contained the fewest (15). Principal component analysis showed that the first six principal components with eigenvalues greater than one contributed 72% of the variability among genotypes. The study unequivocally demonstrated that, even when the genotypes were grouped into a small number of clusters, there was still enough divergence within the clusters to demonstrate the genotypes of amaranth to have a high genetic diversity. These results indicate that there is substantial genetic diversity among Ethiopian amaranth genotypes, which should be safeguarded and may be utilized in breeding in the future.

Composite dyadic models for spatio-temporal data.

Mechanistic statistical models are commonly used to study the flow of biological processes. For example, in landscape genetics, the aim is to infer spatial mechanisms that govern gene flow in populations. Existing statistical approaches in landscape genetics do not account for temporal dependence in the data and may be computationally prohibitive. We infer mechanisms with a Bayesian hierarchical dyadic model that scales well with large data sets and that accounts for spatial and temporal dependence. We construct a fully connected network comprising spatio-temporal data for the dyadic model and use normalized composite likelihoods to account for the dependence structure in space and time. We develop a dyadic model to account for physical mechanisms commonly found in physical-statistical models and apply our methods to ancient human DNA data to infer the mechanisms that affected human movement in Bronze Age Europe.

Field-scale gene flow of fungicide resistance in Pyrenophora teres f. teres and the effect of selection pressure on the population structure.

The effectiveness of fungicides to control foliar fungal crop diseases is being diminished by the increasing spread of resistances to fungicides. One approach that may help to maintain efficacy is remediation of resistant populations by sensitive ones. However, the success of such approaches can be compromised by re-incursion of resistance through aerial spore dispersal; although, knowledge of localized gene flow is lacking. Here, we report on a replicated mark-release-recapture field experiment with several treatments set up to study spore-dispersal-mediated gene flow of a mutated allele that confers demethylase inhibitor resistance in Pyrenophora teres f. teres (Ptt). Artificial inoculation of the host, barley (Hordeum vulgare), was successful across the 12-ha trial, where the introduced sensitive- and resistant-populations were, respectively, 6- and 13-fold the DNA concentration of the native Ptt population. Subsequent disease pressure remained low which hampered spread of the epidemic to such extent that gene flow was not detected at, or beyond 2.5 m from source points. In the absence of gene flow, plots were assessed for treatment effects; fungicide applied to populations that contained 14.3% of allele mutation increased in frequency to 24.5%, whereas sensitive populations had no change in structure. Untreated controls of native Ptt population remained genetically stable, yet untreated controls that were inoculated with sensitive Ptt had half the resistance frequency of the native population structure. The trial demonstrates the potential for management to remediate fungicide resistant pathogen populations, where localized gene flow is minimal; to safeguard chemical crop protection into the future.

Genomic adaptation to small population size and saltwater consumption in the critically endangered Cat Ba langur

Many mammal species have declining populations, but the consequences of small population size on the genomic makeup of species remain largely unknown. We investigated the evolutionary history, genetic load and adaptive potential of the Cat Ba langur (Trachypithecus poliocephalus), a primate species endemic to Vietnam’s famous Ha Long Bay and with less than 100 living individuals one of the most threatened primates in the world. Using high-coverage whole genome data of four wild individuals, we revealed the Cat Ba langur as sister species to its conspecifics of the northern limestone langur clade and found no evidence for extensive secondary gene flow after their initial separation. Compared to other primates and mammals, the Cat Ba langur showed low levels of genetic diversity, long runs of homozygosity, high levels of inbreeding and an excess of deleterious mutations in homozygous state. On the other hand, genetic diversity has been maintained in protein-coding genes and on the gene-rich human chromosome 19 ortholog, suggesting that the Cat Ba langur retained most of its adaptive potential. The Cat Ba langur also exhibits several unique non-synonymous variants that are related to calcium and sodium metabolism, which may have improved adaptation to high calcium intake and saltwater consumption.