The Persian plateau served as hub for Homo sapiens after the main out of Africa dispersal.

A combination of evidence, based on genetic, fossil and archaeological findings, indicates that Homo sapiens spread out of Africa between ~70-60 thousand years ago (kya). However, it appears that once outside of Africa, human populations did not expand across all of Eurasia until ~45 kya. The geographic whereabouts of these early settlers in the timeframe between ~70-60 to 45 kya has been difficult to reconcile. Here we combine genetic evidence and palaeoecological models to infer the geographic location that acted as the Hub for our species during the early phases of colonisation of Eurasia. Leveraging on available genomic evidence we show that populations from the Persian Plateau carry an ancestry component that closely matches the population that settled the Hub outside Africa. With the paleoclimatic data available to date, we built ecological models showing that the Persian Plateau was suitable for human occupation and that it could sustain a larger population compared to other West Asian regions, strengthening this claim.

Aberrantly Expressed Hsa_circ_0060762 and CSE1L as Potential Peripheral Blood Biomarkers for ALS.

Amyotrophic lateral sclerosis (ALS) is a rapidly progressive adult-onset neurodegenerative disease that is often diagnosed with a delay due to initial non-specific symptoms. Therefore, reliable and easy-to-obtain biomarkers are an absolute necessity for earlier and more accurate diagnostics. Circular RNAs (circRNAs) have already been proposed as potential biomarkers for several neurodegenerative diseases. In this study, we further investigated the usefulness of circRNAs as potential biomarkers for ALS. We first performed a microarray analysis of circRNAs on peripheral blood mononuclear cells of a subset of ALS patients and controls. Among the differently expressed circRNA by microarray analysis, we selected only the ones with a host gene that harbors the highest level of conservation and genetic constraints. This selection was based on the hypothesis that genes under selective pressure and genetic constraints could have a major role in determining a trait or disease. Then we performed a linear regression between ALS cases and controls using each circRNA as a predictor variable. With a False Discovery Rate (FDR) threshold of 0.1, only six circRNAs passed the filtering and only one of them remained statistically significant after Bonferroni correction: hsa_circ_0060762 and its host gene CSE1L. Finally, we observed a significant difference in expression levels between larger sets of patients and healthy controls for both hsa_circ_0060762 and CSE1L. CSE1L is a member of the importin ß family and mediates inhibition of TDP-43 aggregation; the central pathogenicity in ALS and hsa_circ_0060762 has binding sites for several miRNAs that have been already proposed as biomarkers for ALS. In addition, receiver operating characteristics curve analysis showed diagnostic potential for CSE1L and hsa_circ_0060762. Hsa_circ_0060762 and CSE1L thus represent novel potential peripheral blood biomarkers and therapeutic targets for ALS.

Insights into gene tissue specificity and protein-protein interactions in the context of purifying selection in humans.

BACKGROUND: How much are natural selection and gene characteristics, such as the number of protein-protein interactions (PPIs), tissue specificity (𝞽), and expression level, connected? METHODS: In order to investigate these relationships, we combined different metrics linked to genetic constraints and analyzed their distribution concerning PPIs, 𝞽 and expression levels. RESULTS: We discovered a positive correlation between genetic constraints, PPIs, and expression levels in all tissues. On the other hand, we obtained a negative correlation between genetic constraints and 𝞽. Furthermore, the fraction of variance in PPI and 𝞽 explained by the constraints metrics is around 6% and 10%, respectively. CONCLUSIONS: We observed that the variance of expression of tissue-specific genes seems not related to their level of selection constraints, which is the opposite of what is found on non-tissue-specific genes. Overall these observations would help to elucidate the relationship between natural selection and gene features.

Impact of cultural and genetic structure on food choices along the Silk Road.

