Identification of driver genes for critical forms of COVID-19 in a deeply phenotyped young patient cohort.

The drivers of critical coronavirus disease 2019 (COVID-19) remain unknown. Given major confounding factors such as age and comorbidities, true mediators of this condition have remained elusive. We used a multi-omics analysis combined with artificial intelligence in a young patient cohort where major comorbidities were excluded at the onset. The cohort included 47 "critical" (in the intensive care unit under mechanical ventilation) and 25 "non-critical" (in a non-critical care ward) patients with COVID-19 and 22 healthy individuals. The analyses included whole-genome sequencing, whole-blood RNA sequencing, plasma and blood mononuclear cell proteomics, cytokine profiling, and high-throughput immunophenotyping. An ensemble of machine learning, deep learning, quantum annealing, and structural causal modeling were used. Patients with critical COVID-19 were characterized by exacerbated inflammation, perturbed lymphoid and myeloid compartments, increased coagulation, and viral cell biology. Among differentially expressed genes, we observed up-regulation of the metalloprotease ADAM9. This gene signature was validated in a second independent cohort of 81 critical and 73 recovered patients with COVID-19 and was further confirmed at the transcriptional and protein level and by proteolytic activity. Ex vivo ADAM9 inhibition decreased severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) uptake and replication in human lung epithelial cells. In conclusion, within a young, otherwise healthy, cohort of individuals with COVID-19, we provide the landscape of biological perturbations in vivo where a unique gene signature differentiated critical from non-critical patients. We further identified ADAM9 as a driver of disease severity and a candidate therapeutic target.

Quantum processor-inspired machine learning in the biomedical sciences.

Recent advances in high-throughput genomic technologies coupled with exponential increases in computer processing and memory have allowed us to interrogate the complex molecular underpinnings of human disease from a genome-wide perspective. While the deluge of genomic information is expected to increase, a bottleneck in conventional high-performance computing is rapidly approaching. Inspired by recent advances in physical quantum processors, we evaluated several unconventional machine-learning (ML) strategies on actual human tumor data, namely "Ising-type" methods, whose objective function is formulated identical to simulated annealing and quantum annealing. We show the efficacy of multiple Ising-type ML algorithms for classification of multi-omics human cancer data from The Cancer Genome Atlas, comparing these classifiers to a variety of standard ML methods. Our results indicate that Ising-type ML offers superior classification performance with smaller training datasets, thus providing compelling empirical evidence for the potential future application of unconventional computing approaches in the biomedical sciences.

Smooth Muscle Cell Reprogramming in Aortic Aneurysms.

The etiology of aortic aneurysms is poorly understood, but it is associated with atherosclerosis, hypercholesterolemia, and abnormal transforming growth factor ß (TGF-ß) signaling in smooth muscle. Here, we investigated the interactions between these different factors in aortic aneurysm development and identified a key role for smooth muscle cell (SMC) reprogramming into a mesenchymal stem cell (MSC)-like state. SMC-specific ablation of TGF-ß signaling in Apoe-/- mice on a hypercholesterolemic diet led to development of aortic aneurysms exhibiting all the features of human disease, which was associated with transdifferentiation of a subset of contractile SMCs into an MSC-like intermediate state that generated osteoblasts, chondrocytes, adipocytes, and macrophages. This combination of medial SMC loss with marked increases in non-SMC aortic cell mass induced exuberant growth and dilation of the aorta, calcification and ossification of the aortic wall, and inflammation, resulting in aneurysm development.

Chronic mTOR activation induces a degradative smooth muscle cell phenotype.

Smooth muscle cell (SMC) proliferation has been thought to limit the progression of thoracic aortic aneurysm and dissection (TAAD) because loss of medial cells associates with advanced disease. We investigated effects of SMC proliferation in the aortic media by conditional disruption of Tsc1, which hyperactivates mTOR complex 1. Consequent SMC hyperplasia led to progressive medial degeneration and TAAD. In addition to diminished contractile and synthetic functions, fate-mapped SMCs displayed increased proteolysis, endocytosis, phagocytosis, and lysosomal clearance of extracellular matrix and apoptotic cells. SMCs acquired a limited repertoire of macrophage markers and functions via biogenesis of degradative organelles through an mTOR/ß-catenin/MITF-dependent pathway, but were distinguishable from conventional macrophages by an absence of hematopoietic lineage markers and certain immune effectors even in the context of hyperlipidemia. Similar mTOR activation and induction of a degradative SMC phenotype in a model of mild TAAD due to Fbn1 mutation greatly worsened disease with near-uniform lethality. The finding of increased lysosomal markers in medial SMCs from clinical TAAD specimens with hyperplasia and matrix degradation further supports the concept that proliferation of degradative SMCs within the media causes aortic disease, thus identifying mTOR-dependent phenotypic modulation as a therapeutic target for combating TAAD.

