Discordant patterns of mtDNA and ethno-linguistic variation in 14 Iranian Ethnic groups

A comprehensive review of HVS-I mitochondrial DNA variation of 19 Iranian populations

Iran is located along the Central Asian corridor, a natural artery that has served as a cross-continental route since the first anatomically modern human populations migrated out of Africa. We compiled and reanalyzed the HVS-I (hypervariable segment-I) of 3840 mitochondrial DNA (mtDNA) sequences from 19 Iranian populations and from 26 groups from adjacent countries to give a comprehensive review of the maternal genetic variation and investigate the impact of historical events and cultural factors on the maternal genetic structure of modern Iranians. We conclude that Iranians have a high level of genetic diversity. Thirty-six haplogroups were observed in Iran's populations, and most of them belong to widespread West-Eurasian haplogroups, such as H, HV, J, N, T, and U. In contrast, the predominant haplogroups observed in most of the adjacent countries studied here are H, M, D, R, U, and C haplogroups. Using principal component analysis, clustering, and genetic distance-based calculations, we estimated moderate genetic relationships between Iranian and other Eurasian groups. Further, analyses of molecular variance and comparing geographic and genetic structures indicate that mtDNA HVS-I sequence diversity does not exhibit any sharp geographic structure in the country. Barring a few from some culturally distinct and naturally separated minorities, most Iranian populations have a homogenous maternal genetic structure.

Genetic Analysis of 27 Y-STR Haplotypes in 11 Iranian Ethnic Groups

BACKGROUND

The study of Y-chromosomal variations provides valuable insights into male susceptibility in certain diseases like cardiovascular disease (CVD). In this study, we analyzed paternal lineage in different Iranian ethnic groups, not only to identify developing medical etiology, but also to pave the way for gender-specific targeted strategies and personalized medicine in medical genetic research studies.

METHODS

The diversity of eleven Iranian ethnic groups was studied using 27 Y-chromosomal short tandem repeat (Y-STR) haplotypes from Y-filer® Plus kit. Analysis of molecular variance (AMOVA) based on pair-wise RST along with multidimensional scaling (MDS) calculation and Network phylogenic analysis was employed to quantify the differences between 503 unrelated individuals from each ethnicity.

RESULTS

Results from AMOVA calculation confirmed that Gilaks and Azeris showed the largest genetic distance (RST=0.35434); however, Sistanis and Lurs had the smallest considerable genetic distance (RST=0.00483) compared to other ethnicities. Although Azeris had a considerable distance from other ethnicities, they were still close to Turkmens. MDS analysis of ethnic groups gave the indication of lack of similarity between different ethnicities. Besides, network phylogenic analysis demonstrated insignificant clustering between samples.

CONCLUSION

The AMOVA analysis results explain that the close distance of Azeris and Turkmens may be the effect of male-dominant expansions across Central Asia that contributed to historical and demographics of populations in the region. Insignificant differences in network analysis could be the consequence of high mutation events that happened in the Y-STR regions over the years. Considering the ethnic group affiliations in medical research, our results provided an understanding and characterization of Iranian male population for future medical and population genetics studies.

Religious and Spiritual Coping Elements in Dealing with Chronic Diseases: A Qualitative Exploration of the Perspectives of Older Iranian Zoroastrians

Religiosity is a desirable alternative coping strategy for many people when facing negative life events including age-related infirmities and stressors. Religious coping mechanisms (RCMs) have been investigated meagerly with regard to religious minorities around the world and, to the best of current knowledge, no study has been conducted on Iranian Zorostrians to explore their religious coping mechanisms in dealing with age-related chronic diseases. This qualitative research, therefore, was aimed to canvas perceptions about RCMs that are utilized by Iranian Zoroastrian older adults to deal with chronic diseases in the city of Yazd, Iran. Semi-structured interviews were conducted with purposefully selected fourteen Zoroastrian older patients and four Zoroastrian priests in 2019. The main extracted themes included performing certain religious behaviors and having sincere religion-based beliefs as employed mechanisms for better coping with their chronic diseases. Prevalent dilemmas/barriers with mitigating impact on the coping capacities in dealing with a persistent illness was another predominant identified theme. Identification of RCMs that religious and ethnic minorities are using to better confront diverse life events, such as chronic diseases, could pave the path to expand new approaches in planing sustainable disease management and proactive quality of life improvement initiatives.

