Using conditional independence tests to elucidate causal links in cell cycle regulation in Escherichia coli.

How cells regulate their cell cycles is a central question for cell biology. Models of cell size homeostasis have been proposed for bacteria, archaea, yeast, plant, and mammalian cells. New experiments bring forth high volumes of data suitable for testing existing models of cell size regulation and proposing new mechanisms. In this paper, we use conditional independence tests in conjunction with data of cell size at key cell cycle events (birth, initiation of DNA replication, and constriction) in the model bacterium Escherichia coli to select between the competing cell cycle models. We find that in all growth conditions that we study, the division event is controlled by the onset of constriction at midcell. In slow growth, we corroborate a model where replication-related processes control the onset of constriction at midcell. In faster growth, we find that the onset of constriction is affected by additional cues beyond DNA replication. Finally, we also find evidence for the presence of additional cues triggering initiations of DNA replication apart from the conventional notion where the mother cells solely determine the initiation event in the daughter cells via an adder per origin model. The use of conditional independence tests is a different approach in the context of understanding cell cycle regulation and it can be used in future studies to further explore the causal links between cell events.

Cell Cycle-Dependent Recruitment of FtsN to the Divisome in Escherichia coli.

Cell division in Escherichia coli starts with the formation of an FtsZ protofilament network at midcell, the Z ring. However, only after a considerable lag period does the cell start to form a midcell constriction. The onset of constriction depends upon the arrival of so-called late divisome proteins, among which, FtsN is the last essential one. The timing and dependency of FtsN arrival to the divisome, along with genetic evidence, suggests it triggers cell division. In this study, we used high-throughput fluorescence microscopy to determine the arrival of FtsN and the early divisome protein ZapA to midcell at a single-cell level during the cell cycle. Our data show while the recruitment of ZapA/FtsZ is gradual in the cell cycle, recruitment of FtsN is rapid and begins at about the onset of constriction. At this time, the fraction of ZapA/FtsZ in the Z ring approaches its peak value. We also find a second increase in FtsN recruitment to the divisome, which begins once the amount of ZapA/FtsZ at midcell starts decreasing. Increasing hypermorphic FtsA* (FtsA R286W), but not FtsA, accelerates FtsN recruitment but not constriction. This finding is consistent with FtsA* recruiting FtsN with some other divisome component being rate-limiting for constriction under these conditions. Finally, our data support the recently proposed idea that ZapA/FtsZ and FtsN are part of physically separate complexes in midcell throughout the whole septation process. IMPORTANCE Cell division in most bacteria starts with the formation of an FtsZ protofilament network at midcell, the Z ring. However, cells only start to constrict after a considerable lag. A factor thought to trigger the onset of constriction in Escherichia coli is FtsN, which is the last essential protein to be recruited to the Z ring. Using a high-throughput quantitative fluorescence microscopy, we determine the cell cycle-dependent recruitment of FtsN to the Z ring. Our data show rapid accumulation of FtsN to the Z ring about a quarter of the cell cycle after the formation of the Z ring. This initial wave is followed by another increase in FtsN recruitment once the FtsZ protofilament network starts to disassemble. The presence of FtsA* accelerates FtsN recruitment to the Z ring but does not lead to earlier constrictions. Our data furthermore suggest FtsZ and FtsN are part of physically separate complexes throughout the division process.

Coupling between DNA replication, segregation, and the onset of constriction in Escherichia coli.

Escherichia coli cell cycle features two critical cell-cycle checkpoints: initiation of replication and the onset of constriction. While the initiation of DNA replication has been extensively studied, it is less clear what triggers the onset of constriction and when exactly it occurs during the cell cycle. Here, using high-throughput fluorescence microscopy in microfluidic devices, we determine the timing for the onset of constriction relative to the replication cycle in different growth rates. Our single-cell data and modeling indicate that the initiation of constriction is coupled to replication-related processes in slow growth conditions. Furthermore, our data suggest that this coupling involves the mid-cell chromosome blocking the onset of constriction via some form of nucleoid occlusion occurring independently of SlmA and the Ter linkage proteins. This work highlights the coupling between replication and division cycles and brings up a new nucleoid mediated control mechanism in E. coli.

Distinguishing different modes of growth using single-cell data.

