Mutational signatures of de-differentiation in functional non-coding regions of melanoma genomes

Much emphasis has been placed on the identification, functional characterization, and therapeutic potential of somatic variants in tumor genomes. However, the majority of somatic variants lie outside coding regions and their role in cancer progression remains to be determined. In order to establish a system to test the functional importance of non-coding somatic variants in cancer, we created a low-passage cell culture of a metastatic melanoma tumor sample. As a foundation for interpreting functional assays, we performed whole-genome sequencing and analysis of this cell culture, the metastatic tumor from which it was derived, and the patient-matched normal genomes. When comparing somatic mutations identified in the cell culture and tissue genomes, we observe concordance at the majority of single nucleotide variants, whereas copy number changes are more variable. To understand the functional impact of non-coding somatic variation, we leveraged functional data generated by the ENCODE Project Consortium. We analyzed regulatory regions derived from multiple different cell types and found that melanocyte-specific regions are among the most depleted for somatic mutation accumulation. Significant depletion in other cell types suggests the metastatic melanoma cells de-differentiated to a more basal regulatory state. Experimental identification of genome-wide regulatory sites in two different melanoma samples supports this observation. Together, these results show that mutation accumulation in metastatic melanoma is nonrandom across the genome and that a de-differentiated regulatory architecture is common among different samples. Our findings enable identification of the underlying genetic components of melanoma and define the differences between a tissue-derived tumor sample and the cell culture created from it. Such information helps establish a broader mechanistic understanding of the linkage between non-coding genomic variations and the cellular evolution of cancer.     

Correlations between coding and contiguous non-coding sequences in isochore families from vertebrate genomes

Many years ago compositional correlations were found to hold between coding and contiguous non-coding sequences. These correlations were essentially studied in whole genomes of mammals, which are characterized by strong compositional heterogeneities. Here we investigated whether these correlations also hold within the much more homogeneous isochore families. This point was checked not only in the case of mammals, but also in that of phylogenetically distant vertebrates, which are characterized by very different compositional patterns. Indeed, these are remarkably different in cold- and warm-blooded vertebrates. Fish genomes, for instance, are much more homogeneous than those of mammals and birds. The compositional correlations between coding sequences and the corresponding introns, or their 5′ and 3′ flanking regions, were studied in the isochore families of the fully sequenced genomes from four fishes (Brachydanio rerio, Oryzias latipes, Gasterosteus aculeatus and Tetraodon nigroviridis), human and chicken.     

Birth-and-death evolution of protein-coding regions and concerted evolution of non-coding regions in the multi-component genomes of nanoviruses

Genomes of the four plant viruses of the genus Nanovirus consist of multiple circular single-stranded DNA components, each of which encodes a single protein. Protein phylogenies supported the hypothesis that faba bean necrotic yellows virus (FBNYV) and milk vetch disease virus (MDV) are sister taxa; that subterranean clover stunt virus (SCSV) branched next; and that banana bunchy top virus (BBTV) is an outgroup to the three other species. The phylogeny of replication (Rep) proteins indicate that this small viral multi-gene family has evolved by a process of duplication and subsequent loss of Rep-encoding genome components, analogous to the "birth-and-death" process of evolution which has been described in eukaryotic multi-gene families. By contrast, repeated recombinational events between components were found to have homogenized the non-coding portions of several components encoding unrelated components. For example, as result of recent recombination a portion of the non-coding region is virtually identical among SCSV components 1, 3, 4, 5, and 7. Thus, there is a process of concerted evolution of non-coding regions of Nanovirus genome components, which raises the possibility that certain non-coding regions are subject to functional constraint.     

Identification of genes with fast-evolving regions in microbial genomes

Complete sequences of multiple strains of the same microbial species provide an invaluable source for studying the evolutionary dynamics between orthologous genes over a relatively short time scale. Usually the intensity of the selection pressure is inferred from a comparison between the nonsynonymous substitution rate and the synonymous substitution rate. In this paper, we propose an alternative method for detecting genes with one or more fast-evolving regions from pairwise comparisons of orthologous genes. Our method looks for regions with overrepresented nonsynonymous mutations along the alignment, and requires a higher nonsynonymous evolution rate in those regions than the neutral evolution rate. It identifies gene targets under intensive selection pressure that are not detected from the conventional rate comparison analysis. For those identified genes with known annotations, most of them have a clear role in processes such as bacterial defense and host–pathogen interactions. Gene sets reported from our method provide a measure of the phenotypic divergence between two closely related genomes.     

