The arrival of European colonists to the Americas, beginning in the 15^th century, contributed to the spread of new viruses amongst Indigenous people. This led to massive outbreaks of disease, and millions of deaths that caused an important Native population to collapse. The exact viruses that caused these outbreaks are unknown, but smallpox, measles, and mumps are all suspected.

During these times, traders and colonists forcibly enslaved and displaced millions of people mainly from the West Coast of Africa to the Americas. The cruel, unsanitary, and overcrowded conditions on ships transporting these people across the Atlantic contributed to the spread of infectious diseases onboard. Once on land, infectious diseases spread quickly, partly due to the poor conditions that enslaved and ndigenous people were made to endure. Native people were also immunologically naïve to the newly introduced pathogens, making them susceptible to severe or fatal outcomes. The new field of paleovirology may help scientists identify the viruses that were circulating in the first years of colonization and trace how viruses arrived in the Americas.

Using next-generation DNA sequencing and other cutting-edge techniques, Guzmán-Solís et al. extracted and enriched viral DNA from skeletal remains dating back to the 16^th century. These remains were found in mass graves that were used to bury epidemic victims at a colonial hospital and chapel in what is now Mexico City. The experiments identified two viruses, human parvovirus B19 and a human hepatitis B virus. These viral genomes were recovered from human remains of first-generation African people in Mexico, as well as an individual who was an Indigenous person.

Although the genetic material of these ancient viruses resembled pathogens that originated in Africa, the study did not determine if the victims died from these viruses or another cause. On the other hand, the results indicate that viruses frequently found in modern Africa were circulating in the Americas during the slave trade period of Mexico. Finally, the results provide evidence that colonists who forcibly brought African people to the Americas participated in the introduction of viruses to Mexico. This constant influx of viruses from the old world, led to dramatic declines in the populations of Indigenous people in the Americas.

European colonization in the Americas resulted in a frequent genetic exchange mainly between Native American populations, Europeans, and Africans (Aguirre-Beltrán, 2005; Rotimi et al., 2016; Salas et al., 2004). Along with human migrations, numerous new species were introduced to the Americas including bacterial and viral pathogens, which played a major role in the dramatic population collapse that afflicted the immunologically naïve Indigenous inhabitants (Acuña-Soto et al., 2004; Lindo et al., 2016). Among these pathogens, viral diseases, such as smallpox, measles, and mumps, have been proposed to be responsible for many of the devastating epidemics during the Colonial period (Acuña-Soto et al., 2004). Remarkably, the pathogen(s) responsible for the deadliest epidemics reported in New Spain (the Spanish viceroyalty that corresponds to Mexico, Central America, and the current US southwest states) remains unknown and is thought to have caused millions of deaths during the 16th century (Acuña-Soto et al., 2004). Indigenous populations were drastically affected by these mysterious epidemics, generically referred to as Cocoliztli (‘pest’ in Nahuatl), followed by Africans and to a lesser extent European people (Acuña-Soto et al., 2004; Malvido and Viesca, 1982; Somolinos d’Árdois, 1982). Accounts of the 1576 Cocoliztli epidemic were described in autopsy reports of victims treated at the ‘Hospital Real de San José de los Naturales’ (HSJN) (Malvido and Viesca, 1982; Wesp, 2017), the first hospital in Mexico dedicated specifically to treat the Indigenous population (Malvido and Viesca, 1982; Wesp, 2017; Figure 1a and b).

Figure 1 with 7 supplements see all

Download asset Open asset

The study of ancient viral genomes has revealed important insights into the evolution of specific viral families (Barquera et al., 2020; Duggan et al., 2016; Düx et al., 2020; Kahila Bar-Gal et al., 2012; Krause-Kyora et al., 2018; Mühlemann et al., 2018a; Mühlemann et al., 2018a; Mühlemann et al., 2018b; Neukamm et al., 2020; Pajer et al., 2017; Patterson Ross et al., 2018; Xiao et al., 2013), as well as their interaction with human populations (Spyrou et al., 2019). To explore the presence of viral pathogens in circulation during epidemic periods in New Spain, we leveraged the vast historical and archeological information available for the Colonial HSJN. These include the skeletal remains of over 600 individuals recovered from mass burials associated with the hospital’s architectural remnants (Figure 1b). Many of these remains were retrieved from burial contexts suggestive of an urgent and simultaneous disposal of the bodies, as in the case of an epidemic (Meza, 2013; Wesp, 2017). Prior bioarcheological research has shown that the remains of a number of individuals in the HSJN collection displayed dental modifications and/or morphological indicators typical of African ancestry (Meza, 2013), consistent with historical and archeological research that documents the presence of a large number of both free and enslaved Africans and their descendants in Colonial Mexico (Aguirre-Beltrán, 2005). Indeed, a recent paleogenomics study reported a sub-Saharan African origin of three individuals from this collection (Barquera et al., 2020).

Here we describe the recovery and characterization of viral pathogens that circulated in New Spain during the Colonial period, using ancient DNA (aDNA) techniques (Figure 1—figure supplement 1). For this work, we sampled skeletal human remains recovered from the HSJN where archeological context suggest victims of epidemics were buried (Meza, 2013) and from ‘La Concepcion’ chapel, one of the first catholic conversion centers in New Spain (Moreno-Cabrera et al., 2015; Figure 1a). We report the reconstruction of ancient hepatitis B virus (HBV) and human parvovirus B19 (B19V) genomes recovered from these remains, providing a direct molecular evidence of human viral pathogens of African origin being introduced to New Spain during the transatlantic slave trade.

We sampled the skeletal remains from two archeological sites, a Colonial Hospital and a Colonial chapel in Mexico City (Figure 1a and b). For the HSJN, 21 dental samples (premolar and molar teeth) were selected based on previous morphometric analyses and dental modifications that suggested an African ancestry (Hernández-Lopez and Negrete, 2012; Karam-Tapia, 2012; Meza, 2013; Ruíz-Albarrán, 2012). The African presence in the Indigenous Hospital might reflect an urgent response to an epidemic outbreak since hospitals treated patients regardless of the origin of the affected individuals during serious public health crises (Meza, 2013). Dental samples of five additional individuals were selected (based on their conservation state) from ‘La Concepción’ chapel (COY), which is located 10 km south of the HSJN in Coyoacán, a Pre-Hispanic Indigenous neighborhood that became the first Spanish settlement in Mexico City after the fall of Tenochtitlan (Moreno-Cabrera et al., 2015). Following strict aDNA protocols, we processed these dental samples to isolate aDNA for next-generation sequencing (NGS) (Figure 1—figure supplement 1, Materials and methods). Tooth roots (which are vascularized) can be a good source of pathogen DNA (Key et al., 2017), especially in the case of viruses that are widespread in the bloodstream during systemic infection. Accordingly, a number of previous studies have successfully recovered ancient viral DNA from tooth roots (Barquera et al., 2020; Krause-Kyora et al., 2018; Mühlemann et al., 2020; Mühlemann et al., 2018a; Mühlemann et al., 2018b).

Metagenomic analyses with MALT (Herbig et al., 2018) (Materials and methods) on the NGS data using the Viral NCBI RefSeq database as a reference (Pruitt et al., 2007) revealed 17 samples containing at least one normalized hit to viral DNA (abundances were normalized to the smallest library since each sample had different number of reads) (Materials and methods), particularly similar to Hepadnaviridae, Herpesviridae, Parvoviridae, and Poxviridae viral families (Figure 1c, Figure 1—figure supplement 2, Materials and methods). These viral hits revealed the potential to recover ancient viral genomes from these samples. We selected 12 samples for further screening (Figure 1c, Figure 1—figure supplement 3) based on the DNA concentration of the NGS library and the quality of the hits to a clinically important virus (HBV, B19V, papillomavirus, smallpox).

To isolate and enrich the viral DNA fraction in the sequencing libraries, biotinylated single-stranded RNA probes designed to capture sequences from diverse human viral pathogens were synthesized (Supplementary file 1A). The selection of the viruses included in the capture design considered the following criteria: (1) DNA viruses previously retrieved from archeological human remains (i.e., hepatitis B virus, human parvovirus B19, variola virus), (2) representative viruses from families capable of integrating into the human genome (i.e., Herpesviridae, Papillomaviridae, Polyomaviridae, Circoviridae), or (3) RNA viruses with a DNA intermediate (i.e., Retroviridae). Given the size constraints of the probe kit, only a couple of genes were selected from some viral families (Materials and methods, Supplementary file 1A). Additionally, a virus-negative aDNA library, which showed no hits to any viral family included in the capture assay (except for a frequent Poxviridae-like region identified as an Alu repeat; Tithi et al., 2018) was captured and sequenced as a negative control (HSJN177) to estimate the efficiency of our capture assay. Only one post-capture library had an ~100-fold increase of Hepadnaviridae-like hits (HBV), while three more libraries had an ~50–200-fold increase of Parvoviridae-like hits (B19V) (Figure 1c, Supplementary file 1B), compared to their corresponding pre-capture libraries (Materials and methods). In contrast, the captured negative control (HSJN177) presented a negligible enrichment of these viral hits (Figure 1c, Supplementary file 1B).

Independently, a metagenomic analysis using Kraken2 (Wood et al., 2019) and Pavian (Breitwieser and Salzberg, 2020) was performed on the non-human (unmapped) reads as part of a different study (Bravo-Lopez et al., 2020). Our samples presented bacterial constituents of human oral and soil microbiota at different proportions between the samples (Figure 1—figure supplements 4–7). Although no lethal bacterial pathogen was retrieved, some ancient dental pathogens (Tannerella forsythia) were reconstructed and described in more detail by Bravo-Lopez et al., 2020 (Figure 1—figure supplements 4–7).

We verified the authenticity of the reads mapped to HBV (BWA) or B19V (BWA/blastn) in the enriched libraries (Materials and methods) by querying the reads against the non-redundant (nr) NCBI database using megaBLAST (Altschul et al., 1990). This step was performed to avoid including in the genome assembly reads that were mapped by BWA or blastn as HBV or B19V, but with a similar identity to a different taxon in the nr database (and absent in DS1; Materials and methods). Therefore, we only retained reads for which the top hit was to either B19V or HBV (Supplementary file 1C). To confirm the ancient origin of these viral reads, we evaluated the misincorporation damage patterns using the program mapDamage 2.0 (Jónsson et al., 2013), which revealed an accumulation of C to T mutations towards their 5′ terminal site with an almost symmetrical G to A pattern on the 3′ end (Figure 2a, Figure 2—figure supplement 1), as expected for aDNA (Briggs et al., 2007). Three ancient B19V genomes were reconstructed (Figure 2b, Supplementary file 1C) with sequence coverages between 92.37% and 99.1%, and average depths of 2.98–15.36× along their single-stranded DNA (ssDNA) coding region, which excludes the double-stranded DNA (dsDNA) hairpin regions at each end of the genome (Luo and Qiu, 2015). These dsDNA inverse terminal repeats (ITRs) displayed considerably higher depth values (<218×) compared to the coding region consistent with the better postmortem preservation of dsDNA compared to ssDNA (Lindahl, 1993; Figure 2b). In addition, we reconstructed one ancient HBV genome (Figure 2c, Supplementary file 1C) at 30.8× average depth and with a sequence coverage of 89.9%, including its ssDNA region at a reduced depth (<10×).