Collection of high-throughput data has become prevalent in biology. Large datasets allow the use of statistical constructs such as binning and linear regression to quantify relationships between variables and hypothesize underlying biological mechanisms based on it. We discuss several such examples in relation to single-cell data and cellular growth. In particular, we show instances where what appears to be ordinary use of these statistical methods leads to incorrect conclusions such as growth being non-exponential as opposed to exponential and vice versa. We propose that the data analysis and its interpretation should be done in the context of a generative model, if possible. In this way, the statistical methods can be validated either analytically or against synthetic data generated via the use of the model, leading to a consistent method for inferring biological mechanisms from data. On applying the validated methods of data analysis to infer cellular growth on our experimental data, we find the growth of length in E. coli to be non-exponential. Our analysis shows that in the later stages of the cell cycle the growth rate is faster than exponential.

All cells ­ from bacteria to humans ­ tightly control their size as they grow and divide. Cells can also change the speed at which they grow, and the pattern of how fast a cell grows with time is called 'mode of growth'. Mode of growth can be 'linear', when cells increase their size at a constant rate, or 'exponential', when cells increase their size at a rate proportional to their current size. A cell's mode of growth influences its inner workings, so identifying how a cell grows can reveal information about how a cell will behave. Scientists can measure the size of cells as they age and identify their mode of growth using single cell imaging techniques. Unfortunately, the statistical methods available to analyze the large amounts of data generated in these experiments can lead to incorrect conclusions. Specifically, Kar et al. found that scientists had been using specific types of plots to analyze growth data that were prone to these errors, and may lead to misinterpreting exponential growth as linear and vice versa. This discrepancy can be resolved by ensuring that the plots used to determine the mode of growth are adequate for this analysis. But how can the adequacy of a plot be tested? One way to do this is to generate synthetic data from a known model, which can have a specific and known mode of growth, and using this data to test the different plots. Kar et al. developed such a 'generative model' to produce synthetic data similar to the experimental data, and used these data to determine which plots are best suited to determine growth mode. Once they had validated the best statistical methods for studying mode of growth, Kar et al. applied these methods to growth data from the bacterium Escherichia coli. This showed that these cells have a form of growth called 'super-exponential growth'. These findings identify a strategy to validate statistical methods used to analyze cell growth data. Furthermore, this strategy ­ the use of generative models to produce synthetic data to test the accuracy of statistical methods ­ could be used in other areas of biology to validate statistical approaches.

Transient Membrane-Linked FtsZ Assemblies Precede Z-Ring Formation in Escherichia coli.

During the early stages of cytokinesis, FtsZ protofilaments form a ring-like structure, the Z-ring, in most bacterial species. This cytoskeletal scaffold recruits downstream proteins essential for septal cell wall synthesis. Despite progress in understanding the dynamic nature of the Z-ring and its role in coordinating septal cell wall synthesis, the early stages of protofilament formation and subsequent assembly into the Z-ring are still not understood. Here we investigate a sequence of assembly steps that lead to the formation of the Z-ring in Escherichia coli using high temporal and spatial resolution imaging. Our data show that formation of the Z-ring is preceded by transient membrane-linked FtsZ assemblies. These assemblies form after attachment of short cytosolic protofilaments, which we estimate to be less than 20 monomers long, to the membrane. The attachments occur at random locations along the length of the cell. The filaments treadmill and show periods of rapid growth and shrinkage. Their dynamic properties imply that protofilaments are bundled in these assemblies. Furthermore, we establish that the size of assemblies is sensitively controlled by the availability of FtsZ molecules and by the presence of ZapA proteins. The latter has been implicated in cross-linking the protofilaments. The likely function of these dynamic FtsZ assemblies is to sample the cell surface for the proper location for the Z-ring.

The effects of polydisperse crowders on the compaction of the Escherichia coli nucleoid.

DNA binding proteins, supercoiling, macromolecular crowders, and transient DNA attachments to the cell membrane have all been implicated in the organization of the bacterial chromosome. However, it is unclear what role these factors play in compacting the bacterial DNA into a distinct organelle-like entity, the nucleoid. By analyzing the effects of osmotic shock and mechanical squeezing on Escherichia coli, we show that macromolecular crowders play a dominant role in the compaction of the DNA into the nucleoid. We find that a 30% increase in the crowder concentration from physiological levels leads to a three-fold decrease in the nucleoid's volume. The compaction is anisotropic, being higher along the long axes of the cell at low crowding levels. At higher crowding levels, the nucleoid becomes spherical, and its compressibility decreases significantly. Furthermore, we find that the compressibility of the nucleoid is not significantly affected by cell growth rates and by prior treatment with rifampicin. The latter results point out that in addition to poly ribosomes, soluble cytoplasmic proteins have a significant contribution in determining the size of the nucleoid. The contribution of poly ribosomes dominates at faster and soluble proteins at slower growth rates.

