Smed454 dataset: unravelling the transcriptome of Schmidtea mediterranea

Background: Freshwater planarians are an attractive model for regeneration and stem cell research and have become a promising tool in the field of regenerative medicine. With the availability of a sequenced planarian genome, the recent application of modern genetic and high-throughput tools has resulted in revitalized interest in these animals, long known for their amazing regenerative capabilities, which enable them to regrow even a new head after decapitation. However, a detailed description of the planarian transcriptome is essential for future investigation into regenerative processes using planarians as a model system. Results: In order to complement and improve existing gene annotations, we used a 454 pyrosequencing approach to analyze the transcriptome of the planarian species Schmidtea mediterranea Altogether, 598,435 454-sequencing reads, with an average length of 327 bp, were assembled together with the ~10,000 sequences of the S. mediterranea UniGene set using different similarity cutoffs. The assembly was then mapped onto the current genome data. Remarkably, our Smed454 dataset contains more than 3 million novel transcribed nucleotides sequenced for the first time. A descriptive analysis of planarian splice sites was conducted on those Smed454 contigs that mapped univocally to the current genome assembly. Sequence analysis allowed us to identify genes encoding putative proteins with defined structural properties, such as transmembrane domains. Moreover, we annotated the Smed454 dataset using Gene Ontology, and identified putative homologues of several gene families that may play a key role during regeneration, such as neurotransmitter and hormone receptors, homeobox-containing genes, and genes related to eye function. Conclusions: We report the first planarian transcript dataset, Smed454, as an open resource tool that can be accessed via a web interface. Smed454 contains significant novel sequence information about most expressed genes of S. mediterranea. Analysis of the annotated data promises to contribute to identification of gene families poorly characterized at a functional level. The Smed454 transcriptome data will assist in the molecular characterization of S. mediterranea as a model organism, which will be useful to a broad scientific community.

Shaping mechanisms of metal specificity in a family of metazoan Metallothioneins: evolutionary differentiation of Mollusc MTs.

Background: The degree of metal binding specificity in metalloproteins such as metallothioneins (MTs) can be crucial for their functional accuracy. Unlike most other animal species, pulmonate molluscs possess homometallic MT isoforms loaded with Cu+ or Cd2+. They have, so far, been obtained as native metal-MT complexes from snail tissues, where they are involved in the metabolism of the metal ion species bound to the respective isoform. However, it has not as yet been discerned if their specific metal occupation is the result of a rigid control of metal availability, or isoform expression programming in the hosting tissues or of structural differences of the respective peptides determining the coordinative options for the different metal ions. In this study, the Roman snail (Helix pomatia) Cu-loaded and Cd-loaded isoforms (HpCuMT and HpCdMT) were used as model molecules in order t o elucidate the biochemical and evolutionary mechanisms permitting pulmonate MTs to achieve specificity for their cognate metal ion. Results: HpCuMT and HpCdMT were recombinantly synthesized in the presence of Cd2+, Zn2+ or Cu2+ and corresponding metal complexes analysed by electrospray mass spectrometry and circular dichroism (CD) and ultra violet-visible (UV-Vis) spectrophotometry. Both MT isoforms were only able to form unique, homometallic and stable complexes (Cd6-HpCdMT and Cu12-HpCuMT) with their cognate metal ions. Yeast complementation assays demonstrated that the two isoforms assumed metal-specific functions, in agreement with their binding preferences, in heterologous eukaryotic environments. In the snail organism, the functional metal specificity of HpCdMT and HpCuMT was contributed by metal-specific transcription programming and cell-specific expression. Sequence elucidation and phylogenetic analysis of MT isoforms from a number of snail species revealed that they possess an unspecific and two metal-specific MT isoforms, whose metal specificity was achieved exclusively by evolutionary modulation of non-cysteine amino acid positions. Conclusion: The Roman snail HpCdMT and HpCuMT isoforms can thus be regarded as prototypes of isoform families that evolved genuine metal-specificity within pulmonate molluscs. Diversification into these isoforms may have been initiated by gene duplication, followed by speciation and selection towards opposite needs for protecting copper-dominated metabolic pathways from nonessential cadmium. The mechanisms enabling these proteins to be metal-specific could also be relevant for other metalloproteins.

Odorant Receptor (Or) genes: polymorphism and divergence in the D. melanogaster and D. pseudoobscura Lineages

