Neutrophil extracellular trap formation and gene programs distinguish TST/IGRA sensitization outcomes among Mycobacterium tuberculosis exposed persons living with HIV.

Persons living with HIV (PLWH) have an increased risk for tuberculosis (TB). After prolonged and repeated exposure, some PLWH never develop TB and show no evidence of immune sensitization to Mycobacterium tuberculosis (Mtb) as defined by persistently negative tuberculin skin tests (TST) and interferon gamma release assays (IGRA). This group has been identified and defined as HIV+ persistently TB, tuberculin and IGRA negative (HITTIN). To investigate potential innate mechanisms unique to individuals with the HITTIN phenotype we compared their neutrophil Mtb infection response to that of PLWH, with no TB history, but who test persistently IGRA positive, and tuberculin positive (HIT). Neutrophil samples from 17 HITTIN (PMNHITTIN) and 11 HIT (PMNHIT) were isolated and infected with Mtb H37Rv for 1h and 6h. RNA was extracted and used for RNAseq analysis. Since there was no significant differential transcriptional response at 1h between infected PMNHITTIN and PMNHIT, we focused on the 6h timepoint. When compared to uninfected PMN, PMNHITTIN displayed 3106 significantly upregulated and 3548 significantly downregulated differentially expressed genes (DEGs) (absolute cutoff of a log2FC of 0.2, FDR < 0.05) whereas PMNHIT demonstrated 3816 significantly upregulated and 3794 significantly downregulated DEGs following 6h Mtb infection. Contrasting the log2FC 6h infection response to Mtb from PMNHITTIN against PMNHIT, 2285 genes showed significant differential response between the two groups. Overall PMNHITTIN had a lower fold change response to Mtb infection compared to PMNHIT. According to pathway enrichment, Apoptosis and NETosis were differentially regulated between HITTIN and HIT PMN responses after 6h Mtb infection. To corroborate the blunted NETosis transcriptional response measured among HITTIN, fluorescence microscopy revealed relatively lower neutrophil extracellular trap formation and cell loss in PMNHITTIN compared to PMNHIT, showing that PMNHITTIN have a distinct response to Mtb.

Challenges of Diagnosing Mendelian Susceptibility to Mycobacterial Diseases in South Africa.

Inborn errors of immunity (IEI) are genetic disorders with extensive clinical presentations. They can range from increased susceptibility to infections to significant immune dysregulation that results in immune impairment. While IEI cases are individually rare, they collectively represent a significant burden of disease, especially in developing countries such as South Africa, where infectious diseases like tuberculosis (TB) are endemic. This is particularly alarming considering that certain high penetrance mutations that cause IEI, such as Mendelian Susceptibility to Mycobacterial Disease (MSMD), put individuals at higher risk for developing TB and other mycobacterial diseases. MSMD patients in South Africa often present with different clinical phenotypes than those from the developed world, therefore complicating the identification of disease-associated variants in this setting with a high burden of infectious diseases. The lack of available data, limited resources, as well as variability in clinical phenotype are the reasons many MSMD cases remain undetected or misdiagnosed. This article highlights the challenges in diagnosing MSMD in South Africa and proposes the use of transcriptomic analysis as a means of potentially identifying dysregulated pathways in affected African populations.

Delineating the Spread and Prevalence of SARS-CoV-2 Omicron Sublineages (BA.1-BA.5) and Deltacron Using Wastewater in the Western Cape, South Africa.

This study was one of the first to detect Omicron sublineages BA.4 and BA.5 in wastewater from South Africa. Spearman rank correlation analysis confirmed a strong positive correlation between severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral RNA in wastewater samples and clinical cases (r = 0.7749, P < .0001). SARS-CoV-2 viral load detected in wastewater, resulting from the Delta-driven third wave, was significantly higher than during the Omicron-driven fourth wave. Whole-genome sequencing confirmed presence of Omicron lineage defining mutations in wastewater with the first occurrence reported 23 November 2021 (BA.1 predominant). The variant spread rapidly, with prevalence of Omicron-positive wastewater samples rising to >80% by 10 January 2022 with BA.2 as the predominant sublineage by 10 March 2022, whilst on 18 April 2022 BA.4 and BA.5 were detected in selected wastewater sites. These findings demonstrate the value of wastewater-based epidemiology to monitor the spatiotemporal spread and potential origin of new Omicron sublineages.

