Analyzing the Diversity of MYB Family Response Strategies to Drought Stress in Different Flax Varieties Based on Transcriptome Data.

The MYB transcription factor family has numerous members, and is involved in biological activities, such as ABA signaling, which plays an important role in a plant's resistance to abiotic stresses such as drought. However, the diversity of MYB members that respond to drought stress and their regulatory mechanisms in different flax varieties were unclear. In this study, we obtained 855.69 Gb of clean data from 120 flax root samples from 20 flax (Linum usitatissimum L.) varieties, assembled 92,861 transcripts, and identified 434 MYB family members in each variety. The expression profiles of the MYB transcription factor family from 20 flax varieties under drought stress were analyzed. The results indicated that there are four strategies by which the MYB family responds to drought stress in these 20 flax varieties, each of which has its own specific processes, such as development, reproduction, and localization processes. The four strategies also include common biological processes, such as stimulus responses, metabolic processes, and biological regulation. The WGCNA method was subsequently employed to identify key members of the MYB family involved in response strategies to drought stress. The results demonstrated that a 1R-MYB subfamily gene co-expression network is significantly related to the gibberellin response and cytokinin-activated signaling pathway processes in the 'Strategy 4' for MYB family response to drought, identifying core genes such as Lus.scaffold70.240. Our results showed a diversity of MYB family responses to drought stress within flax varieties, and these results contribute to deciphering the mechanisms of the MYB family regulation of drought resistance. This will promote the more accurate breeding development of flax to adapt to agricultural production under drought conditions.

A Multi-Omics Analysis Revealed the Diversity of the MYB Transcription Factor Family's Evolution and Drought Resistance Pathways.

The MYB transcription factor family can regulate biological processes such as ABA signal transduction to cope with drought stress, but its evolutionary mechanism and the diverse pathways of response to drought stress in different species are rarely reported. In this study, a total of 4791 MYB family members were identified in 908,757 amino acid sequences from 12 model plants or crops using bioinformatics methods. It was observed that the number of MYB family members had a linear relationship with the chromosome ploidy of species. A phylogenetic analysis showed that the MYB family members evolved in subfamily clusters. In response to drought stress, the pathways of MYB transcription factor families exhibited species-specific diversity, with closely related species demonstrating a higher resemblance. This study provides abundant references for drought resistance research and the breeding of wheat, soybean, and other plants.

Character variation of root space microbial community composition in the response of drought-tolerant spring wheat to drought stress.

Drought is the most prevalent environmental stress in crop production, posing a significant danger to food security. Microorganisms in the crop root zone affect crop growth and development, enhance effective nutrient use, and resist adversity hazards. To analyze the changes and functional differences of root space microbial (endosphere-rhizosphere-bulk soil) communities in spring wheat under drought stress. In this study, the root, rhizosphere, and bulk soil of the drought-tolerant group (DTG, three varieties) and drought-sensitive group (DSG, three varieties) were collected. The control (CK, 25-28%), moderate drought (MD, 15-18%), and severe drought (SD, 9-12%) were analyzed by high-throughput sequencing and bioinformatics. The results showed significant differences in the diversity of Bacteria and Fungi in the root space of spring wheat under drought stress (P < 0.05), with the drought-tolerant group exhibiting higher microbial diversity. The microbial community change in spring wheat root space was mainly determined by the niche differentiation of endosphere, rhizosphere, and bulk soil and declined from endosphere to bulk soil due to drought. The antagonism between microbial and root-space species increased, and the community's complexity and stability deteriorated. Enriching drought-resistant preference groups like Actinobaciota, Variovorax, Streptomyces, and Conocybe altered the structure and function of the microbial community in the root space of spring wheat. Spring wheat's root space Bacteria and Fungi have different strategies to respond to drought.

Drought Sensitivity of Spring Wheat Cultivars Shapes Rhizosphere Microbial Community Patterns in Response to Drought.

