Search tips
Search criteria

Results 1-25 (25)

Clipboard (0)

Select a Filter Below

more »
Year of Publication
more »
1.  ZASP Interacts with the Mechanosensing Protein Ankrd2 and p53 in the Signalling Network of Striated Muscle 
PLoS ONE  2014;9(3):e92259.
ZASP is a cytoskeletal PDZ-LIM protein predominantly expressed in striated muscle. It forms multiprotein complexes and plays a pivotal role in the structural integrity of sarcomeres. Mutations in the ZASP protein are associated with myofibrillar myopathy, left ventricular non-compaction and dilated cardiomyopathy. The ablation of its murine homologue Cypher results in neonatal lethality. ZASP has several alternatively spliced isoforms, in this paper we clarify the nomenclature of its human isoforms as well as their dynamics and expression pattern in striated muscle. Interaction is demonstrated between ZASP and two new binding partners both of which have roles in signalling, regulation of gene expression and muscle differentiation; the mechanosensing protein Ankrd2 and the tumour suppressor protein p53. These proteins and ZASP form a triple complex that appears to facilitate poly-SUMOylation of p53. We also show the importance of two of its functional domains, the ZM-motif and the PDZ domain. The PDZ domain can bind directly to both Ankrd2 and p53 indicating that there is no competition between it and p53 for the same binding site on Ankrd2. However there is competition for this binding site between p53 and a region of the ZASP protein lacking the PDZ domain, but containing the ZM-motif. ZASP is negative regulator of p53 in transactivation experiments with the p53-responsive promoters, MDM2 and BAX. Mutations in the ZASP ZM-motif induce modification in protein turnover. In fact, two mutants, A165V and A171T, were not able to bind Ankrd2 and bound only poorly to alpha-actinin2. This is important since the A165V mutation is responsible for zaspopathy, a well characterized autosomal dominant distal myopathy. Although the mechanism by which this mutant causes disease is still unknown, this is the first indication of how a ZASP disease associated mutant protein differs from that of the wild type ZASP protein.
PMCID: PMC3960238  PMID: 24647531
2.  RNA Sequencing of the Exercise Transcriptome in Equine Athletes 
PLoS ONE  2013;8(12):e83504.
The horse is an optimal model organism for studying the genomic response to exercise-induced stress, due to its natural aptitude for athletic performance and the relative homogeneity of its genetic and environmental backgrounds. Here, we applied RNA-sequencing analysis through the use of SOLiD technology in an experimental framework centered on exercise-induced stress during endurance races in equine athletes. We monitored the transcriptional landscape by comparing gene expression levels between animals at rest and after competition. Overall, we observed a shift from coding to non-coding regions, suggesting that the stress response involves the differential expression of not annotated regions. Notably, we observed significant post-race increases of reads that correspond to repeats, especially the intergenic and intronic L1 and L2 transposable elements. We also observed increased expression of the antisense strands compared to the sense strands in intronic and regulatory regions (1 kb up- and downstream) of the genes, suggesting that antisense transcription could be one of the main mechanisms for transposon regulation in the horse under stress conditions. We identified a large number of transcripts corresponding to intergenic and intronic regions putatively associated with new transcriptional elements. Gene expression and pathway analysis allowed us to identify several biological processes and molecular functions that may be involved with exercise-induced stress. Ontology clustering reflected mechanisms that are already known to be stress activated (e.g., chemokine-type cytokines, Toll-like receptors, and kinases), as well as “nucleic acid binding” and “signal transduction activity” functions. There was also a general and transient decrease in the global rates of protein synthesis, which would be expected after strenuous global stress. In sum, our network analysis points toward the involvement of specific gene clusters in equine exercise-induced stress, including those involved in inflammation, cell signaling, and immune interactions.
PMCID: PMC3877044  PMID: 24391776
3.  The transcriptional landscape of the deep-sea bacterium Photobacterium profundum in both a toxR mutant and its parental strain 
BMC Genomics  2012;13:567.
The deep-sea bacterium Photobacterium profundum is an established model for studying high pressure adaptation. In this paper we analyse the parental strain DB110 and the toxR mutant TW30 by massively parallel cDNA sequencing (RNA-seq). ToxR is a transmembrane DNA-binding protein first discovered in Vibrio cholerae, where it regulates a considerable number of genes involved in environmental adaptation and virulence. In P. profundum the abundance and activity of this protein is influenced by hydrostatic pressure and its role is related to the regulation of genes in a pressure-dependent manner.
To better characterize the ToxR regulon, we compared the expression profiles of wt and toxR strains in response to pressure changes. Our results revealed a complex expression pattern with a group of 22 genes having expression profiles similar to OmpH that is an outer membrane protein transcribed in response to high hydrostatic pressure. Moreover, RNA-seq allowed a deep characterization of the transcriptional landscape that led to the identification of 460 putative small RNA genes and the detection of 298 protein-coding genes previously unknown. We were also able to perform a genome-wide prediction of operon structure, transcription start and termination sites, revealing an unexpected high number of genes (992) with large 5′-UTRs, long enough to harbour cis-regulatory RNA structures, suggesting a correlation between intergenic region size and UTR length.
This work led to a better understanding of high-pressure response in P. profundum. Furthermore, the high-resolution RNA-seq analysis revealed several unexpected features about transcriptional landscape and general mechanisms of controlling bacterial gene expression.
PMCID: PMC3505737  PMID: 23107454
High-pressure adaptation; Deep sea; Extremophile; Transcription; Operon; RNA-seq; UTR; Vibrionaceae; Photobacterium profundum; ToxR
4.  First Survey of the Wheat Chromosome 5A Composition through a Next Generation Sequencing Approach 
PLoS ONE  2011;6(10):e26421.
Wheat is one of the world's most important crops and is characterized by a large polyploid genome. One way to reduce genome complexity is to isolate single chromosomes using flow cytometry. Low coverage DNA sequencing can provide a snapshot of individual chromosomes, allowing a fast characterization of their main features and comparison with other genomes. We used massively parallel 454 pyrosequencing to obtain a 2x coverage of wheat chromosome 5A. The resulting sequence assembly was used to identify TEs, genes and miRNAs, as well as to infer a virtual gene order based on the synteny with other grass genomes. Repetitive elements account for more than 75% of the genome. Gene content was estimated considering non-redundant reads showing at least one match to ESTs or proteins. The results indicate that the coding fraction represents 1.08% and 1.3% of the short and long arm respectively, projecting the number of genes of the whole chromosome to approximately 5,000. 195 candidate miRNA precursors belonging to 16 miRNA families were identified. The 5A genes were used to search for syntenic relationships between grass genomes. The short arm is closely related to Brachypodium chromosome 4, sorghum chromosome 8 and rice chromosome 12; the long arm to regions of Brachypodium chromosomes 4 and 1, sorghum chromosomes 1 and 2 and rice chromosomes 9 and 3. From these similarities it was possible to infer the virtual gene order of 392 (5AS) and 1,480 (5AL) genes of chromosome 5A, which was compared to, and found to be largely congruent with the available physical map of this chromosome.
PMCID: PMC3196578  PMID: 22028874
5.  Multi-Tasking Role of the Mechanosensing Protein Ankrd2 in the Signaling Network of Striated Muscle 
PLoS ONE  2011;6(10):e25519.
Ankrd2 (also known as Arpp) together with Ankrd1/CARP and DARP are members of the MARP mechanosensing proteins that form a complex with titin (N2A)/calpain 3 protease/myopalladin. In muscle, Ankrd2 is located in the I-band of the sarcomere and moves to the nucleus of adjacent myofibers on muscle injury. In myoblasts it is predominantly in the nucleus and on differentiation shifts from the nucleus to the cytoplasm. In agreement with its role as a sensor it interacts both with sarcomeric proteins and transcription factors.
Methodology/Principal Findings
Expression profiling of endogenous Ankrd2 silenced in human myotubes was undertaken to elucidate its role as an intermediary in cell signaling pathways. Silencing Ankrd2 expression altered the expression of genes involved in both intercellular communication (cytokine-cytokine receptor interaction, endocytosis, focal adhesion, tight junction, gap junction and regulation of the actin cytoskeleton) and intracellular communication (calcium, insulin, MAPK, p53, TGF-β and Wnt signaling). The significance of Ankrd2 in cell signaling was strengthened by the fact that we were able to show for the first time that Nkx2.5 and p53 are upstream effectors of the Ankrd2 gene and that Ankrd1/CARP, another MARP member, can modulate the transcriptional ability of MyoD on the Ankrd2 promoter. Another novel finding was the interaction between Ankrd2 and proteins with PDZ and SH3 domains, further supporting its role in signaling. It is noteworthy that we demonstrated that transcription factors PAX6, LHX2, NFIL3 and MECP2, were able to bind both the Ankrd2 protein and its promoter indicating the presence of a regulatory feedback loop mechanism.
In conclusion we demonstrate that Ankrd2 is a potent regulator in muscle cells affecting a multitude of pathways and processes.
PMCID: PMC3189947  PMID: 22016770
6.  Combining ontologies and workflows to design formal protocols for biological laboratories 
Laboratory protocols in life sciences tend to be written in natural language, with negative consequences on repeatability, distribution and automation of scientific experiments. Formalization of knowledge is becoming popular in science. In the case of laboratory protocols two levels of formalization are needed: one for the entities and individuals operations involved in protocols and another one for the procedures, which can be manually or automatically executed. This study aims to combine ontologies and workflows for protocol formalization.
A laboratory domain specific ontology and the COW (Combining Ontologies with Workflows) software tool were developed to formalize workflows built on ontologies. A method was specifically set up to support the design of structured protocols for biological laboratory experiments. The workflows were enhanced with ontological concepts taken from the developed domain specific ontology.
The experimental protocols represented as workflows are saved in two linked files using two standard interchange languages (i.e. XPDL for workflows and OWL for ontologies). A distribution package of COW including installation procedure, ontology and workflow examples, is freely available from
Using COW, a laboratory protocol may be directly defined by wet-lab scientists without writing code, which will keep the resulting protocol's specifications clear and easy to read and maintain.
PMCID: PMC2873243  PMID: 20416048
7.  Large-scale detection and analysis of RNA editing in grape mtDNA by RNA deep-sequencing 
Nucleic Acids Research  2010;38(14):4755-4767.
RNA editing is a widespread post-transcriptional molecular phenomenon that can increase proteomic diversity, by modifying the sequence of completely or partially non-functional primary transcripts, through a variety of mechanistically and evolutionarily unrelated pathways. Editing by base substitution has been investigated in both animals and plants. However, conventional strategies based on directed Sanger sequencing are time-consuming and effectively preclude genome wide identification of RNA editing and assessment of partial and tissue-specific editing sites. In contrast, the high-throughput RNA-Seq approach allows the generation of a comprehensive landscape of RNA editing at the genome level. Short reads from Solexa/Illumina GA and ABI SOLiD platforms have been used to investigate the editing pattern in mitochondria of Vitis vinifera providing significant support for 401 C-to-U conversions in coding regions and an additional 44 modifications in non-coding RNAs. Moreover, 76% of all C-to-U conversions in coding genes represent partial RNA editing events and 28% of them were shown to be significantly tissue specific. Solexa/Illumina and SOLiD platforms showed different characteristics with respect to the specific issue of large-scale editing analysis, and the combined approach presented here reduces the false positive rate of discovery of editing events.
PMCID: PMC2919710  PMID: 20385587
8.  Physical mapping in highly heterozygous genomes: a physical contig map of the Pinot Noir grapevine cultivar 
BMC Genomics  2010;11:204.
Most of the grapevine (Vitis vinifera L.) cultivars grown today are those selected centuries ago, even though grapevine is one of the most important fruit crops in the world. Grapevine has therefore not benefited from the advances in modern plant breeding nor more recently from those in molecular genetics and genomics: genes controlling important agronomic traits are practically unknown. A physical map is essential to positionally clone such genes and instrumental in a genome sequencing project.
We report on the first whole genome physical map of grapevine built using high information content fingerprinting of 49,104 BAC clones from the cultivar Pinot Noir. Pinot Noir, as most grape varieties, is highly heterozygous at the sequence level. This resulted in the two allelic haplotypes sometimes assembling into separate contigs that had to be accommodated in the map framework or in local expansions of contig maps. We performed computer simulations to assess the effects of increasing levels of sequence heterozygosity on BAC fingerprint assembly and showed that the experimental assembly results are in full agreement with the theoretical expectations, given the heterozygosity levels reported for grape. The map is anchored to a dense linkage map consisting of 994 markers. 436 contigs are anchored to the genetic map, covering 342 of the 475 Mb that make up the grape haploid genome.
We have developed a resource that makes it possible to access the grapevine genome, opening the way to a new era both in grape genetics and breeding and in wine making. The effects of heterozygosity on the assembly have been analyzed and characterized by using several complementary approaches which could be easily transferred to the study of other genomes which present the same features.
PMCID: PMC2865496  PMID: 20346114
9.  Correction: High throughput approaches reveal splicing of primary microRNA transcripts and tissue specific expression of mature microRNAs in Vitis vinifera 
BMC Genomics  2010;11:109.
The version of this article published in BMC Genomics 2009, 10:558, contains data in Table 1 which are now known to be unreliable, and an illustration, in Figure 1, of unusual miRNA processing events predicted by these unreliable data. In this full-length correction, new data replace those found to be unreliable, leading to a more straightforward interpretation without altering the principle conclusions of the study. Table 1 and associated methods have been corrected, Figure 1 deleted, supplementary file 1 added, and modifications made to the sections "Deep sequencing of small RNAs from grapevine leaf tissue" and "Microarray analysis of miRNA expression". The editors and authors regret the inconvenience caused to readers by premature publication of the original paper.
MicroRNAs are short (~21 base) single stranded RNAs that, in plants, are generally coded by specific genes and cleaved specifically from hairpin precursors. MicroRNAs are critical for the regulation of multiple developmental, stress related and other physiological processes in plants. The recent annotation of the genome of the grapevine (Vitis vinifera L.) allowed the identification of many putative conserved microRNA precursors, grouped into multiple gene families.
Here we use oligonucleotide arrays to provide the first indication that many of these microRNAs show differential expression patterns between tissues and during the maturation of fruit in the grapevine. Furthermore we demonstrate that whole transcriptome sequencing and deep-sequencing of small RNA fractions can be used both to identify which microRNA precursors are expressed in different tissues and to estimate genomic coordinates and patterns of splicing and alternative splicing for many primary miRNA transcripts.
Our results show that many microRNAs are differentially expressed in different tissues and during fruit maturation in the grapevine. Furthermore, the demonstration that whole transcriptome sequencing can be used to identify candidate splicing events and approximate primary microRNA transcript coordinates represents a significant step towards the large-scale elucidation of mechanisms regulating the expression of microRNAs at the transcriptional and post-transcriptional levels.
PMCID: PMC2831844  PMID: 20152027
10.  High throughput approaches reveal splicing of primary microRNA transcripts and tissue specific expression of mature microRNAs in Vitis vinifera 
BMC Genomics  2009;10:558.
MicroRNAs are short (~21 base) single stranded RNAs that, in plants, are generally coded by specific genes and cleaved specifically from hairpin precursors. MicroRNAs are critical for the regulation of multiple developmental, stress related and other physiological processes in plants. The recent annotation of the genome of the grapevine (Vitis vinifera L.) allowed the identification of many putative conserved microRNA precursors, grouped into multiple gene families.
Here we use oligonucleotide arrays to provide the first indication that many of these microRNAs show differential expression patterns between tissues and during the maturation of fruit in the grapevine. Furthermore we demonstrate that whole transcriptome sequencing and deep-sequencing of small RNA fractions can be used both to identify which microRNA precursors are expressed in different tissues and to estimate genomic coordinates and patterns of splicing and alternative splicing for many primary miRNA transcripts.
Our results show that many microRNAs are differentially expressed in different tissues and during fruit maturation in the grapevine. Furthermore, the demonstration that whole transcriptome sequencing can be used to identify candidate splicing events and approximate primary microRNA transcript coordinates represents a significant step towards the large-scale elucidation of mechanisms regulating the expression of microRNAs at the transcriptional and post-transcriptional levels.
PMCID: PMC2822795  PMID: 19939267
11.  A Class III PDZ Binding Motif in the Myotilin and FATZ Families Binds Enigma Family Proteins: a Common Link for Z-Disc Myopathies▿  
Molecular and Cellular Biology  2008;29(3):822-834.
Interactions between Z-disc proteins regulate muscle functions and disruption of these interactions results in muscle disorders. Mutations in Z-disc components myotilin, ZASP/Cypher, and FATZ-2 (calsarcin-1/myozenin-2) are associated with myopathies. We report here that the myotilin and the FATZ (calsarcin/myozenin) families share high homology at their final C-terminal five amino acids. This C-terminal E[ST][DE][DE]L motif is present almost exclusively in these families and is evolutionary conserved. We show by in vitro and in vivo studies that proteins from the myotilin and FATZ (calsarcin/myozenin) families interact via this novel type of class III PDZ binding motif with the PDZ domains of ZASP/Cypher and other Enigma family members: ALP, CLP-36, and RIL. We show that the interactions can be modulated by phosphorylation. Calmodulin-dependent kinase II phosphorylates the C terminus of FATZ-3 (calsarcin-3/myozenin-3) and myotilin, whereas PKA phosphorylates that of FATZ-1 (calsarcin-2/myozenin-1) and FATZ-2 (calsarcin-1/myozenin-1). This is the first report of a binding motif common to both the myotilin and the FATZ (calsarcin/myozenin) families that is specific for interactions with Enigma family members.
PMCID: PMC2630697  PMID: 19047374
12.  Muscle Research and Gene Ontology: New standards for improved data integration 
The Gene Ontology Project provides structured controlled vocabularies for molecular biology that can be used for the functional annotation of genes and gene products. In a collaboration between the Gene Ontology (GO) Consortium and the muscle biology community, we have made large-scale additions to the GO biological process and cellular component ontologies. The main focus of this ontology development work concerns skeletal muscle, with specific consideration given to the processes of muscle contraction, plasticity, development, and regeneration, and to the sarcomere and membrane-delimited compartments. Our aims were to update the existing structure to reflect current knowledge, and to resolve, in an accommodating manner, the ambiguity in the language used by the community.
The updated muscle terminologies have been incorporated into the GO. There are now 159 new terms covering critical research areas, and 57 existing terms have been improved and reorganized to follow their usage in muscle literature.
The revised GO structure should improve the interpretation of data from high-throughput (e.g. microarray and proteomic) experiments in the area of muscle science and muscle disease. We actively encourage community feedback on, and gene product annotation with these new terms. Please visit the Muscle Community Annotation Wiki .
PMCID: PMC2657163  PMID: 19178689
13.  Annotating genomes with massive-scale RNA sequencing 
Genome Biology  2008;9(12):R175.
A method for de novo genome annotation using high-throughput cDNA sequencing data.
Next generation technologies enable massive-scale cDNA sequencing (so-called RNA-Seq). Mainly because of the difficulty of aligning short reads on exon-exon junctions, no attempts have been made so far to use RNA-Seq for building gene models de novo, that is, in the absence of a set of known genes and/or splicing events. We present G-Mo.R-Se (Gene Modelling using RNA-Seq), an approach aimed at building gene models directly from RNA-Seq and demonstrate its utility on the grapevine genome.
PMCID: PMC2646279  PMID: 19087247
14.  Protein evolution in deep sea bacteria: an analysis of amino acids substitution rates 
Abyssal microorganisms have evolved particular features that enable them to grow in their extreme habitat. Genes belonging to specific functional categories are known to be particularly susceptible to high-pressure; therefore, they should show some evidence of positive selection. To verify this hypothesis we computed the amino acid substitution rates between two deep-sea microorganisms, Photobacterium profundum SS9 and Shewanella benthica KT99, and their respective shallow water relatives.
A statistical analysis of all the orthologs, led to the identification of positive selected (PS) genes, which were then used to evaluate adaptation strategies. We were able to establish "Motility" and "Transport" as two classes significantly enriched with PS genes. The prevalence of transporters led us to analyze variable amino acids (PS sites) by mapping them according to their membrane topology, the results showed a higher frequency of substitutions in the extra-cellular compartment. A similar analysis was performed on soluble proteins, mapping the PS sites on the 3D structure, revealing a prevalence of substitutions on the protein surface. Finally, the presence of some flagellar proteins in the Vibrionaceae PS list confirms the importance of bacterial motility as a SS9 specific adaptation strategy.
The approach presented in this paper is suitable for identifying molecular adaptations to particular environmental conditions. The statistical method takes into account differences in the ratio between non-synonymous to synonymous substitutions, thus allowing the detection of the genes that underwent positive selection. We found that positive selection in deep-sea adapted bacteria targets a wide range of functions, for example solute transport, protein translocation, DNA synthesis and motility. From these data clearly emerges an involvement of the transport and metabolism processes in the deep-sea adaptation strategy of both bathytypes considered, whereas the adaptation of other biological processes seems to be specific to either one or the other. An important role is hypothesized for five PS genes belonging to the transport category that had been previously identified as differentially expressed in microarray experiments. Strikingly, structural mapping of PS sites performed independently on membrane and soluble proteins revealed that residues under positive selection tend to occur in specific protein regions.
PMCID: PMC2600651  PMID: 19014525
15.  Large-Scale Transposon Mutagenesis of Photobacterium profundum SS9 Reveals New Genetic Loci Important for Growth at Low Temperature and High Pressure▿  
Journal of Bacteriology  2007;190(5):1699-1709.
Microorganisms adapted to piezopsychrophilic growth dominate the majority of the biosphere that is at relatively constant low temperatures and high pressures, but the genetic bases for the adaptations are largely unknown. Here we report the use of transposon mutagenesis with the deep-sea bacterium Photobacterium profundum strain SS9 to isolate dozens of mutant strains whose growth is impaired at low temperature and/or whose growth is altered as a function of hydrostatic pressure. In many cases the gene mutation-growth phenotype relationship was verified by complementation analysis. The largest fraction of loci associated with temperature sensitivity were involved in the biosynthesis of the cell envelope, in particular the biosynthesis of extracellular polysaccharide. The largest fraction of loci associated with pressure sensitivity were involved in chromosomal structure and function. Genes for ribosome assembly and function were found to be important for both low-temperature and high-pressure growth. Likewise, both adaptation to temperature and adaptation to pressure were affected by mutations in a number of sensory and regulatory loci, suggesting the importance of signal transduction mechanisms in adaptation to either physical parameter. These analyses were the first global analyses of genes conditionally required for low-temperature or high-pressure growth in a deep-sea microorganism.
PMCID: PMC2258685  PMID: 18156275
16.  CUDA compatible GPU cards as efficient hardware accelerators for Smith-Waterman sequence alignment 
BMC Bioinformatics  2008;9(Suppl 2):S10.
Searching for similarities in protein and DNA databases has become a routine procedure in Molecular Biology. The Smith-Waterman algorithm has been available for more than 25 years. It is based on a dynamic programming approach that explores all the possible alignments between two sequences; as a result it returns the optimal local alignment. Unfortunately, the computational cost is very high, requiring a number of operations proportional to the product of the length of two sequences. Furthermore, the exponential growth of protein and DNA databases makes the Smith-Waterman algorithm unrealistic for searching similarities in large sets of sequences. For these reasons heuristic approaches such as those implemented in FASTA and BLAST tend to be preferred, allowing faster execution times at the cost of reduced sensitivity. The main motivation of our work is to exploit the huge computational power of commonly available graphic cards, to develop high performance solutions for sequence alignment.
In this paper we present what we believe is the fastest solution of the exact Smith-Waterman algorithm running on commodity hardware. It is implemented in the recently released CUDA programming environment by NVidia. CUDA allows direct access to the hardware primitives of the last-generation Graphics Processing Units (GPU) G80. Speeds of more than 3.5 GCUPS (Giga Cell Updates Per Second) are achieved on a workstation running two GeForce 8800 GTX. Exhaustive tests have been done to compare our implementation to SSEARCH and BLAST, running on a 3 GHz Intel Pentium IV processor. Our solution was also compared to a recently published GPU implementation and to a Single Instruction Multiple Data (SIMD) solution. These tests show that our implementation performs from 2 to 30 times faster than any other previous attempt available on commodity hardware.
The results show that graphic cards are now sufficiently advanced to be used as efficient hardware accelerators for sequence alignment. Their performance is better than any alternative available on commodity hardware platforms. The solution presented in this paper allows large scale alignments to be performed at low cost, using the exact Smith-Waterman algorithm instead of the largely adopted heuristic approaches.
PMCID: PMC2323659  PMID: 18387198
17.  Genes involved in TGFβ1-driven epithelial-mesenchymal transition of renal epithelial cells are topologically related in the human interactome map 
BMC Genomics  2007;8:383.
Understanding how mesenchymal cells arise from epithelial cells could have a strong impact in unveiling mechanisms of epithelial cell plasticity underlying kidney regeneration and repair.
In primary human tubular epithelial cells (HUTEC) under different TGFβ1 concentrations we had observed epithelial-to-mesenchymal transition (EMT) but not epithelial-myofibroblast transdifferentiation. We hypothesized that the process triggered by TGFβ1 could be a dedifferentiation event. The purpose of this study is to comprehensively delineate genetic programs associated with TGFβ1-driven EMT in our in vitro model using gene expression profile on large-scale oligonucleotide microarrays.
In HUTEC under TGFβ1 stimulus, 977 genes were found differentially expressed. Thirty genes were identified whose expression depended directly on TGFβ1 concentration. By mapping the differentially expressed genes in the Human Interactome Map using Cytoscape software, we identified a single scale-free network consisting of 2630 interacting proteins and containing 449 differentially expressed proteins. We identified 27 hub proteins in the interactome with more than 29 edges incident on them and encoded by differentially expressed genes. The Gene Ontology analysis showed an excess of up-regulated proteins involved in biological processes, such as "morphogenesis", "cell fate determination" and "regulation of development", and the most up-regulated genes belonged to these categories. In addition, 267 genes were mapped to the KEGG pathways and 14 pathways with more than nine differentially expressed genes were identified. In our model, Smad signaling was not the TGFβ1 action effector; instead, the engagement of RAS/MAPK signaling pathway seems mainly to regulate genes involved in the cell cycle and proliferation/apoptosis.
Our present findings support the hypothesis that context-dependent EMT generated in our model by TGFβ1 might be the outcome of a dedifferentiation. In fact: 1) the principal biological categories involved in the process concern morphogenesis and development; 2) the most up-regulated genes belong to these categories; and, finally, 3) some intracellular pathways are involved, whose engagement during kidney development and nephrogenesis is well known. These long-term effects of TGFβ1 in HUTEC involve genes that are highly interconnected, thereby generating a scale-free network that we named the "TGFβ1 interactome", whose hubs represent proteins that may have a crucial role for HUTEC in response to TGFβ1.
PMCID: PMC2174485  PMID: 17953753
18.  Characterization and Evolution of the Cell Cycle-Associated Mob Domain-Containing Proteins in Eukaryotes 
The MOB family includes a group of cell cycle-associated proteins highly conserved throughout eukaryotes, whose founding members are implicated in mitotic exit and co-ordination of cell cycle progression with cell polarity and morphogenesis. Here we report the characterization and evolution of the MOB domain-containing proteins as inferred from the 43 eukaryotic genomes so far sequenced. We show that genes for Mob-like proteins are present in at least 41 of these genomes, confirming the universal distribution of this protein family and suggesting its prominent biological function. The phylogenetic analysis reveals five distinct MOB domain classes, showing a progressive expansion of this family from unicellular to multicellular organisms, reaching the highest number in mammals. Plant Mob genes appear to have evolved from a single ancestor, most likely after the loss of one or more genes during the early stage of Viridiplantae evolutionary history. Three of the Mob classes are widespread among most of the analyzed organisms. The possible biological and molecular function of Mob proteins and their role in conserved signaling pathways related to cell proliferation, cell death and cell polarity are also presented and critically discussed.
PMCID: PMC2684140  PMID: 19468312
Mob genes; protein structure; phylogenesis; cytokinesis; apoptosis; morphogenesis
19.  A global gene evolution analysis on Vibrionaceae family using phylogenetic profile 
BMC Bioinformatics  2007;8(Suppl 1):S23.
Vibrionaceae represent a significant portion of the cultivable heterotrophic sea bacteria; they strongly affect nutrient cycling and some species are devastating pathogens.
In this work we propose an improved phylogenetic profile analysis on 14 Vibrionaceae genomes, to study the evolution of this family on the basis of gene content.
The phylogenetic profile is based on the observation that genes involved in the same process (e.g. metabolic pathway or structural complex) tend to be concurrently present or absent within different genomes. This allows the prediction of hypothetical functions on the basis of a shared phylogenetic profiles. Moreover this approach is useful to identify putative laterally transferred elements on the basis of their presence on distantly phylogenetically related bacteria.
Vibrionaceae ORFs were aligned against all the available bacterial proteomes. Phylogenetic profile is defined as an array of distances, based on aminoacid substitution matrixes, from single genes to all their orthologues. Final phylogenetic profiles, derived from non-redundant list of all ORFs, was defined as the median of all the profiles belonging to the cluster. The resulting phylogenetic profiles matrix contains gene clusters on the rows and organisms on the columns.
Cluster analysis identified groups of "core genes" with a widespread high similarity across all the organisms and several clusters that contain genes homologous only to a limited set of organisms. On each of these clusters, COG class enrichment has been calculated. The analysis reveals that clusters of core genes have the highest number of enriched classes, while the others are enriched just for few of them like DNA replication, recombination and repair.
We found that mobile elements have heterogeneous profiles not only across the entire set of organisms, but also within Vibrionaceae; this confirms their great influence on bacteria evolution even inside the same family. Furthermore, several hypothetical proteins highly correlate with mobile elements profiles suggesting a possible horizontal transfer mechanism for the evolution of these genes. Finally, we suggested the putative role of some ORFs having an unknown function on the basis of their phylogenetic profile similarity to well characterized genes.
PMCID: PMC1885853  PMID: 17430568
20.  Overview of BITS2005, the Second Annual Meeting of the Italian Bioinformatics Society 
BMC Bioinformatics  2005;6(Suppl 4):S1.
The BITS2005 Conference brought together about 200 Italian scientists working in the field of Bioinformatics, students in Biology, Computer Science and Bioinformatics on March 17–19 2005, in Milan. This Editorial provides a brief overview of the Conference topics and introduces the peer-reviewed manuscripts accepted for publication in this Supplement.
PMCID: PMC1866399
21.  Laterally transferred elements and high pressure adaptation in Photobacterium profundum strains 
BMC Genomics  2005;6:122.
Oceans cover approximately 70% of the Earth's surface with an average depth of 3800 m and a pressure of 38 MPa, thus a large part of the biosphere is occupied by high pressure environments. Piezophilic (pressure-loving) organisms are adapted to deep-sea life and grow optimally at pressures higher than 0.1 MPa. To better understand high pressure adaptation from a genomic point of view three different Photobacterium profundum strains were compared. Using the sequenced piezophile P. profundum strain SS9 as a reference, microarray technology was used to identify the genomic regions missing in two other strains: a pressure adapted strain (named DSJ4) and a pressure-sensitive strain (named 3TCK). Finally, the transcriptome of SS9 grown under different pressure (28 MPa; 45 MPa) and temperature (4°C; 16°C) conditions was analyzed taking into consideration the differentially expressed genes belonging to the flexible gene pool.
These studies indicated the presence of a large flexible gene pool in SS9 characterized by various horizontally acquired elements. This was verified by extensive analysis of GC content, codon usage and genomic signature of the SS9 genome. 171 open reading frames (ORFs) were found to be specifically absent or highly divergent in the piezosensitive strain, but present in the two piezophilic strains. Among these genes, six were found to also be up-regulated by high pressure.
These data provide information on horizontal gene flow in the deep sea, provide additional details of P. profundum genome expression patterns and suggest genes which could perform critical functions for abyssal survival, including perhaps high pressure growth.
PMCID: PMC1239915  PMID: 16162277
22.  Development and production of an oligonucleotide MuscleChip: use for validation of ambiguous ESTs 
BMC Bioinformatics  2002;3:33.
We describe the development, validation, and use of a highly redundant 120,000 oligonucleotide microarray (MuscleChip) containing 4,601 probe sets representing 1,150 known genes expressed in muscle and 2,075 EST clusters from a non-normalized subtracted muscle EST sequencing project (28,074 EST sequences). This set included 369 novel EST clusters showing no match to previously characterized proteins in any database. Each probe set was designed to contain 20–32 25 mer oligonucleotides (10–16 paired perfect match and mismatch probe pairs per gene), with each probe evaluated for hybridization kinetics (Tm) and similarity to other sequences. The 120,000 oligonucleotides were synthesized by photolithography and light-activated chemistry on each microarray.
Hybridization of human muscle cRNAs to this MuscleChip (33 samples) showed a correlation of 0.6 between the number of ESTs sequenced in each cluster and hybridization intensity. Out of 369 novel EST clusters not showing any similarity to previously characterized proteins, we focused on 250 EST clusters that were represented by robust probe sets on the MuscleChip fulfilling all stringent rules. 102 (41%) were found to be consistently "present" by analysis of hybridization to human muscle RNA, of which 40 ESTs (39%) could be genome anchored to potential transcription units in the human genome sequence. 19 ESTs of the 40 ESTs were furthermore computer-predicted as exons by one or more than three gene identification algorithms.
Our analysis found 40 transcriptionally validated, genome-anchored novel EST clusters to be expressed in human muscle. As most of these ESTs were low copy clusters (duplex and triplex) in the original 28,000 EST project, the identification of these as significantly expressed is a robust validation of the transcript units that permits subsequent focus on the novel proteins encoded by these genes.
PMCID: PMC137597  PMID: 12456269
Expression profiling; oligonucleotide microarrays; Affymetrix; muscle; EST
23.  A two-step strategy for constructing specifically self-subtracted cDNA libraries 
Nucleic Acids Research  2002;30(9):e38.
We have developed a new strategy for producing subtracted cDNA libraries that is optimized for connective and epithelial tissues, where a few exceptionally abundant (super-prevalent) RNA species account for a large fraction of the total mRNA mass. Our method consists of a two-step subtraction of the most abundant mRNAs: the first step involves a novel use of oligo-directed RNase H digestion to lower the concentration of tissue-specific, super-prevalent RNAs. In the second step, a highly specific subtraction is achieved through hybridization with probes from a 3′-end ESTs collection. By applying this technique in skeletal muscle, we have constructed subtracted cDNA libraries that are effectively enriched for genes expressed at low levels. We further report on frequent premature termination of transcription in human muscle mitochondria and discuss the importance of this phenomenon in designing subtractive approaches. The tissue-specific collections of cDNA clones generated by our method are particularly well suited for expression profiling.
PMCID: PMC113861  PMID: 11972353
24.  Zasp 
The Journal of Cell Biology  1999;146(2):465-476.
PDZ motifs are modular protein–protein interaction domains, consisting of 80–120 amino acid residues, whose function appears to be the direction of intracellular proteins to multiprotein complexes. In skeletal muscle, there are a few known PDZ-domain proteins, which include neuronal nitric oxide synthase and syntrophin, both of which are components of the dystrophin complex, and actinin-associated LIM protein, which binds to the spectrin-like repeats of α-actinin-2. Here, we report the identification and characterization of a new skeletal muscle protein containing a PDZ domain that binds to the COOH-terminal region of α-actinin-2. This novel 31-kD protein is specifically expressed in heart and skeletal muscle. Using antibodies produced to a fragment of the protein, we can show its location in the sarcomere at the level of the Z-band by immunoelectron microscopy. At least two proteins, 32 kD and 78 kD, can be detected by Western blot analysis of both heart and skeletal muscle, suggesting the existence of alternative forms of the protein. In fact, several forms were found that appear to be the result of alternative splicing. The transcript coding for this Z-band alternatively spliced PDZ motif (ZASP) protein maps on chromosome 10q22.3-10q23.2, near the locus for infantile-onset spinocerebellar ataxia.
PMCID: PMC3206570  PMID: 10427098
skeletal muscle; sarcomeres; muscle proteins; immunoelectron microscopy; alternative splicing
25.  A Septin-based Hierarchy of Proteins Required for Localized Deposition of Chitin in the Saccharomyces cerevisiae Cell Wall  
The Journal of Cell Biology  1997;139(1):75-93.
Just before bud emergence, a Saccharomyces cerevisiae cell forms a ring of chitin in its cell wall; this ring remains at the base of the bud as the bud grows and ultimately forms part of the bud scar marking the division site on the mother cell. The chitin ring seems to be formed largely or entirely by chitin synthase III, one of the three known chitin synthases in S. cerevisiae. The chitin ring does not form normally in temperature-sensitive mutants defective in any of four septins, a family of proteins that are constituents of the “neck filaments” that lie immediately subjacent to the plasma membrane in the mother-bud neck. In addition, a synthetic-lethal interaction was found between cdc12-5, a temperature-sensitive septin mutation, and a mutant allele of CHS4, which encodes an activator of chitin synthase III. Two-hybrid analysis revealed no direct interaction between the septins and Chs4p but identified a novel gene, BNI4, whose product interacts both with Chs4p and Cdc10p and with one of the septins, Cdc10p; this analysis also revealed an interaction between Chs4p and Chs3p, the catalytic subunit of chitin synthase III. Bni4p has no known homologues; it contains a predicted coiled-coil domain, but no other recognizable motifs. Deletion of BNI4 is not lethal, but causes delocalization of chitin deposition and aberrant cellular morphology. Overexpression of Bni4p also causes delocalization of chitin deposition and produces a cellular morphology similar to that of septin mutants. Immunolocalization experiments show that Bni4p localizes to a ring at the mother-bud neck that lies predominantly on the mother-cell side (corresponding to the predominant site of chitin deposition). This localization depends on the septins but not on Chs4p or Chs3p. A GFP-Chs4p fusion protein also localizes to a ring at the mother-bud neck on the mother-cell side. This localization is dependent on the septins, Bni4p, and Chs3p. Chs3p, whose normal localization is similar to that of Chs4p, does not localize properly in bni4, chs4, or septin mutant strains or in strains that accumulate excess Bni4p. In contrast, localization of the septins is essentially normal in bni4, chs4, and chs3 mutant strains and in strains that accumulate excess Bni4p. Taken together, these results suggest that the normal localization of chitin synthase III activity is achieved by assembly of a complex in which Chs3p is linked to the septins via Chs4p and Bni4p.
PMCID: PMC2139831  PMID: 9314530

Results 1-25 (25)