Impacts of parental genomic divergence in non-syntenic regions on cotton heterosis.

INTRODUCTION: Heterosis has revolutionized crop breeding, enhancing global agricultural production. However, the mechanisms underlying heterosis remain obscure. Xiangzamian 2# (XZM2), a super hybrid upland cotton (Gossypium hirsutum L.) characterized by high-yield heterosis, has been developed and extensively planted in China. OBJECTIVES: We conducted a systematic analysis of CRI12 and J8891, two parents of XZM2. We aimed to reveal the precise genetic information and the role of non-syntenic divergence in shaping heterosis, laying a foundation for advancing understanding of heterosis. METHODS: We de novo assembled high-quality genomes of CRI12 and J8891, and further uncovered abundant genetic variations and non-syntenic regions between the parents. Whole-genome comparison, association analysis, transcriptomic analysis and relative identity-by-descent (rIBD) estimation were conducted to identify structural variations (SVs) and introgressions within non-syntenic blocks and to analyze their impacts on promoting heterosis. RESULTS: Parental genetic divergence increased in non-syntenic regions. Furthermore, these regions, accounting for only 16.71% of the total genome, contained more loci with significantly higher heterotic effects, far exceeding the syntenic background. SVs covered 97.26% of non-syntenic sequences and caused widespread gene expression differences in these regions, driving dynamic complementation of gene expression in the hybrid. A set of SVs were responsible for trait improvement and had positive effects on heterosis, contributing larger heritability than short variations. We characterized numerous parental-specific introgressions from G. barbadense. Specifically, a functional introgression segment within non-syntenic blocks introduced an elite haplotype, which significantly increased lint yield and enhanced heterosis. CONCLUSION: Our study clarified non-syntenic regions to harbor more loci with higher heterotic effects, revealed their importance in promoting heterosis and supported the crucial role of genetic complementation in heterosis. SVs and introgressions were identified as key factors responsible for non-syntenic divergence between the parents. They had important effects on gene expression and trait improvement, positively contributing to heterosis.

Dynamic patterns of gene expressional and regulatory variations in cotton heterosis.

Purpose: Although the application of heterosis has significantly increased crop yield over the past century, the mechanisms underlying this phenomenon still remain obscure. Here, we applied transcriptome sequencing to unravel the impacts of parental expression differences and transcriptomic reprogramming in cotton heterosis. Methods: A high-quality transcriptomic atlas covering 15 developmental stages and tissues was constructed for XZM2, an elite hybrid of upland cotton (Gossypium hirsutum L.), and its parental lines, CRI12 and J8891. This atlas allowed us to identify gene expression differences between the parents and to characterize the transcriptomic reprogramming that occurs in the hybrid. Results: Our analysis revealed abundant gene expression differences between the parents, with pronounced tissue specificity; a total of 1,112 genes exhibited single-parent expression in at least one tissue. It also illuminated transcriptomic reprogramming in the hybrid XZM2, which included both additive and non-additive expression patterns. Coexpression networks between parents and hybrid constructed via weighted gene coexpression network analysis identified modules closely associated with fiber development. In particular, key regulatory hub genes involved in fiber development showed high-parent dominant or over dominant patterns in the hybrid, potentially driving the emergence of heterosis. Finally, high-depth resequencing data was generated and allele-specific expression patterns examined in the hybrid, enabling the dissection of cis- and trans-regulation contributions to the observed expression differences. Conclusion: Parental transcriptional differences and transcriptomic reprogramming in the hybrid, especially the non-additive upregulation of key genes, play an important role in shaping heterosis. Collectively, these findings provide new insights into the molecular basis of heterosis in cotton.

Reference genome sequence and population genomic analysis of peas provide insights into the genetic basis of Mendelian and other agronomic traits.