Drought is the most important natural disaster affecting crop growth and development. Crop rhizosphere microorganisms can affect crop growth and development, enhance the effective utilization of nutrients, and resist adversity and hazards. In this paper, six spring wheat varieties were used as research material in the dry farming area of the western foot of the Greater Khingan Mountains, and two kinds of water control treatments were carried out: dry shed rain prevention (DT) and regulated water replenishment (CK). Phenotypic traits, including physiological and biochemical indices, drought resistance gene expression, soil enzyme activity, soil nutrient content, and the responses of potential functional bacteria and fungi under drought stress, were systematically analyzed. The results showed that compared with the control (CK), the leaf wilting, drooping, and yellowing of six spring wheat varieties were enhanced under drought (DT) treatment. The plant height, fresh weight (FW), dry weight (DW), net photosynthetic rate (Pn) and stomatal conductance (Gs), soil total nitrogen (TN), microbial biomass carbon (MBC), microbial biomass nitrogen (MBN), microbial biomass phosphorus (MBP), organic carbon (SOC), and soil alkaline phosphatase (S-ALP) contents were significantly decreased, among which, FW, Gs and MBC decreased by more than 7.84%, 17.43% and 11.31%, respectively. By contrast, the soil total phosphorus (TP), total potassium (TK), and soil catalase (S-CAT) contents were significantly increased (p < 0.05). TaWdreb2 and TaBADHb genes were highly expressed in T.D40, T.L36, and T.L33 and were expressed at low levels in T.N2, T.B12, and T.F5. Among them, the relative expression of the TaWdreb2 gene in T.L36 was significantly increased by 2.683 times compared with CK. Soil TN and TP are the most sensitive to drought stress and can be used as the characteristic values of drought stress. Based on this, a drought-tolerant variety (T.L36) and a drought-sensitive variety (T.B12) were selected to further analyze the changes in rhizosphere microorganisms. Drought treatment and cultivar differences significantly affected the composition of the rhizosphere microbial community. Drought caused a decrease in the complexity of the rhizosphere microbial network, and the structure of bacteria was more complex than that of fungi. The Shannon index and network modular number of bacteria in these varieties (T.L36) increased, with rich small-world network properties. Actinobacteria, Chloroflexi, Firmicutes, Basidiomycota, and Ascomycota were the dominant bacteria under drought treatment. The beneficial bacteria Bacillus, Penicillium, and Blastococcus were enriched in the rhizosphere of T.L36. Brevibacillus and Glycomyce were enriched in the rhizosphere of T.B12. In general, drought can inhibit the growth and development of spring wheat, and spring wheat can resist drought hazards by regulating the expression of drought-related genes, regulating physiological metabolites, and enriching beneficial microorganisms.

Screening and genome-wide analysis of lignocellulose-degrading bacteria from humic soil.

Crop straw contains huge amounts of exploitable energy, and efficient biomass degradation measures have attracted worldwide attention. Mining strains with high yields of cellulose-degrading enzymes is of great significance for developing clean energy and industrial production of related enzymes. In this study, we reported a high-quality genome sequence of Bacillus velezensis SSF6 strain using high-throughput sequencing technology (Illumina PE150 and PacBio) and assessed its lignocellulose degradation potential. The results demonstrated that the genome of B. velezensis SSF6 was 3.89 Mb and contained 4,015 genes, of which 2,972, 3,831 and 158 genes were annotated in the COGs (Clusters of Orthologous Groups), KEGG (Kyoto Encyclopedia of Genes and Genomes) and CAZyme (Carbohydrate-Active enZymes) databases, respectively, and contained a large number of genes related to carbohydrate metabolism. Furthermore, B. velezensis SSF6 has a high cellulose degradation capacity, with a filter paper assay (FPA) and an exoglucanase activity of 64.48 ± 0.28 and 78.59 ± 0.42 U/mL, respectively. Comparative genomic analysis depicted that B. velezensis SSF6 was richer in carbohydrate hydrolase gene. In conclusion, the cellulose-degrading ability of B. velezensis SSF6 was revealed by genome sequencing and the determination of cellulase activity, which laid a foundation for further cellulose degradation and bioconversion.

High-Quality Genome Assembly and Genome-Wide Association Study of Male Sterility Provide Resources for Flax Improvement.

