Zhang, Yuji

11 Background: Currently there are no serum markers of response to AA. We performed a proteome based analysis of serum from prostate cancer patients, followed by mapping of candidate markers to networks based on known molecular interactions in Ingenuity Knowledge Base (IKB). Methods: Isobaric mass tags for relative and absolute quanititation (iTRAQ) analyzed by reverse-phase liquid chromatography coupled to tandem mass spectrometry (LC/MS/MS) was performed in serum of 3 non-localized prostate cancer cohorts. The first cohort included 15 paired untreated hormone-sensitive “pre AA” and 3-month “post AA” specimens; the second included 10 “early AA failure” (median failure time: 11 months) and the third included 10 “late AA failure” specimens (median failure time: 95 months). Differentially expressed candidate proteins were identified by comparisons of (i) paired pre/post AA proteomes and (ii) post AA proteome with the combined AA failure cohorts at a False Discovery Rate of 0.2. To facilitate biological interpretation of multiple candidates identified in the comparisons, IPA was used. Candidate markers implicated in IPA networks with statistical significance were pursued in a second, independent patient cohort for association with time to AA failure using univariate Cox regression analysis. Results: Median PSA for pre/post AA first cohort were 3.15 and 0.29 ng/ml. Median PSA for the second and third cohorts were 27.3 and 4.3 ng/ml. Between post AA and AA failure cohorts, 149 proteins were differentially expressed. Between early and late AA failure 98 proteins were differentially expressed; 47 proteins were common in both comparisons. Network enrichment analysis of the 47 proteins by IPA identified four interaction networks (p<0.01), one of which highlighted a role for 17-β-estradiol (E1). Gas chromatography used for measuring 3-month post AA initiation serum E1, estrone (E2) and testosterone levels (N=38) detected high E1, E2 levels associated with longer time to AA failure (P=0.07 for E1; P=0.08 for E2) in the independent cohort. Conclusions: A global proteomic analysis identified post AA initiation serum E1, E2 levels as potential response markers, which needs validation.

Background— Genetic studies of cardiomyopathy and heart failure have limited throughput in mammalian models. Adult zebrafish have been recently pursued as a vertebrate model with higher throughput, but genetic conservation must be tested. Methods and Results— We conducted transcriptome analysis of zebrafish heart and searched for fish homologues of 51 known human dilated cardiomyopathy–associated genes. We also identified genes with high cardiac expression and genes with differential expression between embryonic and adult stages. Among tested genes, 30 had a single zebrafish orthologue, 14 had 2 homologues, and 5 had ≥3 homologues. By analyzing the expression data on the basis of cardiac abundance and enrichment hypotheses, we identified a single zebrafish gene for 14 of 19 multiple-homologue genes and 2 zebrafish homologues of high priority for ACTC1 . Of note, our data suggested vmhc and vmhcl as functional zebrafish orthologues for human genes MYH6 and MYH7 , respectively, which are established molecular markers for cardiac remodeling. Conclusions— Most known genes for human dilated cardiomyopathy have a corresponding zebrafish orthologue, which supports the use of zebrafish as a conserved vertebrate model. Definition of the cardiac transcriptome and fetal gene program will facilitate systems biology studies of dilated cardiomyopathy in zebrafish.

Abstract Ovarian cancer is a clinically and molecularly heterogeneous disease. The driving forces behind this variability are unknown. Here, we report wide variation in the expression of the DNA cytosine deaminase APOBEC3B, with elevated expression in the majority of ovarian cancer cell lines (three SDs above the mean of normal ovarian surface epithelial cells) and high-grade primary ovarian cancers. APOBEC3B is active in the nucleus of several ovarian cancer cell lines and elicits a biochemical preference for deamination of cytosines in 5′-TC dinucleotides. Importantly, examination of whole-genome sequence from 16 ovarian cancers reveals that APOBEC3B expression correlates with total mutation load as well as elevated levels of transversion mutations. In particular, high APOBEC3B expression correlates with C-to-A and C-to-G transversion mutations within 5′-TC dinucleotide motifs in early-stage high-grade serous ovarian cancer genomes, suggesting that APOBEC3B-catalyzed genomic uracil lesions are further processed by downstream DNA "repair" enzymes including error-prone translesion polymerases. These data identify a potential role for APOBEC3B in serous ovarian cancer genomic instability. Cancer Res; 73(24); 7222–31. ©2013 AACR.

