Biomedical subjects

Jason E Stajich

Publications and source records attributed to Jason E Stajich.

At least 19 recordsLinked to original sources

PubMed2026

Mycodnaviridae is a clade of giant viruses that persistently infect zoosporic fungi.

Giant viruses of the phylum Nucleocytoviricota have emerged as particularly notable due to their increasingly recognized impacts on eukaryotic genome evolution. Their origins are hypothesized to predate or coincide with the diversification of eukaryotes, and they have been detected in hosts that span the eukaryotic tree of life. But surprisingly, such viruses have not been definitively found in Kingdom Fungi, though earlier genomic and metagenomic work suggests putative associations. Here we report both "viral fossils" and active infection by giant viruses in fungi, particularly in the zoosporic phyla Blastocladiomycota and Chytridiomycota. The recovered viral assemblies span up to 350 kb, encode over 300 genes, and form a monophyletic family-level clade within the Nucleocytoviricota related to orders Imitervirales and Algavirales, which we name Mycodnaviridae. We observed variation in infection status among the isolates including apparent active infection and transcriptionally suppressed states, suggesting that viral activation may be constrained to certain life stages of the host. Our experimental findings add to the limited natural virus-host systems available in culture for the study of giant viruses and expand the known host range of Nucleocytoviricota into a new kingdom that contains many model species. Mycodnaviridae have a global distribution, which invites inquiry into the implications of these infections for host traits, host genome evolution, and the metabolic impacts on ecosystems.

Giant Viruses↗

PubMed2026

Deep subsurface organic-rich shale supports abundant, diverse, and novel fungi.

As Earth's principal reservoir of organic carbon and microbial biomass, the deep subsurface hosts microorganisms capable of mobilizing this once-sequestered carbon. Contrary to standard assumptions of eukaryotic scarcity, this study documents abundant fungal communities, ranging from 4.2 × 103 to 6.8 × 103 fungal cells ml-1, across a methane-producing organic-rich shale 247-556 meters below the surface. Although fungal:bacterial cell ratios ranged from 1:7028 to 1:713, application of biomass conversion factors developed for oceanic systems yielded a median fungal:bacterial biomass ratio of 1:4.7. 16S ribosomal ribonucleic acid (rRNA) gene amplicons revealed bacterial and archaeal communities mirroring those found in well-characterized extremophilic, carbon-degrading environments, while sequencing of 18S rRNA gene and internal transcribed spacer rRNA spacer amplicons collectively identified a eukaryotic hotspot with 689 fungal operational taxonomic units across six phyla. The dominant fungal classes, Agaricomycetes and Dothideomycetes, are well-established degraders of recalcitrant carbon compounds at the surface, suggesting they may similarly contribute to organic matter degradation and ecosystem maintenance in the subsurface. Cultivation and isolation efforts yielded 205 fungal strains, including 13 candidate novel taxa, underscoring the deep subsurface as an underexplored eukaryotic habitat. Stable carbon isotopes indicate methane is predominantly generated via microbial conversion of the fossil carbon, while water isotopes suggest in situ geochemical conditions have been relatively stable since the Late Pleistocene, with subglacial recharge as a plausible mechanism for microbial introduction. Collectively, these findings suggest that fungi are underrecognized contributors to organic matter transformation and functional diversity in the deep biosphere, revealing a critical gap in our understanding of deep subsurface ecosystem processes.

Fungi↗

PubMed2006

The evolution of mammalian gene families.

Gene families are groups of homologous genes that are likely to have highly similar functions. Differences in family size due to lineage-specific gene duplication and gene loss may provide clues to the evolutionary forces that have shaped mammalian genomes. Here we analyze the gene families contained within the whole genomes of human, chimpanzee, mouse, rat, and dog. In total we find that more than half of the 9,990 families present in the mammalian common ancestor have either expanded or contracted along at least one lineage. Additionally, we find that a large number of families are completely lost from one or more mammalian genomes, and a similar number of gene families have arisen subsequent to the mammalian common ancestor. Along the lineage leading to modern humans we infer the gain of 689 genes and the loss of 86 genes since the split from chimpanzees, including changes likely driven by adaptive natural selection. Our results imply that humans and chimpanzees differ by at least 6% (1,418 of 22,000 genes) in their complement of genes, which stands in stark contrast to the oft-cited 1.5% difference between orthologous nucleotide sequences. This genomic "revolving door" of gene gain and loss represents a large number of genetic differences separating humans from our closest relatives.

Animals↗

PubMed2006

BASC: an integrated bioinformatics system for Brassica research.