Peas are essential for human nutrition and played a crucial role in the discovery of Mendelian laws of inheritance. In this study, we assembled the genome of the elite vegetable pea cultivar 'Zhewan No. 1' at the chromosome level and analyzed resequencing data from 314 accessions, creating a comprehensive map of genetic variation in peas. We identified 235 candidate loci associated with 57 important agronomic traits through genome-wide association studies. Notably, we pinpointed the causal gene haplotypes responsible for four Mendelian traits: stem length (Le/le), flower color (A/a), cotyledon color (I/i) and seed shape (R/r). Additionally, we discovered the genes controlling pod form (Mendelian P/p) and hilum color. Our study also involved constructing a gene expression atlas across 22 tissues, highlighting key gene modules related to pod and seed development. These findings provide valuable pea genomic information and will facilitate the future genome-informed improvement of pea crops.

Multi-omics analysis unravels chemical roadmap and genetic basis for peach fruit aroma improvement.

Selection of fruits with enhanced health benefits and superior flavor is an important aspect of peach breeding. Understanding the genetic interplay between appearance and flavor chemicals remains a major challenge. We identify the most important volatiles contributing to consumer preferences for peach, thus establishing priorities for improving flavor quality. We quantify volatiles of a peach population consisting of 184 accessions and demonstrate major reductions in the important flavor volatiles linalool and Z-3-hexenyl acetate in red-fleshed accessions. We identify 474 functional gene regulatory networks (GRNs), among which GRN05 plays a crucial role in controlling both red flesh and volatile content through the NAM/ATAF1/2/CUC (NAC) transcription factor PpBL. Overexpressing PpBL results in reduced expression of PpNAC1, a positive regulator for Z-3-hexenyl acetate and linalool synthesis. Additionally, we identify haplotypes for three tandem PpAATs that are significantly correlated with reduced gene expression and ester content. We develop genetic resources for improvement of fruit quality.

Communication between the stem cell niche and an adjacent differentiation niche through miRNA and EGFR signaling orchestrates exit from the stem cell state in the Drosophila ovary.

The signaling environment, or niche, often governs the initial difference in behavior of an adult stem cell and a derivative that initiates a path towards differentiation. The transition between an instructive stem cell niche and differentiation niche must generally have single-cell resolution, suggesting that multiple mechanisms might be necessary to sharpen the transition. Here, we examined the Drosophila ovary and found that Cap cells, which are key constituents of the germline stem cell (GSC) niche, express a conserved microRNA (miR-124). Surprisingly, loss of miR-124 activity in Cap cells leads to a defect in differentiation of GSC derivatives. We present evidence that the direct functional target of miR-124 in Cap cells is the epidermal growth factor receptor (EGFR) and that failure to limit EGFR expression leads to the ectopic expression of a key anti-differentiation BMP signal in neighboring somatic escort cells (ECs), which constitute a differentiation niche. We further found that Notch signaling connects EFGR activity in Cap cells to BMP expression in ECs. We deduce that the stem cell niche communicates with the differentiation niche through a mechanism that begins with the selective expression of a specific microRNA and culminates in the suppression of the major anti-differentiation signal in neighboring cells, with the functionally important overall role of sharpening the spatial distinction between self-renewal and differentiation environments.

Genome-Wide Analysis of the TCP Transcription Factor Gene Family in Pepper (Capsicum annuum L.).

TCP transcription factors play a key role in regulating various developmental processes, particularly in shoot branching, flower development, and leaf development, and these factors are exclusively found in plants. However, comprehensive studies investigating TCP transcription factors in pepper (Capsicum annuum L.) are lacking. In this study, we identified 27 CaTCP members in the pepper genome, which were classified into Class I and Class II through phylogenetic analysis. The motif analysis revealed that CaTCPs in the same class exhibit similar numbers and distributions of motifs. We predicted that 37 previously reported miRNAs target 19 CaTCPs. The expression levels of CaTCPs varied in various tissues and growth stages. Specifically, CaTCP16, a member of Class II (CIN), exhibited significantly high expression in flowers. Class I CaTCPs exhibited high expression levels in leaves, while Class II CaTCPs showed high expression in lateral branches, especially in the CYC/TB1 subclass. The expression profile suggests that CaTCPs play specific roles in the developmental processes of pepper. We provide a theoretical basis that will assist in further functional validation of the CaTCPs.

