Deletion of the Candida albicans TLO gene family using CRISPR-Cas9 mutagenesis allows characterisation of functional differences in α-, ß- and γ- TLO gene function.

The Candida albicans genome contains between ten and fifteen distinct TLO genes that all encode a Med2 subunit of Mediator. In order to investigate the biological role of Med2/Tlo in C. albicans we deleted all fourteen TLO genes using CRISPR-Cas9 mutagenesis. ChIP-seq analysis showed that RNAP II localized to 55% fewer genes in the tloΔ mutant strain compared to the parent, while RNA-seq analysis showed that the tloΔ mutant exhibited differential expression of genes required for carbohydrate metabolism, stress responses, white-opaque switching and filamentous growth. Consequently, the tloΔ mutant grows poorly in glucose- and galactose-containing media, is unable to grow as true hyphae, is more sensitive to oxidative stress and is less virulent in the wax worm infection model. Reintegration of genes representative of the α-, ß- and γ-TLO clades resulted in the complementation of the mutant phenotypes, but to different degrees. TLOα1 could restore phenotypes and gene expression patterns similar to wild-type and was the strongest activator of glycolytic and Tye7-regulated gene expression. In contrast, the two γ-TLO genes examined (i.e., TLOγ5 and TLOγ11) had a far lower impact on complementing phenotypic and transcriptomic changes. Uniquely, expression of TLOß2 in the tloΔ mutant stimulated filamentous growth in YEPD medium and this phenotype was enhanced when Tloß2 expression was increased to levels far in excess of Med3. In contrast, expression of reintegrated TLO genes in a tloΔ/med3Δ double mutant background failed to restore any of the phenotypes tested, suggesting that complementation of these Tlo-regulated processes requires a functional Mediator tail module. Together, these data confirm the importance of Med2/Tlo in a wide range of C. albicans cellular activities and demonstrate functional diversity within the gene family which may contribute to the success of this yeast as a coloniser and pathogen of humans.

Human Genetic Determinants of Viral Diseases.

Much progress has been made in the identification of specific human gene variants that contribute to enhanced susceptibility or resistance to viral diseases. Herein we review multiple discoveries made with genome-wide or candidate gene approaches that have revealed significant insights into virus-host interactions. Genetic factors that have been identified include genes encoding virus receptors, receptor-modifying enzymes, and a wide variety of innate and adaptive immunity-related proteins. We discuss a range of pathogenic viruses, including influenza virus, respiratory syncytial virus, human immunodeficiency virus, human T cell leukemia virus, human papilloma virus, hepatitis B and C viruses, herpes simplex virus, norovirus, rotavirus, parvovirus, and Epstein-Barr virus. Understanding the genetic underpinnings that affect infectious disease outcomes should allow tailored treatment and prevention approaches in the future.

Balancing openness with Indigenous data sovereignty: An opportunity to leave no one behind in the journey to sequence all of life.

The field of genomics has benefited greatly from its "openness" approach to data sharing. However, with the increasing volume of sequence information being created and stored and the growing number of international genomics efforts, the equity of openness is under question. The United Nations Convention of Biodiversity aims to develop and adopt a standard policy on access and benefit-sharing for sequence information across signatory parties. This standardization will have profound implications on genomics research, requiring a new definition of open data sharing. The redefinition of openness is not unwarranted, as its limitations have unintentionally introduced barriers of engagement to some, including Indigenous Peoples. This commentary provides an insight into the key challenges of openness faced by the researchers who aspire to protect and conserve global biodiversity, including Indigenous flora and fauna, and presents immediate, practical solutions that, if implemented, will equip the genomics community with both the diversity and inclusivity required to respectfully protect global biodiversity.

Butyrate Reprograms Expression of Specific Interferon-Stimulated Genes.

