Search

Background Analyses of the genomic variation in the western Mediterranean population are being used to reveal its evolutionary history and to understand the molecular basis of particular diseases. Aim To observe β-thalassemic mutational spectrum in western Andalusia, Spain, in the context of the Mediterranean. In addition, associations between disease and neutral gene variants within the β-globin gene (HBB) were also evaluated. Subjects and methods This study included 63 unrelated individuals diagnosed with β-thalassemia. In addition, 97 unrelated, healthy subjects of the same territory were also analysed as proxies of the normal genetic background. Allele associations and population genetic structure analyses were performed using different methodologies. Results Data have revealed a rather restricted spectrum of β-thalassemia mutations in the analysed sample. Although the detected variants fit well with the Mediterranean pattern, certain singularities support a structure of some specific β-thalassemia alleles. The IVSI-1 (G > A) shows a strong regionalisation. The spatial correlogram revealed a typically narrow wave structure, presumably linked to genetic isolation and genetic drift. Conclusions The long history of endemic malaria in the study territory, the rather high consanguinity rates among its autochthonous population, and other demographic features have been used here to understand the western Andalusian β-thalassemia molecular portrait.

About one-quarter of the euchromatic portion of the male-specific region of the human Y chromosome consists of large duplicated sequences that are organized in eight palindromes (termed P1–P8), which undergo arm-to arm gene conversion, a proposed mechanism for maintaining their sequence integrity. Although the relevance of gene conversion in the evolution of palindromic sequences has been profoundly recognized, the dynamic of this mechanism is still nuanced. To shed light into the evolution of these genomic elements, we performed a high-depth (50×) targeted next-generation sequencing of the palindrome P6 in 157 subjects belonging to the most divergent evolutionary lineages of the Y chromosome. We found 118 new paralogous sequence variants, which were placed into the context of a robust Y chromosome phylogeny based on 7240 SNPs of the X-degenerate region. We mapped along the phylogeny 80 gene conversion events that shaped the diversity of P6 arms during recent human history. In contrast to previous studies, we demonstrated that arm-to-arm gene conversion, which occurs at a rate of 6.01 × 10 −6 conversions/base/year, is not biased toward the retention of the ancestral state of sequences. We also found a significantly lower mutation rate of the arms (6.18 × 10−10 mutations/base/year) compared with the spacer (9.16 × 10−10 mutations/base/year), a finding that may explain the observed higher inter-species conservation of arms, without invoking any bias of conversion. Finally, by formally testing the mutation/conversion balance in P6, we found that the arms of this palindrome reached a steady-state equilibrium between mutation and gene conversion.

The human MSY ampliconic region is mainly composed of large duplicated sequences that are organized in eight palindromes (termed P1–P8), and may undergo arm-to-arm gene conversion. Although the importance of these elements is widely recognized, their evolutionary dynamics are still nuanced. Here, we focused on the P8 palindrome, which shows a complex evolutionary history, being involved in intra- and inter-chromosomal gene conversion. To disclose its evolutionary complexity, we performed a high-depth (50×) targeted next-generation sequencing of this element in 157 subjects belonging to the most divergent lineages of the Y chromosome tree. We found a total of 72 polymorphic paralogous sequence variants that have been exploited to identify 41 Y-Y gene conversion events that occurred during recent human history. Through our analysis, we were able to categorize P8 arms into three portions, whose molecular diversity was modelled by different evolutionary forces. Notably, the outer region of the palindrome is not involved in any gene conversion event and evolves exclusively through the action of mutational pressure. The inner region is affected by Y-Y gene conversion occurring at a rate of 1.52 × 10−5 conversions/base/year, with no bias towards the retention of the ancestral state of the sequence. In this portion, GC-biased gene conversion is counterbalanced by a mutational bias towards AT bases. Finally, the middle region of the arms, in addition to intra-chromosomal gene conversion, is involved in X-to-Y gene conversion (at a rate of 6.013 × 10−8 conversions/base/year) thus being a major force in the evolution of the VCY/VCX gene family.

Objectives The 3' regulatory region of the immunoglobulin heavy chain gene (IGH) includes the HS1.2 enhancer displaying length polymorphism with four known variants. The goal of the research was to provide an overview of this variability and of its evolutionary significance across human populations. Materials and methods We compiled published and original data on HS1.2 polymorphism in 3,100 subjects from 26 human populations. Moreover, we imputed the haplotypic arrangement of the HS1.2 region in the 1000 Genomes Project (1KGP). In this dataset, imputation could also be obtained for the G1m-G3m allotype by virtue of the precise correspondence between serological types and amino acid (and DNA) substitutions in IGHG1 and IGHG3. Results HS1.2 variant frequencies displayed similar patterns of continental partitioning as those reported in the literature for the physically neighboring IGHG1-IGHG3 system. The 1KGP data revealed that linkage disequilibrium (LD) can explain the spread of joint HS1.2-IGHG1-IGHG3 associations across continents and within continental populations, with stronger LD out of Africa and the features of an evolutionarily stable genomic block with differential expression in lymphoblastoid cell lines. Discussion Strong population structuring involves at least the entire 70 kb genomic region here considered, due to the tight LD which maintained HS1.2, IGHG1, and IGHG3 in nonrandom arrangements. This might be key to better understand the evolutionary path of the entire genomic region driven by immune response capabilities, during the formation of continental gene pools.

