Bayesian analysis of gene expression levels: statistical quantification of relative mRNA level across multiple strains or treatments

Light sensing by opsins and fungal ecology: NOP-1 modulates entry into sexual reproduction in response to environmental cues

Understanding the genetic basis of the switch from asexual to sexual lifestyles in response to sometimes rapid environmental changes is one of the major challenges in fungal ecology. Light appears to play a critical role in the asexual-sexual switch-but fungal genomes harbour diverse light sensors. Fungal opsins are homologous to bacterial green-light-sensory rhodopsins, and their organismal functions in fungi have not been well understood. Three of these opsin-like proteins were widely distributed across fungal genomes, but homologs of the Fusarium opsin-like protein CarO were present only in plant-associated fungi. Key amino acids, including potential retinal binding sites, functionally diverged on the phylogeny of opsins. This diversification of opsin-like proteins could be correlated with life history-associated differences among fungi in their expression and function during morphological development. In Neurospora crassa and related species, knockout of the opsin NOP-1 led to a phenotype in the regulation of the asexual-sexual switch, modulating response to both light and oxygen conditions. Sexual development commenced early in ∆nop-1 strains cultured in unsealed plates under constant blue and white light. Furthermore, comparative transcriptomics showed that the expression of nop-1 is light-dependent and that the ∆nop-1 strain abundantly expresses genes involved in oxidative stress response, genes enriched in NAD/NADP binding sites, genes with functions in proton transmembrane movement and catalase activity, and genes involved in the homeostasis of protons. Based on these observations, we contend that light and oxidative stress regulate the switch via light-responsive and ROS pathways in model fungus N. crassa and other fungi.

Reproductive activity triggers accelerated male mortality and decreases lifespan: genetic and gene expression determinants in Drosophila

Reproduction and aging evolved to be intimately associated. Experimental selection for early-life reproduction drives the evolution of decreased longevity in Drosophila whereas experimental selection for increased longevity leads to changes in reproduction. Although life history theory offers hypotheses to explain these relationships, the genetic architecture and molecular mechanisms underlying reproduction-longevity associations remain a matter of debate. Here we show that mating triggers accelerated mortality in males and identify hundreds of genes that are modulated upon mating in the fruit fly Drosophila melanogaster. Interrogation of genome-wide gene expression in virgin and recently mated males revealed coherent responses, with biological processes that are upregulated (testis-specific gene expression) or downregulated (metabolism and mitochondria-related functions) upon mating. Furthermore, using a panel of genotypes from the Drosophila Synthetic Population Resource (DSPR) as a source of naturally occurring genetic perturbation, we uncover abundant variation in longevity and reproduction-induced mortality among genotypes. Genotypes displayed more than fourfold variation in longevity and reproduction-induced mortality that can be traced to variation in specific segments of the genome. The data reveal individual variation in sensitivity to reproduction and physiological processes that are enhanced and suppressed upon mating. These results raise the prospect that variation in longevity and age-related traits could be traced to processes that coordinate germline and somatic function.

The Fast-Evolving phy-2 Gene Modulates Sexual Development in Response to Light in the Model Fungus Neurospora crassa

Rapid responses to changes in incident light are critical to the guidance of behavior and development in most species. Phytochrome light receptors in particular play key roles in bacterial physiology and plant development, but their functions and regulation are less well understood in fungi. Nevertheless, genome-wide expression measurements provide key information that can guide experiments that reveal how genes respond to environmental signals and clarify their role in development. We performed functional genomic and phenotypic analyses of the two phytochromes in Neurospora crassa, a fungal model adapted to a postfire environment that experiences dramatically variable light conditions. Expression of phy-1 and phy-2 was low in early sexual development and in the case of phy-2 increased in late sexual development. Under light stimulation, strains with the phytochromes deleted exhibited increased expression of sexual development-related genes. Moreover, under red light, the phy-2 knockout strain commenced sexual development early. In the evolution of phytochromes within ascomycetes, at least two duplications have occurred, and the faster-evolving phy-2 gene has frequently been lost. Additionally, the three key cysteine sites that are critical for bacterial and plant phytochrome function are not conserved within fungal phy-2 homologs. Through the action of phytochromes, transitions between asexual and sexual reproduction are modulated by light level and light quality, presumably as an adaptation for fast asexual growth and initiation of sexual reproduction of N. crassa in exposed postfire ecosystems.

Environmental signals, including light, play critical roles in regulating fungal growth and pathogenicity, and balance of asexual and sexual reproduction is critical in fungal pathogens’ incidence, virulence, and distribution. Red light sensing by phytochromes is well known to play critical roles in bacterial physiology and plant development. Homologs of phytochromes were first discovered in the fungal model Neurospora crassa and then subsequently in diverse other fungi, including many plant pathogens. Our study investigated the evolution of red light sensors in ascomycetes and confirmed—using the model fungus Neurospora crassa—their roles in modulating the asexual-sexual reproduction balance in fungi. Our findings also provide a key insight into one of the most poorly understood aspects of fungal biology, suggesting that further study of the function of phytochromes in fungi is critical to reveal the genetic basis of the asexual-sexual switch responsible for fungal growth and distribution, including diverse and destructive plant pathogens.

Genomic and Gene-Expression Comparisons among Phage-Resistant Type-IV Pilus Mutants of Pseudomonas syringae pathovar phaseolicola

Pseudomonas syringae pv. phaseolicola (Pph) is a significant bacterial pathogen of agricultural crops, and phage Φ6 and other members of the dsRNA virus family Cystoviridae undergo lytic (virulent) infection of Pph, using the type IV pilus as the initial site of cellular attachment. Despite the popularity of Pph/phage Φ6 as a model system in evolutionary biology, Pph resistance to phage Φ6 remains poorly characterized. To investigate differences between phage Φ6 resistant Pph strains, we examined genomic and gene expression variation among three bacterial genotypes that differ in the number of type IV pili expressed per cell: ordinary (wild-type), non-piliated, and super-piliated. Genome sequencing of non-piliated and super-piliated Pph identified few mutations that separate these genotypes from wild type Pph--and none present in genes known to be directly involved in type IV pilus expression. Expression analysis revealed that 81.1% of gene ontology (GO) terms up-regulated in the non-piliated strain were down-regulated in the super-piliated strain. This differential expression is particularly prevalent in genes associated with respiration--specifically genes in the tricarboxylic acid cycle (TCA) cycle, aerobic respiration, and acetyl-CoA metabolism. The expression patterns of the TCA pathway appear to be generally up and down-regulated, in non-piliated and super-piliated Pph respectively. As pilus retraction is mediated by an ATP motor, loss of retraction ability might lead to a lower energy draw on the bacterial cell, leading to a different energy balance than wild type. The lower metabolic rate of the super-piliated strain is potentially a result of its loss of ability to retract.

Gene Expression Evolves under a House-of-Cards Model of Stabilizing Selection

Divergence in gene regulation is hypothesized to underlie much of phenotypic evolution, but the role of natural selection in shaping the molecular phenotype of gene expression continues to be debated. To resolve the mode of gene expression, evolution requires accessible theoretical predictions for the effect of selection over long timescales. Evolutionary quantitative genetic models of phenotypic evolution can provide such predictions, yet those predictions depend on the underlying hypotheses about the distributions of mutational and selective effects that are notoriously difficult to disentangle. Here, we draw on diverse genomic data sets including expression profiles of natural genetic variation and mutation accumulation lines, empirical estimates of genomic mutation rates, and inferences of genetic architecture to differentiate contrasting hypotheses for the roles of stabilizing selection and mutation in shaping natural expression variation. Our analysis suggests that gene expression evolves in a domain of phenotype space well fit by the House-of-Cards (HC) model. Although the strength of selection inferred is sensitive to the number of loci controlling gene expression, the model is not. The consistency of these results across evolutionary time from budding yeast through fruit fly implies that this model is general and that mutational effects on gene expression are relatively large. Empirical estimates of the genetic architecture of gene expression traits imply that selection provides modest constraints on gene expression levels for most genes, but that the potential for regulatory evolution is high. Our prediction using data from laboratory environments should encourage the collection of additional data sets allowing for more nuanced parameterizations of HC models for gene expression.

Cautions about the reliability of pairwise gene correlations based on expression data

BACKGROUND

