Inherited Germline Variants in Urinary Tract Cancer: A Multicenter Whole-Exome Sequencing Analysis and Correlation With Clinical Features and Tumor Genomics.

PURPOSE: This study investigates a real-world multicenter cohort of patients with urinary tract cancer (UTC), with primary disease sites including the bladder, urethra, and upper tract, who enrolled for research molecular testing of their germline and tumor. The purpose of this study was to evaluate factors that could affect the likelihood of identifying a clinically actionable germline pathogenic variant (PV). METHODS: Patients with UTC were identified from 10 cancer institutes of the Oncology Research Information Exchange Network consortium. The data set comprised abstracted clinical data with germline and tumor genomic data, and comparative analyses were conducted. RESULTS: Clinically actionable germline PVs in cancer predisposition genes were identified in 16 (4.5%) of 354 patients. A higher proportion of patients with the urethra and the upper tract as the primary sites of disease had PVs with a prevalence of 11% (5/45), compared with only 3.6% (11/308) in those with the bladder as the primary site of disease (P = .04). There were no significant differences in markers of genomic instability (such as tumor mutational burden, microsatellite instability [MSI], and loss of heterozygosity, copy number, and chromosomal instability) between those with PVs and those without (P > .05). Of the PVs identified, 10 (62%) were in homologous recombination repair (HRR) genes, three (19%) in mismatch repair (MMR) genes, and three (19%) in genes associated with other pathways. CONCLUSION: Tissue-based assessment of genomic instability, such as MSI, does not reliably indicate germline PV. A comprehensive clinical germline testing approach that includes HRR genes in addition to MMR genes is likely to yield PVs in approximately one of 10 patients with nonbladder primary disease sites such as the upper tract and the urethra.

Zinc-alpha-2-glycoprotein Secreted by Triple-Negative Breast Cancer Promotes Peritumoral Fibrosis.

Obesity is a modifiable predisposition factor for postmenopausal breast cancer. This suggests a localized, reciprocal interaction between breast cancer cells and the surrounding mammary white adipose tissue. To investigate how breast cancer cells alter the composition and function of adipose tissue, we screened the secretomes of 10 human breast cancer cell lines for the ability to modulate the differentiation of adipocyte stem and progenitor cells. The screen identified an adipogenic modulator, zinc-alpha-2-glycoprotein (ZAG/AZGP1) that is secreted by triple-negative breast cancer (TNBC) cells. TNBC-secreted ZAG inhibits adipogenesis and instead induces the expression of fibrotic genes. Accordingly, depletion of ZAG in TNBC cells attenuates fibrosis in white adipose tissue and inhibits tumor growth. Further, high expression of ZAG is linked to poor prognosis in patients with TNBC but not in patients with other clinical subtypes of breast cancer. Our findings suggest a role of TNBC-secreted ZAG in promoting the transdifferentiation of adipocyte stem and progenitor cells into cancer-associated fibroblasts to support tumorigenesis. SIGNIFICANCE: Functional screening of breast cancer secretomes revealed that triple-negative breast cancer promotes fibrosis in the adipose tissue microenvironment by secreting zinc-alpha-2-glycoprotein and promoting the transdifferentiation of adipocyte stem cells into myofibroblasts.

Zinc Alpha-2-Glycoprotein (ZAG/AZGP1) secreted by triple-negative breast cancer promotes tumor microenvironment fibrosis.

Obesity is a predisposition factor for breast cancer, suggesting a localized, reciprocal interaction between breast cancer cells and the surrounding mammary white adipose tissue. To investigate how breast cancer cells alter the composition and function of adipose tissue, we screened the secretomes of ten human breast cancer cell lines for the ability to modulate the differentiation of adipocyte stem and progenitor cells (ASPC). The screen identified a key adipogenic modulator, Zinc Alpha-2-Glycoprotein (ZAG/AZGP1), secreted by triple-negative breast cancer (TNBC) cells. TNBC-secreted ZAG inhibits adipogenesis and instead induces the expression of fibrotic genes. Accordingly, depletion of ZAG in TNBC cells attenuates fibrosis in white adipose tissue and inhibits tumor growth. Further, high expression of ZAG in TNBC patients, but not other clinical subtypes of breast cancer, is linked to poor prognosis. Our findings suggest a role of TNBC-secreted ZAG in promoting the transdifferentiation of ASPCs into cancer-associated fibroblasts to support tumorigenesis.

