Identification and function analysis of GABA branch three gene families in the cotton related to abiotic stresses.

γ -aminobutyric acid (GABA) is closely related to the growth, development and stress resistance of plants. Combined with the previous study of GABA to promote the cotton against abiotic stresses, the characteristics and expression patterns of GABA branch gene family laid the foundation for further explaining its role in cotton stress mechanism. Members of GAD, GAB-T and SSADH (three gene families of GABA branch) were identified from the Gossypium hirsutum, Gossypium barbadense, Gossypium arboreum and Gossypium raimondii genome. The GABA branch genes were 10 GAD genes, 4 GABA-T genes and 2 SSADH genes. The promoter sequences of genes mainly contains response-related elements such as light, hormone and environment.Phylogenetic analysis shows that GAD indicating that even in the same species, the homologous sequences in the family. The GABA-T gene of each cotton genus was in sum the family had gene loss in the process of dicotyledon evolution. SSADH families Gossypium hirsutum, Gossypium barbadense, Gossypium arboreum and Gossypium raimondii were closely related to the dicot plants.GABA gene is involved in the regulation of salt stress and high temperature in Gossypium hirsutum.GABA attenuated part of the abiotic stress damage by increasing leaf protective enzyme activity and reducing reactive oxygen species production.This lays the foundation for a thorough analysis of the mechanism of GABA in cotton stress resistance.

A combination of genome-wide and transcriptome-wide association studies reveals genetic elements leading to male sterility during high temperature stress in cotton.

Global warming has reduced the productivity of many field-grown crops, as the effects of high temperatures can lead to male sterility in such plants. Genetic regulation of the high temperature (HT) response in the major crop cotton is poorly understood. We determined the functionality and transcriptomes of the anthers of 218 cotton accessions grown under HT stress. By analyzing transcriptome divergence and implementing a genome-wide association study (GWAS), we identified three thermal tolerance associated loci which contained 75 protein coding genes and 27 long noncoding RNAs, and provided expression quantitative trait loci (eQTLs) for 13 132 transcripts. A transcriptome-wide association study (TWAS) confirmed six causal elements for the HT response (three genes overlapped with the GWAS results) which are involved in protein kinase activity. The most susceptible gene, GhHRK1, was confirmed to be a previously uncharacterized negative regulator of the HT response in both cotton and Arabidopsis. These functional variants provide a new understanding of the genetic basis for HT tolerance in male reproductive organs.

Genetic diversity and structure of elite cotton germplasm (Gossypium hirsutum L.) using genome-wide SNP data.

Cotton (Gossypium spp.) is the most important natural textile fiber crop, and Gossypium hirsutum L. is responsible for 90% of the annual cotton crop in the world. Information on cotton genetic diversity and population structure is essential for new breeding lines. In this study, we analyzed population structure and genetic diversity of 288 elite Gossypium hirsutum cultivar accessions collected from around the world, and especially from China, using genome-wide single nucleotide polymorphisms (SNP) markers. The average polymorphsim information content (PIC) was 0.25, indicating a relatively low degree of genetic diversity. Population structure analysis revealed extensive admixture and identified three subgroups. Phylogenetic analysis supported the subgroups identified by STRUCTURE. The results from both population structure and phylogenetic analysis were, for the most part, in agreement with pedigree information. Analysis of molecular variance revealed a larger amount of variation was due to diversity within the groups. Establishment of genetic diversity and population structure from this study could be useful for genetic and genomic analysis and systematic utilization of the standing genetic variation in upland cotton.

Genomic and epigenomic insights into the mechanism of cold response in upland cotton (Gossypium hirsutum).