Background: In insects, like in most invertebrates, olfaction is the principal sensory modality, which provides animals with essential information for survival and reproduction. Odorant receptors are involved in this response, mediating interactions between an individual and its environment, as well as between individuals of the same or different species. The adaptive importance of odorant receptors renders them good candidates for having their variation shaped by natural selection. Methodology/Principal Findings: We analyzed nucleotide variation in a subset of eight Or genes located on the 3L chromosomal arm of Drosophila melanogaster in a derived population of this species and also in a population of Drosophila pseudoobscura. Some heterogeneity in the silent polymorphism to divergence ratio was detected in the D. melanogaster/D. simulans comparison, with a single gene (Or67b) contributing ~37% to the test statistic. However, no other signals of a very recent selective event were detected at this gene. In contrast, at the speciation timescale, the MK test uncovered the footprint of positive selection driving the evolution of two of the encoded proteins in both D. melanogaster ¿OR65c and OR67a ¿and D. pseudoobscura ¿OR65b1 and OR67c. Conclusions: The powerful polymorphism/divergence approach provided evidence for adaptive evolution at a rather high proportion of the Or genes studied after relatively recent speciation events. It did not provide, however, clear evidence for very recent selective events in either D. melanogaster or D. pseudoobscura.

Smed454 dataset: unravelling the transcriptome of Schmidtea mediterranea

Background Freshwater planarians are an attractive model for regeneration and stem cell research and have become a promising tool in the field of regenerative medicine. With the availability of a sequenced planarian genome, the recent application of modern genetic and high-throughput tools has resulted in revitalized interest in these animals, long known for their amazing regenerative capabilities, which enable them to regrow even a new head after decapitation. However, a detailed description of the planarian transcriptome is essential for future investigation into regenerative processes using planarians as a model system. Results In order to complement and improve existing gene annotations, we used a 454 pyrosequencing approach to analyze the transcriptome of the planarian species Schmidtea mediterranea Altogether, 598,435 454-sequencing reads, with an average length of 327 bp, were assembled together with the ~10,000 sequences of the S. mediterranea UniGene set using different similarity cutoffs. The assembly was then mapped onto the current genome data. Remarkably, our Smed454 dataset contains more than 3 million novel transcribed nucleotides sequenced for the first time. A descriptive analysis of planarian splice sites was conducted on those Smed454 contigs that mapped univocally to the current genome assembly. Sequence analysis allowed us to identify genes encoding putative proteins with defined structural properties, such as transmembrane domains. Moreover, we annotated the Smed454 dataset using Gene Ontology, and identified putative homologues of several gene families that may play a key role during regeneration, such as neurotransmitter and hormone receptors, homeobox-containing genes, and genes related to eye function. Conclusions We report the first planarian transcript dataset, Smed454, as an open resource tool that can be accessed via a web interface. Smed454 contains significant novel sequence information about most expressed genes of S. mediterranea. Analysis of the annotated data promises to contribute to identification of gene families poorly characterized at a functional level. The Smed454 transcriptome data will assist in the molecular characterization of S. mediterranea as a model organism, which will be useful to a broad scientific community.

The non-LTR retrotransposons in Ciona intestinalis: new insights into the evolution of chordate genomes

Background: Non-long terminal repeat (non-LTR) retrotransposons have contributed to shaping the structure and function of genomes. In silico and experimental approaches have been used to identify the non-LTR elements of the urochordate Ciona intestinalis. Knowledge of the types and abundance of non-LTR elements in urochordates is a key step in understanding their contribution to the structure and function of vertebrate genomes. Results: Consensus elements phylogenetically related to the I, LINE1, LINE2, LOA and R2 elements of the 14 eukaryotic non-LTR clades are described from C. intestinalis. The ascidian elements showed conservation of both the reverse transcriptase coding sequence and the overall structural organization seen in each clade. The apurinic/apyrimidinic endonuclease and nucleic-acid-binding domains encoded upstream of the reverse transcriptase, and the RNase H and the restriction enzyme-like endonuclease motifs encoded downstream of the reverse transcriptase were identified in the corresponding Ciona families. Conclusions: The genome of C. intestinalis harbors representatives of at least five clades of non-LTR retrotransposons. The copy number per haploid genome of each element is low, less than 100, far below the values reported for vertebrate counterparts but within the range for protostomes. Genomic and sequence analysis shows that the ascidian non-LTR elements are unmethylated and flanked by genomic segments with a gene density lower than average for the genome. The analysis provides valuable data for understanding the evolution of early chordate genomes and enlarges the view on the distribution of the non-LTR retrotransposons in eukaryotes.