Hypermedia-based software architecture enables Test-Driven Development.

Objectives: Using agile software development practices, develop and evaluate an architecture and implementation for reliable and user-friendly self-service management of bioinformatic data stored in the cloud. Materials and methods: Comprehensive Oncology Research Environment (CORE) Browser is a new open-source web application for cancer researchers to manage sequencing data organized in a flexible format in Amazon Simple Storage Service (S3) buckets. It has a microservices- and hypermedia-based architecture, which we integrated with Test-Driven Development (TDD), the iterative writing of computable specifications for how software should work prior to development. Relying on repeating patterns found in hypermedia-based architectures, we hypothesized that hypermedia would permit developing test "templates" that can be parameterized and executed for each microservice, maximizing code coverage while minimizing effort. Results: After one-and-a-half years of development, the CORE Browser backend had 121 test templates and 875 custom tests that were parameterized and executed 3031 times, providing 78% code coverage. Discussion: Architecting to permit test reuse through a hypermedia approach was a key success factor for our testing efforts. CORE Browser's application of hypermedia and TDD illustrates one way to integrate software engineering methods into data-intensive networked applications. Separating bioinformatic data management from analysis distinguishes this platform from others in bioinformatics and may provide stable data management while permitting analysis methods to advance more rapidly. Conclusion: Software engineering practices are underutilized in informatics. Similar informatics projects will more likely succeed through application of good architecture and automated testing. Our approach is broadly applicable to data management tools involving cloud data storage.

Replicative Instability Drives Cancer Progression.

In the past decade, defective DNA repair has been increasingly linked with cancer progression. Human tumors with markers of defective DNA repair and increased replication stress exhibit genomic instability and poor survival rates across tumor types. Seminal studies have demonstrated that genomic instability develops following inactivation of BRCA1, BRCA2, or BRCA-related genes. However, it is recognized that many tumors exhibit genomic instability but lack BRCA inactivation. We sought to identify a pan-cancer mechanism that underpins genomic instability and cancer progression in BRCA-wildtype tumors. Methods: Using multi-omics data from two independent consortia, we analyzed data from dozens of tumor types to identify patient cohorts characterized by poor outcomes, genomic instability, and wildtype BRCA genes. We developed several novel metrics to identify the genetic underpinnings of genomic instability in tumors with wildtype BRCA. Associated clinical data was mined to analyze patient responses to standard of care therapies and potential differences in metastatic dissemination. Results: Systematic analysis of the DNA repair landscape revealed that defective single-strand break repair, translesion synthesis, and non-homologous end-joining effectors drive genomic instability in tumors with wildtype BRCA and BRCA-related genes. Importantly, we find that loss of these effectors promotes replication stress, therapy resistance, and increased primary carcinoma to brain metastasis. Conclusions: Our results have defined a new pan-cancer class of tumors characterized by replicative instability (RIN). RIN is defined by the accumulation of intra-chromosomal, gene-level gain and loss events at replication stress sensitive (RSS) genome sites. We find that RIN accelerates cancer progression by driving copy number alterations and transcriptional program rewiring that promote tumor evolution. Clinically, we find that RIN drives therapy resistance and distant metastases across multiple tumor types.

Detection of circulating tumor DNA without a tumor-informed search using next-generation sequencing is a prognostic biomarker in pancreatic ductal adenocarcinoma.

