Mitochondrial DNA control-region and coding-region data highlight geographically structured diversity and post-domestication population dynamics in worldwide donkeys.

Donkeys (Equus asinus) have been used extensively in agriculture and transportations since their domestication, ca. 5000-7000 years ago, but the increased mechanization of the last century has largely spoiled their role as burden animals, particularly in developed countries. Consequently, donkey breeds and population sizes have been declining for decades, and the diversity contributed by autochthonous gene pools has been eroded. Here, we examined coding-region data extracted from 164 complete mitogenomes and 1392 donkey mitochondrial DNA (mtDNA) control-region sequences to (i) assess worldwide diversity, (ii) evaluate geographical patterns of variation, and (iii) provide a new nomenclature of mtDNA haplogroups. The topology of the Maximum Parsimony tree confirmed the two previously identified major clades, i.e. Clades 1 and 2, but also highlighted the occurrence of a deep-diverging lineage within Clade 2 that left a marginal trace in modern donkeys. Thanks to the identification of stable and highly diagnostic coding-region mutational motifs, the two lineages were renamed as haplogroup A and haplogroup B, respectively, to harmonize clade nomenclature with the standard currently adopted for other livestock species. Control-region diversity and population expansion metrics varied considerably between geographical areas but confirmed North-eastern Africa as the likely domestication center. The patterns of geographical distribution of variation analyzed through phylogenetic networks and AMOVA confirmed the co-occurrence of both haplogroups in all sampled populations, while differences at the regional level point to the joint effects of demography, past human migrations and trade following the spread of donkeys out of the domestication center. Despite the strong decline that donkey populations have undergone for decades in many areas of the world, the sizeable mtDNA variability we scored, and the possible identification of a new early radiating lineage further stress the need for an extensive and large-scale characterization of donkey nuclear genome diversity to identify hotspots of variation and aid the conservation of local breeds worldwide.

Mitogenome evidence shows two radiation events and dispersals of matrilineal ancestry from northern coastal China to the Americas and Japan.

Although it is widely recognized that the ancestors of Native Americans (NAs) primarily came from Siberia, the link between mitochondrial DNA (mtDNA) lineage D4h3a (typical of NAs) and D4h3b (found so far only in East China and Thailand) raises the possibility that the ancestral sources for early NAs were more variegated than hypothesized. Here, we analyze 216 contemporary (including 106 newly sequenced) D4h mitogenomes and 39 previously reported ancient D4h data. The results reveal two radiation events of D4h in northern coastal China, one during the Last Glacial Maximum and the other within the last deglaciation, which facilitated the dispersals of D4h sub-branches to different areas including the Americas and the Japanese archipelago. The coastal distributions of the NA (D4h3a) and Japanese lineages (D4h1a and D4h2), in combination with the Paleolithic archaeological similarities among Northern China, the Americas, and Japan, lend support to the coastal dispersal scenario of early NAs.

Genomic evidence for adaptation to tuberculosis in the Andes before European contact.

Most studies focusing on human high-altitude adaptation in the Andean highlands have thus far been focused on Peruvian populations. We present high-coverage whole genomes from Indigenous people living in the Ecuadorian highlands and perform multi-method scans to detect positive natural selection. We identified regions of the genome that show signals of strong selection to both cardiovascular and hypoxia pathways, which are distinct from those uncovered in Peruvian populations. However, the strongest signals of selection were related to regions of the genome that are involved in immune function related to tuberculosis. Given our estimated timing of this selection event, the Indigenous people of Ecuador may have adapted to Mycobacterium tuberculosis thousands of years before the arrival of Europeans. Furthermore, we detect a population collapse that coincides with the arrival of Europeans, which is more severe than other regions of the Andes, suggesting differing effects of contact across high-altitude populations.

A chromosome-level reference genome and pangenome for barn swallow population genomics.

Insights into the evolution of non-model organisms are limited by the lack of reference genomes of high accuracy, completeness, and contiguity. Here, we present a chromosome-level, karyotype-validated reference genome and pangenome for the barn swallow (Hirundo rustica). We complement these resources with a reference-free multialignment of the reference genome with other bird genomes and with the most comprehensive catalog of genetic markers for the barn swallow. We identify potentially conserved and accelerated genes using the multialignment and estimate genome-wide linkage disequilibrium using the catalog. We use the pangenome to infer core and accessory genes and to detect variants using it as a reference. Overall, these resources will foster population genomics studies in the barn swallow, enable detection of candidate genes in comparative genomics studies, and help reduce bias toward a single reference genome.

