Rapidly mining candidate cotton drought resistance genes based on key indicators of drought resistance.

BACKGROUND: Focusing on key indicators of drought resistance is highly important for quickly mining candidate genes related to drought resistance in cotton. RESULTS: In the present study, drought resistance was identified in drought resistance-related RIL populations during the flowering and boll stages, and multiple traits were evaluated; these traits included three key indicators: plant height (PH), single boll weight (SBW) and transpiration rate (Tr). Based on these three key indicators, three groups of extreme mixing pools were constructed for BSA-seq. Based on the mapping interval of each trait, a total of 6.27 Mb QTL intervals were selected on chromosomes A13 (3.2 Mb), A10 (2.45 Mb) and A07 (0.62 Mb) as the focus of this study. Based on the annotation information and qRT‒PCR analysis, three key genes that may be involved in the drought stress response of cotton were screened: GhF6'H1, Gh3AT1 and GhPER55. qRT‒PCR analysis of parental and extreme germplasm materials revealed that the expression of these genes changed significantly under drought stress. Cotton VIGS experiments verified the important impact of key genes on cotton drought resistance. CONCLUSIONS: This study focused on the key indicators of drought resistance, laying the foundation for the rapid mining of drought-resistant candidate genes in cotton and providing genetic resources for directed molecular breeding of drought resistance in cotton.

Genome-Wide and Expression Pattern Analysis of the DVL Gene Family Reveals GhM_A05G1032 Is Involved in Fuzz Development in G. hirsutum.

DVL is one of the small polypeptides which plays an important role in regulating plant growth and development, tissue differentiation, and organ formation in the process of coping with stress conditions. So far, there has been no comprehensive analysis of the expression profile and function of the cotton DVL gene. According to previous studies, a candidate gene related to the development of fuzz was screened, belonging to the DVL family, and was related to the development of trichomes in Arabidopsis thaliana. However, the comprehensive identification and systematic analysis of DVL in cotton have not been conducted. In this study, we employed bioinformatics approaches to conduct a novel analysis of the structural characteristics, phylogenetic tree, gene structure, expression pattern, evolutionary relationship, and selective pressure of the DVL gene family members in four cotton species. A total of 117 DVL genes were identified, including 39 members in G. hirsutum. Based on the phylogenetic analysis, the DVL protein sequences were categorized into five distinct subfamilies. Additionally, we successfully mapped these genes onto chromosomes and visually represented their gene structure information. Furthermore, we predicted the presence of cis-acting elements in DVL genes in G. hirsutum and characterized the repeat types of DVL genes in the four cotton species. Moreover, we computed the Ka/Ks ratio of homologous genes across the four cotton species and elucidated the selective pressure acting on these homologous genes. In addition, we described the expression patterns of the DVL gene family using RNA-seq data, verified the correlation between GhMDVL3 and fuzz development through VIGS technology, and found that some DVL genes may be involved in resistance to biotic and abiotic stress conditions through qRT-PCR technology. Furthermore, a potential interaction network was constructed by WGCNA, and our findings demonstrated the potential of GhM_A05G1032 to interact with numerous genes, thereby playing a crucial role in regulating fuzz development. This research significantly contributed to the comprehension of DVL genes in upland cotton, thereby establishing a solid basis for future investigations into the functional aspects of DVL genes in cotton.

Glyceraldehyde-3-phosphate dehydrogenase Gh_GAPDH9 is associated with drought resistance in Gossypium hirsutum.

Background: Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) is the central enzyme of glycolysis and plays important regulatory roles in plant growth and development and responses to adverse stress conditions. However, studies on the characteristics and functions of cotton GAPDH family genes are still lacking. Methods: In this study, genome-wide identification of the cotton GAPDH gene family was performed, and the phylogeny, gene structures, promoter progenitors and expression profiles of upland cotton GAPDH gene family members were explored by bioinformatics analysis to highlight potential functions. The functions of GhGAPDH9 in response to drought stress were initially validated based on RNA-seq, qRT‒PCR, VIGS techniques and overexpression laying a foundation for further studies on the functions of GAPDH genes. Results: This study is the first systematic analysis of the cotton GAPDH gene family, which contains a total of 84 GAPDH genes, among which upland cotton contains 27 members. Quantitative, phylogenetic and covariance analyses of the genes revealed that the GAPDH gene family has been conserved during the evolution of cotton. Promoter analysis revealed that most cis-acting elements were related to MeJA and ABA. Based on the identified promoter cis-acting elements and RNA-seq data, it was hypothesized that Gh_GAPDH9, Gh_GAPDH11, Gh_GAPDH19 and Gh_GAPDH21 are involved in the response of cotton to abiotic stress. The expression levels of the Gh_GAPDH9 gene in two drought-resistant and two drought-sensitive materials were analyzed by qRT‒PCR and found to be high early in the treatment period in the drought-resistant material. The silencing of Gh_GAPDH9 based on virus-induced gene silencing (VIGS) technology resulted in significant leaf wilting or whole-plant dieback in silenced plants after drought stress compared to the control. The content of-malondialdehyde (MDA) in cotton leaves was significantly increased, and the content of proline (Pro) and chlorophyll (Chl) was reduced. In addition, the leaf wilting and dryness of transgenic lines under drought stress were lower than those of wild-type Arabidopsis, indicating that Gh_GAPDH9 is a positive regulator of drought resistance. In conclusion, our results demonstrate that GAPDH genes play an important role in the response of cotton to abiotic stresses and provide preliminary validation of the function of the Gh_GAPDH9 gene under drought stress. These findings provide an important theoretical basis for further studies on the function of the Gh_GAPDH9 gene and the molecular mechanism of the drought response in cotton.

