The early Miocene rotation of the Corso-Sardinian block. New paleomagnetic constraints for the end of the motion

The paleomagnetic investigations carried out in the 70's on Oligo-Miocene volcanics of Sardinia have demonstrated that the island was turned by 35-30 degrees clockwise from 33 Ma up to 3-1-20.5 Ma and rotated counterclockwise in a few million years [De Jong et al., 1969, 1973; Bobier et Coulon, 1970; Coulon et al., 1974; Manzoni, 1974, 1975; Bellon rr nl.. 1977: Edel et Lortscher, 1977; Edel, 1979, 1980]. Since then, the end of the rotation fixed at 19 Ma by Montigny er al. [1981] was the subject of discussions and several studies associating paleomagnetism and radiometric dating were undertaken [Assorgia er al., 1994: Vigliotti et Langenheim, 1995: Deino et al., 1997; Gattacceca rt Deino, 1999]. This is a contribution to this debate that is hampered by thr important secular variation recorded in the volcanics. The only way to get our of this problem is to sample series of successive flows as completely as possible, and to reduce the effect of secular variation by the calculation of means. Sampling was performed north of Bonorva in 5 pyroclastic flows that belong to the upper ignimbritic series SI2 according to Coulon rr nl. [1974] or LBLS according to Assorgia et al. [1997] (fig. I). Ar-40/Ar-39 dating of biotites from the debris flow (MDF) has yielded an age or 18.35 +/- 0.03 Ma [Dubois, 2000]. Five of the investigated sites are located beneath the debris flow ITV, TVB, TVD, SPM85, SPM86), one site was cured in the matrix of the debris flow (MDF) and one in 4 metric blocks included in the flow (DFC). Another site was sampled in the upper ash flow (PDM) that marks the end of the pyroclastic activity, just before the marine transgression. According to micropaleontological and radiometric dating this transgression has occurred between 18.35 and 17.6 Ma [Dubois, 2000]. After removal of a soft viscous component, the thermal demagnetization generally shows a univectorial behaviour of the remanent magnetization (fig. 2a). The maximum unblocking temperatures of 580-620 degrees (tab. I) and a rapid saturation below 100 mT (fig. 3) indicate that the carrier of the characteristic magnetization is magnetite. The exception comes: from the upper site PDM in which were found two characteristic components, one with a normal polarity and low unblocking temperatures up to 350 degreesC and one with a reversed polarity and maximum unblocking temperatures at 580-600 degreesC of magnetite. After calculation of a mean direction for each flow, the mean << Al >> direction 4 degrees /57 degrees (alpha (95) = 13 degrees) computed with the mean directions for the 5 flows may be considered as weakly affected by secular variation. But the results require a more careful examination. The declinations are N to NNW beneath the debris flow. NNW in the debris flow. and NNE (or SSW) above the debris flow, The elongated distribution of the directions obtained at sites TVB and TVD. scattered from the mean direction of TV to the mean direction of MDF is interpreted as due to partial overprinting during the debris How volcanic episode, The low temperature component PDMa is likely related to the alteration seen on thin sections and is also viewed as an overprint. As NNE/SSW directions occur as well below (mean direction << B >> : 5 degrees /58 degrees) as above the debris flow (PDMb : 200 degrees/-58 degrees). the NNW directions (<< C >> : 337 degrees /64 degrees) associated with the debris flow volcanism may be interpreted as resulting from a magnetic field excursion. According to the polarity scale of Cande and Kent [1992, 1995] and the radiometric age of MDF, the directions with normal polarity (TV, TVB, TVD, SPM85. SPM86a. MDF. DFC) may represent the period 5En. while the directions with reversed polarity PDMb and SPM86b were likely acquired during the period 5Dr. Using the mean << Al >> direction, the mean << B >>, or the PDM direction (tab. I). the deviation in declination with the direction of stable Europe 6.4 degrees /58.7 degrees (alpha (95) = 8 degrees) for a selection of 4 middle Tertiary poles by Besse et Courtillot [1991] or 7 degrees /56 degrees (alpha (95) = 3 degrees) for 19 poles listed by Edel [1980] can be considered as negligible. Using the results from the uppermost ignimbritic layer of Anglona also emplaced around 18.3 Ma [Odin rt al.. 1994]. the mean direction << E >> (3 degrees /51.5 degrees) leads to the same conclusion. On the contrary, when taking into account all dated results available for the period 5En (mean direction << D >> 353 degrees /56 degrees for 45 sites) (tab. II). the deviation 13 degrees is much more significant. As the rotation of Sardinia started around 21-20.5 Ma. the assumption of a constant velocity of rotation and the deviations of the Sardinia directions with respect to the stable Europe direction locate the end of the motion between 18.3 and 17.2 or 16.7 Ma (fig. 4). During the interval 18.35-17.5 Ma, the marine transgression took place. At the same period a NE-SW shortening interpreted as resulting from the collision of Sardinia with Apulia affected different parts of the island [Letouzey et al., 1982]. Consequently, the new paleomagnetic results and the tectono-sedimentary evolution are in favour of an end of the rotation at 17.5-18 Ma.

Pseudogenes in the ENCODE regions: consensus annotation, analysis of transcription, and evolution

