RNA-SEQUENCING APPLICATIONS: GENE EXPRESSION QUANTIFICATION AND METHYLATOR PHENOTYPE IDENTIFICATION

My dissertation focuses on two aspects of RNA sequencing technology. The first is the methodology for modeling the overdispersion inherent in RNA-seq data for differential expression analysis. This aspect is addressed in three sections. The second aspect is the application of RNA-seq data to identify the CpG island methylator phenotype (CIMP) by integrating datasets of mRNA expression level and DNA methylation status. Section 1: The cost of DNA sequencing has reduced dramatically in the past decade. Consequently, genomic research increasingly depends on sequencing technology. However it remains elusive how the sequencing capacity influences the accuracy of mRNA expression measurement. We observe that accuracy improves along with the increasing sequencing depth. To model the overdispersion, we use the beta-binomial distribution with a new parameter indicating the dependency between overdispersion and sequencing depth. Our modified beta-binomial model performs better than the binomial or the pure beta-binomial model with a lower false discovery rate. Section 2: Although a number of methods have been proposed in order to accurately analyze differential RNA expression on the gene level, modeling on the base pair level is required. Here, we find that the overdispersion rate decreases as the sequencing depth increases on the base pair level. Also, we propose four models and compare them with each other. As expected, our beta binomial model with a dynamic overdispersion rate is shown to be superior. Section 3: We investigate biases in RNA-seq by exploring the measurement of the external control, spike-in RNA. This study is based on two datasets with spike-in controls obtained from a recent study. We observe an undiscovered bias in the measurement of the spike-in transcripts that arises from the influence of the sample transcripts in RNA-seq. Also, we find that this influence is related to the local sequence of the random hexamer that is used in priming. We suggest a model of the inequality between samples and to correct this type of bias. Section 4: The expression of a gene can be turned off when its promoter is highly methylated. Several studies have reported that a clear threshold effect exists in gene silencing that is mediated by DNA methylation. It is reasonable to assume the thresholds are specific for each gene. It is also intriguing to investigate genes that are largely controlled by DNA methylation. These genes are called “L-shaped” genes. We develop a method to determine the DNA methylation threshold and identify a new CIMP of BRCA. In conclusion, we provide a detailed understanding of the relationship between the overdispersion rate and sequencing depth. And we reveal a new bias in RNA-seq and provide a detailed understanding of the relationship between this new bias and the local sequence. Also we develop a powerful method to dichotomize methylation status and consequently we identify a new CIMP of breast cancer with a distinct classification of molecular characteristics and clinical features.

BAYESIAN STATISTICAL METHODS IN GENE-ENVIRONMENT AND GENE-GENE INTERACTION STUDIES

