Zhigang Meng

The Single-cell Assay for Transposase-Accessible Chromatin with high throughput sequencing (scATAC-seq) has gained increasing popularity in recent years, allowing for chromatin accessibility to be deciphered and gene regulatory networks (GRNs) to be inferred at single-cell resolution. This cutting-edge technology now enables the genome-wide profiling of chromatin accessibility at the cellular level and the capturing of cell-type-specific cis-regulatory elements (CREs) that are masked by cellular heterogeneity in bulk assays. Additionally, it can also facilitate the identification of rare and new cell types based on differences in chromatin accessibility and the charting of cellular developmental trajectories within lineage-related cell clusters. Due to technical challenges and limitations, the data generated from scATAC-seq exhibit unique features, often characterized by high sparsity and noise, even within the same cell type. To address these challenges, various bioinformatic tools have been developed. Furthermore, the application of scATAC-seq in plant science is still in its infancy, with most research focusing on root tissues and model plant species. In this review, we provide an overview of recent progress in scATAC-seq and its application across various fields. We first conduct scATAC-seq in plant science. Next, we highlight the current challenges of scATAC-seq in plant science and major strategies for cell type annotation. Finally, we outline several future directions to exploit scATAC-seq technologies to address critical challenges in plant science, ranging from plant ENCODE(The Encyclopedia of DNA Elements) project construction to GRN inference, to deepen our understanding of the roles of CREs in plant biology.

Protoplasts, which lack cell walls, are ideal research materials for genetic engineering. They are commonly employed in fusion (they can be used for more distant somatic cell fusion to obtain somatic hybrids), genetic transformation, plant regeneration, and other applications. Cotton is grown throughout the world and is the most economically important crop globally. It is therefore critical to study successful extraction and transformation efficiency of cotton protoplasts. In the present study, a cotton callus protoplast extraction method was tested to optimize the ratio of enzymes (cellulase, pectinase, macerozyme R-10, and hemicellulase) used in the procedure. The optimized ratio significantly increased the quantity and activity of protoplasts extracted. We showed that when enzyme concentrations of 1.5% cellulase and 1.5% pectinase, and either 1.5% or 0.5% macerozyme and 0.5% hemicellulase were used, one can obtain increasingly stable protoplasts. We successfully obtained fluorescent protoplasts by transiently expressing fluorescent proteins in the isolated protoplasts. The protoplasts were determined to be suitable for use in further experimental studies. We also studied the influence of plasmid concentration and transformation time on protoplast transformation efficiency. When the plasmid concentration reaches 16 µg and the transformation time is controlled within 12–16 h, the best transformation efficiency can be obtained. In summary, this study presents efficient extraction and transformation techniques for cotton protoplasts.

Although genetic transformation has opened up a new era for cotton molecular breeding, it still suffers from the limitation problem of long transformation periods, which slows down the generation of new cotton germplasms. In this study, LT gene (SV40 large T antigen), which promotes the transformation efficiency of animal cells, was codon-optimized. Its overexpression vector was transformed into cotton. It was observed that EC (embryogenic callus) formation period was 33% shorter and transformation efficiency was slightly higher in the LT T0 generation than that of control. RNA-seq data of NEC (non-embryonic callus) and EC from LT and control revealed that more DEGs (differential expression genes) in NEC were identified than that of EC, indicating LT mainly functioned in NEC. Further KEGG, GO, and transcription factor analyses showed that DEGs were significantly enriched in brassinosteroid biosynthesis pathways and that bHLH, MYB, and AP2/ERF were the top three gene families, which are involved in EC formation. In addition, the key genes related to the auxin pathway were differentially expressed only in LT overexpression NEC, which caused early response, biosynthesis, and transportation of the hormone, resulting in EC earlier formation. In summary, the results demonstrated that LT can promote somatic embryogenesis in cotton, which provides a new strategy for improving cotton transformation and shortening EC formation time.

