Estimating scale-specific and localized spatial patterns in allele frequency.

Characterizing spatial patterns in allele frequencies is fundamental to evolutionary biology because these patterns contain evidence of underlying processes. However, the spatial scales at which gene flow, changing selection, and drift act are often unknown. Many of these processes can operate inconsistently across space, causing nonstationary patterns. We present a wavelet approach to characterize spatial pattern in allele frequency that helps solve these problems. We show how our approach can characterize spatial patterns in relatedness at multiple spatial scales, i.e. a multilocus wavelet genetic dissimilarity. We also develop wavelet tests of spatial differentiation in allele frequency and quantitative trait loci (QTL). With simulation, we illustrate these methods under different scenarios. We also apply our approach to natural populations of Arabidopsis thaliana to characterize population structure and identify locally adapted loci across scales. We find, for example, that Arabidopsis flowering time QTL show significantly elevated genetic differentiation at 300-1,300 km scales. Wavelet transforms of allele frequencies offer a flexible way to reveal geographic patterns and underlying evolutionary processes.

Whole genome resequencing reveals candidate genes for postaxial polydactyly in Large White pigs.

Polydactyly is a genetic abnormality that affects both pig welfare and industry profits. Despite efforts to explore the genetic basis of pig polydactyly, progress remains limited. In this study, we analyzed a group of Large White pigs with postaxial polydactyly, including 29 cases and 79 controls from 24 families. High-depth sequencing was performed on 20 pigs, while low-depth sequencing was improved through imputation for the remaining pigs. A genome-wide association study (GWAS) and genetic differentiation were conducted using the resequencing dataset, resulting in the identification of 48 significantly associated SNPs and 27 candidate regions. The genetic differentiation regions on chromosomes 5 and 18, which harbored GWAS-identified SNPs, were delineated as confidence regions. The confidence region at Chr18: 1.850-1.925 Mb covers the fifth intron of LMBR1, a gene that contains an important regulatory element for SHH, known as ZRS. Mutations in this ZRS have been found to cause polydactyly in animals and humans. Therefore, we propose LMBR1 as a prospective candidate gene for postaxial polydactyly. These findings emphasize the importance of exploring the role of ZRS within LMBR1 in the pathogenesis of polydactyly in pigs.

Allele frequencies and selection coefficients in locally adapted populations.

Understanding the role of natural selection in driving evolutionary change requires accurate estimates of the strength of selection acting at the genetic level in the wild. This is challenging to achieve but may be easier in the case of populations in migration-selection balance. When two populations are at equilibrium under migration-selection balance, there exist loci whose alleles are selected different ways in the two populations. Such loci can be identified from genome sequencing by their high values of FST. This raises the question of what is the strength of selection on locally-adaptive alleles. To answer this question we analyse a 1-locus 2-allele model of a population distributed between two niches. We show by simulation of selected cases that the outputs from finite-population models are essentially the same as those from deterministic infinite-population models. We then derive theory for the infinite-population model showing the dependence of selection coefficients on equilibrium allele frequencies, migration rates, dominance and relative population sizes in the two niches. An Excel spreadsheet is provided for the calculation of selection coefficients and their approximate standard errors from observed values of population parameters. We illustrate our results with a worked example, with graphs showing the dependence of selection coefficients on equilibrium allele frequencies, and graphs showing how FST depends on the selection coefficients acting on the alleles at a locus. Given the extent of recent progress in ecological genomics, we hope our methods may help those studying migration-selection balance to quantify the advantages conferred by adaptive genes.

FST estimates of 94 populations in China based on STR markers.

