Assessing DNA Degradation through Differential Amplification Efficiency of Total Human and Human Male DNA in a Forensic qPCR Assay.

The assessment of degradation is crucial for the analysis of human DNA samples isolated from forensic specimens. Forensic quantitative PCR (qPCR) assays can include multiple targets of varying amplicon size that display differential amplification efficiency, and thus different concentrations, in the presence of degradation. The possibility of deriving information on DNA degradation was evaluated in a forensic qPCR assay not specifically designed to detect DNA fragmentation, the Plexor HY (Promega), by calculating the ratio between the estimated concentrations of autosomal (99 bp) and Y-chromosomal (133 bp) targets ("[Auto]/[Y]"). The [Auto]/[Y] ratio measured in 57 formalin-fixed, paraffin-embedded samples was compared to a quality score (QS) calculated for corresponding STR profiles using quantitative data (allele peak height). A statistically significant inverse correlation was observed between [Auto]/[Y] and QS (R = -0.65, p < 0.001). The [Auto]/[Y] values were highly correlated (R = 0.75, p < 0.001) with the "[Auto]/[D]" values obtained using the PowerQuant (Promega) assay, expressly designed to detect DNA degradation through simultaneous quantification of a short (Auto) and a long (D) autosomal target. These results indicate that it is possible to estimate DNA degradation in male samples through Plexor HY data and suggest an alternative strategy for laboratories lacking the equipment required for the assessment of DNA integrity through dedicated qPCR assays.

Functional characterization of archaic-specific variants in mitonuclear genes: insights from comparative analysis in S. cerevisiae.

Neanderthal and Denisovan hybridisation with modern humans has generated a non-random genomic distribution of introgressed regions, the result of drift and selection dynamics. Cross-species genomic incompatibility and more efficient removal of slightly deleterious archaic variants have been proposed as selection-based processes involved in the post-hybridisation purge of archaic introgressed regions. Both scenarios require the presence of functionally different alleles across Homo species onto which selection operated differently according to which populations hosted them, but only a few of these variants have been pinpointed so far. In order to identify functionally divergent archaic variants removed in humans, we focused on mitonuclear genes, which are underrepresented in the genomic landscape of archaic humans. We searched for non-synonymous, fixed, archaic-derived variants present in mitonuclear genes, rare or absent in human populations. We then compared the functional impact of archaic and human variants in the model organism Saccharomyces cerevisiae. Notably, a variant within the mitochondrial tyrosyl-tRNA synthetase 2 (YARS2) gene exhibited a significant decrease in respiratory activity and a substantial reduction of Cox2 levels, a proxy for mitochondrial protein biosynthesis, coupled with the accumulation of the YARS2 protein precursor and a lower amount of mature enzyme. Our work suggests that this variant is associated with mitochondrial functionality impairment, thus contributing to the purging of archaic introgression in YARS2. While different molecular mechanisms may have impacted other mitonuclear genes, our approach can be extended to the functional screening of mitonuclear genetic variants present across species and populations.

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.

A fresh perspective on infrared spectroscopy as a prescreening method for molecular and stable isotopes analyses on ancient human bones.

Following the development of modern genome sequencing technologies, the investigation of museum osteological finds is increasingly informative and popular. Viable protocols to help preserve these collections from exceedingly invasive analyses, would allow greater access to the specimens for scientific research. The main aim of this work is to survey skeletal tissues, specifically petrous bones and roots of teeth, using infrared spectroscopy as a prescreening method to assess the bone quality for molecular analyses. This approach could overcome the major problem of identifying useful genetic material in archaeological bone collections without resorting to demanding, time consuming and expensive laboratory studies. A minimally invasive sampling of archaeological bones was developed and bone structural and compositional changes were examined, linking isotopic and genetic data to infrared spectra. The predictive model based on Infrared parameters is effective in determining the occurrence of ancient DNA (aDNA); however, the quality/quantity of aDNA cannot be determined because of the influence of environmental and local factors experienced by the examined bones during the burial period.

Privacy and ethical challenges of the Amelogenin sex test in forensic paternity/kinship analysis: Insights from a 13-year case history.

