VarMixt: efficient variance modelling for the differential analysis of replicated gene expression data

Docosahexaenoic acid mechanisms of action on the bovine oocyte-cumulus complex

Background

Supplementation of bovine oocyte-cumulus complexes during in vitro maturation (IVM) with 1 μM of docosahexaenoic acid (DHA), C22:6 n-3 polyunsaturated fatty acid, was reported to improve in vitro embryo development. The objective of this paper was to decipher the mechanisms of DHA action.

Results

Transcriptomic analysis of 1 μM DHA-treated and control cumulus cells after 4 h IVM showed no significant difference in gene expression. MALDI-TOF mass spectrometry analysis of lipid profiles in DHA-treated and control oocytes and cumulus cells after IVM showed variations of only 3 out of 700 molecular species in oocytes and 7 out of 698 species in cumulus cells (p < 0.01).

We showed expression of free fatty acid receptor FFAR4 in both oocytes and cumulus cells, this receptor is known to be activated by binding to DHA. FFAR4 protein was localized close to the cellular membrane by immunofluorescence. Functional studies demonstrated that supplementation with FFAR4 agonist TUG-891 (1 μM or 5 μM) during IVM led to an increased blastocyst rate (39.5% ± 4.1%, 41.3% ± 4.1%), similar to DHA 1 μM treatment (39.2% ± 4.1%) as compared to control (25.2% ± 3.6%).

FFAR4 activation via TUG-891 led to beneficial effect on oocyte developmental competence and might explain in part similar effects of DHA.

Conclusions

In conclusion, we suggested that low dose of DHA (1 μM) during IVM might activate regulatory mechanisms without evident effect on gene expression and lipid content in oocyte-cumulus complexes, likely through signaling pathways which need to be elucidated in further studies.

Electronic supplementary material

The online version of this article (10.1186/s13048-017-0370-z) contains supplementary material, which is available to authorized users.

A clustering-based approach for efficient identification of microRNA combinatorial biomarkers

BACKGROUND

MicroRNAs (miRNAs) have great potential serving as tumor biomarkers and therapeutic targets. As the rapid development of high-throughput experimental technology, gene expression experiments have become more and more specialized and diversified. The complex data structure has brought great challenge for the identification of biomarkers. In the meantime, current statistical and machine learning methods for detecting biomarkers have the problem of low reliability and biased criteria.

RESULTS

This study aims to select combinatorial miRNA biomarkers, which have higher sensitivity and specificity than single-gene biomarkers. In order to avoid exhaustive search and redundant information, miRNAs are firstly clustered, then the combinations of representative cluster members are assessed as potential biomarkers. Both the criteria for the partition of clusters and selection of representative members are based on Fisher linear discriminant analysis (FDA). The FDA-based criterion has been demonstrated to be superior to three other criteria in selecting representative members, and also good at refining clusters. In the comparison with eight common feature selection methods, this clustering-based method performs the best with regard to the discriminative ability of selected biomarkers.

CONCLUSIONS

Our experimental results demonstrate that the clustering-based method can identify microRNA combinatorial biomarkers with high accuracy and efficiency. Our method and data are available to the public upon request.

The Legionella pneumophila genome evolved to accommodate multiple regulatory mechanisms controlled by the CsrA-system

The carbon storage regulator protein CsrA regulates cellular processes post-transcriptionally by binding to target-RNAs altering translation efficiency and/or their stability. Here we identified and analyzed the direct targets of CsrA in the human pathogen Legionella pneumophila. Genome wide transcriptome, proteome and RNA co-immunoprecipitation followed by deep sequencing of a wild type and a csrA mutant strain identified 479 RNAs with potential CsrA interaction sites located in the untranslated and/or coding regions of mRNAs or of known non-coding sRNAs. Further analyses revealed that CsrA exhibits a dual regulatory role in virulence as it affects the expression of the regulators FleQ, LqsR, LetE and RpoS but it also directly regulates the timely expression of over 40 Dot/Icm substrates. CsrA controls its own expression and the stringent response through a regulatory feedback loop as evidenced by its binding to RelA-mRNA and links it to quorum sensing and motility. CsrA is a central player in the carbon, amino acid, fatty acid metabolism and energy transfer and directly affects the biosynthesis of cofactors, vitamins and secondary metabolites. We describe the first L. pneumophila riboswitch, a thiamine pyrophosphate riboswitch whose regulatory impact is fine-tuned by CsrA, and identified a unique regulatory mode of CsrA, the active stabilization of RNA anti-terminator conformations inside a coding sequence preventing Rho-dependent termination of the gap operon through transcriptional polarity effects. This allows L. pneumophila to regulate the pentose phosphate pathway and the glycolysis combined or individually although they share genes in a single operon. Thus the L. pneumophila genome has evolved to acclimate at least five different modes of regulation by CsrA giving it a truly unique position in its life cycle.

The RNA binding protein CsrA is the master regulator of the bi-phasic life cycle of Legionella pneumophila governing virulence expression in this intracellular pathogen. Here, we have used deep sequencing of RNA enriched by co-immunoprecipitation with epitope-tagged CsrA to identify CsrA-associated transcripts at the genome level. We found 479 mRNAs or non-coding RNAs to be targets of CsrA. Among those major regulators including FleQ, the regulator of flagella expression, LqsR, the regulator of quorum sensing and RpoS implicated in stress response were identified. The expression of over 40 type IV secreted effector proteins important for intracellular survival and virulence are under the control of CsrA. Combined with transcriptomics, whole shotgun proteomics of a wild type and a CsrA mutant strain and functional analyses of several CsrA-targeted RNAs we identified the first riboswitch in L. pneumophila, a thiamine pyrophosphate riboswitch, and discovered a new mode of regulation by CsrA that allows L. pneumophila to regulate the pentose phosphate pathway and the glycolysis combined or individually although they share genes in a single operon. Our results further underline the indispensable role of CsrA in the life cycle of L. pneumophila and provide new insights into its regulatory roles and mechanisms.

A Unique cis-Encoded Small Noncoding RNA Is Regulating Legionella pneumophila Hfq Expression in a Life Cycle-Dependent Manner

Legionella pneumophila is an environmental bacterium that parasitizes protozoa, but it may also infect humans, thereby causing a severe pneumonia called Legionnaires’ disease. To cycle between the environment and a eukaryotic host, L. pneumophila is regulating the expression of virulence factors in a life cycle-dependent manner: replicating bacteria do not express virulence factors, whereas transmissive bacteria are highly motile and infective. Here we show that Hfq is an important regulator in this network. Hfq is highly expressed in transmissive bacteria but is expressed at very low levels in replicating bacteria. A L. pneumophila hfq deletion mutant exhibits reduced abilities to infect and multiply in Acanthamoeba castellanii at environmental temperatures. The life cycle-dependent regulation of Hfq expression depends on a unique cis-encoded small RNA named Anti-hfq that is transcribed antisense of the hfq transcript and overlaps its 5′ untranslated region. The Anti-hfq sRNA is highly expressed only in replicating L. pneumophila where it regulates hfq expression through binding to the complementary regions of the hfq transcripts. This results in reduced Hfq protein levels in exponentially growing cells. Both the small noncoding RNA (sRNA) and hfq mRNA are bound and stabilized by the Hfq protein, likely leading to the cleavage of the RNA duplex by the endoribonuclease RNase III. In contrast, after the switch to transmissive bacteria, the sRNA is not expressed, allowing now an efficient expression of the hfq gene and consequently Hfq. Our results place Hfq and its newly identified sRNA anti-hfq in the center of the regulatory network governing L. pneumophila differentiation from nonvirulent to virulent bacteria.

The abilities of L. pneumophila to replicate intracellularly and to cause disease depend on its capacity to adapt to different extra- and intracellular environmental conditions. Therefore, a timely and fine-tuned expression of virulence factors and adaptation traits is crucial. Yet, the regulatory circuits governing the life cycle of L. pneumophila from replicating to virulent bacteria are only partly uncovered. Here we show that the life cycle-dependent regulation of the RNA chaperone Hfq relies on a small regulatory RNA encoded antisense to the hfq-encoding gene through a base pairing mechanism. Furthermore, Hfq regulates its own expression in an autoregulatory loop. The discovery of this RNA regulatory mechanism in L. pneumophila is an important step forward in the understanding of how the switch from inoffensive, replicating to highly virulent, transmissive L. pneumophila is regulated.

Identifying specific proteins involved in eggshell membrane formation using gene expression analysis and bioinformatics

Background

