Plant mRNAs move into a fungal pathogen via extracellular vesicles to reduce infection.

Cross-kingdom small RNA trafficking between hosts and microbes modulates gene expression in the interacting partners during infection. However, whether other RNAs are also transferred is unclear. Here, we discover that host plant Arabidopsis thaliana delivers mRNAs via extracellular vesicles (EVs) into the fungal pathogen Botrytis cinerea. A fluorescent RNA aptamer reporter Broccoli system reveals host mRNAs in EVs and recipient fungal cells. Using translating ribosome affinity purification profiling and polysome analysis, we observe that delivered host mRNAs are translated in fungal cells. Ectopic expression of two transferred host mRNAs in B. cinerea shows that their proteins are detrimental to infection. Arabidopsis knockout mutants of the genes corresponding to these transferred mRNAs are more susceptible. Thus, plants have a strategy to reduce infection by transporting mRNAs into fungal cells. mRNAs transferred from plants to pathogenic fungi are translated to compromise infection, providing knowledge that helps combat crop diseases.

Extracellular vesicles: cross-organismal RNA trafficking in plants, microbes, and mammalian cells.

Extracellular vesicles (EVs) are membrane-enclosed nanometer-scale particles that transport biological materials such as RNAs, proteins, and metabolites. EVs have been discovered in nearly all kingdoms of life as a form of cellular communication across different cells and between interacting organisms. EV research has primarily focused on EV-mediated intra-organismal transport in mammals, which has led to the characterization of a plethora of EV contents from diverse cell types with distinct and impactful physiological effects. In contrast, research into EV-mediated transport in plants has focused on inter-organismal interactions between plants and interacting microbes. However, the overall molecular content and functions of plant and microbial EVs remain largely unknown. Recent studies into the plant-pathogen interface have demonstrated that plants produce and secrete EVs that transport small RNAs into pathogen cells to silence virulence-related genes. Plant-interacting microbes such as bacteria and fungi also secrete EVs which transport proteins, metabolites, and potentially RNAs into plant cells to enhance their virulence. This review will focus on recent advances in EV-mediated communications in plant-pathogen interactions compared to the current state of knowledge of mammalian EV capabilities and highlight the role of EVs in cross-kingdom RNA interference.

Function and molecular mechanism of a poplar placenta limited MIXTA gene in regulating differentiation of plant epidermal cells.

The placenta in fruits of most plants either desiccate and shrink as the fruits mature or develop further to form the fleshy tissues. In poplars, placental epidermal cells protrude collectively to produce catkin fibers. In this study, three carpel limited MIXTA genes, PdeMIXTA02, PdeMIXTA03, PdeMIXTA04, were find to specifically expressed in carpel immediately after pollination. Heterologous expression of the three genes in Arabidopsis demonstrated that PdeMIXTA04 significantly promoted trichomes density and could restore trichomes in the trichomeless mutant. By contrast, such functions were not observed with PdeMIXTA02, PdeMIXTA03. In situ hybridization revealed that PdeMIXTA04 was explicitly expressed in poplar placental epidermal cells. We also confirmed trichome-specific expression of the PdeMIXTA04 promoter. Multiple experimental proofs have confirmed the interaction between PdeMIXTA04, PdeMYC and PdeWD40, indicating PdeMIXTA04 functioned through the MYB-bHLH-WD40 ternary complex. Our work provided distinctive understanding of the molecular mechanism triggering differentiation of poplar catkins.

Genome-Wide Analysis of the Expansin Gene Family in Populus and Characterization of Expression Changes in Response to Phytohormone (Abscisic Acid) and Abiotic (Low-Temperature) Stresses.

Expansins are a group of cell wall enzyme proteins that help to loosen cell walls by breaking hydrogen bonds between cellulose microfibrils and hemicellulose. Expansins are essential plant proteins that are involved in several key processes, including seed germination, the growth of pollen tubes and root hairs, fruit ripening and abscission processes. Currently, there is a lack of knowledge concerning the role of expansins in woody plants. In this study, we analyzed expansin genes using Populus genome as the study target. Thirty-six members of the expansin gene family were identified in Populus that were divided into four subfamilies (EXPA, EXPB, EXLA and EXLB). We analyzed the molecular structure, chromosome localization, evolutionary relationships and tissue specificity of these genes and investigated expression changes in responses to phytohormone and abiotic stresses of the expansin genes of Populus tremula L. (PtEXs). Molecular structure analysis revealed that each PtEX protein had several conserved motifs and all of the PtEXs genes had multiple exons. Chromosome structure analysis showed that the expansin gene family is distributed on 14 chromosomes. The PtEXs gene family expansion patterns showed segmental duplication. Transcriptome data of Populus revealed that 36 PtEXs genes were differently expressed in different tissues. Cis-element analysis showed that the PtEXs were closely associated with plant development and responses to phytohormone and abiotic stress. Quantitative real-time PCR showed that abscisic acid (ABA) and low-temperature treatment affected the expression of some PtEXs genes, suggesting that these genes are involved in responses to phytohormone and abiotic stress. This study provides a further understanding of the expansin gene family in Populus and forms a basis for future functional research studies.

