Activated dormant stem cells recover spermatogenesis in chemoradiotherapy-induced infertility.

Male infertility is a recognized side effect of chemoradiotherapy. Extant spermatogonial stem cells (SSCs) may act as originators for any subsequent recovery. However, which type of SSCs, the mechanism by which they survive and resist toxicity, and how they act to restart spermatogenesis remain largely unknown. Here, we identify a small population of Set domain-containing protein 4 (Setd4)-expressing SSCs that occur in a relatively dormant state in the mouse seminiferous tubule. Extant beyond high-dose chemoradiotherapy, these cells then activate to recover spermatogenesis. Recovery fails when Setd4+ SSCs are deleted. Confirmed to be of fetal origin, these Setd4+ SSCs are shown to facilitate early testicular development and also contribute to steady-state spermatogenesis in adulthood. Upon activation, chromatin remodeling increases their genome-wide accessibility, enabling Notch1 and Aurora activation with corresponding silencing of p21 and p53. Here, Setd4+ SSCs are presented as the originators of both testicular development and spermatogenesis recovery in chemoradiotherapy-induced infertility.

Automated Cell Lineage Reconstruction using Label-Free 4D Microscopy.

Patterns of lineal descent play a critical role in the development of metazoan embryos. In eutelic organisms that generate a fixed number of somatic cells, invariance in the topology of their cell lineage provides a powerful opportunity to interrogate developmental events with empirical repeatability across individuals. Studies of embryonic development using the nematode Caenorhabditis elegans have been drivers of discovery. These studies have depended heavily on high throughput lineage tracing enabled by 4D fluorescence microscopy and robust computer vision pipelines. For a range of applications, computer-aided yet manual lineage tracing using 4D label-free microscopy remains an essential tool. Deep learning approaches to cell detection and tracking in fluorescence microscopy have advanced significantly in recent years, yet solutions for automating cell detection and tracking in 3D label-free imaging of dense tissues and embryos remain inaccessible. Here we describe embGAN, a deep learning pipeline that addresses the challenge of automated cell detection and tracking in label-free 3D time lapse imaging. embGAN requires no manual data annotation for training, learns robust detections that exhibits a high degree of scale invariance and generalizes well to images acquired in multiple labs on multiple instruments. We characterize embGAN's performance using lineage tracing in the C. elegans embryo as a benchmark. embGAN achieves near state-of-the-art performance in cell detection and tracking, enabling high-throughput studies of cell lineage without the need for fluorescent reporters or transgenics.

Overcoming post-transplant graft failure and adenovirus infection in a patient with FLT3-TKD-mutated mixed-phenotype acute leukemia: A case report.

Mixed-phenotype acute leukemia (MPAL) with FLT3-TKD mutations is a rare and challenging subtype of leukemia. Effective management strategies are crucial for improving patient outcomes. A 31-year-old man with FLT3-TKD-mutated MPAL achieved hematological remission through the JALSG ALL202-O protocol and gilteritinib, followed by cord blood transplantation (CBT). Post-transplant complications included adenovirus-induced hemorrhagic cystitis, managed with bladder irrigation and ribavirin, and engraftment failure, necessitating a second CBT on Day 35. Subsequent adenoviral conjunctivitis resolved with vidarabine. The patient achieved neutrophil engraftment by Day 76 and was discharged on Day 173 without relapse. This case highlights the importance of vigilant supportive care and tailored therapy in managing MPAL with FLT3 mutations, especially in the context of post-transplant complications.

Chimerism and mosaicism shape our physical constitution and impact medical conditions.

ABO blood group discrepancies in healthy individuals were caused by body-wide chimerism and mosaicism. They can be evaluated with new diagnostic options for disease-related cell clones that are typically associated with mosaicism. The observations raise the attention for sporadic mixed-field observations of any blood group antigen. Commentary on: Dauber et al. Body-wide chimerism and mosaicism are predominant causes of naturally occurring ABO discrepancies. Br J Haematol 2024 (Online ahead of print). doi:10.1111/bjh.19618.