The avian eggshell membranes surround the egg white and provide a structural foundation for calcification of the eggshell which is essential for avian reproduction; moreover, it is also a natural biomaterial with many potential industrial and biomedical applications. Due to the insoluble and stable nature of the eggshell membrane fibres, their formation and protein constituents remain poorly characterized. The purpose of this study was to identify genes encoding eggshell membrane proteins, particularly those responsible for its structural features, by analyzing the transcriptome of the white isthmus segment of the oviduct, which is the specialized region responsible for the fabrication of the membrane fibres.

Results

The Del-Mar 14 K chicken microarray was used to investigate up-regulated expression of transcripts in the white isthmus (WI) compared with the adjacent magnum (Ma) and uterine (Ut) segments of the hen oviduct. Analysis revealed 135 clones hybridizing to over-expressed transcripts (WI/Ma + WI/Ut), and corresponding to 107 NCBI annotated non-redundant Gallus gallus gene IDs. This combined analysis revealed that the structural proteins highly over-expressed in the white isthmus include collagen X (COL10A1), fibrillin-1 (FBN1) and cysteine rich eggshell membrane protein (CREMP). These results validate previous proteomics studies which have identified collagen X (α-1) and CREMP in soluble eggshell extracts. Genes encoding collagen-processing enzymes such as lysyl oxidase homologs 1, 2 and 3 (LOXL1, LOXL2 and LOXL3), prolyl 4 hydroxylase subunit α-2 and beta polypeptide (P4HA2 and P4HB) as well as peptidyl-prolyl cis-trans isomerase C (PPIC) were also over-expressed. Additionally, genes encoding proteins known to regulate disulfide cross-linking, including sulfhydryl oxidase (QSOX1) and thioredoxin (TXN), were identified which suggests that coordinated up-regulation of genes in the white isthmus is associated with eggshell membrane fibre formation.

Conclusions

The present study has identified genes associated with the processing of collagen, other structural proteins, and disulfide-mediated cross-linking during eggshell membrane formation in the white isthmus. Identification of these genes will provide new insight into eggshell membrane structure and mechanisms of formation that will assist in the development of selection strategies to improve eggshell quality and food safety of the table egg.

Electronic supplementary material

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

Ectopic lignification in the flax lignified bast fiber1 mutant stem is associated with tissue-specific modifications in gene expression and cell wall composition

Histochemical screening of a flax ethyl methanesulfonate population led to the identification of 93 independent M2 mutant families showing ectopic lignification in the secondary cell wall of stem bast fibers. We named this core collection the Linum usitatissimum (flax) lbf mutants for lignified bast fibers and believe that this population represents a novel biological resource for investigating how bast fiber plants regulate lignin biosynthesis. As a proof of concept, we characterized the lbf1 mutant and showed that the lignin content increased by 350% in outer stem tissues containing bast fibers but was unchanged in inner stem tissues containing xylem. Chemical and NMR analyses indicated that bast fiber ectopic lignin was highly condensed and rich in G-units. Liquid chromatography-mass spectrometry profiling showed large modifications in the oligolignol pool of lbf1 inner- and outer-stem tissues that could be related to ectopic lignification. Immunological and chemical analyses revealed that lbf1 mutants also showed changes to other cell wall polymers. Whole-genome transcriptomics suggested that ectopic lignification of flax bast fibers could be caused by increased transcript accumulation of (1) the cinnamoyl-CoA reductase, cinnamyl alcohol dehydrogenase, and caffeic acid O-methyltransferase monolignol biosynthesis genes, (2) several lignin-associated peroxidase genes, and (3) genes coding for respiratory burst oxidase homolog NADPH-oxidases necessary to increase H2O2 supply.

Gene expression profiling in Entamoeba histolytica identifies key components in iron uptake and metabolism

Entamoeba histolytica is an ameboid parasite that causes colonic dysentery and liver abscesses in humans. The parasite encounters dramatic changes in iron concentration during its invasion of the host, with relatively low levels in the intestinal lumen and then relatively high levels in the blood and liver. The liver notably contains sources of iron; therefore, the parasite's ability to use these sources might be relevant to its survival in the liver and thus the pathogenesis of liver abscesses. The objective of the present study was to identify factors involved in iron uptake, use and storage in E. histolytica. We compared the respective transcriptomes of E. histolytica trophozoites grown in normal medium (containing around 169 µM iron), low-iron medium (around 123 µM iron), iron-deficient medium (around 91 µM iron), and iron-deficient medium replenished with hemoglobin. The differentially expressed genes included those coding for the ATP-binding cassette transporters and major facilitator transporters (which share homology with bacterial siderophores and heme transporters) and genes involved in heme biosynthesis and degradation. Iron deficiency was associated with increased transcription of genes encoding a subset of cell signaling molecules, some of which have previously been linked to adaptation to the intestinal environment and virulence. The present study is the first to have assessed the transcriptome of E. histolytica grown under various iron concentrations. Our results provide insights into the pathways involved in iron uptake and metabolism in this parasite.

The gluconeogenesis pathway is involved in maintenance of enterohaemorrhagic Escherichia coli O157:H7 in bovine intestinal content

Enterohaemorrhagic Escherichia coli (EHEC) are responsible for outbreaks of food- and water-borne illness. The bovine gastrointestinal tract (GIT) is thought to be the principle reservoir of EHEC. Knowledge of the nutrients essential for EHEC growth and survival in the bovine intestine may help in developing strategies to limit their shedding in bovine faeces thus reducing the risk of human illnesses. To identify specific metabolic pathways induced in the animal GIT, the transcriptome profiles of EHEC O157:H7 EDL933 during incubation in bovine small intestine contents (BSIC) and minimal medium supplemented with glucose were compared. The transcriptome analysis revealed that genes responsible for the assimilation of ethanolamine, urea, agmatine and amino acids (Asp, Thr, Gly, Ser and Trp) were strongly up-regulated suggesting that these compounds are the main nitrogen sources for EHEC in BSIC. A central role for the gluconeogenesis pathway and assimilation of gluconeogenic substrates was also pinpointed in EHEC incubated in BSIC. Our results suggested that three amino acids (Asp, Ser and Trp), glycerol, glycerol 3-phosphate, L-lactate and C4-dicarboxylates are important carbon sources for EHEC in BSIC. The ability to use gluconeogenic substrates as nitrogen sources (amino acids) and/or carbon sources (amino acids, glycerol and lactate) may provide a growth advantage to the bacteria in intestinal fluids. Accordingly, aspartate (2.4 mM), serine (1.9 mM), glycerol (5.8 mM) and lactate (3.6 mM) were present in BSIC and may represent the main gluconeogenic substrates potentially used by EHEC. A double mutant of E. coli EDL933 defective for phosphoenolpyruvate synthase (PpsA) and phosphoenolpyruvate carboxykinase (PckA), unable to utilize tricarboxylic acid (TCA) intermediates was constructed. Growth competition experiments between EHEC EDL933 and the isogenic mutant strain in BSIC clearly showed a significant competitive growth advantage of the wild-type strain further illustrating the importance of the gluconeogenesis pathway in maintaining EHEC in the bovine GIT.

Identified members of the Streptomyces lividans AdpA regulon involved in differentiation and secondary metabolism

BACKGROUND

AdpA is a key transcriptional regulator involved in the complex growth cycle of Streptomyces. Streptomyces are Gram-positive bacteria well-known for their production of secondary metabolites and antibiotics. Most work on AdpA has been in S. griseus, and little is known about the pathways it controls in other Streptomyces spp. We recently discovered interplay between ClpP peptidases and AdpA in S. lividans. Here, we report the identification of genes directly regulated by AdpA in S. lividans.

RESULTS

Microarray experiments revealed that the expression of hundreds of genes was affected in a S. lividans adpA mutant during early stationary phase cultures in YEME liquid medium. We studied the expression of the S. lividans AdpA-regulated genes by quantitative real-time PCR analysis after various times of growth. In silico analysis revealed the presence of potential AdpA-binding sites upstream from these genes; electrophoretic mobility shift assays indicated that AdpA binds directly to their promoter regions. This work identifies new pathways directly controlled by AdpA and that are involved in S. lividans development (ramR, SLI7885 also known as hyaS and SLI6586), and primary (SLI0755-SLI0754 encoding CYP105D5 and Fdx4) or secondary (cchA, cchB, and hyaS) metabolism.

CONCLUSIONS

We characterised six S. lividans AdpA-dependent genes whose expression is directly activated by this pleiotropic regulator. Several of these genes are orthologous to bldA-dependent genes in S. coelicolor. Furthermore, in silico analysis suggests that over hundred genes may be directly activated or repressed by S. lividans AdpA, although few have been described as being part of any Streptomyces AdpA regulons. This study increases the number of known AdpA-regulated pathways in Streptomyces spp.