Cell cycle-dependent regulation of FtsZ in Escherichia coli in slow growth conditions.

FtsZ is the key regulator of bacterial cell division. It initiates division by forming a dynamic ring-like structure, the Z-ring, at the mid-cell. What triggers the formation of the Z-ring during the cell cycle is poorly understood. In Escherichia coli, the common view is that FtsZ concentration is constant throughout its doubling time and therefore regulation of assembly is controlled by some yet-to-be-identified protein-protein interactions. Using a newly developed functional, fluorescent FtsZ reporter, we performed a quantitative analysis of the FtsZ concentration throughout the cell cycle under slow growth conditions. In contrast to the common expectation, we show that FtsZ concentrations vary in a cell cycle-dependent manner, and that upregulation of FtsZ synthesis correlates with the formation of the Z-ring. The first half of the cell cycle shows an approximately fourfold upregulation of FtsZ synthesis, followed by its rapid degradation by ClpXP protease in the last 10% of the cell cycle. The initiation of rapid degradation coincides with the dissociation of FtsZ from the septum. Altogether, our data suggest that the Z-ring formation in slow growth conditions in E. coli is partially controlled by a regulatory sequence wherein upregulation of an essential cell cycle factor is followed by its degradation.

Kinetics of large-scale chromosomal movement during asymmetric cell division in Escherichia coli.

Coordination between cell division and chromosome replication is essential for a cell to produce viable progeny. In the commonly accepted view, Escherichia coli realize this coordination via the accurate positioning of its cell division apparatus relative to the nucleoids. However, E. coli lacking proper positioning of its cell division planes can still successfully propagate. Here, we characterize how these cells partition their chromosomes into daughters during such asymmetric divisions. Using quantitative time-lapse imaging, we show that DNA translocase, FtsK, can pump as much as 80% (3.7 Mb) of the chromosome between daughters at an average rate of 1700±800 bp/s. Pauses in DNA translocation are rare, and in no occasions did we observe reversals at experimental time scales of a few minutes. The majority of DNA movement occurs at the latest stages of cell division when the cell division protein ZipA has already dissociated from the septum, and the septum has closed to a narrow channel with a diameter much smaller than the resolution limit of the microscope (~250 nm). Our data suggest that the narrow constriction is necessary for effective translocation of DNA by FtsK.

Lipid Droplets Form from Distinct Regions of the Cell in the Fission Yeast Schizosaccharomyces pombe.

Eukaryotic cells store cholesterol/sterol esters (SEs) and triacylglycerols (TAGs) in lipid droplets, which form from the contiguous endoplasmic reticulum (ER) network. However, it is not known if droplets preferentially form from certain regions of the ER over others. Here, we used fission yeast Schizosaccharomyces pombe cells where the nuclear and cortical/peripheral ER domains are distinguishable by light microscopy to show that SE-enriched lipid droplets form away from the nucleus at the cell tips, whereas TAG-enriched lipid droplets form around the nucleus. Sterols localize to the regions of the cells where droplets enriched in SEs are observed. TAG droplet formation around the nucleus appears to be a strong function of diacylglycerol (DAG) homeostasis with Cpt1p, which coverts DAG into phosphatidylcholine and phosphatidylethanolamine localized exclusively to the nuclear ER. Also, Dgk1p, which converts DAG into phosphatidic acid localized strongly to the nuclear ER over the cortical/peripheral ER. We also show that TAG more readily translocates from the ER to lipid droplets than do SEs. The results augment the standard lipid droplet formation model, which has SEs and TAGs flowing into the same nascent lipid droplet regardless of its biogenesis point in the cell.

The role of MatP, ZapA and ZapB in chromosomal organization and dynamics in Escherichia coli.

