Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Computational prediction”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

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

At least 847 records · Page 47Linked to original sources

Translation initiation of a cardiac voltage-gated potassium channel by internal ribosome entry.

The mammalian Kv1.4 voltage-gated potassium channel mRNA contains an unusually long (1.2 kilobases) 5'-untranslated region (UTR) and includes 18 AUG codons upstream of the authentic site of translation initiation. Computer-predicted secondary structures of this region reveal complex stem-loop structures that would serve as barriers to 5' --> 3' ribosomal scanning. These features suggested that translation initiation in Kv1.4 might occur by the mechanism of internal ribosome entry, a mode of initiation employed by a variety of RNA viruses but only a limited number of vertebrate genes. To test this possibility we introduced the 5'-UTR of mouse Kv1.4 mRNA into the intercistronic region of a bicistronic vector containing two tandem reporter genes, chloramphenicol acetyltransferase and luciferase. The control construct translated only the upstream chloramphenicol cistron in transiently transfected mammalian cells. In contrast, the construct containing the mKv1.4 UTR efficiently translated the luciferase cistron as well, demonstrating the presence of an internal ribosome entry segment. Progressive 5' --> 3' deletions localized the activity to a 3'-proximal 200-nucleotide fragment. Suppression of cap-dependent translation by extracts from poliovirus-infected HeLa cells in an in vitro translation assay eliminated translation of the upstream cistron while allowing translation of the downstream cistron. Our results indicate that the 5'-untranslated region of mKv1.4 contains a functional internal ribosome entry segment that may contribute to unusual and physiologically important modes of translation regulation for this and other potassium channel genes.

Animals↗

Identification of transcription factor binding sites upstream of human genes regulated by the phosphatidylinositol 3-kinase and MEK/ERK signaling pathways.

We have taken an integrated approach in which expression profiling has been combined with the use of small molecule inhibitors and computational analysis of transcription factor binding sites to characterize regulatory sequences of genes that are targets of specific signaling pathways in growth factor-stimulated human cells. T98G cells were stimulated with platelet-derived growth factor (PDGF) and analyzed by DNA microarrays, which identified 74 immediate-early gene transcripts. Cells were then treated with inhibitors to identify subsets of genes that are targets of the phosphatidylinositol 3-kinase (PI3K) and MEK/ERK signaling pathways. Four groups of PDGF-induced genes were defined: independent of PI3K and MEK/ERK signaling, dependent on PI3K signaling, dependent on MEK/ERK signaling, and dependent on both pathways. The upstream regions of all genes in the four groups were scanned using TRANSFAC for putative cis-elements as compared with a background set of non-induced genes. Binding sites for 18 computationally predicted transcription factors were over-represented in the four groups of co-expressed genes compared with the background sequences (p < 0.01). Many of the cis-elements identified were conserved in orthologous mouse genes, and many of the predicted elements and their cognate transcription factors were consistent with previous experimental data. In addition, chromatin immunoprecipitation assays experimentally verified nine predicted SRF binding sites in T98G cells, including a previously unknown SRF site upstream of DUSP5. These results indicate that groups of human genes regulated by discrete intracellular signaling pathways share common cis-regulatory elements.

Animals↗

Detection of protein folding defects caused by BRCA1-BRCT truncation and missense mutations.

Most cancer-associated BRCA1 mutations identified to date result in the premature translational termination of the protein, highlighting a crucial role for the C-terminal, BRCT repeat region in mediating BRCA1 tumor suppressor function. However, the molecular and genetic effects of missense mutations that map to the BRCT region remain largely unknown. Using a protease-based assay, we directly assessed the sensitivity of the folding of the BRCT domain to an extensive set of truncation and single amino acid substitutions derived from breast cancer screening programs. The protein can tolerate truncations of up to 8 amino acids, but further deletion results in drastic BRCT folding defects. This molecular phenotype can be correlated with an increased susceptibility to disease. A cross-validated computational assessment of the BRCT mutation data base suggests that as much as half of all BRCT missense mutations contribute to BRCA1 loss of function and disease through protein-destabilizing effects. The coupled use of proteolytic methods and computational predictive methods to detect mutant BRCA1 conformations at the protein level will augment the efficacy of current BRCA1 screening protocols, especially in the absence of clinical data that can be used to discriminate deleterious BRCT missense mutations from benign polymorphisms.