Functional genome research, including gene transcriptional and posttranslational modifications of histones, can benefit greatly from a high-quality genome assembly. Histone modification plays a significant role in modulating the responses to abiotic stress in plants. However, there are limited reports on the involvement of dynamic changes in histone modification in cold stress response in upland cotton. In this study, the genome of an elite accession, YM11, with considerable cold stress tolerance was de novo assembled, which yielded a genome of 2343.06 Mb with a contig N50 of 88.96 Mb, and a total of 73,821 protein-coding gene models were annotated. Comparisons among YM11 and five Gossypium allopolyploid cotton assemblies highlighted a large amount of structural variations and presence/absence variations. We analyzed transcriptome and metabolome changes in YM11 seedlings subjected to cold stress. Using the CUT&Tag method, genome-wide H3K4me3 and H3K9ac modification patterns and effect of histone changes on gene expression were profiled during cold stress. Significant and consistently changing histone modifications and the gene expressions were screened, of which transcription factors (TFs) were highlighted. Our results suggest a positive correlation between the changes in H3K4me3, H3K9ac modifications and cold stress-responsive gene activation. This genome assembly and comprehensive analysis of genome-wide histone modifications and gene expression provide insights into the genomic variation and epigenetic responses to cold stress in upland cotton.

Driving Principle and Stability Analysis of Vertical Comb-Drive Actuator for Scanning Micromirrors.

We have developed a manufacturing process for micromirrors based on microelectromechanical systems (MEMS) technology. The process involves designing an electrostatic vertically comb-driven actuator and utilizing a self-alignment process to produce a height difference between the movable comb structure and the fixed comb structure of the micromirror. To improve the stability of the micromirror, we propose four instability models in micromirror operation with the quasi-static driving principle and structure of the micromirror considered, which can provide a basic guarantee for the performance of vertical comb actuators. This analysis pinpoints factors leading to instability, including the left and right gap of the movable comb, the torsion beams of the micromirror, and the comb-to-beams distance. Ultimately, the voltages at which device failure occurs can be determined. We successfully fabricated a one-dimensional micromirror featuring a 0.8 mm mirror diameter and a 30 µm device layer thickness. The height difference between the movable and fixed comb structures was 10 µm. The micromirror was able to achieve a static mechanical angle of 2.25° with 60 V@DC. Stable operation was observed at voltages below 60 V, in close agreement with the theoretical calculations and simulations. At the driving voltage of 80 V, we observed the longitudinal displacement movement of the comb fingers. Furthermore, at a voltage of 129 V, comb adhesion occurred, resulting in device failure. This failure voltage corresponds to the lateral torsional failure voltage.

Whole-Genome Bisulfite Sequencing (WGBS) Analysis of Gossypium hirsutum under High-Temperature Stress Conditions.

BACKGROUND: DNA methylation is an important part of epigenetic regulation and plays an important role in the response of plants to adverse stress. METHODS: In this study, whole-genome bisulfite sequencing (WGBS) was performed on the high-temperature-resistant material Xinluzao 36 and the high-temperature-sensitive material Che 61-72 at 0 h and 12 h under high-temperature stress conditions. RESULTS: The results revealed that the Gossypium hirsutum methylation levels of CG and CHG (H = A, C, or T) decreased after the high-temperature stress treatment, and the methylation level of the A subgenome was significantly greater than that of the D subgenome. The methylation level of CHH increased, and the methylation level of CHH in the D subgenome was significantly greater than that in the A subgenome after high-temperature stress treatment. The methylation density of CG is lower than that of CHG and CHH, and the methylation density of the middle region of chromosomes is greater than that of both ends, which is opposite to the distribution density of genes. There were 124 common differentially methylated genes in the CG, CHG, and CHH groups, and 5130 common DEGs and differentially methylated genes were found via joint analysis with RNA-seq; these genes were significantly enriched in the biosynthesis of plant hormones, thiamine metabolism, glutathione metabolism, and tyrosine metabolism pathways. DNA methylation did not affect the expression of many genes (accounting for 85.68% of the differentially methylated genes), DNA methylation-promoted gene expression was located mainly in the downstream region of the gene or gene body, and the expression of inhibitory genes was located mainly in the upstream region of the gene. CONCLUSIONS: This study provides a theoretical basis for further exploration of the gene expression and functional regulatory mechanism of G. hirsutum DNA methylation under high-temperature stress conditions.

RNA-Seq Analysis and Candidate Gene Mining of Gossypium hirsutum Stressed by Verticillium dahliae Cultured at Different Temperatures.