Using association mapping and local interval haplotype association analysis to improve the cotton drought stress response.

Drought stress has a serious impact on the growth and development of cotton. To explore the relevant molecular mechanism of the drought stress response in cotton, gene mapping based on the QTL interval mapped by simplified genome BSA-seq of the drought-resistance-related RIL population was performed. A QTL region spanning 2.02 Mb on chromosome D07 was selected, and 201 resource materials were genotyped using 9 KASP markers in the interval. After local interval haplotype association analysis, the overlap of the 110 kb peak region confirmed the reliability of this region, and at the same time, the role of GhGF14-30, the only gene in the overlapping region, was modeled in the response of cotton to drought stress. qRTPCR analysis of the materials and population parents proved that this gene plays a role in the drought stress response in cotton. Virus-induced gene silencing proved the importance of this gene in drought-sensitive materials, and drought-resistance-related marker genes also proved that the GhGF14-30 gene may play an important role in the ABA and SOS signaling pathways. This study provides a basis for mining drought stress response functional genes in cotton and lays the foundation for the molecular mechanism of the GhGF14-30 gene in response to drought stress in cotton.

Determination of aflatoxin B1 by novel nanofiber-packed solid-phase extraction coupled with a high performance liquid chromatography-fluorescence detector.

A novel analytical proposal based on nanofiber-packed solid-phase extraction coupled with high performance liquid chromatography-fluorescence detector (HPLC-FLD) has been successfully developed for determining aflatoxin B1 (AFB1) in foods. Four types of nanofibers, including polystyrene (PS) nanofibers, polypyrrole (PPY) nanofibers, polystyrene-acrylic resin (PS-AR) nanofibers, and polystyrene-polyvinyl pyrrolidone (PS-PVP) nanofibers, were fabricated by electrospinning and utilized to prepare a home-made extraction device. In this study, the factors of different fibers, namely, fiber dosage, pH of extraction solution, type of salt ion, concentration of salt ion, and volume of the eluent were optimized. Under optimized conditions, the method showed good linearity in the range of 0.1-40 ng mL-1 with a correlation coefficient greater than 0.999 and good inter-day accuracy (90.8-112.7% recovery) and precision (1.8-3.6% intra-day RSDs, 2.6% inter-day RSD), and the limit of detection (LOD) was 0.05 ng mL-1. Due to its cost-effective, time-saving, environmentally friendly, and simple performance, it has the potential to be utilized to determine aflatoxins in complicated matrices.

Weighted Gene Co-Expression Network Analysis Reveals Hub Genes for Fuzz Development in Gossypium hirsutum.

Fuzzless Gossypium hirsutum mutants are ideal materials for investigating cotton fiber initiation and development. In this study, we used the fuzzless G. hirsutum mutant Xinluzao 50 FLM as the research material and combined it with other fuzzless materials for verification by RNA sequencing to explore the gene expression patterns and differences between genes in upland cotton during the fuzz period. A gene ontology (GO) enrichment analysis showed that differentially expressed genes (DEGs) were mainly enriched in the metabolic process, microtubule binding, and other pathways. A weighted gene co-expression network analysis (WGCNA) showed that two modules of Xinluzao 50 and Xinluzao 50 FLM and four modules of CSS386 and Sicala V-2 were highly correlated with fuzz. We selected the hub gene with the highest KME value among the six modules and constructed an interaction network. In addition, we selected some genes with high KME values from the six modules that were highly associated with fuzz in the four materials and found 19 common differential genes produced by the four materials. These 19 genes are likely involved in the formation of fuzz in upland cotton. Several hub genes belong to the arabinogalactan protein and GDSL lipase, which play important roles in fiber development. According to the differences in expression level, 4 genes were selected from the 19 genes and tested for their expression level in some fuzzless materials. The modules, hub genes, and common genes identified in this study can provide new insights into the formation of fiber and fuzz, and provide a reference for molecular design breeding for the genetic improvement of cotton fiber.

