A Comprehensive Analysis of 3 Moroccan Genomes Revealed Contributions From Both African and European Ancestries.

Genetic variations in the human genome represent the differences in DNA sequence within individuals. This highlights the important role of whole human genome sequencing which has become the keystone for precision medicine and disease prediction. Morocco is an important hub for studying human population migration and mixing history. This study presents the analysis of 3 Moroccan genomes; the variant analysis revealed 6 379 606 single nucleotide variants (SNVs) and 1 050 577 small InDels. Of those identified SNVs, 219 152 were novel, with 1233 occurring in coding regions, and 5580 non-synonymous single nucleotide variants (nsSNP) variants were predicted to affect protein functions. The InDels produced 1055 coding variants and 454 non-3n length variants, and their size ranged from -49 and 49 bp. We further analysed the gene pathways of 8 novel coding variants found in the 3 genomes and revealed 5 genes involved in various diseases and biological pathways. We found that the Moroccan genomes share 92.78% of African ancestry, and 92.86% of Non-Finnish European ancestry, according to the gnomAD database. Then, population structure inference, by admixture analysis and network-based approach, revealed that the studied genomes form a mixed population structure, highlighting the increased genetic diversity in Morocco.

Targeted Gene Panel Sequencing Unveiled New Pathogenic Mutations in Patients With Breast Cancer.

The increasing commercialization of new gene panels based on next-generation sequencing for clinical research has significantly improved our understanding of breast cancer genetics and has led to the discovery of new mutation variants. The study included 16 unselected Moroccan breast cancer patients tested with multi-gene panel (HEVA screen panel) using Illumina Miseq, followed by Sanger sequencing to validate the most relevant mutation. Mutational analysis revealed the presence of 13 mutations (11 single-nucleotide polymorphisms [SNPs] and 2 indels), and 6 of 11 identified SNPs were predicted as pathogenic. One of the 6 pathogenic mutations was c.7874G>C, a heterozygous SNP in HD-OB domain of BRCA2 gene, which led to the arginine to threonine change at codon 2625 of the protein. This work describes the first case of a patient with breast cancer harboring this pathogenic variant and analyzes its functional impact using molecular docking and molecular dynamics simulation. Further experimental investigations are needed to validate its pathogenicity and to verify its association with breast cancer.

In Silico Analyses of All STAT3 Missense Variants Leading to Explore Divergent AD-HIES Clinical Phenotypes.

Autosomal dominant hyper-IgE syndrome (AD-HIES) is linked to dominant negative mutations of the STAT3 protein whose molecular basis for dysfunction is unclear and presenting with a variety of clinical manifestations with only supportive treatment. To establish the relationship between the impact of STAT3 mutations in different domains and the severity of the clinical manifestations, 105 STAT3 mutations were analyzed for their impact on protein stability, flexibility, function, and binding affinity using in Silico approaches. Our results showed that 73% of the studied mutations have an impact on the physicochemical properties of the protein, altering the stability, flexibility and function to varying degrees. In particular, mutations affecting the DNA binding domain (DBD) and the Src Homology 2 (SH2) have a significant impact on the protein structure and disrupt its interaction either with DNA or other STAT3 to form a heterodomain complex, leading to severe clinical phenotypes. Collectively, this study suggests that there is a close relationship between the domain involving the mutation, the degree of variation in the properties of the protein and the degree of loss of function ranging from partial loss to complete loss, explaining the variability of clinical manifestations between mild and severe.

Metagenomics Analysis of Breast Microbiome Highlights the Abundance of Rothia Genus in Tumor Tissues.