Complex diseases such as cancer result from multiple genetic changes and environmental exposures. Due to the rapid development of genotyping and sequencing technologies, we are now able to more accurately assess causal effects of many genetic and environmental factors. Genome-wide association studies have been able to localize many causal genetic variants predisposing to certain diseases. However, these studies only explain a small portion of variations in the heritability of diseases. More advanced statistical models are urgently needed to identify and characterize some additional genetic and environmental factors and their interactions, which will enable us to better understand the causes of complex diseases. In the past decade, thanks to the increasing computational capabilities and novel statistical developments, Bayesian methods have been widely applied in the genetics/genomics researches and demonstrating superiority over some regular approaches in certain research areas. Gene-environment and gene-gene interaction studies are among the areas where Bayesian methods may fully exert its functionalities and advantages. This dissertation focuses on developing new Bayesian statistical methods for data analysis with complex gene-environment and gene-gene interactions, as well as extending some existing methods for gene-environment interactions to other related areas. It includes three sections: (1) Deriving the Bayesian variable selection framework for the hierarchical gene-environment and gene-gene interactions; (2) Developing the Bayesian Natural and Orthogonal Interaction (NOIA) models for gene-environment interactions; and (3) extending the applications of two Bayesian statistical methods which were developed for gene-environment interaction studies, to other related types of studies such as adaptive borrowing historical data. We propose a Bayesian hierarchical mixture model framework that allows us to investigate the genetic and environmental effects, gene by gene interactions (epistasis) and gene by environment interactions in the same model. It is well known that, in many practical situations, there exists a natural hierarchical structure between the main effects and interactions in the linear model. Here we propose a model that incorporates this hierarchical structure into the Bayesian mixture model, such that the irrelevant interaction effects can be removed more efficiently, resulting in more robust, parsimonious and powerful models. We evaluate both of the 'strong hierarchical' and 'weak hierarchical' models, which specify that both or one of the main effects between interacting factors must be present for the interactions to be included in the model. The extensive simulation results show that the proposed strong and weak hierarchical mixture models control the proportion of false positive discoveries and yield a powerful approach to identify the predisposing main effects and interactions in the studies with complex gene-environment and gene-gene interactions. We also compare these two models with the 'independent' model that does not impose this hierarchical constraint and observe their superior performances in most of the considered situations. The proposed models are implemented in the real data analysis of gene and environment interactions in the cases of lung cancer and cutaneous melanoma case-control studies. The Bayesian statistical models enjoy the properties of being allowed to incorporate useful prior information in the modeling process. Moreover, the Bayesian mixture model outperforms the multivariate logistic model in terms of the performances on the parameter estimation and variable selection in most cases. Our proposed models hold the hierarchical constraints, that further improve the Bayesian mixture model by reducing the proportion of false positive findings among the identified interactions and successfully identifying the reported associations. This is practically appealing for the study of investigating the causal factors from a moderate number of candidate genetic and environmental factors along with a relatively large number of interactions. The natural and orthogonal interaction (NOIA) models of genetic effects have previously been developed to provide an analysis framework, by which the estimates of effects for a quantitative trait are statistically orthogonal regardless of the existence of Hardy-Weinberg Equilibrium (HWE) within loci. Ma et al. (2012) recently developed a NOIA model for the gene-environment interaction studies and have shown the advantages of using the model for detecting the true main effects and interactions, compared with the usual functional model. In this project, we propose a novel Bayesian statistical model that combines the Bayesian hierarchical mixture model with the NOIA statistical model and the usual functional model. The proposed Bayesian NOIA model demonstrates more power at detecting the non-null effects with higher marginal posterior probabilities. Also, we review two Bayesian statistical models (Bayesian empirical shrinkage-type estimator and Bayesian model averaging), which were developed for the gene-environment interaction studies. Inspired by these Bayesian models, we develop two novel statistical methods that are able to handle the related problems such as borrowing data from historical studies. The proposed methods are analogous to the methods for the gene-environment interactions on behalf of the success on balancing the statistical efficiency and bias in a unified model. By extensive simulation studies, we compare the operating characteristics of the proposed models with the existing models including the hierarchical meta-analysis model. The results show that the proposed approaches adaptively borrow the historical data in a data-driven way. These novel models may have a broad range of statistical applications in both of genetic/genomic and clinical studies.

Patenting genes and gene sequences: Analysis of public health implications and the law