Flax is an economic crop with a long history. It is grown worldwide and is mainly used for edible oil, industry, and textiles. Here, we reported a high-quality genome assembly for "Neiya No. 9", a popular variety widely grown in China. Combining PacBio long reads, Hi-C sequencing, and a genetic map reported previously, a genome assembly of 473.55 Mb was constructed, which covers ~94.7% of the flax genome. These sequences were anchored onto 15 chromosomes. The N50 lengths of the contig and scaffold were 0.91 Mb and 31.72 Mb, respectively. A total of 32,786 protein-coding genes were annotated, and 95.9% of complete BUSCOs were found. Through morphological and cytological observation, the male sterility of flax was considered dominant nuclear sterility. Through GWAS analysis, the gene LUSG00017705 (cysteine synthase gene) was found to be closest to the most significant SNP, and the expression level of this gene was significantly lower in male sterile plants than in fertile plants. Among the significant SNPs identified in the GWAS analysis, only two were located in the coding region, and these two SNPs caused changes in the protein encoded by LUSG00017565 (cysteine protease gene). It was speculated that these two genes may be related to male sterility in flax. This is the first time the molecular mechanism of male sterility in flax has been reported. The high-quality genome assembly and the male sterility genes revealed, provided a solid foundation for flax breeding.

Dynamic changes of soil microorganisms in rotation farmland at the western foot of the Greater Khingan range.

Crop rotation and other tillage systems can affect soil microbial communities and functions. Few studies have reported the response of soil spatial microbial communities to rotation under drought stress. Therefore, the purpose of our study was to explore the dynamic changes of the soil space microbial community under different drought stress-rotation patterns. In this study, two water treatments were set up, control W1 (mass water content 25%-28%), and drought W2 (mass water content 9%-12%). Four crop rotation patterns were set in each water content, spring wheat continuous (R1), spring wheat-potato (R2), spring wheat-potato-rape (R3) and spring wheat-rape (R4), for a total of eight treatments (W1R1, W1R2, W1R3, W1R4, W2R1, W2R2, W2R3, W2R4). Endosphere, rhizosphere and bulk soil of spring wheat in each treatment were collected, and root space microbial community data were generated. The soil microbial community changed under different treatments and their relationship with soil factors were analyzed using a co-occurrence network, mantel test, and other methods. The results revealed that the alpha diversity of microorganisms in the rhizosphere and bulk soil did not differ significantly, but it was significantly greater than in the endosphere. The bacteria community structure was more stable, fungi alpha-diversity significant changes (p < 0.05), that were more sensitive to the response of various treatments than bacteria. The co-occurrence network between fungal species was stable under rotation patterns (R2, R3, R4), while the community stability was poor under continuous cropping pattern (R1), and interactions were strengthened. Soil organic matter (SOM), microbial biomass carbon (MBC), and pH value were the most important factors dominating the bacteria community structural changed in the endosphere, rhizosphere, and bulk soil. The dominant factor that affected the fungal community structural changed in the endosphere, rhizosphere, and bulk soil was SOM. Therefore, we conclude that soil microbial community changes under the drought stress-rotation patterns are mainly influenced by soil SOM and microbial biomass content.

Ubiquitous conservative interaction patterns between post-spliced introns and their mRNAs revealed by genome-wide interspecies comparison.

Introns, as important vectors of biological functions, can influence many stages of mRNA metabolism. However, in recent research, post-spliced introns are rarely considered. In this study, the optimal matched regions between introns and their mRNAs in nine model organism genomes were investigated with improved Smith-Waterman local alignment software. Our results showed that the distributions of mRNA optimal matched frequencies were highly consistent or universal. There are optimal matched frequency peaks in the UTR regions, which are obvious, especially in the 3'-UTR. The matched frequencies are relatively low in the CDS regions of the mRNA. The distributions of the optimal matched frequencies around the functional sites are also remarkably changed. The centers of the GC content distributions for different sequences are different. The matched rate distributions are highly consistent and are located mainly between 60% and 80%. The most probable value of the optimal matched segments is about 20 bp for lower eukaryotes and 30 bp for higher eukaryotes. These results show that there are abundant functional units in the introns, and these functional units are correlated structurally with all kinds of sequences of mRNA. The interaction between the post-spliced introns and their corresponding mRNAs may play a key role in gene expression.

A novel approach to analyze the association characteristics between post-spliced introns and their corresponding mRNA.