High-quality genomes of Bombax ceiba and Ceiba pentandra provide insights into the evolution of Malvaceae species and differences in their natural fiber development.

Members of the Malvaceae family, including Corchorus spp., Gossypium spp., Bombax spp., and Ceiba spp., are important sources of natural fibers. In the past decade, the genomes of several Malvaceae species have been assembled; however, the evolutionary history of Malvaceae species and the differences in their fiber development remain to be clarified. Here, we report the genome assembly and annotation of two natural fiber plants from the Malvaceae, Bombax ceiba and Ceiba pentandra, whose assembled genome sizes are 783.56 Mb and 1575.47 Mb, respectively. Comparative analysis revealed that whole-genome duplication and Gypsy long terminal repeat retroelements have been the major causes of differences in chromosome number (2n = 14 to 2n = 96) and genome size (234 Mb to 2676 Mb) among Malvaceae species. We also used comparative genomic analyses to reconstruct the ancestral Malvaceae karyotype with 11 proto-chromosomes, providing new insights into the evolutionary trajectories of Malvaceae species. MYB-MIXTA-like 3 is relatively conserved among the Malvaceae and functions in fiber cell-fate determination in the epidermis. It appears to perform this function in any tissue where it is expressed, i.e. in fibers on the endocarp of B. ceiba and in ovule fibers of cotton. We identified a structural variation in a cellulose synthase gene and a higher copy number of cellulose synthase-like genes as possible causes of the finer, less spinnable, weaker fibers of B. ceiba. Our study provides two high-quality genomes of natural fiber plants and offers insights into the evolution of Malvaceae species and differences in their natural fiber formation and development through multi-omics analysis.

Gossypium purpurascens genome provides insight into the origin and domestication of upland cotton.

INTRODUCTION: Allotetraploid upland cotton (Gossypium hirsutum L.) is native to the Mesoamerican and Caribbean regions, had been improved in the southern United States by the mid-eighteenth century, was then dispersed worldwide. However, a Hainan Island Native Cotton (HIC) has long been grown extensively on Hainan Island, China. OBJECTIVES: Explore HIC's evolutionary relationship and genomic diversity with other tetraploid cottons, its origin and whether it was used for YAZHOUBU (Yazhou cloth, World Intangible Cultural Heritage) weaving, and the role of structural variations (SVs) in upland cotton domestication. METHODS: We assembled a high-quality genome of one HIC plant. We performed phylogenetic analysis, divergence time estimation, principal component analysis and population differentiation estimation using cotton assemblies and/or resequencing data. SVs were detected by whole-genome comparison. A F2 population was used for linkage analysis and to study effects of SVs. Buoyancy and salt water tolerance tests for seeds were conducted. RESULTS: We found that the HIC belongs to G. purpurascens. G. purpurascens is best classified as a primitive race of G. hirsutum. The potential for long range transoceanic dispersal of G. purpurascens seeds was proved. A set of SVs, selective sweep regions between G. hirsutum races and cultivars, and quantitative trait loci (QTLs) of eleven agronomic traits were obtained. SVs, especially large-scale SVs, were found to have important effects on cotton domestication and improvement. Of them, eight large-scale inversions strongly associated with yield and fiber quality have probably undergone artificial selection in domestication. CONCLUSION: G. purpurascens including HIC is a primitive race of G. hirsutum, probably disperse to Hainan from Central America by floating on ocean currents, may have been partly domesticated, planted and was likely used for YAZHOUBU weaving in Hainan much earlier than the Pre-Columbian period. SV plays an important role in cotton domestication and improvement.

Identification and expression analysis of LEA gene family members in pepper (Capsicum annuum L.).

