Global diversity, recurrent evolution, and recent selection on amylase structural haplotypes in humans

The adoption of agriculture, first documented ~12,000 years ago in the Fertile Crescent, triggered a rapid shift toward starch-rich diets in human populations. Amylase genes facilitate starch digestion and increased salivary amylase copy number has been observed in some modern human populations with high starch intake, though evidence of recent selection is lacking. Here, using 52 long-read diploid assemblies and short read data from ~5,600 contemporary and ancient humans, we resolve the diversity, evolutionary history, and selective impact of structural variation at the amylase locus. We find that both salivary and pancreatic amylase genes have higher copy numbers in populations with agricultural subsistence compared to fishing, hunting, and pastoral groups. We identify 28 distinct amylase structural architectures and demonstrate that identical structures have arisen independently multiple times throughout recent human history. Using a pangenome graph-based approach to infer structural haplotypes across thousands of humans, we identify extensively duplicated haplotypes present at higher frequencies in modern agricultural populations. Leveraging 534 ancient human genomes we find that duplication-containing haplotypes have increased in frequency more than seven-fold over the last 12,000 years providing evidence for recent selection in Eurasians at this locus comparable in magnitude to that at lactase. Together, our study highlights the strong impact of the agricultural revolution on human genomes and the importance of long-read sequencing in identifying signatures of selection at structurally complex loci.

Methods and Insights from Single-Cell Expression Quantitative Trait Loci

Recent advancements in single-cell technologies have enabled expression quantitative trait locus (eQTL) analysis across many individuals at single-cell resolution. Compared with bulk RNA sequencing, which averages gene expression across cell types and cell states, single-cell assays capture the transcriptional states of individual cells, including fine-grained, transient, and difficult-to-isolate populations at unprecedented scale and resolution. Single-cell eQTL (sc-eQTL) mapping can identify context-dependent eQTLs that vary with cell states, including some that colocalize with disease variants identified in genome-wide association studies. By uncovering the precise contexts in which these eQTLs act, single-cell approaches can unveil previously hidden regulatory effects and pinpoint important cell states underlying molecular mechanisms of disease. Here, we present an overview of recently deployed experimental designs in sc-eQTL studies. In the process, we consider the influence of study design choices such as cohort, cell states, and ex vivo perturbations. We then discuss current methodologies, modeling approaches, and technical challenges as well as future opportunities and applications.

Assessing the recovery of Y chromosome microsatellites with population genomic data using Papio and Theropithecus genomes

Y chromosome markers can shed light on male-specific population dynamics but for many species no such markers have been discovered and are available yet, despite the potential for recovering Y-linked loci from available genome sequences. Here, we investigated how effective available bioinformatic tools are in recovering informative Y chromosome microsatellites from whole genome sequence data. In order to do so, we initially explored a large dataset of whole genome sequences comprising individuals at various coverages belonging to different species of baboons (genus: Papio) using Y chromosome references belonging to the same genus and more distantly related species (Macaca mulatta). We then further tested this approach by recovering Y-STRs from available Theropithecus gelada genomes using Papio and Macaca Y chromosome as reference sequences. Identified loci were validated in silico by a) comparing within-species relationships of Y chromosome lineages and b) genotyping male individuals in available pedigrees. Each STR was selected not to extend in its variable region beyond 100 base pairs, so that loci can be developed for PCR-based genotyping of non-invasive DNA samples. In addition to assembling a first set of Papio and Theropithecus Y-specific microsatellite markers, we released TYpeSTeR, an easy-to-use script to identify and genotype Y chromosome STRs using population genomic data which can be modulated according to available male reference genomes and genomic data, making it widely applicable across taxa.

The multifaceted genomic history of Ashaninka from Amazonian Peru

Despite its crucial location, the western side of Amazonia between the Andes and the source(s) of the Amazon River is still understudied from a genomic and archaeogenomic point of view, albeit possibly harboring essential information to clarify the complex genetic history of local Indigenous groups and their interactions with nearby regions,^{1,2,3,4,5,6,7,8} including central America and the Caribbean.^9,10,11,12 Focusing on this key region, we analyzed the genome-wide profiles of 51 Ashaninka individuals from Amazonian Peru, observing an unexpected extent of genomic variation. We identified at least two Ashaninka subgroups with distinctive genomic makeups, which were differentially shaped by the degree and timing of external admixtures, especially with the Indigenous groups from the Andes and the Pacific coast. On a continental scale, Ashaninka ancestors probably derived from a south-north migration of Indigenous groups moving into the Amazonian rainforest from a southeastern area with contributions from the Southern Cone and the Atlantic coast. These ancestral populations diversified in the variegated geographic regions of interior South America, on the eastern side of the Andes, differentially interacting with surrounding coastal groups. In this complex scenario, we also revealed strict connections between the ancestors of present-day Ashaninka, who belong to the Arawakan language family,¹³ and those Indigenous groups that moved further north into the Caribbean, contributing to the early Ceramic (Saladoid) tradition in the islands.^14,15.

Fasting renders immunotherapy effective against low-immunogenic breast cancer while reducing side effects

