Phylogeny and genetic structure of Tunisians and their position within Mediterranean populations

Genetic diversity of variants involved in drug response among Tunisian and Italian populations toward personalized medicine

Adverse drug reactions (ADR) represent a significant contributor to morbidity and mortality, imposing a substantial financial burden. Genetic ancestry plays a crucial role in drug response. The aim of this study is to characterize the genetic variability of selected pharmacogenes involved with ADR in Tunisians and Italians, with a comparative analysis against global populations. A cohort of 135 healthy Tunisians and 737 Italians were genotyped using a SNP array. Variants located in 25 Very Important Pharmacogenes implicated in ADR were extracted from the genotyping data. Distribution analysis of common variants in Tunisian and Italian populations in comparison to 24 publicly available worldwide populations was performed using PLINK and R software. Results from Principle Component and ADMIXTURE analyses showed a high genetic similarity among Mediterranean populations, distinguishing them from Sub-Saharan African and Asian populations. The Fst comparative analysis identified 27 variants exhibiting significant differentiation between the studied populations. Among these variants, four SNPs rs622342, rs3846662, rs7294, rs5215 located in SLC22A1, HMGCR, VKORC1 and KCNJ11 genes respectively, are reported to be associated with ethnic variability in drug responses. In conclusion, correlating the frequencies of genotype risk variants with their associated ADRs would enhance drug outcomes and the implementation of personalized medicine in the studied populations.

Genetic landscape of forensic DNA phenotyping markers among Mediterranean populations

Forensic DNA Phenotyping can reveal the appearance of an unknown individual by predicting the External Visible Characteristics (EVC) from DNA obtained at the crime scene. Our aim is to characterize the genetic landscape of Human identification markers responsible for EVC among Mediterranean populations compared to other worldwide groups. We conducted an exhaustive search for genes involved in EVC variation. Then, variants located on these genes were extracted from public genotypic data of Mediterranean, American, African and East Asiatic populations. The genetic landscape of these Human identification markers, their allelic distribution and admixture analyses, were determined using plink, R and ADMIXTURE softwares. Our results showed that the Mediterranean populations appear close to the Mexican populations and distinguished from sub Saharan African populations living in the USA and from East Asiatic populations. We highlighted a total of 103454 common variants shared between the studied populations and among them, 25 common variants associated with EVC. Interestingly, genotype frequencies results showed that the rs17646946, rs13016869, rs977588, rs1805008 and rs2240751 variants located respectively in the TCHH, PRKCE, OCA2, MC1R and MFSD12 genes are significantly different between the Mediterranean and Asiatic populations. The genotype frequencies of the variants rs977589 and rs7179994 located in the OCA2 gene, and of rs12913832 and rs2240751 located respectively in HERC2 and MFSD12 genes are significantly different between the Mediterranean and American populations. Our work generates a large number of EVC variants that could be a valuable resource for future studies in the forensic field.

The first exome wide association study in Tunisia: identification of candidate loci and pathways with biological relevance for type 2 diabetes

Introduction

Type 2 diabetes (T2D) is a multifactorial disease involving genetic and environmental components. Several genome-wide association studies (GWAS) have been conducted to decipher potential genetic aberrations promoting the onset of this metabolic disorder. These GWAS have identified over 400 associated variants, mostly in the intronic or intergenic regions. Recently, a growing number of exome genotyping or exome sequencing experiments have identified coding variants associated with T2D. Such studies were mainly conducted in European populations, and the few candidate-gene replication studies in North African populations revealed inconsistent results. In the present study, we aimed to discover the coding genetic etiology of T2D in the Tunisian population.

Methods

We carried out a pilot Exome Wide Association Study (EWAS) on 50 Tunisian individuals. Single variant analysis was performed as implemented in PLINK on potentially deleterious coding variants. Subsequently, we applied gene-based and gene-set analyses using MAGMA software to identify genes and pathways associated with T2D. Potential signals were further replicated in an existing large in-silico dataset, involving up to 177116 European individuals.

Results

Our analysis revealed, for the first time, promising associations between T2D and variations in MYORG gene, implicated in the skeletal muscle fiber development. Gene-set analysis identified two candidate pathways having nominal associations with T2D in our study samples, namely the positive regulation of neuron apoptotic process and the regulation of mucus secretion. These two pathways are implicated in the neurogenerative alterations and in the inflammatory mechanisms of metabolic diseases. In addition, replication analysis revealed nominal associations of the regulation of beta-cell development and the regulation of peptidase activity pathways with T2D, both in the Tunisian subjects and in the European in-silico dataset.

Conclusions

The present study is the first EWAS to investigate the impact of single genetic variants and their aggregate effects on T2D risk in Africa. The promising disease markers, revealed by our pilot EWAS, will promote the understanding of the T2D pathophysiology in North Africa as well as the discovery of potential treatments.

