Genome-wide identification and expression-pattern analysis of sulfate transporter (SULTR) gene family in cotton under multiple abiotic stresses and fiber development.

Sulfate transporter (SULTR) proteins are in charge of the transport and absorption on sulfate substances, and have been reported to play vital roles in the biological processes of plant growth and stress response. However, there were few reports of genome-wide identification and expression-pattern analysis of SULTRs in Hibiscus mutabilis. Gossypium genus is a ideal model for studying the allopolyploidy, therefore two diploid species (G. raimondii and G. arboreum) and two tetraploid species (G. hirsutum and G. barbadense) were chosen in this study to perform bioinformatic analyses, identifying 18, 18, 35, and 35 SULTR members, respectively. All the 106 cotton SULTR genes were utilized to construct the phylogenetic tree together with 11 Arabidopsis thaliana, 13 Oryza sativa, and 8 Zea mays ones, which was divided into Group1-Group4. The clustering analyses of gene structures and 10 conserved motifs among the cotton SULTR genes showed the consistent evolutionary relationship with the phylogenetic tree, and the results of gene-duplication identification among the four representative Gossypium species indicated that genome-wide or segment duplication might make main contributions to the expansion of SULTR gene family in cotton. Having conducted the cis-regulatory element analysis in promoter region, we noticed that the existing salicylic acid (SA), jasmonic acid (JA), and abscisic acid (ABA) elements could have influences with expression levels of cotton SULTR genes. The expression patterns of GhSULTR genes were also investigated on the 7 different tissues or organs and the developing ovules and fibers, most of which were highly expressed in root, stem, sepal, receptacel, ovule at 10 DPA, and fiber at 20 and 25 DPA. In addition, more active regulatory were observed in GhSULTR genes responding to multiple abiotic stresses, and 12 highly expressed genes showed the similar expression patterns in the quantitative Real-time PCR experiments under cold, heat, salt, and drought treatments. These findings broaden our insight into the evolutionary relationships and expression patterns of the SULTR gene family in cotton, and provide the valuable information for further screening the vital candidate genes on trait improvement.

Genetic linkage analysis of stable QTLs in Gossypium hirsutum RIL population revealed function of GhCesA4 in fiber development.

INTRODUCTION: Upland cotton is an important allotetrapolyploid crop providing natural fibers for textile industry. Under the present high-level breeding and production conditions, further simultaneous improvement of fiber quality and yield is facing unprecedented challenges due to their complex negative correlations. OBJECTIVES: The study was to adequately identify quantitative trait loci (QTLs) and dissect how they orchestrate the formation of fiber quality and yield. METHODS: A high-density genetic map (HDGM) based on an intraspecific recombinant inbred line (RIL) population consisting of 231 individuals was used to identify QTLs and QTL clusters of fiber quality and yield traits. The weighted gene correlation network analysis (WGCNA) package in R software was utilized to identify WGCNA network and hub genes related to fiber development. Gene functions were verified via virus-induced gene silencing (VIGS) and clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 strategies. RESULTS: An HDGM consisting of 8045 markers was constructed spanning 4943.01 cM of cotton genome. A total of 295 QTLs were identified based on multi-environmental phenotypes. Among 139 stable QTLs, including 35 newly identified ones, seventy five were of fiber quality and 64 yield traits. A total of 33 QTL clusters harboring 74 QTLs were identified. Eleven candidate hub genes were identified via WGCNA using genes in all stable QTLs and QTL clusters. The relative expression profiles of these hub genes revealed their correlations with fiber development. VIGS and CRISPR/Cas9 edition revealed that the hub gene cellulose synthase 4 (GhCesA4, GH_D07G2262) positively regulate fiber length and fiber strength formation and negatively lint percentage. CONCLUSION: Multiple analyses demonstrate that the hub genes harbored in the QTLs orchestrate the fiber development. The hub gene GhCesA4 has opposite pleiotropic effects in regulating trait formation of fiber quality and yield. The results facilitate understanding the genetic basis of negative correlation between cotton fiber quality and yield.