The genetic legacy of the Yaghnobis: A witness of an ancient Eurasian ancestry in the historically reshuffled central Asian gene pool

OBJECTIVES

The Yaghnobis are an ethno-linguistic minority historically settled along the Yaghnob River in the Upper-Zarafshan Valley in Tajikistan. They speak a language of Old Sogdian origin, which is the only present-day witness of the Lingua Franca used along the Silk Road in Late Antiquity. The aim of this study was to reconstruct the genetic history of this community in order to shed light on its isolation and genetic ancestry within the Euro-Asiatic context.

MATERIALS AND METHODS

A total of 100 DNA samples were collected in the Yaghnob and Matcha Valleys during several expeditions and their mitochondrial, Y-chromosome and autosomal genome-wide variation were compared with that from a large set of modern and ancient Euro-Asiatic samples.

RESULTS

Findings from uniparental markers highlighted the long-term isolation of the Yaghnobis. Mitochondrial DNA ancestry traced an ancient link with Middle Eastern populations, whereas Y-chromosome legacy showed more tight relationships with Central Asians. Admixture, outgroup-f3, and D-statistics computed on autosomal variation corroborated Y-chromosome evidence, pointing respectively to low Anatolian Neolithic and high Steppe ancestry proportions in Yaghnobis, and to their closer affinity with Tajiks than to Iranians.

DISCUSSION

Although the Yaghnobis do not show evident signs of recent admixture, they could be considered a modern proxy for the source of gene flow for many Central Asian and Middle Eastern groups. Accordingly, they seem to retain a peculiar genomic ancestry probably ascribable to an ancient gene pool originally wide spread across a vast area and subsequently reshuffled by distinct demographic events occurred in Middle East and Central Asia.

Genotyping of 49-plex autosomal SNP panel in Iranian Turkmens ethnic group

A total of 94 unrelated individuals from Turkmens ethnic group in Iran were typed for forty-nine of the autosomal single nucleotide polymorphisms (SNPs) in the SNPforID 52plex using the SNaPshot assay. Allele frequencies are presents for the 49 SNPs. No deviation from Hardy-Weinberg equilibrium (HWE) was observed in all but one of the 49 SNP systems and no significant linkage disequilibrium was detected for any SNP pairs. F_IS and F_ST were estimated. A statistically significant global F_ST value was obtained when Turkmens ethnic group were compared with other 20 populations in Turkey, Israel, Pakistan, India, China, Taiwan, Japan, Thailand, Siberia, Algeria, Somali, Uganda, Mozambique, Angola, Nigeria, Russia, Slovenia, Sweden, France and Spain. All but 11 pairwise F_ST values were statistically significant. Multidimensional scaling plot drawn based on the pairwise F_ST values showed that the Turkmens ethnic grouped with populations geographically close to Iran and other Middle-Eastern populations. The cumulative values for the match probability using the 49 SNPs was 5.65 × 10^-19 consistent to a combined power of discrimination of >99.99999% and the mean exclusion probability was 99.95%.

The Genetic Legacy of Zoroastrianism in Iran and India: Insights into Population Structure, Gene Flow, and Selection

Zoroastrianism is one of the oldest extant religions in the world, originating in Persia (present-day Iran) during the second millennium BCE. Historical records indicate that migrants from Persia brought Zoroastrianism to India, but there is debate over the timing of these migrations. Here we present genome-wide autosomal, Y chromosome, and mitochondrial DNA data from Iranian and Indian Zoroastrians and neighboring modern-day Indian and Iranian populations and conduct a comprehensive genome-wide genetic analysis in these groups. Using powerful haplotype-based techniques, we find that Zoroastrians in Iran and India have increased genetic homogeneity relative to other sampled groups in their respective countries, consistent with their current practices of endogamy. Despite this, we infer that Indian Zoroastrians (Parsis) intermixed with local groups sometime after their arrival in India, dating this mixture to 690–1390 CE and providing strong evidence that Iranian Zoroastrian ancestry was maintained primarily through the male line. By making use of the rich information in DNA from ancient human remains, we also highlight admixture in the ancestors of Iranian Zoroastrians dated to 570 BCE–746 CE, older than admixture seen in any other sampled Iranian group, consistent with a long-standing isolation of Zoroastrians from outside groups. Finally, we report results, and challenges, from a genome-wide scan to identify genomic regions showing signatures of positive selection in present-day Zoroastrians that might correlate to the prevalence of particular diseases among these communities.