Rapid growth in the availability of genome-wide transcript abundance levels through gene expression microarrays and RNAseq promises to provide deep biological insights into the complex, genome-wide transcriptional behavior of single-celled organisms. However, this promise has not yet been fully realized.

RESULTS

We find that computation of pairwise gene associations (correlation; mutual information) across a set of 2782 total genome-wide expression samples from six diverse bacteria produces unexpectedly large variation in estimates of pairwise gene association-regardless of the metric used, the organism under study, or the number and source of the samples. We pinpoint the cause to sampling bias. In particular, in repositories of expression data (e.g., Gene Expression Omnibus, GEO), many individual genes show small differences in absolute gene expression levels across the set of samples. We demonstrate that these small differences are due mainly to "noise" instead of "signal" attributable to environmental or genetic perturbations. We show that downstream analysis using gene expression levels of genes with small differences yields biased estimates of pairwise association.

CONCLUSIONS

We propose flagging genes with small differences in absolute, RMA-normalized, expression levels (e.g., standard deviation less than 0.5), as potentially yielding biased pairwise association metrics. This strategy has the potential to substantially improve the confidence in genome-wide conclusions about transcriptional behavior in bacterial organisms. Further work is needed to further refine strategies to identify genes with small difference in expression levels prior to computing gene-gene association metrics.

Variation in the X:Autosome Distribution of Male-Biased Genes among Drosophila melanogaster Tissues and Its Relationship with Dosage Compensation

Genes that are expressed differently between males and females (sex-biased genes) often show a nonrandom distribution in their genomic location, particularly with respect to the autosomes and the X chromosome. Previous studies of Drosophila melanogaster found a general paucity of male-biased genes on the X chromosome, although this is mainly limited to comparisons of whole flies or body segments containing the reproductive organs. To better understand the chromosomal distribution of sex-biased genes in various tissues, we used a common analysis framework to analyze microarray and RNA sequence data comparing male and female gene expression in individual tissues (brain, Malpighian tubule, and gonads), composite structures (head and gonadectomized carcass), and whole flies. Although there are relatively few sex-biased genes in the brain, there is a strong and highly significant enrichment of male-biased genes on the X chromosome. A weaker enrichment of X-linked male-biased genes is seen in the head, suggesting that most of this signal comes from the brain. In all other tissues, there is either no departure from the random expectation or a significant paucity of male-biased genes on the X chromosome. The brain and head also differ from other tissues in that their male-biased genes are significantly closer to binding sites of the dosage compensation complex. We propose that the interplay of dosage compensation and sex-specific regulation can explain the observed differences between tissues and reconcile disparate results reported in previous studies.

Social odors conveying dominance and reproductive information induce rapid physiological and neuromolecular changes in a cichlid fish

Background

Social plasticity is a pervasive feature of animal behavior. Animals adjust the expression of their social behavior to the daily changes in social life and to transitions between life-history stages, and this ability has an impact in their Darwinian fitness. This behavioral plasticity may be achieved either by rewiring or by biochemically switching nodes of the neural network underlying social behavior in response to perceived social information. Independent of the proximate mechanisms, at the neuromolecular level social plasticity relies on the regulation of gene expression, such that different neurogenomic states emerge in response to different social stimuli and the switches between states are orchestrated by signaling pathways that interface the social environment and the genotype. Here, we test this hypothesis by characterizing the changes in the brain profile of gene expression in response to social odors in the Mozambique Tilapia, Oreochromis mossambicus. This species has a rich repertoire of social behaviors during which both visual and chemical information are conveyed to conspecifics. Specifically, dominant males increase their urination frequency during agonist encounters and during courtship to convey chemical information reflecting their dominance status.

Results

We recorded electro-olfactograms to test the extent to which the olfactory epithelium can discriminate between olfactory information from dominant and subordinate males as well as from pre- and post-spawning females. We then performed a genome-scale gene expression analysis of the olfactory bulb and the olfactory cortex homolog in order to identify the neuromolecular systems involved in processing these social stimuli.

Conclusions

Our results show that different olfactory stimuli from conspecifics’ have a major impact in the brain transcriptome, with different chemical social cues eliciting specific patterns of gene expression in the brain. These results confirm the role of rapid changes in gene expression in the brain as a genomic mechanism underlying behavioral plasticity and reinforce the idea of an extensive transcriptional plasticity of cichlid genomes, especially in response to rapid changes in their social environment.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1255-4) contains supplementary material, which is available to authorized users.

Identification of differentially-expressed genes in intestinal gastric cancer by microarray analysis

Gastric cancer (GC) is one of the most frequent malignant tumors. In order to systematically characterize the cellular and molecular mechanisms of intestinal GC development, in this study, we used 22K oligonucleotide microarrays and bioinformatics analysis to evaluate the gene expression profiles of GC in 45 tissue samples, including 20 intestinal GC tissue samples, 20 normal appearing tissues (NATs) adjacent to tumors and 5 noncancerous gastric mucosa tissue samples. These profiles allowed us to explore the transcriptional characteristics of GC and determine the change patterns in gene expression that may be of clinical significance. 1519 and 1255 differentially-expressed genes (DEGs) were identified in intestinal GC tissues and NATs, respectively, as determined by Bayesian analysis (P<0.001). These genes were associated with diverse functions such as mucosa secretion, metabolism, proliferation, signaling and development, which occur at different stages of GC development.

Interaction between bisphenol A and dietary sugar affects global gene transcription in Drosophila melanogaster

Human exposure to environmental toxins is a public health issue. The microarray data available in the Gene Expression Omnibus database under accession number GSE55655 and GSE55670GSE55655GSE55670 show the isolated and combined effects of dietary sugar and two organic compounds present in a variety of plastics [bisphenol A (BPA) and Bis(2-ethylhexyl) phthalate (DEHP)] on global gene expression in Drosophila melanogaster. The study was carried out with samples collected from flies exposed to these compounds for a limited period of time (48 h) in the adult stage, or throughout the entire development of the insect. The arrays were normalized using the limma/Bioconductor package. Differential expression was inferred using linear models in limma and BAGEL. The data show that each compound had its unique consequences to gene expression, and that the individual effect of each organic compound is maximized with the joint ingestion of dietary sugar.

Dynamics of transcriptome evolution in the model eukaryote Neurospora

Mounting evidence indicates that changes in the transcriptome contribute significantly to the phenotypic differentiation of closely related species. Nonetheless, further genome-wide studies, spanning a broad range of organisms, are needed to decipher the factors driving transcriptome evolution. The model Neurospora (Ascomycota) comprises a simple system for empirically studying the evolutionary dynamics of the transcriptome. Here, we studied the evolution of gene expression in Neurospora crassa and Neurospora tetrasperma and show that patterns of transcriptome evolution are connected to genome evolution, tissue type and sexual identity (mating types, mat A and mat a) in these eukaryotes. Based on the comparisons of inter- and intraspecies expression divergence, our data reveal that rapid expression divergence is more apt to occur in sexual/female (SF) than vegetative/male (VM) tissues. In addition, interspecies gene expression and protein sequence divergence were strongly correlated for SF, but not VM, tissue. A correlation between transcriptome and protein evolution parallels findings from certain animals, but not yeast, and add support for the theory that expression evolution differs fundamentally among multicellular and unicellular eukaryotes. Finally, we found that sexual identity in these hermaphroditic Neurospora species is connected to interspecies expression divergence in a tissue-dependent manner: rapid divergence occurred for mat A- and mat a-biased genes from SF and VM tissues, respectively. Based on these findings, it is hypothesized that rapid interspecies transcriptome evolution is shifting the mating types of Neurospora towards distinct female and male phenotypes, that is, sexual dimorphism.

Epigenetic inheritance and genome regulation: is DNA methylation linked to ploidy in haplodiploid insects?