Immunotherapy is improving the prognosis and survival of cancer patients, but despite encouraging outcomes in different cancers, the majority of tumors are resistant to it, and the immunotherapy combinations are often accompanied by severe side effects. Here, we show that a periodic fasting-mimicking diet (FMD) can act on the tumor microenvironment and increase the efficacy of immunotherapy (anti-PD-L1 and anti-OX40) against the poorly immunogenic triple-negative breast tumors (TNBCs) by expanding early exhausted effector T cells, switching the cancer metabolism from glycolytic to respiratory, and reducing collagen deposition. Furthermore, FMD reduces the occurrence of immune-related adverse events (irAEs) by preventing the hyperactivation of the immune response. These results indicate that FMD cycles have the potential to enhance the efficacy of anti-cancer immune responses, expand the portion of tumors sensitive to immunotherapy, and reduce its side effects.

Abstract 1712: A single-cell RNA atlas of innate and adaptive intratumoral immunity in triple negative breast cancer during chemo- and immunotherapies

Checkpoint inhibitors (CIs) such as anti-PD-1 and anti PD-L1 have shown combinatorial clinical activity with chemotherapy in triple negative breast cancer (TNBC), but only in a minority of patients and only for a limited period of time. Moreover, it is currently unclear what chemotherapeutic drug is the most appropriate to combine with CIs in TNBC patients. We have investigated at the single cell level the transcriptome and the trajectories of more than 50,000 innate and adaptive intratumoral immune cells in two syngeneic, immune competent, orthotopic murine models of local and metastatic TNBC. Mice injected with 4T1 cells had a predominant lymphoid infiltrate, mice injected with EMT6 cells had a predominant myeloid infiltrate. Mice were treated with CIs and several different types of chemotherapeutics, alone or in combinations. In both models, capecitabine (alone or with CIs) was the less effective drug. Platinum, doxorubicin and taxanes showed synergy with CIs and had superimposable activity. Intermittent, medium dosage cyclophosphamide (CTX) plus vinorelbine and CIs was the most active combinatorial therapy (Falvo et al, Cancer Research 2021). Vinorelbine activated antigen presenting cells and CTX generated new T cell clones including stem cell-like TCF1+ CD8+ T cells. Treatments with most in vivo efficacy were associated to a decrease of regulatory T cells and of gamma delta T cells, which were found to have a pro-tumoral activity in these murine models, likely due to IL-17 expression in the neoplastic microenvironment. An increase of several different clusters of exhausted-like CD8+ T cells was observed in pre-clinical treatments with low efficacy; an opposite trend was found for several clusters of proliferative CD8+ T cells in treatments with high in vivo efficacy. Regarding macrophages, M2-like cells were enriched after treatments with low efficacy, while an opposite behaviour was found in M1-like macrophages. Interestingly, we observed a significant increase of an M1-like cluster with high expression of the Ly6c1/Ly6c2 gene in mice successfully treated with vinorelbine, CTX and CIs. For both cell lines the percentage of plasma B cells increased after in vivo treatments with high efficacy. In particular, the most effective treatment significantly increased the frequency of germinal B cells, which were absent in untreated tumors. These data can lead to new insights on the diagnosis and treatment of TNBC and to possible clinical applications.

