Search

Background In Tunisia a first SARS-CoV-2 confirmed case was reported in March 03, 2020. Since then, an increase of cases number was observed from either imported or local cases. The aim of this preliminary study was to better understand the molecular epidemiology and genetic variability of SARS-CoV-2 viruses circulating in Tunisia and worldwide. Methods Whole genome sequencing was performed using NGS approach on six SARS. CoV-2 highly positive samples detected during the early phase of the outbreak. Results Full genomes sequences of six Tunisian SARS-CoV-2 strains were obtained from imported and locally transmission cases during the COVID-19 outbreak. Reported sequences were non-identical with 0.1% nucleotide divergence rate and clustered into 6 different clades with worldwide sequences. SNPs results favor the distribution of the reported Tunisian sequences into 3 major genotypes. These SNP mutations are critical for diagnosis and vaccine development. Conclusions These results indicate multiple introductions of the virus in Tunisia and add new genomic data on SARS-CoV-2 at the international level.

Background Coronavirus Disease 2019 (COVID-19) is a viral pandemic disease that may induce severe pneumonia in humans. In this paper, we investigated the putative implication of 12 vaccines, including BCG, OPV and MMR in the protection against COVID-19. Sequences of the main antigenic proteins in the investigated vaccines and SARS-CoV-2 proteins were compared to identify similar patterns. The immunogenic effect of identified segments was, then, assessed using a combination of structural and antigenicity prediction tools. Results A total of 14 highly similar segments were identified in the investigated vaccines. Structural and antigenicity prediction analysis showed that, among the identified patterns, three segments in Hepatitis B, Tetanus, and Measles proteins presented antigenic properties that can induce putative protective effect against COVID-19. Conclusions Our results suggest a possible protective effect of HBV, Tetanus and Measles vaccines against COVID-19, which may explain the variation of the disease severity among regions.

SARS-CoV-2 is constantly evolving with lineages emerging and others eclipsing. Some lineages have an important epidemiological impact and are known as variants of interest (VOIs), variants under monitoring (VUMs) or variants of concern (VOCs). Lineage A.27 was first defined as a VUM since it holds mutations of concern. Here, we report additional lineage A.27 data and sequences from five African countries and describe the molecular characteristics, and the genetic history of this lineage worldwide. Based on the new sequences investigated, the most recent ancestor (tMRCA) of lineage A.27 was estimated to be from April 2020 from Niger. It then spread to Europe and other parts of the world with a peak observed between February and April 2021. The detection rate of A.27 then decreased with only a few cases reported during summer 2021. The phylogenetic analysis revealed many sub-lineages. Among them, one was defined by the substitution Q677H in the spike (S) gene, one was defined by the substitution D358N in the nucleoprotein (N) gene and one was defined by the substitution A2143V in the ORF1b gene. This work highlights the importance of molecular characterization and the timely submission of sequences to correctly describe the circulation of particular strains in order to be proactive in monitoring the pandemic.

Documenting the circulation dynamics of SARS-CoV-2 variants in different regions of the world is crucial for monitoring virus transmission worldwide and contributing to global efforts towards combating the pandemic. Tunisia has experienced several waves of COVID-19 with a significant number of infections and deaths. The present study provides genetic information on the different lineages of SARS-CoV-2 that circulated in Tunisia over 17 months. Lineages were assigned for 1359 samples using whole-genome sequencing, partial S gene sequencing and variant-specific real-time RT-PCR tests. Forty-eight different lineages of SARS-CoV-2 were identified, including variants of concern (VOCs), variants of interest (VOIs) and variants under monitoring (VUMs), particularly Alpha, Beta, Delta, A.27, Zeta and Eta. The first wave, limited to imported and import-related cases, was characterized by a small number of positive samples and lineages. During the second wave, a large number of lineages were detected; the third wave was marked by the predominance of the Alpha VOC, and the fourth wave was characterized by the predominance of the Delta VOC. This study adds new genomic data to the global context of COVID-19, particularly from the North African region, and highlights the importance of the timely molecular characterization of circulating strains.

