Disruption of gene expression in hybrids of the fire ants Solenopsis invicta and Solenopsis richteri.

Transcriptome analysis is a powerful tool for unveiling the distribution and magnitude of genetic incompatibilities between hybridizing taxa. The nature of such incompatibilities is closely associated with the evolutionary histories of the parental species and may differ across tissues and between the sexes. In eusocial insects, the presence of castes that experience divergent selection regimes may result in additional distinct patterns of caste-specific hybrid incompatibilities. We analysed levels of expression of >14 000 genes in two life stages of each caste in the fire ants Solenopsis invicta and Solenopsis richteri and in their hybrids. We found strong contributions of both developmental stage and caste to gene expression patterns. In contrast, variability in expression was only weakly associated with taxonomic identity, with hybrid scores falling between those of the two parental species. Hybrid incompatibilities were surprisingly modest, with only 32 genes being mis-expressed, indicating low levels of disruption in gene regulation in hybrids; males and workers each mis-expressed at least seven times as many genes as queens. Interestingly, homologues of many of the mis-expressed genes have been implicated in behavioural variation in Drosophila melanogaster. General expression profiles of hybrids consistently were more similar to those of S. richteri than S. invicta, presumably because S. richteri trans-regulatory elements tend to be dominant and/or because there is an overall bias in the genetic composition of the hybrids towards S. richteri. Altogether, our results suggest that selection acting on each caste may contribute differently to interspecific divergence and speciation in this group of ants.

GLUTAMATE RECEPTOR-LIKE genes mediate leaf-to-leaf wound signalling.

Wounded leaves communicate their damage status to one another through a poorly understood process of long-distance signalling. This stimulates the distal production of jasmonates, potent regulators of defence responses. Using non-invasive electrodes we mapped surface potential changes in Arabidopsis thaliana after wounding leaf eight and found that membrane depolarizations correlated with jasmonate signalling domains in undamaged leaves. Furthermore, current injection elicited jasmonoyl-isoleucine accumulation, resulting in a transcriptome enriched in RNAs encoding key jasmonate signalling regulators. From among 34 screened membrane protein mutant lines, mutations in several clade 3 GLUTAMATE RECEPTOR-LIKE genes (GLRs 3.2, 3.3 and 3.6) attenuated wound-induced surface potential changes. Jasmonate-response gene expression in leaves distal to wounds was reduced in a glr3.3 glr3.6 double mutant. This work provides a genetic basis for investigating mechanisms of long-distance wound signalling in plants and indicates that plant genes related to those important for synaptic activity in animals function in organ-to-organ wound signalling.

Multilayered genetic and omics dissection of mitochondrial activity in a mouse reference population

The manner by which genotype and environment affect complex phenotypes is one of the fundamental questions in biology. In this study, we quantified the transcriptome--a subset of the metabolome--and, using targeted proteomics, quantified a subset of the liver proteome from 40 strains of the BXD mouse genetic reference population on two diverse diets. We discovered dozens of transcript, protein, and metabolite QTLs, several of which linked to metabolic phenotypes. Most prominently, Dhtkd1 was identified as a primary regulator of 2-aminoadipate, explaining variance in fasted glucose and diabetes status in both mice and humans. These integrated molecular profiles also allowed further characterization of complex pathways, particularly the mitochondrial unfolded protein response (UPR(mt)). UPR(mt) shows strikingly variant responses at the transcript and protein level that are remarkably conserved among C. elegans, mice, and humans. Overall, these examples demonstrate the value of an integrated multilayered omics approach to characterize complex metabolic phenotypes.

Comparative and functional genomics provide insights into the pathogenicity of dermatophytic fungi.

ABSTRACT: BACKGROUND: Millions of humans and animals suffer from superficial infections caused by a group of highly specialized filamentous fungi, the dermatophytes, which exclusively infect keratinized host structures. To provide broad insights into the molecular basis of the pathogenicity-associated traits, we report the first genome sequences of two closely phylogenetically related dermatophytes, Arthroderma benhamiae and Trichophyton verrucosum, both of which induce highly inflammatory infections in humans. RESULTS: 97% of the 22.5 megabase genome sequences of A. benhamiae and T. verrucosum are unambiguously alignable and collinear. To unravel dermatophyte-specific virulence-associated traits, we compared sets of potentially pathogenicity-associated proteins, such as secreted proteases and enzymes involved in secondary metabolite production, with those of closely related onygenales (Coccidioides species) and the mould Aspergillus fumigatus. The comparisons revealed expansion of several gene families in dermatophytes and disclosed the peculiarities of the dermatophyte secondary metabolite gene sets. Secretion of proteases and other hydrolytic enzymes by A. benhamiae was proven experimentally by a global secretome analysis during keratin degradation. Molecular insights into the interaction of A. benhamiae with human keratinocytes were obtained for the first time by global transcriptome profiling. Given that A. benhamiae is able to undergo mating, a detailed comparison of the genomes further unraveled the genetic basis of sexual reproduction in this species. CONCLUSIONS: Our results enlighten the genetic basis of fundamental and putatively virulence-related traits of dermatophytes, advancing future research on these medically important pathogens.