Citation Format: Laura Carpen, Paolo Falvo, Stefania Orecchioni, Giulia Mitola, Roman Hillje, Saveria Mazzara, Patrizia Mancuso, Stefano Pileri, Alessandro Raveane, Francesco Bertolini. A single-cell RNA atlas of innate and adaptive intratumoral immunity in triple negative breast cancer during chemo- and immunotherapies [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 1712.

Genomic variation in baboons from central Mozambique unveils complex evolutionary relationships with other Papio species

Background

Gorongosa National Park in Mozambique hosts a large population of baboons, numbering over 200 troops. Gorongosa baboons have been tentatively identified as part of Papio ursinus on the basis of previous limited morphological analysis and a handful of mitochondrial DNA sequences. However, a recent morphological and morphometric analysis of Gorongosa baboons pinpointed the occurrence of several traits intermediate between P. ursinus and P. cynocephalus, leaving open the possibility of past and/or ongoing gene flow in the baboon population of Gorongosa National Park. In order to investigate the evolutionary history of baboons in Gorongosa, we generated high and low coverage whole genome sequence data of Gorongosa baboons and compared it to available Papio genomes.

Results

We confirmed that P. ursinus is the species closest to Gorongosa baboons. However, the Gorongosa baboon genomes share more derived alleles with P. cynocephalus than P. ursinus does, but no recent gene flow between P. ursinus and P. cynocephalus was detected when available Papio genomes were analyzed. Our results, based on the analysis of autosomal, mitochondrial and Y chromosome data, suggest complex, possibly male-biased, gene flow between Gorongosa baboons and P. cynocephalus, hinting to direct or indirect contributions from baboons belonging to the “northern” Papio clade, and signal the presence of population structure within P. ursinus.

Conclusions

The analysis of genome data generated from baboon samples collected in central Mozambique highlighted a complex set of evolutionary relationships with other baboons. Our results provided new insights in the population dynamics that have shaped baboon diversity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-01999-7.

Identification and Characterization of SOG1 (Suppressor of Gamma-Response 1) Homologues in Plants Using Data Mining Resources and Gene Expression Profiling

SOG1 (Suppressor of the Gamma response 1) is the master-regulator of plant DNA damage response (DDR), a highly coordinated network of DNA damage sensors, transducers, mediators, and effectors, with highly coordinated activities. SOG1 transcription factor belongs to the NAC/NAM protein family, containing the well-conserved NAC domain and five serine-glutamine (SQ) motifs, preferential targets for phosphorylation by ATM and ATR. So far, the information gathered for the SOG1 function comes from studies on the model plant Arabidopsis thaliana. To expand the knowledge on plant-specific DDR, it is opportune to gather information on other SOG1 orthologues. The current study identified plants where multiple SOG1 homologues are present and evaluated their functions by leveraging the information contained in publicly available transcriptomics databases. This analysis revealed the presence of multiple SOG1 sequences in thirteen plant species, and four (Medicago truncatula, Glycine max, Kalankoe fedtschenkoi, Populus trichocarpa) were selected for gene expression data mining based on database availability. Additionally, M. truncatula seeds and seedlings exposed to treatments known to activate DDR pathways were used to evaluate the expression profiles of MtSOG1a and MtSOG1b. The experimental workflow confirmed the data retrieved from transcriptomics datasets, suggesting that the SOG1 homologues have redundant functions in different plant species.

A single-cell transcriptomic landscape of innate and adaptive intratumoral immunity in triple negative breast cancer during chemo- and immunotherapies

Breast cancer (BC) constitutes a major health problem worldwide, making it the most common malignancy in women. Current treatment options for BC depend primarily on histological type, molecular markers, clinical aggressiveness and stage of disease. Immunotherapy, such as αPD-1, have shown combinatorial clinical activity with chemotherapy in triple negative breast cancer (TNBC) delineating some therapeutic combinations as more effective than others. However, a clear overview of the main immune cell populations involved in these treatments has never been provided.

Here, an assessment of the immune landscape in the tumor microenvironment (TME) of two TNBC mouse models has been performed using single-cell RNA sequencing technology. Specifically, immune cells were evaluated in untreated conditions and after treatments with chemotherapy or immunotherapy used as single agents or in combination. A decrease of Treg was found in treatments with in vivo efficacy as well as γδ T cells, which have a pro-tumoral activity in mice. Focusing on Cd8 T cells, across all the conditions, a general increase of exhausted-like Cd8 T cells was confirmed in pre-clinical treatments with low efficacy and an opposite trend was found for the proliferative Cd8 T cells. Regarding macrophages, M2-like cells were enriched in treatments with low efficacy while M1-like macrophages followed an opposite trend. For both models, similar proportions of B cells were detected with an increase of proliferative B cells in treatments involving cisplatin in combination with αPD-1. The fine-scale characterization of the immune TME in this work can lead to new insights on the diagnosis and treatment of TNBC.

The genetic origin of Daunians and the Pan-Mediterranean southern Italian Iron Age context

The geographical location and shape of Apulia, a narrow land stretching out in the sea at the South of Italy, made this region a Mediterranean crossroads connecting Western Europe and the Balkans. Such movements culminated at the beginning of the Iron Age with the Iapygian civilization which consisted of three cultures: Peucetians, Messapians, and Daunians. Among them, the Daunians left a peculiar cultural heritage, with one-of-a-kind stelae and pottery, but, despite the extensive archaeological literature, their origin has been lost to time. In order to shed light on this and to provide a genetic picture of Iron Age Southern Italy, we collected and sequenced human remains from three archaeological sites geographically located in Northern Apulia (the area historically inhabited by Daunians) and radiocarbon dated between 1157 and 275 calBCE. We find that Iron Age Apulian samples are still distant from the genetic variability of modern-day Apulians, they show a degree of genetic heterogeneity comparable with the cosmopolitan Republican and Imperial Roman civilization, even though a few kilometers and centuries separate them, and they are well inserted into the Iron Age Pan-Mediterranean genetic landscape. Our study provides for the first time a window on the genetic make-up of pre-Roman Apulia, whose increasing connectivity within the Mediterranean landscape, would have contributed to laying the foundation for modern genetic variability. In this light, the genetic profile of Daunians may be compatible with an at least partial autochthonous origin, with plausible contributions from the Balkan peninsula.

Weaving Mitochondrial DNA and Y-Chromosome Variation in the Panamanian Genetic Canvas

The Isthmus of Panama was a crossroads between North and South America during the continent’s first peopling (and subsequent movements) also playing a pivotal role during European colonization and the African slave trade. Previous analyses of uniparental systems revealed significant sex biases in the genetic history of Panamanians, as testified by the high proportions of Indigenous and sub-Saharan mitochondrial DNAs (mtDNAs) and by the prevalence of Western European/northern African Y chromosomes. Those studies were conducted on the general population without considering any self-reported ethnic affiliations. Here, we compared the mtDNA and Y-chromosome lineages of a new sample collection from 431 individuals (301 males and 130 females) belonging to either the general population, mixed groups, or one of five Indigenous groups currently living in Panama. We found different proportions of paternal and maternal lineages in the Indigenous groups testifying to pre-contact demographic events and genetic inputs (some dated to Pleistocene times) that created genetic structure. Then, while the local mitochondrial gene pool was marginally involved in post-contact admixtures, the Indigenous Y chromosomes were differentially replaced, mostly by lineages of western Eurasian origin. Finally, our new estimates of the sub-Saharan contribution, on a more accurately defined general population, reduce an apparent divergence between genetic and historical data.

The Dual Role of Innate Lymphoid and Natural Killer Cells in Cancer. from Phenotype to Single-Cell Transcriptomics, Functions and Clinical Uses

Simple Summary

Innate lymphoid cells (ILCs), a family of innate immune cells including natural killers (NKs), play a multitude of roles in first-line cancer control, in escape from immunity and in cancer progression. In this review, we summarize preclinical and clinical data on ILCs and NK cells concerning their phenotype, function and clinical applications in cellular therapy trials. We also describe how single-cell transcriptome sequencing has been used and forecast how it will be used to better understand ILC and NK involvement in cancer control and progression as well as their therapeutic potential.

Abstract

The role of innate lymphoid cells (ILCs), including natural killer (NK) cells, is pivotal in inflammatory modulation and cancer. Natural killer cell activity and count have been demonstrated to be regulated by the expression of activating and inhibitory receptors together with and as a consequence of different stimuli. The great majority of NK cell populations have an anti-tumor activity due to their cytotoxicity, and for this reason have been used for cellular therapies in cancer patients. On the other hand, the recently classified helper ILCs are fundamentally involved in inflammation and they can be either helpful or harmful in cancer development and progression. Tissue niche seems to play an important role in modulating ILC function and conversion, as observed at the transcriptional level. In the past, these cell populations have been classified by the presence of specific cellular receptor markers; more recently, due to the advent of single-cell RNA sequencing (scRNA-seq), it has been possible to also explore them at the transcriptomic level. In this article we review studies on ILC (and NK cell) classification, function and their involvement in cancer. We also summarize the potential application of NK cells in cancer therapy and give an overview of the most recent studies involving ILCs and NKs at scRNA-seq, focusing on cancer. Finally, we provide a resource for those who wish to start single-cell transcriptomic analysis on the context of these innate lymphoid cell populations.

Projecting Ancient Ancestry in Modern-Day Arabians and Iranians: A Key Role of the Past Exposed Arabo-Persian Gulf on Human Migrations

The Arabian Peninsula is strategic for investigations centered on the early structuring of modern humans in the wake of the out-of-Africa migration. Despite its poor climatic conditions for the recovery of ancient human DNA evidence, the availability of both genomic data from neighboring ancient specimens and informative statistical tools allow modeling the ancestry of local modern populations. We applied this approach to a data set of 741,000 variants screened in 291 Arabians and 78 Iranians, and obtained insightful evidence. The west-east axis was a strong forcer of population structure in the Peninsula, and, more importantly, there were clear continuums throughout time linking western Arabia with the Levant, and eastern Arabia with Iran and the Caucasus. Eastern Arabians also displayed the highest levels of the basal Eurasian lineage of all tested modern-day populations, a signal that was maintained even after correcting for a possible bias due to a recent sub-Saharan African input in their genomes. Not surprisingly, eastern Arabians were also the ones with highest similarity with Iberomaurusians, who were, so far, the best proxy for the basal Eurasians amongst the known ancient specimens. The basal Eurasian lineage is the signature of ancient non-Africans who diverged from the common European-eastern Asian pool before 50,000 years ago, prior to the later interbred with Neanderthals. Our results appear to indicate that the exposed basin of the Arabo-Persian Gulf was the possible home of basal Eurasians, a scenario to be further investigated by searching ancient Arabian human specimens.

Abstract 1653: A single-cell atlas of the effect of chemotherapeutics over intratumoral immune cells reveals that combining an alkylating agent and a vinca alkaloid can activate antigen presenting cells and increase tcf1+ stem-like CD8+ T-cells, thus improving anti-PD-1 efficacy in triple negative breast cancer and lymphoma

Checkpoint inhibitors (CIs) instigate anticancer immunity, albeit only in a fraction of patients. Triple negative breast cancer (TNBC) is among the most aggressive and lethal BC types, and currently available therapies have an unsatisfactory clinical impact. The association of CIs with chemotherapy had encouraging results in randomized TNBC trials, but so far there has been no clear evidence of what should be considered the best combinatorial therapy. In B-cell non-Hodgkin's lympoma (NHL), CIs' clinical activity has been so far unsatisfactory. We designed the current study to define what could be considered the most efficient combinatorial regimen of chemotherapy plus CI in TNBC and NHL. We investigated different options including platinum (P), doxorubicin (D), taxol (T), vinorelbine (V), and cyclophosphamide ( C), i.e. the most effective chemotherapy drugs for these diseases. Capecitabine was not considered as found not to be effective in combination with CIs. As TNBC patients have a large intratumoral immune cell heterogeneity, we investigated two immune competent TNBC models: 4T1 with a predominant lymphoid intratumoral infiltrate and EMT6 with a predominant myeloid infiltrate. A C57BL/6 Eμ-myc transgenic NHL model, with a predominant lymphatic dissemination pattern, was also studied. We found that 1) Intermittent, medium-dosage C (C140), was more effective than other combinatorial regimens including the CI anti-PD1 plus low-dose C, or regular-dose T, D, or P. The association of V and anti-PD1 further increased the preclinical efficacy of C140; 2) both CD3+CD4+ and CD3+CD8+ T cells (at variance with other subsets of NK and myeloid cells) were needed to control local and metastatic neoplastic growth in mice treated with C140, V and anti-PD1. 3) The combinatorial therapy including C140, V and anti-PD1 was the most efficient in activating antigen presenting cells (APCs), and the investigation of TCR repertoire indicated that this therapy generated a significantly larger clonal expansion of both CD3+CD4+ and CD3+CD8+ T cells. 4) Single-cell transcriptome analysis of >50,000 intratumoural immune cells showed after C140, V and anti-PD1 therapy a gene signature suggestive of a change resulting from exposure to a mitogen, ligand, or an antigen for which it is specific, as well as APC-to-T-cell adhesion. This transcriptional program also increased intratumoural tcf1+ stem-like CD8+ T-cells and altered the balance between terminally and progenitor exhausted T-cells, favoring the latter. These data support the clinical investigation of this therapy in TNBC and NHL, and we are currently investigating the related transcriptional mechanisms inducing APC activation and T cell rewiring.