The proper assessment of DNA evidence in cases of personal identification is a recurring theme in forensics. It is common practice to evaluate the strength of DNA evidence using the likelihood ratio (LR). The accurate use of population allele frequencies is a crucial problem in LR calculation. Allele frequency differences among different populations could be estimated by the FST values. Thus, FST would also affect LR values by correcting the allele frequencies. In this study, Chinese population allele frequency data were selected from population reports published in Chinese and English journals. The population-specific FST values of each population, the overall FST values of each province, each region, and the whole country, and the locus-specific FST values of different loci were calculated. The LRs using different allele frequencies and different FST values were compared based on the combination of simulated genotypes. As a result, the FST values of 94 populations, 19 provinces, 7 regions, and the whole country were obtained. The LR was overestimated using allele frequencies of the combined population containing multiple populations rather than using allele frequencies of a population, and the LRs after FST correction were lower than those without correction. Conclusively, the correction in conjunction with corresponding FST values can make the LRs more accurate and reasonable.

Comparative analyses of copy number variations between swamp buffaloes and river buffaloes.

As an important source of genomic variation, copy number variation (CNV) contributes to environmental adaptation in worldwide buffaloes. Despite this importance, CNV divergence between swamp buffaloes and river buffaloes has not been studied previously. Here, we report 21 152 CNV regions (CNVRs) in 141 buffaloes of 20 breeds detected through multiple CNV calling strategies. Only 248 CNVRs were shared between river buffalo and swamp buffalo, reflecting great variation of CNVRs between the two subspecies. Population structure analysis based on CNVs successfully separated the two buffalo subspecies. We further assessed CNV divergence by calculating FST for genome-wide CNVs. Totally, we identified 110 significantly divergent CNV segments and 44 putatively selected genes between river buffaloes and swamp buffaloes. In particular, LALBA, a key gene controlling milk production in cattle, presented a highly differentiated CNV in the promoter region, which makes it a strong functional candidate gene for differences between swamp buffaloes and river buffaloes in traits related to milk production. Our study provides useful information of CNVs in buffaloes, which may help explain the genetic differences between the two subspecies.

Lymphoid enhancer binding factor 1 is associated with nose color in Yunong black pigs.

Yunong black pig is an indigenous black pig breed being cultivated that has a pure black whole body. However, some individuals appear with a white spot on the nose. We performed case-control association studies and FST approaches in 76 animals with nose color records (26 white-nosed pigs vs. 50 black-nosed pigs) by Illumina Porcine SNP50 BeadChip data. In total, 76 SNPs, which included 2 genome-wide significant SNPs and 18 chromosome-wide suggestive SNPs, were identified by association study. The top-ranked 0.1% windows of FST results as signals under selection and 24 windows were selected. The lymphoid enhancer binding factor 1 was identified as candidate gene with strong signal in analyses of genome-wide association study and FST in black- and white-nosed pigs. Overall, our findings provide evidence that nose color is a heritable trait influenced by many loci. The results contribute to expand our understanding of pigmentation in pigs and provide SNP markers for skin color and related traits selection in Yunong black pigs. Additional research on the genetic link between nose pigmentation is needed.

Exploring genetic diversity of wild and related tetraploid wheat species Triticum turgidum and Triticum timopheevii.

INTRODUCTION: The domestication bottleneck has reduced genetic diversity inwheat, necessitating the use of wild relatives in breeding programs. Wild tetraploid wheat are widely used in the breeding programs but with morphological characters, it is difficult to distinguish these, resulting in misclassification/mislabeling or duplication of accessions in the Gene bank. OBJECTIVES: The study aims to exploreGenotyping by sequencing (GBS) to characterize wild and domesticated tetraploid wheat accessions to generate a core set of accessions to be used in the breeding program. METHODS: TASSEL-GBS pipeline was used for SNP discovery, fastStructure was used to determine the population structure and PowerCore was used to generate a core sets. Nucleotide diversity matrices of Nie's and F-statistics (FST) index were used to determine the center of genetic diversity. RESULTS: We found 65 % and 47 % duplicated accessions in Triticum timopheevii and T. turgidum respectively. Genome-wide nucleotide diversity and FST scan uncovered a lower intra and higher inter-species differentiation. Distinct FST regions were identified in genomic regions belonging to domestication genes: non-brittle rachis (Btr1) and vernalization (VRN-1).Our results suggest that Israel, Jordan, Syria, and Lebanonas the hub of genetic diversity of wild emmer;Turkey, and Georgia for T. durum; and Iraq, Azerbaijan, and Armenia for theT. timopheevii. Identified core set accessions preserved more than 93 % of the available genetic diversity. Genome wide association study (GWAS) indicated the potential chromosomal segment for resistance to leaf rust in T. timopheevii. CONCLUSION: The present study explored the potential of GBS technology in data reduction while maintaining the significant genetic diversity of the species. Wild germplasm showed more differentiation than domesticated accessions, indicating the availability of sufficient diversity for crop improvement. With reduced complexity, the core set preserves the genetic diversity of the gene bank collections and will aid in a more robust characterization of wild germplasm.