Genome-Wide Comparative Analysis of the Fasciclin-like Arabinogalactan Proteins (FLAs) in Salicacea and Identification of Secondary Tissue Development-Related Genes.

Fasciclin-like arabinogalactan proteins (FLAs) are a subclass of arabinogalactan proteins (AGPs) containing both AGP-like glycated domains and fasciclin (FAS) domains, which are involved in plant growth and development and synthesis of the cell wall. However, these proteins have not been identified or analyzed in willow, Salix, the sister genus of Populus. In this study, we performed a whole genome study of the FLA gene family of Salix suchowensis and compared it with the FLA gene family of Populus deltoides. The results showed the presence of 40 and 46 FLA genes in P. deltoides and S. suchowensis, distributed on 17 and 16 chromosomes, respectively. Four pairs of tandem repeat genes were found in willow, while poplar had no tandem repeat genes. Twelve and thirteen pairs of duplicated gene fragments were identified in poplar and willow, respectively. The multispecies phylogenetic tree showed that the FLA gene family could be divided into four groups (I-IV), with Group 1 showing significant expansion in woody plants. A gene expression analysis showed that PdeFLA19/27 in Group I of poplar was highly expressed, specifically during the secondary growth period of the stem and the rapid elongation of seed hairs. In the Group I genes of S. suchowensis, SsuFLA25/26/28 was also highly expressed during the secondary growth period, whereas increased expression of SsuFLA35 was associated with seed hair tissue. These results provide important clues about the differences in the FLA gene family during the evolution of herbs and woody plants, and suggest that the FLA gene family may play an essential role in regulating the secondary growth of woody plants. It also provides a reference for further studies on the regulation of secondary growth and seed hair development by FLA genes in poplar and willow.

Visible gland constantly traces virus-induced gene silencing in cotton.

A virus-induced gene silencing (VIGS) system was established to induce endogenous target gene silencing by post-transcriptional gene silencing (PTGS), which is a powerful tool for gene function analysis in plants. Compared with stable transgenic plant via Agrobacterium-mediated gene transformation, phenotypes after gene knockdown can be obtained rapidly, effectively, and high-throughput through VIGS system. This approach has been successfully applied to explore unknown gene functions involved in plant growth and development, physiological metabolism, and biotic and abiotic stresses in various plants. In this system, GhCLA1 was used as a general control, however, silencing of this gene leads to leaf albino, wilting, and plant death ultimately. As such, it cannot indicate the efficiency of target gene silencing throughout the whole plant growth period. To address this question, in this study, we developed a novel marker gene, Gossypium PIGMENT GLAND FORMATION GENE (GoPGF), as the control to trace the efficiency of gene silencing in the infected tissues. GoPGF has been proved a key gene in gland forming. Suppression of GoPGF does not affect the normal growth and development of cotton. The number of gland altered related to the expression level of GoPGF gene. So it is a good marker that be used to trace the whole growth stages of plant. Moreover, we further developed a method of friction inoculation to enhance and extend the efficiency of VIGS, which facilitates the analysis of gene function in both the vegetative stage and reproductive stage. This improved VIGS technology will be a powerful tool for the rapid functional identification of unknown genes in genomes.

A Selection of Reliable Reference Genes for Gene Expression Analysis in the Female and Male Flowers of Salix suchowensis.

