PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “Biological sequence analysis”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 757 records · Page 42Linked to original sources

Sequences that direct significant levels of frameshifting are frequent in coding regions of Escherichia coli.

It is generally believed that significant ribosomal frameshifting during translation does not occur without a functional purpose. The distribution of two frameshift-prone sequences, A_AAA_AAG and CCC_TGA, in coding regions of Escherichia coli has been analyzed. Although a moderate level of selection against the first sequence is evident, 68 genes contain A_AAA_AAG and 19 contain CCC_TGA. The majority of those tested in their genomic context showed >1% frameshifting. Comparative sequence analysis was employed to assess a potential biological role for frameshifting in decoding these genes. Two new candidates, in pheL and ydaY, for utilized frameshifting have been identified in addition to those previously known in dnaX and nine insertion sequence elements. For the majority of the shift-prone sequences no functional role can be attributed to them, and the frameshifting is likely erroneous. However, none of frameshift sequences is in the 306 most highly expressed genes. The unexpected conclusion is that moderate frameshifting during expression of at least some other genes is not sufficiently harmful for cells to trigger strong negative evolutionary pressure.

Base Sequence↗

Transcriptomic analysis of the cardiac left ventricle in a rodent model of diabetic cardiomyopathy: molecular snapshot of a severe myocardial disease.

Heart disease is the major cause of death in diabetes, a disorder characterized by chronic hyperglycemia and cardiovascular complications. Diabetic cardiomyopathy (DCM) is increasingly recognized as a major contributor to diastolic dysfunction and heart failure in diabetes, but its molecular basis has remained obscure, in part because of its multifactorial origins. Here we employed comparative transcriptomic methods with quantitative verification of selected transcripts by reverse transcriptase quantitative PCR to characterize the molecular basis of DCM in rats with streptozotocin-induced diabetes of 16-wk duration. Diabetes caused left ventricular disease that was accompanied by significant changes in the expression of 1,614 genes, 749 of which had functions assignable by Gene Ontology classification. Genes corresponding to proteins expressed in mitochondria accounted for a disproportionate number of those whose expression was significantly modified in DCM, consistent with the idea that the mitochondrion is a key target of the pathogenic processes that cause myocardial disease in diabetes. Diabetes also induced global perturbations in the expression of genes regulating cardiac fatty acid metabolism, whose dysfunction is likely to play a key role in the promotion of oxidative stress, thereby contributing to the pathogenesis of diabetic myocardial disease. In particular, these data point to impaired regulation of mitochondrial beta-oxidation as central in the mechanisms that generate DCM pathogenesis. This study provides a comprehensive molecular snapshot of the processes leading to myocardial disease in diabetes.

Animals↗

Transcriptional analysis of biofilm formation processes in the anaerobic, hyperthermophilic bacterium Thermotoga maritima.

Thermotoga maritima, a fermentative, anaerobic, hyperthermophilic bacterium, was found to attach to bioreactor glass walls, nylon mesh, and polycarbonate filters during chemostat cultivation on maltose-based media at 80 degrees C. A whole-genome cDNA microarray was used to examine differential expression patterns between biofilm and planktonic populations. Mixed-model statistical analysis revealed differential expression (twofold or more) of 114 open reading frames in sessile cells (6% of the genome), over a third of which were initially annotated as hypothetical proteins in the T. maritima genome. Among the previously annotated genes in the T. maritima genome, which showed expression changes during biofilm growth, were several that corresponded to biofilm formation genes identified in mesophilic bacteria (i.e., Pseudomonas species, Escherichia coli, and Staphylococcus epidermidis). Most notably, T. maritima biofilm-bound cells exhibited increased transcription of genes involved in iron and sulfur transport, as well as in biosynthesis of cysteine, thiamine, NAD, and isoprenoid side chains of quinones. These findings were all consistent with the up-regulation of iron-sulfur cluster assembly and repair functions in biofilm cells. Significant up-regulation of several beta-specific glycosidases was also noted in biofilm cells, despite the fact that maltose was the primary carbon source fed to the chemostat. The reasons for increased beta-glycosidase levels are unclear but are likely related to the processing of biofilm-based polysaccharides. In addition to revealing insights into the phenotype of sessile T. maritima communities, the methodology developed here can be extended to study other anaerobic biofilm formation processes as well as to examine aspects of microbial ecology in hydrothermal environments.