Studies have shown that post-spliced introns promote cell survival when nutrients are scarce, and intron loss/gain can influence many stages of mRNA metabolism. However, few approaches are currently available to study the correlation between intron sequences and their corresponding mature mRNA sequences. Here, based on the results of the improved Smith-Waterman local alignment-based algorithm method (SW method) and binding free energy weighted local alignment algorithm method (BFE method), the optimal matched segments between introns and their corresponding mature mRNAs in Caenorhabditis elegans (C.elegans) and their relative matching frequency (RF) distributions were obtained. The results showed that although the distributions of relative matching frequencies on mRNAs obtained by the BFE method were similar to the SW method, the interaction intensity in 5'and 3'untranslated regions (UTRs) regions was weaker than the SW method. The RF distributions in the exon-exon junction regions were comparable, the effects of long and short introns on mRNA and on the five functional sites with BFE method were similar to the SW method. However, the interaction intensity in 5'and 3'UTR regions with BFE method was weaker than with SW method. Although the matching rate and length distribution shape of the optimal matched fragment were consistent with the SW method, an increase in length was observed. The matching rates and the length of the optimal matched fragments were mainly in the range of 60%-80% and 20-30bp, respectively. Although we found that there were still matching preferences in the 5'and 3'UTR regions of the mRNAs with BFE, the matching intensities were significantly lower than the matching intensities between introns and their corresponding mRNAs with SW method. Overall, our findings suggest that the interaction between introns and mRNAs results from synergism among different types of sequences during the evolutionary process.

Chromosome-scale Genome Assembly of the Yellow Nutsedge (Cyperus esculentus).

The yellow nutsedge (Cyperus esculentus L. 1753) is an unconventional oil plant with oil-rich tubers, and a potential alternative for traditional oil crops. Here, we reported the first high-quality and chromosome-level genome assembly of the yellow nutsedge generated by combining PacBio HiFi long reads, Novaseq short reads, and Hi-C data. The final genome size is 225.6 Mb with an N50 of 4.3 Mb. More than 222.9 Mb scaffolds were anchored to 54 pseudochromosomes with a BUSCO score of 96.0%. We identified 76.5 Mb (33.9%) repetitive sequences across the genome. A total of 23,613 protein-coding genes were predicted in this genome, of which 22,847 (96.8%) were functionally annotated. A whole-genome duplication event was found after the divergence of Carex littledalei and Rhynchospora breviuscula, indicating the rich genetic resources of this species for adaptive evolution. Several significantly enriched GO terms were related to invasiveness of the yellow nutsedge, which may explain its plastic adaptability. In addition, several enriched Kyoto Encyclopedia of Genes and Genomes pathways and expanded gene families were closely related with substances in tubers, partially explaining the genomic basis of characteristics of this oil-rich tuber.

Changes in the Microbial Community in Maize (Zea mays L.) Root Spatial Structure Following Short-Term Nitrogen Application.

The beneficial interactions between crop roots and microbiomes play a key role in crop nutrient availability, growth promotion, and disease suppression. Recent research, however, rarely reported the effects of nitrogen (N) application rate on microbial community composition at different spatial structures in the maize root zone. Therefore, one experiment was conducted to examine the influence of three N-application levels (0, 180, and 360 kg N ha-1) on microbial community composition in three root-associated compartments of maize (bulk soil, rhizoplane, and endosphere). The microbial diversity and community composition differed significantly among the various compartments. The effects of N application on fungal composition decreased in the order bulk soil > rhizosphere > endosphere at different sampling positions. Also, the fungal composition was more sensitive to the N-fertilizer rate in the bulk soil and the rhizosphere than the bacterial community. A total of 14.42, 9.46, and 3.55% of all taxonomic groups were sensitive to N fertilizer, respectively. The keystone species fungal groups were Humicola (bulk soil), Gibberella (rhizosphere soil), and Humicola (endosphere). Together, our results demonstrate that compared with that of the bacterial community, the fungal community composition was more susceptible to different N-application rates. N fertilization affected the distribution of microflora by changing soil physicochemical properties and enzyme activities. There were strong correlations between microbial communities in maize under the N180 treatment. Moreover, the N180 treatment had the maximum fresh yield and biomass at 64.5 and 24.3 kg·ha-1, respectively.

Chemical Constituents of Cyperus esculentus Leaves and the Protective Effect against Agricultural Fungicide-Induced Hepatotoxicity.

