Characterization of human control region sequences of the African American SWGDAM forensic mtDNA data set

Complete genome sequence and bioinformatics analysis of nine Egyptian females with clinical information from different geographic regions in Egypt

Egyptians are at a crossroad between Africa and Eurasia, providing useful genomic resources for analyzing both genetic and environmental factors for future personalized medicine. Two personal Egyptian whole genomes have been published previously by us and here nine female whole genome sequences with clinical information have been added to expand the genomic resource of Egyptian personal genomes. Here we report the analysis of whole genomes of nine Egyptian females from different regions using Illumina short-read sequencers. At 30x sequencing coverage, we identified 12 SNPs that were shared in most of the subjects associated with obesity which are concordant with their clinical diagnosis. Also, we found mtDNA mutation A4282G is common in all the samples and this is associated with chronic progressive external ophthalmoplegia (CPEO). Haplogroup and Admixture analyses revealed that most Egyptian samples are close to the other north Mediterranean, Middle Eastern, and European, respectively, possibly reflecting the into-Africa influx of human migration. In conclusion, we present whole-genome sequences of nine Egyptian females with personal clinical information that cover the diverse regions of Egypt. Although limited in sample size, the whole genomes data provides possible geno-phenotype candidate markers that are relevant to the region's diseases.

Mitochondrial DNA in human identification: a review

Mitochondrial DNA (mtDNA) presents several characteristics useful for forensic studies, especially related to the lack of recombination, to a high copy number, and to matrilineal inheritance. mtDNA typing based on sequences of the control region or full genomic sequences analysis is used to analyze a variety of forensic samples such as old bones, teeth and hair, as well as other biological samples where the DNA content is low. Evaluation and reporting of the results requires careful consideration of biological issues as well as other issues such as nomenclature and reference population databases. In this work we review mitochondrial DNA profiling methods used for human identification and present their use in the main cases of humanidentification focusing on the most relevant issues for forensics.

The association of mitochondrial DNA haplogroups with POAG in African Americans

Mitochondrial dysfunction has been implicated in the pathogenesis of primary open-angle glaucoma (POAG). However, the potential significance of mitochondrial DNA (mtDNA) haplogroups to POAG has not been evaluated in the overaffected African American population. To investigate the association of mtDNA haplogroups with POAG and its phenotypic characteristics, genotyping data from 4081 African American subjects (1919 cases and 2162 controls) was analyzed using 1293 positions on mtDNA. The overall frequency of mtDNA haplogroups in the Primary Open-Angle African American Glaucoma Genetics (POAAGG) study cohort was 37% L3, 29% L2, 21% L1, 4% L0, and 10% non-African haplogroups (non-L). When all haplogroups (L0, L1, L2, and non-L) were compared against theL3 reference group, after adjusting by age and principal component of ancestry, the non-L3 haplogroups showed higher risk of POAG (OR-1.19, p = 0.02), with a particularly strong association among males (OR = 1.41, p = 0.003). More specifically the non-L group was associated with higher POAG risk than the L3 haplogroup (OR = 1.77, p = 0.007, Bonferroni adjusted p = 0.027) and to the L3e (n = 256, OR = 1.92, p = 0.007, Bonferroni adjusted p = 0.029). No significant association was found when genders were analyzed together or in female only analysis. There were no significant differences in various POAG endophenotypes across mtDNA haplogroups. This study expands our knowledge of mitochondrial genetics and mtDNA haplogroup associations in African American POAG. Further work is needed to better understand the functional role of mtDNA polymorphisms and their interactions with nuclear genes that affect POAG.

The mitochondrial landscape of African Americans: An examination of more than 2500 control region haplotypes from 22 U.S. locations

The mitochondrial DNA (mtDNA) control region (16024-576) was Sanger-sequenced for a total of 2563 self-identified African Americans, using automated processing techniques and data review standards exceeding guidelines for forensic applications. Genetic diversity ranged from 0.9952 to 0.9998 in 22 population samples from 20 different states. Haplogroups of African ancestry, found in 82.48% of individuals overall, were most concentrated in the Southeast U.S. and decreased to the north and west. West African and West Central African haplotypes were well-represented in the population samples, especially in the southern U.S. states, while East African haplogroups were observed in low-frequency clusters in a handful of locations across the country. East Asian, Native American, and West Eurasian admixture was present in 3.16%, 2.93%, and 11.43% of samples, respectively. While some geographic substructure was detected across the population samples as clines in admixture frequencies, 20 of the 22 population samples were found to be statistically indistinguishable by pairwise comparisons and AMOVA calculations. Datasets from Hawaii and Idaho, however, were clear outliers. Overall, these more than 2500 control region sequences represent the most comprehensive regional sampling of African American mtDNA diversity to date, and are suitable for use in a forensic mtDNA database. The population data are made available via EMPOP (www.empop.org) and GenBank.