Citation Format: Paolo Falvo, Stefania Orecchioni, Roman Hillje, Alessandro Raveane, Giulia Mitola, Patrizia Mancuso, Chiara Camisaschi, Lucilla Luzi, PierGiuseppe Pelicci, Francesco Bertolini. A single-cell atlas of the effect of chemotherapeutics over intratumoral immune cells reveals that combining an alkylating agent and a vinca alkaloid can activate antigen presenting cells and increase tcf1+ stem-like CD8+ T-cells, thus improving anti-PD-1 efficacy in triple negative breast cancer and lymphoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 1653.

Drug Repurposing in Oncology, an Attractive Opportunity for Novel Combinatorial Regimens

The costs of developing, validating and buying new drugs are dramatically increasing. On the other hand, sobering economies have difficulties in sustaining their healthcare systems, particularly in countries with an elderly population requiring increasing welfare. This conundrum requires immediate action, and a possible option is to study the large, already present arsenal of drugs approved and to use them for innovative therapies. This possibility is particularly interesting in oncology, where the complexity of the cancer genome dictates in most patients a multistep therapeutic approach. In this review, we discuss a) Computational approaches; b) preclinical models; c) currently ongoing or already published clinical trials in the drug repurposing field in oncology; and d) drug repurposing to overcome resistance to previous therapies.