Quantitative trait locus mapping and identification of candidate genes for resistance to Verticillium wilt in four recombinant inbred line populations of Gossypium hirsutum.

Improving resistance to Verticillium wilt is of great significance for achieving high and stable yields of Upland cotton (Gossypium hirsutum). To deeply understand the genetic basis of cotton resistance to Verticillium wilt, Verticillium wilt-resistant Upland Lumianyan 28 and four Verticillium wilt-susceptible Acala cotton cultivars were used to create four recombinant inbred line (RIL) populations of 469 families through nested hybridization. Phenotypic data collected in five stressful environments were used to select resistant and sensitive lines and create a mixed pool of extreme phenotypes for BSA-seq. A total of 8 QTLs associated with Verticillium wilt resistance were identified on 4 chromosomes, of which qVW-A12-5 was detected simultaneously in the RIL populations and in one of the RIL populations and was identified for the first time. According to the sequence comparison and transcriptome analysis of candidate genes in the QTL interval between parents and pools, 4 genes were identified in the qVW-A12-5 interval. qRT-PCR of parental and phenotypically extreme lines revealed that Gh_CPR30 was induced by and may be a candidate gene for resistance to Verticillium wilt in G. hirsutum. Furthermore, VIGS technology revealed that the disease severity index (DSI) of the Gh_CPR30-silenced plants was significantly higher than that of the control. These results indicate that the Gh_CPR30 gene plays an important role in the resistance of G. hirsutum to Verticillium wilt, and the study provides a molecular basis for analyzing the molecular mechanism underlying G. hirsutum resistance to Verticillium wilt.

Genome-wide analysis of serine carboxypeptidase-like protein (SCPL) family and functional validation of Gh_SCPL42 unchromosome conferring cotton Verticillium der Verticillium wilt stress in Gossypium hirsutum.

BACKGROUND: Serine carboxypeptidase-like protein (SCPL) plays an important role in response to stress in plant. However, our knowledge of the function of the SCPL gene family is limited. RESULTS: In this study, a comprehensive and systematic analysis of SCPL gene family was conducted to explore the phylogeny and evolution of the SCPL gene in Gossypium hirsutum. The phenotype and molecular mechanism of silencing of the Gh_SCPL42 under Verticillium wilt stress was also studied. Our results showed that 96 SCPL genes were observed in genome of G. hirsutum, which distributed on 25 chromosomes and most of them were located in the nucleus. The phylogenetic tree analysis showed that members of SCPL gene family can be divided into three subgroups in G. hirsutum, which are relatively conservative in evolution. SCPL gene has a wide range of tissue expression types in G. hirsutum. Promoter analysis showed that the most cis-acting elements related to MeJA and ABA were contained. Through RNA-seq combined with genotyping, it was found that 11 GhSCPL genes not only had significant expression changes during Verticillium wilt stress but also had differential SNPs in the upstream, downstream, exonic or intronic regions. The expression of these 11 genes in the resistant (Zhongzhimian 2) and susceptible (Junmian 1) materials was further analyzed by qRT-PCR, it was found that 6 genes showed significant expression differences in the two materials. Among them, Gh_SCPL42 has the most obvious expression change. Furthermore, virus-induced gene silencing (VIGS) showed necrosis and yellowing of leaves and significantly higher disease severity index (DSI) and disease severity rate (DSR) values in VIGS plants than in control silenced Gh_SCPL42 plants. Moreover, the expression levels of genes related to the SA and JA pathways were significantly downregulated. These results show that Gh_SCPL42 might improve resistance to Verticillium wilt through the SA and JA pathways in G. hirsutum. CONCLUSION: In conclusion, our findings indicated that Gh_SCPL42 gene plays an important role in resistance to Verticillium wilt in cotton. It was provided an important theoretical basis for further research on the function of SCPL gene family and the molecular mechanism of resistance to Verticillium wilt in cotton.

Revealing Genetic Differences in Fiber Elongation between the Offspring of Sea Island Cotton and Upland Cotton Backcross Populations Based on Transcriptome and Weighted Gene Coexpression Networks.

