The mammalian Msx homeobox genes, Msx1 and Msx2, encode transcription factors that control organogenesis and tissue interactions during embryonic development. We observed overlapping expression of these factors in uterine epithelial and stromal compartments of pregnant mice prior to embryo implantation. Conditional ablation of both Msx1 and Msx2 in the uterus resulted in female infertility due to a failure in implantation. In these mutant mice (Msx1/2(d/d)), the uterine epithelium exhibited persistent proliferative activity and failed to attach to the embryos. Gene expression profiling of uterine epithelium and stroma of Msx1/2(d/d) mice revealed an elevated expression of several members of the Wnt gene family in the preimplantation uterus. Increased canonical Wnt signaling in the stromal cells activated β-catenin, stimulating the production of a subset of fibroblast growth factors (FGFs) in these cells. The secreted FGFs acted in a paracrine manner via the FGF receptors in the epithelium to promote epithelial proliferation, thereby preventing differentiation of this tissue and creating a non-receptive uterus refractory to implantation. Collectively, these findings delineate a unique signaling network, involving Msx1/2, Wnts, and FGFs, which operate in the uterus at the time of implantation to control the mesenchymal-epithelial dialogue critical for successful establishment of pregnancy.

Gene duplications within the conserved Hox cluster are rare in animal evolution, but in Lepidoptera an array of divergent Hox-related genes (Shx genes) has been reported between pb and zen. Here, we use genome sequencing of five lepidopteran species (Polygonia c-album, Pararge aegeria, Callimorpha dominula, Cameraria ohridella, Hepialus sylvina) plus a caddisfly outgroup (Glyphotaelius pellucidus) to trace the evolution of the lepidopteran Shx genes. We demonstrate that Shx genes originated by tandem duplication of zen early in the evolution of large clade Ditrysia; Shx are not found in a caddisfly and a member of the basally diverging Hepialidae (swift moths). Four distinct Shx genes were generated early in ditrysian evolution, and were stably retained in all descendent Lepidoptera except the silkmoth which has additional duplications. Despite extensive sequence divergence, molecular modelling indicates that all four Shx genes have the potential to encode stable homeodomains. The four Shx genes have distinct spatiotemporal expression patterns in early development of the Speckled Wood butterfly (Pararge aegeria), with ShxC demarcating the future sites of extraembryonic tissue formation via strikingly localised maternal RNA in the oocyte. All four genes are also expressed in presumptive serosal cells, prior to the onset of zen expression. Lepidopteran Shx genes represent an unusual example of Hox cluster expansion and integration of novel genes into ancient developmental regulatory networks.

Planaria continue to blossom as a model system for understanding all aspects of regeneration. They provide an opportunity to understand how the replacement of missing tissues from preexisting adult tissue is orchestrated at the molecular level. When amputated along any plane, planaria are capable of regenerating all missing tissue and rescaling all structures to the new size of the animal. Recently, rapid progress has been made in understanding the developmental pathways that control planarian regeneration. In particular Wnt/beta-catenin signaling is central in promoting posterior fates and inhibiting anterior identity. Currently the mechanisms that actively promote anterior identity remain unknown. Here, Smed-prep, encoding a TALE class homeodomain, is described as the first gene necessary for correct anterior fate and patterning during planarian regeneration. Smed-prep is expressed at high levels in the anterior portion of whole animals, and Smed-prep(RNAi) leads to loss of the whole brain during anterior regeneration, but not during lateral regeneration or homeostasis in intact worms. Expression of markers of different anterior fated cells are greatly reduced or lost in Smed-prep(RNAi) animals. We find that the ectopic anterior structures induced by abrogation of Wnt signaling also require Smed-prep to form. We use double knockdown experiments with the S. mediterranea ortholog of nou-darake (that when knocked down induces ectopic brain formation) to show that Smed-prep defines an anterior fated compartment within which stem cells are permitted to assume brain fate, but is not required directly for this differentiation process. Smed-prep is the first gene clearly implicated as being necessary for promoting anterior fate and the first homeobox gene implicated in establishing positional identity during regeneration. Together our results suggest that Smed-prep is required in stem cell progeny as they form the anterior regenerative blastema and is required for specifying anterior cell fates and correct patterning.