ATP-Binding Cassette Transporters↗

ISYMOD: a knowledge warehouse for the identification, assembly and analysis of bacterial integrated systems.

MOTIVATION: Complex biological functions emerge from interactions between proteins in stable supra-molecular assemblies and/or through transitory contacts. Most of the time protein partners of the assemblies are composed of one or several domains which exhibit different biochemical functions. Thus the study of cellular process requires the identification of different functional units and their integration in an interaction network; such complexes are referred to as integrated systems. In order to exploit with optimum efficiency the increased release of data, automated bioinformatics strategies are needed to identify, reconstruct and model such systems. For that purpose, we have developed a knowledge warehouse dedicated to the representation and acquisition of bacterial integrated systems involved in the exchange of the bacterial cell with its environment. RESULTS: ISYMOD is a knowledge warehouse that consistently integrates in the same environment the data and the methods used for their acquisition. This is achieved through the construction of (1) a domain knowledge base (DKB) devoted to the storage of the knowledge about the systems, their functional specificities, their partners and how they are related and (2) a methodological knowledge base (MKB) which depicts the task layout used to identify and reconstruct functional integrated systems. Instantiation of the DKB is obtained by solving the tasks of the MKB, whereas some tasks need instances of the DKB to be solved. AROM, an object-based knowledge representation system, has been used to design the DKB, and its task manager, AROMTasks, for developing the MKB. In this study two integrated systems, ABC transporters and two component systems, both involved in adaptation processes of a bacterial cell to its biotope, have been used to evaluate the feasibility of the approach.

ATP-Binding Cassette Transporters↗

In situ analysis of integrin and growth factor receptor signaling pathways in human glioblastomas suggests overlapping relationships with focal adhesion kinase activation.

Deregulated integrin signaling is common in cancers, including glioblastoma. Integrin binding and growth factor receptor signaling activate focal adhesion kinase (FAK) and subsequently up-regulate extracellular regulated kinases (ERK-1/2), leading to cell-cycle progression and cell migration. Most studies of this pathway have used in vitro systems or tumor lysate-based approaches. We examined these pathways primarily in situ using a panel of 30 glioblastomas and gene expression arrays, immunohistochemistry, and fluorescence in situ hybridization, emphasizing the histological distribution of molecular changes. Within individual tumors, increased expression of FAK, p-FAK, paxillin, ERK-1/2, and p-ERK-1/2 occurred in regions of elevated EGFR and/or PDGFRA expression. Moreover, FAK activation levels correlated with EGFR and PDGFRA expression, and p-FAK and EGFR expression co-localized at the single-cell level. In addition, integrin expression was enriched in EGFR/PDGFRA-overexpressing areas but was more regionally confined than FAK, p-FAK, and paxillin. Integrins beta8 and alpha5beta1 were most commonly expressed, often in a perinecrotic or perivascular pattern. Taken together, our data suggest that growth factor receptor overexpression facilitates alterations in the integrin signaling pathway. Thus, FAK may act in glioblastoma as a downstream target of growth factor signaling, with integrins enhancing the impact of such signaling in the tumor microenvironment.

Enzyme Activation↗

Accumulation of genome-specific transcripts, transcription factors and phytohormonal regulators during early stages of fiber cell development in allotetraploid cotton.

