Olaparib and durvalumab in patients with germline BRCA-mutated metastatic breast cancer (MEDIOLA): an open-label, multicentre, phase 1/2, basket study.

BACKGROUND: Poly (ADP-ribose) polymerase inhibitors combined with immunotherapy have shown antitumour activity in preclinical studies. We aimed to assess the safety and activity of olaparib in combination with the PD-L1-inhibitor, durvalumab, in patients with germline BRCA1-mutated or BRCA2-mutated metastatic breast cancer. METHODS: The MEDIOLA trial is a multicentre, open-label, phase 1/2, basket trial of durvalumab and olaparib in solid tumours. Patients were enrolled into four initial cohorts: germline BRCA-mutated, metastatic breast cancer; germline BRCA-mutated, metastatic ovarian cancer; metastatic gastric cancer; and relapsed small-cell lung cancer. Here, we report on the cohort of patients with breast cancer. Patients who were aged 18 years or older (or aged 19 years or older in South Korea) with germline BRCA1-mutated or BRCA2-mutated or both and histologically confirmed, progressive, HER2-negative, metastatic breast cancer were enrolled from 14 health centres in the UK, the USA, Israel, France, Switzerland, and South Korea. Patients should not have received more than two previous lines of chemotherapy for metastatic breast cancer. Patients received 300 mg olaparib in tablet form orally twice daily for 4 weeks and thereafter a combination of olaparib 300 mg twice daily and durvalumab 1·5 g via intravenous infusion every 4 weeks until disease progression. Primary endpoints were safety and tolerability, and 12-week disease control rate. Safety was analysed in patients who received at least one dose of study treatment, and activity analyses were done in the full-analysis set (patients who received at least one dose of study treatment and were not excluded from the study). Recruitment has completed and the study is ongoing. This trial is registered with ClinicalTrials.gov, NCT02734004. FINDINGS: Between June 14, 2016, and May 2, 2017, 34 patients were enrolled and received both study drugs and were included in the safety analysis. 11 (32%) patients experienced grade 3 or worse adverse events, of which the most common were anaemia (four [12%]), neutropenia (three [9%]), and pancreatitis (two [6%]). Three (9%) patients discontinued due to adverse events and four (12%) patients experienced a total of six serious adverse events. There were no treatment-related deaths. 24 (80%; 90% CI 64·3-90·9) of 30 patients eligible for activity analysis had disease control at 12 weeks. INTERPRETATION: Combination of olaparib and durvalumab showed promising antitumour activity and safety similar to that previously observed in olaparib and durvalumab monotherapy studies. Further research in a randomised setting is needed to determine predictors of therapeutic benefit and whether addition of durvalumab improves long-term clinical outcomes compared with olaparib monotherapy. FUNDING: AstraZeneca.

Circadian clock components control daily growth activities by modulating cytokinin levels and cell division-associated gene expression in Populus trees.

Trees are carbon dioxide sinks and major producers of terrestrial biomass with distinct seasonal growth patterns. Circadian clocks enable the coordination of physiological and biochemical temporal activities, optimally regulating multiple traits including growth. To dissect the clock's role in growth, we analysed Populus tremula × P. tremuloides trees with impaired clock function due to down-regulation of central clock components. late elongated hypocotyl (lhy-10) trees, in which expression of LHY1 and LHY2 is reduced by RNAi, have a short free-running period and show disrupted temporal regulation of gene expression and reduced growth, producing 30-40% less biomass than wild-type trees. Genes important in growth regulation were expressed with an earlier phase in lhy-10, and CYCLIN D3 expression was misaligned and arrhythmic. Levels of cytokinins were lower in lhy-10 trees, which also showed a change in the time of peak expression of genes associated with cell division and growth. However, auxin levels were not altered in lhy-10 trees, and the size of the lignification zone in the stem showed a relative increase. The reduced growth rate and anatomical features of lhy-10 trees were mainly caused by misregulation of cell division, which may have resulted from impaired clock function.

Unbiased analysis of potential targets of breast cancer susceptibility loci by Capture Hi-C.

