Strong Bottleneck Research Articles

Identification of potential auxin response candidate genes for soybean rapid canopy coverage through comparative evolution and expression analysis.

Throughout domestication, crop plants have gone through strong genetic bottlenecks, dramatically reducing the genetic diversity in today's available germplasm. This has also reduced the diversity in traits necessary for breeders to develop improved varieties. Many strategies have been developed to improve both genetic and trait diversity in crops, from backcrossing with wild relatives, to chemical/radiation mutagenesis, to genetic engineering. However, even with recent advances in genetic engineering we still face the rate limiting step of identifying which genes and mutations we should target to generate diversity in specific traits. Here, we apply a comparative evolutionary approach, pairing phylogenetic and expression analyses to identify potential candidate genes for diversifying soybean (Glycine max) canopy cover development via the nuclear auxin signaling gene families, while minimizing pleiotropic effects in other tissues. In soybean, rapid canopy cover development is correlated with yield and also suppresses weeds in organic cultivation. We identified genes most specifically expressed during early canopy development from the TIR1/AFB auxin receptor, Aux/IAA auxin co-receptor, and ARF auxin response factor gene families in soybean, using principal component analysis. We defined Arabidopsis thaliana and model legume species orthologs for each soybean gene in these families allowing us to speculate potential soybean phenotypes based on well-characterized mutants in these model species. In future work, we aim to connect genetic and functional diversity in these candidate genes with phenotypic diversity in planta allowing for improvements in soybean rapid canopy cover, yield, and weed suppression. Further development of this and similar algorithms for defining and quantifying tissue- and phenotype-specificity in gene expression may allow expansion of diversity in valuable phenotypes in important crops.

Extreme founder effect associated with hyperglycemia and hyperlipidemia on the island of NIAS/Indonesia

The island of Nias/Indonesia shows an extremely reduced genetic diversity indicating a strong founder effect. As a consequence, the prevalence of some disease genes should significantly differ among populations depending on the gene pool passed on to the founder population and their successive expansion as it has already been documented for several monogenic diseases. Results of the current study based on routine laboratory blood examination give rise to the notion that this might also hold true for polygenic disorders. We observed very high prevalence of hyperglycemia (non-fasting glucose above 200 mg/dL in 14 % Nias population compared to 1.5 % in the population of the neighboring island of Sumatra) accompanied by hypertriglyceridemia, high non-HDL-cholesterol, and low HDL-cholesterol levels. These findings suggest that the Nias population may be disproportionally affected by prediabetes and type 2 diabetes mellitus. By contrast, laboratory parameters potentially indicative of other polygenic disorders such as total plasma cholesterol, electrolytes, creatinine, urea, and uric acid were comparable between the inhabitants of Nias and Sumatra islands. To our knowledge this is the first study suggesting that the extremely strong genetic bottleneck seen in the Nias population translates into the widespread metabolic disease with potentially deleterious influence on public health.

Population genetic structure of the European pond turtle (Emys orbicularis) at the northern edge of the Pannonian Basin

Populations of plants and animals at the edge of a species’ range may exhibit lower genetic diversity due to a combination of historical (e.g., postglacial dispersal) and contemporary (e.g., low effective population size, strong genetic drift or bottlenecks) processes. In this study, we analysed the genetic structure of populations of the European pond turtle (Emys orbicularis) in the northern part of the Pannonian Basin (Slovakia, Hungary, Austria), which represents the edge of the species’ range. Using mitochondrial DNA (mtDNA) and microsatellite markers, we aimed to assign individuals and populations into phylogeographic lineages and to determine genetic variation and population structure. Our study revealed that most individuals from natural populations belong to mtDNA lineage II (widely distributed haplotype IIa, and likely endemic haplotypes IIl and IIo), which is typical for the whole Pannonian region. However, non-native haplotypes were detected in one Slovak population (haplotype Ib) and in the Austrian Donau Auen National Park (haplotype IVa). Microsatellite markers revealed high variability comparable to that in other parts of the range, suggesting that marginal populations may not necessarily have reduced genetic diversity. A genetically mixed population was found in the Donau Auen National Park, consisting of turtles belonging to mtDNA lineage II (taxonomically assessed as subspecies E. orbicularis orbicularis) and lineage IV (taxonomically assessed as E. orbicularis hellenica), which is native to the circum-Adriatic region. The admixture suggests intense hybridization and introgression between the two lineages (subspecies).

