Pritchard, Stephens, and Donnelly on Population Structure

A scalable approach to topic modelling in single-cell data by approximate pseudobulk projection

Probabilistic topic modelling has become essential in many types of single-cell data analysis. Based on probabilistic topic assignments in each cell, we identify the latent representation of cellular states. A dictionary matrix, consisting of topic-specific gene frequency vectors, provides interpretable bases to be compared with known cell type-specific marker genes and other pathway annotations. However, fitting a topic model on a large number of cells would require heavy computational resources-specialized computing units, computing time and memory. Here, we present a scalable approximation method customized for single-cell RNA-seq data analysis, termed ASAP, short for Annotating a Single-cell data matrix by Approximate Pseudobulk estimation. Our approach is more accurate than existing methods but requires orders of magnitude less computing time, leaving much lower memory consumption. We also show that our approach is widely applicable for atlas-scale data analysis; our method seamlessly integrates single-cell and bulk data in joint analysis, not requiring additional preprocessing or feature selection steps.

Genetic patterns reveal geographic drivers of divergence in silvereyes (Zosterops lateralis)

Identifying mechanisms that drive population divergence under varying geographic and ecological scenarios can inform our understanding of evolution and speciation. In particular, analysis of genetic data from island populations with known colonisation timelines allows us to identify potential source populations of diverging island subspecies and current relationships among populations. Silvereyes (Zosterops lateralis) are a small passerine that have served as a valuable study system to investigate evolutionary patterns on both large and small geographic scales. We examined genetic relatedness and diversity of two silvereye subspecies, the mainland Z. l. cornwalli and island Z. l. chlorocephalus, and used 18 077 single nucleotide polymorphisms (SNPs), to compare locations across southeast Queensland, Australia. Although silvereyes are prolific island colonisers our findings revealed population divergence over relatively small spatial scales was strongly influenced by geographic isolation mediated by water barriers. Strong genetic connectivity was displayed between mainland sites, but minimal inter-island connectivity was shown despite comparable sampling distances. Genetic diversity analysis showed little difference in heterozygosity between island and mainland populations, but lower inbreeding scores among the island populations. Our study confirmed the range of the Z. l. chlorocephalus subspecies throughout the southern Great Barrier Reef. Our results show that water barriers and not geographic distance per se are important in driving incipient divergence in island populations. This helps to explain the relatively high number of phenotypically differentiated, but often geographically proximate, island silvereye subspecies compared to a lower number of phenotypically less well-defined Australian continental subspecies.

Contrasting effect of hybridization on genetic differentiation in three rockfish species with similar life history

Hybridization can provide evolutionary benefits (e.g., population resilience to climate change) through the introduction of adaptive alleles and increase of genetic diversity. Nevertheless, management strategies may be designed based only on the parental species within a hybrid zone, without considering the hybrids. This can lead to ineffective spatial management of species, which can directly harm population diversity and negatively impact food webs. Three species of rockfish (Brown Rockfish (Sebastes caurinus), Copper Rockfish (S. auriculatus), and Quillback Rockfish (S. maliger)) are known to hybridize within Puget Sound, Washington, but genetic data from these species are used to infer population structure in the entire genus, including in species that do not hybridize. The goal of this project was to estimate the hybridization rates within the region and determine the effect of hybridization on geographic patterns of genetic structure. We sequenced 290 Brown, Copper, and Quillback rockfish using restriction-site associated DNA sequencing (RADseq) from four regions within and outside Puget Sound, Washington. We show that (i) hybridization within Puget Sound was asymmetrical, not recent, widespread among individuals, and relatively low level within the genome, (ii) hybridization affected population structure in Copper and Brown rockfish, but not in Quillback Rockfish and (iii) after taking hybridization into account we found limited directional dispersal in Brown and Copper rockfish, and evidence for two isolated populations in Quillback Rockfish. Our results suggest that rockfish population structure is species-specific, dependent on the extent of hybridization, and cannot be inferred from one species to another despite similar life history.

Circadian rhythm and redox homeostasis candidate genes showed association with shallow elevation in Norway spruce

The analysis of genetic variation underlying local adaptation in natural populations, together with the response to different external stimuli, is currently a hot topic in forest sciences, with the aim of identifying genetic markers controlling key phenotypic traits of interest for their inclusion in restoration and breeding programs. In Europe, one of the main tree species is Norway spruce (Picea abies (L.) H.Karst.). Using the MassARRAY® platform, 568 trees from North Rhine-Westphalia (Germany) were genotyped with 94 single nucleotide polymorphisms (SNPs) related to circadian and growth rhythms, and to stress response. The association analysis of the selected markers with health status and elevation was performed using three different methods, and those identified by at least two of these were considered as high confidence associated SNPs. While just five markers showed a weak association with health condition, 32 SNPs were correlated with elevation, six of which were considered as high confidence associated SNPs, as indicated by at least two different association methods. Among these genes, thioredoxin and pseudo response regulator 1 (PRR1) are involved in redox homeostasis and ROS detoxification, APETALA2-like 3 (AP2L3), a transcription factor, is involved in seasonal apical growth, and a RPS2-like is a disease resistance gene. The function of some of these genes in controlling light-dependent reactions and metabolic processes suggests signatures of adaptation to local photoperiod and the synchronization of the circadian rhythm. This work provides new insights into the genetic basis of local adaptation over a shallow elevation gradient in Norway spruce.

A Review on Microbiological Source Attribution Methods of Human Salmonellosis: From Subtyping to Whole-Genome Sequencing

Salmonella is one of the main causes of human foodborne illness. It is endemic worldwide, with different animals and animal-based food products as reservoirs and vehicles of infection. Identifying animal reservoirs and potential transmission pathways of Salmonella is essential for prevention and control. There are many approaches for source attribution, each using different statistical models and data streams. Some aim to identify the animal reservoir, while others aim to determine the point at which exposure occurred. With the advance of whole-genome sequencing (WGS) technologies, new source attribution models will greatly benefit from the discriminating power gained with WGS. This review discusses some key source attribution methods and their mathematical and statistical tools. We also highlight recent studies utilizing WGS for source attribution and discuss open questions and challenges in developing new WGS methods. We aim to provide a better understanding of the current state of these methodologies with application to Salmonella and other foodborne pathogens that are common sources of illness in the poultry and human sectors.

A Tibetan group from Ngawa Tibetan and Qiang Autonomous Prefecture, southwest China, is rich in genetic polymorphisms at 36 autosomal STR loci and shares a complex genetic structure with other Chinese populations

The Tibetan people are ancient and populous, constituting the seventh-largest of the fifty-five ethnic minority groups in China. The Ngawa Tibetan and Qiang Autonomous Prefecture (NTQAP), situated on the border of northwest and southwest China, has its distinct group relationships. Short tandem repeat (STR) is extremely polymorphic and extensively used in the application of forensic medicine and population genetics. However, it is not clear the genetic information including linkage disequilibrium (LD) by 36 autosomal STR (A-STR) markers in the Tibetan group from NTQAP. The Tibetan population from NTQAP of southwest China was examined for 36 A-STR loci in the research. Every marker across the 36 A-STR loci was consistent with Hardy-Weinberg equilibrium (HWE). The results of the calculation revealed that the total discrimination power (TDP) is 1-2.2552 × 10-42 and the cumulative probability of exclusion (CPE) is 1-1.3031 × 10-16. Subsequently, a total of 345 alleles with allelic frequencies ranging from 0.00382 to 0.55343 were identified, and the allelic numbers varied from 5 in both the TH01 and TPOX markers to 28 in the SE33 locus. The Ngawa Tibetan population, along with other Chinese populations, exhibited influences from historical factors and regional distribution, as indicated by the results of population genetics analysis. We thus first explored the genetic characteristics and correlated forensic parameters of the 36 A-STR markers in NTQAP to fill the gap in the Tibetan population. It was discovered that these 36 autosomal STR markers supplemented forensic STR databases and offered extremely valuable polymorphisms for Chinese forensic applications, such as parentage testing and personal identification. Moreover, the study would contribute additional information regarding the substructure and diversity in the Chinese population.