Butyrate is an abundant metabolite produced by gut microbiota. While butyrate is a known histone deacetylase inhibitor that activates expression of many genes involved in immune system pathways, its effects on virus infections and on the antiviral type I interferon (IFN) response have not been adequately investigated. We found that butyrate increases cellular infection with viruses relevant to human and animal health, including influenza virus, reovirus, HIV-1, human metapneumovirus, and vesicular stomatitis virus. Mechanistically, butyrate suppresses levels of specific antiviral IFN-stimulated gene (ISG) products, such as RIG-I and IFITM3, in human and mouse cells without inhibiting IFN-induced phosphorylation or nuclear translocation of the STAT1 and STAT2 transcription factors. Accordingly, we discovered that although butyrate globally increases baseline expression of more than 800 cellular genes, it strongly represses IFN-induced expression of 60% of ISGs and upregulates 3% of ISGs. Our findings reveal that there are differences in the IFN responsiveness of major subsets of ISGs depending on the presence of butyrate in the cell environment, and overall, they identify a new mechanism by which butyrate influences virus infection of cells.IMPORTANCE Butyrate is a lipid produced by intestinal bacteria. Here, we newly show that butyrate reprograms the innate antiviral immune response mediated by type I interferons (IFNs). Many of the antiviral genes induced by type I IFNs are repressed in the presence of butyrate, resulting in increased virus infection and replication. Our research demonstrates that metabolites produced by the gut microbiome, such as butyrate, can have complex effects on cellular physiology, including dampening of an inflammatory innate immune pathway resulting in a proviral cellular environment. Our work further suggests that butyrate could be broadly used as a tool to increase growth of virus stocks for research and for the generation of vaccines.

Functional diversification accompanies gene family expansion of MED2 homologs in Candida albicans.

Gene duplication facilitates functional diversification and provides greater phenotypic flexibility to an organism. Expanded gene families arise through repeated gene duplication but the extent of functional divergence that accompanies each paralogous gene is generally unexplored because of the difficulty in isolating the effects of single family members. The telomere-associated (TLO) gene family is a remarkable example of gene family expansion, with 14 members in the more pathogenic Candida albicans relative to two TLO genes in the closely-related species C. dubliniensis. TLO genes encode interchangeable Med2 subunits of the major transcriptional regulatory complex Mediator. To identify biological functions associated with each C. albicans TLO, expression of individual family members was regulated using a Tet-ON system and the strains were assessed across a range of phenotypes involved in growth and virulence traits. All TLOs affected multiple phenotypes and a single phenotype was often affected by multiple TLOs, including simple phenotypes such as cell aggregation and complex phenotypes such as virulence in a Galleria mellonella model of infection. No phenotype was regulated by all TLOs, suggesting neofunctionalization or subfunctionalization of ancestral properties among different family members. Importantly, regulation of three phenotypes could be mapped to individual polymorphic sites among the TLO genes, including an indel correlated with two phenotypes, growth in sucrose and macrophage killing. Different selective pressures have operated on the TLO sequence, with the 5' conserved Med2 domain experiencing purifying selection and the gene/clade-specific 3' end undergoing extensive positive selection that may contribute to the impact of individual TLOs on phenotypic variability. Therefore, expansion of the TLO gene family has conferred unique regulatory properties to each paralog such that it influences a range of phenotypes. We posit that the genetic diversity associated with this expansion contributed to C. albicans success as a commensal and opportunistic pathogen.

Epigenetic control of pheromone MAPK signaling determines sexual fecundity in Candida albicans.

Several pathogenic Candida species are capable of heritable and reversible switching between two epigenetic states, "white" and "opaque." In Candida albicans, white cells are essentially sterile, whereas opaque cells are mating-proficient. Here, we interrogate the mechanism by which the white-opaque switch regulates sexual fecundity and identify four genes in the pheromone MAPK pathway that are expressed at significantly higher levels in opaque cells than in white cells. These genes encode the ß subunit of the G-protein complex (STE4), the pheromone MAPK scaffold (CST5), and the two terminal MAP kinases (CEK1/CEK2). To define the contribution of each factor to mating, C. albicans white cells were reverse-engineered to express elevated, opaque-like levels of these factors, either singly or in combination. We show that white cells co-overexpressing STE4, CST5, and CEK2 undergo mating four orders of magnitude more efficiently than control white cells and at a frequency approaching that of opaque cells. Moreover, engineered white cells recapitulate the transcriptional and morphological responses of opaque cells to pheromone. These results therefore reveal multiple bottlenecks in pheromone MAPK signaling in white cells and that alleviation of these bottlenecks enables efficient mating by these "sterile" cell types. Taken together, our findings establish that differential expression of several MAPK factors underlies the epigenetic control of mating in C. albicans We also discuss how fitness advantages could have driven the evolution of a toggle switch to regulate sexual reproduction in pathogenic Candida species.

Role of Mediator in virulence and antifungal drug resistance in pathogenic fungi.

