Exact mapping of prokaryotic gene starts

It is known that while the programs used to find genes in prokaryotic genomes reliably map protein-coding regions, they often fail in the exact determination of gene starts. This problem is further aggravated by sequencing errors, most notably insertions and deletions leading to frame-shifts. Therefore, the exact mapping of gene starts and identification of frame-shifts are important problems of the computer-assisted functional analysis of newly sequenced genomes. Here we review methods of gene recognition and describe a new algorithm for correction of gene starts and identification of frame-shifts in prokaryotic genomes. The algorithm is based on the comparison of nucleotide and protein sequences of homologous genes from related organisms, using the assumption that the rate of evolutionary changes in protein-coding regions is lower than that in non-coding regions. A dynamic programming algorithm is used to align protein sequences obtained by formal translation of genomic nucleotide sequences. The possibility of frame-shifts is taken into account. The algorithm was tested on several groups of related organisms: gamma-proteobacteria, the Bacillus/Clostridium group, and three Pyrococcus genomes. The testing demonstrated that, dependent or a genome, 1-10 per cent of genes have incorrect starts or contain frame-shifts. The algorithm is implemented in the program package Orthologator-GeneCorrector.     

Identifying common prognostic factors in genomic cancer studies: A novel index for censored outcomes

Background  

With the growing number of public repositories for high-throughput genomic data, it is of great interest to combine the results produced by independent research groups. Such a combination allows the identification of common genomic factors across multiple cancer types and provides new insights into the disease process. In the framework of the proportional hazards model, classical procedures, which consist of ranking genes according to the estimated hazard ratio or the p-value obtained from a test statistic of no association between survival and gene expression level, are not suitable for gene selection across multiple genomic datasets with different sample sizes. We propose a novel index for identifying genes with a common effect across heterogeneous genomic studies designed to remain stable whatever the sample size and which has a straightforward interpretation in terms of the percentage of separability between patients according to their survival times and gene expression measurements.     

Two Distinct Categories of Focal Deletions in Cancer Genomes

One of the key questions about genomic alterations in cancer is whether they are functional in the sense of contributing to the selective advantage of tumor cells. The frequency with which an alteration occurs might reflect its ability to increase cancer cell growth, or alternatively, enhanced instability of a locus may increase the frequency with which it is found to be aberrant in tumors, regardless of oncogenic impact. Here we’ve addressed this on a genome-wide scale for cancer-associated focal deletions, which are known to pinpoint both tumor suppressor genes (tumor suppressors) and unstable loci. Based on DNA copy number analysis of over one-thousand human cancers representing ten different tumor types, we observed five loci with focal deletion frequencies above 5%, including the A2BP1 gene at 16p13.3 and the MACROD2 gene at 20p12.1. However, neither RNA expression nor functional studies support a tumor suppressor role for either gene. Further analyses suggest instead that these are sites of increased genomic instability and that they resemble common fragile sites (CFS). Genome-wide analysis revealed properties of CFS-like recurrent deletions that distinguish them from deletions affecting tumor suppressor genes, including their isolation at specific loci away from other genomic deletion sites, a considerably smaller deletion size, and dispersal throughout the affected locus rather than assembly at a common site of overlap. Additionally, CFS-like deletions have less impact on gene expression and are enriched in cell lines compared to primary tumors. We show that loci affected by CFS-like deletions are often distinct from known common fragile sites. Indeed, we find that each tumor tissue type has its own spectrum of CFS-like deletions, and that colon cancers have many more CFS-like deletions than other tumor types. We present simple rules that can pinpoint focal deletions that are not CFS-like and more likely to affect functional tumor suppressors.     

Comparative genomics and gene expression analysis identifies BBS9, a new Bardet-Biedl syndrome gene