The Amelogenin sex test included in forensic DNA typing kits has the potential to identify congenital conditions such as differences/disorders of sex development (DSD). It can also reveal mismatches between genotypic sex and gender marker in identity documents of transgender persons who obtained legal gender recognition. In a 13-year case history of paternity/kinship tests, involving n = 962 females and n = 1001 males, two mismatches between Amelogenin sex test (male) and gender marker (female), and three cases of chromosomal DSD (Klinefelter syndrome) were observed. The concrete risk of observing Amelogenin anomalies, their potential causes, and the context in which they occur (forensic, i.e. non-medical) mean that laboratory operators are called to strike a complex balance between privacy interests and individual health rights when providing preliminary information and reporting Amelogenin incidental findings. This case history argues for the need of a more responsible approach towards the Amelogenin sex test in the forensic community.

Recombulator-X: A fast and user-friendly tool for estimating X chromosome recombination rates in forensic genetics.

Genetic markers (especially short tandem repeats or STRs) located on the X chromosome are a valuable resource to solve complex kinship cases in forensic genetics in addition or alternatively to autosomal STRs. Groups of tightly linked markers are combined into haplotypes, thus increasing the discriminating power of tests. However, this approach requires precise knowledge of the recombination rates between adjacent markers. The International Society of Forensic Genetics recommends that recombination rate estimation on the X chromosome is performed from pedigree genetic data while taking into account the confounding effect of mutations. However, implementations that satisfy these requirements have several drawbacks: they were never publicly released, they are very slow and/or need cluster-level hardware and strong computational expertise to use. In order to address these key concerns we developed Recombulator-X, a new open-source Python tool. The most challenging issue, namely the running time, was addressed with dynamic programming techniques to greatly reduce the computational complexity of the algorithm. Compared to the previous methods, Recombulator-X reduces the estimation times from weeks or months to less than one hour for typical datasets. Moreover, the estimation process, including preprocessing, has been streamlined and packaged into a simple command-line tool that can be run on a normal PC. Where previous approaches were limited to small panels of STR markers (up to 15), our tool can handle greater numbers (up to 100) of mixed STR and non-STR markers. In conclusion, Recombulator-X makes the estimation process much simpler, faster and accessible to researchers without a computational background, hopefully spurring increased adoption of best practices.

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.

Liver Cancer-Specific Isoform of Serine Protease Inhibitor Kazal for the Detection of Hepatocellular Carcinoma: Results from a Pilot Study in Patients with Dysmetabolic Liver Disease.

Reliable non-invasive biomarkers for the surveillance of patients at risk of hepatocellular carcinoma (HCC) development represent an unmet medical need. Recently, the liver-cancer-specific isoform of serine protease inhibitor Kazal (LC-SPIK) has been proposed as a valuable biomarker for the detection of HCC in patients with chronic liver disease of viral etiology. In the present study, we assessed the diagnostic accuracy of LC-SPIK, alone or in combination with standard serologic biomarkers (i.e., alpha-fetoprotein and protein induced by vitamin K absence or antagonist-II, PIVKA-II), for the detection of HCC among patients with dysmetabolic liver disease. A total of 120 patients with non-alcoholic fatty liver disease (NAFLD), including 62 patients with a diagnosis of HCC and 58 with cirrhosis but without tumor, were retrospectively analyzed. The serum levels of LC-SPIK were measured by enzyme-linked immunosorbent assay (ImCare Biotech, Doylestown, PA). The serum LC-SPIK values were significantly different between patients with HCC (24.3, 17.6−39.8 ng/mL) and those with cirrhosis but without tumor (11.7, 8.7−18.2 ng/mL) (p < 0.001). By receiver operating characteristic curve analysis, we observed an area under the curve (AUC) of 0.841 for the detection of HCC; the combination with PIVKA-II further increased the accuracy to AUC = 0.926 (cross-validation). The promising results observed in the present pilot study foster additional research to investigate the usefulness of LC-SPIK for the stratification of the risk of HCC development in patients with NAFLD and advanced liver disease.

Evaluation of the Effects of Different Sample Collection Strategies on DNA/RNA Co-Analysis of Forensic Stains.