The confounding effects of next-generation sequencing (NGS) noise on detection of low frequency circulating tumor DNA (ctDNA) without a priori knowledge of solid tumor mutations has limited the applications of circulating cell-free DNA (ccfDNA) in clinical oncology. Here, we use a 118 gene panel and leverage ccfDNA technical replicates to eliminate NGS-associated errors while also enhancing detection of ctDNA from pancreatic ductal adenocarcinomas (PDACs). Pre-operative ccfDNA and tumor DNA were acquired from 14 patients with PDAC (78.6% stage II-III). Post-operative ccfDNA was also collected from 11 of the patients within 100 days of surgery. ctDNA detection was restricted to variants corresponding to pathogenic mutations in PDAC present in both replicates. PDAC-associated pathogenic mutations were detected in pre-operative ccfDNA in four genes (KRAS, TP53, SMAD4, ALK) from five patients. Of the nine ctDNA variants detected (variant allele frequency: 0.08%-1.59%), five had a corresponding mutation in tumor DNA. Pre-operative detection of ctDNA was associated with shorter survival (312 vs. 826 days; χ2=5.4, P = 0.021). Guiding ctDNA detection in pre-operative ccfDNA based on mutations present in tumor DNA yielded a similar survival analysis. Detection of ctDNA in the post-operative ccfDNA with or without tumor-informed guidance was not associated with outcomes. Therefore, the detection of PDAC-derived ctDNA during a broad and untargeted survey of ccfDNA with NGS may be a valuable, non-invasive, prognostic biomarker to integrate into the clinical assessment and management of patients prior to surgery.

The stochastic nature of errors in next-generation sequencing of circulating cell-free DNA.

Challenges with distinguishing circulating tumor DNA (ctDNA) from next-generation sequencing (NGS) artifacts limits variant searches to established solid tumor mutations. Here we show early and random PCR errors are a principal source of NGS noise that persist despite duplex molecular barcoding, removal of artifacts due to clonal hematopoiesis of indeterminate potential, and suppression of patterned errors. We also demonstrate sample duplicates are necessary to eliminate the stochastic noise associated with NGS. Integration of sample duplicates into NGS analytics may broaden ctDNA applications by removing NGS-related errors that confound identification of true very low frequency variants during searches for ctDNA without a priori knowledge of specific mutations to target.

Viral miRNA adaptor differentially recruits miRNAs to target mRNAs through alternative base-pairing.

HSUR2 is a viral non-coding RNA (ncRNA) that functions as a microRNA (miRNA) adaptor. HSUR2 inhibits apoptosis in infected cells by recruiting host miRNAs miR-142-3p and miR-16 to mRNAs encoding apoptotic factors. HSUR2's target recognition mechanism is not understood. It is also unknown why HSUR2 utilizes miR-16 to downregulate only a subset of transcripts. We developed a general method for individual-nucleotide resolution RNA-RNA interaction identification by crosslinking and capture (iRICC) to identify sequences mediating interactions between HSUR2 and target mRNAs in vivo. Mutational analyses confirmed identified HSUR2-mRNA interactions and validated iRICC as a method that confidently determines sequences mediating RNA-RNA interactions in vivo. We show that HSUR2 does not display a 'seed' region to base-pair with most target mRNAs, but instead uses different regions to interact with different transcripts. We further demonstrate that this versatile mode of interaction via variable base-pairing provides HSUR2 with a mechanism for differential miRNA recruitment.

Automated size selection for short cell-free DNA fragments enriches for circulating tumor DNA and improves error correction during next generation sequencing.

