Genomic and transcriptomic analysis of the streptomycin-dependent Mycobacterium tuberculosis strain 18b.

The ability of Mycobacterium tuberculosis to establish a latent infection (LTBI) in humans confounds the treatment of tuberculosis. Consequently, there is a need to discover new therapeutic agents that can kill M. tuberculosis both during active disease and LTBI. The streptomycin-dependent strain of M. tuberculosis, 18b, provides a useful tool for this purpose since upon removal of streptomycin (STR) it enters a non-replicating state that mimics latency both in vitro and in animal models. The 4.41 Mb genome sequence of M. tuberculosis 18b was determined and this revealed the strain to belong to clade 3 of the ancient ancestral lineage of the Beijing family. STR-dependence was attributable to insertion of a single cytosine in the 530 loop of the 16S rRNA and to a single amino acid insertion in the N-terminal domain of initiation factor 3. RNA-seq was used to understand the genetic programme activated upon STR-withdrawal and hence to gain insight into LTBI. This revealed reconfiguration of gene expression and metabolic pathways showing strong similarities between non-replicating 18b and M. tuberculosis residing within macrophages, and with the core stationary phase and microaerophilic responses. The findings of this investigation confirm the validity of 18b as a model for LTBI, and provide insight into both the evolution of tubercle bacilli and the functioning of the ribosome.

Nuclear factor I-C links platelet-derived growth factor and transforming growth factor beta1 signaling to skin wound healing progression.

Transforming growth factor beta (TGF-beta) and platelet-derived growth factor A (PDGFAlpha) play a central role in tissue morphogenesis and repair, but their interplay remain poorly understood. The nuclear factor I C (NFI-C) transcription factor has been implicated in TGF-beta signaling, extracellular matrix deposition, and skin appendage pathologies, but a potential role in skin morphogenesis or healing had not been assessed. To evaluate this possibility, we performed a global gene expression analysis in NFI-C(-/-) and wild-type embryonic primary murine fibroblasts. This indicated that NFI-C acts mostly to repress gene expression in response to TGF-beta1. Misregulated genes were prominently overrepresented by regulators of connective tissue inflammation and repair. In vivo skin healing revealed a faster inflammatory stage and wound closure in NFI-C(-/-) mice. Expression of PDGFA and PDGF-receptor alpha were increased in wounds of NFI-C(-/-) mice, explaining the early recruitment of macrophages and fibroblasts. Differentiation of fibroblasts to contractile myofibroblasts was also elevated, providing a rationale for faster wound closure. Taken together with the role of TGF-beta in myofibroblast differentiation, our results imply a central role of NFI-C in the interplay of the two signaling pathways and in regulation of the progression of tissue regeneration.

Side effects of genome structural changes.

The first extensive catalog of structural human variation was recently released. It showed that large stretches of genomic DNA that vary considerably in copy number were extremely abundant. Thus it is conceivable that they play a major role in functional variation. Consistently, genomic insertions and deletions were shown to contribute to phenotypic differences by modifying not only the expression levels of genes within the aneuploid segments but also of normal copy-number neighboring genes. In this report, we review the possible mechanisms behind this latter effect.

Modeling sequencing errors by combining Hidden Markov models.

Among the largest resources for biological sequence data is the large amount of expressed sequence tags (ESTs) available in public and proprietary databases. ESTs provide information on transcripts but for technical reasons they often contain sequencing errors. Therefore, when analyzing EST sequences computationally, such errors must be taken into account. Earlier attempts to model error prone coding regions have shown good performance in detecting and predicting these while correcting sequencing errors using codon usage frequencies. In the research presented here, we improve the detection of translation start and stop sites by integrating a more complex mRNA model with codon usage bias based error correction into one hidden Markov model (HMM), thus generalizing this error correction approach to more complex HMMs. We show that our method maintains the performance in detecting coding sequences.

Results of a collaborative study of the EDNAP group regarding mitochondrial DNA heteroplasmy and segregation in hair shafts.