The complex interplay between genetics, culture, and environment forms an individual's biology, influencing their behavior, choices, and health. However, to what extent information derived from this intertwined network could be quantitatively summarized to provide a glance at an individual's lifestyle is difficult to say. Here, we focused on dietary preferences as cultural proxies and genome-wide data of 543 individuals from six historical Silk Road countries: Georgia, Armenia, Azerbaijan, Uzbekistan, Kazakhstan, and Tajikistan. These lands favored the dispersal of innovations, foods, and DNA halfway across Eurasia, thus representing an ideal subject to explore interactions of cultural factors and genetic ancestry. We used discriminant analysis of principal components to infer cultural clusters, where mixed memberships are allowed. Five different clusters emerged. Of these, clusters 1 and 3, driven by aversion to pork and alcoholic beverages, mirrored genetic admixture patterns with the exception of Azerbaijan, which shares preferences supported by Islamic culture with Eastern countries. Cluster 3 was driven by protein-rich foods, whose preference was significantly related to steppe pastoralist ancestry. Sex and age were secondary clustering factors, with clusters formed by male and young individuals being related to alcohol preference and a reduced liking for vegetables. The soft clustering approach enabled us to model and summarize the individual's dietary information in short and informative vectors, which show meaningful interaction with other nondietary attributes of the studied individuals. Encoding other cultural variables would help summarize an individual's culture quantitatively, thus ultimately supporting its inclusion as a covariate in future association studies.

Ancestry-related distribution of Runs of homozygosity and functional variants in Qatari population.

BACKGROUND: Describing how genetic history shapes the pattern of medically relevant variants could improve the understanding of how specific loci interact with each other and affect diseases and traits prevalence. The Qatari population is characterized by a complex history of admixture and substructure, and the study of its population genomic features would provide valuable insights into the genetic landscape of functional variants. Here, we analyzed the genomic variation of 186 newly-genotyped healthy individuals from the Qatari peninsula. RESULTS: We discovered an intricate genetic structure using ancestry related analyses. In particular, the presence of three different clusters, Cluster 1, Cluster 2 and Cluster 3 (with Near Eastern, South Asian and African ancestry, respectively), was detected with an additional fourth one (Cluster 4) with East Asian ancestry. These subpopulations show differences in the distribution of runs of homozygosity (ROH) and admixture events in the past, ranging from 40 to 5 generations ago. This complex genetic history led to a peculiar pattern of functional markers under positive selection, differentiated in shared signals and private signals. Interestingly we found several signatures of shared selection on SNPs in the FADS2 gene, hinting at a possible common evolutionary link to dietary intake. Among the private signals, we found enrichment for markers associated with HDL and LDL for Cluster 1(Near Eastern ancestry) and Cluster 3 (South Asian ancestry) and height and blood traits for Cluster 2 (African ancestry). The differences in genetic history among these populations also resulted in the different frequency distribution of putative loss of function variants. For example, homozygous carriers for rs2884737, a variant linked to an anticoagulant drug (warfarin) response, are mainly represented by individuals with predominant Bedouin ancestry (risk allele frequency G at 0.48). CONCLUSIONS: We provided a detailed catalogue of the different ancestral pattern in the Qatari population highlighting differences and similarities in the distribution of selected variants and putative loss of functions. Finally, these results would provide useful guidance for assessing genetic risk factors linked to consanguinity and genetic ancestry.

Fine-scale population structure and demographic history of British Pakistanis.

Previous genetic and public health research in the Pakistani population has focused on the role of consanguinity in increasing recessive disease risk, but little is known about its recent population history or the effects of endogamy. Here, we investigate fine-scale population structure, history and consanguinity patterns using genotype chip data from 2,200 British Pakistanis. We reveal strong recent population structure driven by the biraderi social stratification system. We find that all subgroups have had low recent effective population sizes (Ne), with some showing a decrease 15‒20 generations ago that has resulted in extensive identity-by-descent sharing and homozygosity, increasing the risk of recessive disorders. Our results from two orthogonal methods (one using machine learning and the other coalescent-based) suggest that the detailed reporting of parental relatedness for mothers in the cohort under-represents the true levels of consanguinity. These results demonstrate the impact of cultural practices on population structure and genomic diversity in Pakistanis, and have important implications for medical genetic studies.

Poking COVID-19: Insights on Genomic Constraints among Immune-Related Genes between Qatari and Italian Populations.