Kinship recognition between anonymous DNA samples is becoming a relevant issue in forensics, more so with the increasing number of DNA profiles in databanks. Also, NGS-based genotyping is being increasingly used in routine personal identification, to simultaneously type large numbers of markers of different kind. In the present work, we explored computationally and experimentally the performance of the ForenSeq™ DNA Signature Prep Kit in identifying the true relationship between two anonymous samples, distinguishing it from other possible relationships. We analyzed with Familias R series of 10,000 pairs with 9 different simulated relationships, corresponding to different degrees of autosomal sharing. For each pair we obtained likelihood ratios for five kinship hypotheses vs. unrelatedness, and used their ranking to identify the preferred relationship. We also typed 21 subjects from two pedigrees, representing from parent-child to 4th cousins relationships. As expected, the power for identifying the true relationship decays in the order of autosomal sharing. Parent-child and full siblings can be robustly identified against other relationships. For half-siblings the chance of reaching a significant conclusion is already small. For more distant relationships the proportion of cases correctly and significantly identified is 10% or less. Bidirectional errors in kinship attribution include the suggestion of relatedness when this does not exist (10–50%), and the suggestion of independence in pairs of individuals more than 4 generations apart (25–60%). The real cases revealed a relevant effect of genotype miscalling at some loci, which could only be partly avoided by modulating the analysis parameters. In conclusion, with the exception of first degree relatives, the kit can be useful to inform additional investigations, but does not usually provide probatory results.

Several thermal areas, also used for leisure purposes, may represent suitable habitats for free-living amoebae (FLAs), but few studies have been carried out in search for these organisms. The aim of this study was to assess the presence and distribution of FLAs by culture detection and molecular identification, over a one year-round sampling of two sites in Central Italy. Two geothermal springs (Site A and Site B) were investigated for a total of 36 water samples. Four sets of primers were used to amplify FLA DNA from all cultures positive for amoebic growth at both 37 °C and 45 °C. Overall, 33 (91.6%) water samples produced PCR amplification. Eleven taxa were identified. The array of identified species varied over the sampling period, and differed between the two hot springs, Site A harbouring 11 taxa compared to 5 of site B. However, both sites were characterized by the most common species Vermamoeba vermiformis and Naegleria australiensis. Acanthamoeba genotypes T4 and T15 were found at low frequency. Differences in the composition between the two sites could reflect environmental changes in biotic and chemical/physical parameters. From a public health perspective, the detection of potentially pathogenic amoebae could unveil a potential risk for humans.

Background: The compilation of a list of genetic modifiers in Alzheimer’s disease (AD) is an open research field. The GABAergic system is affected in several neurological disorders but its role in AD is largely understudied. Objective/Methods: As an explorative study, we considered variants in genes of GABA catabolism (ABAT, ALDH5A1, AKR7A2), and APOE in 300 Italian patients and 299 controls. We introduce a recent multivariate method to take into account the individual APOE genotype, thus controlling for the effect of the discrepant allele distributions in cases versus controls. We add a genotype-phenotype analysis based on age at onset and the Mini-Mental State Evaluation score. Results: On the background of strongly divergent APOE allele distributions in AD versus controls, two genotypic interactions that represented a subtle but significant peculiarity of the AD cohort emerged. The first is between ABAT and APOE, and the second between some ALDH5A1 genotypes and APOE. Decreased SSADH activity is predicted in AD carriers of APOE ɛ4, representing an additional suggestion for increased oxidative damage. Conclusion: We identified a difference between AD and controls, not in a shift of the allele frequencies at genes of the GABA catabolism pathway, but rather in gene interactions peculiar of the AD cohort. The emerging view is that of a multifactorial contribution to the disease, with a main risk factor (APOE), and additional contributions by the variants here considered. We consider genes of the GABA degradation pathway good candidates as modifiers of AD, contributing to energy impairment in AD brain.

Blastocystis spp. is a common single-celled intestinal symbiont, comprising several genetic subtypes (ST) and transmissible by animal-to-animal, human-to-human, animal-to-human and, possibly, human-to-animal routes. This work was designed to explore the presence of Blastocystis in sympatric domestic and wild suids and their ability to carry zoonotic STs, in a condition of widespread opportunity to come in contact with the microorganism through their shared water and food resources, and other carriers. We sampled 42 and 37 stool samples from wild boars and domestic pigs, respectively. STs were first identified by PCR followed by Sanger sequencing. Sequences represented in double-band PCR products or in Sanger chromatograms displaying multiple peaks, were resolved by next generation sequencing (NGS). Twenty-six (61.9%) wild boar and 26 (70.2%) pig samples were PCR-positive, respectively. ST3, ST5 and ST15 were found in 3.8%, 38.4% and 80.8% of the positive wild boars and 11.5%, 88.5%, 11.5% of the positive pigs, respectively. ST1 was found only in pigs (3.8%). STs 5 and 15 were common in both groups of animals, but in reversed proportions, suggesting preferential colonization. We found significantly different ST distributions among wild boars and domestic pigs. This might indicate that lifestyle differences between the two populations influence their risk for contracting certain subtypes, or that ST5 and ST15 can colonize preferentially wild or domestic animals. Based on the STs described here, wild boars and domestic pigs can act as reservoirs with zoonotic potential. The ability of suids to carry zoonotic STs appears to be higher when using NGS than Sanger sequencing, and resolution of complex sequencing profiles is imperative before excluding the presence of STs of human concern.

