Optimal clone identifier for genomic shotgun libraries: OC Identifier tool

In DNA microarray experiments, the gene fragments that are spotted on the slides are usually obtained by the synthesis of specific oligonucleotides that are able to amplify genes through PCR. Shotgun library sequences are an alternative to synthesis of primers for the study of each gene in the genome. The possibility of putting thousands of gene sequences into a single slide allows the use of shotgun clones in order to proceed with microarray analysis without a completely sequenced genome. We developed an OC Identifier tool (optimal clone identifier for genomic shotgun libraries) for the identification of unique genes in shotgun libraries based on a partially sequenced genome; this allows simultaneous use of clones in projects such as transcriptome and phylogeny studies, using comparative genomic hybridization and genome assembly. The OC Identifier tool allows comparative genome analysis, biological databases, query language in relational databases, and provides bioinformatics tools to identify clones that contain unique genes as alternatives to primer synthesis. The OC Identifier allows analysis of clones during the sequencing phase, making it possible to select genes of interest for construction of a DNA microarray. ©FUNPEC-RP.

An Integrative Genomic and Transcriptomic Analysis Reveals Potential Targets Associated with Cell Proliferation in Uterine Leiomyomas

Background: Uterine Leiomyomas (ULs) are the most common benign tumours affecting women of reproductive age. ULs represent a major problem in public health, as they are the main indication for hysterectomy. Approximately 40-50% of ULs have non-random cytogenetic abnormalities, and half of ULs may have copy number alterations (CNAs). Gene expression microarrays studies have demonstrated that cell proliferation genes act in response to growth factors and steroids. However, only a few genes mapping to CNAs regions were found to be associated with ULs. Methodology: We applied an integrative analysis using genomic and transcriptomic data to identify the pathways and molecular markers associated with ULs. Fifty-one fresh frozen specimens were evaluated by array CGH (JISTIC) and gene expression microarrays (SAM). The CONEXIC algorithm was applied to integrate the data. Principal Findings: The integrated analysis identified the top 30 significant genes (P<0.01), which comprised genes associated with cancer, whereas the protein-protein interaction analysis indicated a strong association between FANCA and BRCA1. Functional in silico analysis revealed target molecules for drugs involved in cell proliferation, including FGFR1 and IGFBP5. Transcriptional and protein analyses showed that FGFR1 (P = 0.006 and P<0.01, respectively) and IGFBP5 (P = 0.0002 and P = 0.006, respectively) were up-regulated in the tumours when compared with the adjacent normal myometrium. Conclusions: The integrative genomic and transcriptomic approach indicated that FGFR1 and IGFBP5 amplification, as well as the consequent up-regulation of the protein products, plays an important role in the aetiology of ULs and thus provides data for potential drug therapies development to target genes associated with cellular proliferation in ULs. © 2013 Cirilo et al.

Genomic and phenotypic profiles of two Brazilian breast cancer cell lines derived from primary human tumors

Breast cancer is the most common type of cancer among women worldwide. Research using breast cancer cell lines derived from primary tumors may provide valuable additional knowledge regarding this type of cancer. Therefore, the aim of this study was to investigate the phenotypic profiles of MACL-1 and MGSO-3, the only Brazilian breast cancer cell lines available for comparative studies. We evaluated the presence of hormone receptors, proliferation, differentiation and stem cell markers, using immunohistochemical staining of the primary tumor, cultured cells and xenografts implanted in immunodeficient mice. We also investigated the ability of the cell lines to form colonies and copy number alterations by array comparative genomic hybridization. Histopathological analysis showed that the invasive primary tumor from which the MACL-1 cell line was derived, was a luminal A subtype carcinoma, while the ductal carcinoma in situ (DCIS) that gave rise to the MGSO-3 cell line was a HER2 subtype tumor, both showing different proliferation levels. The cell lines and the tumor xenografts in mice preserved their high proliferative potential, but did not maintain the expression of the other markers assessed. This shift in expression may be due to the selection of an 'establishment' phenotype in vitro. Whole-genome DNA evaluation showed a large amount of copy number alterations (CNAs) in the two cell lines. These findings render MACL-1 and MGSO-3 the first characterized Brazilian breast cancer cell lines to be potentially used for comparative research. © 2013 Spandidos Publications Ltd. All rights reserved.

