High resolution analysis and phylogenetic network construction using complete mtDNA sequences in sardinian genetic isolates

Helena's Many Daughters: More Mitogenome Diversity behind the Most Common West Eurasian mtDNA Control Region Haplotype in an Extended Italian Population Sample

The high number of matching haplotypes of the most common mitochondrial (mt)DNA lineages are considered to be the greatest limitation for forensic applications. This study investigates the potential to solve this constraint by massively parallel sequencing a large number of mitogenomes that share the most common West Eurasian mtDNA control region (CR) haplotype motif (263G 315.1C 16519C). We augmented a pilot study on 29 to a total of 216 Italian mitogenomes that represents the largest set of the most common CR haplotype compiled from a single country. The extended population sample confirmed and extended the huge coding region diversity behind the most common CR motif. Complete mitogenome sequencing allowed for the detection of 163 distinct haplotypes, raising the power of discrimination from 0 (CR) to 99.6% (mitogenome). The mtDNAs were clustered into 61 named clades of haplogroup H and did not reveal phylogeographic trends within Italy. Rapid individualization approaches for investigative purposes are limited to the most frequent H clades of the dataset, viz. H1, H3, and H7.

Bioarchaeological and palaeogenomic portrait of two Pompeians that died during the eruption of Vesuvius in 79 AD

The archaeological site of Pompeii is one of the 54 UNESCO World Heritage sites in Italy, thanks to its uniqueness: the town was completely destroyed and buried by a Vesuvius’ eruption in 79 AD. In this work, we present a multidisciplinary approach with bioarchaeological and palaeogenomic analyses of two Pompeian human remains from the Casa del Fabbro. We have been able to characterize the genetic profile of the first Pompeian’ genome, which has strong affinities with the surrounding central Italian population from the Roman Imperial Age. Our findings suggest that, despite the extensive connection between Rome and other Mediterranean populations, a noticeable degree of genetic homogeneity exists in the Italian peninsula at that time. Moreover, palaeopathological analyses identified the presence of spinal tuberculosis and we further investigated the presence of ancient DNA from Mycobacterium tuberculosis. In conclusion, our study demonstrates the power of a combined approach to investigate ancient humans and confirms the possibility to retrieve ancient DNA from Pompeii human remains. Our initial findings provide a foundation to promote an intensive and extensive paleogenetic analysis in order to reconstruct the genetic history of population from Pompeii, a unique archaeological site.

From old markers to next generation: reconstructing the history of the peopling of Sardinia

CONTEXT

For many years the Sardinian population has been the object of numerous studies because of its unique genetic structure. Despite the extreme abundance of papers, various aspects of the peopling and genetic structure of Sardinia still remain uncertain and sometimes controversial.

OBJECTIVE

We reviewed what has emerged from different studies, focussing on some still open questions, such as the origin of Sardinians, their relationship with the Corsican population, and the intra-regional genetic heterogeneity.

METHODS

The various issues have been addressed through the analysis of classical markers, molecular markers and, finally, genomic data through next generation sequencing.

RESULTS AND CONCLUSIONS

Although the most ancient human remains date back to the end of the Palaeolithic, Mesolithic populations brought founding lineages that left evident traces in the modern population. Then, with the Neolithic, the island underwent an important demographic expansion. Subsequently, isolation and genetic drift contributed to maintain a significant genetic heterogeneity, but preserving the overall homogeneity on a regional scale. At the same time, isolation and genetic drift contributed to differentiate Sardinia from Corsica, which saw an important gene flow from the mainland. However, the isolation did not prevent gene flow from the neighbouring populations whose contribution are still recognisable in the genome of Sardinians.

Complete mitogenome data for the Serbian population: the contribution to high-quality forensic databases

Mitochondrial genome (mtDNA) is a valuable resource in resolving various human forensic casework. The usage of variability of complete mtDNA genomes increases their discriminatory power to the maximum and enables ultimate resolution of distinct maternal lineages. However, their wider employment in forensic casework is nowadays limited by the lack of appropriate reference database. In order to fill in the gap in the reference data, which, considering Slavic-speaking populations, currently comprises only mitogenomes of East and West Slavs, we present mitogenome data for 226 Serbians, representatives of South Slavs from the Balkan Peninsula. We found 143 (sub)haplogroups among which West Eurasian ones were dominant. The percentage of unique haplotypes was 85%, and the random match probability was as low as 0.53%. We support previous findings on both high levels of genetic diversity in the Serbian population and patterns of genetic differentiation among this and ten studied European populations. However, our high-resolution data supported more pronounced genetic differentiation among Serbians and two Slavic populations (Russians and Poles) as well as expansion of the Serbian population after the Last Glacial Maximum and during the Migration period (fourth to ninth century A.D.), as inferred from the Bayesian skyline analysis. Phylogenetic analysis of haplotypes found in Serbians contributed towards the improvement of the worldwide mtDNA phylogeny, which is essential for the interpretation of the mtDNA casework.

Whole mitochondrial genome diversity in two Hungarian populations

Complete mitochondrial genomics is an effective tool for studying the demographic history of human populations, but there is still a deficit of mitogenomic data in European populations. In this paper, we present results of study of variability of 80 complete mitochondrial genomes in two Hungarian populations from eastern part of Hungary (Szeged and Debrecen areas). The genetic diversity of Hungarian mitogenomes is remarkably high, reaching 99.9% in a combined sample. According to the analysis of molecular variance (AMOVA), European populations showed a low, but statistically significant level of between-population differentiation (Fst = 0.61%, p = 0), and two Hungarian populations demonstrate lack of between-population differences. Phylogeographic analysis allowed us to identify 71 different mtDNA sub-clades in Hungarians, sixteen of which are novel. Analysis of ancestry-informative mtDNA sub-clades revealed a complex genetic structure associated with the genetic impact of populations from different parts of Eurasia, though the contribution from European populations is the most pronounced. At least 8% of ancestry-informative haplotypes found in Hungarians demonstrate similarity with East and West Slavic populations (sub-clades H1c23a, H2a1c1, J2b1a6, T2b25a1, U4a2e, K1c1j, and I1a1c), while the influence of Siberian populations is not so noticeable (sub-clades A12a, C4a1a, and probably U4b1a4).

Mitogenomic diversity in Russians and Poles