Pepper (Capsicum annuum L.) is an economically important crop containing capsaicinoids in the seed and placenta, which has various culinary, medical, and industrial applications. Late embryogenesis abundant (LEA) proteins are a large group of hydrophilic proteins participating in the plant stress response and seed development. However, to date there have been no genome-wide analyses of the LEA gene family in pepper. In the present study, 82 LEA genes were identified in the C. annuum genome and classified into nine subfamilies. Most CaLEA genes contain few introns (≤ 2) and are unevenly distributed across 10 chromosomes. Eight pairs of tandem duplication genes and two pairs of segmental duplication genes were identified in the LEA gene family; these duplicated genes were highly conserved and may have performed similar functions during evolution. Expression profile analysis indicated that CaLEA genes exhibited different tissue expression patterns, especially during embryonic development and stress response, particularly in cold stress. Three out of five CaLEA genes showed induced expression upon cold treatment. In summary, we have comprehensively reviewed the LEA gene family in pepper, offering a new perspective on the evolution of this family.

Integration of eQTL and machine learning to dissect causal genes with pleiotropic effects in genetic regulation networks of seed cotton yield.

The dissection of a gene regulatory network (GRN) that complements the genome-wide association study (GWAS) locus and the crosstalk underlying multiple agronomical traits remains a major challenge. In this study, we generate 558 transcriptional profiles of lint-bearing ovules at one day post-anthesis from a selective core cotton germplasm, from which 12,207 expression quantitative trait loci (eQTLs) are identified. Sixty-six known phenotypic GWAS loci are colocalized with 1,090 eQTLs, forming 38 functional GRNs associated predominantly with seed yield. Of the eGenes, 34 exhibit pleiotropic effects. Combining the eQTLs within the seed yield GRNs significantly increases the portion of narrow-sense heritability. The extreme gradient boosting (XGBoost) machine learning approach is applied to predict seed cotton yield phenotypes on the basis of gene expression. Top-ranking eGenes (NF-YB3, FLA2, and GRDP1) derived with pleiotropic effects on yield traits are validated, along with their potential roles by correlation analysis, domestication selection analysis, and transgenic plants.

The Molecular Mechanism of GhbHLH121 in Response to Iron Deficiency in Cotton Seedlings.

Iron deficiency caused by high pH of saline-alkali soil is a major source of abiotic stress affecting plant growth. However, the molecular mechanism underlying the iron deficiency response in cotton (Gossypium hirsutum) is poorly understood. In this study, we investigated the impacts of iron deficiency at the cotton seedling stage and elucidated the corresponding molecular regulation network, which centered on a hub gene GhbHLH121. Iron deficiency induced the expression of genes with roles in the response to iron deficiency, especially GhbHLH121. The suppression of GhbHLH121 with virus-induced gene silence technology reduced seedlings' tolerance to iron deficiency, with low photosynthetic efficiency and severe damage to the structure of the chloroplast. Contrarily, ectopic expression of GhbHLH121 in Arabidopsis enhanced tolerance to iron deficiency. Further analysis of protein/protein interactions revealed that GhbHLH121 can interact with GhbHLH IVc and GhPYE. In addition, GhbHLH121 can directly activate the expression of GhbHLH38, GhFIT, and GhPYE independent of GhbHLH IVc. All told, GhbHLH121 is a positive regulator of the response to iron deficiency in cotton, directly regulating iron uptake as the upstream gene of GhFIT. Our results provide insight into the complex network of the iron deficiency response in cotton.

Global identification of natural antisense transcripts in Gossypium hirsutum and Gossypium barbadense under chilling stress.