Archaeogenomic distinctiveness of the Isthmo-Colombian area

The recently enriched genomic history of Indigenous groups in the Americas is still meager concerning continental Central America. Here, we report ten pre-Hispanic (plus two early colonial) genomes and 84 genome-wide profiles from seven groups presently living in Panama. Our analyses reveal that pre-Hispanic demographic events contributed to the extensive genetic structure currently seen in the area, which is also characterized by a distinctive Isthmo-Colombian Indigenous component. This component drives these populations on a specific variability axis and derives from the local admixture of different ancestries of northern North American origin(s). Two of these ancestries were differentially associated to Pleistocene Indigenous groups that also moved into South America, leaving heterogenous genetic footprints. An additional Pleistocene ancestry was brought by a still unsampled population of the Isthmus (UPopI) that remained restricted to the Isthmian area, expanded locally during the early Holocene, and left genomic traces up to the present day.

Plant TDP1 (Tyrosyl-DNA Phosphodiesterase 1): A Phylogenetic Perspective and Gene Expression Data Mining

The TDP1 (tyrosyl-DNA phosphodiesterase 1) enzyme removes the non-specific covalent intermediates between topoisomerase I and DNA, thus playing a crucial role in preventing DNA damage. While mammals possess only one TDP1 gene, in plants two genes (TDP1α and TDP1β) are present constituting a small gene subfamily. These display a different domain structure and appear to perform non-overlapping functions in the maintenance of genome integrity. Namely, the HIRAN domain identified in TDP1β is involved in the interaction with DNA during the recognition of stalled replication forks. The availability of transcriptomic databases in a growing variety of experimental systems provides new opportunities to fill the current gaps of knowledge concerning the evolutionary origin and the specialized roles of TDP1 genes in plants. Whereas a phylogenetic approach has been used to track the evolution of plant TDP1 protein, transcriptomic data from a selection of representative lycophyte, eudicots, and monocots have been implemented to explore the transcriptomic dynamics in different tissues and a variety of biotic and abiotic stress conditions. While the phylogenetic analysis indicates that TDP1α is of non-plant origin and TDP1β is plant-specific originating in ancient vascular plants, the gene expression data mining comparative analysis pinpoints for tissue- and stress-specific responses.

Cyclophosphamide and Vinorelbine Activate Stem-Like CD8+ T Cells and Improve Anti-PD-1 Efficacy in Triple-Negative Breast Cancer

Checkpoint inhibitors (CI) instigate anticancer immunity in many neoplastic diseases, albeit only in a fraction of patients. The clinical success of cyclophosphamide (C)-based haploidentical stem-cell transplants indicates that this drug may re-orchestrate the immune system. Using models of triple-negative breast cancer (TNBC) with different intratumoral immune contexture, we demonstrate that a combinatorial therapy of intermittent C, CI, and vinorelbine activates antigen-presenting cells (APC), and abrogates local and metastatic tumor growth by a T-cell-related effect. Single-cell transcriptome analysis of >50,000 intratumoral immune cells after therapy treatment showed a gene signature suggestive of a change resulting from exposure to a mitogen, ligand, or antigen for which it is specific, as well as APC-to-T-cell adhesion. This transcriptional program also increased intratumoral Tcf1⁺ stem-like CD8⁺ T cells and altered the balance between terminally and progenitor-exhausted T cells favoring the latter. Overall, our data support the clinical investigation of this therapy in TNBC. SIGNIFICANCE: A combinatorial therapy in mouse models of breast cancer increases checkpoint inhibition by activating antigen-presenting cells, enhancing intratumoral Tcf1⁺ stem-like CD8⁺ T cells, and increasing progenitor exhausted CD8⁺ T cells.

Circulating endothelial progenitors are increased in COVID-19 patients and correlate with SARS-CoV-2 RNA in severe cases

BACKGROUND

During the course of COVID-19, the disease caused by the new coronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), thrombotic phenomena and/or diffuse vascular damage are frequent, and viral elements have been observed within endothelial cells.

OBJECTIVES