Genome-Wide Identification and Functional Analysis of RF2 Gene Family and the Critical Role of GhRF2-32 in Response to Drought Stress in Cotton.

Cotton is an important natural fiber crop. The RF2 gene family is a member of the bZIP transcription factor superfamily, which plays an important role in plant resistance to environmental stresses. In this paper, the RF2 gene family of four cotton species was analyzed genome-wide, and the key gene RF2-32 was cloned for functional verification. A total of 113 RF2 genes were identified in the four cotton species, and the RF2 family was relatively conserved during the evolution of cotton. Chromosome mapping and collinear analysis indicated that fragment replication was the main expansion mode of RF2 gene family during evolution. Cis-element analysis showed that there were many elements related to light response, hormone response and abiotic stress response in the promoters of RF2 genes. The transcriptome and qRT-PCR analysis of RF2 family genes in upland cotton showed that RF2 family genes responded to salt stress and drought stress. GhRF2-32 protein was localized in the cell nucleus. Silencing the GhRF2-32 gene showed less leaf wilting and increased total antioxidant capacity under drought and salt stress, decreased malondialdehyde content and increased drought and salt tolerance. This study revealed the evolutionary and functional diversity of the RF2 gene family, which laid a foundation for the further study of stress-resistant genes in cotton.

Genome-wide chromatin accessibility landscape and dynamics of transcription factor networks during ovule and fiber development in cotton.

BACKGROUND: The development of cotton fiber is regulated by the orchestrated binding of regulatory proteins to cis-regulatory elements associated with developmental genes. The cis-trans regulatory dynamics occurred throughout the course of cotton fiber development are elusive. Here we generated genome-wide high-resolution DNase I hypersensitive sites (DHSs) maps to understand the regulatory mechanisms of cotton ovule and fiber development. RESULTS: We generated DNase I hypersensitive site (DHS) profiles from cotton ovules at 0 and 3 days post anthesis (DPA) and fibers at 8, 12, 15, and 18 DPA. We obtained a total of 1185 million reads and identified a total of 199,351 DHSs through ~ 30% unique mapping reads. It should be noted that more than half of DNase-seq reads mapped multiple genome locations and were not analyzed in order to achieve a high specificity of peak profile and to avoid bias from repetitive genomic regions. Distinct chromatin accessibilities were observed in the ovules (0 and 3 DPA) compared to the fiber elongation stages (8, 12, 15, and 18 DPA). Besides, the chromatin accessibility during ovules was particularly elevated in genomic regions enriched with transposable elements (TEs) and genes in TE-enriched regions were involved in ovule cell division. We analyzed cis-regulatory modules and revealed the influence of hormones on fiber development from the regulatory divergence of transcription factor (TF) motifs. Finally, we constructed a reliable regulatory network of TFs related to ovule and fiber development based on chromatin accessibility and gene co-expression network. From this network, we discovered a novel TF, WRKY46, which may shape fiber development by regulating the lignin content. CONCLUSIONS: Our results not only reveal the contribution of TEs in fiber development, but also predict and validate the TFs related to fiber development, which will benefit the research of cotton fiber molecular breeding.

Genome-wide characterization of trichome birefringence-like genes provides insights into fiber yield improvement.

Cotton is an important fiber crop. The cotton fiber is an extremely long trichome that develops from the epidermis of an ovule. The trichome is a general and multi-function plant organ, and trichome birefringence-like (TBL) genes are related to trichome development. At the genome-wide scale, we identified TBLs in four cotton species, comprising two cultivated tetraploids (Gossypium hirsutum and G. barbadense) and two ancestral diploids (G. arboreum and G. raimondii). Phylogenetic analysis showed that the TBL genes clustered into six groups. We focused on GH_D02G1759 in group IV because it was located in a lint percentage-related quantitative trait locus. In addition, we used transcriptome profiling to characterize the role of TBLs in group IV in fiber development. The overexpression of GH_D02G1759 in Arabidopsis thaliana resulted in more trichomes on the stems, thereby confirming its function in fiber development. Moreover, the potential interaction network was constructed based on the co-expression network, and it was found that GH_D02G1759 may interact with several genes to regulate fiber development. These findings expand our knowledge of TBL family members and provide new insights for cotton molecular breeding.