A Survey and Comparative Study of Statistical Tests for Identifying Differential Expression from Microarray Data

DNA microarray is a powerful technology that can simultaneously determine the levels of thousands of transcripts (generated, for example, from genes/miRNAs) across different experimental conditions or tissue samples. The motto of differential expression analysis is to identify the transcripts whose expressions change significantly across different types of samples or experimental conditions. A number of statistical testing methods are available for this purpose. In this paper, we provide a comprehensive survey on different parametric and non-parametric testing methodologies for identifying differential expression from microarray data sets. The performances of the different testing methods have been compared based on some real-life miRNA and mRNA expression data sets. For validating the resulting differentially expressed miRNAs, the outcomes of each test are checked with the information available for miRNA in the standard miRNA database PhenomiR 2.0. Subsequently, we have prepared different simulated data sets of different sample sizes (from 10 to 100 per group/population) and thereafter the power of each test have been calculated individually. The comparative simulated study might lead to formulate robust and comprehensive judgements about the performance of each test in the basis of assumption of data distribution. Finally, a list of advantages and limitations of the different statistical tests has been provided, along with indications of some areas where further studies are required.

Identification and characterization of the genetic changes responsible for the characteristic smooth-to-rough morphotype alterations of clinically persistent Mycobacterium abscessus

Mycobacterium abscessus is an emerging pathogen that is increasingly recognized as a relevant cause of human lung infection in cystic fibrosis patients. This highly antibiotic-resistant mycobacterium is an exception within the rapidly growing mycobacteria, which are mainly saprophytic and non-pathogenic organisms. M. abscessus manifests as either a smooth (S) or a rough (R) colony morphotype, which is of clinical importance as R morphotypes are associated with more severe and persistent infections. To better understand the molecular mechanisms behind the S/R alterations, we analysed S and R variants of three isogenic M. abscessus S/R pairs using an unbiased approach involving genome and transcriptome analyses, transcriptional fusions and integrating constructs. This revealed different small insertions, deletions (indels) or single nucleotide polymorphisms within the non-ribosomal peptide synthase gene cluster mps1-mps2-gap or mmpl4b in the three R variants, consistent with the transcriptional differences identified within this genomic locus that is implicated in the synthesis and transport of Glyco-Peptido-Lipids (GPL). In contrast to previous reports, the identification of clearly defined genetic lesions responsible for the loss of GPL-production or transport makes a frequent switching back-and-forth between smooth and rough morphologies in M. abscessus highly unlikely, which is important for our understanding of persistent M. abscessus infections.

Unilateral once daily milking locally induces differential gene expression in both mammary tissue and milk epithelial cells revealing mammary remodeling

Once daily milking reduces milk yield, but the underlying mechanisms are not yet fully understood. Local regulation due to milk stasis in the tissue may contribute to this effect, but such mechanisms have not yet been fully described. To challenge this hypothesis, one udder half of six Holstein dairy cows was milked once a day (ODM), and the other twice a day (TDM). On the 8th day of unilateral ODM, mammary epithelial cells (MEC) were purified from the milk using immunomagnetic separation. Mammary biopsies were harvested from both udder halves. The differences in transcript profiles between biopsies from ODM and TDM udder halves were analyzed by a 22k bovine oligonucleotide array, revealing 490 transcripts that were differentially expressed. The principal category of upregulated transcripts concerned mechanisms involved in cell proliferation and death. We further confirmed remodeling of the mammary tissue by immunohistochemistry, which showed less cell proliferation and more apoptosis in ODM udder halves. Gene expression analyzed by RT-qPCR in MEC purified from milk and mammary biopsies showed a common downregulation of six transcripts (ABCG2, FABP3, NUCB2, RNASE1 and 5, and SLC34A2) but also some discrepancies. First, none of the upregulated transcripts in biopsies varied in milk-purified MEC. Second, only milk-purified MEC showed significant LALBA downregulation, which suggests therefore that they correspond to a mammary epithelial cell subpopulation. Our results, obtained after unilateral milking, suggest that cell remodeling during ODM is due to a local effect, which may be triggered by milk accumulation.

Paralogous metabolism: S-alkyl-cysteine degradation in Bacillus subtilis

Metabolism is prone to produce analogs of essential building blocks in the cell (here named paralogous metabolism). The variants result from lack of absolute accuracy in enzyme-templated reactions as well as from molecular aging. If variants were left to accumulate, the earth would be covered by chemical waste. The way bacteria cope with this situation is essentially unexplored. To gain a comprehensive understanding of Bacillus subtilis sulphur paralogous metabolism, we used expression profiling with DNA arrays to investigate the changes in gene expression in the presence of S-methyl-cysteine (SMeC) and its close analog, methionine, as sole sulphur source. Altogether, more than 200 genes whose relative strength of induction was significantly different depending on the sulphur source used were identified. This allowed us to pinpoint operon ytmItcyJKLMNytmO_ytnIJ_rbfK_ytnLM as controlling the pathway cycling SMeC directly to cysteine, without requiring sulphur oxygenation. Combining genetic and physiological experiments, we deciphered the corresponding pathway that begins with protection of the metabolite by acetylation. Oxygenation of the methyl group then follows, and after deprotection (deacetylation), N-formyl cysteine is produced. This molecule is deformylated by the second deformylase present in B. subtilis DefB, yielding cysteine. This pathway appears to be present in plant-associated microbes.

Absence of cumulus cells during in vitro maturation affects lipid metabolism in bovine oocytes

Cumulus cells (CC) surround the oocyte and are coupled metabolically through regulation of nutrient intake. CC removal before in vitro maturation (IVM) decreases bovine oocyte developmental competence without affecting nuclear meiotic maturation. The objective was to investigate the influence of CC on oocyte cytoplasmic maturation in relation to energy metabolism. IVM with either cumulus-enclosed (CEO) or -denuded (DO) oocytes was performed in serum-free metabolically optimized medium. Transmission electron microscopy revealed different distribution of membrane-bound vesicles and lipid droplets between metaphase II DO and CEO. By Nile Red staining, a significant reduction in total lipid level was evidenced in DO. Global transcriptomic analysis revealed differential expression of genes regulating energy metabolism, transcription, and translation between CEO and DO. By Western blot, fatty acid synthase (FAS) and hormone-sensitive phospholipase (HSL) proteins were detected in oocytes and in CC, indicating a local lipogenesis and lypolysis. FAS protein was significantly less abundant in DO that in CEO and more highly expressed in CC than in the oocytes. On the contrary, HSL protein was more abundant in oocytes than in CC. In addition, active Ser⁵⁶³-phosphorylated HSL was detected in the oocytes only after IVM, and its level was similar in CEO and DO. In conclusion, absence of CC during IVM affected lipid metabolism in the oocyte and led to suboptimal cytoplasmic maturation. Thus, CC may influence the oocyte by orienting the consumption of nutritive storage via regulation of local fatty acid synthesis and lipolysis to provide energy for maturation.

Prevalence of gene expression additivity in genetically stable wheat allohexaploids

The reprogramming of gene expression appears as the major trend in synthetic and natural allopolyploids where expression of an important proportion of genes was shown to deviate from that of the parents or the average of the parents. In this study, we analyzed gene expression changes in previously reported, highly stable synthetic wheat allohexaploids that combine the D genome of Aegilops tauschii and the AB genome extracted from the natural hexaploid wheat Triticum aestivum. A comprehensive genome-wide analysis of transcriptional changes using the Affymetrix GeneChip Wheat Genome Array was conducted. Prevalence of gene expression additivity was observed where expression does not deviate from the average of the parents for 99.3% of 34,820 expressed transcripts. Moreover, nearly similar expression was observed (for 99.5% of genes) when comparing these synthetic and natural wheat allohexaploids. Such near-complete additivity has never been reported for other allopolyploids and, more interestingly, for other synthetic wheat allohexaploids that differ from the ones studied here by having the natural tetraploid Triticum turgidum as the AB genome progenitor. Our study gave insights into the dynamics of additive gene expression in the highly stable wheat allohexaploids.

In vitro maturation of oocytes alters gene expression and signaling pathways in bovine cumulus cells