Endothelial TGF-ß signalling drives vascular inflammation and atherosclerosis.

Atherosclerosis is a progressive vascular disease triggered by interplay between abnormal shear stress and endothelial lipid retention. A combination of these and, potentially, other factors leads to a chronic inflammatory response in the vessel wall, which is thought to be responsible for disease progression characterized by a buildup of atherosclerotic plaques. Yet molecular events responsible for maintenance of plaque inflammation and plaque growth have not been fully defined. Here we show that endothelial TGFß signaling is one of the primary drivers of atherosclerosis-associated vascular inflammation. Inhibition of endothelial TGFß signaling in hyperlipidemic mice reduces vessel wall inflammation and vascular permeability and leads to arrest of disease progression and regression of established lesions. These pro-inflammatory effects of endothelial TGFß signaling are in stark contrast with its effects in other cell types and identify it as an important driver of atherosclerotic plaque growth and show the potential of cell-type specific therapeutic intervention aimed at control of this disease.

Genome Sequence for Candida albicans Clinical Oral Isolate 529L.

The diploid heterozygous yeast Candida albicans is the most common cause of fungal infection. Here, we report the genome sequence assembly of the clinical oral isolate 529L. As this isolate grows as a commensal, this genome will serve as a reference for experimental and genetic studies of mucosal colonization.

Ligand-activated BMP signaling inhibits cell differentiation and death to promote melanoma.

Oncogenomic studies indicate that copy number variation (CNV) alters genes involved in tumor progression; however, identification of specific driver genes affected by CNV has been difficult, as these rearrangements are often contained in large chromosomal intervals among several bystander genes. Here, we addressed this problem and identified a CNV-targeted oncogene by performing comparative oncogenomics of human and zebrafish melanomas. We determined that the gene encoding growth differentiation factor 6 (GDF6), which is the ligand for the BMP family, is recurrently amplified and transcriptionally upregulated in melanoma. GDF6-induced BMP signaling maintained a trunk neural crest gene signature in melanomas. Additionally, GDF6 repressed the melanocyte differentiation gene MITF and the proapoptotic factor SOX9, thereby preventing differentiation, inhibiting cell death, and promoting tumor growth. GDF6 was specifically expressed in melanomas but not melanocytes. Moreover, GDF6 expression levels in melanomas were inversely correlated with patient survival. Our study has identified a fundamental role for GDF6 and BMP signaling in governing an embryonic cell gene signature to promote melanoma progression, thus providing potential opportunities for targeted therapy to treat GDF6-positive cancers.

Brown Fat AKT2 Is a Cold-Induced Kinase that Stimulates ChREBP-Mediated De Novo Lipogenesis to Optimize Fuel Storage and Thermogenesis.

Brown adipose tissue (BAT) is a therapeutic target for metabolic diseases; thus, understanding its metabolic circuitry is clinically important. Many studies of BAT compare rodents mildly cold to those severely cold. Here, we compared BAT remodeling between thermoneutral and mild-cold-adapted mice, conditions more relevant to humans. Although BAT is renowned for catabolic ß-oxidative capacity, we find paradoxically that the anabolic de novo lipogenesis (DNL) genes encoding ACLY, ACSS2, ACC, and FASN were among the most upregulated by mild cold and that, in humans, DNL correlates with Ucp1 expression. The regulation and function of adipocyte DNL and its association with thermogenesis are not understood. We provide evidence suggesting that AKT2 drives DNL in adipocytes by stimulating ChREBPß transcriptional activity and that cold induces the AKT2-ChREBP pathway in BAT to optimize fuel storage and thermogenesis. These data provide insight into adipocyte DNL regulation and function and illustrate the metabolic flexibility of thermogenesis.

