Relationships between genomic dissipation and de novo SNP evolution.

Patterns of single nucleotide polymorphisms (SNPs) in eukaryotic DNA are traditionally attributed to selective pressure, drift, identity descent, or related factors-without accounting for ways in which bias during de novo SNP formation, itself, might contribute. A functional and phenotypic analysis based on evolutionary resilience of DNA points to decreased numbers of non-synonymous SNPs in human and other genomes, with a predominant component of SNP depletion in the human gene pool caused by robust preferences during de novo SNP formation (rather than selective constraint). Ramifications of these findings are broad, belie a number of concepts regarding human evolution, and point to a novel interpretation of evolving DNA across diverse species.

SNP Formation Bias in the Murine Genome Provides Evidence for Parallel Evolution.

In this study, we show novel DNA motifs that promote single nucleotide polymorphism (SNP) formation and are conserved among exons, introns, and intergenic DNA from mice (Sanger Mouse Genomes Project), human genes (1000 Genomes), and tumor-specific somatic mutations (data from TCGA). We further characterize SNPs likely to be very recent in origin (i.e., formed in otherwise congenic mice) and show enrichment for both synonymous and parallel DNA variants occurring under circumstances not attributable to purifying selection. The findings provide insight regarding SNP contextual bias and eukaryotic codon usage as strategies that favor long-term exonic stability. The study also furnishes new information concerning rates of murine genomic evolution and features of DNA mutagenesis (at the time of SNP formation) that should be viewed as "adaptive."

Longevity and plasticity of CFTR provide an argument for noncanonical SNP organization in hominid DNA.

Like many other ancient genes, the cystic fibrosis transmembrane conductance regulator (CFTR) has survived for hundreds of millions of years. In this report, we consider whether such prodigious longevity of an individual gene--as opposed to an entire genome or species--should be considered surprising in the face of eons of relentless DNA replication errors, mutagenesis, and other causes of sequence polymorphism. The conventions that modern human SNP patterns result either from purifying selection or random (neutral) drift were not well supported, since extant models account rather poorly for the known plasticity and function (or the established SNP distributions) found in a multitude of genes such as CFTR. Instead, our analysis can be taken as a polemic indicating that SNPs in CFTR and many other mammalian genes may have been generated--and continue to accrue--in a fundamentally more organized manner than would otherwise have been expected. The resulting viewpoint contradicts earlier claims of 'directional' or 'intelligent design-type' SNP formation, and has important implications regarding the pace of DNA adaptation, the genesis of conserved non-coding DNA, and the extent to which eukaryotic SNP formation should be viewed as adaptive.

Expression-based network biology identifies immune-related functional modules involved in plant defense.

BACKGROUND: Plants respond to diverse environmental cues including microbial perturbations by coordinated regulation of thousands of genes. These intricate transcriptional regulatory interactions depend on the recognition of specific promoter sequences by regulatory transcription factors. The combinatorial and cooperative action of multiple transcription factors defines a regulatory network that enables plant cells to respond to distinct biological signals. The identification of immune-related modules in large-scale transcriptional regulatory networks can reveal the mechanisms by which exposure to a pathogen elicits a precise phenotypic immune response. RESULTS: We have generated a large-scale immune co-expression network using a comprehensive set of Arabidopsis thaliana (hereafter Arabidopsis) transcriptomic data, which consists of a wide spectrum of immune responses to pathogens or pathogen-mimicking stimuli treatments. We employed both linear and non-linear models to generate Arabidopsis immune co-expression regulatory (AICR) network. We computed network topological properties and ascertained that this newly constructed immune network is densely connected, possesses hubs, exhibits high modularity, and displays hallmarks of a "real" biological network. We partitioned the network and identified 156 novel modules related to immune functions. Gene Ontology (GO) enrichment analyses provided insight into the key biological processes involved in determining finely tuned immune responses. We also developed novel software called OCCEAN (One Click Cis-regulatory Elements ANalysis) to discover statistically enriched promoter elements in the upstream regulatory regions of Arabidopsis at a whole genome level. We demonstrated that OCCEAN exhibits higher precision than the existing promoter element discovery tools. In light of known and newly discovered cis-regulatory elements, we evaluated biological significance of two key immune-related functional modules and proposed mechanism(s) to explain how large sets of diverse GO genes coherently function to mount effective immune responses. CONCLUSIONS: We used a network-based, top-down approach to discover immune-related modules from transcriptomic data in Arabidopsis. Detailed analyses of these functional modules reveal new insight into the topological properties of immune co-expression networks and a comprehensive understanding of multifaceted plant defense responses. We present evidence that our newly developed software, OCCEAN, could become a popular tool for the Arabidopsis research community as well as potentially expand to analyze other eukaryotic genomes.