In vitro maturation (IVM) of immature oocytes is widely used in assisted reproduction technologies in cattle, and is increasingly used to treat human infertility. The development competence of IVM oocytes, however, is lower than preovulatory, in vivo-matured oocytes. During maturation, cumulus cells (CC) are metabolically coupled with an oocyte and support the acquisition of its developmental potential. Our objective was to identify genes and pathways that were affected by IVM in bovine CC. Microarray transcriptomic analysis of CC enclosing in vitro- or in vivo-mature oocytes revealed 472 differentially expressed genes, including 28% related to apoptosis, correlating with twofold higher cell death after IVM than in vivo, as detected by TUNEL. Genes overexpressed after IVM were significantly enriched in functions involved in cell movement, focal adhesion, extracellular matrix function, and TGF-beta signaling, whereas under-expressed genes were enriched in regulating gene expression, energy metabolism, stress response, and MAP kinases pathway functions. Differential expression of 15 genes, including PAG11 (increased) and TXNIP (decreased), which were never detected in CC before, was validated by real-time RT-PCR. Moreover, protein quantification confirmed the lower abundance of glutathione S-transferase A1 and prostaglandin G/H synthase 2, and the higher abundance of hyaluronan synthase 2 and SMAD4, a member of TGF-beta pathway, in CC after IVM. Phosphorylation levels of SMAD2, MAPK3/1, and MAPK14, but not MAPK8, were higher after IVM that in vivo. In conclusion, IVM provokes the hyper-activation of TGF-beta and MAPK signaling components, modifies gene expression, leads to increased apoptosis in CC, and thus affects oocyte quality.

Nicotinic acid modulates Legionella pneumophila gene expression and induces virulence traits

In response to environmental fluctuations or stresses, bacteria can activate transcriptional and phenotypic programs to coordinate an adaptive response. The intracellular pathogen Legionella pneumophila converts from a noninfectious replicative form to an infectious transmissive form when the bacterium encounters alterations in either amino acid concentrations or fatty acid biosynthesis. Here, we report that L. pneumophila differentiation is also triggered by nicotinic acid, a precursor of the central metabolite NAD(+). In particular, when replicative L. pneumophila are treated with 5 mM nicotinic acid, the bacteria induce numerous transmissive-phase phenotypes, including motility, cytotoxicity toward macrophages, sodium sensitivity, and lysosome avoidance. Transcriptional profile analysis determined that nicotinic acid induces the expression of a panel of genes characteristic of transmissive-phase L. pneumophila. Moreover, an additional 213 genes specific to nicotinic acid treatment were altered. Although nearly 25% of these genes lack an assigned function, the gene most highly induced by nicotinic acid treatment encodes a putative major facilitator superfamily transporter, Lpg0273. Indeed, lpg0273 protects L. pneumophila from toxic concentrations of nicotinic acid as judged by analyzing the growth of the corresponding mutant. The broad utility of the nicotinic acid pathway to couple central metabolism and cell fate is underscored by this small metabolite's modulation of gene expression by diverse microbes, including Candida glabrata, Bordetella pertussis, Escherichia coli, and L. pneumophila.

New insights into sulfur metabolism in yeasts as revealed by studies of Yarrowia lipolytica

Yarrowia lipolytica, located at the frontier of hemiascomycetous yeasts and fungi, is an excellent candidate for studies of metabolism evolution. This yeast, widely recognized for its technological applications, in particular produces volatile sulfur compounds (VSCs) that fully contribute to the flavor of smear cheese. We report here a relevant global vision of sulfur metabolism in Y. lipolytica based on a comparison between high- and low-sulfur source supplies (sulfate, methionine, or cystine) by combined approaches (transcriptomics, metabolite profiling, and VSC analysis). The strongest repression of the sulfate assimilation pathway was observed in the case of high methionine supply, together with a large accumulation of sulfur intermediates. A high sulfate supply seems to provoke considerable cellular stress via sulfite production, resulting in a decrease of the availability of the glutathione pathway's sulfur intermediates. The most limited effect was observed for the cystine supply, suggesting that the intracellular cysteine level is more controlled than that of methionine and sulfate. Using a combination of metabolomic profiling and genetic experiments, we revealed taurine and hypotaurine metabolism in yeast for the first time. On the basis of a phylogenetic study, we then demonstrated that this pathway was lost by some of the hemiascomycetous yeasts during evolution.

Proteomic analysis of chloroplast-to-chromoplast transition in tomato reveals metabolic shifts coupled with disrupted thylakoid biogenesis machinery and elevated energy-production components

A comparative proteomic approach was performed to identify differentially expressed proteins in plastids at three stages of tomato (Solanum lycopersicum) fruit ripening (mature-green, breaker, red). Stringent curation and processing of the data from three independent replicates identified 1,932 proteins among which 1,529 were quantified by spectral counting. The quantification procedures have been subsequently validated by immunoblot analysis of six proteins representative of distinct metabolic or regulatory pathways. Among the main features of the chloroplast-to-chromoplast transition revealed by the study, chromoplastogenesis appears to be associated with major metabolic shifts: (1) strong decrease in abundance of proteins of light reactions (photosynthesis, Calvin cycle, photorespiration) and carbohydrate metabolism (starch synthesis/degradation), mostly between breaker and red stages and (2) increase in terpenoid biosynthesis (including carotenoids) and stress-response proteins (ascorbate-glutathione cycle, abiotic stress, redox, heat shock). These metabolic shifts are preceded by the accumulation of plastid-encoded acetyl Coenzyme A carboxylase D proteins accounting for the generation of a storage matrix that will accumulate carotenoids. Of particular note is the high abundance of proteins involved in providing energy and in metabolites import. Structural differentiation of the chromoplast is characterized by a sharp and continuous decrease of thylakoid proteins whereas envelope and stroma proteins remain remarkably stable. This is coincident with the disruption of the machinery for thylakoids and photosystem biogenesis (vesicular trafficking, provision of material for thylakoid biosynthesis, photosystems assembly) and the loss of the plastid division machinery. Altogether, the data provide new insights on the chromoplast differentiation process while enriching our knowledge of the plant plastid proteome.

Legionella pneumophila transcriptional response to chlorine treatment

Legionella pneumophila is a ubiquitous environmental microorganism found in freshwater that can cause an acute form of pneumonia known as Legionnaires' disease. Despite widespread use of chlorine to ensure drinking water quality and awareness that L. pneumophila may escape these treatments, little is known about its effects on L. pneumophila. The aim of this study was to investigate the L. pneumophila transcriptional response induced by chlorine treatment. Transcriptome analysis, using DNA arrays, showed that a sublethal dose of chlorine induces a differential expression of 391 genes involved in stress response, virulence, general metabolism, information pathways and transport. Many of the stress response genes were significantly upregulated, whereas a significant number of virulence genes were repressed. In particular, exposure of L. pneumophila to chlorine induced the expression of cellular antioxidant proteins, stress proteins and transcriptional regulators. In addition, glutathione S-transferase specific activity was enhanced following chlorine treatment. Our results clearly indicate that chlorine induces expression of proteins involved in cellular defence mechanisms against oxidative stress that might be involved in adaptation or resistance to chlorine treatment.

A functional genomic approach to the study of the milking ability in dairy sheep

To understand the mechanisms underlying milk ability and more precisely the kinetics of milk emission, we compared teat transcriptome profiles from Lacaune ewes in the tails of the milk flow phenotypic distribution. Two different arrays containing respectively 1896 and 13 168 PCR products selected from several tissue-specific cDNA libraries, including mammary gland, allowed the identification of 73 differentially expressed genes between teats from high and low milk flow ewes. Genes involved in muscle contraction were identified as over-expressed, and genes encoding collagen were found to be under-expressed in teats from low milk flow ewes. We confirmed this underexpression of COL1A1 and COL1A2 in low-milk flow ewes using RT-qPCR. These results suggest that milking ability may be due to the capacity of the teat sphincter to relax during mechanical milking. We propose that an optimal condition for mechanical milking may require proper relaxation of the teats. To our knowledge, this is the first transcriptomic analysis studying milking ability, using udder tissue for gene expression profiling, which demonstrates that mechanical milking ability is not only determined by morphological features but also by tissue composition.

Specific and extensive endometrial deregulation is present before conception in IVF/ICSI repeated implantation failures (IF) or recurrent miscarriages