Unlike many wild grasses, domesticated rice cultivars have uniform culm height and panicle size among tillers and the main shoot, which is an important trait for grain yield. However, the genetic basis of this trait remains unknown. Here, we report that Dwarf Tiller1 (DWT1) controls the developmental uniformity of the main shoot and tillers in rice (Oryza sativa). Most dwt1 mutant plants develop main shoots with normal height and larger panicles, but dwarf tillers bearing smaller panicles compared with those of the wild type. In addition, dwt1 tillers have shorter internodes with fewer and un-elongated cells compared with the wild type, indicating that DWT1 affects cell division and cell elongation. Map-based cloning revealed that DWT1 encodes a Wuschel-related homeobox (WOX) transcription factor homologous to the Arabidopsis WOX8 and WOX9. The DWT1 gene is highly expressed in young panicles, but undetectable in the internodes, suggesting that DWT1 expression is spatially or temporally separated from its effect on the internode growth. Transcriptomic analysis revealed altered expression of genes involved in cell division and cell elongation, cytokinin/gibberellin homeostasis and signaling in dwt1 shorter internodes. Moreover, the non-elongating internodes of dwt1 are insensitive to exogenous gibberellin (GA) treatment, and some of the slender rice1 (slr1) dwt1 double mutant exhibits defective internodes similar to the dwt1 single mutant, suggesting that the DWT1 activity in the internode elongation is directly or indirectly associated with GA signaling. This study reveals a genetic pathway synchronizing the development of tillers and the main shoot, and a new function of WOX genes in balancing branch growth in rice.

The appropriate development of conidia and appressoria is critical in the disease cycle of many fungal pathogens, including Magnaporthe oryzae. A total of eight genes (MoHOX1 to MoHOX8) encoding putative homeobox transcription factors (TFs) were identified from the M. oryzae genome. Knockout mutants for each MoHOX gene were obtained via homology-dependent gene replacement. Two mutants, DeltaMohox3 and DeltaMohox5, exhibited no difference to wild-type in growth, conidiation, conidium size, conidial germination, appressorium formation, and pathogenicity. However, the DeltaMohox1 showed a dramatic reduction in hyphal growth and increase in melanin pigmentation, compared to those in wild-type. DeltaMohox4 and DeltaMohox6 showed significantly reduced conidium size and hyphal growth, respectively. DeltaMohox8 formed normal appressoria, but failed in pathogenicity, probably due to defects in the development of penetration peg and invasive growth. It is most notable that asexual reproduction was completely abolished in DeltaMohox2, in which no conidia formed. DeltaMohox2 was still pathogenic through hypha-driven appressoria in a manner similar to that of the wild-type. However, DeltaMohox7 was unable to form appressoria either on conidial germ tubes, or at hyphal tips, being non-pathogenic. These factors indicate that M. oryzae is able to cause foliar disease via hyphal appressorium-mediated penetration, and MoHOX7 is mutually required to drive appressorium formation from hyphae and germ tubes. Transcriptional analyses suggest that the functioning of M. oryzae homeobox TFs is mediated through the regulation of gene expression and is affected by cAMP and Ca(2+) signaling and/or MAPK pathways. The divergent roles of this gene set may help reveal how the genome and regulatory pathways evolved within the rice blast pathogen and close relatives.