Amino Acid Sequence↗

Structural features of microRNA (miRNA) precursors and their relevance to miRNA biogenesis and small interfering RNA/short hairpin RNA design.

We have established the structures of 10 human microRNA (miRNA) precursors using biochemical methods. Eight of these structures turned out to be different from those that were computer-predicted. The differences localized in the terminal loop region and at the opposite side of the precursor hairpin stem. We have analyzed the features of these structures from the perspectives of miRNA biogenesis and active strand selection. We demonstrated the different thermodynamic stability profiles for pre-miRNA hairpins harboring miRNAs at their 5'- and 3'-sides and discussed their functional implications. Our results showed that miRNA prediction based on predicted precursor structures may give ambiguous results, and the success rate is significantly higher for the experimentally determined structures. On the other hand, the differences between the predicted and experimentally determined structures did not affect the stability of termini produced through "conceptual dicing." This result confirms the value of thermodynamic analysis based on mfold as a predictor of strand section by RNAi-induced silencing complex (RISC).

Base Sequence↗

Quaternary ammonium compounds as water channel blockers. Specificity, potency, and site of action.

Excessive water uptake through Aquaporins (AQP) can be life-threatening and reversible AQP inhibitors are needed. Here, we determined the specificity, potency, and binding site of tetraethylammonium (TEA) to block Aquaporin water permeability. Using oocytes, externally applied TEA blocked AQP1/AQP2/AQP4 with IC50 values of 1.4, 6.2, and 9.8 microM, respectively. Related tetraammonium compounds yielded some (propyl) or no (methyl, butyl, or pentyl) inhibition. TEA inhibition was lost upon a Tyr to Phe amino acid switch in the external water pore of AQP1/AQP2/AQP4, whereas the water permeability of AQP3 and AQP5, which lack a corresponding Tyr, was not blocked by TEA. Consistent with experimental data, multi-nanosecond molecular dynamics simulations showed one stable binding site for TEA, but not tetramethyl (TMA), in AQP1, resulting in a nearly 50% water permeability inhibition, which was reduced in AQP1-Y186F due to effects on the TEA inhibitory binding region. Moreover, in the simulation TEA interacted with charged residues in the C (Asp128) and E (Asp185) loop, and the A(Tyr37-Asn42-Thr44) loop of the neighboring monomer, but not directly with Tyr186. The loss of TEA inhibition in oocytes expressing properly folded AQP1-N42A or -T44A is in line with the computationally predicted binding mode. Our data reveal that the molecular interaction of TEA with AQP1 differs and is about 1000-fold more effective on AQPs than on potassium channels. Moreover, the observed experimental and simulated similarities open the way for rational design and virtual screening for AQP-specific inhibitors, with quaternary ammonium compounds in general, and TEA in particular as a lead compound.

Amino Acid Sequence↗

Polyamine depletion increases cytoplasmic levels of RNA-binding protein HuR leading to stabilization of nucleophosmin and p53 mRNAs.

Polyamines are essential for maintaining normal intestinal epithelial integrity, an effect that relies, at least in part, on their ability to keep low levels of nucleophosmin (NPM) and p53 mRNAs. The RNA-binding protein HuR associates with the p53 mRNA, as reported previously, and with the NPM mRNA, computationally predicted to be a target of HuR. Here, we show that HuR binds the NPM and p53 3'-untranslated regions and stabilizes these mRNAs in polyamine-depleted intestinal epithelial cells. Depletion of cellular polyamines by inhibiting ornithine decarboxylase with alpha-difluoromethylornithine dramatically enhanced the cytoplasmic abundance of HuR, whereas ectopic ornithine decarboxylase overexpression decreased cytoplasmic HuR; neither intervention changed whole-cell HuR levels. HuR was found to specifically bind the 3'-untranslated regions of NPN and p53 mRNAs. HuR silencing rendered the NPM and p53 mRNAs unstable and prevented increases in NPM and p53 mRNA and protein in polyamine-deficient cells. These results indicate that polyamines modulate cytoplasmic HuR levels in intestinal epithelial cells, in turn controlling the stability of the NPM and p53 mRNAs and influencing NPM and p53 protein levels.

3' Untranslated Regions↗

Iron regulation and the cell cycle: identification of an iron-responsive element in the 3'-untranslated region of human cell division cycle 14A mRNA by a refined microarray-based screening strategy.

