MERISTEM-DEFECTIVE regulates the balance between stemness and differentiation in the root meristem through RNA splicing control.

Plants respond to environmental stresses through controlled stem cell maintenance and meristem activity. One level of gene regulation is RNA alternative splicing. However, the mechanistic link between stress, meristem function and RNA splicing is poorly understood. The MERISTEM-DEFECTIVE (MDF) Arabidopsis gene encodes an SR-related family protein, required for meristem function and leaf vascularization, and is the likely orthologue of the human SART1 and yeast Snu66 splicing factors. MDF is required for the correct splicing and expression of key transcripts associated with root meristem function. We identified RSZ33 and ACC1, both known to regulate cell patterning, as splicing targets required for MDF function in the meristem. MDF expression is modulated by osmotic and cold stress, associated with differential splicing and specific isoform accumulation and shuttling between nucleus and cytosol, and acts in part via a splicing target SR34. We propose a model in which MDF controls splicing in the root meristem to promote stemness and to repress stress response, cell differentiation and cell death pathways.

Synergistic interplay of redox homeostasis and polysaccharide synthesis promotes cotton fiber elongation.

Cell polarity is the foundation of cell development and tissue morphogenesis. The investigation of polarized growth provides opportunities to gain profound insights into morphogenesis and tissue functionality in organisms. Currently, there are still many mysteries surrounding the mechanisms that regulate polarized cell growth. Cotton fiber cells serve as an excellent model for studying polarized growth, and provide important clues for unraveling the molecular mechanisms, signaling pathways, and regulatory networks of polarized growth. In this study, we characterized two functional genes, GhMDHAR1AT/DT and GhDHAR2AT/DT with predominant expression during fiber elongation. Loss of function of both genes contributed to a significant increase in fiber length. Transcriptomic data revealed up-regulated expression of antioxidant genes in CRISPR mutant lines, along with delayed expression of secondary wall-related genes and temporally prolonged expression of primary wall-related genes. Experimental evidence demonstrated that the increase in GSH content and glutathione peroxidase (GPX) enzyme activity led to enhanced total antioxidant capacity (T-AOC), resulting in reduced H2O2 levels, which contributed to the extension of fiber elongation stage in CRISPR mutant lines. Moreover, the increased polysaccharide synthesis in CRISPR mutant lines was found to provide an abundant supply of raw materials for fiber cell wall elongation, suggesting that synergistic interplay between redox homeostasis and polysaccharide synthesis in fiber cells may facilitate cell wall remodeling and fiber elongation. This study provides valuable insights for deciphering the mechanisms of cell polarized growth and improving cotton fiber quality.

GhTCE1-GhTCEE1 dimers regulate transcriptional reprogramming during wound-induced callus formation in cotton.

Wounded plant cells can form callus to seal the wound site. Alternatively, wounding can cause adventitious organogenesis or somatic embryogenesis. These distinct developmental pathways require specific cell fate decisions. Here, we identify GhTCE1, a basic helix-loop-helix family transcription factor, and its interacting partners as a central regulatory module of early cell fate transition during in vitro dedifferentiation of cotton (Gossypium hirsutum). RNAi- or CRISPR/Cas9-mediated loss of GhTCE1 function resulted in excessive accumulation of reactive oxygen species (ROS), arrested callus cell elongation, and increased adventitious organogenesis. In contrast, GhTCE1-overexpressing tissues underwent callus cell growth, but organogenesis was repressed. Transcriptome analysis revealed that several pathways depend on proper regulation of GhTCE1 expression, including lipid transfer pathway components, ROS homeostasis, and cell expansion. GhTCE1 bound to the promoters of the target genes GhLTP2 and GhLTP3, activating their expression synergistically, and the heterodimer TCE1-TCEE1 enhances this activity. GhLTP2- and GhLTP3-deficient tissues accumulated ROS and had arrested callus cell elongation, which was restored by ROS scavengers. These results reveal a unique regulatory network involving ROS and lipid transfer proteins, which act as potential ROS scavengers. This network acts as a switch between unorganized callus growth and organized development during in vitro dedifferentiation of cotton cells.

