Next Article in Journal
Insight into the Genetic Population Structure of Wild Red Foxes in Poland Reveals Low Risk of Genetic Introgression from Escaped Farm Red Foxes
Next Article in Special Issue
How a Paleogenomic Approach Can Provide Details on Bioarchaeological Reconstruction: A Case Study from the Globular Amphorae Culture
Previous Article in Journal
Investigating Cellular Trajectories in the Severity of COVID-19 and Their Transcriptional Programs Using Machine Learning Approaches
Previous Article in Special Issue
Maternal Lineages from 10–11th Century Commoner Cemeteries of the Carpathian Basin
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

First Bronze Age Human Mitogenomes from Calabria (Grotta Della Monaca, Southern Italy)

1
Department of Cultural Heritage, University of Bologna, via Degli Ariani 1, 48121 Ravenna, Italy
2
Department of Neuroscience and Rehabilitation, University of Ferrara, Corso Ercole I D’Este 32, 44121 Ferrara, Italy
3
Centro Regionale di Speleologia “Enzo dei Medici”, via Lucania 3, 87070 Roseto Capo Spulico, Italy
4
Department of Biological Geological and Environmental Sciences, University of Bologna, via Selmi 3, 40126 Bologna, Italy
5
Department of Biology, Università degli Studi di Firenze, via del Proconsolo 12, 50125 Firenze, Italy
6
Interdepartmental Centre “Alma Mater Research Institute on Global Challenges and Climate Change (Alma Climate)”, University of Bologna, via Petroni 26, 40126 Bologna, Italy
7
Gruppo di Ricerca Speleo-Archeologica, Università degli Studi di Bari “Aldo Moro”, Piazza Umberto I 1, 70121 Bari, Italy
*
Authors to whom correspondence should be addressed.
These authors equally contributed to the article.
Genes 2021, 12(5), 636; https://doi.org/10.3390/genes12050636
Submission received: 26 February 2021 / Revised: 2 April 2021 / Accepted: 20 April 2021 / Published: 25 April 2021
(This article belongs to the Special Issue Ancient and Archaic Genomes)

Abstract

:
The Italian peninsula was host to a strong history of migration processes that shaped its genomic variability since prehistoric times. During the Metal Age, Sicily and Southern Italy were the protagonists of intense trade networks and settlements along the Mediterranean. Nonetheless, ancient DNA studies in Southern Italy are, at present, still limited to prehistoric and Roman Apulia. Here, we present the first mitogenomes from a Middle Bronze Age cave burial in Calabria to address this knowledge gap. We adopted a hybridization capture approach, which enabled the recovery of one complete and one partial mitochondrial genome. Phylogenetic analysis assigned these two individuals to the H1e and H5 subhaplogroups, respectively. This preliminary phylogenetic analysis supports affinities with coeval Sicilian populations, along with Linearbandkeramik and Bell Beaker cultures maternal lineages from Central Europe and Iberia. Our work represents a starting point which contributes to the comprehension of migrations and population dynamics in Southern Italy, and highlights this knowledge gap yet to be filled by genomic studies.

1. Introduction

The genomic variability within the Italian Peninsula is greater today than in any European country, which may suggest that this area played a pivotal role in the peopling of the Mediterranean in the past [1]. Genetic studies based on autosomal and uniparental markers [2,3,4], genome-wide [1,5,6,7,8,9] and whole genome approaches [10] have dissected the clinal variability of the present Italian population, revealing multilayered patterns of prehistorical and historical processes of migration and admixture that occurred throughout the Peninsula. In particular, these studies provided evidence of an early divergence of Italian groups dating back to the Late Glacial period, with further differentiation attributed to Neolithic and Bronze Age migrations [10].
Modern Southern Italian populations have been extensively studied [4,6,7,10,11], highlighting a long series of migration processes and cultural exchanges that influenced this area. High frequencies of maternal lineages from the Caucasus and the Levant were retrieved, which predate the Neolithic and may support the role of this area as a refugium during the Last Glacial Maximum (LGM) [6,12]. A wide genetic Mediterranean “continuum” was identified, which links Southern Italy with Crete and the Caucasus following Neolithic and Bronze Age migrations from the Near East. In addition, during the Bronze Age, a non-steppe contribution derived from the Caucasus was detected [7,10].
The Metal Age, in particular, deeply characterized the complex population processes of protohistoric Italy, leaving clear signatures in modern populations [13]. Territorialization spread in both Tyrrhenian and Adriatic populations during the Bronze Age, and dynamism was further promoted following novel social structuring [14]. The emergence of elite groups during the second millennium BC [15] resulted in community organization pivoted around kinship and inherited rank in Southern communities [16]. The Cosenza province (central Calabria) played an important role in the landscape of prehistoric studies of Southern Italy due to this structuring. Evidence of consistent mining activities nearby the Sila Plateau allowed protohistoric Calabria to establish itself as one of the fundamental territories for metallurgy in Italy, and to connect this area in a strong exchange relationship within the lower Tyrrhenian coastal communities [17]. The presence of a Protoapennine culture in many Middle Bronze Age (MBA) Calabrian contexts, such as Broglio di Trebisacce [18], was suggested by Ardesia [19] to be the “Calabrian” Rodì-Tindari-Vallelunga (RTV) horizon, thus underlining the strong interactions and shared cultural patterns between the southern areas and the western islands of the peninsula. Nevertheless, a deep analysis of funerary contexts from protohistoric Southern Italy is essential. The presence of numerous cave sites that have not yet been investigated in the area of Tyrrhenian Calabria prompt the need for a deeper knowledge of the ancient communities of Southern Italy. The northwestern sector of Calabria is rich in cave sites, with over two hundred caves present along the Tyrrhenian coast and inland, such as Grotte di Cirella in Diamante, Grotte di Torre Talao in Scalea, Grotta della Madonna in Praia a Mare, Grotta del Romito in Papasidero, Grotte di Sant’Angelo in Cassano allo Ionio and Grotta di Donna Marsilia [20]. Therefore, paleogenetic studies from this area, in particular, are warranted to investigate early population dynamics. To complement genetic studies on modern populations, ancient DNA (aDNA) studies have the potential to provide precise insights into early heritage and migrations, coupled with sociocultural aspects, of past societies. However, aDNA studies in Southern Italy are limited to prehistoric [21,22] and Roman [23] Apulia. Hence, we aimed here to produce mitogenome data on Bronze Age communities from Calabria, in Southern Italy, to provide new information on population dynamics from this understudied area and to contextualize them with archaeological and anthropological evidence. To achieve this, we explored the remains of a MBA cave burial in Tyrrhenian Calabria, Grotta della Monaca (Figure 1).

2. Materials and Methods

2.1. Archaeological Context and Human Remains