Tracing Genetic Exchange and Biogeography of Cryptococcus neoformans var. grubii at the Global Population Level.

Cryptococcus neoformans var. grubii is the causative agent of cryptococcal meningitis, a significant source of mortality in immunocompromised individuals, typically human immunodeficiency virus/AIDS patients from developing countries. Despite the worldwide emergence of this ubiquitous infection, little is known about the global molecular epidemiology of this fungal pathogen. Here we sequence the genomes of 188 diverse isolates and characterize the major subdivisions, their relative diversity, and the level of genetic exchange between them. While most isolates of C. neoformans var. grubii belong to one of three major lineages (VNI, VNII, and VNB), some haploid isolates show hybrid ancestry including some that appear to have recently interbred, based on the detection of large blocks of each ancestry across each chromosome. Many isolates display evidence of aneuploidy, which was detected for all chromosomes. In diploid isolates of C. neoformans var. grubii (serotype AA) and of hybrids with C. neoformans var. neoformans (serotype AD) such aneuploidies have resulted in loss of heterozygosity, where a chromosomal region is represented by the genotype of only one parental isolate. Phylogenetic and population genomic analyses of isolates from Brazil reveal that the previously "African" VNB lineage occurs naturally in the South American environment. This suggests migration of the VNB lineage between Africa and South America prior to its diversification, supported by finding ancestral recombination events between isolates from different lineages and regions. The results provide evidence of substantial population structure, with all lineages showing multi-continental distributions; demonstrating the highly dispersive nature of this pathogen.

Structure of the germline genome of Tetrahymena thermophila and relationship to the massively rearranged somatic genome.

The germline genome of the binucleated ciliate Tetrahymena thermophila undergoes programmed chromosome breakage and massive DNA elimination to generate the somatic genome. Here, we present a complete sequence assembly of the germline genome and analyze multiple features of its structure and its relationship to the somatic genome, shedding light on the mechanisms of genome rearrangement as well as the evolutionary history of this remarkable germline/soma differentiation. Our results strengthen the notion that a complex, dynamic, and ongoing interplay between mobile DNA elements and the host genome have shaped Tetrahymena chromosome structure, locally and globally. Non-standard outcomes of rearrangement events, including the generation of short-lived somatic chromosomes and excision of DNA interrupting protein-coding regions, may represent novel forms of developmental gene regulation. We also compare Tetrahymena's germline/soma differentiation to that of other characterized ciliates, illustrating the wide diversity of adaptations that have occurred within this phylum.

Gld2-catalyzed 3' monoadenylation of miRNAs in the hippocampus has no detectable effect on their stability or on animal behavior.

Gld2, a noncanonical cytoplasmic poly(A) polymerase, interacts with the RNA binding protein CPEB1 to mediate polyadenylation-induced translation in dendrites of cultured hippocampal neurons. Depletion of Gld2 from the hippocampus leads to a deficit in long-term potentiation evoked by theta burst stimulation. At least in mouse liver and human primary fibroblasts, Gld2 also 3' monoadenylates and thereby stabilizes specific miRNAs, which enhance mRNA translational silencing and eventual destruction. These results suggest that Gld2 would be likely to monoadenylate and stabilize miRNAs in the hippocampus, which would produce measurable changes in animal behavior. We now report that using Gld2 knockout mice, there are detectable alterations in specific miRNA monoadenylation in the hippocampus when compared to wild type, but that these modifications produce no detectable effect on miRNA stability. Moreover, we surprisingly find no overt change in animal behavior when comparing Gld2 knockout to wild-type mice. These data indicate that miRNA monoadenylation-mediated stability is cell type-specific and that monoadenylation has no measurable effect on higher cognitive function.

Genome Sequence of Spizellomyces punctatus.

Spizellomyces punctatus is a basally branching chytrid fungus that is found in the Chytridiomycota phylum. Spizellomyces species are common in soil and of importance in terrestrial ecosystems. Here, we report the genome sequence of S. punctatus, which will facilitate the study of this group of early diverging fungi.

Cholesterol-Independent SREBP-1 Maturation Is Linked to ARF1 Inactivation.