Phylogeography and genetic structure of Papaver bracteatum populations in Iran based on genotyping-by-sequencing (GBS)

Papaver bracteatum, known for its high thebaine content and absence of morphine, has emerged as a promising alternative to opium poppy for codeine production. In this study, our objective was to create a diverse panel representing the natural variation of this species in Iran. To achieve this, we employed genotyping-by-sequencing to obtain genome-wide distributed single-nucleotide polymorphisms (SNPs) for phylogeographic analysis, population structure assessment, and evaluation of genetic diversity within P. bracteatum populations. A total of 244 P. bracteatum individuals from 13 distinct populations formed seven genetic groups, along with one highly admixed population. We observed a clear split between the populations inhabiting the Alborz Mts. in the east and Zagros Mts. in the west. In between these mountain ranges, the population of Kachal Mangan exhibited a high degree of genetic admixture between both genetic groups. At or after the end of the last glacial maximum, when climate conditions rapidly changed, all P. bracteatum populations experienced a strong demographic bottleneck reducing the already small effective population sizes further before they increased to their recent strengths. Our results suggest that the ongoing climate change together with human pressure on the species’ habitats and limited seed-dispersal ability are potential factors contributing today to rising genetic isolation of P. bracteatum populations. Our results provide genetic data that can be used for conservation measures to safeguard the species’ genetic diversity as a resource for future breeding approaches in this medicinally important species.

Tracing the fine-scale demographic history and recent admixture in Hmong-Mien speakers.

The linguistic, historical, and subsistent uniqueness of Hmong-Mien (HM) speakers offers a wonderful opportunity to investigate how these factors impact the genetic structure. The genetic differentiation among HM speakers and their population history are not well characterized. Here, we generate genome-wide data from 65 Yao ethnicity samples and analyze them with published data, particularly by leveraging haplotype-based methods. We determined that the fine-scale genetic substructure of HM speakers corresponds better with linguistic classification than with geography. Particularly, parallels between serial founder events and language differentiations can be observed in West Hmongic speakers. Multiple lines of evidence indicate that ~500-year-old GaoHuaHua individuals are most closely related to West Hmongic-speaking Bunu. The strong genetic bottleneck of some HM-speaking groups, especially Bunu, could potentially be associated with their long-term practice of swidden agriculture to some degree. The inferred admixture dates for most of the HM speakers overlap with the reign of the Ming dynasty (1368-1644 CE). Besides a common genetic origin for HM speakers, their genetic ancestry is shared primarily with neighboring Han Chinese and Tai-Kadai speakers in south China. In conclusion, our analyses reveal that recent isolation and admixture events have contributed to the genetic population history of present-day HM speakers.

Pollinator-assisted plant phenotyping, selection, and breeding for crop resilience to abiotic stresses.

Food security is threatened by climate change, with heat and drought being the main stresses affecting crop physiology and ecosystem services, such as plant-pollinator interactions. We hypothesize that tracking and ranking pollinators' preferences for flowers under environmental pressure could be used as a marker of plant quality for agricultural breeding to increase crop stress tolerance. Despite increasing relevance of flowers as the most stress sensitive organs, phenotyping platforms aim at identifying traits of resilience by assessing the plant physiological status through remote sensing-assisted vegetative indexes, but find strong bottlenecks in quantifying flower traits and in accurate genotype-to-phenotype prediction. However, as the transport of photoassimilates from leaves (sources) to flowers (sinks) is reduced in low-resilient plants, flowers are better indicators than leaves of plant well-being. Indeed, the chemical composition and amount of pollen and nectar that flowers produce, which ultimately serve as food resources for pollinators, change in response to environmental cues. Therefore, pollinators' preferences could be used as a measure of functional source-to-sink relationships for breeding decisions. To achieve this challenging goal, we propose to develop a pollinator-assisted phenotyping and selection platform for automated quantification of Genotype × Environment × Pollinator interactions through an insect geo-positioning system. Pollinator-assisted selection can be validated by metabolic, transcriptomic, and ionomic traits, and mapping of candidate genes, linking floral and leaf traits, pollinator preferences, plant resilience, and crop productivity. This radical new approach can change the current paradigm of plant phenotyping and find new paths for crop redomestication and breeding assisted by ecological decisions.

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.

Msaken, Tunisia: A Common Paternal Ancestor Confirmed by Y Chromosome DNA Analysis