At issue is whether or not isolated DNA is patent eligible under the U.S. Patent Law and the implications of that determination on public health. The U.S. Patent and Trademark Office has issued patents on DNA since the 1980s, and scientists and researchers have proceeded under that milieu since that time. Today, genetic research and testing related to the human breast cancer genes BRCA1 and BRCA2 is conducted within the framework of seven patents that were issued to Myriad Genetics and the University of Utah Research Foundation between 1997 and 2000. In 2009, suit was filed on behalf of multiple researchers, professional associations and others to invalidate fifteen of the claims underlying those patents. The Court of Appeals for the Federal Circuit, which hears patent cases, has invalidated claims for analyzing and comparing isolated DNA but has upheld claims to isolated DNA. The specific issue of whether isolated DNA is patent eligible is now before the Supreme Court, which is expected to decide the case by year's end. In this work, a systematic review was performed to determine the effects of DNA patents on various stakeholders and, ultimately, on public health; and to provide a legal analysis of the patent eligibility of isolated DNA and the likely outcome of the Supreme Court's decision. ^ A literature review was conducted to: first, identify principle stakeholders with an interest in patent eligibility of the isolated DNA sequences BRCA1 and BRCA2; and second, determine the effect of the case on those stakeholders. Published reports that addressed gene patents, the Myriad litigation, and implications of gene patents on stakeholders were included. Next, an in-depth legal analysis of the patent eligibility of isolated DNA and methods for analyzing it was performed pursuant to accepted methods of legal research and analysis based on legal briefs, federal law and jurisprudence, scholarly works and standard practice legal analysis. ^ Biotechnology, biomedical and clinical research, access to health care, and personalized medicine were identified as the principle stakeholders and interests herein. Many experts believe that the patent eligibility of isolated DNA will not greatly affect the biotechnology industry insofar as genetic testing is concerned; unlike for therapeutics, genetic testing does not require tremendous resources or lead time. The actual impact on biomedical researchers is uncertain, with greater impact expected for researchers whose work is intended for commercial purposes (versus basic science). The impact on access to health care has been surprisingly difficult to assess; while invalidating gene patents might be expected to decrease the cost of genetic testing and improve access to more laboratories and physicians' offices that provide the test, a 2010 study on the actual impact was inconclusive. As for personalized medicine, many experts believe that the availability of personalized medicine is ultimately a public policy issue for Congress, not the courts. ^ Based on the legal analysis performed in this work, this writer believes the Supreme Court is likely to invalidate patents on isolated DNA whose sequences are found in nature, because these gene sequences are a basic tool of scientific and technologic work and patents on isolated DNA would unduly inhibit their future use. Patents on complementary DNA (cDNA) are expected to stand, however, based on the human intervention required to craft cDNA and the product's distinction from the DNA found in nature. ^ In the end, the solution as to how to address gene patents may lie not in jurisprudence but in a fundamental change in business practices to provide expanded licenses to better address the interests of the several stakeholders. ^

The transcriptional regulation of the {\it neu\/} gene: Examples of activation, repression, and cell-specific expression

The neu oncogene encodes a growth factor receptor-like protein, p185, with an intrinsic tyrosine kinase activity. A single point mutation, an A to T transversion resulting in an amino acid substitution from valine to glutamic acid, in the transmembrane domain of the rat neu gene was found to be responsible for the transforming and tumorigenic phenotype of the cells that carry it. In contrast, the human proto-neu oncogene is frequently amplified in tumors and cell lines derived from tumors and the human neu gene overexpression/amplification in breast and ovarian cancers is known to correlate with poor patient prognosis. Examples of the human neu gene overexpression in the absence of gene amplification have been observed, which may suggest the significant role of the transcriptional and/or post-transcriptional control of the neu gene in the oncogenic process. However, little is known about the transcriptional mechanisms which regulate the neu gene expression. In this study, three examples are presented to demonstrate the positive and negative control of the neu gene expression.^ First, by using band shift assays and methylation interference analyses, I have identified a specific protein-binding sequence, AAGATAAAACC ($-$466 to $-$456), that binds a specific trans-acting factor termed RVF (for EcoRV factor on the neu promoter). The RVF-binding site is required for maximum transcriptional activity of the rat neu promoter. This same sequence is also found in the corresponding regions of both human and mouse neu promoters. Furthermore, this sequence can enhance the CAT activity driven by a minimum promoter of the thymidine kinase gene in an orientation-independent manner, and thus it behaves as an enhancer. In addition, Southwestern (DNA-protein) blot analysis using the RVF-binding site as a probe points to a 60-kDa polypeptide as a potential candidate for RVF.^ Second, it has been reported that the E3 region of adenovirus 5 induces down-regulation of epidermal growth factor (EGF) receptor through endocytosis. I found that the human neu gene product, p185, (an EGF receptor-related protein) is also down-regulated by adenovirus 5, but via a different mechanism. I demonstrate that the adenovirus E1a gene is responsible for the repression of the human neu gene at the transcriptional level.^ Third, a differential expression of the neu gene has been found in two cell model systems: between the mouse fibroblast Swiss-Webster 3T3 (SW3T3) and its variant NR-6 cells; and between the mouse liver tumor cell line, Hep1-a, and the mouse pancreas tumor cell line, 266-6. Both NR-6 and 266-6 cell lines are not able to express the neu gene product, p185. I demonstrate that, in both cases, the transcriptional repression of the neu gene may account for the lack of the p185 expression in these two cell lines. ^