Background: Tetranucleotide Short Tandem Repeats (STRs) for human identification and common use in forensic cases have recently been used to address the population genetics of the North-Eastern Mediterranean area. However, to gain confidence in the inferences made using STRs, this kind of analysis should be challenged with changes in three main aspects of the data, i.e. the sizes of the samples, their distance across space and the genetic background from which they are drawn.Aim: To test the resilience of the gradients previously detected in the North-Eastern Mediterranean to the enlargement of the surveyed area and population set, using revised data.Subjects and methods: STR genotype profiles were obtained from a publicly available database (PopAffilietor databank) and a dataset was assembled including >7000 subjects from the Arabian Peninsula to Scandinavia, genotyped at eight loci. Spatial principal component analysis (sPCA) was applied and the frequency maps of the nine alleles which contributed most strongly to sPC1 were examined in detail.Results: By far the greatest part of diversity was summarised by a single spatial principal component (sPC1), oriented along a SouthEast-to-NorthWest axis. The alleles with the top 5% squared loadings were TH01(9.3), D19S433(14), TH01(6), D19S433(15.2), FGA(20), FGA(24), D3S1358(14), FGA(21) and D2S1338(19). These results confirm a clinal pattern over the whole range for at least four loci (TH01, D19S433, FGA, D3S1358).Conclusions: Four of the eight STR loci (or even alleles) considered here can reproducibly capture continental arrangements of diversity. This would, in principle, allow for the exploitation of forensic data to clarify important aspects in the formation of local gene pools., Italian Ministry of Justice [CUP E81J10001270005]

Naegleria spp. are free-living amoebae belonging to the family Vahlkampfiidae, in the class Heterolobosea. Among the recognized species, Naegleria fowleri causes primary amoebic meningoencephalitis (PAM), while two other species, Naegleria australiensis and Naegleria italica, have been reported as pathogenic in experimental animals. Due to the thermotolerance properties of some species, geothermal water sources including hot springs represent suitable habitats for their proliferation. The main aim of this study was a year-round sampling in two geothermal springs in Central Italy, to investigate the presence of Naegleria spp. using PCR/DNA sequencing based methods. The affinities between the sequences generated here and others reported in the literature were explored by using POY, which implements the concept of dynamic homology. Naegleria australiensis, Naegleria italica, and Naegleria lovaniensis, plus an unassigned Naegleria spp. were detected. Indels in the rDNA ITS1 and ITS2 turned out to be critical to distinguish the three species and confirmed their phylogenetic relationships. This is the first molecular report on the Naegleria spp. occurrence in geothermal waters in Central Italy, coupled with a fine genetic characterization.

Background Today, African pastoralists are found mainly in the Sahel/Savannah belt spanning 6,000 km from west to east, flanked by the Sahara to the north and tropical rainforests to the south. The most significant group among them are the Fulani who not only keep cattle breeds of possible West Eurasian ancestry, but form themselves a gene pool containing some paternally and maternally-transmitted West Eurasian haplogroups. Materials and Methods We generated complete sequences for 33 mitogenomes belonging to haplogroups H1 and U5 (23 and 10, respectively), and genotyped 16 STRs in 65 Y chromosomes belonging to haplogroup R1b-V88. Results We show that age estimates of the maternal lineage H1cb1, occurring almost exclusively in the Fulani, point to the time when the first cattle herders settled the Sahel/Savannah belt. Similar age estimates were obtained for paternal lineage R1b-V88, which occurs today in the Fulani but also in other, mostly pastoral populations. Maternal clade U5b1b1b, reported earlier in the Berbers, shows a shallower age, suggesting another possibly independent input into the Sahelian pastoralist gene pool. Conclusions Despite the fact that animal domestication originated in the Near East ∼ 10 ka, and that it was from there that animals such as sheep, goats as well as cattle were introduced into Northeast Africa soon thereafter, contemporary cattle keepers in the Sahel/Savannah belt show uniparental genetic affinities that suggest the possibility of an ancient contact with an additional ancestral population of western Mediterranean ancestry.

DNA typing and genetic profile data interpretation are among the most relevant topics in forensic science; among other applications, genetic profile's capability to distinguish biogeographic information about population groups, subgroups and affiliations have been largely explored in the last decade. In fact, for investigative and intelligence purposes, it is extremely useful to identify subjects and estimate their biogeographic origins by examining the recovered DNA profiles from evidence on a crime scene. Current approaches for BiogeoGraphic Ancestry (BGA) estimation using STRs profiles are usually based on Bayesian methods, which quantify the evidence in terms of likelihood ratio, supporting or not the hypothesis that a certain profile belongs to a specific ethnic group. The present study provides an alternative approach to the likelihood ratio method that involves multivariate data analysis strategies for the estimation of multiple populations. Starting from the well-known NIST US autosomal STRs dataset involving African-American, Asian, and Caucasian individuals, and moving towards further and more geographically restricted populations (such as Northern Africans vs sub-Saharan Africans, Afghans vs Iraqis and Italians vs Romanians), powerful multivariate techniques such as Sparse and Logistic Principal Component Analysis (SL-PCA), Sparse Partial Least Squares-Discriminant Analysis (sPLS-DA) and Support Vector Machines (SVM) were employed and their discriminating power was also compared. Both sPLS-DA and SVM techniques provided robust classifications, yielding high sensitivity and specificity models capable of discriminating populations on ethnic basis. This application may represent a powerful and dynamic tool for law enforcement agencies whenever a standard autosomal STR profile is obtained from the biological evidence collected at a crime scene or recovered during mass-disaster and missing person investigations.