Iron regulatory proteins (IRPs) 1 and 2 post-transcriptionally control mammalian iron homeostasis by binding to iron-responsive elements (IREs), conserved RNA stem-loop structures located in the 5'- or 3'-untranslated regions of genes involved in iron metabolism (e.g. FTH1, FTL, and TFRC). To identify novel IRE-containing mRNAs, we integrated biochemical, biocomputational, and microarray-based experimental approaches. IRP/IRE messenger ribonucleoproteins were immunoselected, and their mRNA composition was analyzed using an IronChip microarray enriched for genes predicted computationally to contain IRE-like motifs. Among different candidates, this report focuses on a novel IRE located in the 3'-untranslated region of the cell division cycle 14A mRNA. We show that this IRE motif efficiently binds both IRP1 and IRP2. Differential splicing of cell division cycle 14A produces IRE- and non-IRE-containing mRNA isoforms. Interestingly, only the expression of the IRE-containing mRNA isoforms is selectively increased by cellular iron deficiency. This work describes a new experimental strategy to explore the IRE/IRP regulatory network and uncovers a previously unrecognized regulatory link between iron metabolism and the cell cycle.

3' Untranslated Regions↗

The impact of dyskeratosis congenita mutations on the structure and dynamics of the human telomerase RNA pseudoknot domain.

The pseudoknot domain is a functionally crucial part of telomerase RNA and influences the activity and stability of the ribonucleoprotein complex. Autosomal dominant dyskeratosis congenita (DKC) is an inherited disease that is linked to mutations in telomerase RNA and impairs telomerase function. In this paper, we present a computational prediction of the influence of two base DKC mutations on the structure, dynamics, and stability of the pseudoknot domain. We use molecular dynamics simulations, MM-GBSA free energy calculations, static analysis, and melting simulations analysis. Our results show that the DKC mutations stabilize the hairpin form and destabilize the pseudoknot form of telomerase RNA. Moreover, the P3 region of the predicted DKC-mutated pseudoknot structure is unstable and fails to form as a defined helical stem. We directly compare our predictions with experimental observations by calculating the enthalpy of folding and melting profiles for each structure. The enthalpy values are in very good agreement with values determined by thermal denaturation experiments. The melting simulations and simulations at elevated temperatures show the existence of an intermediate structure, which involves the formation of two UU base pairs observed in the hairpin form of the pseudoknot domain.

Base Sequence↗

HLA and H2 class II transgenic mouse models to study susceptibility and protection in autoimmune thyroid disease.

Using single H2 and HLA class II transgenic mice, in the absence of endogenous H2 class II molecules, we have studied the permissiveness of class II molecules for experimental autoimmune thyroiditis (EAT). Resistant strains expressing susceptible class II molecules, H2Ak or HLA-DR3, developed EAT, clearly demonstrating the importance of class II gene inheritance. Polymorphism for HLA-DRB1 was observed, as DR3, but not DR2 or DR4, molecules were permissive for EAT induction with either mouse (m) or human (h) thyroglobulin (Tg). HLA-DQ polymorphism was also detectable, as hTg-induced EAT developed in DQ8+, but not DQ6+, mice. Class II gene interactions leading to reduced EAT severity were observed in H2 transgenic mice, when H2E transgene was expressed in H2A+ mice or H2A molecules were introduced into our novel H2A- E+ transgenic model. Similarly, in DR3+ mice, only the DQ8 transgene reduced EAT severity, depending on both background genes (C57BL/10 or NOD) and Tg species. Based on computer-predicted, class II-binding motifs, potential pathogenic Tg peptides, either unique to hTg (H2A- E+ model) or shared between mTg and hTg (HLA-DR3+ model), were identified. We have also developed a Graves' disease model by immunizing DR3+ mice with TSH receptor DNA. Thus, transgenic models are excellent tools to study human autoimmune thyroid diseases in the context of murine EAT.

Animals↗

Beyond antigen matching: compatibility intelligence theory for transfusion as an emergent biological system.