Complete mtDNA genome sequencing improves molecular resolution for distinguishing variation between individuals and populations, but there is still deficiency of mitogenomic population data. To overcome this limitation, we used Sanger-based protocol to generate complete mtDNA sequences of 376 Russian individuals from six populations of European part of Russia and 100 Polish individuals from northern Poland. Nearly complete resolution of mtDNA haplotypes was achieved - about 97% of haplotypes were unique both in Russians and Poles, and no haplotypes overlapped between them when indels were considered. While European populations showed a low, but statistically significant level of between-population differentiation (Fst=0.66%, p=0), Russians demonstrate lack of between-population differences (Fst=0.22%, p=0.15). Results of the Bayesian skyline analysis of Russian mitogenomes demonstrate not only post-Last Glacial Maximum expansion, but also rapid population growth starting from about 4.3kya (95% CI: 2.9-5.8kya), i.e. in the Bronze Age. This expansion strongly correlates with the Kurgan model established by archaeologists and confirmed by paleogeneticists.

Phylogeography, genetic diversity and demographic history of the Iranian Kurdish groups based on mtDNA sequences

Throughout the history of modern humans, the current Kurdish-inhabited area has served as part of a tricontinental crossroad for major human migrations. Also, a significant body of archaeological evidence points to this area as the site of Neolithic transition. To investigate the phylogeography, origins and demographic history, mtDNA D-loop region of individuals representing four Kurdish groups from Iran were analysed. Our data indicated that most of the Kurds mtDNA lineages belong to branches of the haplogroups with the Western Eurasian origin; with small fractions of the Eastern Eurasian and sub-Saharan African lineages. The low level of mtDNA diversity observed in the Havrami group presented a bias towards isolation or increased drift due to small population size; while in the Kurmanji group it indicated a bias towards drift or mass migration events during the 5-18th century AD. The Mantel test showed strong isolation by distance, and AMOVA results for global and regional scales confirmed that the geography had acted as the main driving force in shaping the current pattern of mtDNA diversity, rather than linguistic similarity. The results of demographic analyses, in agreement with archaeological data, revealed a recent expansion of the Kurds (~9,500 years before present) related to the Neolithic transition from hunting and gathering, to farming and cattle breeding in the Near East. Further, the high frequencies of typical haplogroups for early farmers (H; 37.1%) and hunter-gatherers (U; 13.8%) in the Kurds may testify the earlier hunter-gatherers in the Kurdish-inhabited area that adopted and admixed the Kurds ancestors following the Neolithic transition.

Linking between genetic structure and geographical distance: Study of the maternal gene pool in the Ethiopian population

Background The correlation between genetics and geographical distance has already been examined through the study of the dispersion of human populations, especially in terms of uniparental genetic markers. Aim The present work characterises, at the level of the mitochondrial DNA (mtDNA), two new samples of Amhara and Oromo populations from Ethiopia to evaluate the possible pattern of distribution for mtDNA variation and to test the hypothesis of the Isolation-by-Distance (IBD) model among African, European and Middle-Eastern populations. Subjects and methods This study analysed 173 individuals belonging to two ethnic groups of Ethiopia, Amhara and Oromo, by assaying HVS-I and HVS-II of mtDNA D-loop and informative coding region SNPs of mtDNA. Results The analysis suggests a relationship between genetic and geographic distances, affirming that the mtDNA pool of Africa, Europe and the Middle East might be coherent with the IBD model. Moreover, the mtDNA gene pools of the Sub-Saharan African and Mediterranean populations were very different. Conclusion In this study the pattern of mtDNA distribution, beginning with the Ethiopian plateau, was tested in the IBD model. It could be affirmed that, on a continent scale, the mtDNA pool of Africa, Europe and the Middle East might fall under the IBD model.

Dissecting metabolic syndrome components: data from an epidemiologic survey in a genetic isolate

The metabolic syndrome (MetS) is a large-scale and expanding public-health and clinical threat worldwide. We investigated the determinants of MetS, assessed its prevalence and components and, estimated their genetic contribution, taking advantage of the special characteristics of Sardinian isolated populations. Inhabitants of 10 villages in Ogliastra region participated in a cross-sectional survey in 2002–2008 (n = 9,647). Blood samples, blood pressure (BP), anthropometry and, data from a standardized interview were collected. Prevalence of MetS was estimated by the direct method of standardization. Variables associated with the MetS were identified using multilevel logistic regression. Heritability was determined using variance component models. MetS Prevalence was 19.6% (95% CI 18.9–20.4%) according to NCEP-ATPIII, 24.8% (95% CI 24.0–25.6%) according to IDF and, 29% (95% CI 28.1–29.8%) according to AHA/NHLBI harmonized criteria, ranging from 9 to 26% among villages. The most prevalent combination was BP + HDL-cholesterol (HDL) + triglycerides (TRIG) (19%), followed by BP + HDL + waist circumference (WAIST) (17%) and, BP + HDL + TRIG + WAIST (13.6%). Heritability of MetS was 48% (p = 1.62 × 10⁻²⁵), as the two most common combinations (BP + HDL + TRIG and BP + HDL + WAIST) showed heritability of 53 and 52%, respectively. The larger genetic components of the two most frequent combinations determining MetS deserve greater investigation in order to understand the underlying mechanisms. Besides, further studies are warranted to confirm these findings both in isolated and outbred populations.

Ethnopedology in the study of toponyms connected to the indigenous knowledge on soil resource

In taking an integrated ethnopedological approach, this study aims to investigate the meaning of the distribution of the toponyms used in traditional and recent cartography of Sardinia (southern Italy). It is particularly, but not only, focused on those related to soil resources. Sardinia is particularly interesting in this respect, as its unique history, geography, and linguistic position makes it one of the Italian and Mediterranean regions with the greatest number of toponyms. This research investigated the toponyms belonging to an important sub-region of Sardinia, called Ogliastra (central-eastern Sardinia). The research was conducted through the following integrated approach: i) toponymy research and collection from different sources; ii) database creation and translation of toponyms from the Sardinian language (SL); iii) categorization of toponyms; and iv) graphical, statistical, and cartographic data processing. Distribution and diversity of toponyms were assessed using the compiled database, coupled with a geographical information system (GIS). Of around 7700 toponyms collected, 79% had already been reported in SL, while just 21% were in Italian. Of the toponyms in SL, 77% are of known meaning and 54% of these toponyms were characterized by a meaning directly and/or indirectly connected to specific environmental features. On the whole, morphology would appear to be the primary environmental factor able to explain the complex, articulated presence, distribution, and typology of the investigated toponyms. A least squares regression analysis of pedodiversity vs. topodiversity shows a very closed distribution, with an impressive high correlation index (R² = 0.824). The principal factor analysis (PFA) shows that such a connection may be morphologically based, thereby confirming that pedodiversity and topodiversity are strongly linked each other. Overall, the research shows that an integrated ethnopedological approach, combining indigenous and scientific knowledge may be of great interest in order to mitigate the impressive phenomena of the indigenous knowledge lost.

Traces of forgotten historical events in mountain communities in Central Italy: A genetic insight

OBJECTIVES