Dissecting the multifaceted contribution of the mitochondrial genome to autism spectrum disorder.

Autism spectrum disorder (ASD) is a clinically heterogeneous class of neurodevelopmental conditions with a strong, albeit complex, genetic basis. The genetic architecture of ASD includes different genetic models, from monogenic transmission at one end, to polygenic risk given by thousands of common variants with small effects at the other end. The mitochondrial DNA (mtDNA) was also proposed as a genetic modifier for ASD, mostly focusing on maternal mtDNA, since the paternal mitogenome is not transmitted to offspring. We extensively studied the potential contribution of mtDNA in ASD pathogenesis and risk through deep next generation sequencing and quantitative PCR in a cohort of 98 families. While the maternally-inherited mtDNA did not seem to predispose to ASD, neither for haplogroups nor for the presence of pathogenic mutations, an unexpected influence of paternal mtDNA, apparently centered on haplogroup U, came from the Italian families extrapolated from the test cohort (n = 74) when compared to the control population. However, this result was not replicated in an independent Italian cohort of 127 families and it is likely due to the elevated paternal age at time of conception. In addition, ASD probands showed a reduced mtDNA content when compared to their unaffected siblings. Multivariable regression analyses indicated that variants with 15%-5% heteroplasmy in probands are associated to a greater severity of ASD based on ADOS-2 criteria, whereas paternal super-haplogroups H and JT were associated with milder phenotypes. In conclusion, our results suggest that the mtDNA impacts on ASD, significantly modifying the phenotypic expression in the Italian population. The unexpected finding of protection induced by paternal mitogenome in term of severity may derive from a role of mtDNA in influencing the accumulation of nuclear de novo mutations or epigenetic alterations in fathers' germinal cells, affecting the neurodevelopment in the offspring. This result remains preliminary and needs further confirmation in independent cohorts of larger size. If confirmed, it potentially opens a different perspective on how paternal non-inherited mtDNA may predispose or modulate other complex diseases.

The worldwide spread of Aedes albopictus: New insights from mitogenomes.

The tiger mosquito (Aedes albopictus) is one of the most invasive species in the world and a competent vector for numerous arboviruses, thus the study and monitoring of its fast worldwide spread is crucial for global public health. The small extra-nuclear and maternally-inherited mitochondrial DNA represents a key tool for reconstructing phylogenetic and phylogeographic relationships within a species, especially when analyzed at the mitogenome level. Here the mitogenome variation of 76 tiger mosquitoes, 37 of which new and collected from both wild adventive populations and laboratory strains, was investigated. This analysis significantly improved the global mtDNA phylogeny of Ae. albopictus, uncovering new branches and sub-branches within haplogroup A1, the one involved in its recent worldwide spread. Our phylogeographic approach shows that the current distribution of tiger mosquito mitogenome variation has been strongly affected by clonal and sub-clonal founder events, sometimes involving wide geographic areas, even across continents, thus shedding light on the Asian sources of worldwide adventive populations. In particular, different starting points for the two major clades within A1 are suggested, with A1a spreading mainly along temperate areas from Japanese and Chinese sources, and A1b arising and mainly diffusing in tropical areas from a South Asian source.

Helena's Many Daughters: More Mitogenome Diversity behind the Most Common West Eurasian mtDNA Control Region Haplotype in an Extended Italian Population Sample.

The high number of matching haplotypes of the most common mitochondrial (mt)DNA lineages are considered to be the greatest limitation for forensic applications. This study investigates the potential to solve this constraint by massively parallel sequencing a large number of mitogenomes that share the most common West Eurasian mtDNA control region (CR) haplotype motif (263G 315.1C 16519C). We augmented a pilot study on 29 to a total of 216 Italian mitogenomes that represents the largest set of the most common CR haplotype compiled from a single country. The extended population sample confirmed and extended the huge coding region diversity behind the most common CR motif. Complete mitogenome sequencing allowed for the detection of 163 distinct haplotypes, raising the power of discrimination from 0 (CR) to 99.6% (mitogenome). The mtDNAs were clustered into 61 named clades of haplogroup H and did not reveal phylogeographic trends within Italy. Rapid individualization approaches for investigative purposes are limited to the most frequent H clades of the dataset, viz. H1, H3, and H7.