International audience; We sought to determine the relative value of conventional molecular methods and whole-genome sequencing (WGS) for subtyping Salmonella enterica serovar Enteritidis recovered from 2000 to 2015 in Tunisia and to investigate the genetic diversity of this serotype. A total of 175 Salmonella Enteritidis isolates were recovered from human, animal, and foodborne outbreak samples. Pulsed-field gel electrophoresis (PFGE), multiple locus variable-number tandem repeat analysis (MLVA), and wholegenome sequencing were performed. Eight pulsotypes were detected for all isolates with PFGE (DI = 0.518). Forty-five Salmonella Enteritidis isolates were selected for the MLVA and WGS techniques. Eighteen MLVA profiles were identified and classified into two major clusters (DI = 0.889). Core genome multilocus typing (cgMLST) analysis revealed 16 profiles (DI = 0.785). Whole-genome analysis indicated 660 single-nucleotide polymorphism (SNP) divergences dividing these isolates into 43 haplotypes (DI = 0.997). The phylogenetic tree supported the classification of Salmonella Enteritidis isolates into two distinct lineages subdivided into five clades and seven subclades. Pairwise SNP differences between the isolates ranged between 302 and 350. We observed about 311 SNP differences between the two foodborne outbreaks, while only less or equal to 4 SNP differences within each outbreak. SNP-based WGS typing showed an excellent discriminatory power comparing with the conventional methods such as PFGE and MLVA. Besides, we demonstrate the added value of WGS as a complementary subtyping method to discriminate outbreak from non-outbreak isolates belonging to common subtypes. It is important to continue the survey of Salmonella Enteritidis lineages in Tunisia using WGS.

Effective interventions and treatments for complex diseases have been implemented globally, however, coverage in Africa has been comparatively lower due to lack of capacity, clinical applicability and knowledge on the genetic contribution to disease and treatment. Currently, there is a scarcity of genetic data on African populations, which have enormous genetic diversity. Pharmacogenomics studies have the potential to revolutionise treatment of diseases, therefore, African populations are likely to benefit from these approaches to identify likely responders, reduce adverse side effects and optimise drug dosing. This review discusses clinical pharmacogenetics studies conducted in African populations, focusing on studies that examined drug response in complex diseases relevant to healthcare. Several pharmacogenetics associations have emerged from African studies, as have gaps in knowledge.

International audience; Over the past 25 years, the powerful combination of genome sequencing and bioinformatics analysis has played a crucial role in interpreting information encoded in bacterial genomes. High-throughput sequencing technologies have paved the way towards understanding an increasingly wide range of biological questions. This revolution has enabled advances in areas ranging from genome composition to how proteins interact with nucleic acids. This has created unprecedented opportunities through the integration of genomic data into clinics for the diagnosis of genetic traits associated with disease. Since then, these technologies have continued to evolve, and recently, long-read sequencing has overcome previous limitations in terms of accuracy, thus expanding its applications in genomics, transcriptomics and metagenomics. In this review, we describe a brief history of the bacterial genome sequencing revolution and its application in public health and molecular epidemiology. We present a chronology that encompasses the various technological developments: whole-genome shotgun sequencing, high-throughput sequencing, long-read sequencing. We mainly discuss the application of next-generation sequencing to decipher bacterial genomes. Secondly, we highlight how long-read sequencing technologies go beyond the limitations of traditional short-read sequencing. We intend to provide a description of the guiding principles of the 3rd generation sequencing applications and ongoing improvements in the field of microbial medical research.

H9N2 avian influenza virus (AIV) has been isolated from various species of wild birds and domestic poultry worldwide. It has been reported since the late 1990s, that H9N2 AIV has infected humans as reported in some Asian and North African countries. This subtype has already been circulating and constituting a serious threat to the poultry industry in Tunisia back in 2009. To investigate zoonotic potential and pathogenicity of H9N2 AIV in chickens and mice in Tunisia, five strains have been isolated during the period from 2014 to 2018. Samples were withdrawn from several wild bird species and environment (Lagoon water) of Maamoura and Korba Lagoons as well as Kuriat Island. Phylogenetic analyzes demonstrated that the isolated H9N2 strains belonged to the G1-like sublineage and were close to AIV H9N2 poultry viruses from North Africa, West Africa and the Middle East. All strains carried in their hemagglutinin the residue 226 L, which is an important marker for avian-to-human viral transmission. The hemagglutinin cleavage site has several motifs: PSKSSR/G, PARSSR/G and HARSSR/G. The neuraminidase showed S372A and R403W substitutions that have been previously detected in H3N2 and H2N2 viruses that were reported in human pandemics. Many mutations associated with mammalian infections have been detected in internal proteins. Pathogenicity evaluation in chickens showed that GF/14 replicates effectively in the lungs, tracheas, spleens, kidneys and brains and that it was transmitted among contact chickens. However, GHG/18 replicates poorly in chickens and has not an efficient transmission in contact chickens. GF/14 and GHG/18 could not kill mice though they replicated in their respiratory tract and caused a significant body weight loss (p 0.05). This study highlights the importance of H9N2 AIV monitoring in both migratory birds and the environment to prevent virus transmission to humans.