Lipogenesis requires coordinated expression of genes for fatty acid, phospholipid, and triglyceride synthesis. Transcription factors, such as SREBP-1 (Sterol regulatory element binding protein), may be activated in response to feedback mechanisms linking gene activation to levels of metabolites in the pathways. SREBPs can be regulated in response to membrane cholesterol and we also found that low levels of phosphatidylcholine (a methylated phospholipid) led to SBP-1/SREBP-1 maturation in C. elegans or mammalian models. To identify additional regulatory components, we performed a targeted RNAi screen in C. elegans, finding that both lpin-1/Lipin 1 (which converts phosphatidic acid to diacylglycerol) and arf-1.2/ARF1 (a GTPase regulating Golgi function) were important for low-PC activation of SBP-1/SREBP-1. Mechanistically linking the major hits of our screen, we find that limiting PC synthesis or LPIN1 knockdown in mammalian cells reduces the levels of active GTP-bound ARF1. Thus, changes in distinct lipid ratios may converge on ARF1 to increase SBP-1/SREBP-1 activity.

Endothelial Mitogen-Activated Protein Kinase Kinase Kinase Kinase 4 Is Critical for Lymphatic Vascular Development and Function.

The molecular mechanisms underlying lymphatic vascular development and function are not well understood. Recent studies have suggested a role for endothelial cell (EC) mitogen-activated protein kinase kinase kinase kinase 4 (Map4k4) in developmental angiogenesis and atherosclerosis. Here, we show that constitutive loss of EC Map4k4 in mice causes postnatal lethality due to chylothorax, suggesting that Map4k4 is required for normal lymphatic vascular function. Mice constitutively lacking EC Map4k4 displayed dilated lymphatic capillaries, insufficient lymphatic valves, and impaired lymphatic flow; furthermore, primary ECs derived from these animals displayed enhanced proliferation compared with controls. Yeast 2-hybrid analyses identified the Ras GTPase-activating protein Rasa1, a known regulator of lymphatic development and lymphatic endothelial cell fate, as a direct interacting partner for Map4k4. Map4k4 silencing in ECs enhanced basal Ras and extracellular signal-regulated kinase (Erk) activities, and primary ECs lacking Map4k4 displayed enhanced lymphatic EC marker expression. Taken together, these results reveal that EC Map4k4 is critical for lymphatic vascular development by regulating EC quiescence and lymphatic EC fate.

The Dynamic Genome and Transcriptome of the Human Fungal Pathogen Blastomyces and Close Relative Emmonsia.

Three closely related thermally dimorphic pathogens are causal agents of major fungal diseases affecting humans in the Americas: blastomycosis, histoplasmosis and paracoccidioidomycosis. Here we report the genome sequence and analysis of four strains of the etiological agent of blastomycosis, Blastomyces, and two species of the related genus Emmonsia, typically pathogens of small mammals. Compared to related species, Blastomyces genomes are highly expanded, with long, often sharply demarcated tracts of low GC-content sequence. These GC-poor isochore-like regions are enriched for gypsy elements, are variable in total size between isolates, and are least expanded in the avirulent B. dermatitidis strain ER-3 as compared with the virulent B. gilchristii strain SLH14081. The lack of similar regions in related species suggests these isochore-like regions originated recently in the ancestor of the Blastomyces lineage. While gene content is highly conserved between Blastomyces and related fungi, we identified changes in copy number of genes potentially involved in host interaction, including proteases and characterized antigens. In addition, we studied gene expression changes of B. dermatitidis during the interaction of the infectious yeast form with macrophages and in a mouse model. Both experiments highlight a strong antioxidant defense response in Blastomyces, and upregulation of dioxygenases in vivo suggests that dioxide produced by antioxidants may be further utilized for amino acid metabolism. We identify a number of functional categories upregulated exclusively in vivo, such as secreted proteins, zinc acquisition proteins, and cysteine and tryptophan metabolism, which may include critical virulence factors missed before in in vitro studies. Across the dimorphic fungi, loss of certain zinc acquisition genes and differences in amino acid metabolism suggest unique adaptations of Blastomyces to its host environment. These results reveal the dynamics of genome evolution and of factors contributing to virulence in Blastomyces.

Evolution of Extensively Drug-Resistant Tuberculosis over Four Decades: Whole Genome Sequencing and Dating Analysis of Mycobacterium tuberculosis Isolates from KwaZulu-Natal.

