Reconstructing the genetic history of late Neanderthals.

Although it has previously been shown that Neanderthals contributed DNA to modern humans, not much is known about the genetic diversity of Neanderthals or the relationship between late Neanderthal populations at the time at which their last interactions with early modern humans occurred and before they eventually disappeared. Our ability to retrieve DNA from a larger number of Neanderthal individuals has been limited by poor preservation of endogenous DNA and contamination of Neanderthal skeletal remains by large amounts of microbial and present-day human DNA. Here we use hypochlorite treatment of as little as 9 mg of bone or tooth powder to generate between 1- and 2.7-fold genomic coverage of five Neanderthals who lived around 39,000 to 47,000 years ago (that is, late Neanderthals), thereby doubling the number of Neanderthals for which genome sequences are available. Genetic similarity among late Neanderthals is well predicted by their geographical location, and comparison to the genome of an older Neanderthal from the Caucasus indicates that a population turnover is likely to have occurred, either in the Caucasus or throughout Europe, towards the end of Neanderthal history. We find that the bulk of Neanderthal gene flow into early modern humans originated from one or more source populations that diverged from the Neanderthals that were studied here at least 70,000 years ago, but after they split from a previously sequenced Neanderthal from Siberia around 150,000 years ago. Although four of the Neanderthals studied here post-date the putative arrival of early modern humans into Europe, we do not detect any recent gene flow from early modern humans in their ancestry.

A complete Neandertal mitochondrial genome sequence determined by high-throughput sequencing.

A complete mitochondrial (mt) genome sequence was reconstructed from a 38,000 year-old Neandertal individual with 8341 mtDNA sequences identified among 4.8 Gb of DNA generated from approximately 0.3 g of bone. Analysis of the assembled sequence unequivocally establishes that the Neandertal mtDNA falls outside the variation of extant human mtDNAs, and allows an estimate of the divergence date between the two mtDNA lineages of 660,000 +/- 140,000 years. Of the 13 proteins encoded in the mtDNA, subunit 2 of cytochrome c oxidase of the mitochondrial electron transport chain has experienced the largest number of amino acid substitutions in human ancestors since the separation from Neandertals. There is evidence that purifying selection in the Neandertal mtDNA was reduced compared with other primate lineages, suggesting that the effective population size of Neandertals was small.

Ancient gene flow from early modern humans into Eastern Neanderthals.

It has been shown that Neanderthals contributed genetically to modern humans outside Africa 47,000-65,000 years ago. Here we analyse the genomes of a Neanderthal and a Denisovan from the Altai Mountains in Siberia together with the sequences of chromosome 21 of two Neanderthals from Spain and Croatia. We find that a population that diverged early from other modern humans in Africa contributed genetically to the ancestors of Neanderthals from the Altai Mountains roughly 100,000 years ago. By contrast, we do not detect such a genetic contribution in the Denisovan or the two European Neanderthals. We conclude that in addition to later interbreeding events, the ancestors of Neanderthals from the Altai Mountains and early modern humans met and interbred, possibly in the Near East, many thousands of years earlier than previously thought.

Genetic sub-structuring of Croatian island populations in the Southeastern European context: a meta-analysis.

AIM: To use the method of meta-analysis to assess the influence of island population isolation on the sub-structuring of the Croatian population, as well as the influence of regional population groups on the sub-structuring of the Southeastern European population with regard to basic population genetic statistical parameters calculated by using STR locus analysis. METHODS: Bio-statistical analyses were performed for 2877 unrelated participants of both sexes from Southeastern Europe. Nine autosomal STR loci (D3S1358, vWA, FGA, TH01, TPOX, CSF1PO, D5S818, D13S317, and D7S82) were analyzed by using standard F-statistics and population structure analysis (Structure software). RESULTS: Genetic differentiation of Croatian subpopulations assessed with the FST method was higher at the level of the Croatian population (0.005) than at the level of Southeastern Europe (0.002). The island of Vis showed the most pronounced separation in the Croatian population, and Albanians from Kosovo in the population of Southeast Europe, followed by Croatia, Bosnia and Herzegovina, and Hungary. CONCLUSION: The higher structure of Croatian subpopulations in relation to Southeastern Europe suggest a certain degree of genetic isolation, most likely due to the influence of endogamy within rural island populations.