Resequencing the Yaroslavl cattle genomes reveals signatures of selection and a rare haplotype on BTA28 likely to be related to breed phenotypes.

The genomes of local livestock could shed light on their genetic history, mechanisms of adaptations to environments and unique genetics. Herein we look into the genetics and adaptations of the Russian native dairy Yaroslavl cattle breed using 22 resequenced individuals and comparing them with two related breeds (Russian Kholmogory and Holstein), and to the taurine set of the 1000 Bull Genomes Project (Run 9). HapFLK analysis with Kholmogory and Holstein breeds (using Yakut cattle as outgroup) resulted in 22 regions under selection (q-value < 0.01) on 11 chromosomes assigned to Yaroslavl cattle, including a strong signature of selection in the region of the KIT gene on BTA6. The FST (fixation index) with the 1000 Bull Genomes Dataset showed 48 non-overlapping top (0.1%) FST regions of which three overlapped HapFLK regions. We identified 1982 highly differentiated (FST > 0.40) missense mutations in the Yaroslavl genomes. These genes were enriched in the epidermal growth factor and calcium-binding functional categories. The top FST intervals contained eight genes with allele frequencies quite different between the Yaroslavl and Kholmogory breeds and the rest of the 1000 Bull Genomes Dataset, including KAT6B, which had a nearly Yaroslavl breed-specific deleterious missense mutation with the highest FST in our dataset (0.99). This gene is a part of a long haplotype containing other genes from FST and hapFLK analyses and with a negative association with weight and carcass traits according to the genotyping of 30 phenotyped Yaroslavl cattle individuals. Our work provides the industry with candidate genetic variants to be focused on in breed improvement efforts.

FSTruct: An FST -based tool for measuring ancestry variation in inference of population structure.

In model-based inference of population structure from individual-level genetic data, individuals are assigned membership coefficients in a series of statistical clusters generated by clustering algorithms. Distinct patterns of variability in membership coefficients can be produced for different groups of individuals, for example, representing different predefined populations, sampling sites or time periods. Such variability can be difficult to capture in a single numerical value; membership coefficient vectors are multivariate and potentially incommensurable across predefined groups, as the number of clusters over which individuals are distributed can vary among groups of interest. Further, two groups might share few clusters in common, so that membership coefficient vectors are concentrated on different clusters. We introduce a method for measuring the variability of membership coefficients of individuals in a predefined group, making use of an analogy between variability across individuals in membership coefficient vectors and variation across populations in allele frequency vectors. We show that in a model in which membership coefficient vectors in a population follow a Dirichlet distribution, the measure increases linearly with a parameter describing the variance of a specified component of the membership vector and does not depend on its mean. We apply the approach, which makes use of a normalized FST statistic, to data on inferred population structure in three example scenarios. We also introduce a bootstrap test for equivalence of two or more predefined groups in their level of membership coefficient variability. Our methods are implemented in the r package FSTruct.

Limited evidence of a genetic basis for sex determination in the common creek chub, Semotilus atromaculatus.