A multi-phenotype genome-wide association study of clades causing tuberculosis in a Ghanaian- and South African cohort.

Despite decades of research and advancements in diagnostics and treatment, tuberculosis remains a major public health concern. New computational methods are needed to interrogate the intersection of host- and bacterial genomes. Paired host genotype datum and infecting bacterial isolate information were analysed for associations using a multinomial logistic regression framework implemented in SNPTest. A cohort of 853 admixed South African participants and a Ghanaian cohort of 1359 participants were included. Two directly genotyped variants, namely rs529920 and rs41472447, were identified in the Ghanaian cohort as being statistically significantly associated with risk for infection with strains of different members of the MTBC. Thus, a multinomial logistic regression using paired host-pathogen data may prove valuable for investigating the complex relationships driving infectious disease.

Phenotypic and immune functional profiling of patients with suspected Mendelian Susceptibility to Mycobacterial Disease in South Africa.

BACKGROUND: Mendelian Susceptibility to Mycobacterial Disease (MSMD) is a primary immunodeficiency (PID) characterised by a predisposition to infection by weakly-pathogenic mycobacteria. In countries with a high prevalence of tuberculosis (TB), individuals with MSMD are also prone to infections by Mycobacterium tuberculosis. Several MSMD-associated genes have been described, all resulting in a disruption of IL-12 and IFN-γ cytokine axis, which is essential for control of mycobacterial infections. An accurate molecular diagnosis, confirmed by phenotypic and functional immune investigations, is essential to ensure that the patient receives optimal treatment and prophylaxis for infections. The aim of this study was to implement a set of functional assays to assess the integrity of the IL-12-IFN-γ cytokine pathways in patients presenting with severe, persistent, unusual and/or recurrent TB, mycobacterial infections or other clinical MSMD-defining infections such as Salmonella. METHODS: Blood was collected for subsequent PBMC isolation from 16 participants with MSMD-like clinical phenotypes. A set of flow cytometry (phenotype and signalling integrity) and ELISA-based (cytokine production) functional assays were implemented to assess the integrity of the IL-12-IFN-γ pathway. RESULTS: The combination of the three assays for the assessment of the integrity of the IL-12-IFN-γ pathway was successful in identifying immune deficits in the IL-12-IFN-γ pathway in all of the participants included in this study. CONCLUSIONS: The data presented here emphasise the importance of investigating PID and TB susceptibility in TB endemic regions such as South Africa as MSMD and other previously described PIDs relating to TB susceptibility may present differently in such regions. It is therefore important to have access to in vitro functional investigations to better understand the immune function of these individuals. Although functional assays alone are unlikely to always provide a clear diagnosis, they do give an overview of the integrity of the IL-12-IFN-γ pathway. It would be beneficial to apply these assays routinely to patients with suspected PID relating to mycobacterial susceptibility. A molecular diagnosis with confirmed functional impairment paves the way for targeted treatment and improved disease management options for these patients.

A regulatory variant in the C1Q gene cluster is associated with tuberculosis susceptibility and C1qA plasma levels in a South African population.