Gene expression during the early stages of fiber cell development and in allopolyploid crops is poorly understood. Here we report computational and expression analyses of 32 789 high-quality ESTs derived from Gossypium hirsutum L. Texas Marker-1 (TM-1) immature ovules (GH_TMO). The ESTs were assembled into 8540 unique sequences including 4036 tentative consensus sequences (TCs) and 4504 singletons, representing approximately 15% of the unique sequences in the cotton EST collection. Compared with approximately 178 000 existing ESTs derived from elongating fibers and non-fiber tissues, GH_TMO ESTs showed a significant increase in the percentage of genes encoding putative transcription factors such as MYB and WRKY and genes encoding predicted proteins involved in auxin, brassinosteroid (BR), gibberellic acid (GA), abscisic acid (ABA) and ethylene signaling pathways. Cotton homologs related to MIXTA, MYB5, GL2 and eight genes in the auxin, BR, GA and ethylene pathways were induced during fiber cell initiation but repressed in the naked seed mutant (N1N1) that is impaired in fiber formation. The data agree with the known roles of MYB and WRKY transcription factors in Arabidopsis leaf trichome development and the well-documented phytohormonal effects on fiber cell development in immature cotton ovules cultured in vitro. Moreover, the phytohormonal pathway-related genes were induced prior to the activation of MYB-like genes, suggesting an important role of phytohormones in cell fate determination. Significantly, AA sub-genome ESTs of all functional classifications including cell-cycle control and transcription factor activity were selectively enriched in G. hirsutum L., an allotetraploid derived from polyploidization between AA and DD genome species, a result consistent with the production of long lint fibers in AA genome species. These results suggest general roles for genome-specific, phytohormonal and transcriptional gene regulation during the early stages of fiber cell development in cotton allopolyploids.

Arabidopsis Proteins↗

Global analysis of carbohydrate utilization by Lactobacillus acidophilus using cDNA microarrays.

The transport and catabolic machinery involved in carbohydrate utilization by Lactobacillus acidophilus was characterized genetically by using whole-genome cDNA microarrays. Global transcriptional profiles were determined for growth on glucose, fructose, sucrose, lactose, galactose, trehalose, raffinose, and fructooligosaccharides. Hybridizations were carried out by using a round-robin design, and microarray data were analyzed with a two-stage mixed model ANOVA. Differentially expressed genes were visualized by hierarchical clustering, volcano plots, and contour plots. Overall, only 63 genes (3% of the genome) showed a >4-fold induction. Specifically, transporters of the phosphoenolpyruvate:sugar transferase system were identified for uptake of glucose, fructose, sucrose, and trehalose, whereas ATP-binding cassette transporters were identified for uptake of raffinose and fructooligosaccharides. A member of the LacS subfamily of galactoside-pentose hexuronide translocators was identified for uptake of galactose and lactose. Saccharolytic enzymes likely involved in the metabolism of monosaccharides, disaccharides, and polysaccharides into substrates of glycolysis were also found, including enzymatic machinery of the Leloir pathway. The transcriptome appeared to be regulated by carbon catabolite repression. Although substrate-specific carbohydrate transporters and hydrolases were regulated at the transcriptional level, genes encoding regulatory proteins CcpA, Hpr, HprK/P, and EI were consistently highly expressed. Genes central to glycolysis were among the most highly expressed in the genome. Collectively, microarray data revealed that coordinated and regulated transcription of genes involved in sugar uptake and metabolism is based on the specific carbohydrate provided. L. acidophilus's adaptability to environmental conditions likely contributes to its competitive ability for limited carbohydrate sources available in the human gastrointestinal tract.

ATP-Binding Cassette Transporters↗

An extensive microarray analysis of AAL-toxin-induced cell death in Arabidopsis thaliana brings new insights into the complexity of programmed cell death in plants.

A T-DNA knockout of the Arabidopsis homologue of the tomato disease resistance gene Asc was obtained. The asc gene renders plants sensitive to programmed cell death (PCD) triggered by the fungal AAL toxin. To obtain more insights into the nature of AAL-toxin-induced cell death and to identify genes of potential importance for PCD, we carried out transcription profiling of AAL-toxin-induced cell death in this knockout with an oligonucleotide array representing 21,500 Arabidopsis genes. Genes responsive to reactive oxygen species (ROS) and ethylene were among the earliest to be upregulated, suggesting that an oxidative burst and production of ethylene played a role in the activation of the cell death. This notion was corroborated by induction of several genes encoding ROS-generating proteins, including a respiratory burst oxidase and germin oxalate oxidase. Cytochemical studies confirmed the oxidative burst and, in addition, showed synthesis of callose, a feature of the hypersensitive response. A diverse group of transcription factors was also induced. These events were followed by repression of most of the auxin-regulated genes known to be involved in growth and developmental responses. All photosynthesis-related genes were repressed. Blocking the synthesis of ethylene or NO significantly compromised cell death. In addition, we identified a heterogeneous group of early-induced genes, some of them never before associated with PCD. The group of early-induced genes included a number of proteases that were previously implicated in developmentally regulated types of PCD, suggesting a more principal role for these proteases in the PCD process. These findings provide new insights into the molecular mechanisms of plant PCD.