Vermamoeba vermiformis represents one of the most common free-living amoebae identified in worldwide environmental surveys. We analyzed 56 water samples with varying characteristics, including temperature and the particular settings in which humans may be exposed to water, plus one corneal scraping from a keratitis patient, with the following aims: (i) to investigate the presence of V. vermiformis; (ii) to identify the isolate subtypes; (iii) to place the Italian isolates in the broader picture of the genetic diversity within V. vermiformis. Twenty-two isolates were identified upon culturing and sequencing of > 600 bp in the 18S ribosomal RNA (rRNA) gene sequence, bringing to 27 the number of sequences recovered from Italian sources. By adding deposited sequences, we assembled a dataset of 74 isolates. Three of our isolates were characterized by allelic code 7-5-1-1, never reported before, and two showed 100% identity with an uncultured eukaryote and carried the 719T>C variant. We show that the variable segments E5, E3, F, and G convey most of the information on diversity, enabling the clustering of the isolates in a replicable fashion. The presence of different strains in natural thermal waters and in distribution systems indicated heterogeneity of the amoebic populations. Also, ours and the only other sequence from human infection were mapped in different clades. Overall, we enlarged the repertoire of single nucleotide and indel variants and the list of allelic codes, proceeding one step further in the description of the diversity within the genus.

Throughout the past few years, a lively debate emerged about the timing and magnitude of the human migrations between the Iberian Peninsula and the Maghreb. Several pieces of evidence, including archaeological, anthropological, historical, and genetic data, have pointed to a complex and intermingled evolutionary history in the western Mediterranean area. To study to what extent connections across the Strait of Gibraltar and surrounding areas have shaped the present-day genomic diversity of its populations, we have performed a screening of 2.5 million single-nucleotide polymorphisms in 142 samples from southern Spain, southern Portugal, and Morocco. We built comprehensive data sets of the studied area and we implemented multistep bioinformatic approaches to assess population structure, demographic histories, and admixture dynamics. Both local and global ancestry inference showed an internal substructure in the Iberian Peninsula, mainly linked to a differential African ancestry. Western Iberia, from southern Portugal to Galicia, constituted an independent cluster within Iberia characterized by an enriched African genomic input. Migration time modeling showed recent historic dates for the admixture events occurring both in Iberia and in the North of Africa. However, an integrative vision of both paleogenomic and modern DNA data allowed us to detect chronological transitions and population turnovers that could be the result of transcontinental migrations dating back from Neolithic times. The present contribution aimed to fill the gaps in the modern human genomic record of a key geographic area, where the Mediterranean and the Atlantic come together.

The human cytochrome P4501A1 (CYP1A1) enzyme plays an important role in the metabolism of xenobiotics and endogenous substrates. Because polymorphisms within the CYP1A1 gene have been shown to be associated with various cancer risks and with the predicting clinical efficacy of some chemotherapies in different populations, most studies focus on their clinical significance. We, however, were interested in evaluating whether the polymorphisms could be used to distinguish human populations. Four single nucleotide CYP1A1 polymorphisms (rs4646903/ g.75011641; rs1048943/g.75012985; g.75012235; and rs1799814/ g.75012987) were analysed via PCR-RFLP assay in 1,195 individuals of various human groups from all over the world. In order to gain a more complete view of the genetic variability of the CYP1A1 gene, different statistical analyses were performed upon the populations of the present study and upon the limited data gleaned from previously studied populations. The allele and haplotype frequencies vary among populations: the rs4646903 (C) and rs1048943 (G) have been found to be nearly always linked and were found at the highest frequencies in Native Americans, while the variant associated to the position g.75012235 was only detected in certain African populations. Our work clearly indicates that the CYP1A1 polymorphisms differ among populations and that the prediction of genotypes constitutes an important aspect of precision medicine since some variants were associated with certain cancers and rs1048943 show strong association with optimized chemotherapy. Moreover, the CYP1A1 gene plays an important role in the metabolism of xenobiotics and it is likely that its frequencies could be strongly influenced by environmental factors.

The remote island of Montecristo is renowned for a resident goat population, whose origins are unclear. We describe the complete mtDNA sequence of a male specimen captured on the island. The sequence turned out to belong to haplogroup (Hg) A. Our results exclude that the sequence belongs to a lineage more closely related to bezoars than domestic goats. The lineages most similar to the Montecristo sequence are currently found in Western Europe, favouring partial feralization for the Montecristo isolate. Positioning of the Montecristo mitogenome in the emerging phylogenetic tree for domestic goats reveals a new lineage with multiple derived nucleotide states shared with lineages so far described in Europe as well as in North Asian breeds. This reveals the presence, within Hg A, of lineages with a palearctic distribution, whose descendants are now grown from Eastern Asia to Western Europe.

This work aimed to describe the genetic landscape of the Balkan Peninsula, as revealed by STR markers commonly used in forensics and spatial methods specifically developed for genetic data.We generated and analyzed 16 short tandem repeats (STRs) autosomal genotypes in 287 subjects from ten administrative/geographical regions of Eastern Europe (Romania and the Republic of Moldova). We report estimates of the allele frequencies in these sub-populations, their fixation indexes, and use these results to complement previous spatial analyses of Southern Europe.In seven out of ten analyzed regional samples the heterozygosity, averaged across loci, was lower than expected. The average Fis was 0.011. Among the 16 loci, five returned a significant fixation index Fst. The composite Fst across the 16 loci, among the 10 regional samples, was 0.00417, a figure twice as large as that obtained with the same markers across the entire Northern Mediterranean. The first spatial principal component (sPC1) returned the picture of a Central-European pattern of frequencies for the Carpathians, which extended to the Southern boundary of the Balkan Peninsula. However, the 8 alleles extracted by sPC1 returned a picture of a strong reduction of the migration rate in the Carpathian region, mostly between the inner locations.Our results revealed an unexpected heterogeneity in the area. We believe that populations from some regions will require treatment as distinct entities when considered in forensic applications.