Several genetic studies have implicated genes that encode for components of the innate immune response in tuberculosis (TB) susceptibility. The complement system is an early player in the innate immune response and provides the host with initial protection by promoting phagocytosis of apoptotic or necrotic cells. The C1q molecule is the first component of the classical pathway that leads to the activation of complement by binding to immune complexes and is encoded by the C1Q gene cluster. We investigated variants in this region to determine its association with TB susceptibility. Five single nucleotide polymorphisms (SNPs) (rs12033074, rs631090, rs172378, rs587585, and rs665691) were genotyped using TaqMan® SNP assays in 456 TB cases and 448 healthy controls and analysed by logistic regression models. The rs587585 variant showed a significant additive allelic association where the minor G allele was found more frequently in TB cases than in controls in both the discovery (p = 0.023; OR = 1.30; 95% CI, 1.04-1.64) and validation cohort (p = 0.038; OR = 1.31; 95% CI, 1.22-1.40). In addition, we detected increased C1qA expression when comparing cases and controls (p = 0.037) and linked this to a dosage effect of the G allele, which increased C1qA expression in TB cases. This is the first study to report the association of C1Q gene polymorphisms with progression to tuberculosis.

Identification of a novel WAS mutation in a South African patient presenting with atypical Wiskott-Aldrich syndrome: a case report.

BACKGROUND: The X-linked recessive primary immunodeficiency disease (PIDD) Wiskott-Aldrich syndrome (WAS) is identified by an extreme susceptibility to infections, eczema and thrombocytopenia with microplatelets. The syndrome, the result of mutations in the WAS gene which encodes the Wiskott-Aldrich protein (WASp), has wide clinical phenotype variation, ranging from classical WAS to X-linked thrombocytopaenia and X-linked neutropaenia. In many cases, the diagnosis of WAS in first affected males is delayed, because patients may not present with the classic signs and symptoms, which may intersect with other thrombocytopenia causes. CASE PRESENTATION: Here, we describe a three-year-old HIV negative boy presenting with recurrent infections, skin rashes, features of autoimmunity and atopy. However, platelets were initially reported as normal in numbers and morphology as were baseline immune investigations. An older male sibling had died in infancy from suspected immunodeficiency. Uncertainty of diagnosis and suspected severe PIDD prompted urgent further molecular investigation. Whole exome sequencing identified c. 397 G > A as a novel hemizygous missense mutation located in exon 4 of WAS. CONCLUSION: With definitive molecular diagnosis, we could target treatment and offer genetic counselling and prenatal diagnostic testing to the family. The identification of novel variants is important to confirm phenotype variations of a syndrome.

The X chromosome and sex-specific effects in infectious disease susceptibility.

The X chromosome and X-linked variants have largely been ignored in genome-wide and candidate association studies of infectious diseases due to the complexity of statistical analysis of the X chromosome. This exclusion is significant, since the X chromosome contains a high density of immune-related genes and regulatory elements that are extensively involved in both the innate and adaptive immune responses. Many diseases present with a clear sex bias, and apart from the influence of sex hormones and socioeconomic and behavioural factors, the X chromosome, X-linked genes and X chromosome inactivation mechanisms contribute to this difference. Females are functional mosaics for X-linked genes due to X chromosome inactivation and this, combined with other X chromosome inactivation mechanisms such as genes that escape silencing and skewed inactivation, could contribute to an immunological advantage for females in many infections. In this review, we discuss the involvement of the X chromosome and X inactivation in immunity and address its role in sexual dimorphism of infectious diseases using tuberculosis susceptibility as an example, in which male sex bias is clear, yet not fully explored.

Hyperphosphatasia with mental retardation syndrome type 4 in three unrelated South African patients.

Hyperphosphatasia with mental retardation syndrome (HPMRS) is a rare autosomal recessive disorder caused by pathogenic variants in genes involved in glycosylphosphatidylinositol metabolism that result in a similar phenotype. We describe the first three patients with HPMRS from sub-Saharan Africa. Detection was assisted by Face2Gene phenotype matching and confirmed by the presence of elevated serum alkaline phosphatase. All three patients had severe intellectual disability, absent speech, hypotonia and palatal abnormality (cleft palate in two, very high-arched palate in one), no or minimal brachytelephalangy, and high serum alkaline phosphatase levels. Additional findings included seizures in two, and brain imaging abnormalities in two. In all three patients HPMRS was a top-20 gestalt match using Face2Gene. The overall phenotype is consistent with descriptions in the literature of HPMRS type 4, although not specific to it. Whole exome sequencing in the index patient and his mother detected a candidate variant in a homozygous state in the index patient (PGAP3:c.557G>C, p.Arg186Thr) and heterozygous in the mother. Further variant interpretation indicated pathogenicity. Sanger sequencing of another two patients identified the same homozygous, pathogenic variant, confirming a diagnosis of HPMRS type 4. The shared homozygous variant in apparently unrelated families, and in the absence of consanguinity, suggests the possibility of genetic drift due to a population bottleneck effect, and further research is recommended.