The objective was to examine if IVF/ICSI repeated implantation failures (IF) or recurrent miscarriages (RM) could be related to preconceptional endometrial deregulations. IF was defined as the absence of pregnancy despite the transfer of at least ten IVF/ICSI good quality embryos, and RM as having at least three unexplained miscarriages. Fertile controls (FC) were women who had given birth at least once. Endometrial biopsy was performed in the mild luteal phase of a non-conceptual cycle (five women were selected in each group). Affymetrix chips (GeneChip Human Genome U133 Plus2.0 Array) were used for hybridization. Data were normalized by the gcRMA method, and raw p values adjusted by the Bonferroni procedure (1%). Differential expression of selected genes was analysed using real-time PCR. Gene networks and biological functions were explored using the Ingenuity Pathways Analysis software. Endometrial gene expression profiles at the time of uterine receptivity differ dramatically in the endometrium among FC, RM, and IF patients. Compared to FC, 2126 and 2477 genes are differentially expressed in IF and RM groups, respectively, and 2363 between IF and RM. In both conditions, differential gene expression referred mainly to DNA transcription and expression. Other main cellular functions deregulated in IF conditions correspond to cell morphology, cellular development, cell cycle, and cellular assembly, while in RM conditions, deregulated cellular functions relate to cell signalling (degradation of cyclic AMP and calcium metabolism) and cellular maintenance. In both conditions, there is an over-representation of deregulations related to the haematological system. In the IF condition, cell-mediated immune response and nervous system development and function are highly deregulated, while in RM patients, main deregulations are in organ and tissue development, humoral immune response, and muscular system development and function. Extensive endometrial deregulations are present before conception in patients who experienced IF or RM with both distinct and common deregulation.

Invading basement membrane matrix is sufficient for MDA-MB-231 breast cancer cells to develop a stable in vivo metastatic phenotype

INTRODUCTION

The poor efficacy of various anti-cancer treatments against metastatic cells has focused attention on the role of tumor microenvironment in cancer progression. To understand the contribution of the extracellular matrix (ECM) environment to this phenomenon, we isolated ECM surrogate invading cell populations from MDA-MB-231 breast cancer cells and studied their genotype and malignant phenotype.

METHODS

We isolated invasive subpopulations (INV) from non invasive populations (REF) using a 2D-Matrigel assay, a surrogate of basal membrane passage. INV and REF populations were investigated by microarray assay and for their capacities to adhere, invade and transmigrate in vitro, and to form metastases in nude mice.

RESULTS

REF and INV subpopulations were stable in culture and present different transcriptome profiles. INV cells were characterized by reduced expression of cell adhesion and cell-cell junction genes (44% of down regulated genes) and by a gain in expression of anti-apoptotic and pro-angiogenic gene sets. In line with this observation, in vitro INV cells showed reduced adhesion and increased motility through endothelial monolayers and fibronectin. When injected into the circulation, INV cells induced metastases formation, and reduced injected mice survival by up to 80% as compared to REF cells. In nude mice, INV xenografts grew rapidly inducing vessel formation and displaying resistance to apoptosis.

CONCLUSION

Our findings reveal that the in vitro ECM microenvironment per se was sufficient to select for tumor cells with a stable metastatic phenotype in vivo characterized by loss of adhesion molecules expression and induction of pro-angiogenic and survival factors.

Global regulation of the response to sulfur availability in the cheese-related bacterium Brevibacterium aurantiacum

In this study, we combined metabolic reconstruction, growth assays, and metabolome and transcriptome analyses to obtain a global view of the sulfur metabolic network and of the response to sulfur availability in Brevibacterium aurantiacum. In agreement with the growth of B. aurantiacum in the presence of sulfate and cystine, the metabolic reconstruction showed the presence of a sulfate assimilation pathway, thiolation pathways that produce cysteine (cysE and cysK) or homocysteine (metX and metY) from sulfide, at least one gene of the transsulfuration pathway (aecD), and genes encoding three MetE-type methionine synthases. We also compared the expression profiles of B. aurantiacum ATCC 9175 during sulfur starvation or in the presence of sulfate. Under sulfur starvation, 690 genes, including 21 genes involved in sulfur metabolism and 29 genes encoding amino acids and peptide transporters, were differentially expressed. We also investigated changes in pools of sulfur-containing metabolites and in expression profiles after growth in the presence of sulfate, cystine, or methionine plus cystine. The expression of genes involved in sulfate assimilation and cysteine synthesis was repressed in the presence of cystine, whereas the expression of metX, metY, metE1, metE2, and BL613, encoding a probable cystathionine-γ-synthase, decreased in the presence of methionine. We identified three ABC transporters: two operons encoding transporters were transcribed more strongly during cysteine limitation, and one was transcribed more strongly during methionine depletion. Finally, the expression of genes encoding a methionine γ-lyase (BL929) and a methionine transporter (metPS) was induced in the presence of methionine in conjunction with a significant increase in volatile sulfur compound production.

Overexpression of proinflammatory TLR-2-signalling lipoproteins in hypervirulent mycobacterial variants

Changes in the cell envelope composition of mycobacteria cause major changes in cytokine profiles of infected antigen presenting cells. We describe here the modulation of inflammatory responses by Mycobacterium abscessus, an emerging pathogen in cystic fibrosis. M. abscessus is able to switch from a smooth (S) to a rough (R) morphotype by the loss of a surface glycopeptidolipid. R variants are associated with severe clinical forms and a 'hyper-proinflammatory' response in ex vivo and in vivo models. Using partitioning of cell surface components we found that a complex fraction, more abundant in R variants than in S variants, made a major contribution to the TLR-2-dependent hyper-proinflammatory response induced by R variants. Lipoproteins were the main TLR-2 agonists in this fraction, consistent with the larger amounts of 16 lipoproteins in cell surface extracts from R variants; 15 out of 16 being more strongly induced in R variant than in S variant. Genetic interruption of glycopeptidolipid pathway in wild-type S variant resulted in R phenotype with similar induction of lipoprotein genes. In conclusion, R morphotype in M. abscessus is associated with increased synthesis/exposure at the cell surface of lipoproteins, these changes profoundly modifying the innate immune response through TLR-2-dependent mechanisms.

Two small ncRNAs jointly govern virulence and transmission in Legionella pneumophila

To transit from intra- to extracellular environments, Legionella pneumophila differentiates from a replicative/non-virulent to a transmissive/virulent form using the two-component system LetA/LetS and the global repressor protein CsrA. While investigating how both regulators act co-ordinately we characterized two ncRNAs, RsmY and RsmZ, that link the LetA/LetS and CsrA regulatory networks. We demonstrate that LetA directly regulates their expression and show that RsmY and RsmZ are functional in Escherichia coli and are able to bind CsrA in vitro. Single mutants have no (ΔrsmY) or a little (ΔrsmZ) impact on virulence, but the ΔrsmYZ strain shows a drastic defect in intracellular growth in Acanthamoeba castellanii and THP-1 monocyte-derived macrophages. Analysis of the transcriptional programmes of the ΔletA, ΔletS and ΔrsmYZ strains revealed that the switch to the transmissive phase is partially blocked. One major difference between the ΔletA, ΔletS and ΔrsmYZ strains was that the latter synthesizes flagella. Taken together, LetA activates transcription of RsmY and RsmZ, which sequester CsrA and abolish its post-transcriptional repressive activity. However, the RsmYZ-CsrA pathway appears not to be the main or only regulatory circuit governing flagella synthesis. We suggest that rather RpoS and LetA, by influencing LetE and probably cyclic-di-GMP levels, regulate motility in L. pneumophila.

Transient reduction of placental angiogenesis in PAI-1-deficient mice

Murine placentation is associated with the invasion of maternal endometrium by trophoblasts and an extensive maternal and fetal angiogenesis. Plasminogen activator inhibitor-1 (PAI-1) is transiently produced by spongiotrophoblasts and trophoblast giant cells at 10.5-11.5 days postcoitum (dpc). Knowing the key contribution of PAI-1 in the regulation of angiogenesis, we have now analyzed the consequence of PAI-1 deficiency on murine placentation. Morphological and quantitative computer-assisted image analysis revealed abnormal placental morphology in PAI-1-/- mice at 10.5 and 12.5 dpc. At 10.5 dpc, the genetic ablation of PAI-1 resulted in a transient reduction of both maternal and fetal vascularizations in the placenta and increased trophoblast cell density. This was associated with a poorer development of the labyrinth and an extension of the decidua. A larger spongiotrophoblast layer appeared at 12.5 dpc in PAI-1-deficient mice. Placental morphology was normalized at 14.5 dpc. Microarray analyses performed on laser capture microdissected labyrinths revealed that 46 genes were differentially expressed between the two genotypes at 10.5 dpc. However, only 11 genes were still differently modulated at 14.5 dpc, when normalization of placental morphology had taken place. This transcriptomic profiling highlighted a dysregulation in the expression of placenta-related cathepsin family members. Altogether our data provide evidence for a transient impaired placental morphology in PAI-1-deficient mice that is then normalized, leading to normal embryonic development.

Should we abandon the t-test in the analysis of gene expression microarray data: a comparison of variance modeling strategies