Phylogenomics and pervasive genome-wide phylogenetic discordance among fin whales (Balaenoptera physalus).

Phylogenomics has the power to uncover complex phylogenetic scenarios across the genome. In most cases, no single topology is reflected across the entire genome as the phylogenetic signal differs among genomic regions due to processes, such as introgression and incomplete lineage sorting. Baleen whales are among the largest vertebrates on Earth with a high dispersal potential in a relatively unrestricted habitat, the oceans. The fin whale (Balaenoptera physalus) is one of the most enigmatic baleen whale species, currently divided into four subspecies. It has been a matter of debate whether phylogeographic patterns explain taxonomic variation in fin whales. Here we present a chromosome-level whole genome analysis of the phylogenetic relationships among fin whales from multiple ocean basins. First, we estimated concatenated and consensus phylogenies for both the mitochondrial and nuclear genomes. The consensus phylogenies based upon the autosomal genome uncovered monophyletic clades associated with each ocean basin, aligning with the current understanding of subspecies division. Nevertheless, discordances were detected in the phylogenies based on the Y chromosome, mitochondrial genome, autosomal genome and X chromosome. Furthermore, we detected signs of introgression and pervasive phylogenetic discordance across the autosomal genome. This complex phylogenetic scenario could be explained by a puzzle of introgressive events, not yet documented in fin whales. Similarly, incomplete lineage sorting and low phylogenetic signal could lead to such phylogenetic discordances. Our study reinforces the pitfalls of relying on concatenated or single locus phylogenies to determine taxonomic relationships below the species level by illustrating the underlying nuances which some phylogenetic approaches may fail to capture. We emphasize the significance of accurate taxonomic delineation in fin whales by exploring crucial information revealed through genome-wide assessments.

Identification of Meibomian gland stem cell populations and mechanisms of aging.

Meibomian glands secrete lipid-rich meibum, which prevents tear evaporation. Aging-related Meibomian gland shrinkage may result in part from stem cell exhaustion and is associated with evaporative dry eye disease, a common condition lacking effective treatment. The identities and niche of Meibomian gland stem cells and the signals controlling their activity are poorly defined. Using snRNA-seq, in vivo lineage tracing, ex vivo live imaging, and genetic studies in mice, we identified markers for stem cell populations that maintain distinct regions of the gland and uncovered Hh signaling as a key regulator of stem cell proliferation. Consistent with this, human Meibomian gland carcinoma exhibited increased Hh signaling. Aged glands displayed decreased Hh and EGF signaling, deficient innervation, and loss of collagen I in niche fibroblasts, indicating that alterations in both glandular epithelial cells and their surrounding microenvironment contribute to age-related degeneration. These findings suggest new approaches to treat aging-associated Meibomian gland loss.

The complete mitochondrial genome of Heptathela kimurai (Araneae: Heptathelidae).

Heptathela kimurai (Kishida, 1920) is a spider that belongs to the family Heptathelidae which is a basial lineage of spiders. The molecular information of ancestral species belonging to families like Heptathelidae is comparatively limited when compared to spider species from derived families. Here we present the complete mitochondrial genome sequence (mtDNA) of H. kimurai. The sequence was obtained using massively parallel sequencing technology. The circular genome was 14,224 bp in length, and the AT content was 69.53%. The H. kimurai mitochondrial genome contains 13 protein-coding genes (PCGs), 21 tRNA genes, and 2 rRNA genes. The majority of PCGs were found in the heavy strand.

Baf155 controls hematopoietic differentiation and regeneration through chromatin priming.