Landscape of transcription in human cells.

Eukaryotic cells make many types of primary and processed RNAs that are found either in specific subcellular compartments or throughout the cells. A complete catalogue of these RNAs is not yet available and their characteristic subcellular localizations are also poorly understood. Because RNA represents the direct output of the genetic information encoded by genomes and a significant proportion of a cell's regulatory capabilities are focused on its synthesis, processing, transport, modification and translation, the generation of such a catalogue is crucial for understanding genome function. Here we report evidence that three-quarters of the human genome is capable of being transcribed, as well as observations about the range and levels of expression, localization, processing fates, regulatory regions and modifications of almost all currently annotated and thousands of previously unannotated RNAs. These observations, taken together, prompt a redefinition of the concept of a gene.

The TetR-type MfsR protein of the integrative and conjugative element (ICE) ICEclc controls both a putative efflux system and initiation of ICE transfer.

Integrative and conjugating elements (ICE) are self-transferable DNAs widely present in bacterial genomes, which often carry a variety of auxiliary genes of potential adaptive benefit. One of the model ICE is ICEclc, an element originally found in Pseudomonas knackmussii B13 and known for its propensity to provide its host with the capacity to metabolize chlorocatechols and 2-aminophenol. In this work, we studied the mechanism and target of regulation of MfsR, a TetR-type repressor previously found to exert global control on ICEclc horizontal transfer. By using a combination of ICEclc mutant and transcriptome analysis, gene reporter fusions, and DNA binding assays, we found that MfsR is a repressor of both its own expression and that of a gene cluster putatively coding for a major facilitator superfamily efflux system on ICEclc (named mfsABC). Phylogenetic analysis suggests that mfsR was originally located immediately adjacent to the efflux pump genes but became displaced from its original cis target DNA by a gene insertion. This resulted in divergence of the original bidirectional promoters into two separated individual regulatory units. Deletion of mfsABC did not result in a strong phenotype, and despite screening a large number of compounds and conditions, we were unable to define the precise current function or target of the putative efflux pump. Our data reconstruct how the separation of an ancestor mfsR-mfsABC system led to global control of ICEclc transfer by MfsR.

Gene expression changes in chronic inflammatory demyelinating polyneuropathy skin biopsies.

Chronic-inflammatory demyelinating polyneuropathy (CIDP) is an immune-mediated disease with no known biomarkers for diagnosing the disease or assessing its prognosis. We performed transcriptional profiling microarray analysis on skin punch biopsies from 20 CIDP patients and 17 healthy controls to identify disease-associated gene expression changes. We demonstrate changes in expression of genes involved in immune and chemokine regulation, growth and repair. We also found a combination of two upregulated genes that can be proposed as a novel biomarker of the disorder.

Identification and functional characterization of Rca1, a transcription factor involved in both antifungal susceptibility and host response in Candida albicans.

The identification of novel transcription factors associated with antifungal response may allow the discovery of fungus-specific targets for new therapeutic strategies. A collection of 241 Candida albicans transcriptional regulator mutants was screened for altered susceptibility to fluconazole, caspofungin, amphotericin B, and 5-fluorocytosine. Thirteen of these mutants not yet identified in terms of their role in antifungal response were further investigated, and the function of one of them, a mutant of orf19.6102 (RCA1), was characterized by transcriptome analysis. Strand-specific RNA sequencing and phenotypic tests assigned Rca1 as the regulator of hyphal formation through the cyclic AMP/protein kinase A (cAMP/PKA) signaling pathway and the transcription factor Efg1, but also probably through its interaction with a transcriptional repressor, most likely Tup1. The mechanisms responsible for the high level of resistance to caspofungin and fluconazole observed resulting from RCA1 deletion were investigated. From our observations, we propose that caspofungin resistance was the consequence of the deregulation of cell wall gene expression and that fluconazole resistance was linked to the modulation of the cAMP/PKA signaling pathway activity. In conclusion, our large-scale screening of a C. albicans transcription factor mutant collection allowed the identification of new effectors of the response to antifungals. The functional characterization of Rca1 assigned this transcription factor and its downstream targets as promising candidates for the development of new therapeutic strategies, as Rca1 influences host sensing, hyphal development, and antifungal response.