Organisms show great variation in ploidy level. For example, chromosome copy number varies among cells, individuals and species. One particularly widespread example of ploidy variation is found in haplodiploid taxa, wherein males are typically haploid and females are typically diploid. Despite the prevalence of haplodiploidy, the regulatory consequences of having separate haploid and diploid genomes are poorly understood. In particular, it remains unknown whether epigenetic mechanisms contribute to regulatory compensation for genome dosage. To gain greater insights into the importance of epigenetic information to ploidy compensation, we examined DNA methylation differences among diploid queen, diploid worker, haploid male and diploid male Solenopsis invicta fire ants. Surprisingly, we found that morphologically dissimilar diploid males, queens and workers were more similar to one another in terms of DNA methylation than were morphologically similar haploid and diploid males. Moreover, methylation level was positively associated with gene expression for genes that were differentially methylated in haploid and diploid castes. These data demonstrate that intragenic DNA methylation levels differ among individuals of distinct ploidy and are positively associated with levels of gene expression. Thus, these results suggest that epigenetic information may be linked to ploidy compensation in haplodiploid insects. Overall, this study suggests that epigenetic mechanisms may be important to maintaining appropriate patterns of gene regulation in biological systems that differ in genome copy number.

High intake of dietary sugar enhances bisphenol A (BPA) disruption and reveals ribosome-mediated pathways of toxicity

Bisphenol A (BPA) is an organic compound to which human populations are ubiquitously exposed. Epidemiological data suggest BPA exposure might be associated with higher rates of diabetes and reproductive anomalies. Health concerns also include transgenerational consequences, but these mechanisms are crudely defined. Similarly, little is known about synergistic interactions between BPA and other substances. Here we show that acute and chronic exposure to BPA causes genome-wide modulation of several functionally coherent genetic pathways in the fruit fly Drosophila melanogaster. In particular, BPA exposure causes massive downregulation of testis-specific genes and upregulation of ribosome-associated genes widely expressed across tissues. In addition, it causes the modulation of transposable elements that are specific to the ribosomal DNA loci, suggesting that nucleolar stress might contribute to BPA toxicity. The upregulation of ribosome-associated genes and the impairment of testis-specific gene expression are significantly enhanced upon BPA exposure with a high-sugar diet. Our results suggest that BPA and dietary sugar might functionally interact, with consequences to regulatory programs in both reproductive and somatic tissues.

The TAF9 C-terminal conserved region domain is required for SAGA and TFIID promoter occupancy to promote transcriptional activation

A common function of the TFIID and SAGA complexes, which are recruited by transcriptional activators, is to deliver TBP to promoters to stimulate transcription. Neither the relative contributions of the five shared TBP-associated factor (TAF) subunits in TFIID and SAGA nor the requirement for different domains in shared TAFs for transcriptional activation is well understood. In this study, we uncovered the essential requirement for the highly conserved C-terminal region (CRD) of Taf9, a shared TAF, for transcriptional activation in yeast. Transcriptome profiling performed under Gcn4-activating conditions showed that the Taf9 CRD is required for induced expression of ∼9% of the yeast genome. The CRD was not essential for the Taf9-Taf6 interaction, TFIID or SAGA integrity, or Gcn4 interaction with SAGA in cell extracts. Microarray profiling of a SAGA mutant (spt20Δ) yielded a common set of genes induced by Spt20 and the Taf9 CRD. Chromatin immunoprecipitation (ChIP) assays showed that, although the Taf9 CRD mutation did not impair Gcn4 occupancy, the occupancies of TFIID, SAGA, and the preinitiation complex were severely impaired at several promoters. These results suggest a crucial role for the Taf9 CRD in genome-wide transcription and highlight the importance of conserved domains, other than histone fold domains, as a common determinant for TFIID and SAGA functions.

Genome-wide gene expression effects of sex chromosome imprinting in Drosophila

Imprinting is well-documented in both plant and animal species. In Drosophila, the Y chromosome is differently modified when transmitted through the male and female germlines. Here, we report genome-wide gene expression effects resulting from reversed parent-of-origin of the X and Y chromosomes. We found that hundreds of genes are differentially expressed between adult male Drosophila melanogaster that differ in the maternal and paternal origin of the sex chromosomes. Many of the differentially regulated genes are expressed specifically in testis and midgut cells, suggesting that sex chromosome imprinting might globally impact gene expression in these tissues. In contrast, we observed much fewer Y-linked parent-of-origin effects on genome-wide gene expression in females carrying a Y chromosome, indicating that gene expression in females is less sensitive to sex chromosome parent-of-origin. Genes whose expression differs between females inheriting a maternal or paternal Y chromosome also show sex chromosome parent-of-origin effects in males, but the direction of the effects on gene expression (overexpression or underexpression) differ between the sexes. We suggest that passage of sex chromosome chromatin through male meiosis may be required for wild-type function in F1 progeny, whereas disruption of Y-chromosome function through passage in the female germline likely arises because the chromosome is not adapted to the female germline environment.

Evolutionary dynamics of sex-biased genes in a hermaphrodite fungus

Differential gene expression is believed to largely explain sexually dichotomous phenotypes. This phenomenon is especially significant in hermaphrodites, in which male and female sexual tissues have identical genotypes. Sex differences in transcription have been linked to molecular evolution: genes with higher expression in male compared with female sexual tissues (i.e., male-biased genes) have been associated with rapid gene divergence in various animals and plants, implying that selective differences exist among the sexual structures. In the present investigation, we examined expressed sequence tags, microarrays, and gene sequence data from the hermaphroditic fungus Neurospora crassa and confirmed selective differences of genes with disparate expression among male versus female sexual structures in this organism. The results held across various genotypes and stages of sexual development. Furthermore, our data showed that N. crassa comprises a rare example of an organism where female-biased genes evolve rapidly; they exhibited faster evolution at the protein level and reduced optimal codon usage compared with male-biased genes, sexually unbiased genes, and vegetative genes. Female-biased genes also had a greater portion of sites that experienced positive selection and showed stronger signals of selective sweeps than male-biased genes, suggesting that the rapid evolution is at least partly driven by adaptive evolution. Distinctive aspects of the reproductive biology of N. crassa which might explain the rapid evolution of female-biased genes are discussed, particularly the propensity for female-female competition during mating, as well as the multifunctional nature of male structures. The present findings open new opportunities to test hypotheses about sex-biased gene expression and molecular evolution.

Passage through the mammalian gut triggers a phenotypic switch that promotes Candida albicans commensalism

Among ~5,000,000 fungal species,¹Candida albicans is exceptional in its lifelong association with humans, either within the gastrointestinal microbiome or as an invasive pathogen.² Opportunistic infections are generally ascribed to defective host immunity ³ but may require specific microbial programs. Here, we report that exposure of C. albicans to the mammalian gut triggers a developmental switch, driven by the Wor1 transcription factor, to a commensal cell type. Wor1 expression was previously observed only in rare genetic backgrounds,^4–6 where it controls a white-opaque switch for mating.^4–7 We show that passage of wild-type cells through the murine gastrointestinal tract triggers WOR1 expression and a novel phenotypic switch. The resulting GUT (Gastrointestinally-IndUced Transition) cells differ morphologically and functionally from previously defined cell types, including opaque, and express a transcriptome that is optimized for the digestive tract. The white-GUT switch illuminates how a microorganism utilizes distinct genetic programs to transition between commensalism and invasive pathogenesis.

Sex-linked transcriptional divergence in the hermaphrodite fungus Neurospora tetrasperma

In the filamentous ascomycete Neurospora tetrasperma, a large (approx. 7 Mbp) region of suppressed recombination surrounds the mating-type (mat) locus. While the remainder of the genome is largely homoallelic, this region of recombinational suppression, extending over 1500 genes, is associated with sequence divergence. Here, we used microarrays to examine how the molecular phenotype of gene expression level is linked to this divergent region, and thus to the mating type. Culturing N. tetrasperma on agar media that induce sexual/female or vegetative/male tissue, we found 196 genes significantly differentially expressed between mat A and mat a mating types. Our data show that the genes exhibiting mat-linked expression are enriched in the region genetically linked to mating type, and sequence and expression divergence are positively correlated. Our results indicate that the phenotype of mat A strains is optimized for traits promoting sexual/female development and the phenotype of mat a strains for vegetative/male development. This discovery of differentially expressed genes associated with mating type provides a link between genotypic and phenotypic divergence in this taxon and illustrates a fungal analogue to sexual dimorphism found among animals and plants.

Chromatin-associated proteins HP1 and Mod(mdg4) modify Y-linked regulatory variation in the drosophila testis