Genomic Signatures Predict Poor Outcome in Undifferentiated Pleomorphic Sarcomas and Leiomyosarcomas

Undifferentiated high-grade pleomorphic sarcomas (UPSs) display aggressive clinical behavior and frequently develop local recurrence and distant metastasis. Because these sarcomas often share similar morphological patterns with other tumors, particularly leiomyosarcomas (LMSs), classification by exclusion is frequently used. In this study, array-based comparative genomic hybridization (array CGH) was used to analyze 20 UPS and 17 LMS samples from untreated patients. The LMS samples presented a lower frequency of genomic alterations compared with the UPS samples. The most frequently altered UPS regions involved gains at 20q13.33 and 7q22.1 and losses at 3p26.3. Gains at 8q24.3 and 19q13.12 and losses at 9p21.3 were frequently detected in the LMS samples. Of these regions, gains at 1q21.3, 11q12.2-q12.3, 16p11.2, and 19q13.12 were significantly associated with reduced overall survival times in LMS patients. A multivariate analysis revealed that gains at 1q21.3 were an independent prognostic marker of shorter survival times in LMS patients (HR = 13.76; P = 0.019). Although the copy number profiles of the UPS and LMS samples could not be distinguished using unsupervised hierarchical clustering analysis, one of the three clusters presented cases associated with poor prognostic outcome (P = 0.022). A relative copy number analysis for the ARNT, SLC27A3, and PBXIP1 genes was performed using quantitative real-time PCR in 11 LMS and 16 UPS samples. Gains at 1q21-q22 were observed in both tumor types, particularly in the UPS samples. These findings provide strong evidence for the existence of a genomic signature to predict poor outcome in a subset of UPS and LMS patients. © 2013 Silveira et al.

Recurrent genomic alterations in sequential progressive leukoplakia and oral cancer: drivers of oral tumorigenesis?

A significant proportion (up to 62) of oral squamous cell carcinomas (OSCCs) may arise from oral potential malignant lesions (OPMLs), such as leukoplakia. Patient outcomes may thus be improved through detection of lesions at a risk for malignant transformation, by identifying and categorizing genetic changes in sequential, progressive OPMLs. We conducted array comparative genomic hybridization analysis of 25 sequential, progressive OPMLs and same-site OSCCs from five patients. Recurrent DNA copy number gains were identified on 1p in 20/25 cases (80) with minimal, high-level amplification regions on 1p35 and 1p36. Other regions of gains were frequently observed: 11q13.4 (68), 9q34.13 (64), 21q22.3 (60), 6p21 and 6q25 (56) and 10q24, 19q13.2, 22q12, 5q31.2, 7p13, 10q24 and 14q22 (48). DNA losses were observed in 20 of samples and mainly detected on 5q31.2 (35), 16p13.2 (30), 9q33.1 and 9q33.29 (25) and 17q11.2, 3p26.2, 18q21.1, 4q34.1 and 8p23.2 (20). Such copy number alterations (CNAs) were mapped in all grades of dysplasia that progressed, and their corresponding OSCCs, in 70 of patients, indicating that these CNAs may be associated with disease progression. Amplified genes mapping within recurrent CNAs (KHDRBS1, PARP1, RAB1A, HBEGF, PAIP2, BTBD7) were selected for validation, by quantitative real-time PCR, in an independent set of 32 progressive leukoplakia, 32 OSSCs and 21 non-progressive leukoplakia samples. Amplification of BTBD7, KHDRBS1, PARP1 and RAB1A was exclusively detected in progressive leukoplakia and corresponding OSCC. BTBD7, KHDRBS1, PARP1 and RAB1A may be associated with OSCC progression. Proteinprotein interaction networks were created to identify possible pathways associated with OSCC progression.

Genomic screening of testicular germ cell tumors from monozygotic twins

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Genomic profiling of human penile carcinoma predicts worse prognosis and survival

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Genomic gains in prostatic carcinoma and proliferative inflammatory atrophy in dogs

