Prospective analysis of KRAS, BRAF, and PIK3CA mutations and EGFR copy number in patients (pts) with locally advanced rectal cancer: A translational substudy of a clinical trial (SAKK 41/07) evaluating the effect of neoadjuvant chemoradiation (CRT) with or without panitumumab

Background: Prognostic and predictive markers are of great importance for future study designs and essential for the interpretation of clinical trials incorporating an EGFR-inhibitor. The current study prospectively assessed and validated KRAS, BRAF and PIK3CA mutations in rectal cancer patients screened for the trial SAKK41/07 of concomitant preoperative radio-chemotherapy with or without panitumumab.Methods: Macrodissection was performed on pretreatment formalin fixed paraffin embedded biopsy tissue sections to arrive at a minimum of 50% of tumor cells. DNA was extracted with the Maxwell 16 FFPE Tissue LEV DNA purification kit. After PCR amplification, mutations were identified by pyrosequencing. We prospectively analysed pretreatment biopsy material from 149 rectal cancer pts biopsies for KRAS (exon 2 codon 12 [2-12] and 13 [2-13], exon 3 codon 59 [3-59]) and 61 [3-61], exon 4 codon 117 [4-117] and 146 [4-146]). Sixty-eight pts (KRASwt exon 2, 3 only) were further analysed for BRAF (exon 15 codon 600) and PIK3CA (exon 9 codon 542, 545 and 546, exon 20 codon 1043 [20-1043] and 1047 [20-1047]) mutations, and EGFR copy number by qPCR. For the calculation of the EGFR copy number, we used KRAS copy number as internal reference standard. The calculation was done on the basis of the two standard curves relative quantification method.Results: In 149 screened pts with rectal cancer, the prevalence of KRAS mutations was 36%. Among the 68 pts enrolled in SAKK 41/07 based on initially presumed KRASwt status (exon 2/codons 12+13), 18 pts (26%) had a total of 23 mutations in the RAS/PIK3CA-pathways upon validation analysis. Twelve pts had a KRAS mutation, 7 pts had a PIK3CA mutation, 3 pts had a NRAS mutation, 1 patient a BRAF mutation. Surprisingly, five of these pts had double- mutations, including 4 pts with KRAS plus PIK3CA mutations, and 1 pt with NRAS plus PIK3CA mutations. The median normalized EGFR copy number was 1. Neither mutations of KRAS, BRAF, and PIK3CA, nor EGFR copy number were statistically associated with the primary study endpoint pCR (pathological complete regression).Conclusions: The prevalence of KRAS mutations in rectal and in colon cancer appears to be similar. BRAF mutations are rare; PIK3CA mutations are more common (10%). EGFR copy number is not increased in rectal cancer. A considerable number or KRAS exon 2 wt tumors harbored KRAS exon 3+4 mutations. Further study is needed to determine if KRAS testing should include exons 2-4.

Cell-free reconstitution of vacuole membrane fragmentation reveals regulation of vacuole size and number by TORC1.

Size and copy number of organelles are influenced by an equilibrium of membrane fusion and fission. We studied this equilibrium on vacuoles-the lysosomes of yeast. Vacuole fusion can readily be reconstituted and quantified in vitro, but it had not been possible to study fission of the organelle in a similar way. Here we present a cell-free system that reconstitutes fragmentation of purified yeast vacuoles (lysosomes) into smaller vesicles. Fragmentation in vitro reproduces physiological aspects. It requires the dynamin-like GTPase Vps1p, V-ATPase pump activity, cytosolic proteins, and ATP and GTP hydrolysis. We used the in vitro system to show that the vacuole-associated TOR complex 1 (TORC1) stimulates vacuole fragmentation but not the opposing reaction of vacuole fusion. Under nutrient restriction, TORC1 is inactivated, and the continuing fusion activity then dominates the fusion/fission equilibrium, decreasing the copy number and increasing the volume of the vacuolar compartment. This result can explain why nutrient restriction not only induces autophagy and a massive buildup of vacuolar/lysosomal hydrolases, but also leads to a concomitant increase in volume of the vacuolar compartment by coalescence of the organelles into a single large compartment.