Hereditary hepatic and systemic amyloidosis caused by a new deletion/insertion mutation in the apolipoprotein AI gene.

We report a Spanish family with autosomal-dominant non-neuropathic hereditary amyloidosis with a unique hepatic presentation and death from liver failure, usually by the sixth decade. The disease is caused by a previously unreported deletion/insertion mutation in exon 4 of the apolipoprotein AI (apoAI) gene encoding loss of residues 60-71 of normal mature apoAI and insertion at that position of two new residues, ValThr. Affected individuals are heterozygous for this mutation and have both normal apoAI and variant molecules bearing one extra positive charge, as predicted from the DNA sequence. The amyloid fibrils are composed exclusively of NH2-terminal fragments of the variant, ending mainly at positions corresponding to residues 83 and 92 in the mature wild-type sequence. Amyloid fibrils derived from the other three known amyloidogenic apoAI variants are also composed of similar NH2-terminal fragments. All known amyloidogenic apoAI variants carry one extra positive charge in this region, suggesting that it may be responsible for their enhanced amyloidogenicity. In addition to causing a new phenotype, this is the first deletion mutation to be described in association with hereditary amyloidosis and it significantly extends the value of the apoAI model for investigation of molecular mechanisms of amyloid fibrillogenesis.

Probe-specific mixed-model approach to detect copy number differences using multiplex ligation-dependent probe amplification (MLPA)

Background: MLPA method is a potentially useful semi-quantitative method to detect copy number alterations in targeted regions. In this paper, we propose a method for the normalization procedure based on a non-linear mixed-model, as well as a new approach for determining the statistical significance of altered probes based on linear mixed-model. This method establishes a threshold by using different tolerance intervals that accommodates the specific random error variability observed in each test sample.Results: Through simulation studies we have shown that our proposed method outperforms two existing methods that are based on simple threshold rules or iterative regression. We have illustrated the method using a controlled MLPA assay in which targeted regions are variable in copy number in individuals suffering from different disorders such as Prader-Willi, DiGeorge or Autism showing the best performace.Conclusion: Using the proposed mixed-model, we are able to determine thresholds to decide whether a region is altered. These threholds are specific for each individual, incorporating experimental variability, resulting in improved sensitivity and specificity as the examples with real data have revealed.

PGC-1α induces mitochondrial and myokine transcriptional programs and lipid droplet and glycogen accumulation in cultured human skeletal muscle cells

The transcriptional coactivator peroxisome proliferator-activated receptor-gamma coactivator 1 alpha (PGC-1α) is a chief activator of mitochondrial and metabolic programs and protects against atrophy in skeletal muscle (skm). Here we tested whether PGC-1α overexpression could restructure the transcriptome and metabolism of primary cultured human skm cells, which display a phenotype that resembles the atrophic phenotype. An oligonucleotide microarray analysis was used to reveal the effects of PGC-1α on the whole transcriptome. Fifty-three different genes showed altered expression in response to PGC-1α: 42 upregulated and 11 downregulated. The main gene ontologies (GO) associated with the upregulated genes were mitochondrial components and processes and this was linked with an increase in COX activity, an indicator of mitochondrial content. Furthermore, PGC-1α enhanced mitochondrial oxidation of palmitate and lactate to CO2, but not glucose oxidation. The other most significantly associated GOs for the upregulated genes were chemotaxis and cytokine activity, and several cytokines, including IL-8/CXCL8, CXCL6, CCL5 and CCL8, were within the most highly induced genes. Indeed, PGC-1α highly increased IL-8 cell protein content. The most upregulated gene was PVALB, which is related to calcium signaling. Potential metabolic regulators of fatty acid and glucose storage were among mainly regulated genes. The mRNA and protein level of FITM1/FIT1, which enhances the formation of lipid droplets, was raised by PGC-1α, while in oleate-incubated cells PGC-1α increased the number of smaller lipid droplets and modestly triglyceride levels, compared to controls. CALM1, the calcium-modulated δ subunit of phosphorylase kinase, was downregulated by PGC-1α, while glycogen phosphorylase was inactivated and glycogen storage was increased by PGC-1α. In conclusion, of the metabolic transcriptome deficiencies of cultured skm cells, PGC-1α rescued the expression of genes encoding mitochondrial proteins and FITM1. Several myokine genes, including IL-8 and CCL5, which are known to be constitutively expressed in human skm cells, were induced by PGC-1α.