The BASC system provides tools for the integrated mining and browsing of genetic, genomic and phenotypic data. This public resource hosts information on Brassica species supporting the Multinational Brassica Genome Sequencing Project, and is based upon five distinct modules, ESTDB, Microarray, MarkerQTL, CMap and EnsEMBL. ESTDB hosts expressed gene sequences and related annotation derived from comparison with GenBank, UniRef and the genome sequence of Arabidopsis. The Microarray module hosts gene expression information related to genes annotated within ESTDB. MarkerQTL is the most complex module and integrates information on genetic markers, maps, individuals, genotypes and traits. Two further modules include an Arabidopsis EnsEMBL genome viewer and the CMap comparative genetic map viewer for the visualization and integration of genetic and genomic data. The database is accessible at http://bioinformatics.pbcbasc.latrobe.edu.au.

Arabidopsis↗

PubMed2006

A fungal phylogeny based on 42 complete genomes derived from supertree and combined gene analysis.

BACKGROUND: To date, most fungal phylogenies have been derived from single gene comparisons, or from concatenated alignments of a small number of genes. The increase in fungal genome sequencing presents an opportunity to reconstruct evolutionary events using entire genomes. As a tool for future comparative, phylogenomic and phylogenetic studies, we used both supertrees and concatenated alignments to infer relationships between 42 species of fungi for which complete genome sequences are available. RESULTS: A dataset of 345,829 genes was extracted from 42 publicly available fungal genomes. Supertree methods were employed to derive phylogenies from 4,805 single gene families. We found that the average consensus supertree method may suffer from long-branch attraction artifacts, while matrix representation with parsimony (MRP) appears to be immune from these. A genome phylogeny was also reconstructed from a concatenated alignment of 153 universally distributed orthologs. Our MRP supertree and concatenated phylogeny are highly congruent. Within the Ascomycota, the sub-phyla Pezizomycotina and Saccharomycotina were resolved. Both phylogenies infer that the Leotiomycetes are the closest sister group to the Sordariomycetes. There is some ambiguity regarding the placement of Stagonospora nodurum, the sole member of the class Dothideomycetes present in the dataset. Within the Saccharomycotina, a monophyletic clade containing organisms that translate CTG as serine instead of leucine is evident. There is also strong support for two groups within the CTG clade, one containing the fully sexual species Candida lusitaniae, Candida guilliermondii and Debaryomyces hansenii, and the second group containing Candida albicans, Candida dubliniensis, Candida tropicalis, Candida parapsilosis and Lodderomyces elongisporus. The second major clade within the Saccharomycotina contains species whose genomes have undergone a whole genome duplication (WGD), and their close relatives. We could not confidently resolve whether Candida glabrata or Saccharomyces castellii lies at the base of the WGD clade. CONCLUSION: We have constructed robust phylogenies for fungi based on whole genome analysis. Overall, our phylogenies provide strong support for the classification of phyla, sub-phyla, classes and orders. We have resolved the relationship of the classes Leotiomyctes and Sordariomycetes, and have identified two classes within the CTG clade of the Saccharomycotina that may correlate with sexual status.

Ascomycota↗

PubMed2006

Insights from the genome of the biotrophic fungal plant pathogen Ustilago maydis.

Ustilago maydis is a ubiquitous pathogen of maize and a well-established model organism for the study of plant-microbe interactions. This basidiomycete fungus does not use aggressive virulence strategies to kill its host. U. maydis belongs to the group of biotrophic parasites (the smuts) that depend on living tissue for proliferation and development. Here we report the genome sequence for a member of this economically important group of biotrophic fungi. The 20.5-million-base U. maydis genome assembly contains 6,902 predicted protein-encoding genes and lacks pathogenicity signatures found in the genomes of aggressive pathogenic fungi, for example a battery of cell-wall-degrading enzymes. However, we detected unexpected genomic features responsible for the pathogenicity of this organism. Specifically, we found 12 clusters of genes encoding small secreted proteins with unknown function. A significant fraction of these genes exists in small gene families. Expression analysis showed that most of the genes contained in these clusters are regulated together and induced in infected tissue. Deletion of individual clusters altered the virulence of U. maydis in five cases, ranging from a complete lack of symptoms to hypervirulence. Despite years of research into the mechanism of pathogenicity in U. maydis, no 'true' virulence factors had been previously identified. Thus, the discovery of the secreted protein gene clusters and the functional demonstration of their decisive role in the infection process illuminate previously unknown mechanisms of pathogenicity operating in biotrophic fungi. Genomic analysis is, similarly, likely to open up new avenues for the discovery of virulence determinants in other pathogens.

Fungal Proteins↗

