Fingerprint Finder: Identifying Genomic Fingerprint Sites in Cotton Cohorts for Genetic Analysis and Breeding Advancement.

Genomic data in Gossypium provide numerous data resources for the cotton genomics community. However, to fill the gap between genomic analysis and breeding field work, detecting the featured genomic items of a subset cohort is essential for geneticists. We developed FPFinder v1.0 software to identify a subset of the cohort's fingerprint genomic sites. The FPFinder was developed based on the term frequency-inverse document frequency algorithm. With the short-read sequencing of an elite cotton pedigree, we identified 453 pedigree fingerprint genomic sites and found that these pedigree-featured sites had a role in cotton development. In addition, we applied FPFinder to evaluate the geographical bias of fiber-length-related genomic sites from a modern cotton cohort consisting of 410 accessions. Enriching elite sites in cultivars from the Yangtze River region resulted in the longer fiber length of Yangze River-sourced accessions. Apart from characterizing functional sites, we also identified 12,536 region-specific genomic sites. Combining the transcriptome data of multiple tissues and samples under various abiotic stresses, we found that several region-specific sites contributed to environmental adaptation. In this research, FPFinder revealed the role of the cotton pedigree fingerprint and region-specific sites in cotton development and environmental adaptation, respectively. The FPFinder can be applied broadly in other crops and contribute to genetic breeding in the future.

Cotton pedigree genome reveals restriction of cultivar-driven strategy in cotton breeding.

BACKGROUND: Many elite genes have been identified from the available cotton genomic data, providing various genetic resources for gene-driven breeding. However, backbone cultivar-driven breeding is the most widely applied strategy. Revealing the genetic basis of cultivar-driven strategy's restriction is crucial for transition of cotton breeding strategy. RESULT: CRI12 is a backbone cultivar in cultivar-driven breeding. Here we sequence the pedigree of CRI12 using Nanopore long-read sequencing. We construct a graphical pedigree genome using the high-quality CRI12 genome and 13,138 structural variations within 20 different pedigree members. We find that low hereditary stability of elite segments in backbone cultivars is a drawback of cultivar-driven strategy. We also identify 623 functional segments in CRI12 for multiple agronomic traits in presence and absence variation-based genome-wide association study on three cohorts. We demonstrate that 25 deleterious segments are responsible for the geographical divergence of cotton in pathogen resistance. We also characterize an elite pathogen-resistant gene (GhKHCP) utilized in modern cotton breeding. In addition, we identify 386 pedigree fingerprint segments by comparing the segments of the CRI12 pedigree with those of a large cotton population. CONCLUSION: We characterize the genetic patterns of functional segments in the pedigree of CRI12 using graphical genome method, revealing restrictions of cultivar-driven strategies in cotton breeding. These findings provide theoretical support for transitioning from cultivar-driven to gene-driven strategy in cotton breeding.

Identification and molecular evolution of the GLX genes in 21 plant species: a focus on the Gossypium hirsutum.

BACKGROUND: The glyoxalase system includes glyoxalase I (GLXI), glyoxalase II (GLXII) and glyoxalase III (GLXIII), which are responsible for methylglyoxal (MG) detoxification and involved in abiotic stress responses such as drought, salinity and heavy metal. RESULTS: In this study, a total of 620 GLX family genes were identified from 21 different plant species. The results of evolutionary analysis showed that GLX genes exist in all species from lower plants to higher plants, inferring that GLX genes might be important for plants, and GLXI and GLXII account for the majority. In addition, motif showed an expanding trend in the process of evolution. The analysis of cis-acting elements in 21 different plant species showed that the promoter region of the GLX genes were rich in phytohormones and biotic and abiotic stress-related elements, indicating that GLX genes can participate in a variety of life processes. In cotton, GLXs could be divided into two groups and most GLXIs distributed in group I, GLXIIs and GLXIIIs mainly belonged to group II, indicating that there are more similarities between GLXII and GLXIII in cotton evolution. The transcriptome data analysis and quantitative real-time PCR analysis (qRT-PCR) show that some members of GLX family would respond to high temperature treatment in G.hirsutum. The protein interaction network of GLXs in G.hirsutum implied that most members can participate in various life processes through protein interactions. CONCLUSIONS: The results elucidated the evolutionary history of GLX family genes in plants and lay the foundation for their functions analysis in cotton.