The occurrence and spread of Verticillium dahliae (V. dahliae) in cotton depends on the combined effects of pathogens, host plants, and the environment, among which temperature is one of the most important environmental factors. Studying how temperature impacts the occurrence of V. dahliae in cotton and the mechanisms governing host defense responses is crucial for disease prevention and control. Understanding the dual effects of temperature on both pathogens and hosts can provide valuable insights for developing effective strategies to manage this destructive fungal infection in cotton. This study was based on the deciduous V. dahliae Vd-3. Through cultivation at different temperatures, Vd-3 formed the most microsclerotia and had the largest colony diameter at 25 °C. Endospore toxins were extracted, and 48 h was determined to be the best pathogenic time point for endotoxins to infect cotton leaves through a chlorophyll fluorescence imaging system and phenotypic evaluation. Transcriptome sequencing was performed on cotton leaves infected with Vd-3 endotoxins for 48 h at different culture temperatures. A total of 34,955 differentially expressed genes (DEGs) were identified between each temperature and CK (no pathogen inoculation), including 17,422 common DEGs. The results of the enrichment analysis revealed that all the DEGs were involved mainly in photosynthesis and sugar metabolism. Among the 34,955 DEGs, genes in the biosynthesis and signal transduction pathways of jasmonic acid (JA), salicylic acid (SA), and ethylene (ET) were identified, and their expression patterns were determined. A total of 5652 unique DEGs were clustered into six clusters using the k-means clustering algorithm, and the functions and main transcription factors (TFs) of each cluster were subsequently annotated. In addition, we constructed a gene regulatory network via weighted correlation network analysis (WGCNA) and identified twelve key genes related to cotton defense against V. dahliae at different temperatures, including four genes encoding transcription factors. These findings provide a theoretical foundation for investigating temperature regulation in V. dahliae infecting cotton and introduce novel genetic resources for enhancing resistance to this disease in cotton plants.

Integrating single-cell and bulk RNA sequencing reveals CK19 + cancer stem cells and their specific SPP1 + tumor-associated macrophage niche in HBV-related hepatocellular carcinoma.

PURPOSE: Cytokeratin 19-positive cancer stem cells (CK19 + CSCs) and their tumor-associated macrophages (TAMs) have not been fully explored yet in the hepatitis B virus (HBV)-related hepatocellular carcinoma (HCC). EXPERIMENTAL DESIGN: Single-cell RNA sequencing was performed on the viable cells obtained from 11 treatment-naïve HBV-associated HCC patients, including 8 CK19 + patients, to elucidate their transcriptomic landscape, CK19 + CSC heterogeneity, and immune microenvironment. Two in-house primary HCC cohorts (96 cases-related HBV and 89 cases with recurrence), TCGA external cohort, and in vitro and in vivo experiments were used to validate the results. RESULTS: A total of 64,581 single cells derived from the human HCC and adjacent normal tissues were sequenced, and 11 cell types were identified. The result showed that CK19 + CSCs were phenotypically and transcriptionally heterogeneous, co-expressed multiple hepatics CSC markers, and were positively correlated with worse prognosis. Moreover, the SPP1 + TAMs (TAM_SPP1) with strong M2-like features and worse prognosis were specifically enriched in the CK19 + HCC and promoted tumor invasion and metastasis by activating angiogenesis. Importantly, matrix metalloproteinase 9 (MMP9) derived from TAM_SPP1, as the hub gene of CK19 + HCC, was activated by the VEGFA signal. CONCLUSIONS: This study revealed the heterogeneity and stemness characteristics of CK19 + CSCs and specific immunosuppressive TAM_SPP1 in CK19 + HCC. The VEGFA signal can activate TAM_SPP1-derived MMP9 to promote the invasion and metastasis of CK19 + HCC tumors. This might provide novel insights into the clinical treatment of HCC patients.

Genome-wide identification and expression reveal the involvement of the FCS-like zinc finger (FLZ) gene family in Gossypium hirsutum at low temperature.