BACKGROUND: The continued advance of antibiotic resistance threatens the treatment and control of many infectious diseases. This is exemplified by the largest global outbreak of extensively drug-resistant (XDR) tuberculosis (TB) identified in Tugela Ferry, KwaZulu-Natal, South Africa, in 2005 that continues today. It is unclear whether the emergence of XDR-TB in KwaZulu-Natal was due to recent inadequacies in TB control in conjunction with HIV or other factors. Understanding the origins of drug resistance in this fatal outbreak of XDR will inform the control and prevention of drug-resistant TB in other settings. In this study, we used whole genome sequencing and dating analysis to determine if XDR-TB had emerged recently or had ancient antecedents. METHODS AND FINDINGS: We performed whole genome sequencing and drug susceptibility testing on 337 clinical isolates of Mycobacterium tuberculosis collected in KwaZulu-Natal from 2008 to 2013, in addition to three historical isolates, collected from patients in the same province and including an isolate from the 2005 Tugela Ferry XDR outbreak, a multidrug-resistant (MDR) isolate from 1994, and a pansusceptible isolate from 1995. We utilized an array of whole genome comparative techniques to assess the relatedness among strains, to establish the order of acquisition of drug resistance mutations, including the timing of acquisitions leading to XDR-TB in the LAM4 spoligotype, and to calculate the number of independent evolutionary emergences of MDR and XDR. Our sequencing and analysis revealed a 50-member clone of XDR M. tuberculosis that was highly related to the Tugela Ferry XDR outbreak strain. We estimated that mutations conferring isoniazid and streptomycin resistance in this clone were acquired 50 y prior to the Tugela Ferry outbreak (katG S315T [isoniazid]; gidB 130 bp deletion [streptomycin]; 1957 [95% highest posterior density (HPD): 1937-1971]), with the subsequent emergence of MDR and XDR occurring 20 y (rpoB L452P [rifampicin]; pncA 1 bp insertion [pyrazinamide]; 1984 [95% HPD: 1974-1992]) and 10 y (rpoB D435G [rifampicin]; rrs 1400 [kanamycin]; gyrA A90V [ofloxacin]; 1995 [95% HPD: 1988-1999]) prior to the outbreak, respectively. We observed frequent de novo evolution of MDR and XDR, with 56 and nine independent evolutionary events, respectively. Isoniazid resistance evolved before rifampicin resistance 46 times, whereas rifampicin resistance evolved prior to isoniazid only twice. We identified additional putative compensatory mutations to rifampicin in this dataset. One major limitation of this study is that the conclusions with respect to ordering and timing of acquisition of mutations may not represent universal patterns of drug resistance emergence in other areas of the globe. CONCLUSIONS: In the first whole genome-based analysis of the emergence of drug resistance among clinical isolates of M. tuberculosis, we show that the ancestral precursor of the LAM4 XDR outbreak strain in Tugela Ferry gained mutations to first-line drugs at the beginning of the antibiotic era. Subsequent accumulation of stepwise resistance mutations, occurring over decades and prior to the explosion of HIV in this region, yielded MDR and XDR, permitting the emergence of compensatory mutations. Our results suggest that drug-resistant strains circulating today reflect not only vulnerabilities of current TB control efforts but also those that date back 50 y. In drug-resistant TB, isoniazid resistance was overwhelmingly the initial resistance mutation to be acquired, which would not be detected by current rapid molecular diagnostics employed in South Africa that assess only rifampicin resistance.

Genome Evolution and Innovation across the Four Major Lineages of Cryptococcus gattii.