African genomic medicine and microbiome datasets are usually not well characterized in terms of their origin, making it difficult to find and extract data for specific African ethnic groups or even countries. The Pan-African H3Africa Bioinformatics Network (H3ABioNet) recognized the need for developing data portals for African genomic medicine and African microbiomes to address this and ran a hackathon to initiate their development. The two portals were designed and significant progress was made in their development during the hackathon. All the participants worked in a very synergistic and collaborative atmosphere in order to achieve the hackathon's goals. The participants were divided into content and technical teams and worked over a period of 6 days. In response to one of the survey questions of what the participants liked the most during the hackathon, 55% of the hackathon participants highlighted the familial and friendly atmosphere, the team work and the diversity of team members and their expertise. This paper describes the preparations for the portals hackathon and the interaction between the participants and reflects upon the lessons learned about its impact on successfully developing the two data portals as well as building scientific expertise of younger African researchers. Database URL: The code for developing the two portals was made publicly available in GitHub repositories: [https://github.com/codemeleon/Database; https://github.com/codemeleon/AfricanMicrobiomePortal].

Background Hepatitis C virus (HCV) has a high genetic diversity. Eight genotypes and 90 subtypes are currently described. Genotypes are clinically significant for therapeutic management and their determination is necessary for epidemiological studies. Methods Tunisian patients plasma samples (n = 6) with unassigned HCV-2 subtype using partial sequencing in the NS5B and Core/E1 regions were analyzed by realizing whole-genome sequencing analysis. Phylogenetic analyses were performed to assign subtypes. Results Phylogenetic analysis of the full genome sequences of Tunisian strains shows two subtypes within HCV-2. These later were genetically distinct from all previously established HCV-2 subtypes with nucleotide divergence greater than 15% (20% -31%). These two subtypes are proposed as new subtypes 2v and 2w. Conclusions The discovery of two new HCV-2 subtypes circulating in the Tunisian population confirms the great diversity of HCV-2 viruses and increases the total number of HCV-2 subtypes from 21 to 23.

Orphan diseases (ODs) are progressive genetic disorders, which affect a small number of people. The principal fundamental aspects related to these diseases include insufficient knowledge of mechanisms involved in the physiopathology necessary to access correct diagnosis and to develop appropriate healthcare. Unlike ODs, complex diseases (CDs) have been widely studied due to their high incidence and prevalence allowing to understand the underlying mechanisms controlling their physiopathology. Few studies have focused on the relationship between ODs and CDs to identify potential shared pathways and related molecular mechanisms which would allow improving disease diagnosis, prognosis, and treatment. We have performed a computational approach to studying CDs and ODs relationships through (1) connecting diseases to genes based on genes-diseases associations from public databases, (2) connecting ODs and CDs through binary associations based on common associated genes, and (3) linking ODs and CDs to common enriched pathways. Among the most shared significant pathways between ODs and CDs, we found pathways in cancer, p53 signaling, mismatch repair, mTOR signaling, B cell receptor signaling, and apoptosis pathways. Our findings represent a reliable resource that will contribute to identify the relationships between drugs and disease-pathway networks, enabling to optimise patient diagnosis and disease treatment.