Molecular evidence depicts genetic divergence among Agropyron elongatum and A. cristatum accessions from gene pool of Iran

The genus Agropyron has an important role in soil protection and forage production in rangelands. The investigation utilized 37 ISSR primers, resulting in the detection of 956 loci within the A. elongatum genome and 705 loci within the A. cristatum genome. The findings revealed a high level of polymorphism, with 97% of loci in A. elongatum and 84% of loci in A. cristatum exhibiting variability. Notably, the primer (AC)8GCT emerged as a promising candidate for evaluating genetic diversity due to its ability to amplify numerous loci in both species. Using both the UPGMA algorithm and Bayesian analysis, the examined Agropyron accessions were categorized into two subgroups based on their respective species. The Q values associated with these subgroups suggested that certain accessions, namely "G16," "G19," "G20," "G21," "G22," "G23," "G24," and "G25," displayed potential admixture genomes. An analysis of molecular variance (AMOVA) underscored the significance of within-species variability, which accounted for 69% of the overall diversity, compared to between-species variability at 31%. Various genetic diversity parameters, including Na, Ne, I, He, and the number of private loci, were found to be higher in A. elongatum when compared to A. cristatum. Furthermore, Jaccard similarity coefficients ranged from 0.33 to 0.66 in A. cristatum and from 0.25 to 0.7 in A. elongatum, indicating the extent of genetic relatedness among these species. Intriguingly, the study identified two and three heterotic groups in A. cristatum and A. elongatum, respectively, which could be harnessed in the development of synthetic varieties to exploit heterosis. The results also indicated that a small proportion of ISSR loci pairs (5.2% in A. elongatum and 0.5% in A. cristatum) exhibited significant levels of linkage disequilibrium (LD) (P≤0.05), suggesting the potential utility of LD-based association mapping in Agropyron species. In conclusion, this research sheds light on the genetic diversity of Agropyron species and provides valuable insights into their potential applications in soil protection and forage production, as well as the prospects for enhancing genetic variability and heterosis in these species.

Forensic biogeographical ancestry inference: recent insights and current trends

BACKGROUND

As a powerful complement to the paradigmatic DNA profiling strategy, biogeographical ancestry inference (BGAI) plays a significant part in human forensic investigation especially when a database hit or eyewitness testimony are not available. It indicates one's biogeographical profile based on known population-specific genetic variations, and thus is crucial for guiding authority investigations to find unknown individuals. Forensic biogeographical ancestry testing exploits much of the recent advances in the understanding of human genomic variation and improving of molecular biology.

OBJECTIVE

In this review, recent development of prospective ancestry informative markers (AIMs) and the statistical approaches of inferring biogeographic ancestry from AIMs are elucidated and discussed.

METHODS

We highlight the research progress of three potential AIMs (i.e., single nucleotide polymorphisms, microhaplotypes, and Y or mtDNA uniparental markers) and discuss the prospects and challenges of two methods that are commonly used in BGAI.

CONCLUSION

While BGAI for forensic purposes has been thriving in recent years, important challenges, such as ethics and responsibilities, data completeness, and ununified standards for evaluation, remain for the use of biogeographical ancestry information in human forensic investigations. To address these issues and fully realize the value of BGAI in forensic investigation, efforts should be made not only by labs/institutions around the world independently, but also by inter-lab/institution collaborations.

Consilience Across Multiple, Independent Genomic Data Sets Reveals Species in a Complex with Limited Phenotypic Variation

Species delimitation in the genomic era has focused predominantly on the application of multiple analytical methodologies to a single massive parallel sequencing (MPS) data set, rather than leveraging the unique but complementary insights provided by different classes of MPS data. In this study, we demonstrate how the use of two independent MPS data sets, a sequence capture data set and a single-nucleotide polymorphism (SNP) data set generated via genotyping-by-sequencing, enables the resolution of species in three complexes belonging to the grass genus Ehrharta, whose strong population structure and subtle morphological variation limit the effectiveness of traditional species delimitation approaches. Sequence capture data are used to construct a comprehensive phylogenetic tree of Ehrharta and to resolve population relationships within the focal clades, while SNP data are used to detect patterns of gene pool sharing across populations, using a novel approach that visualizes multiple values of K. Given that the two genomic data sets are independent, the strong congruence in the clusters they resolve provides powerful ratification of species boundaries in all three complexes studied. Our approach is also able to resolve a number of single-population species and a probable hybrid species, both of which would be difficult to detect and characterize using a single MPS data set. Overall, the data reveal the existence of 11 and five species in the E. setacea and E. rehmannii complexes, with the E. ramosa complex requiring further sampling before species limits are finalized. Despite phenotypic differentiation being generally subtle, true crypsis is limited to just a few species pairs and triplets. We conclude that, in the absence of strong morphological differentiation, the use of multiple, independent genomic data sets is necessary in order to provide the cross-data set corroboration that is foundational to an integrative taxonomic approach. [Species delimitation; genotyping-by-sequencing; population structure; integrative taxonomy; cryptic species; Ehrharta (Poaceae).].

The genetic structure and demographic history of Zabelia tyaihyonii, endemic to Korean limestone karst forests, based on genome-wide SNP markers

Similar to the global phenomenon, many plant species endemic to Korean limestone karst forests are at risk of extinction due to human intervention. Zabelia tyaihyonii is a familiar shrub, called "Hardy abelia" and "Fragrant abelia" growing in the karst forests of Korea, where it is one of the most threatened species. We investigated the genetic structure and demographic history of Z. tyaihyonii, which allow us to develop appropriate conservation and management strategies. The genetic structure was evaluated using a total of 187 samples from 14 populations, covering the entire distribution of Z. tyaihyonii in South Korea. We utilized 254 and 1753 SNP loci obtained via MIG-seq (Multiplexed ISSR Genotyping by sequencing) for structure and demographic analyses, respectively. The population demographic modeling was performed with site frequency spectrum. To gain further historical insights, we also employed ENM (Ecological Niche Modeling). We found two distinct clusters (CLI and CLII) of ancient origin (ca. 490 ka). Despite CLII experiencing a more severe bottleneck, both clusters showed similar levels of genetic diversity, indicating mutual historical gene flow. Their historical distribution range seems to have changed very little. We proposed a historical distribution scenario for Z. tyaihyonii, taking into account its intrinsic factors, and emphasized a more complex response to Quaternary climate change beyond simple allopatric speciation models. These findings provide valuable insights for conservation and management strategies for Z. tyaihyonii.

Genetic Differentiation of Aedes aegypti (Diptera: Culicidae) in Areas with High Rates of Infestation in Mid-North Region of Brazil

Aedes aegypti is the principal vector of the arboviruses-yellow fever, dengue virus, chikungunya, and zika virus. Given the epidemiological importance of this mosquito, its capacity to adapt to different habitats, and its resistance to many types of control measures, systematic research into the genetic variability of the populations of this mosquito is one of the most important steps toward a better understanding of its population structure and vector competence. In this context, the present study verified the presence of distinct genetic lineages of Ae. aegypti in areas with high infestation rates, based on the analysis of microsatellite markers. The samples were collected in nine municipalities with high building infestation rates in the Mid-North region of Brazil. Six microsatellite loci were genotyped in the 138 samples, producing a total of 32 alleles, varying from one to nine alleles per locus in each of the different populations. The AMOVA revealed greater within-population genetic differentiation with high fixation rates. The general analysis of population structure, based on a Bayesian approach, revealed K = 2, with two Ae. aegypti lineages that were highly differentiated genetically. These data on the connectivity of the populations and the genetic isolation of the lineages provide important insights for the development of innovative strategies for the control of the populations of this important disease vector.