Whole-exome sequencing reveals novel variants of monogenic diabetes in Tunisia: impact on diagnosis and healthcare management

Introduction: Monogenic diabetes (MD) accounts for 3%-6% of all cases of diabetes. This prevalence is underestimated due to its overlapping clinical features with type 1 and type 2 diabetes. Hence, genetic testing is the most appropriate tool for obtaining an accurate diagnosis. In Tunisia, few cohorts of MD have been investigated until now. The aim of this study is to search for pathogenic variants among 11 patients suspected of having MD in Tunisia using whole-exome sequencing (WES). Materials and methods: WES was performed in 11 diabetic patients recruited from a collaborating medical center. The pathogenicity of genetic variation was assessed using combined filtering and bioinformatics prediction tools. The online ORVAL tool was used to predict the likelihood of combinations of pathogenic variations. Then, Sanger sequencing was carried out to confirm likely pathogenic predicted variants among patients and to check for familial segregation. Finally, for some variants, we performed structural modeling to study their impact on protein function. Results: We identified novel variants related to MD in Tunisia. Pathogenic variants are located in several MODY and non-MODY genes. We highlighted the presence of syndromic forms of diabetes, including the Bardet-Biedl syndrome, Alström syndrome, and severe insulin resistance, as well as the presence of isolated diabetes with significantly reduced penetrance for Wolfram syndrome-related features. Idiopathic type 1 diabetes was also identified in one patient. Conclusion: In this study, we emphasized the importance of genetic screening for MD in patients with a familial history of diabetes, mainly among admixed and under-represented populations living in low- and middle-income countries. An accurate diagnosis with molecular investigation of MD may improve the therapeutic choice for better management of patients and their families. Additional research and rigorous investigations are required to better understand the physiopathological mechanisms of MD and implement efficient therapies that take into account genomic context and other related factors.

Decoding the genetic relationship between Alzheimer's disease and type 2 diabetes: potential risk variants and future direction for North Africa

Introduction

Alzheimer's disease (AD) and Type 2 diabetes (T2D) are both age-associated diseases. Identification of shared genes could help develop early diagnosis and preventive strategies. Although genetic background plays a crucial role in these diseases, we noticed an underrepresentation tendency of North African populations in omics studies.

Materials and methods

First, we conducted a comprehensive review of genes and pathways shared between T2D and AD through PubMed. Then, the function of the identified genes and variants was investigated using annotation tools including PolyPhen2, RegulomeDB, and miRdSNP. Pathways enrichment analyses were performed with g:Profiler and EnrichmentMap. Next, we analyzed variant distributions in 16 worldwide populations using PLINK2, R, and STRUCTURE software. Finally, we performed an inter-ethnic comparison based on the minor allele frequency of T2D-AD common variants.

Results

A total of 59 eligible papers were included in our study. We found 231 variants and 363 genes shared between T2D and AD. Variant annotation revealed six single nucleotide polymorphisms (SNP) with a high pathogenic score, three SNPs with regulatory effects on the brain, and six SNPs with potential effects on miRNA-binding sites. The miRNAs affected were implicated in T2D, insulin signaling pathways, and AD. Moreover, replicated genes were significantly enriched in pathways related to plasma protein binding, positive regulation of amyloid fibril deposition, microglia activation, and cholesterol metabolism. Multidimensional screening performed based on the 363 shared genes showed that main North African populations are clustered together and are divergent from other worldwide populations. Interestingly, our results showed that 49 SNP associated with T2D and AD were present in North African populations. Among them, 11 variants located in DNM3, CFH, PPARG, ROHA, AGER, CLU, BDNF1, CST9, and PLCG1 genes display significant differences in risk allele frequencies between North African and other populations.

Conclusion

Our study highlighted the complexity and the unique molecular architecture of North African populations regarding T2D-AD shared genes. In conclusion, we emphasize the importance of T2D-AD shared genes and ethnicity-specific investigation studies for a better understanding of the link behind these diseases and to develop accurate diagnoses using personalized genetic biomarkers.

Minor histocompatibility antigens HA-8 and PANE1 in the TUNISIAN population

BACKGROUND

Minor histocompatibility antigens (mHAgs) are endogenous immunogenic peptides initially identified due to complications detected in several contexts of HLA geno-identical hematopoietic stem cell transplantation (HSCT). In this study, we chose to examine the molecular polymorphism of the mHAgs HA-8 and PANE1 in the Tunisian population.

MATERIAL AND METHODS

This study was conducted on 150 healthy and unrelated individuals. The DNA extraction and Sequence-Specific Primers PCR (PCR-SSP) methods were used for the molecular genotyping of the selected SNPs: PUM3 (rs2173904) and CENPM (rs5758511).

RESULTS

