Data integration from two microarray platforms identifies bi-allelic genetic inactivation of RIC8A in a breast cancer cell line

Background: Using array comparative genomic hybridization (aCGH), a large number of deleted genomic regions have been identified in human cancers. However, subsequent efforts to identify target genes selected for inactivation in these regions have often been challenging. Methods: We integrated here genome-wide copy number data with gene expression data and non-sense mediated mRNA decay rates in breast cancer cell lines to prioritize gene candidates that are likely to be tumour suppressor genes inactivated by bi-allelic genetic events. The candidates were sequenced to identify potential mutations. Results: This integrated genomic approach led to the identification of RIC8A at 11p15 as a putative candidate target gene for the genomic deletion in the ZR-75-1 breast cancer cell line. We identified a truncating mutation in this cell line, leading to loss of expression and rapid decay of the transcript. We screened 127 breast cancers for RIC8A mutations, but did not find any pathogenic mutations. No promoter hypermethylation in these tumours was detected either. However, analysis of gene expression data from breast tumours identified a small group of aggressive tumours that displayed low levels of RIC8A transcripts. qRT-PCR analysis of 38 breast tumours showed a strong association between low RIC8A expression and the presence of TP53 mutations (P = 0.006). Conclusion: We demonstrate a data integration strategy leading to the identification of RIC8A as a gene undergoing a classical double-hit genetic inactivation in a breast cancer cell line, as well as in vivo evidence of loss of RIC8A expression in a subgroup of aggressive TP53 mutant breast cancers.

Array-based gene expression, CGH and tissue data defines a 12q24 gain in neuroblastic tumors with prognostic implication

Neuroblastoma has successfully served as a model system for the identification of neuroectoderm-derived oncogenes. However, in spite of various efforts, only a few clinically useful prognostic markers have been found. Here, we present a framework, which integrates DNA, RNA and tissue data to identify and prioritize genetic events that represent clinically relevant new therapeutic targets and prognostic biomarkers for neuroblastoma.

Integration of miRNA and mRNA expression with DNA copy number of Ewing sarcoma cell lines

Ewing sarcoma is an aggressive and poorly differentiated malignancy of bone and soft tissue. It primarily affects children, adolescents, and young adults, with a slight male predominance. It is characterized by a translocation between chromosomes 11 and 22 resulting in the EWSR1-FLI1fusion transcription factor. The aim of this study is to identify putative Ewing sarcoma target genes through an integrative analysis of three microarray data sets. Array comparative genomic hybridization is used to measure changes in DNA copy number, and analyzed to detect common chromosomal aberrations. mRNA and miRNA microarrays are used to measure expression of protein-coding and miRNA genes, and these results integrated with the copy number data. Chromosomal aberrations typically contain also bystanders in addition to the driving tumor suppressor and oncogenes, and integration with expression helps to identify the true targets. Correlation between expression of miRNAs and their predicted target mRNAs is also evaluated to assess the results of post-transcriptional miRNA regulation on mRNA levels. The highest frequencies of copy number gains were identified in chromosome 8, 1q, and X. Losses were most frequent in 9p21.3, which also showed an enrichment of copy number breakpoints relative to the rest of the genome. Copy number losses in 9p21.3 were found have a statistically significant effect on the expression of MTAP, but not on CDKN2A, which is a known tumor-suppressor in the same locus. MTAP was also down-regulated in the Ewing sarcoma cell lines compared to mesenchymal stem cells. Genes exhibiting elevated expression in association with copy number gains and up-regulation compared to the reference samples included DCAF7, ENO2, MTCP1, andSTK40. Differentially expressed miRNAs were detected by comparing Ewing sarcoma cell lines against mesenchymal stem cells. 21 up-regulated and 32 down-regulated miRNAs were identified, includingmiR-145, which has been previously linked to Ewing sarcoma. The EWSR1-FLI1 fusion gene represses miR-145, which in turn targets FLI1 forming a mutually repressive feedback loop. In addition higher expression linked to copy number gains and compared to mesenchymal stem cells, STK40 was also found to be a target of four different miRNAs that were all down-regulated in Ewing sarcoma cell lines compared to the reference samples. SLCO5A1 was identified as the only up-regulated gene within a frequently gained region in chromosome 8. This region was gained in over 90 % of the cell lines, and also with a higher frequency than the neighboring regions. In addition, SLCO5A1 was found to be a target of three miRNAs that were down-regulated compared to the mesenchymal stem cells.

Genomotipificación de Camylobacter jejuni mediante microarrays de ADN