The aim of this study was to evaluate the impact of different moistening agents (RNase-free water, absolute anhydrous ethanol, RNAlater®) applied to collection swabs on DNA/RNA retrieval and integrity for capillary electrophoresis applications (STR typing, cell type identification by mRNA profiling). Analyses were conducted on whole blood, luminol-treated diluted blood, saliva, semen, and mock skin stains. The effects of swab storage temperature and the time interval between sample collection and DNA/RNA extraction were also investigated. Water provided significantly higher DNA yields than ethanol in whole blood and semen samples, while ethanol and RNAlater® significantly outperformed water in skin samples, with full STR profiles obtained from over 98% of the skin samples collected with either ethanol or RNAlater®, compared to 71% of those collected with water. A significant difference in mRNA profiling success rates was observed in whole blood samples between swabs treated with either ethanol or RNAlater® (100%) and water (37.5%). Longer swab storage times before processing significantly affected mRNA profiling in saliva stains, with the success rate decreasing from 91.7% after 1 day of storage to 25% after 7 days. These results may contribute to the future development of optimal procedures for the collection of different types of biological traces.

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

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

Genetics and Material Culture Support Repeated Expansions into Paleolithic Eurasia from a Population Hub Out of Africa.

The population dynamics that followed the Out of Africa (OoA) expansion and the whereabouts of the early migrants before the differentiation that ultimately led to the formation of Oceanian, West and East Eurasian macropopulations have long been debated. Shedding light on these events may, in turn, provide clues to better understand the cultural evolution in Eurasia between 50 and 35 ka. Here, we analyze Eurasian Paleolithic DNA evidence to provide a comprehensive population model and validate it in light of available material culture. Leveraging on our integrated approach we propose the existence of a Eurasian population Hub, where Homo sapiens lived between the OoA and the broader colonization of Eurasia, which was characterized by multiple events of expansion and local extinction. A major population wave out of Hub, of which Ust'Ishim, Bacho Kiro, and Tianyuan are unadmixed representatives, is broadly associated with Initial Upper Paleolithic lithics and populated West and East Eurasia before or around 45 ka, before getting largely extinct in Europe. In this light, we suggest a parsimonious placement of Oase1 as an individual related to Bacho Kiro who experienced additional Neanderthal introgression. Another expansion, started before 38 ka, is broadly associated with Upper Paleolithic industries and repopulated Europe with sporadic admixtures with the previous wave (GoyetQ116-1) and more systematic ones, whereas moving through Siberia (Yana, Mal'ta). Before these events, we also confirm Zlatý Kun as the most basal human lineage sequenced to date OoA, potentially representing an earlier wave of expansion out of the Hub.

The Genetic Origin of Daunians and the Pan-Mediterranean Southern Italian Iron Age Context.

The geographical location and shape of Apulia, a narrow land stretching out in the sea at the South of Italy, made this region a Mediterranean crossroads connecting Western Europe and the Balkans. Such movements culminated at the beginning of the Iron Age with the Iapygian civilization which consisted of three cultures: Peucetians, Messapians, and Daunians. Among them, the Daunians left a peculiar cultural heritage, with one-of-a-kind stelae and pottery, but, despite the extensive archaeological literature, their origin has been lost to time. In order to shed light on this and to provide a genetic picture of Iron Age Southern Italy, we collected and sequenced human remains from three archaeological sites geographically located in Northern Apulia (the area historically inhabited by Daunians) and radiocarbon dated between 1157 and 275 calBCE. We find that Iron Age Apulian samples are still distant from the genetic variability of modern-day Apulians, they show a degree of genetic heterogeneity comparable with the cosmopolitan Republican and Imperial Roman civilization, even though a few kilometers and centuries separate them, and they are well inserted into the Iron Age Pan-Mediterranean genetic landscape. Our study provides for the first time a window on the genetic make-up of pre-Roman Apulia, whose increasing connectivity within the Mediterranean landscape, would have contributed to laying the foundation for modern genetic variability. In this light, the genetic profile of Daunians may be compatible with an at least partial autochthonous origin, with plausible contributions from the Balkan peninsula.

Through 40,000 years of human presence in Southern Europe: the Italian case study.

The Italian Peninsula, a natural pier across the Mediterranean Sea, witnessed intricate population events since the very beginning of the human occupation in Europe. In the last few years, an increasing number of modern and ancient genomes from the area have been published by the international research community. This genomic perspective started unveiling the relevance of Italy to understand the post-Last Glacial Maximum (LGM) re-peopling of Europe, the earlier phase of the Neolithic westward migrations, and its linking role between Eastern and Western Mediterranean areas after the Iron Age. However, many open questions are still waiting for more data to be addressed in full. With this review, we summarize the current knowledge emerging from the available ancient Italian individuals and, by re-analysing them all at once, we try to shed light on the avenues future research in the area should cover. In particular, open questions concern (1) the fate of pre-Villabruna Europeans and to what extent their genomic components were absorbed by the post-LGM hunter-gatherers; (2) the role of Sicily and Sardinia before LGM; (3) to what degree the documented genetic structure within the Early Neolithic settlers can be described as two separate migrations; (4) what are the population events behind the marked presence of an Iranian Neolithic-like component in Bronze Age and Iron Age Italian and Southern European samples.

