Intra-tumor genetic heterogeneity in rectal cancer.

Colorectal cancer arises in part from the cumulative effects of multiple gene lesions. Recent studies in selected cancer types have revealed significant intra-tumor genetic heterogeneity and highlighted its potential role in disease progression and resistance to therapy. We hypothesized the existence of significant intra-tumor genetic heterogeneity in rectal cancers involving variations in localized somatic mutations and copy number abnormalities. Two or three spatially disparate regions from each of six rectal tumors were dissected and subjected to the next-generation whole-exome DNA sequencing, Oncoscan SNP arrays, and targeted confirmatory sequencing and analysis. The resulting data were integrated to define subclones using SciClone. Mutant-allele tumor heterogeneity (MATH) scores, mutant allele frequency correlation, and mutation percent concordance were calculated, and copy number analysis including measurement of correlation between samples was performed. Somatic mutations profiles in individual cancers were similar to prior studies, with some variants found in previously reported significantly mutated genes and many patient-specific mutations in each tumor. Significant intra-tumor heterogeneity was identified in the spatially disparate regions of individual cancers. All tumors had some heterogeneity but the degree of heterogeneity was quite variable in the samples studied. We found that 67-97% of exonic somatic mutations were shared among all regions of an individual's tumor. The SciClone computational method identified 2-8 shared and unshared subclones in the spatially disparate areas in each tumor. MATH scores ranged from 7 to 41. Allele frequency correlation scores ranged from R(2)=0.69-0.96. Measurements of correlation between samples for copy number changes varied from R(2)=0.74-0.93. All tumors had some heterogeneity, but the degree was highly variable in the samples studied. The occurrence of significant intra-tumor heterogeneity may allow selected tumors to have a genetic reservoir to draw from in their evolutionary response to therapy and other challenges.

Community transcriptomic assembly reveals microbes that contribute to deep-sea carbon and nitrogen cycling.

The deep ocean is an important component of global biogeochemical cycles because it contains one of the largest pools of reactive carbon and nitrogen on earth. However, the microbial communities that drive deep-sea geochemistry are vastly unexplored. Metatranscriptomics offers new windows into these communities, but it has been hampered by reliance on genome databases for interpretation. We reconstructed the transcriptomes of microbial populations from Guaymas Basin, in the deep Gulf of California, through shotgun sequencing and de novo assembly of total community RNA. Many of the resulting messenger RNA (mRNA) contiguous sequences contain multiple genes, reflecting co-transcription of operons, including those from dominant members. Also prevalent were transcripts with only limited representation (2.8 times coverage) in a corresponding metagenome, including a considerable portion (1.2 Mb total assembled mRNA sequence) with similarity (96%) to a marine heterotroph, Alteromonas macleodii. This Alteromonas and euryarchaeal marine group II populations displayed abundant transcripts from amino-acid transporters, suggesting recycling of organic carbon and nitrogen from amino acids. Also among the most abundant mRNAs were catalytic subunits of the nitrite oxidoreductase complex and electron transfer components involved in nitrite oxidation. These and other novel genes are related to novel Nitrospirae and have limited representation in accompanying metagenomic data. High throughput sequencing of 16S ribosomal RNA (rRNA) genes and rRNA read counts confirmed that Nitrospirae are minor yet widespread members of deep-sea communities. These results implicate a novel bacterial group in deep-sea nitrite oxidation, the second step of nitrification. This study highlights metatranscriptomic assembly as a valuable approach to study microbial communities.

RoBuST: an integrated genomics resource for the root and bulb crop families Apiaceae and Alliaceae.

BACKGROUND: Root and bulb vegetables (RBV) include carrots, celeriac (root celery), parsnips (Apiaceae), onions, garlic, and leek (Alliaceae)--food crops grown globally and consumed worldwide. Few data analysis platforms are currently available where data collection, annotation and integration initiatives are focused on RBV plant groups. Scientists working on RBV include breeders, geneticists, taxonomists, plant pathologists, and plant physiologists who use genomic data for a wide range of activities including the development of molecular genetic maps, delineation of taxonomic relationships, and investigation of molecular aspects of gene expression in biochemical pathways and disease responses. With genomic data coming from such diverse areas of plant science, availability of a community resource focused on these RBV data types would be of great interest to this scientific community. DESCRIPTION: The RoBuST database has been developed to initiate a platform for collecting and organizing genomic information useful for RBV researchers. The current release of RoBuST contains genomics data for 294 Alliaceae and 816 Apiaceae plant species and has the following features: (1) comprehensive sequence annotations of 3663 genes 5959 RNAs, 22,723 ESTs and 11,438 regulatory sequence elements from Apiaceae and Alliaceae plant families; (2) graphical tools for visualization and analysis of sequence data; (3) access to traits, biosynthetic pathways, genetic linkage maps and molecular taxonomy data associated with Alliaceae and Apiaceae plants; and (4) comprehensive plant splice signal repository of 659,369 splice signals collected from 6015 plant species for comparative analysis of plant splicing patterns. CONCLUSIONS: RoBuST, available at http://robust.genome.com, provides an integrated platform for researchers to effortlessly explore and analyze genomic data associated with root and bulb vegetables.

Targeted genome-wide enrichment of functional regions.

