Using evolutionary constraint to define novel candidate driver genes in medulloblastoma.

Current knowledge of cancer genomics remains biased against noncoding mutations. To systematically search for regulatory noncoding mutations, we assessed mutations in conserved positions in the genome under the assumption that these are more likely to be functional than mutations in positions with low conservation. To this end, we use whole-genome sequencing data from the International Cancer Genome Consortium and combined it with evolutionary constraint inferred from 240 mammals, to identify genes enriched in noncoding constraint mutations (NCCMs), mutations likely to be regulatory in nature. We compare medulloblastoma (MB), which is malignant, to pilocytic astrocytoma (PA), a primarily benign tumor, and find highly different NCCM frequencies between the two, in agreement with the fact that malignant cancers tend to have more mutations. In PA, a high NCCM frequency only affects the BRAF locus, which is the most commonly mutated gene in PA. In contrast, in MB, >500 genes have high levels of NCCMs. Intriguingly, several loci with NCCMs in MB are associated with different ages of onset, such as the HOXB cluster in young MB patients. In adult patients, NCCMs occurred in, e.g., the WASF-2/AHDC1/FGR locus. One of these NCCMs led to increased expression of the SRC kinase FGR and augmented responsiveness of MB cells to dasatinib, a SRC kinase inhibitor. Our analysis thus points to different molecular pathways in different patient groups. These newly identified putative candidate driver mutations may aid in patient stratification in MB and could be valuable for future selection of personalized treatment options.

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.

Evolution of pathogenicity and sexual reproduction in eight Candida genomes.

Candida species are the most common cause of opportunistic fungal infection worldwide. Here we report the genome sequences of six Candida species and compare these and related pathogens and non-pathogens. There are significant expansions of cell wall, secreted and transporter gene families in pathogenic species, suggesting adaptations associated with virulence. Large genomic tracts are homozygous in three diploid species, possibly resulting from recent recombination events. Surprisingly, key components of the mating and meiosis pathways are missing from several species. These include major differences at the mating-type loci (MTL); Lodderomyces elongisporus lacks MTL, and components of the a1/2 cell identity determinant were lost in other species, raising questions about how mating and cell types are controlled. Analysis of the CUG leucine-to-serine genetic-code change reveals that 99% of ancestral CUG codons were erased and new ones arose elsewhere. Lastly, we revise the Candida albicans gene catalogue, identifying many new genes.

Identification of the mating-type (MAT) locus that controls sexual reproduction of Blastomyces dermatitidis.

Blastomyces dermatitidis is a dimorphic fungal pathogen that primarily causes blastomycosis in the midwestern and northern United States and Canada. While the genes controlling sexual development have been known for a long time, the genes controlling sexual reproduction of B. dermatitidis (teleomorph, Ajellomyces dermatitidis) are unknown. We identified the mating-type (MAT) locus in the B. dermatitidis genome by comparative genomic approaches. The B. dermatitidis MAT locus resembles those of other dimorphic fungi, containing either an alpha-box (MAT1-1) or an HMG domain (MAT1-2) gene linked to the APN2, SLA2, and COX13 genes. However, in some strains of B. dermatitidis, the MAT locus harbors transposable elements (TEs) that make it unusually large compared to the MAT locus of other dimorphic fungi. Based on the MAT locus sequences of B. dermatitidis, we designed specific primers for PCR determination of the mating type. Two B. dermatitidis isolates of opposite mating types were cocultured on mating medium. Immature sexual structures were observed starting at 3 weeks of coculture, with coiled-hyphae-containing cleistothecia developing over the next 3 to 6 weeks. Genetic recombination was detected in potential progeny by mating-type determination, PCR-restriction fragment length polymorphism (PCR-RFLP), and random amplification of polymorphic DNA (RAPD) analyses, suggesting that a meiotic sexual cycle might have been completed. The F1 progeny were sexually fertile when tested with strains of the opposite mating type. Our studies provide a model for the evolution of the MAT locus in the dimorphic and closely related fungi and open the door to classic genetic analysis and studies on the possible roles of mating and mating type in infection and virulence.