Natural antisense transcripts (NATs) in model plants have been recognized as important regulators of gene expression under abiotic stresses. However, the functional roles of NATs in crops under low temperature are still unclear. Here, we identified 815 and 689 NATs from leaves of Gossypium hirsutum and G. barbadense under chilling stress. Among those, 224 NATs were identified as interspecific homologs between the two species. The correlation coefficients for expression of NATs and their cognate sense genes (CSG) were 0.43 and 0.37 in G. hirsutum and G. barbadense, respectively. Furthermore, expression of interspecific NATs and CSGs alike was highly consistent under chilling stress with correlation coefficients of 0.90-0.91. Four cold-associated NATs were selected for functional validation using virus-induced gene silencing (VIGS). Our results suggest that CAN1 engage in the molecular regulation of chilling stress by regulating SnRK2.8 expression. This highly conserved NAT have valuable potential for applications in breeding cold-tolerant cotton.

Mapping intron retention events contributing to complex traits using splice quantitative trait locus.

BACKGROUND: Alternative splicing (AS) of mRNA plays an important roles in transcriptome diversity, involving regulation of plant growth and stress response. Understanding the variation of AS events underlying GWAS loci in a crop population can provide insight into the molecular mechanisms of complex agronomic traits. To date, genome-wide association studies relating AS events to agronomic traits have rarely been conducted at the population level in crops. RESULTS: Here, a pipeline was constructed to identify candidate AS events related to complex traits. Firstly, ovule transcriptome data were used to characterize intron retention (IR), the predominant type of AS in plants, on a genome-wide scale. This was done in a natural population consisting of 279 upland cotton lines. Secondly, splice quantitative trait locus (sQTL) analysis was carried out, which yielded a total of 2295 sQTLs involving 1607 genes. Of these, 14.25% (n = 427) were cis-regulatory loci. Integration with expression quantitative trait loci (eQTL) revealed that 53 (21.4%) cis-sGenes were regulated by both cis-sQTLs and cis-eQTLs. Finally, co-localization analysis integrated with GWAS loci in this population showed 32 cis-QTLs to be co-located with genetic regulatory loci related to fiber yield and quality traits, indicating that sQTLs are likely to participate in regulating cotton fiber yield and quality. An in-depth evaluation confirmed that differences in the IR rates of sQTL-regulated candidate genes such as GhLRRK1 and GhGC1 are associated with lint percentage (LP), which has potential in breeding applications. CONCLUSION: This study provides a clue that AS of mRNA has an impact on crop yield, along with functional sQTLs are new genetic resources for cotton precision breeding.

An upgraded method of high-throughput chromosome conformation capture (Hi-C 3.0) in cotton (Gossypium spp.).

High-throughput chromosome conformation capture (Hi-C) technology has been applied to explore the chromatin interactions and shed light on the biological functions of three-dimensional genomic features. However, it remains challenging to guarantee the high quality of Hi-C library in plants and hence the reliable capture of chromatin structures, especially loops, due to insufficient fragmentation and low efficiency of proximity ligations. To overcome these deficiencies, we optimized the parameters of the Hi-C protocol, principally the cross-linking agents and endonuclease fragmentation strategy. The double cross-linkers (FA+DSG) and double restriction enzymes (DpnII+DdeI) were utilized. Thus, a systematic in situ Hi-C protocol was designed using plant tissues embedded with comprehensive quality controls to monitor the library construction. This upgraded method, termed Hi-C 3.0, was applied to cotton leaves for trial. In comparison with the conventional Hi-C 2.0, Hi-C 3.0 can obtain more than 50% valid contacts at a given sequencing depth to improve the signal-to-noise ratio. Hi-C 3.0 can furthermore enhance the capturing of loops almost as twice as that of Hi-C 2.0. In addition, Hi-C 3.0 showed higher efficiency of compartment detection and identified compartmentalization more accurately. In general, Hi-C 3.0 contributes to the advancement of the Hi-C method in plants by promoting its capability on decoding the chromatin organization.

Metabolomics Provides New Insights into Host Manipulation Strategies by Asobara japonica (Hymenoptera: Braconidae), a Fruit Fly Parasitoid.