Apoptosis↗

Induction of Egr-1 expression by the retinoid AHPN in human lung carcinoma cells is dependent on activated ERK1/2.

The novel retinoid 6-[3-(1-adamantyl)-4-hydroxyphenyl]-2-naphthalene carboxylic acid (AHPN/CD437) inhibits cell proliferation and is a very effective inducer of apoptosis in a variety of carcinoma cell lines. In order to obtain greater insight into the mechanism of AHPN-induced growth arrest and apoptosis, we began to examine AHPN-induced changes in gene expression by cDNA array screening using human lung carcinoma H460 cells. This analysis identified several AHPN-inducible genes, including the immediate-early genes Egr-1 and Nur77. AHPN was able to increase Egr-1 and Nur77 mRNA expression and protein in a variety of carcinoma cell lines. This induction appeared to be regulated at the transcriptional level and was specific for AHPN since an RAR- and an RXR-selective retinoid were inactive. These results suggest that the induction of Egr-1 and Nur77 by AHPN is independent of nuclear retinoid receptors and involves a novel mechanism. Overexpression of Bcl-2, which inhibits AHPN-induced apoptosis but not growth arrest in human T cell lymphoma Molt-4 cells, did not block the induction of immediate-early gene expression. Treatment of H460 cells with AHPN induced activation of the p38 MAP-kinase but not the ERK1/2 signaling pathway. However, inhibition of the ERK1/2 signaling pathway by PD98059 blocked the induction of Egr-1 and Nur77 mRNA while the p38 MAPK inhibitor PD169316 had little effect. Expression of a dominant-negative ERK1 completely abolished the increase in Egr-1 mRNA. Treatment with MAPK inhibitors or expression of dnERK1 reduced but did not block AHPN-induced apoptosis. Our results suggest that the induction of Egr-1 in H460 by AHPN requires active ERK1/2 and is independent of p38 activation. Egr-1, in cooperation with several other growth-suppressor proteins, is likely involved in AHPN-induced inhibition of cell growth and cell death.

Antineoplastic Agents↗

Inhibitor of DNA binding/differentiation helix-loop-helix proteins mediate bone morphogenetic protein-induced osteoblast differentiation of mesenchymal stem cells.

Bone morphogenetic proteins (BMPs) belong to the TGF-beta superfamily and play an important role in development and in many cellular processes. We have found that BMP-2, BMP-6, and BMP-9 induce the most potent osteogenic differentiation of mesenchymal stem cells. Expression profiling analysis has revealed that the Inhibitors of DNA binding/differentiation (Id)-1, Id-2, and Id-3 are among the most significantly up-regulated genes upon BMP-2, BMP-6, or BMP-9 stimulation. Here, we sought to determine the functional role of these Id proteins in BMP-induced osteoblast differentiation. We demonstrated that the expression of Id-1, Id-2, and Id-3 genes was significantly induced at the early stage of BMP-9 stimulation and returned to basal levels at 3 days after stimulation. RNA interference-mediated knockdown of Id expression significantly diminished the BMP-9-induced osteogenic differentiation of mesenchymal progenitor cells. Surprisingly, a constitutive overexpression of these Id genes also inhibited osteoblast differentiation initiated by BMP-9. Furthermore, we demonstrated that BMP-9-regulated Id expression is Smad4-dependent. Overexpression of the three Id genes was shown to promote cell proliferation that was coupled with an inhibition of osteogenic differentiation. Thus, our findings suggest that the Id helix-loop-helix proteins may play an important role in promoting the proliferation of early osteoblast progenitor cells and that Id expression must be down-regulated during the terminal differentiation of committed osteoblasts, suggesting that a balanced regulation of Id expression may be critical to BMP-induced osteoblast lineage-specific differentiation of mesenchymal stem cells.