Comparative genomic analysis of human fungal pathogens causing paracoccidioidomycosis.

Paracoccidioides is a fungal pathogen and the cause of paracoccidioidomycosis, a health-threatening human systemic mycosis endemic to Latin America. Infection by Paracoccidioides, a dimorphic fungus in the order Onygenales, is coupled with a thermally regulated transition from a soil-dwelling filamentous form to a yeast-like pathogenic form. To better understand the genetic basis of growth and pathogenicity in Paracoccidioides, we sequenced the genomes of two strains of Paracoccidioides brasiliensis (Pb03 and Pb18) and one strain of Paracoccidioides lutzii (Pb01). These genomes range in size from 29.1 Mb to 32.9 Mb and encode 7,610 to 8,130 genes. To enable genetic studies, we mapped 94% of the P. brasiliensis Pb18 assembly onto five chromosomes. We characterized gene family content across Onygenales and related fungi, and within Paracoccidioides we found expansions of the fungal-specific kinase family FunK1. Additionally, the Onygenales have lost many genes involved in carbohydrate metabolism and fewer genes involved in protein metabolism, resulting in a higher ratio of proteases to carbohydrate active enzymes in the Onygenales than their relatives. To determine if gene content correlated with growth on different substrates, we screened the non-pathogenic onygenale Uncinocarpus reesii, which has orthologs for 91% of Paracoccidioides metabolic genes, for growth on 190 carbon sources. U. reesii showed growth on a limited range of carbohydrates, primarily basic plant sugars and cell wall components; this suggests that Onygenales, including dimorphic fungi, can degrade cellulosic plant material in the soil. In addition, U. reesii grew on gelatin and a wide range of dipeptides and amino acids, indicating a preference for proteinaceous growth substrates over carbohydrates, which may enable these fungi to also degrade animal biomass. These capabilities for degrading plant and animal substrates suggest a duality in lifestyle that could enable pathogenic species of Onygenales to transfer from soil to animal hosts.

Obligate biotrophy features unraveled by the genomic analysis of rust fungi.

Rust fungi are some of the most devastating pathogens of crop plants. They are obligate biotrophs, which extract nutrients only from living plant tissues and cannot grow apart from their hosts. Their lifestyle has slowed the dissection of molecular mechanisms underlying host invasion and avoidance or suppression of plant innate immunity. We sequenced the 101-Mb genome of Melampsora larici-populina, the causal agent of poplar leaf rust, and the 89-Mb genome of Puccinia graminis f. sp. tritici, the causal agent of wheat and barley stem rust. We then compared the 16,399 predicted proteins of M. larici-populina with the 17,773 predicted proteins of P. graminis f. sp tritici. Genomic features related to their obligate biotrophic lifestyle include expanded lineage-specific gene families, a large repertoire of effector-like small secreted proteins, impaired nitrogen and sulfur assimilation pathways, and expanded families of amino acid and oligopeptide membrane transporters. The dramatic up-regulation of transcripts coding for small secreted proteins, secreted hydrolytic enzymes, and transporters in planta suggests that they play a role in host infection and nutrient acquisition. Some of these genomic hallmarks are mirrored in the genomes of other microbial eukaryotes that have independently evolved to infect plants, indicating convergent adaptation to a biotrophic existence inside plant cells.

Genomic analysis across 53 canine cancer types reveals novel mutations and high clinical actionability potential.