Our results show that, 94% of Tunisians are carriers of the PANE1R allele (immunogenic variant of the PANE1 mHAg) and 68% of Tunisians are carriers of the HA-8R allele (immunogenic variant of the HA-8 mHAg). Furthermore, this study shows that about 5% of the Tunisians are carrier of the PANE1R antigen and its HLA molecule of presentation (the PANE1R/HLA-A*0301 combination). However, only 2% of Tunisians are carrier of the HA-8R/HLA-A*0201 combination, that is, the HA8 immunogenic variant and its specific HLA molecule of presentation.

CONCLUSION

Our results are close to those reported in Caucasian, Asiatic, and African populations, this may be explained by the historical events experienced by Tunisia for millennia. These results could be used for further clinical and anthropological studies.

Health influenced by genetics: A first comprehensive analysis of breast cancer high and moderate penetrance susceptibility genes in the Tunisian population

Significant advances have been made to understand the genetic basis of breast cancer. High, moderate and low penetrance variants have been identified with inter-ethnic variability in mutation frequency and spectrum. Genome wide association studies (GWAS) are widely used to identify disease-associated SNPs. Understanding the functional impact of these risk-SNPs will help the translation of GWAS findings into clinical interventions. Here we aim to characterize the genetic patterns of high and moderate penetrance breast cancer susceptibility genes and to assess the functional impact of non-coding SNPs. We analyzed BRCA1/2, PTEN, STK11, TP53, ATM, BRIP1, CHEK2 and PALB2 genotype data obtained from 135 healthy participants genotyped using Affymetrix Genome-Wide Human SNP-Array 6.0. Haplotype analysis was performed using Haploview.V4.2 and PHASE.V2.1. Population structure and genetic differentiation were assessed using principal component analysis (PCA) and fixation index (FST). Functional annotation was performed using In Silico web-based tools including RegulomeDB and VARAdb. Haplotype analysis showed distinct LD patterns with high levels of recombination and haplotype blocks of moderate to small size. Our findings revealed also that the Tunisian population tends to have a mixed origin with European, South Asian and Mexican footprints. Functional annotation allowed the selection of 28 putative regulatory variants. Of special interest were BRCA1_ rs8176318 predicted to alter the binding sites of a tumor suppressor miRNA hsa-miR-149 and PALB2_ rs120963 located in tumorigenesis-associated enhancer and predicted to strongly affect the binding of P53. Significant differences in allele frequencies were observed with populations of African and European ancestries for rs8176318 and rs120963 respectively. Our findings will help to better understand the genetic basis of breast cancer by guiding upcoming genome wide studies in the Tunisian population. Putative functional SNPs may be used to develop an efficient polygenic risk score to predict breast cancer risk leading to better disease prevention and management.

Association of HNF1A gene variants and haplotypes with metabolic syndrome: a case-control study in the Tunisian population and a meta-analysis

BACKGROUND

Variants in the Hepatocyte Nuclear Factor 1 Alpha gene (HNF1A) are associated with lipoproteins levels and type 2 diabetes. In this study, we aimed to assess the association of HNF1A gene and haplotypes with the metabolic syndrome (MetS) and its components through an association study in the Tunisian population as well as by a meta-analysis.

METHODS

A total of 594 Tunisian individuals were genotyped for three variants (rs1169288, rs2464196 and rs735396) located in HNF1A gene using KASPar technology. Statistical analyses were performed with R software. The association was furthermore evaluated through a meta-analysis of our results with those obtained in a Moroccan population.

RESULTS

Our results showed no association between HNF1A variants and MetS in the Tunisian population. However, a significant association was observed between the variant rs735396 and a higher waist circumference. The stratified analysis according to the sex highlighted a significant association between the variant rs1169288 and high cholesterol levels only in women. Similarly, Haplotype analysis showed an association between the HNF1A minor haplotype and high total cholesterol mainly in women. Finally, our meta-analysis showed no association between HNF1A variants and MetS.

CONCLUSIONS

Our findings exclude the involvement of the three HNF1A variants rs1169288, rs2464196 and rs735396 in the susceptibility to MetS in our studied Tunisian population but emphasize the role of these variants in the cholesterol homeostasis with sex-specific differences, which may serve to rise clinical consideration to early statin therapy in women carrying these genetic variants.

Correlation between facial measurements and vertical dimension of occlusion among Tunisian populations: An anthropometric study