Human global and population-specific genetic susceptibility to Mycobacterium tuberculosis infection and disease.

PURPOSE OF REVIEW: Multiple lines of evidence support a role of the host genetic component in Mycobacterium tuberculosis infection and disease progression. However, genomic studies of tuberculosis susceptibility have been disappointing compared with that of other complex disorders. Recently the field has explored alternative strategies to facilitate locus discovery. Results emanating from these efforts during the last 18 months are addressed in this review. RECENT FINDINGS: There has been a renewed focus on the refinement of phenotypic definitions of infection and disease as well as on age-related, sex-specific and population-specific effects. Genome-wide association studies have yielded candidate genes but the findings have not always been transferable to all population groups. Candidate gene association studies remain popular as it is used for GWAS replication and is affordable, particularly in lower and middle-income countries. Pharmacogenetic studies involving tuberculosis drugs may locate variants that can be cost-effectively genotyped to identify individuals at risk of developing adverse events during treatment. SUMMARY: Additional GWAS and candidate gene association studies of crudely defined study participants are unlikely to make further important contributions to the TB susceptibility field. Instead refined phenotyping will allow the elucidation of genetic mechanisms contributing to infection and disease in distinct populations and the calculation of polygenic risk scores.

Wild-type and mutant (G2019S) leucine-rich repeat kinase 2 (LRRK2) associate with subunits of the translocase of outer mitochondrial membrane (TOM) complex.

Leucine-rich repeat kinase 2 (LRRK2) is important in various cellular processes including mitochondrial homeostasis and mutations in this gene lead to Parkinson's disease (PD). However, the full spectrum of LRRK2's functions remain to be elucidated. The translocase of outer mitochondrial membrane (TOM) complex is essential for the import of almost all nuclear-encoded mitochondrial proteins and is fundamental for cellular survival. Using co-immunoprecipitation, super-resolution structured illumination microscopy (SR-SIM), and 3D virtual reality (VR) assisted co-localization analysis techniques we show that wild-type and mutant (G2019S) LRRK2 associate and co-localize with subunits of the TOM complex, either under basal (dimethyl sulfoxide, DMSO) or stress-induced (carbonyl cyanide m-chlorophenyl hydrazine, CCCP) conditions. Interestingly, LRRK2 interacted with TOM40 under both DMSO and CCCP conditions, and when the PD causing mutation, G2019S was introduced, the association was not altered. Moreover, overexpression of G2019S LRRK2 resulted in the formation of large, perinuclear aggregates that co-localized with the TOM complex. Taken together, this is the first study to show that both WT and mutant LRRK2 associate with the TOM complex subunits. These findings provide additional evidence for LRRK2's role in mitochondrial function which has important implications for its role in PD pathogenesis.

WD40-repeat 47, a microtubule-associated protein, is essential for brain development and autophagy.