Circulating tumor-derived cell-free DNA (ctDNA) enables non-invasive diagnosis, monitoring, and treatment susceptibility testing in human cancers. However, accurate detection of variant alleles, particularly during untargeted searches, remains a principal obstacle to widespread application of cell-free DNA in clinical oncology. In this study, isolation of short cell-free DNA fragments is shown to enrich for tumor variants and improve correction of PCR- and sequencing-associated errors. Subfractions of the mononucleosome of circulating cell-free DNA (ccfDNA) were isolated from patients with melanoma, pancreatic ductal adenocarcinoma, and colorectal adenocarcinoma using a high-throughput-capable automated gel-extraction platform. Using a 128-gene (128 kb) custom next-generation sequencing panel, variant alleles were on average 2-fold enriched in the short fraction (median insert size: ~142 bp) compared to the original ccfDNA sample, while 0.7-fold reduced in the fraction corresponding to the principal peak of the mononucleosome (median insert size: ~167 bp). Size-selected short fractions compared to the original ccfDNA yielded significantly larger family sizes (i.e., PCR duplicates) during in silico consensus sequence interpretation via unique molecular identifiers. Increments in family size were associated with a progressive reduction of PCR and sequencing errors. Although consensus read depth also decreased at larger family sizes, the variant allele frequency in the short ccfDNA fraction remained consistent, while variant detection in the original ccfDNA was commonly lost at family sizes necessary to minimize errors. These collective findings support the automated extraction of short ccfDNA fragments to enrich for ctDNA while concomitantly reducing false positives through in silico error correction.

Conserved roles of mouse DUX and human DUX4 in activating cleavage-stage genes and MERVL/HERVL retrotransposons.

To better understand transcriptional regulation during human oogenesis and preimplantation development, we defined stage-specific transcription, which highlighted the cleavage stage as being highly distinctive. Here, we present multiple lines of evidence that a eutherian-specific multicopy retrogene, DUX4, encodes a transcription factor that activates hundreds of endogenous genes (for example, ZSCAN4, KDM4E and PRAMEF-family genes) and retroviral elements (MERVL/HERVL family) that define the cleavage-specific transcriptional programs in humans and mice. Remarkably, mouse Dux expression is both necessary and sufficient to convert mouse embryonic stem cells (mESCs) into 2-cell-embryo-like ('2C-like') cells, measured here by the reactivation of '2C' genes and repeat elements, the loss of POU5F1 (also known as OCT4) protein and chromocenters, and the conversion of the chromatin landscape (as assessed by transposase-accessible chromatin using sequencing (ATAC-seq)) to a state strongly resembling that of mouse 2C embryos. Thus, we propose mouse DUX and human DUX4 as major drivers of the cleavage or 2C state.

Adding In Silico Assessment of Potential Splice Aberration to the Integrated Evaluation of BRCA Gene Unclassified Variants.

Clinical mutation screening of the cancer susceptibility genes BRCA1 and BRCA2 generates many unclassified variants (UVs). Most of these UVs are either rare missense substitutions or nucleotide substitutions near the splice junctions of the protein coding exons. Previously, we developed a quantitative method for evaluation of BRCA gene UVs-the "integrated evaluation"-that combines a sequence analysis-based prior probability of pathogenicity with patient and/or tumor observational data to arrive at a posterior probability of pathogenicity. One limitation of the sequence analysis-based prior has been that it evaluates UVs from the perspective of missense substitution severity but not probability to disrupt normal mRNA splicing. Here, we calibrated output from the splice-site fitness program MaxEntScan to generate spliceogenicity-based prior probabilities of pathogenicity for BRCA gene variants; these range from 0.97 for variants with high probability to damage a donor or acceptor to 0.02 for exonic variants that do not impact a splice junction and are unlikely to create a de novo donor. We created a database http://priors.hci.utah.edu/PRIORS/ that provides the combined missense substitution severity and spliceogenicity-based probability of pathogenicity for BRCA gene single-nucleotide substitutions. We also updated the BRCA gene Ex-UV LOVD, available at http://hci-exlovd.hci.utah.edu, with 77 re-evaluable variants.

Potential role for snoRNAs in PKR activation during metabolic stress.