A genomic understanding of the oncogenic processes and individual variability of human cancer has steadily fueled improvement in patient outcomes over the past 20 years. Mutations within tumour tissues are routinely assessed through clinical genomic diagnostic assays by academic and commercial laboratories to facilitate diagnosis, prognosis and effective treatment stratification. The application of genomics has unveiled a wealth of mutation-based biomarkers in canine cancers, suggesting that the transformative principles that have revolutionized human cancer medicine can be brought to bear in veterinary oncology. To advance clinical genomics and genomics-guided medicine in canine oncology, we have developed and validated a canine cancer next-generation sequencing gene panel for the identification of multiple mutation types in clinical specimens. With this panel, we examined the genomic landscapes of 828 tumours from 813 dogs, spanning 53 cancer types. We identified 7856 alterations, encompassing copy number variants, single nucleotide variants, indels and internal tandem duplications. Additionally, we evaluated the clinical utility of these alterations by incorporating a biomarker framework from comprehensive curation of primary canine literature and inferences from human cancer genomic biomarker literature and clinical diagnostics. Remarkably, nearly 90% of the cases exhibited mutations with diagnostic, prognostic or therapeutic implications. Our work represents a thorough assessment of genomic landscapes in a large cohort of canine cancers, the first of its kind for its comprehensive inclusion of multiple mutation types and structured annotation of biomarkers, demonstrating the clinical potential of leveraging mutation-based biomarkers in veterinary oncology.

PIK3CA is recurrently mutated in canine mammary tumors, similarly to in human mammary neoplasia.

Biological features of neoplastic disease affecting mammary gland tissue are shared between canines and humans. Research performed in either species has translational value and early phase clinical trials performed in canines with spontaneous disease could be informative for human trials. The purpose of this study was to investigate the somatic genetic aberrations occurring in canine mammary neoplasia by exome capture and next generation sequencing. Based on 55 tumor-normal pairs we identified the PIK3CA gene as the most commonly mutated gene in canine mammary tumors, with 25% of samples carrying mutations in this gene. A recurrent missense mutation was identified, p.H1047R, which is homologous to the human PIK3CA hotspot mutation found in different types of breast neoplasia. Mutations homologous to other known human mutation hotspots such as the PIK3CA p.E545K and the KRAS p.G12V/D were also identified. We identified copy number aberrations affecting important tumor suppressor and oncogenic pathways including deletions affecting the PTEN tumor suppressor gene. We suggest that activation of the KRAS or PIK3CA oncogenes or loss of the PTEN suppressor gene may be important for mammary tumor development in dogs. This data endorses the conservation of cancer across species and the validity of studying cancer in non-human species.

Genomic tumor analysis provides clinical guidance for the management of diagnostically challenging cancers in dogs.

OBJECTIVE: To evaluate the diagnostic, prognostic, and therapeutic utility of a cancer genomic diagnostic assay (SearchLight DNA; Vidium Animal Health) for diagnostically ambiguous cancer cases. ANIMALS: 69 privately owned dogs with ambiguous cancer diagnoses and for which the genomic assay was performed. PROCEDURES: Genomic assay reports generated between September 28, 2020, and July 31, 2022, for dogs with malignancy or suspected malignancy were reviewed to determine the assay's clinical utility defined as providing diagnostic clarity, prognostic information, and/or therapeutic options. RESULTS: Genomic analysis provided diagnostic clarity in 37 of 69 cases (54%; group 1) and therapeutic and/or prognostic information in 22 of the remaining 32 cases (69%; group 2) for which the diagnosis remained elusive. Overall, the genomic assay was clinically useful in 86% (59/69) of cases. CLINICAL RELEVANCE: To our knowledge, this was the first study to evaluate the multifaceted clinical utility of a single cancer genomic test in veterinary medicine. Study findings supported the use of tumor genomic testing for dogs with cancer, particularly those that are diagnostically ambiguous and therefore inherently challenging to manage. This evidence-driven genomic assay provided diagnostic guidance, prognostic support, and therapeutic options for most patients with an unclear cancer diagnosis that would otherwise have an unsubstantiated clinical plan. Furthermore, 38% (26/69) of samples were easily obtained aspirates. Sample factors (sample type, percentage of tumor cells, and number of mutations) did not influence diagnostic yield. Our study demonstrated the value of genomic testing for the management of canine cancer.

Standing in the canine precision medicine knowledge gap: Improving annotation of canine cancer genomic biomarkers through systematic comparative analysis of human cancer mutations in COSMIC.