Breast cancer is one of the main global priorities in terms of public health. It remains the most frequent cancer in women and is the leading cause of their death. The human microbiome plays various roles in maintaining health by ensuring a dynamic balance with the host or in the appearance of various pathologies including breast cancer. In this study, we performed an analysis of bacterial signature differences between tumor and adjacent tissues of breast cancer patients in Morocco. Using 16S rRNA gene sequencing, we observed that adjacent tissue contained a much higher percentage of the Gammaproteobacteria class (35.7%) while tumor tissue was characterized by a higher percentage of Bacilli and Actinobacteria classes, with about 18.8% and 17.2% average abundance, respectively. Analysis of tumor subtype revealed enrichment of genus Sphingomonodas in TNBC while Sphingomonodas was predominant in HER2. The LEfSe and the genus level heatmap analysis revealed a higher abundance of the Rothia genus in tumor tissues. The identified microbial communities can therefore serve as potential biomarkers for prognosis and diagnosis, while also helping to develop new strategies for the treatment of breast cancer patients.

Draft Genome Sequences of Three Pseudomonas chengduensis Strains Isolated from Desert Soil in Morocco.

We report the draft genome sequences of three Pseudomonas chengduensis strains isolated from the sand dunes of the Merzouga (MDMC17 strain) and Erg Lihoudi (MDMC216 and MDMC224 strains) regions in the Moroccan desert. These bacteria are able to tolerate the harsh environmental conditions of the desert ecosystem.

Draft Genome Sequences of Two Enterobacter hormaechei subsp. xiangfangensis Strains Isolated from the Moroccan Sahara.

We report the draft genome sequences of Enterobacter hormaechei subsp. xiangfangensis strains MDMC82 and MDMC76, which were isolated from the sand dunes of the Merzouga desert in the Moroccan Sahara. These bacteria are able to tolerate the harsh environmental conditions of the Moroccan desert.

Genomic analysis of Paenibacillus sp. MDMC362 from the Merzouga desert leads to the identification of a potentially thermostable catalase.

Microorganisms in hot deserts face heat and other environmental conditions, such as desiccation, UV radiation, or low nutrient availability. Therefore, this hostile environment harbour microorganisms with acquired characteristics related to survival in their habitat, which can be exploited in biotechnology. In this work, the genome of Paenibacillus sp. MDMC362 isolated from the Merzouga desert in Morocco was sequenced to understand its survival strategy's genetic basis; and to evaluate the thermostability of a catalase extracted from genomic annotation files using molecular dynamics. Paenibacillus sp. MDMC362 genome was rich in genetic elements involved in the fight against different stresses, notably temperature stress, UV radiations, osmotic stress, carbon starvation, and oxidative stress. Indeed, we could identify genes of the operons groES-groEL and hrcA-grpE-dnaK and those involved in the different stages of sporulation, which can help the bacteria to survive the high temperatures imposed by a desertic environment. We also observed the genetic components of the UvrABC system and additional mechanisms involved in DNA repair, which help overcome UV radiation damage. Other genes have been identified in the genome, like those coding for ectoine and proline, that aids fight osmotic stress and desiccation. Catalase thermostability investigation using molecular dynamics showed that the protein reached stability and conserved its compactness at temperatures up to 373.15 K. These results suggest a potential thermostability of the enzyme. Since the studied protein is a core protein, thermostability could be conserved among Paenibacillus sp. MDMC362 closely related strains; however, bacteria from harsh environments may have a slight advantage regarding protein stability.

Genomic analysis of two Bacillus safensis isolated from Merzouga desert reveals desert adaptive and potential plant growth-promoting traits.

Deserts represent extreme environments for microorganisms, and conditions such as high soil salinity, nutrient deficiency, and increased levels of UV radiation make desert soil communities of high biotechnological potential. In this study, we isolated, sequenced, and assembled the genomes of Bacillus safensis strains BcP62 and Bcs93, to which we performed comparative genome analyses. Using the DDH and ANI of both strains with the available B. safensis genomes, we identified three potential subspecies within this group. Intra-species core genome phylogenetic analysis did not result in clustering genomes by niche type, with some exceptions. This study also revealed that the genomes of the analyzed strains possessed plant growth-promoting characteristics, most of which were conserved in all B. safensis strains. Furthermore, we highlight the genetic features of B. safensis BcP62 and Bcs93 related to survival in the Merzouga desert in Morocco. These strains could be potentially used in agriculture as PGPB in extreme environments, given their high tolerability to unfavorable conditions.