Network-guided analysis of genes with altered somatic copy number and gene expression reveals pathways commonly perturbed in metastatic melanoma.

Cancer genomes frequently contain somatic copy number alterations (SCNA) that can significantly perturb the expression level of affected genes and thus disrupt pathways controlling normal growth. In melanoma, many studies have focussed on the copy number and gene expression levels of the BRAF, PTEN and MITF genes, but little has been done to identify new genes using these parameters at the genome-wide scale. Using karyotyping, SNP and CGH arrays, and RNA-seq, we have identified SCNA affecting gene expression ('SCNA-genes') in seven human metastatic melanoma cell lines. We showed that the combination of these techniques is useful to identify candidate genes potentially involved in tumorigenesis. Since few of these alterations were recurrent across our samples, we used a protein network-guided approach to determine whether any pathways were enriched in SCNA-genes in one or more samples. From this unbiased genome-wide analysis, we identified 28 significantly enriched pathway modules. Comparison with two large, independent melanoma SCNA datasets showed less than 10% overlap at the individual gene level, but network-guided analysis revealed 66% shared pathways, including all but three of the pathways identified in our data. Frequently altered pathways included WNT, cadherin signalling, angiogenesis and melanogenesis. Additionally, our results emphasize the potential of the EPHA3 and FRS2 gene products, involved in angiogenesis and migration, as possible therapeutic targets in melanoma. Our study demonstrates the utility of network-guided approaches, for both large and small datasets, to identify pathways recurrently perturbed in cancer.

Bidirectional shifts in colony queen number in a socially polymorphic ant population.

The breeding system of social organisms affects many important aspects of social life. Some species vary greatly in the number of breeders per group, but the mechanisms and selective pressures contributing to the maintenance of this polymorphism in social structure remain poorly understood. Here, we take advantage of a genetic dataset that spans 15 years to investigate the dynamics of colony queen number within a socially polymorphic ant species. Our study population of Formica selysi has single- and multiple-queen colonies. We found that the social structure of this species is somewhat flexible: on average, each year 3.2% of the single-queen colonies became polygynous, and conversely 1.4% of the multiple-queen colonies became monogynous. The annualized queen replacement rates were 10.3% and 11.9% for single- and multiple-queen colonies, respectively. New queens were often but not always related to previous colony members. At the population level, the social polymorphism appeared stable. There was no genetic differentiation between single- and multiple-queen colonies at eight microsatellite loci, suggesting ongoing gene flow between social forms. Overall, the regular and bidirectional changes in queen number indicate that social structure is a labile trait in F. selysi, with neither form being favored within a time-frame of 15 years.

Estimating the basic reproductive number from viral sequence data.

Epidemiological processes leave a fingerprint in the pattern of genetic structure of virus populations. Here, we provide a new method to infer epidemiological parameters directly from viral sequence data. The method is based on phylogenetic analysis using a birth-death model (BDM) rather than the commonly used coalescent as the model for the epidemiological transmission of the pathogen. Using the BDM has the advantage that transmission and death rates are estimated independently and therefore enables for the first time the estimation of the basic reproductive number of the pathogen using only sequence data, without further assumptions like the average duration of infection. We apply the method to genetic data of the HIV-1 epidemic in Switzerland.

On the number of independent cultural traits carried by individuals and populations.