CD146 + circulating endothelial cells (CD146 + CECs) and their progenitors (CEPs) are increased in cardiovascular, thrombotic, infectious, and cancer diseases. The present study was designed to investigate their kinetics in novel coronavirus (COVID-19) patients.

METHODS

We used a validated flow cytometry procedure to enumerate viable and apoptotic CD146 + CECs and CEPs in COVID-19 patients during the course of the disease and in patients who recovered.

RESULTS

Viable CEPs per milliliter were significantly increased in COVID-19 patients compared with healthy controls. This increase was observed in patients with mild symptoms and not further augmented in patients with severe symptoms. In patients who recovered, CEPs decreased, but were in a range still significantly higher than normal controls. Regarding mature CD146 + CECs, in COVID-19 patients, their absolute number was similar to those observed in healthy controls, but the viable/apoptotic CD146 + CEC ratio was significantly different. Both mild and severe COVID-19 patients had significantly less apoptotic CD146 + CECs compared with healthy controls. Patients who recovered had significantly less CD146 + CECs per milliliter when compared with controls as well as to mild and severe COVID-19 patients. A positive correlation was found between the copies of SARS-CoV-2 RNA in the cellular fraction and apoptotic CEPs per milliliter in severe COVID-19 patients.

CONCLUSIONS

CD146 + CECs and CEPs might be investigated as candidate biomarkers of endothelial damage in COVID-19 patients.

Abstract 900: Intermittent cyclophosphamide and vinorelbine reshape the immune cell environment, induce T cell clonal replacement and increase the efficacy of PD-1 inhibition in models of triple negative breast cancer

The clinical success of cyclophosphamide (Cy)-based haploidentical stem cell transplants in hematological malignancies indicates that this drug has the potential to reorchestrate the immune system against cancer cells. Along a similar way, Cy is administered before the infusion of CAR-T cells to improve their clinical efficacy. We have previously found that daily, low/medium-dose Cy (20mg/Kg), in association with vinorelbine (V) was able to improve the preclinical efficacy of checkpoint inhibitors (CIs) anti-PD-1 and anti-PD-L1 in preclinical models of breast cancer and lymphoma (Orecchioni et al, 2018). Randomized clinical trials have recently indicated that the addition of CIs to other chemotherapy drugs such as taxanes (T), doxorubicin (D), or platinum (P) might be beneficial in triple negative breast cancer (TNBC) patients. Thus, we designed the present study in two orthotopic, immunocompetent, local and metastatic models of TNBC (CI-resistant 4T1 and EMT6, moderately responsive to CIs) to investigate whether Cy and V, in different schedules and doses, were more effective than T, D, or P. In vivo studies in local and metastatic models using both 4T1 and EMT6 TNBC cells clearly indicated that intermittent (ie every 6 days), medium-dosage (140 mg/Kg) Cy, was more effective than any other combinatorial regimens including CIs plus daily low/medium-dose Cy, or T, D, or P at their regular dosages. The association of V further increased the preclinical efficacy of intermittent Cy plus CIs, and abrogated tumor growth in both 4T1 and EMT6 models. In vivo studies with neutralizing monoclonal antibodies targeting in separate experiments T, B, NK, and myeloid cells demonstrated that CD3+CD4+, CD3+CD8+ T cells as well as CD11b+ monocyte/macrophage dendritic cells were crucial to abrogate tumor growth in TNBC-bearing mice treated with intermittent Cy, V and CIs. Single-cell transcriptome analysis of more than 60,000 intratumoral immune cells and flow cytometry studies in circulating and intratumoral subsets of immune cells indicated that V promoted the generation and maturation of myeloid APC cells, and that intermittent Cy generated new clones of tumor-infiltrating CD3+CD4+ and CD3+CD8+ TCR alpha beta cells. After treatment with intermittent Cy, V and CIs, intratumoral immune cells showed a unique gene signature of 42 genes (CD3D, Ptprc, CD69, Lat, Lck, CD2, CD28, CD27, and CD3E among others) suggestive of a change in morphology and behavior resulting from exposure to a mitogen, cytokine, chemokine, cellular ligand, or an antigen for which it is specific, as well as APC-to-T cell adhesion. Taken together, our data support the hypothesis that APC priming and T cell clonal replacement can significantly improve CI efficacy in vivo in TNBC models. We are further investigating the Cy-related mechanisms inducing intratumoral T-cell clonal replacement and how these T cells cooperate with V-induced APCs.