Asobara japonica (Hymenoptera: Braconidae) is an endoparasitoid wasp that can successfully parasitize a wide range of host species across the Drosophila genus, including the invasive crop pest Drosophila suzukii. Parasitoids are capable of regulating the host metabolism to produce the nutritional metabolites for the survival of their offspring. Here, we intend to investigate the metabolic changes in D. melanogaster hosts after parasitization by A. japonica, using the non-targeted LC-MS (liquid chromatography-mass spectrometry) metabolomics analysis. In total, 3043 metabolites were identified, most of which were not affected by A. japonica parasitization. About 205 metabolites were significantly affected in parasitized hosts in comparison to non-parasitized hosts. The changed metabolites were divided into 10 distinct biochemical groups. Among them, most of the lipid metabolic substances were significantly decreased in parasitized hosts. On the contrary, most of metabolites associated with the metabolism of amino acids and sugars showed a higher abundance of parasitized hosts, and were enriched for a wide range of pathways. In addition, eight neuromodulatory-related substances were upregulated in hosts post A. japonica parasitization. Our results reveal that the metabolites are greatly changed in parasitized hosts, which might help uncover the underlying mechanisms of host manipulation that will advance our understanding of host-parasitoid coevolution.

Biotinylated Tn5 transposase-mediated CUT&Tag efficiently profiles transcription factor-DNA interactions in plants.

In contrast to CUT&Tag approaches for profiling bulk histone modifications, current CUT&Tag methods for analysing specific transcription factor (TF)-DNA interactions remain technically challenging due to TFs having relatively low abundance. Moreover, an efficient CUT&Tag strategy for plant TFs is not yet available. Here, we first applied biotinylated Tn5 transposase-mediated CUT&Tag (B-CUT&Tag) to produce high-quality libraries for interrogating TF-DNA interactions. B-CUT&Tag combines streptavidin-biotin-based DNA purification with routine CUT&Tag, optimizing the removal of large amounts of intact chromatin not targeted by specific TFs. The biotinylated chromatin fragments are then purified for construction of deep sequencing libraries or qPCR analysis. We applied B-CUT&Tag to probe genome-wide DNA targets of Squamosa promoter-binding-like protein 9 (SPL9), a well-established TF in Arabidopsis; the resulting profiles were efficient and consistent in demonstrating its well-established target genes in juvenile-adult transition/flowering, trichome development, flavonoid biosynthesis, wax synthesis and branching. Interestingly, our results indicate functions of AtSPL9 in modulating growth-defence trade-offs. In addition, we established a method for applying qPCR after CUT&Tag (B-CUT&Tag-qPCR) and successfully validated the binding of SPL9 in Arabidopsis and PHR2 in rice. Our study thus provides a convenient and highly efficient CUT&Tag strategy for profiling TF-chromatin interactions that is widely applicable to the annotation of cis-regulatory elements for crop improvement.

HEAT-RESPONSIVE PROTEIN regulates heat stress via fine-tuning ethylene/auxin signaling pathways in cotton.

Plants sense and respond to fluctuating temperature and light conditions during the circadian cycle; however, the molecular mechanism underlying plant adaptability during daytime warm conditions remains poorly understood. In this study, we reveal that the ectopic regulation of a HEAT RESPONSIVE PROTEIN (GhHRP) controls the adaptation and survival of cotton (Gossypium hirsutum) plants in response to warm conditions via modulating phytohormone signaling. Increased ambient temperature promptly enhanced the binding of the phytochrome interacting factor 4 (GhPIF4)/ethylene-insensitive 3 (GhEIN3) complex to the GhHRP promoter to increase its mRNA level. The ectopic expression of GhHRP promoted the temperature-dependent accumulation of GhPIF4 transcripts and hypocotyl elongation by triggering thermoresponsive growth-related genes. Notably, the upregulation of the GhHRP/GhPIF4 complex improved plant growth via modulating the abundance of Arabidopsis thaliana auxin biosynthetic gene YUCCA8 (AtYUC8)/1-aminocyclopropane-1-carboxylate synthase 8 (AtACS8) for fine-tuning the auxin/ethylene interplay, ultimately resulting in decreased ethylene biosynthesis. GhHRP thus protects chloroplasts from photo-oxidative bursts via repressing AtACS8 and AtACS7 and upregulating AtYUC8 and the heat shock transcription factors (HSFA2), heat shock proteins (HSP70 and HSP20). Strikingly, the Δhrp disruption mutant exhibited compromised production of HSP/YUC8 that resulted in an opposite phenotype with the loss of the ability to respond to warm conditions. Our results show that GhHRP is a heat-responsive signaling component that assists plants in confronting the dark phase and modulates auxin signaling to rescue growth under temperature fluctuations.