The accrual of cancer mutation data and related functional and clinical associations have revolutionised human oncology, enabling the advancement of precision medicine and biomarker-guided clinical management. The catalogue of cancer mutations is also growing in canine cancers. However, without direct high-powered functional data in dogs, it remains challenging to interpret and utilise them in research and clinical settings. It is well-recognised that canine and human cancers share genetic, molecular and phenotypic similarities. Therefore, leveraging the massive wealth of human mutation data may help advance canine oncology. Here, we present a structured analysis of sequence conservation and conversion of human mutations to the canine genome through a 'caninisation' process. We applied this analysis to COSMIC, the Catalogue of Somatic Mutations in Cancer, the most prominent human cancer mutation database. For the project's initial phase, we focused on the subset of the COSMIC data corresponding to Cancer Gene Census (CGC) genes. A total of 670 canine orthologs were found for 721 CGC genes. In these genes, 365 K unique mutations across 160 tumour types were converted successfully to canine coordinates. We identified shared putative cancer-driving mutations, including pathogenic and hotspot mutations and mutations bearing similar biomarker associations with diagnostic, prognostic and therapeutic utility. Thus, this structured caninisation of human cancer mutations facilitates the interpretation and annotation of canine mutations and helps bridge the knowledge gap to enable canine precision medicine.

Novel genomic prognostic biomarkers for dogs with cancer.

BACKGROUND: Growing evidence from dogs and humans supports the abundance of mutation-based biomarkers in tumors of dogs. Increasing the use of clinical genomic diagnostic testing now provides another powerful data source for biomarker discovery. HYPOTHESIS: Analyzed clinical outcomes in dogs with cancer profiled using SearchLight DNA, a cancer gene panel for dogs, to identify mutations with prognostic value. ANIMALS: A total of 127 cases of cancer in dogs were analyzed using SearchLight DNA and for which clinical outcome information was available. METHODS: Clinical data points were collected by medical record review. Variables including mutated genes, mutations, signalment, and treatment were fitted using Cox proportional hazard models to identify factors associated with progression-free survival (PFS). The log-rank test was used to compare PFS between patients receiving and not receiving targeted treatment before first progression. RESULTS: Combined genomic and outcomes analysis identified 336 unique mutations in 89 genes across 26 cancer types. Mutations in 6 genes (CCND1, CCND3, SMARCB1, FANCG, CDKN2A/B, and MSH6) were significantly associated with shorter PFS. Dogs that received targeted treatment before first progression (n = 45) experienced significantly longer PFS compared with those that did not (n = 82, P = .01). This significance held true for 29 dogs that received genomically informed targeted treatment compared with those that did not (P = .05). CONCLUSION AND CLINICAL IMPORTANCE: We identified novel mutations with prognostic value and demonstrate the benefit of targeted treatment across multiple cancer types. These results provide clinical evidence of the potential for genomics and precision medicine in dogs with cancer.

Leveraging Base Pair Mammalian Constraint to Understand Genetic Variation and Human Disease.

Although thousands of genomic regions have been associated with heritable human diseases, attempts to elucidate biological mechanisms are impeded by a general inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function that is agnostic to cell type or disease mechanism. Here, single base phyloP scores from the whole genome alignment of 240 placental mammals identified 3.5% of the human genome as significantly constrained, and likely functional. We compared these scores to large-scale genome annotation, genome-wide association studies (GWAS), copy number variation, clinical genetics findings, and cancer data sets. Evolutionarily constrained positions are enriched for variants explaining common disease heritability (more than any other functional annotation). Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

Leveraging base-pair mammalian constraint to understand genetic variation and human disease.

Thousands of genomic regions have been associated with heritable human diseases, but attempts to elucidate biological mechanisms are impeded by an inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function, agnostic to cell type or disease mechanism. Single-base phyloP scores from 240 mammals identified 3.3% of the human genome as significantly constrained and likely functional. We compared phyloP scores to genome annotation, association studies, copy-number variation, clinical genetics findings, and cancer data. Constrained positions are enriched for variants that explain common disease heritability more than other functional annotations. Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

Whole-genome sequencing of multiple isolates of Puccinia triticina reveals asexual lineages evolving by recurrent mutations.