Msaken City (Tunisia) is believed to have been founded around 1360 AD by five related men who migrated from West Asia. The population would have grown with the descendants of these founders and with the arrival of other populations from different regions of Tunisia. In order to elucidate the TMRCA of the founder population and to reveal their geographic origin, 23 males from different families of Msaken were examined, using the services of commercial companies, for 12 to 440 Y chromosome Short Tandem Repeats (STR) and Single Nucleotide Polymorphisms (SNP) markers using NG testing technology. Eight samples were genotyped for SNPs to determine their Haplogroups. In order to refine the phylogeny, traditional Sanger testing was performed on one sample for 300,000 bp in Y Chr (in Walk Through the Y chromosome project). Seven samples were also tested using Next Generation Testing (BigY) covering 20 million bp of Y chromosome overlapping 85% of Gold standard region (chromosome Y positions placed on the phylogenetic tree by the YCC) using NGS instruments, HiSeq 2000 and 2500. A comparison of STR results with data from different sources and databases was made, using SQL scripts and data mining tools, to find matching haplotypes. All the STR results were found to have no more than three mismatches per 12 markers and not more than six mismatches per 67 markers and the SNP results showed that all tested samples belonged to Haplogroup J-M172 inside its subgroup J-L24. Relying on the common STR marker values, we define a Msaken-Haplotype. NG tests for our samples as well as those added to the yfull.com tree allowed us to refine the phylogeny of J-L24 and the samples were all found to belong to J-L271 Haplogroup and share 54 exclusive SNPs. The calculated Time to Most Recent Common Ancestor (TMRCA) based on NG testing, ranges between1500 and 6200 YBP showing a strong bottleneck around 5400 YBP. The variation of the collected results shows a geographic root of J-L192 in East Anatolia, present day Armenia, Azerbaijan, and West Iran. 20 to 30% of random Tunisian STR haplotypes belonging to J-M172 (J2) Haplogroup exhibit the Msaken-Haplotype.

Genomic insights into the evolutionary and demographic histories of the extinct Hoopoe Starling (Fregilupus varius)

The Hoopoe Starling Fregilupus varius is an extinct species of the Sturnidae that was endemic to Réunion Island in the Indian Ocean. The species rapidly disappeared in the middle of the 19th century, primarily because of overexploitation by humans. We generated an approximately 11× coverage genome to reconstruct the demographic history of the Hoopoe Starling and compared these results with the demographic histories of other starlings and myna species. Our analyses confirmed the close affinities of the Hoopoe Starling with the genera Sturnia, Leucopsar and Sturnornis, and revealed that it went through a strong population bottleneck during its evolutionary history, but that its effective population size was not particularly low when compared with other extinct or critically endangered species of birds.

Global genetic diversity and historical demography of the Bull Shark

AbstractAimBiogeographic boundaries and genetic structuring have important effects on the inferences and interpretation of effective population size (Ne) temporal variations, a key genetics parameter. We reconstructed the historical demography and divergence history of a vulnerable coastal high‐trophic shark using population genomics and assessed our ability to detect recent bottleneck events.LocationWestern and Central Indo‐Pacific (IPA), Western Tropical Atlantic (WTA) and Eastern Tropical Pacific (EPA).TaxonCarcharhinus leucas (Müller & Henle, 1839).MethodsA DArTcap™ approach was used to sequence 475 samples and assess global genetic structuring. Three demographic models were tested on each population, using an ABC‐RF framework coupled with coalescent simulations, to investigate within‐cluster structure. Divergence times between clusters were computed, testing multiple scenarios, with fastsimcoal. Ne temporal variations were reconstructed with STAIRWAYPLOT. Coalescent simulations were performed to determine the detectability of recent bottleneck under the estimated historical trend for datasets of this size.ResultsThree genetic clusters corresponding to the IPA, WTA and EPA regions were identified, agreeing with previous studies. The IPA presented the highest genetic diversity and was consistently identified as the oldest. No significant within‐cluster structuring was detected. Ne increased globally, with an earlier onset in the IPA, during the last glacial period. Coalescent simulations showed that weak and recent bottlenecks could not be detected with our dataset, while old and/or strong bottlenecks would erase the observed ancestral expansion.Main ConclusionsThis study further confirms the role of marine biogeographic breaks in shaping the genetic history of large mobile marine predators. Ne historical increases in Ne are potentially linked to extended coastal habitat availability. The limited within‐cluster population structuring suggests that Ne can be monitored over ocean basins. Due to insufficient amount of available genetic data, it cannot be concluded whether overfishing is impacting Bull Shark genetic diversity, calling for whole‐genome sequencing.