BACKGROUND: Despite major advances in serologic testing, extended phenotyping, and blood group genomics, clinically similar transfusion exposures may result in markedly different immune and clinical outcomes. Existing compatibility strategies do not fully explain this biological variability. OBJECTIVES: To examine transfusion compatibility as an emergent donor-recipient biological state and propose a systems-level conceptual framework that integrates established biological determinants into a testable model for future precision transfusion medicine. METHODS: This narrative review critically synthesizes current evidence from blood group genomics, recipient immunobiology, inflammation, disease-specific biology, transfusion medicine, and computational prediction. The proposed framework distinguishes Compatibility Intelligence Theory (CIT) as a biological interpretation from Precision Transfusion Intelligence (PTI) as its potential clinician-supervised translational application. RESULTS: The review argues that transfusion compatibility is shaped by interactions among donor genetics, recipient immune biology, inflammatory physiology, disease context, transfusion history, and longitudinal adaptation rather than by antigen matching alone. CIT provides an organizational framework for integrating these determinants, whereas PTI describes a possible clinician-supervised translation. To address current feasibility, the revised framework separates variables into routinely measurable, contextually available but incompletely standardized, and research-stage domains, and proposes a staged strategy for deriving rather than assuming their quantitative weights. Any clinical implementation would require comparative validation against current serologic, phenotypic, and genotype-based practice. CONCLUSIONS: Compatibility Intelligence Theory offers a testable systems-level framework for understanding transfusion compatibility without replacing established transfusion practices. The framework is not presented as a ready-to-use score: currently measurable variables can be organized for structured risk review, whereas inflammatory, immunogenetic, and multi-omic inputs require prospective standardization and validation. If future studies demonstrate incremental predictive and patient-centered benefit, CIT-informed PTI could support an adaptive, evidence-based extension of current precision transfusion practice.

Humans↗

The initial events in myelin synthesis: orientation of proteolipid protein in the plasma membrane of cultured oligodendrocytes.

Proteolipid protein (PLP) is the most abundant transmembrane protein in myelin of the central nervous system. Conflicting models of PLP topology have been generated by computer predictions based on its primary sequence and experiments with purified myelin. We have examined the initial events in myelin synthesis, including the insertion and orientation of PLP in the plasma membrane, in rat oligodendrocytes which express PLP and the other myelin-specific proteins when cultured without neurons (Dubois-Dalcq, M., T. Behar, L. Hudson, and R. A. Lazzarini. 1986. J. Cell Biol. 102:384-392). These cells, identified by the presence of surface galactocerebroside, the major myelin glycolipid, were stained with six anti-peptide antibodies directed against hydrophilic or short hydrophobic sequences of PLP. Five of these anti-peptide antibodies specifically stained living oligodendrocytes. Staining was only seen approximately 10 d after PLP was first detected in the cytoplasm of fixed and permeabilized cells, suggesting that PLP is slowly transported from the RER to the cell surface. The presence of PLP domains on the extracellular surface was also confirmed by cleavage of such domains with proteases and by antibody-dependent complement-mediated lysis of living oligodendrocytes. Our results indicate that PLP has only two transmembrane domains and that the great majority of the protein, including its amino and carboxy termini, is located on the extracellular face of the oligodendrocyte plasma membrane. This disposition of the PLP molecule suggests that homophilic interactions between PLP molecules of apposed extracellular faces may mediate compaction of adjacent bilayers in the myelin sheath.

Amino Acid Sequence↗

A mathematical model of solute coupled water transport in toad intestine incorporating recirculation of the actively transported solute.

A mathematical model of an absorbing leaky epithelium is developed for analysis of solute coupled water transport. The non-charged driving solute diffuses into cells and is pumped from cells into the lateral intercellular space (lis). All membranes contain water channels with the solute passing those of tight junction and interspace basement membrane by convection-diffusion. With solute permeability of paracellular pathway large relative to paracellular water flow, the paracellular flux ratio of the solute (influx/outflux) is small (2-4) in agreement with experiments. The virtual solute concentration of fluid emerging from lis is then significantly larger than the concentration in lis. Thus, in absence of external driving forces the model generates isotonic transport provided a component of the solute flux emerging downstream lis is taken up by cells through the serosal membrane and pumped back into lis, i.e., the solute would have to be recirculated. With input variables from toad intestine (Nedergaard, S., E.H. Larsen, and H.H. Ussing, J. Membr. Biol. 168:241-251), computations predict that 60-80% of the pumped flux stems from serosal bath in agreement with the experimental estimate of the recirculation flux. Robust solutions are obtained with realistic concentrations and pressures of lis, and with the following features. Rate of fluid absorption is governed by the solute permeability of mucosal membrane. Maximum fluid flow is governed by density of pumps on lis-membranes. Energetic efficiency increases with hydraulic conductance of the pathway carrying water from mucosal solution into lis. Uphill water transport is accomplished, but with high hydraulic conductance of cell membranes strength of transport is obscured by water flow through cells. Anomalous solvent drag occurs when back flux of water through cells exceeds inward water flux between cells. Molecules moving along the paracellular pathway are driven by a translateral flow of water, i.e., the model generates pseudo-solvent drag. The associated flux-ratio equation is derived.