Large-scale metabolome analysis reveals dynamic changes of metabolites during foxtail millet grain filling.

Compared with traditional staple crops, foxtail millet grain is rich in nutrition and beneficial to human health. Foxtail millet is also tolerance to various abiotic stresses, including drought, making it a good plant for growing in barren land. The study on the composition of metabolites and its dynamics changes during grain development is helpful to understand the process of foxtail millet grain formation. In our study, metabolic and transcriptional analysis were used to uncover the metabolic processes that could influence grain filling in foxtail millet. A total of 2104 known metabolites, belonging to 14 categories, were identified during grain filling. Functional analysis of DAMs and DEGs revealed a stage-specific metabolic properties in foxtail millet grain filling. Some important metabolic processes, such as flavonoid biosynthesis, glutathione metabolism, linoleic acid metabolism, starch and sucrose metabolism and valine, leucine and isoleucine biosynthesis were co-mapped for DEGs and DAMs. Thus, we constructed a gene-metabolite regulatory network of these metabolic pathways to explain their potential functions during grain filling. Our study showed the important metabolic processes during grain filling and focused on the dynamic changes of related metabolites and genes at different stages, which provided a reference for us to better understand and improve foxtail millet grain development and yield.

Fine mapping and candidate gene analysis of qFL-A12-5: a fiber length-related QTL introgressed from Gossypium barbadense into Gossypium hirsutum.

KEY MESSAGE: The fiber length-related qFL-A12-5 identified in CSSLs introgressed from Gossypium barbadense into Gossypium hirsutum was fine-mapped to an 18.8 kb region on chromosome A12, leading to the identification of the GhTPR gene as a potential regulator of cotton fiber length. Fiber length is a key determinant of fiber quality in cotton, and it is a key target of artificial selection for breeding and domestication. Although many fiber length-related quantitative trait loci have been identified, there are few reports on their fine mapping or candidate gene validation, thus hampering efforts to understand the mechanistic basis of cotton fiber development. Our previous study identified the qFL-A12-5 associated with superior fiber quality on chromosome A12 in the chromosome segment substitution line (CSSL) MBI7747 (BC4F3:5). A single segment substitution line (CSSL-106) screened from BC6F2 was backcrossed to construct a larger segregation population with its recurrent parent CCRI45, thus enabling the fine mapping of 2852 BC7F2 individuals using denser simple sequence repeat markers to narrow the qFL-A12-5 to an 18.8 kb region of the genome, in which six annotated genes were identified in Gossypium hirsutum. Quantitative real-time PCR and comparative analyses led to the identification of GH_A12G2192 (GhTPR) encoding a tetratricopeptide repeat-like superfamily protein as a promising candidate gene for qFL-A12-5. A comparative analysis of the protein-coding regions of GhTPR among Hai1, MBI7747, and CCRI45 revealed two non-synonymous mutations. The overexpression of GhTPR resulted in longer roots in Arabidopsis, suggesting that GhTPR may regulate cotton fiber development. These results provide a foundation for future efforts to improve cotton fiber length.

Genome-Wide Identification and Characterization of the PPO Gene Family in Cotton (Gossypium) and Their Expression Variations Responding to Verticillium Wilt Infection.