The evolution of the hypotetraploid Catolobus pendulus genome - the poorly known sister species of Capsella

The establishment of Arabidopsis as the most important plant model has also brought other crucifer species into the spotlight of comparative research. While the genus Capsella has become a prominent crucifer model system, its closest relative has been overlooked. The unispecific genus Catolobus is native to temperate Eurasian woodlands, from eastern Europe to the Russian Far East. Here, we analyzed chromosome number, genome structure, intraspecific genetic variation, and habitat suitability of Catolobus pendulus throughout its range. Unexpectedly, all analyzed populations were hypotetraploid (2n = 30, ~330 Mb). Comparative cytogenomic analysis revealed that the Catolobus genome arose by a whole-genome duplication in a diploid genome resembling Ancestral Crucifer Karyotype (ACK, n = 8). In contrast to the much younger Capsella allotetraploid genomes, the presumably autotetraploid Catolobus genome (2n = 32) arose early after the Catolobus/Capsella divergence. Since its origin, the tetraploid Catolobus genome has undergone chromosomal rediploidization, including a reduction in chromosome number from 2n = 32 to 2n = 30. Diploidization occurred through end-to-end chromosome fusion and other chromosomal rearrangements affecting a total of six of 16 ancestral chromosomes. The hypotetraploid Catolobus cytotype expanded toward its present range, accompanied by some longitudinal genetic differentiation. The sister relationship between Catolobus and Capsella allows comparative studies of tetraploid genomes of contrasting ages and different degrees of genome diploidization.

Climate change in the Arctic: Testing the poleward expansion of ticks and tick-borne diseases

Climate change is most strongly felt in the polar regions of the world, with significant impacts on the species that live there. The arrival of parasites and pathogens from more temperate areas may become a significant problem for these populations, but current observations of parasite presence often lack a historical reference of prior absence. Observations in the high Arctic of the seabird tick Ixodes uriae suggested that this species expanded poleward in the last two decades in relation to climate change. As this tick can have a direct impact on the breeding success of its seabird hosts and vectors several pathogens, including Lyme disease spirochaetes, understanding its invasion dynamics is essential for predicting its impact on polar seabird populations. Here, we use population genetic data and host serology to test the hypothesis that I. uriae recently expanded into Svalbard. Both black-legged kittiwakes (Rissa tridactyla) and thick-billed murres (Uria lomvia) were sampled for ticks and blood in Kongsfjorden, Spitsbergen. Ticks were genotyped using microsatellite markers and population genetic analyses were performed using data from 14 reference populations from across the tick's northern distribution. In contrast to predictions, the Spitsbergen population showed high genetic diversity and significant differentiation from reference populations, suggesting long-term isolation. Host serology also demonstrated a high exposure rate to Lyme disease spirochaetes (Bbsl). Targeted PCR and sequencing confirmed the presence of Borrelia garinii in a Spitsbergen tick, demonstrating the presence of Lyme disease bacteria in the high Arctic for the first time. Taken together, results contradict the notion that I. uriae has recently expanded into the high Arctic. Rather, this tick has likely been present for some time, maintaining relatively high population sizes and an endemic transmission cycle of Bbsl. Close future observations of population infestation/infection rates will now be necessary to relate epidemiological changes to ongoing climate modifications.

Exploring genomic regions involved in bread wheat resistance to leaf rust at seedling/adult stages by using GWAS analysis

BACKGROUND

Global wheat productivity is seriously challenged by a range of rust pathogens, especially leaf rust derived from Puccinia triticina. Since the most efficient approach to control leaf rust is genetic resistance, many efforts have been made to uncover resistance genes; however, it demands an ongoing exploration for effective resistance sources because of the advent of novel virulent races. Thus, the current study was focused on detecting leaf rust resistance-related genomic loci against the P. triticina prevalent races by GWAS in a set of Iranian cultivars and landraces.

RESULTS

Evaluation of 320 Iranian bread wheat cultivars and landraces against four prevalent rust pathotypes of P. triticina (LR-99-2, LR-98-12, LR-98-22, and LR-97-12) indicated the diversity in wheat accessions responses to P. triticina. From GWAS results, 80 leaf rust resistance QTLs were located in the surrounding known QTLs/genes on almost chromosomes, except for 1D, 3D, 4D, and 7D. Of these, six MTAs (rs20781/rs20782 associated with resistance to LR-97-12; rs49543/rs52026 for LR-98-22; rs44885/rs44886 for LR-98-22/LR-98-1/LR-99-2) were found on genomic regions where no resistance genes previously reported, suggesting new loci conferring resistance to leaf rust. The GBLUP genomic prediction model appeared better than RR-BLUP and BRR, reflecting that GBLUP is a potent model for genomic selection in wheat accessions.

CONCLUSIONS

Overall, the newly identified MTAs as well as the highly resistant accessions in the recent work provide an opportunity towards improving leaf rust resistance.

Genomics reveals a genetically isolated population of the Pacific white-sided dolphin (Lagenorhynchus obliquidens) distributed in the Sea of Japan

The processes by which animals become genetically isolated in an open environment such as the ocean have not yet been fully elucidated. Morphologically different populations of Pacific white-sided dolphin Lagenorhynchus obliquidens are observed sympatrically in the Sea of Japan. However, genetic studies that exclusively used limited mitochondrial loci or microsatellite DNA have failed to demonstrate the existence of genetically distinct populations. Here, to reveal the population structure, we analysed genome-wide population genetic data using single nucleotide polymorphisms (SNPs) gathered in 2018-2021 from all domestic captive individuals in aquaria, the majority of which originated from the wild, as well as from some stranded individuals, together covering a wide area of coastal water around Japan (n = 123). Multiplexed intersimple sequence repeat genotyping-by-sequencing analysis was performed to obtain the SNP data. Principal coordinate analysis and the clustering method structure indicated that two genetically-distinct populations exist, with little interpopulation gene flow revealed. In addition, the genotypic segregation was reflected in differences in external morphotype. Furthermore, a population demographic analysis based on the whole-genome sequences of an individual from each population indicated that sea-level changes during the Last Glacial Period probably led to allopatric divergence of this species in a limited area of the Sea of Japan, with that group subsequently sharing a distribution area with the other population. These findings yield insights into the formation of genetically isolated sympatric populations in the ocean.

Characterization of experimental cattle populations from Argentina with a low-density SNP genotyping panel

The objectives of the present experiment were to evaluate a low-density SNP array designed for the molecular characterisation of gene banks and to assess the genetic diversity and population structure of beef cattle herds from an Argentinean research station. Forty-nine animals from three breeds (Angus, Hereford, and Argentinean Creole) were genotyped using the multi-species IMAGE001 60-K SNP array (10 K for cattle). Genotypes of other 19 cattle populations from Argentina, other American countries, and Europe were included in the study. Of special interest was the characterization of the Argentinean Creole, the only autochthonous cattle breed in the country. Due to the merging of different datasets, approximately 5 K SNPs were effectively used. Genetic differentiation (F_ST ), principal component analysis, neighbour-joining tree of Reynolds distances and ancestry analysis showed that autochthonous American breeds are clearly differentiated, but all have genetic influences of Iberian cattle. The analysed herds of Argentinean Creole showed no evidence of recent admixture and represent a unique genetic pool within local American breeds. An experimental herd and the local commercial Hereford population have also diverged, probably due to the influence of current selection objectives in the breed. Our results illustrate the utility of using low-cost, low density SNP arrays in the evaluation of animal genetic resources. This type of panels could become a very useful resource in developing countries, where most endangered cattle breeds are located. The results also reinforce the importance of experimental herds as reservoir of genetic diversity, particularly in the case of local breeds under-represented in traditional production systems.