Genome-wide association studies have identified more than 70 common variants that are associated with breast cancer risk. Most of these variants map to non-protein-coding regions and several map to gene deserts, regions of several hundred kilobases lacking protein-coding genes. We hypothesized that gene deserts harbor long-range regulatory elements that can physically interact with target genes to influence their expression. To test this, we developed Capture Hi-C (CHi-C), which, by incorporating a sequence capture step into a Hi-C protocol, allows high-resolution analysis of targeted regions of the genome. We used CHi-C to investigate long-range interactions at three breast cancer gene deserts mapping to 2q35, 8q24.21, and 9q31.2. We identified interaction peaks between putative regulatory elements ("bait fragments") within the captured regions and "targets" that included both protein-coding genes and long noncoding (lnc) RNAs over distances of 6.6 kb to 2.6 Mb. Target protein-coding genes were IGFBP5, KLF4, NSMCE2, and MYC; and target lncRNAs included DIRC3, PVT1, and CCDC26. For one gene desert, we were able to define two SNPs (rs12613955 and rs4442975) that were highly correlated with the published risk variant and that mapped within the bait end of an interaction peak. In vivo ChIP-qPCR data show that one of these, rs4442975, affects the binding of FOXA1 and implicate this SNP as a putative functional variant.

Whole-exome DNA sequence analysis of Brca2- and Trp53-deficient mouse mammary gland tumours.

Germline mutations in the tumour suppressor BRCA2 predispose to breast, ovarian and a number of other human cancers. Brca2-deficient mouse models are used for preclinical studies but the pattern of genomic alterations in these tumours has not yet been described in detail. We have performed whole-exome DNA sequencing analysis of mouse mammary tumours from Blg-Cre Brca2(f/f) Trp53(f/f) animals, a model of BRCA2-deficient human cancer. We also used the sequencing data to estimate DNA copy number alterations in these tumours and identified a recurrent copy number gain in Met, which has been found amplified in other mouse mammary cancer models. Through a comparative genomic analysis, we identified several mouse Blg-Cre Brca2(f/f) Trp53(f/f) mammary tumour somatic mutations in genes that are also mutated in human cancer, but few of these genes have been found frequently mutated in human breast cancer. A more detailed analysis of these somatic mutations revealed a set of genes that are mutated in human BRCA2 mutant breast and ovarian tumours and that are also mutated in mouse Brca2-null, Trp53-null mammary tumours. Finally, a DNA deletion surrounded by microhomology signature found in human BRCA1/2-deficient cancers was not common in the genome of these mouse tumours. Although a useful model, there are some differences in the genomic landscape of tumours arising in Blg-Cre Brca2(f/f) Trp53(f/f) mice compared to human BRCA-mutated breast cancers. Therefore, this needs to be taken into account in the use of this model.

Genome-wide characterization reveals complex interplay between TP53 and TP63 in response to genotoxic stress.

In response to genotoxic stress the TP53 tumour suppressor activates target gene expression to induce cell cycle arrest or apoptosis depending on the extent of DNA damage. These canonical activities can be repressed by TP63 in normal stratifying epithelia to maintain proliferative capacity or drive proliferation of squamous cell carcinomas, where TP63 is frequently overexpressed/amplified. Here we use ChIP-sequencing, integrated with microarray analysis, to define the genome-wide interplay between TP53 and TP63 in response to genotoxic stress in normal cells. We reveal that TP53 and TP63 bind to overlapping, but distinct cistromes of sites through utilization of distinctive consensus motifs and that TP53 is constitutively bound to a number of sites. We demonstrate that cisplatin and adriamycin elicit distinct effects on TP53 and TP63 binding events, through which TP53 can induce or repress transcription of an extensive network of genes by direct binding and/or modulation of TP63 activity. Collectively, this results in a global TP53-dependent repression of cell cycle progression, mitosis and DNA damage repair concomitant with activation of anti-proliferative and pro-apoptotic canonical target genes. Further analyses reveal that in the absence of genotoxic stress TP63 plays an important role in maintaining expression of DNA repair genes, loss of which results in defective repair.

Whole genome sequencing of matched primary and metastatic acral melanomas.

Next generation sequencing has enabled systematic discovery of mutational spectra in cancer samples. Here, we used whole genome sequencing to characterize somatic mutations and structural variation in a primary acral melanoma and its lymph node metastasis. Our data show that the somatic mutational rates in this acral melanoma sample pair were more comparable to the rates reported in cancer genomes not associated with mutagenic exposure than in the genome of a melanoma cell line or the transcriptome of melanoma short-term cultures. Despite the perception that acral skin is sun-protected, the dominant mutational signature in these samples is compatible with damage due to ultraviolet light exposure. A nonsense mutation in ERCC5 discovered in both the primary and metastatic tumors could also have contributed to the mutational signature through accumulation of unrepaired dipyrimidine lesions. However, evidence of transcription-coupled repair was suggested by the lower mutational rate in the transcribed regions and expressed genes. The primary and the metastasis are highly similar at the level of global gene copy number alterations, loss of heterozygosity and single nucleotide variation (SNV). Furthermore, the majority of the SNVs in the primary tumor were propagated in the metastasis and one nonsynonymous coding SNV and one splice site mutation appeared to arise de novo in the metastatic lesion.