Protein kinase RNA-activated (PKR) has long been known to be activated by viral double-stranded RNA (dsRNA) as part of the mammalian immune response. However, in mice PKR is also activated by metabolic stress in the absence of viral infection, and this requires a functional kinase domain, as well as a functional dsRNA-binding domain. The endogenous cellular RNA that potentially leads to PKR activation during metabolic stress is unknown. We investigated this question using mouse embryonic fibroblast cells expressing wild-type PKR (PKRWT) or PKR with a point mutation in each dsRNA-binding motif (PKRRM). Using this system, we identified endogenous RNA that interacts with PKR after induction of metabolic stress by palmitic acid (PA) treatment. Specifically, RIP-Seq analyses showed that the majority of enriched RNAs that interacted with WT PKR (≥twofold, false discovery rate ≤ 5%) were small nucleolar RNAs (snoRNAs). Immunoprecipitation of PKR in extracts of UV-cross-linked cells, followed by RT-qPCR, confirmed that snoRNAs were enriched in PKRWT samples after PA treatment, but not in the PKRRM samples. We also demonstrated that a subset of identified snoRNAs bind and activate PKR in vitro; the presence of a 5'-triphosphate enhanced PKR activity compared with the activity with a 5'-monophosphate, for some, but not all, snoRNAs. Finally, we demonstrated PKR activation in cells upon snoRNA transfection, supporting our hypothesis that endogenous snoRNAs can activate PKR. Our results suggest an unprecedented and unexpected model whereby snoRNAs play a role in the activation of PKR under metabolic stress.

Genome-wide profiling of the C. elegans dsRNAome.

Recent studies hint that endogenous dsRNA plays an unexpected role in cellular signaling. However, a complete understanding of endogenous dsRNA signaling is hindered by an incomplete annotation of dsRNA-producing genes. To identify dsRNAs expressed in Caenorhabditis elegans, we developed a bioinformatics pipeline that identifies dsRNA by detecting clustered RNA editing sites, which are strictly limited to long dsRNA substrates of Adenosine Deaminases that act on RNA (ADAR). We compared two alignment algorithms for mapping both unique and repetitive reads and detected as many as 664 editing-enriched regions (EERs) indicative of dsRNA loci. EERs are visually enriched on the distal arms of autosomes and are predicted to possess strong internal secondary structures as well as sequence complementarity with other EERs, indicative of both intramolecular and intermolecular duplexes. Most EERs were associated with protein-coding genes, with ∼1.7% of all C. elegans mRNAs containing an EER, located primarily in very long introns and in annotated, as well as unannotated, 3' UTRs. In addition to numerous EERs associated with coding genes, we identified a population of prospective noncoding EERs that were distant from protein-coding genes and that had little or no coding potential. Finally, subsets of EERs are differentially expressed during development as well as during starvation and infection with bacterial or fungal pathogens. By combining RNA-seq with freely available bioinformatics tools, our workflow provides an easily accessible approach for the identification of dsRNAs, and more importantly, a catalog of the C. elegans dsRNAome.

Reprogramming the maternal zebrafish genome after fertilization to match the paternal methylation pattern.

Early vertebrate embryos must achieve totipotency and prepare for zygotic genome activation (ZGA). To understand this process, we determined the DNA methylation (DNAme) profiles of zebrafish gametes, embryos at different stages, and somatic muscle and compared them to gene activity and histone modifications. Sperm chromatin patterns are virtually identical to those at ZGA. Unexpectedly, the DNA of many oocyte genes important for germline functions (i.e., piwil1) or early development (i.e., hox genes) is methylated, but the loci are demethylated during zygotic cleavage stages to precisely the state observed in sperm, even in parthenogenetic embryos lacking a replicating paternal genome. Furthermore, this cohort constitutes the genes and loci that acquire DNAme during development (i.e., ZGA to muscle). Finally, DNA methyltransferase inhibition experiments suggest that DNAme silences particular gene and chromatin cohorts at ZGA, preventing their precocious expression. Thus, zebrafish achieve a totipotent chromatin state at ZGA through paternal genome competency and maternal genome DNAme reprogramming.

DNA methylation profiling in zebrafish.

