Whole genome analysis of Rhizopus species causing rhino-cerebral mucormycosis during the COVID-19 pandemic

Introduction

Mucormycosis is an acute invasive fungal disease (IFD) seen mainly in immunocompromised hosts and in patients with uncontrolled diabetes. The incidence of mucormycosis increased exponentially in India during the SARS-CoV-2 (henceforth COVID-19) pandemic. Since there was a lack of data on molecular epidemiology of Mucorales causing IFD during and after the COVID-19 pandemic, whole genome analysis of the Rhizopus spp. isolated during this period was studied along with the detection of mutations that are associated with antifungal drug resistance.

Materials and methods

A total of 50 isolates of Rhizopus spp. were included in this prospective study, which included 28 from patients with active COVID-19 disease, 9 from patients during the recovery phase, and 13 isolates from COVID-19-negative patients. Whole genome sequencing (WGS) was performed for the isolates, and the de novo assembly was done with the Spades assembler. Species identification was done by extracting the ITS gene sequence from each isolate followed by searching Nucleotide BLAST. The phylogenetic trees were made with extracted ITS gene sequences and 12 eukaryotic core marker gene sequences, respectively, to assess the genetic distance between our isolates. Mutations associated with intrinsic drug resistance to fluconazole and voriconazole were analyzed.

Results

All 50 patients presented to the hospital with acute fungal rhinosinusitis. These patients had a mean HbA1c of 11.2%, and a serum ferritin of 546.8 ng/mL. Twenty-five patients had received steroids. By WGS analysis, 62% of the Rhizopus species were identified as R. delemar. Bayesian analysis of population structure (BAPS) clustering categorized these isolates into five different groups, of which 28 belong to group 3, 9 to group 5, and 8 to group 1. Mutational analysis revealed that in the CYP51A gene, 50% of our isolates had frameshift mutations along with 7 synonymous mutations and 46% had only synonymous mutations, whereas in the CYP51B gene, 68% had only synonymous mutations and 26% did not have any mutations.

Conclusion

WGS analysis of Mucorales identified during and after the COVID-19 pandemic gives insight into the molecular epidemiology of these isolates in our community and establishes newer mechanisms for intrinsic azole resistance.

Performance evaluation of multiplex PCR including Aspergillus-not so simple!

Multiplex PCRs have been designed for including species other than Aspergillus fumigatus for the diagnosis of invasive aspergillosis, such as microarrays, liquid-phase array, and electrospray-ionization mass spectrometry (PCR/ESI MS). These methods are based on the selection of multiple primers to amplify different species with the specificity checked by hybridization to a probe or by base composition of the amplicon for the PCR/ESI MS. When testing complex samples such as respiratory specimens, some clinically relevant species can be missed. Indeed, it is impossible to design primers able to amplify all the known fungal species with the same efficiency. Therefore, the best amplified species may not be the most clinically relevant. Multiplex assays have also been proposed to detect A. fumigatus DNA and azole resistance. Since the gene responsible for azole resistance is single copy and the gene used for detection is multicopy, only the high fungal loads can be evaluated. Thus, although interesting for investigating mycobiome, the multiplex assays should be used with cautious for the diagnosis of IA or the detection of resistance. For the diagnosis of invasive aspergillosis, validated quantitative PCRs specifically targeting A. fumigatus or a limited set of species to increase sensitivity is a safer option.