Bardet-Biedl syndrome (BBS) is an autosomal recessive, genetically heterogeneous, pleiotropic human disorder characterized by obesity, retinopathy, polydactyly, renal and cardiac malformations, learning disabilities, and hypogenitalism. Eight BBS genes representing all known mapped loci have been identified. Mutation analysis of the known BBS genes in BBS patients indicate that additional BBS genes exist and/or that unidentified mutations exist in the known genes. To identify new BBS genes, we performed homozygosity mapping of small, consanguineous BBS pedigrees, using moderately dense SNP arrays. A bioinformatics approach combining comparative genomic analysis and gene expression studies of a BBS-knockout mouse model was used to prioritize BBS candidate genes within the newly identified loci for mutation screening. By use of this strategy, parathyroid hormone-responsive gene B1 (B1) was found to be a novel BBS gene (BBS9), supported by the identification of homozygous mutations in BBS patients. The identification of BBS9 illustrates the power of using a combination of comparative genomic analysis, gene expression studies, and homozygosity mapping with SNP arrays in small, consanguineous families for the identification of rare autosomal recessive disorders. We also demonstrate that small, consanguineous families are useful in identifying intragenic deletions. This type of mutation is likely to be underreported because of the difficulty of deletion detection in the heterozygous state by the mutation screening methods that are used in many studies.     

Coupled analysis of gene expression and chromosomal location

Microarray technology can be used to assess simultaneously global changes in expression of mRNA or genomic DNA copy number among thousands of genes in different biological states. In many cases, it is desirable to determine if altered patterns of gene expression correlate with chromosomal abnormalities or assess expression of genes that are contiguous in the genome. We describe a method, differential gene locus mapping (DIGMAP), which aligns the known chromosomal location of a gene to its expression value deduced by microarray analysis. The method partitions microarray data into subsets by chromosomal location for each gene interrogated by an array. Microarray data in an individual subset can then be clustered by physical location of genes at a subchromosomal level based upon ordered alignment in genome sequence. A graphical display is generated by representing each genomic locus with a colored cell that quantitatively reflects its differential expression value. The clustered patterns can be viewed and compared based on their expression signatures as defined by differential values between control and experimental samples. In this study, DIGMAP was tested using previously published studies of breast cancer analyzed by comparative genomic hybridization (CGH) and prostate cancer gene expression profiles assessed by cDNA microarray experiments. Analysis of the breast cancer CGH data demonstrated the ability of DIGMAP to deduce gene amplifications and deletions. Application of the DIGMAP method to the prostate data revealed several carcinoma-related loci, including one at 16q13 with marked differential expression encompassing 19 known genes including 9 encoding metallothionein proteins. We conclude that DIGMAP is a powerful computational tool enabling the coupled analysis of microarray data with genome location.     

Identification of tumor suppressors and oncogenes from genomic and epigenetic features in ovarian cancer

The identification of genetic and epigenetic alterations from primary tumor cells has become a common method to identify genes critical to the development and progression of cancer. We seek to identify those genetic and epigenetic aberrations that have the most impact on gene function within the tumor. First, we perform a bioinformatic analysis of copy number variation (CNV) and DNA methylation covering the genetic landscape of ovarian cancer tumor cells. We separately examined CNV and DNA methylation for 42 primary serous ovarian cancer samples using MOMA-ROMA assays and 379 tumor samples analyzed by The Cancer Genome Atlas. We have identified 346 genes with significant deletions or amplifications among the tumor samples. Utilizing associated gene expression data we predict 156 genes with altered copy number and correlated changes in expression. Among these genes CCNE1, POP4, UQCRB, PHF20L1 and C19orf2 were identified within both data sets. We were specifically interested in copy number variation as our base genomic property in the prediction of tumor suppressors and oncogenes in the altered ovarian tumor. We therefore identify changes in DNA methylation and expression for all amplified and deleted genes. We statistically define tumor suppressor and oncogenic features for these modalities and perform a correlation analysis with expression. We predicted 611 potential oncogenes and tumor suppressors candidates by integrating these data types. Genes with a strong correlation for methylation dependent expression changes exhibited at varying copy number aberrations include CDCA8, ATAD2, CDKN2A, RAB25, AURKA, BOP1 and EIF2C3. We provide copy number variation and DNA methylation analysis for over 11,500 individual genes covering the genetic landscape of ovarian cancer tumors. We show the extent of genomic and epigenetic alterations for known tumor suppressors and oncogenes and also use these defined features to identify potential ovarian cancer gene candidates.     

Integrated genomic and gene expression profiling identifies two major genomic circuits in urothelial carcinoma