The wheat leaf rust fungus, Puccinia triticina Erikss., is a worldwide pathogen of tetraploid durum and hexaploid wheat. Many races of P. triticina differ for virulence to specific leaf rust resistance genes and are found in most wheat-growing regions of the world. Wheat cultivars with effective leaf rust resistance exert selection pressure on P. triticina populations for virulent race types. The objectives of this study were to examine whole-genome sequence data of 121 P. triticina isolates and to gain insight into race evolution. The collection included isolates comprising of many different race phenotypes collected worldwide from common and durum wheat. One isolate from wild wheat relative Aegilops speltoides and two from Ae. cylindrica were also included for comparison. Based on 121,907 informative variants identified relative to the reference Race 1-1 genome, isolates were clustered into 11 major lineages with 100% bootstrap support. The isolates were also grouped based on variation in 1311 predicted secreted protein genes. In gene-coding regions, all groups had high ratios of nonsynonymous to synonymous mutations and nonsense to readthrough mutations. Grouping of isolates based on two main variation principle components for either genome-wide variation or variation just within the secreted protein genes, indicated similar groupings. Variants were distributed across the entire genome, not just within the secreted protein genes. Our results suggest that recurrent mutation and selection play a major role in differentiation within the clonal lineages.

Whole-genome sequencing of glioblastoma reveals enrichment of non-coding constraint mutations in known and novel genes.

BACKGROUND: Glioblastoma (GBM) has one of the worst 5-year survival rates of all cancers. While genomic studies of the disease have been performed, alterations in the non-coding regulatory regions of GBM have largely remained unexplored. We apply whole-genome sequencing (WGS) to identify non-coding mutations, with regulatory potential in GBM, under the hypothesis that regions of evolutionary constraint are likely to be functional, and somatic mutations are likely more damaging than in unconstrained regions. RESULTS: We validate our GBM cohort, finding similar copy number aberrations and mutated genes based on coding mutations as previous studies. Performing analysis on non-coding constraint mutations and their position relative to nearby genes, we find a significant enrichment of non-coding constraint mutations in the neighborhood of 78 genes that have previously been implicated in GBM. Among them, SEMA3C and DYNC1I1 show the highest frequencies of alterations, with multiple mutations overlapping transcription factor binding sites. We find that a non-coding constraint mutation in the SEMA3C promoter reduces the DNA binding capacity of the region. We also identify 1776 other genes enriched for non-coding constraint mutations with likely regulatory potential, providing additional candidate GBM genes. The mutations in the top four genes, DLX5, DLX6, FOXA1, and ISL1, are distributed over promoters, UTRs, and multiple transcription factor binding sites. CONCLUSIONS: These results suggest that non-coding constraint mutations could play an essential role in GBM, underscoring the need to connect non-coding genomic variation to biological function and disease pathology.

Comparative Genomics Reveals Shared Mutational Landscape in Canine Hemangiosarcoma and Human Angiosarcoma.

Angiosarcoma is a highly aggressive cancer of blood vessel-forming cells with few effective treatment options and high patient mortality. It is both rare and heterogenous, making large, well-powered genomic studies nearly impossible. Dogs commonly suffer from a similar cancer, called hemangiosarcoma, with breeds like the golden retriever carrying heritable genetic factors that put them at high risk. If the clinical similarity of canine hemangiosarcoma and human angiosarcoma reflects shared genomic etiology, dogs could be a critically needed model for advancing angiosarcoma research. We assessed the genomic landscape of canine hemangiosarcoma via whole-exome sequencing (47 golden retriever hemangiosarcomas) and RNA sequencing (74 hemangiosarcomas from multiple breeds). Somatic coding mutations occurred most frequently in the tumor suppressor TP53 (59.6% of cases) as well as two genes in the PI3K pathway: the oncogene PIK3CA (29.8%) and its regulatory subunit PIK3R1 (8.5%). The predominant mutational signature was the age-associated deamination of cytosine to thymine. As reported in human angiosarcoma, CDKN2A/B was recurrently deleted and VEGFA, KDR, and KIT recurrently gained. We compared the canine data to human data recently released by The Angiosarcoma Project, and found many of the same genes and pathways significantly enriched for somatic mutations, particularly in breast and visceral angiosarcomas. Canine hemangiosarcoma closely models the genomic landscape of human angiosarcoma of the breast and viscera, and is a powerful tool for investigating the pathogenesis of this devastating disease. IMPLICATIONS: We characterize the genomic landscape of canine hemangiosarcoma and demonstrate its similarity to human angiosarcoma.