Background. Establishing an accurate occlusal vertical dimension (OVD) is a crucial clinical step during full-mouth rehabilitation. Various techniques have been suggested to evaluate OVD, but none of them is practically reliable, and each one has its shortcomings. The correlation between facial proportions and the lower third of the face is a reliable method but needs to be verified in many ethnic groups. Therefore, this study aimed to determine the correlation between OVD and various facial measurements in a Tunisian ethnic group. Methods. A cross-sectional study was conducted between November 2020 to January 2021. The participants were randomly selected from dental students, dental surgeons, and the patients referring to the University Dental Clinic for dental treatments. Seven facial measurements were clinically recorded using a digital caliper. The correlation between OVD and facial measurements was analyzed using Spearman's coefficient and linear regression analysis. Results. A total of 201 dentate participants (134 females and 67 males) were included in the study. The mean OVD in male subjects was higher (67.60±4.49) compared to female subjects (60.72±3.84). The total facial height was positively correlated with OVD in both genders. OVD was statistically correlated with the height of the upper lip. This correlation was highly significant in males while it was weak in the female group. Conclusion. Facial proportions and linear equations are non-invasive, simple, and reliable methods to predict OVD, especially in males.

Spectrum of Genetic Diseases in Tunisia: Current Situation and Main Milestones Achieved

Genetic diseases in Tunisia are a real public health problem given their chronicity and the lack of knowledge concerning their prevalence and etiology, and the high rates of consanguinity. Hence, we performed systematic reviews of the literature in order to provide a more recent spectrum of these disorders and to expose the challenges that still exist to tackle these kinds of diseases. A manual textual data mining was conducted using MeSH and PubMed databases. Collected data were classified according to the CIM-10 classification and the transmission mode. The spectrum of these diseases is estimated to be 589 entities. This suggests remarkable progress through the development of biomedical health research activities and building capacities. Sixty percent of the reported disorders are autosomal recessive, which could be explained by the high prevalence of endogamous mating. Congenital malformations (29.54%) are the major disease group, followed by metabolic diseases (22%). Sixty percent of the genetic diseases have a known molecular etiology. We also reported additional cases of comorbidity that seem to be a common phenomenon in our population. We also noticed that epidemiological data are scarce. Newborn and carrier screening was only limited to pilot projects for a few genetic diseases. Collected data are being integrated into a database under construction that will be a valuable decision-making tool. This study provides the current situation of genetic diseases in Tunisia and highlights their particularities. Early detection of the disease is important to initiate critical intervention and to reduce morbidity and mortality.

Anti-gp210 and anti-Sp100 antibodies in primary biliary cholangitis

BACKGROUND AND STUDY AIMS

To determine the sensitivity and specificity of anti-gp210 and anti-Sp100 autoantibodies in primary biliary cholangitis (PBC) PATIENTS AND METHODS: Sera of 106 PBC patients with positive anti-mitochondrial antibodies and 58 healthy blood donors were analyzed. A line immunoassay was used to evaluate the reactivity of anti-gp210 and anti-Sp100 antibodies.

RESULTS

The frequency of anti-gp210 and anti-Sp100 autoantibodies was 29.2% and 28.3%, respectively. Eight patients had both anti-gp210 and anti-Sp100 antibodies. Of 106 patients, 23 (21.7%) had anti-gp210 antibody, although not anti-Sp100 antibody, and 22 (20.7%) had anti-Sp100, although not anti-gp210 antibodies. Their combination increased the frequency of anti-gp210 and anti-Sp100 antibodies from 29.2% to 50% (P = 0.002) and 28.3% to 50% (P = 0.0012), respectively. In the control group, two subjects had anti-gp210 antibody and none had anti-Sp100 antibody. Thus, the specificity of anti-gp210 and anti-Sp100 antibodies was 96.5% and 100%, respectively. The positive predictive value (PPV) of anti-gp210 antibody was 94%; its negative predictive value (NPV) was 42.7%. The PPV and NPV of anti-Sp100 antibody were 100% and 43.3%, respectively.

CONCLUSION

It is important to combine anti-gp210 and anti-Sp100 antibodies in the immunological exploration of PBC.

A map of copy number variations in the Tunisian population: a valuable tool for medical genomics in North Africa

Copy number variation (CNV) is considered as the most frequent type of structural variation in the human genome. Some CNVs can act on human phenotype diversity, encompassing rare Mendelian diseases and genomic disorders. The North African populations remain underrepresented in public genetic databases in terms of single-nucleotide variants as well as for larger genomic mutations. In this study, we present the first CNV map for a North African population using the Affymetrix Genome-Wide SNP (single-nucleotide polymorphism) array 6.0 array genotyping intensity data to call CNVs in 102 Tunisian healthy individuals. Two softwares, PennCNV and Birdsuite, were used to call CNVs in order to provide reliable data. Subsequent bioinformatic analyses were performed to explore their features and patterns. The CNV map of the Tunisian population includes 1083 CNVs spanning 61.443 Mb of the genome. The CNV length ranged from 1.017 kb to 2.074 Mb with an average of 56.734 kb. Deletions represent 57.43% of the identified CNVs, while duplications and the mixed loci are less represented. One hundred and three genes disrupted by CNVs are reported to cause 155 Mendelian diseases/phenotypes. Drug response genes were also reported to be affected by CNVs. Data on genes overlapped by deletions and duplications segments and the sequence properties in and around them also provided insights into the functional and health impacts of CNVs. These findings represent valuable clues to genetic diversity and personalized medicine in the Tunisian population as well as in the ethnically similar populations from North Africa.