High-throughput post-genomic studies are now routinely and promisingly investigated in biological and biomedical research. The main statistical approach to select genes differentially expressed between two groups is to apply a t-test, which is subject of criticism in the literature. Numerous alternatives have been developed based on different and innovative variance modeling strategies. However, a critical issue is that selecting a different test usually leads to a different gene list. In this context and given the current tendency to apply the t-test, identifying the most efficient approach in practice remains crucial. To provide elements to answer, we conduct a comparison of eight tests representative of variance modeling strategies in gene expression data: Welch's t-test, ANOVA [1], Wilcoxon's test, SAM [2], RVM [3], limma [4], VarMixt [5] and SMVar [6]. Our comparison process relies on four steps (gene list analysis, simulations, spike-in data and re-sampling) to formulate comprehensive and robust conclusions about test performance, in terms of statistical power, false-positive rate, execution time and ease of use. Our results raise concerns about the ability of some methods to control the expected number of false positives at a desirable level. Besides, two tests (limma and VarMixt) show significant improvement compared to the t-test, in particular to deal with small sample sizes. In addition limma presents several practical advantages, so we advocate its application to analyze gene expression data.

Genome-wide gene expression changes in genetically stable synthetic and natural wheat allohexaploids

*The present study aims to understand regulation of gene expression in synthetic and natural wheat (Triticum aestivum) allohexaploids, that combines the AB genome of Triticum turgidum and the D genome of Aegilops tauschii; and which we have recently characterized as genetically stable. *We conducted a comprehensive genome-wide analysis of gene expression that allowed characterization of the effect of variability of the D genome progenitor, the intergenerational stability as well as the comparison with natural wheat allohexaploid. We used the Affymetrix GeneChip Wheat Genome Array, on which 55 049 transcripts are represented. *Additive expression was shown to represent the majority of expression regulation in the synthetic allohexaploids, where expression for more than c. 93% of transcripts was equal to the mid-parent value measured from a mixture of parental RNA. This leaves c. 2000 (c. 7%) transcripts, in which expression was nonadditive. No global gene expression bias or dominance towards any of the progenitor genomes was observed whereas high intergenerational stability and low effect of the D genome progenitor variability were revealed. *Our study suggests that gene expression regulation in wheat allohexaploids is established early upon allohexaploidization and highly conserved over generations, as demonstrated by the high similarity of expression with natural wheat allohexaploids.

The Legionella pneumophila LetA/LetS two-component system exhibits rheostat-like behavior

When confronted with metabolic stress, replicative Legionella pneumophila bacteria convert to resilient, infectious cells equipped for transmission. Differentiation is promoted by the LetA/LetS two-component system, which belongs to a family of signal-transducing proteins that employ a four-step phosphorelay to regulate gene expression. Histidine 307 of LetS was essential to switch on the transmission profile, but a threonine substitution at position 311 (T311M) suggested a rheostat-like function. The letS(T311M) bacteria resembled the wild type (WT) for some traits and letS null mutants for others, whereas they displayed intermediate levels of infectivity, cytotoxicity, and lysosome evasion. Although only 30 to 50% of letS(T311M) mutants became motile, flow cytometry determined that every cell eventually activated the flagellin promoter to WT levels, but expression was delayed. Likewise, letS(T311M) mutants exhibited delayed induction of RsmY and RsmZ, regulatory RNAs that relieve CsrA repression of transmission traits. Transcriptional profile analysis revealed that letS(T311M) mutants expressed the flagellar regulon and multiple other transmissive-phase loci at a higher cell density than the WT. Accordingly, we postulate that the letS(T311M) mutant may relay phosphate less efficiently than the WT LetS sensor protein, leading to sluggish gene expression and a variety of phenotypic profiles. Thus, as first described for BvgA/BvgS, rather than acting as on/off switches, this family of two-component systems exhibit rheostat activity that likely confers versatility as microbes adapt to fluctuating environments.

Whole intact rapeseeds or sunflower oil in high-forage or high-concentrate diets affects milk yield, milk composition, and mammary gene expression profile in goats

This study aimed to ascertain the response of goat mammary metabolic pathways to concentrate and lipid feeding in relation to milk fatty acid (FA) composition and secretion. Sixteen midlactation multiparous goats received diets differing in forage-to-concentrate ratio [high forage (HF) 64:36, and low forage (LF) 43:57] supplemented or not with lipids [HF with 130 g/d of oil from whole intact rapeseeds (RS) and LF with 130 g/d of sunflower oil (SO)] in a 4 x 4 Latin square design. Milk yield, milk composition, FA profile, and FA secretion were measured, as well as the expression profiles of key genes in mammary metabolism and of 8,382 genes, using a bovine oligonucleotide microarray. After 3 wk of treatment, milk, lactose, and protein yields were lower with HF-RS than with the other diets, whereas treatment had no effect on milk protein content. Milk fat content was higher with the HF-RS and LF-SO diets than with the HF and LF diets, and SO supplementation increased milk fat yield compared with the LF diet. Decreasing the forage-to-concentrate ratio from 64:36 to 43:57 had a limited effect on goat milk FA concentrations and secretions. Supplementing the LF diet with SO changed almost all the FA concentrations, including decreases in medium-chain saturated FA and large increases in trans C18:1 and C18:2 isomers (particularly trans-11 C18:1 and cis-9, trans-11 conjugated linoleic acid), without significant changes in C18:0 and cis-9 C18:1, whereas supplementing the HF diet with RS led to a strong decrease in short- and medium-chain saturated FA and a very strong increase in C18:0 and cis-9 C18:1, without significant changes in trans C18:1 and conjugated linoleic acid. Despite the decreases in milk lactose and protein yields observed with HF-RS, and despite the decrease in milk medium-chain FA and the increase in C18 FA secretion with RS or SO supplementation, none of the dietary treatments had any effect on mammary mRNA expression of the key genes involved in lactose (e.g., alpha-lactalbumin), protein (e.g., beta-casein), and lipid metabolism (e.g., lipoprotein lipase) after 3 wk of treatment. In addition, transcriptome analysis did not provide evidence of treatments inducing significant changes in the expression of specific genes in the mammary gland. However, 2-way hierarchical clustering analysis highlighted different global mammary expression profiles between diets, showing that the gene expression profiles corresponding to the same diet were gathered by common groups of genes. This experiment suggests that after 3 wk of dietary treatment, other factors, such as substrate availability for mammary metabolism, could play an important role in contributing to milk FA responses to changes in diet composition in the goat.

Analysis of the Legionella longbeachae genome and transcriptome uncovers unique strategies to cause Legionnaires' disease

Legionella pneumophila and L. longbeachae are two species of a large genus of bacteria that are ubiquitous in nature. L. pneumophila is mainly found in natural and artificial water circuits while L. longbeachae is mainly present in soil. Under the appropriate conditions both species are human pathogens, capable of causing a severe form of pneumonia termed Legionnaires' disease. Here we report the sequencing and analysis of four L. longbeachae genomes, one complete genome sequence of L. longbeachae strain NSW150 serogroup (Sg) 1, and three draft genome sequences another belonging to Sg1 and two to Sg2. The genome organization and gene content of the four L. longbeachae genomes are highly conserved, indicating strong pressure for niche adaptation. Analysis and comparison of L. longbeachae strain NSW150 with L. pneumophila revealed common but also unexpected features specific to this pathogen. The interaction with host cells shows distinct features from L. pneumophila, as L. longbeachae possesses a unique repertoire of putative Dot/Icm type IV secretion system substrates, eukaryotic-like and eukaryotic domain proteins, and encodes additional secretion systems. However, analysis of the ability of a dotA mutant of L. longbeachae NSW150 to replicate in the Acanthamoeba castellanii and in a mouse lung infection model showed that the Dot/Icm type IV secretion system is also essential for the virulence of L. longbeachae. In contrast to L. pneumophila, L. longbeachae does not encode flagella, thereby providing a possible explanation for differences in mouse susceptibility to infection between the two pathogens. Furthermore, transcriptome analysis revealed that L. longbeachae has a less pronounced biphasic life cycle as compared to L. pneumophila, and genome analysis and electron microscopy suggested that L. longbeachae is encapsulated. These species-specific differences may account for the different environmental niches and disease epidemiology of these two Legionella species.

Multiple controls affect arsenite oxidase gene expression in Herminiimonas arsenicoxydans

BACKGROUND

Both the speciation and toxicity of arsenic are affected by bacterial transformations, i.e. oxidation, reduction or methylation. These transformations have a major impact on environmental contamination and more particularly on arsenic contamination of drinking water. Herminiimonas arsenicoxydans has been isolated from an arsenic- contaminated environment and has developed various mechanisms for coping with arsenic, including the oxidation of As(III) to As(V) as a detoxification mechanism.