Mediator complex has recently emerged as an important regulator of gene expression in pathogenic fungi. Mediator is a multi-subunit complex of polypeptides involved in transcriptional activation in eukaryotes, with roles including preinitiation complex (PIC) assembly and chromatin remodeling. Within the last decade, Mediator has been shown to play an integral role in regulating virulence gene expression and drug resistance in human fungal pathogens. In some fungi, specific Mediator subunits have been shown to be required for virulence. In Candida species, duplication and expansion of Mediator subunit encoding genes has occurred on at least three occasions (CgMED15 in C. glabrata and MED2/TLO in C. albicans and C. dubliniensis) suggesting important roles for Mediator in the evolution of these pathogens. This review summarises recent developments in our understanding of Mediator in fungal pathogens and the potential for the development of therapeutic drugs to target Mediator functions.

A Multistate Toggle Switch Defines Fungal Cell Fates and Is Regulated by Synergistic Genetic Cues.

Heritable epigenetic changes underlie the ability of cells to differentiate into distinct cell types. Here, we demonstrate that the fungal pathogen Candida tropicalis exhibits multipotency, undergoing stochastic and reversible switching between three cellular states. The three cell states exhibit unique cellular morphologies, growth rates, and global gene expression profiles. Genetic analysis identified six transcription factors that play key roles in regulating cell differentiation. In particular, we show that forced expression of Wor1 or Efg1 transcription factors can be used to manipulate transitions between all three cell states. A model for tristability is proposed in which Wor1 and Efg1 are self-activating but mutually antagonistic transcription factors, thereby forming a symmetrical self-activating toggle switch. We explicitly test this model and show that ectopic expression of WOR1 can induce white-to-hybrid-to-opaque switching, whereas ectopic expression of EFG1 drives switching in the opposite direction, from opaque-to-hybrid-to-white cell states. We also address the stability of induced cell states and demonstrate that stable differentiation events require ectopic gene expression in combination with chromatin-based cues. These studies therefore experimentally test a model of multistate stability and demonstrate that transcriptional circuits act synergistically with chromatin-based changes to drive cell state transitions. We also establish close mechanistic parallels between phenotypic switching in unicellular fungi and cell fate decisions during stem cell reprogramming.

Genetic and phenotypic intra-species variation in Candida albicans.

Candida albicans is a commensal fungus of the human gastrointestinal tract and a prevalent opportunistic pathogen. To examine diversity within this species, extensive genomic and phenotypic analyses were performed on 21 clinical C. albicans isolates. Genomic variation was evident in the form of polymorphisms, copy number variations, chromosomal inversions, subtelomeric hypervariation, loss of heterozygosity (LOH), and whole or partial chromosome aneuploidies. All 21 strains were diploid, although karyotypic changes were present in eight of the 21 isolates, with multiple strains being trisomic for Chromosome 4 or Chromosome 7. Aneuploid strains exhibited a general fitness defect relative to euploid strains when grown under replete conditions. All strains were also heterozygous, yet multiple, distinct LOH tracts were present in each isolate. Higher overall levels of genome heterozygosity correlated with faster growth rates, consistent with increased overall fitness. Genes with the highest rates of amino acid substitutions included many cell wall proteins, implicating fast evolving changes in cell adhesion and host interactions. One clinical isolate, P94015, presented several striking properties including a novel cellular phenotype, an inability to filament, drug resistance, and decreased virulence. Several of these properties were shown to be due to a homozygous nonsense mutation in the EFG1 gene. Furthermore, loss of EFG1 function resulted in increased fitness of P94015 in a commensal model of infection. Our analysis therefore reveals intra-species genetic and phenotypic differences in C. albicans and delineates a natural mutation that alters the balance between commensalism and pathogenicity.

A chromosome 4 trisomy contributes to increased fluconazole resistance in a clinical isolate of Candida albicans.

Candida albicans is an important opportunistic fungal pathogen capable of causing both mucosal and disseminated disease. Infections are often treated with fluconazole, a front-line antifungal drug that targets the biosynthesis of ergosterol, a major component of the fungal cell membrane. Resistance to fluconazole can arise through a variety of mechanisms, including gain-of-function mutations, loss of heterozygosity events and aneuploidy. The clinical isolate P60002 was found to be highly resistant to azole-class drugs, yet lacked mutations or chromosomal rearrangements known to be associated with azole resistance. Transcription profiling suggested that increased expression of two putative drug efflux pumps, CDR11 and QDR1, might confer azole resistance. However, ectopic expression of the P60002 alleles of these genes in a drug-susceptible strain did not increase fluconazole resistance. We next examined whether the presence of three copies of chromosome 4 (Chr4) or chromosome 6 (Chr6) contributed to azole resistance in P60002. We established that Chr4 trisomy contributes significantly to fluconazole resistance, whereas Chr6 trisomy has no discernible effect on resistance. In contrast, a Chr4 trisomy did not increase fluconazole resistance when present in the standard SC5314 strain background. These results establish a link between Chr4 trisomy and elevated fluconazole resistance, and demonstrate the impact of genetic background on drug resistance phenotypes in C. albicans.

