Differentially methylated genes involved in reproduction and ploidy levels in recent diploidized and tetraploidized Eragrostis curvula genotypes

Epigenetics studies changes in gene activity without changes in the DNA sequence. Methylation is an epigenetic mechanism important in many pathways, such as biotic and abiotic stresses, cell division, and reproduction. Eragrostis curvula is a grass species reproducing by apomixis, a clonal reproduction by seeds. This work employed the MCSeEd technique to identify deferentially methylated positions, regions, and genes in the CG, CHG, and CHH contexts in E. curvula genotypes with similar genomic backgrounds but with different reproductive modes and ploidy levels. In this way, we focused the analysis on the cvs. Tanganyika INTA (4x, apomictic), Victoria (2x, sexual), and Bahiense (4x, apomictic). Victoria was obtained from the diploidization of Tanganyika INTA, while Bahiense was produced from the tetraploidization of Victoria. This study showed that polyploid/apomictic genotypes had more differentially methylated positions and regions than the diploid sexual ones. Interestingly, it was possible to observe fewer differentially methylated positions and regions in CG than in the other contexts, meaning CG methylation is conserved across the genotypes regardless of the ploidy level and reproductive mode. In the comparisons between sexual and apomictic genotypes, we identified differentially methylated genes involved in the reproductive pathways, specifically in meiosis, cell division, and fertilization. Another interesting observation was that several differentially methylated genes between the diploid and the original tetraploid genotype recovered their methylation status after tetraploidization, suggesting that methylation is an important mechanism involved in reproduction and ploidy changes.

Genetic legacy and adaptive signatures: investigating the history, diversity, and selection signatures in Rendena cattle resilient to eighteenth century rinderpest epidemics

BACKGROUND

Rendena is a dual-purpose cattle breed, which is primarily found in the Italian Alps and the eastern areas of the Po valley, and recognized for its longevity, fertility, disease resistance and adaptability to steep Alpine pastures. It is categorized as 'vulnerable to extinction' with only 6057 registered animals in 2022, yet no comprehensive analyses of its molecular diversity have been performed to date. The aim of this study was to analyse the origin, genetic diversity, and genomic signatures of selection in Rendena cattle using data from samples collected in 2000 and 2018, and shed light on the breed's evolution and conservation needs.

RESULTS

Genetic analysis revealed that the Rendena breed shares genetic components with various Alpine and Po valley breeds, with a marked genetic proximity to the Original Braunvieh breed, reflecting historical restocking efforts across the region. The breed shows signatures of selection related to both milk and meat production, environmental adaptation and immune response, the latter being possibly the result of multiple rinderpest epidemics that swept across the Alps in the eighteenth century. An analysis of the Rendena cattle population spanning 18 years showed an increase in the mean level of inbreeding over time, which is confirmed by the mean number of runs of homozygosity per individual, which was larger in the 2018 sample.

CONCLUSIONS

The Rendena breed, while sharing a common origin with Brown Swiss, has developed distinct traits that enable it to thrive in the Alpine environment and make it highly valued by local farmers. Preserving these adaptive features is essential, not only for maintaining genetic diversity and enhancing the ability of this traditional animal husbandry to adapt to changing environments, but also for guaranteeing the resilience and sustainability of both this livestock system and the livelihoods within the Rendena valley.

Global genetic diversity, introgression, and evolutionary adaptation of indicine cattle revealed by whole genome sequencing

Indicine cattle, also referred to as zebu (Bos taurus indicus), play a central role in pastoral communities across a wide range of agro-ecosystems, from extremely hot semiarid regions to hot humid tropical regions. However, their adaptive genetic changes following their dispersal into East Asia from the Indian subcontinent have remained poorly documented. Here, we characterize their global genetic diversity using high-quality whole-genome sequencing data from 354 indicine cattle of 57 breeds/populations, including major indicine phylogeographic groups worldwide. We reveal their probable migration into East Asia was along a coastal route rather than inland routes and we detected introgression from other bovine species. Genomic regions carrying morphology-, immune-, and heat-tolerance-related genes underwent divergent selection according to Asian agro-ecologies. We identify distinct sets of loci that contain promising candidate variants for adaptation to hot semi-arid and hot humid tropical ecosystems. Our results indicate that the rapid and successful adaptation of East Asian indicine cattle to hot humid environments was promoted by localized introgression from banteng and/or gaur. Our findings provide insights into the history and environmental adaptation of indicine cattle.