Despite extensive research over several decades, a comprehensive view of how the Escherichia coli chromosome is organized within the nucleoid, and how two daughter chromosomes segregate has yet to emerge. Here we investigate the role of the MatP, ZapA and ZapB proteins in organizing the replication terminus (Ter) region and in the chromosomal segregation process. Quantitative image analysis of the fluorescently labeled Ter region shows that the replication terminus attaches to the divisome in a single segment along the perimeter of the cell in a MatP, ZapA and ZapB-dependent manner. The attachment does not significantly affect the bulk chromosome segregation in slow growth conditions. With or without the attachment, two chromosomal masses separate from each other at a speed comparable to the cell growth. The separation starts even before the replication terminus region positions itself at the center of the nucleoid. Modeling of the segregation based on conformational entropy correctly predicts the positioning of the replication terminus region within the nucleoid. However, the model produces a distinctly different chromosomal density distribution than the experiment, indicating that the conformational entropy plays a limited role in segregating the chromosomes in the late stages of replication.

Isolation of cellular lipid droplets: two purification techniques starting from yeast cells and human placentas.

Lipid droplets are dynamic organelles that can be found in most eukaryotic and certain prokaryotic cells. Structurally, the droplets consist of a core of neutral lipids surrounded by a phospholipid monolayer. One of the most useful techniques in determining the cellular roles of droplets has been proteomic identification of bound proteins, which can be isolated along with the droplets. Here, two methods are described to isolate lipid droplets and their bound proteins from two wide-ranging eukaryotes: fission yeast and human placental villous cells. Although both techniques have differences, the main method-- density gradient centrifugation--is shared by both preparations. This shows the wide applicability of the presented droplet isolation techniques. In the first protocol, yeast cells are converted into spheroplasts by enzymatic digestion of their cell walls. The resulting spheroplasts are then gently lysed in a loose-fitting homogenizer. Ficoll is added to the lysate to provide a density gradient, and the mixture is centrifuged three times. After the first spin, the lipid droplets are localized to the white-colored floating layer of the centrifuge tubes along with the endoplasmic reticulum (ER), the plasma membrane, and vacuoles. Two subsequent spins are used to remove these other three organelles. The result is a layer that has only droplets and bound proteins. In the second protocol, placental villous cells are isolated from human term placentas by enzymatic digestion with trypsin and DNase I. The cells are homogenized in a loose-fitting homogenizer. Low-speed and medium-speed centrifugation steps are used to remove unbroken cells, cellular debris, nuclei, and mitochondria. Sucrose is added to the homogenate to provide a density gradient and the mixture is centrifuged to separate the lipid droplets from the other cellular fractions. The purity of the lipid droplets in both protocols is confirmed by Western Blot analysis. The droplet fractions from both preps are suitable for subsequent proteomic and lipidomic analysis.

Mid-gestational gene expression profile in placenta and link to pregnancy complications.

Despite the importance of placenta in mediating rapid physiological changes in pregnancy, data on temporal dynamics of placental gene expression are limited. We completed the first transcriptome profiling of human placental gene expression dynamics (GeneChips, Affymetrix®; ~47,000 transcripts) from early to mid-gestation (n = 10; gestational weeks 5-18) and report 154 genes with significant transcriptional changes (ANOVA, FDR P<0.1). TaqMan RT-qPCR analysis (n = 43; gestational weeks 5-41) confirmed a significant (ANOVA and t-test, FDR P<0.05) mid-gestational peak of placental gene expression for BMP5, CCNG2, CDH11, FST, GATM, GPR183, ITGBL1, PLAGL1, SLC16A10 and STC1, followed by sharp decrease in mRNA levels at term (t-test, FDR P<0.05). We hypothesized that normal course of late pregnancy may be affected when genes characteristic to mid-gestation placenta remain highly expressed until term, and analyzed their expression in term placentas from normal and complicated pregnancies [preeclampsia (PE), n = 12; gestational diabetes mellitus (GDM), n = 12; small- and large-for-gestational-age newborns (SGA, LGA), n = 12+12]. STC1 (stanniocalcin 1) exhibited increased mRNA levels in all studied complications, with the most significant effect in PE- and SGA-groups (t-test, FDR P<0.05). In post-partum maternal plasma, the highest STC1 hormone levels (ELISA, n = 129) were found in women who had developed PE and delivered a SGA newborn (median 731 vs 418 pg/ml in controls; ANCOVA, P = 0.00048). Significantly higher expression (t-test, FDR P<0.05) of CCNG2 and LYPD6 accompanied with enhanced immunostaining of the protein was detected in placental sections of PE and GDM cases (n = 15). Our study demonstrates the importance of temporal dynamics of placental transcriptional regulation across three trimesters of gestation. Interestingly, many genes with high expression in mid-gestation placenta have also been implicated in adult complex disease, promoting the discussion on the role of placenta in developmental programming. The discovery of elevated maternal plasma STC1 in pregnancy complications warrants further investigations of its potential as a biomarker.