Maintenance of A/P body regions in planarians by TCEN49, a putative cystine-knot neurotrophin.

In freshwater planarians, the protein TCEN49 has been linked to the regional specification of the central body region, which includes the pharynx.

Identification of Methylated Genes Associated with Aggressive Clinicopathological Features in Mantle Cell Lymphoma

Background: Mantle cell lymphoma (MCL) is genetically characterized by the t(11;14)(q13;q32) translocation and a high number of secondary chromosomal alterations. The contribution of DNA methylation to MCL lymphomagenesis is not well known. We sought to identify epigenetically silenced genes in these tumours that might have clinical relevance. Methodology/Principal Findings: To identify potential methylated genes in MCL we initially investigated seven MCL cell lines treated with epigenetic drugs and gene expression microarray profiling. The methylation status of selected candidate genes was validated by a quantitative assay and subsequently analyzed in a series of primary MCL (n=38). After pharmacological reversion we identified 252 potentially methylated genes. The methylation analysis of a subset of these genes (n=25) in the MCL cell lines and normal B lymphocytes confirmed that 80% of them were methylated in the cell lines but not in normal lymphocytes. The subsequent analysis in primary MCL identified five genes (SOX9,HOXA9,AHR,NR2F2 ,and ROBO1) frequently methylated in these tumours. The gene methylation events tended to occur in the same primary neoplasms and correlated with higher proliferation, increased number of chromosomal abnormalities, and shorter survival of the patients. Conclusions: We have identified a set of genes whose methylation degree and gene expression levels correlate with aggressive clinicopathological features of MCL. Our findings also suggest that a subset of MCL might show a CpG island methylator phenotype (CIMP) that may influence the behaviour of the tumours.

Characterization of human pancreatic orthotopic tumor xenografts suitable for drug screening

Background Efforts to identify novel therapeutic options for human pancreatic ductal adenocarcinoma (PDAC) have failed to result in a clear improvement in patient survival to date. Pancreatic cancer requires efficient therapies that must be designed and assayed in preclinical models with improved predictor ability. Among the available preclinical models, the orthotopic approach fits with this expectation, but its use is still occasional. Methods An in vivo platform of 11 orthotopic tumor xenografts has been generated by direct implantation of fresh surgical material. In addition, a frozen tumorgraft bank has been created, ensuring future model recovery and tumor tissue availability. Results Tissue microarray studies allow showing a high degree of original histology preservation and maintenance of protein expression patterns through passages. The models display stable growth kinetics and characteristic metastatic behavior. Moreover, the molecular diversity may facilitate the identification of tumor subtypes and comparison of drug responses that complement or confirm information obtained with other preclinical models. Conclusions This panel represents a useful preclinical tool for testing new agents and treatment protocols and for further exploration of the biological basis of drug responses.

Characterization of human pancreatic orthotopic tumor xenografts suitable for drug screening

Background Efforts to identify novel therapeutic options for human pancreatic ductal adenocarcinoma (PDAC) have failed to result in a clear improvement in patient survival to date. Pancreatic cancer requires efficient therapies that must be designed and assayed in preclinical models with improved predictor ability. Among the available preclinical models, the orthotopic approach fits with this expectation, but its use is still occasional. Methods An in vivo platform of 11 orthotopic tumor xenografts has been generated by direct implantation of fresh surgical material. In addition, a frozen tumorgraft bank has been created, ensuring future model recovery and tumor tissue availability. Results Tissue microarray studies allow showing a high degree of original histology preservation and maintenance of protein expression patterns through passages. The models display stable growth kinetics and characteristic metastatic behavior. Moreover, the molecular diversity may facilitate the identification of tumor subtypes and comparison of drug responses that complement or confirm information obtained with other preclinical models. Conclusions This panel represents a useful preclinical tool for testing new agents and treatment protocols and for further exploration of the biological basis of drug responses.

