Identification and Functional Analysis of Healing Regulators in
Two major challenges in our understanding of epithelial repair and regeneration is the identification of the signals triggered after injury and the characterization of mechanisms initiated during tissue repair. From a clinical perspective, a key question that remains unanswered is “Why do some wounds fail to heal?” Considering the low genetic redundancy of Drosophila and its high degree of conservation of fundamental functions, the analysis of wound closure in imaginal discs, whose features are comparable to other post-injury events, seems to be a good model. To proceed to genomic studies, we developed a healing-permissive in vitro culture system for discs. Employing this method and microarray analysis, we aimed to identify relevant genes that are involved in healing. We compared cells that were actively involved in healing to those not involved, and identified a set of upregulated or downregulated genes. They were annotated, clustered by expression profiles, chromosomal locations, and presumptive functions. Most importantly, we functionally tested them in a healing assay. This led to the selection of a group of genes whose changes in expression level and functionality are significant for proper tissue repair. Data obtained from these analyses must facilitate the targeting of these genes in gene therapy or pharmacological studies in mammals.
Published in the journal:
. PLoS Genet 11(2): e32767. doi:10.1371/journal.pgen.1004965
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pgen.1004965
Summary
Two major challenges in our understanding of epithelial repair and regeneration is the identification of the signals triggered after injury and the characterization of mechanisms initiated during tissue repair. From a clinical perspective, a key question that remains unanswered is “Why do some wounds fail to heal?” Considering the low genetic redundancy of Drosophila and its high degree of conservation of fundamental functions, the analysis of wound closure in imaginal discs, whose features are comparable to other post-injury events, seems to be a good model. To proceed to genomic studies, we developed a healing-permissive in vitro culture system for discs. Employing this method and microarray analysis, we aimed to identify relevant genes that are involved in healing. We compared cells that were actively involved in healing to those not involved, and identified a set of upregulated or downregulated genes. They were annotated, clustered by expression profiles, chromosomal locations, and presumptive functions. Most importantly, we functionally tested them in a healing assay. This led to the selection of a group of genes whose changes in expression level and functionality are significant for proper tissue repair. Data obtained from these analyses must facilitate the targeting of these genes in gene therapy or pharmacological studies in mammals.
Introduction
Damage to an organism initiates a cascade of events that includes inflammation and the formation and remodeling of new tissue. Multiple studies have revealed significant similarities between how tissues are rebuilt during repair episodes and how they are built during development [1]. Thus, when considering epithelial repair, clear parallels exist at the structural level, as well as in signaling and the control of gene expression with the embryonic dorsal closure or the fusion of imaginal discs in Drosophila, ventral enclosure in C. elegans or eyelid closure in vertebrates [1–5]. Remarkably, co-assembly of actin cables and filopodial protrusions are instrumental in all these processes, with the majority being dependent on signaling by the JNK cascade [3, 6–9].
In invertebrates, the immediate wound healing response involves the formation of a temporary plug that encapsulate invading microbes, along with the activation of melanization and cross-linking enzymes [10]. In Drosophila larvae the outer part of the plug forms a scab within the first few hours of being wounded. The surrounding epidermal cells then orient themselves towards the wound and spread to reestablish a continuous epithelium. In larval epithelia wounds, the inactivation of the JNK pathway inhibits epidermal spreading and reepithelialization [11]. Remarkably, healing of incisional wounds in Drosophila adults also proceeds through a JNK dependent, lamellipodial-directed, epidermal cell spreading and scab formation [12]. In contrast, and similar to vertebrate early embryos, Drosophila embryos wounded by laser beams accumulate an actomyosin cable at the leading edge and display dynamic filopodial protrusions. The same is observed when imaginal discs are subjected to mechanical or genetic injuries [3, 13, 14]. The actin cable and the filopodial protrusions appear to participate in efficient sealing. The cable pulls the wound margins like a purse-string and assembles in response to different signaling events, including the activity of the JNK pathway [3, 13].
While we begin to understand some of the mechanisms of wound healing, many aspects remain unanswered. Genetically amenable systems, such as the healing of Drosophila imaginal discs, do not model all aspects of the repair process, such as inflammation or connective tissue contraction and fibrosis. However, they do allow the less understood aspects of the healing response to be explored; e.g. which cues regulate the migratory and proliferative machinery of reepithelialization? Which signal(s) and transduction pathways are responsible for initiating actin remodeling at the leading edge? How is sealing of the epithelia achieved? What inhibitory cues cancel the healing machinery upon completion? Answering these and other questions first requires the systematic identification and functional analysis of wound specific genes.
Examination of gene expression patterns using whole-genome techniques has been instrumental in the identification of many genes participating in specific developmental processes with potential relevance to disease pathogenesis [15]. To define the transcriptome program and identify pathways and genes critical to wound healing we have carried out genome-wide expression profiling and functional screens for potential healing regulators in Drosophila imaginal discs.
Imaginal discs are epithelial sacs composed by several thousand cells that, upon metamorphosis, evert and differentiate to form the external structures of the adult [16]. Imaginal discs are composed of two opposing layers of epithelial cells; a highly columnar epithelium (CE) that develops into adult structures during metamorphosis, and a layer of flat, squamous cells (peripodial epithelium—PE), which has a leading role in disc eversion and fusion [17, 18]. Shortly after puparium formation, imaginal discs evaginate and assume their positions at the surface of the prepupal animal. About 6 hours after pupariation, the proximal disc regions detach from their connection with the larval epidermis and expand, gradually replacing it. Contralateral and ipsilateral discs adjoin each other and finally meet and fuse. Spreading is driven by filopodia and actin bridges, which mediate the stretching of the imaginal epithelia [7]. This process, like embryonic dorsal closure and the healing of wounds either in embryos, larval and adult epidermis, or imaginal discs, depends on JNK activity and involves the specification of leading edge cells, which express puckered [7]. Wing imaginal discs mutant for the JNK pathway do not spread, or fail to reach the midline and seal the imaginal thorax. Similar abnormalities are observed after overexpressing Puc or dominant-negative JNK transgenes [19].
Injured imaginal discs heal and regenerate after implantation in adult abdomens, or in culture in growth permissive conditions [20, 21]. After cutting, CE and PE cells at the wound edges initiate healing through formation of a supracellular actin cable and emission of filopodial protrusions. Both epithelia expand and the wound edges zip and match to fulfill healing. Subsequently, discs engage in regenerative processes. The activation of the JNK pathway is observed in several rows of cells at the edge of discs wounds and is instrumental in the process promoting actin cable and filopodia formation [13]. These different steps are fully recapitulated in the repair of imaginal tissues after genetic ablation [14].
We have performed a comparative genome-wide profiling of imaginal cells actively involved in healing (JNK-activity positive) vs their non-engaged siblings. In these analyses, we identified a set of potential wound specific genes and defined specific co-expressed gene families participating in the healing process by Gene Ontology (GO), comparisons with published transcriptomic data, and clustering analysis. Importantly, besides identifying potentially relevant genes, we tested the function of many of their gene products by genetic interference and overexpression in two different processes; the fusion of imaginal discs during metamorphosis and wound healing in cultured discs. The fusion of imaginal discs was chosen as a functionally related proxi to discs healing. These non-saturated analyses identified a set of genes whose change in expression level and functionality was significant for proper tissue repair. Amongst those genes identified in the analysis, we explored in detail the role of individual subunits of the TCP1 chaperonin complex. Chaperones recognize proteins in non-native conformations and mediate the acquisition of the native structure [22]. One such chaperone complex is the TCP1 complex, which is known to mediate the folding of actin and tubulins [23]. We found that interference with TCP1 function affects actin cytoskeleton dynamics leading to early healing failure in imaginal discs. This indicates that the controlled expression of TCP1 subunits constitutes a new regulatory element essential to facilitate protein (actin) folding and wound healing.
Studies of the processes leading to wound closure should allow us to elucidate some of the general mechanisms that regulate natural cell movements during embryogenesis and late development. The factors identified in this study may thus be instrumental in the control of morphogenetic epithelial plasticity.
Results
To investigate Drosophila genes that are specifically involved in healing, wounded imaginal discs were cultured in vitro, and one-channel microarrays used to compare the gene expression profiles of healing-engaged cells (those showing activation of the JNK signaling cascade) with cells not participating in healing (silent JNK activity).
Healing of incised wing discs in culture
Firstly, we developed an assay to culture and image wing imaginal discs isolated from synchronized late third instar larvae in vitro. This assay employed a modified medium, which we found to be healing-permissive (see Materials and Methods). We uncovered that the healing of incised wing discs fully resembled that of wounded discs cultured in adult fly abdomens in vivo (Fig. 1A).
Specifically, healing initiation could be morphologically observed after 6 hours in culture. Both disc epithelia (CE and PE) curled towards each other, reducing the wound surface, and establishing heterotypic contacts. This contractile curling seemed to be carried out by microfilaments underlying the CE [24]. Next, the CE initiated wound ‘zippering’ by emitting filopodial extensions from both, the basal and apical surfaces. At this time, actin accumulation was observed at the edges of the wound, initially in the PE. This actin ‘cable’ persisted in healing discs as long as 18 hours. Lastly, the CE cells facing the wound vertex changed their shape elongating towards the wound edge (Fig. 1B). After 12 hours of culture the elongated PE cells initiated wound closure, leading to a remarkable reduction in the wound surface. By this time, homotypic contact had formed within the CE and PE, and the CE completed basolateral zippering (Fig. 1C). Wound closure was completed when the CE and PE met and fused after 18 hours. Afterwards, both, the CE and PE underwent tissue relaxation, leaving no scar behind (Fig. 1D and S1 Movie). All these steps precisely resemble the healing of imaginal discs in vivo [13, 24], although with a slightly slower kinetics (18 instead of 12 hours for completion). A mass of cell debris was frequently seen attached to the wound vertex. However, this was extruded upon edges apposition.
In vitro cultures allow the epithelial sealing process to occur without any influence or assistance from circulating blood cells (haemocytes) and potential inflammatory responses. It also permits the isolation of sibling healing-engaged and healing-silent cells.
Using transgenic flies expressing GFP under the indirect control (Gal4-UAS system) of a specific JNK-responsive puc enhancer (see Materials and Methods), we found that puc is expressed in the disc stalk and at low levels in PE cells in non-wounded discs cultured in vitro. We also found that the JNK signaling becomes activated in cells participating in healing. Four hours after wounding, GFP expression was initiated both, in CE and PE cells along the epithelial wound edges. By 16 hours, high levels of GFP had expanded laterally from the leading edge and could be observed in all cells actively engaged in healing at the time of sealing completion (Fig. 1E). This observed dynamics of puc expression (S1 Movie) supports the previously described correlation between healing and JNK activity.
Expression profiles of wild type and wounded wing imaginal discs
In order to quantitatively isolate puc expressing and non-expressing cells, we cultured wounded and non-wounded discs isolated from synchronized late third instar larvae for a length of time that maximized GFP expression and healing response. To isolate distinct viable cell populations expressing GFP, the discs were then subjected to dissociation and cell sorting by FACS (S1 Fig.).
The genes involved in healing were identified from whole genome expression profiles after four different populations were sampled: (1) JNK-active cells from wild type wing discs (JNK+); (2) JNK-silent cells from wild type wing discs (JNK-); (3) JNK-active cells from wounded wing discs (JNK+ W); and (4) JNK-silent cells from wounded wing discs (JNK- W). Probes prepared from these samples were hybridized to microarrays and their median expression ratios were compared. Using stringent filter settings (see Materials and Methods) different sets of transcripts were identified.
Changes in gene expression between JNK-positive and negative cells were analyzed for both, wounded discs and controls (see Materials and Methods). Firstly, we performed a global comparison in wounded discs (2-fold change, p-value < 0.05) of healing-competent (JNK+ W) cells vs. their non-engaged siblings (JNK- W). This rendered 294 upregulated and 180 downregulated genes. When more relaxed conditions were applied (1.5-fold change, p-value < 0.05), additional sets, comprising 439 upregulated and 568 downregulated genes, were identified (Fig. 2A and S1 Table).
Secondly, we analyzed statistically significant variations between the changes observed in both experimental conditions (wounded vs. controls). Genes whose changes were found to be significant (1.3-fold change, FDR adjusted p-value < 0.05) in any comparison were kept for analysis. Using this method, three subpopulations were identified: 1) W, genes whose expression is modified in wounded discs (JNK+ W vs. JNK- W); 2) NW, genes whose expression is modified in non-wounded discs (JNK+ vs. JNK-) and 3) D, genes that appear differentially modified between JNK-positive and negative cells in wounded vs. unwounded discs. Combining these subpopulations, seven gene sets could be defined according to the significance pattern displayed by the analyzed comparisons (Fig. 2B). In this way, it was possible to identify the genes differentially expressed in wounded discs only (WO—genes present in the W and D but not in the NW subpopulations) or the genes differentially expressed in both, wounded and non-wounded discs, but distinctively in each condition (W/NW/D—genes present in the W, NW and D subpopulations).
The WO subpopulation comprised 430 genes (293 upregulated and 137 downregulated) that were differentially expressed between JNK+ and JNK- cells in wounded discs only (S2 Table). On the other hand, the W/NW/D subpopulation was composed of 610 genes that were differentially expressed between JNK-positive and negative cells both, in wounded (392 upregulated and 218 downregulated) and in non-wounded (425 upregulated and 185 downregulated) discs. In this case the respective changes were significantly different between the two conditions (e.g. CG1225, upregulated in JNK+ W cells in wounded discs, but downregulated in JNK+ cells in non-wounded controls) (S3 Table).
In our analysis, we compared the absolute expression values for each gene for all experimental conditions evaluated. Thus, it was possible to determine if the overall variation in expression was due to changes in wounded discs relative to controls for individual genes; e.g. an apparent differential increase in the expression of a particular gene in JNK+ W cells in wounded discs could be a consequence of either the upregulation of the gene in this population or its specific downregulation in JNK- W siblings (relative to its level of expression in non-wounded controls) (S2 and S3 Figs.). These comparisons allowed us to categorize the genes of the WO set into 4 different subsets: (1) genes that show upregulation in JNK+ W cells; e.g. CG10540 (Fig. 2C); (2) genes that show downregulation in JNK- W cells; e.g. CG9298 (Fig. 2D); (3) genes that show downregulation in JNK+ W cells; e.g. GCR (ich) (Fig. 2E); and (4) genes that show upregulation in JNK- W cells; e.g. CG15313 (Fig. 2F) (S4 Table, see also S2 Fig.). The genes in the W/NW/D set were also categorized into different subsets. Six subsets depict the autonomous upregulation or downregulation of gene expression in JNK+ W cells in the absence of any change in JNK- W cells: (1) genes that show higher upregulation in wounded discs than in non-wounded; e.g. CG13744 (Fig. 2G); (2) genes that show higher upregulation in non-wounded discs than in wounded; e.g. CG13848 (Fig. 2H); (3) genes that shows higher downregulation in wounded discs than in non-wounded; e.g. CG31640 (Fig. 2I); (4) genes that show higher downregulation in non-wounded discs than in wounded discs; e.g. Klar (Fig. 2J); (5) genes upregulated in wounded discs but downregulated in non-wounded discs; e.g. raps (Fig. 2K); and (6) genes downregulated in wounded discs but upregulated in non-wounded discs; e.g. CG12255 (Fig. 2L). Six equivalent subsets illustrate non-autonomous changes in gene expression in JNK- W cells in the absence of any expression change in JNK+ W cells. Last, alternative combinations portraying different simultaneous autonomous and non-autonomous effects (up to 27 subsets in total) were also noted (S5 Table, see also S3 Fig.).
In order to group genes displaying the same pattern of expression changes throughout the four experimental conditions (JNK+ W, JNK- W, JNK+ and JNK-), which could be overlooked in the previous classifications, a further cluster analysis was performed. We used a template matching approach [25] to categorize absolute expression values [high (3), medium (2) and low (1)] for all genes expressed in each condition. These were used to generate different profiles [up to 81 potential combinations—from JNK+ W 1, JNK- W 1, JNK+ 1 JNK- 1 (1.1.1.1) to JNK+ W 3, JNK- W 3, JNK+ 3 JNK- 3 (3.3.3.3)]. In total, it was possible to allocate 288 genes (313 probes) to 34 different clusters (see Materials and Methods, S4 Fig. and S6 Table).
Validation of genomic expression analysis
To validate the datasets obtained from the microarray analyses we compared the level of expression of selected genes. Protein expression was analyzed for selected genes by direct immunostaining of wounded discs during healing. MMP1, a secreted matrix metalloproteinase previously identified as essential for embryonic healing and regulated in response to JNK activity [26] (upregulated 7.43 times in GFP-positive cells) showed ubiquitous staining in the wing discs that was enhanced in those cells surrounding the healing edge (Fig. 3A). An overall but not full correlation between the domain of MMP1 expression enhancement and puc-GFP expressing cells was observed. These differences may result from variations in GFP expression levels (driven by the Puc-Gal4 driver) amongst healing-engaged cells. Similar results were observed for Spaghetti squash (Sqh), a homologue of the regulatory light chain of non-muscle myosin II, whose participation in wound repair has been previously documented [27] (upregulated 1.73 times in GFP-positive cells) (Fig. 3B); Talin [28], a FERM domain protein linking the cytoplasmic tail of integrins to the actin cytoskeleton (upregulated 1.65 times in GFP-positive cells) (Fig. 3C); Paxillin [29], a cytoskeleton scaffolding protein present at focal adhesions (upregulated 2.30 times in GFP-positive cells) (Fig. 3D); Cheerio [30] coding for filamin, an actin binding protein (upregulated 3.12 times in GFP-positive cells) (Fig. 3E) and Windbeutel [31], a resident protein of the endoplasmic reticulum (upregulated 2.98 times in GFP-positive cells) (Fig. 3F). Altogether, these data supported the accuracy of the transcriptome analysis (see Discussion).
Cluster analysis
One potential cause of the co-expression of genes in different cell populations could be their common chromosomal location, which may lead to the hijacking of shared enhancers. To investigate the existence of latent healing-specific enhancers shared by co-expressed genes, we surveyed for the presence of gene clusters at specific chromosomal locations. We defined clusters as group of genes located on the same chromosome and in a close proximity to each other, but not necessarily adjacent to each other that showed the same expression behavior (upregulation or downregulation) during healing (see Materials and Methods and Fig. 4A). For the global comparison we identified 33 upregulated and 19 downregulated clusters (S7 and S8 Tables), while for the W/NW/D set we found no upregulated and 10 downregulated clusters (S9 Table). No clusters were identified for the WO set.
Ontology analysis of healing genes
The differential expression of a gene in healing-engaged cells does not guarantee their significance as a regulator or structural element participating in healing. Their relevance will only be shown after characterizing their functionality. To prioritize “healing” genes for functional assays we performed a GO term enrichment analysis [32]. For the global comparison between JNK+ W and JNK- W cells we identified 5 Cellular Component (CCs), 33 Biological Process (BPs) and 20 Molecular Function (MFs) enriched terms with a p-value < 0.01. In the WO subpopulation, we detected 17 CCs, 26 BPs and 22 MFs enriched terms with a p-value < 0.05; while 10 CCs, 78 BPs and 34 MFs were distinguished for W/NW/D (Fig. 4B and S10 to S12 Tables). These GO term enrichment analyses gave us a complex heterogeneous outline of the potentially relevant gene categories differentially expressed during healing.
Functional analysis of “healing” genes
GO term enrichment was extremely heterogeneous and to select genes for functional analysis we did not take it in consideration. Instead, we utilized three specific criteria: (1) Robust changes of gene expression, (2) availability of functional tools, and (3) relevance of gene annotation (when known) in relation to healing in Drosophila or other models and/or to JNK signaling.
To investigate the potential role of selected genes, we interfered with the expression of upregulated genes, or restore, by overexpression, those that were downregulated. A search for suitable tools led to collect 309 UAS-RNAi lines potentially able to interfere with the expression of 202 upregulated genes (S13 Table). Furthermore, we collected 21 UAS lines designed to direct the expression of 15 downregulated genes (S14 Table). These two collections are obviously non-saturating and future expansion of this analysis should provide additional relevant information on the functionality of untested genes.
In an initial round, we analyzed the functionality of RNAi and overexpression lines in imaginal discs fusion (an amenable model resembling the healing process—see Introduction). We employed a Pnr-Gal4 driver and scored notum fusion phenotypes. Pnr-Gal4 is expressed in a broad domain in the presumptive notum region of the wing disc, reflecting the expression pattern of the pannier (pnr) gene [19]. In summary, 151 RNAi lines (48.86% of those tested) corresponding to 104 genes (51.48% of those tested) displayed distinguishable phenotypes ranging from mild notum suture imperfections to early lethality. It is important to note that these figures probably underscore the real number of potential positive regulators as many of the tested RNAi lines had not been analyzed before and they could eventually fail to abolish the expression of their presumed targets. 16 overexpression constructs (76% of those tested) corresponding to 11 genes (73% of those tested) also led to thorax fusion defects. The phenotypes obtained from this screen were grouped according to their defects [bristle polarity defects (1TC), excess of fusion (2TC), weak notum clefts (3TC), medium notum cleft (4TC), strong notum clefts (5TC), pupal lethal (6TC) and embryonic lethal (7TC)] (S15 Table; See Fig. 5A to 5G for representative phenotypes of RNAi and overexpression lines).
Pupal lethality or strong thorax phenotypes were, in some cases, further analyzed by monitoring the shape of leading cells, actin accumulation or puc (JNK activity) expression. For this purpose we crossed specific UAS-RNAi lines with Pnr-Gal4, pucE69-LacZ recombinant flies. Thoraces from pupae of 5 to 7 h after puparium formation (APF) were dissected and monitored for actin (Phalloidin) and puc expression (LacZ) (See Fig. 5H to 5K). Interestingly, neither a failure to complete thorax closure, nor an unrestrained fusion, necessarily results in the suppression of puc expression. In those cases presented, leading edge cells sustain high levels of JNK activity despite alterations to cell shape, actin accumulation or ectopic cell death. Whether these genes act as effectors of the JNK signaling pathway during healing, or participate in the healing process in a JNK-independent manner, remains to be determined.
A second screen was performed for a collection of genes rendering strong defects in thorax closure. The criteria used to select these genes were their significant upregulation in the microarrays, the magnitude of their defects in thorax fusion, and considerations from GO data about their potential as healing regulators. 53 UAS RNAi lines corresponding to 44 genes and 4 UAS-overexpressing lines for 4 genes were tested (S15 Table). Each UAS line was screened using two different Gal4 lines: Pnr-Gal4 or En-Gal4. For Pnr-Gal4, the expression of several UAS transgenes led to lethality, which prevented us from performing the healing assays. En-Gal4 is a posterior compartment specific Gal4 driver that reflects the expression pattern of the engrailed (en) gene [33]. The phenotypes induced by transgenes expressed under this line are confined to the posterior region of the disc. The anterior compartment can therefore be used as internal control. All crosses were set at 25°C and third instar larval imaginal discs were dissected, wounded in a specific region depending on which Gal4 line was used (in most cases employing the En-Gal4 driver) and cultured in modified medium for 20–24 hours. At least 5 discs were screened for healing defects for each genotype after performing incisive cuts ≈ 30 μm long. The data obtained were extremely robust and almost no variation was observed between individual samples. RNA interference for 35 UAS RNAi lines corresponding to 29 genes (66,0 and 65,9% respectively of those tested) showed healing defects, while all 4 genes tested by overexpression (100%) produced a unique phenotype. These phenotypes (both, for interfered and overexpressed genes) were grouped according to the timing of healing arrest and the qualitative evaluation of defects in cell contacts and actin cytoskeleton organization (see Figs. 6 and 7 and S5 Fig.).
Class 1. Interference with the expression of genes belonging to this group (act42a-CG12051 (Fig. 6A), verprolin (vrp)-CG13503, pvf1-CG7103, reptin-CG9750, CG12007, rhea/Talin-CG6831, TCP1η-CG8351 and TCP1ζ-CG8231) completely prevented healing initiation. The wounded tissue showed no signs of filopodia formation at the leading edge. Extensive apoptotic nuclei could also be observed (See S2 Movie for TCP1ζ).
Class 2. Interference with the expression of jra-Jun related antigen-CG2275 or aats-val-CG4062 (Fig. 6B) caused the healing process to pause as early as 6 hours after wounding. The leading edge of the wounded tissue showed unstructured actin accumulation and no filopodia formation. The cells of the CE align along the edge rounding at the vertex and showed distinct shape changes. Heterotypic contacts were not achieved.
Class 3. Imaginal discs in which the expression of rho1-CG8416 and PDGF- and VEGF-receptor related-CG8222 (Fig. 6C) is abolished, completed most of the initial steps in healing (within 6 hours), including actin accumulation, heterotypic contacts between two membranes and initiation of zippering of the CE. However, they showed no vertex cells rounding, peripodial sealing or obvious filopodia at the leading edge (they are essentially defective in homotypic contact formation).
Class 4. Healing assays on discs defective in capping protein α-CG10540, capping protein β-CG17158 (Fig. 6D), CG15027 and cytC-p-CG17903 yielded incomplete closure and no vertex cells rounding, although heterotypic contacts were present and partial PE sealing was achieved. In the CE, filopodia do not form properly, a great accumulation of actin can be observed and zippering was affected.
Class 5. Imaginal discs in which the expression of CG7296, scab-CG8095 (Fig. 6E), scarface-CG11066, CG15611 or sqh-CG3595 is abolished, or the expression of mirror-CG1943 (Fig. 6F), caup-CG10605 and lamin-CG6944 is increased by overexpression, yielded incomplete healing. Here, the healing process halts after 12 hours of culture, causing both, CE and PE epithelia to proceed up until the closure step, before discontinuing. Gaps are found between the epithelia. Both, the PE and CE established homotypic contacts but they failed to complete closure.
Class 6. Downregulation of fimbrin-CG8649 (Fig. 6G) by UAS-RNAi in healing assays yielded incomplete closure after 18 hours. The PE, and to a greater extent the CE, fail to tightly join. The tissue remains loose and disorganized with frequent apical gaps.
Class 7. Healing assays for arc1-CG12505 (Fig. 6H), serpin55B-CG10913 (spn 6) or CG10200 deficient imaginal discs yielded complete healing of both, CE and PE. However, the disc tissues failed to relax back to an untensed condition, leading to numerous constricted folds. This was most evident at the CE.
The expression of TCP1 chaperonin complex subunits is a regulator of actin folding and wound healing in imaginal discs
The strong phenotypes observed following interference with the TCP1 subunits TCP1η and ζ prompted us to further explore its role in tissue healing. TCP1 as a cytosolic chaperonin that promotes the folding of non-native actins and tubulins in eukaryotic cells in vivo [34]. TCP1η and ζ showed significant upregulation, 1.64- and 2.54-fold respectively, in direct comparisons in healing-engaged cells. Their downregulation by UAS-RNAi during thorax closure resulted in strong defects in imaginal disc fusion (see above). Importantly, strong thoracic clefts were also observed after knocking down every subunit of the complex, with the exception of TCP1γ, which resulted in embryonic lethality (Fig. 8A to 8H). This suggests that the fusion of the imaginal discs is dependent on a systemic function implemented by the TCP1 complex as a whole, i.e., the loss of anyone subunit will reduce the functionality of the complex in this process.
Healing assays with En-Gal4 and RNAi lines specific for the different TCP1 subunits (Fig. 8I and J and S6 Fig.) yielded an open wound phenotype after 20–24 hours of culture. Healing was blocked as early as 6 hours after wounding. This early phenotype was observed for all subunits except TCP1γ, whose inhibition led to embryonic lethality. Both, the CE and PE epithelia failed to close wounds. The wounded tissue showed disorganized actin accumulation while filopodia formation at the leading edge was abolished (Fig. 8I and 8J). On the other hand, microtubules, which in the wild-type condition align along the direction of healing progression (Fig. 8K) reposition themselves to transverse orientations on the epithelial leading cells of interfered discs (Fig. 8L). As for imaginal fusion, all subunits are apparently required for wound closure.
To analyze the potential role of Drosophila TCP1 on the regulation of actin dynamics, we performed RNAi mediated knockdowns of TCP1 subunits in Drosophila S2R+ cells. S2R+ cells were transiently transfected with subunit specific TCP1 dsRNAs and an actin reporter construct, pMT-Act-GFP, carrying an actin-GFP fusion under the control of a metal activated promoter (see Materials and Methods). Four days after transfection, cells were subjected to live imaging, monitoring actin dynamics. TCP1 dsRNA (η subunit) treated S2R+ cells changed their morphology, became rounded, and displayed no filopodial protrusions (S7 Fig.). Further, TCP1 appears to be necessary for the production and polymerization of actin. We analyzed TCP1 ability to influence actin filaments polymerization by employing Latrunculin A (LatA) [35] (see Materials and Methods). Both, control and TCP1 dsRNA (η subunit) transfected cells became fully rounded immediately after LatA treatment. Remarkably, after a recovery period of 48 hours in the absence of the drug, control cells restore their normal morphology and F-actin levels, while dsRNA treated cells remained rounded with no signs of actin polymerization recovery (S7 Fig.). Since imaginal fusion and wound healing are actin-dependent, and the actin cytoskeleton is deregulated following the knockdown of TCP1 subunits, it is reasonable to assume that the regulation of actin synthesis and folding by TCP1 is essential in Drosophila for proper epithelial tissue healing.
Discussion
Wing discs healing transcriptome
We have developed a new in vitro model of wound healing for wing imaginal discs that recapitulates most previously described discs healing stages. We employed this method as a mean to characterize the transcriptomic make up of healing discs.
The transcriptomic profile of healing tissues has been previously explored in Drosophila. Comparison of the RNA profiles of laser wounded vs non-wounded wild type and macrophage-deficient serpent (srp) embryos [36] led to the identification of wound-activated haemocyte genes important in the immune response to wounding. Other analyses were performed in regenerating wing discs manually fragmented and implanted into adult females [37]. This last approach allowed the identification of multiple genes expressed in the early and late stages of regeneration. Still, a lack of discrimination between healing engaged and passive cells reduced its informative capacity.
Based on this current knowledge, we aimed to identify the “healing transcriptome” of those cells activating JNK signaling upon wounding in imaginal discs. In injured imaginal discs, puc expression is detected after just a few hours suggesting that the activation of the JNK pathway is one of the earliest responses to wounding. JNK activity is required for discs healing and early regeneration and cell lineage studies show that puc-expressing cells contribute to the regenerated tissue [13, 38–41]. How the JNK pathway is activated during healing is not known. Neither, the primary responses triggered by the activation of JNK signaling following injury. We employed pucE69-Gal4 I/UAS-GFP transgenic animals, which express GFP under the control of puc enhancers to identify and isolate JNK active cells from wounded and control discs. These cells respond both, to JNK signaling agonists and antagonists, as well as to incisive and ablating wounding [13]. To distinguish the endogenous transcriptional response of JNK active cells [18] from that of the cells that activate JNK in response to wounding, we performed quadruplicated transcriptomic assays. Our analyses tackled the early/intermediate events of the healing process, which predominantly involve cell rearrangements and dynamic changes in cell adhesion and cytoskeleton activity. To single out the transcriptomic healing response, we performed a one to one global comparison between cells involved in healing and neighboring cells in wounded discs and dual comparisons aimed to identify genes with differential expression in wounded discs only, WO (genes specific to wound healing and not to disc eversion and fusion) and with distinct expression changes in wounded and non-wounded discs, W/NW/D (genes important for both, healing and fusion). To identify those cellular and molecular functions augmented amongst the genes whose expression changes during healing we performed GO term enrichment analyses.
Ontology analysis of healing genes
When comparing our results to previously reported expression profiles of imaginal discs during regeneration [37], several similarities arise amongst enriched GO terms. Amongst them, the upregulation of genes involved in stress and immune response and regulation of transcription, or the downregulation of genes related to RNA processing. In contrast, the upregulation of genes involved in cytoskeleton organization, biogenesis and cell proliferation or the downregulation of helix-loop-helix transcription factors in JNK-positive cells in our analyses differs with the overrepresentation of these categories amongst downregulated and upregulated genes respectively during regeneration [37]. These differences are probably due to the fact that each study tackles a different paradigm. Our approach did not evaluate the regeneration response of whole discs [37] but focused on early/intermediate healing events.
Protein and cellular biosynthesis, mitochondrial activity, ATP synthesis, ribosome activity and chromatin remodeling GO terms are enriched amongst those genes whose expression is enhanced in JNK-positive cells in wounded discs only in the WO subpopulation. Cellular localization, secretory pathways, signaling and extracellular attachments are also overrepresented GO terms. In particular, a meaningful swift in the degradative properties of the tissue appears to be necessary to allow healing progression. We found that secreted MMP1 or Timp (Tissue inhibitor of metalloproteases) [26, 42] amongst other extracellular matrix proteins are specifically upregulated in JNK-positive cells upon wounding, while a different set of peptidases is overrepresented amongst downregulated genes.
Otherwise, we found that intracellular communication and regulation of cell motility and adhesion are GO terms enriched amongst the genes upregulated in JNK-positive cells in the W/NW/D subpopulation. In vertebrates, it has been extensively demonstrated the importance of secretory factors in the repair of wounded areas [43]. Also, in flies, multiple growth factors and cytokines, including PDGF, VEGF, TGFβ (Transforming Growth Factor β), FGF (Fibroblast Growth Factor) or Egr (Eiger) have been reported to be necessary for thorax fusion and healing [7, 43–45]. In particular, we found that pvf1 and pvf2, which code for PDGF- and VEGF-related factors are strongly expressed in JNK-positive cells and, remarkably, also expressed at high levels in JNK-negative non-engaged cells during healing just in wounded discs. The meaning of this upregulation is unclear but it could point to a systemic coordination between the cells actively involved in repair and the surrounding ones to accommodate the tissue topology to the healing process. This possibility remains to be tested. On the other hand, organ development, cell metabolism regulation and structural constituents of cuticle GO terms are enriched amongst the downregulated genes in JNK-positive cells in both, wounded and non-wounded discs suggesting that the functions of these genes, e.g. cuticle deposition, may need to be cut short both, for imaginal fusion and proper tissue repair.
Last, some GO terms are overrepresented in both subpopulations, WO and W/NW/D, e.g. cytoskeleton, GTPase activity, JAK/STAT signaling, induction and regulation of cell death and defense and stress responses are enriched amongst upregulated genes in JNK-positive cells. The upregulation of GTPases [2, 46–50], the JAK/STAT cascade [51–53], pro-apoptotic genes [54, 55] and innate immunity genes might be connected to the activation or propagation of JNK activity as a physiological response to stress. GTPases function may be also associated to their known roles as cytoskeleton regulators and linked to the cytoskeletal rearrangements observed in the initial phases of both, healing and discs fusion.
Functional analysis of “healing” genes
Different functional screens to identify genes involved in wound healing have been performed in Drosophila [56–58]. A forward genetic screen by insertional mutagenesis have led to the identification of 30 lethal mutants with defects in embryonic epithelia repair [56]. Further, in a screen of deficiencies and single mutations, several genes have been identified as regulators of the transcription of wound-responsive markers in response to puncture wounds in embryos [58]. Finally, UAS-RNAi transgenes have been overexpressed in larvae to identify genes involved in the control of epithelial migration during postembryonic wounds [57]. These screens are recent, and so their subsequent follow-up has been limited. However, when taken together, they provide solid evidences supporting the conservation of multiple aspects of the mammalian wound response in Drosophila.
To investigate the functionality of the genes identified in our transcriptomic analysis during tissue repair, we performed two reverse functional screens. First, we interfered with the natural process of thorax closure identifying at least 115 genes (approximately 53% of those tested) whose upregulation or downregulation resulted in closure defects. Considering, that many RNAi lines do not efficiently knock down gene expression, these results probably underestimate the real number of regulators in our transcriptome data. Several of those lines reproducibly yielding strong phenotypes were further tested in a disc-healing assay. We found that 33 (69%) of the 48 genes tested display wound healing defects.
The most significant set of relevant genes identified in the functional screenings includes factors involved in the activation or transduction of JNK signaling. Interference in the expression of Jra—DJun (as previously observed in embryo and larvae [56, 57]) and GTPase activity regulators such as loco or JAK/STAT elements (upd) [59, 60] display extreme defects in healing (Phenotypic Class 2). Interestingly, several potential regulators or mediators of JNK activity were also identified for the first time. An element of this class is CG1703, a member of the ABC-F subfamily of ABC proteins [61]. ABC-Fs are ATPases that regulate translation [62]. The downregulation of CG1703 resulted in a phenotype (Class 5TC) strikingly similar to that observed upon downregulation of the JNK pathway [18]. It is tempting to speculate that CG1703 may somehow regulate the activity of the JNK cascade at the level of translation.
Transcriptional activators such as AP-2, mirror, ara, or caup and repressors such as Ssb-c31a or rept constitute a second category of genes highly represented amongst those whose misregulation resulted in strong healing phenotypes. mirror, ara and caup, which code for helix-loop-helix proteins of the Iroquois complex [63, 64] are downregulated during healing. The overexpression in the notum of these genes results in embryonic (Class 7TC) for ara or mirror or pupal lethality (Class 6TC) for caup. Remarkably, the overexpression of mirror or caup in cultured discs resulted in early homotypic contacts between the PE and CE but a failure to proceed further (Phenotypic Class 5). The implication of this complex in wound healing is unexpected and its role remains to be explored.
Interference in genes involved in intercellular communication, such as cytokines and growth factors, which have been previously reported as mediators of healing in vertebrates [43–45], also resulted in discs healing defects. The inhibition of pvf1 or pvr1 resulted in notum malformations (Class 4TC and 5TC respectively) and caused incomplete healing with large wounds (Phenotypic Class 1 and 3 respectively). Pvf and Pvr are necessary in Drosophila for the apical assembly of the actin cytoskeleton during embryogenesis [65] and required for other diverse developmental processes [66–68]. Recently, Pvf1 and Pvr have also been shown to participate in larval wound healing [69].
Regulators of actin cytoskeleton dynamics including Act42a, vrp, TCP1η and ζ cpα and β and fimbrin are also enriched amongst those genes whose function is essential for healing. Downregulation of Act42a, vrp and TCP1η or ζ results in pupal lethality (Class 6TC) and yields an early open wound phenotype (Phenotypic Class 1). Vrp encodes a WH2-containing actin binding protein needed for actin polymerization during polarized growth [70]. Vertebrate Vrps regulate actin dynamics either by binding directly to actin or to proteins of the WASP family [71]. The TCP1 complex is a cytosolic chaperonin involved in the folding of non-native actins and tubulins, amongst other substrates [34]. Milder healing defects are observed upon interference in cpα and β and fimbrin expression. Cpα and Cpβ are regulators of actin filament growth [72] that cap the barbed ends of filaments preventing the addition or release of actin subunits [73]. In Drosophila, Cpα is required for the proper development of the wing blade primordium and its loss results in actin and myosin accumulation, cells extrusion and apoptosis [74]. On the other hand, loss of Cpβ affects border cell migration, reducing the length of actin protrusions. It also causes abnormal actin aggregation on nurse cell membranes [75]. Downregulation of cpα or β resulted in notum malformation (Class 4TC and 5TC respectively) and yielded incomplete healing affecting CE zippering (Phenotypic Class 4). Fimbrin is a cytoskeletal protein associated with microfilaments in microvilli, microspikes, membrane ruffles, and cell-substratum attachment sites [76]. It directs the formation of tightly bundled F-actin assemblies [77] and in Drosophila has been shown to participate in chromosome segregation during meiosis in females [78]. Downregulation of fimbrin results in notum malformation (Class 4TC) and incomplete healing after 18 hours (Phenotypic Class 6). Altogether, these data indicate that proper levels, folding and polymerization of actin appear to be essential for the earliest steps of healing, while the correct regulation of actin filaments growth or bundling seems to be essential for late healing completion.
The TCP1 complex
TCP1 is a cytosolic chaperone complex of about 900 kDa ubiquitously expressed [79]. Each of its eight subunits is coded by individual highly divergent genes [80]. Although many of its substrates have yet to be identified, it is well established that it participates in the folding of actin and tubulin monomers. Clinically, TCP1 has been linked to colorectal cancer progression and to prevention of polyglutamate-containing proteins aggregation in neurodegenerative diseases like Huntington’s [81, 82].
We found that two subunits of the complex, TCP1η and ζ, are upregulated in healing-engaged cells. We also found that interfering in their expression led to strong defects in thorax fusion, and to an early healing arrest. Importantly, although the remainder of TCP1 subunits was not upregulated upon wounding, when their expression was impaired, defects in thorax closure and wound healing were also observed. Fusion and healing arrest in TCP1 compromised discs is accompanied by a failure to organize actin microfilaments and to generate filopodia at the edges. It is also associated with defects in the organization of the microtubular network. This phenotype is very similar to that observed after interfering with Act42a or vrp suggesting that one of the roles of the TCP1 complex during healing might be to fold actin monomers and thus promote the active and rapid rebuilding of the cytoskeleton. Providing a sufficient level of folding capacity by increasing the transcriptional expression of its subunits would synergize with the upregulation of actin and Vrp upon wounding. This mechanism constitutes a new way of improving actin dynamics, working in parallel to those implemented by conventional actin regulators during healing.
It has been shown that specific TCP1 subunits are involved in the recruitment of distinct substrates to the complex [83]. Our findings suggest, however, that, at least during healing, the structural integrity and systemic activity of the complex is essential for its function: i.e. all subunits must be expressed at sufficient levels to modulate the TCP1 folding activity and cannot be compensated by others. Tissue repair is not related to any docking specificity associated to the individual components of the complex. Actin and tubulin have been identified as bona fide general substrates of the TCP1 complex in yeast and vertebrates. Indirectly, this suggests that actin (and maybe tubulin) might be the main substrates targeted by TCP1 during thorax closure and wound healing. Most probably, the upregulation of the TCP1η and ζ subunits responds to the need to reach sufficient expression levels to cope with the exceptional requests imposed by the healing process. How healing-engaged cells sense the need to increase TCP1 subunits mRNA levels and how they do implement it are open questions.
TCP1 in Drosophila appears to share its major substrates, actin and tubulin, with other organisms [79]. We have substantiated this assertion by noticing that interfering with TCP1 expression in Drosophila S2R+ cells results in strong alterations of their morphology and the organization and dynamics of their actin cytoskeleton. Moreover, cells deficient for TCP1 are unable to recover their natural levels of polymerized actin upon a pulse of Latrunculin A. Interestingly, a persistent upregulation of TCP1ζ during fibrotic healing in adult cutaneous and mucosal wounds in rabbits has been reported [84] suggesting that the critical role we uncovered for TCP1 may be conserved between invertebrates and vertebrates.
In summary, in this study, we have uncovered a new set of genes functionally relevant for the control and implementation of healing in Drosophila imaginal discs. We identified 33 genes whose inhibition or overexpression results in defects on healing or its complete absence. Another set of 60 genes show defects in thorax closure and they could also assist in the healing response. In relation to the signals triggering tissue remodeling, we recovered some known (msn, jra, pvf1, pvr, rho1 and scarface) and unknown elements (CG12007 and CG17003) potentially linked to the JNK and Pvf/Pvr cascades. Their study will be instrumental for understanding the initiation of the healing response. With respects to the process of tissue expansion, we, as expected, found multiple regulators of the actin cytoskeleton (act42a, vrp, sqh and cpα and β) and matrix attachments components (scab and rhea/Talin) (some never associated to healing before). We also discovered a novel role for the TCP1 chaperonin (TCP1η and ζ) in the context of healing. We concluded that TCP1 provides a systemic function transcriptionally regulated during healing as a modulator of cytoskeleton dynamics. Further, we uncovered a role for the Iroquois complex of transcriptional regulators (mirror, ara and capu) as a negative factor preventing healing progression. In addition, a series of genes with established (l(1)1Bi, AP-2, Ssb-c31a, rept, aats-val, cytC-p, broad and lamin) or unidentified (CG7296, CG15027 and CG15611) functions were documented as active elements on healing progression for the first time. Their functions in this context remain to be explored. Last, regarding the sealing of the epithelia, we identified five genes fimbrin, aats-trp, arc1, serpin55B (spn 6) and CG10200 whose absence, without affecting healing itself, leads to a failure of final sealing or impedes the tissue relaxation. How do they work? We don’t know, but they constitute genetic entry points that could help us to understand how epithelial closure is achieved or regulated.
Considering the low genetic redundancy of Drosophila, and the high degree of conservation of fundamental cellular functions and signaling pathways when compared to vertebrates, data obtained from our analyses must facilitate investigating and targeting functionally relevant genes in vertebrate healing studies. A variety of different microarray analyses of murine and humans have already described gene expression profiles of cells located at wound sites [85, 86], however, very few studies have addressed the functional role of genes potentially involved in healing regulation. There is promise that our newly identified wound-healing genes will guide further work in the more complex process of mammalian tissue repair.
Materials and Methods
Culture in vitro, wounding, live imaging and immunostaining of imaginal discs
Third instar larvae wing imaginal discs were dissected and incubated in a MM3 medium (Shields and Sang´s) supplied with fly extract. The fly extract was prepared following the protocol suggested by the DGRC, Indiana University [87]. This modified media allows wing imaginal discs to be cultured for more than 24 hours.
Third instar larvae imaginal disc dissections were performed in cultured medium and incisions of precise sizes were carefully executed with a pair of tungsten needles. This was followed by static culture in supplemented MM3 media for 16–18 hours. We found that imaginal discs in culture lose their morphology after 5 to 6 hours, becoming spherical in shape. To flatten the samples we mounted them within a silicon square “sandwich”. This enabled us to perform time-lapse confocal imaging. At least 10 imaginal discs were placed into the silicone square area. When all discs were cut the corresponding slide field was refilled with media mixed with 1 ml of FM4–64 (9 μM) to label the imaginal disc membranes. The imaginal discs were covered with a coverslip and slightly compressed to minimize the free space to avoid imaginal discs spherical deformation (eversion). The ready-made chamber also prevents desiccation. Time-lapse recording was initiated after 5 hours of incubation using inverted Leica SP5 and SP2 or Zeiss LSM700 confocal microscopes and 63 X objectives. Images (Z-stacks of 1 μm thickness) were captured every 10 minutes and in different positions to film at least 3 imaginal discs per slide. Laser intensity was kept at a minimum to avoid photobleaching and to minimize phototoxicity. Image analysis was performed with Leica and Zeiss Confocal Software and ImageJ (NIH Image) was used for mounting time-lapse movies in AVI format.
After 16–18 hours the imaginal discs were recovered using a 200 μl pipette humidified in MM3 media to prevent imaginal discs adhering to the walls of the tips. This was followed by a fast wash in PBS 1X to eliminate MM3 media remnants and fixation in 4% paraformaldehyde for 20 minutes. Antibody incubations were performed following standard procedures. All steps were carried out at room temperature on a shaker. After several washes, imaginal discs were mounted in Vectashield (Vector).
Statistical analysis
The Microarrays datasets were employed for two different types of assessments: a global comparison of healing-competent (JNK-positive) cells to their non-engaged siblings and dual comparisons between JNK-positive and negative cells of both, wounded and non-wounded discs.
For the global comparison, microarray raw intensities were converted to gene expression estimates using a robust multichip average procedure (RMA). Before the statistical analysis, a pre-filtering step was performed to eliminate the genes presenting low signal (genes should have intensity values larger than 40 units, in at least 25% of samples) and/or minimal variability across samples (inter-quartile range should be larger than the 10% percentile). 6722 probe sets, out of 18952 passed this non-specific prefiltering process. The LIMMA package from Bioconductor was used to fit linear models to log-transform expression data. A moderated t-statistic and log-odds of differential expression were computed to assess for statistically significant changes. The corresponding p-values were corrected for multiple hypotheses testing using a False Discovery Rate criterion. Once the statistical analysis was accomplished we selected genes displaying statistically supported evidence of differential change by applying a p-value and a Fold Change cut off. The results obtained with the global analysis just display expression level changes of JNK-positive vs negative cells in wounded discs.
We also performed dual comparisons employing the microarrays generated from wounded and unwounded discs. The numerical raw data for these comparisons were processed through a between-array quantile normalization step, and a final robust summarization of probe-set signals. The statistical significance of differential gene expression signals was assessed using a moderated t-statistic with the LIMMA Bioconductor package. A relatively mild fold change cut-off (1.3) was applied to facilitate the computation of the interaction terms amongst datasets. A Venn diagram (Fig. 2B) was used to represent the interactions of the different identified subpopulations (W, NW and D). Dual comparisons between wounded and unwounded discs not only gave more information about the behavior of the cells involved in healing respect to their neighboring cells, but allowed us to create graphical representations of the expression of JNK-positive and negative cells for all chips in both conditions, wounded and unwounded.
Template matching analysis
In order to perform this analysis, and disregard noisy genes, we picked 5106 genes (out of the 6722 from the pre-filtered set) that displayed an FDR-adjusted significance p-value < 0.05, for at least one of the three following contrasts: changes in expression between JNK-positive and negative cell samples in wounded (W), in non-wounded (NW) conditions, and when significant differences could be established between the both aforementioned expression changes (D).
Gene expression values were then discretized into a three level scale [high (3), medium (2) and low (1)] for all expressed genes in each condition. Different template profiles were generated (1.1.1.1 to 3.3.3.3) and clusters were defined from gene profiles that could be associated (Pearson correlation R2 > 0.95) to respective template profiles.
Gene mapping in the genome
All genes identified through our two analyses (global and dual comparisons) with a fold change higher or lower than 2 and -2 respectively were mapped to the fly genome using the RefSeq track of the UCSC genome browser. Only one transcript per gene was considered (the first one present in the annotations). The file refGene.txt was used to retrieve the coordinates of each gene in RefSeq. The statistics of gene distribution on each chromosome were computed using the 19.670 unique RefSeq genes contained in the file refGene.txt. To measure the statistical significance of the gene distribution, 1.000 random sets of genes with the same size were generated and evaluated by p-value using a hyper-geometric distribution. The sets of up- and downregulated genes in the microarray were then mapped in the genomes using appropriate catalogues.
Localized cluster identification and characterization
A cluster was defined as a group of neighboring genes located in a limited region of the genome, but not necessarily consecutive, that shows the same expression pattern (upregulation or downregulation) in the microarray experiment. Clusters were determined by the physical position in the genome, the number of co-expressed genes and the total number of genes in the genomic fragment. To evaluate the significance of the number of clusters identified in our microarray, 100 random sets with the same size and chromosomal gene distribution were generated from the sets of genes involved in wound healing. The clustering analysis was then performed in such data sets to infer the number of clusters that could be obtained by chance (expected clusters).
Different values were computed (using refGene.txt): the average gene length, the gene length distribution, the number of bidirectional/opposed genes on each cluster according to their strand, and the average intergenic distance. Such values were compared to the whole population of genes annotated in the fruit fly genome, to evaluate their significance.
Functional annotations
A hypergeometric test (GOStats) was performed to assess GO term overrepresentation relative to aleatory assignment in our selected categories (WO and W/NW/D). GO terms and relationships for Drosophila were downloaded (gene_ontology_edit.obo, release 1.2 and gene_association.fb, release 1.94) from the GO web site. Annotation of the different gene sets was done at level 3 of the 'molecular function' ontology tree. We climbed up in the ontology to obtain the corresponding parent term at level 3 for those genes that were annotated at deeper levels in the hierarchical tree. A statistical analysis was performed to test the probability of observing a given GO term significantly enriched in genes belonging to such clusters using a hypergeometric distribution.
Thorax fusion functional screen
The thorax fusion screen was carried out using a collection of transgenic fly lines carrying UAS-RNAi and UAS overexpressing constructs (see above). To drive their expression to the notum, males from each line were crossed with Pnr-Gal4 virgin females and tested at 3 different temperatures 18°C, 25°C and 29°C. Each phenotype was scored in a semi-quantitative manner based on the affected area within the Pnr-Gal4 expression domain. These phenotypes have been further classified based on the severity of phenotypes as described in the text.
Wound healing functional screen
Males carrying UAS-RNAi or overexpressing constructs were crossed with either En-Gal4 or Pnr-Gal4 females and maintained at 25°C. Third instar larval imaginal discs were wounded, cultured at 25°C and evaluated for healing defects after 18 hours of culture.
Further Materials and Methods are described in S1 Text.
Supporting Information
Zdroje
1. Martin P, Parkhurst SM (2004) Parallels between tissue repair and embryo morphogenesis. Development 131: 3021–3034. doi: 10.1242/dev.01253 15197160
2. Harden N (2002) Signaling pathways directing the movement and fusion of epithelial sheets: lessons from dorsal closure in Drosophila. Differentiation 70: 181–203. doi: 10.1046/j.1432-0436.2002.700408.x 12147138
3. Wood W, Jacinto A, Grose R, Woolner S, Gale J, et al. (2002) Wound healing recapitulates morphogenesis in Drosophila embryos. Nature Cell Biology 4: 907–912. doi: 10.1038/ncb875 12402048
4. Martin P (1997) Wound healing--aiming for perfect skin regeneration. Science 276: 75–81. doi: 10.1126/science.276.5309.75 9082989
5. McCluskey J, Martin P (1995) Analysis of the tissue movements of embryonic wound healing--DiI studies in the limb bud stage mouse embryo. Developmental Biology 170: 102–114. doi: 10.1006/dbio.1995.1199 7601301
6. Martin P, Lewis J (1992) Actin cables and epidermal movement in embryonic wound healing. Nature 360: 179–183. doi: 10.1038/360179a0 1436096
7. Martin-Blanco E, Pastor-Pareja JC, Garcia-Bellido A (2000) JNK and decapentaplegic signaling control adhesiveness and cytoskeleton dynamics during thorax closure in Drosophila. Proc Natl Acad Sci U S A 97: 7888–7893. doi: 10.1073/pnas.97.14.7888 10884420
8. Li G, Gustafson-Brown C, Hanks SK, Nason K, Arbeit JM, et al. (2003) c-Jun is essential for organization of the epidermal leading edge. Developmental Cell 4: 865–877. doi: 10.1016/S1534-5807(03)00159-X 12791271
9. Zenz R, Scheuch H, Martin P, Frank C, Eferl R, et al. (2003) c-Jun regulates eyelid closure and skin tumor development through EGFR signaling. Developmental Cell 4: 879–889. doi: 10.1016/S1534-5807(03)00161-8 12791272
10. Mace KA, Pearson JC, McGinnis W (2005) An epidermal barrier wound repair pathway in Drosophila is mediated by grainy head. Science 308: 381–385. doi: 10.1126/science.1107573 15831751
11. Galko MJ, Krasnow MA (2004) Cellular and genetic analysis of wound healing in Drosophila larvae. PLoS biology 2: E239. doi: 10.1371/journal.pbio.0020239 15269788
12. Ramet M, Lanot R, Zachary D, Manfruelli P (2002) JNK signaling pathway is required for efficient wound healing in Drosophila. Developmental Biology 241: 145–156. doi: 10.1006/dbio.2001.0502 11784101
13. Bosch M, Serras F, Martin-Blanco E, Baguna J (2005) JNK signaling pathway required for wound healing in regenerating Drosophila wing imaginal discs. Developmental Biology 280: 73–86. doi: 10.1016/j.ydbio.2005.01.002 15766749
14. Bergantinos C, Corominas M, Serras F (2010) Cell death-induced regeneration in wing imaginal discs requires JNK signalling. Development 137: 1169–1179. doi: 10.1242/dev.045559 20215351
15. Nishimura M, Kumsta C, Kaushik G, Diop SB, Ding Y, et al. (2014) A dual role for integrin-linked kinase and beta1-integrin in modulating cardiac aging. Aging Cell 13: 431–440. doi: 10.1111/acel.12193 24400780
16. Patel U, Myat MM (2013) Receptor guanylyl cyclase Gyc76C is required for invagination, collective migration and lumen shape in the Drosophila embryonic salivary gland. Biol Open 2: 711–717. doi: 10.1242/bio.20134887 23862019
17. Agnes F, Suzanne M, Noselli S (1999) The Drosophila JNK pathway controls the morphogenesis of imaginal discs during metamorphosis. Development 126: 5453–5462. 10556069
18. Pastor-Pareja JC, Grawe F, Martin-Blanco E, Garcia-Bellido A (2004) Invasive cell behavior during Drosophila imaginal disc eversion is mediated by the JNK signaling cascade. Developmental Cell 7: 387–399. doi: 10.1016/j.devcel.2004.07.022 15363413
19. Calleja M, Moreno E, Pelaz S, Morata G (1996) Visualization of gene expression in living adult Drosophila. Science 274: 252–255. doi: 10.1126/science.274.5285.252 8824191
20. Bryant PJ, Fraser SE (1988) Wound healing, cell communication, and DNA synthesis during imaginal disc regeneration in Drosophila. Developmental Biology 127: 197–208. doi: 10.1016/0012-1606(88)90201-1 2452103
21. Schubiger G (1971) Regeneration, duplication and transdetermination in fragments of the leg disc of Drosophila melanogaster. Developmental Biology 26: 277–295. doi: 10.1016/0012-1606(71)90127-8 5003476
22. Frydman J, Nimmesgern E, Ohtsuka K, Hartl FU (1994) Folding of nascent polypeptide chains in a high molecular mass assembly with molecular chaperones. Nature 370: 111–117. doi: 10.1038/370111a0 8022479
23. Brackley KI, Grantham J (2009) Activities of the chaperonin containing TCP-1 (CCT): implications for cell cycle progression and cytoskeletal organisation. Cell Stress Chaperones 14: 23–31. doi: 10.1007/s12192-008-0057-x 18595008
24. Reinhardt CA, Hodgkin NM, Bryant PJ (1977) Wound healing in the imaginal discs of Drosophila. I. Scanning electron microscopy of normal and healing wing discs. Developmental Biology 60: 238–257.
25. Pavlidis P, Noble WS (2001) Analysis of strain and regional variation in gene expression in mouse brain. Genome biology 2: RESEARCH0042. doi: 10.1186/gb-2001-2-10-research0042 11597334
26. Stevens LJ, Page-McCaw A (2012) A secreted MMP is required for reepithelialization during wound healing. Molecular Biology of the Cell 23: 1068–1079. doi: 10.1091/mbc.E11-09-0745 22262460
27. Abreu-Blanco MT, Verboon JM, Parkhurst SM (2011) Cell wound repair in Drosophila occurs through three distinct phases of membrane and cytoskeletal remodeling. Journal of Cell Biology 193: 455–464. doi: 10.1083/jcb.201011018 21518790
28. Brown NH, Gregory SL, Rickoll WL, Fessler LI, Prout M, et al. (2002) Talin is essential for integrin function in Drosophila. Dev Cell 3: 569–579. doi: 10.1016/S1534-5807(02)00290-3 12408808
29. Yagi R, Ishimaru S, Yano H, Gaul U, Hanafusa H, et al. (2001) A novel muscle LIM-only protein is generated from the paxillin gene locus in Drosophila. EMBO Rep 2: 814–820. doi: 10.1093/embo-reports/kve178 11520860
30. Sokol NS, Cooley L (1999) Drosophila filamin encoded by the cheerio locus is a component of ovarian ring canals. Curr Biol 9: 1221–1230. doi: 10.1016/S0960-9822(99)80502-8 10556087
31. Konsolaki M, Schupbach T (1998) windbeutel, a gene required for dorsoventral patterning in Drosophila, encodes a protein that has homologies to vertebrate proteins of the endoplasmic reticulum. Genes Dev 12: 120–131. doi: 10.1101/gad.12.1.120 9420336
32. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, et al. (2000) Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nature genetics 25: 25–29. doi: 10.1038/75556
33. Aza-Blanc P, Ramirez-Weber FA, Laget MP, Schwartz C, Kornberg TB (1997) Proteolysis that is inhibited by hedgehog targets Cubitus interruptus protein to the nucleus and converts it to a repressor. Cell 89: 1043–1053. doi: 10.1016/S0092-8674(00)80292-5 9215627
34. Grantham J, Brackley KI, Willison KR (2006) Substantial CCT activity is required for cell cycle progression and cytoskeletal organization in mammalian cells. Exp Cell Res 312: 2309–2324. doi: 10.1016/j.yexcr.2006.03.028 16765944
35. Coue M, Brenner SL, Spector I, Korn ED (1987) Inhibition of actin polymerization by latrunculin A. FEBS Lett 213: 316–318. doi: 10.1016/0014-5793(87)81513-2 3556584
36. Stramer B, Winfield M, Shaw T, Millard TH, Woolner S, et al. (2008) Gene induction following wounding of wild-type versus macrophage-deficient Drosophila embryos. EMBO Reports 9: 465–471. doi: 10.1038/embor.2008.34 18344972
37. Blanco E, Ruiz-Romero M, Beltran S, Bosch M, Punset A, et al. (2010) Gene expression following induction of regeneration in Drosophila wing imaginal discs. Expression profile of regenerating wing discs. BMC Dev Biol 10: 94.
38. Bosch M, Baguna J, Serras F (2008) Origin and proliferation of blastema cells during regeneration of Drosophila wing imaginal discs. Int J Dev Biol 52: 1043–1050. doi: 10.1387/ijdb.082608mb 18956337
39. Lee N, Maurange C, Ringrose L, Paro R (2005) Suppression of Polycomb group proteins by JNK signalling induces transdetermination in Drosophila imaginal discs. Nature 438: 234–237. doi: 10.1038/nature04120 16281037
40. Mattila J, Omelyanchuk L, Kyttala S, Turunen H, Nokkala S (2005) Role of Jun N-terminal Kinase (JNK) signaling in the wound healing and regeneration of a Drosophila melanogaster wing imaginal disc. Int J Dev Biol 49: 391–399. doi: 10.1387/ijdb.052006jm 15968584
41. Bergantinos C, Vilana X, Corominas M, Serras F (2010) Imaginal discs: Renaissance of a model for regenerative biology. Bioessays 32: 207–217. doi: 10.1002/bies.200900105 20127699
42. McCawley LJ, Matrisian LM (2001) Matrix metalloproteinases: they’re not just for matrix anymore! Curr Opin Cell Biol 13: 534–540. doi: 10.1016/S0955-0674(00)00248-9 11544020
43. Werner S, Grose R (2003) Regulation of wound healing by growth factors and cytokines. Physiol Rev 83: 835–870. 12843410
44. Singer A, Clark R (1999) Cutaneous Wound Healing. New England Journal of Medicine 341: 738–746. doi: 10.1056/NEJM199909023411006 10471461
45. Braddock M (2001) The transcription factor Egr-1: a potential drug in wound healing and tissue repair. Annals of Medicine 33: 313–318. doi: 10.3109/07853890109002083 11491188
46. Brock J, Midwinter K, Lewis J, Martin P (1996) Healing of incisional wounds in the embryonic chick wing bud: characterization of the actin purse-string and demonstration of a requirement for Rho activation. Journal of Cell Biology 135: 1097–1107. doi: 10.1083/jcb.135.4.1097 8922389
47. Su YC, Treisman JE, Skolnik EY (1998) The Drosophila Ste20-related kinase misshapen is required for embryonic dorsal closure and acts through a JNK MAPK module on an evolutionarily conserved signaling pathway. Genes & Development 12: 2371–2380. doi: 10.1101/gad.12.15.2371
48. Stronach B, Perrimon N (2002) Activation of the JNK pathway during dorsal closure in Drosophila requires the mixed lineage kinase, slipper. Genes & Development 16: 377–387. doi: 10.1101/gad.953002
49. Yin J, Lu J, Yu FS (2008) Role of small GTPase Rho in regulating corneal epithelial wound healing. Invest Ophthalmol Vis Sci 49: 900–909. doi: 10.1167/iovs.07-1122 18326710
50. Baek SH, Kwon YC, Lee H, Choe KM (2010) Rho-family small GTPases are required for cell polarization and directional sensing in Drosophila wound healing. Biochem Biophys Res Commun 394: 488–492. doi: 10.1016/j.bbrc.2010.02.124 20184864
51. Zeidler MP, Bach EA, Perrimon N (2000) The roles of the Drosophila JAK/STAT pathway. Oncogene 19: 2598–2606. doi: 10.1038/sj.onc.1203482 10851058
52. Pastor-Pareja JC, Wu M, Xu T (2008) An innate immune response of blood cells to tumors and tissue damage in Drosophila. Dis Model Mech 1: 144–154. doi: 10.1242/dmm.000950 19048077
53. Wu M, Pastor-Pareja JC, Xu T (2010) Interaction between Ras(V12) and scribbled clones induces tumour growth and invasion. Nature 463: 545–548. doi: 10.1038/nature08702 20072127
54. Adachi-Yamada T, Fujimura-Kamada K, Nishida Y, Matsumoto K (1999) Distortion of proximodistal information causes JNK-dependent apoptosis in Drosophila wing. Nature 400: 166–169. doi: 10.1038/22112 10408443
55. Ryoo HD, Gorenc T, Steller H (2004) Apoptotic cells can induce compensatory cell proliferation through the JNK and the Wingless signaling pathways. Developmental Cell 7: 491–501. doi: 10.1016/j.devcel.2004.08.019 15469838
56. Campos I, Geiger J.A., Santos A.C., Carlos V., and Jacinto A. (2010) Genetic Screen in Drosophila melanogaster Uncovers a Novel Set of Genes Required for Embryonic Epithelial Repair. Genetics 184: 129–140. doi: 10.1534/genetics.109.110288 19884309
57. Lesch C, Jo J., Wu Y., Fish G.S., and Galko M.J. (2010) A Targeted UAS-RNAi Screen in Drosophila Larvae Identifies Wound Closure Genes Regulating Distinct Cellular Processes. Genetics 186: 943–957. doi: 10.1534/genetics.110.121822 20813879
58. Juarez MT, Patterson RA, Sandoval-Guillen E, McGinnis W (2011) Duox, Flotillin-2, and Src42A are required to activate or delimit the spread of the transcriptional response to epidermal wounds in Drosophila. PLoS Genetics 7: e1002424. doi: 10.1371/journal.pgen.1002424 22242003
59. Belacortu Y, and Paricio N. (2011) Drosophila as a Model of Wound Healing and Tissue Regeneration in Vertebrates. Developmental Dynamics 240: 2379–2404. doi: 10.1002/dvdy.22753 21953647
60. Schäfer M, and Werner S. (2007) Transcriptional Control of Wound Repair. Annu Rev Cell Dev Biol 23: 69–92. doi: 10.1146/annurev.cellbio.23.090506.123609 17474876
61. Dean M, Rzhetsky A, Allikmets R (2001) The human ATP-binding cassette (ABC) transporter superfamily. Genome Research 11: 1156–1166. doi: 10.1101/gr.GR-1649R 11435397
62. Andersen CB, Becker T, Blau M, Anand M, Halic M, et al. (2006) Structure of eEF3 and the mechanism of transfer RNA release from the E-site. Nature 443: 663–668. doi: 10.1038/nature05126 16929303
63. Gomez-Skarmeta JL, Diez del Corral R, de la Calle-Mustienes E, Ferre-Marco D, Modolell J (1996) Araucan and caupolican, two members of the novel iroquois complex, encode homeoproteins that control proneural and vein-forming genes. Cell 85: 95–105. doi: 10.1016/S0092-8674(00)81085-5 8620542
64. Kankel MW, Hurlbut GD, Upadhyay G, Yajnik V, Yedvobnick B, et al. (2007) Investigating the genetic circuitry of mastermind in Drosophila, a notch signal effector. Genetics 177: 2493–2505. doi: 10.1534/genetics.107.080994 18073442
65. Rosin D, Schejter E, Volk T, Shilo BZ (2004) Apical accumulation of the Drosophila PDGF/VEGF receptor ligands provides a mechanism for triggering localized actin polymerization. Development 131: 1939–1948. doi: 10.1242/dev.01101 15056618
66. Cho NK, Keyes L, Johnson E, Heller J, Ryner L, et al. (2002) Developmental control of blood cell migration by the Drosophila VEGF pathway. Cell 108: 865–876. doi: 10.1016/S0092-8674(02)00676-1 11955438
67. Ishimaru S, Ueda R, Hinohara Y, Ohtani M, Hanafusa H (2004) PVR plays a critical role via JNK activation in thorax closure during Drosophila metamorphosis. EMBO Journal 23: 3984–3994. doi: 10.1038/sj.emboj.7600417 15457211
68. McDonald JA, Pinheiro EM, Montell DJ (2003) PVF1, a PDGF/VEGF homolog, is sufficient to guide border cells and interacts genetically with Taiman. Development 130: 3469–3478. doi: 10.1242/dev.00574 12810594
69. Wu Y, Brock AR, Wang Y, Fujitani K, Ueda R, et al. (2009) A blood-borne PDGF/VEGF-like ligand initiates wound-induced epidermal cell migration in Drosophila larvae. Current Biology 19: 1473–1477. doi: 10.1016/j.cub.2009.07.019 19646875
70. Vaduva G, Martin NC, Hopper AK (1997) Actin-binding verprolin is a polarity development protein required for the morphogenesis and function of the yeast actin cytoskeleton. J Cell Biol 139: 1821–1833. doi: 10.1083/jcb.139.7.1821 9412475
71. Jiang WG, Ye L, Patel G, Harding KG (2010) Expression of WAVEs, the WASP (Wiskott-Aldrich syndrome protein) family of verprolin homologous proteins in human wound tissues and the biological influence on human keratinocytes. Wound Repair Regen 18: 594–604. doi: 10.1111/j.1524-475X.2010.00630.x 20946142
72. Schafer DA, Cooper JA (1995) Control of actin assembly at filament ends. Annu Rev Cell Dev Biol 11: 497–518. doi: 10.1146/annurev.cb.11.110195.002433 8689567
73. Machesky LM, Cooper JA (1999) Cell motility. Bare bones of the cytoskeleton. Nature 401: 542–543.
74. Janody F, Treisman JE (2006) Actin capping protein alpha maintains vestigial-expressing cells within the Drosophila wing disc epithelium. Development 133: 3349–3357. doi: 10.1242/dev.02511 16887822
75. Ogienko AA, Karagodin DA, Lashina VV, Baiborodin SI, Omelina ES, et al. (2013) Capping protein beta is required for actin cytoskeleton organisation and cell migration during Drosophila oogenesis. Cell Biol Int 37: 149–159. doi: 10.1002/cbin.10025 23339103
76. Bretscher A (1981) Fimbrin is a cytoskeletal protein that crosslinks F-actin in vitro. Proc Natl Acad Sci U S A 78: 6849–6853. doi: 10.1073/pnas.78.11.6849 6947259
77. Wu JQ, Bahler J, Pringle JR (2001) Roles of a fimbrin and an alpha-actinin-like protein in fission yeast cell polarization and cytokinesis. Mol Biol Cell 12: 1061–1077. doi: 10.1091/mbc.12.4.1061 11294907
78. Dobie KW, Kennedy CD, Velasco VM, McGrath TL, Weko J, et al. (2001) Identification of chromosome inheritance modifiers in Drosophila melanogaster. Genetics 157: 1623–1637. 11290718
79. Sternlicht H, Farr GW, Sternlicht ML, Driscoll JK, Willison K, et al. (1993) The t-complex polypeptide 1 complex is a chaperonin for tubulin and actin in vivo. Proc Natl Acad Sci U S A 90: 9422–9426. doi: 10.1073/pnas.90.20.9422 8105476
80. Kubota H, Hynes G, Willison K (1995) The chaperonin containing t-complex polypeptide 1 (TCP-1). Multisubunit machinery assisting in protein folding and assembly in the eukaryotic cytosol. Eur J Biochem 230: 3–16.
81. Coghlin C, Carpenter B, Dundas SR, Lawrie LC, Telfer C, et al. (2006) Characterization and over-expression of chaperonin t-complex proteins in colorectal cancer. J Pathol 210: 351–357. doi: 10.1002/path.2056 16981251
82. Kitamura A, Kubota H, Pack CG, Matsumoto G, Hirayama S, et al. (2006) Cytosolic chaperonin prevents polyglutamine toxicity with altering the aggregation state. Nat Cell Biol 8: 1163–1170. doi: 10.1038/ncb1478 16980958
83. Spiess C, Miller EJ, McClellan AJ, Frydman J (2006) Identification of the TRiC/CCT substrate binding sites uncovers the function of subunit diversity in eukaryotic chaperonins. Mol Cell 24: 25–37. doi: 10.1016/j.molcel.2006.09.003 17018290
84. Satish L, Johnson S, Wang JH, Post JC, Ehrlich GD, et al. (2010) Chaperonin containing T-complex polypeptide subunit eta (CCT-eta) is a specific regulator of fibroblast motility and contractility. PLoS One 5: e10063. doi: 10.1371/journal.pone.0010063 20442790
85. Cole J, Tsou R., Wallace K., Gibran N., and Isik F. (2001) Early gene expression profile of human skin to injury using high-density cDNA microarrays. Wound Repair & Regeneration 9: 360–370. doi: 10.1046/j.1524-475x.2001.00360.x
86. Cooper L, Johnson C., Burslem F., and Martin P. (2005) Wound healing and inflammation genes revealed by array analysis of “macrophageless” PU.1 null mice. Genome Biol 6: R5. doi: 10.1186/gb-2004-6-1-r5 15642097
87. Lin G, Zhang X, Ren J, Pang Z, Wang C, et al. (2013) Integrin signaling is required for maintenance and proliferation of intestinal stem cells in Drosophila. Dev Biol 377: 177–187. doi: 10.1016/j.ydbio.2013.01.032 23410794
Štítky
Genetika Reprodukční medicínaČlánek vyšel v časopise
PLOS Genetics
2015 Číslo 2
- Primární hyperoxalurie – aktuální možnosti diagnostiky a léčby
- Srdeční frekvence embrya může být faktorem užitečným v předpovídání výsledku IVF
- Akutní intermitentní porfyrie
- Vztah užívání alkoholu a mužské fertility
- Šanci na úspěšný průběh těhotenství snižují nevhodné hladiny progesteronu vznikající při umělém oplodnění
Nejčtenější v tomto čísle
- Genomic Selection and Association Mapping in Rice (): Effect of Trait Genetic Architecture, Training Population Composition, Marker Number and Statistical Model on Accuracy of Rice Genomic Selection in Elite, Tropical Rice Breeding Lines
- Discovery of Transcription Factors and Regulatory Regions Driving Tumor Development by ATAC-seq and FAIRE-seq Open Chromatin Profiling
- Evolutionary Signatures amongst Disease Genes Permit Novel Methods for Gene Prioritization and Construction of Informative Gene-Based Networks
- Proteotoxic Stress Induces Phosphorylation of p62/SQSTM1 by ULK1 to Regulate Selective Autophagic Clearance of Protein Aggregates