Arising from either retrotransposition or genomic duplication of functional genes, pseudogenes are “genomic fossils” valuable for exploring the dynamics and evolution of genes and genomes. Pseudogene identification is an important problem in computational genomics, and is also critical for obtaining an accurate picture of a genome’s structure and function. However, no consensus computational scheme for defining and detecting pseudogenes has been developed thus far. As part of the ENCyclopedia Of DNA Elements (ENCODE) project, we have compared several distinct pseudogene annotation strategies and found that different approaches and parameters often resulted in rather distinct sets of pseudogenes. We subsequently developed a consensus approach for annotating pseudogenes (derived from protein coding genes) in the ENCODE regions, resulting in 201 pseudogenes, two-thirds of which originated from retrotransposition. A survey of orthologs for these pseudogenes in 28 vertebrate genomes showed that a significant fraction (∼80%) of the processed pseudogenes are primate-specific sequences, highlighting the increasing retrotransposition activity in primates. Analysis of sequence conservation and variation also demonstrated that most pseudogenes evolve neutrally, and processed pseudogenes appear to have lost their coding potential immediately or soon after their emergence. In order to explore the functional implication of pseudogene prevalence, we have extensively examined the transcriptional activity of the ENCODE pseudogenes. We performed systematic series of pseudogene-specific RACE analyses. These, together with complementary evidence derived from tiling microarrays and high throughput sequencing, demonstrated that at least a fifth of the 201 pseudogenes are transcribed in one or more cell lines or tissues.

The DART classification of unannotated transcription within the ENCODE regions: associating transcription with known and novel loci

For the ∼1% of the human genome in the ENCODE regions, only about half of the transcriptionally active regions (TARs) identified with tiling microarrays correspond to annotated exons. Here we categorize this large amount of “unannotated transcription.” We use a number of disparate features to classify the 6988 novel TARs—array expression profiles across cell lines and conditions, sequence composition, phylogenetic profiles (presence/absence of syntenic conservation across 17 species), and locations relative to genes. In the classification, we first filter out TARs with unusual sequence composition and those likely resulting from cross-hybridization. We then associate some of those remaining with proximal exons having correlated expression profiles. Finally, we cluster unclassified TARs into putative novel loci, based on similar expression and phylogenetic profiles. To encapsulate our classification, we construct a Database of Active Regions and Tools (DART.gersteinlab.org). DART has special facilities for rapidly handling and comparing many sets of TARs and their heterogeneous features, synchronizing across builds, and interfacing with other resources. Overall, we find that ∼14% of the novel TARs can be associated with known genes, while ∼21% can be clustered into ∼200 novel loci. We observe that TARs associated with genes are enriched in the potential to form structural RNAs and many novel TAR clusters are associated with nearby promoters. To benchmark our classification, we design a set of experiments for testing the connectivity of novel TARs. Overall, we find that 18 of the 46 connections tested validate by RT-PCR and four of five sequenced PCR products confirm connectivity unambiguously.

Hypothalamus transcriptome profile suggests an anorexia-cachexia syndrome in the anx/anx mouse model

The anx/anx mouse displays poor appetite and lean appearance and is considered a good model for the study of anorexia nervosa. To identify new genes involved in feeding behavior and body weight regulation we performed an expression profiling in the hypothalamus of the anx/anx mice. Using commercial microarrays we detected 156 differentially expressed genes and validated 92 of those using TaqMan low-density arrays. The expression of a set of 87 candidate genes selected based on literature evidences was also quantified by TaqMan low-density arrays. Our results showed enrichment in deregulated genes involved in cell death, cell morphology and cancer as well as an alteration of several signaling circuits involved in energy balance including neuropeptide Y and melanocortin signaling. The expression profile along with the phenotype led us to conclude that anx/anx mice resemble the anorexia-cachexia syndrome typically observed in cancer, infection with human immunodeficiency virus or chronic diseases, rather than starvation, and that anx/anx mice could be considered a good model for the treatment and investigation of this condition.

The quest for a message: budding yeast, a model organism to study the control of pre-mRNA splicing

Removal of introns during pre-mRNA splicing is a critical process in gene expression, and understanding its control at both single-gene and genomic levels is one of the great challenges in Biology. Splicing takes place in a dynamic, large ribonucleoprotein complex known as the spliceosome. Combining Genetics and Biochemistry, Saccharomyces cerevisiae provides insights into its mechanisms, including its regulation by RNA-protein interactions. Recent genome-wide analyses indicate that regulated splicing is broad and biologically relevant even in organisms with a relatively simple intronic structure, such as yeast. Furthermore, the possibility of coordination in splicing regulation at genomic level is becoming clear in this model organism. This should provide a valuable system to approach the complex problem of the role of regulated splicing in genomic expression.

Identification of molecular apocrine breast tumours by microarray analysis.

Previous microarray studies on breast cancer identified multiple tumour classes, of which the most prominent, named luminal and basal, differ in expression of the oestrogen receptor alpha gene (ER). We report here the identification of a group of breast tumours with increased androgen signalling and a 'molecular apocrine' gene expression profile. Tumour samples from 49 patients with large operable or locally advanced breast cancers were tested on Affymetrix U133A gene expression microarrays. Principal components analysis and hierarchical clustering split the tumours into three groups: basal, luminal and a group we call molecular apocrine. All of the molecular apocrine tumours have strong apocrine features on histological examination (P=0.0002). The molecular apocrine group is androgen receptor (AR) positive and contains all of the ER-negative tumours outside the basal group. Kolmogorov-Smirnov testing indicates that oestrogen signalling is most active in the luminal group, and androgen signalling is most active in the molecular apocrine group. ERBB2 amplification is commoner in the molecular apocrine than the other groups. Genes that best split the three groups were identified by Wilcoxon test. Correlation of the average expression profile of these genes in our data with the expression profile of individual tumours in four published breast cancer studies suggest that molecular apocrine tumours represent 8-14% of tumours in these studies. Our data show that it is possible with microarray data to divide mammary tumour cells into three groups based on steroid receptor activity: luminal (ER+ AR+), basal (ER- AR-) and molecular apocrine (ER- AR+).

