A population-level statistic for assessing Mendelian behavior of genotyping-by-sequencing data from highly duplicated genomes

Genome-wide analysis tracks the emergence of intraspecific polyploids in Phragmites australis

Polyploidization plays an important role in plant speciation and adaptation. To address the role of polyploidization in grass diversification, we studied Phragmites australis, an invasive species with intraspecific variation in chromosome numbers ranging from 2n = 36 to 144. We utilized a combined analysis of ploidy estimation, phylogeny, population genetics and model simulations to investigate the evolution of P. australis. Using restriction site-associated DNA sequencing (RAD-seq), we conducted a genome-wide analysis of 88 individuals sourced from diverse populations worldwide, revealing the presence of six distinct intraspecific lineages with extensive genetic admixture. Each lineage was characterized by a specific ploidy level, predominantly tetraploid or octoploid, indicative of multiple independent polyploidization events. The population size of each lineage has declined moderately in history while remaining large, except for the North American native and the US Land types, which experienced constant population size contraction throughout their history. Our investigation did not identify direct association between polyploidization events and grass invasions. Nonetheless, we observed octoploid and hexaploid lineages at contact zones in Romania, Hungary, and South Africa, suggestively due to genomic conflicts arising from allotetraploid parental lineages.

A Public Mid-Density Genotyping Platform for Hexaploid Sweetpotato (Ipomoea batatas [L.] Lam)

Small public breeding programs focused on specialty crops have many barriers to adopting technology, particularly creating and using genetic marker panels for genomic-based decisions in selection. Here, we report the creation of a DArTag panel of 3120 loci distributed across the sweetpotato (Ipomoea batatas [L.] Lam) genome for molecular-marker-assisted breeding and genomic prediction. The creation of this marker panel has the potential to bring cost-effective and rapid genotyping capabilities to sweetpotato breeding programs worldwide. The open access provided by this platform will allow the genetic datasets generated on the marker panel to be compared and joined across projects, institutions, and countries. This genotyping resource has the power to make routine genotyping a reality for any breeder of sweetpotato.

Variant calling in polyploids for population and quantitative genetics

Advancements in genome assembly and sequencing technology have made whole genome sequence (WGS) data and reference genomes accessible to study polyploid species. Compared to popular reduced-representation sequencing approaches, the genome-wide coverage and greater marker density provided by WGS data can greatly improve our understanding of polyploid species and polyploid biology. However, biological features that make polyploid species interesting also pose challenges in read mapping, variant identification, and genotype estimation. Accounting for characteristics in variant calling like allelic dosage uncertainty, homology between subgenomes, and variance in chromosome inheritance mode can reduce errors. Here, I discuss the challenges of variant calling in polyploid WGS data and discuss where potential solutions can be integrated into a standard variant calling pipeline.

Chronosequence of invasion reveals minimal losses of population genomic diversity, niche expansion, and trait divergence in the polyploid, leafy spurge

Rapid evolution may play an important role in the range expansion of invasive species and modify forecasts of invasion, which are the backbone of land management strategies. However, losses of genetic variation associated with colonization bottlenecks may constrain trait and niche divergence at leading range edges, thereby impacting management decisions that anticipate future range expansion. The spatial and temporal scales over which adaptation contributes to invasion dynamics remain unresolved. We leveraged detailed records of the ~130-year invasion history of the invasive polyploid plant, leafy spurge (Euphorbia virgata), across ~500 km in Minnesota, U.S.A. We examined the consequences of range expansion for population genomic diversity, niche breadth, and the evolution of germination behavior. Using genotyping-by-sequencing, we found some population structure in the range core, where introduction occurred, but panmixia among all other populations. Range expansion was accompanied by only modest losses in sequence diversity, with small, isolated populations at the leading edge harboring similar levels of diversity to those in the range core. The climatic niche expanded during most of the range expansion, and the niche of the range core was largely non-overlapping with the invasion front. Ecological niche models indicated that mean temperature of the warmest quarter was the strongest determinant of habitat suitability and that populations at the leading edge had the lowest habitat suitability. Guided by these findings, we tested for rapid evolution in germination behavior over the time course of range expansion using a common garden experiment and temperature manipulations. Germination behavior diverged from the early to late phases of the invasion, with populations from later phases having higher dormancy at lower temperatures. Our results suggest that trait evolution may have contributed to niche expansion during invasion and that distribution models, which inform future management planning, may underestimate invasion potential without accounting for evolution.