Polyphenol oxidases (PPOs) are copper-binding metalloproteinases encoded by nuclear genes, ubiquitously existing in the plastids of microorganisms, plants, and animals. As one of the important defense enzymes, PPOs have been reported to participate in the resistant processes that respond to diseases and insect pests in multiple plant species. However, PPO gene identification and characterization in cotton and their expression patterns under Verticillium wilt (VW) treatment have not been clearly studied. In this study, 7, 8, 14, and 16 PPO genes were separately identified from Gossypium arboreum, G. raimondii, G. hirsutum, and G. barbadense, respectively, which were distributed within 23 chromosomes, though mainly gathered in chromosome 6. The phylogenetic tree manifested that all the PPOs from four cotton species and 14 other plants were divided into seven groups, and the analyses of the conserved motifs and nucleotide sequences showed highly similar characteristics of the gene structure and domains in the cotton PPO genes. The dramatically expressed differences were observed among the different organs at various stages of growth and development or under the diverse stresses referred to in the published RNA-seq data. Quantitative real-time PCR (qRT-PCR) experiments were also performed on the GhPPO genes in the roots, stems, and leaves of VW-resistant MBI8255 and VW-susceptible CCRI36 infected with Verticillium dahliae V991, proving the strong correlation between PPO activity and VW resistance. A comprehensive analysis conducted on cotton PPO genes contributes to the screening of the candidate genes for subsequent biological function studies, which is also of great significance for the in-depth understanding of the molecular genetic basis of cotton resistance to VW.

Genome-wide artificial introgressions of Gossypium barbadense into G. hirsutum reveal superior loci for simultaneous improvement of cotton fiber quality and yield traits.

INTRODUCTION: The simultaneous improvement of fiber quality and yield for cotton is strongly limited by the narrow genetic backgrounds of Gossypium hirsutum (Gh) and the negative genetic correlations among traits. An effective way to overcome the bottlenecks is to introgress the favorable alleles of Gossypium barbadense (Gb) for fiber quality into Gh with high yield. OBJECTIVES: This study was to identify superior loci for the improvement of fiber quality and yield. METHODS: Two sets of chromosome segment substitution lines (CSSLs) were generated by crossing Hai1 (Gb, donor-parent) with cultivar CCRI36 (Gh) and CCRI45 (Gh) as genetic backgrounds, and cultivated in 6 and 8 environments, respectively. The kmer genotyping strategy was improved and applied to the population genetic analysis of 743 genomic sequencing data. A progeny segregating population was constructed to validate genetic effects of the candidate loci. RESULTS: A total of 68,912 and 83,352 genome-wide introgressed kmers were identified in the CCRI36 and CCRI45 populations, respectively. Over 90 % introgressions were homologous exchanges and about 21 % were reverse insertions. In total, 291 major introgressed segments were identified with stable genetic effects, of which 66(22.98 %), 64(21.99 %), 35(12.03 %), 31(10.65 %) and 18(6.19 %) were beneficial for the improvement of fiber length (FL), strength (FS), micronaire, lint-percentage (LP) and boll-weight, respectively. Thirty-nine introgression segments were detected with stable favorable additive effects for simultaneous improvement of 2 or more traits in Gh genetic background, including 6 could increase FL/FS and LP. The pyramiding effects of 3 pleiotropic segments (A07:C45Clu-081, D06:C45Clu-218, D02:C45Clu-193) were further validated in the segregating population. CONCLUSION: The combining of genome-wide introgressions and kmer genotyping strategy showed significant advantages in exploring genetic resources. Through the genome-wide comprehensive mining, a total of 11 clusters (segments) were discovered for the stable simultaneous improvement of FL/FS and LP, which should be paid more attention in the future.

Evolutionary divergence of duplicated genomes in newly described allotetraploid cottons.