Figure 2 with 2 supplements see all

Download asset Open asset

The reconstructed ancient HBV genome shows a 6 nucleotide (nt) insertion in the core gene, which is characteristic of the genotype A (Kramvis, 2014). Further phylogenetic analyses (Materials and methods) revealed that the Colonial HBV genome clustered with modern sequences corresponding to sub-genotype A4 (previously named A6) (Pourkarim et al., 2014; Figure 3a, Figure 3—figure supplement 1). The genotype A (HBV/GtA) has a broad diversity in Africa reflecting its long history in this continent (Kostaki et al., 2018; Kramvis, 2014), while the sub-genotype A4 has been recovered uniquely from African individuals in Belgium (Pourkarim et al., 2010) and has never been found in the Americas. Regarding the three Colonial B19V genomes from individuals HSJN240, COYC4, and HSJNC81 (C81A), these were phylogenetically closer to modern B19V sequences belonging to genotype 3 (Figure 3b, Figure 3—figure supplements 2–3). This B19V genotype is divided into two sub-genotypes: 3a, which is mostly found in Africa, and 3b, which is proposed to have spread outside Africa recently (Hübschen et al., 2009). The viral sequences from the individuals HSJN240 and COYC4 are similar to sub-genotype 3b genomes sampled from immigrants (Morocco, Egypt, and Turkey) in Germany (Schneider et al., 2008; Figure 3b, Figure 3—figure supplements 2–3) while the sequence of the individual HSJNC81 is more similar to a divergent sub-genotype 3a strain (Figure 3b, Figure 3—figure supplements 2–3) retrieved from a child with severe anemia born in France (Nguyen et al., 1999). These observations support the African origin of the reconstructed Colonial viral genomes.

Figure 3 with 7 supplements see all

Download asset Open asset

In order to use our reconstructed viral genomes as molecular fossils to recalibrate each virus phylogeny and perform evolutionary inferences, we first needed to estimate if the phylogenetic relationships among B19V or HBV genomes had a temporal structure (i.e., sufficient genetic change between sampling times to reconstruct a statistical relationship between genetic divergence and time) (Rambaut et al., 2016). In the context of viruses, temporal structure is canonically tested with a root-to-tip distance and date randomization analyses (see Firth et al., 2010; Rieux and Balloux, 2016). Similarly to previous studies (Krause-Kyora et al., 2018; Patterson Ross et al., 2018), we found little or no temporal structure for this HBV phylogeny containing all genotypes (R² = 0.1351; correlation coefficient = 0.3676) (Figure 3—figure supplement 5a-c). The complex evolution of HBV may not be prone to an appropriate genetic dating since multiple inter-genotype recombination and cross-species transmission (Human-Ape) events (Krause-Kyora et al., 2018) occurred throughout its evolution. Since the entire genotype A has been identified as a recombinant genotype before (Mühlemann et al., 2018a), we analyzed it independently and identified a stronger temporal signal within this genotype (R² = 0.722; correlation coefficient = 0.8498) (Figure 3—figure supplement 5d-f). In the case of B19V, we identified a temporal structure when including all three genotypes (R² = 0.3837; correlation coefficient = 0.6194) (Figure 3—figure supplement 6a-c), in agreement with previous studies (Mühlemann et al., 2018b). Furthermore, we corroborated this temporal structure was not an artifact by a set of tip-dated randomized analyses (Rieux and Balloux, 2016), where none of the 95% highest posterior density (HPD) intervals of the clock rate overlapped with the correctly dated dataset (Figure 3—figure supplement 7).

Given its strong temporal structure, we then performed a dated coalescent phylogenetic analysis for B19V (Supplementary file 1D). We inferred a median substitution rate for B19V of 1.03 × 10^–5 (95% HPD: 8.66 × 10^–6–1.21 × 10^–5) s/s/y under a strict clock and a constant population prior, and a substitution rate of 2.62 × 10^–5 (95% HPD: 1.50 × 10^–5–3.98 × 10^–5) s/s/y under a relaxed log normal clock and a constant population prior. The divergence times from the most recent common ancestor of genotypes 1, 2, and 3 under a strict clock were 7.19 (95% HPD: 6.98–7.46), 2.11 (95% HPD: 1.83–2.51), and 3.64 (95% HPD: 3.04–4.33) ka, respectively. The inferred substitution rates and divergence times from the most recent common ancestor for genotypes 1 and 2 were similar to previous estimations (Mühlemann et al., 2018b) that included much older sequences, while the divergence of genotype 3 was subtly older since no other ancient genotype 3 had been reported previously.

Next, we used the shotgun data generated to determine the mitochondrial haplogroup of the hosts, as well as their autosomal genetic ancestry using the 1000 Genomes Project (1000 Genomes Project Consortium, 2015) as a reference panel (Figure 4a, Supplementary file 2A). The nuclear genetic ancestry analysis showed that all three HSJN individuals, from which the reconstructed viral genomes were isolated, fall within African genetic variation in a principal component analysis (PCA) plot (Figure 4a), while their mitochondrial aDNA belong to the L haplogroup, which has high frequency in African populations (Supplementary file 2A, Figure 2—figure supplement 2). Additionally, we performed ⁸⁷Sr/⁸⁶Sr isotopic analysis on two of the HSJN individuals using tooth enamel as well as phalanx (HSJN240) or parietal bone (HSJNC81) to provide insights on the places of birth (adult enamel) and where the last years of life were spent (phalanx/parietal). The ⁸⁷Sr/⁸⁶Sr ratios measured on the enamel of the individual HSJNC81 (0.71098) and HSJN240 (0.71109) are similar to average ⁸⁷Sr/⁸⁶Sr ratios found in soils and rocks from West Africa (average of 0.71044, Figure 4—figure supplement 1, Supplementary file 2), as well as to ⁸⁷Sr/⁸⁶Sr ratios described in first-generation Africans in the Americas (Barquera et al., 2020; Bastos et al., 2016; Fricke et al., 2020; Price et al., 2012; Schroeder et al., 2009). In contrast, the ⁸⁷Sr/⁸⁶Sr ratios on the parietal and phalange bones from the HSJNC81 (0.70672) and HSJN240 (0.70755) show lower values similar to those observed in the Trans Mexican Volcanic Belt where the Mexico City Valley is located (0.70420–0.70550, Figure 4—figure supplement 1, Supplementary file 2). Moreover, radiocarbon dating of HSJN240 (1442–1608 CE, years calibrated for 1σ) and HSJN194 (1472–1625 CE, years calibrated for 1σ) (Supplementary file 2A, Figure 4—figure supplement 2) indicates that these individuals arrived during the first decades of the Colonial period, when the number of enslaved individuals arriving from Africa was particularly high (Aguirre-Beltrán, 2005).

Figure 4 with 3 supplements see all

Download asset Open asset

Strikingly, Colonial individual COYC4, who was also infected with an African B19V strain, clusters with present-day Mexicans and Peruvians from the 1000 Genomes Project (Figure 4a). An ADMIXTURE (Alexander and Lange, 2011) analysis with these data confirmed a Native American genetic component (Figure 4b), as expected for an indigenous individual. The B19V ancient genome from the individual COYC4 is the first genotype 3 genome obtained from a non-African individual and suggests that following the introduction from Africa, the virus (B19V) spread and infected people of different ancestries during the Colonial period.

In this study, we reconstructed one HBV and three B19V ancient genomes from four different individuals using NGS, metagenomics, and in-solution targeted enrichment methods (Figure 2b,c, Figure 1—figure supplement 1). Several lines of evidence support the ancient nature of these viral sequences, in contrast to environmental contamination or a capture artifact. First, our negative control was not enriched for B19V or HBV hits in our capture sequencing (Figure 1c). For those samples that showed an enrichment in viral sequences after capture, the reads covered the reference genomes almost in their entirety and displayed deamination patterns at the terminal ends of the reads, as expected for aDNA (Figure 2a). Moreover, it is important to notice that B19V and HBV are blood-borne human pathogens that are not present in soil or the environment, and that DNA from these viruses had never been extracted before in the aDNA facilities used for this study.

The recovery of aDNA from B19V, which has a ssDNA genome (with dsDNA terminal repeats), in previous studies (Mühlemann et al., 2018b) as well as in our samples is noteworthy considering the NGS libraries were constructed using dsDNA as a template. Therefore, we would not expect to recover the ssDNA from B19V with this library building method. However, it is known that dsDNA intermediates form during the B19V replication cycle (Ganaie and Qiu, 2018), and that throughout the viral infection the replicating genomes are present in both the ssDNA and dsDNA forms. The sequences we retrieved must therefore correspond to the cell-free dsDNA replication intermediates. This is coherent with the peculiar coverage pattern on the B19V genome, where the dsDNA hairpins at its terminal sites and are highly covered, reflecting a better stability of these regions over time (Figure 2b). Similarly, the partially circular dsDNA genome from HBV was poorly covered at the ssDNA region (Figure 2c), which also goes through a dsDNA phase during replication, a similar coverage is reported in three previous ancient HBV genomes (Krause-Kyora et al., 2018). Although HBV and B19V are also capable of integrating into the human host genome (Yuen et al., 2018; Janovitz et al., 2017), the uneven read coverage observed for all reconstructed viruses (higher coverage in dsDNA regions) suggests that these sequences do not correspond to integration events. If the B19V or HBV reads we recovered derived from integrated sequences in the human genome, we would expect an even coverage along the reference viral genome, which is not the case. Further analyses would be needed to determine if the aDNA retrieved in this and other studies comes from systemic circulating virions or from systemic cell-free DNA intermediates (Cheng et al., 2019) produced after viral replication in the bone marrow or liver for B19V and HBV, respectively (Broliden et al., 2006; Yuen et al., 2018).

The ancient B19V genomes were assigned to genotype 3. This genotype is most prevalent in West Africa (Ghana: 100%, n = 11; Burkina Faso: 100%, n = 5) (Candotti et al., 2004; Hübschen et al., 2009; Rinckel et al., 2009) and a potential African origin has been suggested (Candotti et al., 2004). It has also been sporadically found outside of Africa (Jain et al., 2015)^,(Candotti et al., 2004; Rinckel et al., 2009) in Brazil (50%, n = 12) (Freitas et al., 2008; Sanabani et al., 2006), India (15.4%, n = 13) (Jain et al., 2015), France (11.4%, n = 79) (Nguyen et al., 1999; Servant et al., 2002), and the USA (0.85%, n = 117) (Rinckel et al., 2009) as well as in immigrants from Morocco, Egypt, and Turkey in Germany (6.7%, n = 59) (Schneider et al., 2008). Two other genotypes, 1 and 2, exist for this virus. Genotype 1 is the most common and is found worldwide, while the almost extinct genotype 2 is mainly found in elderly people from Northern Europe (Pyöriä et al., 2017). Ancient genomes from genotypes 1 and 2 have been recovered from Eurasian samples, including a genotype 2 B19V genome from a 10th-century Viking burial in Greenland (Mühlemann et al., 2018b). ⁸⁷Sr/⁸⁶Sr isotopes on individuals from this burial revealed that they were immigrants from Iceland (Mühlemann et al., 2018b), suggesting an introduction of the genotype 2 to North America during Viking explorations of Greenland.