The family of WD40-repeat (WDR) proteins is one of the largest in eukaryotes, but little is known about their function in brain development. Among 26 WDR genes assessed, we found 7 displaying a major impact in neuronal morphology when inactivated in mice. Remarkably, all seven genes showed corpus callosum defects, including thicker (Atg16l1, Coro1c, Dmxl2, and Herc1), thinner (Kif21b and Wdr89), or absent corpus callosum (Wdr47), revealing a common role for WDR genes in brain connectivity. We focused on the poorly studied WDR47 protein sharing structural homology with LIS1, which causes lissencephaly. In a dosage-dependent manner, mice lacking Wdr47 showed lethality, extensive fiber defects, microcephaly, thinner cortices, and sensory motor gating abnormalities. We showed that WDR47 shares functional characteristics with LIS1 and participates in key microtubule-mediated processes, including neural stem cell proliferation, radial migration, and growth cone dynamics. In absence of WDR47, the exhaustion of late cortical progenitors and the consequent decrease of neurogenesis together with the impaired survival of late-born neurons are likely yielding to the worsening of the microcephaly phenotype postnatally. Interestingly, the WDR47-specific C-terminal to LisH (CTLH) domain was associated with functions in autophagy described in mammals. Silencing WDR47 in hypothalamic GT1-7 neuronal cells and yeast models independently recapitulated these findings, showing conserved mechanisms. Finally, our data identified superior cervical ganglion-10 (SCG10) as an interacting partner of WDR47. Taken together, these results provide a starting point for studying the implications of WDR proteins in neuronal regulation of microtubules and autophagy.

Primary immunodeficiency diseases in a tuberculosis endemic region: challenges and opportunities.

While individual primary immunodeficiency diseases (PIDs) are rare, collectively they represent a significant burden of disease. Recent estimates show that about one million people in Africa suffer from a PID. However, data from African PID registries reflect only a small percentage of the estimated prevalence. This disparity is partly due to the lack of PID awareness and the masking of PIDs by the endemic pathogens. Over three million tuberculosis (TB) cases were reported in Africa in 2016, with many of these from southern Africa. Despite concerted efforts to address this high burden of disease, the underlying genetic correlates of susceptibility to TB remain poorly understood. High penetrance mutations in immune system genes can cause PIDs that selectively predispose individuals to TB and other mycobacterial diseases. Additionally, the identification of individuals at a heightened risk of developing TB or of presenting with severe or disseminated TB due to their genetic ancestry is crucial to promote a positive treatment outcome. The screening for and identification of PID mutations in TB-endemic regions by next-generation sequencing (NGS) represents a promising approach to improve the understanding of what constitutes an effective immune response to TB, as well as the range of associated PIDs and phenotypes.

Exome sequencing identifies a novel TTC37 mutation in the first reported case of Trichohepatoenteric syndrome (THE-S) in South Africa.

BACKGROUND: Trichohepatoenteric syndrome (THE-S) or phenotypic diarrhoea of infancy is a rare autosomal recessive disorder characterised by severe infantile diarrhoea, facial dysmorphism, immunodeficiency and woolly hair. It was first described in 1982 in two infants with intractable diarrhoea, liver cirrhosis and abnormal hair structure on microscopy. We report on two siblings from a consanguineous family of Somali descent who, despite extensive clinical investigation, remained undiagnosed until their demise. The index patient died of fulminant cytomegalovirus pneumonitis at 3 months of age. METHODS: Whole exome sequencing (WES) was performed on a premortem DNA sample from the index case. Variants in a homozygous recessive state or compound heterozygous state were prioritized as potential candidate variants using TAPER™. Sanger sequencing was done to genotype the parents, unaffected sibling and a deceased sibling for the variant of interest. RESULTS: Exome sequencing identified a novel homozygous mutation (c.4507C > T, rs200067423) in TTC37 which was confirmed by Sanger sequencing in the index case. The identification of this mutation led to the diagnosis of THE-S in the proband and the same homozygous variant was confirmed in a male sibling who died 4 years earlier with severe chronic diarrhoea of infancy. The unaffected parents and sister were heterozygous for the identified variant. CONCLUSIONS: WES permitted definitive genetic diagnosis despite an atypical presentation in the index case and suggests that severe infection, likely secondary to immunodeficiency, may be a presenting feature. In addition definitive molecular diagnosis allows for genetic counseling and future prenatal diagnosis, and demonstrates the value of WES for post-mortem diagnosis of disorders with a non-specific clinical presentation in which a Mendelian cause is suspected.

Ascribing novel functions to the sarcomeric protein, myosin binding protein H (MyBPH) in cardiac sarcomere contraction.