In species subject to individual and social learning, each individual is likely to express a certain number of different cultural traits acquired during its lifetime. If the process of trait innovation and transmission reaches a steady state in the population, the number of different cultural traits carried by an individual converges to some stationary distribution. We call this the trait-number distribution. In this paper, we derive the trait-number distributions for both individuals and populations when cultural traits are independent of each other. Our results suggest that as the number of cultural traits becomes large, the trait-number distributions approach Poisson distributions so that their means characterize cultural diversity in the population. We then analyse how the mean trait number varies at both the individual and population levels as a function of various demographic features, such as population size and subdivision, and social learning rules, such as conformism and anti-conformism. Diversity at the individual and population levels, as well as at the level of cultural homogeneity within groups, depends critically on the details of population demography and the individual and social learning rules.

Segmental copy number variation shapes tissue transcriptomes.

Copy number variation (CNV) is a key source of genetic diversity, but a comprehensive understanding of its phenotypic effect is only beginning to emerge. We have generated a CNV map in wild mice and classical inbred strains. Genome-wide expression data from six major organs show not only that expression of genes within CNVs tend to correlate with copy number changes, but also that CNVs influence the expression of genes in their vicinity, an effect that extends up to half a megabase. Genes within CNVs show lower expression and more specific spatial expression patterns than genes mapping elsewhere. Our analyses reveal differential constraint on copy number changes of genes expressed in different tissues. Dosage alterations of brain-expressed genes are less frequent than those of other genes and are buffered by tighter transcriptional regulation. Our study provides initial evidence that CNVs shape tissue transcriptomes on a global scale.

Copy number variation characteristics in subpopulations of patients with autism spectrum disorders.

Autism spectrum disorders (ASDs) are a heterogeneous group of disorders with a complex genetic etiology. We used high-resolution whole genome array-based comparative genomic hybridization (array-CGH) to screen 223 ASD patients for gene dose alterations associated with susceptibility for autism. Clinically significant copy number variations (CNVs) were identified in 18 individuals (8%), of which 9 cases (4%) had de novo aberrations. In addition, 20 individuals (9%) were shown to have CNVs of unclear clinical relevance. Among these, 13 cases carried rare but inherited CNVs that may increase the risk for developing ASDs, while parental samples were unavailable in the remaining seven cases. Classification of all patients into different phenotypic and inheritance pattern groups indicated the presence of different CNV patterns in different patient groups. Clinically relevant CNVs were more common in syndromic cases compared to non-syndromic cases. Rare inherited CNVs were present in a higher proportion of ASD cases having first- or second-degree relatives with an ASD-related neuropsychiatric phenotype in comparison with cases without reported heredity (P = 0.0096). We conclude that rare CNVs, encompassing potential candidate regions for ASDs, increase the susceptibility for the development of ASDs and related neuropsychiatric disorders giving us further insight into the complex genetics underlying ASDs.

Genome-wide association analysis of copy number variation in recurrent depressive disorder.

