GeneFishing to reconstruct context specific portraits of biological processes.

Rapid advances in genomic technologies have led to a wealth of diverse data, from which novel discoveries can be gleaned through the application of robust statistical and computational methods. Here, we describe GeneFishing, a semisupervised computational approach to reconstruct context-specific portraits of biological processes by leveraging gene-gene coexpression information. GeneFishing incorporates multiple high-dimensional statistical ideas, including dimensionality reduction, clustering, subsampling, and results aggregation, to produce robust results. To illustrate the power of our method, we applied it using 21 genes involved in cholesterol metabolism as "bait" to "fish out" (or identify) genes not previously identified as being connected to cholesterol metabolism. Using simulation and real datasets, we found that the results obtained through GeneFishing were more interesting for our study than those provided by related gene prioritization methods. In particular, application of GeneFishing to the GTEx liver RNA sequencing (RNAseq) data not only reidentified many known cholesterol-related genes, but also pointed to glyoxalase I (GLO1) as a gene implicated in cholesterol metabolism. In a follow-up experiment, we found that GLO1 knockdown in human hepatoma cell lines increased levels of cellular cholesterol ester, validating a role for GLO1 in cholesterol metabolism. In addition, we performed pantissue analysis by applying GeneFishing on various tissues and identified many potential tissue-specific cholesterol metabolism-related genes. GeneFishing appears to be a powerful tool for identifying related components of complex biological systems and may be used across a wide range of applications.

Comparative analysis of metazoan chromatin organization.

Genome function is dynamically regulated in part by chromatin, which consists of the histones, non-histone proteins and RNA molecules that package DNA. Studies in Caenorhabditis elegans and Drosophila melanogaster have contributed substantially to our understanding of molecular mechanisms of genome function in humans, and have revealed conservation of chromatin components and mechanisms. Nevertheless, the three organisms have markedly different genome sizes, chromosome architecture and gene organization. On human and fly chromosomes, for example, pericentric heterochromatin flanks single centromeres, whereas worm chromosomes have dispersed heterochromatin-like regions enriched in the distal chromosomal 'arms', and centromeres distributed along their lengths. To systematically investigate chromatin organization and associated gene regulation across species, we generated and analysed a large collection of genome-wide chromatin data sets from cell lines and developmental stages in worm, fly and human. Here we present over 800 new data sets from our ENCODE and modENCODE consortia, bringing the total to over 1,400. Comparison of combinatorial patterns of histone modifications, nuclear lamina-associated domains, organization of large-scale topological domains, chromatin environment at promoters and enhancers, nucleosome positioning, and DNA replication patterns reveals many conserved features of chromatin organization among the three organisms. We also find notable differences in the composition and locations of repressive chromatin. These data sets and analyses provide a rich resource for comparative and species-specific investigations of chromatin composition, organization and function.

Mouse class III myosins: kinase activity and phosphorylation sites.

As class III unconventional myosins are motor proteins with an N-terminal kinase domain, it seems likely they play a role in both signaling and actin based transport. A growing body of evidence indicates that the motor functions of human class IIIA myosin, which has been implicated in progressive hearing loss, are modulated by intermolecular autophosphorylation. However, the phosphorylation sites have not been identified. We studied the kinase activity and phosphorylation sites of mouse class III myosins, mMyo3A and 3B, which are highly similar to their human orthologs. We demonstrate that the kinase domains of mMyo3A and 3B are active kinases, and that they have similar, if not identical, substrate specificities. We show that the kinase domains of these proteins autophosphorylate, and that they can phosphorylate sites within their myosin and tail domains. Using liquid chromatography-mass spectrometry, we identified phosphorylated sites in the kinase, myosin motor and tail domains of both mMyo3A and 3B. Most of the phosphorylated sites we identified and their consensus phosphorylation motifs are highly conserved among vertebrate class III myosins, including human class III myosins. Our findings are a major step toward understanding how the functions of class III myosins are regulated by phosphorylation.