WITHDRAWN: TNF-Related Apoptosis-Inducing Ligand TRAIL as a Potential Biomarker for Early Pregnancy Complications.

This paper was withdrawn at the request of the editors due to uncertainties inherent in the statistical analysis.

Differential placental expression profile of human Growth Hormone/Chorionic Somatomammotropin genes in pregnancies with pre-eclampsia and gestational diabetes mellitus.

The human GH/CSH cluster consisting of one pituitary-expressed (GH1) and four placenta-expressed loci has been implicated in maternal metabolic adaptation to pregnancy, regulation of intrauterine and postnatal growth. We investigated how the mRNA expression profile of placental GH2, CSH1 and CSH2 genes and their alternative transcripts correlates with maternal pre-eclampsia (PE) and/or gestational diabetes mellitus (GD). The expression of studied genes in PE placentas (n=17) compared to controls (n=17) exhibited a trend for reduced transcript levels. The alternative transcripts retaining intron 4, GH2-2 and CSH1-2 showed significantly reduced expression in PE cases without growth restriction (P=0.007, P=0.008, respectively). In maternal GD (n=23), a tendency of differential expression was detected only for the GH2 gene and in pregnancies with large-for-gestational-age newborns. Our results, together with those reported by others, are consistent with a pleiotropic effect of placental hGH/CSH genes at the maternal-fetal interface relating to the regulation of fetal growth and the risk of affected maternal metabolism.

Differential expression profile of growth hormone/chorionic somatomammotropin genes in placenta of small- and large-for-gestational-age newborns.

CONTEXT: The human growth hormone/chorionic somatomammotropin (hGH/CSH) locus at 17q22-24, consisting of one pituitary-expressed postnatal (GH1) and four placenta-expressed genes (GH2, CSH1, CSH2, and CSHL1), is implicated in regulation of postnatal and intrauterine growth. A positive correlation has been reported between the offspring's birth weight and serum placental GH (coded by GH2) and placental lactogen (coded by CSH1, CSH2) levels in pregnant women. OBJECTIVE: The objective of the study was the investigation of the hypothesis that the mRNA expression profile of placental hGH/CSH genes contributes to the determination of birth weight. DESIGN AND SUBJECTS: We developed a sensitive, fluorescent-labeled semiquantitative RT-PCR assay coupled with gene-specific restriction analysis, capable of distinguishing alternative splice-products of individual placental hGH/CSH genes and quantification of their relative expression levels. The detailed profile of alternative transcripts of GH2, CSH1, CSH2, and CSHL1 genes in placenta from uncomplicated term pregnancies of the REPROMETA sample collection was addressed in association with the birth weight of newborns, grouped as appropriate for gestational age (AGA; n = 23), small for gestational age (SGA; n = 15), and large for gestational age (LGA; n = 34). RESULTS: The majority of pregnancies with SGA newborn showed down-regulation of the entire hGH/CSH cluster in placenta, whereas in the case of LGA, the expression of CSH1-1, CSH2-1, and CSHL1-4 mRNA transcripts in placenta was significantly increased compared with AGA newborns (P < 0.0001, P = 0.009, P = 0.002, respectively). CONCLUSION: The expression profile of placental hGH/CSH genes in placenta is altered in pregnancies accompanied by SGA and LGA compared with AGA newborns, and thus, it may directly affect the circulating fetal and maternal placental GH and placental lactogen levels.

Regeneration of multilineage skin epithelia by differentiated keratinocytes.