Myosin binding protein H (MyBPH) is a protein of unknown function, which shares sequence and structural similarities with myosin binding protein C (cMyBPC), a protein frequently implicated in hypertrophic cardiomyopathy (HCM). Given the similarity between cMyBPC and MyBPH, we proposed that MyBPH, like cMyBPC, could be involved in HCM pathogenesis and we therefore sought to determine its function. We identified MyBPH-interacting proteins by using yeast two-hybrid (Y2H) analysis. The role of MyBPH and cMyBPC in cardiac cell contractility was analysed by measuring the planar cell surface area of differentiated H9c2 rat cardiomyocytes in response to ß-adrenergic stress after siRNA knockdown of MyBPH and cMyBPC. Individual knockdown of either protein had no effect on cardiac contractility, while concurrent knockdowns reduced cardiac contractility. These proteins therefore functionally compensate for one another and are critical for cardiac contractility. We further show that both proteins co-localise with the autophagosomal membrane protein LC3, suggesting that both proteins are involved in autophagosomal membrane maturation processes. The results of this study ascribe novel functions to MyBPH, which may contribute to our understanding of its role in the sarcomere. This study provides evidence for a potential role of MyBPH in HCM, which warrants further investigation.

Polymorphisms in the Pattern Recognition Receptor Mincle Gene (CLEC4E) and Association with Tuberculosis.

The mechanisms involved in interactions between Mycobacterium tuberculosis and host innate immune cells determine outcome. Antigen-presenting cells, including macrophages and dendritic cells, express many pattern recognition receptors to identify pathogen-associated molecular patterns, thereby initiating an immune response. A major mycobacterial virulence factor, trehalose-6',6-dimycolate, is recognised by the macrophage-inducible C-type lectin, Mincle, which leads to the activation of the Syk-Card9 signalling pathway in macrophages. Mincle is encoded by CLEC4E, and we investigated polymorphisms in this gene to assess its role in tuberculosis susceptibility. Four tagging single nucleotide polymorphisms (SNPs) (rs10841845, rs10841847, rs10841856 and rs4620776) were genotyped using TaqMan(®) SNP assays in 416 tuberculosis cases and 405 healthy controls. Logistic regression models were used for analysis. No association was detected with any of the SNPs analysed. This research highlights tuberculosis disease complexity where recognition proteins which specifically bind mycobacterial glycolipids cannot be conclusively associated with the disease in genetic studies.

Mitochondrial impairment observed in fibroblasts from South African Parkinson's disease patients with parkin mutations.

Parkinson's disease (PD), defined as a neurodegenerative disorder, is characterized by the loss of dopaminergic neurons in the substantia nigra in the midbrain. Loss-of-function mutations in the parkin gene are a major cause of autosomal recessive, early-onset PD. Parkin has been implicated in the maintenance of healthy mitochondria, although previous studies show conflicting findings regarding mitochondrial abnormalities in fibroblasts from patients harboring parkin-null mutations. The aim of the present study was to determine whether South African PD patients with parkin mutations exhibit evidence for mitochondrial dysfunction. Fibroblasts were cultured from skin biopsies obtained from three patients with homozygous parkin-null mutations, two heterozygous mutation carriers and two wild-type controls. Muscle biopsies were obtained from two of the patients. The muscle fibers showed subtle abnormalities such as slightly swollen mitochondria in focal areas of the fibers and some folding of the sarcolemma. Although no differences in the degree of mitochondrial network branching were found in the fibroblasts, ultrastructural abnormalities were observed including the presence of electron-dense vacuoles. Moreover, decreased ATP levels which are consistent with mitochondrial dysfunction were observed in the patients' fibroblasts compared to controls. Remarkably, these defects did not manifest in one patient, which may be due to possible compensatory mechanisms. These results suggest that parkin-null patients exhibit features of mitochondrial dysfunction. Involvement of mitochondria as a key role player in PD pathogenesis will have important implications for the design of new and more effective therapies.

Identification and control for the effects of bioinformatic globin depletion on human RNA-seq differential expression analysis.