Grotta della Monaca is a karstic cave located in Sant’Agata di Esaro at 600 m above sea level, close to the Esaro river (Cosenza, northwestern Calabria, Southern Italy) (Figure 1). The cave has a long human history [20,24,25] due to mining activities that began during the Upper Paleolithic period and intensified during the Late Neolithic/Early Chalcolithic period [26]. The cave partially changed its role during the MBA, when the deepest areas were used as burial grounds [27,28]. Archaeological investigations conducted between 2003 and 2010 recovered a mass grave with numerous human skeletal remains in a niche of the area of Cunicoli terminali (hereafter m5v). Collective burials in a natural cave were a common and uninterrupted ritual in protohistoric Southern Italy, particularly in the coastal sites facing the Tyrrhenian Sea [29]. From the Early Bronze Age (EBA), this funerary practice was reserved for high-status individuals [30]. Fragmented pottery and funerary artifacts evidence a link between the burials in Grotta della Monaca and the so-called Apennine culture [31], which was an archaeological complex that was widespread in the Italian Peninsula and Sicily during the second half of the second millennium BC. The remains found in the m5v context represent at least 24 different commingled individuals (File S1 and Table S1), and were in a poor state of preservation due to several taphonomic factors [32]. The osteological remains buried in m5v of Grotta della Monaca were examined by traditional anthropological methodologies at the Laboratory of Archaeo-Anthropology and Forensic Anthropology of the University of Ferrara (See File S1 for details about the methods used and the anthropological data acquired).
Radiocarbon analysis of the human bone samples was conducted at CEDAD (Center of Applied Physics, Dating and Diagnostics, University of Salento, Lecce, Italy). Six out of seven samples were successfully dated and were assigned to the MBA period. Calibrated and uncalibrated dates are summarized in Table S2.

2.2. Experimental Procedure

Ancient DNA analysis was performed on seven individuals (Table S3 and Figure S1) from the m5v mass grave. Due to the poor preservation of the remains, the only petrous bone available (preferred because it is DNA-rich [33,34]) was associated with individual 24. Thus, four shafts of long bones (GdM5, GdM6, GdM7, GdM12) and two cranial fragments (GdM11, GdM14) were sampled along with the petrous bone (GdM24). The samples were selected based on the preservation state, and took into account the bone elements used to calculate the minimum number of individuals (MNI), to make sure we sampled different individuals. Analysis was conducted in the ancient DNA Laboratory (aDNA Lab) of the Department of Cultural Heritage, University of Bologna (Ravenna, Italy). Post-PCR procedures were performed in a separate laboratory, and negative control was carried out alongside each experiment to confirm the absence of intra-laboratory contamination. All materials and instruments used were DNA-free and decontaminated before use with bleach and/or DNA ExitusPlus™ (AppliChem GmbH, Darmstadt, Germany) and researchers followed high-sterility protocols according to ancient DNA authenticity criteria [35,36,37,38].
Specimens were first cleaned with sterile scalpels to remove sediments and calcareous concretions, and 1–2 mm of the outer layer of bones was abraded with circular drilling discs. Bones were decontaminated for 1 h under UV irradiation, and 140 to 340 mg fragments of the bone specimens were pulverized, while 255 mg of bone powder was collected from petrous bone by directly drilling the densest region of the cochlea.

2.3. Ancient DNA Extraction

DNA from the samples was extracted following a two-day silica spin-column protocol [39,40] with slight modifications, as in Cilli et al. [41]. Samples were digested in a rotator at 37 °C overnight in 3000 µL solution composed of 2700 µL EDTA, 37.5 µL Proteinase K and H2O. Supernatant was centrifuged and transferred with PB binding buffer (Qiagen, Hilden, Germany) into silica columns (Roche-High Pure Viral Nucleic Acid Large Volume Kit, Roche, Basel, Switzerland). After two washing steps with PE buffer (Qiagen), DNA was eluted in 50 µL of EB buffer (Qiagen).

2.4. Library Preparation and Enrichment

Double-stranded library preparation was performed according to Carøe et al. [42]. Libraries were indexed with Illumina sequencing adapters and no enzymatic damage repair was carried out in order to preserve patterns of ancient DNA fragments. A qPCR quantification was performed prior to indexing, in order to assign correct amplification cycles to each library and to reduce the formation of heteroduplex structures. A total volume of 30 µL from each library was split in 2/3 aliquots and used for indexing PCR. PCR indexing products were purified with MinElute PCR Purification kit (Qiagen), eluted in EBT buffer (10 mM Tris-Cl, pH 8.5 and 0.05% Tween-20) and analyzed on Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA, USA) to produce qualitative and quantitative estimations of the libraries.
Indexed libraries were equimolarly pooled to a total of maximum 2 µg and enriched for human mitochondrial DNA (mtDNA) for the bead-capture method using long-range PCR products, according to Maricic et al. [43]. The captured libraries were quantified using Agilent 2100 Bioanalyzer and sequenced on Illumina MiSeq platform at the Advanced Genomic Centre of the University of Florence. The sequencing was run as paired-end with 75 × 2 + 8 + 8 cycles.

2.5. Bioinformatic Analysis

Raw reads were first visualized with FastQC [44], and adapter sequences were filtered using AdapterRemoval v2.3.0 [45]. Fragments shorter than 30 bp and with a ‘minquality < 25’ (PHRED values) were discarded, and paired-end reads with overlaps of at least 10 bp were collapsed into single sequences. Filtered reads were mapped against the revised Cambridge Reference Sequence (rCRS, NC_012920.1) [46] with BWA v.0.7.17-r188 [47]. Command ‘aln’ was set with -n number to 0.01, seed-length option deactivated and -o value at 3 for tolerating a higher number of gaps in the alignment. The generated .sai file was aligned against the reference fasta file with ‘samse’ command, and then converted to .bam files and sorted by leftmost coordinate using SAMtools v1.10 [48] ‘view’ and ‘sort’ commands. The Picard tool MarkDuplicates [49] was used on the .bam files to mark and remove PCR duplicates, and reads were locally realigned through Genome Analysis Toolkit tools’ [50] ‘RealignerTargetCreator’ and ‘IndelRealigner’ commands, to minimize the number of mismatches around the indels that may be easily mistaken as SNPs.
Data authenticity and post-mortem damage were investigated through MapDamage v2.0.8 [51] and a new rescaled .bam file was generated with the ‘rescale’ option, that downscales the quality scores at positions likely affected by damage patterns. Damage patterns at the 5′ and 3′ ends of the reads were also computed with contDeam, a tool provided within the Schmutzi package [52], and present-day human contamination estimates were performed with Schmutzi, using a non-redundant database of 256 mitogenomes available in the software package. Results are summarized in Table S3. Mitochondrial genomes were called using Schmutzi, and variants were then investigated using GATK ‘HaplotypeCaller’ and filtered with ‘VariantFiltration’ command, storing only detected variants with quality ‘QUAL20’ (≥20) and depth of coverage ‘COV3’ (≥3). Diagnostic variants were checked and verified with snpToolkit v2.2.3 [53], and only variants with ratio −r 0.9 were used to determine the haplogroup.

2.6. Phylogenetic Analyses