Sexual reproduction is almost universal in vertebrates; therefore, each animal species which uses it must have a mechanism for designating sex as male or female. Fish, especially, have a wide range of sex determining systems. In the present study, we aimed to identify a genetic basis for sex determination in the common creek chub (Semotilus atromaculatus) using genotyping-by-sequencing data. No sex-associated markers were found by RADSex or a GWAS using GEMMA; however, Weir and Cockerham locus-specific FST analysis and discriminant analysis of principal components revealed genetic differentiation between the sexes at several loci. While no explicit sex determination mechanism has been yet discovered in creek chub, these loci are potential candidates for future studies. Incompatible systems are thought to increase reproductive isolation but interspecific hybridization is common among groups such as cyprinid minnows; thus, studies such as ours can provide insight into hybridization and evolutionary diversification of this clade. We also highlight technical challenges involved in studying sex determination in evolutionary groups with extremely variable mechanisms and without heteromorphic sex chromosomes.

Comparative genomics and signatures of selection in North Atlantic eels.

Comparative genomic approaches can identify putative private and shared signatures of selection. We performed a comparative genomic study of North Atlantic eels, European eel (Anguilla Anguilla) and American eel (A. rostrata). The two sister species are nearly undistinguishable at the phenotypic level and despite a wide non-overlapping continental distribution, they spawn in partial sympatry in the Sargasso Sea. Taking advantage of the newly assembled and annotated genome, we used genome wide RAD sequencing data of 359 individuals retrieved from Sequence Nucleotide Archive and state-of-the-art statistic tests to identify putative genomic signatures of selection in North Atlantic eels. First, using the FST and XP-EHH methods, we detected apparent islands of divergence on a total of 7 chromosomes, particularly on chromosomes 6 and 10. Gene ontology analyses suggested candidate genes mainly related to energy production, development and regulation, which could reflect strong selection on traits related to eel migration and larval duration time. Gene effect prediction using SNPeff showed a high number of SNPs in noncoding regions, pointing to a possible regulatory role. Second, using the iHS method we detected shared regions under selection on a total of 11 chromosomes. Several hypotheses might account for the detection of shared islands of selection in North Atlantic eels, including parallel evolution due to adaptation to similar environments and introgression. Future comparative genomic studies will be needed to further clarify the causes and consequences of introgression, including the directionality of these introgression events.

The genetic basis and adult reproductive consequences of developmental thermal plasticity.

Increasing temperature and thermal variability generate profound selection on populations. Given the fast rate of environmental change, understanding the role of plasticity and genetic adaptation in response to increasing temperatures is critical. This may be especially true for thermal effects on reproductive traits in which thermal fertility limits at high temperatures may be lower than for survival traits. Consequences of changing environments during development on adult phenotypes may be particularly problematic for core traits such as reproduction that begin early in development. Here we examine the consequences of developmental thermal plasticity on subsequent adult reproductive traits and its genetic basis. We used a panel of Drosophila melanogaster (the Drosophila Genetic Reference Panel; DGRP) in which male fertility performance was previously defined as either showing relatively little (status = 'high'-performing lines) or substantial ('low'-performing lines) decline when exposed to increasing developmental temperatures. We used a thermal reaction norm approach to quantify variation in the consequences of developmental thermal plasticity on multiple adult reproductive traits, including sex-specific responses, and to identify candidate genes underlying such variation. Developmental thermal stress impacted the means and thermal reaction norms of all reproductive traits except offspring sex ratio. Mating success declined as temperature increased with no difference between high and low lines, whereas increasing temperature resulted in declines for both male and female fertility and productivity but depended on line status. Fertility and offspring number were positively correlated within and between the sexes across lines, but males were more affected than females. We identified 933 SNPs with significant evolved genetic differentiation between high and low lines. In all, 54 of these lie within genomic windows of overall high differentiation, have significant effects of genotype on the male thermal reaction norm for productivity and are associated with 16 genes enriched for phenotypes affecting reproduction, stress responses and autophagy in Drosophila and other organisms. Our results illustrate considerable plasticity in male thermal limits on several reproductive traits following development at high temperature, and we identify differentiated loci with relevant phenotypic effects that may contribute to this population variation. While our work is on a single population, phenotypic results align with an increasing number of studies demonstrating the potential for stronger selection of thermal stress on reproductive traits, particularly in males. Such large fitness costs may have both short- and long-term consequences for the evolution of populations in response to a warming world.