Carriers of human mitochondrial DNA macrohaplogroup M colonized India from southeastern Asia

Background

From a mtDNA dominant perspective, the exit from Africa of modern humans to colonize Eurasia occurred once, around 60 kya, following a southern coastal route across Arabia and India to reach Australia short after. These pioneers carried with them the currently dominant Eurasian lineages M and N. Based also on mtDNA phylogenetic and phylogeographic grounds, some authors have proposed the coeval existence of a northern route across the Levant that brought mtDNA macrohaplogroup N to Australia. To contrast both hypothesis, here we reanalyzed the phylogeography and respective ages of mtDNA haplogroups belonging to macrohaplogroup M in different regions of Eurasia and Australasia.

Results

The macrohaplogroup M has a historical implantation in West Eurasia, including the Arabian Peninsula. Founder ages of M lineages in India are significantly younger than those in East Asia, Southeast Asia and Near Oceania. Moreover, there is a significant positive correlation between the age of the M haplogroups and its longitudinal geographical distribution. These results point to a colonization of the Indian subcontinent by modern humans carrying M lineages from the east instead the west side.

Conclusions

The existence of a northern route, previously proposed for the mtDNA macrohaplogroup N, is confirmed here for the macrohaplogroup M. Both mtDNA macrolineages seem to have differentiated in South East Asia from ancestral L3 lineages. Taking this genetic evidence and those reported by other disciplines we have constructed a new and more conciliatory model to explain the history of modern humans out of Africa.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0816-8) contains supplementary material, which is available to authorized users.

Massive parallel sequencing of human whole mitochondrial genomes with Ion Torrent technology: an optimized workflow for Anthropological and Population Genetics studies

Investigation of human mitochondrial DNA variation patterns and phylogeny has been extensively used in Anthropological and Population Genetics studies and sequencing the whole mitochondrial genome is progressively becoming the gold standard. Among the currently available massive parallel sequencing technologies, Ion Torrent™ semiconductor sequencing represents a promising approach for such studies. Nevertheless, an experimental protocol conceived to enable the achievement of both as high as possible yield and of the most homogeneous sequence coverage through the whole mitochondrial genome is still not available. The present work was thus aimed at improving the overall performance of whole mitochondrial genomes Ion Torrent™ sequencing, with special focus on the capability to obtain robust coverage and highly reliable variants calling. For this purpose, a series of cost-effective modifications in standard laboratory workflows was fine-tuned to optimize them for medium- and large-scale population studies. A total of 54 human samples were thus subjected to sequencing of the whole mitochondrial genome with the Ion Personal Genome Machine™ System in four distinct experiments and using Ion 314 chips. Seven of the selected samples were also characterized by means of conventional Sanger sequencing for the sake of comparison. Obtained results demonstrated that the implemented optimizations had definitely improved sequencing outputs in terms of both variants calling efficiency and coverage uniformity, enabling to setup an effective and accurate protocol for whole mitochondrial genome sequencing and a considerable reduction in experimental time consumption and sequencing costs.

HLA class I variation in Iranian Lur and Kurd populations: high haplotype and allotype diversity with an abundance of KIR ligands