Crop molecular breeding primarily focuses on increasing the trait of plant yield. An elongator-associated protein, KTI12, is closely associated with plant biomass and yield. KTI12 is involved in developmental processes of most organs, including the leaf, root, flower, and seed, through regulating cell division and differentiation. Previous work has shown that in upland cotton (Gossypium hirsutum), GhKTI12 regulates plant height, flowering, and tolerance to salt and drought stress. However, little is known about the molecular regulation mechanism of GhKTI12 in plant developmental processes. In this study, we identified the main GhKTI12 (Gh_D02G144400) gene and transformed it into tobacco (Nicotonia tabacum cv NC89). From seven transgenic lines, we obtained three (OE5, OE6 and OE8) with high expression of GhKTI12; compared with wild type plants, these three lines exhibited larger plant size, later flowering, and higher seed yield. Microscopic observation revealed that the number of leaf epidermal cells and stem parenchyma cells was increased by ~55%. Biochemical analysis showed that chlorophyll content and starch accumulation were significantly increased in younger leaves at the top canopy of transgenic plants, which may contribute to improved photosynthetic rate and, in turn, increased seed yield. To understand the molecular mechanism of GhKTI12 in transgenic plants development, two lines (OE6 and OE8) with higher expression levels of GhKTI12 were used as representative plants to conduct RNA-seq analysis. Through transcriptome analysis of the plant’s shoot apical meristematic tissue of these two lines, we identified 518 upregulated genes and 406 downregulated genes common to both overexpression lines. A large number of cellular component genes associated with cell division and differentiation, such as RD21, TET8, KTN80, AOX1, AOX2, CP1, and KIC, were found to be upregulated, and genes showing the most downregulation included MADS-box genes related to flowering time, such as MADS6, AP1, AP3, AGL8, AGL6, SEP1, and SEP2. Downregulation of these genes caused delayed flowering time and longer vegetative stage during development. Combined with the upregulation of the yield-related gene RD21, the GhKTI12 transgenic plants could produce a higher seed yield. We here show that the overexpression of GhKTI12 could positively improve key agronomic traits in tobacco by regulating cell proliferation, photosynthesis, and organ development, and suggest that homologs of GhKTI12 may also be important in the genetic improvement of other crop plants.

The authors wish to make the following corrections to their paper

Continual reduction in sequencing cost is expanding the accessibility of genome sequencing data for routine clinical applications. However, the lack of methods to construct machine learning-based predictive models using these datasets has become a crucial bottleneck for the application of sequencing technology in clinics. Here, we develop a new algorithm, eTumorMetastasis, which transforms tumor functional mutations into network-based profiles and identifies network operational gene (NOG) signatures. NOG signatures model the tipping point at which a tumor cell shifts from a state that doesn’t favor recurrence to one that does. We show that NOG signatures derived from genomic mutations of tumor founding clones (i.e., the ‘most recent common ancestor’ of the cells within a tumor) significantly distinguish the recurred and non-recurred breast tumors as well as outperform the most popular genomic test (i.e., Oncotype DX). These results imply that mutations of the tumor founding clones are associated with tumor recurrence and can be used to predict clinical outcomes. As such, predictive tools could be used in clinics to guide treatment routes. Finally, the concepts underlying the eTumorMetastasis pave the way for the application of genome sequencing in predictions for other complex genetic diseases. eTumorMetastasis pseudocode and related data used in this study are available at https://github.com/WangEdwinLab/eTumorMetastasis.

Allotetraploid cotton (Gossypium hirsutum and Gossypium barbadense) are cultivated worldwide for its white fiber. For centuries, conventional breeding approaches increase cotton yield at the cost of extensive erosion of natural genetic variability. Sea Island cotton (G. barbadense) is known for its superior fiber quality, but show poor adaptability as compared to Upland cotton. Here, in this study, we use ethylmethanesulfonate (EMS) as a mutagenic agent to induce genome-wide point mutations to improve the current germplasm resources of Sea Island cotton and develop diverse breeding lines with improved adaptability and excellent economic traits. We determined the optimal EMS experimental procedure suitable for construction of cotton mutant library. At M₆ generation, mutant library comprised of lines with distinguished phenotypes of the plant architecture, leaf, flower, boll, and fiber. Genome-wide analysis of SNP distribution and density in yellow leaf mutant reflected the better quality of mutant library. Reduced photosynthetic efficiency and transmission electron microscopy of yellow leaf mutants revealed the effect of induced mutations at physiological and cellular level. Our mutant collection will serve as the valuable resource for basic research on cotton functional genomics, as well as cotton breeding.