Chromatin remodeling is crucial for gene regulation. Remodeling is often mediated through chemical modifications of the DNA template, DNA-associated proteins, and RNA-mediated processes. Y-linked regulatory variation (YRV) refers to the quantitative effects that polymorphic tracts of Y-linked chromatin exert on gene expression of X-linked and autosomal genes. Here we show that naturally occurring polymorphisms in the Drosophila melanogaster Y chromosome contribute disproportionally to gene expression variation in the testis. The variation is dependent on wild-type expression levels of mod(mdg4) as well as Su(var)205; the latter gene codes for heterochromatin protein 1 (HP1) in Drosophila. Testis-specific YRV is abolished in genotypes with heterozygous loss-of-function mutations for mod(mdg4) and Su(var)205 but not in similar experiments with JIL-1. Furthermore, the Y chromosome differentially regulates several ubiquitously expressed genes. The results highlight the requirement for wild-type dosage of Su(var)205 and mod(mdg4) in enabling naturally occurring Y-linked regulatory variation in the testis. The phenotypes that emerge in the context of wild-type levels of the HP1 and Mod(mdg4) proteins might be part of an adaptive response to the environment.

Natural variation of the Y chromosome suppresses sex ratio distortion and modulates testis-specific gene expression in Drosophila simulans

X-linked sex-ratio distorters that disrupt spermatogenesis can cause a deficiency in functional Y-bearing sperm and a female-biased sex ratio. Y-linked modifiers that restore a normal sex ratio might be abundant and favored when a X-linked distorter is present. Here we investigated natural variation of Y-linked suppressors of sex-ratio in the Winters systems and the ability of these chromosomes to modulate gene expression in Drosophila simulans. Seventy-eight Y chromosomes of worldwide origin were assayed for their resistance to the X-linked sex-ratio distorter gene Dox. Y chromosome diversity caused males to sire ∼63% to ∼98% female progeny. Genome-wide gene expression analysis revealed hundreds of genes differentially expressed between isogenic males with sensitive (high sex ratio) and resistant (low sex ratio) Y chromosomes from the same population. Although the expression of about 75% of all testis-specific genes remained unchanged across Y chromosomes, a subset of post-meiotic genes was upregulated by resistant Y chromosomes. Conversely, a set of accessory gland-specific genes and mitochondrial genes were downregulated in males with resistant Y chromosomes. The D. simulans Y chromosome also modulated gene expression in XXY females in which the Y-linked protein-coding genes are not transcribed. The data suggest that the Y chromosome might exert its regulatory functions through epigenetic mechanisms that do not require the expression of protein-coding genes. The gene network that modulates sex ratio distortion by the Y chromosome is poorly understood, other than that it might include interactions with mitochondria and enriched for genes expressed in post-meiotic stages of spermatogenesis.

Comparative transcriptomics of infectious spores from the fungal pathogen Histoplasma capsulatum reveals a core set of transcripts that specify infectious and pathogenic states

Histoplasma capsulatum is a fungal pathogen that infects both healthy and immunocompromised hosts. In regions where it is endemic, H. capsulatum grows in the soil and causes respiratory and systemic disease when inhaled by humans. An interesting aspect of H. capsulatum biology is that it adopts specialized developmental programs in response to its environment. In the soil, it grows as filamentous chains of cells (mycelia) that produce asexual spores (conidia). When the soil is disrupted, conidia aerosolize and are inhaled by mammalian hosts. Inside a host, conidia germinate into yeast-form cells that colonize immune cells and cause disease. Despite the ability of conidia to initiate infection and disease, they have not been explored on a molecular level. We developed methods to purify H. capsulatum conidia, and we show here that these cells germinate into filaments at room temperature and into yeast-form cells at 37°C. Conidia internalized by macrophages germinate into the yeast form and proliferate within macrophages, ultimately lysing the host cells. Similarly, infection of mice with purified conidia is sufficient to establish infection and yield viable yeast-form cells in vivo. To characterize conidia on a molecular level, we performed whole-genome expression profiling of conidia, yeast, and mycelia from two highly divergent H. capsulatum strains. In parallel, we used homology and protein domain analysis to manually annotate the predicted genes of both strains. Analyses of the resultant data defined sets of transcripts that reflect the unique molecular states of H. capsulatum conidia, yeast, and mycelia.

Yeast response to LA virus indicates coadapted global gene expression during mycoviral infection

Viruses that infect fungi have a ubiquitous distribution and play an important role in structuring fungal communities. Most of these viruses have an unusual life history in that they are propagated exclusively via asexual reproduction or fission of fungal cells. This asexual mode of transmission intimately ties viral reproductive success to that of its fungal host and should select for viruses that have minimal deleterious impact on the fitness of their hosts. Accordingly, viral infections of fungi frequently do not measurably impact fungal growth, and in some instances, increase the fitness of the fungal host. Here we determine the impact of the loss of coinfection by LA virus and the virus-like particle M1 upon global gene expression of the fungal host Saccharomyces cerevisiae and provide evidence supporting the idea that coevolution has selected for viral infection minimally impacting host gene expression.

A Bayesian model for pooling gene expression studies that incorporates co-regulation information

Current Bayesian microarray models that pool multiple studies assume gene expression is independent of other genes. However, in prokaryotic organisms, genes are arranged in units that are co-regulated (called operons). Here, we introduce a new Bayesian model for pooling gene expression studies that incorporates operon information into the model. Our Bayesian model borrows information from other genes within the same operon to improve estimation of gene expression. The model produces the gene-specific posterior probability of differential expression, which is the basis for inference. We found in simulations and in biological studies that incorporating co-regulation information improves upon the independence model. We assume that each study contains two experimental conditions: a treatment and control. We note that there exist environmental conditions for which genes that are supposed to be transcribed together lose their operon structure, and that our model is best carried out for known operon structures.

Analyses of expressed sequence tags in Neurospora reveal rapid evolution of genes associated with the early stages of sexual reproduction in fungi

Background

The broadly accepted pattern of rapid evolution of reproductive genes is primarily based on studies of animal systems, although several examples of rapidly evolving genes involved in reproduction are found in diverse additional taxa. In fungi, genes involved in mate recognition have been found to evolve rapidly. However, the examples are too few to draw conclusions on a genome scale.

Results

In this study, we performed microarray hybridizations between RNA from sexual and vegetative tissues of two strains of the heterothallic (self-sterile) filamentous ascomycete Neurospora intermedia, to identify a set of sex-associated genes in this species. We aligned Expressed Sequence Tags (ESTs) from sexual and vegetative tissue of N. intermedia to orthologs from three closely related species: N. crassa, N. discreta and N. tetrasperma. The resulting four-species alignments provided a dataset for molecular evolutionary analyses. Our results confirm a general pattern of rapid evolution of fungal sex-associated genes, compared to control genes with constitutive expression or a high relative expression during vegetative growth. Among the rapidly evolving sex-associated genes, we identified candidates that could be of importance for mating or fruiting-body development. Analyses of five of these candidate genes from additional species of heterothallic Neurospora revealed that three of them evolve under positive selection.

Conclusions

Taken together, our study represents a novel finding of a genome-wide pattern of rapid evolution of sex-associated genes in the fungal kingdom, and provides a list of candidate genes important for reproductive isolation in Neurospora.

Genomic dissection and prioritizing of candidate genes of QTL for regulating spontaneous arthritis on chromosome 1 in mice deficient for interleukin-1 receptor antagonist

Rheumatoid arthritis is a heterogeneous disease with clinical and biological polymorphisms. IL-1RN is a protein that binds to interleukin-1 (IL-1) receptors and inhibits the binding of IL-1-alpha and IL-1-beta. IL-1RN levels are elevated in the blood of patients with a variety of infectious, immune, and traumatic conditions. Balb/c mice deficient in IL-1ra (mouse gene of IL-1RN) develop spontaneous autoimmune arthritis while DBA/1 mice deficient in IL-1ra do not. Previously, we identified a major QTL that regulates the susceptibility to arthritis in Balb/c mice with IL-1ra deficiency. In this study, we found that the QTL may contain two peaks that are regulated by two sets of candidate genes. By haplotype analysis, the total genomic regions of candidate genes were reduced from about 19 Mbp to approximately 9 Mbp. The total number of candidate genes was reduced from 208 to 21.

Diversification of a protein kinase cascade: IME-2 is involved in nonself recognition and programmed cell death in Neurospora crassa