Alterations in the ras oncogene and the p53 tumor suppressor gene in ultraviolet radiation-induced skin carcinogenesis

A rapid increase of the ultraviolet radiation (UVR)-related skin cancer incidence has attracted more and more public attention during the last few decades. Prevention and treatment of UVR-related skin cancer has become an important public health issue in the United States. Recent studies indicate that mutations in ras and/or p53 genes may be involved in UVR-induced skin tumor development but the precise molecular mechanism remains unclear. In this study, alterations of H-ras and p53 genes were investigated in different stages of carcinogenesis in a chronic UVR (solar simulator) exposure-induced Sencar mouse skin carcinogenesis model in order to clarify the role of the alterations of these genes during the skin carcinogenesis process and to further understand the mechanisms by which UVR causes skin cancer.^ Positive ras-p21 staining in cell membranes and cytosol were detected in 18/33 (55%) of squamous cell carcinomas (SCCs), but were not detected in UV-exposed skin, papillomas, or spindle cell tumors (SCTs). Positive staining of the malignant progression marker K13 was found in 17/33 (52%) of SCCs only. A significant positive correlation was observed between the K13 and the ras-p21 expression. Polymerase chain reaction (PCR)-based single strand conformation polymorphism (SSCP) analysis and gene sequencing analysis revealed three point mutations, one (codon 56) in UV-exposed non-tumor bearing skin and the other two (codons 21 and 13) in SCCs. No UV-specific mutation patterns were found.^ Positive p53 nuclear staining was found in 10/37 (27%) of SCCs and 12/24 (50%) of SCTs, but was not detected in normal skin or papillomas. PCR-based SSCP and sequencing analysis revealed eight point mutations in exons 5 and 6 (four in SCTs, two in SCCs, and two in UV-exposed skin) including six C-T or C-A transitions. Four of the mutations occurred at a dipyrimidine (CC) sequence. The pattern of the mutations indicated that the mutagenic lesions were induced by UVR.^ These results indicate that overexpression of ras-p21 in conjunction with aberrant expression of K13 occurred frequently in UVR-induced SCCs in Sencar mouse skin. The point mutation in the H-ras gene appeared to be a rare event in UVR skin carcinogenesis and may not be responsible for overexpression of ras-p21. UVR-induced P53 gene alteration is a frequent event in UVR-induced SCCs and later stage SCT tumors in Sencar mice skin, suggesting the p53 gene mutation plays an important role in skin tumor malignant progression. ^

Functional analysis of SOX9 in chondrogenesis by targeted gene inactivation