Predicting cancer involvement of genes from heterogeneous data

Background: Systematic approaches for identifying proteins involved in different types of cancer are needed. Experimental techniques such as microarrays are being used to characterize cancer, but validating their results can be a laborious task. Computational approaches are used to prioritize between genes putatively involved in cancer, usually based on further analyzing experimental data. Results: We implemented a systematic method using the PIANA software that predicts cancer involvement of genes by integrating heterogeneous datasets. Specifically, we produced lists of genes likely to be involved in cancer by relying on: (i) protein-protein interactions; (ii) differential expression data; and (iii) structural and functional properties of cancer genes. The integrative approach that combines multiple sources of data obtained positive predictive values ranging from 23% (on a list of 811 genes) to 73% (on a list of 22 genes), outperforming the use of any of the data sources alone. We analyze a list of 20 cancer gene predictions, finding that most of them have been recently linked to cancer in literature. Conclusion: Our approach to identifying and prioritizing candidate cancer genes can be used to produce lists of genes likely to be involved in cancer. Our results suggest that differential expression studies yielding high numbers of candidate cancer genes can be filtered using protein interaction networks.

Differential Gene Expression between Fall- and Spring-Run Chinook Salmon Assessed by Long Serial Analysis of Gene Expression

Of all Pacific salmonids, Chinook salmon Oncorhynchus tshawytscha display the greatest variability in return times to freshwater. The molecular mechanisms of these differential return times have not been well described. Current methods, such as long serial analysis of gene expression (LongSAGE) and microarrays, allow gene expression to be analyzed for thousands of genes simultaneously. To investigate whether differential gene expression is observed between fall- and spring-run Chinook salmon from California's Central Valley, LongSAGE libraries were constructed. Three libraries containing between 25,512 and 29,372 sequenced tags (21 base pairs/tag) were generated using messenger RNA from the brains of adult Chinook salmon returning in fall and spring and from one ocean-caught Chinook salmon. Tags were annotated to genes using complementary DNA libraries from Atlantic salmon Salmo salar and rainbow trout O. mykiss. Differentially expressed genes, as estimated by differences in the number of sequence tags, were found in all pairwise comparisons of libraries (freshwater versus saltwater = 40 genes; fall versus spring = 11 genes: and spawning versus nonspawning = 51 genes). The gene for ependymin, an extracellular glycoprotein involved in behavioral plasticity in fish, exhibited the most differential expression among the three groupings. Reverse transcription polymerase chain reaction analysis verified the differential expression of ependymin between the fall- and spring-run samples. These LongSAGE libraries, the first reported for Chinook salmon, provide a window of the transcriptional changes during Chinook salmon return migration to freshwater and spawning and increase the amount of expressed sequence data.

Androgen dependence of hirsutism, acne, and alopecia in women: retrospective analysis of 228 patients investigated for hyperandrogenism.

Hirsutism, acne, alopecia, and oligo-amenorrhea are clinical expressions of hyperandrogenism, one of the most frequent endocrine disorders in women of reproductive age. Women referred to our endocrine clinics for skin symptoms of hyperandrogenism underwent a laboratory workup to evaluate hormone measurements and received antiandrogen therapy. We retrospectively analyzed the outcome of 228 consecutive patients investigated over 6 years.Patients with hirsutism had higher levels of androstenedione, dehydroepiandrosterone sulfate (DHEAS), and salivary testosterone; lower levels of sex hormone-binding globulin (SHBG); and a higher prevalence of oligo-amenorrhea than patients with alopecia, while patients with acne showed intermediate values. Hirsutism score correlated positively with androstenedione, DHEAS, and salivary testosterone, and correlated negatively with SHBG; salivary testosterone showed the highest correlation coefficient. Total testosterone was not significantly different among patients with hirsutism, alopecia, or acne, and did not significantly correlate with hirsutism score. Hirsutism and oligo-amenorrhea were the most sensitive symptoms of hyperandrogenism, and no androgenic parameter alone allowed us to identify all cases of hyperandrogenism.Patients of central European origin sought consultation with milder hirsutism scores than patients of southern European origin. There was, however, no difference in the clinical-biological correlation between these groups, arguing against differences in skin sensitivity to androgens.Polycystic ovary syndrome, defined as hyperandrogenism (hirsutism or elevated androgens) and oligo-amenorrhea, was diagnosed in 63 patients (27.6%), an underestimate compared with other reports that include systematic ovarian ultrasound studies. Neither pelvic ultrasound, used in a limited number of cases, nor the luteinizing hormone/follicle-stimulating hormone ratio helped to distinguish patients with polycystic ovary syndrome from the other diagnostic groups. These included hyperandrogenism (hirsutism or elevated androgens) and eumenorrhea (101 patients; 44.3%); normal androgens (acne or alopecia and eumenorrhea) (51 patients; 22.4%); isolated low SHBG (7 patients; 3.1%); nonclassical congenital adrenal hyperplasia (4 patients; 1.8% of total, 4.9% of patients undergoing cosyntropin stimulation tests); and ovarian tumor (2 patients; 0.9%).Ethinylestradiol and high-dose cyproterone acetate treatment lowered the hirsutism score to 53.5% of baseline at 1 year, and was also effective in treating acne and alopecia. The clinical benefit is ascribed to the peripheral antiandrogenic effect of cyproterone acetate as well as the hormone-suppressive effect of this combination. Salivary testosterone showed the most marked proportional decrease of all the androgens under treatment. Cost-effectiveness and tolerance of ethinylestradiol and high-dose cyproterone acetate compared well with other antiandrogenic drug therapies for hirsutism. The less potent therapy with spironolactone only, a peripheral antiandrogen without hormone-suppressive effect, was effective in treating isolated alopecia in patients with normal androgens.