RESULTS

In the present study, a differential transcriptome analysis was used to identify genes, including arsenite oxidase encoding genes, involved in the response of H. arsenicoxydans to As(III). To get insight into the molecular mechanisms of this enzyme activity, a Tn5 transposon mutagenesis was performed. Transposon insertions resulting in a lack of arsenite oxidase activity disrupted aoxR and aoxS genes, showing that the aox operon transcription is regulated by the AoxRS two-component system. Remarkably, transposon insertions were also identified in rpoN coding for the alternative N sigma factor (sigma54) of RNA polymerase and in dnaJ coding for the Hsp70 co-chaperone. Western blotting with anti-AoxB antibodies and quantitative RT-PCR experiments allowed us to demonstrate that the rpoN and dnaJ gene products are involved in the control of arsenite oxidase gene expression. Finally, the transcriptional start site of the aoxAB operon was determined using rapid amplification of cDNA ends (RACE) and a putative -12/-24 sigma54-dependent promoter motif was identified upstream of aoxAB coding sequences.

CONCLUSION

These results reveal the existence of novel molecular regulatory processes governing arsenite oxidase expression in H. arsenicoxydans. These data are summarized in a model that functionally integrates arsenite oxidation in the adaptive response to As(III) in this microorganism.

Gene expression profiling to identify eggshell proteins involved in physical defense of the chicken egg

BACKGROUND

As uricoletic animals, chickens produce cleidoic eggs, which are self-contained bacteria-resistant biological packages for extra-uterine development of the chick embryo. The eggshell constitutes a natural physical barrier against bacterial penetration if it forms correctly and remains intact. The eggshell's remarkable mechanical properties are due to interactions among mineral components and the organic matrix proteins. The purpose of our study was to identify novel eggshell proteins by examining the transcriptome of the uterus during calcification of the eggshell. An extensive bioinformatic analysis on genes over-expressed in the uterus allowed us to identify novel eggshell proteins that contribute to the egg's natural defenses.

RESULTS

Our 14 K Del-Mar Chicken Integrated Systems microarray was used for transcriptional profiling in the hen's uterus during eggshell deposition. A total of 605 transcripts were over-expressed in the uterus compared with the magnum or white isthmus across a wide range of abundance (1.1- to 79.4-fold difference). The 605 highly-expressed uterine transcripts correspond to 469 unique genes, which encode 437 different proteins. Gene Ontology (GO) analysis was used for interpretation of protein function. The most over-represented GO terms are related to genes encoding ion transport proteins, which provide eggshell mineral precursors. Signal peptide sequence was found for 54 putative proteins secreted by the uterus during eggshell formation. Many functional proteins are involved in calcium binding or biomineralization--prerequisites for interacting with the mineral phase during eggshell fabrication. While another large group of proteins could be involved in proper folding of the eggshell matrix. Many secreted uterine proteins possess antibacterial properties, which would protect the egg against microbial invasion. A final group includes proteases and protease inhibitors that regulate protein activity in the acellular uterine fluid where eggshell formation takes place.

CONCLUSIONS

Our original study provides the first detailed description of the chicken uterus transcriptome during formation of the eggshell. We have discovered a cache of about 600 functional genes and identified a large number of encoded proteins secreted into uterine fluid for fabrication of the eggshell and chemical protection of the egg. Some of these uterine genes could prove useful as biological markers for genetic improvement of phenotypic traits (i.e., egg and eggshell quality).

Control of flagellar gene regulation in Legionella pneumophila and its relation to growth phase

The bacterial pathogen Legionella pneumophila responds to environmental changes by differentiation. At least two forms are well described: replicative bacteria are avirulent; in contrast, transmissive bacteria express virulence traits and flagella. Phenotypic analysis, Western blotting, and electron microscopy of mutants of the regulatory genes encoding RpoN, FleQ, FleR, and FliA demonstrated that flagellin expression is strongly repressed and that the mutants are nonflagellated in the transmissive phase. Transcriptome analyses elucidated that RpoN, together with FleQ, enhances transcription of 14 out of 31 flagellar class II genes, which code for the basal body, hook, and regulatory proteins. Unexpectedly, FleQ independent of RpoN enhances the transcription of fliA encoding sigma 28. Expression analysis of a fliA mutant showed that FliA activates three out of the five remaining flagellar class III genes and the flagellar class IV genes. Surprisingly, FleR does not induce but inhibits expression of at least 14 flagellar class III genes on the transcriptional level. Thus, we propose that flagellar class II genes are controlled by FleQ and RpoN, whereas the transcription of the class III gene fliA is controlled in a FleQ-dependent but RpoN-independent manner. However, RpoN and FleR might influence flagellin synthesis on a posttranscriptional level. In contrast to the commonly accepted view that enhancer-binding proteins such as FleQ always interact with RpoN to fullfill their regulatory functions, our results strongly indicate that FleQ regulates gene expression that is RpoN dependent and RpoN independent. Finally, FliA induces expression of flagellar class III and IV genes leading to the complete synthesis of the flagellum.

Terminal differentiation of goat mammary tissue during pregnancy requires the expression of genes involved in immune functions

Terminal differentiation of mammary tissue into a functional epithelium that synthesizes and secretes milk occurs during pregnancy. The molecular mechanisms underlying this complex process are poorly understood, especially in ruminants. To obtain an overview of the ruminant mammary gland's final differentiation process, we conducted time-course gene expression analysis of five physiological stages: four during pregnancy (P46, P70, P90, and P110) and one after 40 days of lactation (L40). An appropriate loop experimental design was used to follow gene expression profiles. Using three nulliparous (pregnancy) or primiparous (lactation) goats per stage, we performed a comparison starting from nine dye-swaps and using a 22K bovine oligoarray. Statistical analysis revealed that the expression of 1,696 genes varied significantly at least once in the study. These genes fell into 19 clusters based on their expression profiles. Identification of biological functions with Ingenuity Pathway Analysis software revealed several similarities, in keeping with physiological stages described in mice. As in mice, expression of milk protein genes began at midpregnancy, and genes regulating lipid biosynthesis were induced at the onset of lactation. During the first half of pregnancy, the molecular signature of goat mammary tissue was characterized by the expression of genes associated with tissue remodeling and differentiation, while the second half was mainly characterized by the presence of messengers encoding genes involved in cell proliferation. A large number of immune-related genes were also induced, supporting recent speculation that mammary tissue has an original immune function, and the recruitment of migrating hematopoietic cells possibly involved in the branching morphogenesis of the mammary gland. These data hint that the induction of differentiation occurs early in pregnancy, very likely before P46. This period is therefore crucial for obtaining a healthy and productive mammary gland.

Changes in transcriptome after in vivo exposure to ionising radiation reveal a highly specialised liver response

PURPOSE

To identify transcriptional gene-networks involved in the early in vivo response of liver cells to radiation exposure and improve our understanding of the molecular processes responsible for tissue radiosensitivity.

MATERIALS AND METHODS

Transcriptome variations of liver RNA samples were measured 3 hours post-irradiation using microarray technology. The results were confirmed and extended using real-time polymerase-chain-reaction (RT-PCR).

RESULTS

We identified quantitative changes in the expression of 126 genes, most of which were observed for the first time. We show that some modifications, such as the upregulation of the cyclin-dependent kinase inhibitor 1A (Cdkn1A) gene, persisted for at least two months after the initial exposure. Other genes regulated by the transformation-related protein 53 (Trp53/p53) such as Bcl2-associated X protein (Bax) or etoposide-induced-2.4 (Ei24/PIG8) were not upregulated. Grouping differentially expressed genes into functional categories revealed that the primary response of liver cells to radiation exposure was the enhancement of oxidoreductase activity and inhibition of cell proliferation, involving cell cycle progression and apoptosis-related genes.

CONCLUSIONS

The data provides evidence of gene expression modifications associated with the hepatic response to radiation exposure. One of the main differences observed with radiation-sensitive tissues such as the spleen was cell proliferation. The comparison of our data with transcriptome modifications in different biological models enabled the identification of networks of genes that might be co-regulated. Overall, our expression data revealed genes and cellular pathways that might help to improve our understanding of the molecular basis underlying tissue radiosensitivity and to identify possible targets for novel therapeutic strategies.

An evaluation framework for statistical tests on microarray data

