Transcriptome Analysis Reveals Novel Genes Associated with Cartilage Degeneration in Posttraumatic Osteoarthritis Progression.

OBJECTIVE: The current therapeutic strategy for posttraumatic osteoarthritis (PTOA) focuses on early intervention to attenuate disease progression, preserve joint function, and defer joint replacement timing. Sequential transcriptomic changes of articular cartilage in a rat model were investigated to explore the molecular mechanism in early PTOA progression. DESIGN: Anterior cruciate ligament transection and medial meniscectomy (ACLT + MMx)-induced PTOA model was applied on male Wistar rats. Articular cartilages were harvested at time 0 (naïve), 2 week, and 4 weeks after surgery. Affymetrix Rat genome 230 2.0 array was utilized to analyze the gene expression changes of articular cartilages. RESULTS: We identified 849 differentially expressed genes (DEGs) at 2 weeks and 223 DEGs at 4 weeks post-ACLT + MMx surgery compared with time 0 (naïve group). Gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed to gain further insights from these DEGs. 22 novel genes and 1 novel KEGG pathway (axon guidance) in cartilage degeneration of osteoarthritis were identified. Axon guidance molecules-Gnai1, Sema4d, Plxnb1, and Srgap2 commonly dysregulated in PTOA progression. Gnai1 gene showed a concordant change in protein expression by immunohistochemistry staining. CONCLUSIONS: Our study identified 22 novel dysregulated genes and axon guidance pathway associated with articular cartilage degeneration in PTOA progression. These findings provide the potential candidates of biomarkers and therapeutic targets for further investigation.

Comparative functional genomic analysis of Alzheimer's affected and naturally aging brains.

BACKGROUND: Alzheimer's disease (AD) is a prevalent progressive neurodegenerative human disease whose cause remains unclear. Numerous initially highly hopeful anti-AD drugs based on the amyloid-ß (Aß) hypothesis of AD have failed recent late-phase tests. Natural aging (AG) is a high-risk factor for AD. Here, we aim to gain insights in AD that may lead to its novel therapeutic treatment through conducting meta-analyses of gene expression microarray data from AG and AD-affected brain. METHODS: Five sets of gene expression microarray data from different regions of AD (hereafter, ALZ when referring to data)-affected brain, and one set from AG, were analyzed by means of the application of the methods of differentially expressed genes and differentially co-expressed gene pairs for the identification of putatively disrupted biological pathways and associated abnormal molecular contents. RESULTS: Brain-region specificity among ALZ cases and AG-ALZ differences in gene expression and in KEGG pathway disruption were identified. Strong heterogeneity in AD signatures among the five brain regions was observed: HC/PC/SFG showed clear and pronounced AD signatures, MTG moderately so, and EC showed essentially none. There were stark differences between ALZ and AG. OXPHOS and Proteasome were the most disrupted pathways in HC/PC/SFG, while AG showed no OXPHOS disruption and relatively weak Proteasome disruption in AG. Metabolic related pathways including TCA cycle and Pyruvate metabolism were disrupted in ALZ but not in AG. Three pathogenic infection related pathways were disrupted in ALZ. Many cancer and signaling related pathways were shown to be disrupted AG but far less so in ALZ, and not at all in HC. We identified 54 "ALZ-only" differentially expressed genes, all down-regulated and which, when used to augment the gene list of the KEGG AD pathway, made it significantly more AD-specific.

miRDRN-miRNA disease regulatory network: a tool for exploring disease and tissue-specific microRNA regulatory networks.