PubMed2006

Reconstructing the early evolution of Fungi using a six-gene phylogeny.

The ancestors of fungi are believed to be simple aquatic forms with flagellated spores, similar to members of the extant phylum Chytridiomycota (chytrids). Current classifications assume that chytrids form an early-diverging clade within the kingdom Fungi and imply a single loss of the spore flagellum, leading to the diversification of terrestrial fungi. Here we develop phylogenetic hypotheses for Fungi using data from six gene regions and nearly 200 species. Our results indicate that there may have been at least four independent losses of the flagellum in the kingdom Fungi. These losses of swimming spores coincided with the evolution of new mechanisms of spore dispersal, such as aerial dispersal in mycelial groups and polar tube eversion in the microsporidia (unicellular forms that lack mitochondria). The enigmatic microsporidia seem to be derived from an endoparasitic chytrid ancestor similar to Rozella allomycis, on the earliest diverging branch of the fungal phylogenetic tree.

Chytridiomycota↗

PubMed2006

Open source tools and toolkits for bioinformatics: significance, and where are we?

This review summarizes important work in open-source bioinformatics software that has occurred over the past couple of years. The survey is intended to illustrate how programs and toolkits whose source code has been developed or released under an Open Source license have changed informatics-heavy areas of life science research. Rather than creating a comprehensive list of all tools developed over the last 2-3 years, we use a few selected projects encompassing toolkit libraries, analysis tools, data analysis environments and interoperability standards to show how freely available and modifiable open-source software can serve as the foundation for building important applications, analysis workflows and resources.

Algorithms↗

PubMed2006

Phylogenomic analysis of non-ribosomal peptide synthetases in the genus Aspergillus.

Fungi from the genus Aspergillus are important saprophytes and opportunistic human fungal pathogens that contribute in these and other diverse ways to human well-being. Part of their impact on human well-being stems from the production of small molecular weight secondary metabolites, which may contribute to the ability of these fungi to cause invasive fungal infections and allergic diseases. In this study, we identified one group of enzymes responsible for secondary metabolite production in five Aspergillus species, the non-ribosomal peptide synthetases (NRPS). Hidden Markov models were used to search the genome databases of A. fumigatus, A. flavus, A. terreus, A. nidulans, and A. oryzae for domains conserved in NRPS proteins. A genealogy of adenylation domains was utilized to identify orthologous and unique NRPS among the Aspergillus species examined, as well as gain an understanding of the potential evolution of Aspergillus NRPS. mRNA abundance of the 14 NRPS identified in the A. fumigatus genome was analyzed using real-time reverse transcriptase PCR in different environmental conditions to gain a preliminary understanding of the possible functions of the NRPSs' peptide products. Our results suggest that Aspergillus species contain conserved and unique NRPS genes with a complex evolutionary history. This result suggests that the genus Aspergillus produces a substantial diversity of non-ribosomally synthesized peptides. Further analysis of these genes and their peptide products may identify important roles for secondary metabolites produced by NRPS in Aspergillus physiology, ecology, and fungal pathogenicity.

Amino Acid Sequence↗

PubMed2006

Comparative analysis of Saccharomyces cerevisiae WW domains and their interacting proteins.

BACKGROUND: The WW domain is found in a large number of eukaryotic proteins implicated in a variety of cellular processes. WW domains bind proline-rich protein and peptide ligands, but the protein interaction partners of many WW domain-containing proteins in Saccharomyces cerevisiae are largely unknown. RESULTS: We used protein microarray technology to generate a protein interaction map for 12 of the 13 WW domains present in proteins of the yeast S. cerevisiae. We observed 587 interactions between these 12 domains and 207 proteins, most of which have not previously been described. We analyzed the representation of functional annotations within the network, identifying enrichments for proteins with peroxisomal localization, as well as for proteins involved in protein turnover and cofactor biosynthesis. We compared orthologs of the interacting proteins to identify conserved motifs known to mediate WW domain interactions, and found substantial evidence for the structural conservation of such binding motifs throughout the yeast lineages. The comparative approach also revealed that several of the WW domain-containing proteins themselves have evolutionarily conserved WW domain binding sites, suggesting a functional role for inter- or intramolecular association between proteins that harbor WW domains. On the basis of these results, we propose a model for the tuning of interactions between WW domains and their protein interaction partners. CONCLUSION: Protein microarrays provide an appealing alternative to existing techniques for the construction of protein interaction networks. Here we built a network composed of WW domain-protein interactions that illuminates novel features of WW domain-containing proteins and their protein interaction partners.

Amino Acid Motifs↗

PubMed2006