Usefulness of a Hepatitis B Surface Antigen-Based Model for the Prediction of Functional Cure in Patients with Chronic Hepatitis B Virus Infection Treated with Nucleos(t)ide Analogues: A Real-World Study.

In patients with chronic hepatitis B (CHB) under long-term treatment with nucleso(t)ide analogues (NAs), the loss of hepatitis B surface antigen (HBsAg) is a rare event. A growing body of evidence supports the use of quantitative HBsAg for the prediction of functional cure, although these results are mainly derived from studies performed on Asian patients with hepatitis B e antigen (HBeAg)-positive CHB. Here, we investigated the clinical role of quantitative HBsAg in a real-life cohort of CHB patients under treatment with NAs in a tertiary care center from North-West Italy. A total of 101 CHB patients (HBeAg-negative, n = 86) undergoing NAs treatment were retrospectively enrolled. HBsAg was measured at baseline (T0), 6 months (T1), 12 months (T2) and at the last follow-up (FU). Median FU was 5.5 (3.2-8.3) years; at the end of FU, 11 patients lost the HBsAg (annual incidence rate = 1.8%). Baseline HBsAg levels were significantly different between patients with no HBsAg loss and those achieving a functional cure (3.46, 2.91-3.97 vs. 1.11, 0.45-1.98 Log IU/mL, p < 0.001). Similarly, the HBsAg decline (Δ) from T0 to T2 was significantly different between the two groups of patients (0.05, -0.04-0.13, vs. 0.38, 0.11-0.80 Log IU/mL, p = 0.002). By stratified cross-validation analysis, the combination of baseline HBsAg and ΔHBsAg T0-T2 showed an excellent accuracy for the prediction of HBsAg loss (C statistic = 0.966). These results corroborate the usefulness of quantitative HBsAg in Caucasian CHB patients treated with antivirals for the prediction of HBsAg seroclearance.

Biomarkers of Oncogenesis, Adipose Tissue Dysfunction and Systemic Inflammation for the Detection of Hepatocellular Carcinoma in Patients with Nonalcoholic Fatty Liver Disease.

Current surveillance strategy for patients with nonalcoholic fatty liver disease (NAFLD) at risk of hepatocellular carcinoma (HCC) development is unsatisfactory. We aimed to investigate the diagnostic accuracy of alpha-fetoprotein (AFP), protein induced by vitamin K absence or antagonist-II (PIVKA-II), glypican-3 (GPC-3), adiponectin, leptin and interleukin-6 (IL-6), alone or in combination, for the discrimination between NAFLD patients with or without HCC. The biomarkers were investigated in a cohort of 191 NAFLD patients (median age 62, 54-68 years; 121 males and 70 females) with advanced fibrosis/cirrhosis, 72 of whom had a diagnosis of HCC. PIVKA-II showed the best performance for the detection of HCC with an area under the curve (AUC) of 0.853, followed by adiponectin (AUC = 0.770), AFP (AUC = 0.763), GPC-3 (AUC = 0.759) and by IL-6 (AUC = 0.731), while the leptin values were not different between patients with and without HCC. The accuracy of the biomarkers' combination was assessed by a stratified cross-validation approach. The combination of age, gender, PIVKA-II, GPC-3 and adiponectin further improved the diagnostic accuracy (AUC = 0.948); the model correctly identified the 87% of the patients. In conclusion, we developed a model with excellent accuracy for the detection of HCC that may be useful to improve the surveillance of NAFLD patients at risk of tumor development.

Analysis of Italian BRCA1/2 Pathogenic Variants Identifies a Private Spectrum in the Population from the Bergamo Province in Northern Italy.

