A metadata approach for clinical data management in translational genomics studies in breast cancer

Background

Paragangliomas of the head and neck are highly vascular and usually clinically benign tumors arising in the paraganglia of the autonomic nervous system. A significant number of cases (10–50%) are proven to be familial. Multiple genes encoding subunits of the mitochondrial succinate-dehydrogenase (SDH) complex are associated with hereditary paraganglioma: SDHB, SDHC and SDHD. Furthermore, a hereditary paraganglioma family has been identified with linkage to the PGL2 locus on 11q13. No SDH genes are known to be located in the 11q13 region, and the exact gene defect has not yet been identified in this family.

Methods

We have performed a RNA expression microarray study in sporadic, SDHD- and PGL2-linked head and neck paragangliomas in order to identify potential differences in gene expression leading to tumorigenesis in these genetically defined paraganglioma subgroups. We have focused our analysis on pathways and functional gene-groups that are known to be associated with SDH function and paraganglioma tumorigenesis, i.e. metabolism, hypoxia, and angiogenesis related pathways. We also evaluated gene clusters of interest on chromosome 11 (i.e. the PGL2 locus on 11q13 and the imprinted region 11p15).

Results

We found remarkable similarity in overall gene expression profiles of SDHD -linked, PGL2-linked and sporadic paraganglioma. The supervised analysis on pathways implicated in PGL tumor formation also did not reveal significant differences in gene expression between these paraganglioma subgroups. Moreover, we were not able to detect differences in gene-expression of chromosome 11 regions of interest (i.e. 11q23, 11q13, 11p15).

Conclusion

The similarity in gene-expression profiles suggests that PGL2, like SDHD, is involved in the functionality of the SDH complex, and that tumor formation in these subgroups involves the same pathways as in SDH linked paragangliomas. We were not able to clarify the exact identity of PGL2 on 11q13. The lack of differential gene-expression of chromosome 11 genes might indicate that chromosome 11 loss, as demonstrated in SDHD-linked paragangliomas, is an important feature in the formation of paragangliomas regardless of their genetic background.     

Three CRISPR-Cas immune effector complexes coexist in Pyrococcus furiosus

CRISPR-Cas immune systems function to defend prokaryotes against potentially harmful mobile genetic elements including viruses and plasmids. The multiple CRISPR-Cas systems (Types I, II, and III) each target destruction of foreign nucleic acids via structurally and functionally diverse effector complexes (crRNPs). CRISPR-Cas effector complexes are comprised of CRISPR RNAs (crRNAs) that contain sequences homologous to the invading nucleic acids and Cas proteins specific to each immune system type. We have previously characterized a crRNP in Pyrococcus furiosus (Pfu) that contains Cmr (Type III-B) Cas proteins associated with one of two size classes of crRNAs and cleaves complementary target RNAs. Here, we have isolated and characterized two additional native Pfu crRNPs containing either Csa (Type I-A) or Cst (Type I-G) Cas proteins and distinct profiles of associated crRNAs. For each complex, the Cas proteins were identified by mass spectrometry and immunoblotting and the crRNAs by RNA sequencing and Northern blot analysis. The crRNAs associated with both the Csa and Cst complexes originate from all seven Pfu CRISPR loci and contain identical 5′ ends (8-nt repeat-derived 5′ tag sequences) but heterogeneous 3′ ends (containing variable amounts of downstream repeat sequences). These crRNA forms are distinct from Cmr-associated crRNAs, indicating different 3′ end processing pathways following primary cleavage of common pre-crRNAs. Like other previously characterized Type I CRISPR-Cas effector complexes, we predict that the newly identified Pfu Csa and Cst crRNPs each function to target invading DNA, adding an additional layer of protection beyond that afforded by the previously characterized RNA targeting Cmr complex.     

MicroRNAs from the Planarian Schmidtea mediterranea: a model system for stem cell biology