Comparative Transcriptome Analysis Revealed Key Genes Regulating Gossypol Synthesis in Tetraploid Cultivated Cotton.

Tetraploid cultivated cotton (Gossypium spp.) produces cottonseeds rich in protein and oil. Gossypol and related terpenoids, stored in the pigment glands of cottonseeds, are toxic to human beings and monogastric animals. However, a comprehensive understanding of the genetic basis of gossypol and gland formation is still lacking. We performed a comprehensive transcriptome analysis of four glanded versus two glandless tetraploid cultivars distributed in Gossypium hirsutum and Gossypium barbadense. A weighted gene co-expression network analysis (WGCNA) based on 431 common differentially expressed genes (DEGs) uncovered a candidate module that was strongly associated with the reduction in or disappearance of gossypol and pigment glands. Further, the co-expression network helped us to focus on 29 hub genes, which played key roles in the regulation of related genes in the candidate module. The present study contributes to our understanding of the genetic basis of gossypol and gland formation and serves as a rich potential source for breeding cotton cultivars with gossypol-rich plants and gossypol-free cottonseed, which is beneficial for improving food safety, environmental protection, and economic gains of tetraploid cultivated cotton.

Functional characterization of the GhNRT2.1e gene reveals its significant role in improving nitrogen use efficiency in Gossypium hirsutum.

Background: Nitrate is the primary type of nitrogen available to plants, which is absorbed and transported by nitrate transporter 2 (NRT2) at low nitrate conditions. Methods: Genome-wide identification of NRT2 genes in G. hirsutum was performed. Gene expression patterns were revealed using RNA-seq and qRT-PCR. Gene functions were characterized using overexpression in A. thaliana and silencing in G. hirsutum. Protein interactions were verified by yeast two-hybrid and luciferase complementation imaging (LCI) assays. Results: We identified 14, 14, seven, and seven NRT2 proteins in G. hirsutum, G. barbadense, G. raimondii, and G. arboreum. Most NRT2 proteins were predicted in the plasma membrane. The NRT2 genes were classified into four distinct groups through evolutionary relationships, with members of the same group similar in conserved motifs and gene structure. The promoter regions of NRT2 genes included many elements related to growth regulation, phytohormones, and abiotic stresses. Tissue expression pattern results revealed that most GhNRT2 genes were specifically expressed in roots. Under low nitrate conditions, GhNRT2 genes exhibited different expression levels, with GhNRT2.1e being the most up-regulated. Arabidopsis plants overexpressing GhNRT2.1e exhibited increased biomass, nitrogen and nitrate accumulation, nitrogen uptake and utilization efficiency, nitrogen-metabolizing enzyme activity, and amino acid content under low nitrate conditions. In addition, GhNRT2.1e-silenced plants exhibited suppressed nitrate uptake and accumulation, hampered plant growth, affected nitrogen metabolism processes, and reduced tolerance to low nitrate. The results showed that GhNRT2.1e could promote nitrate uptake and transport under low nitrate conditions, thus effectively increasing nitrogen use efficiency (NUE). We found that GhNRT2.1e interacts with GhNAR2.1 by yeast two-hybrid and LCI assays. Discussion: Our research lays the foundation to increase NUE and cultivate new cotton varieties with efficient nitrogen use.

Identification of MYB gene family and functional analysis of GhMYB4 in cotton (Gossypium spp.).