Cyperus esculentus is cultivated as a crop plant due to its edible and oily tubers (tiger nut). However, little is known about the phytochemicals and bioeffects of the leaves. This study was conducted to identify and quantify the chemical constituents of C. esculentus leaves and evaluate their bioactivities. By liquid chromatography-mass spectrometry, 30 compounds including flavan-3-ols, caffeic acid derivatives, and flavones, were identified from the leaves. The quantitative analysis revealed that gallocatechin (8), procyanidin B1 (15), catechin (16), chlorogenic acid (19), orientin (30), and luteolin 7-O-glucuronide (31) are the major chemical constituents of C. esculentus leaves. The contents of these six chemical constituents in the leaves collected in September in Hohhot, China, reached to 1460.85±7.66, 10178.77±302.65, 1048.35±17.37, 1722.15±26.13, 5318.62±277.16, and 1526.54±11.95 µg, respectively, in one gram of the dried leaves. The leaf extract (CELE) showed strong antioxidant activity in vitro, with compounds 8, 15, and 19 contributing the most. CELE showed significant protection against the agricultural fungicide tebuconazole-induced developmental toxicity and hepatotoxicity in zebrafish.

Quantification and purification of procyanidin B1 from food byproducts.

The flavan-3-ols of 10 primarily plant food byproducts, including Muscat Hamburg grape seed, hawthorn sarcocarp, litchi pericarp, cocoa bean, peanut skin, lotus seedpod, Xinyang Maojian green tea, Cinnamomi cortex, Sargentodoxa cuneata stem, and Cyperus esculentus, leaves were analyzed. Ultrahigh-performance liquid chromatography-triple quadrupole mass spectrometry was used for the analysis. Cyperus esculentus leaves contained a high amount of procyanidin B1 (198.9 mg/100 g), second only to Muscat Hamburg grape seed (292.9 mg/100 g). Unlike grape seed that contained several procyanidin B isomers with very similar retention times, C. esculentus leaves contained primarily procyanidin B1 with few isomers. Procyanidin B1 was enriched in the ethyl acetate fraction of a 70% methanol extract of C. esculentus leaves and purified at 95% purity by two runs of open column chromatography. Direct chromatography of the plant extract on octadecylsilane and Sephadex LH20 open columns improved the yield of the resultant leaf procyanidin B1 (95% purity) to 0.21‰. The present research demonstrated that the leaves of C. esculentus, byproducts of tigernut, are ideal plant sources for isolating and providing high-purity procyanidin B1. PRACTICAL APPLICATION: Procyanidin B1 has a broad range of health benefits. Cyperus esculentus is cultivated in many countries with nearly 6190 square hectares (hm2 ) in the Spanish Mediterranean region in 2020-2021 and over 16,700 hm2 in China in recent years, primarily for its tubers. The byproducts, the leaves of C. esculentus, contain high levels of procyanidin B1, with few isomers that interfere with its isolation and purification. Thus, the leaves of this plant provide a viable source for preparing high-purity procyanidin B1.

Genomic regions controlling yield-related traits in spring wheat: a mini review and a case study for rainfed environments in Australia and China.

A genome-wide association study (GWAS) was performed in six environments to identify major or consistent alleles responsible for wheat yield traits in Australia and North China where rainfed farming system is adopted. A panel of 228 spring wheat varieties were genotyped by double digest restriction-site associated DNA genotyping-by-sequencing. A total of 223 significant marker-trait association (MTAs) and 46 candidate genes for large- or consistent-effect MTAs were identified. The results were compared with previous studies based on a mini-review of 23 GWAS analyses on wheat yield. A phenomenon seldom reported in previous studies was that MTAs responsible for the trait tended to cluster together at certain chromosome segments, and many candidate genes were in the form of gene clusters. Although linkage disequilibrium (LD) might contribute to the co-segregation of the regions, it also suggested that marker-assisted selection (MAS) or transgenic method targeting a single gene might not be as effective as MAS targeting a larger genomic region where all the genes or gene clusters underlying play important roles.

Global insights to drought stress perturbed genes in oat (Avena sativa L.) seedlings using RNA sequencing.

Oat (Avena sativa L.) is an important crop in northwestern China. Drought stress is the most significant factor affecting oat yield. In the present study, we explored the changes that occur in oats under drought stress conditions at a global genomic level. RNA sequencing was performed using 15-day-old oat seedlings. The differentially expressed transcripts were identified, and their related functions and pathways were investigated. In total, 1,065 unigenes were differentially expressed in oats under drought stress conditions. Of these, 386 unigenes were upregulated and 679 were downregulated. The perturbed transcripts were closely related to the biosynthesis of secondary metabolites, plant hormone signal transduction, and biosynthesis of antibiotics. DN50483_c0_g1_i3, which was annotated as acetyl-CoA carboxylase, was a significant node in the protein-protein interaction network. Biosynthesis of antibiotics and secondary metabolites may be involved in the drought stress response mechanisms of oats. The perturbed transcripts may provide targets for improving plant stress responses.