Kuwaiti population subgroup of nomadic Bedouin ancestry-Whole genome sequence and analysis

Kuwaiti native population comprises three distinct genetic subgroups of Persian, "city-dwelling" Saudi Arabian tribe, and nomadic "tent-dwelling" Bedouin ancestry. Bedouin subgroup is characterized by presence of 17% African ancestry; it owes it origin to nomadic tribes of the deserts of Arabian Peninsula and North Africa. By sequencing whole genome of a Kuwaiti male from this subgroup at 41X coverage, we report 3,752,878 SNPs, 411,839 indels, and 8451 structural variations. Neighbor-joining tree, based on shared variant positions carrying disease-risk alleles between the Bedouin and other continental genomes, places Bedouin genome at the nexus of African, Asian, and European genomes in concordance with geographical location of Kuwait and Peninsula. In congruence with participant's medical history for morbid obesity and bronchial asthma, risk alleles are seen at deleterious SNPs associated with obesity and asthma. Many of the observed deleterious 'novel' variants lie in genes associated with autosomal recessive disorders characteristic of the region.

Characterization of mitochondrial haplogroups in a large population-based sample from the United States

Mitochondrial DNA (mtDNA) haplogroups are valuable for investigations in forensic science, molecular anthropology, and human genetics. In this study, we developed a custom panel of 61 mtDNA markers for high-throughput classification of European, African, and Native American/Asian mitochondrial haplogroup lineages. Using these mtDNA markers, we constructed a mitochondrial haplogroup classification tree and classified 18,832 participants from the National Health and Nutrition Examination Surveys (NHANES). To our knowledge, this is the largest study to date characterizing mitochondrial haplogroups in a population-based sample from the United States, and the first study characterizing mitochondrial haplogroup distributions in self-identified Mexican Americans separately from Hispanic Americans of other descent. We observed clear differences in the distribution of maternal genetic ancestry consistent with proposed admixture models for these subpopulations, underscoring the genetic heterogeneity of the United States Hispanic population. The mitochondrial haplogroup distributions in the other self-identified racial/ethnic groups within NHANES were largely comparable to previous studies. Mitochondrial haplogroup classification was highly concordant with self-identified race/ethnicity (SIRE) in non-Hispanic whites (94.8 %), but was considerably lower in admixed populations including non-Hispanic blacks (88.3 %), Mexican Americans (81.8 %), and other Hispanics (61.6 %), suggesting SIRE does not accurately reflect maternal genetic ancestry, particularly in populations with greater proportions of admixture. Thus, it is important to consider inconsistencies between SIRE and genetic ancestry when performing genetic association studies. The mitochondrial haplogroup data that we have generated, coupled with the epidemiologic variables in NHANES, is a valuable resource for future studies investigating the contribution of mtDNA variation to human health and disease.

Recurrent tissue-specific mtDNA mutations are common in humans

Mitochondrial DNA (mtDNA) variation can affect phenotypic variation; therefore, knowing its distribution within and among individuals is of importance to understanding many human diseases. Intra-individual mtDNA variation (heteroplasmy) has been generally assumed to be random. We used massively parallel sequencing to assess heteroplasmy across ten tissues and demonstrate that in unrelated individuals there are tissue-specific, recurrent mutations. Certain tissues, notably kidney, liver and skeletal muscle, displayed the identical recurrent mutations that were undetectable in other tissues in the same individuals. Using RFLP analyses we validated one of the tissue-specific mutations in the two sequenced individuals and replicated the patterns in two additional individuals. These recurrent mutations all occur within or in very close proximity to sites that regulate mtDNA replication, strongly implying that these variations alter the replication dynamics of the mutated mtDNA genome. These recurrent variants are all independent of each other and do not occur in the mtDNA coding regions. The most parsimonious explanation of the data is that these frequently repeated mutations experience tissue-specific positive selection, probably through replication advantage.