Germline pathogenic variants (PVs) in the BRCA1 or BRCA2 genes cause high breast cancer risk. Recurrent or founder PVs have been described worldwide including some in the Bergamo province in Northern Italy. The aim of this study was to compare the BRCA1/2 PV spectra of the Bergamo and of the general Italian populations. We retrospectively identified at five Italian centers 1019 BRCA1/2 PVs carrier individuals affected with breast cancer and representative of the heterogeneous national population. Each individual was assigned to the Bergamo or non-Bergamo cohort based on self-reported birthplace. Our data indicate that the Bergamo BRCA1/2 PV spectrum shows less heterogeneity with fewer different variants and an average higher frequency compared to that of the rest of Italy. Consistently, four PVs explained about 60% of all carriers. The majority of the Bergamo PVs originated locally with only two PVs clearly imported. The Bergamo BRCA1/2 PV spectrum appears to be private. Hence, the Bergamo population would be ideal to study the disease risk associated with local PVs in breast cancer and other disease-causing genes. Finally, our data suggest that the Bergamo population is a genetic isolate and further analyses are warranted to prove this notion.

Functional and clinical implications of genetic structure in 1686 Italian exomes.

To reconstruct the phenotypical and clinical implications of the Italian genetic structure, we thoroughly analyzed a whole-exome sequencing data set comprised of 1686 healthy Italian individuals. We found six previously unreported variants with remarkable frequency differences between Northern and Southern Italy in the HERC2, OR52R1, ADH1B, and THBS4 genes. We reported 36 clinically relevant variants (submitted as pathogenic, risk factors, or drug response in ClinVar) with significant frequency differences between Italy and Europe. We then explored putatively pathogenic variants in the Italian exome. On average, our Italian individuals carried 16.6 protein-truncating variants (PTVs), with 2.5% of the population having a PTV in one of the 59 American College of Medical Genetics (ACMG) actionable genes. Lastly, we looked for PTVs that are likely to cause Mendelian diseases. We found four heterozygous PTVs in haploinsufficient genes (KAT6A, PTCH1, and STXBP1) and three homozygous PTVs in genes causing recessive diseases (DPYD, FLG, and PYGM). Comparing frequencies from our data set to other public databases, like gnomAD, we showed the importance of population-specific databases for a more accurate assessment of variant pathogenicity. For this reason, we made aggregated frequencies from our data set publicly available as a tool for both clinicians and researchers (http://nigdb.cineca.it; NIG-ExIT).

Small Non-Coding RNA Profiling in Plasma Extracellular Vesicles of Bladder Cancer Patients by Next-Generation Sequencing: Expression Levels of miR-126-3p and piR-5936 Increase with Higher Histologic Grades.

Bladder cancer (BC) is the tenth most frequent cancer worldwide. Due to the need for recurrent cystoscopies and the lack of non-invasive biomarkers, BC is associated with a high management burden. In this respect, small non-coding RNAs (sncRNAs) have been investigated in urine as possible biomarkers for BC, but in plasma their potential has not yet been defined. The expression levels of sncRNAs contained in plasma extracellular vesicles (EVs) from 47 men with BC and 46 healthy controls were assessed by next-generation sequencing. The sncRNA profiles were compared with urinary profiles from the same subjects. miR-4508 resulted downregulated in plasma EVs of muscle-invasive BC patients, compared to controls (adj-p = 0.04). In World Health Organization (WHO) grade 3 (G3) BC, miR-126-3p was upregulated both in plasma EVs and urine, when compared to controls (for both, adj-p < 0.05). Interestingly, two sncRNAs were associated with the risk class: miR-4508 with a downward trend going from controls to high risk BC, and piR-hsa-5936 with an upward trend (adj-p = 0.04 and adj-p = 0.05, respectively). Additionally, BC cases with low expression of miR-185-5p and miR-106a-5p or high expression of miR-10b-5p showed shorter survival (adj-p = 0.0013, adj-p = 0.039 and adj-p = 0.047, respectively). SncRNAs from plasma EVs could be diagnostic biomarkers for BC, especially in advanced grade.

Ancient Rome: A genetic crossroads of Europe and the Mediterranean.

Ancient Rome was the capital of an empire of ~70 million inhabitants, but little is known about the genetics of ancient Romans. Here we present 127 genomes from 29 archaeological sites in and around Rome, spanning the past 12,000 years. We observe two major prehistoric ancestry transitions: one with the introduction of farming and another prior to the Iron Age. By the founding of Rome, the genetic composition of the region approximated that of modern Mediterranean populations. During the Imperial period, Rome's population received net immigration from the Near East, followed by an increase in genetic contributions from Europe. These ancestry shifts mirrored the geopolitical affiliations of Rome and were accompanied by marked interindividual diversity, reflecting gene flow from across the Mediterranean, Europe, and North Africa.