While serological evidence indicates that B19V currently circulates in Mexico, only the presence of genotype 1 has been formally described using molecular analyses (Valencia Pacheco et al., 2017). Taken together, our results are consistent with an introduction of the genotype 3 to New Spain as a consequence of the transatlantic slave trade imposed by European colonization. This hypothesis is supported by the ⁸⁷Sr/⁸⁶Sr isotopic analysis, which suggests that the individuals from the HSJN with B19V (HSJN240, HSJNC81) were born in West Africa and spent their last years of life in New Spain (Figure 4—figure supplement 1). Furthermore, the radiocarbon analysis for individuals HSJN240 and HSJN194 (Figure 4—figure supplement 2) support this notion as they correspond to the Early Colonial period, during which the number of enslaved Africans arriving was higher compared to later periods (Aguirre-Beltrán, 2005). Remarkably, a B19V genome belonging to the genotype 3 was recovered from an individual (COYC4) with 100% Indigenous ancestry (Figure 4b). COY4 was excavated in an independent archeological site 10 km south of the HSJN (Figure 1a), supporting the notion that viral transmissions between African individuals and Native Americans occurred during the Colonial period in Mexico City.

The HBV genotype A is highly diverse in Africa, reflecting its long evolutionary history, and likely originated somewhere between Africa, the Middle East, and Central Asia (Kostaki et al., 2018). The introduction of the genotype A from Africa to the Americas has been proposed based on phylogenetic analysis of modern strains from Brazil (Freitas et al., 2008; Kostaki et al., 2018) and Mexico (Roman et al., 2010), and more precisely of the sub-genotype A1 using sequences from Martinique (Brichler et al., 2013), Venezuela (Quintero et al., 2002), Haiti (Andernach et al., 2009), and Colombia (Alvarado-Mora et al., 2012). Recently, a similar introduction pattern was proposed for the quasi genotype A3 based on an ancient HBV genome recovered from an ancient African individual sampled in Mexico (Barquera et al., 2020). The origin of the sub-genotype A4 is controversial since the apparent African origin is based on modern sequences recovered from African immigrants living in Europe (Pourkarim et al., 2010). The Colonial ancient HBV genome reconstructed in our work represents the first ancient A4 linked to the transatlantic slave trade (Figure 3a, Figure 3—figure supplement 1), and the only report of this sub-genotype in the Americas, further supporting its African origin. The introduction of pathogens from Africa to the Americas has been proposed for other human-infecting viruses such as smallpox (Mandujano-Sánchez et al., 1982; Somolinos d’Árdois, 1982), based on historical records; or yellow fever virus (Bryant et al., 2007), HTLMV-1 (Gadelha et al., 2014), hepatitis C virus (genotype 2) (Markov et al., 2009), and human herpes simplex virus (Forni et al., 2020) based on phylogenetic analysis of modern strains from Afro-descendant or admixed human populations.

Although we cannot assert where exactly the African-born individuals in this study contracted B19V or HBV (Africa, America, or the Middle Passage) nor if the cause of their deaths can be attributed to such infections, the identification of ancient B19V and HBV in contexts associated with Colonial epidemics in Mexico City is still relevant in light of their paleopathological marks and the clinical information available for the closest sequences in the phylogenetic analyses. Notably, individual HSJNC81 displayed cribra orbitalia in the eye sockets and porotic hyperostosis on the cranial vault (Figure 4—figure supplement 3). The reconstructed ancient B19V genome from this individual is closest to the V9 strain, which was isolated from an infant with severe anemia and G6PD deficiency (Nguyen et al., 1999; GenBank: AJ249437; Figure 3b). The HSJN skeletal collection has a notably higher rate of cribra orbitalia and porotic hyperostosis compared to other Colonial archeological sites, marks that were proposed to be caused by an unknown infectious disease (Castillo-Chavez, 2000). These skeletal indicators are caused by irregular hematopoiesis in the bone marrow and are typically associated with genetic anemias such as thalassemia and sickle cell anemia (Angel, 1966), as well as to nutritional stress or parasitic infections (Walker et al., 2009). It is acknowledged that B19V infection can cause severe or even fatal anemia due to the low level of hemoglobin in individuals with other blood disorders, such as thalassemia, sickle cell anemia, malaria, or iron deficiency (Broliden et al., 2006; Heegaard and Brown, 2002). Therefore, since B19V infects precursors of the erythroid lineage (Broliden et al., 2006), it is possible that the morphological changes found in HSJNC81 might be the result of a severe anemia caused or enhanced by a B19V infection. With our data we cannot discard the simultaneous presence of a genetic disease since the loci for thalassemia, sickle cell anemia, and G6PD deficiency were not covered with our human-mapped NGS data. Nevertheless, the identification of ancient B19V in a Colonial context is noteworthy considering several recent reports reveal that measles-like cases were actually attributable to B19V (De Los Ángeles Ribas et al., 2019; Rezaei et al., 2016) or rubella (Anderson et al., 1985; Davidkin et al., 1998; De Los Ángeles Ribas et al., 2019; Rezaei et al., 2016), which produce a similar kind of rash and fever. Therefore, it is possible that B19V might have been responsible for some of the numerous cases attributed to measles that were described in early 16th-century Mexico (Acuña-Soto et al., 2004; Mandujano-Sánchez et al., 1982; Wesp, 2017), in particular historical records that document the treatment of an outbreak of measles at the HSJN in 1531 CE (Meza, 2013). Our study, however, does not reject the notorious role that measles played during the Colonial outbreaks (as it is strongly supported by historical records), but provides evidence of the presence of B19V during the Colonial period in Mexico City to facilitate discussions about the paradigmatic etiology of the supposed measles epidemics reported in historical records (Malvido and Viesca, 1982; Mandujano-Sánchez et al., 1982; Somolinos d’Árdois, 1982). This hypothesis requires additional comprehensive studies aimed at characterizing the presence of measles and rubella viruses from ancient remains, a task that currently poses difficult technical challenges given that RNA is known to degrade rapidly. In fact, most ancient viral RNA genomes have been recovered only from formalin-fixed tissue (Düx et al., 2020; Xiao et al., 2013).

Furthermore, historical records of the autopsies of the victims of the 1576 CE Cocoliztli epidemic treated at the HSJN describe the presence of enlarged hard liver and jaundice (Acuña-Soto et al., 2002; Acuña-Soto et al., 2004; Malvido and Viesca, 1982; Marr and Kiracofe, 2000; Somolinos d’Árdois, 1982) as well as a black spleen and lungs and heart with yellow liquid and black blood (Acuña-Soto et al., 2000; Malvido and Viesca, 1982; Somolinos d’Árdois, 1982). This is noteworthy given that both HBV and B19V viruses proliferate in the liver and are associated with hepatitis and jaundice (Broliden et al., 2006; Yuen et al., 2018). The radiocarbon dating of individuals HSJN194 (HBV) and HSJN240 (B19V) suggests that these individuals died between 1472–1625 CE and 1442–1608 CE (years calibrated for 1σ), respectively (Figure 4—figure supplement 2), which overlaps with the period of time when the hepatitis symptoms were reported in the autopsies after the 1576 Cocoliztli epidemic at the HSJN (Acuña-Soto et al., 2004; Marr and Kiracofe, 2000; Somolinos d’Árdois, 1982). However, additional analyses are needed before being able to establish a link between these viruses and the wide array of symptoms described for Cocoliztli. Currently, technological limitations prevent the direct identification of ancient RNA viruses in bone or dental remains. However, future studies, with larger sample sizes from different contexts associated with the outbreak, should explore a wider range of pathogens previously suggested as potential causative agents, like arthropod-borne pathogens (malaria, yellow fever virus, and dengue virus) (Marr and Kiracofe, 2000) or hemorrhagic fever RNA viruses (Acuña-Soto et al., 2004).

Furthermore, it is important to acknowledge that both viruses have also been previously identified in aDNA datasets not necessarily associated with disease or epidemic contexts (Kahila Bar-Gal et al., 2012; Krause-Kyora et al., 2018; Mühlemann et al., 2018a; Patterson Ross et al., 2018). Additionally, our data is not sufficient to elucidate the age when the individuals acquired the viruses or if it is related to their cause of death.

In the case of HSJN194, we cannot establish if he acquired HBV vertically or horizontally, nor if this individual presented an acute or chronic infection. Finally, although our data does not allow us to associate these viruses to a specific epidemic outbreak, the identification of HBV and B19V in Post-Contact remains opens up new opportunities for investigating their presence in similar contexts and expand our knowledge on their evolution and potential link to disease in Colonial Mexico. This type of research is particularly relevant when considering previous hypotheses favoring the synergistic action of different types of pathogens in these devastating Colonial epidemics (Somolinos d’Árdois, 1982).

It is important to emphasize that our findings should be interpreted with careful consideration of the historical and social context of the transatlantic slave trade. This cruel episode in history involved the forced displacement of millions of individuals to the Americas (ca. 250,000 to New Spain; Aguirre-Beltrán, 2005) under inhumane, unsanitary, and overcrowded conditions that, with no doubt, favored the spread of infectious diseases (Mandujano-Sánchez et al., 1982). Therefore, the introduction of these and other pathogens from Africa to the Americas should be attributed to the brutal and harsh conditions of the Middle passage that enslaved Africans were subjected to by traders and colonizers, and not to the African peoples themselves. Moreover, the adverse life conditions for enslaved Africans and Native Americans, especially during the first decades after colonization, surely favored the spread of diseases and emergence of epidemics (Mandujano-Sánchez et al., 1982). Integrative and multidisciplinary approaches are thus needed to understand this phenomenon in full.

In summary, our study provides direct aDNA evidence of HBV and B19V introduced to the Americas from Africa during the transatlantic slave trade. The isolation and characterization of these ancient HBV and B19V genomes represent an important contribution to the ancient viral genomes reported in the Americas (Barquera et al., 2020; Duggan et al., 2020; Schrick et al., 2017). Our results expand our knowledge on the viral agents that were in circulation during Colonial epidemics like Cocoliztli, some of which resulted in the catastrophic collapse of the immunologically naïve Indigenous population. Although we cannot assign a direct causality link between HBV and B19V and Cocoliztli, our findings confirm that these potentially harmful viruses were indeed circulating in individuals found in archeological contexts associated with this epidemic outbreak. Further analyses from different sites and samples will help understand the possible role of these and other pathogens in Colonial epidemics, as well as the full spectrum of pathogens that were introduced to the Americas during European colonization.