The differentiation efficiency of adult stem cells undergoes a significant decline in aged animals, which is closely related to the decline in organ function and age-associated diseases. However, the underlying mechanisms that ultimately lead to this observed decline of the differentiation efficiency of stem cells remain largely unclear. This study investigated Drosophila midguts and identified an obvious upregulation of caudal (cad), which encodes a homeobox transcription factor. This factor is traditionally known as a central regulator of embryonic anterior-posterior body axis patterning. This study reports that depletion of cad in intestinal stem/progenitor cells promotes quiescent intestinal stem cells (ISCs) to become activate and produce enterocytes in the midgut under normal gut homeostasis conditions. However, overexpression of cad results in the failure of ISC differentiation and intestinal epithelial regeneration after injury. Moreover, this study suggests that cad prevents intestinal stem/progenitor cell differentiation by modulating the Janus kinase/signal transducers and activators of the transcription pathway and Sox21a-GATAe signaling cascade. Importantly, the reduction of cad expression in intestinal stem/progenitor cells restrained age-associated gut hyperplasia in Drosophila. This study identified a function of the homeobox gene cad in the modulation of adult stem cell differentiation and suggested a potential gene target for the treatment of age-related diseases induced by age-related stem cell dysfunction.

Neural tube defects (NTDs) is a general term for central nervous system malformations secondary to a failure of closure or development of the neural tube. The resulting pathologies may involve the brain, spinal cord and/or vertebral column, in addition to associated structures such as soft tissue or skin. The condition is reported among the more common birth defects in humans, leading to significant infant morbidity and mortality. The etiology remains poorly understood but genetic, nutritional, environmental factors, or a combination of these, are known to play a role in the development of NTDs. The variable conditions associated with NTDs occur naturally in dogs, and have been previously reported in the Weimaraner breed. Taking advantage of the strong linkage-disequilibrium within dog breeds we performed genome-wide association analysis and mapped a genomic region for spinal dysraphism, a presumed NTD, using 4 affected and 96 unaffected Weimaraners. The associated region on canine chromosome 8 (pgenome  =3.0 × 10(-5)), after 100,000 permutations, encodes 18 genes, including NKX2-8, a homeobox gene which is expressed in the developing neural tube. Sequencing NKX2-8 in affected Weimaraners revealed a G to AA frameshift mutation within exon 2 of the gene, resulting in a premature stop codon that is predicted to produce a truncated protein. The exons of NKX2-8 were sequenced in human patients with spina bifida and rare variants (rs61755040 and rs10135525) were found to be significantly over-represented (p=0.036). This is the first documentation of a potential role for NKX2-8 in the etiology of NTDs, made possible by investigating the molecular basis of naturally occurring mutations in dogs.

Previously, we identified an adolescent idiopathic scoliosis susceptibility locus near human ladybird homeobox 1 (LBX1) and FLJ41350 by a genome-wide association study. Here, we characterized the associated non-coding variant and investigated the function of these genes. A chromosome conformation capture assay revealed that the genome region with the most significantly associated single nucleotide polymorphism (rs11190870) physically interacted with the promoter region of LBX1-FLJ41350. The promoter in the direction of LBX1, combined with a 590-bp region including rs11190870, had higher transcriptional activity with the risk allele than that with the non-risk allele in HEK 293T cells. The ubiquitous overexpression of human LBX1 or either of the zebrafish lbx genes (lbx1a, lbx1b, and lbx2), but not FLJ41350, in zebrafish embryos caused body curvature followed by death prior to vertebral column formation. Such body axis deformation was not observed in transcription activator-like effector nucleases mediated knockout zebrafish of lbx1b or lbx2. Mosaic expression of lbx1b driven by the GATA2 minimal promoter and the lbx1b enhancer in zebrafish significantly alleviated the embryonic lethal phenotype to allow observation of the later onset of the spinal curvature with or without vertebral malformation. Deformation of the embryonic body axis by lbx1b overexpression was associated with defects in convergent extension, which is a component of the main axis-elongation machinery in gastrulating embryos. In embryos overexpressing lbx1b, wnt5b, a ligand of the non-canonical Wnt/planar cell polarity (PCP) pathway, was significantly downregulated. Injection of mRNA for wnt5b or RhoA, a key downstream effector of Wnt/PCP signaling, rescued the defective convergent extension phenotype and attenuated the lbx1b-induced curvature of the body axis. Thus, our study presents a novel pathological feature of LBX1 and its zebrafish homologs in body axis deformation at various stages of embryonic and subsequent growth in zebrafish.