Salix is a dioecious plant. Research on the molecular regulation mechanism of male and female inflorescence differentiation and development is necessary to analyze sex differentiation in the willow and the underlying mechanisms of unisexual flower development. However, at present, there are no reference genes suitable for stable expression in the process of willow inflorescence development. In this study, Salix suchowensis was used as the research material, nine candidate reference genes (α-TUB1, α-TUB2, ACT, H2A, DnaJ, CDC2, GAPDH, TIP41, ß-TUB) were selected, and qRT-PCR technology was used to detect the expression of each candidate reference gene in female and male flowers at different developmental stages and using five algorithms (geNorm, Normfinder, Delta Ct, BestKeeper, and RefFinder) to comprehensively evaluate the stability of candidate reference genes. The results showed that ACT and DnaJ were stably expressed in all samples and could be used as reference genes. In addition, the reliability of the screening results was further verified via an expression pattern analysis of the CFS gene that encodes flower specific transcription factor in different samples. The stable reference genes selected in this study provide the basis for future research on the expression analysis of functional genes related to the development of male and female flowers of S. suchowensis.

Genome-Wide Comparative Analysis of R2R3 MYB Gene Family in Populus and Salix and Identification of Male Flower Bud Development-Related Genes.

The MYB transcription factor (TF) family is one of the largest plant transcription factor gene family playing vital roles in plant growth and development, including defense, cell differentiation, secondary metabolism, and responses to biotic and abiotic stresses. As a model tree species of woody plants, in recent years, the identification and functional prediction of certain MYB family members in the poplar genome have been reported. However, to date, the characterization of the gene family in the genome of the poplar's sister species willow has not been done, nor are the differences and similarities between the poplar and willow genomes understood. In this study, we conducted the first genome-wide investigation of the R2R3 MYB subfamily in the willow, identifying 216 R2R3 MYB gene members, and combined with the poplar R2R3 MYB genes, performed the first comparative analysis of R2R3 MYB genes between the poplar and willow. We identified 81 and 86 pairs of R2R3 MYB paralogs in the poplar and willow, respectively. There were 17 pairs of tandem repeat genes in the willow, indicating active duplication of willow R2R3 MYB genes. A further 166 pairs of poplar and willow orthologs were identified by collinear and synonymous analysis. The findings support the duplication of R2R3 MYB genes in the ancestral species, with most of the R2R3 MYB genes being retained during the evolutionary process. The phylogenetic trees of the R2R3 MYB genes of 10 different species were drawn. The functions of the poplar and willow R2R3 MYB genes were predicted using reported functional groupings and clustering by OrthoFinder. Identified 5 subgroups in general expanded in woody species, three subgroups were predicted to be related to lignin synthesis, and we further speculate that the other two subgroups also play a role in wood formation. We analyzed the expression patterns of the GAMYB gene of subgroup 18 (S18) related to pollen development in the male flower buds of poplar and willow at different developmental stages by qRT-PCR. The results showed that the GAMYB gene was specifically expressed in the male flower bud from pollen formation to maturity, and that the expression first increased and then decreased. Both the specificity of tissue expression specificity and conservation indicated that GAMYB played an important role in pollen development in both poplar and willow and was an ideal candidate gene for the analysis of male flower development-related functions of the two species.

Evidences for a role of two Y-specific genes in sex determination in Populus deltoides.

Almost all plants in the genus Populus are dioecious (i.e. trees are either male or female), but it is unknown whether dioecy evolved in a common ancestor or independently in different subgenera. Here, we sequence the small peritelomeric X- and Y-linked regions of P. deltoides chromosome XIX. Two genes are present only in the Y-linked region. One is a duplication of a non-Y-linked, female-specifically expressed response regulator, which produces siRNAs that block this gene's expression, repressing femaleness. The other is an LTR/Gypsy transposable element family member, which generates long non-coding RNAs. Overexpression of this gene in A. thaliana promotes androecium development. We also find both genes in the sex-determining region of P. simonii, a different poplar subgenus, which suggests that they are both stable components of poplar sex-determining systems. By contrast, only the duplicated response regulator gene is present in the sex-linked regions of P. davidiana and P. tremula. Therefore, findings in our study suggest dioecy may have evolved independently in different poplar subgenera.

Efficient CRISPR/Cas9-Mediated Gene Editing in an Interspecific Hybrid Poplar With a Highly Heterozygous Genome.