SETD2 Is Recurrently Mutated in Whole-Exome Sequenced Canine Osteosarcoma.

Osteosarcoma is a debilitating bone cancer that affects humans, especially children and adolescents. A homologous form of osteosarcoma spontaneously occurs in dogs, and its differential incidence observed across breeds allows for the investigation of tumor mutations in the context of multiple genetic backgrounds. Using whole-exome sequencing and dogs from three susceptible breeds (22 golden retrievers, 21 Rottweilers, and 23 greyhounds), we found that osteosarcoma tumors show a high frequency of somatic copy-number alterations (SCNA), affecting key oncogenes and tumor-suppressor genes. The across-breed results are similar to what has been observed for human osteosarcoma, but the disease frequency and somatic mutation counts vary in the three breeds. For all breeds, three mutational signatures (one of which has not been previously reported) and 11 significantly mutated genes were identified. TP53 was the most frequently altered gene (83% of dogs have either mutations or SCNA in TP53), recapitulating observations in human osteosarcoma. The second most frequently mutated gene, histone methyltransferase SETD2, has known roles in multiple cancers, but has not previously been strongly implicated in osteosarcoma. This study points to the likely importance of histone modifications in osteosarcoma and highlights the strong genetic similarities between human and dog osteosarcoma, suggesting that canine osteosarcoma may serve as an excellent model for developing treatment strategies in both species.Significance: Canine osteosarcoma genomics identify SETD2 as a possible oncogenic driver of osteosarcoma, and findings establish the canine model as a useful comparative model for the corresponding human disease. Cancer Res; 78(13); 3421-31. ©2018 AACR.

Comparative Genomics Integrated with Association Analysis Identifies Candidate Effector Genes Corresponding to Lr20 in Phenotype-Paired Puccinia triticina Isolates from Australia.

Leaf rust is one of the most common and damaging diseases of wheat, and is caused by an obligate biotrophic basidiomycete, Puccinia triticina (Pt). In the present study, 20 Pt isolates from Australia, comprising 10 phenotype-matched pairs with contrasting pathogenicity for Lr20, were analyzed using whole genome sequencing. Compared to the reference genome of the American Pt isolate 1-1 BBBD Race 1, an average of 404,690 single nucleotide polymorphisms (SNPs) per isolate was found and the proportion of heterozygous SNPs was above 87% in the majority of the isolates, demonstrating a high level of polymorphism and a high rate of heterozygosity. From the genome-wide SNPs, a phylogenetic tree was inferred, which consisted of a large clade of 15 isolates representing diverse presumed clonal lineages including 14 closely related isolates and the more diverged isolate 670028, and a small clade of five isolates characterized by lower heterozygosity level. Principle component analysis detected three distinct clusters, corresponding exactly to the two major subsets of the small clade and the large clade comprising all 15 isolates without further separation of isolate 670028. While genome-wide association analysis identified 302 genes harboring at least one SNP associated with Lr20 virulence (p < 0.05), a Wilcoxon rank sum test revealed that 36 and 68 genes had significant (p < 0.05) and marginally significant (p < 0.1) differences in the counts of non-synonymous mutations between Lr20 avirulent and virulent groups, respectively. Twenty of these genes were predicted to have a signal peptide without a transmembrane segment, and hence identified as candidate effector genes corresponding to Lr20. SNP analysis also implicated the potential involvement of epigenetics and small RNA in Pt pathogenicity. Future studies are thus warranted to investigate the biological functions of the candidate effectors as well as the gene regulation mechanisms at epigenetic and post-transcription levels. Our study is the first to integrate phenotype-genotype association with effector prediction in Pt genomes, an approach that may circumvent some of the technical difficulties in working with obligate rust fungi and accelerate avirulence gene identification.