Sox9 is a Sry-related HMG-domain containing transcription factor. Lines of evidence suggest that Sox9 has a potential role in skeletal development. During mouse development, Sox9 is predominantly expressed in all chondroprogenitors and differentiated chondrocytes, throughout the deposition of cartilage matrix. Mutations in one allele of SOX9 in humans result in campomelic dysplasia (CD), a skeletal dysplasia. syndrome characterized by the bowing of long bones. Moreover, Sox9 binds to and activates chondrocyte-specific enhancers in Col2a1 and Col11a2 genes. To further investigate the function of Sox9 in chondrogenesis, we analyzed chimeras derived from Sox9 heterozygous and homozygous null embryonic stem (ES) cells. In mouse chimeras, Sox9 −/− cells were excluded from all cartilages and did not express chondrocyte-specific genes. The segregation occurred during mesenchymal condensation. No cartilages developed in teratocarcinomas derived from Sox9 −/− ES cells. Mice heterozygous for the Sox9 mutation died neonatally and exhibited skeletal abnormalities resembling those of the CD patients. The Sox9 +/− mutants had a cleft palate and hypoplasia of scapula, pelvis and other skeletal structures derived by endochondral ossification. Bending of the radius, ulna and tibia cartilage was prominent at embryonic day 14.5 (E14.5). At E12.5 many pre-cartilaginous condensations were already defective. Advanced ossification was observed and the hypertrophic zone was enlarged in the growth plates, suggesting that Sox9 also regulates hypertrophic chondrocyte differentiation. Our results identify Sox9 as the first essential regulator of chondrocyte differentiation, which plays multiple roles in chondrogenesis. ^

Gene expression and physiological changes of the Antarctic krill Euphausia superba under different hypoxia intensities

Antarctic krill (Euphausia superba) from South Georgia comprise one of the most northern and abundant krill stocks. South Georgia waters are undergoing rapid warming, as a result of climate change, which in turn could alter the oxygen concentration of the water. We investigated gene expression in Antarctic krill related to aerobic metabolism, antioxidant defence, and heat-shock response under severe (2.5% O2 saturation or 0.6 kPa) and threshold (20% O2 saturation or 4 kPa) hypoxia exposure compared to in situ levels (normoxic; 100% O2 saturation or 21 kPa). Biochemical metabolic and oxidative stress indicators complemented the genic expression analysis to detect in vivo signs of stress during the hypoxia treatments. Expression levels of the genes citrate synthase (CS), mitochondrial manganese superoxide dismutase (SODMn-m) and one heat-shock protein isoform (E) were higher in euphausiids incubated 6 h at 20% O2 saturation than in animals exposed to control (normoxic) conditions. All biochemical antioxidant defence parameters remained unchanged among treatments. Levels of lipid peroxidation were raised after 6 h of severe hypoxia. Overall, short-term exposure to hypoxia altered mitochondrial metabolic and antioxidant capacity, but did not induce anaerobic metabolism. Antarctic krill are swarming organisms and may experience short periods of hypoxia when present in dense swarms. A future, warmer Southern ocean, where oxygen saturation levels are decreased, may result in smaller, less dense swarms as they act to avoid greater levels of hypoxia.

Gene expression profiling in gills of the great spider crab Hyas araneus in respond to ocean acidification and warming, supplementary data

Hypercapnia and elevated temperatures resulting from climate change may have adverse consequences for many marine organisms. While diverse physiological and ecological effects have been identified, changes in those molecular mechanisms, which shape the physiological phenotype of a species and limit its capacity to compensate, remain poorly understood. Here, we use global gene expression profiling through RNA-Sequencing to study the transcriptional responses to ocean acidification and warming in gills of the boreal spider crab Hyas araneus exposed medium-term (10 weeks) to intermediate (1,120 µatm) and high (1,960 µatm) PCO2 at different temperatures (5°C and 10°C). The analyses reveal shifts in steady state gene expression from control to intermediate and from intermediate to high CO2 exposures. At 5°C acid-base, energy metabolism and stress response related genes were upregulated at intermediate PCO2, whereas high PCO2 induced a relative reduction in expression to levels closer to controls. A similar pattern was found at elevated temperature (10°C). There was a strong coordination between acid-base, metabolic and stress-related processes. Hemolymph parameters at intermediate PCO2 indicate enhanced capacity in acid-base compensation potentially supported by upregulation of a V-ATPase. The likely enhanced energy demand might be met by the upregulation of the electron transport system (ETS), but may lead to increased oxidative stress reflected in upregulated antioxidant defense transcripts. These mechanisms were attenuated by high PCO2, possibly as a result of limited acid-base compensation and metabolic down-regulation. Our findings indicate a PCO2 dependent threshold beyond which compensation by acclimation fails progressively. They also indicate a limited ability of this stenoecious crustacean to compensate for the effects of ocean acidification with and without concomitant warming.