Identification of the CNGC Gene Family in Rice and Mining of Alleles for Application in Rice Improvement.

Cyclic nucleotide-gated ion channel (CNGC) gene regulation plays important roles in plant immune and abiotic stress response. Here, we identified 16 CNGC genes in rice (Oryza sativa). Then, we analyzed their chromosomal location, physicochemical properties, subcellular localization, gene functional interaction network, cis-acting elements, phylogenetic relationships, collinearity, expression in tissues under normal conditions and abiotic stresses, and geng-cds-haplotype (gcHap) diversity in 3010 gcHaps. As a result, OsCNGC3 (Os06g0527300) was identified as a gene different from previous report, and OsCNGC genes were found to play important roles in rice population differentiation and rice improvement. Our results revealed their very strong differentiation between subspecies and populations, important roles in response to abiotic stresses, as well as strong genetic bottleneck effects and artificial selection of gcHap diversity in the modern breeding process of Xian (indica) and Geng (japonica) populations. The results also suggested that natural variations in most rice CNGC loci are potentially valuable for improving rice productivity and tolerance to abiotic stresses. The favorable alleles at the CNGC loci should be explored to facilitate their application in future rice improvement.

A source of hidden diversity: soil seed bank and aboveground populations of a common herb contain similar levels of genetic variation.

In many landscapes, successful re-establisment of plant populations depends on the presence of diaspores, either near or directly beneath sites to be restored. The soil seed bank is, therefore, an important part of ecosystem resilience and a vital pillar for regeneration of genetic diversity in many plant populations. However, regeneration from the soil seed bank and the siubsequent restoration can only be considered successful when genetic diversity of restored populations is not eroded nor genetic differentiation inflated. We compared genetic variation within and among soil seed bank and aboveground populations of Origanum vulgare, to test whether genetically variable populations can be restored from the soil seed bank. We explored levels of genetic diversity within aboveground populations and the corresponding soil seed banks. Furthermore, we assessed the extent to which the soil seed bank differs genetically from the aboveground population. Levels of genetic diversity were to generally similar in aboveground populations and the corresponding soil seed banks. Only levels of inbreeding were slightly higher in the lower layer of the soil seed bank compared to the aboveground populations, probably because of selection processes acting against homozygotes accumulating in the seed bank. Furthermore, significant genetic differentiation between the aboveground population and the corresponding seed banks was completely absent. Across all sites, genetic differentiation between the soil seed bank was similar to that between aboveground populations, probably due to the absence of severe climate conditions, strong bottlenecks or disturbance events. Our conclusions support the possibility of successful re-establishment of healthy, genetically variable plant populations after aboveground destruction or following soil re-allocation from persistent seed banks.

Genomic analyses indicate the North American Ap-ha variant of the tick-vectored bacterium Anaplasma phagocytophilum was introduced from Europe

BackgroundAnaplasma phagocytophilum is a tick-vectored, obligately intracellular bacterium that infects a diversity of vertebrate hosts. In North America, the Ap-ha variant of A. phagocytophilum can cause dangerous infections in humans, whereas symptomatic human infections in Europe are rare. Conversely, the European host-generalist ecotype of A. phagocytophilum frequently causes illness in domestic ruminants while no comparable infections have been recorded from North America. Despite these differences in pathogenicity, the Ap-ha variant is closely aligned phylogenetically with the European host-generalist ecotype. Furthermore, North American populations of A. phagocytophilum are less genetically diverse than those in Europe. Taken together, these observations suggest that the North American Ap-ha variant may represent an introduced population of this bacterium.Methods Data from publicly available whole genomes of A. phagocytophilum were used to compare phylogeographic patterns and the extent of genetic divergence between the North American Ap-ha variant and the European host-generalist ecotype.ResultsThe results confirm that North American Ap-ha samples are phylogenetically nested within the diversity of the European host-generalist ecotype, and that Ap-ha likely radiated within the last 100 years. As expected, the Ap-ha variant also exhibited relatively low genetic diversity levels compared to the European host-generalist ecotype. Finally, North American Ap-ha harbored significantly more derived alleles than the European host-generalist A. phagocytophilum population.ConclusionsCollectively, these results support the hypothesis that the Ap-ha variant was recently introduced to North America from Europe and underwent a strong genetic bottleneck during this process (i.e. a ‘founder event’). Adaptation to novel vectors may have also played a role in shaping genetic diversity and divergence patterns in these pathogenic bacteria. These findings have implications for future studies aimed at understanding evolutionary patterns and pathogenicity variation within A. phagocytophilum.Graphical