Genetic diversity in global populations of the critically endangered addax (Addax nasomaculatus) and its implications for conservation

Threatened species are frequently patchily distributed across small wild populations, ex situ populations managed with varying levels of intensity and reintroduced populations. Best practice advocates for integrated management across in situ and ex situ populations. Wild addax (Addax nasomaculatus) now number fewer than 100 individuals, yet 1000 of addax remain in ex situ populations, which can provide addax for reintroductions, as has been the case in Tunisia since the mid-1980s. However, integrated management requires genetic data to ascertain the relationships between wild and ex situ populations that have incomplete knowledge of founder origins, management histories, and pedigrees. We undertook a global assessment of genetic diversity across wild, ex situ and reintroduced populations in Tunisia to assist conservation planning for this Critically Endangered species. We show that the remnant wild populations retain more mitochondrial haplotypes that are more diverse than the entirety of the ex situ populations across Europe, North America and the United Arab Emirates, and the reintroduced Tunisian population. Additionally, 1704 SNPs revealed that whilst population structure within the ex situ population is minimal, each population carries unique diversity. Finally, we show that careful selection of founders and subsequent genetic management is vital to ensure genetic diversity is provided to, and minimize drift and inbreeding within reintroductions. Our results highlight a vital need to conserve the last remaining wild addax population, and we provide a genetic foundation for determining integrated conservation strategies to prevent extinction and optimize future reintroductions.

A box on the river: The phylogenetics and phylogeography of Eucalyptus baueriana (Eucalyptus sect. Adnataria ser. Heterophloiae)

We present a phylogeographic study of the tree species Eucalyptus baueriana Schauer, which occurs in disjunct areas on the near coastal plains and ranges of the south-east Australian mainland. DArTseq data are used to build a phylogeny including E. baueriana and closely related taxa to test its monophyly, test the genetic distinctness of the three subspecies of E. baueriana, and investigate relationships between its disjunct populations. Additionally, we use population structure analysis to investigate the genetic distinctness of populations, and MaxEnt to investigate the environmental factors potentially influencing the species' distribution. We show E. baueriana is monophyletic and most closely related to three other Blue Box eucalypt species: E. conica H.Deane & Maiden, E. dalveenica T.L.Collins, R.L.Andrew & J.J.Bruhl and E. magnificata L.A.S.Johnson & K.D.Hill, with some evidence for genetic introgression between these taxa. Within E. baueriana, the deepest genetic breaks do not correspond with the subspecies classification as the two geographically restricted subspecies, together with samples of the more widespread E. baueriana subsp. baueriana from west of the Gippsland lowlands, form a south-western clade with that is sister to other populations of subsp. baueriana. The oldest genetic break in the species occurs in far eastern Gippsland (Victoria), corresponding to one of the shortest geographic disjunctions in the species' distribution. Genetic breaks in other species have been observed in this region which is broadly referred to as the southern transition zone. Both total annual rainfall and the seasonality of this rainfall are hypothesised to affect the species' distribution; gaps in its distribution are in areas of higher rainfall that support closed forest and in regions with more winter dominated rainfall.

Maximizing genetic gain through unlocking genetic variation in different ecotypes of kalmegh (Andrographis paniculata (Burm. f.) Nee)

Andrographis paniculata, commonly known as kalmegh is among the most popular medicinal herbs in Southeast Asia. It is widely cultivated for medicinal purposes. The bioactive molecule, Andrographolide accumulated in herb leaves has immense therapeutic and economic potential. However, comprehensive information regarding genetic diversity is very limited in this species. The present study assessed genetic diversity between and within the six populations (ecotypes) of twenty-four kalmegh accessions using multiple datasets (agro-morphological traits, phytochemical traits, and genic markers). This is the established report where EST-SSR (Expressed sequence tags-Simple Sequence Repeat) markers have been used to unlock genetic variation in kalmegh. Here, we identified and developed ninety-one metabolic pathway-specific EST-SSR markers. Finally, 32 random EST-SSR primer pairs were selected for genetic diversity assessment. Multivariate analysis to unveil the agro-morphological, phytochemical and genotypic variability was helpful in discriminating various germplasms studied in the present study. Among all the morphological discriptors used in present study, days to fifty percent flowering and dry herb yield were found as potential selection index for AP genetic improvement. Hierarchical cluster analysis built with agro-morphological data identified three major groups. However, corresponding analysis with phytochemical and molecular data generated two clear-cut groups among the studied individuals. Moreover, the grouping of individuals into different clusters using multiple datasets was geographically independent, and also showed inconsistency in grouping among agromorphological, phytochemical and molecular dataset based clusters. However, joint analysis using agro-morphological, phytochemical and genotypic information generated two genetic groups, which could be a valuable resource for identifying complementary crossing panels in the kalmegh breeding program. The accessions AP7, AP13, AP5, AP3 belong to cluster I and accessions AP17, AP18 belong to cluster II could be utilized as potential donors for high dry herb yield and andrographolide content, respectively in different selective breeding programs of AP. Thus, our results provided useful information about the overall genetic diversity and variation in economic traits useful for initiating selective breeding programs for contrasting traits of interest and maximizing genetic gain in kalmegh.

Exploring the genetic makeup and population structure among Capsicum accessions for crop improvement and breeding curriculum insights

Background

Capsicum or chilli is an important crop in India which exhibits immense structural and genetic variations reflecting their intra- and inter-specific relationships. The aim of this study was to establish relationships amongst 54 Capsicum accessions through analysis of genetic and population structure using ISSR markers.

Results

Out of 19, successful DNA amplifications were shown by 7 ISSR primers and a total of 80 bands were identified ranging between 8 and 14 with an average of 11.43 bands/primer. A significant degree of polymorphic information content (PIC), discriminating power (DP), resolving power (RP), effective multiplex ratio (EMR), and marker index (MI) were identified as 0.39, 0.70, 6.40, 5.88, and 2.30, respectively, using ISSR markers in chillies. The cross-transferability ranged from 8.0 to 72.15% with an average of 52.63% among chillies. Amongst genetic information, grand mean values were 0.264, 0.180, 0.376, 0.296, and 0.180, which correspond to Shannon’s information index (I), expected heterozygosity (He), Nei’s gene diversity, total diversity among species (Ht), diversity within species (Hs), respectively. Further, the coefficients of gene differentiation (Gst) and gene flow (Nm) were 0.393 and 0.773, representing higher genetic variation among the population which was confirmed by analysis of molecular variance (AMOVA).

Conclusion

ISSR markers represented a potent system for the estimation of relationships or variation studies and generated information useful for planning crop management and improvement strategies in chilli breeding.

Population Dynamics of Wide Compatibility System and Evaluation of Intersubspecific Hybrids by indica-japonica Hybridization in Rice