The Cryptococcus neoformans catalase gene family and its role in antioxidant defense.

In the present study, we sought to elucidate the contribution of the Cryptococcus neoformans catalase gene family to antioxidant defense. We employed bioinformatics techniques to identify four members of the C. neoformans catalase gene family and created mutants lacking single or multiple catalase genes. Based on a phylogenetic analysis, CAT1 and CAT3 encode putative spore-specific catalases, CAT2 encodes a putative peroxisomal catalase, and CAT4 encodes a putative cytosolic catalase. Only Cat1 exhibited detectable biochemical activity in vitro, and Cat1 activity was constitutive in the yeast form of this organism. Although they were predicted to be important in spores, neither CAT1 nor CAT3 was essential for mating or spore viability. Consistent with previous studies of Saccharomyces cerevisiae, the single (cat1, cat2, cat3, and cat4) and quadruple (cat1 cat2 cat3 cat4) catalase mutant strains exhibited no oxidative-stress phenotypes under conditions in which either exogenous or endogenous levels of reactive oxygen species were elevated. In addition, there were no significant differences in the mean times to mortality between groups of mice infected with C. neoformans catalase mutant strains (the cat1 and cat1 cat2 cat3 cat4 mutants) and those infected with wild-type strain H99. We conclude from the results of this study that C. neoformans possesses a robust antioxidant system, composed of functionally overlapping and compensatory components that provide protection against endogenous and exogenous oxidative stresses.

Animals↗

PubMed2006

Evidence of mRNA-mediated intron loss in the human-pathogenic fungus Cryptococcus neoformans.

Introns are a defining feature of eukaryotic genomes, though the mechanism of intron gain or loss is not well understood. Reverse transcription of mRNA followed by homologous recombination with the genome has been posited as a mechanism of intron loss, though little direct evidence of recent loss events has been described to support this model. We find supporting evidence for an mRNA-mediated mechanism of loss through comparative genome analyses that revealed a recent loss of 10 adjacent introns in a 22-exon gene in the human-pathogenic fungus Cryptococcus neoformans. We surveyed the gene structures of the entire genomes of Cryptococcus gattii, which diverged from the C. neoformans lineage 37 million years ago (Mya), and C. neoformans var. grubii and var. neoformans, which diverged 18 Mya. Our comparison revealed greater than 99.9% intron conservation, with evidence from 20 genes showing evidence of intron loss, but no convincing evidence of intron gain. Our findings confirm that Cryptococcus introns have been quite stable over recent evolutionary time, with occasional mRNA-mediated intron loss events.

Biological Evolution↗

PubMed2006

Comparative genomics in C. elegans, C. briggsae, and other Caenorhabditis species.

The genome of the nematode Caenorhabditis elegans was the first animal genome sequenced. Subsequent sequencing of the Caenorhabditis briggsae genome enabled a comparison of the genomes of two nematode species. In this chapter, we describe the methods that we used to compare the C. elegans genome to that of C. briggsae. We discuss how these methods could be developed to compare the C. elegans and C. briggsae genomes to those of Caenorhabditis remanei, C. n. sp. represented by strains PB2801 and CB5161, among others (1), and Caenorhabditis japonica, which are currently being sequenced.

Animals↗

PubMed2005

Same-sex mating and the origin of the Vancouver Island Cryptococcus gattii outbreak.

Genealogy can illuminate the evolutionary path of important human pathogens. In some microbes, strict clonal reproduction predominates, as with the worldwide dissemination of Mycobacterium leprae, the cause of leprosy. In other pathogens, sexual reproduction yields clones with novel attributes, for example, enabling the efficient, oral transmission of the parasite Toxoplasma gondii. However, the roles of clonal or sexual propagation in the origins of many other microbial pathogen outbreaks remain unknown, like the recent fungal meningoencephalitis outbreak on Vancouver Island, Canada, caused by Cryptococcus gattii. Here we show that the C. gattii outbreak isolates comprise two distinct genotypes. The majority of isolates are hypervirulent and have an identical genotype that is unique to the Pacific Northwest. A minority of the isolates are significantly less virulent and share an identical genotype with fertile isolates from an Australian recombining population. Genotypic analysis reveals evidence of sexual reproduction, in which the majority genotype is the predicted offspring. However, instead of the classic a-alpha sexual cycle, the majority outbreak clone appears to have descended from two alpha mating-type parents. Analysis of nuclear content revealed a diploid environmental isolate homozygous for the major genotype, an intermediate produced during same-sex mating. These studies demonstrate how cryptic same-sex reproduction can enable expansion of a human pathogen to a new geographical niche and contribute to the ongoing production of infectious spores. This has implications for the emergence of other microbial pathogens and inbreeding in host range expansion in the fungal and other kingdoms.