Bacillus thuringiensis (Bt) is a gram negative soil bacterium. This bacterium secretes various proteins during different growth phases with an insecticidal potential against many economically important crop pests. One of the important families of Bt proteins is vegetative insecticidal proteins (Vip), which are secreted into the growth medium during vegetative growth. There are three subfamilies of Vip proteins. Vip1 and Vip2 heterodimer toxins have an insecticidal activity against many Coleopteran and Hemipteran pests. Vip3, the most extensively studied family of Vip toxins, is effective against Lepidopteron. Vip proteins do not share homology in sequence and binding sites with Cry proteins, but share similarities at some points in their mechanism of action. Vip3 proteins are expressed as pyramids alongside Cry proteins in crops like maize and cotton, so as to control resistant pests and delay the evolution of resistance. Biotechnological- and in silico-based analyses are promising for the generation of mutant Vip proteins with an enhanced insecticidal activity and broader spectrum of target insects.

Cotton is an important economic crop affected by different abiotic stresses at different developmental stages. Salinity limits the growth and productivity of crops worldwide. Na⁺/H⁺ antiporters play a key role during the plant development and in its tolerance to salt stress. The aim of the present study was a genome-wide characterization and expression pattern analysis under the salinity stress of the sodium-proton antiporter (NHX) of Gossypium barbadense in comparison with Gossypium hirsutum. In G. barbadense, 25 NHX genes were identified on the basis of the Na⁺_H⁺ exchanger domain. All except one of the G. barbadense NHX transporters have an Amiloride motif that is a known inhibitor of Na⁺ ions in plants. A phylogenetic analysis inferred three classes of GbNHX genes—viz., Vac (GbNHX1, 2 and 4), Endo (GbNHX6), and PM (GbNHX7). A high number of the stress-related cis-acting elements observed in promoters show their role in tolerance against abiotic stresses. The Ka/Ks values show that the majority of GbNHX genes are subjected to strong purifying selection under the course of evolution. To study the functional divergence of G. barbadense NHX transporters, the real-time gene expression was analyzed under salt stress in the root, stem, and leaf tissues. In G. barbadense, the expression was higher in the stem, while in G. hirsutum the leaf and root showed a high expression. Moreover, our results revealed that NHX2 homologues in both species have a high expression under salinity stress at higher time intervals, followed by NHX7. The protein-protein prediction study revealed that GbNHX7 is involved in the CBL-CIPK protein interaction pathway. Our study also provided valuable information explaining the molecular mechanism of Na⁺ transport for the further functional study of Gossypium NHX genes.

Members of the casein kinase 1 (CK1) family are key regulators in numerous cellular signal transduction pathways and in order to prevent the development of certain diseases, CK1 kinase activity needs to be tightly regulated. Modulation of kinase activity by site-specific phosphorylation within the C-terminal regulatory domain of CK1δ has already been shown for several cellular kinases. By using biochemical methods, we now identified residues T161, T174, T176, and S181 within the kinase domain of CK1δ as target sites for checkpoint kinase 1 (Chk1). At least residues T176 and S181 show full conservation among CK1δ orthologues from different eukaryotic species. Enzyme kinetic analysis furthermore led to the hypothesis that site-specific phosphorylation within the kinase domain finally contributes to fine-tuning of CK1δ kinase activity. These data provide a basis for the extension of our knowledge about the role of site-specific phosphorylation for regulation of CK1δ and associated signal transduction pathways.