Allotetraploid cotton (Gossypium) species represents a model system for the study of plant polyploidy, molecular evolution, and domestication. Here, chromosome-scale genome sequences were obtained and assembled for two recently described wild species of tetraploid cotton, Gossypium ekmanianum [(AD)6, Ge] and Gossypium stephensii [(AD)7, Gs], and one early form of domesticated Gossypium hirsutum, race punctatum [(AD)1, Ghp]. Based on phylogenomic analysis, we provide a dated whole-genome level perspective for the evolution of the tetraploid Gossypium clade and resolved the evolutionary relationships of Gs, Ge, and domesticated G. hirsutum. We describe genomic structural variation that arose during Gossypium evolution and describe its correlates-including phenotypic differentiation, genetic isolation, and genetic convergence-that contributed to cotton biodiversity and cotton domestication. Presence/absence variation is prominent in causing cotton genomic structural variations. A presence/absence variation-derived gene encoding a phosphopeptide-binding protein is implicated in increasing fiber length during cotton domestication. The relatively unimproved Ghp offers the potential for gene discovery related to adaptation to environmental challenges. Expanded gene families enoyl-CoA δ isomerase 3 and RAP2-7 may have contributed to abiotic stress tolerance, possibly by targeting plant hormone-associated biochemical pathways. Our results generate a genomic context for a better understanding of cotton evolution and for agriculture.

Linkage and association analyses reveal that hub genes in energy-flow and lipid biosynthesis pathways form a cluster in upland cotton.

Upland cotton is an important allotetraploid crop that provides both natural fiber for the textile industry and edible vegetable oil for the food or feed industry. To better understand the genetic mechanism that regulates the biosynthesis of storage oil in cottonseed, we identified the genes harbored in the major quantitative trait loci/nucleotides (QTLs/QTNs) of kernel oil content (KOC) in cottonseed via both multiple linkage analyses and genome-wide association studies (GWAS). In 'CCRI70' RILs, six stable QTLs were simultaneously identified by linkage analysis of CHIP and SLAF-seq strategies. In '0-153' RILs, eight stable QTLs were detected by consensus linkage analysis integrating multiple strategies. In the natural panel, thirteen and eight loci were associated across multiple environments with two algorithms of GWAS. Within the confidence interval of a major common QTL on chromosome 3, six genes were identified as participating in the interaction network highly correlated with cottonseed KOC. Further observations of gene differential expression showed that four of the genes, LtnD, PGK, LPLAT1, and PAH2, formed hub genes and two of them, FER and RAV1, formed the key genes in the interaction network. Sequence variations in the coding regions of LtnD, FER, PGK, LPLAT1, and PAH2 genes may support their regulatory effects on oil accumulation in mature cottonseed. Taken together, clustering of the hub genes in the lipid biosynthesis interaction network provides new insights to understanding the mechanism of fatty acid biosynthesis and TAG assembly and to further genetic improvement projects for the KOC in cottonseeds.

Protoplast Dissociation and Transcriptome Analysis Provides Insights to Salt Stress Response in Cotton.

As one of the pioneer crops widely planted in saline-alkaline areas, Gossypium provides daily necessities, including natural fiber, vegetable proteins, and edible oils. However, cotton fiber yield and quality are highly influenced by salt stress. Therefore, elucidating the molecular mechanisms of cotton in response to salinity stress is importance to breed new cultivars with high tolerance. In this study, we first developed a method for single-cell RNA-seq based on isolating protoplast from cotton root tips; then, we studied the impact of salinity stress on gene expression profiling and their dynamic changes using the developed high-efficiency method for protoplast dissociation suitable for single-cell RNA-seq. A total of 3391 and 2826 differentially expressed genes (DEGs) were identified in salt-treated samples before and after protoplast dissociation, respectively, which were enriched into several molecular components, including response to stimulus, response to stress, and cellular macromolecule metabolic process by gene ontology (GO) analysis. Plant hormone signal transduction, phenylpropanoid biosynthesis, and MAPK signaling pathway were found to be enriched via Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Twenty-two and nine salinity-responsive DEGs participated in plant hormone signaling and MAPK signaling in roots, before and after protoplast dissociation, respectively; six upregulated DEGs were involved in ABA signaling transduction, namely, Ga04G2111, Ga07G0142, Ga09G2061, Ga10G0262, Ga01G0063, and Ga08G1915 which indicates their potential functions on plants adapting to salt stress. Additionally, 384 and 257 transcription factors (TFs) were differentially expressed in salt-stress roots before and after protoplast dissociation, respectively, of which significantly up-regulated TFs mainly belonged to the AP2/ERF-ERF family, which implied their potential roles responding to salt stress. These results not only provide novel insights to reveal the regulatory networks in plant's root response to salt stress, but also lay the solid foundation for further exploration on cellular heterogeneity by single-cell transcriptome sequencing.