A collaborative exercise was carried out by the European DNA Profiling Group (EDNAP) in order to evaluate the distribution of mitochondrial DNA (mtDNA) heteroplasmy amongst the hairs of an individual who displays point heteroplasmy in blood and buccal cells. A second aim of the exercise was to study reproducibility of mtDNA sequencing of hairs between laboratories using differing chemistries, further to the first mtDNA reproducibility study carried out by the EDNAP group. Laboratories were asked to type 2 sections from each of 10 hairs, such that each hair was typed by at least two laboratories. Ten laboratories participated in the study, and a total of 55 hairs were typed. The results showed that the C/T point heteroplasmy observed in blood and buccal cells at position 16234 segregated differentially between hairs, such that some hairs showed only C, others only T and the remainder, C/T heteroplasmy at varying ratios. Additionally, differential segregation of heteroplasmic variants was confirmed in independent extracts at positions 16093 and the poly(C) tract at 302-309, whilst a complete A-G transition was confirmed at position 16129 in one hair. Heteroplasmy was observed at position 16195 on both strands of a single extract from one hair segment, but was not observed in the extracts from any other segment of the same hair. Similarly, heteroplasmy at position 16304 was observed on both strands of a single extract from one hair. Additional variants at positions 73, 249 and the HVII poly(C) region were reported by one laboratory; as these were not confirmed in independent extracts, the possibility of contamination cannot be excluded. Additionally, the electrophoresis and detection equipment used by this laboratory was different to those of the other laboratories, and the discrepancies at position 249 and the HVII poly(C) region appear to be due to reading errors that may be associated with this technology. The results, and their implications for forensic mtDNA typing, are discussed in the light of the biology of hair formation.

Biodiversity of amoebae and amoebae-resisting bacteria in a drinking water treatment plant

The complex ecology of free-living amoebae (FLA) and their role in spreading pathogenic microorganisms through water systems have recently raised considerable interest. In this study, we investigated the presence of FLA and amoebae-resisting bacteria (ARB) at various stages of a drinking water plant fed with river water. We isolated various amoebal species from the river and from several points within the plant, mostly at early steps of water treatment. Echinamoeba- and Hartmannella-related amoebae were mainly recovered in the drinking water plant whereas Acanthamoeba- and Naegleria-related amoebae were recovered from the river water and the sand filtration units. Some FLA isolates were recovered immediately after the ozonation step, thus suggesting resistance of these microorganisms to this disinfection procedure. A bacterial isolate related to Mycobacterium mucogenicum was recovered from an Echinamoeba-related amoeba isolated from ozone-treated water. Various other ARB were recovered using co-culture with axenic Acanthamoeba castellanii, including mycobacteria, legionella, Chlamydia-like organisms and various proteobacteria. Noteworthy, a new Parachlamydia acanthamoebae strain was recovered from river water and from granular activated carbon (GAC) biofilm. As amoebae mainly multiply in sand and GAC filters, optimization of filter backwash procedures probably offers a possibility to better control these protists and the risk associated with their intracellular hosts

High throughput sequencing and proteomics to identify immunogenic proteins of a new pathogen: the dirty genome approach.

With the availability of new generation sequencing technologies, bacterial genome projects have undergone a major boost. Still, chromosome completion needs a costly and time-consuming gap closure, especially when containing highly repetitive elements. However, incomplete genome data may be sufficiently informative to derive the pursued information. For emerging pathogens, i.e. newly identified pathogens, lack of release of genome data during gap closure stage is clearly medically counterproductive. We thus investigated the feasibility of a dirty genome approach, i.e. the release of unfinished genome sequences to develop serological diagnostic tools. We showed that almost the whole genome sequence of the emerging pathogen Parachlamydia acanthamoebae was retrieved even with relatively short reads from Genome Sequencer 20 and Solexa. The bacterial proteome was analyzed to select immunogenic proteins, which were then expressed and used to elaborate the first steps of an ELISA. This work constitutes the proof of principle for a dirty genome approach, i.e. the use of unfinished genome sequences of pathogenic bacteria, coupled with proteomics to rapidly identify new immunogenic proteins useful to develop in the future specific diagnostic tests such as ELISA, immunohistochemistry and direct antigen detection. Although applied here to an emerging pathogen, this combined dirty genome sequencing/proteomic approach may be used for any pathogen for which better diagnostics are needed. These genome sequences may also be very useful to develop DNA based diagnostic tests. All these diagnostic tools will allow further evaluations of the pathogenic potential of this obligate intracellular bacterium.