Gene Expression Profiling in Limb-Girdle Muscular Dystrophy 2A

Limb-girdle muscular dystrophy type 2A (LGMD2A) is a recessive genetic disorder caused by mutations in calpain 3 (CAPN3). Calpain 3 plays different roles in muscular cells, but little is known about its functions or in vivo substrates. The aim of this study was to identify the genes showing an altered expression in LGMD2A patients and the possible pathways they are implicated in. Ten muscle samples from LGMD2A patients with in which molecular diagnosis was ascertained were investigated using array technology to analyze gene expression profiling as compared to ten normal muscle samples. Upregulated genes were mostly those related to extracellular matrix (different collagens), cell adhesion (fibronectin), muscle development (myosins and melusin) and signal transduction. It is therefore suggested that different proteins located or participating in the costameric region are implicated in processes regulated by calpain 3 during skeletal muscle development. Genes participating in the ubiquitin proteasome degradation pathway were found to be deregulated in LGMD2A patients, suggesting that regulation of this pathway may be under the control of calpain 3 activity. As frizzled-related protein (FRZB) is upregulated in LGMD2A muscle samples, it could be hypothesized that β-catenin regulation is also altered at the Wnt signaling pathway, leading to an incorrect myogenesis. Conversely, expression of most transcription factor genes was downregulated (MYC, FOS and EGR1). Finally, the upregulation of IL-32 and immunoglobulin genes may induce the eosinophil chemoattraction explaining the inflammatory findings observed in presymptomatic stages. The obtained results try to shed some light on identification of novel therapeutic targets for limb-girdle muscular dystrophies

Fusarium proliferatum isolated from garlic in Spain: Identification, toxigenic potential and pathogenicity on related Allium species

Fusarium proliferatum has been reported on garlic in the Northwest USA, Spain and Serbia, causing water-soaked tan-colored lesions on cloves. In this work, Fusarium proliferatum was isolated from 300 symptomatic garlic bulbs. Morphological identification of Fusarium was confirmed using species-specific PCR assays and EF-1α sequencing. Confirmation of pathogenicity was conducted with eighteen isolates. Six randomly selected F. proliferatum isolates from garlic were tested for specific pathogenicity and screened for fusaric acid production. Additionally, pathogenicity of each F. proliferatum isolate was tested on healthy seedlings of onion (Allium cepa), leek (A. porrum), scallions (A. fistulosum), chives (A. schoenoprasum) and garlic (A. sativum). A disease severity index (DSI) was calculated as the mean severity on three plants of each species with four test replicates. Symptoms on onion and garlic plants were observed three weeks after inoculation. All isolates tested produced symptoms on all varieties inoculated. Inoculation of F. proliferatum isolates from diseased garlic onto other Allium species provided new information on host range and pathogenicity. The results demonstrated differences in susceptibility with respect to host species and cultivar. The F. proliferatum isolates tested all produced fusaric acid (FA); correlations between FA production and isolate pathogenicity are discussed. Additionally, all isolates showed the presence of the FUM1 gene suggesting the ability of Spanish isolates to produce fumonisins.

MedVir: 3D visual interface applied to gene profile analysis