Animals↗

Betatron electron beam characterisation for dosimetry calculations.

Parameters have been specified for electron beams with energies 5--45 MeV from a Brown Boveri betatron for use in computerised dosimetry calculations. A semi-empirical equation is given for the dose at any point at various depths in water. This equation is a modification to Kawachi's predictive model which was based on solutions of a general age--diffusion equation. The depth doses and isodose curves are predicted as a function of the practical range, source skin distance (SSD) and field size. Depth dose accuracy requirements of +/-2% above 50% depth dose and +/-5% at lower doses, relative to maximum dose, have been set and achieved. Also, the shape of the isodose curves with the constrictions at higher doses and bulging at lower values are accurately predicted. Computer calculated beams have been used to generate summed isodose distributions for certain clinical situations.

Computers↗

Wide variation in serum anion gap measurements by chemistry analyzers.

The traditional anion gap [AG = Na-Cl-(total CO2)] mean value of 12 mEq/L was established during the 1970s with analyzer methods that are no longer used widely. No studies have systematically compared mean AG values from analyzers in current use. We used data from healthy subjects obtained from 27 clinical laboratories, 5 manufacturers, and 8 publications to compute mean AG values from 1970s analyzers and 8 current analyzers. We also compared mean AG values by evaluating Na, Cl, and total CO2 data from the College of American Pathologists Chemistry Surveys (1990-1996). Data from healthy subjects showed that overall mean AG values of the 9 analyzers ranged from 5.9 to 12.4 mEq/L. The pooled (i.e., average) AG SD was 2.3 mEq/L. We then used the data of the Surveys and the mean value from 1 analyzer to compute predicted mean values for the other 7 current analyzers. Almost all mean AG values predicted from the Surveys agreed (within 1.5 mEq/L) with mean values from healthy subjects. These results show that mean values of analyzers vary widely, indicating that analytic bias strongly influences the AG. The results should be a useful guide for the AG measurements that can be expected from different analyzers.

Acid-Base Equilibrium↗

Syntactic recognition of regulatory regions in Escherichia coli.

MOTIVATION: One of the most common methodologies to identify cis-regulatory sites in regulatory regions in the DNA is that of weight matrices, as testified by several articles in this issue. An alternative to strengthen the computational predictions in regulatory regions is to develop methods that incorporate more biological properties present in such DNA regions. The grammatical implementation presented in this paper provides a concrete example in this direction. RESULTS: On the basis of the analysis of an exhaustive collection of regulatory regions in Escherichia coli, a grammatical model for the regulatory regions of sigma 70 promoters has been developed. The terminal symbols of the grammar represent individual sites for the binding of activator and repressor proteins, and include the precise position of sites in relation to transcription initiation. Combining these symbols, the grammar generates a large number of different sentences, each of which can be searched for matching against a collection of regulatory regions by means of weight matrices specific for each set of sites for individual proteins. On the basis of this grammatical model, a Prolog syntactic recognizer is presented here. Specific subgrammars for ArgR, LexA and TyrR were implemented. When parsing a collection of 128 sigma 70 promoter regions, the syntactic recognizer produces a much lower number of false-positive sites than the standard search using weight matrices.

Algorithms↗

Comparative interactomics analysis of protein family interaction networks using PSIMAP (protein structural interactome map).

MOTIVATION: Many genomes have been completely sequenced. However, detecting and analyzing their protein-protein interactions by experimental methods such as co-immunoprecipitation, tandem affinity purification and Y2H is not as fast as genome sequencing. Therefore, a computational prediction method based on the known protein structural interactions will be useful to analyze large-scale protein-protein interaction rules within and among complete genomes. RESULTS: We confirmed that all the predicted protein family interactomes (the full set of protein family interactions within a proteome) of 146 species are scale-free networks, and they share a small core network comprising 36 protein families related to indispensable cellular functions. We found two fundamental differences among prokaryotic and eukaryotic interactomes: (1) eukarya had significantly more hub families than archaea and bacteria and (2) certain special hub families determined the topology of the eukaryotic interactomes. Our comparative analysis suggests that a very small number of expansive protein families led to the evolution of interactomes and seemed to have played a key role in species diversification. SUPPLEMENTARY INFORMATION: http://interactomics.org.