Although homeostasis of rapidly renewing tissues like skin epithelia is maintained by stem cells, the committed progeny of stem cells in the basal layer of epidermis retain regenerative potential and are capable of forming epidermis in response to environmental cues. It is not clear, however, at what point within the epidermal lineage keratinocytes lose this regenerative potential. In this study, we examined the extent of tissue formation by post-mitotic differentiated keratinocytes. We show that cultures of mouse keratinocytes that were, by all measures, differentiated were able to reform a self-renewing, hair-bearing skin when transplanted onto suitable sites in vivo. Genetic labeling and lineage-tracing studies in combination with an involucrin-driven Cre/lox reporter system confirmed that transplanted differentiated keratinocytes were indeed the source of the regenerated skin. More importantly, analysis of early stages of skin regeneration showed hallmarks of dedifferentiation of transplanted differentiated keratinocytes. These data indicate that commitment to differentiation does not prohibit cells from re-entering the cell cycle, de-differentiating, and acquiring "stemness". These findings suggest that epidermis can use different strategies for homeostasis and tissue regeneration.

The epidermal differentiation-associated Grainyhead gene Get1/Grhl3 also regulates urothelial differentiation.

Skin and bladder epithelia form effective permeability barriers through the activation of distinct differentiation gene programs. Using a genome-wide gene-expression study, we identified transcriptional regulators whose expression correlates highly with that of differentiation markers in both the bladder and skin, including the Grainyhead factor Get1/Grhl3, which is already known to be important for epidermal barrier formation. In the bladder, Get1 is most highly expressed in the differentiated umbrella cells and its mutation in mice leads to a defective bladder epithelial barrier formation due to the failure of apical membrane specialization. Genes encoding components of the specialized urothelial membrane, the uroplakins, were downregulated in Get1(-/-) mice. At least one of these genes, uroplakin II, is a direct target of Get1. The urothelial-specific activation of the uroplakin II gene is due to selective binding of Get1 to the uroplakin II promoter in urothelial cells, which is most likely regulated by histone modifications. These results show a crucial role for Get1 in urothelial differentiation and barrier formation.

Grainyhead-like factor Get1/Grhl3 regulates formation of the epidermal leading edge during eyelid closure.

Grainyhead transcription factors play an evolutionarily conserved role in regulating epidermal terminal differentiation. One such factor, the mammalian Grainyhead-like epithelial transactivator (Get1/Grhl3), is important for epidermal barrier formation. In addition to a role in barrier formation, Grainyhead genes play roles in closure of several structures such as the mouse neural tube and Drosophila wounds. Consistent with these observations, we found that Get1 knockout mice have an eye-open at birth phenotype. The failure of eyelid closure appears to be due to critical functions of Get1 in promoting F-actin polymerization, filopodia formation, and the cell shape changes that are required for migration of the keratinocytes at the leading edge during eyelid closure. The expression of TGFalpha, a known regulator of leading edge formation, is decreased in the eyelid tip of Get1(-/-) mice. Levels of phospho-EGFR and phospho-ERK are also decreased at the leading edge tip. Furthermore, in an organ culture model, TGFalpha can increase levels of phospho-EGFR and promote cell shape changes as well as leading edge formation in Get1(-/-) eyelids, indicating that in eyelid closure Get1 acts upstream of TGFalpha in the EGFR/ERK pathway.

Co-factors of LIM domains (Clims/Ldb/Nli) regulate corneal homeostasis and maintenance of hair follicle stem cells.

The homeostasis of both cornea and hair follicles depends on a constant supply of progeny cells produced by populations of keratin (K) 14-expressing stem cells localized in specific niches. To investigate the potential role of Co-factors of LIM domains (Clims) in epithelial tissues, we generated transgenic mice expressing a dominant-negative Clim molecule (DN-Clim) under the control of the K14 promoter. As expected, the K14 promoter directed high level expression of the transgene to the basal cells of cornea and epidermis, as well as the outer root sheath of hair follicles. In corneal epithelium, the transgene expression causes decreased expression of adhesion molecule BP180 and defective hemidesmosomes, leading to detachment of corneal epithelium from the underlying stroma, which in turn causes blisters, wounds and an inflammatory response. After a period of epithelial thinning, the corneal epithelium undergoes differentiation to an epidermis-like structure. The K14-DN-Clim mice also develop progressive hair loss due to dysfunctional hair follicles that fail to generate hair shafts. The number of hair follicle stem cells is decreased by at least 60% in K14-DN-Clim mice, indicating that Clims are required for hair follicle stem cell maintenance. In addition, Clim2 interacts with Lhx2 in vivo, suggesting that Clim2 is an essential co-factor for the LIM homeodomain factor Lhx2, which was previously shown to play a role in hair follicle stem cell maintenance. Together, these data indicate that Clim proteins play important roles in the homeostasis of corneal epithelium and hair follicles.