Dissecting the biological relationship between TCGA miRNA and mRNA sequencing data using MMiRNA-Viewer.

BACKGROUND: MicroRNAs (miRNA) are short nucleotides that interact with their target genes through 3' untranslated regions (UTRs). The Cancer Genome Atlas (TCGA) harbors an increasing amount of cancer genome data for both tumor and normal samples. However, there are few visualization tools focusing on concurrently displaying important relationships and attributes between miRNAs and mRNAs of both cancer tumor and normal samples. Moreover, a deep investigation of miRNA-mRNA target and biological relationships across multiple cancer types by integrating web-based analysis has not been thoroughly conducted. RESULTS: We developed an interactive visualization tool called MMiRNA-Viewer that can concurrently present the co-relationships of expression between miRNA-mRNA pairs of both tumor and normal samples into a single graph. The input file of MMiRNA-Viewer contains the expression information including fold changes between normal and tumor samples for mRNAs and miRNAs, the correlation between mRNA and miRNA, and the predicted target relationship by a number of databases. Users can also load their own input data into MMiRNA-Viewer and visualize and compare detailed information about cancer-related gene expression changes, and also changes in the expression of transcription-regulating miRNAs. To validate the MMiRNA-Viewer, eight types of TCGA cancer datasets with both normal and control samples were selected in this study and three filter steps were applied subsequently. We performed Gene Ontology (GO) analysis for genes available in final selected 238 pairs and also for genes in the top 5 % (95 percentile) for each of eight cancer types to report a significant number of genes involved in various biological functions and pathways. We also calculated various centrality measurement matrices for the largest connected component(s) in each of eight cancers and reported top genes and miRNAs with high centrality measurements. CONCLUSIONS: With its user-friendly interface, dynamic visualization and advanced queries, we also believe MMiRNA-Viewer offers an intuitive approach for visualizing and elucidating co-relationships between miRNAs and mRNAs of both tumor and normal samples. We suggest that miRNA and mRNA pairs with opposite fold changes of their expression and with inverted correlation values between tumor and normal samples might be most relevant for explaining the decoupling of mRNAs and their targeting miRNAs in tumor samples for certain cancer types.

Correlating bladder cancer risk genes with their targeting microRNAs using MMiRNA-Tar.