Silencing is noisy: population and cell level noise in telomere-adjacent genes is dependent on telomere position and sir2.

Cell-to-cell gene expression noise is thought to be an important mechanism for generating phenotypic diversity. Furthermore, telomeric regions are major sites for gene amplification, which is thought to drive genetic diversity. Here we found that individual subtelomeric TLO genes exhibit increased variation in transcript and protein levels at both the cell-to-cell level as well as at the population-level. The cell-to-cell variation, termed Telomere-Adjacent Gene Expression Noise (TAGEN) was largely intrinsic noise and was dependent upon genome position: noise was reduced when a TLO gene was expressed at an ectopic internal locus and noise was elevated when a non-telomeric gene was expressed at a telomere-adjacent locus. This position-dependent TAGEN also was dependent on Sir2p, an NAD+-dependent histone deacetylase. Finally, we found that telomere silencing and TAGEN are tightly linked and regulated in cis: selection for either silencing or activation of a TLO-adjacent URA3 gene resulted in reduced noise at the neighboring TLO but not at other TLO genes. This provides experimental support to computational predictions that the ability to shift between silent and active chromatin states has a major effect on cell-to-cell noise. Furthermore, it demonstrates that these shifts affect the degree of expression variation at each telomere individually.

Telomeric ORFs (TLOs) in Candida spp. Encode mediator subunits that regulate distinct virulence traits.

The TLO genes are a family of telomere-associated ORFs in the fungal pathogens Candida albicans and C. dubliniensis that encode a subunit of the Mediator complex with homology to Med2. The more virulent pathogen C. albicans has 15 copies of the gene whereas the less pathogenic species C. dubliniensis has only two (CdTLO1 and CdTLO2). In this study we used C. dubliniensis as a model to investigate the role of TLO genes in regulating virulence and also to determine whether TLO paralogs have evolved to regulate distinct functions. A C. dubliniensis tlo1Δ/tlo2Δ mutant is unable to form true hyphae, has longer doubling times in galactose broth, is more susceptible to oxidative stress and forms increased levels of biofilm. Transcript profiling of the tlo1Δ/tlo2Δ mutant revealed increased expression of starvation responses in rich medium and retarded expression of hypha-induced transcripts in serum. ChIP studies indicated that Tlo1 binds to many ORFs including genes that exhibit high and low expression levels under the conditions analyzed. The altered expression of these genes in the tlo1Δ/tlo2Δ null mutant indicates roles for Tlo proteins in transcriptional activation and repression. Complementation of the tlo1Δ/tlo2Δ mutant with TLO1, but not TLO2, restored wild-type filamentous growth, whereas only TLO2 fully suppressed biofilm growth. Complementation with TLO1 also had a greater effect on doubling times in galactose broth. The different abilities of TLO1 and TLO2 to restore wild-type functions was supported by transcript profiling studies that showed that only TLO1 restored expression of hypha-specific genes (UME6, SOD5) and galactose utilisation genes (GAL1 and GAL10), whereas TLO2 restored repression of starvation-induced gene transcription. Thus, Tlo/Med2 paralogs encoding Mediator subunits regulate different virulence properties in Candida spp. and their expansion may account for the increased adaptability of C. albicans relative to other Candida species.

Candida albicans' inorganic phosphate transport and evolutionary adaptation to phosphate scarcity.