Kinase cascades and the modification of proteins by phosphorylation are major mechanisms for cell signaling and communication, and evolution of these signaling pathways can contribute to new developmental or environmental response pathways. The Saccharomyces cerevisiae kinase Ime2 has been well characterized for its role in meiosis. However, recent studies have revealed alternative functions for Ime2 in both S. cerevisiae and other fungi. In the filamentous fungus Neurospora crassa, the IME2 homolog (ime-2) is not required for meiosis. Here we determine that ime-2 interacts genetically with a transcription factor vib-1 during nonself recognition and programmed cell death (PCD). Mutations in vib-1 (Δvib-1) suppress PCD due to nonself recognition events; however, a Δvib-1 Δime-2 mutant restored wild-type levels of cell death. A role for ime-2 in the post-translational processing and localization of a mitochondrial matrix protein was identified, which may implicate mitochondria in N. crassa nonself recognition and PCD. Further, Δvib-1 strains do not produce extracellular proteases, but protease secretion reverted to near wild-type levels in a Δvib-1 Δime-2 strain. Mass spectrometry analysis revealed that the VIB-1 protein is phosphorylated at several sites, including a site that matches the IME-2 consensus. The genetic and biochemical data for ime-2 and vib-1 indicate that IME-2 is a negative regulator of VIB-1 and suggest parallel negative regulation by IME-2 of a cell death pathway in N. crassa that functions in concert with the VIB-1 cell death pathway. Thus, IME2 kinase function has evolved following the divergence of S. cerevisiae and N. crassa and provides insight into the evolution of kinases and their regulatory targets.

Sex-specific gene expression during asexual development of Neurospora crassa

The impact of loci that determine sexual identity upon the asexual, dominant stage of fungal life history has been well studied. To investigate their impact, expression differences between strains of different mating type during asexual development were assayed, with RNA sampled from otherwise largely isogenic mat A and mat a strains of Neurospora crassa at early, middle, and late clonal stages of development. We observed significant differences in overall gene expression between mating types across clonal development, especially at late development stages. The expression levels of mating-type genes and pheromone genes were assayed by reverse transcription and quantitative PCR, revealing expression of pheromone and receptor genes in strains of both mating types in all development stages, and revealing that mating type (mat) genes were increasingly expressed over the course of asexual development. Interestingly, among differentially expressed genes, the mat A genotype more frequently exhibited a higher expression level than mat a, and demonstrated greater transcriptional regulatory dynamism. Significant up-regulation of expression was observed for many late light-responsive genes at late asexual development stages. Further investigation of the impact of light and the roles of light response genes in asexual development of both mating types are warranted.

Demasculinization of the Anopheles gambiae X chromosome

Background

In a number of organisms sex-biased genes are non-randomly distributed between autosomes and the shared sex chromosome X (or Z). Studies on Anopheles gambiae have produced conflicting results regarding the underrepresentation of male-biased genes on the X chromosome and it is unclear to what extent sexual antagonism, dosage compensation or X-inactivation in the male germline, the evolutionary forces that have been suggested to affect the chromosomal distribution of sex-biased genes, are operational in Anopheles.

Results

We performed a meta-analysis of sex-biased gene expression in Anopheles gambiae which provides evidence for a general underrepresentation of male-biased genes on the X-chromosome that increased in significance with the observed degree of sex-bias. A phylogenomic comparison between Drosophila melanogaster, Aedes aegypti and Culex quinquefasciatus also indicates that the Anopheles X chromosome strongly disfavours the evolutionary conservation of male-biased expression and that novel male-biased genes are more likely to arise on autosomes. Finally, we demonstrate experimentally that transgenes situated on the Anopheles gambiae X chromosome are transcriptionally silenced in the male germline.

Conclusion

The data presented here support the hypothesis that the observed demasculinization of the Anopheles X chromosome is driven by X-chromosome inactivation in the male germline and by sexual antagonism. The demasculinization appears to be the consequence of a loss of male-biased expression, rather than a failure in the establishment or the extinction of male-biased genes.

Unravelling the molecular basis for light modulated cellulase gene expression - the role of photoreceptors in Neurospora crassa

BACKGROUND

Light represents an important environmental cue, which exerts considerable influence on the metabolism of fungi. Studies with the biotechnological fungal workhorse Trichoderma reesei (Hypocrea jecorina) have revealed an interconnection between transcriptional regulation of cellulolytic enzymes and the light response. Neurospora crassa has been used as a model organism to study light and circadian rhythm biology. We therefore investigated whether light also regulates transcriptional regulation of cellulolytic enzymes in N. crassa.

RESULTS

We show that the N. crassa photoreceptor genes wc-1, wc-2 and vvd are involved in regulation of cellulase gene expression, indicating that this phenomenon is conserved among filamentous fungi. The negative effect of VVD on production of cellulolytic enzymes is thereby accomplished by its role in photoadaptation and hence its function in White collar complex (WCC) formation. In contrast, the induction of vvd expression by the WCC does not seem to be crucial in this process. Additionally, we found that WC-1 and WC-2 not only act as a complex, but also have individual functions upon growth on cellulose.

CONCLUSIONS

Genome wide transcriptome analysis of photoreceptor mutants and evaluation of results by analysis of mutant strains identified several candidate genes likely to play a role in light modulated cellulase gene expression. Genes with functions in amino acid metabolism, glycogen metabolism, energy supply and protein folding are enriched among genes with decreased expression levels in the wc-1 and wc-2 mutants. The ability to properly respond to amino acid starvation, i. e. up-regulation of the cross pathway control protein cpc-1, was found to be beneficial for cellulase gene expression. Our results further suggest a contribution of oxidative depolymerization of cellulose to plant cell wall degradation in N. crassa.

Differential impact of nutrition on developmental and metabolic gene expression during fruiting body development in Neurospora crassa

Fungal fruiting body size and form are influenced by the ecology of the species, including diverse environmental stimuli. Accordingly, nutritional resources available to the fungus during development can be vital to successful production of fruiting bodies. To investigate the effect of nutrition, perithecial development of Neurospora crassa was induced on two different media, a chemically sparsely nutritive Synthetic Crossing Medium (SCM) and a natural Carrot Agar (CA). Protoperithecia were collected before crossing, and perithecia were collected at 2, 24, 48, 72, 96, 120, and at full maturity 144 h after crossing. No differences in fruiting body morphology were observed between the two media at any time point. A circuit of microarray hybridizations comparing cDNA from all neighboring stages was performed. For a majority of differentially expressed genes, expression was higher in SCM than in CA, and expression of core metabolic genes was particularly affected. Effects of nutrition were highest in magnitude before crossing, lowering in magnitude during early perithecial development. Interestingly, metabolic effects of the media were also large in magnitude during late perithecial development, at which stage the lower expression in CA presumably reflected the continued intake of diverse complex initial compounds, diminishing the need for expression of anabolic pathways. However, for genes with key regulatory roles in sexual development, including pheromone precursor ccg-4 and poi2, expression patterns were similar between treatments. When possible, a common nutritional environment is ideal for comparing transcriptional profiles between different fungi. Nevertheless, the observed consistency of the developmental program across media, despite considerable metabolic differentiation is reassuring. This result facilitates comparative studies that will require different nutritional resources for sexual development in different fungi.

Genome features of "Dark-fly", a Drosophila line reared long-term in a dark environment

Organisms are remarkably adapted to diverse environments by specialized metabolisms, morphology, or behaviors. To address the molecular mechanisms underlying environmental adaptation, we have utilized a Drosophila melanogaster line, termed “Dark-fly”, which has been maintained in constant dark conditions for 57 years (1400 generations). We found that Dark-fly exhibited higher fecundity in dark than in light conditions, indicating that Dark-fly possesses some traits advantageous in darkness. Using next-generation sequencing technology, we determined the whole genome sequence of Dark-fly and identified approximately 220,000 single nucleotide polymorphisms (SNPs) and 4,700 insertions or deletions (InDels) in the Dark-fly genome compared to the genome of the Oregon-R-S strain, a control strain. 1.8% of SNPs were classified as non-synonymous SNPs (nsSNPs: i.e., they alter the amino acid sequence of gene products). Among them, we detected 28 nonsense mutations (i.e., they produce a stop codon in the protein sequence) in the Dark-fly genome. These included genes encoding an olfactory receptor and a light receptor. We also searched runs of homozygosity (ROH) regions as putative regions selected during the population history, and found 21 ROH regions in the Dark-fly genome. We identified 241 genes carrying nsSNPs or InDels in the ROH regions. These include a cluster of alpha-esterase genes that are involved in detoxification processes. Furthermore, analysis of structural variants in the Dark-fly genome showed the deletion of a gene related to fatty acid metabolism. Our results revealed unique features of the Dark-fly genome and provided a list of potential candidate genes involved in environmental adaptation.