Although the CRISPR/Cas9 system has been widely used for crop breeding, its application for the genetic improvement of trees has been limited, partly because of the outcrossing nature and substantial genomic heterozygosity of trees. Shanxin yang (Populus davidiana × P. bolleana), is a commercially important poplar clone that is widely grown in northern China. An established transformation protocol for this interspecific hybrid enables researchers to simultaneously investigate the efficiency and specificity of the CRISPR/Cas9-mediated manipulation of a highly heterozygous genome. Using the phytoene desaturase gene (PDS) as an example, we revealed that the CRISPR/Cas9 system could efficiently edit the Shanxin yang genome. Two sgRNAs were designed and incorporated into a single binary vector containing the Cas9 expression cassette. Among 62 independent transgenic lines, 85.5% exhibited an exclusively albino phenotype, revealing the total loss of PDS function. The Illumina sequencing results confirmed the targeted mutation of PdbPDS homologs induced by CRISPR/Cas9, and small insertions/deletions were the most common mutations. Biallelic and homozygous knockout mutations were detected at both target sites of the T0 transformants. Off-target activity was detected for sgRNA2 with a frequency of 3.2%. Additionally, the SNP interference of targeting specificity was assessed based on the sequence variation among PdbPDS homologs. A single mismatch at 19- or 10-bp from the PAM was tolerated by the CRISPR/Cas9 system. Therefore, multiple homologous genes were simultaneously edited despite the presence of a mismatch between the sgRNA and the target site. The establishment of a viable CRISPR/Cas9-based strategy for editing the Shanxin yang genome will not only accelerate the breeding process, but may also be relevant for other economically or scientifically important non-model plants species.

Fine mapping of the sex locus in Salix triandra confirms a consistent sex determination mechanism in genus Salix.

Salix triandra belongs to section Amygdalinae in genus Salix, which is in a different section from the willow species in which sex determination has been well studied. Studying sex determination in distantly related willow species will help to clarify whether the sexes of different willows arise through a common sex determination system. For this purpose, we generated an intraspecific full-sib F1 population for S. triandra and constructed high-density genetic linkage maps for the crossing parents using restriction site-associated DNA sequencing and following a two-way pseudo-testcross strategy. With the established maps, the sex locus was positioned in linkage group XV only in the maternal map, and no sex linkage was detected in the paternal map. Consistent with previous findings in other willow species, our study showed that chromosome XV was the incipient sex chromosome and that females were the heterogametic sex in S. triandra. Therefore, sex in this willow species is also determined through a ZW sex determination system. We further performed fine mapping in the vicinity of the sex locus with SSR markers. By comparing the physical and genetic distances for the target interval encompassing the sex determination gene confined by SSRs, severe recombination repression was revealed in the sex determination region in the female map. The recombination rate in the confined interval encompassing the sex locus was approximately eight-fold lower than the genome-wide average. This study provides critical information relevant to sex determination in S. triandra.

The transcription factor MML4_D12 regulates fiber development through interplay with the WD40-repeat protein WDR in cotton.

In planta, a vital regulatory complex, MYB-basic helix-loop-helix (bHLH)-WD40 (MBW), is involved in trichome development and synthesis of anthocyanin and proanthocyanin in Arabidopsis. Usually, WD40 proteins provide a scaffold for protein-protein interaction between MYB and bHLH proteins. Members of subgroup 9 of the R2R3 MYB transcription factors, which includes MYBMIXTA-Like (MML) genes important for plant cell differentiation, are unable to interact with bHLH. In this study, we report that a cotton (Gossypium hirsutum) seed trichome or lint fiber-related GhMML factor, GhMML4_D12, interacts with a diverged WD40 protein (GhWDR) in a process similar to but different from that of the MBW ternary complex involved in Arabidopsis trichome development. Amino acids 250-267 of GhMML4_D12 and the first and third WD40 repeat domains of GhWDR determine their interaction. GhWDR could rescue Arabidopsis ttg1 to its wild type, confirming its orthologous function in trichome development. Our findings shed more light towards understanding the key role of the MML and WD40 families in plants and in the improvement of cotton fiber production.

Identification of Genes Underlying the Resistance to Melampsora larici-populina in an R Gene Supercluster of the Populus deltoides Genome.

Identification of the particular genes in an R genes supercluster underlying resistance to the rust fungus Melampsora larici-populina in poplar genome remains challenging. Based on the de novo assembly of the Populus deltoides genome, all of the detected major genetic loci conferring resistance to M. larici-populina were confined to a 3.5-Mb region on chromosome 19. The transcriptomes of the resistant and susceptible genotypes were sequenced for a timespan from 0 to 168 hours postinoculation. By mapping the differentially expressed genes to the target genomic region, we identified two constitutive expression R genes and one inducible expression R gene that might confer resistance to M. larici-populina. Nucleotide variations were predicted based on the reconstructed haplotypes for each allele of the candidate genes. We also confirmed that salicylic acid was the phytohormone mediating signal transduction pathways, and PR-1 was identified as a key gene inhibiting rust reproduction. Finally, quantitative reverse transcription PCR assay revealed consistent expressions with the RNA-sequencing data for the detected key genes. This study presents an efficient approach for the identification of particular genes underlying phenotype of interest by the combination of genetic mapping, transcriptome profiling, and candidate gene sequences dissection. The identified key genes would be useful for host resistance diagnosis and for molecular breeding of elite poplar cultivars exhibiting resistance to M. larici-populina infection. The detected R genes are also valuable for testing whether the combination of individual R genes can induce durable quantitative resistance.