Blood pressure loci identified with a gene-centric array.

Raised blood pressure (BP) is a major risk factor for cardiovascular disease. Previous studies have identified 47 distinct genetic variants robustly associated with BP, but collectively these explain only a few percent of the heritability for BP phenotypes. To find additional BP loci, we used a bespoke gene-centric array to genotype an independent discovery sample of 25,118 individuals that combined hypertensive case-control and general population samples. We followed up four SNPs associated with BP at our p < 8.56 × 10(-7) study-specific significance threshold and six suggestively associated SNPs in a further 59,349 individuals. We identified and replicated a SNP at LSP1/TNNT3, a SNP at MTHFR-NPPB independent (r(2) = 0.33) of previous reports, and replicated SNPs at AGT and ATP2B1 reported previously. An analysis of combined discovery and follow-up data identified SNPs significantly associated with BP at p < 8.56 × 10(-7) at four further loci (NPR3, HFE, NOS3, and SOX6). The high number of discoveries made with modest genotyping effort can be attributed to using a large-scale yet targeted genotyping array and to the development of a weighting scheme that maximized power when meta-analyzing results from samples ascertained with extreme phenotypes, in combination with results from nonascertained or population samples. Chromatin immunoprecipitation and transcript expression data highlight potential gene regulatory mechanisms at the MTHFR and NOS3 loci. These results provide candidates for further study to help dissect mechanisms affecting BP and highlight the utility of studying SNPs and samples that are independent of those studied previously even when the sample size is smaller than that in previous studies.

Skin fibroblast model to study an impaired glutathione synthesis: consequences of a genetic polymorphism on the proteome.

An impaired glutathione (GSH) synthesis was observed in several multifactorial diseases, including schizophrenia and myocardial infarction. Genetic studies revealed an association between schizophrenia and a GAG trinucleotide repeat (TNR) polymorphism in the catalytic subunit (GCLC) of the glutamate cysteine ligase (GCL). Disease-associated genotypes of this polymorphism correlated with a decrease in GCLC protein expression, GCL activity and GSH content. To clarify consequences of a decreased GCL activity at the proteome level, three schizophrenia patients and three controls have been selected based on the GCLC GAG TNR polymorphism. Fibroblast cultures were obtained by skin biopsy and were challenged with tert-butylhydroquinone (t-BHQ), a substance known to induce oxidative stress. Proteome changes were analyzed by two dimensional gel electrophoresis (2-DE) and results revealed 10 spots that were upregulated in patients following t-BHQ treatment, but not in controls. Nine corresponding proteins could be identified by MALDI mass spectrometry and these proteins are involved in various cellular functions, including energy metabolism, oxidative stress response, and cytoskeletal reorganization. In conclusion, skin fibroblasts of subjects with an impaired GSH synthesis showed an altered proteome reaction in response to oxidative stress. Furthermore, the study corroborates the use of fibroblasts as an additional mean to study vulnerability factors of psychiatric diseases.

Recurrent deletions and reciprocal duplications of 10q11.21q11.23 including CHAT and SLC18A3 are likely mediated by complex low-copy repeats.