Analysis of human genetic variation in mountain communities can shed light on the peopling of mountainous regions, perhaps revealing whether the remote geographic location spared them from outside invasion and preserved their gene pool from admixture. In this study, we created a model to assess genetic traces of historical events by reconstructing the paternal and maternal genetic history of seven small mountain villages in inland valleys of Central Italy.

METHODS

The communities were selected for their geographic isolation, attested biodemographic stability, and documented history prior to the Roman conquest. We studied the genetic structure by analyzing two hypervariable segments (HVS-I and HVS-II) of the mtDNA D-loop and several informative single nucleotide polymorphisms (SNPs) of the mtDNA coding region in 346 individuals, in addition to 17 short tandem repeats (STRs) and Y-chromosome SNPs in 237 male individuals.

RESULTS

For both uniparental markers, most of the haplogroups originated in Western Europe while some Near Eastern haplogroups were identified at low frequencies. However, there was an evident genetic similarity between the Central Italian samples and Near Eastern populations mainly in the male genetic pool.

CONCLUSIONS

The samples highlight an overall European genetic pattern both for mtDNA and Y chromosome. Notwithstanding this scenario, Y chromosome haplogroup Q, a common paternal lineage in Central/Western Asia but almost Europe-wide absent, was found, suggesting that Central Italy could have hosted a settlement from Anatolia that might be supported by cultural, topographic and genetic evidence.

Helena, the hidden beauty: Resolving the most common West Eurasian mtDNA control region haplotype by massively parallel sequencing an Italian population sample

The analysis of mitochondrial (mt)DNA is a powerful tool in forensic genetics when nuclear markers fail to give results or maternal relatedness is investigated. The mtDNA control region (CR) contains highly condensed variation and is therefore routinely typed. Some samples exhibit an identical haplotype in this restricted range. Thus, they convey only weak evidence in forensic queries and limited phylogenetic information. However, a CR match does not imply that also the mtDNA coding regions are identical or samples belong to the same phylogenetic lineage. This is especially the case for the most frequent West Eurasian CR haplotype 263G 315.1C 16519C, which is observed in various clades within haplogroup H and occurs at a frequency of 3-4% in many European populations. In this study, we investigated the power of massively parallel complete mtGenome sequencing in 29 Italian samples displaying the most common West Eurasian CR haplotype - and found an unexpected high diversity. Twenty-eight different haplotypes falling into 19 described sub-clades of haplogroup H were revealed in the samples with identical CR sequences. This study demonstrates the benefit of complete mtGenome sequencing for forensic applications to enforce maximum discrimination, more comprehensive heteroplasmy detection, as well as highest phylogenetic resolution.

The episode of genetic drift defining the migration of humans out of Africa is derived from a large east African population size

Human genetic variation particularly in Africa is still poorly understood. This is despite a consensus on the large African effective population size compared to populations from other continents. Based on sequencing of the mitochondrial Cytochrome C Oxidase subunit II (MT-CO2), and genome wide microsatellite data we observe evidence suggesting the effective size (Ne) of humans to be larger than the current estimates, with a foci of increased genetic diversity in east Africa, and a population size of east Africans being at least 2-6 fold larger than other populations. Both phylogenetic and network analysis indicate that east Africans possess more ancestral lineages in comparison to various continental populations placing them at the root of the human evolutionary tree. Our results also affirm east Africa as the likely spot from which migration towards Asia has taken place. The study reflects the spectacular level of sequence variation within east Africans in comparison to the global sample, and appeals for further studies that may contribute towards filling the existing gaps in the database. The implication of these data to current genomic research, as well as the need to carry out defined studies of human genetic variation that includes more African populations; particularly east Africans is paramount.

Sardinians genetic background explained by runs of homozygosity and genomic regions under positive selection

The peculiar position of Sardinia in the Mediterranean sea has rendered its population an interesting biogeographical isolate. The aim of this study was to investigate the genetic population structure, as well as to estimate Runs of Homozygosity and regions under positive selection, using about 1.2 million single nucleotide polymorphisms genotyped in 1077 Sardinian individuals. Using four different methods--fixation index, inflation factor, principal component analysis and ancestry estimation--we were able to highlight, as expected for a genetic isolate, the high internal homogeneity of the island. Sardinians showed a higher percentage of genome covered by RoHs>0.5 Mb (F(RoH%0.5)) when compared to peninsular Italians, with the only exception of the area surrounding Alghero. We furthermore identified 9 genomic regions showing signs of positive selection and, we re-captured many previously inferred signals. Other regions harbor novel candidate genes for positive selection, like TMEM252, or regions containing long non coding RNA. With the present study we confirmed the high genetic homogeneity of Sardinia that may be explained by the shared ancestry combined with the action of evolutionary forces.

Demographic histories, isolation and social factors as determinants of the genetic structure of Alpine linguistic groups

Great European mountain ranges have acted as barriers to gene flow for resident populations since prehistory and have offered a place for the settlement of small, and sometimes culturally diverse, communities. Therefore, the human groups that have settled in these areas are worth exploring as an important potential source of diversity in the genetic structure of European populations. In this study, we present new high resolution data concerning Y chromosomal variation in three distinct Alpine ethno-linguistic groups, Italian, Ladin and German. Combining unpublished and literature data on Y chromosome and mitochondrial variation, we were able to detect different genetic patterns. In fact, within and among population diversity values observed vary across linguistic groups, with German and Italian speakers at the two extremes, and seem to reflect their different demographic histories. Using simulations we inferred that the joint effect of continued genetic isolation and reduced founding group size may explain the apportionment of genetic diversity observed in all groups. Extending the analysis to other continental populations, we observed that the genetic differentiation of Ladins and German speakers from Europeans is comparable or even greater to that observed for well known outliers like Sardinian and Basques. Finally, we found that in south Tyroleans, the social practice of Geschlossener Hof, a hereditary norm which might have favored male dispersal, coincides with a significant intra-group diversity for mtDNA but not for Y chromosome, a genetic pattern which is opposite to those expected among patrilocal populations. Together with previous evidence regarding the possible effects of “local ethnicity” on the genetic structure of German speakers that have settled in the eastern Italian Alps, this finding suggests that taking socio-cultural factors into account together with geographical variables and linguistic diversity may help unveil some yet to be understood aspects of the genetic structure of European populations.

Complete mitochondrial DNA diversity in Iranians