Multi-environment Evaluations Across Ecological Regions Reveal That the Kernel Oil Content of Cottonseed Is Equally Determined by Genotype and Environment.

Cotton is the fifth-largest oil crop in the world. A high kernel oil content (KOC) and high stability are important cottonseed attributes for food security. In this study, the phenotype of KOC and the genotype-by-environment interaction factors were collectively dissected using 250 recombinant inbred lines, their parental cultivars sGK156 and 901-001, and CCRI70 across multi-environments. ANOVA and correlation analysis showed that both genotype and environment contributed significantly to KOC accumulation. Analyses of additive main effect multiplicative interaction and genotype-by-environment interaction biplot models presented the effects of genotype, environment, and genotype by environment on KOC performance and the stability of the experimental materials. Interaction network analysis revealed that meteorological and geographical factors explained 38% of the total KOC variance, with average daily rainfall contributing the largest positive impact and cumulative rainfall having the largest negative impact on KOC accumulation. This study provides insight into KOC accumulation and could direct selection strategies for improved KOC and field management of cottonseed in the future.

QTL mapping and candidate gene prediction for fiber yield and quality traits in a high-generation cotton chromosome substitution line with Gossypium barbadense segments.

Gossypium provides the foremost natural fiber for supporting the rapid development of the textile industry. Quantitative trait locus (QTL) mapping of fiber yield and quality traits is, thus, of great significance for providing a foundation for the genetic improvement of key target traits in cotton production. In this study, a superior chromosome segment substitution line (CSSL), MBI8255, with high yield and premium fiber quality characteristics was cultivated from the BC5F3:5 lineage derived from G. barbadense Hai1 and G. hirsutum CCRI36, and was chosen to construct a segregation population containing 123 F2 individuals with CCRI36. A total of 71 polymorphic SSR (simple sequence repeat) markers were identified based on a previous high-density linkage map, and 17 QTLs distributed on five chromosomes were detected, of which 10 QTLs for cotton yield explained 0.26-15.41% of phenotypic variations, while 7 QTLs for fiber quality explained 0.84-9.38% of phenotypic variations, separately containing four and one stable QTLs detected from over two environments. Among three identified QTL clusters, only the Chr19 QTL cluster harbored two stable and one unstable QTL for three different traits, and hence this significant region, which included 1546 genes, was subjected to functional enrichment and transcriptome expression analyses, ultimately screening eight candidate genes relevant to fiber development. This study not only provides useful information for the further fine-mapping and functional verification of candidate genes, but also offers a solid foundation for revealing the molecular mechanisms of fiber formation.

Genome-wide analysis of the CalS gene family in cotton reveals their potential roles in fiber development and responses to stress.