Chromatin priming promotes cell-type-specific gene expression, lineage differentiation, and development. The mechanism of chromatin priming has not been fully understood. Here, we report that mouse hematopoietic stem and progenitor cells (HSPCs) lacking the Baf155 subunit of the BAF (BRG1/BRM-associated factor) chromatin remodeling complex produce a significantly reduced number of mature blood cells, leading to a failure of hematopoietic regeneration upon transplantation and 5-fluorouracil (5-FU) injury. Baf155-deficient HSPCs generate particularly fewer neutrophils, B cells, and CD8+ T cells at homeostasis, supporting a more immune-suppressive tumor microenvironment and enhanced tumor growth. Single-nucleus multiomics analysis reveals that Baf155-deficient HSPCs fail to establish accessible chromatin in selected regions that are enriched for putative enhancers and binding motifs of hematopoietic lineage transcription factors. Our study provides a fundamental mechanistic understanding of the role of Baf155 in hematopoietic lineage chromatin priming and the functional consequences of Baf155 deficiency in regeneration and tumor immunity.

Hepatic stellate cells in zone 1 engage in capillarization rather than myofibroblast formation in murine liver fibrosis.

The combination of lineage tracing and immunohistochemistry has helped to identify subpopulations and fate of hepatic stellate cells (HSC) in murine liver. HSC are sinusoidal pericytes that act as myofibroblast precursors after liver injury. Single cell RNA sequencing approaches have recently helped to differentiate central and portal HSC. A specific Cre line to lineage trace portal HSC has not yet been described. We used three Cre lines (Lrat-Cre, PDGFRß-CreERT2 and SMMHC-CreERT2) known to label mesenchymal cells including HSC in combination with a tdTomato-expressing reporter. All three Cre lines labeled populations of HSC as well as smooth muscle cells (SMC). Using the SMMHC-CreERT2, we identified a subtype of HSC in the periportal area of the hepatic lobule (termed zone 1-HSC). We lineage traced tdTomato-expressing zone 1-HSC over 1 year, described fibrotic behavior in two fibrosis models and investigated their possible role during fibrosis. This HSC subtype resides in zone 1 under healthy conditions; however, zonation is disrupted in preclinical models of liver fibrosis (CCl4 and MASH). Zone 1-HSC do not transform into αSMA-expressing myofibroblasts. Rather, they participate in sinusoidal capillarization. We describe a novel subtype of HSC restricted to zone 1 under physiological conditions and its possible function after liver injury. In contrast to the accepted notion, this HSC subtype does not transform into αSMA-positive myofibroblasts; rather, zone 1-HSC adopt properties of capillary pericytes, thereby participating in sinusoidal capillarization.

Anatomy, Histology, and Embryonic Origin of Adipose Tissue: Insights to Understand Adipose Tissue Homofunctionality in Regeneration and Therapies.

Preadipocytes are formed during the 14th and 16th weeks of gestation. White adipose tissue, in particular, is generated in specific areas and thereby assembles after birth, rapidly increasing following the propagation of adipoblasts, which are considered the preadipocyte cell precursors. The second trimester of gestation is a fundamental phase of adipogenesis, and in the third trimester, adipocytes, albeit small may be present within the main deposition areas. In the course of late gestation, adipose tissue develops in the foetus and promotes the synthesis of large amounts of uncoupling protein 1, in similar quantities relative to differentiated brown adipose tissue. In mammals, differentiation occurs in two functionally different types of adipose cells: white adipose cells resulting from lipid storage and brown adipose cells from increased metabolic energy consumption. During skeletogenesis, synovial joints develop through the condensation of mesenchymal cells, which forms an insertional layer of flattened cells that umlaut skeletal elements, by sharing the same origin in the development of synovium. Peri-articular fat pads possess structural similarity with body subcutaneous white adipose tissue; however, they exhibit a distinct metabolic function due to the micro-environmental cues in which they are embedded. Fat pads are an important component of the synovial joint and play a key role in the maintenance of joint homeostasis. They are also implicated in pathological states such as osteoarthritis.In this paper we explore the therapeutic potential of adipocyte tissue mesenchymal precursor-based stem cell therapy linking it back to the anatomic origin of adipose tissue.

