On the genes, genealogies, and geographies of Quebec.

Population genetic models only provide coarse representations of real-world ancestry. We used a pedigree compiled from 4 million parish records and genotype data from 2276 French and 20,451 French Canadian individuals to finely model and trace French Canadian ancestry through space and time. The loss of ancestral French population structure and the appearance of spatial and regional structure highlights a wide range of population expansion models. Geographic features shaped migrations, and we find enrichments for migration, genetic, and genealogical relatedness patterns within river networks across regions of Quebec. Finally, we provide a freely accessible simulated whole-genome sequence dataset with spatiotemporal metadata for 1,426,749 individuals reflecting intricate French Canadian population structure. Such realistic population-scale simulations provide opportunities to investigate population genetics at an unprecedented resolution.

Spatially mapping the immune landscape of melanoma using imaging mass cytometry.

Melanoma is an immunogenic cancer with a high response rate to immune checkpoint inhibitors (ICIs). It harbors a high mutation burden compared with other cancers and, as a result, has abundant tumor-infiltrating lymphocytes (TILs) within its microenvironment. However, understanding the complex interplay between the stroma, tumor cells, and distinct TIL subsets remains a substantial challenge in immune oncology. To properly study this interplay, quantifying spatial relationships of multiple cell types within the tumor microenvironment is crucial. To address this, we used cytometry time-of-flight (CyTOF) imaging mass cytometry (IMC) to simultaneously quantify the expression of 35 protein markers, characterizing the microenvironment of 5 benign nevi and 67 melanomas. We profiled more than 220,000 individual cells to identify melanoma, lymphocyte subsets, macrophage/monocyte, and stromal cell populations, allowing for in-depth spatial quantification of the melanoma microenvironment. We found that within pretreatment melanomas, the abundance of proliferating antigen-experienced cytotoxic T cells (CD8+CD45RO+Ki67+) and the proximity of antigen-experienced cytotoxic T cells to melanoma cells were associated with positive response to ICIs. Our study highlights the potential of multiplexed single-cell technology to quantify spatial cell-cell interactions within the tumor microenvironment to understand immune therapy responses.

A review of UMAP in population genetics.

Uniform manifold approximation and projection (UMAP) has been rapidly adopted by the population genetics community to study population structure. It has become common in visualizing the ancestral composition of human genetic datasets, as well as searching for unique clusters of data, and for identifying geographic patterns. Here we give an overview of applications of UMAP in population genetics, provide recommendations for best practices, and offer insights on optimal uses for the technique.

Legacy Data Confound Genomics Studies.

Recent reports have identified differences in the mutational spectra across human populations. Although some of these reports have been replicated in other cohorts, most have been reported only in the 1000 Genomes Project (1kGP) data. While investigating an intriguing putative population stratification within the Japanese population, we identified a previously unreported batch effect leading to spurious mutation calls in the 1kGP data and to the apparent population stratification. Because the 1kGP data are used extensively, we find that the batch effects also lead to incorrect imputation by leading imputation servers and a small number of suspicious GWAS associations. Lower quality data from the early phases of the 1kGP thus continue to contaminate modern studies in hidden ways. It may be time to retire or upgrade such legacy sequencing data.

UMAP reveals cryptic population structure and phenotype heterogeneity in large genomic cohorts.

Human populations feature both discrete and continuous patterns of variation. Current analysis approaches struggle to jointly identify these patterns because of modelling assumptions, mathematical constraints, or numerical challenges. Here we apply uniform manifold approximation and projection (UMAP), a non-linear dimension reduction tool, to three well-studied genotype datasets and discover overlooked subpopulations within the American Hispanic population, fine-scale relationships between geography, genotypes, and phenotypes in the UK population, and cryptic structure in the Thousand Genomes Project data. This approach is well-suited to the influx of large and diverse data and opens new lines of inquiry in population-scale datasets.

Sulfur isotope fractionation during the evolutionary adaptation of a sulfate-reducing bacterium.

Dissimilatory sulfate reduction is a microbial catabolic pathway that preferentially processes less massive sulfur isotopes relative to their heavier counterparts. This sulfur isotope fractionation is recorded in ancient sedimentary rocks and generally is considered to reflect a phenotypic response to environmental variations rather than to evolutionary adaptation. Modern sulfate-reducing microorganisms isolated from similar environments can exhibit a wide range of sulfur isotope fractionations, suggesting that adaptive processes influence the sulfur isotope phenotype. To date, the relationship between evolutionary adaptation and isotopic phenotypes has not been explored. We addressed this by studying the covariation of fitness, sulfur isotope fractionation, and growth characteristics in Desulfovibrio vulgaris Hildenborough in a microbial evolution experiment. After 560 generations, the mean fitness of the evolved lineages relative to the starting isogenic population had increased by ∼ 17%. After 927 generations, the mean fitness relative to the initial ancestral population had increased by ∼ 20%. Growth rate in exponential phase increased during the course of the experiment, suggesting that this was a primary influence behind the fitness increases. Consistent changes were observed within different selection intervals between fractionation and fitness. Fitness changes were associated with changes in exponential growth rate but changes in fractionation were not. Instead, they appeared to be a response to changes in the parameters that govern growth rate: yield and cell-specific sulfate respiration rate. We hypothesize that cell-specific sulfate respiration rate, in particular, provides a bridge that allows physiological controls on fractionation to cross over to the adaptive realm.