Mutation and selection on the anticodon of tRNA genes in vertebrate mitochondrial genomes

The H-strand of vertebrate mitochondrial DNA is left single-stranded for hours during the slow DNA replication. This facilitates C-->U mutations on the H-strand (and consequently G-->A mutations on the L-strand) via spontaneous deamination which occurs much more frequently on single-stranded than on double-stranded DNA. For the 12 coding sequences (CDS) collinear with the L-strand, NNY synonymous codon families (where N stands for any of the four nucleotides and Y stands for either C or U) end mostly with C, and NNR and NNN codon families (where R stands for either A or G) end mostly with A. For the lone ND6 gene on the other strand, the codon bias is the opposite, with NNY codon families ending mostly with U and NNR and NNN codon families ending mostly with G. These patterns are consistent with the strand-specific mutation bias. The codon usage biased towards C-ending and A-ending in the 12 CDS sequences affects the codon-anticodon adaptation. The wobble site of the anticodon is always G for NNY codon families dominated by C-ending codons and U for NNR and NNN codon families dominated by A-ending codons. The only, but consistent, exception is the anticodon of tRNA-Met which consistently has a 5'-CAU-3' anticodon base-pairing with the AUG codon (the translation initiation codon) instead of the more frequent AUA. The observed CAU anticodon (matching AUG) would increase the rate of translation initiation but would reduce the rate of peptide elongation because most methionine codons are AUA, whereas the unobserved UAU anticodon (matching AUA) would increase the elongation rate at the cost of translation initiation rate. The consistent CAU anticodon in tRNA-Met suggests the importance of maximizing the rate of translation initiation.     

On the phylogeny of Mustelidae subfamilies: analysis of seventeen nuclear non-coding loci and mitochondrial complete genomes

Background  

Mustelidae, as the largest and most-diverse family of order Carnivora, comprises eight subfamilies. Phylogenetic relationships among these Mustelidae subfamilies remain argumentative subjects in recent years. One of the main reasons is that the mustelids represent a typical example of rapid evolutionary radiation and recent speciation event. Prior investigation has been concentrated on the application of different mitochondrial (mt) sequence and nuclear protein-coding data, herein we employ 17 nuclear non-coding loci (>15 kb), in conjunction with mt complete genome data (>16 kb), to clarify these enigmatic problems.     

Identification of transposon-like elements in non-coding regions of tomato ACC oxidase genes

1-aminocyclopropane-1-carboxylate (ACC) oxidase, which catalyses the terminal step in ethylene biosynthesis, is encoded by a small multigene family in tomato that is differentially expressed in response to developmental and environmental cues. In this study we report the isolation and sequencing of approximately 2 kb of 5′-flanking sequence of three tomato ACC oxidase genes (LEACO1, LEACO2, LEACO3) and the occurrence of class I and class II mobile element-like insertions in promoter and intron regions of two of them. The LEACO1 upstream region contains a 420-bp direct repeat which is present in multiple copies in the tomato genome and is very similar to sequences in the promoters of the tomato E4 and 2A11 genes. The region covering the repeats resembles the remnant of a retrotransposon. Two copies of a small transposable element, belonging to the Stowaway inverted repeat element family, have been found in the 5′-flanking sequence and the third intron of LEACO3. Received: 8 August 1996 / Accepted: 4 November 1996     

Identification of a conserved sequence in the non-coding regions of many human genes.

We have analyzed a sequence of approximately 70 base pairs (bp) that shows a high degree of similarity to sequences present in the non-coding regions of a number of human and other mammalian genes. The sequence was discovered in a fragment of human genomic DNA adjacent to an integrated hepatitis B virus genome in cells derived from human hepatocellular carcinoma tissue. When one of the viral flanking sequences was compared to nucleotide sequences in GenBank, more than thirty human genes were identified that contained a similar sequence in their non-coding regions. The sequence element was usually found once or twice in a gene, either in an intron or in the 5' or 3' flanking regions. It did not share any similarities with known short interspersed nucleotide elements (SINEs) or presently known gene regulatory elements. This element was highly conserved at the same position within the corresponding human and mouse genes for myoglobin and N-myc, indicating evolutionary conservation and possible functional importance. Preliminary DNase I footprinting data suggested that the element or its adjacent sequences may bind nuclear factors to generate specific DNase I hypersensitive sites. The size, structure, and evolutionary conservation of this sequence indicates that it is distinct from other types of short interspersed repetitive elements. It is possible that the element may have a cis-acting functional role in the genome.     