Microevolution of Serial Clinical Isolates of Cryptococcus neoformans var. grubii and C. gattii.

The pathogenic species of Cryptococcus are a major cause of mortality owing to severe infections in immunocompromised as well as immunocompetent individuals. Although antifungal treatment is usually effective, many patients relapse after treatment, and in such cases, comparative analyses of the genomes of incident and relapse isolates may reveal evidence of determinative, microevolutionary changes within the host. Here, we analyzed serial isolates cultured from cerebrospinal fluid specimens of 18 South African patients with recurrent cryptococcal meningitis. The time between collection of the incident isolates and collection of the relapse isolates ranged from 124 days to 290 days, and the analyses revealed that, during this period within the patients, the isolates underwent several genetic and phenotypic changes. Considering the vast genetic diversity of cryptococcal isolates in sub-Saharan Africa, it was not surprising to find that the relapse isolates had acquired different genetic and correlative phenotypic changes. They exhibited various mechanisms for enhancing virulence, such as growth at 39°C, adaptation to stress, and capsule production; a remarkable amplification of ERG11 at the native and unlinked locus may provide stable resistance to fluconazole. Our data provide a deeper understanding of the microevolution of Cryptococcus species under pressure from antifungal chemotherapy and host immune responses. This investigation clearly suggests a promising strategy to identify novel targets for improved diagnosis, therapy, and prognosis.IMPORTANCE Opportunistic infections caused by species of the pathogenic yeast Cryptococcus lead to chronic meningoencephalitis and continue to ravage thousands of patients with HIV/AIDS. Despite receiving antifungal treatment, over 10% of patients develop recurrent disease. In this study, we collected isolates of Cryptococcus from cerebrospinal fluid specimens of 18 patients at the time of their diagnosis and when they relapsed several months later. We then sequenced and compared the genomic DNAs of each pair of initial and relapse isolates. We also tested the isolates for several key properties related to cryptococcal virulence as well as for their susceptibility to the antifungal drug fluconazole. These analyses revealed that the relapsing isolates manifested multiple genetic and chromosomal changes that affected a variety of genes implicated in the pathogenicity of Cryptococcus or resistance to fluconazole. This application of comparative genomics to serial clinical isolates provides a blueprint for identifying the mechanisms whereby pathogenic microbes adapt within patients to prolong disease.

Comparative Analysis Highlights Variable Genome Content of Wheat Rusts and Divergence of the Mating Loci.

Three members of the Puccinia genus, Pucciniatriticina (Pt), Pstriiformis f.sp. tritici (Pst), and Pgraminis f.sp. tritici (Pgt), cause the most common and often most significant foliar diseases of wheat. While similar in biology and life cycle, each species is uniquely adapted and specialized. The genomes of Pt and Pst were sequenced and compared to that of Pgt to identify common and distinguishing gene content, to determine gene variation among wheat rust pathogens, other rust fungi, and basidiomycetes, and to identify genes of significance for infection. Pt had the largest genome of the three, estimated at 135 Mb with expansion due to mobile elements and repeats encompassing 50.9% of contig bases; in comparison, repeats occupy 31.5% for Pst and 36.5% for Pgt We find all three genomes are highly heterozygous, with Pst [5.97 single nucleotide polymorphisms (SNPs)/kb] nearly twice the level detected in Pt (2.57 SNPs/kb) and that previously reported for Pgt Of 1358 predicted effectors in Pt, 784 were found expressed across diverse life cycle stages including the sexual stage. Comparison to related fungi highlighted the expansion of gene families involved in transcriptional regulation and nucleotide binding, protein modification, and carbohydrate degradation enzymes. Two allelic homeodomain pairs, HD1 and HD2, were identified in each dikaryotic Puccinia species along with three pheromone receptor (STE3) mating-type genes, two of which are likely representing allelic specificities. The HD proteins were active in a heterologous Ustilago maydis mating assay and host-induced gene silencing (HIGS) of the HD and STE3 alleles reduced wheat host infection.