Bioarchaeological study of ancient Teotihuacans based on complete mitochondrial genome sequences and diet isotopes.

BACKGROUND: The Teotihuacan civilisation was the largest one in ancient Mesoamerica. The Teotihuacan city was born in the north-eastern Basin of Mexico around the second century BC, reached its peak in the fourth century AD, and had cultural influence throughout Mesoamerica. At its peak, the size of the city reached more than 20 km2, and the total population is estimated to have increased from 100,000 to 200,000. However, knowledge of the genetic background of the Teotihuacan people is still limited. AIM: We aimed to determine the mitogenome sequences of the Teotihuacan human remains and compare the ancient and present Mesoamericans. In addition, we aimed to identify the food habits of ancient Teotihuacans. SUBJECTS AND METHODS: We determined the mitogenome sequences of human remains dated to 250-636 cal AD using target enrichment-coupled next generation sequencing. We also performed stable isotope analysis. RESULTS: We successfully obtained nearly full-length sequences newly unearthed from a civilian dwelling in the Teotihuacan site. Teotihuacan mitochondrial DNA was classified into the haplogroups in present and ancient Mesoamericans. In addition, Teotihuacan individuals had a diet dependent on C4 plants such as maize. CONCLUSION: Genetic diversity varied among the Teotihuacans.

Diversity in matrilineages among the Jomon individuals of Japan.

BACKGROUND: The Jomon period of Japan is characterised by a unique combination of sedentary and hunting/gathering lifestyles, spanning for more than 10,000 years from the final Pleistocene to the Holocene. The transition from the preceding Palaeolithic period to the Jomon period is known to have begun with the appearance of pottery usage. However, knowledge of the genetic background of the Jomon people is still limited. AIM: We aimed to determine the population-scale complete mitogenome sequences of the Initial Jomon human remains and compare the occurrence of mitochondrial haplogroups in the Jomon period from temporal and regional perspectives. SUBJECTS AND METHODS: For human remains dated to 8200-8600 cal BP, we determined their complete mitogenome sequences using target enrichment-coupled next-generation sequencing. RESULTS: We successfully obtained the complete mitogenome sequences with high depth of coverage and high concordance on consensus sequences. These sequences differed by more than three bases each, except for two individuals having completely identical sequences. Co-existence of individuals with haplogroups N9b and M7a was first observed at the same archaeological site from the Initial Jomon period. CONCLUSION: The genetic diversity within the population was not found to be low even in the Initial Jomon period.

Population dynamics in the Japanese Archipelago since the Pleistocene revealed by the complete mitochondrial genome sequences.

The Japanese Archipelago is widely covered with acidic soil made of volcanic ash, an environment which is detrimental to the preservation of ancient biomolecules. More than 10,000 Palaeolithic and Neolithic sites have been discovered nationwide, but few skeletal remains exist and preservation of DNA is poor. Despite these challenging circumstances, we succeeded in obtaining a complete mitogenome (mitochondrial genome) sequence from Palaeolithic human remains. We also obtained those of Neolithic (the hunting-gathering Jomon and the farming Yayoi cultures) remains, and over 2,000 present-day Japanese. The Palaeolithic mitogenome sequence was not found to be a direct ancestor of any of Jomon, Yayoi, and present-day Japanese people. However, it was an ancestral type of haplogroup M, a basal group of the haplogroup M. Therefore, our results indicate continuity in the maternal gene pool from the Palaeolithic to present-day Japanese. We also found that a vast increase of population size happened and has continued since the Yayoi period, characterized with paddy rice farming. It means that the cultural transition, i.e. rice agriculture, had significant impact on the demographic history of Japanese population.

The number of SNPs required for distinguishing Japanese from other East Asians.

In some cases, it is necessary to estimate the national origin of an unknown subject in forensic medicine. The use of single nucleotide polymorphism (SNP) markers appears to be very effective for this purpose, since genome-wide SNP genotype data of many human populations are publicly available. In this study, we examined the number of SNPs that could objectively and accurately distinguish Japanese subjects (1KG-JPT) from the other East Asians (1KG-CDX, -CHB, -CHS, and -KHV) using the combination of principal component analysis and hierarchical cluster analysis. A computer simulation showed that approximately 3000 randomly selected SNPs were enough for the discrimination. Our results suggest that at least a 0.024% coverage is needed in the next generation sequencing experiment to objectively determine whether an unknown person is Japanese or not if the amount of DNA sample from him/her is insufficient or the quality is low.