We report 24 unrelated individuals with deletions and 17 additional cases with duplications at 10q11.21q21.1 identified by chromosomal microarray analysis. The rearrangements range in size from 0.3 to 12 Mb. Nineteen of the deletions and eight duplications are flanked by large, directly oriented segmental duplications of >98% sequence identity, suggesting that nonallelic homologous recombination (NAHR) caused these genomic rearrangements. Nine individuals with deletions and five with duplications have additional copy number changes. Detailed clinical evaluation of 20 patients with deletions revealed variable clinical features, with developmental delay (DD) and/or intellectual disability (ID) as the only features common to a majority of individuals. We suggest that some of the other features present in more than one patient with deletion, including hypotonia, sleep apnea, chronic constipation, gastroesophageal and vesicoureteral refluxes, epilepsy, ataxia, dysphagia, nystagmus, and ptosis may result from deletion of the CHAT gene, encoding choline acetyltransferase, and the SLC18A3 gene, mapping in the first intron of CHAT and encoding vesicular acetylcholine transporter. The phenotypic diversity and presence of the deletion in apparently normal carrier parents suggest that subjects carrying 10q11.21q11.23 deletions may exhibit variable phenotypic expressivity and incomplete penetrance influenced by additional genetic and nongenetic modifiers.

Uncoupling phosphate deficiency from its major effects on growth and transcriptome via PHO1 expression in Arabidopsis.

Inorganic phosphate (Pi) is one of the most limiting nutrients for plant growth in both natural and agricultural contexts. Pi-deficiency leads to a strong decrease in shoot growth, and triggers extensive changes at the developmental, biochemical and gene expression levels that are presumably aimed at improving the acquisition of this nutrient and sustaining growth. The Arabidopsis thaliana PHO1 gene has previously been shown to participate in the transport of Pi from roots to shoots, and the null pho1 mutant has all the hallmarks associated with shoot Pi deficiency. We show here that A. thaliana plants with a reduced expression of PHO1 in roots have shoot growth similar to Pi-sufficient plants, despite leaves being strongly Pi deficient. Furthermore, the gene expression profile normally triggered by Pi deficiency is suppressed in plants with low PHO1 expression. At comparable levels of shoot Pi supply, the wild type reduces shoot growth but maintains adequate shoot vacuolar Pi content, whereas the PHO1 underexpressor maintains maximal growth with strongly depleted Pi reserves. Expression of the Oryza sativa (rice) PHO1 ortholog in the pho1 null mutant also leads to plants that maintain normal growth and suppression of the Pi-deficiency response, despite the low shoot Pi. These data show that it is possible to unlink low shoot Pi content with the responses normally associated with Pi deficiency through the modulation of PHO1 expression or activity. These data also show that reduced shoot growth is not a direct consequence of Pi deficiency, but is more likely to be a result of extensive gene expression reprogramming triggered by Pi deficiency.

The G0/G1 switch gene 2 is a novel PPAR target gene.

PPARs (peroxisome-proliferator-activated receptors) alpha, beta/delta and gamma are a group of transcription factors that are involved in numerous processes, including lipid metabolism and adipogenesis. By comparing liver mRNAs of wild-type and PPARalpha-null mice using microarrays, a novel putative target gene of PPARalpha, G0S2 (G0/G1 switch gene 2), was identified. Hepatic expression of G0S2 was up-regulated by fasting and by the PPARalpha agonist Wy14643 in a PPARalpha-dependent manner. Surprisingly, the G0S2 mRNA level was highest in brown and white adipose tissue and was greatly up-regulated during mouse 3T3-L1 and human SGBS (Simpson-Golabi-Behmel syndrome) adipogenesis. Transactivation, gel shift and chromatin immunoprecipitation assays indicated that G0S2 is a direct PPARgamma and probable PPARalpha target gene with a functional PPRE (PPAR-responsive element) in its promoter. Up-regulation of G0S2 mRNA seemed to be specific for adipogenesis, and was not observed during osteogenesis or myogenesis. In 3T3-L1 fibroblasts, expression of G0S2 was associated with growth arrest, which is required for 3T3-L1 adipogenesis. Together, these data indicate that G0S2 is a novel target gene of PPARs that may be involved in adipocyte differentiation.