DNA mutations are expected to be formed randomly, thus any reproducible pattern of DNA somatic mutations across multiple individuals or even across organs within each individual is highly unexpected. Using next generation sequencing of multiple tissues from the same individuals we found several somatic mutations in mitochondrial DNA that appear in a heteroplasmic state in all individuals examined, but only in particular tissues. These mutations were only found in known regions of replication control for the mitochondrial DNA. These data imply the presence of tissue-specific positive selection for these variants.

Mitochondrial DNA diversity in the African American population

Genetic polymorphism along mitochondrial DNA (mtDNA) defines population-specific signatures called mtDNA haplogroups. Estimation of mtDNA haplogroup distribution may be prone to errors, notably if the study sample is not drawn from a multicenter cohort. Here, we report on mtDNA diversity in a sample of African American individuals (n = 343) enrolled in a multicenter cohort. Sequencing of the hypervariable regions I and II of the D-loop control region showed that the most common mitochondrial variants are 73G, 146C, 150T, 152C, 189G, 16278T, and 16311C. In agreement with the published data, we observed 17 common mtDNA haplogroups: L0, L1, L1b, L1c, L2, L2a, L2b, L2c, L2e, L3, L3b, L3d, L3e, L3f, L3h, L3x, and L4. The most commonly observed haplogroup is L2a (19.8%), followed by L1b (10.2%). Overall, the observed mtDNA haplogroup distribution in our study is similar to those published for the African American and the African populations.

Evaluating self-declared ancestry of U.S. Americans with autosomal, Y-chromosomal and mitochondrial DNA

The current U.S. population represents an amalgam of individuals originating mainly from four continental regions (Africa, Europe, Asia and America). To study the genetic ancestry and compare with self-declared ancestry we have analyzed paternally, maternally and bi-parentally inherited DNA markers sensitive for indicating continental genetic ancestry in all four major U.S. American groups. We found that self-declared U.S. Hispanics and U.S. African Americans tend to show variable degrees of continental genetic admixture among the three genetic systems, with evidence for a marked sex-biased admixture history. Moreover, for these two groups we observed significant regional variation across the country in genetic admixture. In contrast, self-declared U.S. European and U.S. Asian Americans were genetically more homogeneous at the continental ancestry level. Two autosomal ancestry-sensitive markers located in skin pigmentation candidate genes showed significant differences in self-declared U.S. African Americans or U.S. European Americans, relative to their assumed parental populations from Africa or Europe. This provides genetic support for the importance of skin color in the complex process of ancestry identification.

The imprint of the Slave Trade in an African American population: mitochondrial DNA, Y chromosome and HTLV-1 analysis in the Noir Marron of French Guiana

Background

Retracing the genetic histories of the descendant populations of the Slave Trade (16^th-19^th centuries) is particularly challenging due to the diversity of African ethnic groups involved and the different hybridisation processes with Europeans and Amerindians, which have blurred their original genetic inheritances. The Noir Marron in French Guiana are the direct descendants of maroons who escaped from Dutch plantations in the current day Surinam. They represent an original ethnic group with a highly blended culture. Uniparental markers (mtDNA and NRY) coupled with HTLV-1 sequences (env and LTR) were studied to establish the genetic relationships linking them to African American and African populations.

Results

All genetic systems presented a high conservation of the African gene pool (African ancestry: mtDNA = 99.3%; NRY = 97.6%; HTLV-1 env = 20/23; HTLV-1 LTR = 6/8). Neither founder effect nor genetic drift was detected and the genetic diversity is within a range commonly observed in Africa. Higher genetic similarities were observed with the populations inhabiting the Bight of Benin (from Ivory Coast to Benin). Other ancestries were identified but they presented an interesting sex-bias. Whilst male origins spread throughout the north of the bight (from Benin to Senegal), female origins were spread throughout the south (from the Ivory Coast to Angola).

Conclusions

The Noir Marron are unique in having conserved their African genetic ancestry, despite major cultural exchanges with Amerindians and Europeans through inhabiting the same region for four centuries. Their maroon identity and the important number of slaves deported in this region have maintained the original African diversity. All these characteristics permit to identify a major origin located in the former region of the Gold Coast and the Bight of Benin; regions highly impacted by slavery, from which goes a sex-biased longitudinal gradient of ancestry.