Loss of West Nile virus genetic diversity during mosquito infection due to species-dependent population bottlenecks

Vector competence (VC) refers to the efficiency of pathogen transmission by vectors. Each step in the infection of a mosquito vector constitutes a barrier to transmission that may impose bottlenecks on virus populations. West Nile virus (WNV) is maintained by multiple mosquito species with varying VC. However, the extent to which bottlenecks and VC are linked is poorly understood. Similarly, quantitative analyses of mosquito-imposed bottlenecks on virus populations are limited. We used molecularly barcoded WNV to quantify tissue-associated population bottlenecks in three variably competent WNV vectors. Our results confirm strong population bottlenecks during mosquito infection that are capable of dramatically reshaping virus population structure in a non-selective manner. In addition, we found that mosquitoes with differing VC uniquely shape WNV population structure: highly competent vectors are more likely to contribute to the maintenance of rare viral genotypes. These findings have important implications for arbovirus emergence and evolution.

Functional characterization and allelic mining of OsGLR genes for potential uses in rice improvement.

Glutamate-like receptor (GLR) genes are a group of regulatory genes involved in many physiological processes of plants. With 26 members in the rice genome, the functionalities of most rice GLR genes remain unknown. To facilitate their potential uses in rice improvement, an integrated strategy involving CRISPR-Cas9 mediated knockouts, deep mining and analyses of transcriptomic responses to different abiotic stresses/hormone treatments and gene CDS haplotype (gcHap) diversity in 3,010 rice genomes was taken to understand the functionalities of the 26 rice GLR genes, which led us to two conclusions. First, the expansion of rice GLR genes into a large gene family during evolution had gone through repeated gene duplication events occurred primarily in two large GLR gene clusters on rice chromosomes 9 and 6, which was accompanied with considerable functional differentiation. Secondly, except for two extremely conserved ones (OsGLR6.2 and OsGLR6.3), rich gcHap diversity exists at the remaining GLR genes which played important roles in rice population differentiation and rice improvement, evidenced by their very strong sub-specific and population differentiation, by their differentiated responses to day-length and different abiotic stresses, by the large phenotypic effects of five GLR gene knockout mutants on rice yield traits, by the significant association of major gcHaps at most GLR loci with yield traits, and by the strong genetic bottleneck effects and artificial selection on the gcHap diversity in populations Xian (indica) and Geng (japonica) during modern breeding. Our results suggest the potential values of the natural variation at most rice GLR loci for improving the productivity and tolerances to abiotic stresses. Additional efforts are needed to determine the phenotypic effects of major gcHaps at these GLR loci in order to identify 'favorable' alleles at specific GLR loci specific target traits in specific environments to facilitate their application to rice improvement in future.

Population demographic history and evolutionary rescue: Influence of a bottleneck event.

Rapid environmental change presents a significant challenge to the persistence of natural populations. Rapid adaptation that increases population growth, enabling populations that declined following severe environmental change to grow and avoid extinction, is called evolutionary rescue. Numerous studies have shown that evolutionary rescue can indeed prevent extinction. Here, we extend those results by considering the demographic history of populations. To evaluate how demographic history influences evolutionary rescue, we created 80 populations of red flour beetle, Tribolium castaneum, with three classes of demographic history: diverse populations that did not experience a bottleneck, and populations that experienced either an intermediate or a strong bottleneck. We subjected these populations to a new and challenging environment for six discrete generations and tracked extinction and population size. Populations that did not experience a bottleneck in their demographic history avoided extinction entirely, while more than 20% of populations that experienced an intermediate or strong bottleneck went extinct. Similarly, among the extant populations at the end of the experiment, adaptation increased the growth rate in the novel environment the most for populations that had not experienced a bottleneck in their history. Taken together, these results highlight the importance of considering the demographic history of populations to make useful and effective conservation decisions and management strategies for populations experiencing environmental change that pushes them toward extinction.