DNA methylation on cytosine in vertebrates such as zebrafish serves to silence gene expression by interfering with the binding of certain transcription factors and through the recruitment of repressive chromatin machinery. Cytosine DNA methylation is chemically stable and heritable through the germline - but also reversible through many modes, making it a useful and dynamic epigenetic modification. Virtually all of the enzymes and factors involved in the deposition, binding, and removal of cytosine methylation are conserved in zebrafish, and therefore the organism an excellent model for understanding the use of DNA methylation in the control of gene regulation and other processes. Here, we discuss the main approaches to quantifying DNA methylation levels genome-wide in zebrafish: one is an established method for revealing regional methylation (methylated DNA immunoprecipitation (MeDIP)), and the other is an emerging method that reveals DNA methylation at base-pair resolution (shotgun bisulphite sequencing). We also introduce some of the analytical methods that are useful for identifying regions of hypo- or hyper-methylation, and ways to identify differentially methylated regions.

Genome-wide analysis identifies changes in histone retention and epigenetic modifications at developmental and imprinted gene loci in the sperm of infertile men.

BACKGROUND: The sperm chromatin of fertile men retains a small number of nucleosomes that are enriched at developmental gene promoters and imprinted gene loci. This unique chromatin packaging at certain gene promoters provides these genomic loci the ability to convey instructive epigenetic information to the zygote, potentially expanding the role and significance of the sperm epigenome in embryogenesis. We hypothesize that changes in chromatin packaging may be associated with poor reproductive outcome. METHODS: Seven patients with reproductive dysfunction were recruited: three had unexplained poor embryogenesis during IVF and four were diagnosed with male infertility and previously shown to have altered protamination. Genome-wide analysis of the location of histones and histone modifications was analyzed by isolation and purification of DNA bound to histones and protamines. The histone-bound fraction of DNA was analyzed using high-throughput sequencing, both initially and following chromatin immunoprecipitation. The protamine-bound fraction was hybridized to agilent arrays. DNA methylation was examined using bisulfite sequencing. RESULTS: Unlike fertile men, five of seven infertile men had non-programmatic (randomly distributed) histone retention genome-wide. Interestingly, in contrast to the total histone pool, the localization of H3 Lysine 4 methylation (H3K4me) or H3 Lysine 27 methylation (H3K27me) was highly similar in the gametes of infertile men compared with fertile men. However, there was a reduction in the amount of H3K4me or H3K27me retained at developmental transcription factors and certain imprinted genes. Finally, the methylation status of candidate developmental promoters and imprinted loci were altered in a subset of the infertile men. CONCLUSIONS: This initial genome-wide analysis of epigenetic markings in the sperm of infertile men demonstrates differences in composition and epigenetic markings compared with fertile men, especially at certain imprinted and developmental loci. Although no single locus displays a complete change in chromatin packaging or DNA modification, the data suggest that moderate changes throughout the genome exist and may have a cumulative detrimental effect on fecundity.

Genome-wide RNA-seq analysis of human and mouse platelet transcriptomes.

Inbred mice are a useful tool for studying the in vivo functions of platelets. Nonetheless, the mRNA signature of mouse platelets is not known. Here, we use paired-end next-generation RNA sequencing (RNA-seq) to characterize the polyadenylated transcriptomes of human and mouse platelets. We report that RNA-seq provides unprecedented resolution of mRNAs that are expressed across the entire human and mouse genomes. Transcript expression and abundance are often conserved between the 2 species. Several mRNAs, however, are differentially expressed in human and mouse platelets. Moreover, previously described functional disparities between mouse and human platelets are reflected in differences at the transcript level, including protease activated receptor-1, protease activated receptor-3, platelet activating factor receptor, and factor V. This suggests that RNA-seq is a useful tool for predicting differences in platelet function between mice and humans. Our next-generation sequencing analysis provides new insights into the human and murine platelet transcriptomes. The sequencing dataset will be useful in the design of mouse models of hemostasis and a catalyst for discovery of new functions of platelets. Access to the dataset is found in the "Introduction."

ENCODE tiling array analysis identifies differentially expressed annotated and novel 5' capped RNAs in hepatitis C infected liver.