The exploitation of heterosis through intersubspecific hybridisation between indica and japonica has been a major breeding target in rice, but is marred by the cross incompatibility between the genomes. Wide compatibility (WC) is a triallelic system at the S5 locus on chromosome 6 that ensures the specificity of hybridisation within and between indica and japonica. The S5n allele that favours intercrossing is sparsely distributed in the rice gene pool and therefore warrants identification of diverse WC sources to develop superior intersubspecific hybrids. In this study, we have identified several novel WC sources through the marker-assisted screening of a large set of 950 rice genotypes. Seventeen percent of the genotypes carried S5n, which fell into two subpopulations. The WC genotypes showed wide phenotypic and genotypic variability, including both indica and japonica lines. Based on phenotypic performance, the WC varieties were grouped into three clusters. A subset of 41 WC varieties was used to develop 164 hybrids, of which WC/japonica hybrids showed relative superiority over WC/indica hybrids. The multilocation evaluation of hybrids indicated that hybrids derived from WC varieties, such as IRG137, IRG143, OYR128, and IRGC10658, were higher yielding across all the three different locations. Most of the hybrids showed the stability of performance across locations. The identified diverse set of wide compatible varieties (WCVs) can be used in the development of intersubspecific hybrids and also for parental line development in hybrid rice breeding.

Population and Landscape Genetics Provide Insights Into Species Conservation of Two Evergreen Oaks in Qinghai-Tibet Plateau and Adjacent Regions

The combination of population and landscape genetics can facilitate the understanding of conservation strategy under the changing climate. Here, we focused on the two most diverse and ecologically important evergreen oaks: Quercus aquifolioides and Quercus spinosa in Qinghai-Tibetan Plateau (QTP), which is considered as world's biodiversity hotspot. We genotyped 1,657 individuals of 106 populations at 15 nuclear microsatellite loci throughout the species distribution range. Spatial patterns of genetic diversity were identified by mapping the allelic richness (AR) and locally common alleles (LCA) according to the circular neighborhood methodology. Migration routes from QTP were detected by historical gene flow estimation. The response pattern of genetic variation to environmental gradient was assessed by the genotype-environment association (GEA) analysis. The overall genetic structure showed a high level of intra-species genetic divergence of a strong west-east pattern. The West-to-East migration route indicated the complex demographic history of two oak species. We found evidence of isolation by the environment in Q. aqu-East and Q. spi-West lineage but not in Q. aqu-West and Q. spi-East lineage. Furthermore, priority for conservation should be given to populations that retain higher spatial genetic diversity or isolated at the edge of the distribution range. Our findings indicate that knowledge of spatial diversity and migration route can provide valuable information for the conservation of existing populations. This study provides an important guide for species conservation for two oak species by the integration of population and landscape genetic methods.

Posterior probabilities of membership of repertoires in acoustic clades

Recordings of calls may be used to assess population structure for acoustic species. This can be particularly effective if there are identity calls, produced nearly exclusively by just one population segment. The identity call method, IDcall, classifies calls into types using contaminated mixture models, and then clusters repertoires of calls into identity clades (potential population segments) using identity calls that are characteristic of the repertoires in each identity clade. We show how to calculate the Bayesian posterior probabilities that each repertoire is a member of each identity clade, and display this information as a stacked bar graph. This methodology (IDcallPP) is introduced using the output of IDcall but could easily be adapted to estimate posterior probabilities of clade membership when acoustic clades are delineated using other methods. This output is similar to that of the STRUCTURE software which uses molecular genetic data to assess population structure and has become a standard in conservation genetics. The technique introduced here should be a valuable asset to those who use acoustic data to address evolution, ecology, or conservation, and creates a methodological and conceptual bridge between geneticists and acousticians who aim to assess population structure.

Re-situations of scientific knowledge: a case study of a skirmish over clusters vs clines in human population genomics

We track and analyze the re-situation of scientific knowledge in the field of human population genomics ancestry studies. We understand re-situation as a process of accommodating the direct or indirect transfer of objects of knowledge from one site/situation to (one or many) other sites/situations. Our take on the concept borrows from Mary S. Morgan's work on facts traveling while expanding it to include other objects of knowledge such as models, data, software, findings, and visualizations. We structure a specific case study by tracking the re-situation of these objects between three research projects studying human population diversity reported in three articles in Science, Genome Research and PLoS Genetics between 2002 and 2005. We characterize these three engagements as a unit of analysis, a "skirmish," in order to compare: (a) the divergence of interests in how life-scientists answer similar research questions and (b) to track the challenging transformation of workflows in research laboratories as these scientific objects are re-situated individually or in bundles. Our analysis of the case study shows that an accurate understanding of re-situation requires tracking the whole bundle of objects in a project because they interact in particular key ways. The absence or dismissal of these interactions opens the door to unforeseen trade-offs, misunderstandings and misrepresentations about research design(s) and workflow(s) and what these say about the questions asked and the findings produced.

Genomic evidence of an ancient Inland Temperate Rainforest in the Pacific Northwest of North America

The disjunct temperate rainforests of the Pacific Northwest of North America (PNW) are characterized by late‐successional dominant tree species Thuja plicata (western redcedar) and Tsuga heterophylla (western hemlock). The demographic histories of these species, along with the PNW rainforest ecosystem in its entirety, have been heavily impacted by geological and climatic changes the PNW has experienced over the last 5 million years, including mountain orogeny and repeated Pleistocene glaciations. These environmental events have ultimately shaped the history of these species, with inland populations potentially being extirpated during the Pleistocene glaciations. Here, we collect genomic data for both species across their ranges to test multiple demographic models, each reflecting a different phylogeographical hypothesis on how the ecosystem‐dominating species may have responded to dramatic climatic change. Our results indicate that inland and coastal populations in both species diverged ~2.5 million years ago in the early Pleistocene and experienced decreases in population size during glacial cycles, with subsequent population expansion. Importantly, we found evidence for gene flow between coastal and inland populations during the mid‐Holocene. It is likely that intermittent migration in these species during this time has prevented allopatric speciation via genetic drift alone. In conclusion, our results from combining genomic data and demographic inference procedures establish that populations of the ecosystem dominants Thuja plicata and Tsuga heterophylla persisted in refugia located in both the coastal and inland regions of the PNW throughout the Pleistocene, with populations expanding and contracting in response to glacial cycles with occasional gene flow.

Global flyway evolution in red knots Calidris canutus and genetic evidence for a Nearctic refugium

Present‐day ecology and population structure are the legacies of past climate and habitat perturbations, and this is particularly true for species that are widely distributed at high latitudes. The red knot, Calidris canutus, is an arctic‐breeding, long‐distance migratory shorebird with six recognized subspecies defined by differences in morphology, migration behavior, and annual cycle phenology, in a global distribution thought to have arisen just since the last glacial maximum (LGM). We used nextRAD sequencing of 10,881 single‐nucleotide polymorphisms (SNPs) to assess the neutral genetic structure and phylogeographic history of 172 red knots representing all known global breeding populations. Using population genetics approaches, including model‐based scenario‐testing in an approximate Bayesian computation (ABC) framework, we infer that red knots derive from two main lineages that diverged ca. 34,000 years ago, and thus most probably persisted at the LGM in both Palearctic and Nearctic refugia, followed by at least two instances of secondary contact and admixture. Within two Beringian subspecies (C. c. roselaari and rogersi), we detected previously unknown genetic structure among sub‐populations sharing a migratory flyway, reflecting additional complexity in the phylogeographic history of the region. Conversely, we found very weak genetic differentiation between two Nearctic populations (rufa and islandica) with clearly divergent migratory phenotypes and little or no apparent contact throughout the annual cycle. Together, these results suggest that relative gene flow among migratory populations reflects a complex interplay of historical, geographical, and ecological factors.

Review of the Forensic Applicability of Biostatistical Methods for Inferring Ancestry from Autosomal Genetic Markers