Phosphorus is essential in all cells' structural, metabolic and regulatory functions. For fungal cells that import inorganic phosphate (Pi) up a steep concentration gradient, surface Pi transporters are critical capacitators of growth. Fungi must deploy Pi transporters that enable optimal Pi uptake in pH and Pi concentration ranges prevalent in their environments. Single, triple and quadruple mutants were used to characterize the four Pi transporters we identified for the human fungal pathogen Candida albicans, which must adapt to alkaline conditions during invasion of the host bloodstream and deep organs. A high-affinity Pi transporter, Pho84, was most efficient across the widest pH range while another, Pho89, showed high-affinity characteristics only within one pH unit of neutral. Two low-affinity Pi transporters, Pho87 and Fgr2, were active only in acidic conditions. Only Pho84 among the Pi transporters was clearly required in previously identified Pi-related functions including Target of Rapamycin Complex 1 signaling and hyphal growth. We used in vitro evolution and whole genome sequencing as an unbiased forward genetic approach to probe adaptation to prolonged Pi scarcity of two quadruple mutant lineages lacking all 4 Pi transporters. Lineage-specific genomic changes corresponded to divergent success of the two lineages in fitness recovery during Pi limitation. In this process, initial, large-scale genomic alterations like aneuploidies and loss of heterozygosity were eventually lost as populations presumably gained small-scale mutations. Severity of some phenotypes linked to Pi starvation, like cell wall stress hypersensitivity, decreased in parallel to evolving populations' fitness recovery in Pi scarcity, while that of others like membrane stress responses diverged from these fitness phenotypes. C. albicans therefore has diverse options to reconfigure Pi management during prolonged scarcity. Since Pi homeostasis differs substantially between fungi and humans, adaptive processes to Pi deprivation may harbor small-molecule targets that impact fungal growth and virulence.

The three clades of the telomere-associated TLO gene family of Candida albicans have different splicing, localization, and expression features.

Candida albicans grows within a wide range of host niches, and this adaptability enhances its success as a commensal and as a pathogen. The telomere-associated TLO gene family underwent a recent expansion from one or two copies in other CUG clade members to 14 expressed copies in C. albicans. This correlates with increased virulence and clinical prevalence relative to those of other Candida clade species. The 14 expressed TLO gene family members have a conserved Med2 domain at the N terminus, suggesting a role in general transcription. The C-terminal half is more divergent, distinguishing three clades: clade α and clade ß have no introns and encode proteins that localize primarily to the nucleus; clade γ sometimes undergoes splicing, and the gene products localize within the mitochondria as well as the nuclei. Additionally, TLOα genes are generally expressed at much higher levels than are TLOγ genes. We propose that expansion of the TLO gene family and the predicted role of Tlo proteins in transcription regulation provide C. albicans with the ability to adapt rapidly to the broad range of different environmental niches within the human host.

Microbiome ownership for Indigenous peoples.

Several studies have reported increased microbial diversity, or distinct microbial community compositions, in the microbiomes of Indigenous peoples around the world. However, there is a widespread failure to include Indigenous cultures and perspectives in microbiome research programmes, and ethical issues pertaining to microbiome research involving Indigenous participants have not received enough attention. We discuss the benefits and risks arising from microbiome research involving Indigenous peoples and analyse microbiome ownership as an ethical concept in this context. We argue that microbiome ownership represents an opportunity for Indigenous peoples to steward and protect their resident microbial communities at every stage of research.

A relational framework for microbiome research with Indigenous communities.

Ethical practices in human microbiome research have failed to keep pace with scientific advances in the field. Researchers seeking to 'preserve' microbial species associated with Indigenous groups, but absent from industrialized populations, have largely failed to include Indigenous people in knowledge co-production or benefit, perpetuating a legacy of intellectual and material extraction. We propose a framework centred on relationality among Indigenous peoples, researchers and microbes, to guide ethical microbiome research. Our framework centres accountability to flatten historical power imbalances that favour researcher perspectives and interests to provide space for Indigenous worldviews in pursuit of Indigenous research sovereignty. Ethical inclusion of Indigenous communities in microbiome research can provide health benefits for all populations and reinforce mutually beneficial partnerships between researchers and the public.

The Many Hosts of Mycobacteria 9 (MHM9): A conference report.

The Many Hosts of Mycobacteria (MHM) meeting series brings together basic scientists, clinicians and veterinarians to promote robust discussion and dissemination of recent advances in our knowledge of numerous mycobacterial diseases, including human and bovine tuberculosis (TB), nontuberculous mycobacteria (NTM) infection, Hansen's disease (leprosy), Buruli ulcer and Johne's disease. The 9th MHM conference (MHM9) was held in July 2022 at The Ohio State University (OSU) and centered around the theme of "Confounders of Mycobacterial Disease." Confounders can and often do drive the transmission of mycobacterial diseases, as well as impact surveillance and treatment outcomes. Various confounders were presented and discussed at MHM9 including those that originate from the host (comorbidities and coinfections) as well as those arising from the environment (e.g., zoonotic exposures), economic inequality (e.g. healthcare disparities), stigma (a confounder of leprosy and TB for millennia), and historical neglect (a confounder in Native American Nations). This conference report summarizes select talks given at MHM9 highlighting recent research advances, as well as talks regarding the historic and ongoing impact of TB and other infectious diseases on Native American Nations, including those in Southwestern Alaska where the regional TB incidence rate is among the highest in the Western hemisphere.