Targeted next-generation sequencing of Mycobacterium tuberculosis from patient samples: lessons learned from high drug-resistant burden clinical settings in Bangladesh.

AbstractLack of appropriate early diagnostic tools for drug-resistant tuberculosis (DR-TB) and their incomplete drug susceptibility testing (DST) profiling is concerning for TB disease control. Existing methods, such as phenotypic DST (pDST), are time-consuming, while Xpert MTB/RIF (Xpert) and line probe assay (LPA) are limited to detecting resistance to few drugs. Targeted next-generation sequencing (tNGS) has been recently approved by WHO as an alternative approach for rapid and comprehensive DST. We aimed to investigate the performance and feasibility of tNGS for detecting DR-TB directly from clinical samples in Bangladesh. pDST, LPA and tNGS were performed among 264 sputum samples, either rifampicin-resistant (RR) or rifampicin-sensitive (RS) TB cases confirmed by Xpert assay. Resistotypes of tNGS were compared with pDST, LPA and composite reference standard (CRS, resistant if either pDST or LPA showed a resistant result). tNGS results revealed higher sensitivities for rifampicin (RIF) (99.3%), isoniazid (INH) (96.3%), fluoroquinolones (FQs) (94.4%), and aminoglycosides (AMGs) (100%) but comparatively lower for ethambutol (76.6%), streptomycin (68.7%), ethionamide (56.0%) and pyrazinamide (50.7%) when compared with pDST. The sensitivities of tNGS for INH, RIF, FQs and AMGs were 93.0%, 96.6%, 90.9%, and 100%, respectively and the specificities ranged from 91.3% to 100% when compared with CRS. This proof of concept study, conducted in a high-burden setting demonstrated that tNGS is a valuable tool for identifying DR-TB directly from the clinical specimens. Its feasibility in our laboratory suggests potential implementation and moving tNGS from research settings into clinical settings.

The cell fates of intermediate cell population in prostate development.

Organ development, regeneration and cancer initiation are typically influenced by the proliferation and lineage plasticity of tissue-specific stem cells. Prostate intermediate cells, which exhibit characteristics of both basal and luminal cells, are prevalent in pathological states and during organ development. However, the identity, fate and function of these intermediate cells in prostate development are not well understood. Through single-cell RNA-seq analysis on neonatal urogenital sinus tissue, we identified intermediate cells exhibiting stem cell potential. A notable decline in the population of intermediate cells was observed during prostate development. Prostate intermediate cells were specifically labeled in early and late postnatal development by the enhanced dual-recombinase-mediated genetic tracing systems. Our findings revealed that these cells possess significant stem cell capabilities as demonstrated in organoid formation and cell fate mapping assays. These intermediate cells also exhibited intrinsic bipotential properties, enabling them to differentiate into both basal and luminal cells. Additionally, we discovered a novel transition from intermediate cell expressing neuroendocrine markers to neuroendocrine cell during prostate development. This study highlights intermediate cells as a crucial stem cell population and enhances our understanding of their role in prostate development and the plasticity of prostate cancer lineage.

Structural changes shaping the Drosophila ellipsoid body ER-neurons during development and aging.