Host genomic information, specifically genomic variations, may characterize susceptibility to disease and identify people with a higher risk of harm, leading to better targeting of care and vaccination. Italy was the epicentre for the spread of COVID-19 in Europe, the first country to go into a national lockdown and has one of the highest COVID-19 associated mortality rates. Qatar, on the other hand has a very low mortality rate. In this study, we compared whole-genome sequencing data of 14398 adults and Qatari-national to 925 Italian individuals. We also included in the comparison whole-exome sequence data from 189 Italian laboratory-confirmed COVID-19 cases. We focused our study on a curated list of 3619 candidate genes involved in innate immunity and host-pathogen interaction. Two population-gene metric scores, the Delta Singleton-Cohort variant score (DSC) and Sum Singleton-Cohort variant score (SSC), were applied to estimate the presence of selective constraints in the Qatari population and in the Italian cohorts. Results based on DSC and SSC metrics demonstrated a different selective pressure on three genes (MUC5AC, ABCA7, FLNA) between Qatari and Italian populations. This study highlighted the genetic differences between Qatari and Italian populations and identified a subset of genes involved in innate immunity and host-pathogen interaction.

Natural human knockouts and Mendelian disorders: deep phenotyping in Italian isolates.

Whole genome sequencing (WGS) allows the identification of human knockouts (HKOs), individuals in whom loss of function (LoF) variants disrupt both alleles of a given gene. HKOs are a valuable model for understanding the consequences of genes function loss. Naturally occurring biallelic LoF variants tend to be significantly enriched in "genetic isolates," making these populations specifically suited for HKO studies. In this work, a meticulous WGS data analysis combined with an in-depth phenotypic assessment of 947 individuals from three Italian genetic isolates led to the identification of ten biallelic LoF variants in ten OMIM genes associated with known autosomal recessive diseases. Notably, only a minority of the identified HKOs (C7, F12, and GPR68 genes) displayed the expected phenotype. For most of the genes, instead, (ACADSB, FANCL, GRK1, LGI4, MPO, PGAM2, and RP1L1), the carriers showed none or few of the signs and symptoms typically associated with the related diseases. Of particular interest is a case presenting with a FANCL biallelic LoF variant and a positive diepoxybutane test but lacking a full Fanconi anemia phenotypic spectrum. Identifying KO subjects displaying expected phenotypes suggests that the lack of correct genetic diagnoses may lead to inappropriate and delayed treatment. In contrast, the presence of HKOs with phenotypes deviating from the expected patterns underlines how LoF variants may be responsible for broader phenotypic spectra. Overall, these results highlight the importance of in-depth phenotypical characterization to understand the role of LoF variants and the advantage of studying these variants in genetic isolates.

Runs of homozygosity are associated with staging of periodontitis in isolated populations.

Periodontitis is a common inflammatory disease characterized by a complex etiology, which is the result of a combination of genetic and environmental factors. Genetic variants linked to the periodontitis disease were already investigated, however, little was known regarding the severity of this disease. Recently, long runs of homozygosity (ROH) were associated with several multifactorial diseases. Therefore, in our work, we tried to assess the role of ROH and periodontitis status. We found an association between the excess of homozygosity owing to ROH and staging of periodontitis. More in detail, the total amount of homozygosity owing to ROH is positively associated with an increased severity of periodontitis (P = 0.0001). Regression tree analysis showed the impact of ROH burden in discriminating individuals with mild periodontitis stages I and II and periodontitis stages III and IV (P < 0.001). Furthermore, ROH mapping highlights several regions associated with a severe status of periodontitis (odds ratio > 1). Among them, we found a total of 33 genes. Interestingly, some of these genes were previously associated with granulocyte or platelet measures, both linked to the onset and the progression of periodontal disease. Our results suggest the not only single variants association test could help to risk assessment but even individual genomic features; furthermore, our ROH mapping highlighted the possible role of multiple genes in periodontal development.

Glucagon-like peptide-1 receptor and sarcoglycan delta genetic variants can affect cardiovascular risk in chronic kidney disease patients under hemodialysis.