In order to improve the phylogeography of the male-specific genetic traces of Greek and Phoenician colonizations on the Northern coasts of the Mediterranean, we performed a geographically structured sampling of seven subclades of haplogroup J in Turkey, Greece and Italy. We resequenced 4.4 Mb of Y-chromosome in 58 subjects, obtaining 1079 high quality variants. We did not find a preferential coalescence of Turkish samples to ancestral nodes, contradicting the simplistic idea of a dispersal and radiation of Hg J as a whole from the Middle East. Upon calibration with an ancient Hg J chromosome, we confirmed that signs of Holocenic Hg J radiations are subtle and date mainly to the Bronze Age. We pinpointed seven variants which could potentially unveil star clusters of sequences, indicative of local expansions. By directly genotyping these variants in Hg J carriers and complementing with published resequenced chromosomes (893 subjects), we provide strong temporal and distributional evidence for markers of the Greek settlement of Magna Graecia (J2a-L397) and Phoenician migrations (rs760148062). Our work generated a minimal but robust list of evolutionarily stable markers to elucidate the demographic dynamics and spatial domains of male-mediated movements across and around the Mediterranean, in the last 6,000 years.

Migration is one of several marine vertebrate traits increasingly affected by human encroachment. The Adriatic Sea is an important foraging and wintering site for loggerhead turtles (Caretta caretta), and one of the Mediterranean regions where they are most heavily impacted, particularly by fisheries bycatch. Conservation measures concern foraging and wintering areas and nesting grounds, and must, therefore, be informed by the natal origin of individuals in these mixed aggregates. Genetic diversity was investigated among 488 loggerheads found stranded or incidentally captured in fishing gear across the Adriatic and the origin of individuals was assessed through mixed-stock analysis based on comparison of 755 bp of mitochondrial DNA control region sequence. In addition, we highlight the importance of the Gulf of Manfredonia (Apulia, Italy), and of the northeastern Adriatic--a previously genetically undescribed foraging ground. No significant genetic divergence was recorded among sampling areas, between turtles sampled in cold and warm months or between adults and juveniles. The distribution of turtles across the Adriatic Sea appeared not to depend on individual origin. Rookeries in western Greece and Crete provided the most important contributions to the Adriatic mixed stock. In particular, the Greek populations were the most abundant locally in the Gulf of Manfredonia, so they are likely to experience an even greater impact than previously thought because of the severe fishery bycatch levels in this area. This study also provides an example of how substantial increases in sample sizes permit a relatively comprehensive testing of genetic structure across groupings in foraging aggregations.

Background The structure of haplogroup H reveals significant differences between the western and eastern edges of the Mediterranean, as well as between the northern and southern regions. Human populations along the westernmost Mediterranean coasts, which were settled by individuals from two continents separated by a relatively narrow body of water, show the highest frequencies of mitochondrial haplogroup H. These characteristics permit the analysis of ancient migrations between both shores, which may have occurred via primitive sea crafts and early seafaring. We collected a sample of 750 autochthonous people from the southern Iberian Peninsula (Andalusians from Huelva and Granada provinces). We performed a high-resolution analysis of haplogroup H by control region sequencing and coding SNP screening of the 337 individuals harboring this maternal marker. Our results were compared with those of a wide panel of populations, including individuals from Iberia, the Maghreb, and other regions around the Mediterranean, collected from the literature. Results Both Andalusian subpopulations showed a typical western European profile for the internal composition of clade H, but eastern Andalusians from Granada also revealed interesting traces from the eastern Mediterranean. The basal nodes of the most frequent H sub-haplogroups, H1 and H3, harbored many individuals of Iberian and Maghrebian origins. Derived haplotypes were found in both regions; haplotypes were shared far more frequently between Andalusia and Morocco than between Andalusia and the rest of the Maghreb. These and previous results indicate intense, ancient and sustained contact among populations on both sides of the Mediterranean. Conclusions Our genetic data on mtDNA diversity, combined with corresponding archaeological similarities, provide support for arguments favoring prehistoric bonds with a genetic legacy traceable in extant populations. Furthermore, the results presented here indicate that the Strait of Gibraltar and the adjacent Alboran Sea, which have often been assumed to be an insurmountable geographic barrier in prehistory, served as a frequently traveled route between continents. Electronic supplementary material The online version of this article (doi:10.1186/s12863-017-0514-6) contains supplementary material, which is available to authorized users.

Sequence diversity and the ages of the deepest nodes of the MSY phylogeny remain largely unexplored due to the severely biased collection of SNPs available for study. We characterized 68 worldwide Y chromosomes by high-coverage next-generation sequencing, including 18 deep-rooting ones, and identified 2386 SNPs, 80% of which were novel. Many aspects of this pool of variants resembled the pattern observed among genome-wide de novo events, suggesting that in the MSY, a large proportion of newly arisen alleles has survived in the phylogeny. Some degree of purifying selection emerged in the form of an excess of private missense variants. Our tree recapitulated the previously known topology, but the relative lengths of major branches were drastically modified and the associated node ages were remarkably older. We found significantly different branch lengths when comparing the rare deep-rooted A1b African lineage with the rest of the tree. Our dating results and phylogeography led to the following main conclusions: (1) Patrilineal lineages with ages approaching those of early AMH fossils survive today only in central-western Africa; (2) only a few evolutionarily successful MSY lineages survived between 160 and 115 kya; and (3) an early exit out of Africa (before 70 kya), which fits recent western Asian archaeological evidence, should be considered. Our experimental design produced an unbiased resource of new MSY markers informative for the initial formation of the anatomically modern human gene pool, i.e., a period of our evolution that had been previously considered to be poorly accessible with paternally inherited markers.