Callose deposition occurs during plant growth and development, as well as when plants are under biotic and abiotic stress. Callose synthase is a key enzyme for the synthesis of callose. In this study, 27, 28, 16, and 15 callose synthase family members were identified in Gossypium hirsutum, Gossypium barbadense, Gossypium raimondii, and Gossypium arboreum using the sequence of Arabidopsis callose synthase. The CalSs were divided into five groups by phylogenetic, gene structure, and conservative motif analysis. The conserved motifs and gene structures of CalSs in each group were highly similar. Based on the analysis of cis-acting elements, it is inferred that GhCalSs were regulated by abiotic stress. WGD/Segmental duplication promoted the amplification of the CalS gene in cotton, and purification selection had an important function in the CalS family. The transcriptome data and qRT-PCR under cold, heat, salt, and PEG treatments showed that GhCalSs were involved in abiotic stress. The expression patterns of GhCalSs were different in various tissues. We predicted that GhCalS4, which was highly expressed in fibers, had an important effect on fiber elongation. Hence, these results help us understand the role of GhCalSs in fiber development and stress response.

Identification of Candidate Cotton Genes Associated With Fiber Length Through Quantitative Trait Loci Mapping and RNA-Sequencing Using a Chromosome Segment Substitution Line.

Fiber length is an important determinant of fiber quality, and it is a quantitative multi-genic trait. Identifying genes associated with fiber length is of great importance for efforts to improve fiber quality in the context of cotton breeding. Integrating transcriptomic information and details regarding candidate gene regions can aid in candidate gene identification. In the present study, the CCRI45 line and a chromosome segment substitution line (CSSL) with a significantly higher fiber length (MBI7747) were utilized to establish F2 and F2:3 populations. Using a high-density genetic map published previously, six quantitative trait loci (QTLs) associated with fiber length and two QTLs associated with fiber strength were identified on four chromosomes. Within these QTLs, qFL-A07-1, qFL-A12-2, qFL-A12-5, and qFL-D02-1 were identified in two or three environments and confirmed by a meta-analysis. By integrating transcriptomic data from the two parental lines and through qPCR analyses, four genes associated with these QTLs including Cellulose synthase-like protein D3 (CSLD3, GH_A12G2259 for qFL-A12-2), expansin-A1 (EXPA1, GH_A12G1972 for qFL-A12-5), plasmodesmata callose-binding protein 3 (PDCB3, GH_A12G2014 for qFL-A12-5), and Polygalacturonase (At1g48100, GH_D02G0616 for qFL-D02-1) were identified as promising candidate genes associated with fiber length. Overall, these results offer a robust foundation for further studies regarding the molecular basis for fiber length and for efforts to improve cotton fiber quality.

DNA Methylation and RNA-Sequencing Analysis Show Epigenetic Function During Grain Filling in Foxtail Millet (Setaria italica L.).

Grain filling is a crucial process for crop yield and quality. Certain studies already gained insight into the molecular mechanism of grain filling. However, it is unclear whether epigenetic modifications are associated with grain filling in foxtail millet. Global DNA methylation and transcriptome analysis were conducted in foxtail millet spikelets during different stages to interpret the epigenetic effects of the grain filling process. The study employed the whole-genome bisulfite deep sequencing and advanced bioinformatics to sequence and identify all DNA methylation during foxtail millet grain filling; the DNA methylation-mediated gene expression profiles and their involved gene network and biological pathway were systematically studied. One context of DNA methylation, namely, CHH methylation, was accounted for the largest percentage, and it was gradually increased during grain filling. Among all developmental stages, the methylation levels were lowest at T2, followed by T4, which mainly occurred in CHG. The distribution of differentially methylated regions (DMR) was varied in the different genetic regions for three contexts. In addition, gene expression was negatively associated with DNA methylation. Evaluation of the interconnection of the DNA methylome and transcriptome identified some stage-specific differentially expressed genes associated with the DMR at different stages compared with the T1 developmental stage, indicating the potential function of epigenetics on the expression regulation of genes related to the specific pathway at different stages of grain development. The results demonstrated that the dynamic change of DNA methylation plays a crucial function in gene regulation, revealing the potential function of epigenetics in grain development in foxtail millet.

Genome-wide identification of the MIOX gene family and their expression profile in cotton development and response to abiotic stress.