Phylogenetic analyses were performed on the complete mitogenome of individual GdM24, while GdM7 was only analyzed for haplogroup assignment. The other samples were not included due to inconsistencies in the data (Table S3). The mtDNA haplogroup assignment was predicted using both Haplofind [54] and HaploGrep [55] software, and checked manually according to PhyloTree build 17 [56]. The consensus sequence obtained for GdM24 was then compared to a reference dataset of 97 ancient individuals extracted from the literature (Table S4). In addition, 13 published mitogenomes from ancient African and European samples that supply basal lineage information for phylogenetic reconstruction were included in the analyses (Table S5). The RSRS and rCRS reference sequences were also added for comparison purposes in the Network analysis.
Sequence alignment was performed with the DNA Alignment software [57] and was visually confirmed. The poly-C and AC-indels at 303–315, 515–524, 573–576 and 16,180–16,193 positions, as well as the nucleotide position 16,519, were excluded from phylogenetic analyses [58]. A Median Joining Network was calculated by using Network v.10.2 [59] with no pre- or post-processing steps. Furthermore, a phylogenetic tree, based on the same multi-alignment dataset after retaining unique sequences, was reconstructed with BEAST v1.8.0 [60] by using reference ancient dated samples as calibration points (Table S5). The best substitution model for the dataset was tested with Mega 5.2 [61], resulting in the Hasegawa–Kishino–Yano model with fixed fraction of invariable sites and gamma distributed rates (HKY + I + G). We ran BEAST analysis with 200,000,000 MCMC generations, and sampled every 2000 iterations by testing different combinations of clock models (i.e., Strict clock vs. Uncorrelated Relaxed Lognormal-ULN clock) and tree models (i.e., Constant Population Size vs. Bayesian Skyline). Chain convergence was assessed with Tracer v 1.6, resulting in Effective Sampling Size (ESS) values higher than 200 for all the parameters and in all the tested model combinations. Evaluating the maximum likelihood estimates (MLE) for the four combinations of clock and tree models further revealed the best support for the Relaxed ULN clock model and Bayesian Skyline tree model; therefore, this combination was chosen for the final step of analysis. The Maximum Clade Credibility tree was calculated using TreeAnnotator v1.8.0 by discarding the first 20% of iterations as burn-in. The resulting phylogenetic tree was graphically represented with FigTree [62].

3. Results

Six human bone fragments and one petrous bone from the m5v burial area of Grotta della Monaca were analyzed for ancient mtDNA on the Illumina MiSeq platform. Five out of seven samples contained low concentrations of DNA (Table S3) and were not further analyzed. For the sample GdM24, we obtained 10,471 mapping reads with quality ≥30, which resulted in a mean coverage of 34X, and 100% of positions covered at least 3-fold. GdM7 contained 550 reads mapped with quality ≥30, resulting in a mean coverage of 2X, and 91% of positions covered at least 1-fold. Analysis of the read datasets with MapDamage v2.0.8 revealed miscoding lesion distribution patterns typical of ancient DNA in both samples (Figure S2). The presence of positions’ specific substitutions increases at the ends of the reads, with an average frequency of 28.67% (GdM7) and 45.46% (GdM24) of C-to-T substitution at the first base at 5′, and similar average frequency of 33.80% (GdM7) and 42.85% (GdM24) G-to-A substitution at 3′ ends. The average fragment length is around 56.49 for GdM24 and 77.33 for GdM7. These data are consistent with the antiquity of the samples and the environmental conditions of the site.
The GdM24 mitochondrial genome was assigned with 85.63% accuracy (Haplogrep) and 1.0 score (Haplofind) to the H1e mitochondrial subclade, based on the relevance of the eleven diagnostic variants detected (263G, 1438G, 3010A, 4769G, 5460A, 6779G, 8860G, 15326G, 16172C, 16224C, 16519C). Regarding GdM7, haplogroup assignment was inferred on the partial sequence, evaluating the presence of diagnostic variants (263G, 1438G, 2626C, 4769G, 8860G, 16304C, 16310A). This appears to belong to the mitochondrial subclade H5, with 79.95% accuracy on Haplogrep and 1.0 quality score on Haplofind, respectively. However, due to the low coverage and incompleteness of the mitogenome, GdM7 was not included in the subsequent phylogenetic analysis.
Phylogenetic reconstruction performed with both Median Joining Network (MJN) and Bayesian Evolutionary Analysis of Sampling Trees (BEAST) confirmed the attribution of GdM24 to the mtDNA haplogroup H1e. Accordingly, the sample indeed clusters with the other H1e samples extracted from the literature, within a clade that splits up after H1 and the other basal-considered lineages (Figure 2 and Figure 3). A star-like pattern characterized both H1 and H1e haplotypes in the Network analysis (Figure 2). In particular, the GdM24 sample branches, along with two Sicilian Bronze Age samples (I3876 and I7774 [63]) from a median vector, originated from the H1e basal node which contains samples from Middle/Late Neolithic and Bronze Age cultures from Central Europe, and a Bell Beaker sample from Southern Italy (I4930 [63]). The same clustering pattern was consistently confirmed in the BEAST reconstructed phylogenetic tree (Figure 3), with the branch leading to the GdM24 sequence, particularly, dating at 4125 years BP (95% HPD: 2301–6332).

4. Discussion