MicroRNAs (miRNAs) are approximately 22-nt RNA molecules that typically bind to the 3' untranslated regions of target mRNAs and function to either induce mRNA degradation or repress translation. miRNAs have been shown to play important roles in the function of stem cells and cell lineage decisions in a variety of organisms, including humans. Planarians are bilaterally symmetric metazoans that have the unique ability to completely regenerate lost tissues or organs. This regenerative capacity is facilitated by a population of stem cells known as neoblasts. Planarians are therefore an excellent model system for studying many aspects of stem cell biology. Here we report the cloning and initial characterization of 71 miRNAs from the planarian Schmidtea mediterranea. While several of the S. mediterranea miRNAs are members of miRNA families identified in other species, we also identified a number of planarian-specific miRNAs. This work lays the foundation for functional studies aimed at addressing the role of these miRNAs in regeneration, cell lineage decisions, and basic stem cell biology.     

Genomic imprinting absent in Drosophila melanogaster adult females

Genomic imprinting occurs when expression of an allele differs based on the sex of the parent that transmitted the allele. In D. melanogaster, imprinting can occur, but its impact on allelic expression genome-wide is unclear. Here, we search for imprinted genes in D. melanogaster using RNA-seq to compare allele-specific expression between pools of 7- to 10-day-old adult female progeny from reciprocal crosses. We identified 119 genes with allelic expression consistent with imprinting, and these genes showed significant clustering within the genome. Surprisingly, additional analysis of several of these genes showed that either genomic heterogeneity or high levels of intrinsic noise caused imprinting-like allelic expression. Consequently, our data provide no convincing evidence of imprinting for D. melanogaster genes in their native genomic context. Elucidating sources of false-positive signals for imprinting in allele-specific RNA-seq data, as done here, is critical given the growing popularity of this method for identifying imprinted genes.     

A Statistical Model of Protein Sequence Similarity and Function Similarity Reveals Overly-Specific Function Predictions

Background

Predicting protein function from primary sequence is an important open problem in modern biology. Not only are there many thousands of proteins of unknown function, current approaches for predicting function must be improved upon. One problem in particular is overly-specific function predictions which we address here with a new statistical model of the relationship between protein sequence similarity and protein function similarity.

Methodology

Our statistical model is based on sets of proteins with experimentally validated functions and numeric measures of function specificity and function similarity derived from the Gene Ontology. The model predicts the similarity of function between two proteins given their amino acid sequence similarity measured by statistics from the BLAST sequence alignment algorithm. A novel aspect of our model is that it predicts the degree of function similarity shared between two proteins over a continuous range of sequence similarity, facilitating prediction of function with an appropriate level of specificity.

Significance

Our model shows nearly exact function similarity for proteins with high sequence similarity (bit score >244.7, e-value >1e⁻⁶², non-redundant NCBI protein database (NRDB)) and only small likelihood of specific function match for proteins with low sequence similarity (bit score <54.6, e-value <1e⁻⁰⁵, NRDB). For sequence similarity ranges in between our annotation model shows an increasing relationship between function similarity and sequence similarity, but with considerable variability. We applied the model to a large set of proteins of unknown function, and predicted functions for thousands of these proteins ranging from general to very specific. We also applied the model to a data set of proteins with previously assigned, specific functions that were electronically based. We show that, on average, these prior function predictions are more specific (quite possibly overly-specific) compared to predictions from our model that is based on proteins with experimentally determined function.     

Cyclic nucleotide content of tobacco BY-2 cells

The cyclic nucleotide content of cultured tobacco bright yellow-2 (BY-2) cells was determined, after freeze-killing, perchlorate extraction and sequential chromatography, by radioimmunoassay. The identities of the putative cyclic nucleotides, adenosine 3',5'-cyclic monophosphate (cyclic AMP), guanosine 3',5'-cyclic monophosphate (cyclic GMP) and cytidine 3',5'-cyclic monophosphate (cyclic CMP) were unambiguously confirmed by tandem mass spectrometry. The potential of BY-2 cell cultures as a model system for future investigations of cyclic nucleotide function in higher plants is discussed.     

The External RNA Controls Consortium: a progress report

Standard controls and best practice guidelines advance acceptance of data from research, preclinical and clinical laboratories by providing a means for evaluating data quality. The External RNA Controls Consortium (ERCC) is developing commonly agreed-upon and tested controls for use in expression assays, a true industry-wide standard control.