Citation Format: Paolo Falvo, Stefania Orecchioni, Roman Hillje, Alessandro Raveane, Patrizia Mancuso, Chiara Camisaschi, Lucilla Luzi, Francesco Bertolini. Intermittent cyclophosphamide and vinorelbine reshape the immune cell environment, induce T cell clonal replacement and increase the efficacy of PD-1 inhibition in models of triple negative breast cancer [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 900.

The mitogenome portrait of Umbria in Central Italy as depicted by contemporary inhabitants and pre-Roman remains

Umbria is located in Central Italy and took the name from its ancient inhabitants, the Umbri, whose origins are still debated. Here, we investigated the mitochondrial DNA (mtDNA) variation of 545 present-day Umbrians (with 198 entire mitogenomes) and 28 pre-Roman individuals (obtaining 19 ancient mtDNAs) excavated from the necropolis of Plestia. We found a rather homogeneous distribution of western Eurasian lineages across the region, with few notable exceptions. Contemporary inhabitants of the eastern part, delimited by the Tiber River and the Apennine Mountains, manifest a peculiar mitochondrial proximity to central-eastern Europeans, mainly due to haplogroups U4 and U5a, and an overrepresentation of J (30%) similar to the pre-Roman remains, also excavated in East Umbria. Local genetic continuities are further attested to by six terminal branches (H1e1, J1c3, J2b1, U2e2a, U8b1b1 and K1a4a) shared between ancient and modern mitogenomes. Eventually, we identified multiple inputs from various population sources that likely shaped the mitochondrial gene pool of ancient Umbri over time, since early Neolithic, including gene flows with central-eastern Europe. This diachronic mtDNA portrait of Umbria fits well with the genome-wide population structure identified on the entire peninsula and with historical sources that list the Umbri among the most ancient Italic populations.

Y-chromosome and Surname Analyses for Reconstructing Past Population Structures: The Sardinian Population as a Test Case

Many anthropological, linguistic, genetic and genomic analyses have been carried out to evaluate the potential impact that evolutionary forces had in shaping the present-day Sardinian gene pool, the main outlier in the genetic landscape of Europe. However, due to the homogenizing effect of internal movements, which have intensified over the past fifty years, only partial information has been obtained about the main demographic events. To overcome this limitation, we analyzed the male-specific region of the Y chromosome in three population samples obtained by reallocating a large number of Sardinian subjects to the place of origin of their monophyletic surnames, which are paternally transmitted through generations in most of the populations, much like the Y chromosome. Three Y-chromosome founding lineages, G2-L91, I2-M26 and R1b-V88, were identified as strongly contributing to the definition of the outlying position of Sardinians in the European genetic context and marking a significant differentiation within the island. The present distribution of these lineages does not always mirror that detected in ancient DNAs. Our results show that the analysis of the Y-chromosome gene pool coupled with a sampling method based on the origin of the family name, is an efficient approach to unravelling past heterogeneity, often hidden by recent movements, in the gene pool of modern populations. Furthermore, the reconstruction and comparison of past genetic isolates represent a starting point to better assess the genetic information deriving from the increasing number of available ancient DNA samples.

The Genomic Impact of European Colonization of the Americas

The human genetic diversity of the Americas has been affected by several events of gene flow that have continued since the colonial era and the Atlantic slave trade. Moreover, multiple waves of migration followed by local admixture occurred in the last two centuries, the impact of which has been largely unexplored. Here, we compiled a genome-wide dataset of ∼12,000 individuals from twelve American countries and ∼6,000 individuals from worldwide populations and applied haplotype-based methods to investigate how historical movements from outside the New World affected (1) the genetic structure, (2) the admixture profile, (3) the demographic history, and (4) sex-biased gene-flow dynamics of the Americas. We revealed a high degree of complexity underlying the genetic contribution of European and African populations in North and South America, from both geographic and temporal perspectives, identifying previously unreported sources related to Italy, the Middle East, and to specific regions of Africa.

Genome-wide analysis of Corsican population reveals a close affinity with Northern and Central Italy

Despite being the fourth largest island in the Mediterranean basin, the genetic variation of Corsica has not been explored as exhaustively as Sardinia, which is situated only 11 km South. However, it is likely that the populations of the two islands shared, at least in part, similar demographic histories. Moreover, the relative small size of the Corsica may have caused genetic isolation, which, in turn, might be relevant under medical and translational perspectives. Here we analysed genome wide data of 16 Corsicans, and integrated with newly (33 individuals) and previously generated samples from West Eurasia and North Africa. Allele frequency, haplotype-based, and ancient genome analyses suggest that although Sardinia and Corsica may have witnessed similar isolation and migration events, the latter is genetically closer to populations from continental Europe, such as Northern and Central Italians.

Population structure of modern-day Italians reveals patterns of ancient and archaic ancestries in Southern Europe

European populations display low genetic differentiation as the result of long-term blending of their ancient founding ancestries. However, it is unclear how the combination of ancient ancestries related to early foragers, Neolithic farmers, and Bronze Age nomadic pastoralists can explain the distribution of genetic variation across Europe. Populations in natural crossroads like the Italian peninsula are expected to recapitulate the continental diversity, but have been systematically understudied. Here, we characterize the ancestry profiles of Italian populations using a genome-wide dataset representative of modern and ancient samples from across Italy, Europe, and the rest of the world. Italian genomes capture several ancient signatures, including a non-steppe contribution derived ultimately from the Caucasus. Differences in ancestry composition, as the result of migration and admixture, have generated in Italy the largest degree of population structure detected so far in the continent, as well as shaping the amount of Neanderthal DNA in modern-day populations.

Analysis of the human Y-chromosome haplogroup Q characterizes ancient population movements in Eurasia and the Americas

BACKGROUND

Recent genome studies of modern and ancient samples have proposed that Native Americans derive from a subset of the Eurasian gene pool carried to America by an ancestral Beringian population, from which two well-differentiated components originated and subsequently mixed in different proportion during their spread in the Americas. To assess the timing, places of origin and extent of admixture between these components, we performed an analysis of the Y-chromosome haplogroup Q, which is the only Pan-American haplogroup and accounts for virtually all Native American Y chromosomes in Mesoamerica and South America.

RESULTS

Our analyses of 1.5 Mb of 152 Y chromosomes, 34 re-sequenced in this work, support a "coastal and inland routes scenario" for the first entrance of modern humans in North America. We show a major phase of male population growth in the Americas after 15 thousand years ago (kya), followed by a period of constant population size from 8 to 3 kya, after which a secondary sign of growth was registered. The estimated dates of the first expansion in Mesoamerica and the Isthmo-Colombian Area, mainly revealed by haplogroup Q-Z780, suggest an entrance in South America prior to 15 kya. During the global constant population size phase, local South American hints of growth were registered by different Q-M848 sub-clades. These expansion events, which started during the Holocene with the improvement of climatic conditions, can be ascribed to multiple cultural changes rather than a steady population growth and a single cohesive culture diffusion as it occurred in Europe.

CONCLUSIONS

We established and dated a detailed haplogroup Q phylogeny that provides new insights into the geographic distribution of its Eurasian and American branches in modern and ancient samples.

Reconstructing the genetic history of Italians: new insights from a male (Y-chromosome) perspective

BACKGROUND

Due to its central and strategic position in Europe and in the Mediterranean Basin, the Italian Peninsula played a pivotal role in the first peopling of the European continent and has been a crossroad of peoples and cultures since then.

AIM

This study aims to gain more information on the genetic structure of modern Italian populations and to shed light on the migration/expansion events that led to their formation.

SUBJECTS AND METHODS

High resolution Y-chromosome variation analysis in 817 unrelated males from 10 informative areas of Italy was performed. Haplogroup frequencies and microsatellite haplotypes were used, together with available data from the literature, to evaluate Mediterranean and European inputs and date their arrivals.

RESULTS

Fifty-three distinct Y-chromosome lineages were identified. Their distribution is in general agreement with geography, southern populations being more differentiated than northern ones.

CONCLUSIONS

A complex genetic structure reflecting the multifaceted peopling pattern of the Peninsula emerged: southern populations show high similarity with those from the Middle East and Southern Balkans, while those from Northern Italy are close to populations of North-Western Europe and the Northern Balkans. Interestingly, the population of Volterra, an ancient town of Etruscan origin in Tuscany, displays a unique Y-chromosomal genetic structure.

Phylogeography of human Y-chromosome haplogroup Q3-L275 from an academic/citizen science collaboration

Background

The Y-chromosome haplogroup Q has three major branches: Q1, Q2, and Q3. Q1 is found in both Asia and the Americas where it accounts for about 90% of indigenous Native American Y-chromosomes; Q2 is found in North and Central Asia; but little is known about the third branch, Q3, also named Q1b-L275. Here, we combined the efforts of population geneticists and genetic genealogists to use the potential of full Y-chromosome sequencing for reconstructing haplogroup Q3 phylogeography and suggest possible linkages to events in population history.

Results

We analyzed 47 fully sequenced Y-chromosomes and reconstructed the haplogroup Q3 phylogenetic tree in detail. Haplogroup Q3-L275, derived from the oldest known split within Eurasian/American haplogroup Q, most likely occurred in West or Central Asia in the Upper Paleolithic period. During the Mesolithic and Neolithic epochs, Q3 remained a minor component of the West Asian Y-chromosome pool and gave rise to five branches (Q3a to Q3e), which spread across West, Central and parts of South Asia. Around 3–4 millennia ago (Bronze Age), the Q3a branch underwent a rapid expansion, splitting into seven branches, some of which entered Europe. One of these branches, Q3a1, was acquired by a population ancestral to Ashkenazi Jews and grew within this population during the 1st millennium AD, reaching up to 5% in present day Ashkenazi.

Conclusions

This study dataset was generated by a massive Y-chromosome genotyping effort in the genetic genealogy community, and phylogeographic patterns were revealed by a collaboration of population geneticists and genetic genealogists. This positive experience of collaboration between academic and citizen science provides a model for further joint projects. Merging data and skills of academic and citizen science promises to combine, respectively, quality and quantity, generalization and specialization, and achieve a well-balanced and careful interpretation of the paternal-side history of human populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0870-2) contains supplementary material, which is available to authorized users.