HLA-A, -B and -C alleles of 285 individuals, representing three Iranian Lur populations and one Iranian Kurd population were sequenced completely, yielding human leukocyte antigen (HLA) class I genotypes at high resolution and filling four fields of the official HLA nomenclature. Each population has 87-99 alleles, evenly distributed between the three HLA class I genes, 145 alleles being identified in total. These alleles were already known, named and deposited in the HLA database. The alleles form 316 different HLA A-B-C haplotypes, with each population having between 80 and 112 haplotypes. The four Iranian populations form a related group that is distinguished from other populations, including other Iranians. All four KIR ligands - the A3/11, Bw4, C1 and C2 epitopes - are well represented, particularly Bw4, which is carried by three high-frequency allotypes: HLA-A*24:02, HLA-A*32:01 and HLA-B*51:01. In the Lur and Kurd populations, between 82% and 94% of individuals have the Bw4 epitope, the ligand for KIR3DL1. HLA-B*51:01 is likely of Neandertal origin and associated with Behcet's disease, also known as the Silk Road disease. The Lordegan Lur have the highest frequency of HLA-B*51:01 in the world. This allele is present on 46 Lur and Kurd haplotypes. Present at lower frequency is HLA-B*51:08, which is also associated with Behcet's disease. In the four Iranian populations, 31 haplotypes encode both Bw4(+) HLA-A and Bw4(+) HLA-B, a dual combination of Bw4 epitopes that is relatively rare in other populations, worldwide. This study both demonstrates and emphasizes the value of studying HLA class I polymorphism at highest resolution in anthropologically well-defined populations.

Fine Dissection of Human Mitochondrial DNA Haplogroup HV Lineages Reveals Paleolithic Signatures from European Glacial Refugia

Genetic signatures from the Paleolithic inhabitants of Eurasia can be traced from the early divergent mitochondrial DNA lineages still present in contemporary human populations. Previous studies already suggested a pre-Neolithic diffusion of mitochondrial haplogroup HV*(xH,V) lineages, a relatively rare class of mtDNA types that includes parallel branches mainly distributed across Europe and West Asia with a certain degree of structure. Up till now, variation within haplogroup HV was addressed mainly by analyzing sequence data from the mtDNA control region, except for specific sub-branches, such as HV4 or the widely distributed haplogroups H and V. In this study, we present a revised HV topology based on full mtDNA genome data, and we include a comprehensive dataset consisting of 316 complete mtDNA sequences including 60 new samples from the Italian peninsula, a previously underrepresented geographic area. We highlight points of instability in the particular topology of this haplogroup, reconstructed with BEAST-generated trees and networks. We also confirm a major lineage expansion that probably followed the Late Glacial Maximum and preceded Neolithic population movements. We finally observe that Italy harbors a reservoir of mtDNA diversity, with deep-rooting HV lineages often related to sequences present in the Caucasus and the Middle East. The resulting hypothesis of a glacial refugium in Southern Italy has implications for the understanding of late Paleolithic population movements and is discussed within the archaeological cultural shifts occurred over the entire continent.

High prevalence of type 2 diabetes and pre-diabetes in adult Zoroastrians in Yazd, Iran: a cross-sectional study

Background:

The prevalence of type 2 diabetes mellitus (T2DM) varies among ethnic groups. We aimed to estimate the prevalence of diagnosed and undiagnosed diabetes mellitus, impaired fasting glucose (IFG), and impaired glucose tolerance (IGT) for the first time in an ethnic population, specifically Zoroastrian citizens in Yazd, Iran whose ages were 30 or older.

Methods:

In a cross-sectional study, participants aged≥30 years were selected using systematic random sampling. An inventory, including socio-demographic data, was completed. Weight, height, body mass index (BMI), and blood pressure (BP) were measured using standard methods. Also, blood levels of glucose, triglycerides (TG), total cholesterol (TC), high density lipoprotein (HDL), low density lipoprotein (LDL), urea, creatinine (Cr), and uric acid were measured. The latest criteria established by the American Diabetes Association (ADA) were used to diagnose DM.

Results:

The mean age of the participants (n=403) was 56.9±12.8 years. The total prevalence of diabetes, including previously diagnosed and undiagnosed diabetes, IFG, and IGT was 26.1%, 18.6%, 7.5%, 34.7% and 25.8%, respectively. Participants with diabetes had higher fasting blood sugar (FBS) (P<0.001), oral glucose tolerance test (OGTT) (P<0.001), urea (P=0.019), BMI (P=0.001), systolic blood pressures (P<0.001), TG (P=0.007) and lower HDL (P=0.034) than patients with IFG, IGT, and normoglycemic subjects.