The dog can spontaneously develop prostate cancer and consequently can be used as an experimental model for prostatic diseases associated with aging, including benign prostate hyperplasia (BPH) and prostate carcinoma (PCa). DNA copy number variations (CNVs) have been used to identify genes associated with cancer development and progression. DNA microarray based comparative genomic hybridization (aCGH) is a technique that allows to identify copy number of thousands of genes throughout the genome. aCGH was used to identify genomic regions with significantly different DNA copy number in three benign prostatic hyperplasia (BPH), four proliferative inflammatory atrophy (PIA), and 14 canine prostate carcinoma (PCa). Five histologically normal prostate tissue were used as reference. Genomic DNA was extracted from formalin fixed and paraffin embedded samples and CNVs data was evaluated in Canine Genome CGH Microarray 4x44K (G2519F, Design ID021193, Agilent). Data analysis was performed using Genomic Workbench Standard Edition 5.0.14 (Agilent). PCa showed higher number of altered genes related to canonical diseases process, cellular functions and molecular pathways as well as greater inter-relationship between genes, compared with PIA and BPH. In conclusion, PCa showed a more complex genotype, being losses the most frequent genomic changes. Some discrepancies between genomic alterations in human and canine carcinomas may indicate the different clinical behavior of these tumors in these two species. In addition, it was observed was an ascending pattern of genomic complexity in BPH, PIA and CA consistent with a model of multistep tumor progression.

Tumor microenvironmental genomic alterations in juvenile nasopharyngeal angiofibroma

Background To better characterize the pathophysiology of juvenile nasopharyngeal angiofibroma (JNA), endothelial and stromal cells were evaluated by genomic imbalances in association with transcript expression levels of genes mapped on these altered regions. Methods. High-resolution comparative genomic hybridization (HR-CGH) was used in laser-captured endothelial and stromal cells from 9 JNAs. Ten genes were evaluated by quantitative real-timereverse transcription polymerase chain reaction (qRT-PCR) in 15 cases. Results. Although gains were more frequently detected in endothelial cells, 57% of chromosomal alterations were common by both components. Gene expression analyses revealed a positive correlation between endothelial and stromal components for ASPM, CDH1, CTNNB1, FGF18, and SUPT16H. A significant difference was found for FGF18 and AURKB overexpression in stromal cells and AR down-expression in endothelial cells. Conclusions. A similar pattern of gene expression and chromosomal imbalances in both exponents would suggest a common mechanism of functional regulation. AURKB, FGF18, and SUPT16H were identified as potential molecular markers in JNA. (C) 2011 Wiley Periodicals, Inc. Head Neck 34: 485-492, 2012

Xylella fastidiosa comparative genomic database is an information resource to explore the annotation, genomic features, and biology of different strains

The Xylella fastidiosa comparative genomic database is a scientific resource with the aim to provide a user-friendly interface for accessing high-quality manually curated genomic annotation and comparative sequence analysis, as well as for identifying and mapping prophage-like elements, a marked feature of Xylella genomes. Here we describe a database and tools for exploring the biology of this important plant pathogen. The hallmarks of this database are the high quality genomic annotation, the functional and comparative genomic analysis and the identification and mapping of prophage-like elements. It is available from web site http://www.xylella.lncc.br.

Comparative genomic and transcriptome analyses of pathotypes of Xanthomonas citri subsp. citri provide insights into mechanisms of bacterial virulence and host range