Major genomic regions responsible for wheat yield and its components as revealed by meta-QTL and genotype-phenotype association analyses.

MAIN CONCLUSION: Meta-QTL (MQTL) analysis was done for yield-related traits in wheat. Candidate genes were identified within the refined MQTL and further validated by genotype-phenotype association analysis. Extensive studies have been undertaken on quantitative trait locus/loci (QTL) for wheat yield and its component traits. This study conducted a meta-analysis of 381 QTL related to wheat yield under various environments, including irrigated, drought- and/or heat-stressed conditions. Markers flanking meta-QTL (MQTL) were mapped on the wheat reference genome for their physical positions. Putative candidate genes were examined for MQTL with a physical interval of less than 20 Mbp. A total of 86 MQTL were identified as responsible for yield, of which 34 were for irrigated environments, 39 for drought-stressed environments, 36 for heat-stressed environments, and 23 for both drought- and heat-stressed environments. The high-confidence genes within the physical positions of the MQTL flanking markers were screened in the reference genome RefSeq V1.0, which identified 210 putative candidate genes. The phenotypic data for 14 contrasting genotypes with either high or low yield performance-according to the Australian National Variety Trials-were associated with their genotypic data obtained through ddRAD sequencing, which validated 18 genes or gene clusters associated with MQTL that had important roles for wheat yield. The detected and refined MQTL and candidate genes will be useful for marker-assisted selection of high yield in wheat breeding.

Effects of Trace Irrigation at Different Depths on Transcriptome Expression Pattern in Cotton (G. hirsutum L.) Leaves.

Drought is a limiting factor for cotton productivity and quality. Irrigation could increase cotton yield. This study is aimed at formulating a proper irrigation depth for cotton at China' Inner Mongolia and at investigating the molecular mechanism underlying the difference induced by irrigation. Transcriptomic analysis was carried out to reveal the global transcriptome profiles on the leaves of cotton seedlings (G. hirsutum L. cv. "Zhongmian 92") with trace irrigation tapes at 30 cm (D30) and 50 cm (D50) underground. The differentially expressed genes (DEGs) were identified and clustered by functional enrichment analysis. The results showed that no significant differences were found in the lint percentage. The yields of unpinned and lint cotton were increased by the D30 regime but decreased by the D50 regime. Transcriptomic analysis showed that 4,549 nonoverlapped DEGs were identified by comparative analysis. Transcription factors, including bZIP, WARK, Myb, and NAC, were altered between D50 and D30. The D50 regime induced more DEGs compared with the D30 regime, which was associated with plant tolerance to abiotic stresses and drought. In conclusion, trace irrigation at 30 cm underground was suitable for cotton irrigation at China's Inner Mongolia, while the D50 irrigation regime influenced the cotton yield via drought stress in cotton plants.

Nitrogen Fertilization Increases Root Growth and Coordinates the Root-Shoot Relationship in Cotton.

The root system plays an important role in the growth and development of cotton, and root growth is closely related to shoot growth, both of which are affected by N availability in the soil. However, it is unknown how N affects root growth and the root-shoot relationship under various N rates in the Yellow River Basin, China. Thus, the aim of this study was to assess the impacts of the application rate of N on root growth and the root-shoot relationship, to provide insight into the N regulation of root and shoot growth and N efficiency from the perspective of the root system. A field experiment conducted in 2014 and 2015 was used to determine the effects of N rates (0, 120, 240, and 480 kg ha-1) on root morphology, root distribution, the root-shoot relationship, and cotton yield. A moderate N fertilization rate (240 kg ha-1) increased root length, root surface area, and root biomass in most soil layers and significantly increased total root growth and total root biomass by more than 36.06% compared to the 0 kg ha-1 treatment. In addition, roots in the surface soil layers were more strongly affected by N fertilization than roots distributed in the deeper soil layers. Total root length, total root surface area, and root biomass in the 0-15 cm layer were significantly correlated with shoot biomass and boll biomass. In the 60-75 cm layer, total root length, total root surface area, and root length were significantly positively correlated with seed cotton yield. The application of a moderate level of N markedly increased total shoot biomass, boll biomass, and seed cotton yield. Our results show that increased shoot and boll biomasses were correlated with a significant increase in the root system especially the shallow roots in the moderate N treatment (240 kg ha-1), leading to an increase in cotton seed yield.