Visualization and quality assessment of de novo genome assemblies.

Recent technological progress has greatly facilitated de novo genome sequencing. However, de novo assemblies consist in many pieces of contiguous sequence (contigs) arranged in thousands of scaffolds instead of small numbers of chromosomes. Confirming and improving the quality of such assemblies is critical for subsequent analysis. We present a method to evaluate genome scaffolding by aligning independently obtained transcriptome sequences to the genome and visually summarizing the alignments using the Cytoscape software. Applying this method to the genome of the red fire ant Solenopsis invicta allowed us to identify inconsistencies in 7%, confirm contig order in 20% and extend 16% of scaffolds.Scripts that generate tables for visualization in Cytoscape from FASTA sequence and scaffolding information files are publicly available at https://github.com/ksanao/TGNet.

Test of four colon cancer risk-scores in formalin fixed paraffin embedded microarray gene expression data.

BACKGROUND: Prognosis prediction for resected primary colon cancer is based on the T-stage Node Metastasis (TNM) staging system. We investigated if four well-documented gene expression risk scores can improve patient stratification. METHODS: Microarray-based versions of risk-scores were applied to a large independent cohort of 688 stage II/III tumors from the PETACC-3 trial. Prognostic value for relapse-free survival (RFS), survival after relapse (SAR), and overall survival (OS) was assessed by regression analysis. To assess improvement over a reference, prognostic model was assessed with the area under curve (AUC) of receiver operating characteristic (ROC) curves. All statistical tests were two-sided, except the AUC increase. RESULTS: All four risk scores (RSs) showed a statistically significant association (single-test, P < .0167) with OS or RFS in univariate models, but with HRs below 1.38 per interquartile range. Three scores were predictors of shorter RFS, one of shorter SAR. Each RS could only marginally improve an RFS or OS model with the known factors T-stage, N-stage, and microsatellite instability (MSI) status (AUC gains < 0.025 units). The pairwise interscore discordance was never high (maximal Spearman correlation = 0.563) A combined score showed a trend to higher prognostic value and higher AUC increase for OS (HR = 1.74, 95% confidence interval [CI] = 1.44 to 2.10, P < .001, AUC from 0.6918 to 0.7321) and RFS (HR = 1.56, 95% CI = 1.33 to 1.84, P < .001, AUC from 0.6723 to 0.6945) than any single score. CONCLUSIONS: The four tested gene expression-based risk scores provide prognostic information but contribute only marginally to improving models based on established risk factors. A combination of the risk scores might provide more robust information. Predictors of RFS and SAR might need to be different.

LRF-mediated Dll4 repression in erythroblasts is necessary for hematopoietic stem cell maintenance.

Hematopoietic stem cells (HSCs) are the most primitive cells in the hematopoietic system and are under tight regulation for self-renewal and differentiation. Notch signals are essential for the emergence of definitive hematopoiesis in mouse embryos and are critical regulators of lymphoid lineage fate determination. However, it remains unclear how Notch regulates the balance between HSC self-renewal and differentiation in the adult bone marrow (BM). Here we report a novel mechanism that prevents HSCs from undergoing premature lymphoid differentiation in BM. Using a series of in vivo mouse models and functional HSC assays, we show that leukemia/lymphoma related factor (LRF) is necessary for HSC maintenance by functioning as an erythroid-specific repressor of Delta-like 4 (Dll4) expression. Lrf deletion in erythroblasts promoted up-regulation of Dll4 in erythroblasts, sensitizing HSCs to T-cell instructive signals in the BM. Our study reveals novel cross-talk between HSCs and erythroblasts, and sheds a new light on the regulatory mechanisms regulating the balance between HSC self-renewal and differentiation.

Mice carrying ubiquitin-specific protease 2 (Usp2) gene inactivation maintain normal sodium balance and blood pressure.