Characterization of the genomic features and expressed fusion genes in micropapillary carcinomas of the breast.

Micropapillary carcinoma (MPC) is a rare histological special type of breast cancer, characterized by an aggressive clinical behaviour and a pattern of copy number aberrations (CNAs) distinct from that of grade- and oestrogen receptor (ER)-matched invasive carcinomas of no special type (IC-NSTs). The aims of this study were to determine whether MPCs are underpinned by a recurrent fusion gene(s) or mutations in 273 genes recurrently mutated in breast cancer. Sixteen MPCs were subjected to microarray-based comparative genomic hybridization (aCGH) analysis and Sequenom OncoCarta mutation analysis. Eight and five MPCs were subjected to targeted capture and RNA sequencing, respectively. aCGH analysis confirmed our previous observations about the repertoire of CNAs of MPCs. Sequencing analysis revealed a spectrum of mutations similar to those of luminal B IC-NSTs, and recurrent mutations affecting mitogen-activated protein kinase family genes and NBPF10. RNA-sequencing analysis identified 17 high-confidence fusion genes, eight of which were validated and two of which were in-frame. No recurrent fusions were identified in an independent series of MPCs and IC-NSTs. Forced expression of in-frame fusion genes (SLC2A1-FAF1 and BCAS4-AURKA) resulted in increased viability of breast cancer cells. In addition, genomic disruption of CDK12 caused by out-of-frame rearrangements was found in one MPC and in 13% of HER2-positive breast cancers, identified through a re-analysis of publicly available massively parallel sequencing data. In vitro analyses revealed that CDK12 gene disruption results in sensitivity to PARP inhibition, and forced expression of wild-type CDK12 in a CDK12-null cell line model resulted in relative resistance to PARP inhibition. Our findings demonstrate that MPCs are neither defined by highly recurrent mutations in the 273 genes tested, nor underpinned by a recurrent fusion gene. Although seemingly private genetic events, some of the fusion transcripts found in MPCs may play a role in maintenance of a malignant phenotype and potentially offer therapeutic opportunities.

Genome-wide functional screen identifies a compendium of genes affecting sensitivity to tamoxifen.

Therapies that target estrogen signaling have made a very considerable contribution to reducing mortality from breast cancer. However, resistance to tamoxifen remains a major clinical problem. Here we have used a genome-wide functional profiling approach to identify multiple genes that confer resistance or sensitivity to tamoxifen. Combining whole-genome shRNA screening with massively parallel sequencing, we have profiled the impact of more than 56,670 RNA interference reagents targeting 16,487 genes on the cellular response to tamoxifen. This screen, along with subsequent validation experiments, identifies a compendium of genes whose silencing causes tamoxifen resistance (including BAP1, CLPP, GPRC5D, NAE1, NF1, NIPBL, NSD1, RAD21, RARG, SMC3, and UBA3) and also a set of genes whose silencing causes sensitivity to this endocrine agent (C10orf72, C15orf55/NUT, EDF1, ING5, KRAS, NOC3L, PPP1R15B, RRAS2, TMPRSS2, and TPM4). Multiple individual genes, including NF1, a regulator of RAS signaling, also correlate with clinical outcome after tamoxifen treatment.

Intraclonal heterogeneity and distinct molecular mechanisms characterize the development of t(4;14) and t(11;14) myeloma.

We have used whole exome sequencing to compare a group of presentation t(4;14) with t(11;14) cases of myeloma to define the mutational landscape. Each case was characterized by a median of 24.5 exonic nonsynonymous single-nucleotide variations, and there was a consistently higher number of mutations in the t(4;14) group, but this number did not reach statistical significance. We show that the transition and transversion rates in the 2 subgroups are similar, suggesting that there was no specific mechanism leading to mutation differentiating the 2 groups. Only 3% of mutations were seen in both groups, and recurrently mutated genes include NRAS, KRAS, BRAF, and DIS3 as well as DNAH5, a member of the axonemal dynein family. The pattern of mutation in each group was distinct, with the t(4;14) group being characterized by deregulation of chromatin organization, actin filament, and microfilament movement. Recurrent RAS pathway mutations identified subclonal heterogeneity at a mutational level in both groups, with mutations being present as either dominant or minor subclones. The presence of subclonal diversity was confirmed at a single-cell level using other tumor-acquired mutations. These results are consistent with a distinct molecular pathogenesis underlying each subgroup and have important impacts on targeted treatment strategies. The Medical Research Council Myeloma IX trial is registered under ISRCTN68454111.