Compositions in life science data

The aim of this talk is to convince the reader that there are a lot of interesting statisticalproblems in presentday life science data analysis which seem ultimately connected withcompositional statistics.Key words: SAGE, cDNA microarrays, (1D-)NMR, virus quasispecies

Mononuclear cell transcriptome response after sustained virgin olive oil consumption in humans: an exploratory nutrigenomics study

Virgin olive oil (VOO) is considered to be one of the main components responsible for the health benefits of the Mediterranean diet, particularly against atherosclerosis where peripheral blood mononuclear cells (PBMNCs) play a crucial role in atherosclerosis development and progression. The objective of this article was to identify the PBMNC genes that respond to VOO consumption in order to ascertain the molecular mechanisms underlying the beneficial action of VOO in the prevention of atherosclerosis. Gene expression profiles of PBMNCs from healthy individuals were examined in pooled RNA samples by microarrays after 3 weeks of moderate and regular consumption of VOO, as the main fat source in a diet controlled for antioxidant content. Gene expression was verified by qPCR. The response to VOO consumption was confirmed for individual samples (n = 10) by qPCR for 10 upregulated genes (ADAM17, ALDH1A1, BIRC1, ERCC5, LIAS, OGT, PPARBP, TNFSF10, USP48, and XRCC5). Their putative role in the molecular mechanisms involved in atherosclerosis development and progression is discussed, focusing on a possible relation with VOO consumption. Our data support the hypothesis that 3 weeks of nutritional intervention with VOO supplementation, at doses common in the Mediterranean diet, can alter the expression of genes related to atherosclerosis development and progression.

VARIABILITY OF PHENOTYPIC, PROTEOMIC AND TRANSCRIPTOMIC EXPRESSION OF STAPHYLOCOCCUS AUREUS SURFACE ADHESINS

Staphylococcus aureus is a highly successful pathogen responsible of a wide variety of diseases, from minor skin infection to life-threatening sepsis or infective endocarditis, as well as food poisoning and toxic shock syndrome. This heterogeneity of infections and the ability of S. aureus to develop antibiotic-resistance to virtually any available drugs reflect its extraordinary capacity to adapt and survive in a great variety of environments. The pathogenesis of S. aureus infection involves a wide range of cell wall-associated adhesins and extracellular toxins that promote host colonization and invasion. In addition, S. aureus is extremely well equipped with regulatory systems that sense environmental conditions and respond by fine tuning the expression of metabolic and virulence determinants. Surface adhesins referred to MSCRAMMs - for Microbial Surface Component Recognizing Adherence Matrix Molecules - mediate binding to the host extracellular matrix or serum components, including fibrinogen, fibronectin, collagen and elastin, and promote tissue colonization and invasion. Major MSCRAMMs include a family of surface-attached proteins covalently bound to the cell wall peptidoglycan via a conserved LPXTG motif. Genomic analyses indicate that S. aureus contain up to 22 LPXTG surface proteins, which could potentially act individually or in synergy to promote infection. In the first part of this study we determined the range of adherence phenotypes to fibrinogen and fibronectin among 30 carriage isolates of S. aureus and compared it to the adherence phenotypes of 30 infective endocarditis and 30 blood culture isolates. Overall there were great variations in in vitro adherence, but no differences were observed between carriage and infection strains. We further determined the relation between in vitro adherence and in vivo infectivity in a rat model of experimental endocarditis, using 4 isolates that displayed either extremely low or high adherence phenotypes. Unexpectedly, no differences were observed between the in vivo infectivity of isolates that were poorly and highly adherent in vitro. We concluded that the natural variability of in vitro adherence to fibrinogen and fibronectin did not correlate with in vivo infectivity, and thus that pathogenic differences between various strains might only be expressed in in vivo conditions, but not in vitro. Therefore, considering the importance of adhesins expression for infection, direct measurement of those adhesins present on the bacterial surface were made by proteomic approach. 5 In the second series of experiments we assessed the physical presence of the LPXTG species at the staphylococcal surface, as measured at various time points during growth in different culture media. S. aureus Newman was grown in either tryptic soy broth (TSB) or in Roswell Park Memorial Institute (RPMI) culture medium, and samples were removed from early exponential growth phase to late stationary phase. Experiments were performed with mutants in the global accessory-gene regulator (agr), surface protein A (Spa) and clumping factor A (ClfA). Peptides of surface proteins were recovered by "trypsin-shaving" of live bacteria, and semi-quantitative proteomic analysis was performed by tandem liquid-chromatography and mass-spectrometry (LC-MS). We also determined in parallel the mRNA expression by microarrays analysis, as well as the phenotypic adherence of the bacteria to fibrinogen in vitro. The surface proteome was highly complex and contained numerous proteins theoretically not belonging to the bacterial envelope, including ribosomal proteins and metabolic enzymes. Sixteen of the 21 known LPXTG species were detected, but were differentially expressed. As expected, 9 known agr-regulated proteins (e.g. including Spa, FnBPA, ClfA, IsdA, IsdB, SasH, SasD, SasG and FmtB) increased up to the late exponential growth phase, and were abrogated in agr-negative mutants. However, only Spa and SasH modified their proteomic and mRNA profiles in parallel in the parent and its agr negative mutant, while all other LPXTG proteins modified their proteomic profiles independently of their mRNA. Moreover, ClfA became highly transcribed and active in in vitro fibrinogen adherence tests during late growth (24h), whereas it remained poorly detected by proteomics. Differential expression was also detected in iron-rich TSB versus iron-poor RPMI. Proteins from the iron-regulated surface determinant (isd) system, including IsdA, IsdB and IsdH were barely expressed in iron-rich TSB, whereas they increased their expression by >10 time in iron-poor RPMI. We conclude that semi-quantitative proteomic analysis of specific protein species is feasible in S. aureus and that proteomic, transcriptomic and adherence phenotypes demonstrated differential profiles in S. aureus. Furthermore, peptide signatures released by trypsin shaving suggested differential protein domain exposures in various environments, which might be relevant for antiadhesins vaccines. A comprehensive understanding of the S. aureus physiology should integrate all these approaches.

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.