Only a small fraction of large genomes such as that of the human contains the functional regions such as the exons, promoters, and polyA sites. A platform technique for selective enrichment of functional genomic regions will enable several next-generation sequencing applications that include the discovery of causal mutations for disease and drug response. Here, we describe a powerful platform technique, termed "functional genomic fingerprinting" (FGF), for the multiplexed genomewide isolation and analysis of targeted regions such as the exome, promoterome, or exon splice enhancers. The technique employs a fixed part of a uniquely designed Fixed-Randomized primer, while the randomized part contains all the possible sequence permutations. The Fixed-Randomized primers bind with full sequence complementarity at multiple sites where the fixed sequence (such as the splice signals) occurs within the genome, and multiplex amplify many regions bounded by the fixed sequences (e.g., exons). Notably, validation of this technique using cardiac myosin binding protein-C (MYBPC3) gene as an example strongly supports the application and efficacy of this method. Further, assisted by genomewide computational analyses of such sequences, the FGF technique may provide a unique platform for high-throughput sample production and analysis of targeted genomic regions by the next-generation sequencing techniques, with powerful applications in discovering disease and drug response genes.

AspAlt: A tool for inter-database, inter-genomic and user-specific comparative analysis of alternative transcription and alternative splicing in 46 eukaryotes.

We have developed AspAlt-a web-based comparative analytical platform for exploring the variations in alternative transcription (AT) events and alternative splicing (AS) events in eukaryotes. AspAlt provides integrated access to 2.1 million AT-AS annotations from 1,58,876 multi-isoform genes and has the following user-friendly analytical features: (1) advanced graphical display to visualize and analyze AT-AS events in 46 eukaryotic genomes; (2) compare and identify the differences in AT-AS patterns among a group of genes specified by the user or among homologous gene groups; (3) inter-database comparative viewer to analyze the differences in the AT-AS annotations for the same gene among Ensembl, RefSeq and AceView databases; (4) dynamically classify and generate graphical plots of AT-AS events from mRNA annotations submitted by the user; and (5) download genomic AT-AS annotations of 46 eukaryotes in XML and tab-delimited formats. The AspAlt resource is available at http://66.170.16.154/AspAlt.

ExDom: an integrated database for comparative analysis of the exon-intron structures of protein domains in eukaryotes.

We have developed ExDom, a unique database for the comparative analysis of the exon-intron structures of 96 680 protein domains from seven eukaryotic organisms (Homo sapiens, Mus musculus, Bos taurus, Rattus norvegicus, Danio rerio, Gallus gallus and Arabidopsis thaliana). ExDom provides integrated access to exon-domain data through a sophisticated web interface which has the following analytical capabilities: (i) intergenomic and intragenomic comparative analysis of exon-intron structure of domains; (ii) color-coded graphical display of the domain architecture of proteins correlated with their corresponding exon-intron structures; (iii) graphical analysis of multiple sequence alignments of amino acid and coding nucleotide sequences of homologous protein domains from seven organisms; (iv) comparative graphical display of exon distributions within the tertiary structures of protein domains; and (v) visualization of exon-intron structures of alternative transcripts of a gene correlated to variations in the domain architecture of corresponding protein isoforms. These novel analytical features are highly suited for detailed investigations on the exon-intron structure of domains and make ExDom a powerful tool for exploring several key questions concerning the function, origin and evolution of genes and proteins. ExDom database is freely accessible at: http://66.170.16.154/ExDom/.

Origination of the split structure of spliceosomal genes from random genetic sequences.

The mechanism by which protein-coding portions of eukaryotic genes came to be separated by long non-coding stretches of DNA, and the purpose for this perplexing arrangement, have remained unresolved fundamental biological problems for three decades. We report here a plausible solution to this problem based on analysis of open reading frame (ORF) length constraints in the genomes of nine diverse species. If primordial nucleic acid sequences were random in sequence, functional proteins that are innately long would not be encoded due to the frequent occurrence of stop codons. The best possible way that a long protein-coding sequence could have been derived was by evolving a split-structure from the random DNA (or RNA) sequence. Results of the systematic analyses of nine complete genome sequences presented here suggests that perhaps the major underlying structural features of split-genes have evolved due to the indigenous occurrence of split protein-coding genes in primordial random nucleotide sequence. The results also suggest that intron-rich genes containing short exons may have been the original form of genes intrinsically occurring in random DNA, and that intron-poor genes containing long exons were perhaps derived from the original intron-rich genes.

EuSplice: a unified resource for the analysis of splice signals and alternative splicing in eukaryotic genes.

MOTIVATION: Despite increased availability of genome annotation data, a comprehensive resource for in-depth analysis of splice signal distributions and alternative splicing (AS) patterns in eukaryote genomes is still lacking. To meet this need, we have developed EuSplice--a unique splice-centric database which provides reliable splice signal and AS information for 23 eukaryotes. RESULTS: The EuSplice database contains 95,822 AS events and 2.1 million splice signals associated with over 270,000 protein-coding genes. The intuitive, user-friendly EuSplice web interface has powerful data mining and graphics capabilities for inter-genomic comparative analysis of splice signals, putative cryptic splice sites and AS events. Moreover, the seamless integration of splicing data to extensive gene-specific annotations, such as homolog annotations, functional information, mutations and sequence details makes EuSplice a powerful one-stop information resource for investigating the molecular mechanisms of complex splicing events, disease associations and the evolution of splicing in eukaryotes. AVAILABILITY: http://66.170.16.154/EuSplice. SUPPLEMENTARY INFORMATION: Supplementary tables and figures at Bioinfo online.