BACKGROUND: MicroRNA (miRNA) regulates cellular processes by acting on specific target genes, and cellular processes proceed through multiple interactions often organized into pathways among genes and gene products. Hundreds of miRNAs and their target genes have been identified, as are many miRNA-disease associations. These, together with huge amounts of data on gene annotation, biological pathways, and protein-protein interactions are available in public databases. Here, using such data we built a database and web service platform, miRNA disease regulatory network (miRDRN), for users to construct disease and tissue-specific miRNA-protein regulatory networks, with which they may explore disease related molecular and pathway associations, or find new ones, and possibly discover new modes of drug action. METHODS: Data on disease-miRNA association, miRNA-target association and validation, gene-tissue association, gene-tumor association, biological pathways, human protein interaction, gene ID, gene ontology, gene annotation, and product were collected from publicly available databases and integrated. A large set of miRNA target-specific regulatory sub-pathways (RSPs) having the form (T, G 1, G 2) was built from the integrated data and stored, where T is a miRNA-associated target gene, G 1 (G 2) is a gene/protein interacting with T (G 1). Each sequence (T, G 1, G 2) was assigned a p-value weighted by the participation of the three genes in molecular interactions and reaction pathways. RESULTS: A web service platform, miRDRN (http://mirdrn.ncu.edu.tw/mirdrn/), was built. The database part of miRDRN currently stores 6,973,875 p-valued RSPs associated with 116 diseases in 78 tissue types built from 207 diseases-associated miRNA regulating 389 genes. miRDRN also provides facilities for the user to construct disease and tissue-specific miRNA regulatory networks from RSPs it stores, and to download and/or visualize parts or all of the product. User may use miRDRN to explore a single disease, or a disease-pair to gain insights on comorbidity. As demonstrations, miRDRN was applied: to explore the single disease colorectal cancer (CRC), in which 26 novel potential CRC target genes were identified; to study the comorbidity of the disease-pair Alzheimer's disease-Type 2 diabetes, in which 18 novel potential comorbid genes were identified; and, to explore possible causes that may shed light on recent failures of late-phase trials of anti-AD, BACE1 inhibitor drugs, in which genes downstream to BACE1 whose suppression may affect signal transduction were identified.

Genotypes of cancer stem cells characterized by epithelial-to-mesenchymal transition and proliferation related functions.

Cancer stem cells (CSCs), or cancer cells with stem cell-like properties, generally exhibit drug resistance and have highly potent cancer inducing capabilities. Genome-wide expression data collected at public repositories over the last few years provide excellent material for studies that can lead to insights concerning the molecular and functional characteristics of CSCs. Here, we conducted functional genomic studies of CSC based on fourteen PCA-screened high quality public CSC whole genome gene expression datasets and, as control, four high quality non-stem-like cancer cell and non-cancerous stem cell datasets from the Gene Expression Omnibus database. A total of 6,002 molecular signatures were taken from the Molecular Signatures Database and used to characterize the datasets, which, under two-way hierarchical clustering, formed three genotypes. Type 1, consisting of mainly glia CSCs, had significantly enhanced proliferation, and significantly suppressed epithelial-mesenchymal transition (EMT), related functions. Type 2, mainly breast CSCs, had significantly enhanced EMT, but not proliferation, related functions. Type 3, composed of ovarian, prostate, and colon CSCs, had significantly suppressed proliferation related functions and mixed expressions on EMT related functions.

An in vivo molecular response analysis of colorectal cancer treated with Astragalus membranaceus extract.

The fact that many chemotherapeutic drugs cause chemoresistance and side effects during the course of colorectal cancer treatment necessitates development of novel cytotoxic agents aiming to attenuate new molecular targets. Here, we show that Astragalus membranaceus (Fischer) Bge. var. mongolicus (Bge.) Hsiao (AM), a traditional Chinese medicine, can inhibit tumor growth in vivo and elucidate the underlying molecular mechanisms. The antitumor effect of AM was assessed on the subcutaneous tumors of human colorectal cancer cell line HCT116 grafted into nude mice. The mice were treated with either water or 500 mg/kg AM once per day, before being sacrificed for extraction of tumors, which were then subjected to microarray expression profiling. The gene expression of the extraction was then profiled using microarray analysis. The identified genes differentially expressed between treated mice and controls reveal that administration of AM suppresses chromosome organization, histone modification, and regulation of macromolecule metabolic process. A separate analysis focused on differentially expressed microRNAs revealing involvement of macromolecule metabolism, and intracellular transport, as well as several cancer signaling pathways. For validation, the input of the identified genes to The Library of Integrated Network-based Cellular Signatures led to many chemopreventive agents of natural origin that produce similar gene expression profiles to that of AM. The demonstrated effectiveness of AM suggests a potential therapeutic drug for colorectal cancer.

Non-Catalytic RISCs and Kinetics Determine Mammalian siRNA Sub-Cellular Localization.

Small interfering RNAs (siRNAs) are fundamental to the regulation of cell function. Much is known about its gene interfering mechanism, but a kinetic description of it is still lacking. Here, we derived a set of reaction-diffusion equations for multiple RNA-induced silencing complex (RISC) pathways that give quantitative temporal and spatial descriptions of the siRNA process in mammalian cell, and are able to correctly describe all salient experimentally observed patterns of sub-cellular siRNA localization, including those that, at first glance, appear irreconcilable. These results suggest siRNA sub-cellular localization mainly concerns the non-catalytic RISC-target complex, and is caused by the selectiveness of RISC-target interaction and the permeability of the nuclear membrane to siRNA strands but not to RISC-target complexes. Our method is expected to be useful in devising RNAi based cell regulation strategies.

Gene-Set Local Hierarchical Clustering (GSLHC)--A Gene Set-Based Approach for Characterizing Bioactive Compounds in Terms of Biological Functional Groups.

Gene-set-based analysis (GSA), which uses the relative importance of functional gene-sets, or molecular signatures, as units for analysis of genome-wide gene expression data, has exhibited major advantages with respect to greater accuracy, robustness, and biological relevance, over individual gene analysis (IGA), which uses log-ratios of individual genes for analysis. Yet IGA remains the dominant mode of analysis of gene expression data. The Connectivity Map (CMap), an extensive database on genomic profiles of effects of drugs and small molecules and widely used for studies related to repurposed drug discovery, has been mostly employed in IGA mode. Here, we constructed a GSA-based version of CMap, Gene-Set Connectivity Map (GSCMap), in which all the genomic profiles in CMap are converted, using gene-sets from the Molecular Signatures Database, to functional profiles. We showed that GSCMap essentially eliminated cell-type dependence, a weakness of CMap in IGA mode, and yielded significantly better performance on sample clustering and drug-target association. As a first application of GSCMap we constructed the platform Gene-Set Local Hierarchical Clustering (GSLHC) for discovering insights on coordinated actions of biological functions and facilitating classification of heterogeneous subtypes on drug-driven responses. GSLHC was shown to tightly clustered drugs of known similar properties. We used GSLHC to identify the therapeutic properties and putative targets of 18 compounds of previously unknown characteristics listed in CMap, eight of which suggest anti-cancer activities. The GSLHC website http://cloudr.ncu.edu.tw/gslhc/ contains 1,857 local hierarchical clusters accessible by querying 555 of the 1,309 drugs and small molecules listed in CMap. We expect GSCMap and GSLHC to be widely useful in providing new insights in the biological effect of bioactive compounds, in drug repurposing, and in function-based classification of complex diseases.

Hylan G-F 20 attenuates posttraumatic osteoarthritis progression: Association with upregulated expression of the circardian gene NPAS2.

AIMS: The study was to examine the effect of Hylan G-F 20 on the progression of posttraumatic osteoarthritis (PTOA) and the expression of the circadian genes neuronal PAS domain protein 2 (NPAS2) and period 2 (Per2). MAIN METHODS: We used the anterior cruciate ligament transaction and medial menisectomy (ACLT+MMx) model in Wistar rats. The rats were divided into three groups, the sham-operated group, the Hylan G-F 20-treated group, and the saline-treated group. Rats which underwent ACLT + MMx surgery were injected intraarticularly with, respectively, Hylan G-F 20 or saline once a week for 3 consecutive weeks, starting 7days after surgery. The gross morphology and histopathology of the experimental knee joints were evaluated at the end of week 6. Expression of the NPAS2 and Per2 genes was measured by real-time PCR. KEY FINDINGS: Hylan G-F 20 suppressed the articular cartilage destruction and synovitis compared to the saline-treated group. Compared to the sham-operated group, the Hylan G-F 20-treated group showed significantly upregulated expression of NPAS2 in cartilage (2.53±0.08-fold higher; p<0.05) and a non-significant increase in Per2 expression (2.35±1.26-fold higher p=0.28), while the saline-treated group showed significant downregulation of NPAS2 expression and a non-significant decrease in Per2 expression. SIGNIFICANCE: Our data suggested that early intraarticular injection of Hylan G-F 20 attenuates the progression of PTOA and significantly upregulates NPAS2 expression. These findings provide a new direction for studying associations between the use of a pharmacological agent, the degenerative process, and circadian gene expression.

Zfra activates memory Hyal-2+ CD3- CD19- spleen cells to block cancer growth, stemness, and metastasis in vivo.

Zfra is a 31-amino-acid zinc finger-like protein, which participates in the tumor necrosis factor signaling. Here, we determined that when nude mice and BALB/c mice were pre-injected with nanogram levels of a synthetic Zfra1-31 or truncated Zfra4-10 peptide via tail veins, these mice became resistant to the growth, metastasis and stemness of melanoma cells, and many malignant cancer cells. The synthetic peptides underwent self-polymerization in phosphate-buffered saline. Alteration of the Ser8 phosphorylation site to Gly8 abolished Zfra aggregation and its-mediated cancer suppression in vivo. Injected Zfra peptide autofluoresced due to polymerization and was trapped mainly in the spleen. Transfer of Zfra-stimulated spleen cells to naïve mice conferred resistance to cancer growth. Zfra-binding cells, designated Hyal-2+ CD3- CD19- Z cells, are approximately 25-30% in the normal spleen, but are significantly downregulated (near 0-3%) in tumor-growing mice. Zfra prevented the loss of Z cells caused by tumors. In vitro stimulation or education of naïve spleen cells with Zfra allowed generation of activated Z cells to confer a memory anticancer response in naïve or cancer-growing mice. In particular, Z cells are abundant in nude and NOD-SCID mice, and can be readily activated by Zfra to mount against cancer growth.

miR-524-5p suppresses the growth of oncogenic BRAF melanoma by targeting BRAF and ERK2.

It has been well documented that miRNAs can modulate the effectiveness of cancer-associated signaling pathways. Mitogen-activated protein kinase (MAPK/ERK) signaling plays an essential role in the progression of many cancers, including melanoma and colon cancers. However, no single miRNA is reported to directly target multiple components of the MAPK/ERK pathway. We performed a miRNA PCR array screening with various MAPK/ERK signaling activities. The miRNA array data revealed that the expression of miR-524-5p was decreased in cells with an active MAPK/ERK pathway and confirmed that the expression of miR-524-5p is inversely associated with the activity of the MAPK/ERK pathway. We demonstrated that miR-524-5p directly binds to the 3'-untranslated regions of both BRAFandERK2 and suppresses the expression of these proteins. Because BRAF and ERK2 are the main components of MAPK signaling, the overexpression of miR-524-5p effectively inhibits MAPK/ERK signaling, tumor proliferation, and melanoma cell migration. Moreover, tumors overexpressing miR-524-5p were significantly smaller than those of the negative control mice. Our findings provide new insight into the role of miR-524-5p as an important miRNA that negatively regulates the MAPK/ERK signaling pathway, suggesting that miR-524-5p could be a potent therapeutic candidate for melanoma treatment.

Functional Module Connectivity Map (FMCM): a framework for searching repurposed drug compounds for systems treatment of cancer and an application to colorectal adenocarcinoma.

Drug repurposing has become an increasingly attractive approach to drug development owing to the ever-growing cost of new drug discovery and frequent withdrawal of successful drugs caused by side effect issues. Here, we devised Functional Module Connectivity Map (FMCM) for the discovery of repurposed drug compounds for systems treatment of complex diseases, and applied it to colorectal adenocarcinoma. FMCM used multiple functional gene modules to query the Connectivity Map (CMap). The functional modules were built around hub genes identified, through a gene selection by trend-of-disease-progression (GSToP) procedure, from condition-specific gene-gene interaction networks constructed from sets of cohort gene expression microarrays. The candidate drug compounds were restricted to drugs exhibiting predicted minimal intracellular harmful side effects. We tested FMCM against the common practice of selecting drugs using a genomic signature represented by a single set of individual genes to query CMap (IGCM), and found FMCM to have higher robustness, accuracy, specificity, and reproducibility in identifying known anti-cancer agents. Among the 46 drug candidates selected by FMCM for colorectal adenocarcinoma treatment, 65% had literature support for association with anti-cancer activities, and 60% of the drugs predicted to have harmful effects on cancer had been reported to be associated with carcinogens/immune suppressors. Compounds were formed from the selected drug candidates where in each compound the component drugs collectively were beneficial to all the functional modules while no single component drug was harmful to any of the modules. In cell viability tests, we identified four candidate drugs: GW-8510, etacrynic acid, ginkgolide A, and 6-azathymine, as having high inhibitory activities against cancer cells. Through microarray experiments we confirmed the novel functional links predicted for three candidate drugs: phenoxybenzamine (broad effects), GW-8510 (cell cycle), and imipenem (immune system). We believe FMCM can be usefully applied to repurposed drug discovery for systems treatment of other types of cancer and other complex diseases.

Cystatin C in cerebrospinal fluid is upregulated in elderly patients with chronic osteoarthritis pain and modulated through matrix metalloproteinase 9-specific pathway.

OBJECTIVES: Although high abundant cystatin c (CysC) in cerebrospinal fluid (CSF) is well known, its ambiguous role associated with pain still remains unclear. This study evaluated the effects of intrathecal CysC content from chronic pain caused by osteoarthritis (OA) and the novel relationship with matrix metalloproteinases 2 and 9 (MMP2 and MMP9) in CSF. METHODS: Samples of CSF were obtained from 8 elderly patients (65 y and above) with OA with lower limb pain for at least 6 months (OA group) and 8 sex-matched and age-matched relatively healthy elderly individuals without any pain problems (control group). The intrathecal CysC, MMP2, and MMP9 were examined by Western blotting. The analysis of CysC cleavage under different conditions was performed through silver staining and using mass-spectroscopy (SELDI-TOF) on 2 groups. RESULTS: Expression of full-length CysC and pro-MMP2 proteins showed statistically significant upregulation (P=0.0004 vs. 0.03), and expression of MMP9 protein showed downregulation (P=0.007) in the OA group. Both MMP9 and MMP2 initiated the mechanism for full-length CysC cleavage only in the presence of CSF. However, MMP9 showed greater ability than MMP2 for CysC cleavage in control and OA groups in sliver staining. Incubation of CSF with the MMP9 inhibitor led to the suppression of CysC cleavage in SELDI-TOF. DISCUSSION: These findings provide the first in vivo evidence on a relationship between CysC and gelatinases (MMP2 and MMP9), and could facilitate further investigation of novel interactions among these proteins within the proteomics field, especially protein-protein interactions involved in pain.

ToP: a trend-of-disease-progression procedure works well for identifying cancer genes from multi-state cohort gene expression data for human colorectal cancer.

Significantly expressed genes extracted from microarray gene expression data have proved very useful for identifying genetic biomarkers of diseases, including cancer. However, deriving a disease related inference from a list of differentially expressed genes has proven less than straightforward. In a systems disease such as cancer, how genes interact with each other should matter just as much as the level of gene expression. Here, in a novel approach, we used the network and disease progression properties of individual genes in state-specific gene-gene interaction networks (GGINs) to select cancer genes for human colorectal cancer (CRC) and obtain a much higher hit rate of known cancer genes when compared with methods not based on network theory. We constructed GGINs by integrating gene expression microarray data from multiple states--healthy control (Nor), adenoma (Ade), inflammatory bowel disease (IBD) and CRC--with protein-protein interaction database and Gene Ontology. We tracked changes in the network degrees and clustering coefficients of individual genes in the GGINs as the disease state changed from one to another. From these we inferred the state sequences Nor-Ade-CRC and Nor-IBD-CRC both exhibited a trend of (disease) progression (ToP) toward CRC, and devised a ToP procedure for selecting cancer genes for CRC. Of the 141 candidates selected using ToP, ∼50% had literature support as cancer genes, compared to hit rates of 20% to 30% for standard methods using only gene expression data. Among the 16 candidate cancer genes that encoded transcription factors, 13 were known to be tumorigenic and three were novel: CDK1, SNRPF, and ILF2. We identified 13 of the 141 predicted cancer genes as candidate markers for early detection of CRC, 11 and 2 at the Ade and IBD states, respectively.

A fast universal immobilization of immunoglobulin G at 4 °C for the development of array-based immunoassays.

To maintain the antibody activity and enhance performance of array-based immunoassays, protein G was used to allow a shorter duration of immunoglobulin G immobilization at 4 °C, with the antibody placed in the appropriate orientation. The multiplexed detection of six pain-related message molecules (PRMMs) was used as examples for the development of array-based immunoassays: substance P, calcitonin gene-related peptide, nerve growth factor, brain-derived neurotrophic factor, tumor necrosis factor-α, and ß-endorphin. Protein G- and non-protein G-coated slides were tested. Compared to non-protein G immunoassays, protein G shortened the antibody immobilization time at 4 °C from overnight to 2 hours. Only protein G-facilitated immunoassays succeeded in simultaneously detecting all six PRMMs with high specificity. Dose-response curves showed that the limits of detection of the protein G-multiplexed immunoassays for the PRMMs was approximately 164, 167, 120, 60, 80, and 92 pg/ml, respectively. Thus, protein G effectively shortens the duration of antibody immobilization at 4 °C, allowing the use of sensitive array-based immunoassays for the simultaneous detection of PRMMs.

Universal global imprints of genome growth and evolution--equivalent length and cumulative mutation density.

BACKGROUND: Segmental duplication is widely held to be an important mode of genome growth and evolution. Yet how this would affect the global structure of genomes has been little discussed. METHODS/PRINCIPAL FINDINGS: Here, we show that equivalent length, or L(e), a quantity determined by the variance of fluctuating part of the distribution of the k-mer frequencies in a genome, characterizes the latter's global structure. We computed the L(e)s of 865 complete chromosomes and found that they have nearly universal but (k-dependent) values. The differences among the L(e) of a chromosome and those of its coding and non-coding parts were found to be slight. CONCLUSIONS: We verified that these non-trivial results are natural consequences of a genome growth model characterized by random segmental duplication and random point mutation, but not of any model whose dominant growth mechanism is not segmental duplication. Our study also indicates that genomes have a nearly universal cumulative "point" mutation density of about 0.73 mutations per site that is compatible with the relatively low mutation rates of (1-5) x 10(-3)/site/Mya previously determined by sequence comparison for the human and E. coli genomes.

Inverse symmetry in complete genomes and whole-genome inverse duplication.

The cause of symmetry is usually subtle, and its study often leads to a deeper understanding of the bearer of the symmetry. To gain insight into the dynamics driving the growth and evolution of genomes, we conducted a comprehensive study of textual symmetries in 786 complete chromosomes. We focused on symmetry based on our belief that, in spite of their extreme diversity, genomes must share common dynamical principles and mechanisms that drive their growth and evolution, and that the most robust footprints of such dynamics are symmetry related. We found that while complement and reverse symmetries are essentially absent in genomic sequences, inverse-complement plus reverse-symmetry is prevalent in complex patterns in most chromosomes, a vast majority of which have near maximum global inverse symmetry. We also discovered relations that can quantitatively account for the long observed but unexplained phenomenon of -mer skews in genomes. Our results suggest segmental and whole-genome inverse duplications are important mechanisms in genome growth and evolution, probably because they are efficient means by which the genome can exploit its double-stranded structure to enrich its code-inventory.

Layer IV of the primary somatosensory cortex has the highest complexity under anesthesia and cortical complexity is modulated by specific thalamic inputs.

The system complexity, as calculated from correlation dimension, embedded in each layer and its modulation by specific inputs and general excitatory state are not yet known. The aims of present study were to estimate the system complexity across the cortical layers by analyzing intracortical EEG signals using a nonlinear analytical method, and to identify how layer-related complexity varies with the alteration of thalamic input and brain state. Male Sprague-Dawley rats were anesthetized under l% halothane. Sixteen channels of evoked or spontaneous EEG signals were recorded simultaneously across the six cortical layers in the somatosensory cortex with a single Michigan probe. The system complexity was assessed by computing correlation dimension, D(2), based on the Nonlinear Time Series Analysis data analysis program. Cortical layer IV exhibited a D(2) value, 3.24, that was significantly higher than that of the other cortical layers. The D(2) values in layers IV and II/III were significantly reduced after reversible deactivation of the ventral posterior lateral thalamic nucleus. D(2) decreased with increases in administered halothane concentration from 0.75% to 2.0%, particularly in layer IV. The present findings suggest that cortical layer IV maintains a higher complexity than the other layers and that the complexity of the mid-cortical layers is subject to regulation from specific thalamic inputs and more sensitive to changes in the general state of brain excitation.

Evaluation of the current models for the evolution of bacterial DNA uptake signal sequences.

Current opinion considers two main hypotheses for the evolutionary origin of uptake signal sequences in bacteria: one model regards the uptake signal sequence (USS) as the result of biased gene conversion, whereas the second model views the USS as a molecular tag that evolved as an adaptation. In this article, we present various computational models that implement specific versions of those hypotheses. Those models show that the two hypothesis are not necessarily as opposed to each other as may appear at first glance.

Evolution of DNA uptake signal sequences.

The DNA of some naturally competent species of bacteria contains a large number of evenly distributed copies of a short sequence. This highly overrepresented sequence is believed to be an uptake signal sequence (USS) that helps bacteria to take up DNA selectively from (dead) members of their own species. For some time it has been assumed that the USS evolved in order to enable bacteria to distinguish between conspecific and nonconspecific DNA fragments (the preference-first hypothesis). Recently, Redfield suggested that this hypothesis is not in fact realistic, as it would require biologically implausible group selection. In this article we present a model designed to demonstrate the emergence of similar USSs in a population of simulated evolving agents. We use this model to examine the conditions under which a USS will emerge in a preference-first scenario.

Shannon information in complete genomes.

Shannon information in the genomes of all completely sequenced prokaryotes and eukaryotes are measured in word lengths of two to ten letters. It is found that in a scale-dependent way, the Shannon information in complete genomes are much greater than that in matching random sequences--thousands of times greater in the case of short words. Furthermore, with the exception of the 14 chromosomes of Plasmodium falciparum, the Shannon information in all available complete genomes belong to a universality class given by an extremely simple formula. The data are consistent with a model for genome growth composed of two main ingredients: random segmental duplications that increase the Shannon information in a scale-independent way, and random point mutations that preferentially reduces the larger-scale Shannon information. The inference drawn from the present study is that the large-scale and coarse-grained growth of genomes was selectively neutral and this suggests an independent corroboration of Kimura's neutral theory of evolution.