Influence of genetic substructuring of statistical forensic parameters on genetic short tandem repeat markers in the populations of Southeastern Europe.

AIM: To investigate the influence of specific intrapopulation genetic structures on interpopulation relationships. Special focus was the influence of island population isolation on the substructuring of the Croatian population, and the influence of regional population groups on the substructuring of Southeast European populations. METHODS: Autosomal short tandem repeat (STR) loci were analyzed by using four forensic parameters: matching probability (PM), power of discrimination (PD), power of exclusion (PE), and polymorphic information content (PIC) on a sample of 2877 unrelated participants of both sexes. A sample set comprising 590 participants was analyzed for the first time, and 2287 participants were included from previous studies. The analysis was performed with PowerStats v. 1.2. RESULTS: The analysis of forensic parameters for all nine loci in the Croatian subpopulations showed the largest deviations in the populations of the islands of Korcula and Hvar. The smallest deviations were found in the mainland population. As for Southeast European populations, the largest deviations were found in the population of North Macedonia, followed by Romania, Albanians from Kosovo, and Montenegro, while the smallest deviations were found in the population of Hungary. CONCLUSION: The comparison of forensic parameters between different subpopulations of Croatia and Southeast Europe indicates that the isolation of individual Croatian subpopulations and rare alleles in their gene pool affect the values of forensic parameters. Specific features of (sub)populations should be taken into account for appropriate sampling of the total population when creating a DNA database of STR markers.

Patterns of coding variation in the complete exomes of three Neandertals.

We present the DNA sequence of 17,367 protein-coding genes in two Neandertals from Spain and Croatia and analyze them together with the genome sequence recently determined from a Neandertal from southern Siberia. Comparisons with present-day humans from Africa, Europe, and Asia reveal that genetic diversity among Neandertals was remarkably low, and that they carried a higher proportion of amino acid-changing (nonsynonymous) alleles inferred to alter protein structure or function than present-day humans. Thus, Neandertals across Eurasia had a smaller long-term effective population than present-day humans. We also identify amino acid substitutions in Neandertals and present-day humans that may underlie phenotypic differences between the two groups. We find that genes involved in skeletal morphology have changed more in the lineage leading to Neandertals than in the ancestral lineage common to archaic and modern humans, whereas genes involved in behavior and pigmentation have changed more on the modern human lineage.

Genetic heritage of Croatians in the Southeastern European gene pool-Y chromosome analysis of the Croatian continental and Island population.

OBJECTIVES: The research objective of this study is to enlarge and deepen the Y chromosome research on the Croatian population and enable additional insights into the population diversity and historic events that shaped the current genetic landscape of Croatia and Southeastern Europe (SEE). MATERIALS AND METHODS: A high-resolution phylogenetic and phylogeographic analysis of 66 biallelic (SNPs) and 17 microsatellite (STRs) markers of the Y chromosome was performed using 720 Croatian samples. The obtained results were placed in a wider European context by comparison with ∼4450 samples from a number of other European populations. RESULTS: A high diversity of haplogroups was observed in the overall Croatian sample, and all typical European Y chromosome haplogroups with corresponding clinal patterns were observed. Three distinct genetic signals were identifiable in the Croatian paternal gene pool - I2a1b-M423, R1a1a1b1a*-M558, and E1b1b1a1b1a-V13 haplogroups. DISCUSSION: The analyses of the dominant and autochthonous I2a1b-M423 lineage (>30%) suggest that SEE had a significant role in the Upper Paleolithic, the R1a1a1b1a*-M558 lineage (19%) represents a signal from present day Slavic populations of Central Europe in the Croatian population, and the phylogeography of the E1b1b1a1b1a-V13 clade (around 9%) implies cultural diffusion of agriculture into Europe via the Balkan Peninsula. Am. J. Hum. Biol., 2016. © 2016 Wiley Periodicals, Inc. Am. J. Hum. Biol. 28:837-845, 2016. © 2016Wiley Periodicals, Inc.

Variability and distribution of COL1A2 (type I collagen) polymorphisms in the central-eastern Mediterranean Basin.