Due to its pivotal geographical location and proximity to transcontinental migratory routes, Iran has played a key role in subsequent migrations, both prehistoric and historic, between Africa, Asia and Europe. To shed light on the genetic structure of the Iranian population as well as on the expansion patterns and population movements which affected this region, the complete mitochondrial genomes of 352 Iranians were obtained. All Iranian populations studied here exhibit similarly high diversity values comparable to the other groups from the Caucasus, Anatolia and Europe. The results of AMOVA and MDS analyses did not associate any regional and/or linguistic group of populations in the Anatolia/Caucasus and Iran region pointing to close genetic positions of Persians and Qashqais to each other and to Armenians, and Azeris from Iran to Georgians. By reconstructing the complete mtDNA phylogeny of haplogroups R2, N3, U1, U3, U5a1g, U7, H13, HV2, HV12, M5a and C5c we have found a previously unexplored genetic connection between the studied Iranian populations and the Arabian Peninsula, India, Near East and Europe, likely the result of both ancient and recent gene flow. Our results for Persians and Qashqais point to a continuous increase of the population sizes from ∼24 kya to the present, although the phase between 14-24 kya is thought to be hyperarid according to the Gulf Oasis model. Since this would have affected hunter-gatherer ranges and mobility patterns and forced them to increasingly rely on coastal resources, this transition can explain the human expansion across the Persian Gulf region.

MT-CYB mutations in hypertrophic cardiomyopathy

Mitochondrial dysfunction is a characteristic of heart failure. Mutations in mitochondrial DNA, particularly in MT-CYB coding for cytochrome B in complex III (CIII), have been associated with isolated hypertrophic cardiomyopathy (HCM). We hypothesized that MT-CYB mutations might play an important causal or modifying role in HCM. The MT-CYB gene was sequenced from DNA isolated from blood from 91 Danish HCM probands. Nonsynonymous variants were analyzed by bioinformatics, molecular modeling and simulation. Two germline-inherited, putative disease-causing, nonsynonymous variants: m.15024G>A; p.C93Y and m.15482T>C; p.S246P were identified. Modeling showed that the p.C93Y mutation leads to disruption of the tertiary structure of Cytb by helix displacement, interfering with protein–heme interaction. The p.S246P mutation induces a diproline structure, which alters local secondary structure and induces a kink in the protein backbone, interfering with macromolecular interactions. These molecular effects are compatible with a leaky phenotype, that is, limited but progressive mitochondrial dysfunction. In conclusion, we find that rare, putative leaky mtDNA variants in MT-CYB can be identified in a cohort of HCM patients. We propose that further patients with HCM should be examined for mutations in MT-CYB in order to clarify the role of these variants.

Environmental and genetic contribution to hypertension prevalence: data from an epidemiological survey on Sardinian genetic isolates

Background and Objectives

Hypertension represents a major cause of cardiovascular morbidity and mortality worldwide but its prevalence has been shown to vary in different countries. The reasons for such differences are still matter of debate, the relative contributions given by environmental and genetic factors being still poorly defined. We estimated the current prevalence, distribution and determinants of hypertension in isolated Sardinian populations and also investigated the environmental and genetic contribution to hypertension prevalence taking advantage of the characteristics of such populations.

Methods and Results

An epidemiological survey with cross-sectional design was carried out measuring blood pressure in 9845 inhabitants of 10 villages of Ogliastra region between 2002 and 2008. Regression analysis for assessing blood pressure determinants and variance component models for estimating heritability were performed. Overall 38.8% of this population had hypertension, its prevalence varying significantly by age, sex and among villages taking into account age and sex structure of their population. About 50% of hypertensives had prior cardiovascular disease. High blood pressure was independently associated with age, obesity related factors, heart rate, total cholesterol, alcohol consumption, low education and smoking status, all these factors contributing more in women than in men. Heritability was 27% for diastolic and 36% for systolic blood pressure, its contribution being significantly higher in men (57%) than in women (46%). Finally, the genetic correlation between systolic and diastolic blood pressure was 0.74, indicating incomplete pleiotropy.

Conclusion

Genetic factors involved in the expression of blood pressure traits account for about 30% of the phenotypic variance, but seem to play a larger role in men; comorbidities and environmental factors remain of predominant importance, but seem to contribute much more in women.

Divorcing the Late Upper Palaeolithic demographic histories of mtDNA haplogroups M1 and U6 in Africa

Background

A Southwest Asian origin and dispersal to North Africa in the Early Upper Palaeolithic era has been inferred in previous studies for mtDNA haplogroups M1 and U6. Both haplogroups have been proposed to show similar geographic patterns and shared demographic histories.

Results

We report here 24 M1 and 33 U6 new complete mtDNA sequences that allow us to refine the existing phylogeny of these haplogroups. The resulting phylogenetic information was used to genotype a further 131 M1 and 91 U6 samples to determine the geographic spread of their sub-clades. No southwest Asian specific clades for M1 or U6 were discovered. U6 and M1 frequencies in North Africa, the Middle East and Europe do not follow similar patterns, and their sub-clade divisions do not appear to be compatible with their shared history reaching back to the Early Upper Palaeolithic. The Bayesian Skyline Plots testify to non-overlapping phases of expansion, and the haplogroups’ phylogenies suggest that there are U6 sub-clades that expanded earlier than those in M1. Some M1 and U6 sub-clades could be linked with certain events. For example, U6a1 and M1b, with their coalescent ages of ~20,000–22,000 years ago and earliest inferred expansion in northwest Africa, could coincide with the flourishing of the Iberomaurusian industry, whilst U6b and M1b1 appeared at the time of the Capsian culture.

Conclusions

Our high-resolution phylogenetic dissection of both haplogroups and coalescent time assessments suggest that the extant main branching pattern of both haplogroups arose and diversified in the mid-later Upper Palaeolithic, with some sub-clades concomitantly with the expansion of the Iberomaurusian industry. Carriers of these maternal lineages have been later absorbed into and diversified further during the spread of Afro-Asiatic languages in North and East Africa.

Results of Mitochondrial DNA Sequence Analysis in Patients with Clinically Diagnosed Leber's Hereditary Optic Neuropathy

OBJECTIVE

To investigate possible mitochondrial DNA (mtDNA) mutations in patients with Leber's hereditary optic neuropathy (LHON) in order to provide a precise diagnosis and genetic counseling.

MATERIAL AND METHODS

Between 1982 and 2007, ten patients were clinically diagnosed with LHON and six of these patients agreed to be involved in this study. Six healthy individuals were also included as a control group. mtDNA was isolated from peripheral blood samples and polymerase chain reaction and mtDNA sequence analysis were performed.

RESULTS

In one of the six patients, a homoplasmic mutant m.11778G>A mutation was detected. All of the clinically diagnosed LHON patients and the control groups had the m.14212C>T and m.14580G>A single nucleotide polymorphisms (SNPs). The m.11719A>G SNP was detected in three of six patients and four of the controls. Two of the six patients had the m.3197T>C SNP and, in addition, the m.14258G>A SNP was found in one of these two patients, while neither of these mutations were present in the control group.