egr-4, a target of EGFR signaling, is required for the formation of the brian primordial and head regeneration in planarians

During the regeneration of freshwater planarians, polarity and patterning programs play essential roles in determining whether a head or a tail regenerates at anterior or posterior-facing wounds. This decision is made very soon after amputation. The pivotal role of the Wnt/β-catenin and Hh signaling pathways in re-establishing anterior-posterior (AP) polarity has been well documented. However, the mechanisms that control the growth and differentiation of the blastema in accordance with its AP identity are less well understood. Previous studies have described a role of Smed-egfr-3, a planarian epidermal growth factor receptor, in blastema growth and differentiation. Here, we identify Smed-egr-4, a zinc-finger transcription factor belonging to the early growth response gene family, as a putative downstream target of Smed-egfr-3. Smed-egr-4 is mainly expressed in the central nervous system and its silencing inhibits anterior regeneration without affecting the regeneration of posterior regions. Single and combinatorial RNA interference to target different elements of the Wnt/β-catenin pathway, together with expression analysis of brain- and anterior-specific markers, revealed that Smed-egr-4: (1) is expressed in two phases - an early Smed-egfr-3-independent phase and a late Smed-egfr-3-dependent phase; (2) is necessary for the differentiation of the brain primordia in the early stages of regeneration; and (3) that it appears to antagonize the activity of the Wnt/β-catenin pathway to allow head regeneration. These results suggest that a conserved EGFR/egr pathway plays an important role in cell differentiation during planarian regeneration and indicate an association between early brain differentiation and the proper progression of head regeneration.

Association of MMP-2 polymorphisms with severe and very severe COPD: a case control study of MMPs-1, 9 and 12 in a European population.

BACKGROUND: Genetic factors play a role in chronic obstructive pulmonary disease (COPD) but are poorly understood. A number of candidate genes have been proposed on the basis of the pathogenesis of COPD. These include the matrix metalloproteinase (MMP) genes which play a role in tissue remodelling and fit in with the protease--antiprotease imbalance theory for the cause of COPD. Previous genetic studies of MMPs in COPD have had inadequate coverage of the genes, and have reported conflicting associations of both single nucleotide polymorphisms (SNPs) and SNP haplotypes, plausibly due to under-powered studies. METHODS: To address these issues we genotyped 26 SNPs, providing comprehensive coverage of reported SNP variation, in MMPs- 1, 9 and 12 from 977 COPD patients and 876 non-diseased smokers of European descent and evaluated their association with disease singly and in haplotype combinations. We used logistic regression to adjust for age, gender, centre and smoking history. RESULTS: Haplotypes of two SNPs in MMP-12 (rs652438 and rs2276109), showed an association with severe/very severe disease, corresponding to GOLD Stages III and IV. CONCLUSIONS: Those with the common A-A haplotype for these two SNPs were at greater risk of developing severe/very severe disease (p = 0.0039) while possession of the minor G variants at either SNP locus had a protective effect (adjusted odds ratio of 0.76; 95% CI 0.61 - 0.94). The A-A haplotype was also associated with significantly lower predicted FEV1 (42.62% versus 44.79%; p = 0.0129). This implicates haplotypes of MMP-12 as modifiers of disease severity.

Vanishing chromosomal inversion clines in Drosophila subobscura from Chile: Is behavioral thermoregulation to blame?

Chromosomal inversion clines paralleling the long-standing ones in native Palearctic populations of Drosophila subobscura evolved swiftly after this species invaded the Americas in the late 1970s and early 1980s. However, the new clines did not consistently continue to converge on the Old World baseline. Our recent survey of Chilean populations of D. subobscura shows that inversion clines have faded or even changed sign with latitude. Here, we investigate the hypothesis that this fading of inversion clines might be due to the Bogert effect, namely, that flies' thermoregulatory behavior has eventually compensated for environmental variation in temperature, thus buffering selection on thermal-related traits. We show that latitudinal divergence in thermal preference (T-p) has evolved in Chile for females, with higher-latitude flies having a lower mean T-p. Plastic responses in T-p also lessen latitudinal thermal variation because flies developed at colder temperatures prefer warmer microclimates. Our results are consistent with the idea that active behavioral thermoregulation might buffer environmental variation and reduce the potential effect of thermal selection on other traits as chromosomal arrangements.