INTEGRATIVE BIOMARKER IDENTIFICATION AND CLASSIFICATION USING HIGH THROUGHPUT ASSAYS

It is well accepted that tumorigenesis is a multi-step procedure involving aberrant functioning of genes regulating cell proliferation, differentiation, apoptosis, genome stability, angiogenesis and motility. To obtain a full understanding of tumorigenesis, it is necessary to collect information on all aspects of cell activity. Recent advances in high throughput technologies allow biologists to generate massive amounts of data, more than might have been imagined decades ago. These advances have made it possible to launch comprehensive projects such as (TCGA) and (ICGC) which systematically characterize the molecular fingerprints of cancer cells using gene expression, methylation, copy number, microRNA and SNP microarrays as well as next generation sequencing assays interrogating somatic mutation, insertion, deletion, translocation and structural rearrangements. Given the massive amount of data, a major challenge is to integrate information from multiple sources and formulate testable hypotheses. This thesis focuses on developing methodologies for integrative analyses of genomic assays profiled on the same set of samples. We have developed several novel methods for integrative biomarker identification and cancer classification. We introduce a regression-based approach to identify biomarkers predictive to therapy response or survival by integrating multiple assays including gene expression, methylation and copy number data through penalized regression. To identify key cancer-specific genes accounting for multiple mechanisms of regulation, we have developed the integIRTy software that provides robust and reliable inferences about gene alteration by automatically adjusting for sample heterogeneity as well as technical artifacts using Item Response Theory. To cope with the increasing need for accurate cancer diagnosis and individualized therapy, we have developed a robust and powerful algorithm called SIBER to systematically identify bimodally expressed genes using next generation RNAseq data. We have shown that prediction models built from these bimodal genes have the same accuracy as models built from all genes. Further, prediction models with dichotomized gene expression measurements based on their bimodal shapes still perform well. The effectiveness of outcome prediction using discretized signals paves the road for more accurate and interpretable cancer classification by integrating signals from multiple sources.

Single-copy transgenic mice with chosen-site integration.

We describe a general way of introducing transgenes into the mouse germ line for comparing different sequences without the complications of variation in copy number and insertion site. The method uses homologous recombination in embryonic stem (ES) cells to generate mice having a single copy of a transgene integrated into a chosen location in the genome. To test the method, a single copy murine bcl-2 cDNA driven by either a chicken beta-actin promoter or a human beta-actin promoter has been inserted immediately 5' to the X-linked hypoxanthine phosphoribosyltransferase locus by a directly selectable homologous recombination event. The level of expression of the targeted bcl-2 transgene in ES cells is identical in independently isolated homologous recombinants having the same promoter yet varies between the different promoters. In contrast, the expression of bcl-2 transgenes having the same (chicken beta-actin) promoter varies drastically when they are independently integrated at random insertion sites. Both promoters direct broad expression of the single-copy transgene in mice derived from the respective targeted ES cells. In vitro and in vivo, the human beta-actin promoter consistently directed a higher level of transgene expression than the chicken beta-actin promoter.

Copy-up mutants of the plasmid RK2 replication initiation protein are defective in coupling RK2 replication origins.

The broad host range plasmid RK2 replicates and regulates its copy number in a wide range of Gram-negative bacteria. The plasmid-encoded trans-acting replication protein TrfA and the origin of replication oriV are sufficient for controlled replication of the plasmid in all Gram-negative bacteria tested. The TrfA protein binds specifically to direct repeat sequences (iterons) at the origin of replication. A replication control model, designated handcuffing or coupling, has been proposed whereby the formation of coupled TrfA-oriV complexes between plasmid molecules results in hindrance of origin activity and, consequently, a shut-down of plasmid replication under conditions of higher than normal copy number. Therefore, according to this model, the coupling activity of an initiation protein is essential for copy number control and a copy-up initiation protein mutant should have reduced ability to form coupled complexes. To test this model for plasmid RK2, two previously characterized copy-up TrfA mutations, trfA-254D and trfA-267L, were combined and the resulting copy-up double mutant TFrfA protein TrfA-254D/267L was characterized. Despite initiating runaway (uncontrolled) replication in vivo, the copy-up double-mutant TrfA protein exhibited replication kinetics similar to the wild-type protein in vitro. Purified TrfA-254D, TrfA-267L, and TrfA-254D/267L proteins were then examined for binding to the iterons and for coupling activity using an in vitro ligase-catalyzed multimerization assay. It was found that both single and double TrfA mutant proteins exhibited substantially reduced (single mutants) or barely detectable (double mutant) levels of coupling activity while not being diminished in their capacity to bind to the origin of replication. These observations provide direct evidence in support of the coupling model of replication control.