CONCLUSION

The clinical diagnosis of LHON could be supported by molecular genetics only in one patient by the detection of one mutation. The m.3197T>C and m.14258G>A SNPs should be considered as potential mtDNA mutations due to the fact that they were detected in the patient group. These mutations should be investigated further in large case groups for suspected gene loci that could lead to optic neuropathy.

Evidence of high genetic variation among linguistically diverse populations on a micro-geographic scale: a case study of the Italian Alps

Although essential for the fine-scale reconstruction of genetic structure, only a few micro-geographic studies have been carried out in European populations. This study analyzes mitochondrial variation (651 bp of the hypervariable region plus 17 single-nucleotide polymorphisms) in 393 samples from nine populations from Trentino (Eastern Italian Alps), a small area characterized by a complex geography and high linguistic diversity. A high level of genetic variation, comparable to geographically dispersed European groups, was observed. We found a difference in the intensity of peopling processes between two longitudinal areas, as populations from the west-central part of the region show stronger signatures of expansion, whereas those from the eastern area are closer to the expectations of a stationary demographic state. This may be explained by geomorphological factors and is also supported by archeological data. Finally, our results reveal a striking difference in the way in which the two linguistically isolated populations are genetically related to the neighboring groups. The Ladin speakers were found to be genetically close to the Italian-speaking populations and differentiated from the other Dolomitic Ladins, whereas the German-speaking Cimbri behave as an outlier, showing signatures of founder effects and low growth rate.

Complete mitochondrial genomes reveal neolithic expansion into Europe

The Neolithic transition from hunting and gathering to farming and cattle breeding marks one of the most drastic cultural changes in European prehistory. Short stretches of ancient mitochondrial DNA (mtDNA) from skeletons of pre-Neolithic hunter-gatherers as well as early Neolithic farmers support the demic diffusion model where a migration of early farmers from the Near East and a replacement of pre-Neolithic hunter-gatherers are largely responsible for cultural innovation and changes in subsistence strategies during the Neolithic revolution in Europe. In order to test if a signal of population expansion is still present in modern European mitochondrial DNA, we analyzed a comprehensive dataset of 1,151 complete mtDNAs from present-day Europeans. Relying upon ancient DNA data from previous investigations, we identified mtDNA haplogroups that are typical for early farmers and hunter-gatherers, namely H and U respectively. Bayesian skyline coalescence estimates were then used on subsets of complete mtDNAs from modern populations to look for signals of past population expansions. Our analyses revealed a population expansion between 15,000 and 10,000 years before present (YBP) in mtDNAs typical for hunters and gatherers, with a decline between 10,000 and 5,000 YBP. These corresponded to an analogous population increase approximately 9,000 YBP for mtDNAs typical of early farmers. The observed changes over time suggest that the spread of agriculture in Europe involved the expansion of farming populations into Europe followed by the eventual assimilation of resident hunter-gatherers. Our data show that contemporary mtDNA datasets can be used to study ancient population history if only limited ancient genetic data is available.

Sampling strategies in a linguistic isolate: results from mtDNA analysis

OBJECTIVES

Sampling strategies are crucial issues in population genetics and anthropological studies. The sampling choice is related to the research question and the type of markers used. In this research, we compared two different sampling strategies in the Sardinian linguistic isolate of Carloforte (Italy).

METHODS

A first sampling (N = 49) was carried out through grandparents criterion: individuals selected for the study were born and resident in Carloforte, and unrelated for at least three generations. A second sampling (N = 50) was based on founders surnames (FS): selected participants were proved to be descendants of the village founders, and to have no ancestors in common, at least up to the grandparental generation.

RESULTS

The group selected through FS showed a greater gene diversity, which was confirmed by both network and haplogroup analysis. Among the shared haplogroups, we find clear differences in their frequencies. Sampling through grandparents criterion showed essentially the same haplogroups found in Sardinia, and with similar frequencies. Interesting results came from genetic tree. The FS sampling clustered with Northern African populations and it is located very far from Italian and Sardinian populations, whereas the grandparents criterion sampling clustered with Italian populations and it is located close to the other Sardinian populations.

CONCLUSIONS

Results showed that different sampling strategies can lead to contrasting results. As sampling through grandparents criterion is influenced by recent gene flow, we hypothesize that the difference observed with the two sampling strategies is due to the merging of Carloforte with Sardinian populations.

Mitochondrial Cytochrome c Oxidase subunit 1 Sequence Variation in Prostate Cancer

PURPOSE

Mitochondrial DNA (mtDNA) gene mutations have been described in nearly every adult solid neoplasm including prostate cancer. There are marked racial differences in specific inherited mutations within the cytochrome c oxidase subunit 1 (COI) gene in individuals with prostate cancer (PCa). The purpose of this study was to identify the variation in COI gene sequence in prostate cancer patients and to compare the mutations in African and Caucasian Americans.

MATERIALS AND METHODS

We sequenced the COI gene in DNA derived from peripheral blood in 482 prostate cancer patients and 189 controls. All bases that differed from the revised Cambridge Reference Sequence (rCRS) were classified as either silent (non-amino acid altering) or missense (amino acid altering) and the compiled alterations were then compared between races and published reports of mutations in this gene in both Caucasian and African-Americans.

RESULTS AND CONCLUSIONS

We found inherited mtDNA COI missense variants in 8.8% of Caucasian prostate cancer patients (vs. 0.0% controls) and 72.8 % of African-American prostate cancer patients (vs. 64.3% controls) A total of 144 COI variants were identified, of which 30 were missense mutations. Of 482 PCa patients, 116 (24.1%) had one or more missense mutations. Further evaluation of this gene and these mutations may allow for the identification of genetically at-risk populations. The high rate of COI mutations in African-Americans may account for some of the racial disparity observed in prostate cancer.

Mendelian breeding units versus standard sampling strategies: Mitochondrial DNA variation in southwest Sardinia

We report a sampling strategy based on Mendelian Breeding Units (MBUs), representing an interbreeding group of individuals sharing a common gene pool. The identification of MBUs is crucial for case-control experimental design in association studies. The aim of this work was to evaluate the possible existence of bias in terms of genetic variability and haplogroup frequencies in the MBU sample, due to severe sample selection. In order to reach this goal, the MBU sampling strategy was compared to a standard selection of individuals according to their surname and place of birth. We analysed mitochondrial DNA variation (first hypervariable segment and coding region) in unrelated healthy subjects from two different areas of Sardinia: the area around the town of Cabras and the western Campidano area. No statistically significant differences were observed when the two sampling methods were compared, indicating that the stringent sample selection needed to establish a MBU does not alter original genetic variability and haplogroup distribution. Therefore, the MBU sampling strategy can be considered a useful tool in association studies of complex traits.