Molecular phylogeny and taxonomy of the Hydrangeaserrata complex (Hydrangeaceae) in western Japan, including a new subspecies of H.acuminata from Yakushima.

According to the contemporary classification of Hydrangea native to Japan, H.serrata is a polymorphic species including six varieties. We discovered a plant identified as H.serrata, but morphologically distinct from previously known varieties, in Yakushima island where approximately 50 endemic species are known. To determine the relationship of this plant with previously known varieties, we examined morphology and constructed a highly resolved phylogeny of H.serrata and its relatives using three chloroplast genomic regions, rbcL, trnL intron, psbA-trnH, and two nuclear genomic regions, ITS1 and ITS2, and Multiplex ISSR genotyping by sequencing (MIG-seq). Based on these morphological and phylogenetic observations, we describe Hydrangeaacuminatasubsp.yakushimensissubsp. nov. as a newly discovered lineage in Yakushima, Japan and propose Hydrangeaminamitanii stat. nov. and Hydrangeaacuminatasubsp.australisstat. nov. which were previously treated as varieties of H.serrata.

Genetic differentiation of grain, fodder and pod vegetable type cowpeas (Vigna unguiculata L.) identified through single nucleotide polymorphisms from genotyping-by-sequencing.

The species Vigna unguiculata L. (Walp), commonly known as cowpea, is a multi-purpose legume that has been selected into three subspecies that are divided into grain, fodder and pod (yardlong bean) types. However, genetic bases for distinctions are not well understood. The purpose of this study was to apply genotyping-by-sequencing (GBS) and current reference genome for V. unguiculata to distinguish three subspecies and identify signatures of divergence. The collection of 130 accessions included 128 cultivated from: 1) ssp. cylindrica, fodder type; 2) ssp. sesquipedalis, pod vegetable type; and 3) ssp. unguiculata, grain type. Two wilds genotypes from spp. dekindtiana and spp. pubescens, were used to anchor phylogeny. A total of 11,083 highly informative single nucleotide polymorphisms (SNPs) were discovered. Wild accessions showed distinct genetic fingerprints and were separated from cultivated subspecies. Principal component analysis showed closer relationship between ssp. unguiculata and ssp. cylindrica compared to ssp. sesquipedalis. Relative differentiation of cultivated subspecies (with Fixation Index, FST) indicated the existence of discrete signatures of selection. This work clarifies the population structure, phylogeny, and domestication of cultivated cowpeas. Furthermore, significant genetic differences between grain and pod vegetable types can provide valuable information for future breeding in three cowpea groups.

Ethnic Markers without Ethnic Conflict : Why do Interdependent Masikoro, Mikea, and Vezo of Madagascar Signal their Ethnic Differences?

People often signal their membership in groups through their clothes, hairstyle, posture, and dialect. Most existing evolutionary models argue that markers label group members so individuals can preferentially interact with those in their group. Here we ask why people mark ethnic differences when interethnic interaction is routine, necessary, and peaceful. We asked research participants from three ethnic groups in southwestern Madagascar to sort photos of unfamiliar people by ethnicity, and by with whom they would prefer or not prefer to cooperate, in a wage labor vignette. Results indicate that southwestern Malagasy reliably send and detect ethnic signals; they signal less in the marketplace, a primary site of interethnic coordination and cooperation; and they do not prefer co-ethnics as cooperation partners in novel circumstances. Results from a cultural knowledge survey and calculations of cultural FST suggest that these ethnic groups have relatively little cultural differentiation. We concur with Moya and Boyd (Human Nature 26:1-27, 2015) that ethnicity is unlikely to be a singular social phenomenon. The current functions of ethnic divisions and marking may be different from those at the moment of ethnogenesis. Group identities may persist without group conflict or differentiation.