Similar to other malignancies, urothelial carcinoma (UC) is characterized by specific recurrent chromosomal aberrations and gene mutations. However, the interconnection between specific genomic alterations, and how patterns of chromosomal alterations adhere to different molecular subgroups of UC, is less clear. We applied tiling resolution array CGH to 146 cases of UC and identified a number of regions harboring recurrent focal genomic amplifications and deletions. Several potential oncogenes were included in the amplified regions, including known oncogenes like E2F3, CCND1, and CCNE1, as well as new candidate genes, such as SETDB1 (1q21), and BCL2L1 (20q11). We next combined genome profiling with global gene expression, gene mutation, and protein expression data and identified two major genomic circuits operating in urothelial carcinoma. The first circuit was characterized by FGFR3 alterations, overexpression of CCND1, and 9q and CDKN2A deletions. The second circuit was defined by E3F3 amplifications and RB1 deletions, as well as gains of 5p, deletions at PTEN and 2q36, 16q, 20q, and elevated CDKN2A levels. TP53/MDM2 alterations were common for advanced tumors within the two circuits. Our data also suggest a possible RAS/RAF circuit. The tumors with worst prognosis showed a gene expression profile that indicated a keratinized phenotype. Taken together, our integrative approach revealed at least two separate networks of genomic alterations linked to the molecular diversity seen in UC, and that these circuits may reflect distinct pathways of tumor development.     

Novel reference gene RPN1 for normalization of quantitative data in lung and kidney cancer

Quantitative methods of gene expression analysis in tumors require accurate data normalization, which allows comparison of different mRNA/cDNA samples with unknown concentration. For this purpose reference genes with stable expression level (such as GAPDH, ACTB, HPRT1, TBP) are used. The choice of appropriate reference genes is still actual because well-known reference genes are not suitable for certain cancer types frequently and their unreasonable use without additional tests lead to wrong conclusions. We have developed the bioinformatic approach and selected a new potential reference gene RPN1 for lung and kidney tumors. This gene is located at the long arm of chromosome 3. Our method includes mining of the dbEST and Oncomine databases and functional analysis of genes. The RPN1 was selected from 1500 candidate housekeeping genes. Using comparative genomic hybridization with NotI-microarrays we found no methylation, deletions and/or amplifications at the RPN1-containing locus in 56 non-small cell lung and 42 clear cell renal cancer samples. Using RT-qPCR we showed low variability of RPN1 mRNA level comparable to those of reference genes GAPDH and GUSB in lung and kidney cancer. The mRNA levels of two target genes coding hyalouronidases--HYAL1 and HYAL2--were estimated and normalized relative to pair RPN1--GAPDH genes for lung cancer and RPN1--GUSB for kidney cancer. These combinations were shown to be optimal for obtaining accurate and reproducible data. All obtained results allow us to suggest RPN1 as novel reference gene for quantitative data normalization in gene expression studies for lung and kidney cancers.     

The FUR1 gene of Saccharomyces cerevisiae: cloning, structure and expression of wild-type and mutant alleles

The FUR1 gene of Saccharomyces cerevisiae encodes uracil phosphoribosyltransferase (UPRTase) which catalyses the conversion of uracil into uridine 5'-monophosphate (UMP) in the pyrimidine salvage pathway. The FUR1 gene is included in a 2.1 kb genomic segment of DNA and is transcribed into a 1 kb poly(A)+mRNA. Sequencing has determined a 753 bp open reading frame capable of encoding a protein of 251 amino acids. The FUR1 genes for three recessive fur1 alleles, having different sensibilities to 5-fluorouridine (5-FUR) but identical levels of resistance to 5-fluorouracil (5-FU), were cloned and sequenced. Single bp changes located in different regions of the gene were found in each mutant. Two in vitro-constructed deletions of the FUR1 gene have been integrated at the chromosomal locus, giving strains with 5-FURR and 5-FURR mutant phenotype. Assays of UPRTase, uridine kinase, uridine ribohydrolase and uridine 5'-monophosphate nucleotidase enzymatic activities, in extracts of strains where the FUR1 gene is overexpressed or deleted, indicate that the FUR1 encoded protein possesses only UPRTase activity.     

Many ribosomal protein genes are cancer genes in zebrafish