Berbers and Arabs: Tracing the genetic diversity and history of Southern Tunisia through genome wide analysis

OBJECTIVES

Tunisia has been a crossroads for people from Africa, Europe, and the Middle East since prehistoric times. At present, it is inhabited by two main ethnic groups, Arabs and Berbers, and several minorities. This study aims to advance knowledge regarding their genetic structure using new population samplings and a genome-wide approach.

MATERIALS AND METHODS

We investigated genomic variation, estimated ancestry components and dated admixture events in three Berber and two Arab populations from Southern Tunisia, mining a dataset including Middle Eastern, sub-Saharan, and European populations.

RESULTS

Differences in the proportion of North African, Arabian, and European ancestries and the varying impact of admixture and isolation determined significant heterogeneity in the genetic structure of Southern Tunisian populations. Admixture time estimates show a multilayer pattern of admixture events, involving both ethno-linguistic groups, which started around the mid XI century and lasted for nearly five centuries.

DISCUSSION

Our study provides evidence that the relationships between genetic and cultural diversity of old and new inhabitants of North Africa in southern Tunisia follow different patterns. The Berbers seem to have preserved a significant part of their common genomic heritage despite Islamization, Arab cultural influence, and linguistic diversity. Compared to Morocco and Algeria, southern Tunisian Arabs have retained a higher level of Arabian ancestry. This is more evident in the semi-nomad R'Baya, who have kept their original Bedouin lifestyle, than in the population from Douz, who have undergone multiple events of stratification and admixture.

Association of rs662799 variant and APOA5 gene haplotypes with metabolic syndrome and its components: a meta-analysis in North Africa

Apolipoprotein A5 (APOA5) has been linked to metabolic syndrome (MetS) in several populations. In North Africa, only the Tunisian and Moroccan populations were investigated. Our aim is to assess the association between APOA5 gene variant (rs662799) and haplotypes with MetS in Tunisian population and to perform a meta-analysis in North Africa. A total of 594 Tunisian participants were genotyped for polymorphism rs662799 using KASPar technology. Two polymorphisms rs3135506 and rs651821 in APOA5 gene genotyped in our previous study, were used in addition to rs662799 to assess the haplotype association with MetS. The genotype of 875 participants was used for the meta-analysis. Statistical analyses were performed with R software.

The rs662799 increases the risk of MetS under the dominant (P=0.018) and the additive models (P=0.028) in the Tunisian population. After stratification of the cohort following the sex and the geographic origin, a positive association of rs662799 with MetS was found for participant from the Northern region and for the women group. Only the haplotype AGT showed a significant association with MetS by decreasing the risk of the disease.

The meta-analysis reported a significant association of rs662799 and rs3135506 with MetS.

Our results showed a significant association between the APOA5 gene variants rs662799 and haplotypes with MetS and its traits in Tunisia. An impact of the sex and the geographic origin on the genotype distribution was highlighted. Our funding emphasizes the role of APOA5 in the development of MetS in North Africa.

PORTABILITY OF GWAS RESULTS BETWEEN ETHNIC POPULATIONS: GENETIC MARKERS FOR POLYCYSTIC OVARY SYNDROME (PCOS) IN MEDITERRANEAN AREA

Genome Wide Association Studies (GWAS) are excellent opportunities to define culprit genes in complex disorders such as the polycystic ovary syndrome (PCOS). PCOS is a prevalent disorder characterized by anovulation, hyperandrogenism and polycystic ovaries, which benefitted from several GWASs in Asians and Europeans revealing more than 20 potential culprit genes near associated single nucleotide variations (SNV). Translation of these findings into the clinical practice raises difficulties since positive hits are surrogate SNVs linked with causative mutations by linkage disequilibrium (LD). Studies in Mediterranean populations (e.g. Southern Europe and North Africa) raise supplementary problems because of a different LD-pattern, which may disrupt the link with causative mutations. Our experience in MEDIGENE program between Tunisia and France enforces the necessity of genetic anthropology studies before translating GWAS data. Tunisians are a heterogeneous population with ancestral Berbers, European, Arab and Sub-Saharan African components while South Europeans display a high level of genetic diversity, partially explained by gene flow from North Africa. Human diversity studies require sampling from Middle East and North Africa (MENA) region that will help to understand genetic factors in complex diseases.

A genome wide SNP genotyping study in the Tunisian population: specific reporting on a subset of common breast cancer risk loci

Background