Fiber length is an important indicator of cotton fiber quality, and the time and rate of cotton fiber cell elongation are key factors in determining the fiber length of mature cotton. To gain insight into the differences in fiber elongation mechanisms in the offspring of backcross populations of Sea Island cotton Xinhai 16 and land cotton Line 9, we selected two groups with significant differences in fiber length (long-fiber group L and short-fiber group S) at different fiber development stages 0, 5, 10 and 15 days post-anthesis (DPA) for transcriptome comparison. A total of 171.74 Gb of clean data was obtained by RNA-seq, and eight genes were randomly selected for qPCR validation. Data analysis identified 6055 differentially expressed genes (DEGs) between two groups of fibers, L and S, in four developmental periods, and gene ontology (GO) term analysis revealed that these DEGs were associated mainly with microtubule driving, reactive oxygen species, plant cell wall biosynthesis, and glycosyl compound hydrolase activity. Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis indicated that plant hormone signaling, mitogen-activated protein kinase (MAPK) signaling, and starch and sucrose metabolism pathways were associated with fiber elongation. Subsequently, a sustained upregulation expression pattern, profile 19, was identified and analyzed using short time-series expression miner (STEM). An analysis of the weighted gene coexpression network module uncovered 21 genes closely related to fiber development, mainly involved in functions such as cell wall relaxation, microtubule formation, and cytoskeletal structure of the cell wall. This study helps to enhance the understanding of the Sea Island-Upland backcross population and identifies key genes for cotton fiber development, and these findings will provide a basis for future research on the molecular mechanisms of fiber length formation in cotton populations.

Genome-wide identification of the DUF668 gene family in cotton and expression profiling analysis of GhDUF668 in Gossypium hirsutum under adverse stress.

BACKGROUND: Domain of unknown function 668 (DUF668) may play a crucial role in the plant growth and developmental response to adverse stress. However, our knowledge of the function of the DUF668 gene family is limited. RESULTS: Our study was conducted based on the DUF668 gene family identified from cotton genome sequencing. Phylogenetic analysis showed that the DUF668 family genes can be classified into four subgroups in cotton. We identified 32 DUF668 genes, which are distributed on 17 chromosomes and most of them located in the nucleus of Gossypium hirsutum. Gene structure and motif analyses revealed that the members of the DUF668 gene family can be clustered in G. hirsutum into two broad groups, which are relatively evolutionarily conserved. Transcriptome data analysis showed that the GhDUF668 genes are differentially expressed in different tissues under various stresses (cold, heat, drought, salt, and Verticillium dahliae), and expression is generally increased in roots and stems. Promoter and expression analyses indicated that Gh_DUF668-05, Gh_DUF668-08, Gh_DUF668-11, Gh_DUF668-23 and Gh_DUF668-28 in G. hirsutum might have evolved resistance to adverse stress. Additionally, qRT-PCR revealed that these 5 genes in four cotton lines, KK1543 (drought resistant), Xinluzao 26 (drought sensitive), Zhongzhimian 2 (disease resistant) and Simian 3 (susceptible), under drought and Verticillium wilt stress were all significantly induced. Roots had the highest expression of these 5 genes before and after the treatment. Among them, the expression levels of Gh_DUF668-08 and Gh_DUF668-23 increased sharply at 6 h and reached a maximum at 12 h under biotic and abiotic stress, which showed that they might be involved in the process of adverse stress resistance in cotton. CONCLUSION: The significant changes in GhDUF668 expression in the roots after adverse stress indicate that GhDUF668 is likely to increase plant resistance to stress. This study provides an important theoretical basis for further research on the function of the DUF668 gene family and the molecular mechanism of adverse stress resistance in cotton.

Effect of long-distance inter-basin water transfer on the bioavailability of Cu for the receiving water.

Long-distance inter-basin water transfer solves the problem of unbalanced water resources in different regions. However, it also changes the natural water chemistry characteristics as well as the bioavailability of heavy metals in the receiving water. In this study, taking the South-to-North Water Transfer Project in China as an example, the basic physicochemical characteristics of the source water (in the Danjiangkou (DJK) Reservoir) and receiving water (in the Beitang (BT) Reservoir) were studied. The BLM (biotic ligand model) was used to study the effect of long-distance inter-basin water transfer on the bioavailability of Cu in receiving waters. The results showed that the TOC (total organic carbon) and TDS (total dissolved solids) in the BT Reservoir water were 10 times and 4.6 times greater than those of the DJK Reservoir water, respectively. The ions in the BT Reservoir were mainly (K++Na+)-(SO42-+CI-), while the ions in the DJK Reservoir were mainly (Ca2++Mg2+)-HCO3-. The results from the BLM showed that the main species of Cu in the water was total organic Cu (Torg Cu), which accounted for 98.69% and 99.77% of the Cu in the DJK Reservoir and BT Reservoir, respectively. The LC50 of Cu for Daphnia magna was 1203.40 ± 57.70 µg/l in the BT Reservoir and only 101.93 ± 7.60 µg/l in the DJK Reservoir. The criteria maximum concentration value of the BT Reservoir was 13.75 times that of the DJK Reservoir, while the criteria continuous concentration value of the BT Reservoir was 13.76 times that of the DJK Reservoir. These results showed that the heavy metals content in water bodies should not be used as the only consideration for water ecological security in the inter-basin water transfer process, and that differences in water quality criteria values caused by differing water environmental qualities in the river basins must be taken into consideration.