When profiling blood samples by RNA-sequencing (RNA-seq), RNA from haemoglobin (Hgb) can account for up to 70% of the transcriptome. Due to considerations of sequencing depth and power to detect biological variation, Hgb RNA is typically depleted prior to sequencing by hybridisation-based methods; an alternative approach is to deplete reads arising from Hgb RNA bioinformatically. In the present study, we compared the impact of these two approaches on the outcome of differential gene expression analysis performed using RNA-seq data from 58 human tuberculosis (TB) patient or contact whole blood samples-29 globin kit-depleted and 29 matched non-depleted-a subset of which were taken at TB diagnosis and at six months post-TB treatment from the same patient. Bioinformatic depletion of Hgb genes from the non-depleted samples (bioinformatic-depleted) substantially reduced library sizes (median = 57.24%) and fewer long non-coding, micro, small nuclear and small nucleolar RNAs were captured in these libraries. Profiling published TB gene signatures across all samples revealed inferior correlation between kit-depleted and bioinformatic-depleted pairs when the proportion of reads mapping to Hgb genes was higher in the non-depleted sample, particularly at the TB diagnosis time point. A set of putative "globin-fingerprint" genes were identified by directly comparing kit-depleted and bioinformatic-depleted samples at each timepoint. Two TB treatment response signatures were also shown to have decreased differential performance when comparing samples at TB diagnosis to six months post-TB treatment when profiled on the bioinformatic-depleted samples compared with their kit-depleted counterparts. These results demonstrate that failure to deplete Hgb RNA prior to sequencing has a negative impact on the sensitivity to detect disease-relevant gene expression changes even when bioinformatic removal is performed.

Molecular epidemiology of SARS-CoV-2 in Northern South Africa: wastewater surveillance from January 2021 to May 2022.

Introduction: Wastewater-based genomic surveillance of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) provides a comprehensive approach to characterize evolutionary patterns and distribution of viral types in a population. This study documents the molecular epidemiology of SARS-CoV-2, in Northern South Africa, from January 2021 to May 2022. Methodology: A total of 487 wastewater samples were collected from the influent of eight wastewater treatment facilities and tested for SARS-CoV-2 RNA using quantitative reverse transcriptase polymerase chain reaction (qRT-PCR). SARS-CoV-2 positive samples with genome copies/mL ≥1,500 were subjected to allele-specific genotyping (ASG) targeting the Spike protein; 75 SARS-CoV-2 positive samples were subjected to whole genome sequencing (WGS) on the ATOPlex platform. Variants of concern (VoC) and lineages were assigned using the Nextclade and PangoLIN Software. Concordance for VoC between ASG and WGS analyses was determined. Sequence relationship was determined by phylogenetic analysis. Results: Seventy-five percent (365/487) of the influent samples were positive for SARS-CoV-2 RNA. Delta and Omicron VoC were more predominant at a prevalence of 45 and 32%, respectively, and they were detected as early as January and February 2021, while Beta VoC was least detected at a prevalence of 5%. A total of 11/60 (18%) sequences were assigned lineages and clades only, but not a specific VoC name. Phylogenetic analysis was used to investigate the relationship of these sequences to other study sequences, and further characterize them. Concordance in variant assignment between ASG and WGS was seen in 51.2% of the study sequences. There was more intra-variant diversity among Beta VoC sequences; mutation E484K was absent. Three previously undescribed mutations (A361S, V327I, D427Y) were seen in Delta VoC. Discussion and Conclusion: The detection of Delta and Omicron VoCs in study sites earlier in the outbreak than has been reported in other regions of South Africa highlights the importance of population-based approaches over individual sample-based approaches in genomic surveillance. Inclusion of non-Spike protein targets could improve the specificity of ASG, since all VoCs share similar Spike protein mutations. Finally, continuous molecular epidemiology with the application of sensitive technologies such as next generation sequencing (NGS) is necessary for the documentation of mutations whose implications when further investigated could enhance diagnostics, and vaccine development efforts.