RéSUMé: Les gens montrent souvent leur appartenance à un groupe à travers leurs modes vestimentaires, leur style de coiffure, leur posture et surtout leur dialecte. La plupart des modèles évolutifs existants soutiennent que les marqueurs caractérisent les membres du groupe afin que les individus puissent interagir de manière préférentielle avec les membres de leur groupe. Nous nous demandons ici, pourquoi les gens marquent les différences ethniques lorsque l'interaction interethnique est routinière, nécessaire et pacifique. Nous avons alors demandé à des participants issus de trois groupes ethniques du Sud-Ouest de Madagascar de trier des photos de personnes inconnues en fonction de leur appartenance ethnique, et en fonction des personnes avec lesquelles ils préféreraient ou non coopérer, dans une vignette hypothétique. Les résultats recueilli indiquent clairement que les Malgaches du Sud-Ouest émettent et détectent de manière fiable les indicateurs ethniques; ils émettent moins de signaux indicatifs sur la place du marché, dans un site primaire de coordination et de coopération interethnique; et ils ne préfèrent pas les co-ethnies comme partenaires de coopération dans des circonstances nouvelles. Basé sur les résultats obtenus d'une enquête réalisée sur les connaissances culturelles et les calculs du FST culturel suggèrent que ces groupes ethniques présentent une différenciation culturelle relativement faible. Nous partageons l'opinion de Moya et Boyd (Human Nature 26:1­27, 2015) pour dire qu'il est peu probable que l'ethnicité soit un phénomène social singulier. Les fonctions actuelles des divisions et du marquage ethniques peuvent être différentes de celles du moment de l'ethnogenèse. Les identités de groupe peuvent persister sans qu'il y ait conflit ou différenciation de groupe.

Comparing divergence landscapes from reduced-representation and whole genome resequencing in the yellow-rumped warbler (Setophaga coronata) species complex.

Researchers seeking to generate genomic data for non-model organisms are faced with a number of trade-offs when deciding which method to use. The selection of reduced representation approaches versus whole genome resequencing will ultimately affect the marker density, sequencing depth, and the number of individuals that can multiplexed. These factors can affect researchers' ability to accurately characterize certain genomic features, such as landscapes of divergence-how FST varies across the genomes. To provide insight into the effect of sequencing method on the estimation of divergence landscapes, we applied an identical bioinformatic pipeline to three generations of sequencing data (GBS, ddRAD, and WGS) produced for the same system, the yellow-rumped warbler species complex. We compare divergence landscapes generated using each method for the myrtle warbler (Setophaga coronata coronata) and the Audubon's warbler (S. c. auduboni), and for Audubon's warblers with deeply divergent mtDNA resulting from mitochondrial introgression. We found that most high-FST peaks were not detected in the ddRAD data set, and that while both GBS and WGS were able to identify the presence of large peaks, WGS was superior at a finer scale. Comparing Audubon's warblers with divergent mitochondrial haplotypes, only WGS allowed us to identify small (10-20 kb) regions of elevated differentiation, one of which contained the nuclear-encoded mitochondrial gene NDUFAF3. We calculated the cost per base pair for each method and found it was comparable between GBS and WGS, but significantly higher for ddRAD. These comparisons highlight the advantages of WGS over reduced representation methods when characterizing landscapes of divergence.

Reduced Global Genetic Differentiation of Exploited Marine Fish Species.

Knowledge on genetic structure is key to understand species connectivity patterns and to define the spatiotemporal scales over which conservation management plans should be designed and implemented. The distribution of genetic diversity (within and among populations) greatly influences species ability to cope and adapt to environmental changes, ultimately determining their long-term resilience to ecological disturbances. Yet, the drivers shaping connectivity and structure in marine fish populations remain elusive, as are the effects of fishing activities on genetic subdivision. To investigate these questions, we conducted a meta-analysis and compiled genetic differentiation data (FST/ΦST estimates) for more than 170 fish species from over 200 published studies globally distributed. We modeled the effects of multiple life-history traits, distance metrics, and methodological factors on observed population differentiation indices and specifically tested whether any signal arising from different exposure to fishing exploitation could be detected. Although the myriad of variables shaping genetic structure makes it challenging to isolate the influence of single drivers, results showed a significant correlation between commercial importance and genetic structure, with widespread lower population differentiation in commercially exploited species. Moreover, models indicate that variables commonly used as proxy for connectivity, such as larval pelagic duration, might be insufficient, and suggest that deep-sea species may disperse further. Overall, these results contribute to the growing body of knowledge on marine genetic connectivity and suggest a potential effect of commercial fisheries on the homogenization of genetic diversity, highlighting the need for additional research focused on dispersal ecology to ensure long-term sustainability of exploited marine species.