Animals↗

PubMed2005

Likelihoods from summary statistics: recent divergence between species.

We describe an importance-sampling method for approximating likelihoods of population parameters based on multiple summary statistics. In this first application, we address the demographic history of closely related members of the Drosophila pseudoobscura group. We base the maximum-likelihood estimation of the time since speciation and the effective population sizes of the extant and ancestral populations on the pattern of nucleotide variation at DPS2002, a noncoding region tightly linked to a paracentric inversion that strongly contributes to reproductive isolation. Consideration of summary statistics rather than entire nucleotide sequences permits a compact description of the genealogy of the sample. We use importance sampling first to propose a genealogical and mutational history consistent with the observed array of summary statistics and then to correct the likelihood with the exact probability of the history determined from a system of recursions. Analysis of a subset of the data, for which recursive computation of the exact likelihood was feasible, indicated close agreement between the approximate and exact likelihoods. Our results for the complete data set also compare well with those obtained through Metropolis-Hastings sampling of fully resolved genealogies of entire nucleotide sequences.

Animals↗

PubMed2005

Investigating hookworm genomes by comparative analysis of two Ancylostoma species.

BACKGROUND: Hookworms, infecting over one billion people, are the mostly closely related major human parasites to the model nematode Caenorhabditis elegans. Applying genomics techniques to these species, we analyzed 3,840 and 3,149 genes from Ancylostoma caninum and A. ceylanicum. RESULTS: Transcripts originated from libraries representing infective L3 larva, stimulated L3, arrested L3, and adults. Most genes are represented in single stages including abundant transcripts like hsp-20 in infective L3 and vit-3 in adults. Over 80% of the genes have homologs in C. elegans, and nearly 30% of these were with observable RNA interference phenotypes. Homologies were identified to nematode-specific and clade V specific gene families. To study the evolution of hookworm genes, 574 A. caninum/A. ceylanicum orthologs were identified, all of which were found to be under purifying selection with distribution ratios of nonsynonymous to synonymous amino acid substitutions similar to that reported for C. elegans/C. briggsae orthologs. The phylogenetic distance between A. caninum and A. ceylanicum is almost identical to that for C. elegans/C. briggsae. CONCLUSION: The genes discovered should substantially accelerate research toward better understanding of the parasites' basic biology as well as new therapies including vaccines and novel anthelmintics.

Ancylostoma↗

PubMed2005

Estimating the tempo and mode of gene family evolution from comparative genomic data.

Comparison of whole genomes has revealed that changes in the size of gene families among organisms is quite common. However, there are as yet no models of gene family evolution that make it possible to estimate ancestral states or to infer upon which lineages gene families have contracted or expanded. In addition, large differences in family size have generally been attributed to the effects of natural selection, without a strong statistical basis for these conclusions. Here we use a model of stochastic birth and death for gene family evolution and show that it can be efficiently applied to multispecies genome comparisons. This model takes into account the lengths of branches on phylogenetic trees, as well as duplication and deletion rates, and hence provides expectations for divergence in gene family size among lineages. The model offers both the opportunity to identify large-scale patterns in genome evolution and the ability to make stronger inferences regarding the role of natural selection in gene family expansion or contraction. We apply our method to data from the genomes of five yeast species to show its applicability.

Evolution, Molecular↗

PubMed2004

Disentangling the effects of demography and selection in human history.

Demographic events affect all genes in a genome, whereas natural selection has only local effects. Using publicly available data from 151 loci sequenced in both European-American and African-American populations, we attempt to distinguish the effects of demography and selection. To analyze large sets of population genetic data such as this one, we introduce "Perlymorphism," a Unix-based suite of analysis tools. Our analyses show that the demographic histories of human populations can account for a large proportion of effects on the level and frequency of variation across the genome. The African-American population shows both a higher level of nucleotide diversity and more negative values of Tajima's D statistic than does a European-American population. Using coalescent simulations, we show that the significantly negative values of the D statistic in African-Americans and the positive values in European-Americans are well explained by relatively simple models of population admixture and bottleneck, respectively. Working within these nonequilibrium frameworks, we are still able to show deviations from neutral expectations at a number of loci, including ABO and TRPV6. In addition, we show that the frequency spectrum of mutations--corrected for levels of polymorphism--is correlated with recombination rate only in European-Americans. These results are consistent with repeated selective sweeps in non-African populations, in agreement with recent reports using microsatellite data.

ABO Blood-Group System↗