354 p. (Bibliogr. 271-303) - Correo electrónico de la autora: andrea.guridi@gmail.com

Prioritization of candidate cancer genes-an aid to oncogenomic studies

The development of techniques for oncogenomic analyses such as array comparative genomic hybridization, messenger RNA expression arrays and mutational screens have come to the fore in modern cancer research. Studies utilizing these techniques are able to highlight panels of genes that are altered in cancer. However, these candidate cancer genes must then be scrutinized to reveal whether they contribute to oncogenesis or are coincidental and non-causative. We present a computational method for the prioritization of candidate (i) proto-oncogenes and (ii) tumour suppressor genes from oncogenomic experiments. We constructed computational classifiers using different combinations of sequence and functional data including sequence conservation, protein domains and interactions, and regulatory data. We found that these classifiers are able to distinguish between known cancer genes and other human genes. Furthermore, the classifiers also discriminate candidate cancer genes from a recent mutational screen from other human genes. We provide a web-based facility through which cancer biologists may access our results and we propose computational cancer gene classification as a useful method of prioritizing candidate cancer genes identified in oncogenomic studies.

Pseudohypoparathyroidism Type Ib Associated with Novel Duplications in the GNAS Locus

Context Pseudohypoparathyroidism type 1b (PHP-Ib) is characterized by renal resistance to PTH (and, sometimes, a mild resistance to TSH) and absence of any features of Albright's hereditary osteodystrophy. Patients with PHP-Ib suffer of defects in the methylation pattern of the complex GNAS locus. PHP-Ib can be either sporadic or inherited in an autosomal dominant pattern. Whereas familial PHP-Ib is well characterized at the molecular level, the genetic cause of sporadic PHP-Ib cases remains elusive, although some molecular mechanisms have been associated with this subtype. Objective The aim of the study was to investigate the molecular and imprinting defects in the GNAS locus in two unrelated patients with PHP-Ib. Design We have analyzed the GNAS locus by direct sequencing, Methylation-Specific Multiplex Ligation-dependent Probe Amplification, microsatellites, Quantitative Multiplex PCR of Short Fluorescent fragments and array-Comparative Genomic Hybridization studies in order to characterize two unrelated families with clinical features of PHP-Ib. Results We identified two duplications in the GNAS region in two patients with PHP-Ib: one of them, comprising similar to 320 kb, occurred 'de novo' in the patient, whereas the other one, of similar to 179 kb in length, was inherited from the maternal allele. In both cases, no other known genetic cause was observed. Conclusion In this article, we describe the to-our-knowledge biggest duplications reported so far in the GNAS region. Both are associated to PHP-Ib, one of them occurring 'de novo' and the other one being maternally inherited.

Reciprocal chromosome painting among human, aardvark, and elephant (superorder Afrotheria) reveals the likely eutherian ancestral karyotype

The Afrotheria, a supraordinal grouping of mammals whose radiation is rooted in Africa, is strongly supported by DNA sequence data but not by their disparate anatomical features. We have used flow-sorted human, aardvark, and African elephant chromosome painting probes and applied reciprocal painting schemes to representatives of two of the Afrotherian orders, the Tubulidentata (aardvark) and Proboscidea (elephants), in an attempt to shed additional light on the evolutionary affinities of this enigmatic group of mammals. Although we have not yet found any unique cytogenetic signatures that support the monophyly of the Afrotheria, embedded within the aardvark genome we find the strongest evidence yet of a mammalian ancestral karyotype comprising 2n = 44. This karyotype includes nine chromosomes that show complete conserved synteny to those of man, six that show conservation as single chromosome arms or blocks in the human karyotype but that occur on two different chromosomes in the ancestor, and seven neighbor-joining combinations (i.e., the synteny is maintained in the majority of species of the orders studied so far, but which corresponds to two chromosomes in humans). The comparative chromosome maps presented between human and these Afrotherian species provide further insight into mammalian genome organization and comparative genomic data for the Afrotheria, one of the four major evolutionary clades postulated for the Eutheria.

Genome 10K: A Proposal to Obtain Whole-Genome Sequence for 10 000 Vertebrate Species

The human genome project has been recently complemented by whole-genome assessment sequence of 32 mammals and 24 nonmammalian vertebrate species suitable for comparative genomic analyses. Here we anticipate a precipitous drop in costs and increase in sequ

Molecular characterization of three Rana grylio virus (RGV) isolates and Paralichthys olivaceus lymphocystis disease virus (LCDV-C) in iridoviruses