Global patterns of population genetic differentiation in seed plants.

Evaluating the factors that drive patterns of population differentiation in plants is critical for understanding several biological processes such as local adaptation and incipient speciation. Previous studies have given conflicting results regarding the significance of pollination mode, seed dispersal mode, mating system, growth form and latitudinal region in shaping patterns of genetic structure, as estimated by FST values, and no study to date has tested their relative importance together across a broad scale. Here, we assembled a 337-species data set for seed plants from publications with data on FST from nuclear markers and species traits, including variables pertaining to the sampling scheme of each study. We used species traits, while accounting for sampling variables, to perform phylogenetic multiple regressions. Results demonstrated that FST values were higher for tropical, mixed-mating, non-woody species pollinated by small insects, indicating greater population differentiation, and lower for temperate, outcrossing trees pollinated by wind. Among the factors we tested, latitudinal region explained the largest portion of variance, followed by pollination mode, mating system and growth form, while seed dispersal mode did not significantly relate to FST . Our analyses provide the most robust and comprehensive evaluation to date of the main ecological factors predicted to drive population differentiation in seed plants, with important implications for understanding the basis of their genetic divergence. Our study supports previous findings showing greater population differentiation in tropical regions and is the first that we are aware of to robustly demonstrate greater population differentiation in species pollinated by small insects.

Widespread high-risk clones of multidrug-resistant extended-spectrum ß-lactamase-producing Escherichia coli B2-ST131 and F-ST648 in public aquatic environments.

Aquatic environments are considered a reservoir for the dissemination of multidrug-resistant (MDR) bacteria, principally Escherichia coli, with the consequent spread of acquired antimicrobial resistance genes (ARGs). Widespread high-risk clones of MDR E. coli are responsible for human infections worldwide. This study aimed to characterise, through whole-genome sequencing (WGS), isolates of MDR E. coli harbouring ARGs obtained from public aquatic environments in Brazil. MDR E. coli isolates were obtained from rivers, streams and lakes that presented different Water Quality Index records and were submitted to WGS. The resistome, mobilome and virulome showed a great diversity of ARGs, plasmids and virulence genes, respectively. In addition, mutations in the quinolone resistance-determining regions of GyrA, ParC and ParE as well as several metal resistance genes (MRGs) and antibacterial biocide resistance genes (ABGs) were detected. Typing and subtyping of MDR E. coli revealed different lineages, with two belonging to widespread high-risk clones (i.e. B2-ST131-fimH30 and F-ST648-fimH27), which are grouped by core genome multilocus sequence typing (cgMLST) in clusters with E. coli lineages obtained from different sources distributed worldwide. MDR bacteria carrying MRGs and ABGs have emerged as a global human and environmental health problem. Detection of widespread high-risk clones calls for attention to the dissemination of fluoroquinolone-resistant QnrS1- and CTX-M-producing E. coli lineages associated with human infections in public aquatic environments.

Disagreement in FST estimators: A case study from sex chromosomes.

Sewall Wright developed FST for describing population differentiation and it has since been extended to many novel applications, including the detection of homomorphic sex chromosomes. However, there has been confusion regarding the expected estimate of FST for a fixed difference between the X- and Y-chromosome when comparing males and females. Here, we attempt to resolve this confusion by contrasting two common FST estimators and explain why they yield different estimates when applied to the case of sex chromosomes. We show that this difference is true for many allele frequencies, but the situation characterized by fixed differences between the X- and Y-chromosome is among the most extreme. To avoid additional confusion, we recommend that all authors using FST clearly state which estimator of FST their work uses.