The use of data mining techniques for the gene profile discovery of diseases, such as cancer, is becoming usual in many researches. These techniques do not usually analyze the relationships between genes in depth, depending on the different variety of manifestations of the disease (related to patients). This kind of analysis takes a considerable amount of time and is not always the focus of the research. However, it is crucial in order to generate personalized treatments to fight the disease. Thus, this research focuses on finding a mechanism for gene profile analysis to be used by the medical and biologist experts. Results: In this research, the MedVir framework is proposed. It is an intuitive mechanism based on the visualization of medical data such as gene profiles, patients, clinical data, etc. MedVir, which is based on an Evolutionary Optimization technique, is a Dimensionality Reduction (DR) approach that presents the data in a three dimensional space. Furthermore, thanks to Virtual Reality technology, MedVir allows the expert to interact with the data in order to tailor it to the experience and knowledge of the expert.

Arabidopsis thaliana bZIP44: a transcription factor affecting seed germination and expression of the mannanase-encoding gene AtMAN7.

Endo-β-mannanases (MAN; EC. 3.2.1.78) catalyze the cleavage of β1[RIGHTWARDS ARROW]4 bonds in mannan polymers and have been associated with the process of weakening the tissues surrounding the embryo during seed germination. In germinating Arabidopsis thaliana seeds, the most highly expressed MAN gene is AtMAN7 and its transcripts are restricted to the micropylar endosperm and to the radicle tip just before radicle emergence. Mutants with a T-DNA insertion in AtMAN7 have a slower germination than the wild type. To gain insight into the transcriptional regulation of the AtMAN7 gene, a bioinformatic search for conserved non-coding cis-elements (phylogenetic shadowing) within the Brassicaceae MAN7 gene promoters has been done, and these conserved motifs have been used as bait to look for their interacting transcription factors (TFs), using as a prey an arrayed yeast library from A. thaliana. The basic-leucine zipper TF AtbZIP44, but not the closely related AtbZIP11, has thus been identified and its transcriptional activation upon AtMAN7 has been validated at the molecular level. In the knock-out lines of AtbZIP44, not only is the expression of the AtMAN7 gene drastically reduced, but these mutants have a significantly slower germination than the wild type, being affected in the two phases of the germination process, both in the rupture of the seed coat and in the breakage of the micropylar endosperm cell walls. In the over-expression lines the opposite phenotype is observed.

CsRAV1: a candidate gene for improving lignocellulosic biomass yield

Fast-growing tree species of Populus spp.,Salix spp. and Eucalyptus spp. are cultivated to produce wood in a short time. Poplars are cultivated with cycles of 15-18 years to obtain saw timber and peeler logs, but when grown as short -rotation coppice(SRC) to produce biomass, planting density increases and rotation is considerably reduced (3-5 years). In this regard, research efforts are focused in the identification of traits and loci that allow the generation of improved SRC biomass-yielding genotypes. Biomass yield is a highly complex trait as it is the combined outcome of many other complex traits, each under separate polygenic control. Among profitable biomass yield-related traits are the amount of sylleptic branching and the length of winter dormancy. In poplar and in a few other Salicaceae species some lateral buds grow out sylleptically, the same season in which they form without the need of an intervening rest period. Sylleptic branching in poplar increases branch number, leaf area and general growth of the tree in its early years, and is a reasonable predictor of coppice yield. On the other hand, the length of winter dormancy determines the extent of the growth period. Our group has characterized the RAV1 gene of Castanea sativa (CsRAV1), encoding a transcription factor of the subfamily RAV (Related to ABI3/VP1). CsRAV1 expression shows a marked seasonal pattern, being higher in autumn and winter both in stems and buds. We generated transgenic lines of the hybrid clone Populus tremulax P. alba INRA 717 1B4 constitutively expressing CsRAV 1. These CsRAV1-expressing poplars develop sylleptic branches only a few weeks after potting. In addition to the sylleptic branching phenotype, these trees show phenological features that could give rise to an extended growth period. We are currently assessing the phenotype and behavior of these transgenic trees in a field trial, and ultimately, we will evaluate the impact on lignocellulosic biomass quality and production.