African mitochondrial DNA subhaplogroups and peripheral neuropathy during antiretroviral therapy

Susceptibility to peripheral neuropathy during antiretroviral therapy with nucleoside reverse-transcriptase inhibitors was previously associated with a European mitochondrial DNA (mtDNA) haplogroup among non-Hispanic white persons. To determine whether nucleoside reverse-transcriptase inhibitor-associated peripheral neuropathy was related to mtDNA variation in non-Hispanic black persons, we sequenced mtDNA of participants from AIDS Clinical Trials Group study 384. Of 156 non-Hispanic black persons with genomic data, 51 (33%) developed peripheral neuropathy. In a multivariate model, African mtDNA subhaplogroup L1c was an independent predictor of peripheral neuropathy (odds ratio, 3.7 [95% confidence interval, 1.1-12.0]). An African mtDNA subhaplogroup is for the first time implicated in susceptibility to nucleoside reverse-transcriptase inhibitor-associated toxicity.

Mitochondrial control region sequences from an African American population sample

Entire mitochondrial control region data were generated for 248 African American individuals, which had been previously typed for 15 autosomal STRs [J.M. Butler, R. Schoske, P.M. Vallone, J.W. Redman, M.C. Kline, Allele frequencies for 15 autosomal STR loci on U.S. Caucasian, African American, and Hispanic populations, J. Forensic Sci. 48 (2003) 908-911].

Evaluation of group genetic ancestry of populations from Philadelphia and Dakar in the context of sex-biased admixture in the Americas

BACKGROUND

Population history can be reflected in group genetic ancestry, where genomic variation captured by the mitochondrial DNA (mtDNA) and non-recombining portion of the Y chromosome (NRY) can separate female- and male-specific admixture processes. Genetic ancestry may influence genetic association studies due to differences in individual admixture within recently admixed populations like African Americans.

PRINCIPAL FINDINGS

We evaluated the genetic ancestry of Senegalese as well as European Americans and African Americans from Philadelphia. Senegalese mtDNA consisted of approximately 12% U haplotypes (U6 and U5b1b haplotypes, common in North Africa) while the NRY haplotypes belonged solely to haplogroup E. In Philadelphia, we observed varying degrees of admixture. While African Americans have 9-10% mtDNAs and approximately 31% NRYs of European origin, these results are not mirrored in the mtDNA/NRY pools of European Americans: they have less than 7% mtDNAs and less than 2% NRYs from non-European sources. Additionally, there is <2% Native American contribution to Philadelphian African American ancestry and the admixture from combined mtDNA/NRY estimates is consistent with the admixture derived from autosomal genetic data. To further dissect these estimates, we have analyzed our samples in the context of different demographic groups in the Americas.

CONCLUSIONS

We found that sex-biased admixture in African-derived populations is present throughout the Americas, with continual influence of European males, while Native American females contribute mainly to populations of the Caribbean and South America. The high non-European female contribution to the pool of European-derived populations is consistently characteristic of Iberian colonization. These data suggest that genomic data correlate well with historical records of colonization in the Americas.

The roles of ethnicity and antiretrovirals in HIV-associated polyneuropathy: a pilot study

BACKGROUND

In the pre-highly active antiretroviral therapy (HAART) era, distal sensory polyneuropathy (DSP) was associated with markers of advanced HIV infection and the use of neurotoxic antiretrovirals (ARVs). As HAART became widespread, and the AIDS epidemic shifted into minority populations, the risk factors for DSP became less clear. We explore the roles of ethnicity and ARV in the development of DSP in an HAART era cohort.

METHODS

Data from 336 HIV-positive adults were obtained from the Manhattan HIV Brain Bank. One hundred four participants had no DSP at entry visit; at least 1 follow-up visit; and a self-identified ethnicity of non-Hispanic white, Hispanic, or African American.

RESULTS

Fifty percent of participants developed DSP; of those, 67% were symptomatic. Participants who developed DSP were older (P = 0.02) and had higher CD4 counts (P = 0.001). ARV-DSP was more common in Hispanics (P = 0.02) and intravenous drug users (P = 0.02). There was a trend for higher pain scores in Hispanics with symptomatic DSP (P = 0.08).