Neofunctionalization following gene duplication is thought to be one of the key drivers in generating evolutionary novelty. A gene duplication in a common ancestor of land plants produced two classes of KNOTTED-like TALE homeobox genes, class I (KNOX1) and class II (KNOX2). KNOX1 genes are linked to tissue proliferation and maintenance of meristematic potentials of flowering plant and moss sporophytes, and modulation of KNOX1 activity is implicated in contributing to leaf shape diversity of flowering plants. While KNOX2 function has been shown to repress the gametophytic (haploid) developmental program during moss sporophyte (diploid) development, little is known about KNOX2 function in flowering plants, hindering syntheses regarding the relationship between two classes of KNOX genes in the context of land plant evolution. Arabidopsis plants harboring loss-of-function KNOX2 alleles exhibit impaired differentiation of all aerial organs and have highly complex leaves, phenocopying gain-of-function KNOX1 alleles. Conversely, gain-of-function KNOX2 alleles in conjunction with a presumptive heterodimeric BELL TALE homeobox partner suppressed SAM activity in Arabidopsis and reduced leaf complexity in the Arabidopsis relative Cardamine hirsuta, reminiscent of loss-of-function KNOX1 alleles. Little evidence was found indicative of epistasis or mutual repression between KNOX1 and KNOX2 genes. KNOX proteins heterodimerize with BELL TALE homeobox proteins to form functional complexes, and contrary to earlier reports based on in vitro and heterologous expression, we find high selectivity between KNOX and BELL partners in vivo. Thus, KNOX2 genes confer opposing activities rather than redundant roles with KNOX1 genes, and together they act to direct the development of all above-ground organs of the Arabidopsis sporophyte. We infer that following the KNOX1/KNOX2 gene duplication in an ancestor of land plants, neofunctionalization led to evolution of antagonistic biochemical activity thereby facilitating the evolution of more complex sporophyte transcriptional networks, providing plasticity for the morphological evolution of land plant body plans.

The Rhox cluster on the mouse X chromosome contains reproduction-related homeobox genes expressed in a sexually dimorphic manner. We report that two members of the Rhox cluster, Rhox6 and 9, are regulated by de-methylation of histone H3 at lysine 27 by KDM6A, a histone demethylase with female-biased expression. Consistent with other homeobox genes, Rhox6 and 9 are in bivalent domains prior to embryonic stem cell differentiation and thus poised for activation. In female mouse ES cells, KDM6A is specifically recruited to Rhox6 and 9 for gene activation, a process inhibited by Kdm6a knockdown in a dose-dependent manner. In contrast, KDM6A occupancy at Rhox6 and 9 is low in male ES cells and knockdown has no effect on expression. In mouse ovary where Rhox6 and 9 remain highly expressed, KDM6A occupancy strongly correlates with expression. Our study implicates Kdm6a, a gene that escapes X inactivation, in the regulation of genes important in reproduction, suggesting that KDM6A may play a role in the etiology of developmental and reproduction-related effects of X chromosome anomalies.

Cleft palate (CP) is one of the most commonly occurring craniofacial birth defects in humans. In order to study cleft palate in a naturally occurring model system, we utilized the Nova Scotia Duck Tolling Retriever (NSDTR) dog breed. Micro-computed tomography analysis of CP NSDTR craniofacial structures revealed that these dogs exhibit defects similar to those observed in a recognizable subgroup of humans with CP: Pierre Robin Sequence (PRS). We refer to this phenotype in NSDTRs as CP1. Individuals with PRS have a triad of birth defects: shortened mandible, posteriorly placed tongue, and cleft palate. A genome-wide association study in 14 CP NSDTRs and 72 unaffected NSDTRs identified a significantly associated region on canine chromosome 14 (24.2 Mb-29.3 Mb; p(raw )= 4.64 × 10(-15)). Sequencing of two regional candidate homeobox genes in NSDTRs, distal-less homeobox 5 (DLX5) and distal-less homeobox 6 (DLX6), identified a 2.1 kb LINE-1 insertion within DLX6 in CP1 NSDTRs. The LINE-1 insertion is predicted to insert a premature stop codon within the homeodomain of DLX6. This prompted the sequencing of DLX5 and DLX6 in a human cohort with CP, where a missense mutation within the highly conserved DLX5 homeobox of a patient with PRS was identified. This suggests the involvement of DLX5 in the development of PRS. These results demonstrate the power of the canine animal model as a genetically tractable approach to understanding naturally occurring craniofacial birth defects in humans.