Comparative genomic hybridization analysis of craniopharyngiomas

Object. Craniopharyngioma is the most common childhood brain tumor and is thought to arise from embryonic remnants of the Rathke pouch. Some craniopharyngiomas are monoclonal in origin and hence presumably harbor somatic genetic alterations, although the precise molecular mechanisms involved in craniopharyngioma development are unknown. The goal of this study was to identify genetic alterations in craniopharyngiomas. Methods. To gain insight into the molecular mechanisms involved in development of these tumors, the authors analyzed nine adamantinomatous craniopharyngiomas by using comparative genomic hybridization. Six tumors (67%) displayed at least one genomic alteration, and three had six or more alterations. Only two tumors displayed a decrease in DNA copy number, and in all others an increase in DNA copy number was noted. Conclusions. The authors conclude that a subset of craniopharyngiomas consists of monoclonal tumors arising from activation of oncogenes located at specific chromosomal loci.

Integrated molecular analysis of endometrial carcinoma: translation of biomarker profiles into the clinical practice

The Cancer Genome Atlas (TCGA) collaborative project identified four distinct prognostic groups of endometrial carcinoma (EC) based on molecular alterations: (i) the ultramutated subtype that encompassed POLE mutated (POLE) cases; (ii) the hypermutated subtype, characterized by MisMatch Repair deficiency (MMRd); (iii) the copy-number high subtype, with p53 abnormal/mutated features (p53abn); (iv) the copy-number low subtype, known as No Specific Molecular Profile (NSMP). Although the prognostic value of TCGA molecular classification, NSMP tumors present a wide variability in molecular alterations and biological aggressiveness. This study aims to investigate the impact of ARID1A and CTNNB1/β-catenin alterations by targeted Next-generation sequencing (NGS) and immunohistochemistry (IHC) in a consecutive series of 125 molecularly classified ECs. NGS and IHC were used to assign surrogate TCGA groups and to identify molecular alterations of multiple target genes including POLE, PTEN, ARID1A, CTNNB1, TP53. Associations with clinicopathologic parameters, molecular subtypes, and outcomes identified NSMP category as the most heterogeneous group in terms of clinicopathologic features and outcome. Integration of surrogate TCGA molecular classification with ARID1A and β-catenin analysis showed NSMP cases with ARID1A mutation characterized by the worst outcome with early recurrence, while NSMP tumors with ARID1A wild-type and β-catenin alteration had indolent clinicopathologic features and no recurrence. This study indicates how the identification of ARID1A and β-catenin alterations in EC represents a simple and effective way to characterize NSMP tumor aggressiveness and metastatic potential.

Identification of molecular biomarkers in esophageal adenocarcinoma

Esophageal adenocarcinoma (EAC) is a severe cancer that has been on the rise in Western nations over the past few decades. It has a high mortality rate and the 5-year survival rate is only 35%–45%. EAC has been included in a group of tumors with one of the highest rates of copy number alterations (CNAs), somatic structural rearrangements, high mutation frequency, with different mutational signatures, and with epigenetic mechanisms. The vast heterogeneity of EAC mutations makes it challenging to comprehend the biology that underlies tumor onset and development, identify prognostic biomarkers, and define a molecular classification to stratify patients. The only way to resolve the current disagreements is through an exhaustive molecular analysis of EAC. We examined the genetic profile of 164 patients' esophageal adenocarcinoma samples (without chemo-radiotherapy). The included patients did not receive neoadjuvant therapies, which can change the genetic and molecular composition of the tumor. Using next-generation sequencing technologies (NGS) at high coverage, we examined a custom panel of 26 cancer-related genes. Over the entire cohort, 337 variants were found, with the TP53 gene showing the most frequent alteration (67.27%). Poorer cancer-specific survival was associated with missense mutations in the TP53 gene (Log Rank P=0.0197). We discovered HNF1alpha gene disruptive mutations in 7 cases that were also affected by other gene changes. We started to investigate its role in EAC cell lines by silencing HNF1alpha to mimic our EAC cohort and we use Seahorse technique to analyze its role in the metabolism in esophageal cell. No significant changes were found in transfected cell lines. We conclude by finding that a particular class of TP53 mutations (missense changes) adversely impacted cancer-specific survival in EAC. HNF1alpha, a new EAC-mutated gene, was found, but more research is required to fully understand its function as a tumor suppressor gene.