Adenoviridae↗

Rosiglitazone induction of Insig-1 in white adipose tissue reveals a novel interplay of peroxisome proliferator-activated receptor gamma and sterol regulatory element-binding protein in the regulation of adipogenesis.

Insulin-induced gene 1 (INSIG-1) is a key regulator in the processing of the sterol regulatory element-binding proteins (SREBPs). We demonstrated that Insig-1 is regulated by peroxisome proliferator-activated receptor gamma (PPARgamma) providing a link between insulin sensitization/glucose homeostasis and lipid homeostasis. Insig-1 was identified as a PPARgamma target gene using microarray analysis of mRNA from the white adipose tissue of diabetic (db/db) animals treated with PPARgamma agonists. Insig-1 was induced in subcutaneous (9-fold) and epididymal (4-fold) fat pads from db/db mice treated for 8 days with the PPARgamma agonist rosiglitazone (30 mg/kg/day). This in vivo response was confirmed in differentiated C3H10T1/2 adipocytes treated with rosiglitazone. To elucidate the molecular mechanisms regulating INSIG-1 expression, we cloned and characterized the human INSIG-1 promoter. Co-expression of PPARgamma and RXRalpha transactivated the INSIG-1 promoter in the presence of PPARgamma agonists. This induction was attenuated when a dominant negative PPARgamma construct was transfected into cells. Furthermore, a PPARgamma antagonist repressed the transactivation of the INSIG-1 promoter-reporter construct. Truncations of the promoter resulted in the identification of a PPAR response element that mediated the regulation of the promoter. We demonstrated with recombinant proteins that the PPARgamma/RXRalpha heterodimer binds directly to this PPAR response element. In addition to regulation by PPARgamma/RXRalpha, we demonstrated that the INSIG-1 promoter is regulated by transcriptionally active SREBP. The sterol response element was identified 380 base pairs upstream of the transcriptional start site. These findings suggest that the regulation of Insig-1 by PPARgamma agonists could in turn regulate SREBP processing and thus couple insulin sensitizers with the regulation of lipid homeostasis.

Adipocytes↗

SiMCAL 1 algorithm for analysis of gene expression data related to the phosphatidylserine receptor.

OBJECTIVE: SiMCAL 1 (simple multilevel clustering and linking, version 1) is a novel clustering algorithm for time-series microarray data, presented here with an application to a specific data set. The purpose of the algorithm is to present a complete feature set not found in either Jarvis-Patrick clustering, from which it is derived, or in other popular clustering methods such as hierarchical and k-means. The data concern the activity of the phosphatidylserine receptor (PSR) which is believed to be a crucial molecular switch in the mediation of inflammatory response in apoptosis and lysis. By analyzing the behavior of PSR-related genes in mouse macrophages, we hope to elucidate the mechanisms involved in this important biological process. METHODS AND MATERIALS: SiMCAL 1 is implemented in the Python programming language using the Numerical Python extensions, and the data are stored using the MySQL database management system. The data are derived from exposures of multiple Affymetrix mouse gene microarray chips to elevated levels of PSR antibody and control conditions. Code and data are available at (accessed: 17 January 2005). RESULTS: The algorithm meets its objectives: it is simple, in that it is computationally inexpensive; it is multilevel, in that it provides a small number of clearly defined hierarchical levels of clusters; and it offers linking between clusters at the same level in each hierarchy. Clustering and linking results indicate previously unknown co-regulation for genes expressing PGH synthase (COX2) and PGE2, appear to confirm increased production of proteins for clearance of apoptotic cells in the presence of PSR antibody, and correspond to other findings regarding the temporal relationship between PGE2 production and B cell proliferation and differentiation. These results are promising but should be taken as highly preliminary. CONCLUSION: Both the algorithm and its application to this problem show great potential for future development. We plan to improve and extend the SiMCAL family of algorithms, and to obtain new data so that the algorithm(s) may be further applied to this and other problems of interest.