Excluding Sicily and Sardinia that experienced different processes of migration and admixture and were covered by several paleogenomic studies [63,64,65,66,67,68,69], Bronze Age Southern Italy has not been studied with an aDNA approach. In this study we attempted analysis of the mitochondrial DNA variation of seven MBA samples from Grotta della Monaca cave, located in Calabria, Southern Italy, with the aim of filling the gap in the contribution of ancient genetic data to the population dynamics of this peculiar area. This is of particular interest and importance considering the key role that the Italian Peninsula played in past demographic processes, as evidenced by the wide variability in genomic diversity in present-day Italian populations [1].
Despite the poor preservation of the remains, the hybridization capture method allowed us to retrieve and reconstruct the complete mitochondrial genome of individual GdM24, and partially reconstruct the mitochondrial genome of individual GdM7. These are the first attested ancient molecular data from the Calabria region, as highlighted by a recent review of the state-of-the-art European paleogenomic data [70]. Unsurprisingly, GdM24 was the only sample for which the petrous bone was available from this burial area. Previous studies suggest that this skeletal element preserves the endogenous DNA better than other bones or tissues [33,71]. Thus, this type of sample is confirmed as fundamental in the recovery of ancient endogenous DNA from particularly difficult contexts, such as those at low latitudes.
The two samples for which we obtained complete mitogenomes are both assigned to mitochondrial macrohaplogroup H. Nowadays, this is the most widespread in Europe, encompassing over 40% of the total mtDNA variation, and could be considered as the quintessential Eurasiatic marker. Previous studies have shown how haplogroup H arrived in Europe from the Near East before LGM (∼22,000 BP), survived in south western glacial refugia and then colonized Central and Northern Europe [72,73]. In particular, GdM24 was assigned to mitochondrial sublineage H1e, and GdM7 to H5.
Among the subclades, H1 is the most frequent in modern-day Europeans, followed by H3 and H5 [74]. When compared to other present-day European and Middle Eastern populations, haplogroup H1 shows frequency peaks in the Franco-Cantabrian region, declining from the west towards Eastern and Southern Europe [75]. Taking into account both the distribution and the coalescence age of this subhaplogroup (9.3–12.8 kya; [76]), it has been hypothesized that it could be correlated with a Late Glacial re-expansion of populations from the Franco-Iberian refugial areas at the end of the Ice Age, from ∼15 kya, or at the end of the Younger Dryas, ∼11.5 kya [77,78,79,80,81]. Along with H3, H1 subclade was also proposed as a marker for a Bell Beaker complex expansion originating in Iberia [72]; however, a recent genome-wide study on ancient genomes does not support migration as an important mechanism of the spread of this culture between Iberia and Central Europe [66]. Ancient DNA evidence also highlighted the earliest presence of the H1e subhaplogroup in Neolithic samples from Hungary (∼7000 BP) [82] and Germany (∼6100 BP) [72]. In Italy, the H1e subhaplogroup was retrieved in Sicily from two Middle Neolithic samples (∼6800 BP), an Eneolithic (∼4700 BP) and an EBA sample (∼4000 BP) [63]. In the present study, phylogenetic analyses of GdM24 highlighted the proximity of this sample, among others, with two Sicilian Bronze Age samples (I3876 and I7774 [63]) recovered in Marcita and Contrada Paolina necropolises, respectively.
Regarding H5, this subhaplogroup shows a coalescence time of ∼10.7–17.1 kya [76]. There is no agreement about the place of origin of H5; hypotheses have been made about eastern Mediterranean [74,80] or south-west European origin, where it has expanded after the Ice Age, with several possible dispersal routes [76]. However, the dynamics that led to the diffusion of this haplogroup are not yet clearly understood. The oldest evidence of H5 lineage was retrieved in Anatolia (∼8200 BP) [83] and Bulgaria (∼7600 BP) [84]. A migration with a center of origin in the eastern Mediterranean, that would have carried this lineage into Italy and, to a lesser extent, the western Mediterranean, until northern and western Iberia, has been hypothesized [80].
The mitochondrial haplogroups found in present-day Calabria ([7], personal data) include all of the haplogroups that have been found in the ancient individuals analyzed here (H1e and H5). None of these are unique to the ancient Southern Italian populations and, therefore, cannot be taken as indisputable proof of local population continuity. Nonetheless, they deserve to be more deeply analyzed at the nuclear level to obtain more information about the ancestry and migration patterns of this area.
In the neighboring Sicily, past population dynamics are better understood due to sampling conducted in two projects [63,66] that detected ancestry typical of early European farmers as a mixture of Anatolia Neolithic and Western Hunter-Gatherer, and no evidence of steppe ancestry in Middle Neolithic samples. They found evidence of steppe ancestry in the EBA, by around 2200 BC, with the forming of a clade with EBA Mallorca, suggesting that the population may harbor ancestry most plausibly from Iberia, with a scenario of west-to-east gene flow. Fernandes et al. [63] retrieved Iranian-related ancestry in Sicily, from the MBA (1800–1500 BC) and Late Bronze Age (LBA), which was widespread among the Aegean by the MBA, along with the Mycenaean cultural expansion or earlier.
The active role of Tyrrhenian Calabria in cultural connections between eastern Sicily, the Aeolian Islands and the Adriatic environments in protohistoric times have been highlighted [85,86,87,88,89]. These connections are demonstrated by the chrono-typological study of pottery, and reflect strong relations between central-southern Calabria, the Aeolian Islands and the northern part of Sicily [90]. However, the so-called Palma Campania facies that identifies the MBA contexts of southwestern Italy still generates uncertainties, especially in Tyrrhenian Calabria, due to the lack of in-depth knowledge of the archaeological contexts, and the absence of ancient genomic data. In light of this evidence, our data represents a baseline for a deeper comprehension of the population pattern of prehistoric and protohistoric Southern Italy.
In the area of Tyrrhenian Calabria, the presence of over two hundred cave sites with attested past human frequentation, not yet investigated, prompts the need for a deeper knowledge of the ancient communities of Southern Italy. Future studies, focused on genome-wide data from this peculiar area, could clarify migratory and demographic processes that took place in the prehistory and protohistory of the Southern European continent, provide information about the ancestry of individuals through time and allow the study of intrinsic differences in migratory behavior related to sex-biased processes, adaptation and admixture.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/genes12050636/s1, File S1: Archaeological context and anthropological analyses, Figure S1: specimens used for molecular analyses, Figure S2: MapDamage plot, Table S1: Summary of biological profile and pathologies in individuals from Grotta della Monaca involved in genetic analyses, Table S2: Radiocarbon dating and calibration executed on human samples from Grotta della Monaca, Table S3: Summary of the genomic analyses, Table S4: Ancient samples belonging to H1 and H1e haplogroups used in phylogenetic analyses; Table S5: Dataset of reference ancient sequences used as tip calibration points in the BEAST analyses.

Author Contributions

Conceptualization, D.L., E.G.-R. and F.L.; methodology E.C., F.F., and A.M.; software F.F., S.S., and A.J.A.; validation E.C., and F.F.; formal analysis E.C., F.F., F.A., A.L., A.J.A., and S.D.F.; investigation E.C., F.F., and F.A.; resources D.L.; data curation F.F., and S.S.; writing—original drafting preparation, F.F., E.C. and F.A.; writing—review and editing F.F., E.C., F.A., S.D.F., S.S., A.M., A.J.A., A.L., P.A., D.C., M.L., F.L., E.G.-R., and D.L.; visualization S.S. and F.F.; supervision D.L.; project administration D.L.; funding acquisition D.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the MIUR-PRIN 20177PJ9XF grant to D.L. and D.C. A.J.A. is supported by SeaChanges—European Union’s Horizon 2020 programme, under the Marie Skłodowska-Curie grant agreement No 813383.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Analyzed data for GdM24 sample are available on NCBI GenBank (Accession number: MW853672). Raw data for all samples have been uploaded on ENA database (Study accession number: PRJEB44125).

Acknowledgments