A comparison of methods for differential expression analysis of RNA-seq data.

BACKGROUND: Finding genes that are differentially expressed between conditions is an integral part of understanding the molecular basis of phenotypic variation. In the past decades, DNA microarrays have been used extensively to quantify the abundance of mRNA corresponding to different genes, and more recently high-throughput sequencing of cDNA (RNA-seq) has emerged as a powerful competitor. As the cost of sequencing decreases, it is conceivable that the use of RNA-seq for differential expression analysis will increase rapidly. To exploit the possibilities and address the challenges posed by this relatively new type of data, a number of software packages have been developed especially for differential expression analysis of RNA-seq data. RESULTS: We conducted an extensive comparison of eleven methods for differential expression analysis of RNA-seq data. All methods are freely available within the R framework and take as input a matrix of counts, i.e. the number of reads mapping to each genomic feature of interest in each of a number of samples. We evaluate the methods based on both simulated data and real RNA-seq data. CONCLUSIONS: Very small sample sizes, which are still common in RNA-seq experiments, impose problems for all evaluated methods and any results obtained under such conditions should be interpreted with caution. For larger sample sizes, the methods combining a variance-stabilizing transformation with the 'limma' method for differential expression analysis perform well under many different conditions, as does the nonparametric SAMseq method.

Pollution stress on aquatic microbial communities