Analyzing and Characterizing the Chloroplast Genome of Salix wilsonii.

Salix wilsonii is an important ornamental willow tree widely distributed in China. In this study, an integrated circular chloroplast genome was reconstructed for S. wilsonii based on the chloroplast reads screened from the whole-genome sequencing data generated with the PacBio RSII platform. The obtained pseudomolecule was 155,750 bp long and had a typical quadripartite structure, comprising a large single copy region (LSC, 84,638 bp) and a small single copy region (SSC, 16,282 bp) separated by two inverted repeat regions (IR, 27,415 bp). The S. wilsonii chloroplast genome encoded 115 unique genes, including four rRNA genes, 30 tRNA genes, 78 protein-coding genes, and three pseudogenes. Repetitive sequence analysis identified 32 tandem repeats, 22 forward repeats, two reverse repeats, and five palindromic repeats. Additionally, a total of 118 perfect microsatellites were detected, with mononucleotide repeats being the most common (89.83%). By comparing the S. wilsonii chloroplast genome with those of other rosid plant species, significant contractions or expansions were identified at the IR-LSC/SSC borders. Phylogenetic analysis of 17 willow species confirmed that S. wilsonii was most closely related to S. chaenomeloides and revealed the monophyly of the genus Salix. The complete S. wilsonii chloroplast genome provides an additional sequence-based resource for studying the evolution of organelle genomes in woody plants.

Intrinsic checkpoint deficiency during cell cycle re-entry from quiescence.

To maintain tissue homeostasis, cells transition between cell cycle quiescence and proliferation. An essential G1 process is minichromosome maintenance complex (MCM) loading at DNA replication origins to prepare for S phase, known as origin licensing. A p53-dependent origin licensing checkpoint normally ensures sufficient MCM loading before S phase entry. We used quantitative flow cytometry and live cell imaging to compare MCM loading during the long first G1 upon cell cycle entry and the shorter G1 phases in the second and subsequent cycles. We discovered that despite the longer G1 phase, the first G1 after cell cycle re-entry is significantly underlicensed. Consequently, the first S phase cells are hypersensitive to replication stress. This underlicensing results from a combination of slow MCM loading with a severely compromised origin licensing checkpoint. The hypersensitivity to replication stress increases over repeated rounds of quiescence. Thus, underlicensing after cell cycle re-entry from quiescence distinguishes a higher-risk first cell cycle that likely promotes genome instability.

Gossypium barbadense and Gossypium hirsutum genomes provide insights into the origin and evolution of allotetraploid cotton.

Allotetraploid cotton is an economically important natural-fiber-producing crop worldwide. After polyploidization, Gossypium hirsutum L. evolved to produce a higher fiber yield and to better survive harsh environments than Gossypium barbadense, which produces superior-quality fibers. The global genetic and molecular bases for these interspecies divergences were unknown. Here we report high-quality de novo-assembled genomes for these two cultivated allotetraploid species with pronounced improvement in repetitive-DNA-enriched centromeric regions. Whole-genome comparative analyses revealed that species-specific alterations in gene expression, structural variations and expanded gene families were responsible for speciation and the evolutionary history of these species. These findings help to elucidate the evolution of cotton genomes and their domestication history. The information generated not only should enable breeders to improve fiber quality and resilience to ever-changing environmental conditions but also can be translated to other crops for better understanding of their domestication history and use in improvement.

Genetics and evolution of MIXTA genes regulating cotton lint fiber development.