Cloning and distribution of myosin 3B in the mouse retina: differential distribution in cone outer segments.

Class III myosins are important for the function and survival of photoreceptors and ciliary hair cells. Although vertebrates possess two class III myosin genes, myo3A and myo3B, recent studies have focused on Myo3A because mutations in the human gene are implicated in progressive hearing loss. Myo3B may compensate for defects in Myo3A, yet little is known about its distribution and function. This study focuses on Myo3B expression in the mouse retina. We cloned two variants of myo3B from mouse retina and determined that they are expressed early in retinal development. In this study we show for the first time in a mammal that both Myo3B and Myo3A proteins are present in inner segments of all photoreceptors. Myo3B is also present in outer segments of S opsin-immunoreactive cones but not M opsin dominant cones. Myo3B is also detected in rare cells of the inner nuclear layer and some ganglion cells. Myo3B may have diverse roles in retinal neurons. In photoreceptor inner segments Myo3B is positioned appropriately to prevent photoreceptor loss of function caused by Myo3A defects.

A new compartment at stereocilia tips defined by spatial and temporal patterns of myosin IIIa expression.

Class III myosins are motor proteins that contain an N-terminal kinase domain and a C-terminal actin-binding domain. We show that myosin IIIa, which has been implicated in nonsyndromic progressive hearing loss, is localized at stereocilia tips. Myosin IIIa progressively accumulates during stereocilia maturation in a thimble-like pattern around the stereocilia tip, distinct from the cap-like localization of myosin XVa and the shaft localization of myosin Ic. Overexpression of deletion mutants for functional domains of green fluorescent protein (GFP)-myosin IIIa shows that the motor domain, but not the actin-binding tail domain, is required for stereocilia tip localization. Deletion of the kinase domain produces stereocilia elongation and bulging of the stereocilia tips. The thimble-like localization and the influence myosin IIIa has on stereocilia shape reveal a previously unrecognized molecular compartment at the distal end of stereocilia, the site of actin polymerization as well as operation of the mechanoelectrical transduction apparatus.

Localization of a class III myosin to filopodia tips in transfected HeLa cells requires an actin-binding site in its tail domain.

Bass Myo3A, a class III myosin, was expressed in HeLa cells as a GFP fusion in order to study its cellular localization. GFP-Myo3A localized to the cytoplasm and to the tips of F-actin bundles in filopodia, a localization that is consistent with the observed concentration toward the distal ends of F-actin bundles in photoreceptor cells. A mutation in the motor active site resulted in a loss of filopodia localization, suggesting that Myo3A motor activity is required for filopodial tip localization. Deletion analyses showed that the NH2-terminal kinase domain is not required but the CO2H-terminal 22 amino acids of the Myo3A tail are required for filopodial localization. Expression of this tail fragment alone produced fluorescence associated with F-actin throughout the cytoplasm and filopodia and a recombinant tail fragment bound to F-actin in vitro. An actin-binding motif was identified within this tail fragment, and a mutation within this motif abolished both filopodia localization by Myo3A and F-actin binding by the tail fragment alone. Calmodulin localized to filopodial tips when coexpressed with Myo3A but not in the absence of Myo3A, an observation consistent with the previous proposal that class III myosins bind calmodulin and thereby localize it in certain cell types.

Myo3A, one of two class III myosin genes expressed in vertebrate retina, is localized to the calycal processes of rod and cone photoreceptors and is expressed in the sacculus.