Myeloblastosis (MYB) transcription factors (TFs) form a large gene family involved in a variety of biological processes in plants. Little is known about their roles in the development of cotton pigment glands. In this study, 646 MYB members were identified in Gossypium hirsutum genome and phylogenetic classification was analyzed. Evolution analysis revealed assymetric evolution of GhMYBs during polyploidization and sequence divergence of MYBs in G. hirustum was preferentially happend in D sub-genome. WGCNA (weighted gene co-expression network analysis) showed that four modules had potential relationship with gland development or gossypol biosynthesis in cotton. Eight differentially expressed GhMYB genes were identified by screening transcriptome data of three pairs of glanded and glandless cotton lines. Of these, four were selected as candidate genes for cotton pigment gland formation or gossypol biosynthesis by qRT-PCR assay. Silencing of GH_A11G1361 (GhMYB4) downregulated expression of multiple genes in gossypol biosynthesis pathway, indicating it could be involved in gossypol biosynthesis. The potential protein interaction network suggests that several MYBs may have indirect interaction with GhMYC2-like, a key regulator of pigment gland formation. Our study was the systematic analysis of MYB genes in cotton pigment gland development, providing candidate genes for further study on the roles of cotton MYB genes in pigment gland formation, gossypol biosynthesis and future crop plant improvement.

Genome-wide analysis elucidates the roles of GhHMA genes in different abiotic stresses and fiber development in upland cotton.

The heavy metal-binding domain is involved in heavy metal transporting and plays a significant role in plant detoxification. However, the functions of HMAs are less well known in cotton. In this study, a total of 143 GhHMAs (heavy metal-binding domain) were detected by genome-wide identification in G. hirsutum L. All the GhHMAs were classified into four groups via phylogenetic analysis. The exon/intron structure and protein motifs indicated that each branch of the GhHMA genes was highly conserved. 212 paralogous GhHMA gene pairs were identified, and the segmental duplications were the main role to the expansion of GhHMAs. The Ka/Ks values suggested that the GhHMA gene family has undergone purifying selection during the long-term evolutionary process. GhHMA3 and GhHMA75 were located in the plasma membrane, while GhHMA26, GhHMA117 and GhHMA121 were located in the nucleus, respectively. Transcriptomic data and qRT-PCR showed that GhHMA26 exhibited different expression patterns in each tissue and during fiber development or under different abiotic stresses. Overexpressing GhHMA26 significantly promoted the elongation of leaf trichomes and also improved the tolerance to salt stress. Therefore, GhHMA26 may positively regulate fiber elongation and abiotic stress. Yeast two-hybrid assays indicated that GhHMA26 and GhHMA75 participated in multiple biological functions. Our results suggest some genes in the GhHMAs might be associated with fiber development and the abiotic stress response, which could promote further research involving functional analysis of GhHMA genes in cotton.

Infection of SARS-CoV-2 causes severe pathological changes in mouse testis.

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has affected more than 600 million people worldwide. Several organs including lung, intestine, and brain are infected by SARS-CoV-2. It has been reported that SARS-CoV-2 receptor angiotensin-converting enzyme-2 (ACE2) is expressed in human testis. However, whether testis is also affected by SARS-CoV-2 is still unclear. In this study, we generate a human ACE2 (hACE2) transgenic mouse model in which the expression of hACE2 gene is regulated by hACE2 promoter. Sertoli and Leydig cells from hACE2 transgenic mice can be infected by SARS-CoV-2 pseudovirus in vitro, and severe pathological changes are observed after injecting the SARS-CoV-2 pseudovirus into the seminiferous tubules. Further studies reveal that Sertoli and Leydig cells from hACE2 transgenic mice are also infected by authentic SARS-CoV-2 virus in vitro. After testis interstitium injection, authentic SARS-CoV-2 viruses are first disseminated to the interstitial cells, and then detected inside the seminiferous tubules which in turn cause germ cell loss and disruption of seminiferous tubules. Our study demonstrates that testis is most likely a target of SARS-CoV-2 virus. Attention should be paid to the reproductive function in SARS-CoV-2 patients.

Development of a Lung Vacancy Mouse Model through CRISPR/Cas9-Mediated Deletion of Thyroid Transcription Factor 1 Exon 2.