Distal 22q11.2 Microduplication Combined With Typical 22q11.2 Proximal Deletion: A Case Report.

The 22q11 chromosomal region contains low copy repeats (LCRs) sequences that mediate non-allelic homologous recombination, which predisposes to copy number variations (CNVs) at this locus. Hemizygous deletions of the proximal 22q11.2 region result in the 22q11.2 deletion syndrome (22q11.2 DS). In addition, 22q11.2 duplications involving the distal LCR22s have been reported. This article describes a patient presenting a 2.5-Mb de novo deletion at proximal 22q11.21 region (between LCRs A-D), combined with a 1.3-Mb maternally inherited duplication at distal 22q11.23 region (between LCRs F-H). The presence of concomitant chromosomal imbalances found in this patient has not been reported previously. Clinical and molecular data were compared with literature, in order to contribute to genotype-phenotype correlation. These findings exemplify the complexity and genetic heterogeneity observed in 22q11.2 deletion syndrome and highlights the difficulty to make genetic counseling and predict phenotypic consequences in these situations.

Evaluation of muscle strength and motor abilities in children with type II and III spinal muscle atrophy treated with valproic acid

Background: Spinal muscular atrophy (SMA) is an autosomal recessive disorder that affects the motoneurons of the spinal anterior horn, resulting in hypotonia and muscle weakness. The disease is caused by deletion or mutation in the telomeric copy of SMN gene (SMN1) and clinical severity is in part determined by the copy number of the centromeric copy of the SMN gene (SMN2). The SMN2 mRNA lacks exon 7, resulting in a production of lower amounts of the full-length SMN protein. Knowledge of the molecular mechanism of diseases has led to the discovery of drugs capable of increasing SMN protein level through activation of SMN2 gene. One of these drugs is the valproic acid (VPA), a histone deacetylase inhibitor. Methods: Twenty-two patients with type II and III SMA, aged between 2 and 18 years, were treated with VPA and were evaluated five times during a one-year period using the Manual Muscle Test (Medical Research Council scale-MRC), the Hammersmith Functional Motor Scale (HFMS), and the Barthel Index. Results: After 12 months of therapy, the patients did not gain muscle strength. The group of children with SMA type II presented a significant gain in HFMS scores during the treatment. This improvement was not observed in the group of type III patients. The analysis of the HFMS scores during the treatment period in the groups of patients younger and older than 6 years of age did not show any significant result. There was an improvement of the daily activities at the end of the VPA treatment period. Conclusion: Treatment of SMA patients with VPA may be a potential alternative to alleviate the progression of the disease.

Cuff-induced vascular intima thickening is influenced by titration of the Ace gene in mice

Lacchini S, Heimann AS, Evangelista FS, Cardoso L, Silva GJ, Krieger JE. Cuff-induced vascular intima thickening is influenced by titration of the Ace gene in mice. Physiol Genomics 37: 225-230, 2009. First published March 3, 2009; doi:10.1152/physiolgenomics.90288.2008.-We tested the hypothesis that small changes in angiotensin I-converting enzyme (ACE) expression can alter the vascular response to injury. Male mice containing one, two, three, and four copies of the Ace gene with no detectable vascular abnormality or changes in blood pressure were submitted to cuff-induced femoral artery injury. Femoral thickening was higher in 3- and 4-copy mice (42.4 +/- 4.3% and 45.7 +/- 6.5%, respectively) compared with 1- and 2-copy mice (8.3 +/- 1.3% and 8.5 +/- 0.9%, respectively). Femoral ACE levels from control and injured vessels were assessed in 1- and 3-copy Ace mice, which represent the extremes of the observed response. ACE vascular activity was higher in 3- vs. 1-copy Ace mice (2.4-fold, P < 0.05) in the control uninjured vessel. Upon injury, ACE activity significantly increased in both groups [2.41-fold and 2.14-fold (P < 0.05) for 1- and 3-copy groups, respectively] but reached higher levels in 3- vs. 1-copy Ace mice (P < 0.05). Pharmacological interventions were then used as a counterproof and to indirectly assess the role of angiotensin II (ANG II) on this response. Interestingly, ACE inhibition (enalapril) and ANG II AT(1) receptor blocker (losartan) reduced intima thickening in 3-copy mice to 1-copy mouse values (P < 0.05) while ANG II treatment significantly increased intima thickening in 1-copy mice to 3-copy mouse levels (P < 0.05). Together, these data indicate that small physiologically relevant changes in ACE, not associated with basal vascular abnormalities or blood pressure levels, do influence the magnitude of cuff-induced neointima thickening in mice.