In migratory species female- and male-mediated gene flow are important for defining relevant Management Units, and for evaluating connectivity between these and their respective foraging grounds. The stock composition at five Mediterranean foraging areas was investigated by analysing variation in the mitochondrial D-loop and six microsatellite loci in a sample of 268 loggerhead turtles (Caretta caretta) stranded or accidentally caught by fisheries. This involved a comprehensive Mixed Stock Analysis which considers also recent data from major rookeries in Libya and Turkey, and the generation of a standardized nomenclature of allele sizes at the microsatellite loci. The results indicate: that the north Adriatic, the Tunisian continental shelf, the waters around Malta and the Italian Ionian Sea represent important areas for the conservation of rookeries in Greece, Libya and Turkey, respectively; that waters off the Italian peninsula and the islands of Lampedusa and Malta are mainly inhabited by individuals of Mediterranean origin, with a major contribution from the nearest and largest colonies, while Atlantic turtles are restricted to the western areas; that specific migratory routes exist from rookeries to foraging grounds; a poor bi-parental genetic structuring, which suggests a high male-mediated gene flow in the Mediterranean; mixing of small turtles in waters distant from natal rookeries, and recovery of structuring for large-sized individuals; and that uncommon mtDNA haplotypes are more powerful markers than microsatellite alleles in assessing an individual's origin, owing to their higher geographic specificity. Copyright © 2013 John Wiley & Sons, Ltd.

The African Sahel is conducive to studies of divergence/admixture genetic events as a result of its population history being so closely related with past climatic changes. Today, it is a place of the co-existence of two differing food-producing subsistence systems, i.e., that of sedentary farmers and nomadic pastoralists, whose populations have likely been formed from several dispersed indigenous hunter-gatherer groups. Using new methodology, we show here that the male gene pool of the extant populations of the African Sahel harbors signatures of multiple and differentiated contributions from different genetic sources. We also show that even if the Fulani pastoralists and their neighboring farmers share high frequencies of four Y chromosome subhaplogroups of E, they have drawn on molecularly differentiated subgroups at different times. These findings, based on combinations of SNP and STR polymorphisms, add to our previous knowledge and highlight the role of differences in the demographic history and displacements of the Sahelian populations as a major factor in the segregation of the Y chromosome lineages in Africa. Interestingly, within the Fulani pastoralist population as a whole, a differentiation of the groups from Niger is characterized by their high presence of R1b-M343 and E1b1b1-M35. Moreover, the R1b-M343 is represented in our dataset exclusively in the Fulani group and our analyses infer a north-to-south African migration route during a recent past.

We sequenced coding portions (1.6 kb) of the mtDNA in 170 loggerhead (Caretta caretta) turtles sampled in the central Mediterranean. The sequences spanned the entire ND1 and ND3 genes, the tRNAGly and tRNAArg, plus the 3' and 5' termini of COXIII and ND4L genes, respectively. Based on our sequencing results and published complete mitogenomes, we constructed a maximum parsimony phylogeny of C. caretta matrilines that sheds new light on the evolutionary relationships within the collection of lineages found in the Mediterranean and so far recognized by D-loop haplotypes only. We show that the new variants are useful to understand the ancestry of extant haplotypes, to improve genetically based studies on the philopatry and migratory behavior of the species, and for conservation purposes. In order to better understand the biological significance of the observed variation, we addressed intraspecific nonsynonymous substitutions in the context of the three-dimensional modeled structures of ND1 and ND3. The positions of variant amino acids within the folded subunits are consistent with a coadaptation with the restructuring of membrane thickness, fluidity, and lipid composition, a well-known response mechanism to thermal conditions. The pattern of amino acid substitutions departs from neutrality, suggesting local adaptation and/or polymorphism-based local selection.

Human forensic STRs used for individual identification have been reported to have little power for inter-population analyses. Several methods have been developed which incorporate information on the spatial distribution of individuals to arrive at a description of the arrangement of diversity. We genotyped at 16 forensic STRs a large population sample obtained from many locations in Italy, Greece and Turkey, i.e. three countries crucial to the understanding of discontinuities at the European/Asian junction and the genetic legacy of ancient migrations, but seldom represented together in previous studies. Using spatial PCA on the full dataset, we detected patterns of population affinities in the area. Additionally, we devised objective criteria to reduce the overall complexity into reduced datasets. Independent spatially explicit methods applied to these latter datasets converged in showing that the extraction of information on long- to medium-range geographical trends and structuring from the overall diversity is possible. All analyses returned the picture of a background clinal variation, with regional discontinuities captured by each of the reduced datasets. Several aspects of our results are confirmed on external STR datasets and replicate those of genome-wide SNP typings. High levels of gene flow were inferred within the main continental areas by coalescent simulations. These results are promising from a microevolutionary perspective, in view of the fast pace at which forensic data are being accumulated for many locales. It is foreseeable that this will allow the exploitation of an invaluable genotypic resource, assembled for other (forensic) purposes, to clarify important aspects in the formation of local gene pools.