A biomolecular anthropological investigation of William Adams, the first SAMURAI from England.

William Adams (Miura Anjin) was an English navigator who sailed with a Dutch trading fleet to the far East and landed in Japan in 1600. He became a vassal under the Shogun, Tokugawa Ieyasu, was bestowed with a title, lands and swords, and became the first SAMURAI from England. "Miura" comes from the name of the territory given to him and "Anjin" means "pilot". He lived out the rest of his life in Japan and died in Hirado, Nagasaki Prefecture, in 1620, where he was reportedly laid to rest. Shortly after his death, graveyards designated for foreigners were destroyed during a period of Christian repression, but Miura Anjin's bones were supposedly taken, protected, and reburied. Archaeological investigations in 1931 uncovered human skeletal remains and it was proposed that they were those of Miura Anjin. However, this could not be confirmed from the evidence at the time and the remains were reburied. In 2017, excavations found skeletal remains matching the description of those reinterred in 1931. We analyzed these remains from various aspects, including genetic background, dietary habits, and burial style, utilizing modern scientific techniques to investigate whether they do indeed belong to the first English SAMURAI.

A study of 8,300-year-old Jomon human remains in Japan using complete mitogenome sequences obtained by next-generation sequencing.

Ancient human remains have been assigned to their mitochondrial DNA (mtDNA) haplogroups. To obtain efficiently deep and reliable nucleotide sequences of ancient DNA of interest, we achieved target enrichment followed by next-generation sequencing (NGS). Complete mitochondrial genome (mitogenome) sequences were obtained for three human remains from the Iyai rock-shelter site of the Initial Jomon Period in Japan. All the Jomon mitogenomes belong to haplogroup N9b, but no sequences among them were identical. High genetic diversity was clarified even among the Jomon human remains belonging to haplogroup N9b, which has been described as a haplogroup representing the Jomon people.

MitoIMP: A Computational Framework for Imputation of Missing Data in Low-Coverage Human Mitochondrial Genome.

The incompleteness of partial human mitochondrial genome sequences makes it difficult to perform relevant comparisons among multiple resources. To deal with this issue, we propose a computational framework for deducing missing nucleotides in the human mitochondrial genome. We applied it to worldwide mitochondrial haplogroup lineages and assessed its performance. Our approach can deduce the missing nucleotides with a precision of 0.99 or higher in most human mitochondrial DNA lineages. Furthermore, although low-coverage mitochondrial genome sequences often lead to a blurred relationship in the multidimensional scaling analysis, our approach can correct this positional arrangement according to the corresponding mitochondrial DNA lineages. Therefore, our framework will provide a practical solution to compensate for the lack of genome coverage in partial and fragmented human mitochondrial genome sequences. In this study, we developed an open-source computer program, MitoIMP, implementing our imputation procedure. MitoIMP is freely available from https://github.com/omics-tools/mitoimp.

Characterization of complete mitochondrial genomes of indigenous Mayans in Mexico.

BACKGROUND: The authors have previously published the complete mitochondrial genome (mitogenome) sequences of two indigenous Mesoamerican populations, Mazahua (n = 25) and Zapotec (n = 88). METHODS: This study determined the complete mitogenome sequences of nine unrelated individuals from the indigenous Maya population living in Mexico. RESULTS: Their mitogenome sequences could be classified into either of the haplogroups A2 and C1. Surprisingly, there were no mitogenome sequences (haplotypes) that the Maya, Mazahua, and Zapotec people share in common. CONCLUSIONS: This indicates that no genetic exchange, at least matrilineally, has occurred among them.

Imputation approach for deducing a complete mitogenome sequence from low-depth-coverage next-generation sequencing data: application to ancient remains from the Moon Pyramid, Mexico.

It is considered that more than 15 depths of coverage are necessary for next-generation sequencing (NGS) data to obtain reliable complete nucleotide sequences of the mitogenome. However, it is difficult to satisfy this requirement for all nucleotide positions because of problems obtaining a uniform depth of coverage for poorly preserved materials. Thus, we propose an imputation approach that allows a complete mitogenome sequence to be deduced from low-depth-coverage NGS data. We used different types of mitogenome data files as panels for imputation: a worldwide panel comprising all the major haplogroups, a worldwide panel comprising sequences belonging to the estimated haplogroup alone, a panel comprising sequences from the population most closely related to an individual under investigation, and a panel comprising sequences belonging to the estimated haplogroup from the population most closely related to an individual under investigation. The number of missing nucleotides was drastically reduced in all the panels, but the contents obtained by imputation were quite different among the panels. The efficiency of the imputation method differed according to the panels used. The missing nucleotides were most credibly imputed using sequences of the estimated haplogroup from the population most closely related to the individual under investigation as a panel.