Abstract Background Citrus bacterial canker is a disease that has severe economic impact on citrus industries worldwide and is caused by a few species and pathotypes of Xanthomonas. X. citri subsp. citri strain 306 (XccA306) is a type A (Asiatic) strain with a wide host range, whereas its variant X. citri subsp. citri strain Aw12879 (Xcaw12879, Wellington strain) is restricted to Mexican lime. Results To characterize the mechanism for the differences in host range of XccA and Xcaw, the genome of Xcaw12879 that was completed recently was compared with XccA306 genome. Effectors xopAF and avrGf1 are present in Xcaw12879, but were absent in XccA306. AvrGf1 was shown previously for Xcaw to cause hypersensitive response in Duncan grapefruit. Mutation analysis of xopAF indicates that the gene contributes to Xcaw growth in Mexican lime but does not contribute to the limited host range of Xcaw. RNA-Seq analysis was conducted to compare the expression profiles of Xcaw12879 and XccA306 in Nutrient Broth (NB) medium and XVM2 medium, which induces hrp gene expression. Two hundred ninety two and 281 genes showed differential expression in XVM2 compared to in NB for XccA306 and Xcaw12879, respectively. Twenty-five type 3 secretion system genes were up-regulated in XVM2 for both XccA and Xcaw. Among the 4,370 common genes of Xcaw12879 compared to XccA306, 603 genes in NB and 450 genes in XVM2 conditions were differentially regulated. Xcaw12879 showed higher protease activity than XccA306 whereas Xcaw12879 showed lower pectate lyase activity in comparison to XccA306. Conclusions Comparative genomic analysis of XccA306 and Xcaw12879 identified strain specific genes. Our study indicated that AvrGf1 contributes to the host range limitation of Xcaw12879 whereas XopAF contributes to virulence. Transcriptome analyses of XccA306 and Xcaw12879 presented insights into the expression of the two closely related strains of X. citri subsp. citri. Virulence genes including genes encoding T3SS components and effectors are induced in XVM2 medium. Numerous genes with differential expression in Xcaw12879 and XccA306 were identified. This study provided the foundation to further characterize the mechanisms for virulence and host range of pathotypes of X. citri subsp. citri.

Comparative genomic sequencing analysis of a region in human chromosome 11p15.3, mouse chromosome 7 and analysis of the novel STK33, Stk33 gene

The comparative genomic sequence analysis of a region in human chromosome 11p15.3 and its homologous segment in mouse chromosome 7 between ST5 and LMO1 genes has been performed. 158,201 bases were sequenced in the mouse and compared with the syntenic region in human, partially available in the public databases. The analysed region exhibits the typical eukaryotic genomic structure and compared with the close neighbouring regions, strikingly reflexes the mosaic pattern distribution of (G+C) and repeats content despites its relative short size. Within this region the novel gene STK33 was discovered (Stk33 in the mouse), that codes for a serine/threonine kinase. The finding of this gene constitutes an excellent example of the strength of the comparative sequencing approach. Poor gene-predictions in the mouse genomic sequence were corrected and improved by the comparison with the unordered data from the human genomic sequence publicly available. Phylogenetical analysis suggests that STK33 belongs to the calcium/calmodulin-dependent protein kinases group and seems to be a novelty in the chordate lineage. The gene, as a whole, seems to evolve under purifying selection whereas some regions appear to be under strong positive selection. Both human and mouse versions of serine/threonine kinase 33, consists of seventeen exons highly conserved in the coding regions, particularly in those coding for the core protein kinase domain. Also the exon/intron structure in the coding regions of the gene is conserved between human and mouse. The existence and functionality of the gene is supported by the presence of entries in the EST databases and was in vivo fully confirmed by isolating specific transcripts from human uterus total RNA and from several mouse tissues. Strong evidence for alternative splicing was found, which may result in tissue-specific starting points of transcription and in some extent, different protein N-termini. RT-PCR and hybridisation experiments suggest that STK33/Stk33 is differentially expressed in a few tissues and in relative low levels. STK33 has been shown to be reproducibly down-regulated in tumor tissues, particularly in ovarian tumors. RNA in-situ hybridisation experiments using mouse Stk33-specific probes showed expression in dividing cells from lung and germinal epithelium and possibly also in macrophages from kidney and lungs. Preliminary experimentation with antibodies designed in this work, performed in parallel to the preparation of this manuscript, seems to confirm this expression pattern. The fact that the chromosomal region 11p15 in which STK33 is located may be associated with several human diseases including tumor development, suggest further investigation is necessary to establish the role of STK33 in human health.

Genetic analysis using genomic representations