BACKGROUND: Chronic kidney disease (CKD) patients under hemodialysis show a higher risk of cardiovascular (CV) mortality and morbidity than the general population. This study aims to identify genetic markers that could explain the increased CV risk in hemodialysis. METHODS: A total of 245 CKD patients under hemodialysis were recruited and followed up for 5 years to record CV events. Genetic analysis was performed using single-nucleotide polymorphisms (SNPs) genotyping by Infinium Expanded Multi-Ethnic Genotyping Array (Illumina, San Diego, CA, USA) comparing patients with and without a history of CV events [161 cardiovascular diseases (CVDs) and 84 no CVDs]. The fixation index (Fst) measure was used to identify the most differentiated SNPs, and gene ontology analysis [Protein Analysis THrough Evolutionary Relationships (PANTHER) and Ingenuity Pathway Analysis (IPA)] was applied to define the biological/pathological roles of the associated SNPs. Partitioning tree analysis interrogated the genotype-phenotype relationship between discovered genetic variants and CV phenotypes. Cox regression analysis measured the effect of these SNPs on new CV events during the follow-up (FU). RESULTS: Fst analysis identified 3218 SNPs that were significantly different between CVD and no CVD. Gene ontology analysis identified two of these SNPs as involved in cardiovascular disease pathways (Ingenuity Pathway) and heart development (Panther) and belonging to 2 different genes: Glucagon-like peptide-1 receptor (GLP1R) and Sarcoglycan delta (SGCD). The phenotype-genotype analysis found a higher percentage of CVD patients carrying the GLP1R rs10305445 allele A (P = 0.03) and lower percentages of CVD patients carrying the SGCD rs145292439 allele A (P = 0.038). Moreover, SGCD rs145292439 was associated with higher levels of high-density lipoprotein (P = 0.015). Cox analysis confirmed the increased frequency of CV events during the 5-year FU in patients carrying GLP1R rs1035445 allele A but it did not show any significant association with SGCD rs145292439. CONCLUSIONS: This study identified GLP1R rs10305445 and SCGD rs145292439 as potential genetic markers that may explain the higher risk of CVD in hemodialysis patients.

A population-based approach for gene prioritization in understanding complex traits.

Gene prioritization is the process of determining which variants and genes identified in genetic analyses are likely to cause a disease or a variation in a phenotype. For many genes, neither in vitro nor in vivo testing is available, thus assessing their pathogenic role could be challenging, leading to false-positive or false-negative results. In this paper, we propose an innovative score of gene prioritization based on the population of interest. We introduce the concept of singleton-cohort variants (SC variant), a variant that has allele count equal to one in the cohort under study. The difference between the normalized count of SC variants in the coding region and the normalized count of SC variants in the non-coding region should give a hint regarding the level of constraints for that gene in a specific population. This scoring system is negative when there are constraints that allow the presence of SC variants only in the non-coding region; on the contrary, it is positive when there are no constraints. A complimentary score is the sum of SC variants normalized count in both coding and non-coding regions, which could be used as a proxy of positive or strong purifying selection in a specific population. Our methodology showed a high level of constraining for genes such as USP34 in all subpopulations tested (1000 G dataset). In contrast, some genes showed a high negative score only in specific populations, e.g., MYT1L in Europeans, UBR5 in East Asians, and FBXO11 in Africans.

Author Correction: Establishment and equilibrium levels of deleterious mutations in large populations.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

A bird's-eye view of Italian genomic variation through whole-genome sequencing.

The genomic variation of the Italian peninsula populations is currently under characterised: the only Italian whole-genome reference is represented by the Tuscans from the 1000 Genome Project. To address this issue, we sequenced a total of 947 Italian samples from three different geographical areas. First, we defined a new Italian Genome Reference Panel (IGRP1.0) for imputation, which improved imputation accuracy, especially for rare variants, and we tested it by GWAS analysis on red blood traits. Furthermore, we extended the catalogue of genetic variation investigating the level of population structure, the pattern of natural selection, the distribution of deleterious variants and occurrence of human knockouts (HKOs). Overall the results demonstrate a high level of genomic differentiation between cohorts, different signatures of natural selection and a distinctive distribution of deleterious variants and HKOs, confirming the necessity of distinct genome references for the Italian population.

A transcriptomic atlas of mammalian olfactory mucosae reveals an evolutionary influence on food odor detection in humans.