Secondary mutations in BRCA2 associated with clinical resistance to a PARP inhibitor.

PARP inhibitors (PARPi) for the treatment of BRCA1 or BRCA2 deficient tumours are currently the focus of seminal clinical trials exploiting the concept of synthetic lethality. Although clinical resistance to PARPi has been described, the mechanism underlying this has not been elucidated. Here, we investigate tumour material from patients who had developed resistance to the PARPi olaparib, subsequent to showing an initial clinical response. Massively parallel DNA sequencing of treatment-naive and post-olaparib treatment biopsies identified tumour-specific BRCA2 secondary mutations in olaparib-resistant metastases. These secondary mutations restored full-length BRCA2 protein, and most likely cause olaparib resistance by re-establishing BRCA2 function in the tumour cells.

The genomic landscape of oesophagogastric junctional adenocarcinoma.

The incidence of oesophagogastric junctional (OGJ) adenocarcinoma is rising rapidly in western countries, in contrast to the declining frequency of distal gastric carcinoma. Treatment options for adenocarcinomas involving the oesophagogastric junction are limited and the overall prognosis is extremely poor. To determine the genomic landscape of OGJ adenocarcinoma, exomes of eight tumours and matched germline DNA were subjected to massively parallel DNA sequencing. Microsatellite instability was observed in three tumours which coincided with an elevated number of somatic mutations. In total, 117 genes were identified that had predicted coding alterations in more than one tumour. Potentially actionable coding mutations were identified in 67 of these genes, including those in CR2, HGF , FGFR4, and ESRRB. Twenty-nine genes harbouring somatic coding mutations and copy number changes in the MSS OGJ dataset are also known to be altered with similar predicted functional consequence in other tumour types. Compared with the published mutational profile of gastric cancers, 49% (57/117) of recurrently mutated genes were unique to OGJ tumours. TP53, SYNE1, and ARID1A were amongst the most frequently mutated genes in a larger OGJ cohort. Our study provides an insight into the mutational landscape of OGJ adenocarcinomas and confirms that this is a highly mutated and heterogeneous disease. Furthermore, we have uncovered somatic mutations in therapeutically relevant genes which may represent candidate drug targets.

Genome-wide analysis of p63 binding sites identifies AP-2 factors as co-regulators of epidermal differentiation.

The p63 transcription factor (TP63) is critical in development, growth and differentiation of stratifying epithelia. This is highlighted by the severity of congenital abnormalities caused by TP63 mutations in humans, the dramatic phenotypes in knockout mice and de-regulation of TP63 expression in neoplasia altering the tumour suppressive roles of the TP53 family. In order to define the normal role played by TP63 and provide the basis for better understanding how this network is perturbed in disease, we used chromatin immunoprecipitation combined with massively parallel sequencing (ChIP-seq) to identify >7500 high-confidence TP63-binding regions across the entire genome, in primary human neonatal foreskin keratinocytes (HFKs). Using integrative strategies, we demonstrate that only a subset of these sites are bound by TP53 in response to DNA damage. We identify a role for TP63 in transcriptional regulation of multiple genes genetically linked to cleft palate and identify AP-2alpha (TFAP2A) as a co-regulator of a subset of these genes. We further demonstrate that AP-2gamma (TFAP2C) can bind a subset of these regions and that acute depletion of either TFAP2A or TFAP2C alone is sufficient to reduce terminal differentiation of organotypic epidermal skin equivalents, indicating overlapping physiological functions with TP63.

Target-enrichment strategies for next-generation sequencing.

We have not yet reached a point at which routine sequencing of large numbers of whole eukaryotic genomes is feasible, and so it is often necessary to select genomic regions of interest and to enrich these regions before sequencing. There are several enrichment approaches, each with unique advantages and disadvantages. Here we describe our experiences with the leading target-enrichment technologies, the optimizations that we have performed and typical results that can be obtained using each. We also provide detailed protocols for each technology so that end users can find the best compromise between sensitivity, specificity and uniformity for their particular project.