CONCLUSIONS

This study suggests that there are ethnic disparities in the clinical manifestations of HIV-related neuropathies including pain and the susceptibility to ARV-DSP. Further studies of larger cohorts are indicated to explore the etiology of these differences.

Two common mitochondrial DNA polymorphisms are highly associated with migraine headache and cyclic vomiting syndrome

Mitochondrial dysfunction is a hypothesized component in the multifactorial pathogenesis of migraine without aura (MoA, 'common migraine') and the related condition of cyclic vomiting syndrome (CVS). In this study, the entire mitochondrial genome was sequenced in 20 haplogroup-H CVS patients, a subject group studied because of greater genotypic and phenotypic homogeneity. Sequences were compared against haplogroup-H controls. Polymorphisms of interest were tested in 10 additional CVS subjects and in 112 haplogroup-H adults with MoA. The 16519C-->T polymorphism was found to be highly disease associated: 21/30 CVS subjects [70%, odds ratio (OR) 6.2] and 58/112 migraineurs (52%, OR 3.6) vs. 63/231 controls (27%). A second polymorphism, 3010G-->A, was found to be highly disease associated in those subjects with 16519T: 6/21 CVS subjects (29%, OR 17) and 15/58 migraineurs (26%, OR 15) vs. 1/63 controls (1.6%). Our data suggest that these polymorphisms constitute a substantial proportion of the genetic factor in migraine pathogenesis, and strengthen the hypothesis that there is a component of mitochondrial dysfunction in migraine.

Mitochondrial DNA control region polymorphism in the population of Alagoas state, north-eastern Brazil

The sequences of the two hypervariable (HV) segments of the mitochondrial DNA (mtDNA) control region were determined in 167 randomly selected, unrelated individuals living in the state of Alagoas, north-eastern Brazil. One hundred and forty-five different haplotypes, associated with 139 variable positions, were determined. More than 95% of the mtDNA sequences could be allocated to specific mtDNA haplogroups according to the mutational motifs. Length heteroplasmy in the C-stretch HV1 and HV2 regions was observed in 22 and 11%, respectively, of the population sample. The genetic diversity was estimated to be 0.9975 and the probability of two random individuals presenting identical mtDNA haplotypes was 0.0084. The most frequent haplotype was shared by six individuals. All sequences showed high-quality values and phantom mutations were not detected. The diversity revealed in the mitochondrial control region indicates the importance of this locus for forensic casework and population studies within Alagoas, Brazil.

Characterization of human control region sequences for Spanish individuals in a forensic mtDNA data set

Population data on the hypervariable regions of the mitochondrial DNA (mtDNA) genome are used to convey the relative rarity of mtDNA profiles obtained from evidence samples and of profiles used to identify missing persons. In this study, mtDNA profiles of Spanish individuals (n=312) were analyzed to describe haplogroup distributions and to determine relevant single nucleotide polymorphisms (SNPs) of those haplogroups. All nine common European haplogroups were observed in the sample, and these were divided into subgroups when possible. Haplogroup H was the most common haplogroup. The haplogroups U, J, T, and V were the next most frequent groups, each occurring at a frequency of 6.4% or greater. In addition, African and Asian sequences were present though rare in the samples. The data were compared with and found to be similar to other published data sets. There were 109 SNPs observed in the data set, including 10 positions not previously reported. The most variable sites are consistent with other studies.

Application of a quasi-median network analysis for the visualization of character conflicts to a population sample of mitochondrial DNA control region sequences from southern Germany (Ulm)