Massive microRNA sequence conservation and prevalence in human and chimpanzee introns.

Human and chimpanzee introns contain numerous sequences strongly related to known microRNA hairpin structures. The relative frequency is precisely maintained across all chromosomes, suggesting the possible co-evolution of gene networks dependent upon microRNA regulation and with origins corresponding to the advent of primate transposable elements (TEs). While the motifs are known to be derived from transposable elements, the most common are far more numerous than expected from the number of TEs and their paralogous sequences, and exhibit striking conservation in comparison to the surrounding TE sequence context. Several of these motifs also exhibit structural complimentarity to each other, suggesting a pairing function at the level of DNA or RNA. These "pseudomicroRNAs," in semblance to pseudogenes, include hundreds of thousands of vestigial paralogs of primate microRNAs, many of which may have functioned historically or remain active today.

The non-random distribution of intronless human genes across molecular function categories.

It has been proposed that one function of introns is to coordinate expression across/within networks of related genes. This same hypothesis also predicts that intronless genes will not be uniformly distributed among the functional categories of human genes, but will be found most frequently in those categories that have less need to communicate changes in their expression. Statistical analysis demonstrates significant clustering of single exon genes among those with binding or signal transduction/receptor activities and fewer than expected among those with catalytic function.

Micro-RNA-like effects of complete intronic sequences.

MicroRNAs (miRNAs) have been suggested as suppressors of numerous target genes in human cells. In this report, we present gene chip array data indicating that in the absence of miRNA sequences, complete human introns are similarly capable of coordinating expression of large numbers of gene products at spatially diverse sites within the genome. The expression of selected intronic sequences (6a, 14b and 23) derived from the cystic fibrosis transmembrane conductance regulator (CFTR) gene caused extensive and specific transcriptional changes in epithelial cells (HeLa) that do not normally express this gene product. Each intron initiated a distinctive pattern of gene transcription. Affected genes such as FOXF1, sucrase-isomaltase, collagen, interferon, complement and thrombospondin 1 have previously been linked to CFTR function or are known to contribute to the related processes of epithelial differentiation and repair. A possible regulatory function of this nature has not been demonstrated previously for non-coding sequences within eukaryotic DNA. The results are consistent with the observation that splicesomal introns are found only in eukaryotic organisms and that the number of introns increases with phylogenetic complexity.

Comparative gene expression profile analysis of GLI and c-MYC in an epithelial model of malignant transformation.

Transcription factor oncogenes such as GLI and c-MYC are central to the pathogenesis of human tumors. GLI encodes a zinc finger protein that is activated by Sonic Hedgehog signaling. Mutations in this pathway induce GLI expression in basal cell carcinoma, and expression of GLI in mice is sufficient to induce these skin tumors. We used microarrays to identify transcripts regulated by GLI or c-MYC after retroviral transduction and short-term culture of epithelial RK3E cells. Although each of these oncogenes induces malignant transformation of RK3E, two distinct sets of genes were identified. Of approximately 17,500 transcripts represented on the microarrays, GLI up-regulated the expression of 158 and repressed the expression of 52. In contrast, transcripts regulated by c-MYC were mainly repressed (424 of 682 regulated transcripts). Transcripts induced by the GLI transgene are likewise expressed in association with endogenous GLI in Ptch-deficient murine fibroblasts or in human skin tumors, but are not up-regulated in RK3E cells transformed by c-MYC, KLF4, or HRAS1. Unlike these other oncogenes, GLI induced the expression of mesenchymal cell markers including Snail, a zinc finger protein implicated in epithelial-mesenchymal transition in development and during tumor progression. A novel GLI-estrogen receptor fusion protein rapidly induced Snail mRNA expression in a manner like Ptch, a known direct transcriptional target gene. Induction of Snail expression and epithelial-mesenchymal transition by GLI may account for certain histopathological features of basal cell carcinoma, such as the absence of a well-defined, intraepithelial precursor lesion. In addition, consistent expression of the newly identified GLI-induced transcripts within GLI-expressing tumors in vivo indicates that oncogene-specific transcriptional profiles may be useful diagnostic tools for analysis of human tumors.