Small-scale genetic structure of populations of the bulb mite Rhizoglyphus robini.

Bulb mites are an economically significant pest of subterranean parts of plants and a versatile laboratory animal. However, the genetic structure of their populations remains unknown. To fill this gap in our knowledge of their biology, we set up a field experiment in which we allowed mites to colonize onion bulbs, and then determined the genetic structure of colonisers based on a panel of microsatellite loci. We found moderate but significant population structure among sites separated by ca. 20m (FST range 0.03-0.21), with 7% of genetic variance distributed among sites. Allelic richness within some bulbs was nearly as high as that in the total population, suggesting that colonisation of bulbs was not associated with strong population bottlenecks. The significant genetic structure we observed over small spatial scales seems to reflect limited dispersal of mites in soil.

High-yield preparation of coal tar pitch based porous carbon via low melting point fire retardant carbonation strategy for supercapacitor

Currently, the preparation of porous carbon materials with low-cost, high yield and excellent capacitive performance that carbonized carbon precursors in air is suffering from strong technical bottlenecks. In this work, we reported an efficient and simple method for the preparation of high yield porous carbon by carbonizing coal tar pitch (CTP) with the assistance of fire retardants in muffle furnace in air. We investigated in detail the effect of potassium salt fire retardants with different melting points on the yield and morphology of porous carbon, and the fire retarding mechanism in air. Fire retardants could not only prevent CTP from burning completely, but also used as heteroatomic source and morphology regulator to get the heteroatom doped honeycomb structure carbon. It was found that the lower the melting point of the fire retardants, the higher the carbon yield. The yield of N and F doped porous carbon (CKKF700-3) prepared by the optimized K2CO3-KCl-KF molten salt fire retardant was 40.2 wt%, and exhibited excellent supercapacitor performance (379.8F g−1). The energy density of the symmetrical supercapacitor device composed of CKKF700-3 and K2CO3 based deep eutectic solvent electrolyte was 43.7 Wh kg−1, and it could work properly at extreme temperature from −25 to 105 °C

Genomic insights into the postintroduction failure of the Asian icefish Protosalanx chinensis in China.

Biological introductions provide a natural ecological experiment unfolding in a recent historical timeframe to elucidate how evolutionary processes (such as founder effects, genetic diversity and adaptation) shape the genomic landscape of populations postintroduction. The Asian icefish, Protosalanx chinensis, is an economically important fishery resource, deliberately introduced into dozens of provinces across China for decades. However, while invading and disturbing the local ecosystem, many introduced populations declined, disappearing mysteriously in a very short time. The way in which various evolutionary forces integrate to result in invasion failure of an introduced population remains unknown. Here, we performed whole-genome sequencing of 10 species from the Salangidae family and 70 Asian icefish (Protosalanx chinensis) individuals from 7 geographic populations in China, aiming to characterize the evolutionary fate of introduced populations. Our results show that compared to other Salangidae species, P. chinensis has low genetic diversity, potentially due to the long-lasting decline in population size. In a recently introducted population, Lugu lake, severe sampling effects and a strong bottleneck further deteriorated the genomic landscape. Although the introduced population showed signs of reduced genetic load, the purging selection efficiency was low. Our selective sweep analysis revealed site frequency changes in candidate genes, including gata1a and hoxd4b, which could be associated with a decrease in dissolved oxygen in the deep-water plateau lake. These findings caution against the widespread introduction of P. chinensis in China and lay the groundwork for future use of this economically species.

ANKH: A Generalized O(N) Interpolated Ewald Strategy for Molecular Dynamics Simulations

To evaluate electrostatics interactions, molecular dynamics (MD) simulations rely on Particle Mesh Ewald (PME), an algorithm that uses Fast Fourier Transforms (FFTs) or, alternatively, on Fast Multipole Methods (FMM) approaches. However, the FFTs low scalability remains a strong bottleneck for large-scale PME simulations on supercomputers. On the opposite, FFT-free FMM techniques are able to deal efficiently with such systems but they fail to reach PME performances for small- to medium-size systems, limiting their real-life applicability. We propose ANKH, a strategy grounded on interpolated Ewald summations and designed to remain efficient/scalable for any size of systems. The method is generalized for distributed point multipoles, and so for induced dipoles, which makes it suitable for high performance simulations using new generation polarizable force fields toward exascale computing.