Microarray analysis has become a popular and routine method in functional genomics. It is typical for such experiments to involve a small number of replicates, which causes unreliable estimates of the sample variance. Microarrays have fostered the development of new statistical methods to analyze data resulting from experiments with small sample sizes. In this study, we tackle the problem of evaluating the performance of statistical tests for generating ranked gene lists from two-channel direct comparisons. We propose an evaluation method based on a oligonucleotide microarray with a large number of replicate spots yielding a maximum of 400 replicates per gene. We apply Spearman's rank correlation coefficient to ranked gene-lists generated by eight widely used microarray specific test statistics, which are applied to small random samples. We could show that variance stabilizing methods such as Cyber-T, SAM, and LIMMA can be beneficial for very small sample sizes and that SAM and the t-test provide stronger control of the type I error rate than the other methods. Specifically, we report that for four replicates all methods reach a high to very high correlation with our reference standard.

The global regulator CodY regulates toxin gene expression in Bacillus anthracis and is required for full virulence

In gram-positive bacteria, CodY is an important regulator of genes whose expression changes upon nutrient limitation and acts as a repressor of virulence gene expression in some pathogenic species. Here, we report the role of CodY in Bacillus anthracis, the etiologic agent of anthrax. Disruption of codY completely abolished virulence in a toxinogenic, noncapsulated strain, indicating that the activity of CodY is required for full virulence of B. anthracis. Global transcriptome analysis of a codY mutant and the parental strain revealed extensive differences. These differences could reflect direct control for some genes, as suggested by the presence of CodY binding sequences in their promoter regions, or indirect effects via the CodY-dependent control of other regulatory proteins or metabolic rearrangements in the codY mutant strain. The differences included reduced expression of the anthrax toxin genes in the mutant strain, which was confirmed by lacZ reporter fusions and immunoblotting. The accumulation of the global virulence regulator AtxA protein was strongly reduced in the mutant strain. However, in agreement with the microarray data, expression of atxA, as measured using an atxA-lacZ transcriptional fusion and by assaying atxA mRNA, was not significantly affected in the codY mutant. An atxA-lacZ translational fusion was also unaffected. Overexpression of atxA restored toxin component synthesis in the codY mutant strain. These results suggest that CodY controls toxin gene expression by regulating AtxA accumulation posttranslationally.

Gene expression profiles of bovine caruncular and intercaruncular endometrium at implantation

At implantation the endometrium undergoes modifications necessary for its physical interactions with the trophoblast as well as the development of the conceptus. We aim to identify endometrial factors and pathways essential for a successful implantation in the caruncular (C) and the intercaruncular (IC) areas in cattle. Using a 13,257-element bovine oligonucleotide array, we established expression profiles at day 20 of the estrous cycle or pregnancy (implantation), revealing 446 and 1,295 differentially expressed genes (DEG) in C and IC areas, respectively (false discovery rate = 0.08). The impact of the conceptus was higher on the immune response function in C but more prominent on the regulation of metabolism function in IC. The C vs. IC direct comparison revealed 1,177 and 453 DEG in cyclic and pregnant animals respectively (false discovery rate = 0.05), with a major impact of the conceptus on metabolism and cell adhesion. We selected 15 genes including C11ORF34, CXCL12, CXCR4, PLAC8, SCARA5, and NPY and confirmed their differential expression by quantitative RT-PCR. The cellular localization was analyzed by in situ hybridization and, upon pregnancy, showed gene-specific patterns of cell distribution, including a high level of expression in the luminal epithelium for C11ORF34 and MX1. Using primary cultures of bovine endometrial cells, we identified PTN, PLAC8, and CXCL12 as interferon-tau (IFNT) target genes and MSX1 and CXCR7 as IFNT-regulated genes, whereas C11ORF34 was not an IFNT-regulated gene. Our transcriptomic data provide novel molecular insights accounting for the biological functions related to the C or IC endometrial areas and may contribute to the identification of potential biomarkers for normal and perturbed early pregnancy.

Variance estimation in the analysis of microarray data

Microarrays are one of the most widely used high throughput technologies. One of the main problems in the area is that conventional estimates of the variances that are required in the t-statistic and other statistics are unreliable owing to the small number of replications. Various methods have been proposed in the literature to overcome this lack of degrees of freedom problem. In this context, it is commonly observed that the variance increases proportionally with the intensity level, which has led many researchers to assume that the variance is a function of the mean. Here we concentrate on estimation of the variance as a function of an unknown mean in two models: the constant coefficient of variation model and the quadratic variance-mean model. Because the means are unknown and estimated with few degrees of freedom, naive methods that use the sample mean in place of the true mean are generally biased because of the errors-in-variables phenomenon. We propose three methods for overcoming this bias. The first two are variations on the theme of the so-called heteroscedastic simulation-extrapolation estimator, modified to estimate the variance function consistently. The third class of estimators is entirely different, being based on semiparametric information calculations. Simulations show the power of our methods and their lack of bias compared with the naive method that ignores the measurement error. The methodology is illustrated by using microarray data from leukaemia patients.

Postgenomic analysis of streptococcus thermophilus cocultivated in milk with Lactobacillus delbrueckii subsp. bulgaricus: involvement of nitrogen, purine, and iron metabolism

Streptococcus thermophilus is one of the most widely used lactic acid bacteria in the dairy industry, in particular in yoghurt manufacture, where it is associated with Lactobacillus delbrueckii subsp. bulgaricus. This bacterial association, known as a proto-cooperation, is poorly documented at the molecular and regulatory levels. We thus investigate the kinetics of the transcriptomic and proteomic modifications of S. thermophilus LMG 18311 in response to the presence of L. delbrueckii subsp. bulgaricus ATCC 11842 during growth in milk at two growth stages. Seventy-seven different genes or proteins (4.1% of total coding sequences), implicated mainly in the metabolism of nitrogen (24%), nucleotide base (21%), and iron (20%), varied specifically in coculture. One of the most unpredicted results was a significant decrease of most of the transcripts and enzymes involved in purine biosynthesis. Interestingly, the expression of nearly all genes potentially encoding iron transporters of S. thermophilus decreased, whereas that of iron-chelating dpr as well as that of the fur (perR) regulator genes increased, suggesting a reduction in the intracellular iron concentration, probably in response to H(2)O(2) production by L. bulgaricus. The present study reveals undocumented nutritional exchanges and regulatory relationships between the two yoghurt bacteria, which provide new molecular clues for the understanding of their associative behavior.

Endometrium as an early sensor of in vitro embryo manipulation technologies

Implantation is crucial for placental development that will subsequently impact fetal growth and pregnancy success with consequences on postnatal health. We postulated that the pattern of genes expressed by the endometrium when the embryo becomes attached to the mother uterus could account for the final outcome of a pregnancy. As a model, we used the bovine species where the embryo becomes progressively and permanently attached to the endometrium from day 20 of gestation onwards. At that stage, we compared the endometrial genes profiles in the presence of an in vivo fertilized embryo (AI) with the endometrial patterns obtained in the presence of nuclear transfer (SCNT) or in vitro fertilized embryos (IVF), both displaying lower and different potentials for term development. Our data provide evidence that the endometrium can be considered as a biological sensor able to fine-tune its physiology in response to the presence of embryos whose development will become altered much later after the implantation process. Compared with AI, numerous biological functions and several canonical pathways with a major impact on metabolism and immune function were found to be significantly altered in the endometrium of SCNT pregnancies at implantation, whereas the differences were less pronounced with IVF embryos. Determining the limits of the endometrial plasticity at the onset of implantation should bring new insights on the contribution of the maternal environment to the development of an embryo and the success of pregnancy.

Kerfdr: a semi-parametric kernel-based approach to local false discovery rate estimation

Background

The use of current high-throughput genetic, genomic and post-genomic data leads to the simultaneous evaluation of a large number of statistical hypothesis and, at the same time, to the multiple-testing problem. As an alternative to the too conservative Family-Wise Error-Rate (FWER), the False Discovery Rate (FDR) has appeared for the last ten years as more appropriate to handle this problem. However one drawback of FDR is related to a given rejection region for the considered statistics, attributing the same value to those that are close to the boundary and those that are not. As a result, the local FDR has been recently proposed to quantify the specific probability for a given null hypothesis to be true.

Results

In this context we present a semi-parametric approach based on kernel estimators which is applied to different high-throughput biological data such as patterns in DNA sequences, genes expression and genome-wide association studies.

Conclusion

The proposed method has the practical advantages, over existing approaches, to consider complex heterogeneities in the alternative hypothesis, to take into account prior information (from an expert judgment or previous studies) by allowing a semi-supervised mode, and to deal with truncated distributions such as those obtained in Monte-Carlo simulations. This method has been implemented and is available through the R package kerfdr via the CRAN or at .