We have generated several hundred lines of zebrafish (Danio rerio), each heterozygous for a recessive embryonic lethal mutation. Since many tumor suppressor genes are recessive lethals, we screened our colony for lines that display early mortality and/or gross evidence of tumors. We identified 12 lines with elevated cancer incidence. Fish from these lines develop malignant peripheral nerve sheath tumors, and in some cases also other tumor types, with moderate to very high frequencies. Surprisingly, 11 of the 12 lines were each heterozygous for a mutation in a different ribosomal protein (RP) gene, while one line was heterozygous for a mutation in a zebrafish paralog of the human and mouse tumor suppressor gene, neurofibromatosis type 2. Our findings suggest that many RP genes may act as haploinsufficient tumor suppressors in fish. Many RP genes might also be cancer genes in humans, where their role in tumorigenesis could easily have escaped detection up to now.     

Copy number variations due to large genomic deletion in X-linked chronic granulomatous disease

Mutations in genes for any of the six subunits of NADPH oxidase cause chronic granulomatous disease (CGD), but almost 2/3 of CGD cases are caused by mutations in the X-linked CYBB gene, also known as NAD (P) H oxidase 2. Approximately 260 patients with CGD have been reported in Japan, of whom 92 were shown to have mutations of the CYBB gene and 16 to have chromosomal deletions. However, there has been very little detailed analysis of the range of the deletion or close understanding of the disease based on this. We therefore analyzed genomic rearrangements in X-linked CGD using array comparative genomic hybridization analysis, revealing the extent and the types of the deletion genes. The subjects were five Japanese X-linked CGD patients estimated to have large base deletions of 1 kb or more in the CYBB gene (four male patients, one female patient) and the mothers of four of those patients. The five Japanese patients were found to range from a patient exhibiting deletions only of the CYBB gene to a female patient exhibiting an extensive DNA deletion and the DMD and CGD phenotype manifested. Of the other three patients, two exhibited CYBB, XK, and DYNLT3 gene deletions. The remaining patient exhibited both a deletion encompassing DNA subsequent to the CYBB region following intron 2 and the DYNLT3 gene and a complex copy number variation involving the insertion of an inverted duplication of a region from the centromere side of DYNLT3 into the deleted region.     

A hypervariable 130-kilobase genomic region of Magnetospirillum gryphiswaldense comprises a magnetosome island which undergoes frequent rearrangements during stationary growth

Genes involved in magnetite biomineralization are clustered in the genome of the magnetotactic bacterium Magnetospirillum gryphiswaldense. We analyzed a 482-kb genomic fragment, in which we identified an approximately 130-kb region representing a putative genomic "magnetosome island" (MAI). In addition to all known magnetosome genes, the MAI contains genes putatively involved in magnetosome biomineralization and numerous genes with unknown functions, as well as pseudogenes, and it is particularly rich in insertion elements. Substantial sequence polymorphism of clones from different subcultures indicated that this region undergoes frequent rearrangements during serial subcultivation in the laboratory. Spontaneous mutants affected in magnetosome formation arise at a frequency of up to 10(-2) after prolonged storage of cells at 4 degrees C or exposure to oxidative stress. All nonmagnetic mutants exhibited extended and multiple deletions in the MAI and had lost either parts of or the entire mms and mam gene clusters encoding magnetosome proteins. The mutations were polymorphic with respect to the sites and extents of deletions, but all mutations were found to be associated with the loss of various copies of insertion elements, as revealed by Southern hybridization and PCR analysis. Insertions and deletions in the MAI were also found in different magnetosome-producing clones, indicating that parts of this region are not essential for the magnetic phenotype. Our data suggest that the genomic MAI undergoes frequent transposition events, which lead to subsequent deletion by homologous recombination under physiological stress conditions. This can be interpreted in terms of adaptation to physiological stress and might contribute to the genetic plasticity and mobilization of the magnetosome island.     

Chromosomal deletions and inversions mediated by TALENs and CRISPR/Cas in zebrafish