A central issue of myogenesis is the acquisition of identity by individual muscles. In Drosophila, at the time muscle progenitors are singled out, they already express unique combinations of muscle identity genes. This muscle code results from the integration of positional and temporal signalling inputs. Here we identify, by means of loss-of-function and ectopic expression approaches, the Iroquois Complex homeobox genes araucan and caupolican as novel muscle identity genes that confer lateral transverse muscle identity. The acquisition of this fate requires that Araucan/Caupolican repress other muscle identity genes such as slouch and vestigial. In addition, we show that Caupolican-dependent slouch expression depends on the activation state of the Ras/Mitogen Activated Protein Kinase cascade. This provides a comprehensive insight into the way Iroquois genes integrate in muscle progenitors, signalling inputs that modulate gene expression and protein activity.

Gliomas, the most common malignant tumors of the nervous system, frequently harbor mutations that activate the epidermal growth factor receptor (EGFR) and phosphatidylinositol-3 kinase (PI3K) signaling pathways. To investigate the genetic basis of this disease, we developed a glioma model in Drosophila. We found that constitutive coactivation of EGFR-Ras and PI3K pathways in Drosophila glia and glial precursors gives rise to neoplastic, invasive glial cells that create transplantable tumor-like growths, mimicking human glioma. Our model represents a robust organotypic and cell-type-specific Drosophila cancer model in which malignant cells are created by mutations in signature genes and pathways thought to be driving forces in a homologous human cancer. Genetic analyses demonstrated that EGFR and PI3K initiate malignant neoplastic transformation via a combinatorial genetic network composed primarily of other pathways commonly mutated or activated in human glioma, including the Tor, Myc, G1 Cyclins-Cdks, and Rb-E2F pathways. This network acts synergistically to coordinately stimulate cell cycle entry and progression, protein translation, and inappropriate cellular growth and migration. In particular, we found that the fly orthologs of CyclinE, Cdc25, and Myc are key rate-limiting genes required for glial neoplasia. Moreover, orthologs of Sin1, Rictor, and Cdk4 are genes required only for abnormal neoplastic glial proliferation but not for glial development. These and other genes within this network may represent important therapeutic targets in human glioma.

CRISPR-based homing gene drives can be designed to disrupt essential genes whilst biasing their own inheritance, leading to suppression of mosquito populations in the laboratory. This class of gene drives relies on CRISPR-Cas9 cleavage of a target sequence and copying ('homing') therein of the gene drive element from the homologous chromosome. However, target site mutations that are resistant to cleavage yet maintain the function of the essential gene are expected to be strongly selected for. Targeting functionally constrained regions where mutations are not easily tolerated should lower the probability of resistance. Evolutionary conservation at the sequence level is often a reliable indicator of functional constraint, though the actual level of underlying constraint between one conserved sequence and another can vary widely. Here we generated a novel adult lethal gene drive (ALGD) in the malaria vector Anopheles gambiae, targeting an ultra-conserved target site in a haplosufficient essential gene (AGAP029113) required during mosquito development, which fulfils many of the criteria for the target of a population suppression gene drive. We then designed a selection regime to experimentally assess the likelihood of generation and subsequent selection of gene drive resistant mutations at its target site. We simulated, in a caged population, a scenario where the gene drive was approaching fixation, where selection for resistance is expected to be strongest. Continuous sampling of the target locus revealed that a single, restorative, in-frame nucleotide substitution was selected. Our findings show that ultra-conservation alone need not be predictive of a site that is refractory to target site resistance. Our strategy to evaluate resistance in vivo could help to validate candidate gene drive targets for their resilience to resistance and help to improve predictions of the invasion dynamics of gene drives in field populations.