We would like to thank Thomas Mignani, Andrea De Giovanni, Irene Boreggio and Isabella Cavaliere who assisted during the experimental procedures.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sazzini, M.; Ruscone, G.A.G.; Giuliani, C.; Sarno, S.; Quagliariello, A.; De Fanti, S.; Boattini, A.; Gentilini, D.; Fiorito, G.; Catanoso, M.; et al. Complex interplay between neutral and adaptive evolution shaped differential genomic background and disease susceptibility along the Italian peninsula. Sci. Rep. 2016, 6, 32513. [Google Scholar] [CrossRef] [Green Version]
  2. Capelli, C.; Brisighelli, F.; Scarnicci, F.; Arredi, B.; Caglia’, A.; Vetrugno, G.; Tofanelli, S.; Onofri, V.; Tagliabracci, A.; Paoli, G.; et al. Y chromosome genetic variation in the Italian peninsula is clinal and supports an admixture model for the Mesolithic–Neolithic encounter. Mol. Phylogenet. Evol. 2007, 44, 228–239. [Google Scholar] [CrossRef]
  3. Brisighelli, F.; Álvarez-Iglesias, V.; Fondevila, M.; Blanco-Verea, A.; Carracedo, Á.; Pascali, V.L.; Capelli, C.; Salas, A. Uniparental Markers of Contemporary Italian Population Reveals Details on Its Pre-Roman Heritage. PLoS ONE 2012, 7, e50794. [Google Scholar] [CrossRef] [Green Version]
  4. Boattini, A.; Martinez-Cruz, B.; Sarno, S.; Harmant, C.; Useli, A.; Sanz, P.; Yang-Yao, D.; Manry, J.; Ciani, G.; Luiselli, D.; et al. Uniparental Markers in Italy Reveal a Sex-Biased Genetic Structure and Different Historical Strata. PLoS ONE 2013, 8, e65441. [Google Scholar] [CrossRef] [Green Version]
  5. Di Gaetano, C.; Voglino, F.; Guarrera, S.; Fiorito, G.; Rosa, F.; Di Blasio, A.M.; Manzini, P.; Dianzani, I.; Betti, M.; Cusi, D.; et al. An Overview of the Genetic Structure within the Italian Population from Genome-Wide Data. PLoS ONE 2012, 7, e43759. [Google Scholar] [CrossRef]
  6. Sarno, S.; Boattini, A.; Carta, M.; Ferri, G.; Alù, M.; Yao, D.Y.; Ciani, G.; Pettener, D.; Luiselli, D. An Ancient Mediterranean Melting Pot: Investigating the Uniparental Genetic Structure and Population History of Sicily and Southern Italy. PLoS ONE 2014, 9, e96074. [Google Scholar] [CrossRef] [Green Version]
  7. Sarno, S.; Boattini, A.; Pagani, L.; Sazzini, M.; De Fanti, S.; Quagliariello, A.; Ruscone, G.A.G.; Guichard, E.; Ciani, G.; Bortolini, E.; et al. Ancient and recent admixture layers in Sicily and Southern Italy trace multiple migration routes along the Mediterranean. Sci. Rep. 2017, 7, 1984. [Google Scholar] [CrossRef]
  8. Fiorito, G.; Di Gaetano, C.; Guarrera, S.; Rosa, F.; Feldman, M.W.; Piazza, A.; Matullo, G. The Italian genome reflects the history of Europe and the Mediterranean basin. Eur. J. Hum. Genet. 2016, 24, 1056–1062. [Google Scholar] [CrossRef] [Green Version]
  9. Raveane, A.; Aneli, S.; Montinaro, F.; Athanasiadis, G.; Barlera, S.; Birolo, G.; Boncoraglio, G.; Di Blasio, A.M.; Di Gaetano, C.; Pagani, L.; et al. Population structure of modern-day Italians reveals patterns of ancient and archaic ancestries in Southern Europe. Sci. Adv. 2019, 5, eaaw3492. [Google Scholar] [CrossRef] [Green Version]
  10. Sazzini, M.; Abondio, P.; Sarno, S.; Gnecchi-Ruscone, G.A.; Ragno, M.; Giuliani, C.; De Fanti, S.; Ojeda-Granados, C.; Boattini, A.; Marquis, J.; et al. Genomic history of the Italian population recapitulates key evolutionary dynamics of both Continental and Southern Europeans. BMC Biol. 2020, 18, 51. [Google Scholar] [CrossRef]
  11. Sarno, S.; Petrilli, R.; Abondio, P.; De Giovanni, A.; Boattini, A.; Sazzini, M.; De Fanti, S.; Cilli, E.; Ciani, G.; Gentilini, D.; et al. Genetic history of Calabrian Greeks reveals ancient events and long term isolation in the Aspromonte area of Southern Italy. Sci. Rep. 2021, 11, 3045. [Google Scholar] [CrossRef]
  12. De Fanti, S.; Barbieri, C.; Sarno, S.; Sevini, F.; Vianello, D.; Tamm, E.; Metspalu, E.; Van Oven, M.; Hübner, A.; Sazzini, M.; et al. Fine Dissection of Human Mitochondrial DNA Haplogroup HV Lineages Reveals Paleolithic Signatures from European Glacial Refugia. PLoS ONE 2015, 10, e0144391. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Capocasa, M.; Anagnostou, P.; Bachis, V.; Battaggia, C.; Bertoncini, S.; Biondi, G.; Boattini, A.; Boschi, I.; Brisighelli, F.; Caló, C.M.; et al. Linguistic, geographic and genetic isolation: A collaborative study of Italian populations. J. Anthropol. Sci 2014, 92, 201–231. [Google Scholar] [PubMed]
  14. Cavazzuti, C.; Skeates, R.; Millard, A.R.; Nowell, G.; Peterkin, J.; Bernabò Brea, M.; Cardarelli, A.; Salzani, L. Flows of people in villages and large centres in Bronze Age Italy through strontium and oxygen isotopes. PLoS ONE 2019, 14, e0209693. [Google Scholar] [CrossRef]
  15. Cavazzuti, C.; Arena, A. The Bioarchaeology of Social Stratification in Bronze Age Italy. Arheo 2020, 37, 69–106. [Google Scholar]
  16. Iacono, F. Archaeology of Late Bronze Age Interaction and Mobility at the Gates of Europe: People, Things… and Networks around the Southern Adriatic Sea; Bloomsbury Academic: London, UK; New York, NY, USA, 2019; ISBN 978-1-350-03614-7. [Google Scholar]
  17. Vanzetti, A.; Tinè, V. La Calabria dal Neolitico all’età del ferro. In Museo dei Brettii e Degli Enotri. Catalogo Dell’esposizione; Rubbettino: Soveria Mannelli, Italy, 2014; pp. 40–43. ISBN 978-88-498-4233-3. (In Italian) [Google Scholar]
  18. Guzzo, P.G.; Peroni, R.; Bergonzi, G.; Cardarelli, A.; Vagnetti, L. Ricerche Sulla Protostoria Della SIBARITIDE, 1; Publications du Centre Jean Bérard: Napoli, Italy, 1982; ISBN 978-2-903189-17-4. (In Italian) [Google Scholar]
  19. Ardesia, V. La cultura di Rodì-Tindari-Vallelunga in Sicilia: Origini, diffusione e cronologia alla luce dei recenti studi. Parte 1. IpoTESI Preist. 2014, 6, 35–98. (In Italian) [Google Scholar] [CrossRef]
  20. Larocca, F. Grotta Della Monaca (Calabria, Italia Meridionale). Una Miniera Neolitica per l’estrazione Dell’ocra. Rubricatum Rev. Mus. Gavà 2012, 5, 249–256. (In Italian) [Google Scholar]
  21. Fu, Q.; Posth, C.; Hajdinjak, M.; Petr, M.; Mallick, S.; Fernandes, D.; Furtwängler, A.; Haak, W.; Meyer, M.; Mittnik, A.; et al. The genetic history of Ice Age Europe. Nature 2016, 534, 200–205. [Google Scholar] [CrossRef] [Green Version]
  22. Posth, C.; Renaud, G.; Mittnik, A.; Drucker, D.G.; Rougier, H.; Cupillard, C.; Valentin, F.; Thevenet, C.; Furtwängler, A.; Wißing, C.; et al. Pleistocene Mitochondrial Genomes Suggest a Single Major Dispersal of Non-Africans and a Late Glacial Population Turnover in Europe. Curr. Biol. 2016, 26, 827–833. [Google Scholar] [CrossRef] [Green Version]
  23. Emery, M.V.; Duggan, A.T.; Murchie, T.J.; Stark, R.J.; Klunk, J.; Hider, J.; Eaton, K.; Karpinski, E.; Schwarcz, H.P.; Poinar, H.N.; et al. Ancient Roman mitochondrial genomes and isotopes reveal relationships and geographic origins at the local and pan-Mediterranean scales. J. Archaeol. Sci. Rep. 2018, 20, 200–209. [Google Scholar] [CrossRef]
  24. Quarta, G.; LaRocca, F.; D’Elia, M.; Gaballo, V.; Macchia, M.; Palestra, G.; Calcagnile, L. Radiocarbon Dating the Exploitation Phases of the Grotta Della Monaca Cave in Calabria, Southern Italy: A Prehistoric Mine for the Extraction of Iron and Copper. Radiocarbon 2013, 55, 1246–1251. [Google Scholar] [CrossRef]
  25. Levato, C.; Larocca, F. The Prehistoric Iron Mine of Grotta Della Monaca (Calabria, Italy). Anthropol. Praehist. 2015, 126, 25–37. [Google Scholar]
  26. Caricola, I.; Breglia, F.; LaRocca, F.; Hamon, C.; Lemorini, C.; Giligny, F. Prehistoric exploitation of minerals resources. Experimentation and use-wear analysis of grooved stone tools from Grotta della Monaca (Calabria, Italy). Archaeol. Anthropol. Sci. 2020, 12, 259. [Google Scholar] [CrossRef]
  27. Larocca, F. La Miniera Pre-Protostorica di Grotta Della Monaca (Sant’Agata di Esaro—Cosenza); Centro Regionale di Speleologia “Enzo dei Medici”: Roseto Capo Spulico, Italy, 2005. (In Italian) [Google Scholar]
  28. Arena, F.; Larocca, F.; Onisto, N.; Gualdi-Russo, E. Il Sepolcreto Protostorico Di Grotta Della Monaca in Calabria. Aspetti Antropologici; Annali Dell’Università di Ferrara: Ferrara, Italy, 2013. (In Italian) [Google Scholar]
  29. Bietti Sestieri, A.M. L’Italia Nell’età del Bronzo e del Ferro. Dalle Palafitte a ROMOLO (2200-700 a.C.); Carocci Editore: Rome, Italy, 2010; ISBN 88-430-5207-1. (In Italian) [Google Scholar]
  30. Cocchi Genick, D. Criteri di Nomenclatura e di Terminologia Inerente Alla Definizione Delle Forme Vascolari del Neolitico/Eneolitico e del Bronzo/Ferro: Atti del Congresso di Lido di Camaiore, 26–29 Marzo 1998. Vol. 1; Octavo: Firenze, Italy, 1999; ISBN 978-88-8030-190-5. (In Italian) [Google Scholar]
  31. Larocca, F. Dal Pollino all’Orsomarso. L’uso funerario delle cavità naturali in età pre-protostorica. In Il Pollino. Barriera Naturale e Crocevia di Culture; Giornate Internazionali di Archeologia (San Lorenzo Bellizzi, 16–17 Aprile 2016); Colelli, C.E., Larocca, A., Eds.; Collana del Dipartimento di Studi Umanistici, Sezione Archeologia, Università della Calabria: Calabria, Italy, 2018; Volume Ricerche XII, pp. 13–27. (In Italian) [Google Scholar]
  32. Arena, F.; Gualdi-Russo, E. Taphonomy and Post-Depositional Movements of a Bronze Age Mass Grave in the Archaeological Site of Grotta Della Monaca; Annali dell’Università di Ferrara: Ferrara, Italy, 2014. (In Italian) [Google Scholar]
  33. Gamba, C.; Jones, E.R.; Teasdale, M.D.; McLaughlin, R.L.; Gonzalez-Fortes, G.; Mattiangeli, V.; Domboróczki, L.; Kővári, I.; Pap, I.; Anders, A.; et al. Genome flux and stasis in a five millennium transect of European prehistory. Nat. Commun. 2014, 5, 5257. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Hansen, H.B.; Damgaard, P.B.; Margaryan, A.; Stenderup, J.; Lynnerup, N.; Willerslev, E.; Allentoft, M.E. Comparing Ancient DNA Preservation in Petrous Bone and Tooth Cementum. PLoS ONE 2017, 12, e0170940. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Cooper, A.; Poinar, H.N. Ancient DNA: Do It Right or Not at All. Science 2000, 289, 1139. [Google Scholar] [CrossRef]
  36. Gilbert, M.T.P.; Bandelt, H.-J.; Hofreiter, M.; Barnes, I. Assessing ancient DNA studies. Trends Ecol. Evol. 2005, 20, 541–544. [Google Scholar] [CrossRef]
  37. Fulton, T.L. Setting up an Ancient DNA Laboratory. Methods Mol. Biol. 2012, 840, 1–11. [Google Scholar] [CrossRef]
  38. Fulton, T.L.; Shapiro, B. Setting up an Ancient DNA Laboratory. Methods Mol. Biol. 2019, 1963, 1–13. [Google Scholar] [CrossRef]
  39. Dabney, J.; Knapp, M.; Glocke, I.; Gansauge, M.-T.; Weihmann, A.; Nickel, B.; Valdiosera, C.; García, N.; Pääbo, S.; Arsuaga, J.-L.; et al. Complete mitochondrial genome sequence of a Middle Pleistocene cave bear reconstructed from ultrashort DNA fragments. Proc. Natl. Acad. Sci. USA 2013, 110, 15758–15763. [Google Scholar] [CrossRef] [Green Version]
  40. Damgaard, P.B.; Margaryan, A.; Schroeder, H.; Orlando, L.; Willerslev, E.; Allentoft, M.E. Improving access to endogenous DNA in ancient bones and teeth. Sci. Rep. 2015, 5, 11184. [Google Scholar] [CrossRef] [Green Version]
  41. Cilli, E.; Gabanini, G.; Ciucani, M.M.; De Fanti, S.; Serventi, P.; Bazaj, A.; Sarno, S.; Ferri, G.; Fregnani, A.; Cornaglia, G.; et al. A multifaceted approach towards investigating childbirth deaths in double burials: Anthropology, paleopathology and ancient DNA. J. Archaeol. Sci. 2020, 122, 105219. [Google Scholar] [CrossRef]
  42. Carøe, C.; Gopalakrishnan, S.; Vinner, L.; Mak, S.S.T.; Sinding, M.H.S.; Samaniego, J.A.; Wales, N.; Sicheritz-Pontén, T.; Gilbert, M.T.P. Single-tube library preparation for degraded DNA. Methods Ecol. Evol. 2018, 9, 410–419. [Google Scholar] [CrossRef]
  43. Maricic, T.; Whitten, M.; Pääbo, S. Multiplexed DNA Sequence Capture of Mitochondrial Genomes Using PCR Products. PLoS ONE 2010, 5, e14004. [Google Scholar] [CrossRef] [PubMed]
  44. Andrews, S. FastQC. A Quality Control Tool for High Throughput Sequence Data. Available online: https://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (accessed on 22 April 2021).
  45. Schubert, M.; Lindgreen, S.; Orlando, L. AdapterRemoval v2: Rapid adapter trimming, identification, and read merging. BMC Res. Notes 2016, 9, 88. [Google Scholar] [CrossRef] [Green Version]
  46. Andrews, R.M.; Kubacka, I.; Chinnery, P.F.; Lightowlers, R.N.; Turnbull, D.M.; Howell, N. Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA. Nat. Genet. 1999, 23, 147. [Google Scholar] [CrossRef]
  47. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R.; 1000 Genome Project Data Processing Subgroup. The Sequence Alignment/Map Format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Picard Toolkit. Available online: https://broadinstitute.github.io/picard/ (accessed on 22 April 2021).
  50. McKenna, A.; Hanna, M.; Banks, E.; Sivachenko, A.; Cibulskis, K.; Kernytsky, A.; Garimella, K.; Altshuler, D.; Gabriel, S.; Daly, M.; et al. The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010, 20, 1297–1303. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Jónsson, H.; Ginolhac, A.; Schubert, M.; Johnson, P.L.F.; Orlando, L. mapDamage2.0: Fast approximate Bayesian estimates of ancient DNA damage parameters. Bioinformatics 2013, 29, 1682–1684. [Google Scholar] [CrossRef] [PubMed]
  52. Renaud, G.; Slon, V.; Duggan, A.T.; Kelso, J. Schmutzi: Estimation of contamination and endogenous mitochondrial consensus calling for ancient DNA. Genome Biol. 2015, 16, 224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Namouchi, A. SnpToolkit. Available online: https://github.com/Amine-Namouchi/snpToolkit (accessed on 22 April 2021).
  54. Vianello, D.; Sevini, F.; Castellani, G.; Lomartire, L.; Capri, M.; Franceschi, C. HAPLOFIND: A New Method for High-Throughput mtDNA Haplogroup Assignment. Hum. Mutat. 2013, 34, 1189–1194. [Google Scholar] [CrossRef] [PubMed]
  55. Weissensteiner, H.; Pacher, D.; Kloss-Brandstätter, A.; Forer, L.; Specht, G.; Bandelt, H.-J.; Kronenberg, F.; Salas, A.; Schönherr, S. HaploGrep 2: Mitochondrial haplogroup classification in the era of high-throughput sequencing. Nucleic Acids Res. 2016, 44, W58–W63. [Google Scholar] [CrossRef] [PubMed]
  56. Van Oven, M. PhyloTree Build 17: Growing the human mitochondrial DNA tree. Forensic Sci. Int. Genet. Suppl. Ser. 2015, 5, e392–e394. [Google Scholar] [CrossRef] [Green Version]
  57. DNA Alignment Software. Available online: https://www.fluxus-engineering.com/align.htm (accessed on 22 April 2021).
  58. Van Oven, M.; Kayser, M. Updated comprehensive phylogenetic tree of global human mitochondrial DNA variation. Hum. Mutat. 2009, 30, E386–E394. [Google Scholar] [CrossRef]
  59. Network. Available online: https://www.fluxus-engineering.com/sharenet.htm (accessed on 22 April 2021).
  60. Drummond, A.J.; Suchard, M.A.; Xie, D.; Rambaut, A. Bayesian Phylogenetics with BEAUti and the BEAST 1.7. Mol. Biol. Evol. 2012, 29, 1969–1973. [Google Scholar] [CrossRef] [Green Version]
  61. Tamura, K.; Peterson, D.; Peterson, N.; Stecher, G.; Nei, M.; Kumar, S. MEGA5: Molecular Evolutionary Genetics Analysis Using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Mol. Biol. Evol. 2011, 28, 2731–2739. [Google Scholar] [CrossRef] [Green Version]
  62. Rambaut, A. Figtree v1.4.2, A Graphical Viewer of Phylogenetic Trees. Available online: tree.bio.ed.ac.uk/software/figtree/ (accessed on 22 April 2021).
  63. Fernandes, D.M.; Mittnik, A.; Olalde, I.; Lazaridis, I.; Cheronet, O.; Rohland, N.; Mallick, S.; Bernardos, R.; Broomandkhoshbacht, N.; Carlsson, J.; et al. The spread of steppe and Iranian-related ancestry in the islands of the western Mediterranean. Nat. Ecol. Evol. 2020, 4, 334–345. [Google Scholar] [CrossRef]
  64. Modi, A.; Tassi, F.; Susca, R.R.; Vai, S.; Rizzi, E.; Bellis, G.D.; Lugliè, C.; Gonzalez Fortes, G.; Lari, M.; Barbujani, G.; et al. Complete mitochondrial sequences from Mesolithic Sardinia. Sci. Rep. 2017, 7, 42869. [Google Scholar] [CrossRef] [Green Version]
  65. Olivieri, A.; Sidore, C.; Achilli, A.; Angius, A.; Posth, C.; Furtwängler, A.; Brandini, S.; Capodiferro, M.R.; Gandini, F.; Zoledziewska, M.; et al. Mitogenome Diversity in Sardinians: A Genetic Window onto an Island’s Past. Mol. Biol. Evol. 2017, 34, 1230–1239. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Olalde, I.; Brace, S.; Allentoft, M.E.; Armit, I.; Kristiansen, K.; Booth, T.; Rohland, N.; Mallick, S.; Szécsényi-Nagy, A.; Mittnik, A.; et al. The Beaker phenomenon and the genomic transformation of northwest Europe. Nature 2018, 555, 190–196. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Catalano, G.; Lo Vetro, D.; Fabbri, P.F.; Mallick, S.; Reich, D.; Rohland, N.; Sineo, L.; Mathieson, I.; Martini, F. Late Upper Palaeolithic hunter-gatherers in the Central Mediterranean: New archaeological and genetic data from the Late Epigravettian burial Oriente C (Favignana, Sicily). Quat. Int. 2020, 537, 24–32. [Google Scholar] [CrossRef]
  68. Marcus, J.H.; Posth, C.; Ringbauer, H.; Lai, L.; Skeates, R.; Sidore, C.; Beckett, J.; Furtwängler, A.; Olivieri, A.; Chiang, C.W.K.; et al. Genetic history from the Middle Neolithic to present on the Mediterranean island of Sardinia. Nat. Commun. 2020, 11, 939. [Google Scholar] [CrossRef] [PubMed]
  69. Modi, A.; Catalano, G.; D’Amore, G.; Di Marco, S.; Lari, M.; Sineo, L.; Caramelli, D. Paleogenetic and morphometric analysis of a Mesolithic individual from Grotta d’Oriente: An oldest genetic legacy for the first modern humans in Sicily. Quat. Sci. Rev. 2020, 248, 106603. [Google Scholar] [CrossRef]
  70. Olalde, I.; Posth, C. Latest trends in archaeogenetic research of west Eurasians. Curr. Opin. Genet. Dev. 2020, 62, 36–43. [Google Scholar] [CrossRef]
  71. Pinhasi, R.; Fernandes, D.; Sirak, K.; Novak, M.; Connell, S.; Alpaslan-Roodenberg, S.; Gerritsen, F.; Moiseyev, V.; Gromov, A.; Raczky, P.; et al. Optimal Ancient DNA Yields from the Inner Ear Part of the Human Petrous Bone. PLoS ONE 2015, 10, e0129102. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Brotherton, P.; Haak, W.; Templeton, J.; Brandt, G.; Soubrier, J.; Jane Adler, C.; Richards, S.M.; Sarkissian, C.D.; Ganslmeier, R.; Friederich, S.; et al. Neolithic mitochondrial haplogroup H genomes and the genetic origins of Europeans. Nat. Commun. 2013, 4, 1764. [Google Scholar] [CrossRef] [Green Version]
  73. Brandt, G.; Haak, W.; Adler, C.J.; Roth, C.; Szécsényi-Nagy, A.; Karimnia, S.; Möller-Rieker, S.; Meller, H.; Ganslmeier, R.; Friederich, S.; et al. Ancient DNA Reveals Key Stages in the Formation of Central European Mitochondrial Genetic Diversity. Science 2013, 342, 257–261. [Google Scholar] [CrossRef] [Green Version]
  74. Hernández, C.L.; Dugoujon, J.M.; Novelletto, A.; Rodríguez, J.N.; Cuesta, P.; Calderón, R. The distribution of mitochondrial DNA haplogroup H in southern Iberia indicates ancient human genetic exchanges along the western edge of the Mediterranean. BMC Genet. 2017, 18, 46. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Álvarez-Iglesias, V.; Mosquera-Miguel, A.; Cerezo, M.; Quintáns, B.; Zarrabeitia, M.T.; Cuscó, I.; Lareu, M.V.; García, O.; Pérez-Jurado, L.; Carracedo, Á.; et al. New Population and Phylogenetic Features of the Internal Variation within Mitochondrial DNA Macro-Haplogroup R0. PLoS ONE 2009, 4, e5112. [Google Scholar] [CrossRef] [PubMed]
  76. Soares, P.; Achilli, A.; Semino, O.; Davies, W.; Macaulay, V.; Bandelt, H.-J.; Torroni, A.; Richards, M.B. The Archaeogenetics of Europe. Curr. Biol. 2010, 20, R174–R183. [Google Scholar] [CrossRef] [Green Version]
  77. Torroni, A.; Bandelt, H.-J.; D’Urbano, L.; Lahermo, P.; Moral, P.; Sellitto, D.; Rengo, C.; Forster, P.; Savontaus, M.-L.; Bonné-Tamir, B.; et al. mtDNA Analysis Reveals a Major Late Paleolithic Population Expansion from Southwestern to Northeastern Europe. Am. J. Hum. Genet. 1998, 62, 1137–1152. [Google Scholar] [CrossRef] [Green Version]
  78. Torroni, A.; Bandelt, H.-J.; Macaulay, V.; Richards, M.; Cruciani, F.; Rengo, C.; Martinez-Cabrera, V.; Villems, R.; Kivisild, T.; Metspalu, E.; et al. A Signal, from Human mtDNA, of Postglacial Recolonization in Europe. Am. J. Hum. Genet. 2001, 69, 844–852. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Achilli, A.; Rengo, C.; Magri, C.; Battaglia, V.; Olivieri, A.; Scozzari, R.; Cruciani, F.; Zeviani, M.; Briem, E.; Carelli, V.; et al. The Molecular Dissection of mtDNA Haplogroup H Confirms That the Franco-Cantabrian Glacial Refuge Was a Major Source for the European Gene Pool. Am. J. Hum. Genet. 2004, 75, 910–918. [Google Scholar] [CrossRef]
  80. Pereira, L. High-resolution mtDNA evidence for the late-glacial resettlement of Europe from an Iberian refugium. Genome Res. 2005, 15, 19–24. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Soares, P.; Ermini, L.; Thomson, N.; Mormina, M.; Rito, T.; Röhl, A.; Salas, A.; Oppenheimer, S.; Macaulay, V.; Richards, M.B. Correcting for Purifying Selection: An Improved Human Mitochondrial Molecular Clock. Am. J. Hum. Genet. 2009, 84, 740–759. [Google Scholar] [CrossRef] [Green Version]
  82. Lipson, M.; Szécsényi-Nagy, A.; Mallick, S.; Pósa, A.; Stégmár, B.; Keerl, V.; Rohland, N.; Stewardson, K.; Ferry, M.; Michel, M.; et al. Parallel palaeogenomic transects reveal complex genetic history of early European farmers. Nature 2017, 551, 368–372. [Google Scholar] [CrossRef]
  83. Mathieson, I.; Lazaridis, I.; Rohland, N.; Mallick, S.; Patterson, N.; Roodenberg, S.A.; Harney, E.; Stewardson, K.; Fernandes, D.; Novak, M.; et al. Genome-wide patterns of selection in 230 ancient Eurasians. Nature 2015, 528, 499–503. [Google Scholar] [CrossRef] [Green Version]
  84. Mathieson, I.; Alpaslan-Roodenberg, S.; Posth, C.; Szécsényi-Nagy, A.; Rohland, N.; Mallick, S.; Olalde, I.; Broomandkhoshbacht, N.; Candilio, F.; Cheronet, O.; et al. The genomic history of southeastern Europe. Nature 2018, 555, 197–203. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Bernabò Brea, L. La Sicilia Prima dei Greci; Il Saggiatore: Milano, Italy, 1961. (In Italian) [Google Scholar]
  86. Bietti Sestieri, A.M. Implicazioni Del Concetto Di Territorio in Situazioni Culturali Complesse: Le Isole Eolie Nell’età Del Bronzo. Dialoghi Archeol. 1982, 2, 39–60. (In Italian) [Google Scholar]
  87. Cazzella, A.; Recchia, G. Malta, Sicily, Aeolian Islands and southern Italy during the Bronze Age: The meaning of a changing relationship. In Exchange Networks and Local Transformation: Interaction and Local Change in Europe and the Mediterranean from the Bronze Age to the Iron Age; Oxbow: Oxford, UK, 2013; pp. 80–91. [Google Scholar]
  88. Pacciarelli, M.; Scarano, T.; Crispino, A. The Transition between the Copper and Bronze Ages in Southern Italy and Sicily. In Proceedings of the 7th Archaeological Conference of Central Germany, Halle, Germany, 23–26 October 2014; Landesmuseums für Vorgeschichte Halle: Halle, Germany, 2015. [Google Scholar]
  89. Gori, M.; Revello Lami, M.; Pintucci, A. Editorial: Practices, Representations and Meanings of Human Mobility in Archaeology. Ex Novo 2018, 3, 1–6. [Google Scholar] [CrossRef]
  90. Marino, S. Across the Strait. New evidence on cultural interconnections and exchanges between Calabria and Sicily during the early Bronze Age. J. Archaeol. Sci. Rep. 2020, 29, 102164. [Google Scholar] [CrossRef]