The striped bass has two retina-expressed class III myosin genes, each composed of a kinase, motor, and tail domain. We report the cloning, sequence analysis, and expression patterns of the long (Myo3A) and short (Myo3B) class III myosins, as well as cellular localization and biochemical characterization of the long isoform, Myo3A. Myo3A (209 kDa) is expressed in the retina, brain, testis, and sacculus, and Myo3B (155 kDa) is expressed in the retina, intestine, and testis. The tails of these two isoforms contain two highly conserved domains, 3THDI and 3THDII. Whereas Myo3B has three IQ motifs, Myo3A has nine IQ motifs, four in its neck and five in its tail domain. Myo3A localizes to actin filament bundles of photoreceptors and is concentrated in the calycal processes. An anti-Myo3A antibody decorates the actin cytoskeleton of rod inner/outer segments, and this labeling is reduced by the presence of ATP. The ATP-sensitive actin association is a feature characteristic of myosin motors. The numerous IQ motifs may play a structural or signaling role in the Myo3A, and its localization to calycal processes indicates that this myosin mediates a local function at this site in vertebrate photoreceptors.

RP1-13D10.2 Is a Novel Modulator of Statin-Induced Changes in Cholesterol.

BACKGROUND: Numerous genetic contributors to cardiovascular disease risk have been identified through genome-wide association studies; however, identifying the molecular mechanism underlying these associations is not straightforward. The Justification for the Use of Statins in Primary Prevention: An Intervention Trial Evaluating Rosuvastatin (JUPITER) trial of rosuvastatin users identified a sub-genome-wide association of rs6924995, a single-nucleotide polymorphism ≈10 kb downstream of myosin regulatory light chain interacting protein (MYLIP, aka IDOL and inducible degrader of low-density lipoprotein receptor [LDLR]), with LDL cholesterol statin response. Interestingly, although this signal was initially attributed to MYLIP, rs6924995 lies within RP1-13D10.2, an uncharacterized long noncoding RNA. METHODS AND RESULTS: Using simvastatin and sham incubated lymphoblastoid cell lines from participants of the Cholesterol and Pharmacogenetics (CAP) simvastatin clinical trial, we found that statin-induced change in RP1-13D10.2 levels differed between cell lines from the tails of the white and black low-density lipoprotein cholesterol response distributions, whereas no difference in MYLIP was observed. RP1-13D10.2 overexpression in Huh7 and HepG2 increased LDLR transcript levels, increased LDL uptake, and decreased media levels of apolipoprotein B. In addition, we found a trend of slight differences in the effects of RP1-13D10.2 overexpression on LDLR transcript levels between hepatoma cells transfected with the rs6924995 A versus G allele and a suggestion of an association between rs6924995 and RP1-10D13.2 expression levels in the CAP lymphoblastoid cell lines. Finally, RP1-13D10.2 expression levels seem to be sterol regulated, consistent with its potential role as a novel lipid regulator. CONCLUSIONS: RP1-13D10.2 is a long noncoding RNA that regulates LDLR and may contribute to low-density lipoprotein cholesterol response to statin treatment. These findings highlight the potential role of noncoding RNAs as determinants of interindividual variation in drug response.

Parathyroid hormone regulation of hypoxia-inducible factor signaling in osteoblastic cells.

Osteoblasts perceive and respond to changes in their pericellular environment, including biophysical signals and oxygen availability, to elicit an anabolic or catabolic response. Parathyroid hormone (PTH) affects each arm of skeletal remodeling, with net anabolic or catabolic effects dependent upon duration of exposure. Similarly, the capacity of osteoblastic cells to perceive pericellular oxygen has a profound effect on skeletal mass and architecture, as mice expressing stable hypoxia-inducible factor (HIF)-1α and -2α demonstrate age-dependent increases in bone volume per tissue volume and osteoblast number. Further, HIF levels and signaling can be influenced in an oxygen-independent manner. Because the cellular mechanisms involved in PTH regulation of the skeleton remain vague, we sought whether PTH could influence HIF-1α expression and HIF-α-driven luciferase activity independently of altered oxygen availability. Using UMR106.01 mature osteoblasts, we observed that 100nM hPTH(1-34) decreased HIF-1α and HIF-responsive luciferase activity in a process involving heat shock protein 90 (Hsp90) and cyclic AMP but not intracellular calcium. Altering activity of the small GTPase RhoA and its effector kinase ROCK altered HIF-α-driven luciferase activity in the absence and presence of PTH. Taken together, these data introduce PTH as a regulator of oxygen-independent HIF-1α levels through a mechanism involving cyclic AMP, Hsp90, and the cytoskeleton.