BACKGROUND: The most abundant of the collagen protein family, type I collagen is encoded by the COL1A2 gene. The COL1A2 restriction fragment length polymorphisms (RFLPs) EcoRI, RsaI and MspI in samples from several different central-eastern Mediterranean populations were analysed and found to be potentially informative anthropogenetic markers. AIM: The objective was to define the genetic variability of COL1A2 in the central-eastern Mediterranean and to shed light on its genetic distribution in human groups over a wide geographic area. SUBJECTS AND METHODS: PCR-RFLP analysis of EcoRI, RsaI and MspI polymorphisms of the COL1A2 gene was performed on oral swab and blood samples from 308 individuals from the central-eastern Mediterranean Basin. The genetic similarities among these groups and other populations described in the literature were investigated through correspondence analysis. RESULTS: Single-marker data and haplotype frequencies seemed to suggest a genetic homogeneity within the European populations, whereas a certain degree of differentiation was noted for the Egyptians and the Turks. CONCLUSIONS: The genetic variability in the central-eastern Mediterranean area is probably a result of the geographical barrier of the Mediterranean Sea, which separated European and African populations over time.

Genetic influences on plasma CFH and CFHR1 concentrations and their role in susceptibility to age-related macular degeneration.

It is a longstanding puzzle why non-coding variants in the complement factor H (CFH) gene are more strongly associated with age-related macular degeneration (AMD) than functional coding variants that directly influence the alternative complement pathway. The situation is complicated by tight genetic associations across the region, including the adjacent CFH-related genes CFHR3 and CFHR1, which may themselves influence the alternative complement pathway and are contained within a common deletion (CNP147) which is associated with protection against AMD. It is unclear whether this association is mediated through a protective effect of low plasma CFHR1 concentrations, high plasma CFH or both. We examined the triangular relationships of CFH/CFHR3/CFHR1 genotype, plasma CFH or CFHR1 concentrations and AMD susceptibility in combined case-control (1256 cases, 1020 controls) and cross-sectional population (n = 1004) studies and carried out genome-wide association studies of plasma CFH and CFHR1 concentrations. A non-coding CFH SNP (rs6677604) and the CNP147 deletion were strongly correlated both with each other and with plasma CFH and CFHR1 concentrations. The plasma CFH-raising rs6677604 allele and raised plasma CFH concentration were each associated with AMD protection. In contrast, the protective association of the CNP147 deletion with AMD was not mediated by low plasma CFHR1, since AMD-free controls showed increased plasma CFHR1 compared with cases, but it may be mediated by the association of CNP147 with raised plasma CFH concentration. The results are most consistent with a regulatory locus within a 32 kb region of the CFH gene, with a major effect on plasma CFH concentration and AMD susceptibility.

Fish and Shellfish Intake and Diabetes in a Costal Population of the Adriatic.

The objective of the study was to examine the association between fish and shellfish intake and diabetes in an island population, and the design of the study was Cross-sectional. Two independent population-based field surveys were conducted in Hvar Island of the eastern Adriatic coast of Croatia in May 2007 and May 2008, with a total of 1,379 adult participants. In multivariable logistic regression models, total fish intake was positively associated with diabetes prevalence in the total population (OR(Q4 vs. Q1) = 1.64; 95% CI = 1.01-2.66; p-trend = 0.09). Oily fish intake also exhibited a positive association with diabetes prevalence in the total population (OR(Q4 vs. Q1) = 2.22; 95% CI = 1.35-3.64; p-trend = 0.01) and in analyses stratified by body mass index, males and those with a high waist circumference. The study suggests an association between oily fish intake and diabetes in the population of the Hvar Island in Croatia. Longitudinal studies incorporating measures of persistent organic pollutants and local cooking practices are warranted to identify factors in fatty fish that may influence the development or persistence of diabetes.

Maternal genetic heritage of Southeastern Europe reveals a new Croatian isolate and a novel, local sub-branching in the x2 haplogroup.