Thousands of chemical compounds enter the natural environment but many have unknown effects and consequences, in particular at low concentrations. This thesis work contributes to our understanding of pollution effects by using bacteria as test organisms. Bacteria are important for this question because some of them degrade and transform pollutants into less harmful compounds, but secondly because they themselves can be inhibited in their reproduction by exposure to toxic compounds. When inhibitory effects occur this may change the composition of the microbial com¬munity in the long run, leading to altered or diminished ecosystem services by those communities. As a result chemicals of anthropogenic origin may accumulate and per¬sist in the environment, and finally, affect higher organisms as well. In addition to acquiring basic understanding of pollutant effects at low concentrations on bacterial communities an applied goal of this thesis work was to develop bacteria-based tests to screen new organic chemicals for toxicity and biodégradation. In the first part of this work we developed a flow cytometry-based assay on SYT09 plus ethidium-bromide or propidium-iodide stained cells of Pseudomonas ûuorescens exposed or not to a variety of pollutants under oligotrophic growth conditions. Flow cytometry (FC) allows fast and accurate counting of bacterial cells under simul¬taneous assessment of their physiological state, in particular in combination with different fluorescent dyes. Here we employed FC and fluorescent dyes to monitor the effect that pollutants may exert on Pseudomonas ûuorescens SV3. First we designed an oligotrophic growth test, which enabled us to follow population growth at low densities (104 - 10 7 cells per ml) using 0.1 mM sodium acetate as carbon source. Cells in the oligotrophic milieu were then exposed or not to a variety of common pollutants, such as 2-chlorobiphenyl (2CBP), naphthalene (NAH), 4-chlorophenol (4CP), tetradecane (TD), mercury chloride (HgCl2) or benzene, in different dosages. Exposed culture samples were stained with SYT09 (green fluorescent dye binding nucleic acids, generally staining all cells) in combination with propidium iodide (PI) or ethidium bromide (EB), both dyes being membrane integrity indicators. We ob- served that most of the tested compounds decreased population growth in a dosage- dependent manner. SYT09/PI or SYT09/EB staining then revealed that chemical exposure led to arisal of subpopulations of live and injured or dead cells. By modeling population growth on the total cell numbers in population or only the subpopulation of live cells we inferred that even in stressed populations live cells multiply at rates no different to unexposed controls. The net decrease in population growth would thus be a consequence of more and more cells being not able to multiply at all, rather than all cells multiplying at slower rates. In addition, the proportion of injured cells correlated to the compound dosage. We concluded that the oligotrophic test may be useful to asses toxicity of unknown chemicals on a variety of model bacteria. Mul¬tiple tests can be run in parallel and effects are rapidly measured within a period of 8 hours. Interestingly, in the same exposure tests with P. fluorescens SV3 we observed that some chemicals which did not lead to a reduction of net population growth rates did cause measurable effects on live cells. This was mainly observed in cells within the live subpopulation as an increase of the EB fluorescence signal. We showed that SYT09/EB is a more useful combination of dyes than SYT09/PI because PI fluorescence tend to increase only when cells are effectively dead, but not so much in live cells (less then twofold). In contrast, EB geometric mean fluorescence in live cells increased up to eightfold after exposure to toxic compounds. All compounds even at the lowest concentration caused a measurable increase in EB geometric mean fluorescence especially after 2 h incubation time. This effect was found to be transient for cells exposed to 2CBP and 4CP, but chronic for cells incubated with TD and NAH (ultimately leading to cell death). In order to understand the mechanism underlying the observed effects we used known membrane or energy uncouplers. The pattern of EB signal increase in chemical-exposed populations resembled mostly that of EDTA, although EB fluorescence in EDTA-treated or pasteurized cells was even higher than after exposure to the four test chemicals. We conclude that the ability of cells to efflux EB under equilibrium conditions is an appropriate measure for the potential of a chemical to exert toxicity. Since most bacterial species possess efflux systems for EB that all require cellular energy, our test should be more widely relevant to infer toxicity effects of chemical exposure on the physiological status of the bacterial cell. To better understand the effect of toxicant exposure on efflux defense systems, we studied 2-hydroxybiphenyl toxicity to Pseudomonas azeiaica HBP1. We showed that 2-HBP exerts toxicity even to P. azelaica HBP1, but only at concentrations higher than 0.5 mM. Above this concentration transient loss of membrane polarization and integrity occurred, which we conclude from staining of growing cells with fluorescent dyes. Cells finally recover and resume growth on 2HBP. The high resistance of P. azelaica HBP1 to 2-HBP was found to be the result of an efficient MexABOprM- type efflux pump system counteracting passive influx of this compound into the membrane and cellular interior. Mutants with disrupted mexA, mexB and oprM genes did no longer grow on 2-HBP at concentrations above 100 μΜ, whereas below this concentration we found 2-HBP-concentration dependent decrease of growth rate. The MexAB-OprM system in P. azeiaica HBP1 is indeed an efflux pump for ethidium bromide as well. By introducing gfp reporter fusions responsive to intracellular 2- HBP concentrations into HBP1 wild-type or the mutants we demonstrated that 2HBP enters into the cells in a similar way. In contrast, the reporter system in the wild-type cells does not react to 2-HBP at an outside concentration of 2.4 μΜ, whereas in mutant cells it does. This suggests that wild-type cells pump 2-HBP to the outside very effectively preventing accumulation of 2-HBP. 2HBP metabolism, therefore, is not efficient enough to lower the intracellular concentration and prevent toxicity. We conclude that P. azelaica HBP1 resistance to 2-HBP is mainly due to an efficient efflux system and that 2HBP in high concentrations exerts narcotic effects on the bacterial membrane. In the part of this thesis, we investigated the possibilities of bacteria to degrade pollutants at low concentrations (1 mg per L and below). As test components we used 2-hydroxybiphenyl, antibiotics and a variety of fragrances, many of which are known to be difficult to biodegrade. By using accurate counting of low numbers of bacterial cells we could demonstrate that specific growth on these compounds is possible. We demonstrated the accuracy of FC counting at low cell numbers (down to 103 bacterial cells per ml). Then we tested whether bacterial population growth could be specifically monitored at the expense of low substrate concentrations, us¬ing P. azelaica HBP1. A perfect relationship was found between growth rate, yield and 2-HBP concentrations in the range of 0.1 up to 5 mg per L. Mixing P. azelaica within sludge, however, suggested that growth yields in a mixed community can be much lower than in pure culture, perhaps because of loss of metabolic intermediates. We then isolated new strains from activated sludge using 2-HBP or antibiotics (Nal, AMP, SMX) at low concentrations (0.1-1 mg per L) as sole carbon and energy sub¬strate and PAO microdishes. The purified strains were then examined for growth on their respective substrate, which interestingly, showed that all strains can not with¬stand higher than 1 or 10 mg per L concentrations of target substrate. Thus, bacteria must exist that contribute to compound degradation at low pollutant concentrations but are inhibited at higher concentrations. Finally we tested whether specific biomass growth (in number of cells) at the expense of pollutants can also be detected with communities as starting material. Hereto, we focused on a number of fragrance chemicals and measured community biomass increase by flow cytometry cell counting on two distinct starter communities: (i) diluted Lake Geneva water, and dilute activated sludge from a wastewater treatment plant. We observed that most of the test compounds indeed resulted in significant biomass increase in the starter community compared to a no-carbon added control, but activated sludge and lake Geneva water strongly differed (almost mutually ex¬clusive) in their capacity to degrade the test chemicals. In two cases for activated sludge the same type of microbial community developed upon compound exposure, as concluded from transcription fragment length polymorphism analysis on community purified and PCR amplified 16S rRNA gene fragments. To properly test compound biodegradability it is thus important to use starter communities of different origin. We conclude that FC counting can be a valuable tool to screen chemicals for their biodegradability and toxicity. - Des milliers de produits chimiques sont libérés dans l'environnement mais beaucoup ont des effets inconnus, en particulier à basses concentrations. Ce travail de thèse contribue à notre comprehension des effets de la pollution en utilisant des bacteries comme des organismes-tests. Les bacteries sont importantes pour etudier cette ques¬tion car certaines d'entre elles peuvent degrader ou transformer les polluants, mais également parce qu'elles-mmes peuvent tre inhibees dans leur reproduction après avoit ete exposees à ces composes toxiques. Quand des effets inhibiteurs ont lieu, la composition de la communauté microbienne peut tre changee à long terme, ce qui mène à une reduction du service d'ecosystème offert par ces communautés. En consequence, après leur liberation dans l'environnement, les produits chimiques d'origine anthropogenique peuvent soit s'y accumuler et per¬sister, exerant ainsi des effets encore inconnus sur les organismes vivants. En plus d'acquérir des connaissances de base sur les effets des polluants à basses concentra¬tions sur les communautés microbiennes, un but applique de cette thèse était de développer des tests bases sur les bacteries afin d'identifier de nouveau composes pour leur toxicité ou leur biodégradation. Dans la première partie de ce travail, nous avons developpe un test base sur la cytometrie de flux (FC) sur des cellules de Pseudomonas fluorescens colorees par du bromure d'ethidium ou de l'iodure de propidium et exposees ou non à une palette de polluants sous des conditions de croissance oligotrophique. La cytometrie de flux est une technique qui connaît de nombreuses applications dans la microbiologie environ¬nementale. Cela est principalement du au fait qu'elle permet un comptage rapide et precis ainsi que l'évaluation de l'état physiologique, en particulier lorsqu'elle est combinée h des colorations fluorescentes. Ici, nous avons utilise la technique FC et des colorants fluorescents afin de mesurer l'effet que peuvent exercer certains pollu¬ants sur Pseudomonas ûuorescens SV3 . D'abord nous avons conu des tests oligo- trophiques qui nous permettent de suivre la croissance complète de cellules en culture h des densites faibles (104 -10 7 cellules par ml), sur de l'acetate de sodium à 0.1 mM, en presence ou absence de produits chimiques (2-chlorobiphenyl (2CBP), naphthalène (NAH), 4-chlorophenol (4CP), tetradecane (TD), chlorure de mercure(II) (HgCl2)) à différentes concentrations. Afin de montrer le devenir des bacteries tant au niveau de la cellule individuelle que celui de la population globale, après exposition à des series de composes chimiques, nous avons compte les cellules colorees avec du SYT09 (col¬orant fluorescent vert des acides nucléiques pour la discrimination des cellules par rapport au bruit de fond) en combinaison avec l'iodure de propidium (PI) ou le bromure d'ethidium (EB), indicateurs de l'intégrité de la membrane cellulaire avec FC. Nous avons observe que de nombreux composes testes avaient un effet sur la croissance bacterienne, resultant en une baisse du taux de reproduction de la pop¬ulation. En outre, la double coloration que nous avons utilisee dans cette etude SYT09/PI ou SYT09/EB a montre que les produits chimiques testes induisaient une reponse heterogène des cellules dans la population, divisant celle-ci en sous- populations "saine", "endommagee" ou "morte". Les nombres de cellules à partir du comptage et de la proportion de celles "saines" et "endommagees/mortes" ont ensuite ete utilises pour modeliser la croissance de P. ûuorescens SV3 exposee aux produits chimiques. La reduction nette dans la croissance de population est une consequence du fait que de plus en plus de cellules sont incapables de se reproduire, plutt que du fait d'une croissance plus lente de l'ensemble de la population. De plus, la proportion de cellules endommagees est correllee au dosage du compose chimique. Les résultats obtenus nous ont permis de conclure que le test oligotrophique que nous avons developpe peut tre utilise pour l'évaluation de la toxicité de produits chimiques sur différents modèles bacteriens. Des tests multiples peuvent tre lances en parallèle et les effets sont mesures en l'espace de huit heures. Par ailleurs, nous en déduisons que les produits chimiques exercént un effet sur la croissance des cellules de P. ûuorescens SV3, qui est heterogène parmi les cellules dans la population et depend du produit chimique. Il est intéressant de noter que dans les mmes tests d'exposition avec P. ûuorescens SV3, nous avons observe que certains composes qui n'ont pas conduit à une reduction du taux de la croissance nette de la population, ont cause des effets mesurables sur les cellule saines. Ceci a ete essentiellement observe dans la portion "saine" des cellules en tant qu'augmentation du signal de la fluorescence de 1ΈΒ. D'abord nous avons montre que SYT09/EB était une com¬binaison de colorants plus utile que celle de SYT09/PI parce que la fluorescence du PI a tendance à augmenter uniquement lorsque les cellules sont effectivement mortes, et non pas dans les cellules saines (moins de deux fois plus). Par opposi¬tion, la fluorescence moyenne de l'EB dans les cellules saines augmente jusqu'à huit fois plus après exposition aux composes toxiques. Tous les composes, mme aux plus basses concentrations, induisent une augmentation mesurable de la fluorescence moy¬enne de 1ΈΒ, plus particulièrement après deux heures d'incubation. Cet effet s'est revele tre transitoire pour les cellules exposees aux 2CNP et 4CP, mais est chro¬nique pour les cellules incubees avec le TD et le NAH (entranant la mort cellulaire). Afin de comprendre les mécanismes qui sous-tendent les effets observes, nous avons utilise des decoupleurs d'energie ou de membrane. L'augmentation du signal EB dans les populations causee par des produits chimiques ressemblait à celle exerce par le chelateur des ions divalents EDTA. Cependant, les intensités du signal EB des cellules exposees aux produits chimiques testees n'ont jamais atteint les valeurs des cellules traitees avec l'EDTA ou pasteurises. Nous en concluons que le test oli- gotrophique utilisant la coloration (SYT09/)EB des cellules exposees ou non à un produit chimique est utile afin d'evaluer l'effet toxique exerce par les polluants sur la physiologie bacterienne. Afin de mieux comprendre la reaction d'un système de defense par pompe à efflux après exposition à une toxine, nous avons étudié la toxicité du 2-hydroxybiphenyl (2-HBP) sur Pseudomonas azeiaica HBP1. Nous avons montre que le 2-HBP exerce une toxicité mme sur HBP1, mais uniquement à des concentrations supérieures à 0.5 mM. Au-dessus de cette concentration, des pertes transitoires d'intégrité et de polarization membranaire ont lieu, comme cela nous a ete montre par coloration des cellules en croissance. Les cellules sont finalement capables de se rétablir et de reprendre leur croissance sur 2-HBP. La forte resistance de P. azeiaica HBP1 h 2-HBP physiologie bacterienne s'est revele tre le résultat d'un système de pompe h efflux de type MexABOprM qui contre-balance l'influx passif de ce compose h travers la membrane. Nous avons montre, en construisant des mutants avec des insertions dans les gènes mexA, mexB and oprM et des fusions avec le gène rapporteur gfp, que l'altération de n'importe quelle partie du système d'efflux conduisait à accroître l'accumulation de 2-HBP dans la cellule, en comparaison avec la souche sauvage HBP1, provoquant une diminution de la resistance au 2-HBP ainsi qu'une baisse du taux de reproduction des cellules. Des systèmes d'efflux similaires sont répandus chez de nombreuses espèces bactériennes. Ils seraient responsables de la resistance aux produits chimiques tels que les colorants fluorescents (bromure d'ethidium) et des antibiotiques. Nous concluons que la resistance de P. azelaica HBP1 à 2-HBP est principalement due à un système d'efflux efficace et que 2-HBP, à des concentrations elevees, exerce un effet deletère sur la membrane bacterienne. En se basant sur le comptage des cellules avec la FC, nous avons developpe ensuite une methode pour evaluer la biodegradabilite de polluants tels que le 2-HBP ainsi que les antibiotiques (acide nalidixique (Nal), ampicilline (AMP) ou sulfamethoxazole (SMX)) à de faibles concentrations lmg par L et moins), par le suivi de la croissance spécifique sur le compose de cultures microbiennes pures et mixtes. En utilisant un comptage precis de faibles quantités de cellules nous avons pu demontrer que la croissance spécifique sur ces composes est possible. Nous avons pu illustrer la precision du comptage par cytometrie de flux à faible quantité de cellules (jusqu'à 10 3 cellules par ml). Ensuite, nous avons teste s'il était possible de suivre dynamiquement la croissance de la population de cellules sur faibles concentrations de substrats, en utilisant P. azelaica HBP1. Une relation parfaite a ete trouvee entre le taux de croissance, le rendement et les concentrations de 2-HBP (entre 0.1 et 5 mg par L). En mélangeant HBP1 à de la boue active, nous avons pu montrer que le rendement en communauté mixtes pouvait tre bien inférieur qu'en culture pure. Ceci étant peut tre le résultat d'une perte d'intermédiaires métaboliques. Nous avons ensuite isole de nouvelles souches à partir de la boue active en utilisant le 2-HBP ou des antibiotiques (Nal, AMP, SMX) h basses concentrations (0.1-1 mg par L) comme seules sources de carbone et d'energie. En combinaison avec ceci, nous avons également utilise des microplaques PAO. Les souches purifiees ont ensuite ete examinees pour leurs croissances sur leurs substrats respectifs. De faon intéressante, toutes ces souches ont montre qu'elles ne pouvaient pas survivre à des concentrations de substrats supérieures à 1 ou 10 mg par L. Ainsi, il existe des bacteries qui contribuent à la degradation de composes à basses concentrations de polluant mais sont inhibes lorsque ces concentrations deviennent plus hautes. Finalement, nous avons cherche à savoir s'il est possible de detecter une croissance spécifique à une biomasse au depend d'un polluant, en partant d'une communauté microbienne. Ainsi, nous nous sommes concentre sur certains composes et avons mesure l'augmentation de la biomasse d'une communauté grce à la cytometrie de flux. Nous avons compte deux communautés de depart distinctes: (i) une dilution d'eau du Lac Léman, et une dilution de boue active d'une station d'épuration. Nous avons observe que la plupart des composes testes ont entrane une augmentation de la biomasse de depart par rapport au control sans addition de source de carbone. Néanmoins, les échantillons du lac Léman et de la station d'épuration différaient largement (s'excluant mutuellement l'un l'autre) dans leur capacité à degrader les composes chimiques. Dans deux cas provenant de la station d'épuration, le mme type de communauté microbienne s'est developpe après exposition aux composes, comme l'a démontré l'analyse TRFLP sur les fragments d'ARN 16S purifie de la communauté et amplifie par PCR. Afin de tester correctement la biodegradabilite d'un compose, il est donc important d'utiliser des communautés de depart de différentes origines Nous en concluons que le comptage par cytometrie de flux peut tre un outil de grande utilité pour mettre en valeur la biodegradabillite et la toxicité des composes chimiques.