The ellipsoid body (EB) of the insect brain performs pivotal functions in controlling navigation. Input and output of the EB is provided by multiple classes of R-neurons (now referred to as ER-neurons) and columnar neurons which interact with each other in a stereotypical and spatially highly ordered manner. The developmental mechanisms that control the connectivity and topography of EB neurons are largely unknown. One indispensable prerequisite to unravel these mechanisms is to document in detail the sequence of events that shape EB neurons during their development. In this study, we analyzed the development of the Drosophila EB. In addition to globally following the ER-neuron and columnar neuron (sub)classes in the spatial context of their changing environment we performed a single cell analysis using the multi-color flip out (MCFO) system to analyze the developmental trajectory of ER-neurons at different pupal stages, young adults (4d) and aged adults (∼60d). We show that the EB develops as a merger of two distinct elements, a posterior and anterior EB primordium (prEBp and prEBa, respectively. ER-neurons belonging to different subclasses form growth cones and filopodia that associate with the prEBp and prEBa in a pattern that, from early pupal stages onward, foreshadows their mature structure. Filopodia of all ER-subclasses are initially much longer than the dendritic and terminal axonal branches they give rise to, and are pruned back during late pupal stages. Interestingly, extraneous branches, particularly significant in the dendritic domain, are a hallmark of ER-neuron structure in aged brains. Aging is also associated with a decline in synaptic connectivity from columnar neurons, as well as upregulation of presynaptic protein (Brp) in ER-neurons. Our findings advance the EB (and ER-neurons) as a favorable system to visualize and quantify the development and age-related decline of a complex neuronal circuitry.

Prognostic value of a novel myeloid-to-lymphoid ratio biomarker in advanced gastric cancer.

BACKGROUND: Currently, immune checkpoint inhibitors (ICIs) have excellent performance in the clinical treatment of advanced gastric cancer (AGC). However, precisely selecting AGC patients who can benefit from immunotherapy is an urgent difficulty. In this study, we investigated the immunoprognostic role of myeloid-to-lymphocyte ratio (M:L) in AGC patients. METHODS: We collected information on 268 AGC patients who were hospitalized in the Department of Medical Oncology of PLA General Hospital from December 2014 to May 2021. The patients were divided into low M: L group (< 3.76) and high M:L group (≥ 3.76). Survival differences between different M: L level groups at baseline and after treatment were analyzed by methods such as Kaplan-Meier, Cox or Logistic regression model. RESULTS: Progression free survival (PFS) (5.8 months vs. 3.4 months, p = 0.001) and overall survival (OS) (14.1 months vs. 9.0 months, p = 0.001) were significantly longer in the low M:L group than in the high M:L group. After analyses of Cox regression modeling it was concluded that M:L was an independent prognostic factor for PFS (HR 1.371 95%CI 1.057-1.777 p = 0.017) and OS (HR 1.352 95%CI 1.003-1.824 p = 0.048), respectively. Subsequent subgroup analyses performed across immunotherapy lines, regimens, PD-1 inhibitor agents, and age groups revealed a poorer prognosis in the high M:L group. Notably, an increase in the value of M:L after treatment significantly increased the risk of poor prognosis. CONCLUSIONS: M:L ≥ 3.76 is associated with poor prognostic outcomes in AGC patients receiving immunotherapy and may be a predictive biomarker of prognosis. This result needs to be confirmed by larger prospective studies.

Pan-lineage Mycobacterium tuberculosis reference genome for enhanced molecular diagnosis.

In Mycobacterium tuberculosis (MTB) control, whole genome sequencing-based molecular drug susceptibility testing (molDST-WGS) has emerged as a pivotal tool. However, the current reliance on a single-strain reference limits molDST-WGS's true potential. To address this, we introduce a new pan-lineage reference genome, 'MtbRf'. We assembled 'unmapped' reads from 3,614 MTB genomes (751 L1; 881 L2; 1,700 L3; and 282 L4) into 35 shared, annotated contigs (54 coding sequences [CDSs]). We constructed MtbRf through: (1) searching for contig homologues among genome database that precipitate results uniquely within Mycobacteria genus; (2) comparing genomes with H37Rv ('lift-over') to define 18 insertions; and (3) filling gaps in H37Rv with insertions. MtbRf adds 1.18% sequences to H37rv, salvaging >60% of previously unmapped reads. Transcriptomics confirmed gene expression of new CDSs. The new variants provided a moderate DST predictive value (AUROC 0.60-0.75). MtbRf thus unveils previously hidden genomic information and lays the foundation for lineage-specific molDST-WGS.

Dynamic regulation of tissue fluidity controls skin repair during wound healing.

During wound healing, different pools of stem cells (SCs) contribute to skin repair. However, how SCs become activated and drive the tissue remodeling essential for skin repair is still poorly understood. Here, by developing a mouse model allowing lineage tracing and basal cell lineage ablation, we monitor SC fate and tissue dynamics during regeneration using confocal and intravital imaging. Analysis of basal cell rearrangements shows dynamic transitions from a solid-like homeostatic state to a fluid-like state allowing tissue remodeling during repair, as predicted by a minimal mathematical modeling of the spatiotemporal dynamics and fate behavior of basal cells. The basal cell layer progressively returns to a solid-like state with re-epithelialization. Bulk, single-cell RNA, and epigenetic profiling of SCs, together with functional experiments, uncover a common regenerative state regulated by the EGFR/AP1 axis activated during tissue fluidization that is essential for skin SC activation and tissue repair.

LineageVAE: Reconstructing Historical Cell States and Transcriptomes toward Unobserved Progenitors.

MOTIVATION: Single-cell RNA sequencing (scRNA-seq) enables comprehensive characterization of the cell state. However, its destructive nature prohibits measuring gene expression changes during dynamic processes such as embryogenesis. Although recent studies integrating scRNA-seq with lineage tracing have provided clonal insights between progenitor and mature cells, challenges remain. Because of their experimental nature, observations are sparse, and cells observed in the early state are not the exact progenitors of cells observed at later time points. To overcome these limitations, we developed LineageVAE, a novel computational methodology that utilizes deep learning based on the property that cells sharing barcodes have identical progenitors. RESULTS: LineageVAE is a deep generative model that transforms scRNA-seq observations with identical lineage barcodes into sequential trajectories toward a common progenitor in a latent cell state space. This method enables the reconstruction of unobservable cell state transitions, historical transcriptomes, and regulatory dynamics at a single-cell resolution. Applied to hematopoiesis and reprogrammed fibroblast datasets, LineageVAE demonstrated its ability to restore backward cell state transitions and infer progenitor heterogeneity and transcription factor activity along differentiation trajectories. AVAILABILITY AND IMPLEMENTATION: The LineageVAE model was implemented in Python using the PyTorch deep learning library. The code is available on GitHub at https://github.com/LzrRacer/LineageVAE/. SUPPLEMENTARY INFORMATION: Available at Bioinformatics online.

Comparative analysis of the RVA VP7 and VP4 antigenic epitopes circulating in Iran and the Rotarix and RotaTeq vaccines.

Analyzing the lineages and detecting antigenic variation in immunogenic motifs of Group A Rotavirus (RVA) variants is crucial because it can impact vaccine efficacy. This study investigated the circulating lineages of VP4 and VP7 proteins of human RVA isolates and their phylogeny in ≤24-month-old symptomatic, rotavirus-positive children with transudative diarrhea within 48 h of admission to Mofid Children's Hospital between December 2020 and March 2022 in Tehran, Iran. Antigen detection was performed by ELISA, RNA extraction, and semi-nested multiplex PCR for G/P genotypes, followed by sequencing and bioinformatic analysis using multiple sequence alignments in MEGA and phylogenetic analysis by Geneious Prime. The similarity of VP7 and VP4 amino acids with the RotaTeq and Rotarix vaccine strains for cytotoxic T cell and antigenic epitopes was evaluated using the UCSF Chimera Molecular Modeling System. Overall, 27.3 % of the samples were RVA positive, showing untypeable (2.5 %), single (76.9 %), and mixed (20.5 %) genotypic characteristics. The strains clustered in the G1/II, G2/IV, G3/I, G4/I, G9/III, P (Kachooei et al., 2023) [8]/III, P (Howley et al., 2020) [4]/V, and P (Wahyuni et al., 2021) [6]/I lineages. Comparative analysis of VP7 antigenic epitopes showed that the G1/II strains were completely conserved, while the G2/IV, G3/I, G4/I, G6, G9/III strains contained 2, 3-5, 2, 4 and 9 amino acid substitutions, respectively. The P (Kachooei et al., 2023) [8]/III genotypes differed by 3 amino acids, while the P (Wahyuni et al., 2021) [6]/I genotype had the most substitutions. CTL epitopes were completely conserved in G3/I strains, but other genotypes differed by 1-4 amino acids compared to the vaccine strains. Given the diversity of circulating RVA genotypes and the observed mutations in neutralizing and CTL epitopes, immune escape by some of the strains is likely in Iran. This finding underscores the importance of evaluating the effect of rotavirus vaccines on local genotypes and related lineages before implementing a vaccination program.

H3K27me3-mediated epigenetic regulation in pluripotency maintenance and lineage differentiation.

Cell fate determination is an intricate process which is orchestrated by multiple regulatory layers including signal pathways, transcriptional factors, epigenetic modifications, and metabolic rewiring. Among the sophisticated epigenetic modulations, the repressive mark H3K27me3, deposited by PRC2 (polycomb repressive complex 2) and removed by demethylase KDM6, plays a pivotal role in mediating the cellular identity transition through its dynamic and precise alterations. Herein, we overview and discuss how H3K27me3 and its modifiers regulate pluripotency maintenance and early lineage differentiation. We primarily highlight the following four aspects: 1) the two subcomplexes PRC2.1 and PRC2.2 and the distribution of genomic H3K27 methylation; 2) PRC2 as a critical regulator in pluripotency maintenance and exit; 3) the emerging role of the eraser KDM6 in early differentiation; 4) newly identified additional factors influencing H3K27me3. We present a comprehensive insight into the molecular principles of the dynamic regulation of H3K27me3, as well as how this epigenetic mark participates in pluripotent stem cell-centered cell fate determination.