Breast cancer is the most common cancer in women worldwide. Around 50% of breast cancer familial risk has been so far explained by known susceptibility alleles with variable levels of risk and prevalence. The vast majority of these breast cancer associated variations reported to date are from populations of European ancestry. In spite of its heterogeneity and genetic wealth, North-African populations have not been studied by the HapMap and the 1000Genomes projects. Thus, very little is known about the genetic architecture of these populations.

Methods

This study aimed to investigate a subset of common breast cancer loci in the general Tunisian population and to compare their genetic composition to those of other ethnic groups. We undertook a genome-wide haplotype study by genotyping 135 Tunisian subjects using the Affymetrix 6.0-Array. We compared Tunisian allele frequencies and linkage disequilibrium patterns to those of HapMap populations and we performed a comprehensive assessment of the functional effects of several selected variants.

Results

Haplotype analyses showed that at risk haplotypes on 2p24, 4q21, 6q25, 9q31, 10q26, 11p15, 11q13 and 14q32 loci are considerably frequent in the Tunisian population (> 20%). Allele frequency comparison showed that the frequency of rs13329835 is significantly different between Tunisian and all other HapMap populations. LD-blocks and Principle Component Analysis revealed that the genetic characteristics of breast cancer variants in the Tunisian, and so probably the North-African populations, are more similar to those of Europeans than Africans.

Using eQTl analysis, we characterized rs9911630 as the most strongly expression-associated SNP that seems to affect the expression levels of BRCA1 and two long non coding RNAs (NBR2 and LINC008854). Additional in-silico analysis also suggested a potential functional significance of this variant.

Conclusions

We illustrated the utility of combining haplotype analysis in diverse ethnic groups with functional analysis to explore breast cancer genetic architecture in Tunisia. Results presented in this study provide the first report on a large number of common breast cancer genetic polymorphisms in the Tunisian population which may establish a baseline database to guide future association studies in North Africa.

Electronic supplementary material

The online version of this article (10.1186/s12885-018-5133-8) contains supplementary material, which is available to authorized users.

Clinical and genetic study of Tunisian families with genetic generalized epilepsy: contribution of CACNA1H and MAST4 genes

Genetic generalized epilepsies (GGE) (childhood absence epilepsy (CAE), juvenile myoclonic epilepsy (JME) and epilepsy with generalized tonic-clonic seizures (GTCS)) are mainly determined by genetic factors. Since few mutations were identified in rare families with autosomal dominant GGE, a polygenic inheritance was suspected in most patients. Recent studies on large American or European cohorts of sporadic cases showed that susceptibility genes were numerous although their variants were rare, making their identification difficult. Here, we reported clinical and genetic characteristics of 30 Tunisian GGE families, including 71 GGE patients. The phenotype was close to that in sporadic cases. Nineteen pedigrees had a homogeneous type of GGE (JME-CAE-CGTS), and 11 combined these epileptic syndromes. Rare non-synonymous variants were selected in probands using a targeted panel of 30 candidate genes and their segregation was determined in families. Molecular studies incriminated different genes, mainly CACNA1H and MAST4. The segregation of at least two variants in different genes in some pedigrees was compatible with the hypothesis of an oligogenic inheritance, which was in accordance with the relatively low frequency of consanguineous probands. Since at least 2 susceptibility genes were likely shared by different populations, genetic factors involved in the majority of Tunisian GGE families remain to be discovered. Their identification should be easier in families with a homogeneous type of GGE, in which an intra-familial genetic homogeneity could be suspected.

Pharmacogenetic landscape of Metabolic Syndrome components drug response in Tunisia and comparison with worldwide populations

Genetic variation is an important determinant affecting either drug response or susceptibility to adverse drug reactions. Several studies have highlighted the importance of ethnicity in influencing drug response variability that should be considered during drug development. Our objective is to characterize the genetic variability of some pharmacogenes involved in the response to drugs used for the treatment of Metabolic Syndrome (MetS) in Tunisia and to compare our results to the worldwide populations. A set of 135 Tunisians was genotyped using the Affymetrix Chip 6.0 genotyping array. Variants located in 24 Very Important Pharmacogenes (VIP) involved in MetS drug response were extracted from the genotyping data. Analysis of variant distribution in Tunisian population compared to 20 worldwide populations publicly available was performed using R software packages. Common variants between Tunisians and the 20 investigated populations were extracted from genotyping data. Multidimensional screening showed that Tunisian population is clustered with North African and European populations. The greatest divergence was observed with the African and Asian population. In addition, we performed Inter-ethnic comparison based on the genotype frequencies of five VIP biomarkers. The genotype frequencies of the biomarkers rs3846662, rs1045642, rs7294 and rs12255372 located respectively in HMGCR, ABCB1, VKORC1 and TCF7L2 are similar between Tunisian, Tuscan (TSI) and European (CEU). The genotype frequency of the variant rs776746 located in CYP3A5 gene is similar between Tunisian and African populations and different from CEU and TSI. The present study shows that the genetic make up of the Tunisian population is relatively complex in regard to pharmacogenes and reflects previous historical events. It is important to consider this ethnic difference in drug prescription in order to optimize drug response to avoid serious adverse drug reactions. Taking into account similarities with other neighboring populations, our study has an impact not only on the Tunisian population but also on North African population which are underrepresented in pharmacogenomic studies.