Microarray studies of chronic hepatitis C infection have provided valuable information regarding the host response to viral infection. However, recent studies of the human transcriptome indicate pervasive transcription in previously unannotated regions of the genome and that many RNA transcripts have short or lack 3' poly(A) ends. We hypothesized that using ENCODE tiling arrays (1% of the genome) in combination with affinity purifying Pol II RNAs by their unique 5' m7GpppN cap would identify previously undescribed annotated and unannotated genes that are differentially expressed in liver during hepatitis C virus (HCV) infection. Both 5'-capped and poly(A)+ populations of RNA were analyzed using ENCODE tiling arrays. Sixty-four annotated genes were significantly increased in HCV cirrhotic as compared to control liver; twenty-seven (42%) of these genes were identified only by analyzing 5' capped RNA. Thirty-one annotated genes were significantly decreased; sixteen (50%) of these were identified only by analyzing 5' capped RNA. Bioinformatic analysis showed that capped RNA produced more consistent results, provided a more extensive expression profile of intronic regions and identified upregulated Pol II transcriptionally active regions in unannotated areas of the genome in HCV cirrhotic liver. Two of these regions were verified by PCR and RACE analysis. qPCR analysis of liver biopsy specimens demonstrated that these unannotated transcripts, as well as IRF1, TRIM22 and MET, were also upregulated in hepatitis C with mild inflammation and no fibrosis. The analysis of 5' capped RNA in combination with ENCODE tiling arrays provides additional gene expression information and identifies novel upregulated Pol II transcripts not previously described in HCV infected liver. This approach, particularly when combined with new RNA sequencing technologies, should also be useful in further defining Pol II transcripts differentially regulated in specific disease states and in studying RNAs regulated by changes in pre-mRNA splicing or 3' polyadenylation status.

Next generation tools for genomic data generation, distribution, and visualization.

BACKGROUND: With the rapidly falling cost and availability of high throughput sequencing and microarray technologies, the bottleneck for effectively using genomic analysis in the laboratory and clinic is shifting to one of effectively managing, analyzing, and sharing genomic data. RESULTS: Here we present three open-source, platform independent, software tools for generating, analyzing, distributing, and visualizing genomic data. These include a next generation sequencing/microarray LIMS and analysis project center (GNomEx); an application for annotating and programmatically distributing genomic data using the community vetted DAS/2 data exchange protocol (GenoPub); and a standalone Java Swing application (GWrap) that makes cutting edge command line analysis tools available to those who prefer graphical user interfaces. Both GNomEx and GenoPub use the rich client Flex/Flash web browser interface to interact with Java classes and a relational database on a remote server. Both employ a public-private user-group security model enabling controlled distribution of patient and unpublished data alongside public resources. As such, they function as genomic data repositories that can be accessed manually or programmatically through DAS/2-enabled client applications such as the Integrated Genome Browser. CONCLUSIONS: These tools have gained wide use in our core facilities, research laboratories and clinics and are freely available for non-profit use. See http://sourceforge.net/projects/gnomex/, http://sourceforge.net/projects/genoviz/, and http://sourceforge.net/projects/useq.

Dicer's helicase domain is required for accumulation of some, but not all, C. elegans endogenous siRNAs.

Years after the discovery that Dicer is a key enzyme in gene silencing, the role of its helicase domain remains enigmatic. Here we show that this domain is critical for accumulation of certain endogenous small interfering RNAs (endo-siRNAs) in Caenorhabditis elegans. The domain is required for the production of the direct products of Dicer, or primary endo-siRNAs, and consequently affects levels of downstream intermediates, the secondary endo-siRNAs. Consistent with the role of endo-siRNAs in silencing, their loss correlates with an increase in cognate mRNA levels. We find that the helicase domain of Dicer is not necessary for microRNA (miRNA) processing, or RNA interference following exposure to exogenous double-stranded RNA. Comparisons of wild-type and helicase-defective strains using deep-sequencing analyses show that the helicase domain is required by a subset of annotated endo-siRNAs, in particular, those associated with the slightly longer 26-nucleotide small RNA species containing a 5' guanosine.