FCS-like zinc finger (FLZ) is a plant-specific gene family that plays an important regulatory role in plant growth and development and its response to stress. However, studies on the characteristics and functions of cotton FLZ family genes are still lacking. This study systematically identified members of the cotton FLZ gene family based on cotton genome data. The cotton FLZ family genes were systematically analyzed by bioinformatics, and their expression patterns in different tissues and under low-temperature stress were analyzed by transcriptome and qRT-PCR. The G. hirsutum genome contains 56 FLZ genes distributed on 20 chromosomes, and most of them are located in the nucleus. According to the number and evolution analysis of FLZ family genes, FLZ family genes can be divided into five subgroups in cotton. The G. hirsutum FLZ gene has a wide range of tissue expression types, among which the expression is generally higher in roots, stems, leaves, receptacles and calyx. Through promoter analysis, it was found that it contained the most cis-acting elements related to methyl jasmonate (MeJA) and abscisic acid (ABA). Combined with the promoter and qRT-PCR results, it was speculated that GhFLZ11, GhFLZ25, GhFLZ44 and GhFLZ55 were involved in the response of cotton to low-temperature stress. Taken together, our findings suggest an important role for the FLZ gene family in the cotton response to cold stress. This study provides an important theoretical basis for further research on the function of the FLZ gene family and the molecular mechanism of the cotton response to low temperature.

Experimental Investigation of Vibration Isolator for Large Aperture Electromagnetic MEMS Micromirror.

The Micro-Electro-Mechanical-System (MEMS) micromirror has shown great advantages in Light Detection and Ranging (LiDAR) for autonomous vehicles. The equipment on vehicles is usually exposed to environmental vibration that may degrade or even destroy the flexure of the micromirror for its delicate structure. In this work, a mechanical low-pass filter (LPF) acting as a vibration isolator for a micromirror is proposed. The research starts with the evaluation of vibration influences on the micromirror by theoretical calculation and simulation. The results illustrate that mechanical load concentrates at the slow flexure of the micromirror as it is excited to resonate in second-order mode (named piston mode) in Z-direction vibration. A specific LPF for the micromirror is designed to attenuate the response to high-frequency vibration, especially around piston mode. The material of the LPF is a beryllium-copper alloy, chosen for its outstanding properties of elasticity, ductility, and fatigue resistance. To measure the mechanical load on the micromirror in practical, the on-chip piezoresistive sensor is utilized and a relevant test setup is built to validate the effect of the LPF. Micromirrors with or without the LPF are both tested under 10 g vibration in the Z-direction. The sensor output of the device with the LPF is 35.9 mV in piston mode, while the device without the LPF is 70.42 mV. The attenuation ratio is 0.51. This result demonstrates that the LPF structure can effectively reduce the stress caused by piston mode vibration.

Whole Transcriptome Sequencing Reveals Drought Resistance-Related Genes in Upland Cotton.

China, particularly the cotton-growing province of Xinjiang, is experiencing acute agricultural water shortages, stifling the expansion of the cotton sector. Discovering drought resistance genes in cotton and generating high-quality, drought-resistant cotton varieties through molecular breeding procedures are therefore critical to the cotton industry's success. The drought-resistant cotton variety Xinluzhong No. 82 and the drought-sensitive cotton variety Kexin No. 1 were utilised in this study to uncover a batch of drought-resistant candidate genes using whole transcriptome sequencing. The following are the key research findings: A competing endogenous RNA network (ceRNA) was built using complete transcriptional sequencing to screen the core genes in the core pathway, and two drought-related candidate genes were discovered. It was found that γ-aminobutyric acid aminotransferase (GhGABA-T, Gohir.A11G156000) was upregulated at 0 h vs. 12 h and downregulated at 12 h vs. 24 h. L-Aspartate oxidase (GhAO, Gohir.A07G220600) was downregulated at 0 h vs. 12 h and upregulated at 12 h vs. 24 h. GABA-T is analogous to a pyridoxal phosphate-dependent transferase superfamily protein (POP2) in Arabidopsis thaliana and influences plant drought resistance by controlling γ-aminobutyric acid (GABA) concentration. The analogue of GhAO in A. thaliana is involved in the early steps of nicotinamide adenine dinucleotide (NAD) production as well as in plant antioxidant responses. This study revealed that gene expression regulatory networks can be used for rapid screening of reliable drought resistance genes and then utilised to validate gene function.

Genome-Wide Comparative Analysis of Heat Shock Transcription Factors Provides Novel Insights for Evolutionary History and Expression Characterization in Cotton Diploid and Tetraploid Genomes.