Ancient and recent Middle Eastern maternal genetic contribution to North Africa as viewed by mtDNA diversity in Tunisian Arab populations

OBJECTIVES

Through previous mitochondrial DNA studies, the Middle Eastern maternal genetic contribution to Tunisian populations appears limited. In fact, most of the studied communities were cosmopolitan, or of Berber or Andalusian origin. To provide genetic evidence for the actual contribution of Middle Eastern mtDNA lineages to Tunisia, we focused on two Arab speaking populations from Kairouan and Wesletia known to belong to an Arab genealogical lineage.

MATERIALS AND METHODS

A total of 114 samples were sequenced for the mtDNA HVS-I and HVS-II regions. Using these data, we evaluated the distribution of Middle Eastern haplogroups in the study populations, constructed interpolation maps, and established phylogenetic networks allowing estimation of the coalescence time for three specific Middle Eastern subclades (R0a, J1b, and T1).

RESULTS

Both studied populations displayed North African genetic structure and Middle Eastern lineages with a frequency of 12% and 28.12% in Kairouan and Wesletia, respectively. TMRCA estimates for haplogroups T1a, R0a, and J1b in Tunisian Arabian samples were around 15 000 YBP, 9000 to 5000 YBP, and 960 to 600 YBP, respectively.

CONCLUSIONS

The Middle Eastern maternal genetic contribution to Tunisian populations, as to other North African populations, occurred mostly in deep prehistory. They were brought in different migration waves during the Upper Paleolithic, probably with the expansion of Iberomaurusian culture, and during Epipaleolithic and Early Neolithic periods, which are concomitant with the Capsian civilization. Middle Eastern lineages also came to Tunisia during the recent Islamic expansion of the 7th CE and the subsequent massive Bedouin migration during the 11th CE.

Lactase persistence in Tunisia as a result of admixture with other Mediterranean populations

Background

The ability to digest lactose after weaning, namely, lactase persistence (LP), is encoded by polymorphisms in the MCM6 gene and varies widely in frequency among different human populations. Although, evolution of LP-related genetic variants was investigated in many groups of Sub-Saharan African, Middle Eastern, and European ancestry, only few studies have focused on populations from North Africa and no data are especially available from the Tunisian one. For this reason, there is an urgent need to investigate the frequency patterns at these loci in Tunisia since this adaptive trait is implicated in health.

Methods

Forty SNPs covering the LCT/MCM6 genes and including the two functional variants − 13,910 C > T and − 22,018 G > A were genotyped in 117 Tunisian individuals using the Sequenom Mass Array technology. The observed nucleotide and haplotype patterns of variation were then compared with those of several African, European, and Mediterranean human groups for which comparable data were publicly available. Admixture analysis on a 5 Mb genomic region surrounding the LCT/MCM6 loci was also performed by extracting genotypes from a previously generated genome-wide dataset in order to deepen the reconstruction of the evolutionary history of these loci.

Results

We found that lactase non-persistence (LNP)-related alleles and haplotypes were predominantly present in the examined population. A clear differentiation between Tunisian, African, and North European/North Italian samples was found, while the Tunisian population showed more genetic affinity to Central and South Italian groups.

Conclusions

Our study provided a first report of LP-associated alleles and haplotypes in the Tunisian population. We highlighted a gradient followed by LP diffusion from Europe to North Africa. Based on the rich historic background of Tunisia, we suggest that this adaptive trait was introduced in that geographic region by a relatively recent gene flow.

Electronic supplementary material

The online version of this article (doi:10.1186/s12263-017-0573-3) contains supplementary material, which is available to authorized users.

Mitochondrial DNA analysis of Tunisians reveals a mosaic genetic structure with recent population expansion

Tunisia is a country of great interest for human population genetics due to its strategic geographic position and rich human settlement history. These factors significantly contributed to the genetic makeup of present-day Tunisians harbouring components of diverse geographic origins. Here, we investigated the genetic structure of Tunisians by performing a mitochondrial DNA (mtDNA) comparison of 15 Tunisian population groups, in order to explore their complex genetic landscape. All Tunisian data were also analysed against 40 worldwide populations. Statistical results (Tajima's D and Fu's F_S tests) suggested recent population expansion for the majority of studied populations, as well as showed (AMOVA test) that all populations were significantly different from each other, which is evidence of population structure even if it is not guided by geographic and ethnic effects. Gene flow analysis revealed the assignment of Tunisians to multiple ancestries, which agrees with their genetic heterogeneity. The resulting picture for the mtDNA pool confirms the evidence of a recent expansion of the Tunisian population and is in accordance with a mosaic structure, composed by North African, Middle Easterner, European and Sub-Saharan lineages, resulting from a complex settlement history.