Application of a new method for GWAS in a related case/control sample with known pedigree structure: identification of new loci for nephrolithiasis

In contrast to large GWA studies based on thousands of individuals and large meta-analyses combining GWAS results, we analyzed a small case/control sample for uric acid nephrolithiasis. Our cohort of closely related individuals is derived from a small, genetically isolated village in Sardinia, with well-characterized genealogical data linking the extant population up to the 16(th) century. It is expected that the number of risk alleles involved in complex disorders is smaller in isolated founder populations than in more diverse populations, and the power to detect association with complex traits may be increased when related, homogeneous affected individuals are selected, as they are more likely to be enriched with and share specific risk variants than are unrelated, affected individuals from the general population. When related individuals are included in an association study, correlations among relatives must be accurately taken into account to ensure validity of the results. A recently proposed association method uses an empirical genotypic covariance matrix estimated from genome-screen data to allow for additional population structure and cryptic relatedness that may not be captured by the genealogical data. We apply the method to our data, and we also investigate the properties of the method, as well as other association methods, in our highly inbred population, as previous applications were to outbred samples. The more promising regions identified in our initial study in the genetic isolate were then further investigated in an independent sample collected from the Italian population. Among the loci that showed association in this study, we observed evidence of a possible involvement of the region encompassing the gene LRRC16A, already associated to serum uric acid levels in a large meta-analysis of 14 GWAS, suggesting that this locus might lead a pathway for uric acid metabolism that may be involved in gout as well as in nephrolithiasis.

Population expansion in the North African late Pleistocene signalled by mitochondrial DNA haplogroup U6

BACKGROUND

The archaeology of North Africa remains enigmatic, with questions of population continuity versus discontinuity taking centre-stage. Debates have focused on population transitions between the bearers of the Middle Palaeolithic Aterian industry and the later Upper Palaeolithic populations of the Maghreb, as well as between the late Pleistocene and Holocene.

RESULTS

Improved resolution of the mitochondrial DNA (mtDNA) haplogroup U6 phylogeny, by the screening of 39 new complete sequences, has enabled us to infer a signal of moderate population expansion using Bayesian coalescent methods. To ascertain the time for this expansion, we applied both a mutation rate accounting for purifying selection and one with an internal calibration based on four approximate archaeological dates: the settlement of the Canary Islands, the settlement of Sardinia and its internal population re-expansion, and the split between haplogroups U5 and U6 around the time of the first modern human settlement of the Near East.

CONCLUSIONS

A Bayesian skyline plot placed the main expansion in the time frame of the Late Pleistocene, around 20 ka, and spatial smoothing techniques suggested that the most probable geographic region for this demographic event was to the west of North Africa. A comparison with U6's European sister clade, U5, revealed a stronger population expansion at around this time in Europe. Also in contrast with U5, a weak signal of a recent population expansion in the last 5,000 years was observed in North Africa, pointing to a moderate impact of the late Neolithic on the local population size of the southern Mediterranean coast.

History, geography and population structure influence the distribution and heritability of blood and anthropometric quantitative traits in nine Sardinian genetic isolates

Isolated founder populations which exhibit great genetic and environmental homogeneity provide an attractive setting for the study of quantitative traits (QTs). Geneticists have repeatedly turned to population isolates and the past successes have prompted increased interest among medical researchers. We studied nine small isolated villages of a secluded area of Sardinia (Ogliastra), all of them characterized by a few founders, high endogamy rates, slow population expansion and a distinct genetic makeup. Anthropometric and blood parameters, 43 QTs in all, were analysed in about 9000 voluntary subjects for whom extended genealogical information was available. We explored the distribution and examined mean differences of each trait among villages by analysis of variance (ANOVA). A heritability analysis with the variance component (VC) method was performed. Results show significant differences in the distribution of most traits between groups of villages located in two distinct geographical areas already identified by a previous population structure analysis, thus supporting the existence of differentiation among sub-populations in the same region. Heritability estimates range between 30 and 89%, demonstrating that genetic effects substantially contribute to phenotypic variation of all investigated traits and that this population provides excellent research conditions for gene-mapping projects. Results suggest that history, geographic location and population structure may have influenced the genetic and phenotypic features of these isolates. Our findings may be useful for the ongoing linkage and association studies in these isolates and suggest that a thorough characterization of population is valuable to better identify genes or variants that may be rare in the population at large and peculiar to single villages.

Craniofacial morphometric variation and the biological history of the peopling of Sardinia

The aim of this work is to explore the pattern of craniofacial morphometric variation and the relationships among five prehistoric Sardinian groups dated from Late Neolithic to the Nuragic Period (Middle and Late Bronze Age), in order to formulate hypotheses on the peopling history of Sardinia. Biological relationships with coeval populations of central peninsular Italy were also analysed to detect influences from and towards extra-Sardinian sources. Furthermore, comparison with samples of contemporary populations from Sardinia and from continental Italy provided an indication of the trend leading to the final part of the peopling history. Finally, Upper Palaeolithic and Mesolithic samples were included in the analyses to compare the prehistoric Sardinians with some of their potential continental ancestors. The analysis is based on multivariate techniques including Mahalanobis D(2) distance, non-parametric multidimensional scaling (MDS) and principal component analysis (PCA). The results showed the tendency to progressive differentiation between Sardinian groups and peninsular Italian groups, with the possible exception of a discontinuity showed by the Bonnànaro (Early Bronze Age) Sardinian sample. Several aspects of the morphological results were found to agree with the current genetic evidence available for the present-day Sardinian population and a Nuragic sample: (1) biological divergence between the Sardinian and peninsular Italian populations; (2) similarity/continuity among Neolithic, Bronze Age and recent Sardinians; (3) biological separation between the Nuragic and Etruscan populations; (4) contribution of a Palaeo-Mesolithic gene pool to the genetic structure of current Sardinians.

Analysis of 12,517 inhabitants of a Sardinian geographic isolate reveals that predispositions to thrombocytopenia and thrombocytosis are inherited traits

BACKGROUND

Thrombocytopenia is a common finding in several diseases but almost nothing is known about the prevalence of thrombocytopenia in the general population. We examined the prevalence of thrombocytopenia and determinants of platelet count in a healthy population with a wide age range.

DESIGN AND METHODS

We performed a cross-sectional study on 12,517 inhabitants of ten villages (80% of residents) in a secluded area of Sardinia (Ogliastra). Participants underwent a complete blood count evaluation and a structured questionnaire, used to collect epidemiological data.

RESULTS