Architectural groups of a subtelomeric gene family evolve along distinct paths in Candida albicans.

Subtelomeres are dynamic genomic regions shaped by elevated rates of recombination, mutation, and gene birth/death. These processes contribute to formation of lineage-specific gene family expansions that commonly occupy subtelomeres across eukaryotes. Investigating the evolution of subtelomeric gene families is complicated by the presence of repetitive DNA and high sequence similarity among gene family members that prevents accurate assembly from whole genome sequences. Here, we investigated the evolution of the telomere-associated (TLO) gene family in Candida albicans using 189 complete coding sequences retrieved from 23 genetically diverse strains across the species. Tlo genes conformed to the 3 major architectural groups (α/ß/γ) previously defined in the genome reference strain but significantly differed in the degree of within-group diversity. One group, Tloß, was always found at the same chromosome arm with strong sequence similarity among all strains. In contrast, diverse Tloα sequences have proliferated among chromosome arms. Tloγ genes formed 7 primary clades that included each of the previously identified Tloγ genes from the genome reference strain with 3 Tloγ genes always found on the same chromosome arm among strains. Architectural groups displayed regions of high conservation that resolved newly identified functional motifs, providing insight into potential regulatory mechanisms that distinguish groups. Thus, by resolving intraspecies subtelomeric gene variation, it is possible to identify previously unknown gene family complexity that may underpin adaptive functional variation.

Parasexuality of Candida Species.

While most fungi have the ability to reproduce sexually, multiple independent lineages have lost meiosis and developed parasexual cycles in its place. Emergence of parasexual cycles is particularly prominent in medically relevant fungi from the CUG paraphyletic group of Candida species. Since the discovery of parasex in C. albicans roughly two decades ago, it has served as the model for Candida species. Importantly, parasex in C. albicans retains hallmarks of meiosis including genetic recombination and chromosome segregation, making it a potential driver of genetic diversity. Furthermore, key meiotic genes play similar roles in C. albicans parasex and highlights parallels between these processes. Yet, the evolutionary role of parasex in Candida adaptation and the extent of resulting genotypic and phenotypic diversity remain as key knowledge gaps in this facultative reproductive program. Here, we present our current understanding of parasex, the mechanisms governing its regulation, and its relevance to Candida biology.

Intraspecies Transcriptional Profiling Reveals Key Regulators of Candida albicans Pathogenic Traits.

The human commensal and opportunistic fungal pathogen Candida albicans displays extensive genetic and phenotypic variation across clinical isolates. Here, we performed RNA sequencing on 21 well-characterized isolates to examine how genetic variation contributes to gene expression differences and to link these differences to phenotypic traits. C. albicans adapts primarily through clonal evolution, and yet hierarchical clustering of gene expression profiles in this set of isolates did not reproduce their phylogenetic relationship. Strikingly, strain-specific gene expression was prevalent in some strain backgrounds. Association of gene expression with phenotypic data by differential analysis, linear correlation, and assembly of gene networks connected both previously characterized and novel genes with 23 C. albicans traits. Construction of de novo gene modules produced a gene atlas incorporating 67% of C. albicans genes and revealed correlations between expression modules and important phenotypes such as systemic virulence. Furthermore, targeted investigation of two modules that have novel roles in growth and filamentation supported our bioinformatic predictions. Together, these studies reveal widespread transcriptional variation across C. albicans isolates and identify genetic and epigenetic links to phenotypic variation based on coexpression network analysis.IMPORTANCE Infectious fungal species are often treated uniformly despite clear evidence of genotypic and phenotypic heterogeneity being widespread across strains. Identifying the genetic basis for this phenotypic diversity is extremely challenging because of the tens or hundreds of thousands of variants that may distinguish two strains. Here, we use transcriptional profiling to determine differences in gene expression that can be linked to phenotypic variation among a set of 21 Candida albicans isolates. Analysis of this transcriptional data set uncovered clear trends in gene expression characteristics for this species and new genes and pathways that were associated with variation in pathogenic processes. Direct investigation confirmed functional predictions for a number of new regulators associated with growth and filamentation, demonstrating the utility of these approaches in linking genes to important phenotypes.