Algorithms↗

A High-Resolution Stereo-Seq Spatial Transcriptomic Resource for Adult Holstein Cattle Liver.

The bovine liver is a highly compartmentalized organ that plays essential roles in continuous gluconeogenesis and nitrogen recycling; however, its spatial molecular architecture has remained largely uncharacterized due to the limitations of traditional bulk and single-cell approaches. To address this gap, Spatial Enhanced Resolution Omics-sequencing (Stereo-seq) was utilized to generate a subcellular-resolution (500 nm) transcriptomic map of an adult Holstein cattle liver, and a refined reference-guided workflow was implemented to overcome standard annotation limitations in livestock. Raw sequencing data were processed using the Stereo-seq Analysis Workflow and analyzed with Stereopy, Seurat, SingleR, and reference-guided workflows. Spatial aggregation was evaluated at Bin20, Bin50, Bin100, Bin150, and Bin200. Increasing bin size increased molecular identifier counts and detected-gene complexity while progressively reducing spatial granularity. Bin50, corresponding to 50 × 50 DNA nanoballs and an approximate nominal footprint of 25 × 25 µm, was therefore selected as a practical intermediate aggregation level for the primary analyses. Quality-control assessment, Leiden clustering, UMAP visualization, reference-based cell-type annotation, cluster-marker analysis, and spatial mapping of canonical hepatic genes demonstrated preservation of biologically interpretable liver transcriptional organization. Raw sequencing data processed spatial matrices, annotated objects, and analysis code are publicly available to support reanalysis and computational benchmarking. In summary, we present a Stereo-seq spatial transcriptomic resource generated from liver tissue of an adult Holstein cow. This initial resource provides a valuable foundation for future studies of bovine liver biology, comparative genomics, and the spatial basis of livestock health and production traits.

Animals↗

Gene expression analysis of EpiDerm following exposure to SLS using cDNA microarrays.

There is a need to investigate the mechanistic basis of the human skin irritation response if relevant in vitro test systems for the predictive identification of skin irritation hazards are to be developed. Recent progress in genomics technologies mean that tools for the identification and investigation of important biochemical events in the processes of skin irritation are now available. The aim of this work was to identify genes (for further mechanistic investigation) which may be regulated in response to skin irritation, following exposure of the EpiDerm skin model to the known skin irritant sodium lauryl sulphate (SLS). EpiDerm cultures were treated in triplicate with a non-cytotoxic dose of SLS (0.1 mg/ml, as determined by the MTT assay and histological examination) for 15 min, 30 min, 1 h, 2 h, 3 h, 4 h and 24 h. Total RNA was extracted from the pooled EpiDerm cultures and used to probe Atlas human arrays (Clontech) covering approximately 3600 genes. Preliminary data indicated an up-regulation at early time points (15-30 min) of a number of genes involved in transportation (e.g. the sodium and chloride dependent taurine transporter) and receptors (e.g. ZAP70 and protocadherin 42 precursor). The gene encoding the UV excision repair protein and other DNA repair genes (e.g. DNA-directed RNA polymerase II) were up-regulated after 1-3 h, along with TGF beta 3 and other tumour suppressors, which play a role in cellular development and wound healing. At the later time points of 4-24 h, genes involved in protein translation (e.g. Cathepsin D receptor) and metabolism (e.g. CYP27A) were up-regulated. In addition, a number of genes were down-regulated in response to treatment with SLS, although these followed less of a time dependent pattern. These results indicate the differential regulation of a number of genes in response to treatment with SLS, some of which may provide additional clues to the molecular events underpinning the irritation response to this particular surfactant and possibly to other chemical irritants.

Animal Testing Alternatives↗

Identification of mixups among DNA sequencing plates.