We observed a platelet count lower than 150 × 10⁹/L in 3.2% (2.8%-3.6%) of females and 4.8% (4.3%-5.4%) of males, with a value of 3.9% (3.6%-4.3%) in the entire population. Thrombocytopenia was mild (platelet count: 100 × 10⁹/L-150 × 10⁹/L), asymptomatic and not associated with other cytopenias or overt disorders in most cases. Its standardized prevalence was quite different in different villages, with values ranging from 1.5% to 6.8%, and was negatively correlated with the prevalence of a mild form of thrombocytosis, which ranged from 0.9% to 4.5%. Analysis of platelet counts across classes of age revealed that platelet number decreased progressively with aging. As a consequence, thrombocytopenia was nearly absent in young people and its prevalence increased regularly during lifetime. The opposite occurred for thrombocytosis.

CONCLUSIONS

Given the high genetic differentiation among Ogliastra villages with "high" and "low" platelet counts and the substantial heritability of this quantitative trait (54%), we concluded that the propensity to present mild and transient thrombocytosis in youth and to acquire mild thrombocytopenia during aging are new genetic traits.

Near eastern neolithic genetic input in a small oasis of the Egyptian Western Desert

The Egyptian Western Desert lies on an important geographic intersection between Africa and Asia. Genetic diversity of this region has been shaped, in part, by climatic changes in the Late Pleistocene and Holocene epochs marked by oscillating humid and arid periods. We present here a whole genome analysis of mitochondrial DNA (mtDNA) and high-resolution molecular analysis of nonrecombining Y-chromosomal (NRY) gene pools of a demographically small but autochthonous population from the Egyptian Western Desert oasis el-Hayez. Notwithstanding signs of expected genetic drift, we still found clear genetic evidence of a strong Near Eastern input that can be dated into the Neolithic. This is revealed by high frequencies and high internal variability of several mtDNA lineages from haplogroup T. The whole genome sequencing strategy and molecular dating allowed us to detect the accumulation of local mtDNA diversity to 5,138 +/- 3,633 YBP. Similarly, theY-chromosome gene pool reveals high frequencies of the Near Eastern J1 and the North African E1b1b1b lineages, both generally known to have expanded within North Africa during the Neolithic. These results provide another piece of evidence of the relatively young population history of North Africa.

Mitochondrial haplogroup U5b3: a distant echo of the epipaleolithic in Italy and the legacy of the early Sardinians

There are extensive data indicating that some glacial refuge zones of southern Europe (Franco-Cantabria, Balkans, and Ukraine) were major genetic sources for the human recolonization of the continent at the beginning of the Holocene. Intriguingly, there is no genetic evidence that the refuge area located in the Italian Peninsula contributed to this process. Here we show, through phylogeographic analyses of mitochondrial DNA (mtDNA) variation performed at the highest level of molecular resolution (52 entire mitochondrial genomes), that the most likely homeland for U5b3-a haplogroup present at a very low frequency across Europe-was the Italian Peninsula. In contrast to mtDNA haplogroups that expanded from other refugia, the Holocene expansion of haplogroup U5b3 toward the North was restricted by the Alps and occurred only along the Mediterranean coasts, mainly toward nearby Provence (southern France). From there, approximately 7,000-9,000 years ago, a subclade of this haplogroup moved to Sardinia, possibly as a result of the obsidian trade that linked the two regions, leaving a distinctive signature in the modern people of the island. This scenario strikingly matches the age, distribution, and postulated geographic source of a Sardinian Y chromosome haplogroup (I2a2-M26), a paradigmatic case in the European context of a founder event marking both female and male lineages.

Correcting for purifying selection: an improved human mitochondrial molecular clock

There is currently no calibration available for the whole human mtDNA genome, incorporating both coding and control regions. Furthermore, as several authors have pointed out recently, linear molecular clocks that incorporate selectable characters are in any case problematic. We here confirm a modest effect of purifying selection on the mtDNA coding region and propose an improved molecular clock for dating human mtDNA, based on a worldwide phylogeny of > 2000 complete mtDNA genomes and calibrating against recent evidence for the divergence time of humans and chimpanzees. We focus on a time-dependent mutation rate based on the entire mtDNA genome and supported by a neutral clock based on synonymous mutations alone. We show that the corrected rate is further corroborated by archaeological dating for the settlement of the Canary Islands and Remote Oceania and also, given certain phylogeographic assumptions, by the timing of the first modern human settlement of Europe and resettlement after the Last Glacial Maximum. The corrected rate yields an age of modern human expansion in the Americas at approximately 15 kya that-unlike the uncorrected clock-matches the archaeological evidence, but continues to indicate an out-of-Africa dispersal at around 55-70 kya, 5-20 ky before any clear archaeological record, suggesting the need for archaeological research efforts focusing on this time window. We also present improved rates for the mtDNA control region, and the first comprehensive estimates of positional mutation rates for human mtDNA, which are essential for defining mutation models in phylogenetic analyses.

Exaggerated status of "novel" and "pathogenic" mtDNA sequence variants due to inadequate database searches

Given its relative ease, screening the entire mitochondrial DNA (mtDNA) for heteroplasmic or novel homoplasmic mutations has become part of the routine diagnostic workup for the molecular geneticist confronted with a disease case exhibiting clinical and biochemical features of mitochondrial dysfunction. "Novelty" of a given mtDNA variant is most often equated with nonregistration in the extensive MITOMAP database (www.mitomap.org). This practice has led to a number of spurious findings and wrong conclusions concerning the pathogenic status of specific mtDNA mutations, especially in the absence of proper evaluation and pathogenicity scoring. We demonstrate by way of real cases targeting the mt-tRNA(Cys) (MT-TC) gene and a stretch within the MT-ND3 gene, that a straightforward Google search can identify twice as many previously observed mutations than any MITOMAP query could achieve. Further, we reassess the recent rediscovery of m.15287T>C by listing all known occurrences and, where possible, providing the haplogroup context, shedding new light on the potential pathogenicity status of m.15287T>C.

High differentiation among eight villages in a secluded area of Sardinia revealed by genome-wide high density SNPs analysis

To better design association studies for complex traits in isolated populations it's important to understand how history and isolation moulded the genetic features of different communities. Population isolates should not “a priori” be considered homogeneous, even if the communities are not distant and part of a small region. We studied a particular area of Sardinia called Ogliastra, characterized by the presence of several distinct villages that display different history, immigration events and population size. Cultural and geographic isolation characterized the history of these communities. We determined LD parameters in 8 villages and defined population structure through high density SNPs (about 360 K) on 360 unrelated people (45 selected samples from each village). These isolates showed differences in LD values and LD map length. Five of these villages show high LD values probably due to their reduced population size and extreme isolation. High genetic differentiation among villages was detected. Moreover population structure analysis revealed a high correlation between genetic and geographic distances. Our study indicates that history, geography and biodemography have influenced the genetic features of Ogliastra communities producing differences in LD and population structure. All these data demonstrate that we can consider each village an isolate with specific characteristics. We suggest that, in order to optimize the study design of complex traits, a thorough characterization of genetic features is useful to identify the presence of sub-populations and stratification within genetic isolates.