Reconstructed genomes from this study are available in Genbank under accession number MT108214, MT108215, MT108216, MT108217. NGS reads used to reconstruct ancient viral genomes reported in this study are available in Dryad (https://doi.org/10.5061/dryad.5x69p8d2s).

1. 1. Avila-Acros MC

   (2021) Dryad Digital Repository

   Data from: Ancient viral genomes reveal introduction of human pathogenic viruses into Mexico during the transatlantic slave trade.

   https://doi.org/10.5061/dryad.5x69p8d2s

1. 1. 1000 Genomes Project Consortium

   (2015) A global reference for human genetic variation

   Nature 526:68–74.

   https://doi.org/10.1038/nature15393

   • PubMed
   • Google Scholar
2. 1. Acuña-Soto R
   2. Romero LC
   3. Maguire JH

   (2000) Large epidemics of hemorrhagic fevers in Mexico 1545-1815

   The American Journal of Tropical Medicine and Hygiene 62:733–739.

   https://doi.org/10.4269/ajtmh.2000.62.733

   • PubMed
   • Google Scholar
3. 1. Acuña-Soto R
   2. Stahle DW
   3. Cleaveland MK
   4. Therrell MD

   (2002) Megadrought and Megadeath in 16th Century Mexico

   Emerging Infectious Diseases 8:360–362.

   https://doi.org/10.3201/eid0804.010175

   • PubMed
   • Google Scholar
4. 1. Acuña-Soto R
   2. Stahle DW
   3. Therrell MD
   4. Griffin RD
   5. Cleaveland MK

   (2004) When half of the population died: The epidemic of hemorrhagic fevers of 1576 in Mexico

   FEMS Microbiology Letters 240:1–5.

   https://doi.org/10.1016/j.femsle.2004.09.011

   • PubMed
   • Google Scholar
5. 1. Aguirre-Beltrán G

   (2005)

   La presencia del negro en México

   Revista Del CESLA 7:351–367.

   • Google Scholar
6. 1. Alexander DH
   2. Lange K

   (2011) Enhancements to the ADMIXTURE algorithm for individual ancestry estimation

   BMC Bioinformatics 12:246.

   https://doi.org/10.1186/1471-2105-12-246

   • Google Scholar
7. 1. Altschul SF
   2. Gish W
   3. Miller W
   4. Myers EW
   5. Lipman DJ

   (1990) Basic local alignment search tool

   Journal of Molecular Biology 215:403–410.

   https://doi.org/10.1016/S0022-2836(05)80360-2

   • PubMed
   • Google Scholar
8. 1. Alvarado-Mora M
   2. Romano CM
   3. Gomes-Gouvêa MS
   4. Gutierrez MF
   5. Carrilho FJ
   6. Pinho JRR

   (2012) Phylogenetic analysis of complete genome sequences of hepatitis B virus from an Afro-Colombian community: Presence of HBV F3/A1 recombinant strain

   Virology Journal 9:244.

   https://doi.org/10.1186/1743-422X-9-244

   • PubMed
   • Google Scholar
9. 1. Andernach IE
   2. Nolte C
   3. Pape JW
   4. Muller CP

   (2009) Slave trade and hepatitis B virus genotypes and subgenotypes in Haiti and Africa

   Emerging Infectious Diseases 15:1222–1228.

   https://doi.org/10.3201/eid1508.081642

   • PubMed
   • Google Scholar
10. 1. Anderson MJ
    2. Kidd IM
    3. Morgan-Capner P

    (1985) Human parvovirus and rubella-like illness

    The Lancet 326:663.

    https://doi.org/10.1016/S0140-6736(85)90025-X

    • Google Scholar
11. 1. Angel JL

    (1966) Porotic hyperostosis, anemias, malarias, and marshes in the prehistoric Eastern Mediterranean

    Science 153:760–763.

    https://doi.org/10.1126/science.153.3737.760

    • PubMed
    • Google Scholar
12. 1. Barquera R
    2. Lamnidis TC
    3. Lankapalli AK
    4. Kocher A
    5. Hernández-Zaragoza DI
    6. Nelson EA
    7. Zamora-Herrera AC
    8. Ramallo P
    9. Bernal-Felipe N
    10. Immel A
    11. Bos K
    12. Acuña-Alonzo V
    13. Barbieri C
    14. Roberts P
    15. Herbig A
    16. Kühnert D
    17. Márquez-Morfín L
    18. Krause J

    (2020) Origin and health status of first-generation Africans from early colonial Mexico

    Current Biology 30:2078–2091.

    https://doi.org/10.1016/j.cub.2020.04.002

    • PubMed
    • Google Scholar
13. 1. Bastos MQR
    2. de Souza SMF
    3. Rodrigues-Carvalho C
    4. Tykot RH
    5. Cook DC
    6. Santos R

    (2016) Isotopic study of geographic origins and diet of enslaved Africans buried in two Brazilian cemeteries

    Journal of Archaeological Science 70:82–90.

    https://doi.org/10.1016/j.jas.2016.04.020

    • Google Scholar
14. 1. Bernard-Griffiths J
    2. Peucat JJ
    3. Fourcade S
    4. Kienast JR
    5. Ouzegane K

    (1988) Origin and evolution of 2 Ga old carbonatite complex (Ihouhaouene, Ahaggar, Algeria): Nd and Sr isotopic evidence

    Contributions to Mineralogy and Petrology 100:339–348.

    https://doi.org/10.1007/BF00379743

    • Google Scholar
15. 1. Bravo-Lopez M
    2. Villa-Islas V
    3. Rocha Arriaga C
    4. Villaseñor-Altamirano AB
    5. Guzmán-Solís A
    6. Sandoval-Velasco M
    7. Wesp JK
    8. Alcantara K
    9. López-Corral A
    10. Gómez-Valdés J
    11. Mejía E
    12. Herrera A
    13. Meraz-Moreno A
    14. Moreno-Cabrera M de la L
    15. Moreno-Estrada A
    16. Nieves-Colón MA
    17. Olvera J
    18. Pérez-Pérez J
    19. Iversen KH
    20. Rasmussen S
    21. Sandoval K
    22. Zepeda G
    23. Ávila-Arcos MC

    (2020) Paleogenomic insights into the red complex bacteria tannerella forsythia in pre-hispanic and colonial individuals from Mexico

    Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences 375:20190580.

    https://doi.org/10.1098/rstb.2019.0580

    • PubMed
    • Google Scholar
16. 1. Breitwieser FP
    2. Salzberg SL

    (2020) Pavian: Interactive analysis of metagenomics data for microbiome studies and pathogen identification

    Bioinformatics 36:1303–1304.

    https://doi.org/10.1093/bioinformatics/btz715

    • PubMed
    • Google Scholar
17. 1. Brichler S
    2. Lagathu G
    3. Chekaraou MA
    4. Le Gal F
    5. Edouard A
    6. Dény P
    7. Césaire R
    8. Gordien E

    (2013) African, Amerindian and European hepatitis B virus strains circulate on the Caribbean Island of Martinique

    The Journal of General Virology 94:2318–2329.

    https://doi.org/10.1099/vir.0.055459-0

    • PubMed
    • Google Scholar
18. 1. Briggs AW
    2. Stenzel U
    3. Johnson PLF
    4. Green RE
    5. Kelso J
    6. Prüfer K
    7. Meyer M
    8. Krause J
    9. Ronan MT
    10. Lachmann M
    11. Pääbo S

    (2007) Patterns of damage in genomic DNA sequences from a Neandertal

    PNAS 104:14616–14621.

    https://doi.org/10.1073/pnas.0704665104

    • PubMed
    • Google Scholar
19. 1. Broliden K
    2. Tolfvenstam T
    3. Norbeck O

    (2006) Clinical aspects of parvovirus B19 infection

    Journal of Internal Medicine 260:285–304.

    https://doi.org/10.1111/j.1365-2796.2006.01697.x

    • PubMed
    • Google Scholar
20. 1. Bronk Ramsey C

    (2013) Recent and Planned Developments of the Program OxCal

    Radiocarbon 55:3–4.

    https://doi.org/10.2458/azu_js_rc.55.16215

    • Google Scholar
21. 1. Bryant JE
    2. Holmes EC
    3. Barrett ADT

    (2007) Out of Africa: A Molecular Perspective on the Introduction of Yellow Fever Virus into the Americas

    PLOS Pathogens 3:e0030075.

    https://doi.org/10.1371/journal.ppat.0030075

    • Google Scholar
22. Book

    1. Cabrera-Torres JJ
    2. García-Martínez M

    (1998)

    Utilización, Modificación y Reuso de Los Espacios Del Edificio Sede Del Hospital Real de San José de Los Naturales

    Escuela Nacional de Antropología e Historia.

    • Google Scholar
23. 1. Candotti D
    2. Etiz N
    3. Parsyan A
    4. Allain JP

    (2004) Identification and Characterization of Persistent Human Erythrovirus Infection in Blood Donor Samples

    Journal of Virology 78:12169–12178.

    https://doi.org/10.1128/JVI.78.22.12169-12178.2004

    • PubMed
    • Google Scholar
24. Book

    1. Castillo-Chavez O

    (2000)

    Condiciones de Vida y Salud de Una Muestra Poblacional En La Ciudad de México En La Época Colonia

    Escuela Nacional de Antropología e Historia.

    • Google Scholar
25. 1. Cheng AP
    2. Burnham P
    3. Lee JR
    4. Cheng MP
    5. Suthanthiran M
    6. Dadhania D
    7. Vlaminck D

    (2019) A cell-free DNA metagenomic sequencing assay that integrates the host injury response to infection

    PNAS 116:18738–18744.

    https://doi.org/10.1073/pnas.1906320116

    • Google Scholar
26. Software

    1. Dabdoub S
    2. Connell J
    3. Clarke E
    4. Borry M

    (2018) kraken-biom, version 46b32df

    GitHub.

    https://github.com/smdabdoub/kraken-biom
27. 1. Dabney J
    2. Knapp M
    3. Glocke I
    4. Gansauge MT
    5. Weihmann A
    6. Nickel B
    7. Valdiosera C
    8. García N
    9. Pääbo S
    10. Arsuaga JL
    11. Meyer M

    (2013) Complete mitochondrial genome sequence of a middle pleistocene cave bear reconstructed from ultrashort dna fragments

    PNAS 110:15758–15763.

    https://doi.org/10.1073/pnas.1314445110

    • PubMed
    • Google Scholar
28. 1. Darriba D
    2. Taboada GL
    3. Doallo R
    4. Posada D

    (2012) jModelTest 2: More models, new heuristics and high-performance computing

    Nature Methods 9:772.

    https://doi.org/10.1038/nmeth.2109

    • Google Scholar
29. 1. Davidkin I
    2. Valle M
    3. Peltola H
    4. Hovi T
    5. Paunio M
    6. Roivainen M
    7. Linnavuori K
    8. Jokinen S
    9. Leinikki P

    (1998) Etiology of Measles‐ and Rubella‐like Illnesses in Measles, Mumps, and Rubella–Vaccinated Children

    The Journal of Infectious Diseases 178:1567–1570.

    https://doi.org/10.1086/314513

    • PubMed
    • Google Scholar
30. 1. De Los Ángeles Ribas M
    2. Tejero Y
    3. Cordero Y
    4. Pérez D
    5. Sausy A
    6. Muller CP
    7. Hübschen JM

    (2019) Identification of human parvovirus B19 among measles and rubella suspected patients from Cuba

    Journal of Medical Virology 91:1351–1354.

    https://doi.org/10.1002/jmv.25444

    • PubMed
    • Google Scholar
31. 1. Drexler JF
    2. Geipel A
    3. Konig A
    4. Corman VM
    5. van Riel D
    6. Leijten LM
    7. Bremer CM
    8. Rasche A
    9. Cottontail VM
    10. Maganga GD
    11. Schlegel M
    12. Muller MA
    13. Adam A
    14. Klose SM
    15. Borges Carneiro AJ
    16. Stocker A
    17. Franke CR
    18. Gloza-Rausch F
    19. Geyer J
    20. Drosten C

    (2013) Bats carry pathogenic hepadnaviruses antigenically related to hepatitis B virus and capable of infecting human hepatocytes

    PNAS 110:16151–16156.

    https://doi.org/10.1073/pnas.1308049110

    • Google Scholar
32. 1. Drummond AJ
    2. Suchard MA
    3. Xie D
    4. Rambaut A

    (2012) Bayesian Phylogenetics with BEAUti and the BEAST 1.7

    Molecular Biology and Evolution 29:1969–1973.

    https://doi.org/10.1093/molbev/mss075

    • PubMed
    • Google Scholar
33. 1. Duggan AT
    2. Perdomo MF
    3. Piombino-Mascali D
    4. Marciniak S
    5. Poinar D
    6. Emery M
    7. Buchmann JP
    8. Duchêne S
    9. Jankauskas R
    10. Humphreys M
    11. Golding GB
    12. Southon J
    13. Devault A
    14. Rouillard JM
    15. Sahl JW
    16. Dutour O
    17. Hedman K
    18. Sajantila A
    19. Smith GL
    20. Poinar HN

    (2016) 17th Century Variola Virus Reveals the Recent History of Smallpox

    Current Biology: CB 26:3407–3412.

    https://doi.org/10.1016/j.cub.2016.10.061

    • PubMed
    • Google Scholar
34. 1. Duggan AT
    2. Klunk J
    3. Porter AF
    4. Dhody AN
    5. Hicks R
    6. Smith GL
    7. Humphreys M
    8. McCollum AM
    9. Davidson WB
    10. Wilkins K
    11. Li Y
    12. Burke A
    13. Polasky H
    14. Flanders L
    15. Poinar D
    16. Raphenya AR
    17. Lau TTY
    18. Alcock B
    19. McArthur AG
    20. Golding GB
    21. Holmes EC
    22. Poinar HN

    (2020) The origins and genomic diversity of American Civil War Era smallpox vaccine strains

    Genome Biology 21:175.

    https://doi.org/10.1186/s13059-020-02079-z

    • PubMed
    • Google Scholar
35. 1. Düx A
    2. Lequime S
    3. Patrono L
    4. Vrancken B
    5. Boral S
    6. Gogarten JF
    7. Hilbig A
    8. Horst D
    9. Merkel K
    10. Prepoint B
    11. Santibanez S
    12. Schlotterbeck J
    13. Suchard MA
    14. Ulrich M
    15. Widulin N
    16. Mankertz A
    17. Leendertz FH
    18. Harper K
    19. Schnalke T
    20. Calvignac-Spencer S

    (2020) Measles virus and rinderpest virus divergence dated to the sixth century BCE

    Science 368:1367–1370.

    https://doi.org/10.1126/science.aba9411

    • PubMed
    • Google Scholar
36. 1. Edgar RC

    (2004) MUSCLE: Multiple sequence alignment with high accuracy and high throughput

    Nucleic Acids Research 32:1792–1797.

    https://doi.org/10.1093/nar/gkh340

    • PubMed
    • Google Scholar
37. 1. Fairbairn HW
    2. Moorbath S
    3. Ramo AO
    4. Pinson WH
    5. Hurley PM

    (1967) Rb-Sr age of granitic rocks of southeastern massachusetts and the age of the lower cambrian at Hoppin Hill

    Earth and Planetary Science Letters 2:321–328.

    https://doi.org/10.1016/0012-821X(67)90149-5

    • Google Scholar
38. 1. Feng Q
    2. Lu D
    3. Xu S

    (2018) AncestryPainter: A Graphic Program for Displaying Ancestry Composition of Populations and Individuals

    Genomics, Proteomics & Bioinformatics 16:382–385.

    https://doi.org/10.1016/j.gpb.2018.05.002

    • Google Scholar
39. 1. Ferrière L
    2. Koeberl C
    3. Thöni M
    4. Liang C

    (2010) Single crystal U–Pb zircon age and Sr–Nd isotopic composition of impactites from the Bosumtwi impact structure, Ghana: Comparison with country rocks and Ivory Coast tektites

    Chemical Geology 275:254–261.

    https://doi.org/10.1016/j.chemgeo.2010.05.016

    • PubMed
    • Google Scholar
40. 1. Firth C
    2. Kitchen A
    3. Shapiro B
    4. Suchard MA
    5. Holmes EC
    6. Rambaut A

    (2010) Using time-structured data to estimate evolutionary rates of double-stranded DNA viruses

    Molecular Biology and Evolution 27:2038–2051.

    https://doi.org/10.1093/molbev/msq088

    • PubMed
    • Google Scholar
41. 1. Forni D
    2. Pontremoli C
    3. Clerici M
    4. Pozzoli U
    5. Cagliani R
    6. Sironi M

    (2020) Recent out-of-africa migration of human Herpes simplex viruses

    Molecular Biology and Evolution 37:1259–1271.

    https://doi.org/10.1093/molbev/msaa001

    • PubMed
    • Google Scholar
42. 1. Freitas RB
    2. Melo FL
    3. Oliveira DS
    4. Romano CM
    5. Freitas MRC
    6. Macêdo O
    7. Linhares AC
    8. Durigon EL

    (2008) Molecular characterization of human erythrovirus B19 strains obtained from patients with several clinical presentations in the Amazon region of Brazil

    Journal of Clinical Virology: The Official Publication of the Pan American Society for Clinical Virology 43:60–65.

    https://doi.org/10.1016/j.jcv.2008.03.033

    • PubMed
    • Google Scholar
43. 1. Fricke F
    2. Laffoon J
    3. Victorina A
    4. Haviser J

    (2020) Delayed physical development in a first generation enslaved African woman from Pietermaai, Curaçao

    International Journal of Osteoarchaeology 30:43–52.

    https://doi.org/10.1002/oa.2829

    • Google Scholar
44. 1. Fullgraf T
    2. Ndiaye PM
    3. Blein O
    4. Buscail F
    5. Lahondère D
    6. Métour JL
    7. Sergeev S
    8. Tegyey M

    (2013) Silurian magmatism in eastern Senegal and its significance for the Paleozoic evolution of NW-Gondwana

    Journal of African Earth Sciences 78:66–85.

    https://doi.org/10.1016/j.jafrearsci.2012.08.003

    • Google Scholar
45. 1. Gadelha SR
    2. Aleluia MM
    3. Mello MA
    4. Rego FF
    5. Pereira LS
    6. Galvão-Castro B
    7. Gonçalves M
    8. Sousa SM
    9. Alcântara LC

    (2014) The origin of HTLV-1 in the South Bahia by phylogenetic, mitochondrial DNA and β-globin analysis

    Retrovirology 11:P49.

    https://doi.org/10.1186/1742-4690-11-S1-P49

    • Google Scholar
46. 1. Ganaie SS
    2. Qiu J

    (2018) Recent Advances in Replication and Infection of Human Parvovirus B19

    Frontiers in Cellular and Infection Microbiology 8:166.

    https://doi.org/10.3389/fcimb.2018.00166

    • Google Scholar
47. 1. Heegaard ED
    2. Brown KE

    (2002) Human Parvovirus B19

    Clinical Microbiology Reviews 15:485–505.

    https://doi.org/10.1128/CMR.15.3.485-505.2002

    • PubMed
    • Google Scholar
48. 1. Herbig A
    2. Campana MG
    3. García NMR
    4. Warinner C
    5. Sabin S
    6. Spyrou MA
    7. Valtueña AA
    8. Huson D
    9. Tuross N
    10. Bos K
    11. Krause J

    (2018) Salmonella enterica genomes from victims of a major sixteenth-century epidemic in Mexico

    Nature Ecology & Evolution 2:520.

    https://doi.org/10.1038/s41559-017-0446-6

    • Google Scholar
49. Book

    1. Hernández-Lopez PE
    2. Negrete S

    (2012)

    Realmente Eran Indios? Afinidad Biológica Entre Las Personas Atendidas En El Hospital Real San Jose de Los Naturales, Siglos XVI - XVIII

    Escuela Nacional de Antropología e Historia.

    • Google Scholar
50. 1. Hübschen JM
    2. Mihneva Z
    3. Mentis AF
    4. Schneider F
    5. Aboudy Y
    6. Grossman Z
    7. Rudich H
    8. Kasymbekova K
    9. Sarv I
    10. Nedeljkovic J
    11. Tahita MC
    12. Tarnagda Z
    13. Ouedraogo JB
    14. Gerasimova AG
    15. Moskaleva TN
    16. Tikhonova NT
    17. Chitadze N
    18. Forbi JC
    19. Faneye AO
    20. Muller CP

    (2009) Phylogenetic analysis of human parvovirus b19 sequences from eleven different countries confirms the predominance of genotype 1 and suggests the spread of genotype 3b

    Journal of Clinical Microbiology 47:3735–3738.

    https://doi.org/10.1128/JCM.01201-09

    • PubMed
    • Google Scholar
51. 1. Huson DH
    2. Beier S
    3. Flade I
    4. Górska A
    5. El-Hadidi M
    6. Mitra S
    7. Ruscheweyh HJ
    8. Tappu R

    (2016) MEGAN Community Edition—Interactive Exploration and Analysis of Large-Scale Microbiome Sequencing Data

    PLOS Computational Biology 12:e1004957.

    https://doi.org/10.1371/journal.pcbi.1004957

    • Google Scholar
52. 1. Jain P
    2. Jain A
    3. Prakash S
    4. Khan DN
    5. Singh DD
    6. Kumar A
    7. Moulik NR
    8. Chandra T

    (2015) Prevalence and genotypic characterization of human parvovirus B19 in children with hemato-oncological disorders in North India

    Journal of Medical Virology 87:303–309.

    https://doi.org/10.1002/jmv.24028

    • PubMed
    • Google Scholar
53. 1. Janovitz T
    2. Wong S
    3. Young NS
    4. Oliveira T
    5. Falck-Pedersen E

    (2017) Parvovirus B19 integration into human CD36+ erythroid progenitor cells

    Virology 511:40–48.

    https://doi.org/10.1016/j.virol.2017.08.011

    • PubMed
    • Google Scholar
54. 1. Jónsson H
    2. Ginolhac A
    3. Schubert M
    4. Johnson PLF
    5. Orlando L

    (2013) mapDamage2.0: Fast approximate Bayesian estimates of ancient DNA damage parameters

    Bioinformatics (Oxford, England) 29:1682–1684.

    https://doi.org/10.1093/bioinformatics/btt193

    • PubMed
    • Google Scholar
55. 1. Kahila Bar-Gal G
    2. Kim MJ
    3. Klein A
    4. Shin DH
    5. Oh CS
    6. Kim JW
    7. Kim TH
    8. Kim SB
    9. Grant PR
    10. Pappo O
    11. Spigelman M
    12. Shouval D

    (2012) Tracing hepatitis B virus to the 16th century in a Korean mummy

    Hepatology 56:1671–1680.

    https://doi.org/10.1002/hep.25852

    • PubMed
    • Google Scholar
56. Book

    1. Karam-Tapia CE

    (2012)

    Estimación Del Mestizaje Mediante La Morfología Dental En La Ciudad de México (Siglo XVI al XIX

    Escuela Nacional de Antropología e Historia.

    • Google Scholar
57. 1. Kearse M
    2. Moir R
    3. Wilson A
    4. Stones-Havas S
    5. Cheung M
    6. Sturrock S
    7. Buxton S
    8. Cooper A
    9. Markowitz S
    10. Duran C
    11. Thierer T
    12. Ashton B
    13. Meintjes P
    14. Drummond A

    (2012) Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data

    Bioinformatics 28:1647–1649.

    https://doi.org/10.1093/bioinformatics/bts199

    • PubMed
    • Google Scholar
58. 1. Key FM
    2. Posth C
    3. Krause J
    4. Herbig A
    5. Bos K

    (2017) Mining Metagenomic Data Sets for Ancient DNA: Recommended Protocols for Authentication

    Trends in Genetics: TIG 33:508–520.

    https://doi.org/10.1016/j.tig.2017.05.005

    • PubMed
    • Google Scholar
59. 1. Kloss-Brandstätter A
    2. Pacher D
    3. Schönherr S
    4. Weissensteiner H
    5. Binna R
    6. Specht G
    7. Kronenberg F

    (2011) HaploGrep: A fast and reliable algorithm for automatic classification of mitochondrial DNA haplogroups

    Human Mutation 32:25–32.

    https://doi.org/10.1002/humu.21382

    • PubMed
    • Google Scholar
60. 1. Kostaki EG
    2. Karamitros T
    3. Stefanou G
    4. Mamais I
    5. Angelis K
    6. Hatzakis A
    7. Kramvis A
    8. Paraskevis D

    (2018) Unravelling the history of hepatitis B virus genotypes A and D infection using a full-genome phylogenetic and phylogeographic approach

    eLife 7:e36709.

    https://doi.org/10.7554/eLife.36709

    • Google Scholar
61. 1. Kramvis A

    (2014) Genotypes and Genetic Variability of Hepatitis B Virus

    Intervirology 57:141–150.

    https://doi.org/10.1159/000360947

    • PubMed
    • Google Scholar
62. 1. Krause-Kyora B
    2. Susat J
    3. Key FM
    4. Kühnert D
    5. Bosse E
    6. Immel A
    7. Rinne C
    8. Kornell SC
    9. Yepes D
    10. Franzenburg S
    11. Heyne HO
    12. Meier T
    13. Lösch S
    14. Meller H
    15. Friederich S
    16. Nicklisch N
    17. Alt KW
    18. Schreiber S
    19. Tholey A
    20. Krause J

    (2018) Neolithic and medieval virus genomes reveal complex evolution of hepatitis B

    eLife 7:e36666.

    https://doi.org/10.7554/eLife.36666

    • PubMed
    • Google Scholar
63. 1. Larsson A

    (2014) AliView: A fast and lightweight alignment viewer and editor for large datasets

    Bioinformatics 30:3276–3278.

    https://doi.org/10.1093/bioinformatics/btu531

    • PubMed
    • Google Scholar
64. 1. Li H
    2. Durbin R

    (2009) Fast and accurate short read alignment with Burrows–Wheeler transform

    Bioinformatics 25:1754–1760.

    https://doi.org/10.1093/bioinformatics/btp324

    • PubMed
    • Google Scholar
65. 1. Li H
    2. Handsaker B
    3. Wysoker A
    4. Fennell T
    5. Ruan J
    6. Homer N
    7. Marth G
    8. Abecasis G
    9. Durbin R
    10. 1000 Genome Project Data Processing Subgroup

    (2009) The sequence Alignment/map format and samtools

    Bioinformatics 25:2078–2079.

    https://doi.org/10.1093/bioinformatics/btp352

    • PubMed
    • Google Scholar
66. 1. Liégeois JP
    2. Sauvage JF
    3. Black R

    (1991) The Permo-Jurassic alkaline province of Tadhak, Mali: Geology, geochronology and tectonic significance

    Lithos 27:95–105.

    https://doi.org/10.1016/0024-4937(91)90022-D

    • Google Scholar
67. 1. Lindahl T

    (1993) Instability and decay of the primary structure of DNA

    Nature 362:709–715.

    https://doi.org/10.1038/362709a0

    • PubMed
    • Google Scholar
68. 1. Lindo J
    2. Huerta-Sánchez E
    3. Nakagome S
    4. Rasmussen M
    5. Petzelt B
    6. Mitchell J
    7. Cybulski JS
    8. Willerslev E
    9. DeGiorgio M
    10. Malhi RS

    (2016) A time transect of exomes from a Native American population before and after European contact

    Nature Communications 7:13175.

    https://doi.org/10.1038/ncomms13175

    • PubMed
    • Google Scholar
69. 1. Luo Y
    2. Qiu J

    (2015) Human parvovirus B19: A mechanistic overview of infection and DNA replication

    Future Virology 10:155–167.

    https://doi.org/10.2217/fvl.14.103

    • PubMed
    • Google Scholar
70. Book

    1. Malvido E
    2. Viesca C

    (1982)

    La epidemia de cocoliztli de 1576

    In: Florescano E, Malvido E, editors. Ensayos Sobre La Historia de Las Epidemias En México. Instituto Mexicano del Seguro Social. pp. 27–32.

    • Google Scholar
71. Book

    1. Mandujano-Sánchez A
    2. Solache LC
    3. Mandujano MA

    (1982)

    Historia de las Epidemias en el México Antiguo: Algunos Aspectos Biológicos y Sociales

    In: Florescano E, editors. (EM Id= Id="ba2a0096-2ce4-4952-A88b-Ea42e65865c5">’’c95e7f0d-715a-46b5-B60d-476e6d922acb’’)(EmID=’’ba2a00(/Em)EEnssobre-A88laeHistDe5c5’’)IcLasepidemias. Instituto Mexicano del Seguro Social. pp. 2086–2096.

    • Google Scholar
72. 1. Markov P
    2. Pepin J
    3. Frost E
    4. Deslandes S
    5. Labbé AC
    6. Pybus OG

    (2009) Phylogeography and molecular epidemiology of hepatitis C virus genotype 2 in Africa

    The Journal of General Virology 90:2086–2096.

    https://doi.org/10.1099/vir.0.011569-0

    • PubMed
    • Google Scholar
73. 1. Marr JS
    2. Kiracofe JB

    (2000)

    Was the huey cocoliztli a haemorrhagic fever

    Medical History 44:341–362.

    • Google Scholar
74. 1. Meyer M
    2. Kircher M

    (2010) Illumina Sequencing Library Preparation for Highly Multiplexed Target Capture and Sequencing

    Cold Spring Harbor Protocols 2010:1–10.

    https://doi.org/10.1101/pdb.prot5448

    • Google Scholar
75. 1. Meza A

    (2013)

    Presencia africana en el cementerio del Hospital Real de San Jóse de los Naturales

    Arqueología Mexicana 119:40–44.

    • Google Scholar
76. Book

    1. Moreno-Cabrera M
    2. Meraz-Moreno A
    3. Cervantes-Rosado J

    (2015)

    Relleno aligerado con vasijas cerámicas en el templo de la Inmaculada Concepción, en Coyoacán

    In: Moreno-Cabrera M, editors. Boletín de Monumentos Históricos. Instituto Nacional de Antropología e Historia. pp. 121–134.

    • Google Scholar
77. 1. Mühlemann B
    2. Jones TC
    3. Damgaard P
    4. Allentoft ME
    5. Shevnina I
    6. Logvin A
    7. Usmanova E
    8. Panyushkina IP
    9. Boldgiv B
    10. Bazartseren T
    11. Tashbaeva K
    12. Merz V
    13. Lau N
    14. Smrčka V
    15. Voyakin D
    16. Kitov E
    17. Epimakhov A
    18. Pokutta D
    19. Vicze M
    20. Willerslev E

    (2018a) Ancient hepatitis B viruses from the Bronze Age to the Medieval period

    Nature 557:418–423.

    https://doi.org/10.1038/s41586-018-0097-z

    • PubMed
    • Google Scholar
78. 1. Mühlemann B
    2. Margaryan A
    3. Damgaard P
    4. Allentoft ME
    5. Vinner L
    6. Hansen AJ
    7. Weber A
    8. Bazaliiskii VI
    9. Molak M
    10. Arneborg J
    11. Bogdanowicz W
    12. Falys C
    13. Sablin M
    14. Smrčka V
    15. Sten S
    16. Tashbaeva K
    17. Lynnerup N
    18. Sikora M
    19. Smith DJ
    20. Fouchier RAM
    21. Drosten C
    22. Sjögren K-G
    23. Kristiansen K
    24. Willerslev E
    25. Jones TC

    (2018b) Ancient human parvovirus b19 in Eurasia reveals its long-term association with humans

    PNAS 115:7557–7562.

    https://doi.org/10.1073/pnas.1804921115

    • PubMed
    • Google Scholar
79. 1. Mühlemann B
    2. Vinner L
    3. Margaryan A
    4. Wilhelmson H
    5. de la Fuente Castro C
    6. Allentoft ME
    7. de Barros Damgaard P
    8. Hansen AJ
    9. Holtsmark Nielsen S
    10. Strand LM
    11. Bill J
    12. Buzhilova A
    13. Pushkina T
    14. Falys C
    15. Khartanovich V
    16. Moiseyev V
    17. Jørkov MLS
    18. Østergaard Sørensen P
    19. Magnusson Y
    20. Sikora M

    (2020) Diverse variola virus (smallpox) strains were widespread in northern Europe in the Viking Age

    Science 369:6502.

    https://doi.org/10.1126/science.aaw8977

    • PubMed
    • Google Scholar
80. 1. Neukamm J
    2. Pfrengle S
    3. Molak M
    4. Seitz A
    5. Francken M
    6. Eppenberger P
    7. Avanzi C
    8. Reiter E
    9. Urban C
    10. Welte B
    11. Stockhammer PW
    12. Teßmann B
    13. Herbig A
    14. Harvati K
    15. Nieselt K
    16. Krause J
    17. Schuenemann VJ

    (2020) 2000-year-old pathogen genomes reconstructed from metagenomic analysis of egyptian mummified individuals

    BMC Biology 18:108.

    https://doi.org/10.1186/s12915-020-00839-8

    • PubMed
    • Google Scholar
81. 1. Nguyen QT
    2. Sifer C
    3. Schneider V
    4. Allaume X
    5. Servant A
    6. Bernaudin F
    7. Auguste V
    8. Garbarg-Chenon A

    (1999) Novel Human Erythrovirus Associated with Transient Aplastic Anemia

    Journal of Clinical Microbiology 37:2483–2487.

    https://doi.org/10.1128/JCM.37.8.2483-2487.1999

    • PubMed
    • Google Scholar
82. 1. Ondov BD
    2. Bergman NH
    3. Phillippy AM

    (2011) Interactive metagenomic visualization in a Web browser

    BMC Bioinformatics 12:385.

    https://doi.org/10.1186/1471-2105-12-385

    • Google Scholar
83. 1. Pajer P
    2. Dresler J
    3. Kabíckova H
    4. Písa L
    5. Aganov P
    6. Fucik K
    7. Elleder D
    8. Hron T
    9. Kuzelka V
    10. Velemínsky P
    11. Klimentova J
    12. Fucikova A
    13. Pejchal J
    14. Hrabakova R
    15. Benes V
    16. Rausch T
    17. Dundr P
    18. Pilin A
    19. Cabala R
    20. Hubalek M
    21. Stríbrny J
    22. Antwerpen MH
    23. Meyer H

    (2017) Characterization of two historic smallpox specimens from a czech museum

    Viruses 9:v9080200.

    https://doi.org/10.3390/v9080200

    • PubMed
    • Google Scholar
84. 1. Paraskevis D
    2. Angelis K
    3. Magiorkinis G
    4. Kostaki E
    5. Ho SYW
    6. Hatzakis A

    (2015) Dating the origin of hepatitis B virus reveals higher substitution rate and adaptation on the branch leading to F/H genotypes

    Molecular Phylogenetics and Evolution 93:44–54.

    https://doi.org/10.1016/j.ympev.2015.07.010

    • PubMed
    • Google Scholar
85. 1. Patterson N
    2. Price AL
    3. Reich D

    (2006) Population Structure and Eigenanalysis

    PLOS Genetics 2:e190.

    https://doi.org/10.1371/journal.pgen.0020190

    • PubMed
    • Google Scholar
86. 1. Patterson Ross Z
    2. Klunk J
    3. Fornaciari G
    4. Giuffra V
    5. Duchêne S
    6. Duggan AT
    7. Poinar D
    8. Douglas MW
    9. Eden JS
    10. Holmes EC
    11. Poinar HN

    (2018) The paradox of HBV evolution as revealed from a 16th century mummy

    PLOS Pathogens 14:e1006750.

    https://doi.org/10.1371/journal.ppat.1006750

    • PubMed
    • Google Scholar
87. 1. Peucat JJ
    2. Capdevila R
    3. Drareni A
    4. Mahdjoub Y
    5. Kahoui M

    (2005) The Eglab massif in the West African Craton (Algeria), an original segment of the Eburnean orogenic belt: Petrology, geochemistry and geochronology

    Precambrian Research 136:309–352.

    https://doi.org/10.1016/j.precamres.2004.12.002

    • Google Scholar
88. 1. Pourkarim MR
    2. Lemey P
    3. Amini-Bavil-Olyaee S
    4. Maes P
    5. Van Ranst M

    (2010) Novel hepatitis B virus subgenotype A6 in African-Belgian patients

    Journal of Clinical Virology: The Official Publication of the Pan American Society for Clinical Virology 47:93–96.

    https://doi.org/10.1016/j.jcv.2009.09.032

    • PubMed
    • Google Scholar
89. 1. Pourkarim MR
    2. Amini-Bavil-Olyaee S
    3. Kurbanov F
    4. Van Ranst M
    5. Tacke F

    (2014) Molecular identification of hepatitis B virus genotypes/subgenotypes: Revised classification hurdles and updated resolutions

    World Journal of Gastroenterology 20:7152–7168.

    https://doi.org/10.3748/wjg.v20.i23.7152

    • PubMed
    • Google Scholar
90. 1. Price AL
    2. Patterson NJ
    3. Plenge RM
    4. Weinblatt ME
    5. Shadick NA
    6. Reich D

    (2006) Principal components analysis corrects for stratification in genome-wide association studies

    Nature Genetics 38:904–909.

    https://doi.org/10.1038/ng1847

    • PubMed
    • Google Scholar
91. 1. Price TD
    2. Burton JH
    3. Cucina A
    4. Zabala P
    5. Frei R
    6. Tykot RH
    7. Tiesler V

    (2012) Isotopic Studies of Human Skeletal Remains from a Sixteenth to Seventeenth Century AD Churchyard in Campeche, Mexico: Diet, Place of Origin, and Age

    Current Anthropology 53:396–433.

    https://doi.org/10.1086/666492

    • Google Scholar
92. 1. Pruitt KD
    2. Tatusova T
    3. Maglott DR

    (2007) NCBI reference sequences (RefSeq): A curated non-redundant sequence database of genomes, transcripts and proteins

    Nucleic Acids Research 35:D61–D65.

    https://doi.org/10.1093/nar/gkl842

    • PubMed
    • Google Scholar
93. 1. Purcell S
    2. Neale B
    3. Todd-Brown K
    4. Thomas L
    5. Ferreira MAR
    6. Bender D
    7. Maller J
    8. Sklar P
    9. de Bakker PIW
    10. Daly MJ
    11. Sham PC

    (2007) PLINK: A Tool Set for Whole-Genome Association and Population-Based Linkage Analyses

    The American Journal of Human Genetics 81:559–575.

    https://doi.org/10.1086/519795

    • Google Scholar
94. 1. Pyöriä L
    2. Toppinen M
    3. Mäntylä E
    4. Hedman L
    5. Aaltonen LM
    6. Vihinen-Ranta M
    7. Ilmarinen T
    8. Söderlund-Venermo M
    9. Hedman K
    10. Perdomo MF

    (2017) Extinct type of human parvovirus b19 persists in tonsillar b cells

    Nature Communications 8:14930.

    https://doi.org/10.1038/ncomms14930

    • PubMed
    • Google Scholar
95. 1. Quintero A
    2. Martínez D
    3. Alarcón De Noya B
    4. Costagliola A
    5. Urbina L
    6. González N
    7. Liprandi F
    8. Castro De Guerra D
    9. Pujol FH

    (2002) Molecular epidemiology of hepatitis B virus in Afro-Venezuelan populations

    Archives of Virology 147:1829–1836.

    https://doi.org/10.1007/s00705-002-0842-2

    • PubMed
    • Google Scholar
96. 1. Rambaut A
    2. Lam TT
    3. Max Carvalho L
    4. Pybus OG

    (2016) Exploring the temporal structure of heterochronous sequences using TempEst (formerly Path-O-Gen)

    Virus Evolution 2:vew007.

    https://doi.org/10.1093/ve/vew007

    • Google Scholar
97. 1. Rambaut A
    2. Drummond AJ
    3. Xie D
    4. Baele G
    5. Suchard MA

    (2018) Posterior Summarization in Bayesian Phylogenetics Using Tracer 1.7

    Systematic Biology 67:901–904.

    https://doi.org/10.1093/sysbio/syy032

    • PubMed
    • Google Scholar
98. 1. Reimer PJ
    2. Bard E
    3. Bayliss A
    4. Beck JW
    5. Blackwell PG
    6. Ramsey CB
    7. Buck CE
    8. Cheng H
    9. Edwards RL
    10. Friedrich M
    11. Grootes PM
    12. Guilderson TP
    13. Haflidason H
    14. Hajdas I
    15. Hatté C
    16. Heaton TJ
    17. Hoffmann DL
    18. Hogg AG
    19. Hughen KA
    20. van der Plicht J

    (2013) IntCal13 and Marine13 Radiocarbon Age Calibration Curves 0–50,000 Years cal BP

    Radiocarbon 55:1869–1887.

    https://doi.org/10.2458/azu_js_rc.55.16947

    • Google Scholar
99. 1. Renaud G
    2. Slon V
    3. Duggan AT
    4. Kelso J

    (2015) Schmutzi: Estimation of contamination and endogenous mitochondrial consensus calling for ancient DNA

    Genome Biology 16:224.

    https://doi.org/10.1186/s13059-015-0776-0

    • Google Scholar
100. 1. Rezaei F
     2. Sarshari B
     3. Ghavami N
     4. Meysami P
     5. Shadab A
     6. Salimi H
     7. Mokhtari-Azad T

     (2016) Prevalence and genotypic characterization of Human Parvovirus B19 in children with measles- and rubella-like illness in Iran

     Journal of Medical Virology 88:947–953.

     https://doi.org/10.1002/jmv.24425

     • PubMed
     • Google Scholar
101. 1. Rieux A
     2. Balloux F

     (2016) Inferences from tip-calibrated phylogenies: A review and a practical guide

     Molecular Ecology 25:1911–1924.

     https://doi.org/10.1111/mec.13586

     • PubMed
     • Google Scholar
102. 1. Rieux A
     2. Khatchikian CE

     (2017) Tipdatingbeast an r package to assist the implementation of phylogenetic tip-dating tests using beast

     Molecular Ecology Resources 17:608–613.

     https://doi.org/10.1111/1755-0998.12603

     • PubMed
     • Google Scholar
103. 1. Rinckel LA
     2. Buno BR
     3. Gierman TM
     4. Lee DC

     (2009) Discovery and analysis of a novel parvovirus B19 Genotype 3 isolate in the United States

     Transfusion 49:1488–1492.

     https://doi.org/10.1111/j.1537-2995.2009.02160.x

     • PubMed
     • Google Scholar
104. 1. Rohland N
     2. Hofreiter M

     (2007) Ancient DNA extraction from bones and teeth

     Nature Protocols 2:1756–1762.

     https://doi.org/10.1038/nprot.2007.247

     • PubMed
     • Google Scholar
105. 1. Roman S
     2. Tanaka Y
     3. Khan A
     4. Kurbanov F
     5. Kato H
     6. Mizokami M
     7. Panduro A

     (2010) Occult hepatitis B in the genotype H-infected Nahuas and Huichol native Mexican population

     Journal of Medical Virology 82:1527–1536.

     https://doi.org/10.1002/jmv.21846

     • PubMed
     • Google Scholar
106. 1. Rotimi CN
     2. Tekola-Ayele F
     3. Baker JL
     4. Shriner D

     (2016) The African Diaspora: History, Adaptation and Health

     Current Opinion in Genetics & Development 41:77–84.

     https://doi.org/10.1016/j.gde.2016.08.005

     • Google Scholar
107. Book

     1. Ruíz-Albarrán P

     (2012)

     Estudio de Variabilidad Biológica En La Colección Esquelética Hospital Real de San José de Los Naturales. Un Acercamiento A Tráves de La Técnica de Morfometría Geométrica

     Escuela Nacional de Antropología e Historia.

     • Google Scholar
108. 1. Sakyi PA
     2. Anum S
     3. Su BX
     4. Nude PM
     5. Su BC
     6. Asiedu DK
     7. Nyame F
     8. Kwayisi D

     (2018) Geochemical and Sr-Nd isotopic records of Paleoproterozoic metavolcanics and mafic intrusive rocks from the West African Craton: Evidence for petrogenesis and tectonic setting

     Geological Journal 53:725–741.

     https://doi.org/10.1002/gj.2923

     • Google Scholar
109. 1. Salas A
     2. Richards M
     3. Lareu MV
     4. Scozzari R
     5. Coppa A
     6. Torroni A
     7. Macaulay V
     8. Carracedo Á

     (2004) The African Diaspora: Mitochondrial DNA and the Atlantic Slave Trade

     American Journal of Human Genetics 74:454–465.

     https://doi.org/10.1086/382194

     • PubMed
     • Google Scholar
110. 1. Sanabani S
     2. Neto WK
     3. Pereira J
     4. Sabino EC

     (2006) Sequence Variability of Human Erythroviruses Present in Bone Marrow of Brazilian Patients with Various Parvovirus B19-Related Hematological Symptoms

     Journal of Clinical Microbiology 44:604–606.

     https://doi.org/10.1128/JCM.44.2.604-606.2006

     • PubMed
     • Google Scholar
111. 1. Schneider B
     2. Höne A
     3. Tolba RH
     4. Fischer H-P
     5. Blümel J
     6. Eis-Hübinger AM

     (2008) Simultaneous persistence of multiple genome variants of human Parvovirus b19

     Journal of General Virology 89:164–176.

     https://doi.org/10.1099/vir.0.83053-0

     • Google Scholar
112. 1. Schrick L
     2. Tausch SH
     3. Dabrowski PW
     4. Damaso CR
     5. Esparza J
     6. Nitsche A

     (2017) An Early American Smallpox Vaccine Based on Horsepox

     New England Journal of Medicine 377:1491–1492.

     https://doi.org/10.1056/NEJMc1707600

     • Google Scholar
113. 1. Schroeder H
     2. O’Connell TC
     3. Evans JA
     4. Shuler KA
     5. Hedges REM

     (2009) Trans-Atlantic slavery: Isotopic evidence for forced migration to Barbados

     American Journal of Physical Anthropology 139:547–557.

     https://doi.org/10.1002/ajpa.21019

     • PubMed
     • Google Scholar
114. 1. Schubert M
     2. Ginolhac A
     3. Lindgreen S
     4. Thompson JF
     5. Al-Rasheid KAS
     6. Willerslev E
     7. Krogh A
     8. Orlando L

     (2012) Improving ancient DNA read mapping against modern reference genomes

     BMC Genomics 13:178.

     https://doi.org/10.1186/1471-2164-13-178

     • PubMed
     • Google Scholar
115. 1. Schubert M
     2. Lindgreen S
     3. Orlando L

     (2016) AdapterRemoval v2: Rapid adapter trimming, identification, and read merging

     BMC Research Notes 9:88.

     https://doi.org/10.1186/s13104-016-1900-2

     • PubMed
     • Google Scholar
116. 1. Servant A
     2. Laperche S
     3. Lallemand F
     4. Marinho V
     5. De Saint Maur G
     6. Meritet JF
     7. Garbarg-Chenon A

     (2002) Genetic Diversity within Human Erythroviruses: Identification of Three Genotypes

     Journal of Virology 76:9124–9134.

     https://doi.org/10.1128/jvi.76.18.9124-9134.2002

     • PubMed
     • Google Scholar
117. 1. Skoglund P
     2. Storå J
     3. Götherström A
     4. Jakobsson M

     (2013) Accurate sex identification of ancient human remains using DNA shotgun sequencing

     Journal of Archaeological Science 40:4477–4482.

     https://doi.org/10.1016/j.jas.2013.07.004

     • Google Scholar
118. Book

     1. Solís-Pichardo G
     2. Schaaf P
     3. Hernández-Treviño T
     4. Lailson B
     5. Manzanilla LR
     6. Horn P

     (2017)

     Migration in Teopancazco, Teotihuacan: Evidence from Sr isotopic studies

     In: Solís-Pichardo G, editors. Multiethnicity and Migration at Teopancazco. Investigations of a Teotihuacan Neighborhood Center. University Florida Press. pp. 143–163.

     • Google Scholar
119. Book

     1. Somolinos d’Árdois G

     (1982)

     Las epidemias en México Durante El Siglo XVI

     In: Florescano E, Malvido E, editors. Ensayos Sobre La Historia de Las Epidemias En México. Instituto Mexicano del Seguro Social. pp. 138–143.

     • Google Scholar
120. 1. Spyrou MA
     2. Bos K
     3. Herbig A
     4. Krause J

     (2019) Ancient pathogen genomics as an emerging tool for infectious disease research

     Nature Reviews. Genetics 20:323–340.

     https://doi.org/10.1038/s41576-019-0119-1

     • Google Scholar
121. 1. Stamatakis A

     (2014) RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies

     Bioinformatics 30:1312–1313.

     https://doi.org/10.1093/bioinformatics/btu033

     • PubMed
     • Google Scholar
122. 1. Tithi SS
     2. Aylward FO
     3. Jensen R
     4. Zhang L

     (2018) FastViromeExplorer: A pipeline for virus and phage identification and abundance profiling in metagenomics data

     PeerJ 6:e4227.

     https://doi.org/10.7717/peerj.4227

     • Google Scholar
123. 1. Valencia Pacheco G
     2. Nakazawa Ueji YE
     3. Rodríguez Dzul EA
     4. Angulo Ramírez A
     5. López Villanueva RF
     6. Quintal Ortíz IG
     7. Rosado Paredes EP

     (2017) Serological and molecular analysis of parvovirus B19 infection in Mayan women with systemic lupus erythematosus in Mexico

     Colombia Médica v48:105–112.

     https://doi.org/10.25100/cm.v48i3.2981

     • Google Scholar
124. 1. Walker PL
     2. Bathurst RR
     3. Richman R
     4. Gjerdrum T
     5. Andrushko VA

     (2009) The causes of porotic hyperostosis and cribra orbitalia: A reappraisal of the iron-deficiency-anemia hypothesis

     American Journal of Physical Anthropology 139:109–125.

     https://doi.org/10.1002/ajpa.21031

     • PubMed
     • Google Scholar
125. 1. Warinner C
     2. Herbig A
     3. Mann A
     4. Fellows Yates JA
     5. Weiß CL
     6. Burbano HA
     7. Orlando L
     8. Krause J

     (2017) A Robust Framework for Microbial Archaeology

     Annual Review of Genomics and Human Genetics 18:321–356.

     https://doi.org/10.1146/annurev-genom-091416-035526

     • PubMed
     • Google Scholar
126. 1. Weissensteiner H
     2. Pacher D
     3. Kloss-Brandstätter A
     4. Forer L
     5. Specht G
     6. Bandelt HJ
     7. Kronenberg F
     8. Salas A
     9. Schönherr S

     (2016) HaploGrep 2: Mitochondrial haplogroup classification in the era of high-throughput sequencing

     Nucleic Acids Research 44:W58.

     https://doi.org/10.1093/nar/gkw233

     • Google Scholar
127. Book

     1. Wesp JK

     (2017) Caring for Bodies or Simply Saving Souls: The Emergence of Institutional Care in Spanish Colonial America

     In: Tilley L, Schrenk AA, editors. New Developments in the Bioarchaeology of Care: Further Case Studies and Expanded Theory. Springer International Publishing. pp. 253–276.

     https://doi.org/10.1007/978-3-319-39901-0_13

     • Google Scholar
128. 1. Wood DE
     2. Lu J
     3. Langmead B

     (2019) Improved metagenomic analysis with Kraken 2

     Genome Biology 20:257.

     https://doi.org/10.1186/s13059-019-1891-0

     • PubMed
     • Google Scholar
129. 1. Xiao YL
     2. Kash JC
     3. Beres SB
     4. Sheng ZM
     5. Musser JM
     6. Taubenberger JK

     (2013) High-throughput RNA sequencing of a formalin-fixed, paraffin-embedded autopsy lung tissue sample from the 1918 influenza pandemic

     The Journal of Pathology 229:535–545.

     https://doi.org/10.1002/path.4145

     • PubMed
     • Google Scholar
130. 1. Yuen MF
     2. Chen DS
     3. Dusheiko GM
     4. Janssen HLA
     5. Lau DTY
     6. Locarnini SA
     7. Peters MG
     8. Lai CL

     (2018) Hepatitis B virus infection

     Nature Reviews. Disease Primers 4:18035.

     https://doi.org/10.1038/nrdp.2018.35

     • Google Scholar

Author details

1. Axel A Guzmán-Solís

   Laboratorio Internacional de Investigación sobre el Genoma Humano, Universidad Nacional Autónoma de México, Querétaro, Mexico

   Contribution

   Data curation, Formal analysis, Investigation, Methodology, Visualization, Writing – original draft, Writing – review and editing

   Competing interests

   None

   "This ORCID iD identifies the author of this article:" 0000-0001-6878-206X

2. Viridiana Villa-Islas

   Laboratorio Internacional de Investigación sobre el Genoma Humano, Universidad Nacional Autónoma de México, Querétaro, Mexico

   Contribution

   Data curation, Formal analysis, Methodology, Visualization, Writing – review and editing

   Competing interests

   None

3. Miriam J Bravo-López

   Laboratorio Internacional de Investigación sobre el Genoma Humano, Universidad Nacional Autónoma de México, Querétaro, Mexico

   Contribution

   Data curation, Methodology, Visualization, Writing – review and editing

   Competing interests

   none

4. Marcela Sandoval-Velasco

   Section for Evolutionary Genomics, The Globe Institute, Faculty of Health, University of Copenhagen, Copenhagen, Denmark

   Contribution

   Conceptualization, Funding acquisition, Resources, Writing – review and editing

   Competing interests

   None

5. Julie K Wesp

   Department of Sociology and Anthropology, North Carolina State University, Raleigh, United States

   Contribution

   Conceptualization, Funding acquisition, Investigation, Resources, Supervision, Writing – review and editing

   Competing interests

   None

6. Jorge A Gómez-Valdés

   Escuela Nacional de Antropología e Historia, Mexico City, Mexico

   Contribution

   Resources, Supervision

   Competing interests

   None

   "This ORCID iD identifies the author of this article:" 0000-0001-6996-2732

7. María de la Luz Moreno-Cabrera

   Instituto Nacional de Antropología e Historia, Mexico City, Mexico

   Contribution

   Resources

   Competing interests

   None

8. Alejandro Meraz

   Instituto Nacional de Antropología e Historia, Mexico City, Mexico

   Contribution

   Resources

   Competing interests

   None

9. Gabriela Solís-Pichardo

   Laboratorio Universitario de Geoquímica Isotópica (LUGIS), Instituto de Geología, Universidad Nacional Autónoma de México, Mexico City, Mexico

   Contribution

   Formal analysis, Methodology, Writing – review and editing

   Competing interests

   None

10. Peter Schaaf

    LUGIS, Instituto de Geofísica, Universidad Nacional Autónoma de México, Mexico City, Mexico

    Contribution

    Formal analysis, Methodology, Writing – review and editing

    Competing interests

    None

11. Benjamin R TenOever

    Department of Microbiology, Icahn School of Medicine at Mount Sinai, New York, United States

    Contribution

    Resources, Writing – review and editing

    Competing interests

    None

    "This ORCID iD identifies the author of this article:" 0000-0003-0324-3078

12. Daniel Blanco-Melo

    1. Department of Microbiology, Icahn School of Medicine at Mount Sinai, New York, United States
    2. Vaccine and Infectious Disease Division, Fred Hutchinson Cancer Research Center, Seattle, WA, United States

    Contribution

    Conceptualization, Data curation, Funding acquisition, Supervision, Visualization, Writing – original draft, Writing – review and editing

    For correspondence

    dblancom@fredhutch.org

    Competing interests

    None

    "This ORCID iD identifies the author of this article:" 0000-0002-0256-5019

13. María C Ávila Arcos

    Laboratorio Internacional de Investigación sobre el Genoma Humano, Universidad Nacional Autónoma de México, Querétaro, Mexico

    Contribution

    Conceptualization, Funding acquisition, Investigation, Project administration, Resources, Supervision, Writing – original draft, Writing – review and editing

    For correspondence

    mavila@liigh.unam.mx

    Competing interests

    None

    "This ORCID iD identifies the author of this article:" 0000-0003-1691-1696

• 6,672

  views

• 821

  downloads

• 30

  citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.