The search of miRNA genes in Bombyx mori nuclear polyhedrosis virus genomes regions complementary to the latest genes

The search of miRNA genes in Bombyx mori nuclear polyhedrosis virus genome region complementary to very late genes has been carried out. The search miRNA algorithm in silico was developed by us. It was shown that NPV B. mori genome region containing orf4 gene complementary to ph gene encodes the potential miRNA. NPV B. mori genome region containing p74 gene complementary to p10 gene encodes mature miRNA and potential miRNA. The genome region containing orf1629 encodes two small non-coding RNAs complementary to orf 5'-end of polyhedrin miRNA. From obtained results it is proposed that two small noncoding RNAs complementary to regions of polyhedrin miRNA are included in polyhedra.     

Orthologs of key vertebrate neural genes are expressed during neurogenesis in the annelid Platynereis dumerilii

SUMMARY The molecular mechanisms underlying the formation and patterning of the nervous system are relatively poorly understood for lophotrochozoans (like annelids) as compared with ecdysozoans (especially Drosophila ) and deuterostomes (especially vertebrates). Therefore, we have undertaken a candidate gene approach to study aspects of neurogenesis in a polychaete annelid Platynereis dumerilii . We determined the spatiotemporal expression for Platynereis orthologs of four genes ( SoxB, Churchill, prospero / Prox , and SoxC) known to play key roles in vertebrate neurogenesis. During Platynereis development, SoxB is expressed in the neuroectoderm and its expression switches off when committed neural precursors are formed. Subsequently, Prox is expressed in all differentiating neural precursors in the central and peripheral nervous systems. Finally, SoxC and Churchill are transcribed in patterns consistent with their involvement in neural differentiation. The expression patterns of Platynereis SoxB and Prox closely resemble those in Drosophila and vertebrates—this suggests that orthologs of these genes play similar neurogenic roles in all bilaterians. Whereas Platynereis SoxC , like its vertebrate orthologs, plays a role in neural cell differentiation, related genes in Drosophila do not appear to be involved in neurogenesis. Finally, conversely to Churchill in Platynereis , vertebrate orthologs of this gene are expressed during neuroectoderm formation, but not later during nerve cell differentiation; in the insect lineage, homologs of these genes have been secondarily lost. In spite of such instances of functional divergence or loss, the present study shows conspicuous similarities in the genetic control of neurogenesis among bilaterians. These commonalities suggest that key features of the genetic program for neurogenesis are ancestral to bilaterians.     

The spectra of point mutations in vertebrate genomes

In spite of the importance of point mutations for evolution and human diseases, their natural spectrum of incidence in different species is not known. Here I propose to determine these spectra by comparing consecutive sequence periods in stretches of repetitive DNA. The article presents the analysis of more than 51,000 such point mutations identified by this approach in the genomes of human, chimpanzee, rat, mouse, pufferfish, zebrafish, and sea squirt. I propose to explain the observed spectra by auto‐mutagenic mechanisms of genome variation involving the inter‐conversions of nucleotides, single base‐pair inversions and their combinations.     

Chaperonin genes on the rise: new divergent classes and intense duplication in human and other vertebrate genomes

Background  

Chaperonin proteins are well known for the critical role they play in protein folding and in disease. However, the recent identification of three diverged chaperonin paralogs associated with the human Bardet-Biedl and McKusick-Kaufman Syndromes (BBS and MKKS, respectively) indicates that the eukaryotic chaperonin-gene family is larger and more differentiated than previously thought. The availability of complete genome sequences makes possible a definitive characterization of the complete set of chaperonin sequences in human and other species.     

Identification and analysis of genes and pseudogenes within duplicated regions in the human and mouse genomes

The identification and classification of genes and pseudogenes in duplicated regions still constitutes a challenge for standard automated genome annotation procedures. Using an integrated homology and orthology analysis independent of current gene annotation, we have identified 9,484 and 9,017 gene duplicates in human and mouse, respectively. On the basis of the integrity of their coding regions, we have classified them into functional and inactive duplicates, allowing us to define the first consistent and comprehensive collection of 1,811 human and 1,581 mouse unprocessed pseudogenes. Furthermore, of the total of 14,172 human and mouse duplicates predicted to be functional genes, as many as 420 are not included in current reference gene databases and therefore correspond to likely novel mammalian genes. Some of these correspond to partial duplicates with less than half of the length of the original source genes, yet they are conserved and syntenic among different mammalian lineages. The genes and unprocessed pseudogenes obtained here will enable further studies on the mechanisms involved in gene duplication as well as of the fate of duplicated genes.