Myosin IIIB uses an actin-binding motif in its espin-1 cargo to reach the tips of actin protrusions.

Myosin IIIA (MYO3A) targets actin protrusion tips using a motility mechanism dependent on both motor and tail actin-binding activity [1]. We show that myosin IIIB (MYO3B) lacks tail actin-binding activity and is unable to target COS7 cell filopodia tips, yet is somehow able to target stereocilia tips. Strikingly, when MYO3B is coexpressed with espin-1 (ESPN1), a MYO3A cargo protein endogenously expressed in stereocilia [2], MYO3B targets and carries ESPN1 to COS7 filopodia tips. We show that this tip localization is lost when we remove the ESPN1 C terminus actin-binding site. We also demonstrate that, like MYO3A [2], MYO3B can elongate filopodia by transporting ESPN1 to the polymerizing end of actin filaments. The mutual dependence of MYO3B and ESPN1 for tip localization reveals a novel mechanism for the cell to regulate myosin tip localization via a reciprocal relationship with cargo that directly participates in actin binding for motility. Our results are consistent with a novel form of motility for class III myosins that requires both motor and tail domain actin-binding activity and show that the actin-binding tail can be replaced by actin-binding cargo. This study also provides a framework to better understand the late-onset hearing loss phenotype in patients with MYO3A mutations.

Integrative analysis of the Caenorhabditis elegans genome by the modENCODE project.

We systematically generated large-scale data sets to improve genome annotation for the nematode Caenorhabditis elegans, a key model organism. These data sets include transcriptome profiling across a developmental time course, genome-wide identification of transcription factor-binding sites, and maps of chromatin organization. From this, we created more complete and accurate gene models, including alternative splice forms and candidate noncoding RNAs. We constructed hierarchical networks of transcription factor-binding and microRNA interactions and discovered chromosomal locations bound by an unusually large number of transcription factors. Different patterns of chromatin composition and histone modification were revealed between chromosome arms and centers, with similarly prominent differences between autosomes and the X chromosome. Integrating data types, we built statistical models relating chromatin, transcription factor binding, and gene expression. Overall, our analyses ascribed putative functions to most of the conserved genome.

Myosin IIIa boosts elongation of stereocilia by transporting espin 1 to the plus ends of actin filaments.

Two proteins implicated in inherited deafness, myosin IIIa, a plus-end-directed motor, and espin, an actin-bundling protein containing the actin-monomer-binding motif WH2, have been shown to influence the length of mechanosensory stereocilia. Here we report that espin 1, an ankyrin repeat-containing isoform of espin, colocalizes with myosin IIIa at stereocilia tips and interacts with a unique conserved domain of myosin IIIa. We show that combined overexpression of these proteins causes greater elongation of stereocilia, compared with overexpression of either myosin IIIa alone or espin 1 alone. When these two proteins were co-expressed in the fibroblast-like COS-7 cell line they induced a tenfold elongation of filopodia. This extraordinary filopodia elongation results from the transport of espin 1 to the plus ends of F-actin by myosin IIIa and depends on espin 1 WH2 activity. This study provides the basis for understanding the role of myosin IIIa and espin 1 in regulating stereocilia length, and presents a physiological example where myosins can boost elongation of actin protrusions by transporting actin regulatory factors to the plus ends of actin filaments.

The kinase domain alters the kinetic properties of the myosin IIIA motor.