On the origin of Iberomaurusians: new data based on ancient mitochondrial DNA and phylogenetic analysis of Afalou and Taforalt populations

The Western North African population was characterized by the presence of Iberomaurusian civilization at the Epiplaeolithic period (around 20,000 years before present (YBP) to 10,000 YBP). The origin of this population is still not clear: they may come from Europe, Near East, sub-Saharan Africa or they could have evolved in situ in North Africa. With the aim to contribute to a better knowledge of the settlement of North Africa we analysed the mitochondrial DNA extracted from Iberomaurusian skeletons exhumed from the archaeological site of Afalou (AFA) (15,000-11,000 YBP) in Algeria and from the archaeological site of Taforalt (TAF) (23,000-10,800 YBP) in Morocco. Then, we carried out a phylogenetic analysis relating these Iberomaurusians to 61 current Mediterranean populations. The genetic structure of TAF and AFA specimens contains only North African and Eurasian maternal lineages. These finding demonstrate the presence of these haplotypes in North Africa from at least 20,000 YBP. The very low contribution of a Sub-Saharan African haplotype in the Iberomaurusian samples is confirmed. We also highlighted the existence of genetic flows between Southern and Northern coast of the Mediterranean.

The first complete mitochondrial genome sequence of Uraspis secunda (Perciformes: Carangidae) and its phylogenetic relationship

In this study, the first complete mitochondrial genome sequence of Uraspis secunda was determined by using Illumina and Sanger sequencing technology. This genome was 16 554 bp long, consisting of 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and a putative control region. Twenty-eight genes were located on heavy strand, while nine genes were located on light strand. Totally two types of start codons and four kinds of stop codons were observed from 13 protein-coding genes. The overall base composition of this genome was 28.17% for A, 16.23% for G, 25.81% for T, and 29.78% for C, with a slightly higher A + T content of 53.99%. The highest A + T content (63.18%) was found in the control region of the mitochondrial genome. A phylogenetic tree was constructed based on 13 concatenated protein-coding genes, which indicated that U. secunda was closest to Parastromateus niger, and these two fish species clustered a 1. monophyletic group with Alectis ciliaris, Carangoides malabaricus and C. armatus. This study will be helpful for researches on molecular phylogeny and genetic diversity of U. secunda and the related fish species.

Characterization of the complete mitochondrial genome and phylogenetic relationship of Caranx tille (Perciformes: Carangidae)

In the present study, we determined the complete mitochondrial genome sequence of Caranx tille by using Illumina and Sanger sequencing technology that was 16 593 bp in size. This mitochondrial genome was made up of 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and a putative control region. Of these 37 genes, 28 were encoded by heavy strand, while nine were encoded by light strand. Three types of start codons and four kinds of stop codons were separately detected in 13 protein-coding genes. The overall base composition of this genome was 28.95% for A, 15.84% for G, 26.26% for T, and 28.95% for C, with a slight higher A + T content of 55.21%. From the phylogenetic tree constructed based on 13 concatenated protein-coding genes, we can see that C. tille, C. melampygus and C. ignobilis grouped together and formed a single cluster. In addition, C. tille was genetically closest to C. melampygus. This work should facilitate the studies on molecular evolution and population genetic structure for fishes under family Carangidae.

Study of the T16189C variant and mitochondrial lineages in Tunisian and overall Mediterranean region

The mitochondrial DNA (mtDNA) variant T16189C has been investigated in several metabolic diseases. In this study, we aimed to estimate the frequency of the T16189C variant in Tunisian and other Mediterranean populations and to evaluate the impact of this variant on the phylogeny of Mediterranean populations. Blood sample of 240 unrelated Tunisian subjects were recruited from several Tunisian localities. The hypervariable region 1 of the mtDNA were amplified and sequenced. Additional sequences (N = 4921) from Mediterranean populations were compiled from previous studies. The average frequency of T16189C variant in Tunisia (29%) is similar to that observed in North African and Near Eastern populations. Our findings showed positive correlation of the T16189C variant with Sub-Saharan and North African lineages, while a negative correlation was found with the Eurasian haplogroups, reaching its maximum with the Eurasian haplogroup H. The principal component analyses showed a high internal heterogeneity between Tunisian localities. At the Mediterranean scale, Tunisians are closer to North African (Algerian and Moroccan) and Near Eastern populations (Syrians and Palestinians) than to Europeans.