Human forensic STRs are used for individual identification but have been reported to have little power for inter-population analyses. Several methods have been developed which incorporate information on the spatial distribution of individuals to arrive at a description of the arrangement of diversity. We genotyped at 16 forensic STRs a large population sample obtained from many locations in Italy, Greece and Turkey, i.e. three countries seldom represented together in previous studies. Using spatial PCA on the full dataset, we detected patterns of population affinities in the area similar to those of genome-wide SNP and STR studies. Additionally, we devised objective criteria to reduce the overall complexity into reduced datasets. Independent spatially explicit methods applied to these latter datasets converged in showing that the extraction of information on long-to medium-range geographical trends and structuring from the overall diversity is possible. All analyses returned the picture of a background clinal variation, with regional discontinuities captured by each of the reduced datasets. These coincided with the main bodies of water, i.e. the Adriatic/Ionian and the Aegean Seas. High levels of gene flow were inferred within the main continental areas by coalescent simulations. These results are promising in a microevolutionary perspective, in view of the fast pace at which forensic data are being accumulated for many locales. It is foreseeable that this will allow the exploitation of an invaluable genotypic resource, assembled for other (forensic) purposes, to clarify important aspects in the formation of local gene pools.

Huntington's disease is a neurodegenerative disorder caused by an expanded CAG trinucleotide repeat whose length is the major determinant of age at onset but remaining variation appears to be due in part to the effect of genetic modifiers. GRIK2, which encodes GluR6, a mediator of excitatory neurotransmission in the brain, has been suggested in several studies to be a modifier gene based upon a 3′ untranslated region TAA trinucleotide repeat polymorphism. Prior to investing in detailed studies of the functional impact of this polymorphism, we sought to confirm its effect on age at onset in a much larger dataset than in previous investigations. We genotyped the HD CAG repeat and the GRIK2 TAA repeat in DNA samples from 2,911 Huntington's disease subjects with known age at onset, and tested for a potential modifier effect of GRIK2 using a variety of statistical approaches. Unlike previous reports, we detected no evidence of an influence of the GRIK2 TAA repeat polymorphism on age at motor onset. Similarly, the GRIK2 polymorphism did not show significant modifier effect on psychiatric and cognitive age at onset in HD. Comprehensive analytical methods applied to a much larger sample than in previous studies do not support a role for GRIK2 as a genetic modifier of age at onset of clinical symptoms in Huntington's disease.