Mitochondrial DNA haplogroup H structure in North Africa

BACKGROUND

The Strait of Gibraltar separating the Iberian Peninsula from North Africa is thought to be a stronger barrier to gene flow for male than for female lineages. However, the recent subdivision of the haplogroup H at mitochondrial DNA (mtDNA) level has revealed greater genetic differentiation among geographic regions than previously detected. The dissection of the mtDNA haplogroup H in North Africa, and its comparison with the Iberian Peninsula and Near-East profiles would help clarify the relative affinities among these regions.

RESULTS

Like the Iberian Peninsula, the dominant mtDNA haplogroup H subgroups in North Africa are H1 (42%) and H3 (13%). The similarity between these regions is stronger in the North-West edge affecting mainly Moroccan Arabs, West Saharans and Mauritanians, and decreases eastwards probably due to gene flow from Near East as attested for the higher frequencies of H4, H5, H7, H8 and H11 subgroups. Moroccan Berbers show stronger affinities with Tunisian and Tunisian Berbers than with Moroccan Arabs. Coalescence ages for H1 (11 +/- 2 ky) and H3 (11 +/- 4 ky) in North Africa point to the possibility of a late Palaeolithic settlement for these lineages similar to those found for other mtDNA haplogroups. Total and partial mtDNA genomic sequencing unveiled stronger mtDNA differentiation among regions than previously found using HVSI mtDNA based analysis.

CONCLUSION

The subdivision of the mtDNA haplogroup H in North Africa has confirmed that the genetic differentiation found among Western and Eastern populations is mainly due to geographical rather than cultural barriers. It also shows that the historical Arabian role on the region had more a cultural than a demic effect. Whole mtDNA sequencing of identical H haplotypes based on HVSI and RFLP information has unveiled additional mtDNA differences between North African and Iberian Peninsula lineages, pointing to an older mtDNA genetic flow between regions than previously thought. Based on this new information, it seems that the Strait of Gibraltar barrier affected both male and female gene flow in a similar fashion.

Data from complete mtDNA sequencing of Tunisian centenarians: testing haplogroup association and the "golden mean" to longevity

Since the mitochondrial theory of ageing was proposed, mitochondrial DNA (mtDNA) diversity has been largely studied in old people, however complete genomes are still rare, being limited to Japanese and UK/US samples. In this work, we evaluated possible longevity associated polymorphisms/haplogroups in an African population, from Tunisia, by performing complete mtDNA sequencing. This population has a mixed Eurasian/sub-Saharan mtDNA gene pool, which could potentially facilitate the evaluation of association for sub-Saharan lineages. Sub-Saharan haplogroups were shown to be significantly less represented in centenarians (9.5%) than in controls (54.5%), but it is not possible to rule out an influence of population structure, which is high in these populations. No recurrent polymorphism were more frequent in centenarians than in controls, and although the Tunisian centenarians presented less synonymous and replacement polymorphisms than controls, this difference was not statistically significant. So far, it does not seem that centenarians have significantly less mildly deleterious substitutions, not only in Tunisia but also in Japanese and UK/US samples, as tested here, not favouring a "golden mean" to longevity.

Complete mitochondrial genome sequence of the Tyrolean Iceman

The Tyrolean Iceman was a witness to the Neolithic-Copper Age transition in Central Europe 5350-5100 years ago, and his mummified corpse was recovered from an Alpine glacier on the Austro-Italian border in 1991 [1]. Using a mixed sequencing procedure based on PCR amplification and 454 sequencing of pooled amplification products, we have retrieved the first complete mitochondrial-genome sequence of a prehistoric European. We have then compared it with 115 related extant lineages from mitochondrial haplogroup K. We found that the Iceman belonged to a branch of mitochondrial haplogroup K1 that has not yet been identified in modern European populations. This is the oldest complete Homo sapiens mtDNA genome generated to date. The results point to the potential significance of complete-ancient-mtDNA studies in addressing questions concerning the genetic history of human populations that the phylogeography of modern lineages is unable to tackle.

Direct maximum parsimony phylogeny reconstruction from genotype data

Background

Maximum parsimony phylogenetic tree reconstruction from genetic variation data is a fundamental problem in computational genetics with many practical applications in population genetics, whole genome analysis, and the search for genetic predictors of disease. Efficient methods are available for reconstruction of maximum parsimony trees from haplotype data, but such data are difficult to determine directly for autosomal DNA. Data more commonly is available in the form of genotypes, which consist of conflated combinations of pairs of haplotypes from homologous chromosomes. Currently, there are no general algorithms for the direct reconstruction of maximum parsimony phylogenies from genotype data. Hence phylogenetic applications for autosomal data must therefore rely on other methods for first computationally inferring haplotypes from genotypes.

Results

In this work, we develop the first practical method for computing maximum parsimony phylogenies directly from genotype data. We show that the standard practice of first inferring haplotypes from genotypes and then reconstructing a phylogeny on the haplotypes often substantially overestimates phylogeny size. As an immediate application, our method can be used to determine the minimum number of mutations required to explain a given set of observed genotypes.

Conclusion

Phylogeny reconstruction directly from unphased data is computationally feasible for moderate-sized problem instances and can lead to substantially more accurate tree size inferences than the standard practice of treating phasing and phylogeny construction as two separate analysis stages. The difference between the approaches is particularly important for downstream applications that require a lower-bound on the number of mutations that the genetic region has undergone.

What is a 'novel' mtDNA mutation--and does 'novelty' really matter?

The hunt for pathogenic mitochondrial DNA (mtDNA) mutations is often fueled by the seeming novelty of mutations that are either nonsynonymous or affect the protein synthesis machinery in patients. In order to determine the novelty of a detected mutation, the working geneticist nearly always consults MITOMAP--often exclusively. By reanalyzing some case studies of refractory anemia with ring sideroblasts, prostate cancer, and hearing impairment, we demonstrate that the practice of solely relying on MITOMAP can be most misleading. A notorious example is the T1243C mutation, which was assessed to be novel and deemed to be associated with some (rare) disease simply because researchers did not realize that T1243C defines a deep branch in the Eurasian mtDNA phylogeny. The majority of 'novel' mutations suspected of being pathogenic are in actual fact known (and presumably neutral) polymorphisms (although unknown to MITOMAP), and this becomes glaringly evident when proper database searches and straightforward Internet queries are carried out.