The inference of ancestry has become a part of the services many forensic genetic laboratories provide. Interest in ancestry may be to provide investigative leads or identify the region of origin in cases of unidentified missing persons. There exist many biostatistical methods developed for the study of population structure in the area of population genetics. However, the challenges and questions are slightly different in the context of forensic genetics, where the origin of a specific sample is of interest compared to the understanding of population histories and genealogies. In this paper, the methodologies for modelling population admixture and inferring ancestral populations are reviewed with a focus on their strengths and weaknesses in relation to ancestry inference in the forensic context.

OptM: estimating the optimal number of migration edges on population trees using Treemix

The software Treemix has become extensively used to estimate the number of migration events, or edges (m), on population trees from genome-wide allele frequency data. However, the appropriate number of edges to include remains unclear. Here, I show that an optimal value of m can be inferred from the second-order rate of change in likelihood (Δm) across incremental values of m. Repurposed from its original use to estimate the number of population clusters in the software Structure (ΔK), I show using simulated populations that Δm performs equally as well as current recommendations for Treemix. A demonstration of an empirical dataset from domestic dogs indicates that this method may be preferable in large, complex population histories and can prioritize migration events for subsequent investigation. The method has been implemented in a freely available R package called “OptM” and as a web application (https://rfitak.shinyapps.io/OptM/) to interface directly with the output files of Treemix.

The impact of estimator choice: Disagreement in clustering solutions across K estimators for Bayesian analysis of population genetic structure across a wide range of empirical data sets

The software program STRUCTURE is one of the most cited tools for determining population structure. To infer the optimal number of clusters from STRUCTURE output, the ΔK method is often applied. However, a recent study relying on simulated microsatellite data suggested that this method has a downward bias in its estimation of K and is sensitive to uneven sampling. If this finding holds for empirical data sets, conclusions about the scale of gene flow may have to be revised for a large number of studies. To determine the impact of method choice, we applied recently described estimators of K to re-estimate genetic structure in 41 empirical microsatellite data sets; 15 from a broad range of taxa and 26 from one phylogenetic group, coral. We compared alternative estimates of K (Puechmaille statistics) with traditional (ΔK and posterior probability) estimates and found widespread disagreement of estimators across data sets. Thus, one estimator alone is insufficient for determining the optimal number of clusters; this was regardless of study organism or evenness of sampling scheme. Subsequent analysis of molecular variance (AMOVA) did not necessarily clarify which clustering solution was best. To better infer population structure, we suggest a combination of visual inspection of STRUCTURE plots and calculation of the alternative estimators at various thresholds in addition to ΔK. Disagreement between traditional and recent estimators may have important biological implications, such as previously unrecognized population structure, as was the case for many studies reanalysed here.

Population genetics, invasion pathways and public health risks of the raccoon and its roundworm Baylisascaris procyonis in northwestern Europe

The geographic range of the zoonotic raccoon roundworm (Baylisascaris procyonis) is expanding together with the range of its host, the raccoon (Procyon lotor). This creates a new public health risk in parts of Europe where this parasite was previously absent. In the Netherlands, a raccoon population is becoming established and incidental findings of B. procyonis have been reported. To assess the risk to public health, the prevalence of B. procyonis was determined in the province of Limburg, where currently the largest Dutch raccoon population is present, as well as in the adjoining region of southern Belgium. Furthermore, genetic methods were employed to assess invasion pathways of both the raccoon and B. procyonis to aid in the development of control measures. Macroscopic analysis of intestinal content and testing of faecal samples were performed to detect B. procyonis adults and eggs. The population genetics of both B. procyonis and its raccoon host were analysed using samples from central and northwestern Europe. B. procyonis was found in 14/23 (61%, 95% CI: 41%-78%) raccoons from Limburg, but was not detected in 50 Belgian raccoons. Genetic analyses showed that the majority of the Dutch raccoons and their roundworms were introduced through ex-captive individuals. As long as free-living raccoon populations originate from captivity, population control methods may be pursued. However, natural dispersal from the border regions will complicate prolonged population control. To reduce the public health risk posed by B. procyonis, public education to increase awareness and adapt behaviour towards raccoons is key.

Historical effective population size of North American hoary bat (Lasiurus cinereus) and challenges to estimating trends in contemporary effective breeding population size from archived samples

Background

Hoary bats (Lasiurus cinereus) are among the bat species most commonly killed by wind turbine strikes in the midwestern United States. The impact of this mortality on species census size is not understood, due in part to the difficulty of estimating population size for this highly migratory and elusive species. Genetic effective population size (Ne) could provide an index of changing census population size if other factors affecting Ne are stable.

Methods

We used the NeEstimator package to derive effective breeding population size (Nb) estimates for two temporally spaced cohorts: 93 hoary bats collected in 2009-2010 and an additional 93 collected in 2017-2018. We sequenced restriction-site associated polymorphisms and generated a de novo genome assembly to guide the removal of sex-linked and multi-copy loci, as well as identify physically linked markers.

Results

Analysis of the reference genome with psmc suggested at least a doubling of Ne in the last 100,000 years, likely exceeding Ne = 10,000 in the Holocene. Allele and genotype frequency analyses confirmed that the two cohorts were comparable, although some samples had unusually high or low observed heterozygosities. Additionally, the older cohort had lower mean coverage and greater variability in coverage, and batch effects of sampling locality were observed that were consistent with sample degradation. We therefore excluded samples with low coverage or outlier heterozygosity, as well as loci with sequence coverage far from the mode value, from the final data set. Prior to excluding these outliers, contemporary Nb estimates were significantly higher in the more recent cohort, but this finding was driven by high values for the 2018 sample year and low values for all other years. In the reduced data set, Nb did not differ significantly between cohorts. We found base substitutions to be strongly biased toward cytosine to thymine or the complement, and further partitioning loci by substitution type had a strong effect on Nb estimates. Minor allele frequency and base quality bias thresholds also had strong effects on Nb estimates. Instability of Nb with respect to common data filtering parameters and empirically identified factors prevented robust comparison of the two cohorts. Given that confidence intervals frequently included infinity as the stringency of data filtering increased, contemporary trends in Nb of North American hoary bats may not be tractable with the linkage disequilibrium method, at least using the protocol employed here.

Transcriptome sequencing and microsatellite marker discovery in Ailanthus altissima (Mill.) Swingle (Simaroubaceae)

Ailanthus altissima Swingle, is a tree species native to East Asia and has a great potential in decorative, bioenergy and industrial applications in many countries. To date, despite its commercial importance, the genomic and genetic resources available for this species are still insufficient. In this study, we characterized the transcriptome of A. altissima and developed thirteen EST-SSRs (expressed sequence tag-simple sequence repeats) based on Illumina paired-end RNA sequencing (RNA-seq). Besides, we developed ten polymorphic chloroplast microsatellite (cpSSR) markers using the available chloroplast genome of A. altissima. The transcriptome data produced 87,797 unigenes, of which 64,891 (73.91%) unigenes were successfully annotated in at least one protein database. For cpSSR markers the number of detected alleles (N) per marker varied from three at cpSSR12 to twelve at cpSSR8, the unbiased haploid diversity indices (uh) varied from 0.111 to 0.485, and haploid diversity indices (h) ranged from 0.101 to 0.444 with an average unbiased haploid diversity index (uh) of 0.274. Overall, a total of 65 different cpSSR alleles were identified at the ten loci among 165 individuals of A. altissima. The allele number per locus for EST-SSRs varied from 2.143 to 9.357, and the values of observed and expected heterozygosity ranged from 0.312 to 1.000 and 0.505 to 0.826, respectively. The molecular markers developed in this study will facilitate future genetic diversity, population structure, long distance-gene transfer and pollen-based gene flow analyses of A. altissima populations from its known distribution ranges in China focusing on planted and natural forest stands.

Exploring Population Structure with Admixture Models and Principal Component Analysis

Population structure is a commonplace feature of genetic variation data, and it has importance in numerous application areas, including evolutionary genetics, conservation genetics, and human genetics. Understanding the structure in a sample is necessary before more sophisticated analyses are undertaken. Here we provide a protocol for running principal component analysis (PCA) and admixture proportion inference-two of the most commonly used approaches in describing population structure. Along with hands-on examples with CEPH-Human Genome Diversity Panel and pragmatic caveats, readers will learn to analyze and visualize population structure on their own data.

Quantifying and contextualizing the impact of bioRxiv preprints through automated social media audience segmentation

Engagement with scientific manuscripts is frequently facilitated by Twitter and other social media platforms. As such, the demographics of a paper's social media audience provide a wealth of information about how scholarly research is transmitted, consumed, and interpreted by online communities. By paying attention to public perceptions of their publications, scientists can learn whether their research is stimulating positive scholarly and public thought. They can also become aware of potentially negative patterns of interest from groups that misinterpret their work in harmful ways, either willfully or unintentionally, and devise strategies for altering their messaging to mitigate these impacts. In this study, we collected 331,696 Twitter posts referencing 1,800 highly tweeted bioRxiv preprints and leveraged topic modeling to infer the characteristics of various communities engaging with each preprint on Twitter. We agnostically learned the characteristics of these audience sectors from keywords each user's followers provide in their Twitter biographies. We estimate that 96% of the preprints analyzed are dominated by academic audiences on Twitter, suggesting that social media attention does not always correspond to greater public exposure. We further demonstrate how our audience segmentation method can quantify the level of interest from nonspecialist audience sectors such as mental health advocates, dog lovers, video game developers, vegans, bitcoin investors, conspiracy theorists, journalists, religious groups, and political constituencies. Surprisingly, we also found that 10% of the preprints analyzed have sizable (>5%) audience sectors that are associated with right-wing white nationalist communities. Although none of these preprints appear to intentionally espouse any right-wing extremist messages, cases exist in which extremist appropriation comprises more than 50% of the tweets referencing a given preprint. These results present unique opportunities for improving and contextualizing the public discourse surrounding scientific research.

Mycorrhiza: genotype assignment using phylogenetic networks

MOTIVATION

The genotype assignment problem consists of predicting, from the genotype of an individual, which of a known set of populations it originated from. The problem arises in a variety of contexts, including wildlife forensics, invasive species detection and biodiversity monitoring. Existing approaches perform well under ideal conditions but are sensitive to a variety of common violations of the assumptions they rely on.

RESULTS

In this article, we introduce Mycorrhiza, a machine learning approach for the genotype assignment problem. Our algorithm makes use of phylogenetic networks to engineer features that encode the evolutionary relationships among samples. Those features are then used as input to a Random Forests classifier. The classification accuracy was assessed on multiple published empirical SNP, microsatellite or consensus sequence datasets with wide ranges of size, geographical distribution and population structure and on simulated datasets. It compared favorably against widely used assessment tests or mixture analysis methods such as STRUCTURE and Admixture, and against another machine-learning based approach using principal component analysis for dimensionality reduction. Mycorrhiza yields particularly significant gains on datasets with a large average fixation index (FST) or deviation from the Hardy-Weinberg equilibrium. Moreover, the phylogenetic network approach estimates mixture proportions with good accuracy.

AVAILABILITY AND IMPLEMENTATION

Mycorrhiza is released as an easy to use open-source python package at github.com/jgeofil/mycorrhiza.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Evaluation of model fit of inferred admixture proportions

Model based methods for genetic clustering of individuals, such as those implemented in structure or ADMIXTURE, allow the user to infer individual ancestries and study population structure. The underlying model makes several assumptions about the demographic history that shaped the analysed genetic data. One assumption is that all individuals are a result of K homogeneous ancestral populations that are all well represented in the data, while another assumption is that no drift happened after the admixture event. The histories of many real world populations do not conform to that model, and in that case taking the inferred admixture proportions at face value might be misleading. We propose a method to evaluate the fit of admixture models based on estimating the correlation of the residual difference between the true genotypes and the genotypes predicted by the model. When the model assumptions are not violated, the residuals from a pair of individuals are not correlated. In the case of a bad fitting admixture model, individuals with similar demographic histories have a positive correlation of their residuals. Using simulated and real data, we show how the method is able to detect a bad fit of inferred admixture proportions due to using an insufficient number of clusters K or to demographic histories that deviate significantly from the admixture model assumptions, such as admixture from ghost populations, drift after admixture events and nondiscrete ancestral populations. We have implemented the method as an open source software that can be applied to both unphased genotypes and low depth sequencing data.

Genomic analysis of natural intra-specific hybrids among Ethiopian isolates of Leishmania donovani

Parasites of the genus Leishmania (Kinetoplastida: Trypanosomatidae) cause widespread and devastating human diseases. Visceral leishmaniasis due to Leishmania donovani is endemic in Ethiopia where it has also been responsible for major epidemics. The presence of hybrid genotypes has been widely reported in surveys of natural populations, genetic variation reported in a number of Leishmania species, and the extant capacity for genetic exchange demonstrated in laboratory experiments. However, patterns of recombination and the evolutionary history of admixture that produced these hybrid populations remain unclear. Here, we use whole-genome sequence data to investigate Ethiopian L. donovani isolates previously characterized as hybrids by microsatellite and multi-locus sequencing. To date there is only one previous study on a natural population of Leishmania hybrids based on whole-genome sequences. We propose that these hybrids originate from recombination between two different lineages of Ethiopian L. donovani occurring in the same region. Patterns of inheritance are more complex than previously reported with multiple, apparently independent, origins from similar parents that include backcrossing with parental types. Analysis indicates that hybrids are representative of at least three different histories. Furthermore, isolates were highly polysomic at the level of chromosomes with differences between parasites recovered from a recrudescent infection from a previously treated individual. The results demonstrate that recombination is a significant feature of natural populations and contributes to the growing body of data that shows how recombination, and gene flow, shape natural populations of Leishmania.

Genetic structure and diversity analysis of tall fescue populations by EST-SSR and ISSR markers

Tall fescue is a perennial cool-season grass with economic importance especially in temperate regions of the northern hemisphere. This study was done to assess the genetic diversity and population structure of 90 tall fescue populations and cultivars using ISSR and EST-SSR markers in order to categorize valuable populations for breeding programs and to construct the core collection of tall fescue collection in Iran. The 10 EST-SSR primer pairs amplified 92 alleles. The allele numbers varied from 4 to 13 alleles per locus with an average of 9.2 alleles, of which 84 (90.6%) were polymorphic with an average of 8.4 polymorphic bands per primer. The 39 ISSR primers totally produced 387 scorable bands, of which 335 (86.6%) were polymorphic with an average of 8.6 polymorphic bands per primer. The amplified markers by ISSR primers varied from 6 to 14 markers per primer with an average of 9.92 markers per primer. The 90 tall fescue populations using both EST-SSR and ISSR data were classified into two clusters by UPGMA method that was coincide with PCA and structure analysis results. The turf-type and forage-type populations were clearly separated. Based on the results, the Iranian populations provide a valuable and novel germplasm to employ in tall fescue varietal improvement programs for both forage and turf-type applications. This progression is an important step to introduce this collection for development of a core collection of tall fescue germplasm in Iran.

Assessment of genetic diversity and population structure in cultured Australian Pacific oysters

The recent development of Pacific oyster (Crassostrea gigas) SNP genotyping arrays has allowed detailed characterisation of genetic diversity and population structure within and between oyster populations. It also raises the potential of harnessing genomic selection for genetic improvement in oyster breeding programmes. The aim of this study was to characterise a breeding population of Australian oysters through genotyping and analysis of 18 027 SNPs, followed by comparison with genotypes of oyster sampled from Europe and Asia. This revealed that the Australian populations had similar population diversity (H_E ) to oysters from New Zealand, the British Isles, France and Japan. Population divergence was assessed using PCA of genetic distance and revealed that Australian oysters were distinct from all other populations tested. Australian Pacific oysters originate from planned introductions sourced from three Japanese populations. Approximately 95% of these introductions were from geographically, and potentially genetically, distinct populations from the Nagasaki oysters assessed in this study. Finally, in preparation for the application of genomic selection in oyster breeding programmes, the strength of LD was evaluated and subsets of loci were tested for their ability to accurately infer relationships. Weak LD was observed on average; however, SNP subsets were shown to accurately reconstitute a genomic relationship matrix constructed using all loci. This suggests that low-density SNP panels may have utility in the Australian population tested, and the findings represent an important first step towards the design and implementation of genomic approaches for applied breeding in Pacific oysters.

Inference of Population Structure from Time-Series Genotype Data

Sequencing ancient DNA can offer direct probing of population history. Yet, such data are commonly analyzed with standard tools that assume DNA samples are all contemporary. We present DyStruct, a model and inference algorithm for inferring shared ancestry from temporally sampled genotype data. DyStruct explicitly incorporates temporal dynamics by modeling individuals as mixtures of unobserved populations whose allele frequencies drift over time. We develop an efficient inference algorithm for our model using stochastic variational inference. On simulated data, we show that DyStruct outperforms the current state of the art when individuals are sampled over time. Using a dataset of 296 modern and 80 ancient samples, we demonstrate DyStruct is able to capture a well-supported admixture event of steppe ancestry into modern Europe. We further apply DyStruct to a genome-wide dataset of 2,067 modern and 262 ancient samples used to study the origin of farming in the Near East. We show that DyStruct provides new insight into population history when compared with alternate approaches, within feasible run time.

Potential limits to the benefits of admixture during biological invasion

Species introductions often bring together genetically divergent source populations, resulting in genetic admixture. This geographic reshuffling of diversity has the potential to generate favourable new genetic combinations, facilitating the establishment and invasive spread of introduced populations. Observational support for the superior performance of admixed introductions has been mixed, however, and the broad importance of admixture to invasion questioned. Under most underlying mechanisms, admixture's benefits should be expected to increase with greater divergence among and lower genetic diversity within source populations, though these effects have not been quantified in invaders. We experimentally crossed source populations differing in divergence in the invasive plant Centaurea solstitialis. Crosses resulted in many positive (heterotic) interactions, but fitness benefits declined and were ultimately negative at high source divergence, with patterns suggesting cytonuclear epistasis. We explored the literature to assess whether such negative epistatic interactions might be impeding admixture at high source population divergence. Admixed introductions reported for plants came from sources with a wide range of genetic variation, but were disproportionately absent where there was high genetic divergence among native populations. We conclude that while admixture is common in species introductions and often happens under conditions expected to be beneficial to invaders, these conditions may be constrained by predictable negative genetic interactions, potentially explaining conflicting evidence for admixture's benefits to invasion.

Interspecific Gene Exchange as a Driver of Adaptive Evolution in Fungi

Throughout evolutionary history in the kingdom Fungi, taxa have exchanged genetic information among species, as revealed in particular by analyses of genome sequences. In fungi, hybridization can occur by sexual mating or by fusion of vegetative structures giving rise to new species or leaving traces of introgression in the genome. Furthermore, gene exchange can occur by horizontal gene transfer between species and can even include organisms outside the kingdom Fungi. In several cases, interspecific gene exchange has been instrumental in rapid adaptive evolution of fungal species and has notably played a role in the emergence of new pathogens. Here we summarize mechanisms and examples of gene exchange in fungi with a particular focus on the genomic context. We emphasize the need for and potential of applying population genetic approaches to better understand the processes and the impact of interspecific gene exchange in rapid adaptive evolution and species diversification. The broad occurrence of gene exchange among fungal species challenges our species concepts in the kingdom Fungi.

A tutorial on how not to over-interpret STRUCTURE and ADMIXTURE bar plots

Genetic clustering algorithms, implemented in programs such as STRUCTURE and ADMIXTURE, have been used extensively in the characterisation of individuals and populations based on genetic data. A successful example is the reconstruction of the genetic history of African Americans as a product of recent admixture between highly differentiated populations. Histories can also be reconstructed using the same procedure for groups that do not have admixture in their recent history, where recent genetic drift is strong or that deviate in other ways from the underlying inference model. Unfortunately, such histories can be misleading. We have implemented an approach, badMIXTURE, to assess the goodness of fit of the model using the ancestry "palettes" estimated by CHROMOPAINTER and apply it to both simulated data and real case studies. Combining these complementary analyses with additional methods that are designed to test specific hypotheses allows a richer and more robust analysis of recent demographic history.

Integrative approaches to guide conservation decisions: Using genomics to define conservation units and functional corridors

Climate change and increasing habitat loss greatly impact species survival, requiring range shifts, phenotypic plasticity and/or evolutionary change for long-term persistence, which may not readily occur unaided in threatened species. Therefore, defining conservation actions requires a detailed assessment of evolutionary factors. Existing genetic diversity needs to be thoroughly evaluated and spatially mapped to define conservation units (CUs) in an evolutionary context, and we address that here. We also propose a multidisciplinary approach to determine corridors and functional connectivity between CUs by including genetic diversity in the modelling while controlling for isolation by distance and phylogeographic history. We evaluate our approach on a Near Threatened Iberian endemic rodent by analysing genotyping-by-sequencing (GBS) genomic data from 107 Cabrera voles (Microtus cabrerae), screening the entire species distribution to define categories of CUs and their connectivity: We defined six management units (MUs) which can be grouped into four evolutionarily significant units (ESUs) and three (putatively) adaptive units (AUs). We demonstrate that the three different categories of CU can be objectively defined using genomic data, and their characteristics and connectivity can inform conservation decision-making. In particular, we show that connectivity of the Cabrera vole is very limited in eastern Iberia and that the pre-Pyrenean and part of the Betic geographic nuclei contribute the most to the species genetic diversity. We argue that a multidisciplinary framework for CU definition is essential and that this framework needs a strong evolutionary basis.

Fine population structure analysis method for genomes of many

Fine population structure can be examined through the clustering of individuals into subpopulations. The clustering of individuals in large sequence datasets into subpopulations makes the calculation of subpopulation specific allele frequency possible, which may shed light on selection of candidate variants for rare diseases. However, as the magnitude of the data increases, computational burden becomes a challenge in fine population structure analysis. To address this issue, we propose fine population structure analysis (FIPSA), which is an individual-based non-parametric method for dissecting fine population structure. FIPSA maximizes the likelihood ratio of the contingency table of the allele counts multiplied by the group. We demonstrated that its speed and accuracy were superior to existing non-parametric methods when the simulated sample size was up to 5,000 individuals. When applied to real data, the method showed high resolution on the Human Genome Diversity Project (HGDP) East Asian dataset. FIPSA was independently validated on 11,257 human genomes. The group assignment given by FIPSA was 99.1% similar to those assigned based on supervised learning. Thus, FIPSA provides high resolution and is compatible with a real dataset of more than ten thousand individuals.