A whole-genome massively parallel sequencing analysis of BRCA1 mutant oestrogen receptor-negative and -positive breast cancers.

BRCA1 encodes a tumour suppressor protein that plays pivotal roles in homologous recombination (HR) DNA repair, cell-cycle checkpoints, and transcriptional regulation. BRCA1 germline mutations confer a high risk of early-onset breast and ovarian cancer. In more than 80% of cases, tumours arising in BRCA1 germline mutation carriers are oestrogen receptor (ER)-negative; however, up to 15% are ER-positive. It has been suggested that BRCA1 ER-positive breast cancers constitute sporadic cancers arising in the context of a BRCA1 germline mutation rather than being causally related to BRCA1 loss-of-function. Whole-genome massively parallel sequencing of ER-positive and ER-negative BRCA1 breast cancers, and their respective germline DNAs, was used to characterize the genetic landscape of BRCA1 cancers at base-pair resolution. Only BRCA1 germline mutations, somatic loss of the wild-type allele, and TP53 somatic mutations were recurrently found in the index cases. BRCA1 breast cancers displayed a mutational signature consistent with that caused by lack of HR DNA repair in both ER-positive and ER-negative cases. Sequencing analysis of independent cohorts of hereditary BRCA1 and sporadic non-BRCA1 breast cancers for the presence of recurrent pathogenic mutations and/or homozygous deletions found in the index cases revealed that DAPK3, TMEM135, KIAA1797, PDE4D, and GATA4 are potential additional drivers of breast cancers. This study demonstrates that BRCA1 pathogenic germline mutations coupled with somatic loss of the wild-type allele are not sufficient for hereditary breast cancers to display an ER-negative phenotype, and has led to the identification of three potential novel breast cancer genes (ie DAPK3, TMEM135, and GATA4).

High altitude adaptation in Daghestani populations from the Caucasus.

We have surveyed 15 high-altitude adaptation candidate genes for signals of positive selection in North Caucasian highlanders using targeted re-sequencing. A total of 49 unrelated Daghestani from three ethnic groups (Avars, Kubachians, and Laks) living in ancient villages located at around 2,000 m above sea level were chosen as the study population. Caucasian (Adygei living at sea level, N = 20) and CEU (CEPH Utah residents with ancestry from northern and western Europe; N = 20) were used as controls. Candidate genes were compared with 20 putatively neutral control regions resequenced in the same individuals. The regions of interest were amplified by long-PCR, pooled according to individual, indexed by adding an eight-nucleotide tag, and sequenced using the Illumina GAII platform. 1,066 SNPs were called using false discovery and false negative thresholds of ~6%. The neutral regions provided an empirical null distribution to compare with the candidate genes for signals of selection. Two genes stood out. In Laks, a non-synonymous variant within HIF1A already known to be associated with improvement in oxygen metabolism was rediscovered, and in Kubachians a cluster of 13 SNPs located in a conserved intronic region within EGLN1 showing high population differentiation was found. These variants illustrate both the common pathways of adaptation to high altitude in different populations and features specific to the Daghestani populations, showing how even a mildly hypoxic environment can lead to genetic adaptation.

Whole genome in vivo RNAi screening identifies the leukemia inhibitory factor receptor as a novel breast tumor suppressor.

Cancer is caused by mutations in oncogenes and tumor suppressor genes, resulting in the deregulation of processes fundamental to the normal behavior of cells. The identification and characterization of oncogenes and tumor suppressors has led to new treatment strategies that have significantly improved cancer outcome. The advent of next generation sequencing has allowed the elucidation of the fine structure of cancer genomes, however, the identification of pathogenic changes is complicated by the inherent genomic instability of cancer cells. Therefore, functional approaches for the identification of novel genes involved in the initiation and development of tumors are critical. Here we report the first whole human genome in vivo RNA interference screen to identify functionally important tumor suppressor genes. Using our novel approach, we identify previously validated tumor suppressor genes including TP53 and MNT, as well as several novel candidate tumor suppressor genes including leukemia inhibitory factor receptor (LIFR). We show that LIFR is a key novel tumor suppressor, whose deregulation may drive the transformation of a significant proportion of human breast cancers. These results demonstrate the power of genome wide in vivo RNAi screens as a method for identifying novel genes regulating tumorigenesis.