The Cancer Genome Atlas (TCGA) (http://cancergenome.nih.gov) is a valuable data resource focused on an increasing number of well-characterized cancer genomes. In part, TCGA provides detailed information about cancer-dependent gene expression changes, including changes in the expression of transcription-regulating microRNAs. We developed a web interface tool MMiRNA-Tar (http://bioinf1.indstate.edu/MMiRNA-Tar) that can calculate and plot the correlation of expression for mRNA-microRNA pairs across samples or over a time course for a list of pairs under different prediction confidence cutoff criteria. Prediction confidence was established by requiring that the proposed mRNA-microRNA pair appears in at least one of three target prediction databases: TargetProfiler, TargetScan, or miRanda. We have tested our MMiRNA-Tar tool through analyzing 53 tumor and 11 normal samples of bladder urothelial carcinoma (BLCA) datasets obtained from TCGA and identified 204 microRNAs. These microRNAs were correlated with the mRNAs of five previously-reported bladder cancer risk genes and these selected pairs exhibited correlations in opposite direction between the tumor and normal samples based on the customized cutoff criterion of prediction. Furthermore, we have identified additional 496 genes (830 pairs) potentially targeted by 79 significant microRNAs out of 204 using three cutoff criteria, i.e., false discovery rate (FDR)<0.1, opposite correlation coefficient between the tumor and normal samples, and predicted by at least one of three target prediction databases. Therefore, MMiRNA-Tar provides researchers a convenient tool to visualize the co-relationship between microRNAs and mRNAs and to predict their targeting relationship. We believe that correlating expression profiles for microRNAs and mRNAs offers a complementary approach for elucidating their interactions.

Insights into the evolution of longevity from the bowhead whale genome.

The bowhead whale (Balaena mysticetus) is estimated to live over 200 years and is possibly the longest-living mammal. These animals should possess protective molecular adaptations relevant to age-related diseases, particularly cancer. Here, we report the sequencing and comparative analysis of the bowhead whale genome and two transcriptomes from different populations. Our analysis identifies genes under positive selection and bowhead-specific mutations in genes linked to cancer and aging. In addition, we identify gene gain and loss involving genes associated with DNA repair, cell-cycle regulation, cancer, and aging. Our results expand our understanding of the evolution of mammalian longevity and suggest possible players involved in adaptive genetic changes conferring cancer resistance. We also found potentially relevant changes in genes related to additional processes, including thermoregulation, sensory perception, dietary adaptations, and immune response. Our data are made available online (http://www.bowhead-whale.org) to facilitate research in this long-lived species.

A study on the distribution of 37 well conserved families of C2H2 zinc finger genes in eukaryotes.

BACKGROUND: The C2H2 zinc-finger (ZNF) containing gene family is one of the largest and most complex gene families in metazoan genomes. These genes are known to exist in almost all eukaryotes, and they constitute a major subset of eukaryotic transcription factors. The genes of this family usually occur as clusters in genomes and are thought to have undergone a massive expansion in vertebrates by multiple tandem duplication events (BMC Evol Biol 8:176, 2008). RESULTS: In this study, we combined two popular approaches for homolog detection, Reciprocal Best Hit (RBH) (Proc Natl Acad Sci USA 95:6239-6244, 1998) and Hidden-Markov model (HMM) profiles search (Bioinformatics 14:755-763, 1998), on a diverse set of complete genomes of 124 eukaryotic species ranging from excavates to humans to identify all detectable members of 37 C2H2 ZNF gene families. We succeeded in identifying 3,890 genes as distinct members of 37 C2H2 gene families. These 37 families are distributed among the eukaryotes as progressive additions of gene blocks with increasing complexity of the organisms. The first block featuring the protists had 7 families, the second block featuring plants had 2 families, the third block featuring the fungi had 2 families (one of which was also present in plants) and the final block consisted of metazoans with 25 families. Among the metazoans, the simpler unicellular metazoans had just 15 of the 25 families while most of the bilaterians had all 25 families making up a total of 37 families. Multiple potential examples of lineage-specific gene duplications and gene losses were also observed. CONCLUSIONS: Our hybrid approach combines features of the both RBH and HMM methods for homolog detection. This largely automated technique is much faster than manual methods and is able to detect homologs accurately and efficiently among a diverse set of organisms. Our analysis of the 37 evolutionarily conserved C2H2 ZNF gene families revealed a stepwise appearance of ZNF families, agreeing well with the phylogenetic relationship of the organisms compared and their presumed stepwise increase in complexity (Science 300:1694, 2003).

PolyLens: software for map-based visualisation and analysis of genome-scale polymorphism data.

Software tools for the flexible examination of genomic sequence information derived from populations of organisms in a geospatial context are few in number, closely tied to Web-based resources, generally focused on one or a few loci or haplotypes, and typically produce a global phylogeny as a summary of relatedness. We sought instead to produce a portable, self-contained analysis tool that is efficiently focused on a geospatial display of specifically chosen polymorphism frequencies or combination frequencies from very large data sets of genome-scale sequence from multiple individuals. PolyLens is a Java-based, integral visual analytical toolkit which can systematically process population genomic data, visualise geographic distributions of genealogical lineages, and display allele distribution patterns. PolyLens is designed for users to visualise specific DNA sequences within each individual and its related location information in the existing data set.

Whole genome phylogeny for 21 Drosophila species using predicted 2b-RAD fragments.

Type IIB restriction endonucleases are site-specific endonucleases that cut both strands of double-stranded DNA upstream and downstream of their recognition sequences. These restriction enzymes have recognition sequences that are generally interrupted and range from 5 to 7 bases long. They produce DNA fragments which are uniformly small, ranging from 21 to 33 base pairs in length (without cohesive ends). The fragments are generated from throughout the entire length of a genomic DNA providing an excellent fractional representation of the genome. In this study we simulated restriction enzyme digestions on 21 sequenced genomes of various Drosophila species using the predicted targets of 16 Type IIB restriction enzymes to effectively produce a large and arbitrary selection of loci from these genomes. The fragments were then used to compare organisms and to calculate the distance between genomes in pair-wise combination by counting the number of shared fragments between the two genomes. Phylogenetic trees were then generated for each enzyme using this distance measure and the consensus was calculated. The consensus tree obtained agrees well with the currently accepted tree for the Drosophila species. We conclude that multi-locus sub-genomic representation combined with next generation sequencing, especially for individuals and species without previous genome characterization, can accelerate studies of comparative genomics and the building of accurate phylogenetic trees.

Whole genome phylogenies for multiple Drosophila species.

BACKGROUND: Reconstructing the evolutionary history of organisms using traditional phylogenetic methods may suffer from inaccurate sequence alignment. An alternative approach, particularly effective when whole genome sequences are available, is to employ methods that don't use explicit sequence alignments. We extend a novel phylogenetic method based on Singular Value Decomposition (SVD) to reconstruct the phylogeny of 12 sequenced Drosophila species. SVD analysis provides accurate comparisons for a high fraction of sequences within whole genomes without the prior identification of orthologs or homologous sites. With this method all protein sequences are converted to peptide frequency vectors within a matrix that is decomposed to provide simplified vector representations for each protein of the genome in a reduced dimensional space. These vectors are summed together to provide a vector representation for each species, and the angle between these vectors provides distance measures that are used to construct species trees. RESULTS: An unfiltered whole genome analysis (193,622 predicted proteins) strongly supports the currently accepted phylogeny for 12 Drosophila species at higher dimensions except for the generally accepted but difficult to discern sister relationship between D. erecta and D. yakuba. Also, in accordance with previous studies, many sequences appear to support alternative phylogenies. In this case, we observed grouping of D. erecta with D. sechellia when approximately 55% to 95% of the proteins were removed using a filter based on projection values or by reducing resolution by using fewer dimensions. Similar results were obtained when just the melanogaster subgroup was analyzed. CONCLUSIONS: These results indicate that using our novel phylogenetic method, it is possible to consult and interpret all predicted protein sequences within multiple whole genomes to produce accurate phylogenetic estimations of relatedness between Drosophila species. Furthermore, protein filtering can be effectively applied to reduce incongruence in the dataset as well as to generate alternative phylogenies.

Integrating HIV prevention in reproductive health settings.

CONTEXT: This article describes results of a process evaluation of a cooperative agreement between the Centers for Disease Control and Prevention's Division of Reproductive Health and 10 regional training centers to increase the number of reproductive health (RH) settings that integrate human immunodeficiency virus (HIV) prevention services at an appropriate level into routine care. OBJECTIVE: Our goal was to learn about the process of integrating HIV prevention into RH settings. DESIGN: We conducted a retrospective evaluation, using qualitative methods. SETTING: The clinics were from 10 US Department of Health and Human Services regions. PARTICIPANTS: We interviewed 16 key informants from 10 selected model clinics. MAIN OUTCOME MEASURES: The main outcome was organization change. RESULTS: The most common obstacles to integration were staff issues, logistics barriers, inadequate clinic structure to support integration, and staff training barriers. Using the transtheoretical model (TTM) applied to organizations, we documented organizational change as informants described their clinics' progression to integration and overcoming obstacles. All model clinics began in the contemplation stage of transtheoretical model. Every clinic exhibited at least 1 process of change for every stage. In the contemplation stage, most informants discussed fears about not changing, stated that the integration was consistent with the agency's mission, and described thinking about commitment to the change. In the preparation stage, all informants described building teams that supported integration of HIV prevention. During the action stage, informants talked about assessments of facilities, staff and protocols, commitments through grants or agreements, and then using training to support new behaviors and adopting new cognitions. In the maintenance stage, all reported changing policies, procedures, or protocols, most promoted helping relationships among the staff, and nearly all reported rewards for the new ways of working. CONCLUSIONS: RH settings were able to integrate HIV prevention services by employing a systematic process.

A survey of well conserved families of C2H2 zinc-finger genes in Daphnia.

BACKGROUND: A recent comparative genomic analysis tentatively identified roughly 40 orthologous groups of C2H2 Zinc-finger proteins that are well conserved in "bilaterians" (i.e. worms, flies, and humans). Here we extend that analysis to include a second arthropod genome from the crustacean, Daphnia pulex. RESULTS: Most of the 40 orthologous groups of C2H2 zinc-finger proteins are represented by just one or two proteins within each of the previously surveyed species. Likewise, Daphnia were found to possess a similar number of orthologs for all of these small orthology groups. In contrast, the number of Sp/KLF homologs tends to be greater and to vary between species. Like the corresponding mammalian Sp/KLF proteins, most of the Drosophila and Daphnia homologs can be placed into one of three sub-groups: Class I-III. Daphnia were found to have three Class I proteins that roughly correspond to their Drosophila counterparts, dSP1, btd, CG5669, and three Class II proteins that roughly correspond to Luna, CG12029, CG9895. However, Daphnia have four additional KLF-Class II proteins that are most similar to the vertebrate KLF1/2/4 proteins, a subset not found in Drosophila. Two of these four proteins are encoded by genes linked in tandem. Daphnia also have three KLF-Class III members, one more than Drosophila. One of these is a likely Bteb2 homolog, while the other two correspond to Cabot and KLF13, a vertebrate homolog of Cabot. CONCLUSION: Consistent with their likely roles as fundamental determinants of bilaterian form and function, most of the 40 groups of C2H2 zinc-finger proteins are conserved in kind and number in Daphnia. However, the KLF family includes several additional genes that are most similar to genes present in vertebrates but missing in Drosophila.

An SVD-based comparison of nine whole eukaryotic genomes supports a coelomate rather than ecdysozoan lineage.

BACKGROUND: Eukaryotic whole genome sequences are accumulating at an impressive rate. Effective methods for comparing multiple whole eukaryotic genomes on a large scale are needed. Most attempted solutions involve the production of large scale alignments, and many of these require a high stringency pre-screen for putative orthologs in order to reduce the effective size of the dataset and provide a reasonably high but unknown fraction of correctly aligned homologous sites for comparison. As an alternative, highly efficient methods that do not require the pre-alignment of operationally defined orthologs are also being explored. RESULTS: A non-alignment method based on the Singular Value Decomposition (SVD) was used to compare the predicted protein complement of nine whole eukaryotic genomes ranging from yeast to man. This analysis resulted in the simultaneous identification and definition of a large number of well conserved motifs and gene families, and produced a species tree supporting one of two conflicting hypotheses of metazoan relationships. CONCLUSIONS: Our SVD-based analysis of the entire protein complement of nine whole eukaryotic genomes suggests that highly conserved motifs and gene families can be identified and effectively compared in a single coherent definition space for the easy extraction of gene and species trees. While this occurs without the explicit definition of orthologs or homologous sites, the analysis can provide a basis for these definitions.

Transgene manipulation in zebrafish by using recombinases.

Although much remains to be done, our results to date suggest that efficient and precise genome engineering in zebrafish will be possible in the future by using Cre recombinase and SB transposase in combination with their respective target sites. In this study, we provide the first evidence that Cre recombinase can mediate effective site-specific deletion of transgenes in zebrafish. We found that the efficiency of target site utilization could approach 100%, independent of whether the target site was provided transiently by injection or stably within an integrated transgene. Microinjection of Cre mRNA appeared to be slightly more effective for this purpose than microinjection of Cre-expressing plasmid DNA. Our work has not yet progressed to the point where SB-mediated mobilization of our transgene constructs would be observed. However, a recent report has demonstrated that SB can enhance transgenesis rates sixfold over conventional methods by efficiently mediating multiple single-copy insertion of transgenes into the zebrafish genome (Davidson et al., 2003). Therefore, it seems likely that a combined system should eventually allow both SB-mediated transgene mobilization and Cre-mediated transgene modification. Our goal is to validate methods for the precise reengineering of the zebrafish genome by using a combination of Cre-loxP and SB transposon systems. These methods can be used to delete, replace, or mobilize large pieces of DNA or to modify the genome only when and where required by the investigator. For example, it should be possible to deliver particular RNAi genes to well-expressed chromosomal loci and then exchange them easily with alternative RNAi genes for the specific suppression of alternative targets. As a nonviral vector for gene therapy, the transposon component allows for the possibility of highly efficient integration, whereas the Cre-loxP component can target the integration and/or exchange of foreign DNA into specific sites within the genome. The specificity and efficiency of this system also make it ideal for applications in which precise genome modifications are required (e.g., stock improvement). Future work should establish whether alternative recombination systems (e.g., phiC31 integrase) can improve the utility of this system. After the fish system is fully established, it would be interesting to explore its application to genome engineering in other organisms.

A comprehensive whole genome bacterial phylogeny using correlated peptide motifs defined in a high dimensional vector space.

As whole genome sequences continue to expand in number and complexity, effective methods for comparing and categorizing both genes and species represented within extremely large datasets are required. Methods introduced to date have generally utilized incomplete and likely insufficient subsets of the available data. We have developed an accurate and efficient method for producing robust gene and species phylogenies using very large whole genome protein datasets. This method relies on multidimensional protein vector definitions supplied by the singular value decomposition (SVD) of a large sparse data matrix in which each protein is uniquely represented as a vector of overlapping tetrapeptide frequencies. Quantitative pairwise estimates of species similarity were obtained by summing the protein vectors to form species vectors, then determining the cosines of the angles between species vectors. Evolutionary trees produced using this method confirmed many accepted prokaryotic relationships. However, several unconventional relationships were also noted. In addition, we demonstrate that many of the SVD-derived right basis vectors represent particular conserved protein families, while many of the corresponding left basis vectors describe conserved motifs within these families as sets of correlated peptides (copeps). This analysis represents the most detailed simultaneous comparison of prokaryotic genes and species available to date.

A comprehensive vertebrate phylogeny using vector representations of protein sequences from whole genomes.

We recently developed a method for producing comprehensive gene and species phylogenies from unaligned whole genome data using singular value decomposition (SVD) to analyze character string frequencies. This work provides an integrated gene and species phylogeny for 64 vertebrate mitochondrial genomes composed of 832 total proteins. In addition, to provide a theoretical basis for the method, we present a graphical interpretation of both the original frequency matrix and the SVD-derived matrix. These large matrices describe high-dimensional Euclidean spaces within which biomolecular sequences can be uniquely represented as vectors. In particular, the SVD-derived vector space describes each protein relative to a restricted set of newly defined, independent axes, each of which represents a novel form of conserved motif, termed a correlated peptide motif. A quantitative comparison of the relative orientations of protein vectors in this space provides accurate and straightforward estimates of sequence similarity, which can in turn be used to produce comprehensive gene trees. Alternatively, the vector representations of genes from individual species can be summed, allowing species trees to be produced.

Integrated gene and species phylogenies from unaligned whole genome protein sequences.

MOTIVATION: Most molecular phylogenies are based on sequence alignments. Consequently, they fail to account for modes of sequence evolution that involve frequent insertions or deletions. Here we present a method for generating accurate gene and species phylogenies from whole genome sequence that makes use of short character string matches not placed within explicit alignments. In this work, the singular value decomposition of a sparse tetrapeptide frequency matrix is used to represent the proteins of organisms uniquely and precisely as vectors in a high-dimensional space. Vectors of this kind can be used to calculate pairwise distance values based on the angle separating the vectors, and the resulting distance values can be used to generate phylogenetic trees. Protein trees so derived can be examined directly for homologous sequences. Alternatively, vectors defining each of the proteins within an organism can be summed to provide a vector representation of the organism, which is then used to generate species trees. RESULTS: Using a large mitochondrial genome dataset, we have produced species trees that are largely in agreement with previously published trees based on the analysis of identical datasets using different methods. These trees also agree well with currently accepted phylogenetic theory. In principle, our method could be used to compare much larger bacterial or nuclear genomes in full molecular detail, ultimately allowing accurate gene and species relationships to be derived from a comprehensive comparison of complete genomes. In contrast to phylogenetic methods based on alignments, sequences that evolve by relative insertion or deletion would tend to remain recognizably similar.