Algorithms↗

Elevated intron retention implicates neuroinflammation in brains of individuals with alcohol use disorder.

Intron retention, a form of alternative RNA splicing, can occur as part of normal gene regulation or result from disruption of the splicing machinery. Retained introns can potentially form double-stranded RNA, activating innate immune sensors and inflammation. This mechanism has been implicated in cancer but has not been studied in neuropsychiatric diseases like alcohol use disorder. We systematically analysed transcriptome-wide intron retention events in post-mortem brain tissue from 142 individuals (66 with alcohol use disorder and 76 controls), encompassing 320 region-specific samples from the superior frontal cortex, nucleus accumbens, central nucleus and basolateral amygdala. Analyses were adjusted for demographic, technical and biological covariates. Validation was performed in alcohol-preferring (P) rats using long-read sequencing. In complementary experiments, immunofluorescent staining was used to detect double-stranded RNA in rat brain tissue, while single-cell RNA-sequencing was performed to test activation of double-stranded RNA-sensing pathways in human brains. Brains from individuals with alcohol use disorder showed significantly higher total intron retention compared with controls, independent of age, with females showing greater increases than males. A total of 368 introns were positively associated with alcohol use disorder, and these introns were significantly longer and had weaker splice acceptor sites compared with non-associated introns. Genes harbouring these intron retention events were enriched in Purkinje neurons, visual cortex neurons and oligodendrocytes. Computational predictions indicated these long introns could form duplex RNA structures. Increased double-stranded RNA was confirmed experimentally in multiple brain regions of alcohol-consuming rats, where it co-localized primarily with neuronal nuclei and dendrites. In individuals with alcohol use disorder, we found that multiple pathways including double-stranded RNA responses, neuroinflammation, interferon and NF-&#x3ba;B signalling, adaptive immunity and apoptosis were activated. In addition, NeuN-positive neuronal counts significantly decreased in both the prefrontal and visual cortices. Furthermore, single-cell analysis demonstrated upregulation of TICAM1, the target of double-stranded RNA sensor TLR3, in oligodendrocytes, as well as widespread activation of downstream inflammatory pathways across glial and neuronal cell types. These findings provide the first evidence that chronic alcohol consumption promotes an overall increase of intron retention in the brain and is associated with the presence of double-stranded RNA. Furthermore, the double-stranded RNA may contribute to neuronal loss and brain pathology by activating a neuroinflammatory response.

alcohol use disorder↗

Complete genome structure of the thermophilic cyanobacterium Thermosynechococcus elongatus BP-1.

The entire genome of a thermophilic unicellular cyanobacterium, Thermosynechococcus elongatus BP-1, was sequenced. The genome consisted of a circular chromosome 2,593,857 bp long, and no plasmid was detected. A total of 2475 potential protein-encoding genes, one set of rRNA genes, 42 tRNA genes representing 42 tRNA species and 4 genes for small structural RNAs were assigned to the chromosome by similarity search and computer prediction. The translated products of 56% of the potential protein-encoding genes showed sequence similarity to experimentally identified and predicted proteins of known function, and the products of 34% of these genes showed sequence similarity to the translated products of hypothetical genes. The remaining 10% lacked significant similarity to genes for predicted proteins in the public DNA databases. Sixty-three percent of the T. elongatus genes showed significant sequence similarity to those of both Synechocystis sp. PCC 6803 and Anabaena sp. PCC 7120, while 22% of the genes were unique to this species, indicating a high degree of divergence of the gene information among cyanobacterial strains. The lack of genes for typical fatty acid desaturases and the presence of more genes for heat-shock proteins in comparison with other mesophilic cyanobacteria may be genomic features of thermophilic strains. A remarkable feature of the genome is the presence of 28 copies of group II introns, 8 of which contained a presumptive gene for maturase/reverse transcriptase. A trace of genome rearrangement mediated by the group II introns was also observed.

Cyanobacteria↗