Assessing temporal and geographic contacts across the Adriatic Sea through the analysis of genome-wide data from Southern Italy.

Southern Italy was characterised by a complex prehistory that started with different Palaeolithic cultures, later followed by the Neolithization and the demic dispersal from the Pontic-Caspian Steppe during the Bronze Age. Archaeological and historical evidences point to a link between Southern Italians and the Balkans still present in modern times. To shed light on these dynamics, we analysed around 700 South Mediterranean genomes combined with informative ancient DNAs. Our findings revealed high affinities of South-Eastern Italians with modern Eastern Peloponnesians, and a closer affinity of ancient Greek genomes with those from specific regions of South Italy than modern Greek genomes. The higher similarity could be associated with a Bronze Age component ultimately originating from the Caucasus with high Iranian and Anatolian Neolithic ancestries. Furthermore, extremely differentiated allele frequencies among Northern and Southern Italy revealed putatively adapted SNPs in genes involved in alcohol metabolism, nevi features and immunological traits.

The Mitogenome Relationships and Phylogeography of Barn Swallows (Hirundo rustica).

The barn swallow (Hirundo rustica) poses a number of fascinating scientific questions, including the taxonomic status of postulated subspecies. Here, we obtained and assessed the sequence variation of 411 complete mitogenomes, mainly from the European H. r. rustica, but other subspecies as well. In almost every case, we observed subspecies-specific haplogroups, which we employed together with estimated radiation times to postulate a model for the geographical and temporal worldwide spread of the species. The female barn swallow carrying the Hirundo rustica ancestral mitogenome left Africa (or its vicinity) around 280 thousand years ago (kya), and her descendants expanded first into Eurasia and then, at least 51 kya, into the Americas, from where a relatively recent (<20 kya) back migration to Asia took place. The exception to the haplogroup subspecies specificity is represented by the sedentary Levantine H. r. transitiva that extensively shares haplogroup A with the migratory European H. r. rustica and, to a lesser extent, haplogroup B with the Egyptian H. r. savignii. Our data indicate that rustica and transitiva most likely derive from a sedentary Levantine population source that split at the end of the Younger Dryas (YD) (11.7 kya). Since then, however, transitiva received genetic inputs from and admixed with both the closely related rustica and the adjacent savignii. Demographic analyses confirm this species' strong link with climate fluctuations and human activities making it an excellent indicator for monitoring and assessing the impact of current global changes on wildlife.

Overview of the Americas' First Peopling from a Patrilineal Perspective: New Evidence from the Southern Continent.

Uniparental genetic systems are unique sex indicators and complement the study of autosomal diversity by providing landmarks of human migrations that repeatedly shaped the structure of extant populations. Our knowledge of the variation of the male-specific region of the Y chromosome in Native Americans is still rather scarce and scattered, but by merging sequence information from modern and ancient individuals, we here provide a comprehensive and updated phylogeny of the distinctive Native American branches of haplogroups C and Q. Our analyses confirm C-MPB373, C-P39, Q-Z780, Q-M848, and Q-Y4276 as the main founding haplogroups and identify traces of unsuccessful (pre-Q-F1096) or extinct (C-L1373*, Q-YP4010*) Y-chromosome lineages, indicating that haplogroup diversity of the founder populations that first entered the Americas was greater than that observed in the Indigenous component of modern populations. In addition, through a diachronic and phylogeographic dissection of newly identified Q-M848 branches, we provide the first Y-chromosome insights into the early peopling of the South American hinterland (Q-BY104773 and Q-BY15730) and on overlying inland migrations (Q-BY139813).

The role of mtDNA haplogroups on metabolic features in narcolepsy type 1.