MOTIVATION: During the process of high-throughput genome sequencing there are opportunities for mixups of reagents and data associated with particular projects. The sequencing templates or sequence data generated for an assembly may become contaminated with reagents or sequences from another project, resulting in poorer quality and inaccurate assemblies. RESULTS: We have developed a system to assess sequence assemblies and monitor for laboratory mixups. We describe several methods for testing the consistency of assemblies and resolving mixed ones. We use statistical tests to evaluate the distribution of sequencing reads from different plates into contigs, and a graph-based approach to resolve situations where data has been inappropriately combined. While these methods have been designed for use in a high-throughput DNA sequencing environment processing thousands of clones, they can be applied in any situation where distinct sequencing projects are performed at redundant coverage.

Algorithms↗

Molecular cloning and characterization of IL-15R alpha gene in rainbow trout (Oncorhynchus mykiss).

The alpha-chain of IL-15 receptor complex serves as a high-affinity, specific subunit for IL-15 binding. It shares functional and structural features with IL-2 receptor alpha-chain. Here we report for the first time that the molecular cloning, characterization and sequence analysis of the teleostean IL-15R alpha from Oncorhynchus mykiss. The full-length cDNA comprises of a 5' untranslated region (UTR) of 167 bp, an open reading frame of 732 bp and a large 3'UTR of 630 bp, constitutively expressed in all the tissues examined. Another two variants are found derived from alternative splicing. Two copies of rainbow trout IL-15R alpha (rtIL-15R alpha) were detected in the genome by Southern blot hybridization. Moreover, EST evidence is also traced from other fishes. The deduced rtIL-15R alpha of 243 amino acids contains a 17 aa signal peptide at N-terminus and a transmembrane region at C-terminus, as well as a typical sushi domain included in the extracellular region. Phylogenetic tree analysis groups rtIL-15R alpha with other IL-15R alphas but separates it from IL-2R alphas. In addition, a putative sequence was found in the sushi domain which is conserved and versatile among all species and this might relate to cytokine recognition.

Amino Acid Sequence↗

Comparative mapping of sequence-based and structure-based protein domains.

BACKGROUND: Protein domains have long been an ill-defined concept in biology. They are generally described as autonomous folding units with evolutionary and functional independence. Both structure-based and sequence-based domain definitions have been widely used. But whether these types of models alone can capture all essential features of domains is still an open question. METHODS: Here we provide insight on domain definitions through comparative mapping of two domain classification databases, one sequence-based (Pfam) and the other structure-based (SCOP). A mapping score is defined to indicate the significance of the mapping, and the properties of the mapping matrices are studied. RESULTS: The mapping results show a general agreement between the two databases, as well as many interesting areas of disagreement. In the cases of disagreement, the functional and evolutionary characteristics of the domains are examined to determine which domain definition is biologically more informative.

Algorithms↗

A mathematical and computational framework for quantitative comparison and integration of large-scale gene expression data.

Analysis of large-scale gene expression studies usually begins with gene clustering. A ubiquitous problem is that different algorithms applied to the same data inevitably give different results, and the differences are often substantial, involving a quarter or more of the genes analyzed. This raises a series of important but nettlesome questions: How are different clustering results related to each other and to the underlying data structure? Is one clustering objectively superior to another? Which differences, if any, are likely candidates to be biologically important? A systematic and quantitative way to address these questions is needed, together with an effective way to integrate and leverage expression results with other kinds of large-scale data and annotations. We developed a mathematical and computational framework to help quantify, compare, visualize and interactively mine clusterings. We show that by coupling confusion matrices with appropriate metrics (linear assignment and normalized mutual information scores), one can quantify and map differences between clusterings. A version of receiver operator characteristic analysis proved effective for quantifying and visualizing cluster quality and overlap. These methods, plus a flexible library of clustering algorithms, can be called from a new expandable set of software tools called CompClust 1.0 (http://woldlab.caltech.edu/compClust/). CompClust also makes it possible to relate expression clustering patterns to DNA sequence motif occurrences, protein-DNA interaction measurements and various kinds of functional annotations. Test analyses used yeast cell cycle data and revealed data structure not obvious under all algorithms. These results were then integrated with transcription motif and global protein-DNA interaction data to identify G1 regulatory modules.

Algorithms↗