Evaluation of DNA Methylation Profiles of LINE-1, Alu and Ribosomal DNA Repeats in Human Cell Lines Exposed to Radiofrequency Radiation.

A large body of evidence indicates that environmental agents can induce alterations in DNA methylation (DNAm) profiles. Radiofrequency electromagnetic fields (RF-EMFs) are radiations emitted by everyday devices, which have been classified as "possibly carcinogenic"; however, their biological effects are unclear. As aberrant DNAm of genomic repetitive elements (REs) may promote genomic instability, here, we sought to determine whether exposure to RF-EMFs could affect DNAm of different classes of REs, such as long interspersed nuclear elements-1 (LINE-1), Alu short interspersed nuclear elements and ribosomal repeats. To this purpose, we analysed DNAm profiles of cervical cancer and neuroblastoma cell lines (HeLa, BE(2)C and SH-SY5Y) exposed to 900 MHz GSM-modulated RF-EMF through an Illumina-based targeted deep bisulfite sequencing approach. Our findings showed that radiofrequency exposure did not affect the DNAm of Alu elements in any of the cell lines analysed. Conversely, it influenced DNAm of LINE-1 and ribosomal repeats in terms of both average profiles and organisation of methylated and unmethylated CpG sites, in different ways in each of the three cell lines studied.

A Targeted Epigenetic Clock for the Prediction of Biological Age.

