Abstract

Introduction. β-Thalassaemia is the most common genetic disorder and is considered as a major public health concern in Iran. Different countrywide studies have shown a heterogeneous mutational basis of β-thalassaemia with different frequencies in each area. This study is aimed at investigating the common and rare mutations in Mazandaran and Golestan, northern provinces of Iran. Methods. 5425 microcytic and hypochromic individuals were investigated from Mazandaran and Golestan provinces. From these, 1323 beta carrier or affected individuals were selected where 938 persons were from Mazandaran and 385 people were from Golestan province, respectively. Result. 53 different mutations were identified, IVSII-1 (G>A) was the most common (59.14%) followed by Cd 22/23/24 (-7 bp) (5.34%), Cd 8 (-AA) (4.93%), Cd30 (G>A) (4.00%), and IVSI-5 (G>C) (3.70%) with a total of 77.11% in Mazandaran Province, respectively. In Golestan Province, IVSI-5 (G>C) was the most frequent (44.62%) followed by IVSII-1 (G>A) (27.18%), Cd 15 (TGG>TAG) (4.36%), Fr 8/9(+G) (3.85%), and Cd 8(-AA) (2.05%) with a total of 82.06%, respectively. From the 53 different mutations, 22 numbers have been observed in both provinces. Two deletions of the beta gene named Sicilian and Asian-Indian have been detected in Mazandaran with a frequency of 0.72% each. Conclusion. The 53 different mutations identified in this study were the most ever reported mutations in the country. Due to diversity of different ethnic groups, there are many varieties of mutation in beta globin gene in Iran. It could be assumed that both founder effect and natural selection caused by migration from neighboring areas have complemented each other to produce the high frequency of unique alleles within each region.

1. Introduction