Exploring the Y Chromosomal Ancestry of Modern Panamanians

Geologically, Panama belongs to the Central American land-bridge between North and South America crossed by Homo sapiens >14 ka ago. Archaeologically, it belongs to a wider Isthmo-Colombian Area. Today, seven indigenous ethnic groups account for 12.3% of Panama’s population. Five speak Chibchan languages and are characterized by low genetic diversity and a high level of differentiation. In addition, no evidence of differential structuring between maternally and paternally inherited genes has been reported in isthmian Chibchan cultural groups. Recent data have shown that 83% of the Panamanian general population harbour mitochondrial DNAs (mtDNAs) of Native American ancestry. Considering differential male/female mortality at European contact and multiple degrees of geographical and genetic isolation over the subsequent five centuries, the Y-chromosome Native American component is expected to vary across different geographic regions and communities in Panama. To address this issue, we investigated Y-chromosome variation in 408 modern males from the nine provinces of Panama and one indigenous territory (the comarca of Kuna Yala). In contrast to mtDNA data, the Y-chromosome Native American component (haplogroup Q) exceeds 50% only in three populations facing the Caribbean Sea: the comarca of Kuna Yala and Bocas del Toro province where Chibchan languages are spoken by the majority, and the province of Colón where many Kuna and people of mixed indigenous-African-and-European descent live. Elsewhere the Old World component is dominant and mostly represented by western Eurasian haplogroups, which signal the strong male genetic impact of invaders. Sub-Saharan African input accounts for 5.9% of male haplotypes. This reflects the consequences of the colonial Atlantic slave trade and more recent influxes of West Indians of African heritage. Overall, our findings reveal a local evolution of the male Native American ancestral gene pool, and a strong but geographically differentiated unidirectional sex bias in the formation of local modern Panamanian populations.