Deciphering transcriptional regulatory mechanisms associated with hemicellulose degradation in Neurospora crassa

Hemicellulose, the second most abundant plant biomass fraction after cellulose, is widely viewed as a potential substrate for the production of liquid fuels and other value-added materials. Degradation of hemicellulose by filamentous fungi requires production of many different enzymes, which are induced by biopolymers or its derivatives and regulated mainly at the transcriptional level through transcription factors (TFs). Neurospora crassa, a model filamentous fungus, expresses and secretes enzymes required for plant cell wall deconstruction. To better understand genes specifically associated with degradation of hemicellulose, we applied secretome and transcriptome analysis to N. crassa grown on beechwood xylan. We identified 34 secreted proteins and 353 genes with elevated transcription on xylan. The xylanolytic phenotype of strains with deletions in genes identified from the secretome and transcriptome analysis of the wild type was assessed, revealing functions for known and unknown proteins associated with hemicellulose degradation. By evaluating phenotypes of strains containing deletions of predicted TF genes in N. crassa, we identified a TF (XLR-1; xylan degradation regulator 1) essential for hemicellulose degradation that is an ortholog to XlnR/XYR1 in Aspergillus and Trichoderma species, respectively, a major transcriptional regulator of genes encoding both cellulases and hemicellulases. Deletion of xlr-1 in N. crassa abolished growth on xylan and xylose, but growth on cellulose and cellulolytic activity were only slightly affected. To determine the regulatory mechanisms for hemicellulose degradation, we explored the transcriptional regulon of XLR-1 under xylose, xylanolytic, and cellulolytic conditions. XLR-1 regulated only some predicted hemicellulase genes in N. crassa and was required for a full induction of several cellulase genes. Hemicellulase gene expression was induced by a combination of release from carbon catabolite repression (CCR) and induction. This systematic analysis illustrates the similarities and differences in regulation of hemicellulose degradation among filamentous fungi.

Transcriptome profiling of sexual maturation and mating in the Mediterranean fruit fly, Ceratitis capitata

Sexual maturation and mating in insects are generally accompanied by major physiological and behavioural changes. Many of these changes are related to the need to locate a mate and subsequently, in the case of females, to switch from mate searching to oviposition behaviour. The prodigious reproductive capacity of the Mediterranean fruit fly, Ceratitis capitata, is one of the factors that has led to its success as an invasive pest species. To identify the molecular changes related to maturation and mating status in male and female medfly, a microarray-based gene expression approach was used to compare the head transcriptomes of sexually immature, mature virgin, and mated individuals. Attention was focused on the changes in abundance of transcripts related to reproduction, behaviour, sensory perception of chemical stimulus, and immune system processes. Broad transcriptional changes were recorded during female maturation, while post-mating transcriptional changes in females were, by contrast, modest. In male medfly, transcriptional changes were consistent both during maturation and as a consequence of mating. Of particular note was the lack of the mating-induced immune responses that have been recorded for Drosophila melanogaster, that may be due to the different reproductive strategies of these species. This study, in addition to increasing our understanding of the molecular machinery behind maturation and mating in the medfly, has identified important gene targets that might be useful in the future management of this pest.

Copy-number variation: the balance between gene dosage and expression in Drosophila melanogaster

Copy-number variants (CNVs) reshape gene structure, modulate gene expression, and contribute to significant phenotypic variation. Previous studies have revealed CNV patterns in natural populations of Drosophila melanogaster and suggested that selection and mutational bias shape genomic patterns of CNV. Although previous CNV studies focused on heterogeneous strains, here, we established a number of second-chromosome substitution lines to uncover CNV characteristics when homozygous. The percentage of genes harboring CNVs is higher than found in previous studies. More CNVs are detected in homozygous than heterozygous substitution strains, suggesting the comparative genomic hybridization arrays underestimate CNV owing to heterozygous masking. We incorporated previous gene expression data collected from some of the same substitution lines to investigate relationships between CNV gene dosage and expression. Most genes present in CNVs show no evidence of increased or diminished transcription, and the fraction of such dosage-insensitive CNVs is greater in heterozygotes. More than 70% of the dosage-sensitive CNVs are recessive with undetectable effects on transcription in heterozygotes. A deficiency of singletons in recessive dosage-sensitive CNVs supports the hypothesis that most CNVs are subject to negative selection. On the other hand, relaxed purifying selection might account for the higher number of protein–protein interactions in dosage-insensitive CNVs than in dosage-sensitive CNVs. Dosage-sensitive CNVs that are upregulated and downregulated coincide with copy-number increases and decreases. Our results help clarify the relation between CNV dosage and gene expression in the D. melanogaster genome.

Discovering transcription factor regulatory targets using gene expression and binding data

MOTIVATION

Identifying the target genes regulated by transcription factors (TFs) is the most basic step in understanding gene regulation. Recent advances in high-throughput sequencing technology, together with chromatin immunoprecipitation (ChIP), enable mapping TF binding sites genome wide, but it is not possible to infer function from binding alone. This is especially true in mammalian systems, where regulation often occurs through long-range enhancers in gene-rich neighborhoods, rather than proximal promoters, preventing straightforward assignment of a binding site to a target gene.

RESULTS

We present EMBER (Expectation Maximization of Binding and Expression pRofiles), a method that integrates high-throughput binding data (e.g. ChIP-chip or ChIP-seq) with gene expression data (e.g. DNA microarray) via an unsupervised machine learning algorithm for inferring the gene targets of sets of TF binding sites. Genes selected are those that match overrepresented expression patterns, which can be used to provide information about multiple TF regulatory modes. We apply the method to genome-wide human breast cancer data and demonstrate that EMBER confirms a role for the TFs estrogen receptor alpha, retinoic acid receptors alpha and gamma in breast cancer development, whereas the conventional approach of assigning regulatory targets based on proximity does not. Additionally, we compare several predicted target genes from EMBER to interactions inferred previously, examine combinatorial effects of TFs on gene regulation and illustrate the ability of EMBER to discover multiple modes of regulation.

AVAILABILITY

All code used for this work is available at http://dinner-group.uchicago.edu/downloads.html.

Array CGH phylogeny: how accurate are comparative genomic hybridization-based trees?

Background

Array-based Comparative Genomic Hybridization (CGH) data have been used to infer phylogenetic relationships. However, the reliability of array CGH analysis to determine evolutionary relationships has not been well established. In most CGH work, all species and strains are compared to a single reference species, whose genome was used to design the array. In the accompanying work, we critically evaluated CGH-based phylogeny using simulated competitive hybridization data. This work showed that a limited number of conditions, principally the tree topology and placement of the reference taxon in the tree, had a strong effect on the ability to recover the correct tree topology. Here, we add to our simulation study by testing the use of CGH as a phylogenetic tool with experimental CGH data from competitive hybridizations between N. crassa and other Neurospora species. In the discussion, we add to our empirical study of Neurospora by reanalyzing of data from a previous CGH phylogenetic analysis of the yeast sensu stricto complex.

Results

Array ratio data for Neurospora and related species were normalized with loess, robust spline, and linear ratio based methods, and then used to construct Neighbor-Joining and parsimony trees. These trees were compared to published phylogenetic analyses for Neurospora based on multilocus sequence analysis (MLSA). For the Neurospora dataset, the best combination of methods resulted in recovery of the MLSA tree topology less than half the time. Our reanalysis of a yeast dataset found that trees identical to established phylogeny were recovered only by pruning taxa - including the reference taxon - from the analysis.

Conclusion

Our results indicate that CGH data can be problematic for phylogenetic analysis. Success fluctuates based on the methods utilized to construct the tree and the taxa included. Selective pruning of the taxa improves the results - an impractical approach for normal phylogenetic analysis. From the more successful methods we make suggestions on the normalization and post-normalization methods that work best in estimating genetic distance between taxa.