Narcolepsy type 1 (NT1) is due to selective loss of hypocretin (hcrt)-producing-neurons. Hcrt is a neuropeptide regulating the sleep/wake cycle, as well as feeding behavior. A subset of NT1 patients become overweight/obese, with a dysmetabolic phenotype. We hypothesized that mitochondrial DNA (mtDNA) sequence variation might contribute to the metabolic features in NT1 and we undertook an exploratory survey of mtDNA haplogroups in a cohort of well-characterized patients. We studied 246 NT1 Italian patients, fully defined for their metabolic features, including obesity, hypertension, low HDL, hypertriglyceridemia and hyperglycemia. For haplogroup assignment, the mtDNA control region was sequenced in combination with an assessment of diagnostic markers in the coding region. NT1 patients displayed the same mtDNA haplogroups (H, HV, J, K, T, U) frequency as those reported in the general Italian population. The majority of NT1 patients (64%) were overweight: amongst these, 35% were obese, 48% had low HDL cholesterol levels, and 31% had hypertriglyceridemia. We identified an association between haplogroups J, K and hypertriglyceridemia (P = 0.03, 61.5% and 61.5%, respectively vs. 31.3% of the whole sample) and after correction for age and sex, we observed a reduction of these associations (OR = 3.65, 95%CI = 0.76-17.5, p = 0.106 and 1.73, 0.52-5.69, p = 0.368, respectively). The low HDL level showed a trend for association with haplogroup J (P = 0.09, 83.3% vs. 47.4% of the whole sample) and after correction we observed an OR = 6.73, 95%CI = 0.65-69.9, p = 0.110. Our study provides the first indication that mtDNA haplogroups J and K can modulate metabolic features of NT1 patients, linking mtDNA variation to the dysmetabolic phenotype in NT1.

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

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

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

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

Biomolecular insights into North African-related ancestry, mobility and diet in eleventh-century Al-Andalus.

Historical records document medieval immigration from North Africa to Iberia to create Islamic al-Andalus. Here, we present a low-coverage genome of an eleventh century CE man buried in an Islamic necropolis in Segorbe, near Valencia, Spain. Uniparental lineages indicate North African ancestry, but at the autosomal level he displays a mosaic of North African and European-like ancestries, distinct from any present-day population. Altogether, the genome-wide evidence, stable isotope results and the age of the burial indicate that his ancestry was ultimately a result of admixture between recently arrived Amazigh people (Berbers) and the population inhabiting the Peninsula prior to the Islamic conquest. We detect differences between our sample and a previously published group of contemporary individuals from Valencia, exemplifying how detailed, small-scale aDNA studies can illuminate fine-grained regional and temporal differences. His genome demonstrates how ancient DNA studies can capture portraits of past genetic variation that have been erased by later demographic shifts-in this case, most likely the seventeenth century CE expulsion of formerly Islamic communities as tolerance dissipated following the Reconquista by the Catholic kingdoms of the north.

Archaeogenomic distinctiveness of the Isthmo-Colombian area.

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

Cattle mitogenome variation reveals a post-glacial expansion of haplogroup P and an early incorporation into northeast Asian domestic herds.

Surveys of mitochondrial DNA (mtDNA) variation have shown that worldwide domestic cattle are characterized by just a few major haplogroups. Two, T and I, are common and characterize Bos taurus and Bos indicus, respectively, while the other three, P, Q and R, are rare and are found only in taurine breeds. Haplogroup P is typical of extinct European aurochs, while intriguingly modern P mtDNAs have only been found in northeast Asian cattle. These Asian P mtDNAs are extremely rare with the exception of the Japanese Shorthorn breed, where they reach a frequency of 45.9%. To shed light on the origin of this haplogroup in northeast Asian cattle, we completely sequenced 14 Japanese Shorthorn mitogenomes belonging to haplogroup P. Phylogenetic and Bayesian analyses revealed: (1) a post-glacial expansion of aurochs carrying haplogroup P from Europe to Asia; (2) that all Asian P mtDNAs belong to a single sub-haplogroup (P1a), so far never detected in either European or Asian aurochs remains, which was incorporated into domestic cattle of continental northeastern Asia possibly ~ 3700 years ago; and (3) that haplogroup P1a mtDNAs found in the Japanese Shorthorn breed probably reached Japan about 650 years ago from Mongolia/Russia, in agreement with historical evidence.

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

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

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

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

The Genomic Impact of European Colonization of the Americas.

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