Analysis of the genetic changes in human tumors is often problematical because of the presence of normal stroma and the limited availability of pure tumor DNA. However, large amounts of highly reproducible “representations” of tumor and normal genomes can be made by PCR from nanogram amounts of restriction endonuclease cleaved DNA that has been ligated to oligonucleotide adaptors. We show here that representations are useful for many types of genetic analyses, including measuring relative gene copy number, loss of heterozygosity, and comparative genomic hybridization. Representations may be prepared even from sorted nuclei from fixed and archived tumor biopsies.

Analysis of genomic alterations in benign, atypical, and anaplastic meningiomas: Toward a genetic model of meningioma progression

Nineteen benign [World Health Organization (WHO) grade I; MI], 21 atypical (WHO grade II; MII), and 19 anaplastic (WHO grade III; MIII) sporadic meningiomas were screened for chromosomal imbalances by comparative genomic hybridization (CGH). These data were supplemented by molecular genetic analyses of selected chromosomal regions and genes. With increasing malignancy grade, a marked accumulation of genomic aberrations was observed; i.e., the numbers (mean ± SEM) of total alterations detected per tumor were 2.9 ± 0.7 for MI, 9.2 ± 1.2 for MII, and 13.3 ± 1.9 for MIII. The most frequent alteration detected in MI was loss on 22q (58%). In MII, aberrations most commonly identified were losses on 1p (76%), 22q (71%), 14q (43%), 18q (43%), 10 (38%), and 6q (33%), as well as gains on 20q (48%), 12q (43%), 15q (43%), 1q (33%), 9q (33%), and 17q (33%). In MIII, most of these alterations were found at similar frequencies. However, an increase in losses on 6q (53%), 10 (68%), and 14q (63%) was observed. In addition, 32% of MIII demonstrated loss on 9p. Homozygous deletions in the CDKN2A gene at 9p21 were found in 4 of 16 MIII (25%). Highly amplified DNA sequences were mapped to 12q13–q15 by CGH in 1 MII. Southern blot analysis of this tumor revealed amplification of CDK4 and MDM2. By CGH, DNA sequences from 17q were found to be amplified in 1 MII and 8 MIII, involving 17q23 in all cases. Despite the high frequency of chromosomal aberrations in the MII and MIII investigated, none of these tumors showed mutations in exons 5–8 of the TP53 gene. On the basis of the most common aberrations identified in the various malignancy grades, a model for the genomic alterations associated with meningioma progression is proposed.

Unlocking pathology archives for molecular genetic studies: a reliable method to generate probes for chromogenic and fluorescent in situ hybridization

Chromogenic (CISH) and fluorescent ( FISH) in situ hybridization have emerged as reliable techniques to identify amplifications and chromosomal translocations. CISH provides a spatial distribution of gene copy number changes in tumour tissue and allows a direct correlation between copy number changes and the morphological features of neoplastic cells. However, the limited number of commercially available gene probes has hindered the use of this technique. We have devised a protocol to generate probes for CISH that can be applied to formalin-fixed, paraffin-embedded tissue sections (FFPETS). Bacterial artificial chromosomes ( BACs) containing fragments of human DNA which map to specific genomic regions of interest are amplified with phi 29 polymerase and random primer labelled with biotin. The genomic location of these can be readily confirmed by BAC end pair sequencing and FISH mapping on normal lymphocyte metaphase spreads. To demonstrate the reliability of the probes generated with this protocol, four strategies were employed: (i) probes mapping to cyclin D1 (CCND1) were generated and their performance was compared with that of a commercially available probe for the same gene in a series of 10 FFPETS of breast cancer samples of which five harboured CCND1 amplification; (ii) probes targeting cyclin-dependent kinase 4 were used to validate an amplification identified by microarray-based comparative genomic hybridization (aCGH) in a pleomorphic adenoma; (iii) probes targeting fibroblast growth factor receptor 1 and CCND1 were used to validate amplifications mapping to these regions, as defined by aCGH, in an invasive lobular breast carcinoma with FISH and CISH; and (iv) gene-specific probes for ETV6 and NTRK3 were used to demonstrate the presence of t(12; 15)(p12; q25) translocation in a case of breast secretory carcinoma with dual colour FISH. In summary, this protocol enables the generation of probes mapping to any gene of interest that can be applied to FFPETS, allowing correlation of morphological features with gene copy number.