A developmental niche vacancy in host embryos is necessary for stem cell complementation-based organ regeneration (SCOG). Thyroid transcription factor 1 (TTF-1) is a tissue-specific transcription factor that regulates the embryonic development and differentiation of the thyroid and, more importantly, lungs; thus, it has been considered as a master gene to knockout in order to develop a lung vacancy host. TTF-1 knockout mice were originally produced by inserting a stop codon in Exon 3 of the gene (E3stop) through embryonic stem cell-based homologous recombination. The main problems of utilizing E3stop host embryos for lung SCOG are that these animals all have a tracheoesophageal fistula (TEF), which cannot be corrected by donor stem cells, and most of them have monolateral sac-like lungs. To improve the mouse model towards achieving SCOG-based lung generation, in this project, we used the CRISPR/Cas9 tool to remove Exon 2 of the gene by zygote microinjection and successfully produced TTF-1 knockout (E2del) mice. Similar to E3stop, E2del mice are birth-lethal due to retarded lung development with sac-like lungs and only a rudimentary bronchial tree, increased basal cells but without alveolar type II cells and blood vessels, and abnormal thyroid development. Unlike E3stop, 57% of the E2del embryos presented type I tracheal agenesis (TA, a kind of human congenital malformation) with a shortened trachea and clear separations of the trachea and esophagus, while the remaining 43% had TEF. Furthermore, all the E2del mice had bilateral sac-like lungs. Both TA and bilateral sac-like lungs are preferred in SCOG. Our work presents a new strategy for producing SCOG host embryos that may be useful for lung regeneration.

Assessment of the predictive ability of the Johns Hopkins Fall Risk Assessment Tool (Chinese Version) in inpatient settings.

AIMS: This study was to assess the predictive ability of the Johns Hopkins Fall Risk Assessment Tool (Chinese Version) in inpatient settings. DESIGN: A case-control study. METHODS: This study was conducted in a tertiary hospital based on 2019 data. With a case-control design in a 1:2 ratio, the predictive ability of the Johns Hopkins Fall Risk Assessment Tool (Chinese Version) was determined by ROC curve. The best cut point was identified based on sensitivity, specificity, positive predict value and negative predict value. Conditional logistical regression analysis was conducted to test the predictive ability of each indicator. RESULTS: The study included 309 patients, with 103 in the case group and 206 in the control groups. Generally, the predictive ability was acceptable with the area under ROC curve value at 0.73 (95% CI: 0.67-0.79). Positive predict value and negative predict value performed best at the cut point of 13. Sensitivity at cut point 6 was much higher than that at cut point 13, though specificity was lower. Except for age, all indicators in the Johns Hopkins Fall Risk Assessment Tool (Chinese Version) demonstrated significant predictive ability as to occurrence of fall. CONCLUSION: The Johns Hopkins Fall Risk Assessment Tool (Chinese Version) is a reliable assessment instrument in the inpatient settings. IMPACT: This is the first study that evaluated the predictive ability of the Johns Hopkins Fall Risk Assessment Tool (Chinese version) in the inpatient settings, and proved that the instrument is reliable for assessing inpatient fall risks. Further studies could be carried out to assess the predict ability of Johns Hopkins Fall Risk Assessment Tool (Chinese version) among specific populations.

Interpretation of convolutional neural networks reveals crucial sequence features involving in transcription during fiber development.

BACKGROUND: Upland cotton provides the most natural fiber in the world. During fiber development, the quality and yield of fiber were influenced by gene transcription. Revealing sequence features related to transcription has a profound impact on cotton molecular breeding. We applied convolutional neural networks to predict gene expression status based on the sequences of gene transcription start regions. After that, a gradient-based interpretation and an N-adjusted kernel transformation were implemented to extract sequence features contributing to transcription. RESULTS: Our models had approximate 80% accuracies, and the area under the receiver operating characteristic curve reached over 0.85. Gradient-based interpretation revealed 5' untranslated region contributed to gene transcription. Furthermore, 6 DOF binding motifs and 4 transcription activator binding motifs were obtained by N-adjusted kernel-motif transformation from models in three developmental stages. Apart from 10 general motifs, 3 DOF5.1 genes were also detected. In silico analysis about these motifs' binding proteins implied their potential functions in fiber formation. Besides, we also found some novel motifs in plants as important sequence features for transcription. CONCLUSIONS: In conclusion, the N-adjusted kernel transformation method could interpret convolutional neural networks and reveal important sequence features related to transcription during fiber development. Potential functions of motifs interpreted from convolutional neural networks could be validated by further wet-lab experiments and applied in cotton molecular breeding.

Identification of the CesA Subfamily and Functional Analysis of GhMCesA35 in Gossypium hirsutum L.