Detection of ß-agonists in pork tissue with novel electrospun nanofibers-based solid-phase extraction followed ultra-high performance liquid chromatography/tandem mass spectrometry.

A selective analytical method based on packed-fiber solid-phase extraction and ultra-high performance liquid chromatography-tandem mass spectrometry (PFSPE-UPLC-MS/MS) has been developed for determination of six ß-agonists (clorprenaline, bambuterol, clenbuterol, brombuterol, mabuterol, and penbuterol) in pork tissue. Polystyrene-polymeric crown ether (PS-PCE) composite nanofibers were fabricated by electrospinning and utilized to prepare the homemade extraction columns. With optimal conditions, all analytes were separated very well and the blank pork did not disturb the determination, and the linearity is good in a range of 5.0µg/kg-25.0µg/kg. The recoveries were 79.3-110.1%. RSDs for intra-day were in the range of 1.5-10.5% and RSDs for inter-day were 4.7-11.8%. Above all, only 5mg of sorbent and 200µL of elution solvent were favorable to directly extract all analytes in a complex matrix. The method is simple and cost-effective, and has the potential to be applied to quantitatively analyze the concentrations of polar species in food samples containing complex matrix.

Heavy metals in aquatic organisms of different trophic levels and their potential human health risk in Bohai Bay, China.

Fourteen aquatic organism samples were collected from Bohai Bay, and concentrations of five heavy metals were measured to evaluate the pollution levels in aquatic organisms and the potential risk to human health. The concentrations of Zn and Cu were much higher than those of Cd, Cr, and Pb in all the organisms. In general, the heavy metal concentration levels were in the order phytoplankton < zooplankton < fish < shrimp < shellfish. Heavy metal concentrations in higher trophic-level aquatic organisms in Bohai Bay were compared to those in the organisms from other worldwide coastal waters. The concentration levels of most heavy metals were higher than the 75th percentile, except that Pb concentration was between the 25th and 50th percentiles. The calculated bioconcentration factors (BCF) of Cr, Cu, and Pb for phytoplankton were less than 100, indicating no accumulation in primary producers. The bioaccumulation factor (BAF) of Pb for zooplankton was the highest, indicating significant Pb accumulation in zooplankton. For higher trophic-level aquatic organisms, the order of BAF values was fish < shrimp < shellfish for most metals except for Pb. The human health risk assessment suggests that strict abatement measures of heavy metals must be taken to decrease the health risk caused by consuming aquatic products.

[Comparison and application of biological indices of macroinvertebrates in river health assessment].

The different biological indices usually result in different results in the river health assessment. It is imperative and valuable to identify the correlation among different indices and their applicability for assessing stream health. In this study, totally five biological indices were selected and compared in the investigation of macroinvertebrate communities in the Taizi river. The results showed significant correlations among the five indices. However, due to the difference in health rating criteria for each biological index, different results of health ratings were obtained when different indices were used. The responding sensitivities to disturbance caused by different types of human activities were studied for each index to determine their applicability in assessment of river health. The data indicated that the BI index had significant correlations with land use and dissolved oxygen and was a good indicator for these two types of disturbance. The FBI index could well reflect the acid and ammonia contamination of the investigated stream. Strong negative correlation was found between the ASPT index and several water quality parameters concerning oxygen consumption. The B-IBI index had a significant negative correlation with the total nitrogen concentration, being a good indicator for nitrogen contamination. Besides, the B-IBI index was also significantly correlated to disturbance caused by other types of human activities and can be used as an indicator for both land use and aquatic pollution. To be concluded, the BI index and ASPT index can be individually used to assess the land use of a riverine and the impact of hydrochemical index on the ecosystems, whereas the B-IBI index could be a suitable indicator for evaluating the stream health correlated with various human activities.