BREVIPEDICELLUS (BP or KNAT1), a class-I KNOTTED1-like homeobox (KNOX) transcription factor in Arabidopsis thaliana, contributes to shaping the normal inflorescence architecture through negatively regulating other two class-I KNOX genes, KNAT2 and KNAT6. However, the molecular mechanism of BP-mediated transcription regulation remains unclear. In this study, we showed that BP directly interacts with the SWI2/SNF2 chromatin remodeling ATPase BRAHMA (BRM) both in vitro and in vivo. Loss-of-function BRM mutants displayed inflorescence architecture defects, with clustered inflorescences, horizontally orientated pedicels, and short pedicels and internodes, a phenotype similar to the bp mutants. Furthermore, the transcript levels of KNAT2 and KNAT6 were elevated in brm-3, bp-9 and brm-3 bp-9 double mutants. Increased histone H3 lysine 4 tri-methylation (H3K4me3) levels were detected in brm-3, bp-9 and brm-3 bp-9 double mutants. Moreover, BRM and BP co-target to KNAT2 and KNAT6 genes, and BP is required for the binding of BRM to KNAT2 and KNAT6. Taken together, our results indicate that BP interacts with the chromatin remodeling factor BRM to regulate the expression of KNAT2 and KNAT6 in control of inflorescence architecture.

We have characterized the biological functions of the chromatin remodeling protein Bptf (Bromodomain PHD-finger Transcription Factor), the largest subunit of NURF (Nucleosome Remodeling Factor) in a mammal. Bptf mutants manifest growth defects at the post-implantation stage and are reabsorbed by E8.5. Histological analyses of lineage markers show that Bptf(-/-) embryos implant but fail to establish a functional distal visceral endoderm. Microarray analysis at early stages of differentiation has identified Bptf-dependent gene targets including homeobox transcriptions factors and genes essential for the development of ectoderm, mesoderm, and both definitive and visceral endoderm. Differentiation of Bptf(-/-) embryonic stem cell lines into embryoid bodies revealed its requirement for development of mesoderm, endoderm, and ectoderm tissue lineages, and uncovered many genes whose activation or repression are Bptf-dependent. We also provide functional and physical links between the Bptf-containing NURF complex and the Smad transcription factors. These results suggest that Bptf may co-regulate some gene targets of this pathway, which is essential for establishment of the visceral endoderm. We conclude that Bptf likely regulates genes and signaling pathways essential for the development of key tissues of the early mouse embryo.

The expression of specific transcription factors determines the differentiated features of postmitotic neurons. However, the mechanism by which specific molecules determine neuronal cell fate and the extent to which the functions of transcription factors are conserved in evolution are not fully understood. In C. elegans, the cholinergic and peptidergic SMB sensory/inter/motor neurons innervate muscle quadrants in the head and control the amplitude of sinusoidal movement. Here we show that the LIM homeobox protein LIM-4 determines neuronal characteristics of the SMB neurons. In lim-4 mutant animals, expression of terminal differentiation genes, such as the cholinergic gene battery and the flp-12 neuropeptide gene, is completely abolished and thus the function of the SMB neurons is compromised. LIM-4 activity promotes SMB identity by directly regulating the expression of the SMB marker genes via a distinct cis-regulatory motif. Two human LIM-4 orthologs, LHX6 and LHX8, functionally substitute for LIM-4 in C. elegans. Furthermore, C. elegans LIM-4 or human LHX6 can induce cholinergic and peptidergic characteristics in the human neuronal cell lines. Our results indicate that the evolutionarily conserved LIM-4/LHX6 homeodomain proteins function in generation of precise neuronal subtypes.