Aims Characterization of quinolone-resistant Salmonella Kentucky and Typhimurium isolates in Tunisia from various sources, detection of some plasmid-mediated quinolone resistance genes and the genetic relatedness. Methods A total of 1404 isolates of S. Kentucky (n = 1059)/S. Typhimurium (n = 345) from various sources from all over Tunisia were tested for quinolone resistance by disk diffusion method. Minimum inhibitory concentrations of nalidixic acid, ciprofloxacin and ofloxacin were determined. Quinolone-resistant isolates were screened for plasmid-mediated quinolone-resistance genes (qnrA,qnrB,qnrS, aac(6')-Ib-cr and qepA) by polymerase chain reaction (PCR). Mutations in the quinolone-resistance-determining regions of the gyrA and parC genes were detected by PCR and DNA sequencing. Pulsed-field gel electrophoresis and multilocus sequence typing were accomplished for isolates harbouring plasmid-mediated quinolone-resistance genes. Results According to our selection criteria (NAL = resistance phenotype; CIP = resistant with diameter 0, or intermediate), only 63 S. Kentucky/41 S. Typhimurium isolates were investigated: 49% (5/104) were multidrug resistant. Two S. Typhimurium isolates harboured qnrB19 with different PFGE profiles. A mutation was detected in the gyrA gene for each of these two isolates. MLST revealed the presence of ST313 and ST34, an endemic sequence type. Conclusion Our study highlights the presence of quinolone multidrug-resistant Salmonella in humans and animals in Tunisia. This is the first report of S. Typhimurium ST34 in Africa and qnrB19 in Tunisia. Significance and impact of the study This is the first report that describes not only the current epidemiological situation of the quinolone resistance in S. Kentucky and Typhimurium isolated from various sources and regions in Tunisia, but also, the genetic resistance determinants associated with phenotypic antibiotic resistance and the molecular mechanisms of their quinolone-resistance. Also, we provide the first report of S. Typhimurium ST34 in Africa, and the first report of qnrB19 in Salmonella in Tunisia.

The spread of COVID-19 caused by the SARS-CoV-2 outbreak has been growing since its first identification in December 2019. The publishing of the first SARS-CoV-2 genome made a valuable source of data to study the details about its phylogeny, evolution, and interaction with the host. Protein-protein binding assays have confirmed that Angiotensin-converting enzyme 2 (ACE2) is more likely to be the cell receptor through which the virus invades the host cell. In the present work, we provide an insight into the interaction of the viral spike Receptor Binding Domain (RBD) from different coronavirus isolates with host ACE2 protein. By calculating the binding energy score between RBD and ACE2, we highlighted the putative jump in the affinity from a progenitor form of SARS-CoV-2 to the current virus responsible for COVID-19 outbreak. Our result was consistent with previously reported phylogenetic analysis and corroborates the opinion that the interface segment of the spike protein RBD might be acquired by SARS-CoV-2 via a complex evolutionary process rather than a progressive accumulation of mutations. We also highlighted the relevance of Q493 and P499 amino acid residues of SARS-CoV-2 RBD for binding to human ACE2 and maintaining the stability of the interface. Moreover, we show from the structural analysis that it is unlikely for the interface residues to be the result of genetic engineering. Finally, we studied the impact of eight different variants located at the interaction surface of ACE2, on the complex formation with SARS-CoV-2 RBD. We found that none of them is likely to disrupt the interaction with the viral RBD of SARS-CoV-2., Highlights • We noticed a large difference in the binding energy of RBD from the spike protein of SARS-CoV-2 with ACE2, compared to the closest isolates to the proposed progenitor forms. • Currently, available sequence data do not explain the transition from low-affinity RBD forms to high-affinity RBD that are capable to infect the Humans. • We suggest that residues N493Q and not N501T as was reported previously is responsible for a higher affinity toward ACE2 when we compare SARS-CoV-2 to SARS-CoV. • The mutation T499P is responsible for stabilizing the interface of RBD interacting with ACE2 when we compare SARS-CoV-2 to SARS-CoV. • It is unlikely that the interface of RBD from SARS-CoV-2 is a result of human intervention aiming to increase the affinity toward ACE2. • The role of receptors other than ACE2 is shown, and these may have a more critical function in crossing the species barrier. • Variants that corresponding to residues on the ACE2 interface are unlikely to be associated with resistance or sensibility forms compared to the reference allele.

International audience; Molecular diagnosis of rare inherited palmoplantar keratoderma (PPK) is still challenging. We investigated at the clinical and genetic level a consanguineous Tunisian family presenting an autosomal dominant atypical form of transgrediens and progrediens PPK to better characterize this ultrarare disease and to identify its molecular etiology. Whole-exome sequencing (WES), filtering strategies, and bioinformatics analysis have been achieved. Clinical investigation and follow up over 13 years of this Tunisian family with three siblings formerly diagnosed as an autosomal recessive form of Mal de Melela-like conducted us to reconsider its initial phenotype. Indeed, the three patients presented clinical features that overlap both Mal de Meleda and progressive symmetric erythrokeratoderma (PSEK). The mode of inheritance was also reconsidered, since the mother, initially classified as unaffected, exhibited a similar expression of the disease. WES analysis showed the absence of potentially functional rare variants in known PPKs or PSEK-related genes. Results revealed a novel heterozygous nonsynonymous variant in cadherin-12 gene (CDH12, NM_004061, c.1655C > A, p.Thr552Asn) in all affected family members. This variant is absent in dbSNP and in 50 in-house control exomes. In addition, in silico analysis of the mutated 3D domain structure predicted that this variant would result in cadherin-12 protein destabilization and thermal instability. Functional annotation and biological network construction data provide further supporting evidence for the potential role of CDH12 in the maintenance of skin integrity. Taken together, these results suggest that CDH12 gene is a potential candidate gene for an atypical presentation of an autosomal dominant form of transgrediens and progrediens PPK.

Host-pathogen molecular cross-talks are critical in determining the pathophysiology of a specific infection. Most of these cross-talks are mediated via protein-protein interactions between the host and the pathogen (HP-PPI). Thus, it is essential to know how some pathogens interact with their hosts to understand the mechanism of infections. Malaria is a life-threatening disease caused by an obligate intracellular parasite belonging to the Plasmodium genus, of which P. falciparum is the most prevalent. Several previous studies predicted human-plasmodium protein-protein interactions using computational methods have demonstrated their utility, accuracy, and efficiency to identify the interacting partners and therefore complementing experimental efforts to characterize host-pathogen interaction networks. To predict potential putative HP-PPIs, we use an integrative computational approach based on the combination of multiple OMICS-based methods including human red blood cells (RBC) and Plasmodium falciparum 3D7 strain expressed proteins, domain-domain based PPI, similarity of gene ontology terms, structure similarity method homology identification, and machine learning prediction. Our results reported a set of 716 protein interactions involving 302 human proteins and 130 Plasmodium proteins. This work provides a list of potential human-Plasmodium interacting proteins. These findings will contribute to better understand the mechanisms underlying the molecular determinism of malaria disease and potentially to identify candidate pharmacological targets.

Rhinoviruses (RV) are a major cause of Severe Acute Respiratory Infection (SARI) in children, with high genotypic diversity in different regions. However, RV type diversity remains unknown in several regions of the world. In this study, the genetic variability of the frequently circulating RV types in Northern Tunisia was investigated, using phylogenetic and phylogeographic analyses with a specific focus on the most frequent RV types: RV-A101 and RV-C45. This study concerned 13 RV types frequently circulating in Northern Tunisia. They were obtained from respiratory samples collected in 271 pediatric SARI cases, between September 2015 and November 2017. A total of 37 RV VP4-VP2 sequences, selected among a total of 49 generated sequences, was compared to 359 sequences from different regions of the world. Evolutionary analysis of RV-A101 and RV-C45 showed high genetic relationship between different Tunisian strains and Malaysian strains. RV-A101 and C45 progenitor viruses’ dates were estimated in 1981 and 1995, respectively. Since the early 2000s, the two types had a wide spread throughout the world. Phylogenetic analyses of other frequently circulating strains showed significant homology of Tunisian strains from the same epidemic period, in contrast with earlier strains. The genetic relatedness of RV-A101 and RV-C45 might result from an introduction of viruses from different clades followed by local dissemination rather than a local persistence of an endemic clades along seasons. International traffic may play a key role in the spread of RV-A101, RV-C45, and other RVs.

Computational approaches have been developed to prioritize candidate genes in disease gene identification. They are based on different pieces of evidences associating each gene with the given disease. In this study, 648 genes underlying genodermatoses have been compared to 1808 genes involved in other genetic diseases using a bioinformatic approach. These genes were studied at the structural, evolutionary and functional levels. Results show that genes underlying genodermatoses present longer CDS and have more exons. Significant differences were observed in nucleotide motif and amino-acid compositions. Evolutionary conservation analysis revealed that genodermatoses genes have less paralogs, more orthologs in Mouse and Dog and are less conserved. Functional analysis revealed that genodermatosis genes seem to be involved in immune system and skin layers. The Bayesian network model returned a rate of good classification of around 80%. This computational approach could help investigators working in the field of dermatology by prioritizing positional candidate genes for mutation screening.

Sickle cell disease (SCD) is one of the most common monogenic diseases in humans with multiple phenotypic expressions that can manifest as both acute and chronic complications. Although described more than a century ago, challenges in comprehensive disease management and collaborative research on this disease are compounded by the complex molecular and clinical phenotypes of SCD, environmental and psychosocial factors, limited therapeutic options and ambiguous terminology. This ambiguous terminology has hampered the integration and interoperability of existing SCD knowledge, and SCD research translation. The SCD Ontology (SCDO), which is a community-driven integrative and universal knowledge representation system for SCD, overcomes this issue by providing a controlled vocabulary developed by a group of experts in both SCD and ontology design. SCDO is the first and most comprehensive standardized human- and machine-readable resource that unambiguously represents terminology and concepts about SCD for researchers, patients and clinicians. It is built around the central concept ‘hemoglobinopathy’, allowing inclusion of non-SCD haemoglobinopathies, such as thalassaemias, which may interfere with or influence SCD phenotypic manifestations. This collaboratively developed ontology constitutes a comprehensive knowledge management system and standardized terminology of various SCD-related factors. The SCDO will promote interoperability of different research datasets, facilitate seamless data sharing and collaborations, including meta-analyses within the SCD community, and support the development and curation of data-basing and clinical informatics in SCD.

Hereditary breast cancer accounts for 5–10% of all breast cancer cases. So far, known genetic risk factors account for only 50% of the breast cancer genetic component and almost a quarter of hereditary cases are carriers of pathogenic mutations in BRCA1/2 genes. Hence, the genetic basis for a significant fraction of familial cases remains unsolved. This missing heritability may be explained in part by Copy Number Variations (CNVs). We herein aimed to evaluate the contribution of CNVs to hereditary breast cancer in Tunisia. Whole exome sequencing was performed for 9 BRCA negative cases with a strong family history of breast cancer and 10 matched controls. CNVs were called using the ExomeDepth R-package and investigated by pathway analysis and web-based bioinformatic tools. Overall, 483 CNVs have been identified in breast cancer patients. Rare CNVs affecting cancer genes were detected, of special interest were those disrupting APC2, POU5F1, DOCK8, KANSL1, TMTC3 and the mismatch repair gene PMS2. In addition, common CNVs known to be associated with breast cancer risk have also been identified including CNVs on APOBECA/B, UGT2B17 and GSTT1 genes. Whereas those disrupting SULT1A1 and UGT2B15 seem to correlate with good clinical response to tamoxifen. Our study revealed new insights regarding CNVs and breast cancer risk in the Tunisian population. These findings suggest that rare and common CNVs may contribute to disease susceptibility. Those affecting mismatch repair genes are of interest and require additional attention since it may help to select candidates for immunotherapy leading to better outcomes.

Typhimurium is one of the main Salmonella serovar responsible for non-typhoidal gastro-enteritis in Tunisia. Here, we aimed to assess the genetic diversity of 88 clinical Salmonella Typhimurium strains recovered during 14 years from 2000 to 2013. Phage typing, CRISPR polymorphisms (CRISPOL), pulsed-field gel electrophoresis (PFGE), multi-locus variable-number tandem repeat analysis (MLVA) and Whole genome sequencing (WGS) were used to study the relatedness and spatio-temporal evolution of Salmonella Typhimurium populations (Typhimurium (n = 81), monophasic (n = 3) and nonmotile (n = 4) variants). Seven-locus MLST from whole genome assemblies showed that all isolates, except one, belonged to ST19. The isolates were divided into 10 definitive phage (DT) types, dominated by DT104-L (39.8%), DT41 (14.8%), DT116 (11.4%) and DT120 (5.7%). Fifty-seven MLVA patterns (DI, 0.978) were obtained compared to 11 different CRISPOL types and 15 PFGE types (DI,0.845). For cgMLST analysis, 20 profiles were found. A total of 3056 SNPs were identified from the whole genome of the 88 Salmonella Typhimurium isolates. These SNPs resolved these isolates into 86 SNP haplotypes. The phylogeny result allocated most Salmonella Typhimurium isolates into four distinct clades and seven subclades. Genetic diversity between the four clades ranged in the order of 249 to 720 nucleotide changes. The prevalent phage type DT104L formed a major clade on the phylogenetic tree. Pairwise SNP differences between the strains of this clade ranged between 0 and 59. SNP-based WGS typing seems to be the most valuable molecular markers for studying the evolutionary relationships of homogeneous serovar Typhimurium isolates.

International audience; A chronic inflammatory state to a large extent explains sickle cell disease (SCD) pathophysi-ology. Nonetheless, the principal dysregulated factors affecting this major pathway and their mechanisms of action still have to be fully identified and elucidated. Integrating gene expression and genome-wide association study (GWAS) data analysis represents a novel approach to refining the identification of key mediators and functions in complex diseases. Here, we performed gene expression meta-analysis of five independent publicly available microarray datasets related to homozygous SS patients with SCD to identify a consensus SCD transcriptomic profile. The meta-analysis conducted using the MetaDE R package based on combining p values (maxP approach) identified 335 differentially expressed genes (DEGs; 224 upregulated and 111 downregulated). Functional gene set enrichment revealed the importance of several metabolic pathways, of innate immune responses, erythrocyte development, and hemostasis pathways. Advanced analyses of GWAS data generated within the framework of this study by means of the atSNP R package and SIFT tool identified 60 regulatory single-nucleotide polymorphisms (rSNPs) occurring in the promoter of 20 DEGs and a deleterious SNP, affecting CAMKK2 protein function. This novel database of candidate genes, transcription factors, and rSNPs associated with SCD provides new markers that may help to identify new therapeutic targets.

Background A family history of breast cancer has long been thought to indicate the presence of inherited genetic events that predispose to this disease. In North Africa, many specific epidemio-genetic characteristics have been observed in breast cancer families when compared to Western populations. Despite these specificities, the majority of breast cancer genetics studies performed in North Africa remain restricted to the investigation of the BRCA1 and BRCA2 genes. Thus, comprehensive data at a whole exome or whole genome level from local patients are lacking. Methods A whole exome sequencing (WES) of seven breast cancer Tunisian families have been performed using a family-based approach. We focused our analysis on BC-TN-F001 family that included two affected members that have been sequenced using WES. Relevant variants identified in BC-TN-F001 have been confirmed using Sanger sequencing. Then, we conducted an integrative analysis by combining our results with those from other WES studies in order to figure out the genetic transmission model of the newly identified genes. Biological network construction and protein–protein interactions analyses have been performed to decipher the molecular mechanisms likely accounting for the role of these genes in breast cancer risk. Results Sequencing, filtering strategies, and validation analysis have been achieved. For BC-TN-F001, no deleterious mutations have been identified on known breast cancer genes. However, 373 heterozygous, exonic and rare variants have been identified on other candidate genes. After applying several filters, 12 relevant high-risk variants have been selected. Our results showed that these variants seem to be inherited in a family specific model. This hypothesis has been confirmed following a thorough analysis of the reported WES studies. Enriched biological process and protein–protein interaction networks resulted in the identification of four novel breast cancer candidate genes namely MMS19, DNAH3, POLK and KATB6. Conclusions In this first WES application on Tunisian breast cancer patients, we highlighted the impact of next generation sequencing technologies in the identification of novel breast cancer candidate genes which may bring new insights into the biological mechanisms of breast carcinogenesis. Our findings showed that the breast cancer predisposition in non-BRCA families may be ethnic and/or family specific. Electronic supplementary material The online version of this article (10.1186/s12967-018-1504-9) contains supplementary material, which is available to authorized users.

International audience; The application of genomics technologies to medicine and biomedical research is increasing in popularity, made possible by new high-throughput genotyping and sequencing technologies and improved data analysis capabilities. Some of the greatest genetic diversity among humans, animals, plants, and microbiota occurs in Africa, yet genomic research outputs from the continent are limited. The Human Heredity and Health in Africa (H3Africa) initiative was established to drive the development of genomic research for human health in Africa, and through recognition of the critical role of bioinformatics in this process, spurred the establishment of H3ABioNet, a pan-African bioinformatics network for H3Africa. The limitations in bioinformatics capacity on the continent have been a major contributory factor to the lack of notable outputs in high-throughput biology research. Although pockets of high-quality bioinformatics teams have existed previously, the majority of research institutions lack experienced faculty who can train and supervise bioinformatics students. H3ABioNet aims to address this dire need, specifically in the area of human genetics and genomics, but knock-on effects are ensuring this extends to other areas of bioinformatics. Here, we describe the emergence of genomics research and the development of bioinformatics in Africa through H3ABioNet.

HCV genotype 2 (HCV-2) has a worldwide distribution with prevalence rates that vary from country to country. High genetic diversity and long-term endemicity were suggested in West African countries. A global dispersal of HCV-2 would have occurred during the 20th century, especially in European countries. In Tunisia, genotype 2 was the second prevalent genotype after genotype 1 and most isolates belong to subtypes 2c and 2k. In this study, phylogenetic analyses based on the NS5B genomic sequences of 113 Tunisian HCV isolates from subtypes 2c and 2k were carried out. A Bayesian coalescent-based framework was used to estimate the origin and the spread of these subtypes circulating in Tunisia. Phylogenetic analyses of HCV-2c sequences suggest the absence of country-specific or time-specific variants. In contrast, the phylogenetic grouping of HCV-2k sequences shows the existence of two major genetic clusters that may represent two distinct circulating variants. Coalescent analysis indicated a most recent common ancestor (tMRCA) of Tunisian HCV-2c around 1886 (1869-1902) before the introduction of HCV-2k in 1901 (1867-1931). Our findings suggest that the introduction of HCV-2c in Tunisia is possibly a result of population movements between Tunisia and European population following the French colonization.

The present study describes the in silico prediction of the regulatory network of Leishmania infected human macrophages. The construction of the gene regulatory network requires the identification of Transcription Factor Binding Sites (TFBSs) in the regulatory regions (promoters, enhancers) of genes that are regulated upon Leishmania infection. The promoters of human, mouse, rat, dog and chimpanzee genes were first identified in the whole genomes using available experimental data on full length cDNA sequences or deep CAGE tag data (DBTSS, FANTOM3, FANTOM4), mRNA models (ENSEMBL), or using hand annotated data (EPD, TRANSFAC). A phylogenetic footprinting analysis and a MATCH analysis of the promoter sequences were then performed to predict TFBS. Then, an SQL database that integrates all results of promoter analysis as well as other genome annotation information obtained from ENSEMBL, TRANSFAC, TRED (Transcription Regulatory Element Database), ORegAnno and the ENCODE project, was established. Finally publicly available expression data from human Leishmania infected macrophages were analyzed using the genome-wide information on predicted TFBS with the computer system ExPlain™. The gene regulatory network was constructed and activated signal transduction pathways were revealed. The Irak1 pathway was identified as a key pathway regulating gene expression changes in Leishmania infected macrophages.

Background Leishmania (L.) are intracellular protozoan parasites able to survive and replicate in the hostile phagolysosomal environment of infected macrophages. They cause leishmaniasis, a heterogeneous group of worldwide-distributed affections, representing a paradigm of neglected diseases that are mainly embedded in impoverished populations. To establish successful infection and ensure their own survival, Leishmania have developed sophisticated strategies to subvert the host macrophage responses. Despite a wealth of gained crucial information, these strategies still remain poorly understood. MicroRNAs (miRNAs), an evolutionarily conserved class of endogenous 22-nucleotide non-coding RNAs, are described to participate in the regulation of almost every cellular process investigated so far. They regulate the expression of target genes both at the levels of mRNA stability and translation; changes in their expression have a profound effect on their target transcripts. Methodology/Principal Findings We report in this study a comprehensive analysis of miRNA expression profiles in L. major-infected human primary macrophages of three healthy donors assessed at different time-points post-infection (three to 24 h). We show that expression of 64 out of 365 analyzed miRNAs was consistently deregulated upon infection with the same trends in all donors. Among these, several are known to be induced by TLR-dependent responses. GO enrichment analysis of experimentally validated miRNA-targeted genes revealed that several pathways and molecular functions were disturbed upon parasite infection. Finally, following parasite infection, miR-210 abundance was enhanced in HIF-1α-dependent manner, though it did not contribute to inhibiting anti-apoptotic pathways through pro-apoptotic caspase-3 regulation. Conclusions/Significance Our data suggest that alteration in miRNA levels likely plays an important role in regulating macrophage functions following L. major infection. These results could contribute to better understanding of the dynamics of gene expression in host cells during leishmaniasis., Author Summary Leishmania parasites belong to different species, each one characterized by specific vectors and reservoirs, and causes cutaneous or visceral disease(s) of variable clinical presentation and severity. In its mammalian host, the parasite is an obligate intracellular pathogen infecting the monocyte/macrophage lineage. Leishmania have developed ambiguous relationships with macrophages. Indeed, these cells are the shelter of invading parasites, where they will grow and eventually will reside in a silent state for life. But macrophages are also the cells that participate, through the induction of several pro-inflammatory mediators and antigen presentation, to shape the host immune response and ultimately kill the invader. To subvert these anti-parasite responses, Leishmania manipulate the host machinery for their own differentiation and survival. We aimed to evaluate the impact of L. major (the causative agent of zoonotic cutaneous leishmaniasis) infection on deregulation of non-coding miRNAs, a class of important regulators of gene expression. Our results revealed the implication of several miRNAs on macrophage fate upon parasite infection through regulation of different pathways, including cell death. Our findings provided a new insight for understanding mechanisms governing this miRNA deregulation by parasite infection and will help to provide clues for the development of control strategies for this disease.

The discipline of bioinformatics has developed rapidly since the complete sequencing of the first genomes in the 1990s. The development of many high-throughput techniques during the last decades has ensured that bioinformatics has grown into a discipline that overlaps with, and is required for, the modern practice of virtually every field in the life sciences. This has placed a scientific premium on the availability of skilled bioinformaticians, a qualification that is extremely scarce on the African continent. The reasons for this are numerous, although the absence of a skilled bioinformatician at academic institutions to initiate a training process and build sustained capacity seems to be a common African shortcoming. This dearth of bioinformatics expertise has had a knock-on effect on the establishment of many modern high-throughput projects at African institutes, including the comprehensive and systematic analysis of genomes from African populations, which are among the most genetically diverse anywhere on the planet. Recent funding initiatives from the National Institutes of Health and the Wellcome Trust are aimed at ameliorating this shortcoming. In this paper, we discuss the problems that have limited the establishment of the bioinformatics field in Africa, as well as propose specific actions that will help with the education and training of bioinformaticians on the continent. This is an absolute requirement in anticipation of a boom in high-throughput approaches to human health issues unique to data from African populations.

Recent efforts have reported numerous variants that influence severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral characteristics, including pathogenicity, transmission rate, and detectability by molecular tests. Whole-genome sequencing based on next-generation sequencing technologies is the method of choice to identify all viral variants; however, the resources needed to use these techniques for a representative number of specimens remain limited in many low- and middle-income countries. To decrease sequencing costs, we developed a primer set allowing partial sequences to be generated in the viral S gene, enabling rapid detection of numerous variants of concern (VOCs) and variants of interest (VOIs); whole-genome sequencing is then performed on a selection of viruses based on partial sequencing results. Two hundred one nasopharyngeal specimens collected during the decreasing phase of a high-transmission COVID-19 wave in Tunisia were analyzed. The results reveal high genetic variability within the sequenced fragment and allow the detection of first introductions in the country of already-known VOCs and VOIs, as well as other variants that have interesting genomic mutations and need to be kept under surveillance. IMPORTANCE The method of choice for SARS-CoV-2 variant detection is whole-genome sequencing using next-generation sequencing (NGS) technologies. Resources for this technology remain limited in many low- and middle-income countries, where it is not possible to perform whole-genome sequencing for representative numbers of SARS-CoV-2-positive cases. In the present work, we developed a novel strategy based on a first partial Sanger screening in the S gene, which includes key mutations of the already known VOCs and VOIs, for rapid identification of these VOCs and VOIs and to help better select specimens that need to be sequenced by NGS technologies. The second step consists of whole-genome sequencing to allow a holistic view of all variants within the selected viral strains and confirm the initial classification of the strains based on partial S gene sequencing.