β-Thalassaemia is the most common genetic disorder and is considered as a major public health concern in Iran. It is estimated that up to 7% of the world population are carriers of hemoglobin diseases, such as β- or α-zero thalassemia (α0-thal), or hemoglobin S, C, D Punjab, or E [1]. The prevalence of β-thalassemia is high in the Middle East, Southeast Asia, and Mediterranean Region and increasing in nonendemic regions including North America and Western Europe due to migration pattern [2]. The gene frequency of β-thalassaemia is high in Iran and varies considerably from area to area. It is believed that the highest rate is around the Caspian Sea with up to 8% [35]. Nearly 2800 patients are register in health care services in Mazandaran and about 402 patients in Golestan Province according to Thalassemia Registry in Mazandaran Province (http://thr.mazums.ac.ir). According to previous studies, the frequency of carriers is between 2 and 9.5% in different provinces, of which Mazandaran and Golestan are among the highest number of carriers with about 8% [6, 7]. One study in Golestan Province showed a frequency of 9.8% for beta-thalassemia carriers [8]. The high prevalence of β-thalassaemia in this region is believed to be related to endemic malaria [9]. National thalassemia screening program that is a premarriage screening started since 1995, which initially uses a cell blood count (CBC) and hemoglobin electrophoresis tests from the couple before marriage, to determine beta-thalassemia carriers as well as anemic people due to hemoglobin defects such as sickle cell anemia, hemoglobin E, and alpha thalassemia [10, 11].

Prenatal diagnosis (PND) requires a genetic test to identify defects or mutations in parents followed by mutation detection in fetus after pregnancy [12]. Performing PND for the fetus and medical abortion considerably reduced the incidence of thalassemia in Iran, including Mazandaran and Golestan provinces [6, 7, 11, 13].

A long list including 950 different mutations is reported in hemoglobin database HbVar or “A Database of Human Hemoglobin Variants and Thalassemias” from different countries to date (http://globin.bx.psu.edu/hbvar/menu.html). Also, several countrywide studies have shown that the mutational basis of β-thalassaemia in Iran is heterogeneous, and more than 50 different mutations are detected across different geographic regions of Iran [1419]. The spectrum of these mutations is variable in each area and depends on geographic regions, ethnic group, and immigration. Few reports are available regarding the spectrum and frequency of mutations causing β-thalassaemia in Mazandaran and Golestan provinces to date [2022]. This study is aimed at investigating the common mutations in the region. This data will help us to update our knowledge about common as well as rare mutations in HBB gene in northern provinces of the country. It also will be helpful for β-thalassaemia carriers screening, genetic counseling, and prenatal diagnosis of the disease in these area as well as in the country.

2. Materials and Methods

2.1. Study Population

As part of the national thalassemia screening program, high-risk couples or couples who were already detected as beta-thalassemia carrier referred to a designated local laboratory for genetic tests. It is a retrospective and cross-sectional study, with nonprobability and consecutive sampling between 2008 and 2021 among the individuals who were referred to a diagnostic genetic laboratory. People with beta globin gene mutation are included, and those with alpha globin gene mutation or without mutation were excluded from this study.

People were first screened by CBC (cell blood count) and hemoglobin electrophoresis. Those who had reduced MCV ( fl) or MCH ( pg/cell) [23] were subjected to further study and referred to a local genetic laboratory as a routine national services [10]. Also, for the HbA2 level, concentration more than 3.5% was considered as beta-thalassemia carrier [10, 24], and HbF hemoglobin normal level in adults is considered from 0.5% to 1.5% [24, 25]. From 5425 microcytic and hypochromic people referred to “Novin Genetic Diagnostic Laboratory” during 2008 to 2021, 1323 persons were beta carrier or affected individuals, who were referred from Mazandaran or Golestan provinces (Table 1). All individuals were informed about this research, and written consent was obtained from each patient; meanwhile, name and personal data of the participants will remain confidential.

2.2. Genomic DNA Extraction

A 5 to 10 ml venous blood was obtained from each subject and stored in Na-EDTA tube at -25°C until processing. Lymphocytic genomic DNA was extracted by Nucleon BACCII method [26] followed by DNA concentration measurement using NanoDrop2000/2000c Spectrophotometers (Thermo Fisher Scientific, USA).

2.3. ARMS-PCR Amplification for Common Mutations

The specific primers were used to amplify each specific mutation separately. The ARMS-PCR was applied for the detection of 10 common mutations previously reported including IVSII-1 (G>A), IVSI-5 (G>C), Cd 22/23/24 (del GAAGTTGGT), Cd 30 (G>C), Cd 8 (-AA), Fr 8/9 (+G), Cd 44 (-C), IVSI-I (G>A), Cd 36/37 (-T), and IVSI-25bp del [19, 27]. Briefly, three specific primers were used to genotype the mentioned mutations, two specific forward primers for normal and mutant alleles, respectively, and one common reverse primer. For each sample, two tubes were used as M and N. A 25 μl master mix PCR containing 11 μl “Taq DNA Polymerase 2x Master Mix RED” (Amplicon, Denmark), 0.6 μl of forward and 0.6 μl reverse primer, and 2 μl DNA template up to 25 μl dH2O were applied in each tube. The PCR reactions were carried out with the following conditions: 93°C, 40 seconds; annealing temperature for 40 seconds; and 72°C 40 seconds for 30 cycles. PCR products were visualized on 1% agarose gel. Results were evaluated as normal when only N tube amplified, homozygous mutant (MM) when only M tube amplified, and heterozygous normal/mutant (NM), when both tubes amplified (Figure 1(a)).

2.4. DNA Sequencing

Samples that showed no mutations using conventional ARMS-PCR methods were subjected to further investigations using Sanger DNA sequencing method (Applied Biosystems, Waltham, USA). In brief, DNA sequencing method was used to screen entire coding and noncoding regions of the HBB gene using specific primers (Table 2) to find any point mutations or short deletion/insertions. A DNA sequence analysis software, Gene Runner (http://www.generunner.com), was applied along with using reference sequences from GenBank database. Finch TV, a DNA sequence chromatogram viewer software (Geospiza, Inc., USA), was also applied to view nucleotide changes.

2.5. Gap-PCR and Multiplex Ligation-Dependent Probe Amplification (MLPA)

Samples that showed no mutations using different methods were subjected for further study using Gap-PCR and the MLPA method (MRC-Holland, Netherland) to identify any potential large deletion or insertion [28]. Some frequent deletions including Hb Lepore, Hb Sicilian, and Hb Asian-Indian were detected using previously described Gap-PCR [29]. Subsequently, in those that mutation was not identified, MLPA assay using SALSA MLPA Probemix P102 HBB (MRC-Holland, Amsterdam, Netherlands) was applied searching for deletions or gene fusion within the HBB cluster [12, 28, 30].

3. Result

5425 microcytic and hypochromic people who were referred to the Novin Genetic Diagnostic Laboratory during 2008 to 2021 were investigated. In total, 4500 persons were from Mazandaran and 925 individuals were from Golestan Province, respectively. From these, 1323 beta carrier or affected individuals were selected where 938 persons were from Mazandaran (902 carrier and 36 affected) and 385 persons were from Golestan Province (380 carrier and 5 affected), respectively (Table 1). 53 different mutations were identified in this study (Table 3). From the 974 mutant chromosomes investigated in Mazandaran Province, IVSII-1 (G>A) was the most common mutation with a frequency of 59.14% followed by four mutations including Cd 22/23/24 (-7 bp) (5.34%), Cd 8 (-AA) (4.93%), C30 (G>A) (4.00%), and IVSI-5 (G>C) (3.70%), respectively. These 5 mutations represented 77.11% of the total mutations in Mazandaran Province (Table 3). From the 4500 individuals investigated from Mazandaran Province, 28 (0.62%) individuals were HbS patients including 27 heterozygous carrier and 1 homozygous and 2 sickle-thal affected, respectively. We also found 49 individuals (1.08%) with HbD Punjab (Cd 121 GAA>CAA) variant including two HbD homozygous in Mazandaran. Also, 5 individuals (0.11%) were identified with HbE (Cd 26 GAG>TAG) variant in Mazandaran population among the 4500 individuals.

In Golestan Province, IVSI-5 G>C was the most frequent mutation with the frequency of 44.62% followed by IVSII-1 (G>A) (27.18%), Cd 15 (TGG>TAG) (4.36%), Fr 8/9 (+G) (3.85%), and Cd 8 (-AA) (2.05%), respectively. These 5 mutations represented 82.06% of the total mutations in Golestan Province (Table 3). Also, among the 925 individuals, there were 2 (0.22%) carrier individuals for HbS hemoglobin and 28 individuals (3.03%) including 27 heterozygous and 1 homozygous for HbD Punjab (Cd 121 GAA>CAA) hemoglobin. Also, 4 beta-thalassemia carrier had HbD hemoglobin as well. Also, 3 individuals (0.32%) with HbE (Cd26 GAG>TAG) variant were identified in Golestan Province, respectively.

Among the beta-thalassemia carriers that were studied in this research, some showed normal HbA2 level between 1 and 3.5%. Table 4 shows 92 persons from Mazandaran and 20 from Golestan Province with normal level of HbA2 hemoglobin. Although some previously known mutations like -101, UTR+22, are IVSI-6 are listed among these mutations, some beta-zero mutations like IVSII-1 and IVSI-5 can also be seen among this list (Table 4).

4. Discussion

Beta-thalassemia is a heterogeneous and diverse disease from molecular point of view. In HbVar hemoglobin database (http://globin.bx.psu.edu/hbvar/menu.html), 950 variant alleles have been reported in beta globin gene to date [31]. A small number of thalassemia mutations are predominant in most parts of the world, and the most common ones tend to be those that are geographically the most widespread [32]. For instance, six common mutations including IVSII-745 (C>G), IVSI-110 (G>A), IVSI-6 (T>C), IVSI-1 (G>A), IVSII-1 (G>A), and C39 (C>T) represent about 90% of all mutations in the Mediterranean area [33], and four alleles account for about 91% of the mutations in China and Southeast Asia [34]. Iran with a population of approximately 85 million represents a highly heterogeneous gene pool and mutation spectrum due to geographical, cultural, and ethnical diversity [14, 19, 35].

Beta-thalassemia gene mutation distribution and frequency are well studied and reported in different provinces in Iran. For example, a study in Hormozgan Province (south of Iran) reported IVSI-5 (G>C) (69%) as a most frequent, followed by IVSII-1 (G>A) (9.6%) [36]. Other study in Azerbaijan Province (northwest of the country) showed that IVSII-1 (G>A) (21%) and IVSI-110 (G>A) (18%) were more frequent [37]. Kurdistan Province in the west also showed IVSII-1 (G>A) with 31% as the most frequent mutation [38]. In the southwest province of Khuzestan, IVSII-1 (G>A) with 34% was the most frequent, followed by Fr 8/9 with 17.3% [39]. A recent study on the central province of Isfahan showed that IVSII-1 (G>A) (27.9%) and Fr 36/37(-T) (19.7%) were the most frequent mutations, respectively [40]. IVSII-1 (G>A) is a common mutation in the Mediterranean area and was reported as the most predominant in three northern (Gilan, Mazandaran, and Golestan) and northwestern (Azerbaijan and Ardebil) provinces of Iran [20, 37].

In this study, 53 different mutations were identified. From these, 22 mutations have been observed in both provinces, including IVSII-1 (G>A), IVSI-5 (G>C), Fr8/9 (+G), IVSI-25 del, C8 (-AA), IVSI-130 (G>C), IVS-110 (G>A), and IVSI-1 (G>A) which were more frequent, respectively (Table 3). Except of these common, about 31 remaining are detected in either both or one province with a frequency of less than 1% (Table 3). While IVSII-1 (G>A) was the most common mutation in Mazandaran (59.14%), IVSI-5 (G>C) was the most prevalent in Golestan Province (44.62%) followed by IVSII-1 (G>A) (27.18%) instead. One new or unreported mutation also is identified in a carrier person from Golestan Province in this study, a base pair deletion in codon 80, c.243 C>-, and C80 (AAC>AA-). From the 5 most common mutations detected in the two provinces, 3 of them are the same (IVSII-1, IVSI-5, and Cd 8), but with different frequencies (Table 3).

In HBB gene, point mutations are prevalent, but still, there are some gene deletions which are detected and may be different from area to area. For instance, two gross deletions named Sicilian and Asian-Indian have been detected in Mazandaran Province with frequency of 0.72% and 0.72% in this study (Table 3) [29]. In addition, two other unknown deletions identified by MLPA method, deletion of exon 3 of beta globin gene, and a large deletion with about 132 kb length (from 5, A-gamma gene, beta gene, and delta gene toward TRIM68 gene in 3) were recognized in this study and need to be fully characterized later (Figure 2). The carrier person with large gene deletion had the following hematologic index: RBC: 6.35, Hbg: 17.1, MCV: 80.3, MCH: 26.9, HbA1: 79.4, HbA2: 4.4, and HbF: 16.2, respectively.

Different investigations revealed IVSI-5 (C>G) as the most prevalent mutation in the Middle East, India, and South and Southeast Asia before [32, 34, 41]. The IVSI-5 is the most common mutation among the people with Sistani and Baluch ethnic origin, from Sistan-Baluchestan Province, a southeastern part of Iran which has border with Pakistan [42]. Earlier study in this province reported IVSI-5 (G>A) (44.8%) as the most common mutation, similar to the frequency in Pakistan (53%) [43]. Because there was internal migration or relocation from Sistan and Baluchestan provinces to Golestan in recent decades and they are made up to 30% of the total population in Golestan Province now [44], so as a consequence, IVSI-5 is one of the frequent mutations in this province. Surprisingly, it is only found 3.7% in Mazandaran Province in the present study. Additionally, our results revealed the frequency of 3 hemoglobin variants, HbS, HbD, and HbE. In Mazandaran Province, frequency of HbS was 0.62%, HbD 1.08%, and HbE about 0.32%, respectively, in a total of 4500 people studied. In Golestan Province, frequency of HbS was found 0.22%, HbD 3.03%, and HbE about 0.32% in 925 individuals, respectively.

CBC index detection is the first line of laboratory examination to find carrier individual in Iranian national thalassemia program [23]. In the present study, carrier individuals with different mutations generally showed the mean MCV below 70 fl and MCH below 20 pg per cell, respectively (Table 5), but for the mutations such as -88, -101, 5-UTR+22, CD27, CD126, IVSI-6, and IVSI-128, these values achieved were usually above 70 and 20, respectively. Moreover, for carriers of two deletion mutations, Sicilian and Asian-Indian, the MCV and MCH values were higher than 70 and 20 (Table 5). Furthermore, the HbA2 level in beta-thalassemia carrier is more than the normal level of 3.5% [23], but in this study, there were some carrier individuals with normal HbA2 level. Table 4 lists a number of mutations including -101, 5-UTR+22, CD27, CD126, IVSI-6, and IVSI-128. The Sicilian and Asian-Indian deletion mutations also have normal HbA2 level and increased HbF with averages of 6.85% and 6.4%, respectively (Table 4). Some beta-thalassemia carriers with common mutations such as IVSII-1 and IVSI-5 also show normal HbA2 level with unknown reason (Table 4). There were some limitations in this study. For instance, some families were not fully cooperative for additional investigations, so the blood samples from other family members were not available. Also because of the limitation in the budget, it was not possible to investigate and determine the exact length of the deletion in individuals who had new gene deletion.

5. Conclusion

The 53 different heterogeneous mutations identified among carriers and patients in this study were the most ever reported mutations both in the number and diversity compared with the previous reports in the country. Two previous studies in Mazandaran reported 17 and 25 different mutations in beta-thalassemia carrier and patients, respectively [20, 21]. Although the diversity and number of mutations reported here are higher, the frequency of common mutations was similar with the previous studies (Table 3). Five mutations were more common in both provinces in which three of them, IVSII-1, IVSI-5, and Cd 8, were similar between two provinces.

Due to presence and diversity of different ethnic groups in the country, there are many varieties in the molecular genetics of beta globin gene in Iran. The presence of different highly frequent mutants’ alleles, along with some similarities in different geographical regions, is supporting a role of nonisolating genetically areas. It could be assumed that both founder effect and natural selection caused by migration from neighboring areas have complemented each other to produce the high frequency of unique alleles within each region.

Data Availability

The data used to support the findings of this study are included within the article. The data were supplied by a diagnostic laboratory and so cannot be made freely available. Requests for access to these data should be made to the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors are extremely grateful to all health centers and health care staff from Mazandaran and Golestan provinces as well as staff of the Novin Genetic Laboratory, Mazandaran Medical Sciences, Sari, for their assistance and cooperation. This project is supported by Novin Genetic Diagnostic Laboratory.