In animals, circadian rhythms are driven by oscillations in transcription, translation, and proteasomal degradation of highly conserved genes, resulting in diel cycles in the expression of numerous clock-regulated genes. Transcription is largely regulated through the binding of transcription factors to cis-regulatory elements within accessible regions of the chromatin. Chromatin remodeling is linked to circadian regulation in mammals, but it is unknown whether cycles in chromatin accessibility are a general feature of clock-regulated genes throughout evolution. To assess this, we applied an ATAC-seq approach using Nematostella vectensis, grown under two separate light regimes (light:dark (LD) and constant darkness (DD)). Based on previously identified N. vectensis circadian genes, our results show the coupling of chromatin accessibility and circadian transcription rhythmicity under LD conditions. Out of 180 known circadian genes, we were able to list 139 gene promoters that were highly accessible compared to common promoters. Furthermore, under LD conditions, we identified 259 active enhancers as opposed to 333 active enhancers under DD conditions, with 171 enhancers shared between the two treatments. The development of a highly reproducible ATAC-seq protocol integrated with published RNA-seq and ChIP-seq databases revealed the enrichment of transcription factor binding sites (such as C/EBP, homeobox, and MYB), which have not been previously associated with circadian signaling in cnidarians. These results provide new insight into the regulation of cnidarian circadian machinery. Broadly speaking, this supports the notion that the association between chromatin remodeling and circadian regulation arose early in animal evolution as reflected in this non-bilaterian lineage.

Microarrays enable comparative analyses of gene expression on a genomic scale, however these experiments frequently identify an abundance of differentially expressed genes such that it may be difficult to identify discrete functional networks that are hidden within large microarray datasets. Microarray analyses in which mutant organisms are compared to nonmutant siblings can be especially problematic when the gene of interest is expressed in relatively few cells. Here, we describe the use of laser microdissection microarray to perform transcriptional profiling of the maize shoot apical meristem (SAM), a ~100-microm pillar of organogenic cells that is required for leaf initiation. Microarray analyses compared differential gene expression within the SAM and incipient leaf primordium of nonmutant and narrow sheath mutant plants, which harbored mutations in the duplicate genes narrow sheath1 (ns1) and narrow sheath2 (ns2). Expressed in eight to ten cells within the SAM, ns1 and ns2 encode paralogous WUSCHEL1-like homeobox (WOX) transcription factors required for recruitment of leaf initials that give rise to a large lateral domain within maize leaves. The data illustrate the utility of laser microdissection-microarray analyses to identify a relatively small number of genes that are differentially expressed within the SAM. Moreover, these analyses reveal potentially conserved WOX gene functions and implicate specific hormonal and signaling pathways during early events in maize leaf development.

Rice (Oryza sativa) has long and narrow leaves with parallel veins, similar to other grasses. Relative to Arabidopsis thaliana which has oval-shaped leaves, our understanding of the mechanism of leaf development is insufficient in grasses. In this study, we show that OsWOX4, a member of the WUSCHEL-RELATED HOMEOBOX gene family, plays important roles in early leaf development in rice. Inducible downregulation of OsWOX4 resulted in severe defects in leaf development, such as an arrest of vascular differentiation, a partial defect in the early cell proliferation required for midrib formation, and a failure to maintain cellular activity in general parenchyma cells. In situ analysis showed that knockdown of OsWOX4 reduced the expression of two LONELY GUY genes, which function in the synthesis of active cytokinin, in developing vascular bundles. Consistent with this, cytokinin levels were downregulated by OsWOX4 knockdown. Transcriptome analysis further showed that OsWOX4 regulates multiple genes, including those responsible for cell cycle progression and hormone action, consistent with the effects of OsWOX4 downregulation on leaf phenotypes. Collectively, these results suggest that OsWOX4 acts as a key regulator at an early stage of leaf development. Our previous work revealed that OsWOX4 is involved in the maintenance of shoot apical meristem in rice, whereas AtWOX4 is specifically associated with the maintenance of vascular stem cells in Arabidopsis. Thus, the function of the two orthologous genes seems to be diversified between rice and Arabidopsis.