The enzyme myo-inositol oxygenase (MIOX) catalyzes the myo-inositol into glucuronic acid. In this study, 6 MIOX genes were identified from all of the three diploid cotton species (Gossypium arboretum, Gossypium herbaceum and Gossypium raimondii) and Gossypioides kirkii, 12 MIOX genes were identified from two domesticated tetraploid cottons Gossypium hirsutum, Gossypium barbadense, and 11 MIOX genes were identified from three wild tetraploid cottons Gossypium tomentosum, Gossypium mustelinum and Gossypium darwinii. The number of MIOX genes in tetraploid cotton genome is roughly twice that of diploid cotton genome. Members of MIOX family were classified into six groups based on the phylogenetic analysis. Integrated analysis of collinearity events and chromosome locations suggested that both whole genome duplication and segmental duplication events contributed to the expansion of MIOX genes during cotton evolution. The ratios of non-synonymous (Ka) and synonymous (Ks) substitution rates revealed that purifying selection was the main force driving the evolution of MIOX genes. Numerous cis-acting elements related to light responsive element, defense and stress responsive element were identified in the promoter of the MIOX genes. Expression analyses of MIOX genes based on RNA-seq data and quantitative real time PCR showed that MIOX genes within the same group shared similar expression patterns with each other. All of these results provide the foundation for further study of the biological functions of MIOX genes in cotton environmental adaptability.

Genome-wide quantitative trait loci mapping on Verticillium wilt resistance in 300 chromosome segment substitution lines from Gossypium hirsutum × Gossypium barbadense.

Cotton Verticillium wilt (VW) is a devastating disease seriously affecting fiber yield and quality, and the most effective and economical prevention measure at present is selection and extension of Gossypium varieties harboring high resistance to VW. However, multiple attempts to improve the VW resistance of the most widely cultivated upland cottons have made little significant progress. The introduction of chromosome segment substitution lines (CSSLs) provide the practical solutions for merging the superior genes related with high yield and wide adaptation from Gossypium hirsutum and VW resistance and the excellent fiber quality from Gossypium barbadense. In this study, 300 CSSLs were chosen from the developed BC5F3:5 CSSLs constructed from CCRI36 (G. hirsutum) and Hai1 (G. barbadense) to conduct quantitative trait locus (QTL) mapping of VW resistance, and a total of 40 QTL relevant to VW disease index (DI) were identified. Phenotypic data were obtained from a 2-year investigation in two fields with two replications per year. All the QTL were distributed on 21 chromosomes, with phenotypic variation of 1.05%-10.52%, and 21 stable QTL were consistent in at least two environments. Based on a meta-analysis, 34 novel QTL were identified, while 6 loci were consistent with previously identified QTL. Meanwhile, 70 QTL hotspot regions were detected, including 44 novel regions. This study concentrates on QTL identification and screening for hotspot regions related with VW in the 300 CSSLs, and the results lay a solid foundation not only for revealing the genetic and molecular mechanisms of VW resistance but also for further fine mapping, gene cloning and molecular designing in breeding programs for resistant cotton varieties.

Micro-FTIR combined with curve fitting method to study cellulose crystallinity of developing cotton fibers.

This study aimed to use micro-FTIR with transmission mode to investigate cellulose crystallinity of developing cotton fibers. Compared with ATR-FTIR method, we found that micro-FTIR can obtain more information of cellulose inside of the developing cotton fibers, especially in high wavenumber of 2800-3000 cm-1 region. Combined with curve fitting method, a new IR crystallinity index (CI) method named wax crystallinity index (WCI) was introduced to evaluate the cellulose crystallinity in the development of cotton fibers based on the peak and area ratios of 2900 cm-1/2850 cm-1 and 2900 cm-1/2920 cm-1. The obtained WCI values demonstrated an excellent coefficient of determination with X-ray diffraction (XRD) CI method with the value up to 0.99. This study suggested that micro-FTIR was an effective technique to qualitatively analyze the crystallinity in developing cotton fibers combined with curve fitting method.