Modulation of fungal sensitivity to staurosporine by targeting proteins identified by transcriptional profiling

An analysis of the time-dependent genetic response to the death-inducer staurosporine was performed in Neurospora crassa by transcriptional profiling. Staurosporine induced two major genes encoding an ABC transporter and a protein with similarity to regulatory subunits of potassium channels. The transcriptional response is dependent on the activity of a novel transcription factor. Deletion mutants in differentially expressed genes displayed altered sensitivity to staurosporine, underscoring significant proteins involved in the response to the drug. A null-mutant of the ABC transporter (abc3) is extremely sensitive to staurosporine, accumulates more staurosporine than the wild type strain and is defective in energy-dependent export of the drug, indicating that the ABC3 protein is the first described staurosporine transporter. It was located in the plasma membrane by immunofluorescence microscopy. The combination of inhibitors of ABC transporters or of potassium channels with staurosporine leads to an enhanced activity against N. crassa and pathogenic fungi paving the way to the development of more potent and specific antifungals. Our results highlight the general use of transcriptional profiling for the identification of novel proteins involved in cell death and their potential use as drug targets.

Identification of the CRE-1 cellulolytic regulon in Neurospora crassa

BACKGROUND

In filamentous ascomycete fungi, the utilization of alternate carbon sources is influenced by the zinc finger transcription factor CreA/CRE-1, which encodes a carbon catabolite repressor protein homologous to Mig1 from Saccharomyces cerevisiae. In Neurospora crassa, deletion of cre-1 results in increased secretion of amylase and β-galactosidase.

METHODOLOGY/PRINCIPAL FINDINGS

Here we show that a strain carrying a deletion of cre-1 has increased cellulolytic activity and increased expression of cellulolytic genes during growth on crystalline cellulose (Avicel). Constitutive expression of cre-1 complements the phenotype of a N. crassa Δcre-1 strain grown on Avicel, and also results in stronger repression of cellulolytic protein secretion and enzyme activity. We determined the CRE-1 regulon by investigating the secretome and transcriptome of a Δcre-1 strain as compared to wild type when grown on Avicel versus minimal medium. Chromatin immunoprecipitation-PCR of putative target genes showed that CRE-1 binds to only some adjacent 5'-SYGGRG-3' motifs, consistent with previous findings in other fungi, and suggests that unidentified additional regulatory factors affect CRE-1 binding to promoter regions. Characterization of 30 mutants containing deletions in genes whose expression level increased in a Δcre-1 strain under cellulolytic conditions identified novel genes that affect cellulase activity and protein secretion.

CONCLUSIONS/SIGNIFICANCE

Our data provide comprehensive information on the CRE-1 regulon in N. crassa and contribute to deciphering the global role of carbon catabolite repression in filamentous ascomycete fungi during plant cell wall deconstruction.

Relaxed selection is a precursor to the evolution of phenotypic plasticity

Phenotypic plasticity allows organisms to produce alternative phenotypes under different conditions and represents one of the most important ways by which organisms adaptively respond to the environment. However, the relationship between phenotypic plasticity and molecular evolution remains poorly understood. We addressed this issue by investigating the evolution of genes associated with phenotypically plastic castes, sexes, and developmental stages of the fire ant Solenopsis invicta. We first determined if genes associated with phenotypic plasticity in S. invicta evolved at a rapid rate, as predicted under theoretical models. We found that genes differentially expressed between S. invicta castes, sexes, and developmental stages all exhibited elevated rates of evolution compared with ubiquitously expressed genes. We next investigated the evolutionary history of genes associated with the production of castes. Surprisingly, we found that orthologs of caste-biased genes in S. invicta and the social bee Apis mellifera evolved rapidly in lineages without castes. Thus, in contrast to some theoretical predictions, our results suggest that rapid rates of molecular evolution may not arise primarily as a consequence of phenotypic plasticity. Instead, genes evolving under relaxed purifying selection may more readily adopt new forms of biased expression during the evolution of alternate phenotypes. These results suggest that relaxed selective constraint on protein-coding genes is an important and underappreciated element in the evolutionary origin of phenotypic plasticity.

Genomic investigation of year-long and multigenerational exposures of fathead minnow to the munitions compound RDX

We assessed the impacts of exposure to an environmentally representative concentration (0.83 mg/L) of the explosive cyclotrimethylenetrinitramine (RDX) on fathead minnows (Pimephales promelas) in one-year and multigenerational bioassays. In the one-year bioassay, impacts were assessed by statistical comparisons of females from breeding groups reared in control or RDX-exposure conditions. The RDX had no significant effect on gonadosomatic index or condition factor assayed at 1 d and at one, three, six, nine, and 12 months. The liver-somatic index was significantly increased versus controls only at the 12-month timepoint. RDX had no significant effect on live-prey capture rates, egg production, or fertilization. RDX caused minimal differential-transcript expression with no consistent discernable effect on gene-functional categories for either brain or liver tissues in the one-year exposure. In the multigenerational assay, the effects of acute (96 h) exposure to RDX were compared in fish reared to the F(2) generation in either control or RDX-exposure conditions. Enrichment of gene functions including neuroexcitatory glutamate metabolism, sensory signaling, and neurological development were observed comparing control-reared and RDX-reared fish. Our results indicated that exposure to RDX at a concentration representing the highest levels observed in the environment (0.83 mg/L) had limited impacts on genomic, individual, and population-level endpoints in fathead minnows in a one-year exposure. However, multigenerational exposures altered transcript expression related to neural development and function. Environ.

Conserved toxic responses across divergent phylogenetic lineages: a meta-analysis of the neurotoxic effects of RDX among multiple species using toxicogenomics

At military training sites, a variety of pollutants such as hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX), may contaminate the area originating from used munitions. Studies investigating the mechanism of toxicity of RDX have shown that it affects the central nervous system causing seizures in humans and animals. Environmental pollutants such as RDX have the potential to affect many different species, therefore it is important to establish how phylogenetically distant species may respond to these types of emerging pollutants. In this paper, we have used a transcriptional network approach to compare and contrast the neurotoxic effects of RDX among five phylogenetically disparate species: rat (Sprague-Dawley), Northern bobwhite quail (Colinus virginianus), fathead minnow (Pimephales promelas), earthworm (Eisenia fetida), and coral (Acropora formosa). Pathway enrichment analysis indicated a conservation of RDX impacts on pathways related to neuronal function in rat, Northern bobwhite quail, fathead minnows and earthworm, but not in coral. As evolutionary distance increased common responses decreased with impacts on energy and metabolism dominating effects in coral. A neurotransmission related transcriptional network based on whole rat brain responses to RDX exposure was used to identify functionally related modules of genes, components of which were conserved across species depending upon evolutionary distance. Overall, the meta-analysis using genomic data of the effects of RDX on several species suggested a common and conserved mode of action of the chemical throughout phylogenetically remote organisms.

Ribosomal DNA deletions modulate genome-wide gene expression: "rDNA-sensitive" genes and natural variation

The ribosomal rDNA gene array is an epigenetically-regulated repeated gene locus. While rDNA copy number varies widely between and within species, the functional consequences of subtle copy number polymorphisms have been largely unknown. Deletions in the Drosophila Y-linked rDNA modifies heterochromatin-induced position effect variegation (PEV), but it has been unknown if the euchromatic component of the genome is affected by rDNA copy number. Polymorphisms of naturally occurring Y chromosomes affect both euchromatin and heterochromatin, although the elements responsible for these effects are unknown. Here we show that copy number of the Y-linked rDNA array is a source of genome-wide variation in gene expression. Induced deletions in the rDNA affect the expression of hundreds to thousands of euchromatic genes throughout the genome of males and females. Although the affected genes are not physically clustered, we observed functional enrichments for genes whose protein products are located in the mitochondria and are involved in electron transport. The affected genes significantly overlap with genes affected by natural polymorphisms on Y chromosomes, suggesting that polymorphic rDNA copy number is an important determinant of gene expression diversity in natural populations. Altogether, our results indicate that subtle changes to rDNA copy number between individuals may contribute to biologically relevant phenotypic variation.

Exploring the bZIP transcription factor regulatory network in Neurospora crassa