UNLABELLED: Cryptococcus gattii is a fungal pathogen of humans, causing pulmonary infections in otherwise healthy hosts. To characterize genomic variation among the four major lineages of C. gattii (VGI, -II, -III, and -IV), we generated, annotated, and compared 16 de novo genome assemblies, including the first for the rarely isolated lineages VGIII and VGIV. By identifying syntenic regions across assemblies, we found 15 structural rearrangements, which were almost exclusive to the VGI-III-IV lineages. Using synteny to inform orthology prediction, we identified a core set of 87% of C. gattii genes present as single copies in all four lineages. Remarkably, 737 genes are variably inherited across lineages and are overrepresented for response to oxidative stress, mitochondrial import, and metal binding and transport. Specifically, VGI has an expanded set of iron-binding genes thought to be important to the virulence of Cryptococcus, while VGII has expansions in the stress-related heat shock proteins relative to the other lineages. We also characterized genes uniquely absent in each lineage, including a copper transporter absent from VGIV, which influences Cryptococcus survival during pulmonary infection and the onset of meningoencephalitis. Through inclusion of population-level data for an additional 37 isolates, we identified a new transcontinental clonal group that we name VGIIx, mitochondrial recombination between VGII and VGIII, and positive selection of multidrug transporters and the iron-sulfur protein aconitase along multiple branches of the phylogenetic tree. Our results suggest that gene expansion or contraction and positive selection have introduced substantial variation with links to mechanisms of pathogenicity across this species complex. IMPORTANCE: The genetic differences between phenotypically different pathogens provide clues to the underlying mechanisms of those traits and can lead to new drug targets and improved treatments for those diseases. In this paper, we compare 16 genomes belonging to four highly differentiated lineages of Cryptococcus gattii, which cause pulmonary infections in otherwise healthy humans and other animals. Half of these lineages have not had their genomes previously assembled and annotated. We identified 15 ancestral rearrangements in the genome and over 700 genes that are unique to one or more lineages, many of which are associated with virulence. In addition, we found evidence for recent transcontinental spread, mitochondrial genetic exchange, and positive selection in multidrug transporters. Our results suggest that gene expansion/contraction and positive selection are diversifying the mechanisms of pathogenicity across this species complex.

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.

Analysis of the genome and transcriptome of Cryptococcus neoformans var. grubii reveals complex RNA expression and microevolution leading to virulence attenuation.

Cryptococcus neoformans is a pathogenic basidiomycetous yeast responsible for more than 600,000 deaths each year. It occurs as two serotypes (A and D) representing two varieties (i.e. grubii and neoformans, respectively). Here, we sequenced the genome and performed an RNA-Seq-based analysis of the C. neoformans var. grubii transcriptome structure. We determined the chromosomal locations, analyzed the sequence/structural features of the centromeres, and identified origins of replication. The genome was annotated based on automated and manual curation. More than 40,000 introns populating more than 99% of the expressed genes were identified. Although most of these introns are located in the coding DNA sequences (CDS), over 2,000 introns in the untranslated regions (UTRs) were also identified. Poly(A)-containing reads were employed to locate the polyadenylation sites of more than 80% of the genes. Examination of the sequences around these sites revealed a new poly(A)-site-associated motif (AUGHAH). In addition, 1,197 miscRNAs were identified. These miscRNAs can be spliced and/or polyadenylated, but do not appear to have obvious coding capacities. Finally, this genome sequence enabled a comparative analysis of strain H99 variants obtained after laboratory passage. The spectrum of mutations identified provides insights into the genetics underlying the micro-evolution of a laboratory strain, and identifies mutations involved in stress responses, mating efficiency, and virulence.

Comparative genomic and transcriptomic analysis of wangiella dermatitidis, a major cause of phaeohyphomycosis and a model black yeast human pathogen.

Black or dark brown (phaeoid) fungi cause cutaneous, subcutaneous, and systemic infections in humans. Black fungi thrive in stressful conditions such as intense light, high radiation, and very low pH. Wangiella (Exophiala) dermatitidis is arguably the most studied phaeoid fungal pathogen of humans. Here, we report our comparative analysis of the genome of W. dermatitidis and the transcriptional response to low pH stress. This revealed that W. dermatitidis has lost the ability to synthesize alpha-glucan, a cell wall compound many pathogenic fungi use to evade the host immune system. In contrast, W. dermatitidis contains a similar profile of chitin synthase genes as related fungi and strongly induces genes involved in cell wall synthesis in response to pH stress. The large portfolio of transporters may provide W. dermatitidis with an enhanced ability to remove harmful products as well as to survive on diverse nutrient sources. The genome encodes three independent pathways for producing melanin, an ability linked to pathogenesis; these are active during pH stress, potentially to produce a barrier to accumulated oxidative damage that might occur under stress conditions. In addition, a full set of fungal light-sensing genes is present, including as part of a carotenoid biosynthesis gene cluster. Finally, we identify a two-gene cluster involved in nucleotide sugar metabolism conserved with a subset of fungi and characterize a horizontal transfer event of this cluster between fungi and algal viruses. This work reveals how W. dermatitidis has adapted to stress and survives in diverse environments, including during human infections.