The cellulose synthase genes control the biosynthesis of cellulose in plants. Nonetheless, the gene family members of CesA have not been identified in the newly assembled genome of Gossypiumhirsutum (AD1, HEBAU_NDM8). We identified 38 CesA genes in G. hirsutum (NDM8) and found that the protein sequence of GhMCesA35 is 100% identical to CelA1 in a previous study. It is already known that CelA1 is involved in cellulose biosynthesis in vitro. However, the function of this gene in vivo has not been validated. In this study, we verified the function of GhMCesA35 in vivo based on overexpressed Arabidopsis thaliana. In addition, we found that it interacted with GhCesA7 through the yeast two-hybrid assay. This study provides new insights for studying the biological functions of CesA genes in G. hirsutum, thereby improving cotton fiber quality and yield.

Genome-wide identification and expression analysis of GL2-interacting-repressor (GIR) genes during cotton fiber and fuzz development.

MAIN CONCLUSION: GL2-interacting-repressor (GIR) family members may contribute to fiber/fuzz formation via a newly discovered unique pathway in Gossypium arboreum. There are similarities between cotton fiber development and the formation of trichomes and root hairs. The GL2-interacting-repressors (GIRs) are crucial regulators of root hair and trichome formation. The GaFzl gene, annotated as GaGIR1, is negatively associated with trichome development and fuzz initiation. However, there is relatively little available information regarding the other GIR genes in cotton, especially regarding their effects on cotton fiber development. In this study, 21 GIR family genes were identified in the diploid cotton species Gossypium arboreum; these genes were divided into three groups. The GIR genes were characterized in terms of their phylogenetic relationships, structures, chromosomal distribution and evolutionary dynamics. These GIR genes were revealed to be unequally distributed on 12 chromosomes in the diploid cotton genome, with no GIR gene detected on Ga06. The cis-acting elements in the promoter regions were predicted to be responsive to light, phytohormones, defense activities and stress. The transcriptomic data and qRT-PCR results revealed that most GIR genes were not differentially expressed between the wild-type control and the fuzzless mutant line. Moreover, 14 of 21 family genes were expressed at high levels, indicating these genes may play important roles during fiber development and fuzz formation. Furthermore, Ga01G0231 was predominantly expressed in root samples, suggestive of a role in root hair formation rather than in fuzz initiation and development. The results of this study have enhanced our understanding of the GIR genes and their potential utility for improving cotton fiber through breeding.

Genome-Wide Identification and Characterization of GhCOMT Gene Family during Fiber Development and Verticillium Wilt Resistance in Cotton.

Caffeic acid O-methyltransferases (COMTs) play an essential role in lignin synthesis procession, especially in the plant's phenylalanine metabolic pathway. The content of COMT genes in cotton and the relationship between their expression patterns have not been studied clearly in cotton. In this study, we have identified 190 COMT genes in cotton, which were classified into three groups (I, II and III), and mapped on the cotton chromosomes. In addition, we found that 135 of the 190 COMT genes result from dispersed duplication (DSD) and whole-genome duplication (WGD), indicating that DSD and WGD were the main forces driving COMT gene expansion. The Ka/Ks analysis showed that GhCOMT43 and GhCOMT41 evolved from GaCOMT27 and GrCOMT14 through positive selection. The results of qRT-PCR showed that GhCOMT13, GhCOMT28, GhCOMT39 and GhCOMT55 were related to lignin content during the cotton fiber development. GhCOMT28, GhCOMT39, GhCOMT55, GhCOMT56 and GhCOMT57 responded to Verticillium Wilt (VW) and maybe related to VW resistance through lignin synthesis. Conclusively, this study found that GhCOMTs were highly expressed in the secondary wall thickening stage and VW. These results provide a clue for studying the functions of GhCOMTs in the development of cotton fiber and VW resistance and could lay a foundation for breeding cotton cultivates with higher quantity and high resistance to VW.

Genome-Wide Analysis of AAT Genes and Their Expression Profiling during Fiber Development in Cotton.