Drought response revealed by chromatin organization variation and transcriptional regulation in cotton.

BACKGROUND: Cotton is a major world cash crop and an important source of natural fiber, oil, and protein. Drought stress is becoming a restrictive factor affecting cotton production. To facilitate the development of drought-tolerant cotton varieties, it is necessary to study the molecular mechanism of drought stress response by exploring key drought-resistant genes and related regulatory factors. RESULTS: In this study, two cotton varieties, ZY007 (drought-sensitive) and ZY168 (drought-tolerant), showing obvious phenotypic differences under drought stress, were selected. A total of 25,898 drought-induced genes were identified, exhibiting significant enrichment in pathways related to plant stress responses. Under drought induction, At subgenome expression bias was observed at the whole-genome level, which may be due to stronger inhibition of Dt subgenome expression. A gene co-expression module that was significantly associated with drought resistance was identified. About 90% of topologically associating domain (TAD) boundaries were stable, and 6613 TAD variation events were identified between the two varieties under drought. We identified 92 genes in ZY007 and 98 in ZY168 related to chromatin 3D structural variation and induced by drought stress. These genes are closely linked to the cotton response to drought stress through canonical hormone-responsive pathways, modulation of kinase and phosphatase activities, facilitation of calcium ion transport, and other related molecular mechanisms. CONCLUSIONS: These results lay a foundation for elucidating the molecular mechanism of the cotton drought response and provide important regulatory locus and gene resources for the future molecular breeding of drought-resistant cotton varieties.

High-resolution sequencing of nine elite upland cotton cultivars uncovers genic variations and breeding improvement targets.

Structural variations (SVs) are critical factors affecting genome evolution and important traits. However, identification results and functional analyses of SVs in upland cotton are rare. Here, based on the genetic relationships, breeding history and cumulative planting area of upland cotton in China, nine predominant cultivars from the past 60 years (1950s-2010s) were selected for long read sequencing to uncover genic variations and breeding improvement targets for this crop. Based on the ZM24 reference genome, 0.88-1.47 × 104 SVs per cultivar were identified, and an SV set was constructed. SVs affected the expression of a large number of genes during fiber elongation, and a transposable element insertion resulted in the glandless phenotype in upland cotton. Six widespread inversions were identified based on nine draft genomes and high-throughput chromosome conformation capture data. Multiple haplotype blocks that were always associated with aggregated SVs were demonstrated to play a pivotal role in the agronomic traits of upland cotton and drove its adaptation to the northern planting region. Exotic introgression was the source of these haplotype blocks and increased the genetic diversity of upland cotton. Our results enrich the genome resources of upland cotton, and the identified SVs will promote genetic and breeding research in cotton.

The elicitor VP2 from Verticillium dahliae triggers defence response in cotton.

Verticillium dahliae is a widespread and destructive soilborne vascular pathogenic fungus that causes serious diseases in dicot plants. Here, comparative transcriptome analysis showed that the number of genes upregulated in defoliating pathotype V991 was significantly higher than in the non-defoliating pathotype 1cd3-2 during the early response of cotton. Combined with analysis of the secretome during the V991-cotton interaction, an elicitor VP2 was identified, which was highly upregulated at the early stage of V991 invasion, but was barely expressed during the 1cd3-2-cotton interaction. Full-length VP2 could induce cell death in several plant species, and which was dependent on NbBAK1 but not on NbSOBIR1 in N. benthamiana. Knock-out of VP2 attenuated the pathogenicity of V991. Furthermore, overexpression of VP2 in cotton enhanced resistance to V. dahliae without causing abnormal plant growth and development. Several genes involved in JA, SA and lignin synthesis were significantly upregulated in VP2-overexpressing cotton. The contents of JA, SA, and lignin were also significantly higher than in the wild-type control. In summary, the identified elicitor VP2, recognized by the receptor in the plant membrane, triggers the cotton immune response and enhances disease resistance.