Three Rana grylio virus (RGV) isolates and lymphocystis disease virus (LCDV-C) were molecularly characterized by antigenicity comparison, Western blot detection of viral polypeptides, restriction fragment length polymorphism analysis of viral genomes, and MCP sequence analysis. Significant antigenicity differences existed among the three RGV isolates and LCDV-C. Western blot detection indicated that the viral polypeptides of three RGV isolates could be recognized by the anti-RGV9807 serum, whereas no bands were observed in the LCDV-C, and significant differences exist among the band patterns of three RGV isolates. Restriction fragment length polymorphism (RFLP) analysis was performed by digesting genomic DNA of the four iridovirus isolates with restriction endonucleases HindIII, KpnI, XbaI and BamHI. On the whole, obvious discrepancies existed between LCDV-C and RGV isolates, and some significant band pattern differences were also revealed between RGV9808 and RGV9506 (or RGV9807) in the profiles of restriction endonucleases Xbal, Kpn I and BamHI. PCR amplification and sequence analysis of MCP gene sequence further revealed their phylogenetic relationship among the three RGV isolates, LCDV-C and other iridoviruses. RGV9506, RGV9807 and RGV9808 are clustered together with other ranaviruses, such as FV3, BIV, TFV and ENHV, although the RGV9808 is more close to EHNV than to other ranaviruses. Additionally, LCDV-C is clustered with LCDV-1, the type species of genus Lymphocystisvirus. The current study provides clear evidence that significant genetic difference exists among the three RGV isolates. Therefore, further work on comparative genomic studies will contribute significantly to understanding of their taxonomic position and pathological mechanism. (C) 2005 Elsevier B.V. All rights reserved.

Construction and characterization of a fosmid library for pathogenic bacterium Vibrio anguillarum

Vibrio anguillarum is a common bacterial pathogen in fish. However, little is known about its pathogenic mechanism, in part, because the entire genome has not been completely sequenced. We constructed a fosmid library for V. anguillarum containing 960 clones with an average insert size of 37.7 kb and 8.6-fold genome coverage. We characterized the library by end-sequencing 50 randomly selected clones. This generated 93 sequences with a total length of 57 485 by covering 1.4% of the whole genome. Of these sequences, 58 (62.4%) were homologous to known genes, 30 (32.3%) were genes with hypothetical functions, and the remaining 5 (5.3%) were unknown genes. We demonstrated the utility of this library by PCR screening of 10 genes. This resulted in an average of 6.2 fosmid clones per screening. This fosmid library offers a new tool for gene screening and cloning of V. anguillarum, and for comparative genomic studies among Vibrio species.

蓝藻信号转导系统的比较基因组学分析及重要基因功能的验证

蓝藻经过漫长的进化，在地球上具有极其广范的分布，而且对环境胁迫具有极强的耐受性，这其中信号转导系统起了至关重要的作用。蓝藻信号转导系统主要包括二元信号转导系统和丝氨酸/苏氨酸激酶。蓝藻全基因组测序工作的完成为蓝藻信号转导系统的研究奠定了基础。目前国外对蓝藻信号转导系统的研究主要集中在对单个蓝藻的单个信号系统的功能研究。本论文首先采用比较基因组学的方法对蓝藻中的两大信号转导系统从分布、结构及进化等角度进行了系统的分析，进而对其中重要信号转导基因的功能进行了验证。 通过比较基因组学分析发现蓝藻中的信号转导系统的分布、结构特征与物种的生理、生态特征之间关系密切；蓝藻中丝氨酸/苏氨酸激酶具有与真核型激酶类似的催化方式；在信号转导系统的进化过程中存在着基因复制、缺失、附属功能域的获得及随机重排等复杂的进化现象；且二元信号转导系统中相互作用的组氨酸激酶与反应调控蛋白之间并不是完全同步的进化关系，可能有着不同的进化历程。该研究建立了蓝藻信号转导系统中的基因-结构-功能的框架结构，为其功能的研究奠定了基础。 聚球藻PCC7942和集胞藻PCC6803是单细胞的淡水蓝藻，具有天然的外源DNA转化系统，是蓝藻分子遗传学研究的模式生物。通过基因突变的方法对这两株蓝藻中二元信号转导系统的部分重要反应调控蛋白同源基因rre28和syn7942_0095的功能进行了验证，发现高度保守的同源序列在这两株蓝藻中起着不同的作用，且一个组氨酸激酶可能与不同的反应调控蛋白相互作用，说明同源序列的功能在进化的过程中已经发生了分化。 目前对蓝藻信号转导系统中丝氨酸/苏氨酸激酶的功能研究较少。通过基因突变及表达差异分析发现集胞藻PCC6803中的丝氨酸/苏氨酸激酶SpkG参与高盐胁迫的信号传递。通过高盐胁迫条件下SpkG对整个转录图谱影响的研究，发现了60个差异表达基因，涉及转运、能量传递、蛋白加工修饰和信号转导等多个生理过程。该研究首次通过实验验证了蓝藻中的丝氨酸/苏氨酸激酶对环境胁迫的响应，并首次发现由二元信号转导系统和丝氨酸/苏氨酸激酶共同参与的对逆境胁迫的调控网络。 蓝藻兼具细菌和植物的特点，蓝藻成熟的转化体系也为真核生物基因功能的研究提供了新的模式宿主。蓝藻信号转导系统中激酶功能的研究为我们进一步研究真核生物激酶的功能提供借鉴，且对蓝藻信号转导的研究将对植物的抗逆胁迫研究提供重要的理论依据。