The mammalian olfactory system displays species-specific adaptations to different ecological niches. To investigate the evolutionary dynamics of olfactory sensory neuron (OSN) subtypes across mammalian evolution, we applied RNA sequencing of whole olfactory mucosa samples from mouse, rat, dog, marmoset, macaque, and human. We find that OSN subtypes, representative of all known mouse chemosensory receptor gene families, are present in all analyzed species. Further, we show that OSN subtypes expressing canonical olfactory receptors are distributed across a large dynamic range and that homologous subtypes can be either highly abundant across all species or species/order specific. Highly abundant mouse and human OSN subtypes detect odorants with similar sensory profiles and sense ecologically relevant odorants, such as mouse semiochemicals or human key food odorants. Together, our results allow for a better understanding of the evolution of mammalian olfaction in mammals and provide insights into the possible functions of highly abundant OSN subtypes.

Next generation sequencing study in a cohort of Italian patients with syndromic hearing loss.

Hearing loss (HL), one of the most common congenital disorder, affects about one child in 1000. Among the genetic forms of HL, ∼30% of the cases are associated with other signs or symptoms, leading to Syndromic Hearing Loss (SHL) with about 700 different forms described so far. In this report, we refer the clinical and molecular data of 38 Italian SHL unrelated patients, and their relatives, affected by the most common syndromes associated with HL (i.e., Usher, Pendred, Charge, Waardenburg, Alport, Stickler, Branchiootorenal and Microdeletions syndromes). Patients have been analysed using next-generation sequencing (NGS) and High Density (HD)-SNP array technologies. Data analysis led to the identification of nine novel and 27 known causative mutations in 12 genes and two microdeletions in chromosomes 1 and 10, respectively. In particular, as regards to Usher syndrome, that affects 32% of our patients, we were able to reach a molecular diagnosis in 83% of the cases and to identify in Northern Eastern Italy a very common USH2A gene mutation (39%) (c.11864G > A, p.(Trp3955*) which can be defined "Central-Eastern European allele." As regards to Alport syndrome, we were able to potentially reclassify a pathogenic allele in the COL4A3 gene, previously associated only with benign familial hematuria. In all the other cases, the genomic analysis allowed us to confirm the role of known causative genes and to identify several novel and known alleles. Overall, our results highlight the effectiveness of combining an accurate clinical characterization with the use of genomic technologies (NGS and SNP arrays) for the molecular diagnosis of SHL, with a clear positive impact in the management and treatment of all the patients.

Mutations in PLS1, encoding fimbrin, cause autosomal dominant nonsyndromic hearing loss.

Nonsyndromic hearing loss (NSHL), a common sensory disorder, is characterized by high clinical and genetic heterogeneity (i.e., approximately 115 genes and 170 loci so far identified). Nevertheless, almost half of patients submitted for genetic testing fail to receive a conclusive molecular diagnosis. We used next-generation sequencing to identify causal variants in PLS1 (c.805G>A, p.[E269K]; c.713G>T, p.[L238R], and c.383T>C, p.[F128S]) in three unrelated families of European ancestry with autosomal dominant NSHL. PLS1 encodes Plastin 1 (also called fimbrin), one of the most abundant actin-bundling proteins of the stereocilia. In silico protein modeling suggests that all variants destabilize the structure of the actin-binding domain 1, likely reducing the protein's ability to bind F actin. The role of PLS1 gene in hearing function is further supported by the recent demonstration that Pls1-/- mice show a hearing loss phenotype similar to that of our patients. In summary, we report PLS1 as a novel gene for autosomal dominant NSHL, suggesting that this gene is required for normal hearing in humans and mice.

Establishment and equilibrium levels of deleterious mutations in large populations.

Analytical and statistical stochastic approaches are used to model the dispersion of monogenic variants through large populations. These approaches are used to quantify the magnitude of the selective advantage of a monogenic heterozygous variant in the presence of a homozygous disadvantage. Dunbar's results regarding the cognitive upper limit of the number of stable social relationships that humans can maintain are used to determine a realistic effective community size from which an individual can select mates. By envisaging human community structure as a network where social proximity rather than physical geography predominates, a significant simplification is achieved, implicitly accounting for the effects of migration and consanguinity, and with population structure and genetic drift becoming emergent features of the model. Effective community size has a dramatic effect on the probability of establishing beneficial alleles. It also affects the eventual equilibrium values that are reached in the case of variants conferring a heterozygous selective advantage, but a homozygous disadvantage, as in the case of cystic fibrosis and sickle cell disease. The magnitude of this selective advantage can then be estimated based on observed occurrence levels of a specific allele in a population, without requiring prior information regarding its phenotypic manifestation.

Investigation of the link between PROP taste perception and vegetables consumption using FAOSTAT data.

In this work we investigated, in populations located in Central Asia, the relationship between PROP taste perception and vegetables liking and consumption using FAOSTAT dataset. Collected data were analysed using distance matrices, Mantel test and Pearson correlation. Populations showing similar ability in tasting PROP bitterness are more similar as respect to vegetable consumption (r = 0.63, p-value = .05). Moreover, a significant negative correlation was found between the percentage of Non Taster (NT) in different countries and the percentage of vegetable consumption (r = -0.87, p-value = .02), while a significant positive correlation emerged between the percentage of Super Taster (ST) and the percentage of vegetable liking (r = 0.87, p-value = .02). In our work we showed that differences in bitter perception among populations contributes to differences in vegetable liking and vegetable consumption. More in detail, populations with higher percentage of ST consume more vegetables than population where the majority of individuals are NT.

Genetic Landscape of Slovenians: Past Admixture and Natural Selection Pattern.

The Slovenian territory played a crucial role in the past serving as gateway for several human migrations. Previous studies used Slovenians as a source population to interpret different demographic events happened in Europe but not much is known about the genetic background and the demographic history of this population. Here, we analyzed genome-wide data from 96 individuals to shed light on the genetic role and history of the Slovenian population. Y chromosome diversity splits into two major haplogroups R1b and R1a with the latter suggesting a genetic contribution from the steppe. Slovenian individuals are more closely related to Northern and Eastern European populations than Southern European populations even though they are geographically closer. This pattern is confirmed by an admixture and clustering analysis. We also identified a single stream of admixture events between the Slovenians with Sardinians and Russians around ∼2630 BCE (2149-3112). Using ancient samples, we found a significant admixture in Slovenians using Yamnaya and the early Neolithic Hungarians as sources, dated around ∼1762 BCE (1099-2426) suggesting a strong contribution from the steppe to the foundation of the observed modern genetic diversity. Finally, we looked for signals of selection in candidate variants and we found significant hits in HERC2 and FADS responsible for blue eye color and synthesis of long-chain unsaturated fatty acids, respectively, when Slovenians were compared to Southern Europeans. While the comparison was done with Eastern Europeans, we identified significant signals in PKD2L1 and IL6R which are genes associated with taste and coronary artery disease, respectively.

Whole-genome sequencing reveals new insights into age-related hearing loss: cumulative effects, pleiotropy and the role of selection.

Age-related hearing loss (ARHL) is the most common sensory disorder in the elderly. Although not directly life threatening, it contributes to loss of autonomy and is associated with anxiety, depression and cognitive decline. To search for genetic risk factors underlying ARHL, a large whole-genome sequencing (WGS) approach has been carried out in a cohort of 212 cases and controls, both older than 50 years to select genes characterized by a burden of variants specific to cases or controls. Accordingly, the total variation load per gene was compared and two groups were detected: 375 genes more variable in cases and 371 more variable in controls. In both cases, Gene Ontology analysis showed that the largest enrichment for biological processes (fold > 5, p-value = 0.042) was the "sensory perception of sound", suggesting cumulative genetic effects were involved. Replication confirmed 141 genes, while additional analysis based on natural selection led to a prioritization of 21 genes. The majority of them (20 out of 21) showed positive expression in mouse cochlea cDNA and were associated with two functional pathways. Among them, two genes were previously associated with hearing (CSMD1 and PTRPD) and re-sequenced in a large Italian cohort of ARHL patients (N = 389). Results led to the identification of six coding variants not detected in cases so far, suggesting a possible protective role, which requires investigation. In conclusion, we show that this multistep strategy (WGS, selection, expression, pathway analysis and targeted re-sequencing) can provide major insights into the molecular characterization of complex diseases such as ARHL.