Heat shock transcription factors (HSFs) are involved in environmental stress response and plant development, such as heat stress and flowering development. According to the structural characteristics of the HSF gene family, HSF genes were classified into three major types (HSFA, HSFB, and HSFC) in plants. Using conserved domains of HSF genes, we identified 621 HSF genes among 13 cotton genomes, consisting of eight diploid and five tetraploid genomes. Phylogenetic analysis indicated that HSF genes among 13 cotton genomes were grouped into two different clusters: one cluster contained all HSF genes of HSFA and HSFC, and the other cluster contained all HSF genes of HSFB. Comparative analysis of HSF genes in Arabidopsis thaliana, Gossypium herbaceum (A1), Gossypium arboreum (A2), Gossypium raimondii (D5), and Gossypium hirsutum (AD1) genomes demonstrated that four HSF genes were inherited from a common ancestor, A0, of all existing cotton A genomes. Members of the HSF gene family in G. herbaceum (A1) genome indicated a significant loss compared with those in G. arboretum (A2) and G. hirsutum (AD1) A genomes. However, HSF genes in G. raimondii (D5) showed relative loss compared with those in G. hirsutum (AD1) D genome. Analysis of tandem duplication (TD) events of HSF genes revealed that protein-coding genes among different cotton genomes have experienced TD events, but only the two-gene tandem array was detected in Gossypium thurberi (D1) genome. The expression analysis of HSF genes in G. hirsutum (AD1) and Gossypium barbadense (AD2) genomes indicated that the expressed HSF genes were divided into two different groups, respectively, and the expressed HSF orthologous genes between the two genomes showed totally different expression patterns despite the implementation of the same abiotic stresses. This work will provide novel insights for the study of evolutionary history and expression characterization of HSF genes in different cotton genomes and a widespread application model for the study of HSF gene families in plants.

Nanopore-Based Comparative Transcriptome Analysis Reveals the Potential Mechanism of High-Temperature Tolerance in Cotton (Gossypium hirsutum L.).

Extreme high temperatures are threatening cotton production around the world due to the intensification of global warming. To cope with high-temperature stress, heat-tolerant cotton cultivars have been bred, but the heat-tolerant mechanism remains unclear. This study selected heat-tolerant ('Xinluzao36') and heat-sensitive ('Che61-72') cultivars of cotton treated with high-temperature stress as plant materials and performed comparative nanopore sequencing transcriptome analysis to reveal the potential heat-tolerant mechanism of cotton. Results showed that 120,605 nonredundant sequences were generated from the raw reads, and 78,601 genes were annotated. Differentially expressed gene (DEG) analysis showed that a total of 19,600 DEGs were screened; the DEGs involved in the ribosome, heat shock proteins, auxin and ethylene signaling transduction, and photosynthesis pathways may be attributed to the heat tolerance of the heat-tolerant cotton cultivar. This study also predicted a total of 5118 long non-coding RNAs (lncRNAs)and 24,462 corresponding target genes. Analysis of the target genes revealed that the expression of some ribosomal, heat shock, auxin and ethylene signaling transduction-related and photosynthetic proteins may be regulated by lncRNAs and further participate in the heat tolerance of cotton. This study deepens our understandings of the heat tolerance of cotton.

Genome-Wide Association Analysis of Salt-Tolerant Traits in Terrestrial Cotton at Seedling Stage.

Soil salinization is the main abiotic stress factor affecting agricultural production worldwide, and salt stress has a significant impact on plant growth and development. Cotton is one of the most salt-tolerant crops. Therefore, the selection and utilization of salt-tolerant germplasm resources and the excavation of salt resistance genes play important roles in improving cotton production in saline-alkali soils. In this study, we analysed the population structure and genetic diversity of a total 149 cotton plant materials including 137 elite Gossypium hirsutum cultivar accessions collected from China and 12 elite Gossypium hirsutum cultivar accessions collected from around the world. Illumina Cotton SNP 70 K was used to obtain genome-wide single-nucleotide polymorphism (SNP) data for 149 elite Gossypium hirsutum cultivar accessions, and 18,430 highly consistent SNP loci were obtained by filtering. It was assessed by using PCA principal component analysis so that the 149 elite Gossypium hirsutum cultivar accessions could be divided into two subgroups, including subgroup 1 with 78 materials and subgroup 2 with 71 materials. Using the obtained SNP and other marker genotype test results, under salt stress, the salt tolerance traits 3d Germination potential, 3d Radicle length drop rate, 7d Germination rate, 7d Radicle length drop rate, 7d Germination weight, 3d Radicle length, 7d Radicle length, Relative Germination potential, Relative Germination rate, 7d Radicle weight drop rate, Salt tolerance index 3d Germination potential index, 3d Radicle length index, 7d Radicle length index, 7d Radicle weight index and 7d Germination rate index were evaluated by GWAS (genome-wide association analysis). A total of 27 SNP markers closely related to the salt tolerance traits and 15 SNP markers closely related to the salt tolerance index were detected. At the SNP locus associated with phenotyping, Gh_D01G0943, Gh_D01G0945, Gh_A01G0906, Gh_A01G0908, Gh_D08G1308 and Gh_D08G1309 related to plant salt tolerance were detected, and they were found to be involved in intracellular transport, sucrose synthesis, osmotic pressure balance, transmembrane transport, N-glycosylation, auxin response and cell amplification. This study provides a theoretical basis for the selection and breeding of salt-tolerant upland cotton varieties.