New lineages of RNA viruses from clinical isolates of Rhizopus microsporus revealed by fragmented and primer-ligated dsRNA sequencing (FLDS) analysis.

Rhizopus microsporus is a species in the order Mucorales that is known to cause mucormycosis, but it is poorly understood as a host of viruses. Here, we examined 25 clinical strains of R. microsporus for viral infection with a conventional double-stranded RNA (dsRNA) assay using agarose gel electrophoresis (AGE) and the recently established fragmented and primer-ligated dsRNA sequencing (FLDS) protocol. By AGE, five virus-infected strains were detected. Then, full-length genomic sequences of 12 novel RNA viruses were revealed by FLDS, which were related to the families Mitoviridae, Narnaviridae, and Endornaviridae, ill-defined groups of single-stranded RNA (ssRNA) viruses with similarity to the established families Virgaviridae and Phasmaviridae, and the proposed family "Ambiguiviridae." All the characterized viruses, except a potential phasmavirid with a negative-sense RNA genome, had positive-sense RNA genomes. One virus belonged to a previously established species within the family Mitoviridae, whereas the other 11 viruses represented new species or even new genera. These results show that the fungal pathogen R. microsporus harbors diverse RNA viruses and extend our understanding of the diversity of RNA viruses in the fungal order Mucorales, division Mucoromycota. Identifying RNA viruses from clinical isolates of R. microsporus may expand the repertoire of natural therapeutic agents for mucormycosis in the future.IMPORTANCEThe diversity of mycoviruses in fungal hosts in the division Mucoromycota has been underestimated, mainly within the species Rhizopus microsporus. Only five positive-sense RNA genomes had previously been discovered in this species. Because current sequencing methods poorly complete the termini of genomes, we used fragmented and primer-ligated double-stranded RNA sequencing to acquire the full-length genomes. Eleven novel mycoviruses were detected in this study, including the first negative-sense RNA genome reported in R. microsporus. Our findings extend the understanding of the viral diversity in clinical strains of Mucoromycota, may provide insights into the pathogenesis and ecology of this fungus, and may offer therapeutic options.