Huntington’s disease (HD) is an inherited neurodegenerative disorder characterized by motor, cognitive and behavioral disturbances, caused by the expansion of a CAG trinucleotide repeat in the HD gene. The CAG allele size is the major determinant of age at onset (AO) of motor symptoms, although the remaining variance in AO is highly heritable. The rs7665116 SNP in PPARGC1A, encoding the mitochondrial regulator PGC-1α, has been reported to be a significant modifier of AO in three European HD cohorts, perhaps due to affected cases from Italy. We attempted to replicate these findings in a large collection of (1,727) HD patient DNA samples of European origin. In the entire cohort, rs7665116 showed a significant effect in the dominant model (p value = 0.008) and the additive model (p value = 0.009). However, when examined by origin, cases of Southern European origin had an increased rs7665116 minor allele frequency (MAF), consistent with this being an ancestry-tagging SNP. The Southern European cases, despite similar mean CAG allele size, had a significantly older mean AO (p

Age at the onset of motor symptoms in Huntington disease (HD) is determined largely by the length of a CAG repeat expansion in HTT but is also influenced by other genetic factors. We tested whether common genetic variation near the mutation site is associated with differences in the distribution of expanded CAG alleles or age at the onset of motor symptoms. To define disease-associated single-nucleotide polymorphisms (SNPs), we compared 4p16.3 SNPs in HD subjects with population controls in a case:control strategy, which revealed that the strongest signals occurred at a great distance from the HD mutation as a result of “synthetic association” with SNP alleles that are of low frequency in population controls. Detailed analysis delineated a prominent ancestral haplotype that accounted for ∼50% of HD chromosomes and extended to at least 938 kb on about half of these. Together, the seven most abundant haplotypes accounted for ∼83% of HD chromosomes. Neither the extended shared haplotype nor the individual local HTT haplotypes were associated with altered CAG-repeat length distribution or residual age at the onset of motor symptoms, arguing against modification of these disease features by common cis-regulatory elements. Similarly, the 11 most frequent control haplotypes showed no trans-modifier effect on age at the onset of motor symptoms. Our results argue against common local regulatory variation as a factor influencing HD pathogenesis, suggesting that genetic modifiers be sought elsewhere in the genome. They also indicate that genome-wide association analysis with a small number of cases can be effective for regional localization of genetic defects, even when a founder effect accounts for only a fraction of the disorder.

Seventeen Y-chromosomal short tandem repeats (STRs) were analyzed in 347 healthy, unrelated, autochthonous males from the Andalusian provinces of Huelva ( N =167) and Granada ( N =180). AmpFlSTR Y-filer PCR Amplification kit (Applied Biosystems) was used to type the Y-STR markers. A total of 156 and 166 different haplotypes for the 17 Y-STR set were detected in Huelva, and Granada, respectively. The same haplotype diversity was found for both samples (0.998±0.001), and the overall discrimination capacity was 0.904. The most common minimal haplotype (DYS19, DYS389 I, DYS389 II, DYS390, DYS391, DYS392, DYS393) in both subpopulations was 14-13-16-24-11-13-13, which is also the most frequent haplotype among Atlantic European populations. Comparison analysis using pairwise R ST values and Analysis of Molecular Variance (AMOVA) revealed a significant genetic distance between our Andalusian samples and other ones from the northern Iberian fringe (including Basque and Pyrenean populations). However, results from the multi-dimensional scaling analysis (MDS) yielded a well-defined group of Iberian populations separated from the other Mediterranean clusters observed.

Objectives This work aimed to describe the genetic landscape of the Balkan Peninsula, as revealed by STR markers commonly used in forensics and spatial methods specifically developed for genetic data. Methods We generated and analyzed 16 short tandem repeats (STRs) autosomal genotypes in 287 subjects from ten administrative/geographical regions of Eastern Europe (Romania and the Republic of Moldova). We report estimates of the allele frequencies in these sub-populations, their fixation indexes, and use these results to complement previous spatial analyses of Southern Europe. Results In seven out of ten analyzed regional samples the heterozygosity, averaged across loci, was lower than expected. The average Fis was 0.011. Among the 16 loci, five returned a significant fixation index Fst. The composite Fst across the 16 loci, among the 10 regional samples, was 0.00417, a figure twice as large as that obtained with the same markers across the entire Northern Mediterranean. The first spatial principal component (sPC1) returned the picture of a Central-European pattern of frequencies for the Carpathians, which extended to the Southern boundary of the Balkan Peninsula. However, the 8 alleles extracted by sPC1 returned a picture of a strong reduction of the migration rate in the Carpathian region, mostly between the inner locations. Conclusions Our results revealed an unexpected heterogeneity in the area. We believe that populations from some regions will require treatment as distinct entities when considered in forensic applications.

Udgivelsesdato: 2008-Feb To investigate the genetic contribution to familial similarity in longevity, we set up a novel experimental design where cousin-pairs born from siblings who were concordant or discordant for the longevity trait were analyzed. To check this design, two chromosomal regions already known to encompass longevity-related genes were examined: 6p21.3 (genes TNFalpha, TNFbeta, HSP70.1) and 11p15.5 (genes SIRT3, HRAS1, IGF2, INS, TH). Population pools of 1.6, 2.3 and 2.0 million inhabitants were screened, respectively, in Denmark, France and Italy to identify families matching the design requirements. A total of 234 trios composed by one centenarian, his/her child and a child of his/her concordant or discordant sib were collected. By using population-specific allele frequencies, we reconstructed haplotype phase and estimated the likelihood of Identical By Descent (IBD) haplotype sharing in cousin-pairs born from concordant and discordant siblings. In addition, we analyzed haplotype transmission from centenarians to offspring, and a statistically significant Transmission Ratio Distortion (TRD) was observed for both chromosomal regions in the discordant families (P=0.007 for 6p21.3 and P=0.015 for 11p15.5). In concordant families, a marginally significant TRD was observed at 6p21.3 only (P=0.06). Although no significant difference emerged between the two groups of cousin-pairs, our study gave new insights on the hindrances to recruiting a suitable sample to obtain significant IBD data on longevity-related chromosomal regions. This will allow to dimension future sampling campaigns to study-genetic basis of human longevity.European Journal of Human Genetics (2008) 16, 236-242; doi:10.1038/sj.ejhg.5201950; published online 7 November 2007.

In questa rassegna vengono discusse soprattutto quelle caratteristiche genetiche che accomunano o differenziano l’uomo da orango, gorilla e scimpanze. Sono discussi i concetti di inquadramento gerarchico della diversita e di coalescenza e viene valutato il significato dell’evoluzione delle regioni codificanti per l’interpretazione dei dati sperimentali accumulatisi negli ultimi anni anche come risultato dei progetti di sequenziamento del genoma. I risultati sperimentali citati coprono tre aspetti della comparazione uomo-grandi scimmie: genetica di popolazioni, genomica comparata e profili di espressione genica. Gli studi di genetica delle popolazioni concordano nell’indicare, per le diverse sottospecie di orango, gorilla, e scimpanze, un grado di diversita genetica assai superiore a quello di qualsiasi popolazione umana, indicativo sia di un isolamento protrattosi per tempi lunghi che di una dimensione della popolazione piu grande di quella che ha dato origine alle popolazioni continentali umane. Il sequenziamento completo del genoma dello scimpanze ha rivelato una divergenza dell’ 1,23% rispetto all’uomo, per le regioni ortologhe, ovvero che possono essere direttamente confrontate. A questo si deve aggiungere il 3% relativo a regioni presenti o assenti in una delle due specie. Vengono illustrati alcuni esempi di geni per i quali e stato proposto un ruolo preminente per la comparsa di caratteristiche peculiari umane. Complessivamente, i dati indicano che nella storia delle linee evolutive umana e dello scimpanze, la maggioranza delle novita genetiche ha comportato cambiamenti soprattutto di carattere quantitativo di strutture e modalita di sviluppo gia esistenti. Infine, i profili di espressione genica sono stati analizzati per saggiare l’ipotesi di un’evoluzione principalmente regolativa. Anche in questo caso i risultati indicano la prevalenza di una selezione naturale contraria a cambiamenti drastici (purificante). Complessivamente questi dati indicano che l’evoluzione delle grandi scimmie e dell’uomo si e svolta secondo modalita e con velocita comuni anche agli altri mammiferi. Infine si discute come, considerando dati interspecifici e intra-specifici, diversi autori hanno messo in discussione l’identificazione di numerose specie nella linea evolutiva che ha condotto all’uomo an atomicamente moderno.