IDbSV: An Open-Access Repository for Monitoring SARS-CoV-2 Variations and Evolution.

Ending COVID-19 pandemic requires a collaborative understanding of SARS-CoV-2 and COVID-19 mechanisms. Yet, the evolving nature of coronaviruses results in a continuous emergence of new variants of the virus. Central to this is the need for a continuous monitoring system able to detect potentially harmful variants of the virus in real-time. In this manuscript, we present the International Database of SARS-CoV-2 Variations (IDbSV), the result of ongoing efforts in curating, analyzing, and sharing comprehensive interpretation of SARS-CoV-2's genetic variations and variants. Through user-friendly interactive data visualizations, we aim to provide a novel surveillance tool to the scientific and public health communities. The database is regularly updated with new records through a 4-step workflow (1-Quality control of curated sequences, 2-Call of variations, 3-Functional annotation, and 4-Metadata association). To the best of our knowledge, IDbSV provides access to the largest repository of SARS-CoV-2 variations and the largest analysis of SARS-CoV-2 genomes with over 60 thousand annotated variations curated from the 1,808,613 genomes alongside their functional annotations, first known appearance, and associated genetic lineages, enabling a robust interpretation tool for SARS-CoV-2 variations to help understanding SARS-CoV-2 dynamics across the world.

Variants of human CLDN9 cause mild to profound hearing loss.

Hereditary deafness is clinically and genetically heterogeneous. We investigated deafness segregating as a recessive trait in two families. Audiological examinations revealed an asymmetric mild to profound hearing loss with childhood or adolescent onset. Exome sequencing of probands identified a homozygous c.475G>A;p.(Glu159Lys) variant of CLDN9 (NM_020982.4) in one family and a homozygous c.370_372dupATC;p.(Ile124dup) CLDN9 variant in an affected individual of a second family. Claudin 9 (CLDN9) is an integral membrane protein and constituent of epithelial bicellular tight junctions (TJs) that form semipermeable, paracellular barriers between inner ear perilymphatic and endolymphatic compartments. Computational structural modeling predicts that substitution of a lysine for glutamic acid p.(Glu159Lys) alters one of two cis-interactions between CLDN9 protomers. The p.(Ile124dup) variant is predicted to locally misfold CLDN9 and mCherry tagged p.(Ile124dup) CLDN9 is not targeted to the HeLa cell membrane. In situ hybridization shows that mouse Cldn9 expression increases from embryonic to postnatal development and persists in adult inner ears coinciding with prominent CLDN9 immunoreactivity in TJs of epithelia outlining the scala media. Together with the Cldn9 deaf mouse and a homozygous frameshift of CLDN9 previously associated with deafness, the two bi-allelic variants of CLDN9 described here point to CLDN9 as a bona fide human deafness gene.

Large scale genomic analysis of 3067 SARS-CoV-2 genomes reveals a clonal geo-distribution and a rich genetic variations of hotspots mutations.