Transcription factors (TFs) are key nodes of regulatory networks in eukaryotic organisms, including filamentous fungi such as Neurospora crassa. The 178 predicted DNA-binding TFs in N. crassa are distributed primarily among six gene families, which represent an ancient expansion in filamentous ascomycete genomes; 98 TF genes show detectable expression levels during vegetative growth of N. crassa, including 35 that show a significant difference in expression level between hyphae at the periphery versus hyphae in the interior of a colony. Regulatory networks within a species genome include paralogous TFs and their respective target genes (TF regulon). To investigate TF network evolution in N. crassa, we focused on the basic leucine zipper (bZIP) TF family, which contains nine members. We performed baseline transcriptional profiling during vegetative growth of the wild-type and seven isogenic, viable bZIP deletion mutants. We further characterized the regulatory network of one member of the bZIP family, NCU03905. NCU03905 encodes an Ap1-like protein (NcAp-1), which is involved in resistance to multiple stress responses, including oxidative and heavy metal stress. Relocalization of NcAp-1 from the cytoplasm to the nucleus was associated with exposure to stress. A comparison of the NcAp-1 regulon with Ap1-like regulons in Saccharomyces cerevisiae, Schizosaccharomyces pombe, Candida albicans and Aspergillus fumigatus showed both conservation and divergence. These data indicate how N. crassa responds to stress and provide information on pathway evolution.

CSP and takeout genes modulate the switch between attraction and repulsion during behavioral phase change in the migratory locust

Behavioral plasticity is the most striking trait in locust phase transition. However, the genetic basis for behavioral plasticity in locusts is largely unknown. To unravel the molecular mechanisms underlying the behavioral phase change in the migratory locust Locusta migratoria, the gene expression patterns over the time courses of solitarization and gregarization were compared by oligonucleotide microarray analysis. Data analysis revealed that several gene categories relevant to peripheral olfactory perception are strongly regulated in a total of 1,444 differentially expressed genes during both time courses. Among these candidate genes, several CSP (chemosensory protein) genes and one takeout gene, LmigTO1, showed higher expression in gregarious and solitarious locusts, respectively, and displayed opposite expression trends during solitarization and gregarization. qRT-PCR experiments revealed that most CSP members and LmigTO1 exhibited antenna-rich expressions. RNA interference combined with olfactory behavioral experiments confirmed that the CSP gene family and one takeout gene, LmigTO1, are involved in the shift from repulsion to attraction between individuals during gregarization and in the reverse transition during solitarization. These findings suggest that the response to locust-emitted olfactory cues regulated by CSP and takeout genes is involved in the behavioral phase change in the migratory locust and provide a previously undescribed molecular mechanism linked to the formation of locust aggregations.

The migratory locust, Locusta migratoria, is a worldwide agricultural pest whose outbreaks can result in plagues during which locusts mass migrate in marching bands and flying swarms. They exhibit striking phenotypic plasticity depending on the population density. They can transform between the barely-visible solitarious phase and the swarm-forming gregarious phase. A key step in the formation of large aggregation is the initial shift from the strong mutual aversion in solitarious locusts to the attraction to each other in gregarious locusts. Previous genomics studies have revealed 532 differentially expressed genes and many regulating small RNAs between the two phases of the migratory locust. Here, we developed a large-scale oligonucleotide microarray to quantify the expression of 9,154 genes during the time courses of isolating gregarious locusts and crowding solitarious locusts. We found that several olfaction-related genes, several CSPs and one takeout, are strongly regulated during both processes and display opposite expression trends in response to population density change. We then identified that these genes mediate the rapid switch of attraction/repulsion behaviors in the migratory locust. Our findings are particularly significant for understanding the process of locust aggregation and may provide new targets to manipulate locust behavior as part of novel management strategies.

Population transcriptomics of Drosophila melanogaster females

BACKGROUND

Variation at the level of gene expression is abundant in natural populations and is thought to contribute to the adaptive divergence of populations and species. Gene expression also differs considerably between males and females. Here we report a microarray analysis of gene expression variation among females of 16 Drosophila melanogaster strains derived from natural populations, including eight strains from the putative ancestral range in sub-Saharan Africa and eight strains from Europe. Gene expression variation among males of the same strains was reported previously.

RESULTS

We detected relatively low levels of expression polymorphism within populations, but much higher expression divergence between populations. A total of 569 genes showed a significant expression difference between the African and European populations at a false discovery rate of 5%. Genes with significant over-expression in Europe included the insecticide resistance gene Cyp6g1, as well as genes involved in proteolysis and olfaction. Genes with functions in carbohydrate metabolism and vision were significantly over-expressed in the African population. There was little overlap between genes expressed differently between populations in females and males.

CONCLUSIONS

Our results suggest that adaptive changes in gene expression have accompanied the out-of-Africa migration of D. melanogaster. Comparison of female and male expression data indicates that the vast majority of genes differing in expression between populations do so in only one sex and suggests that most regulatory adaptation has been sex-specific.

Investigations of transcript expression in fathead minnow (Pimephales promelas) brain tissue reveal toxicological impacts of RDX exposure

Production, usage and disposal of the munitions constituent (MC) cyclotrimethylenetrinitramine (RDX) has led to environmental releases on military facilities. The chemical attributes of RDX are conducive for leaching to surface water which may put aquatic organisms at risk of exposure. Because RDX has been observed to cause aberrant neuromuscular effects across a wide range of animal phyla, we assessed the effects of RDX on central nervous system (CNS) functions in the representative aquatic ecotoxicological model species, fathead minnow (Pimephales promelas). We developed a fathead minnow brain-tissue cDNA library enriched for transcripts differentially expressed in response to RDX and trinitrotoluene (TNT) exposure. All 4,128 cDNAs were sequenced, quality filtered and assembled yielding 2230 unique sequences and 945 significant blastx matches (E ≤10(-5)). The cDNA library was leveraged to create custom-spotted microarrays for use in transcript expression assays. The impact of RDX on transcript expression in brain tissue was examined in fathead minnows exposed to RDX at 0.625, 2.5, 5, 10mg/L or an acetone-spike control for 10 days. Overt toxicity of RDX in fathead minnow occurred only at the highest exposure concentration resulting in 50% mortality and weight loss. Conversely, Bayesian analysis of microarray data indicated significant changes in transcript expression at concentrations as low as 0.625 mg/L. In total, 154 cDNAs representing 44 unique transcripts were differentially expressed in RDX exposures, the majority of which were validated by reverse transcriptase-quantitative PCR (RT-qPCR). Investigation of molecular pathways, gene ontology (GO) and individual gene functions affected by RDX exposures indicated changes in metabolic processes involved in: oxygen transport, neurological function, calcium binding/signaling, energy metabolism, cell growth/division, oxidative stress and ubiquitination. In total, our study indicated that RDX exposure affected molecular processes critical to CNS function in fathead minnow.

Towards a more nuanced understanding of the relationship between sex-biased gene expression and rates of protein-coding sequence evolution

Genes that are differentially expressed between the sexes (sex-biased genes) are among the fastest evolving genes in animal genomes. The majority of sex-biased expression is attributable to genes that are primarily expressed in sex-limited reproductive tissues, and these reproductive genes are often rapidly evolving because of intra- and intersexual selection pressures. Additionally, studies of multiple taxa have revealed that genes with sex-biased expression are also expressed in a limited number of tissues. This is worth noting because narrowly expressed genes are known to evolve faster than broadly expressed genes. Therefore, it is not clear whether sex-biased genes are rapidly evolving because they have sexually dimorphic expression, because they are expressed in sex-limited reproductive tissues, or because they are narrowly expressed. To determine the extend to which other confounding variables can explain the rapid evolution of sex-biased genes, I analyzed the rates of evolution of sex-biased genes in Drosophila melanogaster and Mus musculus in light of tissue-specific measures of expression. I find that genes with sex-biased expression in somatic tissues shared by both sexes are often evolving faster than non-sex-biased genes, but this is best explained by the narrow expression profiles of sex-biased genes. Sex-biased genes in sex-limited tissues in D. melanogaster, however, evolve faster than other narrowly expressed genes. Therefore, the rapid evolution of sex-biased genes is limited only to those genes primarily expressed in sex-limited reproductive tissues.