High mtDNA variation in Southeastern Europe (SEE) is a reflection of the turbulent and complex demographic history of this area, influenced by gene flow from various parts of Eurasia and a long history of intermixing. Our results of 1035 samples (488 from Croatia, 239 from Bosnia and 130 from Herzegovina, reported earlier, and 97 Slovenians and 81 individuals from Zumberak, reported here for the first time) show that the SEE maternal genetic diversity fits within a broader European maternal genetic landscape. The study also shows that the population of Zumberak, located in the continental part of Croatia, developed some unique mtDNA haplotypes and elevated haplogroup frequencies due to distinctive demographic history and can be considered a moderate genetic isolate. We also report seven samples from the Bosnian population and one Herzegovinian sample designated as X2* individuals that could not be assigned to any of its sublineages (X2a'o) according to the existing X2 phylogeny. In an attempt to clarify the phylogeny of our X2 samples, their mitochondrial DNA has been completely sequenced. We suppose that these lineages are signs of local microdifferentiation processes that occurred in the recent demographic past in this area and could possibly be marked as SEE-specific X2 sublineages.

Modeling of environmental effects in genome-wide association studies identifies SLC2A2 and HP as novel loci influencing serum cholesterol levels.

Genome-wide association studies (GWAS) have identified 38 larger genetic regions affecting classical blood lipid levels without adjusting for important environmental influences. We modeled diet and physical activity in a GWAS in order to identify novel loci affecting total cholesterol, LDL cholesterol, HDL cholesterol, and triglyceride levels. The Swedish (SE) EUROSPAN cohort (N(SE) = 656) was screened for candidate genes and the non-Swedish (NS) EUROSPAN cohorts (N(NS) = 3,282) were used for replication. In total, 3 SNPs were associated in the Swedish sample and were replicated in the non-Swedish cohorts. While SNP rs1532624 was a replication of the previously published association between CETP and HDL cholesterol, the other two were novel findings. For the latter SNPs, the p-value for association was substantially improved by inclusion of environmental covariates: SNP rs5400 (p(SE,unadjusted) = 3.6 x 10(-5), p(SE,adjusted) = 2.2 x 10(-6), p(NS,unadjusted) = 0.047) in the SLC2A2 (Glucose transporter type 2) and rs2000999 (p(SE,unadjusted) = 1.1 x 10(-3), p(SE,adjusted) = 3.8 x 10(-4), p(NS,unadjusted) = 0.035) in the HP gene (Haptoglobin-related protein precursor). Both showed evidence of association with total cholesterol. These results demonstrate that inclusion of important environmental factors in the analysis model can reveal new genetic susceptibility loci.

Dietary patterns in adults from an Adriatic Island of Croatia and their associations with metabolic syndrome and its components.

Adriatic islanders have a high prevalence of metabolic syndrome (MetS) although they have traditionally practiced an active lifestyle and adhered to a Mediterranean diet. We performed a cross-sectional study to identify dietary patterns in a sample of 1442 adults from the island of Hvar, and determined whether MetS and its components: waist-circumference, serum triglycerides, fasting plasma glucose, HDL-cholesterol, and blood pressure, were related to an altered pattern of the traditional Mediterranean diet. Dietary intake was assessed by a food frequency questionnaire. MetS was defined using the International Diabetes Federation criteria. Our study showed that dietary patterns in this population have diversified from the traditional diet. Principal component analysis identified three major patterns. The meat, alcohol, and fish pattern (MAFp), sweets, grains, and fats pattern (SGFp), and an olive-oil, vegetables, and fruits pattern (OVFp) explained 30.6% of total dietary variance. The MAFp associated significantly with MetS (p = 0.027) and high plasma glucose (p = 0.006).

Refining the Global Phylogeny of Mitochondrial N1a, X, and HV2 Haplogroups Based on Rare Mitogenomes from Croatian Isolates.

Mitochondrial DNA (mtDNA) has been used for decades as a predominant tool in population genetics and as a valuable addition to forensic genetic research, owing to its unique maternal inheritance pattern that enables the tracing of individuals along the maternal lineage across numerous generations. The dynamic interplay between evolutionary forces, primarily genetic drift, bottlenecks, and the founder effect, can exert significant influence on genetic profiles. Consequently, the Adriatic islands have accumulated a subset of lineages that exhibits remarkable absence or rarity within other European populations. This distinctive genetic composition underscores the islands' potential as a significant resource in phylogenetic research, with implications reaching beyond regional boundaries to contribute to a global understanding. In the initial attempt to expand the mitochondrial forensic database of the Croatian population with haplotypes from small isolated communities, we sequenced mitogenomes of rare haplogroups from different Croatian island and mainland populations using next-generation sequencing (NGS). In the next step and based on the obtained results, we refined the global phylogeny of haplogroup N1a, HV2, and X by analyzing rare haplotypes, which are absent from the current phylogenetic tree. The trees were based on 16 novel and 52 previously published samples, revealing completely novel branches in the X and HV2 haplogroups and a new European cluster in the ancestral N1a variant, previously believed to be an exclusively African-Asian haplogroup. The research emphasizes the importance of investigating geographically isolated populations and their unique characteristics within a global context.

Genome-wide association of serum uric acid concentration: replication of sequence variants in an island population of the Adriatic coast of Croatia.

A genome-wide association study of serum uric acid (SUA) laevels was performed in a relatively isolated population of European descent from an island of the Adriatic coast of Croatia. The study sample included 532 unrelated and 768 related individuals from 235 pedigrees. Inflation due to relatedness was controlled by using genomic control. Genetic association was assessed with 2,241,249 single nucleotide polymorphisms (SNPs) in 1300 samples after adjusting for age and gender. Our study replicated four previously reported SUA loci (SLC2A9, ABCG2, RREB1, and SLC22A12). The strongest association was found with a SNP in SLC2A9 (rs13129697, P=2.33×10(-19)), which exhibited significant gender-specific effects, 35.76 µmol/L (P=2.11×10(-19)) in females and 19.58 µmol/L (P=5.40×10(-5)) in males. Within this region of high linkage disequilibrium, we also detected a strong association with a nonsynonymous SNP, rs16890979 (P=2.24×10(-17)), a putative causal variant for SUA variation. In addition, we identified several novel loci suggestive of association with uric acid levels (SEMA5A, TMEM18, SLC28A2, and ODZ2), although the P-values (P<5×10(-6)) did not reach the threshold of genome-wide significance. Together, these findings provide further confirmation of previously reported uric-acid-related genetic variants and highlight suggestive new loci for additional investigation.

Influence of evolutionary forces and demographic processes on the genetic structure of three Croatian populations: a maternal perspective.

BACKGROUND: Many Croatian islands are examples of genetic isolates, with low level of heterozygosity and high level of inbreeding, due to practice of endogamy. AIM: The aim was to study the genetic structure of two insular and one mainland population through high-resolution phylogenetic analysis of mitochondrial DNA (mtDNA). SUBJECTS AND METHODS: MtDNA polymorphisms were explored in 300 unrelated individuals from Mljet, Lastovo and the coastal city of Dubrovnik, based on SNP polymorphisms. RESULTS: All mtDNA haplogroups found in the sample were of typical European origin. However, the frequency distribution of their subclades differed significantly from other Croatian and European populations. MtDNA haplotype analysis revealed only two possible founder lineages on Mljet and six on Lastovo, accounting for almost half of the sample on both islands. The island of Mljet also has the lowest reported haplotype and nucleotide diversity among Croatian isolates and the island of Lastovo, a new sublineage of a usually quite rare U1b clade. CONCLUSION: The results can be explained by the effect evolutionary forces have on genetic structure, which is in line with the specific demographic histories of the islands. An additional research value of these two island isolates is the appearance of certain Mendelian disorders, highlighting their importance in epidemiological studies.

Maternal genetic legacy of the eastern Adriatic island of Krk--an interplay of evolutionary forces and island's historical events in shaping the genetic structure of contemporary island population.

This study presents genetic diversity and structure of contemporary Krk islanders revealed by high-resolution mitochondrial DNA analysis on a sample of 132 unrelated autochthonous adults from seven different settlements and regions of the island. Relatively high level of haplogroup and haplotype diversity in the overall island sample is an indicator of numerous migrations and gene flows throughout the history. Expectedly, the results show the highest frequency of haplogroup H (33.3%), yet this value is much lower compared to different Croatian and other European mainland populations. An interesting finding refers to highly elevated frequencies of some haplogroups, otherwise rare in Croatia and most of the Europe, such as I (11.3%) and W (7.6%) in Krk population, especially pronounced in some settlements. At the level of settlements, many of the major European haplogroups were found to be absent from their mtDNA gene pools, whereas several others show a pronounced deviation from an average. Overall, our results suggest a tangled interplay of different evolutionary forces, such as founder effects and a few strong bottlenecks, presumably due to epidemics, which have occurred in various periods of the island's history. Cultural customs, such as frequent endogamy in some regions of the island during past centuries, have additionally shaped its genetic structure into the observed present-day diversity patterns.

Genetic landscape of the Hvar Island - highlight of a 50-year long bioanthropological research.

The paper presents an overview of the 50-year long bioanthropological research of the Hvar islanders and depicts the maternal and paternal genetic landscape of the Hvar population (mtDNA and NRY lineages) in more detail. MtDNA haplogroups were determined in 169 and NRY haplogroups in 407 autochthonous individuals from the Hvar Island. The relatively high level of diversity of mtDNA and NRY lineages has been observed, however with interesting deviations from both the maternal (F1b1 lineage) and paternal (Q2a1a lineage) perspective. Additionally, population substructuring revealed differences between Hvar communities (east-west substructuring), in line with the ethnohistoric background and observed migration patterns on the island. Genetic analysis of the Hvar islanders presents a highlight of the 50-year long anthropological research on this island and offers insight into the current genetic structure of Dalmatia, Croatia, shaped by dynamic and diverse population movements throughout history.

Stable methylation loci are associated with systolic blood pressure in a Croatian island population.

Background: The objective was to identify stable and dynamic DNA methylation loci associated with cardiometabolic traits among an adult population from the Croatian island of Hvar. Materials & methods: An epigenome-wide association study was conducted using peripheral blood longitudinally collected at two time points 10 years apart via Infinium MethylationEPIC beadarray (n = 112). Stable and dynamic loci were identified using linear mixed models. Associations between cardiometabolic traits and loci were assessed using linear models. Results: 22 CpG loci were significantly associated with systolic blood pressure. Twenty were stable and two were dynamic. Conclusion: Multiple genes may be involved in the determination of systolic blood pressure level via stable epigenetic programming, potentially established earlier in life.

Cardiovascular disease is the leading cause of death worldwide. Previous studies have found that genetics incompletely explain susceptibility to cardiovascular disease. To find new potential risk factors, the authors investigated the possible contribution of DNA methylation (modifications to DNA that can affect gene expression but do not alter the underlying genetic code) in an adult population on the Croatian island of Hvar, which has a high number of people with cardiovascular and metabolic disease. By examining DNA methylation in blood collected at two time points, 10 years apart, the authors were able to identify DNA methylation that either stayed the same over time (stable) or changed the most over time (dynamic). These were then compared with clinical test results related to cardiovascular or metabolic diseases to determine if they are associated. Twenty-two methylation sites were found to be associated with systolic blood pressure. Of those, 20 were considered stable and two were dynamic. Additionally, there was one stable methylation site associated with serum calcium and one with C-reactive protein. These findings suggest that systolic blood pressure may be regulated through stable DNA methylation that is potentially established earlier in life.

Runs of homozygosity in European populations.

Estimating individual genome-wide autozygosity is important both in the identification of recessive disease variants via homozygosity mapping and in the investigation of the effects of genome-wide homozygosity on traits of biomedical importance. Approaches have tended to involve either single-point estimates or rather complex multipoint methods of inferring individual autozygosity, all on the basis of limited marker data. Now, with the availability of high-density genome scans, a multipoint, observational method of estimating individual autozygosity is possible. Using data from a 300,000 SNP panel in 2618 individuals from two isolated and two more-cosmopolitan populations of European origin, we explore the potential of estimating individual autozygosity from data on runs of homozygosity (ROHs). Termed F(roh), this is defined as the proportion of the autosomal genome in runs of homozygosity above a specified length. Mean F(roh) distinguishes clearly between subpopulations classified in terms of grandparental endogamy and population size. With the use of good pedigree data for one of the populations (Orkney), F(roh) was found to correlate strongly with the inbreeding coefficient estimated from pedigrees (r = 0.86). Using pedigrees to identify individuals with no shared maternal and paternal ancestors in five, and probably at least ten, generations, we show that ROHs measuring up to 4 Mb are common in demonstrably outbred individuals. Given the stochastic variation in ROH number, length, and location and the fact that ROHs are important whether ancient or recent in origin, approaches such as this will provide a more useful description of genomic autozygosity than has hitherto been possible.