The complex hexaploid oil-Camellia genome traces back its phylogenomic history and multi-omics analysis of Camellia oil biosynthesis.

Oil-Camellia (Camellia oleifera), belonging to the Theaceae family Camellia, is an important woody edible oil tree species. The Camellia oil in its mature seed kernels, mainly consists of more than 90% unsaturated fatty acids, tea polyphenols, flavonoids, squalene and other active substances, which is one of the best quality edible vegetable oils in the world. However, genetic research and molecular breeding on oil-Camellia are challenging due to its complex genetic background. Here, we successfully report a chromosome-scale genome assembly for a hexaploid oil-Camellia cultivar Changlin40. This assembly contains 8.80 Gb genomic sequences with scaffold N50 of 180.0 Mb and 45 pseudochromosomes comprising 15 homologous groups with three members each, which contain 135 868 genes with an average length of 3936 bp. Referring to the diploid genome, intragenomic and intergenomic comparisons of synteny indicate homologous chromosomal similarity and changes. Moreover, comparative and evolutionary analyses reveal three rounds of whole-genome duplication (WGD) events, as well as the possible diversification of hexaploid Changlin40 with diploid occurred approximately 9.06 million years ago (MYA). Furthermore, through the combination of genomics, transcriptomics and metabolomics approaches, a complex regulatory network was constructed and allows to identify potential key structural genes (SAD, FAD2 and FAD3) and transcription factors (AP2 and C2H2) that regulate the metabolism of Camellia oil, especially for unsaturated fatty acids biosynthesis. Overall, the genomic resource generated from this study has great potential to accelerate the research for the molecular biology and genetic improvement of hexaploid oil-Camellia, as well as to understand polyploid genome evolution.

The chromosome-scale reference genome of mirid bugs (Adelphocoris suturalis) genome provides insights into omnivory, insecticide resistance, and survival adaptation.

BACKGROUND: Adelphocoris suturalis (Hemiptera: Miridae) is a notorious agricultural pest, which causes serious economic losses to a diverse range of agricultural crops around the world. The poor understanding of its genomic characteristics has seriously hindered the establishment of sustainable and environment-friendly agricultural pest management through biotechnology and biological insecticides. RESULTS: Here, we report a chromosome-level assembled genome of A. suturalis by integrating Illumina short reads, PacBio, 10x Chromium, and Hi-C mapping technologies. The resulting 1.29 Gb assembly contains twelve chromosomal pseudomolecules with an N50 of 1.4 and 120.6 Mb for the contigs and scaffolds, respectively, and carries 20,010 protein-coding genes. The considerable size of the A. suturalis genome is predominantly attributed to a high amount of retrotransposons, especially long interspersed nuclear elements (LINEs). Transcriptomic and phylogenetic analyses suggest that A. suturalis-specific candidate effectors, and expansion and expression of gene families associated with omnivory, insecticide resistance and reproductive characteristics, such as digestion, detoxification, chemosensory receptors and long-distance migration likely contribute to its strong environmental adaptability and ability to damage crops. Additionally, 19 highly credible effector candidates were identified and transiently overexpressed in Nicotiana benthamiana for functional assays and potential targeting for insect resistance genetic engineering. CONCLUSIONS: The high-quality genome of A. suturalis provides an important genomic landscape for further investigations into the mechanisms of omnivory, insecticide resistance and survival adaptation, and for the development of integrated management strategies.

N6-methyladenosine RNA modification regulates cotton drought response in a Ca2+ and ABA-dependent manner.

N6 -methyladenosine (m6 A) is the most prevalent internal modification present in mRNAs, and is considered to participate in a range of developmental and biological processes. Drought response is highly regulated at the genomic, transcriptional and post-transcriptional levels. However, the biological function and regulatory mechanism of m6 A modification in the drought stress response is still poorly understood. We generated a transcriptome-wide m6 A map using drought-resistant and drought-sensitive varieties of cotton under different water deficient conditions to uncover patterns of m6 A methylation in cotton response to drought stress. The results reveal that m6 A represents a common modification and exhibit dramatic changes in distribution during drought stress. More 5'UTR m6 A was deposited in the drought-resistant variety and was associated with a positive effect on drought resistance by regulating mRNA abundance. Interestingly, we observed that increased m6 A abundance was associated with increased mRNA abundance under drought, contributing to drought resistance, and vice versa. The demethylase GhALKBH10B was found to decrease m6 A levels, facilitating the mRNA decay of ABA signal-related genes (GhZEP, GhNCED4 and GhPP2CA) and Ca2+ signal-related genes (GhECA1, GhCNGC4, GhANN1 and GhCML13), and mutation of GhALKBH10B enhanced drought resistance at seedling stage in cotton. Virus-induced gene silencing (VIGS) of two Ca2+ -related genes, GhECA1 and GhCNGC4, reduced drought resistance with the decreased m6 A enrichment on silenced genes in cotton. Collectively, we reveal a novel mechanism of post-transcriptional modification involved in affecting drought response in cotton, by mediating m6 A methylation on targeted transcripts in the ABA and Ca2+ signalling transduction pathways.

GhWRKY41 forms a positive feedback regulation loop and increases cotton defence response against Verticillium dahliae by regulating phenylpropanoid metabolism.

Despite the established significance of WRKY proteins and phenylpropanoid metabolism in plant immunity, how WRKY proteins modulate aspects of the phenylpropanoid pathway remains undetermined. To understand better the role of WRKY proteins in plant defence, we identified a cotton (Gossypium hirsutum) protein, GhWRKY41, that is, universally and rapidly induced in three disease-resistant cotton cultivars following inoculation with the plant pathogenic fungus, Verticillium dahliae. We show that overexpression of GhWRKY41 in transgenic cotton and Arabidopsis enhances resistance to V. dahliae, while knock-down increases cotton more susceptibility to the fungus. GhWRKY41 physically interacts with itself and directly activates its own transcription. A genome-wide chromatin immunoprecipitation and high-throughput sequencing (ChIP-seq), in combination with RNA sequencing (RNA-seq) analyses, revealed that 43.1% of GhWRKY41-binding genes were up-regulated in cotton upon inoculation with V. dahliae, including several phenylpropanoid metabolism master switches, receptor kinases, and disease resistance-related proteins. We also show that GhWRKY41 homodimer directly activates the expression of GhC4H and Gh4CL, thereby modulating the accumulation of lignin and flavonoids. This finding expands our understanding of WRKY-WRKY protein interactions and provides important insights into the regulation of the phenylpropanoid pathway in plant immune responses by a WRKY protein.

The complex genome and adaptive evolution of polyploid Chinese pepper (Zanthoxylum armatum and Zanthoxylum bungeanum).

Zanthoxylum armatum and Zanthoxylum bungeanum, known as 'Chinese pepper', are distinguished by their extraordinary complex genomes, phenotypic innovation of adaptive evolution and species-special metabolites. Here, we report reference-grade genomes of Z. armatum and Z. bungeanum. Using high coverage sequence data and comprehensive assembly strategies, we derived 66 pseudochromosomes comprising 33 homologous phased groups of two subgenomes, including autotetraploid Z. armatum. The genomic rearrangements and two whole-genome duplications created large (~4.5 Gb) complex genomes with a high ratio of repetitive sequences (>82%) and high chromosome number (2n = 4x = 132). Further analysis of the high-quality genomes shed lights on the genomic basis of involutional reproduction, allomones biosynthesis and adaptive evolution in Chinese pepper, revealing a high consistent relationship between genomic evolution, environmental factors and phenotypic innovation. Our study provides genomic resources and new insights for investigating diversification and phenotypic innovation in Chinese pepper, with broader implications for the protection of plants under severe environmental changes.

Development of an efficient and precise adenine base editor (ABE) with expanded target range in allotetraploid cotton (Gossypium hirsutum).

BACKGROUND: Base editors (BEs) display diverse applications in a variety of plant species such as Arabidopsis, rice, wheat, maize, soybean, and cotton, where they have been used to mediate precise base pair conversions without the collateral generation of undesirable double-stranded breaks (DSB). Studies of single-nucleotide polymorphisms (SNPs) underpinning plant traits are still challenging, particularly in polyploidy species where such SNPs are present in multiple copies, and simultaneous modification of all alleles would be required for functional analysis. Allotetraploid cotton has a number of homoeologous gene pairs located in the A and D sub-genomes with considerable SNPs, and it is desirable to develop adenine base editors (ABEs) for efficient and precise A-to-G single-base editing without DSB in such complex genome. RESULTS: We established various ABE vectors based on different engineered adenosine deaminase (TadA) proteins fused to Cas9 variants (dCas9, nCas9), enabling efficient A to G editing up to 64% efficiency on-target sites of the allotetraploid cotton genome. Comprehensive analysis showed that GhABE7.10n exhibited the highest editing efficiency, with the main editing sites specifically located at the position A5 (counting the PAM as positions 21-23). Furthermore, DNA and RNA off-target analysis of cotton plants edited with GhABE7.10n and GhABE7.10d by whole genome and whole-transcriptome sequencing revealed no DNA off-target mutations, while very low-level RNA off-target mutations were detected. A new base editor, namely GhABE7.10dCpf1 (7.10TadA + dCpf1), that recognizes a T-rich PAM, was developed for the first time. Targeted A-to-G substitutions generated a single amino acid change in the cotton phosphatidyl ethanolamine-binding protein (GhPEBP), leading to a compact cotton plant architecture, an ideotype for mechanized harvesting of modern cotton production. CONCLUSIONS: Our data illustrate the robustness of adenine base editing in plant species with complex genomes, which provides efficient and precise toolkit for cotton functional genomics and precise molecular breeding.

Single-cell RNA-seq reveals fate determination control of an individual fibre cell initiation in cotton (Gossypium hirsutum).

Cotton fibre is a unicellular seed trichome, and lint fibre initials per seed as a factor determines fibre yield. However, the mechanisms controlling fibre initiation from ovule epidermis are not understood well enough. Here, with single-cell RNA sequencing (scRNA-seq), a total of 14 535 cells were identified from cotton ovule outer integument of Xu142_LF line at four developmental stages (1.5, 1, 0.5 days before anthesis and the day of anthesis). Three major cell types, fibre, non-fibre epidermis and outer pigment layer were identified and then verified by RNA in situ hybridization. A comparative analysis on scRNA-seq data between Xu142 and its fibreless mutant Xu142 fl further confirmed fibre cluster definition. The developmental trajectory of fibre cell was reconstructed, and fibre cell was identified differentiated at 1 day before anthesis. Gene regulatory networks at four stages revealed the spatiotemporal pattern of core transcription factors, and MYB25-like and HOX3 were demonstrated played key roles as commanders in fibre differentiation and tip-biased diffuse growth respectively. A model for early development of a single fibre cell was proposed here, which sheds light on further deciphering mechanism of plant trichome and the improvement of cotton fibre yield.

Degradation of de-esterified pctin/homogalacturonan by the polygalacturonase GhNSP is necessary for pollen exine formation and male fertility in cotton.

The pollen wall exine provides a protective layer for the male gametophyte and is largely composed of sporopollenin, which comprises fatty acid derivatives and phenolics. However, the biochemical nature of the external exine is poorly understood. Here, we show that the male sterile line 1355A of cotton mutated in NO SPINE POLLEN (GhNSP) leads to defective exine formation. The GhNSP locus was identified through map-based cloning and confirmed by genetic analysis (co-segregation test and allele prediction using the CRISPR/Cas9 system). In situ hybridization showed that GhNSP is highly expressed in tapetum. GhNSP encodes a polygalacturonase protein homologous to AtQRT3, which suggests a function for polygalacturonase in pollen exine formation. These results indicate that GhNSP is functionally different from AtQRT3, the latter has the function of microspore separation. Biochemical analysis showed that the percentage of de-esterified pectin was significantly increased in the 1355A anthers at developmental stage 8. Furthermore, immunofluorescence studies using antibodies to the de-esterified and esterified homogalacturonan (JIM5 and JIM7) showed that the Ghnsp mutant exhibits abundant of de-esterified homogalacturonan in the tapetum and exine, coupled with defective exine formation. The characterization of GhNSP provides new understanding of the role of polygalacturonase and de-esterified homogalacturonan in pollen exine formation.

Gibberellin signaling mediates lateral root inhibition in response to K+-deprivation.

The potassium ion (K+) is vital for plant growth and development, and K+-deprivation leads to reduced crop yields. Here we describe phenotypic, transcriptomic, and mutant analyses to investigate the signaling mechanisms mediating root architectural changes in Arabidopsis (Arabidopsis thaliana) Columbia. We showed effects on root architecture are mediated through a reduction in cell division in the lateral root (LR) meristems, the rate of LR initiation is reduced but LR density is unaffected, and primary root growth is reduced only slightly. This was primarily regulated through gibberellic acid (GA) signaling, which leads to the accumulation of growth-inhibitory DELLA proteins. The short LR phenotype was rescued by exogenous application of GA but not of auxin or by the inhibition of ethylene signaling. RNA-seq analysis showed upregulation by K+-deprivation of the transcription factors JUNGBRUNNEN1 (JUB1) and the C-repeat-binding factor (CBF)/dehydration-responsive element-binding factor 1 regulon, which are known to regulate GA signaling and levels that regulate DELLAs. Transgenic overexpression of JUB1 and CBF1 enhanced responses to K+ stress. Attenuation of the reduced LR growth response occurred in mutants of the CBF1 target gene SFR6, implicating a role for JUB1, CBF1, and SFR6 in the regulation of LR growth in response to K+-deprivation via DELLAs. We propose this represents a mechanism to limit horizontal root growth in conditions where K+ is available deeper in the soil.

Silencing of a LIM gene in cotton exhibits enhanced resistance against Apolygus lucorum.

Plant bugs (Miridae species) have become major agricultural pests that cause increasing and severe economic damage. Plant-mediated RNA interference (RNAi) is emerging as an eco-friendly, efficient, and reliable strategy for pest management. In this study, we isolated and characterized a lethal gene of Apolygus lucorum and named it Apolygus lucorum LIM (AlLIM), which produced A. lucorum mortality rates ranging from 38% to 81%. Downregulation of the AlLIM gene expression in A. lucorum by injection of a double-stranded RNA (dsRNA) led to muscle structural disorganization that resulted in metamorphosis deficiency and increased mortality. Then we constructed a plant expression vector that enabled transgenic cotton to highly and stably express dsRNA of AlLIM (dsAlLIM) by Agrobacterium-mediated genetic transformation. In the field bioassay, dsAlLIM transgenic cotton was protected from A. lucorum damage with high efficiency, with almost no detectable yield loss. Therefore, our study successfully provides a promising genetically modified strategy to overpower A. lucorum attack.

Epidermal expression of a sterol biosynthesis gene regulates root growth by a non-cell-autonomous mechanism in Arabidopsis.