Phylogenomic analysis of cytochrome P450 multigene family and its differential expression analysis in pepper (Capsicum annuum L.).

Plant cytochrome P450 is a multifamily enzyme widely involved in biochemical reactions for the synthesis of antioxidants, pigments, structural polymers, and defense-related compounds. Pepper (Capsicum annuum L.) is an economically important plant. A comprehensive identification and characterization of P450 genes would provide valuable information on the evolutionary relationships of genes and their functional characteristics. In this study, we identified P450 genes in pepper with the aid of bioinformatics methods to investigate the phylogenetic relation, gene structure, chromosomal localization, duplicated events, and collinearity among Solanaceae species. We identified and classified 478 genes of P450 from the pepper genome into two major clades and nine subfamilies through phylogenetic analysis. Massive duplication events were found in the P450 gene family, which may explain the expansion of the P450 gene family. In addition, we also found that these duplication genes may have undergone strict purification selection during evolution. Gene expression analysis showed that some P450 genes that belong to clan 71 in pepper may play an important role in placenta and pericarp development. Through quantitative real-time polymerase chain reaction and transcriptome analysis, we also found that many P450 genes were related to defensive and phytohormone response in pepper. These findings provide insight for further studies to identify the biological functions of the P450 genes in pepper.

Development of the engineered "glanded plant and glandless seed" cotton.

After fiber, cottonseed is the second most important by-product of cotton production. However, high concentrations of toxic free gossypol deposited in the glands of the cottonseed greatly hamper its effective usage as food or feed. Here, we developed a cotton line with edible cottonseed by specifically silencing the endogenous expression of GoPGF in the seeds, which led to a glandless phenotype with an ultra-low gossypol content in the seeds and nearly normal gossypol in other parts of the plants. This engineered cotton maintains normal resistance to insect pests, but the gossypol content in the seeds dropped by 98%, and thus, it can be consumed directly as food. The trait of a low gossypol content in the cottonseeds was stable and heritable, while the protein, oil content, and fiber yield or quality were nearly unchanged compared to the transgenic receptor W0. In addition, comparative transcriptome analysis showed that down-regulated genes in the ovules of the glandless cotton were enriched in terpenoid biosynthesis, indicating the underlying relationship between gland formation and gossypol biosynthesis. These results pave the way for the comprehensive utilization of cotton as a fiber, oil, and feed crop in the future.

Profiling of H3K4me3 Modification in Plants using Cleavage under Targets and Tagmentation.

Epigenomic regulation at the chromatic level, including DNA and histone modifications, behaviors of transcription factors, and non-coding RNAs with their recruited proteins, lead to temporal and spatial control of gene expression. Cleavage under targets and tagmentation (CUT&Tag) is an enzyme-tethering method in which the specific chromatin protein is firstly recognized by its specific antibody, and then the antibody tethers a protein A-transposase (pA-Tn5) fusion protein, which cleaves the targeted chromatin in situ by the activation of magnesium ions. Here, we provide our previously published CUT&Tag protocol using intact nuclei isolated from allortetraploid cotton leaves with modification. This step-by-step protocol can be used for epigenomic research in plants. In addition, substantial modifications for plant nuclei isolation are provided with critical comments.