In late December 2019, an emerging viral infection COVID-19 was identified in Wuhan, China, and became a global pandemic. Characterization of the genetic variants of SARS-CoV-2 is crucial in following and evaluating it spread across countries. In this study, we collected and analyzed 3,067 SARS-CoV-2 genomes isolated from 55 countries during the first three months after the onset of this virus. Using comparative genomics analysis, we traced the profiles of the whole-genome mutations and compared the frequency of each mutation in the studied population. The accumulation of mutations during the epidemic period with their geographic locations was also monitored. The results showed 782 variants sites, of which 512 (65.47%) had a non-synonymous effect. Frequencies of mutated alleles revealed the presence of 68 recurrent mutations, including ten hotspot non-synonymous mutations with a prevalence higher than 0.10 in this population and distributed in six SARS-CoV-2 genes. The distribution of these recurrent mutations on the world map revealed that certain genotypes are specific to geographic locations. We also identified co-occurring mutations resulting in the presence of several haplotypes. Moreover, evolution over time has shown a mechanism of mutation co-accumulation which might affect the severity and spread of the SARS-CoV-2. The phylogentic analysis identified two major Clades C1 and C2 harboring mutations L3606F and G614D, respectively and both emerging for the first time in China. On the other hand, analysis of the selective pressure revealed the presence of negatively selected residues that could be taken into considerations as therapeutic targets. We have also created an inclusive unified database (http://covid-19.medbiotech.ma) that lists all of the genetic variants of the SARS-CoV-2 genomes found in this study with phylogeographic analysis around the world.

Draft Genome Sequences of Six Moroccan Helicobacter pylori Isolates Belonging to the hspWAfrica Group.

Helicobacter pylori affects up to 50% of people worldwide. Here, we present the draft genome sequences of six H. pylori strains isolated from Moroccan patients with different gastric diseases. Multilocus sequence typing analysis showed that all of the H. pylori isolates belonged to the hspWAfrica group.

Genomic Diversity and Hotspot Mutations in 30,983 SARS-CoV-2 Genomes: Moving Toward a Universal Vaccine for the "Confined Virus"?

The COVID-19 pandemic has been ongoing since its onset in late November 2019 in Wuhan, China. Understanding and monitoring the genetic evolution of the virus, its geographical characteristics, and its stability are particularly important for controlling the spread of the disease and especially for the development of a universal vaccine covering all circulating strains. From this perspective, we analyzed 30,983 complete SARS-CoV-2 genomes from 79 countries located in the six continents and collected from 24 December 2019, to 13 May 2020, according to the GISAID database. Our analysis revealed the presence of 3206 variant sites, with a uniform distribution of mutation types in different geographic areas. Remarkably, a low frequency of recurrent mutations has been observed; only 169 mutations (5.27%) had a prevalence greater than 1% of genomes. Nevertheless, fourteen non-synonymous hotspot mutations (>10%) have been identified at different locations along the viral genome; eight in ORF1ab polyprotein (in nsp2, nsp3, transmembrane domain, RdRp, helicase, exonuclease, and endoribonuclease), three in nucleocapsid protein, and one in each of three proteins: Spike, ORF3a, and ORF8. Moreover, 36 non-synonymous mutations were identified in the receptor-binding domain (RBD) of the spike protein with a low prevalence (<1%) across all genomes, of which only four could potentially enhance the binding of the SARS-CoV-2 spike protein to the human ACE2 receptor. These results along with intra-genomic divergence of SARS-CoV-2 could indicate that unlike the influenza virus or HIV viruses, SARS-CoV-2 has a low mutation rate which makes the development of an effective global vaccine very likely.

Draft Genome Sequence of Stenotrophomonas maltophilia MDMC339, Isolated from Soil of Merzouga Desert in Morocco.

Here, we report the draft genome sequence of Stenotrophomonas maltophilia MDMC339, a strain able to survive in the difficult conditions imposed by the Merzouga desert. The analyzed genome contains 4,788,525 bp with 4,262 genes coding for proteins, including several genes related to stress.

Genome Sequences of Six SARS-CoV-2 Strains Isolated in Morocco, Obtained Using Oxford Nanopore MinION Technology.

Here, we report the draft genome sequences of six severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) strains. SARS-CoV-2 is responsible for the COVID-19 pandemic, which started at the end of 2019 in Wuhan, China. The isolates were obtained from nasopharyngeal swabs from Moroccan patients with COVID-19. Mutation analysis revealed the presence of the spike D614G mutation in all six genomes, which is widely present in several genomes around the world.

First Whole-Genome Sequences of Two Multidrug-Resistant Acinetobacter baumannii Strains Isolated from a Moroccan Hospital Floor.

This report describes the whole-genome shotgun sequences of two multidrug-resistant Acinetobacter baumannii strains, ABE8_07 and ABE12_M, isolated from a Moroccan hospital floor. These two genome sequences will initiate the study and characterization of the Acinetobacter baumannii genome in Morocco.