Phenotypic analysis of genes whose mRNA accumulation is dependent on calcineurin in Aspergillus fumigatus

Calcineurin plays an important role in the control of cell morphology and virulence in fungi. Calcineurin is a serine/threonine-specific protein phosphatase heterodimer consisting of a catalytic subunit A and a regulatory subunit B. A mutant of Aspergillus fumigatus lacking the calcineurin A (calA) catalytic subunit exhibited defective hyphal morphology related to apical extension and branching growth, which resulted in drastically decreased filamentation. Here, we investigated which pathways are influenced by A. fumigatus calcineurin during proliferation by comparatively determining the transcriptional profile of A. fumigatus wild type and Delta calA mutant strains. Our results showed that the mitochondrial copy number is reduced in the Delta calA mutant strain, and the mutant has increased alternative oxidase (aoxA) mRNA accumulation and activity. Furthermore, we identified four genes that encode transcription factors that have increased mRNA expression in the Delta calA mutant. Deletion mutants for these transcription factors had reduced susceptibility to itraconazole, caspofungin, and sodium dodecyl sulfate (SDS). (C) 2009 Elsevier Inc. All rights reserved.

Efficient transformation and regeneration of diverse cultivars of peanut (Arachis hypogaea L.) by particle bombardment into embryogenic callus produced from mature seeds

Peanut, one of the world's most important oilseed crops, has a narrow germplasm base and lacks sources of resistance to several major diseases. The species is considered recalcitrant to transformation, with few confirmed transgenic plants upon particle bombardment or Agrobacterium treatment. Reported transformation methods are limited by low efficiency, cultivar specificity, chimeric or infertile transformants, or availability of explants. Here we present a method to efficiently transform cultivars in both botanical types of peanut, by (1) particle bombardment into embryogenic callus derived from mature seeds, (2) escape-free (not stepwise) selection for hygromycin B resistance, (3) brief osmotic desiccation followed by sequential incubation on charcoal and cytokinin-containing media; resulting in efficient conversion of transformed somatic embryos into fertile, non-chimeric, transgenic plants. The method produces three to six independent transformants per bombardment of 10 cm(2) embryogenic callus. Potted, transgenic plant lines can be regenerated within 9 months of callus initiation, or 6 months after bombardment. Transgene copy number ranged from one to 20 with multiple integration sites. There was ca. 50% coexpression of hph and luc or uidA genes coprecipitated on separate plasmids. Reporter gene (luc) expression was confirmed in T-1 progeny from each of six tested independent transformants. Insufficient seeds were produced under containment conditions to determine segregation ratios. The practicality of the technique for efficient cotransformation with selected and unselected genes is demonstrated using major commercial peanut varieties in Australia (cv. NC-7, a virginia market type) and Indonesia (cv. Gajah, a spanish market type).

EDD, the human orthologue of the hyperplastic discs tumour suppressor gene, is amplified and overexpressed in cancer

EDD (E3 isolated by differential display), located at chromosome 8q22.3, is the human orthologue of the Drosophila melanogaster tumour suppressor gene 'hyperplastic discs' and encodes a HECT domain E3 ubiquitin protein-ligase. To investigate the possible involvement of EDD in human cancer, several cancers from diverse tissue sites were analysed for allelic gain or loss (allelic imbalance, AI) at the EDD locus using an EDD-specific microsatellite, CEDD, and other polymorphic microsatellites mapped in the vicinity of the 8q22.3 locus. Of 143 cancers studied, 38 had AI at CEDD (42% of 90 informative cases). In 14 of these cases, discrete regions of imbalance encompassing 8q22.3 were present, while the remainder had more extensive 8q aberrations. AI of CEDD was most frequent in ovarian cancer (22/47 informative cases, 47%), particularly in the serous subtype (16/22, 73%), but was rare in benign and borderline ovarian tumours. AI was also common in breast cancer (31%), hepatocellular carcinoma (46%), squamous cell carcinoma of the tongue (50%) and metastatic melanoma (18%). AI is likely to represent amplification of the EDD gene locus rather than loss of heterozygosity, as quantitative RT-PCR and immunohistochemistry showed that EDD mRNA and protein are frequently overexpressed in breast and ovarian cancers, while among breast cancer cell lines EDD overexpression and increased gene copy number were correlated. These results demonstrate that AI at the EDD locus is common in a diversity of carcinomas and that the EDD gene is frequently overexpressed in breast and ovarian cancer, implying a potential role in cancer progression.