Conclusions:

The current study showed a high prevalence of T2DM in the Zoroastrian population of Yazd, Iran. One-third of the total cases with diabetes were undiagnosed.

Evaluation of glutathione S-transferase T1 deletion polymorphism on type 2 diabetes mellitus risk in Zoroastrian females in Yazd, Iran

BACKGROUND

There has been much interest in the role of free radicals and oxidative stress in the pathogenesis of diabetes mellitus (DM). The aim of this study was to assess the possible association between genetic polymorphisms of the glutathione S-transferase-Theta (GSTT1) and the risk of the development of DM in Zoroastrian females in Yazd, Iran.

MATERIALS AND METHODS

This was a case-control study in which GSTT1 polymorphism was genotyped in 51 randomly selected DM patients and 50 randomly selected healthy controls among Zoroastrian females whose ages ranged from 40 to 70.

RESULTS

The frequencies of GSTT1 null genotype and GSTT1 present were 72% and 28%, respectively, in control samples, while in patients with type 2 diabetes (T2DM), the frequencies of GSTT1 null genotype and GSTT1 present were 27.5% and 72.5%, respectively. There were higher levels of triglyceride (TG), fasting blood sugar (FBS), total cholesterol (TC), low-density lipoprotein (LDL), Urea, and high-density lipoprotein (HDL) in cases of GSTT1 null genotype compared to the GSTT1 present genotype in controls.

CONCLUSIONS

Our results indicated that healthy subjects had a higher frequency of the GSTT1 null genotype than patients with T2DM. However, we observed no significant association between the GSTT1 null genotype and T2DM in the current study.

Increased efficiency in geographic ancestry assignment and human identification by combining lineage profiles: The case of the Iranians

OBJECTIVES

This research is a first empirical attempt to quantify the increase of the among-groups variance and the probative value of a DNA evidence when combining profiles based on markers with uniparental inheritance.

METHODS

Yfiler and HVS-I panels of loci were analyzed in 130 healthy unrelated males from six Iranian native groups.

RESULTS

A separate analysis of DNA profiles at the two lineage markers failed to detect a population substructure, whereas maximum levels of genetic diversity (HD = 1) and discrimination capacity (DC = 1) were obtained by combining the two profiles.

CONCLUSIONS

When combined, the forensic efficiency of routinely used panels of lineage markers can be largely sufficient to resolve cases of geographic ancestry and human identification even in genetically homogeneous 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.

The dual origin of Tati-speakers from Dagestan as written in the genealogy of uniparental variants

OBJECTIVES

Tat language is classified in an Iranian subbranch of the Indo-European family. It is spoken in the Caucasus and in the West Caspian region by populations with heterogeneous cultural traditions and religion whose ancestry is unknown. The aim of this study is to get a first insight about the genetic history of this peculiar linguistic group.

METHODS

We investigated the uniparental gene pools, defined by NRY and mtDNA high-resolution markers, in two Tati-speaking communities from Dagestan: Mountain Jews or Juhur, who speak the Judeo-Tat dialect, and the Tats, who speak the Muslim-Tat dialect. The samples have been collected in monoethnic rural villages and selected on the basis of genealogical relationships. A novel approach aimed at resolving cryptic cases in the recent history of human populations, which combines the properties of uniparental genetic markers with the potential of "forward-in-time" computer simulations, is presented.

RESULTS

Judeo-Tats emerged as a group with tight matrilineal genetic legacy who separated early from other Jewish communities. Tats exhibited genetic signals of a much longer in situ evolution, which appear as substantially unlinked with other Indo-Iranian enclaves in the Caucasus.

CONCLUSIONS

The independent demographic histories of the two samples, with mutually reversed profiles at paternally and maternally transmitted genetic systems, suggest that geographic proximity and linguistic assimilation of Tati-speakers from Dagestan do not reflect a common ancestry.