mtDNA diversity of the Zapotec in Mexico suggests a population decline long before the first contact with Europeans.

The New World is the last continent colonized by anatomically modern humans, Homo sapiens. The first migrants entered the New World from Asia through Beringia. It is suggested that there were three streams of Asian gene flow, one major and two additional minor gene flows. The first major migrants took a Pacific coastal route and began spreading to the American continent before the opening of the ice-free corridor. We investigated the diversity of full-length mitochondrial DNA genomes of the Zapotec population, residing in the Mesoamerican region, and reconstructed their demographic history using Bayesian Skyline Plots. We estimated the initial date of gene flow into the New World by Zapotec ancestors at around 17 000­19 000 years ago,which is highly concordant with previous studies. We also show a population decline after the initial expansion. This decline started 4000 years ago, long before European contact with Native Americans. This indicates that other factors including climatec hange should be considered to explain the observed demographic pattern.

Complete mitogenome analysis of indigenous populations in Mexico: its relevance for the origin of Mesoamericans.

Mesoamerica has an important role in the expansion of Paleoamericans as the route to South America. In this study, we determined complete mitogenome sequences of 113 unrelated individuals from two indigenous populations of Mesoamerica, Mazahua and Zapotec. All newly sequenced mitogenomes could be classified into haplogroups A2, B2, C1 and D1, but one sequence in Mazahua was D4h3a, a subclade of haplogroup D4. This haplogroup has been mostly found in South America along the Pacific coast. Haplogroup X2a was not found in either population. Genetic similarity obtained using phylogenetic tree construction and principal component analysis showed that these two populations are distantly related to each other. Actually, the Mazahua and the Zapotec shared no sequences (haplotypes) in common, while each also showed a number of unique subclades. Surprisingly, Zapotec formed a cluster with indigenous populations living in an area from central Mesoamerica to Central America. By contrast, the Mazahua formed a group with indigenous populations living in external areas, including southwestern North America and South America. This intriguing genetic relationship suggests the presence of two paleo-Mesoamerican groups, invoking a scenario in which one group had expanded into South America and the other resided in Mesoamerica.

Emulsion PCR-coupled target enrichment: an effective fishing method for high-throughput sequencing of poorly preserved ancient DNA.

Due to the difficulties in deep sequencing, high-throughput sequencing of ancient DNA has been limited to exceptionally well-preserved ancient materials. The primary factor is microbial attack popularly observed in the buried materials, and it causes drastic increase in relative ratio of microbial DNA in the extracted DNA. We present a unified strategy in which emulsion PCR is coupled with target enrichment followed by next-generation sequencing. The method made it possible to obtain efficiently non-duplicated reads mapped to target sequences of interest, and this can achieve deep and reliable sequencing of ancient DNA from typical materials, even though poorly preserved.

Comparative aspects of polyglutamine binding domain in PQBP-1 among Vertebrata.

We investigated the evolutionary conservation of polyglutamine binding protein-1 (PQBP-1) among Vertebrata. PQBP-1s were highly conserved and shared the same domain features including a WW domain, a polar amino acid rich domain (PRD), a nuclear localization signal (NLS), and a C-terminal domain (CTD) among Eutheria, but not always among Vertebrata. PQBP-1s of Vertebrata contained a variable region in the middle portion corresponding to the position of PRD. The full form of PRD including both 7aa and DR/ER repeats was specific to Eutheria. PRD of non-eutherian Amniota was minimal. Amphibia had no PRD. The DR/ER repeat was solo in fishes. Agnatha PRD was also rich in polar amino acids, but contained no repetitive sequence. We investigated 3 polyQ-containing proteins known to interact with PQBP-1: BRN-2, Huntingtin, and ATAXIN-1, and showed a diverse nature of protein-protein interaction in Vertebrata. There appears to be no interaction between PQBP-1 and BRN-2, Huntingtin, or ATAXIN-1 in Amphibia, while the interaction between PQBP-1 and BRN-2 is expected to be conserved among Mammalia, and the interaction between PQBP-1 and Huntingtin or ATAXIN-1 depends on the lineage in Eutheria.