The genetic history of the Southern Andes from present-day Mapuche ancestry

The southernmost regions of South America harbor some of the earliest evidence of human presence in the Americas. However, connections with the rest of the continent and the contextualization of present-day indigenous ancestries remain poorly resolved. In this study, we analyze the genetic ancestry of one of the largest indigenous groups in South America: the Mapuche. We generate genome-wide data from 64 participants from three Mapuche populations in Southern Chile: Pehuenche, Lafkenche, and Huilliche. Broadly, we describe three main ancestry blocks with a common origin, which characterize the Southern Cone, the Central Andes, and Amazonia. Within the Southern Cone, ancestors of the Mapuche lineages differentiated from those of the Far South during the Middle Holocene and did not experience further migration waves from the north. We find that the deep genetic split between the Central and Southern Andes is followed by instances of gene flow, which may have accompanied the southward spread of cultural traits from the Central Andes, including crops and loanwords from Quechua into Mapudungun (the language of the Mapuche). Finally, we report close genetic relatedness between the three populations analyzed, with the Huilliche characterized additionally by intense recent exchanges with the Far South. Our findings add new perspectives on the genetic (pre)history of South America, from the first settlement through to the present-day indigenous presence. Follow-up fieldwork took these results back to the indigenous communities to contextualize the genetic narrative alongside indigenous knowledge and perspectives.

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.

The multifaceted genomic history of Ashaninka from Amazonian Peru

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

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.

Identification of genomic regions associated with total and progressive sperm motility in Italian Holstein bulls

Sperm motility is directly related to the ability of sperm to move through the female reproductive tract to reach the ovum. Sperm motility is a complex trait that is influenced by environmental and genetic factors and is associated with male fertility, oocyte penetration rate, and reproductive success of cattle. In this study we carried out a GWAS in Italian Holstein bulls to identify candidate regions and genes associated with variations in progressive and total motility (PM and TM, respectively). After quality control, the final data set consisted of 5,960 records from 949 bulls having semen collected in 10 artificial insemination stations and genotyped at 412,737 SNPs (call rate >95%; minor allele frequency >5%). (Co)variance components were estimated using single trait mixed models, and associations between SNPs and phenotypes were assessed using a genomic BLUP approach. Ten windows that explained the greatest percentage of genetic variance were located on Bos taurus autosomes 1, 2, 4, 6, 7, 23, and 26 for TM and Bos taurus autosomes 1, 2, 4, 6, 8, 16, 23, and 26 for PM. A total of 150 genes for TM and 72 genes for PM were identified within these genomic regions. Gene Ontology enrichment analyses identified significant Gene Ontology terms involved with energy homeostasis, membrane functions, sperm-egg interactions, protection against oxidative stress, olfactory receptors, and immune system. There was significant enrichment of quantitative trait loci for fertility, calving ease, immune response, feed intake, and carcass weight within the candidate windows. These results contribute to understanding the architecture of the genetic control of sperm motility and may aid in the development of strategies to identify subfertile bulls and improve reproductive success.

Genetic structure and differentiation from early bronze age in the mediterranean island of sicily: Insights from ancient mitochondrial genomes