Origin of avian genome size and structure in non-avian dinosaurs

Chris L. Organ, Andrew M. Shedlock, Andrew Meade, Mark Pagel and Scott V. Edwards (2007). Origin of avian genome size and structure in non-avian dinosaurs. Nature, 46(7132), 180-184. RAE2008

Conservation, duplication, and loss of the Tor signaling pathway in the fungal kingdom.

BACKGROUND: The nutrient-sensing Tor pathway governs cell growth and is conserved in nearly all eukaryotic organisms from unicellular yeasts to multicellular organisms, including humans. Tor is the target of the immunosuppressive drug rapamycin, which in complex with the prolyl isomerase FKBP12 inhibits Tor functions. Rapamycin is a gold standard drug for organ transplant recipients that was approved by the FDA in 1999 and is finding additional clinical indications as a chemotherapeutic and antiproliferative agent. Capitalizing on the plethora of recently sequenced genomes we have conducted comparative genomic studies to annotate the Tor pathway throughout the fungal kingdom and related unicellular opisthokonts, including Monosiga brevicollis, Salpingoeca rosetta, and Capsaspora owczarzaki. RESULTS: Interestingly, the Tor signaling cascade is absent in three microsporidian species with available genome sequences, the only known instance of a eukaryotic group lacking this conserved pathway. The microsporidia are obligate intracellular pathogens with highly reduced genomes, and we hypothesize that they lost the Tor pathway as they adapted and streamlined their genomes for intracellular growth in a nutrient-rich environment. Two TOR paralogs are present in several fungal species as a result of either a whole genome duplication or independent gene/segmental duplication events. One such event was identified in the amphibian pathogen Batrachochytrium dendrobatidis, a chytrid responsible for worldwide global amphibian declines and extinctions. CONCLUSIONS: The repeated independent duplications of the TOR gene in the fungal kingdom might reflect selective pressure acting upon this kinase that populates two proteinaceous complexes with different cellular roles. These comparative genomic analyses illustrate the evolutionary trajectory of a central nutrient-sensing cascade that enables diverse eukaryotic organisms to respond to their natural environments.

Expression in aneuploid Drosophila S2 cells.

Extensive departures from balanced gene dose in aneuploids are highly deleterious. However, we know very little about the relationship between gene copy number and expression in aneuploid cells. We determined copy number and transcript abundance (expression) genome-wide in Drosophila S2 cells by DNA-Seq and RNA-Seq. We found that S2 cells are aneuploid for >43 Mb of the genome, primarily in the range of one to five copies, and show a male genotype ( approximately two X chromosomes and four sets of autosomes, or 2X;4A). Both X chromosomes and autosomes showed expression dosage compensation. X chromosome expression was elevated in a fixed-fold manner regardless of actual gene dose. In engineering terms, the system "anticipates" the perturbation caused by X dose, rather than responding to an error caused by the perturbation. This feed-forward regulation resulted in precise dosage compensation only when X dose was half of the autosome dose. Insufficient compensation occurred at lower X chromosome dose and excessive expression occurred at higher doses. RNAi knockdown of the Male Specific Lethal complex abolished feed-forward regulation. Both autosome and X chromosome genes show Male Specific Lethal-independent compensation that fits a first order dose-response curve. Our data indicate that expression dosage compensation dampens the effect of altered DNA copy number genome-wide. For the X chromosome, compensation includes fixed and dose-dependent components.