Myosin IIIA is unique among myosin proteins in that it contains an N-terminal kinase domain capable of autophosphorylating sites on the motor domain. A construct of myosin IIIA lacking the kinase domain localizes more efficiently to the stereocilia tips and alters the morphology of the tips in inner ear hair cells. Therefore, we performed a kinetic analysis of myosin IIIA without the kinase domain (MIII DeltaK) and compared these results with our reported analysis of myosin IIIA containing the kinase domain (MIII). The steady-state kinetic properties of MIII DeltaK indicate that it has a 2-fold higher maximum actin-activated ATPase rate (kcat = 1.5 +/- 0.1 s-1) and a 5-fold tighter actin affinity (KATPase = 6.0 +/- 1.4 microM, and KActin = 1.4 +/- 0.4 microM) compared to MIII. The rate of ATP binding to the motor domain is enhanced in MIII DeltaK (K1k+2 approximately 0.10 +/- 0.01 microM-1.s-1) to a level similar to the rate of binding to MIII in the presence of actin. The rate of ATP hydrolysis in the absence of actin is slow and may be rate limiting. Actin-activated phosphate release is identical with and without the kinase domain. The transition between actomyosin.ADP states, which is rate limiting in MIII, is enhanced in MIII DeltaK. MIII DeltaK accumulates more efficiently at the tips of filopodia in HeLa cells. Our results suggest a model in which the activity and concentration of myosin IIIA localized to the tips of actin bundles mediates the morphology of the tips in sensory cells.

Kinetic mechanism of human myosin IIIA.

Myosin IIIA is specifically expressed in photoreceptors and cochlea and is important for the phototransduction and hearing processes. In addition, myosin IIIA contains a unique N-terminal kinase domain and C-terminal tail actin-binding motif. We examined the kinetic properties of baculovirus expressed human myosin IIIA containing the kinase, motor, and two IQ domains. The maximum actin-activated ATPase rate is relatively slow (k(cat) = 0.77 +/- 0.08 s(-1)), and high actin concentrations are required to fully activate the ATPase rate (K(ATPase) = 34 +/- 11 microm). However, actin co-sedimentation assays suggest that myosin III has a relatively high steady-state affinity for actin in the presence of ATP (K(actin) approximately 7 microm). The rate of ATP binding to the motor domain is quite slow both in the presence and absence of actin (K(1)k(+2) = 0.020 and 0.001 microm(-1).s(-1), respectively). The rate of actin-activated phosphate release is more than 100-fold faster (85 s(-1)) than the k(cat), whereas ADP release in the presence of actin follows a two-step mechanism (7.0 and 0.6 s(-1)). Thus, our data suggest a transition between two actomyosin-ADP states is the rate-limiting step in the actomyosin III ATPase cycle. Our data also suggest the myosin III motor spends a large fraction of its cycle in an actomyosin ADP state that has an intermediate affinity for actin (K(d) approximately 5 microm). The long lived actomyosin-ADP state may be important for the ability of myosin III to function as a cellular transporter and actin cross-linker in the actin bundles of sensory cells.

A class III myosin expressed in the retina is a potential candidate for Bardet-Biedl syndrome.

Class III myosins are actin-based motors with amino-terminal kinase domains. Expression of these motors is highly enhanced in retinal photoreceptors. As mutations in the gene encoding NINAC, a Drosophila melanogaster class III myosin, cause retinal degeneration, human homologs of this gene are potential candidates for human retinal disease. We have recently reported the cloning of MYO3A, a human myosin III expressed predominantly in the retina and retinal pigmented epithelium [1]. The map locus of MYO3A is close to, but does not overlap, that of human Usher's 1F [2]. Here we introduce a shorter class III myosin isoform, MYO3B, which is expressed in the retina, kidney, and testis. We describe the cDNA sequence, genomic organization, and splice variants of MYO3B expressed in the human retina. A product of 36 exons, MYO3B has several splice variants containing either one or two calmodulin binding (IQ) motifs in the neck domain and one of three predominant tail variations: a short tail ending just past the second IQ motif, or two alternatively spliced longer tails. MYO3B maps to 2q31.1-q31.2, a region that overlaps the locus for a Bardet-Biedl syndrome (BBS5) linked to markers at 2q31 [3].