Large, rare copy number variants (CNVs) have been implicated in a variety of psychiatric disorders, but the role of CNVs in recurrent depression is unclear. We performed a genome-wide analysis of large, rare CNVs in 3106 cases of recurrent depression, 459 controls screened for lifetime-absence of psychiatric disorder and 5619 unscreened controls from phase 2 of the Wellcome Trust Case Control Consortium (WTCCC2). We compared the frequency of cases with CNVs against the frequency observed in each control group, analysing CNVs over the whole genome, genic, intergenic, intronic and exonic regions. We found that deletion CNVs were associated with recurrent depression, whereas duplications were not. The effect was significant when comparing cases with WTCCC2 controls (P=7.7 × 10(-6), odds ratio (OR) =1.25 (95% confidence interval (CI) 1.13-1.37)) and to screened controls (P=5.6 × 10(-4), OR=1.52 (95% CI 1.20-1.93). Further analysis showed that CNVs deleting protein coding regions were largely responsible for the association. Within an analysis of regions previously implicated in schizophrenia, we found an overall enrichment of CNVs in our cases when compared with screened controls (P=0.019). We observe an ordered increase of samples with deletion CNVs, with the lowest proportion seen in screened controls, the next highest in unscreened controls and the highest in cases. This may suggest that the absence of deletion CNVs, especially in genes, is associated with resilience to recurrent depression.

Integrated Analysis of High-Resolution Gene Expression and Copy Number Profiling Identified Biallelic Deletion of CDKN2A/2B Tumor Suppressor Locus As the Most Frequent and Unique Genomic Abnormality in Diffuse Large B-Cell Lymphoma (DLBCL) with Strong Prognostic Value in Both GCB and ABC Subtypes and Not Overcome by a Dose-Intensive Immunochemotherapy Regimen Plus Rituximab. Results of a Prospective GELA Clinical Trial Program

Background and aim of the study: Genomic gains and losses play a crucial role in the development and progression of DLBCL and are closely related to gene expression profiles (GEP), including the germinal center B-cell like (GCB) and activated B-cell like (ABC) cell of origin (COO) molecular signatures. To identify new oncogenes or tumor suppressor genes (TSG) involved in DLBCL pathogenesis and to determine their prognostic values, an integrated analysis of high-resolution gene expression and copy number profiling was performed. Patients and methods: Two hundred and eight adult patients with de novo CD20+ DLBCL enrolled in the prospective multicentric randomized LNH-03 GELA trials (LNH03-1B, -2B, -3B, 39B, -5B, -6B, -7B) with available frozen tumour samples, centralized reviewing and adequate DNA/RNA quality were selected. 116 patients were treated by Rituximab(R)-CHOP/R-miniCHOP and 92 patients were treated by the high dose (R)-ACVBP regimen dedicated to patients younger than 60 years (y) in frontline. Tumour samples were simultaneously analysed by high resolution comparative genomic hybridization (CGH, Agilent, 144K) and gene expression arrays (Affymetrix, U133+2). Minimal common regions (MCR), as defined by segments that affect the same chromosomal region in different cases, were delineated. Gene expression and MCR data sets were merged using Gene expression and dosage integrator algorithm (GEDI, Lenz et al. PNAS 2008) to identify new potential driver genes. Results: A total of 1363 recurrent (defined by a penetrance > 5%) MCRs within the DLBCL data set, ranging in size from 386 bp, affecting a single gene, to more than 24 Mb were identified by CGH. Of these MCRs, 756 (55%) showed a significant association with gene expression: 396 (59%) gains, 354 (52%) single-copy deletions, and 6 (67%) homozygous deletions. By this integrated approach, in addition to previously reported genes (CDKN2A/2B, PTEN, DLEU2, TNFAIP3, B2M, CD58, TNFRSF14, FOXP1, REL...), several genes targeted by gene copy abnormalities with a dosage effect and potential physiopathological impact were identified, including genes with TSG activity involved in cell cycle (HACE1, CDKN2C) immune response (CD68, CD177, CD70, TNFSF9, IRAK2), DNA integrity (XRCC2, BRCA1, NCOR1, NF1, FHIT) or oncogenic functions (CD79b, PTPRT, MALT1, AUTS2, MCL1, PTTG1...) with distinct distribution according to COO signature. The CDKN2A/2B tumor suppressor locus (9p21) was deleted homozygously in 27% of cases and hemizygously in 9% of cases. Biallelic loss was observed in 49% of ABC DLBCL and in 10% of GCB DLBCL. This deletion was strongly correlated to age and associated to a limited number of additional genetic abnormalities including trisomy 3, 18 and short gains/losses of Chr. 1, 2, 19 regions (FDR < 0.01), allowing to identify genes that may have synergistic effects with CDKN2A/2B inactivation. With a median follow-up of 42.9 months, only CDKN2A/2B biallelic deletion strongly correlates (FDR p.value < 0.01) to a poor outcome in the entire cohort (4y PFS = 44% [32-61] respectively vs. 74% [66-82] for patients in germline configuration; 4y OS = 53% [39-72] vs 83% [76-90]). In a Cox proportional hazard prediction of the PFS, CDKN2A/2B deletion remains predictive (HR = 1.9 [1.1-3.2], p = 0.02) when combined with IPI (HR = 2.4 [1.4-4.1], p = 0.001) and GCB status (HR = 1.3 [0.8-2.3], p = 0.31). This difference remains predictive in the subgroup of patients treated by R-CHOP (4y PFS = 43% [29-63] vs. 66% [55-78], p=0.02), in patients treated by R-ACVBP (4y PFS = 49% [28-84] vs. 83% [74-92], p=0.003), and in GCB (4y PFS = 50% [27-93] vs. 81% [73-90], p=0.02), or ABC/unclassified (5y PFS = 42% [28-61] vs. 67% [55-82] p = 0.009) molecular subtypes (Figure 1). Conclusion: We report for the first time an integrated genetic analysis of a large cohort of DLBCL patients included in a prospective multicentric clinical trial program allowing identifying new potential driver genes with pathogenic impact. However CDKN2A/2B deletion constitutes the strongest and unique prognostic factor of chemoresistance to R-CHOP, regardless the COO signature, which is not overcome by a more intensified immunochemotherapy. Patients displaying this frequent genomic abnormality warrant new and dedicated therapeutic approaches.

A step beyond--the relevance of depressed mood and mastery in the interplay between the number of social roles and alcohol use.

OBJECTIVES: The present study examines whether depressed mood and external control mediate or moderate the relationship between the number of social roles and alcohol use. PARTICIPANTS: The analysis was based on a national representative sample of 25- to 45-year-old male and female drinkers in Switzerland. METHOD: The influence of depressed mood and external control on the relationship between the number of social roles (parenthood, partnership, employment) and alcohol use was examined in linear structural equation models (mediation) and in multiple regressions (moderation) stratified by gender. All analyses were adjusted for age and education level. RESULTS: Holding more roles was associated with lower alcohol use, lower external control and lower depressed mood. The study did not find evidence of depressed mood or external control mediating the social roles-alcohol relationship. A moderation effect was identified among women only, whereby a protective effect of having more roles could not be found among those who scored high on external control. In general, a stronger link was observed between roles and alcohol use, while depressed mood and external control acted independently on drinking. With the exception of women with high external control, the study found no link between a higher number of social roles and greater alcohol use. CONCLUSION: Our results indicate that drinking behaviours are more strongly linked to external control and depressed mood than they are to the number of social roles. The study also demonstrates that in any effective alcohol prevention policy, societal actions that enable individuals to combine more social roles play a central role.

Prominent use of distal 5' transcription start sites and discovery of a large number of additional exons in ENCODE regions.

This report presents systematic empirical annotation of transcript products from 399 annotated protein-coding loci across the 1% of the human genome targeted by the Encyclopedia of DNA elements (ENCODE) pilot project using a combination of 5' rapid amplification of cDNA ends (RACE) and high-density resolution tiling arrays. We identified previously unannotated and often tissue- or cell-line-specific transcribed fragments (RACEfrags), both 5' distal to the annotated 5' terminus and internal to the annotated gene bounds for the vast majority (81.5%) of the tested genes. Half of the distal RACEfrags span large segments of genomic sequences away from the main portion of the coding transcript and often overlap with the upstream-annotated gene(s). Notably, at least 20% of the resultant novel transcripts have changes in their open reading frames (ORFs), most of them fusing ORFs of adjacent transcripts. A significant fraction of distal RACEfrags show expression levels comparable to those of known exons of the same locus, suggesting that they are not part of very minority splice forms. These results have significant implications concerning (1) our current understanding of the architecture of protein-coding genes; (2) our views on locations of regulatory regions in the genome; and (3) the interpretation of sequence polymorphisms mapping to regions hitherto considered to be "noncoding," ultimately relating to the identification of disease-related sequence alterations.