Figure 1. Location of Grotta della Monaca. The entrance of the cave is marked by a yellow circle on the bottom right image. The scheme of the cave with the m5v area is represented at the top right.
Figure 1. Location of Grotta della Monaca. The entrance of the cave is marked by a yellow circle on the bottom right image. The scheme of the cave with the m5v area is represented at the top right.
Genes 12 00636 g001
Figure 2. Median Joining Network based on the complete GdM24 mitogenome (shown in red) compared to a dataset of ancient reference (in white), H1 (grey) and H1e (black) mitogenomes (see Table S4 for dataset).
Figure 2. Median Joining Network based on the complete GdM24 mitogenome (shown in red) compared to a dataset of ancient reference (in white), H1 (grey) and H1e (black) mitogenomes (see Table S4 for dataset).
Genes 12 00636 g002
Figure 3. Phylogenetic tree reconstructed with BEAST using dated ancient reference samples as calibration points (Tables S4 and S5).
Figure 3. Phylogenetic tree reconstructed with BEAST using dated ancient reference samples as calibration points (Tables S4 and S5).
Genes 12 00636 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fontani, F.; Cilli, E.; Arena, F.; Sarno, S.; Modi, A.; De Fanti, S.; Andrews, A.J.; Latorre, A.; Abondio, P.; Larocca, F.; et al. First Bronze Age Human Mitogenomes from Calabria (Grotta Della Monaca, Southern Italy). Genes 2021, 12, 636. https://doi.org/10.3390/genes12050636