Abstract Background: Given the high technical reproducibility and orders of magnitude greater resolution than microarrays, next-generation sequencing of mRNA (RNA-Seq) is quickly becoming the de facto standard for measuring levels of gene expression in biological experiments. Two important questions must be taken into consideration when designing a particular experiment, namely, 1) how deep does one need to sequence? and, 2) how many biological replicates are necessary to observe a significant change in expression? Results: Based on the gene expression distributions from 127 RNA-Seq experiments, we find evidence that 91% ± 4% of all annotated genes are sequenced at a frequency of 0.1 times per million bases mapped, regardless of sample source. Based on this observation, and combining this information with other parameters such as biological variation and technical variation that we empirically estimate from our large datasets, we developed a model to estimate the statistical power needed to identify differentially expressed genes from RNA-Seq experiments. Conclusions: Our results provide a needed reference for ensuring RNA-Seq gene expression studies are conducted with the optimally sample size, power, and sequencing depth. We also make available both R code and an Excel worksheet for investigators to calculate for their own experiments.

AbstractThe biological role of vitamin D receptors (VDR), which are abundantly expressed in developing zebrafish (Danio rerio) as early as 48 h after fertilization, and before the development of a mineralized skeleton and mature intestine and kidney, is unknown. We probed the role of VDR in developing zebrafish biology by examining changes in expression of RNA by whole transcriptome shotgun sequencing (RNA-seq) in fish treated with picomolar concentrations of the VDR ligand and hormonal form of vitamin D3, 1α,25-dihydroxyvitamin D3 [1α,25(OH)2D3)].We observed significant changes in RNAs of transcription factors, leptin, peptide hormones, and RNAs encoding proteins of fatty acid, amino acid, xenobiotic metabolism, receptor-activator of NFκB ligand (RANKL), and calcitonin-like ligand receptor pathways. Early highly restricted, and subsequent massive changes in more than 10% of expressed cellular RNA were observed. At days post fertilization (dpf) 2 [24 h 1α,25(OH)2D3-treatment], only four RNAs were differentially expressed (hormone vs. vehicle). On dpf 4 (72 h treatment), 77 RNAs; on dpf 6 (120 h treatment) 1039 RNAs; and on dpf 7 (144 h treatment), 2407 RNAs were differentially expressed in response to 1α,25(OH)2D3. Fewer RNAs (n = 481) were altered in dpf 7 larvae treated for 24 h with 1α,25(OH)2D3vs. those treated with hormone for 144 h. At dpf 7, in 1α,25(OH)2D3-treated larvae, pharyngeal cartilage was larger and mineralization was greater. Changes in expression of RNAs for transcription factors, peptide hormones, and RNAs encoding proteins integral to fatty acid, amino acid, leptin, calcitonin-like ligand receptor, RANKL, and xenobiotic metabolism pathways, demonstrate heretofore unrecognized mechanisms by which 1α,25(OH)2D3 functions in vivo in developing eukaryotes.

Abstract Summary: Reduced representation bisulfite sequencing (RRBS) is a cost-effective approach for genome-wide methylation pattern profiling. Analyzing RRBS sequencing data is challenging and specialized alignment/mapping programs are needed. Although such programs have been developed, a comprehensive solution that provides researchers with good quality and analyzable data is still lacking. To address this need, we have developed a Streamlined Analysis and Annotation Pipeline for RRBS data (SAAP-RRBS) that integrates read quality assessment/clean-up, alignment, methylation data extraction, annotation, reporting and visualization. This package facilitates a rapid transition from sequencing reads to a fully annotated CpG methylation report to biological interpretation. Availability and implementation: SAAP-RRBS is freely available to non-commercial users at the web site http://ndc.mayo.edu/mayo/research/biostat/stand-alone-packages.cfm. Contact: baheti.saurabh@mayo.edu or sun.zhifu@mayo.edu Supplementary Information: Supplementary data are available at Bioinformatics online.