Amino acid transporters (AATs) are a kind of membrane proteins that mediate the transport of amino acids across cell membranes in higher plants. The AAT proteins are involved in regulating plant cell growth and various developmental processes. However, the biological function of this gene family in cotton fiber development is not clear. In this study, 190, 190, 101, and 94 full-length AAT genes were identified from Gossypiumhirsutum, G. barbadense, G. arboreum, and G. raimondii. A total of 575 AAT genes from the four cotton species were divided into two subfamilies and 12 clades based on phylogenetic analysis. The AAT genes in the four cotton species were distributed on all the chromosomes. All GhAAT genes contain multiple exons, and each GhAAT protein has multiple conserved motifs. Transcriptional profiling and RT qPCR analysis showed that four GhATT genes tend to express specifically at the fiber initiation stage. Eight genes tend to express specifically at the fiber elongation and maturity stage, and four genes tend to express specifically at the fiber initiation and elongation stages. Our results provide a solid basis for further elucidating the biological function of AAT genes related to cotton fiber development and offer valuable genetic resources for crop improvement in the future.

Identification of the Group III WRKY Subfamily and the Functional Analysis of GhWRKY53 in Gossypium hirsutum L.

WRKY transcription factors had multiple functions in plant secondary metabolism, leaf senescence, fruit ripening, adaptation to biotic and abiotic stress, and plant growth and development. However, knowledge of the group III WRKY subfamily in fiber development in upland cotton (Gossypium hirsutum L.) is largely absent. Previous studies have shown that there were 21 putative group III WRKY members in G. hirsutum L. These putative amino acid sequences from the III WRKY group were phylogenetically clustered into three clades. Multiple alignment, conservative motif analysis, and gene structure analysis showed that the members clustered together in the phylogenetic tree had similar motifs and gene structures. Expression pattern analysis revealed that variation in the expression levels of these genes in different tissues and fiber development stages. To better understand the functions of putative group III WRKY genes in G. hirsutum L., we selected the cotton fiber initiation-related gene GhWRKY53 for cloning and functional identification. The subcellular localization experiment of GhWRKY53 in Nicotiana tabacum leaves showed that it was located in the nucleus. The heterologous expression of GhWRKY53 in Arabidopsis thaliana could significantly increase the density of trichomes. Twelve proteins that interacted with GhWRKY53 were screened from the cotton fiber cDNA library by yeast two-hybrid experiment. This study findings lay a foundation for further research on the role of the GhWRKY53 during cotton fiber development and provide a new insight for further studying putative group III WRKY genes in G. hirsutum L. Our research results also provide vital information for the genetic mechanism of high-quality cotton fiber formation and essential genetic resources for cotton fiber quality improvement.

Weighted Gene Co-Expression Network Analysis Reveals Hub Genes Contributing to Fuzz Development in Gossypium arboreum.

Fuzzless mutants are ideal materials to decipher the regulatory network and mechanism underlying fuzz initiation and formation. In this study, we utilized two Gossypium arboreum accessions differing in fuzz characteristics to explore expression pattern differences and discriminate genes involved in fuzz development using RNA sequencing. Gene ontology (GO) analysis was conducted and found that DEGs were mainly enriched in the regulation of transcription, metabolic processes and oxidation-reduction-related processes. Weighted gene co-expression network analysis discerned the MEmagenta module highly associated with a fuzz/fuzzless trait, which included a total of 50 hub genes differentially expressed between two materials. GaFZ, which negatively regulates trichome and fuzz formation, was found involved in MEmagenta cluster1. In addition, twenty-eight hub genes in MEmagenta cluster1 were significantly up-regulated and expressed in fuzzless mutant DPL972. It is noteworthy that Ga04G1219 and Ga04G1240, which, respectively, encode Fasciclin-like arabinogalactan protein 18(FLA18) and transport protein, showed remarkable differences of expression level and implied that they may be involved in protein glycosylation to regulate fuzz formation and development. This module and hub genes identified in this study will provide new insights on fiber and fuzz formation and be useful for the molecular design breeding of cotton genetic improvement.

Gene Expression Correlation Analysis Reveals MYC-NAC Regulatory Network in Cotton Pigment Gland Development.