AMA Style

Fontani F, Cilli E, Arena F, Sarno S, Modi A, De Fanti S, Andrews AJ, Latorre A, Abondio P, Larocca F, et al. First Bronze Age Human Mitogenomes from Calabria (Grotta Della Monaca, Southern Italy). Genes. 2021; 12(5):636. https://doi.org/10.3390/genes12050636

Chicago/Turabian Style

Fontani, Francesco, Elisabetta Cilli, Fabiola Arena, Stefania Sarno, Alessandra Modi, Sara De Fanti, Adam Jon Andrews, Adriana Latorre, Paolo Abondio, Felice Larocca, and et al. 2021. "First Bronze Age Human Mitogenomes from Calabria (Grotta Della Monaca, Southern Italy)" Genes 12, no. 5: 636. https://doi.org/10.3390/genes12050636

APA Style

Fontani, F., Cilli, E., Arena, F., Sarno, S., Modi, A., De Fanti, S., Andrews, A. J., Latorre, A., Abondio, P., Larocca, F., Lari, M., Caramelli, D., Gualdi-Russo, E., & Luiselli, D. (2021). First Bronze Age Human Mitogenomes from Calabria (Grotta Della Monaca, Southern Italy). Genes, 12(5), 636. https://doi.org/10.3390/genes12050636

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop