The maize single-nucleus transcriptome comprehensively describes signaling networks governing movement and development of grass stomata.

The unique morphology of grass stomata enables rapid responses to environmental changes. Deciphering the basis for these responses is critical for improving food security. We have developed a planta platform of single-nucleus RNA-sequencing by combined fluorescence-activated nuclei flow sorting, and used it to identify cell types in mature and developing stomata from 33,098 nuclei of the maize epidermis-enriched tissues. Guard cells (GCs) and subsidiary cells (SCs) displayed differential expression of genes, besides those encoding transporters, involved in the abscisic acid, CO2, Ca2+, starch metabolism, and blue light signaling pathways, implicating coordinated signal integration in speedy stomatal responses, and of genes affecting cell wall plasticity, implying a more sophisticated relationship between GCs and SCs in stomatal development and dumbbell-shaped guard cell formation. The trajectory of stomatal development identified in young tissues, and by comparison to the bulk RNA-seq data of the MUTE defective mutant in stomatal development, confirmed known features, and shed light on key participants in stomatal development. Our study provides a valuable, comprehensive, and fundamental foundation for further insights into grass stomatal function.

Genome-wide characterization of the WAK gene family and expression analysis under plant hormone treatment in cotton.

BACKGROUND: Wall-associated kinases (WAK), one of the receptor-like kinases (RLK), function directly in the connection and communication between the plant cell wall and the cytoplasm. WAK genes are highly conserved and have been identified in plants, such as rice, but there is little research on the WAK gene family in cotton. RESULTS: In the present study, we identified 29 GhWAK genes in Gossypium hirsutum. Phylogenetic analysis showed that cotton WAK proteins can be divided into five clades. The results of synteny and Ka/Ks analysis showed that the GhWAK genes mainly originated from whole genome duplication (WGD) and were then mainly under purifying selection. Transcriptome data and real-time PCR showed that 97% of GhWAK genes highly expressed in cotton fibers and ovules. ß-glucuronidase (GUS) staining assays showed that GhWAK5 and GhWAK16 expressed in Arabidopsis leaf trichomes. Fourteen GhWAK genes were found to possess putative gibberellin (GA) response elements in the promoter regions, 13 of which were significantly induced by GA treatment. Ten GhWAK genes contained auxin (IAA) response elements and the expression level of nine GhWAKs significantly increased under auxin treatment. CONCLUSIONS: We provide a preliminary analysis of the WAK gene family in G. hirsutum, which sheds light on the potantial roles of GhWAK genes in cotton fiber cell development. Our data also provides a useful resource for future studies on the functional roles of GhWAK genes.

GhARF16-1 modulates leaf development by transcriptionally regulating the GhKNOX2-1 gene in cotton.

The leaf is a crucial organ evolved with remarkable morphological diversity to maximize plant photosynthesis. The leaf shape is a key trait that affects photosynthesis, flowering rates, disease resistance and yield. Although many genes regulating leaf development have been identified in the past years, the precise regulatory architecture underlying the generation of diverse leaf shapes remains to be elucidated. We used cotton as a reference model to probe the genetic framework underlying divergent leaf forms. Comparative transcriptome analysis revealed that the GhARF16-1 and GhKNOX2-1 genes might be potential regulators of leaf shape. We functionally characterized the auxin-responsive factor ARF16-1 acting upstream of GhKNOX2-1 to determine leaf morphology in cotton. The transcription of GhARF16-1 was significantly higher in lobed-leaved cotton than in smooth-leaved cotton. Furthermore, the overexpression of GhARF16-1 led to the up-regulation of GhKNOX2-1 and resulted in more and deeper serrations in cotton leaves, similar to the leaf shape of cotton plants overexpressing GhKNOX2-1. We found that GhARF16-1 specifically bound to the promoter of GhKNOX2-1 to induce its expression. The heterologous expression of GhARF16-1 and GhKNOX2-1 in Arabidopsis led to lobed and curly leaves, and a genetic analysis revealed that GhKNOX2-1 is epistatic to GhARF16-1 in Arabidopsis, suggesting that the GhARF16-1 and GhKNOX2-1 interaction paradigm also functions to regulate leaf shape in Arabidopsis. To our knowledge, our results uncover a novel mechanism by which auxin, through the key component ARF16-1 and its downstream-activated gene KNOX2-1, determines leaf morphology in eudicots.

Mutations in a subfamily of abscisic acid receptor genes promote rice growth and productivity.

Abscisic acid (ABA) is a key phytohormone that controls plant growth and stress responses. It is sensed by the pyrabactin resistance 1 (PYR1)/PYR1-like (PYL)/regulatory components of the ABA receptor (RCAR) family of proteins. Here, we utilized CRISPR/Cas9 technology to edit group I (PYL1-PYL6 and PYL12) and group II (PYL7-PYL11 and PYL13) PYL genes in rice. Characterization of the combinatorial mutants suggested that genes in group I have more important roles in stomatal movement, seed dormancy, and growth regulation than those in group II. Among all of the single pyl mutants, only pyl1 and pyl12 exhibited significant defects in seed dormancy. Interestingly, high-order group I mutants, but not any group II mutants, displayed enhanced growth. Among group I mutants, pyl1/4/6 exhibited the best growth and improved grain productivity in natural paddy field conditions, while maintaining nearly normal seed dormancy. Our results suggest that a subfamily of rice PYLs has evolved to have particularly important roles in regulating plant growth and reveal a genetic strategy to improve rice productivity.

The plasma-membrane polyamine transporter PUT3 is regulated by the Na+ /H+ antiporter SOS1 and protein kinase SOS2.

In Arabidopsis, the plasma membrane transporter PUT3 is important to maintain the cellular homeostasis of polyamines and plays a role in stabilizing mRNAs of some heat-inducible genes. The plasma membrane Na+ /H+ transporter SOS1 and the protein kinase SOS2 are two salt-tolerance determinants crucial for maintaining intracellular Na+ and K+ homeostasis. Here, we report that PUT3 genetically and physically interacts with SOS1 and SOS2, and these interactions modulate PUT3 transport activity. Overexpression of PUT3 (PUT3OE) results in hypersensitivity of the transgenic plants to polyamine and paraquat. The hypersensitivity of PUT3OE is inhibited by the sos1 and sos2 mutations, which indicates that SOS1 and SOS2 are required for PUT3 transport activity. A protein interaction assay revealed that PUT3 physically interacts with SOS1 and SOS2 in yeast and plant cells. SOS2 phosphorylates PUT3 both in vitro and in vivo. SOS1 and SOS2 synergistically activate the polyamine transport activity of PUT3, and PUT3 also modulates SOS1 activity by activating SOS2 in yeast cells. Overall, our findings suggest that both plasma-membrane proteins PUT3 and SOS1 could form a complex with the protein kinase SOS2 in response to stress conditions and modulate the transport activity of each other through protein interactions and phosphorylation.

Genome-wide identification, phylogeny, and expression analysis of pectin methylesterases reveal their major role in cotton fiber development.

BACKGROUND: Pectin methylesterase (PME, EC 3.1.1.11) is a hydrolytic enzyme that utilizes pectin as substrates, and plays a significant role in regulating pectin reconstruction thereby regulating plant growth. Pectin is one of the important components of the plant cell wall, which forms the main structural material of cotton fiber. In this research, cotton genome information was used to identify PMEs. RESULTS: We identified 80 (GaPME01-GaPME80) PME genes from diploid G. arboreum (A genome), 78 (GrPME01-GrPME78) PME genes from G. raimondii (D genome), and 135 (GhPME001-GhPME135) PME genes from tetraploid cotton G. hirsutum (AD genome). We further analyzed their gene structure, conserved domain, gene expression, and systematic evolution to lay the foundation for deeper research on the function of PMEs. Phylogenetic data indicated that members from the same species demonstrated relatively high sequence identities and genetic similarities. Analysis of gene structures showed that most of the PMEs genes had 2-3 exons, with a few having a variable number of exons from 4 to 6. There are nearly no differences in the gene structure of PMEs among the three (two diploid and one tetraploid) cotton species. Selective pressure analysis showed that the Ka/Ks value for each of the three cotton species PME families was less than one. CONCLUSION: Conserved domain analysis showed that PMEs members had a relatively conserved C-terminal pectinesterase domain (PME) while the N-terminus was less conserved. Moreover, some of the family members contained a pectin methylesterase inhibitor (PMEI) domain. The Ka/Ks ratios suggested that the duplicated PMEs underwent purifying selection after the duplication events. This study provided an important basis for further research on the functions of cotton PMEs. Results from qRT-PCR indicated that the expression level of different PMEs at various fiber developmental stages was different. Moreover, some of the PMEs showed fiber predominant expression in secondary wall thickening indicating tissue-specific expression patterns.

Fine mapping and candidate gene analysis of the dominant glandless gene Gl 2 (e) in cotton (Gossypium spp.).

KEY MESSAGE: Dominant glandless gene Gl 2 (e) was fine-mapped to a 15 kb region containing one candidate gene encoding an MYC transcription factor, sequence and expression level of the gene were analyzed. Cottonseed product is an excellent source of oil and protein. However, this nutrition source is greatly limited in utilization by the toxic gossypol in pigment glands. It is reported that the Gl 2 (e) gene could effectively inhibit the formation of the pigment glands. Here, three F2 populations were constructed using two pairs of near isogenic lines (NILs), which differ nearly only by the gland trait, for fine mapping of Gl 2 (e) . DNA markers were identified from recently developed cotton genome sequence. The Gl 2 (e) gene was located within a 15-kb genomic interval between two markers CS2 and CS4 on chromosome 12. Only one gene was identified in the genomic interval as the candidate for Gl 2 (e) which encodes a family member of MYC transcription factor with 475-amino acids. Unexpectedly, the results of expression analysis indicated that the MYC gene expresses in glanded lines while almost does not express in glandless lines. These results suggest that the MYC gene probably serves as a vital positive regulator in the organogenesis pathway of pigment gland, and low expression of this gene will not launch the downstream pathway of pigment gland formation. This is the first pigment gland-related gene identification in cotton and will facilitate the research on glandless trait, cotton MYC proteins and low-gossypol cotton breeding.

Development of chromosome-specific markers with high polymorphism for allotetraploid cotton based on genome-wide characterization of simple sequence repeats in diploid cottons (Gossypium arboreum L. and Gossypium raimondii Ulbrich).

BACKGROUND: Tetraploid cotton contains two sets of homologous chromosomes, the At- and Dt-subgenomes. Consequently, many markers in cotton were mapped to multiple positions during linkage genetic map construction, posing a challenge to anchoring linkage groups and mapping economically-important genes to particular chromosomes. Chromosome-specific markers could solve this problem. Recently, the genomes of two diploid species were sequenced whose progenitors were putative contributors of the At- and Dt-subgenomes to tetraploid cotton. These sequences provide a powerful tool for developing chromosome-specific markers given the high level of synteny among tetraploid and diploid cotton genomes. In this study, simple sequence repeats (SSRs) on each chromosome in the two diploid genomes were characterized. Chromosome-specific SSRs were developed by comparative analysis and proved to distinguish chromosomes. RESULTS: A total of 200,744 and 142,409 SSRs were detected on the 13 chromosomes of Gossypium arboreum L. and Gossypium raimondii Ulbrich, respectively. Chromosome-specific SSRs were obtained by comparing SSR flanking sequences from each chromosome with those from the other 25 chromosomes. The average was 7,996 per chromosome. To confirm their chromosome specificity, these SSRs were used to distinguish two homologous chromosomes in tetraploid cotton through linkage group construction. The chromosome-specific SSRs and previously-reported chromosome markers were grouped together, and no marker mapped to another homologous chromosome, proving that the chromosome-specific SSRs were unique and could distinguish homologous chromosomes in tetraploid cotton. Because longer dinucleotide AT-rich repeats were the most polymorphic in previous reports, the SSRs on each chromosome were sorted by motif type and repeat length for convenient selection. The primer sequences of all chromosome-specific SSRs were also made publicly available. CONCLUSION: Chromosome-specific SSRs are efficient tools for chromosome identification by anchoring linkage groups to particular chromosomes during genetic mapping and are especially useful in mapping of qualitative-trait genes or quantitative trait loci with just a few markers. The SSRs reported here will facilitate a number of genetic and genomic studies in cotton, including construction of high-density genetic maps, positional gene cloning, fingerprinting, and genetic diversity and comparative evolutionary analyses among Gossypium species.

[Recent progress in gene mapping through high-throughput sequencing technology and forward genetic approaches].

Traditional gene mapping using forward genetic approaches is conducted primarily through construction of a genetic linkage map, the process of which is tedious and time-consuming, and often results in low accuracy of mapping and large mapping intervals. With the rapid development of high-throughput sequencing technology and decreasing cost of sequencing, a variety of simple and quick methods of gene mapping through sequencing have been developed, including direct sequencing of the mutant genome, sequencing of selective mutant DNA pooling, genetic map construction through sequencing of individuals in population, as well as sequencing of transcriptome and partial genome. These methods can be used to identify mutations at the nucleotide level and has been applied in complex genetic background. Recent reports have shown that sequencing mapping could be even done without the reference of genome sequence, hybridization, and genetic linkage information, which made it possible to perform forward genetic study in many non-model species. In this review, we summarized these new technologies and their application in gene mapping.

The role of WRKY transcription factors in plant abiotic stresses.

The WRKY gene family has been suggested to play important roles in the regulation of transcriptional reprogramming associated with plant stress responses. Modification of the expression patterns of WRKY genes and/or changes in their activity contribute to the elaboration of various signaling pathways and regulatory networks. Furthermore, a single WRKY gene often responds to several stress factors, and then their proteins may participate in the regulation of several seemingly disparate processes as negative or positive regulators. WRKY proteins also function via protein-protein interaction and autoregulation or cross-regulation is extensively recorded among WRKY genes, which help us understand the complex mechanisms of signaling and transcriptional reprogramming controlled by WRKY proteins. Here, we review recent progress made in starting to reveal the role of WRKY transcription factors in plant abiotic stresses. This article is part of a Special Issue entitled: Plant gene regulation in response to abiotic stress.

Analyses of the NAC transcription factor gene family in Gossypium raimondii Ulbr.: chromosomal location, structure, phylogeny, and expression patterns.

NAC domain proteins are plant-specific transcription factors known to play diverse roles in various plant developmental processes. In the present study, we performed the first comprehensive study of the NAC gene family in Gossypium raimondii Ulbr., incorporating phylogenetic, chromosomal location, gene structure, conserved motif, and expression profiling analyses. We identified 145 NAC transcription factor (NAC-TF) genes that were phylogenetically clustered into 18 distinct subfamilies. Of these, 127 NAC-TF genes were distributed across the 13 chromosomes, 80 (55%) were preferentially retained duplicates located in both duplicated regions and six were located in triplicated chromosomal regions. The majority of NAC-TF genes showed temporal-, spatial-, and tissue-specific expression patterns based on transcriptomic and qRT-PCR analyses. However, the expression patterns of several duplicate genes were partially redundant, suggesting the occurrence of sub-functionalization during their evolution. Based on their genomic organization, we concluded that genomic duplications contributed significantly to the expansion of the NAC-TF gene family in G. raimondii. Comprehensive analysis of their expression profiles could provide novel insights into the functional divergence among members of the NAC gene family in G. raimondii.

Genome-wide analysis of the Sus gene family in cotton.

Sucrose synthase (Sus) is a key enzyme in plant sucrose metabolism. In cotton, Sus (EC 2.4.1.13) is the main enzyme that degrades sucrose imported into cotton fibers from the phloem of the seed coat. This study demonstrated that the genomes of Gossypium arboreum L., G. raimondii Ulbr., and G. hirsutum L., contained 8, 8, and 15 Sus genes, respectively. Their structural organizations, phylogenetic relationships, and expression profiles were characterized. Comparisons of genomic and coding sequences identified multiple introns, the number and positions of which were highly conserved between diploid and allotetraploid cotton species. Most of the phylogenetic clades contained sequences from all three species, suggesting that the Sus genes of tetraploid G. hirsutum derived from those of its diploid ancestors. One Sus group (Sus I) underwent expansion during cotton evolution. Expression analyses indicated that most Sus genes were differentially expressed in various tissues and had development-dependent expression profiles in cotton fiber cells. Members of the same orthologous group had very similar expression patterns in all three species. These results provide new insights into the evolution of the cotton Sus gene family, and insight into its members' physiological functions during fiber growth and development.

Introgressing the Aegilops tauschii genome into wheat as a basis for cereal improvement.

Increasing crop production is necessary to feed the world's expanding population, and crop breeders often utilize genetic variations to improve crop yield and quality. However, the narrow diversity of the wheat D genome seriously restricts its selective breeding. A practical solution is to exploit the genomic variations of Aegilops tauschii via introgression. Here, we established a rapid introgression platform for transferring the overall genetic variations of A. tauschii to elite wheats, thereby enriching the wheat germplasm pool. To accelerate the process, we assembled four new reference genomes, resequenced 278 accessions of A. tauschii and constructed the variation landscape of this wheat progenitor species. Genome comparisons highlighted diverse functional genes or novel haplotypes with potential applications in wheat improvement. We constructed the core germplasm of A. tauschii, including 85 accessions covering more than 99% of the species' overall genetic variations. This was crossed with elite wheat cultivars to generate an A. tauschii-wheat synthetic octoploid wheat (A-WSOW) pool. Laboratory and field analysis with two examples of the introgression lines confirmed its great potential for wheat breeding. Our high-quality reference genomes, genomic variation landscape of A. tauschii and the A-WSOW pool provide valuable resources to facilitate gene discovery and breeding in wheat.

Male gametophyte-specific WRKY34 transcription factor mediates cold sensitivity of mature pollen in Arabidopsis.

Mature pollen is very sensitive to cold stress in chilling-sensitive plants. Plant WRKY DNA-binding transcription factors are key regulators in plant responses to abiotic and biotic stresses. Previous studies have suggested that WRKY34 (At4g26440) gene might be involved in pollen viability, although the mechanism involved is unclear. In this study, it is shown that cold treatment increased WRKY34 expression in the wild type, and promoter-GUS analysis revealed that WRKY34 expression is pollen-specific. Enhanced green fluorescent protein-tagged WRKY34 was localized in the nuclei. Pollen harbouring the wrky34 allele showed higher viability than pollen with the WRKY34 allele after cold treatment. Further functional analysis indicated that the WRKY34 transcription factor was involved in pollen development regulated by the pollen-specific MIKC* class of MADS-domain transcription factors under cold stress, and cold-insensitivity of mature wrky34 pollen might be partly attributable to the enhanced expression of transcriptional activator CBFs in the mutants. Thus, the WRKY34 transcription factor negatively mediated cold sensitivity of mature Arabidopsis pollen and might be involved in the CBF signal cascade in mature pollen.

Genome-wide comparative analysis of RNA-binding Glycine-rich protein family genes between Gossypium arboreum and Gossypium raimondii.

RB-GRP (RNA-binding Glycine-rich protein gene) family belongs to the fourth subfamily of the GRP (Glycine-rich protein gene) superfamily, which plays a great role in plant growth and development, as well as in abiotic stresses response, while it has not been identified in cotton. Here, we identified 37 and 32 RB-GRPs from two cotton species (Gossypium arboreum and Gossypium raimondii, respectively), which were divided into four distinct subfamilies based on the presence of additional motifs and the arrangement of the glycine repeats. The distribution of RB-GRPs was nonrandom and uneven among the chromosomes both in two cotton species. The expansion of RB-GRP gene family between two cultivars was mainly attributed to segmental and tandem duplication events indicated by synteny analysis, and the tandem duplicated genes were mapped into homologous collinear blocks, indicated that they shared a common ancestral gene in both species. Furthermore, most RB-GRPs in two cotton species undergone stronger negative selective pressure by evolutionary analysis of RB-GRP orthologous genes. Meanwhile, RB-GRPs participated in different abiotic stresses (Abscisic acid, salt and Polyethylene glycol) responses and tissues at different developmental stages between two cotton species were showed by gene expression analysis. This research would provide insight into the evolution and function of the RB-GRPs in Gossypium species.

The genome of broomcorn millet.

Broomcorn millet (Panicum miliaceum L.) is the most water-efficient cereal and one of the earliest domesticated plants. Here we report its high-quality, chromosome-scale genome assembly using a combination of short-read sequencing, single-molecule real-time sequencing, Hi-C, and a high-density genetic map. Phylogenetic analyses reveal two sets of homologous chromosomes that may have merged ~5.6 million years ago, both of which exhibit strong synteny with other grass species. Broomcorn millet contains 55,930 protein-coding genes and 339 microRNA genes. We find Paniceae-specific expansion in several subfamilies of the BTB (broad complex/tramtrack/bric-a-brac) subunit of ubiquitin E3 ligases, suggesting enhanced regulation of protein dynamics may have contributed to the evolution of broomcorn millet. In addition, we identify the coexistence of all three C4 subtypes of carbon fixation candidate genes. The genome sequence is a valuable resource for breeders and will provide the foundation for studying the exceptional stress tolerance as well as C4 biology.

Phylogenetic analysis on the bacteria producing non-volatile fungistatic substances.

This study characterized the soil bacteria producing non-volatile fungistatic substances. Among the 2,100 colonies of soil bacteria randomly isolated from seven agricultural soil samples, 518 isolates (24.67% of total) showed fungistatic activity toward nematophagous fungi Paecilomyces lilacinus and Trichoderma viride by producing non-volatile substances. A phylogenetic analysis based on amplified ribosomal DNA restriction analysis (ARDRA) and 16S rDNA sequence placed the 518 bacteria in three groups of the domain Bacteria: Actinomycetales, Bacillales, and Gammaproteobacteria. Three genera, Arthrobacter, Bacillus, and Pseudomonas, were the most frequently encountered groups.

Phylogenetic diversity of bacteria in an earth-cave in Guizhou province, southwest of China.

The objective of this study was to analyze the phylogenetic composition of bacterial community in the soil of an earth-cave (Niu Cave) using a culture-independent molecular approach. 16S rRNA genes were amplified directly from soil DNA with universally conserved and Bacteria-specific rRNA gene primers and cloned. The clone library was screened by restriction fragment length polymorphism (RFLP), and representative rRNA gene sequences were determined. A total of 115 bacterial sequence types were found in 190 analyzed clones. Phylogenetic sequence analyses revealed novel 16S rRNA gene sequence types and a high diversity of putative bacterial community. Members of these bacteria included Proteobacteria (42.6%), Acidobacteria (18.6%), Planctomycetes (9.0%), Chloroflexi (Green nonsulfur bacteria, 7.5%), Bacteroidetes (2.1%), Gemmatimonadetes (2.7%), Nitrospirae (8.0%), Actinobacteria (High G+C Gram-positive bacteria, 6.4%) and candidate divisions (including the OP3, GN08, and SBR1093, 3.2%). Thirty-five clones were affiliated with bacteria that were related to nitrogen, sulfur, iron or manganese cycles. The comparison of the present data with the data obtained previously from caves based on 16S rRNA gene analysis revealed similarities in the bacterial community components, especially in the high abundance of Proteobacteria and Acidobacteria. Furthermore, this study provided the novel evidence for presence of Gemmatimonadetes, Nitrosomonadales, Oceanospirillales, and Rubrobacterales in a karstic hypogean environment.

A high-quality genome assembly of quinoa provides insights into the molecular basis of salt bladder-based salinity tolerance and the exceptional nutritional value.

Chenopodium quinoa is a halophytic pseudocereal crop that is being cultivated in an ever-growing number of countries. Because quinoa is highly resistant to multiple abiotic stresses and its seed has a better nutritional value than any other major cereals, it is regarded as a future crop to ensure global food security. We generated a high-quality genome draft using an inbred line of the quinoa cultivar Real. The quinoa genome experienced one recent genome duplication about 4.3 million years ago, likely reflecting the genome fusion of two Chenopodium parents, in addition to the γ paleohexaploidization reported for most eudicots. The genome is highly repetitive (64.5% repeat content) and contains 54 438 protein-coding genes and 192 microRNA genes, with more than 99.3% having orthologous genes from glycophylic species. Stress tolerance in quinoa is associated with the expansion of genes involved in ion and nutrient transport, ABA homeostasis and signaling, and enhanced basal-level ABA responses. Epidermal salt bladder cells exhibit similar characteristics as trichomes, with a significantly higher expression of genes related to energy import and ABA biosynthesis compared with the leaf lamina. The quinoa genome sequence provides insights into its exceptional nutritional value and the evolution of halophytes, enabling the identification of genes involved in salinity tolerance, and providing the basis for molecular breeding in quinoa.