The genomic basis of geographic differentiation and fiber improvement in cultivated cotton.

Large-scale genomic surveys of crop germplasm are important for understanding the genetic architecture of favorable traits. The genomic basis of geographic differentiation and fiber improvement in cultivated cotton is poorly understood. Here, we analyzed 3,248 tetraploid cotton genomes and confirmed that the extensive chromosome inversions on chromosomes A06 and A08 underlies the geographic differentiation in cultivated Gossypium hirsutum. We further revealed that the haplotypic diversity originated from landraces, which might be essential for understanding adaptative evolution in cultivated cotton. Introgression and association analyses identified new fiber quality-related loci and demonstrated that the introgressed alleles from two diploid cottons had a large effect on fiber quality improvement. These loci provided the potential power to overcome the bottleneck in fiber quality improvement. Our study uncovered several critical genomic signatures generated by historical breeding effects in cotton and a wealth of data that enrich genomic resources for the research community.

Heterologous overexpression of the GbTCP5 gene increased root hair length, root hair and stem trichome density, and lignin content in transgenic Arabidopsis.

Teosinte branched1/cycloidea/proliferating cell factor1 (TCP) is a plant-specific protein family member involved in plant growth and development. However, the functions of most members of the cotton TCP family are unknown. In this study, the GbTCP5 gene encodes a sea-island cotton class II TCP CIN subclass transcription factor. The GbTCP5 transcription factor is located in the nucleus, has transcriptional activation activity, and can bind to TCP II cis-acting elements. GbTCP5 was widely expressed in tissues with the highest transcript level in the calyx. GbTCP5 is expressed at different developmental stages of the fiber and has significantly high transcriptional level expression in the fibers at 20, 30 and 35 days post anthesis (DPA). Heterologous overexpression of the GbTCP5 gene increased root hair length, root hair and stem trichome density, and stem lignin content in transgenic Arabidopsis compared to the wild type (WT). GbTCP5 binds the promoters of the GL3, EGL3, CPC, MYB46, LBD30, CesA4, VND7, CCOMT1, and CAD5 genes to upregulate their expression. Moreover, the homologous genes of these genes are expressed in the fibers of different developmental stages of the sea-island cotton fiber. These results indicate that GbTCP5 regulates root hair development and secondary wall formation in Arabidopsis and may be a candidate gene for improving cotton fiber quality.

Resequencing a core collection of upland cotton identifies genomic variation and loci influencing fiber quality and yield.

Upland cotton is the most important natural-fiber crop. The genomic variation of diverse germplasms and alleles underpinning fiber quality and yield should be extensively explored. Here, we resequenced a core collection comprising 419 accessions with 6.55-fold coverage depth and identified approximately 3.66 million SNPs for evaluating the genomic variation. We performed phenotyping across 12 environments and conducted genome-wide association study of 13 fiber-related traits. 7,383 unique SNPs were significantly associated with these traits and were located within or near 4,820 genes; more associated loci were detected for fiber quality than fiber yield, and more fiber genes were detected in the D than the A subgenome. Several previously undescribed causal genes for days to flowering, fiber length, and fiber strength were identified. Phenotypic selection for these traits increased the frequency of elite alleles during domestication and breeding. These results provide targets for molecular selection and genetic manipulation in cotton improvement.