Amplification-free Illumina sequencing-library preparation facilitates improved mapping and assembly of (G+C)-biased genomes.

Amplification artifacts introduced during library preparation for the Illumina Genome Analyzer increase the likelihood that an appreciable proportion of these sequences will be duplicates and cause an uneven distribution of read coverage across the targeted sequencing regions. As a consequence, these unfavorable features result in difficulties in genome assembly and variation analysis from the short reads, particularly when the sequences are from genomes with base compositions at the extremes of high or low G+C content. Here we present an amplification-free method of library preparation, in which the cluster amplification step, rather than the PCR, enriches for fully ligated template strands, reducing the incidence of duplicate sequences, improving read mapping and single nucleotide polymorphism calling and aiding de novo assembly. We illustrate this by generating and analyzing DNA sequences from extremely (G+C)-poor (Plasmodium falciparum), (G+C)-neutral (Escherichia coli) and (G+C)-rich (Bordetella pertussis) genomes.

ZEITLUPE Promotes ABA-Induced Stomatal Closure in Arabidopsis and Populus.

Plants balance water availability with gas exchange and photosynthesis by controlling stomatal aperture. This control is regulated in part by the circadian clock, but it remains unclear how signalling pathways of daily rhythms are integrated into stress responses. The serine/threonine protein kinase OPEN STOMATA 1 (OST1) contributes to the regulation of stomatal closure via activation of S-type anion channels. OST1 also mediates gene regulation in response to ABA/drought stress. We show that ZEITLUPE (ZTL), a blue light photoreceptor and clock component, also regulates ABA-induced stomatal closure in Arabidopsis thaliana, establishing a link between clock and ABA-signalling pathways. ZTL sustains expression of OST1 and ABA-signalling genes. Stomatal closure in response to ABA is reduced in ztl mutants, which maintain wider stomatal apertures and show higher rates of gas exchange and water loss than wild-type plants. Detached rosette leaf assays revealed a stronger water loss phenotype in ztl-3, ost1-3 double mutants, indicating that ZTL and OST1 contributed synergistically to the control of stomatal aperture. Experimental studies of Populus sp., revealed that ZTL regulated the circadian clock and stomata, indicating ZTL function was similar in these trees and Arabidopsis. PSEUDO-RESPONSE REGULATOR 5 (PRR5), a known target of ZTL, affects ABA-induced responses, including stomatal regulation. Like ZTL, PRR5 interacted physically with OST1 and contributed to the integration of ABA responses with circadian clock signalling. This suggests a novel mechanism whereby the PRR proteins-which are expressed from dawn to dusk-interact with OST1 to mediate ABA-dependent plant responses to reduce water loss in time of stress.

A large genome center's improvements to the Illumina sequencing system.

The Wellcome Trust Sanger Institute is one of the world's largest genome centers, and a substantial amount of our sequencing is performed with 'next-generation' massively parallel sequencing technologies: in June 2008 the quantity of purity-filtered sequence data generated by our Genome Analyzer (Illumina) platforms reached 1 terabase, and our average weekly Illumina production output is currently 64 gigabases. Here we describe a set of improvements we have made to the standard Illumina protocols to make the library preparation more reliable in a high-throughput environment, to reduce bias, tighten insert size distribution and reliably obtain high yields of data.

Circadian clock components regulate entry and affect exit of seasonal dormancy as well as winter hardiness in Populus trees.

This study addresses the role of the circadian clock in the seasonal growth cycle of trees: growth cessation, bud set, freezing tolerance, and bud burst. Populus tremula x Populus tremuloides (Ptt) LATE ELONGATED HYPOCOTYL1 (PttLHY1), PttLHY2, and TIMING OF CAB EXPRESSION1 constitute regulatory clock components because down-regulation by RNA interference of these genes leads to altered phase and period of clock-controlled gene expression as compared to the wild type. Also, both RNA interference lines show about 1-h-shorter critical daylength for growth cessation as compared to the wild type, extending their period of growth. During winter dormancy, when the diurnal variation in clock gene expression stops altogether, down-regulation of PttLHY1 and PttLHY2 expression compromises freezing tolerance and the expression of C-REPEAT BINDING FACTOR1, suggesting a role of these genes in cold hardiness. Moreover, down-regulation of PttLHY1 and PttLHY2 causes a delay in bud burst. This evidence shows that in addition to a role in daylength-controlled processes, PttLHY plays a role in the temperature-dependent processes of dormancy in Populus such as cold hardiness and bud burst.