Cotton, with cellulose-enriched mature fibers, is the largest source of natural textiles. Through a map-based cloning strategy, we isolated an industrially important lint fiber development gene (Li3 ) that encodes an MYB-MIXTA-like transcription factor (MML) on chromosome D12 (GhMML4_D12). Virus-induced gene silencing or decreasing the expression of the GhMML4_D12 gene in n2 NSM plants resulted in a significant reduction in epidermal cell prominence and lint fiber production. GhMML4_D12 is arranged in tandem with GhMML3, another MIXTA gene responsible for fuzz fiber development. These two very closely related MIXTA genes direct fiber initiation production in two specialized cell forms: lint and fuzz fibers. They may control the same metabolic pathways in different cell types. The MIXTAs expanded in Malvaceae during their evolution and produced a Malvaceae-specific family that regulates epidermal cell differentiation, different from the gene family that regulates leaf hair trichome development. Cotton has developed a unique transcriptional regulatory network for fiber development. Characterization of target genes regulating fiber production has provided insights into the molecular mechanisms underlying cotton fiber development and has allowed the use of genetic engineering to increase lint yield by inducing more epidermal cells to develop into lint rather than fuzz fibers.

Rapid mapping and cloning of the virescent-1 gene in cotton by bulked segregant analysis-next generation sequencing and virus-induced gene silencing strategies.

Map-based gene cloning is a vital strategy for the identification of the quantitative trait loci or genes underlying important agronomic traits. The conventional map-based cloning method is powerful but generally time-consuming and labor-intensive. In this context, we introduce an improved bulked segregant analysis method in combination with a virus-induced gene silencing (VIGS) strategy for rapid and reliable gene mapping, identification and functional verification. This method was applied to a multiple recessive marker line of upland cotton, Texas 582 (T582), and identified unique genomic positions harboring mutant loci, showing the reliability and efficacy of this method. The v1 locus was further fine-mapped. Only one gene, GhCHLI, which encodes one of the subunits of Mg chelatase, was differentially down-regulated in T582 compared with TM-1. A point mutation occurred in the AAA+ conserved region of GhCHLI and led to an amino acid substitution. Suppression of its expression by VIGS in TM-1 resulted in a yellow blade phenotype that was similar to T582. This integrated approach provides a paradigm for the rapid mapping and identification of the candidate genes underlying the genetic traits in plants with large and complex genomes in the future.

Genomic insights into divergence and dual domestication of cultivated allotetraploid cottons.

BACKGROUND: Cotton has been cultivated and used to make fabrics for at least 7000 years. Two allotetraploid species of great commercial importance, Gossypium hirsutum and Gossypium barbadense, were domesticated after polyploidization and are cultivated worldwide. Although the overall genetic diversity between these two cultivated species has been studied with limited accessions, their population structure and genetic variations remain largely unknown. RESULTS: We resequence the genomes of 147 cotton accessions, including diverse wild relatives, landraces, and modern cultivars, and construct a comprehensive variation map to provide genomic insights into the divergence and dual domestication of these two important cultivated tetraploid cotton species. Phylogenetic analysis shows two divergent groups for G. hirsutum and G. barbadense, suggesting a dual domestication processes in tetraploid cottons. In spite of the strong genetic divergence, a small number of interspecific reciprocal introgression events are found between these species and the introgression pattern is significantly biased towards the gene flow from G. hirsutum into G. barbadense. We identify selective sweeps, some of which are associated with relatively highly expressed genes for fiber development and seed germination. CONCLUSIONS: We report a comprehensive analysis of the evolution and domestication history of allotetraploid cottons based on the whole genomic variation between G. hirsutum and G. barbadense and between wild accessions and modern cultivars. These results provide genomic bases for improving cotton production and for further evolution analysis of polyploid crops.

Insights into Interspecific Hybridization Events in Allotetraploid Cotton Formation from Characterization of a Gene-Regulating Leaf Shape.

The morphology of cotton leaves varies considerably. Phenotypes, including okra, sea-island, super-okra, and broad leaf, are controlled by a multiple allele locus, L2 Okra leaf (L2°) is an incomplete mutation that alters leaf shape by increasing the length of lobes with deeper sinuses. Using a map-based cloning strategy, we cloned the L2 locus gene, which encodes a LATE MERISTEM IDENTITY 1 (LMI1)-like transcription factor (GhOKRA). Silencing GhOKRA leads to a change in phenotype from okra to broad leaf. Overexpression of GhOKRA in Arabidopsis thaliana greatly increases the degree of the leaf lobes and changes the leaf shape. Premature termination of translation in GhOKRA results in the production of broad leaves. The sequences of OKRA from diploid progenitor D-genome species, and wild races and domesticated allotetraploid cottons in Gossypium hirsutum show that a premature termination mutation occurred before and after the formation of tetraploid cotton, respectively. This study provides genomic insights into the two interspecific hybridization events: one produced the present broad leaf and another formed okra leaf phenotype with complete OKRA, that occurred during allotetraploid cotton formation.