The epidermis is hypothesized to play a signalling role during plant development. One class of mutants showing defects in signal transduction and radial patterning are those in sterol biosynthesis. The expectation is that living cells require sterols, but it is not clear that all cell types express sterol biosynthesis genes. The HYDRA1 (HYD1) gene of Arabidopsis encodes sterol Δ8-Δ7 isomerase, and although hyd1 seedlings are defective in radial patterning across several tissues, we show that the HYD1 gene is expressed most strongly in the root epidermis. Transgenic activation of HYD1 transcription in the epidermis of hyd1 null mutants reveals a major role in root patterning and growth. HYD1 expression in the vascular tissues and root meristem, though not endodermis or pericycle, also leads to some phenotypic rescue. Phenotypic rescue is associated with rescued patterning of the PIN1 and PIN2 auxin efflux carriers. The importance of the epidermis in controlling root growth and development is proposed to be, in part, due to its role as a site for sterol biosynthesis, and auxin is a candidate for the non-cell-autonomous signal.

Plant 3D genomics: the exploration and application of chromatin organization.

Eukaryotic genomes are highly folded for packing into higher-order chromatin structures in the nucleus. With the emergence of state-of-the-art chromosome conformation capture methods and microscopic imaging techniques, the spatial organization of chromatin and its functional implications have been interrogated. Our knowledge of 3D chromatin organization in plants has improved dramatically in the past few years, building on the early advances in animal systems. Here, we review recent advances in 3D genome mapping approaches, our understanding of the sophisticated organization of spatial structures, and the application of 3D genomic principles in plants. We also discuss directions for future developments in 3D genomics in plants.

Root growth responses to mechanical impedance are regulated by a network of ROS, ethylene and auxin signalling in Arabidopsis.

The growth and development of root systems is influenced by mechanical properties of the substrate in which the plants grow. Mechanical impedance, such as by compacted soil, can reduce root elongation and limit crop productivity. To understand better the mechanisms involved in plant root responses to mechanical impedance stress, we investigated changes in the root transcriptome and hormone signalling responses of Arabidopsis to artificial root barrier systems in vitro. We demonstrate that upon encountering a barrier, reduced Arabidopsis root growth and a characteristic 'step-like' growth pattern is due to a reduction in cell elongation associated with changes in signalling gene expression. Data from RNA-sequencing combined with reporter line and mutant studies identified essential roles for reactive oxygen species, ethylene and auxin signalling during the barrier response. We propose a model in which early responses to mechanical impedance include reactive oxygen signalling integrated with ethylene and auxin responses to mediate root growth changes. Inhibition of ethylene responses allows improved growth in response to root impedance, an observation that may inform future crop breeding programmes.

The Arabidopsis R-SNARE VAMP714 is essential for polarisation of PIN proteins and auxin responses.

The plant hormone auxin and its directional intercellular transport play a major role in diverse aspects of plant growth and development. The establishment of auxin gradients requires the asymmetric distribution of members of the auxin efflux carrier PIN-FORMED (PIN) protein family to the plasma membrane. An endocytic pathway regulates the recycling of PIN proteins between the plasma membrane and endosomes, providing a mechanism for dynamic localisation. N-Ethylmaleimide-sensitive factor adaptor protein receptors (SNAP receptors, SNAREs) mediate fusion between vesicles and target membranes and are classed as Q- or R-SNAREs based on their sequence. We analysed gain- and loss-of-function mutants, dominant-negative transgenics and localisation of the Arabidopsis R-SNARE VAMP714 protein to understand its function. We demonstrate that VAMP714 is essential for the insertion of PINs into the plasma membrane, for polar auxin transport, root gravitropism and morphogenesis. VAMP714 gene expression is upregulated by auxin, and the VAMP714 protein co-localises with endoplasmic reticulum and Golgi vesicles and with PIN proteins at the plasma membrane. It is proposed that VAMP714 mediates the delivery of PIN-carrying vesicles to the plasma membrane, and that this forms part of a positive regulatory loop in which auxin activates a VAMP714-dependent PIN/auxin transport system to control development.