Epigenetic clocks were initially developed to track chronological age, but accumulating evidence indicates that they can also predict biological age. They are usually based on the analysis of DNA methylation by genome-wide methods, but targeted approaches, based on the assessment of a small number of CpG sites, are advisable in several settings. In this study, we developed a targeted epigenetic clock purposely optimized for the measurement of biological age. The clock includes six genomic regions mapping in ELOVL2, NHLRC1, AIM2, EDARADD, SIRT7 and TFAP2E genes, selected from a re-analysis of existing microarray data, whose DNA methylation is measured by EpiTYPER assay. In healthy subjects (n = 278), epigenetic age calculated using the targeted clock was highly correlated with chronological age (Spearman correlation = 0.89). Most importantly, and in agreement with previous results from genome-wide clocks, epigenetic age was significantly higher and lower than expected in models of increased (persons with Down syndrome, n = 62) and decreased (centenarians, n = 106; centenarians' offspring, n = 143; nutritional intervention in elderly, n = 233) biological age, respectively. These results support the potential of our targeted epigenetic clock as a new marker of biological age and open its evaluation in large cohorts to further promote the assessment of biological age in healthcare practice.

DNA Methylation Analysis of Ribosomal DNA in Adults With Down Syndrome.

Control of ribosome biogenesis is a critical aspect of the regulation of cell metabolism. As ribosomal genes (rDNA) are organized in repeated clusters on chromosomes 13, 14, 15, 21, and 22, trisomy of chromosome 21 confers an excess of rDNA copies to persons with Down syndrome (DS). Previous studies showed an alteration of ribosome biogenesis in children with DS, but the epigenetic regulation of rDNA genes has not been investigated in adults with DS so far. In this study, we used a targeted deep-sequencing approach to measure DNA methylation (DNAm) of rDNA units in whole blood from 69 adults with DS and 95 euploid controls. We further evaluated the expression of the precursor of ribosomal RNAs (RNA45S) in peripheral blood mononuclear cells (PBMCs) from the same subjects. We found that the rDNA promoter tends to be hypermethylated in DS concerning the control group. The analysis of epihaplotypes (the combination of methylated and unmethylated CpG sites along the same DNA molecule) showed a significantly lower intra-individual diversity in the DS group, which at the same time was characterized by a higher interindividual variability. Finally, we showed that RNA45S expression is lower in adults with DS. Collectively, our results suggest a rearrangement of the epigenetic profile of rDNA in DS, possibly to compensate for the extranumerary rDNA copies. Future studies should assess whether the regulation of ribosome biogenesis can contribute to the pathogenesis of DS and explain the clinical heterogeneity characteristic of the syndrome.

Dietary, Cultural, and Pathogens-Related Selective Pressures Shaped Differential Adaptive Evolution among Native Mexican Populations.

Native American genetic ancestry has been remarkably implicated with increased risk of diverse health issues in several Mexican populations, especially in relation to the dramatic changes in environmental, dietary, and cultural settings they have recently undergone. In particular, the effects of these ecological transitions and Westernization of lifestyles have been investigated so far predominantly on Mestizo individuals. Nevertheless, indigenous groups, rather than admixed Mexicans, have plausibly retained the highest proportions of genetic components shaped by natural selection in response to the ancient milieu experienced by Mexican ancestors during their pre-Columbian evolutionary history. These formerly adaptive variants have the potential to represent the genetic determinants of some biological traits that are peculiar to Mexican people, as well as a reservoir of loci with possible biomedical relevance. To test such a hypothesis, we used genome-wide genotype data to infer the unique adaptive evolution of Native Mexican groups selected as reasonable descendants of the main pre-Columbian Mexican civilizations. A combination of haplotype-based and gene-network analyses enabled us to detect genomic signatures ascribable to polygenic adaptive traits plausibly evolved by the main genetic clusters of Mexican indigenous populations to cope with local environmental and/or cultural conditions. Some of these adaptations were found to play a role in modulating the susceptibility/resistance of these groups to certain pathological conditions, thus providing new evidence that diverse selective pressures have contributed to shape the current biological and disease-risk patterns of present-day Native and Mestizo Mexican populations.

Association of rs3027178 polymorphism in the circadian clock gene PER1 with susceptibility to Alzheimer's disease and longevity in an Italian population.

Many physiological processes in the human body follow a 24-h circadian rhythm controlled by the circadian clock system. Light, sensed by retina, is the predominant "zeitgeber" able to synchronize the circadian rhythms to the light-dark cycles. Circadian rhythm dysfunction and sleep disorders have been associated with aging and neurodegenerative diseases including mild cognitive impairment (MCI) and Alzheimer's disease (AD). In the present study, we aimed at investigating the genetic variability of clock genes in AD patients compared to healthy controls from Italy. We also included a group of Italian centenarians, considered as super-controls in association studies given their extreme phenotype of successful aging. We analyzed the exon sequences of eighty-four genes related to circadian rhythms, and the most significant variants identified in this first discovery phase were further assessed in a larger independent cohort of AD patients by matrix assisted laser desorption/ionization-time of flight mass spectrometry. The results identified a significant association between the rs3027178 polymorphism in the PER1 circadian gene with AD, the G allele being protective for AD. Interestingly, rs3027178 showed similar genotypic frequencies among AD patients and centenarians. These results collectively underline the relevance of circadian dysfunction in the predisposition to AD and contribute to the discussion on the role of the relationship between the genetics of age-related diseases and of longevity.

CNA Profiling of Single CTCs in Locally Advanced Esophageal Cancer Patients during Therapy Highlights Unexplored Molecular Pathways.

BACKGROUND: Here, we monitored the evolution of CTCs spread in 11 patients affected by locally advanced EC who were undergoing therapy. METHODS: In this perspective study, we designed multiple blood biopsies from individual patients: before and after neoadjuvant chemo-radio therapy and after surgery. We developed a multi-target array, named Grab-all assay, to estimate CTCs for their epithelial (EpCAM/E-Cadherin/Cytokeratins) and mesenchymal/stem (N-Cadherin/CD44v6/ABCG2) phenotypes. Identified CTCs were isolated as single cells by DEPArray, subjected to whole genome amplification, and copy number aberration (CNA) profiles were determined. Through bioinformatic analysis, we assessed the genomic imbalance of single CTCs, investigated specific focal copy number changes previously reported in EC and aberrant pathways using enrichment analysis. RESULTS: Longitudinal monitoring allowed the identification of CTCs in at least one time-point per patient. Through single cell CNA analysis, we revealed that CTCs showed significantly dynamic genomic imbalance during treatment. Individual CTCs from relapsed patients displayed a higher degree of genomic imbalance relative to disease-free patients' groups. Genomic aberrations previously reported in EC occurred mostly in post-neoadjuvant therapy CTCs. In-depth analysis showed that networks enrichment in all time-point CTCs were inherent to innate immune system. Transcription/gene regulation, post-transcriptional and epigenetic modifications were uniquely affected in CTCs of relapsed patients. CONCLUSIONS: Our data add clues to the comprehension of the role of CTCs in EC aggressiveness: chromosomal aberrations on genes related to innate immune system behave as relevant to the onset of CTC-status, whilst pathways of transcription/gene regulation, post-transcriptional and epigenetic modifications seem linked to patients' outcome.

Internal validation and improvement of mitochondrial genome sequencing using the Precision ID mtDNA Whole Genome Panel.

With the recent advances in next-generation sequencing (NGS), mitochondrial whole-genome sequencing has begun to be applied to the field of the forensic biology as an alternative to the traditional Sanger-type sequencing (STS). However, experimental workflows, commercial solutions, and output data analysis must be strictly validated before being implemented into the forensic laboratory. In this study, we performed an internal validation for an NGS-based typing of the entire mitochondrial genome using the Precision ID mtDNA Whole Genome Panel (Thermo Fisher Scientific) on the Ion S5 sequencer (Thermo Fisher Scientific). Concordance, repeatability, reproducibility, sensitivity, and heteroplasmy detection analyses were assessed using the 2800 M and 9947A standard control DNA as well as typical casework specimens, and results were compared with conventional Sanger sequencing and another NGS sequencer in a different laboratory. We discuss the strengths and limitations of this approach, highlighting some issues regarding noise thresholds and heteroplasmy detection, and suggesting solutions to mitigate these effects and improve overall data interpretation. Results confirmed that the Precision ID Whole mtDNA Genome Panel is highly reproducible and sensitive, yielding useful full mitochondrial DNA sequences also from challenging DNA specimens, thus providing further support for its use in forensic practice.

Insights into Punic genetic signatures in the southern necropolis of Tharros (Sardinia).

BACKGROUND: Phoenician and Punic expansions have been protagonists of intense trade networks and settlements in the Mediterranean Sea. AIMS: The maternal genetic variability of ancient Punic samples from the Sardinian necropolis of Tharros was analysed, with the aim to explore genetic interactions and signatures of past population events. SUBJECTS AND METHODS: The mtDNA HVS-I and coding region SNPs were analysed in 14 Punic samples and 74 modern individuals from Cabras and Belvì (for which the HVS-II region was also analysed). The results were compared with 5,590 modern Euro-Mediterranean sequences and 127 ancient samples. RESULTS: While contemporary groups fall within the genetic variability of other modern Sardinians, our Punic samples reveal proximity to present-day North-African and Iberian populations. Furthermore, Cabras and Belvì cluster mainly with pre-Phoenician groups, while samples from Tharros project with other Punic Sardinian individuals. CONCLUSION: This study provides the first preliminary insights into the population dynamics of the Punic site of Tharros. While the number of currently available samples does not allow definitive investigation of the connection with indigenous Sardinian groups, our results seem to confirm internal migratory phenomena in the central-western Mediterranean and female participation in the Punic mobility.

Uniparental Lineages from the Oldest Indigenous Population of Ecuador: The Tsachilas.

Together with Cayapas, the Tsachilas constitute the oldest population in the country of Ecuador and, according to some historians, they are the last descendants of the ancient Yumbos. Several anthropological issues underlie the interest towards this peculiar population: the uncertainty of their origin, their belonging to the Barbacoan linguistic family, which is still at the center of an intense linguistic debate, and the relations of their Yumbo ancestors with the Inca invaders who occupied their ancient territory. Our contribution to the knowledge of their complex past was the reconstruction of their genetic maternal and paternal inheritance through the sequencing of 70 entire mitochondrial genomes and the characterization of the non-recombinant region of the Y chromosome in 26 males. For both markers, we built comprehensive datasets of various populations from the surrounding geographical area, northwestern South America, NW, with a known linguistic affiliation, and we could then compare our sample against the overall variability to infer relationships with other Barbacoan people and with other NW natives. We found contrasting patterns of genetic diversity for the two markers, but generally, our results indicated a possible common origin between the Tsachilas, the Chachi, and other Ecuadorian and Colombian Barbacoans and are suggestive of an interesting ancient linkage to the Inca invaders in Yumbo country.

Detecting Variants in the NBN Gene While Testing for Hereditary Breast Cancer: What to Do Next?

The NBN gene has been included in breast cancer (BC) multigene panels based on early studies suggesting an increased BC risk for carriers, though not confirmed by recent research. To evaluate the impact of NBN analysis, we assessed the results of NBN sequencing in 116 BRCA-negative BC patients and reviewed the literature. Three patients (2.6%) carried potentially relevant variants: two, apparently unrelated, carried the frameshift variant c.156_157delTT and another one the c.628G>T variant. The latter was subsequently found in 4/1390 (0.3%) BC cases and 8/1580 (0.5%) controls in an independent sample, which, together with in silico predictions, provided evidence against its pathogenicity. Conversely, the rare c.156_157delTT variant was absent in the case-control set; moreover, a 50% reduction of NBN expression was demonstrated in one carrier. However, in one family it failed to co-segregate with BC, while the other carrier was found to harbor also a probably pathogenic TP53 variant that may explain her phenotype. Therefore, the c.156_157delTT, although functionally deleterious, was not supported as a cancer-predisposing defect. Pathogenic/likely pathogenic NBN variants were detected by multigene panels in 31/12314 (0.25%) patients included in 15 studies. The risk of misinterpretation of such findings is substantial and supports the exclusion of NBN from multigene panels.

Y-chromosome variability and genetic history of Commons from Northern Italy.

OBJECTIVES: Genetic drift and admixture are driving forces in human evolution, but their concerted impact to population evolution in historical times and at a micro-geographic scale is poorly assessed. In this study we test a demographic model encompassing both admixture and drift to the case of social-cultural isolates such as the so-called "Commons." MATERIALS AND METHODS: Commons are peculiar institutions of medieval origins whose key feature is the tight relationship between population and territory, mediated by the collective property of shared resources. Here, we analyze the Y-chromosomal genetic structure of four Commons (for a total of 366 samples) from the Central and Eastern Padana plain in Northern Italy. RESULTS: Our results reveal that all these groups exhibit patterns of significant diversity reduction, peripheral/outlier position within the Italian/European genetic space and high frequency of Common-specific haplogroups. By explicitly testing different drift-admixture models, we show that a drift-only model is more probable for Central Padana Commons, while additional admixture (~20%) from external population around the same time of their foundation cannot be excluded for the Eastern ones. DISCUSSION: Building on these results, we suggest central Middle Ages as the most probable age of foundation for three of the considered Commons, the remaining one pointing to late antiquity. We conclude that an admixture-drift model is particularly useful for interpreting the genetic structure and recent demographic history of small-scale populations in which social-cultural features play a significant role.

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

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

Genomic adaptations to cereal-based diets contribute to mitigate metabolic risk in some human populations of East Asian ancestry.

Adoption of diets based on some cereals, especially on rice, signified an iconic change in nutritional habits for many Asian populations and a relevant challenge for their capability to maintain glucose homeostasis. Indeed, rice shows the highest carbohydrates content and glycemic index among the domesticated cereals and its usual ingestion represents a potential risk factor for developing insulin resistance and related metabolic diseases. Nevertheless, type 2 diabetes and obesity epidemiological patterns differ among Asian populations that rely on rice as a staple food, with higher diabetes prevalence and increased levels of central adiposity observed in people of South Asian ancestry rather than in East Asians. This may be at least partly due to the fact that populations from East Asian regions where wild rice or other cereals such as millet have been already consumed before their cultivation and/or were early domesticated have relied on these nutritional resources for a period long enough to have possibly evolved biological adaptations that counteract their detrimental side effects. To test such a hypothesis, we compared adaptive evolution of these populations with that of control groups from regions where the adoption of cereal-based diets occurred many thousand years later and which were identified from a genome-wide dataset including 2,379 individuals from 124 East Asian and South Asian populations. This revealed selective sweeps and polygenic adaptive mechanisms affecting functional pathways involved in fatty acids metabolism, cholesterol/triglycerides biosynthesis from carbohydrates, regulation of glucose homeostasis, and production of retinoic acid in Chinese Han and Tujia ethnic groups, as well as in people of Korean and Japanese ancestry. Accordingly, long-standing rice- and/or millet-based diets have possibly contributed to trigger the evolution of such biological adaptations, which might represent one of the factors that play a role in mitigating the metabolic risk of these East Asian populations.

Genetic history of Calabrian Greeks reveals ancient events and long term isolation in the Aspromonte area of Southern Italy.

Calabrian Greeks are an enigmatic population that have preserved and evolved a unique variety of language, Greco, survived in the isolated Aspromonte mountain area of Southern Italy. To understand their genetic ancestry and explore possible effects of geographic and cultural isolation, we genome-wide genotyped a large set of South Italian samples including both communities that still speak Greco nowadays and those that lost the use of this language earlier in time. Comparisons with modern and ancient populations highlighted ancient, long-lasting genetic links with Eastern Mediterranean and Caucasian/Near-Eastern groups as ancestral sources of Southern Italians. Our results suggest that the Aspromonte communities might be interpreted as genetically drifted remnants that departed from such ancient genetic background as a consequence of long-term isolation. Specific patterns of population structuring and higher levels of genetic drift were indeed observed in these populations, reflecting geographic isolation amplified by cultural differences in the groups that still conserve the Greco language. Isolation and drift also affected the current genetic differentiation at specific gene pathways, prompting for future genome-wide association studies aimed at exploring trait-related loci that have drifted up in frequency in these isolated groups.

Haplotype data and forensic evaluation of 23 Y-STR and 12 X-STR loci in eight ethnic groups from Eritrea.

Eritrea is a multi-ethnic country of over 3 million of people consisting of different ethnic groups, having each its own language and cultural tradition. Due to the lack of population genetic data for markers of forensic interest, in this study, we analyzed the genetic polymorphisms of 23 Y-chromosome STR loci and of 12 X-chromosome STR loci in a sample of 255 unrelated individuals from 8 Eritrean ethnic groups, with the aim to generate a reference haplotype database for anthropological and forensic applications. X- and Y-chromosomes markers may indeed offer information especially in personal identification and kinship testing, when relying on the availability of large local population data to derive sufficiently accurate frequency estimates. The population genetic analyses in the Eritrean sample for both the two set of Y- and X-STR markers showed high power of discrimination both at country-based and population levels. Comparison population results highlight the importance of considering the ethnic composition within the analyzed country and the necessity of increasing available data especially when referring to heterogeneous populations such as the African ones.

Population Dynamics in Italian Canids between the Late Pleistocene and Bronze Age.

Dog domestication is still largely unresolved due to time-gaps in the sampling of regions. Ancient Italian canids are particularly understudied, currently represented by only a few specimens. In the present study, we sampled 27 canid remains from Northern Italy dated between the Late Pleistocene and Bronze Age to assess their genetic variability, and thus add context to dog domestication dynamics. They were targeted at four DNA fragments of the hypervariable region 1 of mitochondrial DNA. A total of 11 samples had good DNA preservation and were used for phylogenetic analyses. The dog samples were assigned to dog haplogroups A, C and D, and a Late Pleistocene wolf was set into wolf haplogroup 2. We present our data in the landscape of ancient and modern dog genetic variability, with a particular focus on the ancient Italian samples published thus far. Our results suggest there is high genetic variability within ancient Italian canids, where close relationships were evident between both a ~24,700 years old Italian canid, and Iberian and Bulgarian ancient dogs. These findings emphasize that disentangling dog domestication dynamics benefits from the analysis of specimens from Southern European regions.

Genomic history of the Italian population recapitulates key evolutionary dynamics of both Continental and Southern Europeans.

BACKGROUND: The cline of human genetic diversity observable across Europe is recapitulated at a micro-geographic scale by variation within the Italian population. Besides resulting from extensive gene flow, this might be ascribable also to local adaptations to diverse ecological contexts evolved by people who anciently spread along the Italian Peninsula. Dissecting the evolutionary history of the ancestors of present-day Italians may thus improve the understanding of demographic and biological processes that contributed to shape the gene pool of European populations. However, previous SNP array-based studies failed to investigate the full spectrum of Italian variation, generally neglecting low-frequency genetic variants and examining a limited set of small effect size alleles, which may represent important determinants of population structure and complex adaptive traits. To overcome these issues, we analyzed 38 high-coverage whole-genome sequences representative of population clusters at the opposite ends of the cline of Italian variation, along with a large panel of modern and ancient Euro-Mediterranean genomes. RESULTS: We provided evidence for the early divergence of Italian groups dating back to the Late Glacial and for Neolithic and distinct Bronze Age migrations having further differentiated their gene pools. We inferred adaptive evolution at insulin-related loci in people from Italian regions with a temperate climate, while possible adaptations to pathogens and ultraviolet radiation were observed in Mediterranean Italians. Some of these adaptive events may also have secondarily modulated population disease or longevity predisposition. CONCLUSIONS: We disentangled the contribution of multiple migratory and adaptive events in shaping the heterogeneous Italian genomic background, which exemplify population dynamics and gene-environment interactions that played significant roles also in the formation of the Continental and Southern European genomic landscapes.

Gene Panel Analysis in a Large Cohort of Patients With Autosomal Dominant Polycystic Kidney Disease Allows the Identification of 80 Potentially Causative Novel Variants and the Characterization of a Complex Genetic Architecture in a Subset of Families.

Introduction: Autosomal dominant polycystic kidney disease (ADPKD) is one of the most common inherited disorders in humans and the majority of patients carry a variant in either PKD1 or PKD2. Genetic testing is increasingly required for diagnosis, prognosis, and treatment decision, but it is challenging due to segmental duplications of PKD1, genetic and allelic heterogeneity, and the presence of many variants hypomorphic or of uncertain significance. We propose an NGS-based testing strategy for molecular analysis of ADPKD and its phenocopies, validated in a diagnostic setting. Materials and Methods: Our protocol is based on high-throughput simultaneous sequencing of PKD1 and PKD2 after long range PCR of coding regions, followed by a masked reference genome alignment, and MLPA analysis. A further screening of additional 14 cystogenes was performed in negative cases. We applied this strategy to analyze 212 patients with a clinical suspicion of ADPKD. Results and Discussion: We detected causative variants (interpreted as pathogenic/likely pathogenic) in 61.3% of our index patients, and variants of uncertain clinical significance in 12.5%. The majority (88%) of genetic variants was identified in PKD1, 12% in PKD2. Among 158 distinct variants, 80 (50.6%) were previously unreported, confirming broad allelic heterogeneity. Eleven patients showed more than one variant. Segregation analysis indicated biallelic disease in five patients, digenic in one, de novo variant with unknown phase in two. Furthermore, our NGS protocol allowed the identification of two patients with somatic mosaicism, which was undetectable with Sanger sequencing. Among patients without PKD1/PKD2 variants, we identified three with possible alternative diagnosis: a patient with biallelic mutations in PKHD1, confirming the overlap between recessive and dominant PKD, and two patients with variants in ALG8 and PRKCSH, respectively. Genotype-phenotype correlations showed that patients with PKD1 variants predicted to truncate (T) the protein experienced end-stage renal disease 9 years earlier than patients with PKD1 non-truncating (NT) mutations and >13 years earlier than patients with PKD2 mutations. ADPKD-PKD1 T cases showed a disease onset significantly earlier than ADPKD-PKD1 NT and ADPK-PKD2, as well as a significant earlier diagnosis. These data emphasize the need to combine clinical information with genetic data to achieve useful prognostic predictions.

Gut microbiota composition in Himalayan and Andean populations and its relationship with diet, lifestyle and adaptation to the high-altitude environment.

Human populations living at high altitude evolved a number of biological adjustments to cope with a challenging environment characterised especially by reduced oxygen availability and limited nutritional resources. This condition may also affect their gut microbiota composition. Here, we explored the impact of exposure to such selective pressures on human gut microbiota by considering different ethnic groups living at variable degrees of altitude: the high-altitude Sherpa and low-altitude Tamang populations from Nepal, the high-altitude Aymara population from Bolivia, as well as a low-altitude cohort of European ancestry, used as control. We thus observed microbial profiles common to the Sherpa and Aymara, but absent in the low-altitude cohorts, which may contribute to the achievement of adaptation to high-altitude lifestyle and nutritional conditions. The collected evidences suggest that microbial signatures associated to these rural populations may enhance metabolic functions able to supply essential compounds useful for the host to cope with high altitude-related physiological changes and energy demand. Therefore, these results add another valuable piece of the puzzle to the understanding of the beneficial effects of symbiosis between microbes and their human host even from an evolutionary perspective.

Estimating Y-Str Mutation Rates and Tmrca Through Deep-Rooting Italian Pedigrees.

In the population genomics era, the study of Y-chromosome variability is still of the greatest interest for several fields ranging from molecular anthropology to forensics and genetic genealogy. In particular, mutation rates of Y-chromosomal Short Tandem Repeats markers (Y-STRs) are key parameters for different interdisciplinary applications. Among them, testing the patrilineal relatedness between individuals and calculating their Time of Most Recent Common Ancestors (TMRCAs) are of the utmost importance. To provide new valuable estimates and to address these issues, we typed 47 Y-STRs (comprising Yfiler, PowerPlex23 and YfilerPlus loci, the recently defined Rapidly Mutating [RM] panel and 11 additional markers often used in genetic genealogical applications) in 135 individuals belonging to 66 deep-rooting paternal genealogies from Northern Italy. Our results confirmed that the genealogy approach is an effective way to obtain reliable Y-STR mutation rate estimates even with a limited number of samples. Moreover, they showed that the impact of multi-step mutations and backmutations is negligible within the temporal scale usually adopted by forensic and genetic genealogy analyses. We then detected a significant association between the number of mutations within genealogies and observed TMRCAs. Therefore, we compared observed and expected TMRCAs by implementing a Bayesian procedure originally designed by Walsh (2001) and showed that the method yields a good performance (up to 96.72%), especially when using the Infinite Alleles Model (IAM).