The Conservation Genetics of Iris lacustris (Dwarf Lake Iris), a Great Lakes Endemic

Iris lacustris, a northern Great Lakes endemic, is a rare species known from 165 occurrences across Lakes Michigan and Huron in the United States and Canada. Due to multiple factors, including habitat loss, lack of seed dispersal, patterns of reproduction, and forest succession, the species is threatened. Early population genetic studies using isozymes and allozymes recovered no to limited genetic variation within the species. To better explore genetic variation across the geographic range of I. lacustris and to identify units for conservation, we used tunable Genotyping-by-Sequencing (tGBS) with 171 individuals across 24 populations from Michigan and Wisconsin, and because the species is polyploid, we filtered the single nucleotide polymorphism (SNP) matrices using polyRAD to recognize diploid and tetraploid loci. Based on multiple population genetic approaches, we resolved three to four population clusters that are geographically structured across the range of the species. The species migrated from west to east across its geographic range, and minimal genetic exchange has occurred among populations. Four units for conservation are recognized, but nine adaptive units were identified, providing evidence for local adaptation across the geographic range of the species. Population genetic analyses with all, diploid, and tetraploid loci recovered similar results, which suggests that methods may be robust to variation in ploidy level.

Computation of the distribution of model accuracy statistics in machine learning: Comparison between analytically derived distributions and simulation-based methods

Background and Aims

All fields have seen an increase in machine-learning techniques. To accurately evaluate the efficacy of novel modeling methods, it is necessary to conduct a critical evaluation of the utilized model metrics, such as sensitivity, specificity, and area under the receiver operator characteristic curve (AUROC). For commonly used model metrics, we proposed the use of analytically derived distributions (ADDs) and compared it with simulation-based approaches.

Methods

A retrospective cohort study was conducted using the England National Health Services Heart Disease Prediction Cohort. Four machine learning models (XGBoost, Random Forest, Artificial Neural Network, and Adaptive Boost) were used. The distribution of the model metrics and covariate gain statistics were empirically derived using boot-strap simulation (N = 10,000). The ADDs were created from analytic formulas from the covariates to describe the distribution of the model metrics and compared with those of bootstrap simulation.

Results

XGBoost had the most optimal model having the highest AUROC and the highest aggregate score considering six other model metrics. Based on the Anderson-Darling test, the distribution of the model metrics created from bootstrap did not significantly deviate from a normal distribution. The variance created from the ADD led to smaller SDs than those derived from bootstrap simulation, whereas the rest of the distribution remained not statistically significantly different.

Conclusions

ADD allows for cross study comparison of model metrics, which is usually done with bootstrapping that rely on simulations, which cannot be replicated by the reader.

Extensive polyploid clonality was a successful strategy for seagrass to expand into a newly submerged environment

Polyploidy has the potential to allow organisms to outcompete their diploid progenitor(s) and occupy new environments. Shark Bay, Western Australia, is a World Heritage Area dominated by temperate seagrass meadows including Poseidon's ribbon weed, Posidonia australis. This seagrass is at the northern extent of its natural geographic range and experiences extremes in temperature and salinity. Our genomic and cytogenetic assessments of 10 meadows identified geographically restricted, diploid clones (2n = 20) in a single location, and a single widespread, high-heterozygosity, polyploid clone (2n = 40) in all other locations. The polyploid clone spanned at least 180 km, making it the largest known example of a clone in any environment on earth. Whole-genome duplication through polyploidy, combined with clonality, may have provided the mechanism for P. australis to expand into new habitats and adapt to new environments that became increasingly stressful for its diploid progenitor(s). The new polyploid clone probably formed in shallow waters after the inundation of Shark Bay less than 8500 years ago and subsequently expanded via vegetative growth into newly submerged habitats.