Ubiquitylation plays an important role in the control of Na⁺ homeostasis by the kidney. It is well established that the epithelial Na⁺ channel ENaC is regulated by the ubiquitin-protein ligase NEDD4-2, limiting ENaC cell surface expression and activity. Ubiquitylation can be reversed by the action of deubiquitylating enzymes (DUBs). One such DUB, USP2-45, was identified previously as an aldosterone-induced protein in the kidney and is also a circadian output gene. In heterologous expression systems, USP2-45 binds to ENaC, deubiquitylates it, and enhances channel density and activity at the cell surface. Because the role of USP2-45 in renal Na⁺ transport had not been studied in vivo, we investigated here the effect of Usp2 gene inactivation in this process. We demonstrate first that USP2-45 protein has a rhythmic expression with a peak at ZT12. Usp2-KO mice did not show any differences from wild-type littermates with respect to the diurnal control of Na⁺ or K⁺ urinary excretion and plasma levels either on a standard diet or after acute and chronic changes to low- and high-Na⁺ diets, respectively. Moreover, they had similar aldosterone levels on either a low- or high-Na⁺ diet. Blood pressure measurements using telemetry did not reveal variations compared with control mice. Usp2-KO mice did not display alterations in expression of genes involved in sodium homeostasis or the ubiquitin system, as evidenced by transcriptome analysis in the kidney. Our data suggest that USP2 does not play a primary role in the control of Na⁺ balance or blood pressure.

The evolution of gene expression levels in mammalian organs.

Changes in gene expression are thought to underlie many of the phenotypic differences between species. However, large-scale analyses of gene expression evolution were until recently prevented by technological limitations. Here we report the sequencing of polyadenylated RNA from six organs across ten species that represent all major mammalian lineages (placentals, marsupials and monotremes) and birds (the evolutionary outgroup), with the goal of understanding the dynamics of mammalian transcriptome evolution. We show that the rate of gene expression evolution varies among organs, lineages and chromosomes, owing to differences in selective pressures: transcriptome change was slow in nervous tissues and rapid in testes, slower in rodents than in apes and monotremes, and rapid for the X chromosome right after its formation. Although gene expression evolution in mammals was strongly shaped by purifying selection, we identify numerous potentially selectively driven expression switches, which occurred at different rates across lineages and tissues and which probably contributed to the specific organ biology of various mammals.

Phenotypic switching in Pseudomonas brassicacearum involves GacS- and GacA-dependent Rsm small RNAs.

The plant-beneficial bacterium Pseudomonas brassicacearum forms phenotypic variants in vitro as well as in planta during root colonization under natural conditions. Transcriptome analysis of typical phenotypic variants using microarrays containing coding as well as noncoding DNA fragments showed differential expression of several genes relevant to secondary metabolism and of the small RNA (sRNA) genes rsmX, rsmY, and rsmZ. Naturally occurring mutations in the gacS-gacA system accounted for phenotypic switching, which was characterized by downregulation of antifungal secondary metabolites (2,4-diacetylphloroglucinol and cyanide), indoleacetate, exoenzymes (lipase and protease), and three different N-acyl-homoserine lactone molecules. Moreover, in addition to abrogating these biocontrol traits, gacS and gacA mutations resulted in reduced expression of the type VI secretion machinery, alginate biosynthesis, and biofilm formation. In a gacA mutant, the expression of rsmX was completely abolished, unlike that of rsmY and rsmZ. Overexpression of any of the three sRNAs in the gacA mutant overruled the pleiotropic changes and restored the wild-type phenotypes, suggesting functional redundancy of these sRNAs. In conclusion, our data show that phenotypic switching in P. brassicacearum results from mutations in the gacS-gacA system.

KAP1 regulates gene networks controlling T-cell development and responsiveness.

Chromatin remodeling at specific genomic loci controls lymphoid differentiation. Here, we investigated the role played in this process by Kruppel-associated box (KRAB)-associated protein 1 (KAP1), the universal cofactor of KRAB-zinc finger proteins (ZFPs), a tetrapod-restricted family of transcriptional repressors. T-cell-specific Kap1-deleted mice displayed a significant expansion of immature thymocytes, imbalances in CD4(+)/CD8(+) cell ratios, and altered responses to TCR and TGFβ stimulation when compared to littermate KAP1 control mice. Transcriptome and chromatin studies revealed that KAP1 binds T-cell-specific cis-acting regulatory elements marked by the H3K9me3 repressive mark and enriched in Ikaros/NuRD complexes. Also, KAP1 directly controls the expression of several genes involved in TCR and cytokine signaling. Among these, regulation of FoxO1 seems to play a major role in this system. Likely responsible for tethering KAP1 to at least part of its genomic targets, a small number of KRAB-ZFPs are selectively expressed in T-lymphoid cells. These results reveal the so far unsuspected yet important role of KAP1-mediated epigenetic regulation in T-lymphocyte differentiation and activation.