Entire mtDNA control region sequences from 100 individuals in a west Eurasian population sample from southern Germany (around the city of Ulm) were generated and analyzed. The control region was amplified in one piece and sequenced with ten different sequencing primers. Sequence evaluation was performed independently. Phylogenetic analyses were used for quality assurance purposes and for the determination of the haplogroup affiliation of the samples. The sequences were scrutinized performing a quasi-median network analysis. To visualize character conflicts, frequent mutations were filtered, and the reduced data were represented by the torso of their quasi-median network. Character incompatibilities were found to be based on real biological patterns of homoplasy. The population data will be incorporated in the EMPOP database ( http://www.empop.org ).

Phylogeographic investigations: the role of trees in forensic genetics

The human mitochondrial DNA (mtDNA) genome is commonly analyzed in various disciplines, such as population, medical, and forensic genetics, but conceptual and scientific exchange between them is still limited. Here we review several aspects of the mtDNA phylogeny that are particularly--but not exclusively--of interest to the forensic community. Among the issues that arise, we emphasize the importance of integrating evolutionary concepts into the forensic routine. We also discuss topics such as mtDNA mutation-rate heterogeneity and the weight of evidence, ethnic affiliations of mtDNA profiles, and the abuse of reference databases. Finally, we show the usefulness of coding-region variation in a forensic context.

Quality assessment of DNA sequence data: autopsy of a mis-sequenced mtDNA population sample

Published DNA data sets constitute a body of sequencing results resting in silico that are supposed to reflect the variation of (once) living cells. In cases where the DNA variation reported is suspected to be fraught with artefacts, an autopsy of the full body of data is needed to clarify the amount and causes of mis-sequencing. In this paper we elaborate on strategies that allow a clear-cut identification of the problems in severely flawed mtDNA data. This approach is applied, by way of example, to a data set of HVS-I sequences from the Caucasus, published by Nasidze & Stoneking in 2001. These data bear numerous ambiguous nucleotide positions and suffer from an even higher number of phantom mutations, indicating that severe biochemical problems adversely influenced those sequencing results at the time. Furthermore, systematic omission of sequences with a long C-stretch (incurred by a transition at position 16189) must have severely biased the data set. Since no complete correction of these data has appeared to date, this example of mis-sequencing necessitates circumstantial evidence that is bullet-proof.

Evaluation of Variation in Control Region Sequences for Hispanic Individuals in the SWGDAM mtDNA Data Set

The Scientific Working Group on DNA Analysis Methods (SWGDAM) Hispanic data set was analyzed to determine the diversity, phylogeny, and relevant single nucleotide polymorphisms (SNPs) that describe haplogroup patterns for Hispanic Americans (N=686), and to assess the degree of admixture regarding mitochondrial DNA (mtDNA). The largest component of admixture based on mtDNA analysis derives from the four major haplogroups previously observed in Native American ancestry, including A (29.3%), B (15.7%), C (20.6%), and D (4.8%). European (17.8%) and African (11.8%) haplogroups also were observed within this data set. Hispanic SWGDAM samples from the southwest, compared with other SWGDAM Hispanic samples, were observed to have a greater percent of Native American haplogroups present (79.9%), and fewer African American haplogroups (4.5%). A total of 234 SNPs were observed in the data set, including 36 newly reported variable positions. These SWGDAM Hispanic data set SNPs ranged from having 1 to 31 changes (Length=L) on the phylogenetic tree, with site 16519 being the most variable. On average, there were 3.9 character changes for each variable position on the tree. The most variable sites (with 13 or more changes each listed from fastest to slowest) observed were 16519 (L=31), 16189 (L=23), 152 (L=23), 16311 (L=19), 146 (L=17), 195 (L=17), 16093 (L=15), 16362 (L=14), 16129 (L=13), 150 (L=13), and 153 (L=13). These sites are consistent with other reports on highly variable positions. A total of 27 SNPs were chosen to identify all clusters containing 1% (N=7) or more individuals in the SWGDAM Hispanic data set. The descriptive analyses revealed that the SWGDAM Hispanic data set is similar to published Native American and Hispanic data sets.

A practical guide to mitochondrial DNA error prevention in clinical, forensic, and population genetics

Several suggestions have been made for avoiding errors in mitochondrial DNA (mtDNA) sequencing and documentation. Unfortunately, the current clinical, forensic, and population genetic literature on mtDNA still delivers a large number of studies with flawed sequence data, which, in extreme cases, damage the whole message of a study. The phylogenetic approach has been shown to be useful for pinpointing most of the errors. However, many geneticists, especially in the forensic and medical fields, are not familiar with either effective search strategies or the evolutionary terminology. We here provide a manual that should help prevent errors at any stage by re-examining data fresh from the sequencer in the light of previously published data. A fictitious case study of a European mtDNA data set (albeit composed from the literature) then demonstrates the steps one has to go through in order to assess the quality of sequencing and documentation.