Customized TALENs and Cas9/gRNAs have been used for targeted mutagenesis in zebrafish to induce indels into protein-coding genes. However, indels are usually not sufficient to disrupt the function of non-coding genes, gene clusters or regulatory sequences, whereas large genomic deletions or inversions are more desirable for this purpose. By injecting two pairs of TALEN mRNAs or two gRNAs together with Cas9 mRNA targeting distal DNA sites of the same chromosome, we obtained predictable genomic deletions or inversions with sizes ranging from several hundred bases to nearly 1 Mb. We have successfully achieved this type of modifications for 11 chromosomal loci by TALENs and 2 by Cas9/gRNAs with different combinations of gRNA pairs, including clusters of miRNA and protein-coding genes. Seven of eight TALEN-targeted lines transmitted the deletions and one transmitted the inversion through germ line. Our findings indicate that both TALENs and Cas9/gRNAs can be used as an efficient tool to engineer genomes to achieve large deletions or inversions, including fragments covering multiple genes and non-coding sequences. To facilitate the analyses and application of existing ZFN, TALEN and CRISPR/Cas data, we have updated our EENdb database to provide a chromosomal view of all reported engineered endonucleases targeting human and zebrafish genomes.     

Localization of the Fanconi anemia complementation group D gene to a 200-kb region on chromosome 3p25.3

Fanconi anemia (FA) is a rare autosomal recessive disease manifested by bone-marrow failure and an elevated incidence of cancer. Cells taken from patients exhibit spontaneous chromosomal breaks and rearrangements. These breaks and rearrangements are greatly elevated by treatment of FA cells with the use of DNA cross-linking agents. The FA complementation group D gene (FANCD) has previously been localized to chromosome 3p22-26, by use of microcell-mediated chromosome transfer. Here we describe the use of noncomplemented microcell hybrids to identify small overlapping deletions that narrow the FANCD critical region. A 1.2-Mb bacterial-artificial-chromosome (BAC)/P1 contig was constructed, bounded by the marker D3S3691 distally and by the gene ATP2B2 proximally. The contig contains at least 36 genes, including the oxytocin receptor (OXTR), hOGG1, the von Hippel-Lindau tumor-suppressor gene (VHL), and IRAK-2. Both hOGG1 and IRAK-2 were excluded as candidates for FANCD. BACs were then used as probes for FISH analyses, to map the extent of the deletions in four of the noncomplemented microcell hybrid cell lines. A narrow region of common overlapping deletions limits the FANCD critical region to approximately 200 kb. The three candidate genes in this region are TIGR-A004X28, SGC34603, and AA609512.     

Estimating genomic instability mediated by Alu retroelements in breast cancer

Alu-PCR is a relatively simple technique that can be used to investigate genomic instability in cancer. This technique allows identification of the loss, gain or amplification of gene sequences based on the analysis of segments between two Alu elements coupled with quantitative and qualitative analyses of the profiles obtained from tumor samples, surgical margins and blood. In this work, we used Alu-PCR to identify gene alterations in ten patients with invasive ductal breast cancer. Several deletions and insertions were identified, indicating genomic instability in the tumor and adjacent normal tissue. Although not associated with specific genes, the alterations, which involved chromosomal bands 1p36.23, 1q41, 11q14.3, 13q14.2, occurred in areas of well-known genomic instability in breast and other types of cancer. These results indicate the potential usefulness of Alu-PCR in identifying altered gene sequences in breast cancer. However, caution is required in its application since the Alu primer can produce non-specific amplification.     

DNA Methylation Differences Associated with Tumor Tissues Identified by Genome Scanning Analysis

Most investigations on the role of DNA methylation in cancer have focused on epigenetic changes associated with known tumor suppressor genes. This may have led to an underestimation of the number of CpG islands altered by DNA methylation, since it is possible that a subset of unknown genes relevant to cancer development may preferentially be affected by epigenetic rather than genetic means and would not be identified as familial deletions, mutations, or loss of heterozygosity. We used a recently developed screening procedure (methylation-sensitive arbitrarily primed-polymerase chain reaction to scan genomic DNA for CpG islands methylated in white blood cells (WBCs) and in tumor tissues. DNA methylation pattern analysis showed little interindividual differences in the WBCs and normal epithelium (adjacent to colon, bladder, and prostate cancer cells), but with some tissue-specific differences. Cancer cells showed marked methylation changes that varied considerably between different tumors, suggesting variable penetrance of the methylation phenotype in patients. Direct sequencing of 8 of 45 bands altered in these cancers showed that several of them were CpG islands, and 2 of these sequences were identified in GenBank. Surprisingly, three of the bands studied corresponded to transcribed regions of genes. Thus, hypermethylation of CpG islands in cancer cells is not confined to the promoters of growth regulatory genes but is also found in actively transcribed regions.