Plant NAC (NAM, ATAF1/2, and CUC2) family is involved in various development processes including Programmed Cell Death (PCD) associated development. However, the relationship between NAC family and PCD-associated cotton pigment gland development is largely unknown. In this study, we identified 150, 153 and 299 NAC genes in newly updated genome sequences of G. arboreum, G. raimondii and G. hirsutum, respectively. All NAC genes were divided into 8 groups by the phylogenetic analysis and most of them were conserved during cotton evolution. Using the vital regulator of gland formation GhMYC2-like as bait, expression correlation analysis screened out 6 NAC genes which were low-expressed in glandless cotton and high-expressed in glanded cotton. These 6 NAC genes acted downstream of GhMYC2-like and were induced by MeJA. Silencing CGF1(Cotton Gland Formation1), another MYC-coding gene, caused almost glandless phenotype and down-regulated expression of GhMYC2-like and the 6 NAC genes, indicating a MYC-NAC regulatory network in gland development. In addition, predicted regulatory mechanism showed that the 6 NAC genes were possibly regulated by light, various phytohormones and transcription factors as well as miRNAs. The interaction network and DNA binding sites of the 6 NAC transcription factors were also predicted. These results laid the foundation for further study of gland-related genes and gland development regulatory network.

Molecular cloning and characterization of GhERF105, a gene contributing to the regulation of gland formation in upland cotton (Gossypium hirsutum L.).

BACKGROUND: Gossypium hirsutum L. (cotton) is one of the most economically important crops in the world due to its significant source of fiber, feed, foodstuff, oil and biofuel products. However, the utilization of cottonseed was limited due to the presence of small and darkly pigmented glands that contain large amounts of gossypol, which is toxic to human beings and non-ruminant animals. To date, some progress has been made in the pigment gland formation, but the underlying molecular mechanism of its formation was still unclear. RESULTS: In this study, we identified an AP2/ERF transcription factor named GhERF105 (GH_A12G2166), which was involved in the regulation of gland pigmentation by the comparative transcriptome analysis of the leaf of glanded and glandless plants. It encoded an ERF protein containing a converved AP2 domain which was localized in the nucleus with transcriptional activity, and showed the high expression in glanded cotton accessions that contained much gossypol. Virus-induced gene silencing (VIGS) against GhERF105 caused the dramatic reduction in the number of glands and significantly lowered levels of gossypol in cotton leaves. GhERF105 showed the patterns of spatiotemporal and inducible expression in the glanded plants. CONCLUSIONS: These results suggest that GhERF105 contributes to the pigment gland formation and gossypol biosynthesis in partial organs of glanded plant. It also provides a potential molecular basis to generate 'glandless-seed' and 'glanded-plant' cotton cultivar.

Genome-wide identification of the expansin gene family reveals that expansin genes are involved in fibre cell growth in cotton.

BACKGROUND: Expansins (EXPs), a group of proteins that loosen plant cell walls and cellulosic materials, are involved in regulating cell growth and diverse developmental processes in plants. However, the biological functions of this gene family in cotton are still unknown. RESULTS: In this paper, we identified a total of 93 expansin genes in Gossypium hirsutum. These genes were classified into four subfamilies, including 67 GhEXPAs, 8 GhEXPBs, 6 GhEXLAs, and 12 GhEXLBs, and divided into 15 subgroups. The 93 expansin genes are distributed over 24 chromosomes, excluding Ghir_A02 and Ghir_D06. All GhEXP genes contain multiple exons, and each GhEXP protein has multiple conserved motifs. Transcript profiling and qPCR analysis revealed that the expansin genes have distinct expression patterns among different stages of cotton fibre development. Among them, 3 genes (GhEXPA4o, GhEXPA1A, and GhEXPA8h) were highly expressed in the initiation stage, 9 genes (GhEXPA4a, GhEXPA13a, GhEXPA4f, GhEXPA4q, GhEXPA8f, GhEXPA2, GhEXPA8g, GhEXPA8a, and GhEXPA4n) had high expression during the fast elongation stage, and GhEXLA1c and GhEXLA1f were preferentially expressed in the transition stage of fibre development. CONCLUSIONS: Our results provide a solid basis for further elucidation of the biological functions of expansin genes in relation to cotton fibre development and valuable genetic resources for future crop improvement.