Sequence analysis, chromosomal distribution and long-range organization show that rapid turnover of new and old pBuM satellite DNA repeats leads to different patterns of variation in seven species of the Drosophila buzzatii cluster

We aimed to study patterns of variation and factors influencing the evolutionary dynamics of a satellite DNA, pBuM, in all seven Drosophila species from the buzzatii cluster (repleta group). We analyzed 117 alpha pBuM-1 (monomer length 190 bp) and 119 composite alpha/beta (370 bp) pBuM-2 repeats and determined the chromosome location and long-range organization on DNA fibers of major sequence variants. Such combined methodologies in the study of satDNAs have been used in very few organisms. In most species, concerted evolution is linked to high copy number of pBuM repeats. Species presenting low-abundance and scattered distributed pBuM repeats did not undergo concerted evolution and maintained part of the ancestral inter-repeat variability. The alpha and alpha/beta repeats colocalized in heterochromatic regions and were distributed on multiple chromosomes, with notable differences between species. High-resolution FISH revealed array sizes of a few kilobases to over 0.7 Mb and mutual arrangements of alpha and alpha/beta repeats along the same DNA fibers, but with considerable changes in the amount of each variant across species. From sequence, chromosomal and phylogenetic data, we could infer that homogenization and amplification events involved both new and ancestral pBuM variants. Altogether, the data on the structure and organization of the pBuM satDNA give insights into genome evolution including mechanisms that contribute to concerted evolution and diversification.

Two distinct regions in 2q24.2-q24.3 associated with idiopathic epilepsy

P>Approximately 50% of all carriers of 2q21-q31 deletions present epileptic seizures. The band 2q24 constitutes the smallest commonly deleted segment in these patients, and contains the voltage-gated sodium channel genes SCN1A and SCN2A, associated with Dravet syndrome and benign familial neonatal-infantile seizures, respectively. A further putative locus involving epilepsy in the region was previously identified through disruption of the SLC4A10 gene by translocation. In the course of performing high-resolution DNA copy number analyses on syndromic mentally impaired individuals, we encountered three patients with overlapping deletions in chromosome region 2q24. Two of these patients exhibited epileptic seizures in addition to mental deficiency. The deletion in one of the epileptic patients did not include the SCN cluster, demonstrating that a less severe form of epilepsy maps to an adjacent genomic region. This second region comprises about 3 Mb and contains the candidate gene SLC4A10, providing further support for the potential role of this gene in epilepsy.

CYP1A2*1C, CYP2E1*5B, and GSTM1 polymorphisms are predictors of risk and poor outcome in head and neck squamous cell carcinoma patients

Head and neck squamous cell carcinoma (HNSCC) is associated with environmental factors, especially tobacco and alcohol consumption. Most of the carcinogens present in tobacco smoke are converted into DNA-reactive metabolites by cytochrome P450 (CYPs) enzymes and detoxification of these substances is performed by glutathione S-transferases (GSTs). It has been suggested that genetic alterations, such as polymorphisms, play an important role in tumorigenesis and HNSCC progression. The aim of this study was to investigate CYP1A1, CYP1A2, CYP2E1, GSTM1, and GSTT1 polymorphisms as risk factors in HNSCC and their association with clinicopathologic data. The patients comprised 153 individuals with HNSCC (cases) and 145 with no current or previous diagnosis of cancer (controls). Genotyping of the single nucleotide polymorphisms (SNPs) of the CYP1A1, CYP1A2, and CYP2E1 genes was performed by PCR-RFLP and the GSTM1 and GSTT1 copy number polymorphisms (CNPs) were analyzed by PCR-multiplex. As expected, a significant difference was detected for tobacco and alcohol consumption between cases and controls (P < 0.001). It was observed that the CYP1A2*1D (OR = 16.24) variant and GSTM1 null alleles (OR = 0.02) confer increased risk of HNSCC development (P < 0.001). In addition, head and neck cancer alcohol consumers were more frequently associated with the CYP2E1*5B variant allele than control alcohol users (P < 0.0001, OR = 190.6). The CYP1A2*1C polymorphism was associated with tumor recurrence (log-rank test, P = 0.0161). The CYP2E1*5B and GSTM1 null alleles were significantly associated with advanced clinical stages (T3 + T4; P = 0.022 and P = 0.028, respectively). Overall, the findings suggested that the genetic polymorphisms studied are predictors of risk and are also associated with tumor recurrence, since they are important for determining the parameters associated with tumor progression and poor outcomes in HNSCC. (C) 2009 Elsevier Ltd. All rights reserved.