Cryptococcus neoformans ( C. neoformans var. grubii ) is an environmentally acquired pathogen causing 181,000 HIV-associated deaths each year. We sequenced 699 isolates, primarily C. neoformans from HIV-infected patients, from 5 countries in Asia and Africa. The phylogeny of C. neoformans reveals a recent exponential population expansion, consistent with the increase in the number of susceptible hosts. In our study population, this expansion has been driven by three sub-clades of the C. neoformans VNIa lineage; VNIa-4, VNIa-5 and VNIa-93. These three sub-clades account for 91% of clinical isolates sequenced in our study. Combining the genome data with clinical information, we find that the VNIa-93 sub-clade, the most common sub-clade in Uganda and Malawi, was associated with better outcomes than VNIa-4 and VNIa-5, which predominate in Southeast Asia. This study lays the foundation for further work investigating the dominance of VNIa-4, VNIa-5 and VNIa-93 and the association between lineage and clinical phenotype.

The vast burden of cryptococcal meningitis occurs in immunosuppressed patients, driven by HIV, and is caused by Cryptococcus neoformans var. grubii. We previously reported cryptococcal meningitis in Vietnam arising atypically in HIV uninfected, apparently immunocompetent patients, caused by a single amplified fragment length polymorphism (AFLP) cluster of C. neoformans var. grubii (VNIγ). This variant was less common in HIV infected individuals; it remains unclear why this lineage is associated with apparently immunocompetent patients. To study this host tropism we aimed to further our understanding of clinical phenotype and genomic variation within Vietnamese C. neoformans var. grubii. After performing MLST on C. neoformans clinical isolates we identified 14 sequence types (STs); ST5 correlated with the VNIγ cluster. We next compared clinical phenotype by lineage and found HIV infected patients with cryptococcal meningitis caused by ST5 organisms were significantly more likely to have lymphadenopathy (11% vs. 4%, p = 0.05 Fisher’s exact test) and higher blood lymphocyte count (median 0.76 versus 0.55 X109 cells/L, p = 0.001, Kruskal-Wallis test). Furthermore, survivors of ST5 infections had evidence of worse disability outcomes at 70 days (72.7% (40/55) in ST5 infections versus 57.1% (52/91) non-ST5 infections (OR 2.11, 95%CI 1.01 to 4.41), p = 0.046). To further investigate the relationship between strain and disease phenotype we performed genome sequencing on eight Vietnamese C. neoformans var. grubii. Eight genome assemblies exhibited >99% nucleotide sequence identity and we identified 165 kbp of lineage specific to Vietnamese isolates. ST5 genomes harbored several strain specific regions, incorporating 19 annotated coding sequences and eight hypothetical proteins. These regions included a phenolic acid decarboxylase, a DEAD-box ATP-dependent RNA helicase 26, oxoprolinases, a taurine catabolism dioxygenase, a zinc finger protein, membrane transport proteins and various drug transporters. Our work outlines the complexity of genomic pathogenicity in cryptococcal infections and identifies a number of gene candidates that may aid the disaggregation of the pathways associated with the pathogenesis of Cryptococcus neoformans var. grubii.

Cryptococcosis causes approximately 180 000 deaths each year in patients with human immunodeficiency virus (HIV). Patients with other forms of immunosuppression are also at risk, and disease is increasingly recognized in apparently immunocompetent individuals. Cryptococcus neoformans var. grubii , responsible for the majority of cases, is distributed globally. We used the consensus ISHAM Multilocus sequence typing (MLST) scheme to define the population structure of clinical C. neoformans var. grubii isolates from Laos ( n = 81), which we placed into the global context using published MLST data from other countries (total N = 1047), including a reanalysis of 136 Vietnamese isolates previously reported. We observed a phylogeographical relationship in which the Laotian population was similar to its neighbor Thailand, being dominated (83%) by Sequence Types (ST) 4 and 6. This phylogeographical structure changed moving eastwards, with Vietnam's population consisting of an admixture of isolates dominated by the ST4/ST6 (35%) and ST5 (48%) lineages. The ST5 lineage is the predominant ST reported from China and East Asia, where it accounts for >90% of isolates. Analysis of genetic distance ( Fst ) between different populations of C. neoformans var. grubii supports this intermediate structure of the Vietnamese population. The pathogen and host diversity reported from Vietnam provide the strongest epidemiological evidence of the association between ST5 and HIV-uninfected patients. Regional anthropological genetic distances suggest diversity in the C. neoformans var. grubii population across Southeast Asia is driven by ecological rather than human host factors. Where the ST5 lineage is present, disease in HIV-uninfected patients is to be expected.

We previously observed a substantial burden of cryptococcal meningitis in Vietnam atypically arising in individuals who are uninfected with human immunodeficiency virus (HIV). This disease was associated with a single genotype of Cryptococcus neoformans (sequence type [ST]5), which was significantly less common in HIV-infected individuals. Aiming to compare the phenotypic characteristics of ST5 and non-ST5 C. neoformans, we selected 30 representative Vietnamese isolates and compared their in vitro pathogenic potential and in vivo virulence. ST5 and non-ST5 organisms exhibited comparable characteristics with respect to in vitro virulence markers including melanin production, replication at 37°C, and growth in cerebrospinal fluid. However, the ST5 isolates had significantly increased variability in cellular and capsular sizing compared with non-ST5 organisms (P < .001). Counterintuitively, mice infected with ST5 isolates had significantly longer survival with lower fungal burdens at day 7 than non-ST5 isolates. Notably, ST5 isolates induced significantly greater initial inflammatory responses than non-ST5 strains, measured by TNF-α concentrations (P < .001). Despite being generally less virulent in the mouse model, we hypothesize that the significant within strain variation seen in ST5 isolates in the tested phenotypes may represent an evolutionary advantage enabling adaptation to novel niches including apparently immunocompetent human hosts.

27C. neoformans var. grubii (C. neoformans) is an environmentally acquired pathogen causing 181 000 28 HIV-associated deaths each year. We used whole genome sequencing (WGS) to characterise 699 29 isolates, primarily C. neoformans from HIV-infected patients, from 5 countries in Asia and Africa. We 30 found that 91% of our clinical isolates belonged to one of three highly clonal sub-clades of VNIa, 31 which we have termed VNIa-4, VNIa-5 and VNIa-93. Parsimony analysis revealed frequent, long 32 distance transmissions of C. neoformans; international transmissions took place on 13% of VNIa-4 33 branches, and intercontinental transmissions on 7% of VNIa-93 branches. The median length of 34 within sub-clade internal branches was 3-6 SNPs, while terminal branches were 44.5-77.5 SNPs. The 35 short median internal branches were partly driven by the large number (12-15% of internal 36 branches) of polytomies in the within-sub-clade trees. To simultaneously explain our observation of 37 no apparent molecular clock, short internal branches and frequent polytomies we hypothesise that 38 C. neoformans VNIa spends much of its time in the environment in a quiescent state, while, when it 39 is sampled, it has almost always undergone an extended period of growth. Infections with VNIa-93 40 were associated with a significantly reduced risk of death by 10 weeks compared with infections 41 with VNIa-4 (Hazard Ratio = 0.45, p = 0.003). We detected a recombination in the mitochondrial 42 sequence of VNIa-5, suggesting that mitochondria could be involved in the propensity of this sub-43 clade to infect HIV-uninfected patients. These data highlight the insight into the biology and 44 epidemiology of pathogenic fungi which can be gained from WGS data. 45 46 3