Sicily is one of the main islands of the Mediterranean Sea, and it is characterized by a variety of archaeological records, material culture and traditions, reflecting the history of migrations and populations’ interaction since its first colonization, during the Paleolithic. These deep and complex demographic and cultural dynamics should have affected the genomic landscape of Sicily at different levels; however, the relative impact of these migrations on the genomic structure and differentiation within the island remains largely unknown. The available Sicilian modern genetic data gave a picture of the current genetic structure, but the paucity of ancient data did not allow so far to make predictions about the level of historical variation. In this work, we sequenced and analyzed the complete mitochondrial genomes of 36 individuals from five different locations in Sicily, spanning from Early Bronze Age to Iron Age, and with different cultural backgrounds. The comparison with coeval groups from the Mediterranean Basin highlighted structured genetic variation in Sicily since Early Bronze Age, thus supporting a demic impact of the cultural transitions within the Island. Explicit model testing through Approximate Bayesian Computation allowed us to make predictions about the origin of Sicanians, one of the three indigenous peoples of Sicily, whose foreign origin from Spain, historically attributed, was not confirmed by our analysis of genetic data. Sicilian modern mitochondrial data show a different, more homogeneous, genetic composition, calling for a recent genetic replacement in the Island of pre-Iron Age populations, that should be further investigated.

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.

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 (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.

Mitochondrial DNA Footprints from Western Eurasia in Modern Mongolia

Mongolia is located in a strategic position at the eastern edge of the Eurasian Steppe. Nomadic populations moved across this wide area for millennia before developing more sedentary communities, extended empires, and complex trading networks, which connected western Eurasia and eastern Asia until the late Medieval period. We provided a fine-grained portrait of the mitochondrial DNA (mtDNA) variation observed in present-day Mongolians and capable of revealing gene flows and other demographic processes that took place in Inner Asia, as well as in western Eurasia. The analyses of a novel dataset (N = 2,420) of mtDNAs highlighted a clear matrilineal differentiation within the country due to a mixture of haplotypes with eastern Asian (EAs) and western Eurasian (WEu) origins, which were differentially lost and preserved. In a wider genetic context, the prevalent EAs contribution, larger in eastern and central Mongolian regions, revealed continuous connections with neighboring Asian populations until recent times, as attested by the geographically restricted haplotype-sharing likely facilitated by the Genghis Khan’s so-called Pax Mongolica. The genetic history beyond the WEu haplogroups, notably detectable on both sides of Mongolia, was more difficult to explain. For this reason, we moved to the analysis of entire mitogenomes (N = 147). Although it was not completely possible to identify specific lineages that evolved in situ, two major changes in the effective (female) population size were reconstructed. The more recent one, which began during the late Pleistocene glacial period and became steeper in the early Holocene, was probably the outcome of demographic events connected to western Eurasia. The Neolithic growth could be easily explained by the diffusion of dairy pastoralism, as already proposed, while the late glacial increase indicates, for the first time, a genetic connection with western Eurasian refuges, as supported by the unusual high frequency and internal sub-structure in Mongolia of haplogroup H1, a well-known post-glacial marker in Europe. Bronze Age events, without a significant demographic impact, might explain the age of some mtDNA haplogroups. Finally, a diachronic comparison with available ancient mtDNAs made it possible to link six mitochondrial lineages of present-day Mongolians to the timeframe and geographic path of the Silk Route.

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.

Evaluating the Impact of Sex-Biased Genetic Admixture in the Americas through the Analysis of Haplotype Data

A general imbalance in the proportion of disembarked males and females in the Americas has been documented during the Trans-Atlantic Slave Trade and the Colonial Era and, although less prominent, more recently. This imbalance may have left a signature on the genomes of modern-day populations characterised by high levels of admixture. The analysis of the uniparental systems and the evaluation of continental proportion ratio of autosomal and X chromosomes revealed a general sex imbalance towards males for European and females for African and Indigenous American ancestries. However, the consistency and degree of this imbalance are variable, suggesting that other factors, such as cultural and social practices, may have played a role in shaping it. Moreover, very few investigations have evaluated the sex imbalance using haplotype data, containing more critical information than genotypes. Here, we analysed genome-wide data for more than 5000 admixed American individuals to assess the presence, direction and magnitude of sex-biased admixture in the Americas. For this purpose, we applied two haplotype-based approaches, ELAI and NNLS, and we compared them with a genotype-based method, ADMIXTURE. In doing so, besides a general agreement between methods, we unravelled that the post-colonial admixture dynamics show higher complexity than previously described.

The Mitochondrial DNA Landscape of Modern Mexico

Mexico is a rich source for anthropological and population genetic studies with high diversity in ethnic and linguistic groups. The country witnessed the rise and fall of major civilizations, including the Maya and Aztec, but resulting from European colonization, the population landscape has dramatically changed. Today, the majority of Mexicans do not identify themselves as Indigenous but as admixed, and appear to have very little in common with their pre-Columbian predecessors. However, when the maternally inherited mitochondrial (mt)DNA is investigated in the modern Mexican population, this is not the case. Control region sequences of 2021 samples deriving from all over the country revealed an overwhelming Indigenous American legacy, with almost 90% of mtDNAs belonging to the four major pan-American haplogroups A2, B2, C1, and D1. This finding supports a very low European contribution to the Mexican gene pool by female colonizers and confirms the effectiveness of employing uniparental markers as a tool to reconstruct a country’s history. In addition, the distinct frequency and dispersal patterns of Indigenous American and West Eurasian clades highlight the benefit such large and country-wide databases provide for studying the impact of colonialism from a female perspective and population stratification. The importance of geographical database subsets not only for forensic application is clearly demonstrated.

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.

Complete vertebrate mitogenomes reveal widespread repeats and gene duplications

BACKGROUND

Modern sequencing technologies should make the assembly of the relatively small mitochondrial genomes an easy undertaking. However, few tools exist that address mitochondrial assembly directly.

RESULTS

As part of the Vertebrate Genomes Project (VGP) we develop mitoVGP, a fully automated pipeline for similarity-based identification of mitochondrial reads and de novo assembly of mitochondrial genomes that incorporates both long (> 10 kbp, PacBio or Nanopore) and short (100-300 bp, Illumina) reads. Our pipeline leads to successful complete mitogenome assemblies of 100 vertebrate species of the VGP. We observe that tissue type and library size selection have considerable impact on mitogenome sequencing and assembly. Comparing our assemblies to purportedly complete reference mitogenomes based on short-read sequencing, we identify errors, missing sequences, and incomplete genes in those references, particularly in repetitive regions. Our assemblies also identify novel gene region duplications. The presence of repeats and duplications in over half of the species herein assembled indicates that their occurrence is a principle of mitochondrial structure rather than an exception, shedding new light on mitochondrial genome evolution and organization.

CONCLUSIONS

Our results indicate that even in the "simple" case of vertebrate mitogenomes the completeness of many currently available reference sequences can be further improved, and caution should be exercised before claiming the complete assembly of a mitogenome, particularly from short reads alone.

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.

Ancient genomes reveal tropical bovid species in the Tibetan Plateau contributed to the prevalence of hunting game until the late Neolithic

Local wild bovids have been determined to be important prey on the northeastern Tibetan Plateau (NETP), where hunting game was a major subsistence strategy until the late Neolithic, when farming lifestyles dominated in the neighboring Loess Plateau. However, the species affiliation and population ecology of these prehistoric wild bovids in the prehistoric NETP remain unknown. Ancient DNA (aDNA) analysis is highly informative in decoding this puzzle. Here, we applied aDNA analysis to fragmented bovid and rhinoceros specimens dating ∼5,200 y B.P. from the Neolithic site of Shannashuzha located in the marginal area of the NETP. Utilizing both whole genomes and mitochondrial DNA, our results demonstrate that the range of the present-day tropical gaur (Bos gaurus) extended as far north as the margins of the NETP during the late Neolithic from ∼29°N to ∼34°N. Furthermore, comparative analysis with zooarchaeological and paleoclimatic evidence indicated that a high summer temperature in the late Neolithic might have facilitated the northward expansion of tropical animals (at least gaur and Sumatran-like rhinoceros) to the NETP. This enriched the diversity of wildlife, thus providing abundant hunting resources for humans and facilitating the exploration of the Tibetan Plateau as one of the last habitats for hunting game in East Asia.

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.

Genomic analyses reveal distinct genetic architectures and selective pressures in buffaloes

BACKGROUND

The domestic buffalo (Bubalus bubalis) is an essential farm animal in tropical and subtropical regions, whose genomic diversity is yet to be fully discovered.

RESULTS

In this study, we describe the demographic events and selective pressures of buffalo by analyzing 121 whole genomes (98 newly reported) from 25 swamp and river buffalo breeds. Both uniparental and biparental markers were investigated to provide the final scenario. The ancestors of swamp and river buffalo diverged ∼0.23 million years ago and then experienced independent demographic histories. They were domesticated in different regions, the swamp buffalo at the border between southwest China and southeast Asia, while the river buffalo in south Asia. The domestic stocks migrated to other regions and further differentiated, as testified by (at least) 2 ancestral components identified in each subspecies. Different signals of selective pressures were also detected in these 2 types of buffalo. The swamp buffalo, historically used as a draft animal, shows selection signatures in genes associated with the nervous system, while in river dairy breeds, genes under selection are related to heat stress and immunity.

CONCLUSIONS

Our findings substantially expand the catalogue of genetic variants in buffalo and reveal new insights into the evolutionary history and distinct selective pressures in river and swamp buffalo.