Genome-wide identification of GhRLCK-VII subfamily genes in Gossypium hirsutum and investigation of their functions in resistance to Verticillium wilt.

BACKGROUND: The receptor-like cytoplasmic kinases subfamily VII (RLCK-VII) is critical in regulating plant growth, development, and pattern-triggered immunity. However, a comprehensive exploration of these genes in the allotetraploid Gossypium hirsutum is still lacking. This study aimed to identify RLCK-VII genes in G. hirsutum and investigate their evolutionary history, structural features, expression patterns, and role in plant defense. RESULTS: Seventy-two RLCK-VII genes in the G. hirsutum genome were unveiled and classified into nine groups following their phylogenetic analysis with Arabidopsis thaliana. Group VII-1 was the largest, accounting for 28%, while Groups VII-2 and VII-3 had only one member each. The analysis using MCScanX revealed that these 72 genes formed 166 collinear gene pairs and were resided on 26 chromosomes of G. hirsutum, suggesting that they were derived from whole genome segmental duplication events. Their calculated Ka/Ks values were below one, implying the occurrence of purification selection during the evolution and inhibition of gene function differentiation/loss. All members of the RLCK-VII subfamily possessed two conserved domains, PKinase-Tyr and PKinase, and several conserved PBS1 kinase subdomains, individually included in one of the ten motifs identified using MEME. The RNA-Seq results showed that RLCK-VII genes exhibited different spatiotemporal expression, indicating their involvement in cotton growth, development, and defense responses to Verticillium dahliae. The transcription patterns of RLCK-VII genes found by RNA-Seq were further validated using qRT-PCR assays after inoculating "20B12" (cotton cultivar) with "V991" (V. dahliae). The virus-induced gene silencing (VIGS) assays uncovered that two RLCK-VII genes (Gohir.A13G227248 and Gohir.A10G219900) were essential to G. hirsutum resistance to Verticillium wilt. CONCLUSIONS: These observations offer valuable insight into the attributes and roles of RLCK-VII genes in G. hirsutum, potentially enable the breeding of new cotton cultivars with enhanced resistance to Verticillium wilt.

Identification of small effect quantitative trait loci of plant architectural, flowering, and early maturity traits in reciprocal interspecific introgression population in cotton.

Plant architecture, flowering time and maturity traits are important determinants of yield and fiber quality of cotton. Genetic dissection of loci determining these yield and quality components is complicated by numerous loci with alleles conferring small differences. Therefore, mapping populations segregating for smaller numbers and sizes of introgressed segments is expected to facilitate dissection of these complex quantitative traits. At an advanced stage in the development of reciprocal advanced backcross populations from crosses between elite Gossypium hirsutum cultivar 'Acala Maxxa' (GH) and G. barbadense 'Pima S6' (GB), we undertook mapping of plant architectural traits, flowering time and maturity. A total of 284 BC4F1 and BC4F2 progeny rows, 120 in GH and 164 in GB background, were evaluated for phenotype, with only 4 and 3 (of 7) traits showing significant differences among progenies. Genotyping by sequencing yielded 3,186 and 3,026 SNPs, respectively, that revealed a total of 27 QTLs in GH background and 22 in GB, for plant height, days to flowering, residual flowering at maturity and maturity. More than of 90% QTLs identified in both backgrounds had small effects (%PV < 10), supporting the merit of this population structure to reduce background noise and small effect QTLs. Germplasm developed in this study may serve as potential pre-breeding material to develop improved cotton cultivars.

Comparative Transcriptome Analysis Reveals Genes Associated with the Gossypol Synthesis and Gland Morphogenesis in Gossypium hirsutum.

Gossypium hirsutum is an important source of natural textile fibers. Gossypol, which is a sesquiterpenoid compound mainly existing in the cotton pigment glands, can facilitate resistance to the stress from diseases and pests. The level of gossypol in the cotton is positively correlated to the quantity of pigment glands. However, the underlying regulatory mechanisms of gossypol synthesis and gland morphogenesis are still poorly understood, especially from a transcriptional perspective. The transcripts of young leaves and ovules at 30 DPA of the glanded plants and glandless plants were studied by RNA-Seq and 865 million clean reads were obtained. A total of 34,426 differentially expressed genes (DEGs) were identified through comparative transcriptome analysis. Genes related to gossypol synthesis or gland morphogenesis displayed significant differential expression between the two cultivars. Functional annotation revealed that the candidate genes related to catalytic activity, the biosynthesis of secondary metabolites, and biomolecular decomposition processes. Our work herein unveiled several potential candidate genes related to gossypol synthesis or gland morphogenesis and may provide useful clues for a breeding program of cotton cultivars with low cottonseed gossypol contents.

Genome-wide identification and characterization of terpene synthase genes in Gossypium hirsutum.

Terpenoids are widely distributed in plants and play important roles in the regulation of plant growth and development and in the interactions between plants and both the environment and other organisms. However, terpene synthase (TPS) genes have not been systematically investigated in the tetraploid Gossypium hirsutum. In this study, whole genome identification and characterization of the TPS family from G. hirsutum were carried out. Eighty-five TPS genes, including 47 previously unidentified genes, were identified in the G. hirsutum genome and classified into 5 subfamilies according to protein sequence similarities, as follows: 43 GhTPS-a, 29 GhTPS-b, 4 GhTPS-c, 7 GhTPS-e/f, and 2 GhTPS-g members. These 85 TPS genes were mapped onto 19 chromosomes of the G. hirsutum genome. Segmental duplications and tandem duplications contributed greatly to the expansion of TPS genes in G. hirsutum and were followed by intense purifying selection during evolution. Indentification of cis-acting regulatory elements suggest that the expression of TPS genes is regulated by a variety of hormones. RNA sequencing (RNA-seq) expression profile analysis revealed that the TPS genes had distinct spatiotemporal expression patterns, and several genes were highly and preferentially expressed in the leaves of cotton with gossypol glands (glanded cotton) versus a glandless strain. Virus-induced gene silencing (VIGS) of three TPS genes yielded plants characterized by fewer, smaller, and lighter gossypol glands, which indicated that these three genes were responsible for gland activity. Taken together, our results provide a solid basis for further elucidation of the biological functions of TPS genes in relation to gland activity and gossypol biosynthesis to develop cotton cultivars with low cottonseed gossypol contents.

Genetic Analysis of Gossypium Fiber Quality Traits in Reciprocal Advanced Backcross Populations.

In mapping populations segregating for many loci, the large amount of variation among genotypes often masks small-effect quantitative trait loci (QTL). This problem can be reduced by development of populations with fewer chromosome segments segregating. Here, we report early QTL detection in reciprocal advanced backcross populations from crosses between elite Gossypium hirsutum L. 'Acala Maxxa' (GH) and G. barbadense L. 'Pima S6' (GB). A total of 297 BCF and BCF progeny rows-127 segregating for GB chromosome segments in GH background and 170 segregating for GH chromosome segments in GB background-were evaluated in three environments. Totals of 3186 and 3026 polymorphic single-nucleotide polymorphisms (SNPs) in GH and GB backgrounds, respectively, were identified and used for trait mapping. Small-effect QTL (<10% variance explained) made up 87 and 100% of QTL in GH and GB backgrounds, respectively. In both species, favorable alleles were found with effects being masked or neutralized by unfavorable alleles, with greater scope for improvement of GH than GB by introgressive breeding. A total of three stable QTL-two in GH background for fiber elongation (ELO) and micronaire (MIC) and one in GB background for upper-half mean length (UHM)-were identified in two out of three environments. Curiously, only four QTL-three for UHM and one for ELO-showed the expected opposite effects in reciprocal backgrounds, perhaps reflecting the combined consequences of epistasis, small phenotypic effects, and low coverage of some genomic regions. Along with new information for marker-assisted breeding, this study adds to knowledge that can be used to unravel complex genetic networks governing fiber quality traits.

A global survey of alternative splicing in allopolyploid cotton: landscape, complexity and regulation.

Alternative splicing (AS) is a crucial regulatory mechanism in eukaryotes, which acts by greatly increasing transcriptome diversity. The extent and complexity of AS has been revealed in model plants using high-throughput next-generation sequencing. However, this technique is less effective in accurately identifying transcript isoforms in polyploid species because of the high sequence similarity between coexisting subgenomes. Here we characterize AS in the polyploid species cotton. Using Pacific Biosciences single-molecule long-read isoform sequencing (Iso-Seq), we developed an integrated pipeline for Iso-Seq transcriptome data analysis (https://github.com/Nextomics/pipeline-for-isoseq). We identified 176 849 full-length transcript isoforms from 44 968 gene models and updated gene annotation. These data led us to identify 15 102 fibre-specific AS events and estimate that c. 51.4% of homoeologous genes produce divergent isoforms in each subgenome. We reveal that AS allows differential regulation of the same gene by miRNAs at the isoform level. We also show that nucleosome occupancy and DNA methylation play a role in defining exons at the chromatin level. This study provides new insights into the complexity and regulation of AS, and will enhance our understanding of AS in polyploid species. Our methodology for Iso-Seq data analysis will be a useful reference for the study of AS in other species.

Genome-wide identification and analysis of the aldehyde dehydrogenase (ALDH) gene superfamily of Gossypium raimondii.

BACKGROUND: Aldehyde dehydrogenases (ALDHs) are members of the NAD(P)(+)-dependent protein superfamily which catalyzes aliphatic and aromatic aldehyde oxidation to non-toxic carboxylic acids. ALDH genes may offer promise for improving plant adaptation to environmental stress. Recently, elucidated genome sequences of Gossypium raimondii provide a foundation for systematic identification and analysis of ALDH genes. To date, this has been accomplished for many plant species except G. raimondii. RESULTS: In this study, thirty unique ALDH sequences that code for 10 ALDH families were identified in the G. raimondii genome. Phylogenetic analysis revealed that ALDHs were split into six clades in G. raimondii, and ALDH proteins from the same families were clustered together. Phylogenetic relationships of ALDHs from 11 plant species suggest that ALDHs in G. raimondii shared the highest protein homology with ALDHs from poplar. Members within ALDH families possessed homologous exon-intron structures. Chromosomal distribution of ALDH did not occur evenly in the G. raimondii genome and many ALDH genes were involved in the syntenic region as documented by identification of physical locations among single chromosomes. In addition, syntenic analysis revealed that homologues of many G. raimondii ALDHs appeared in corresponding Arabidopsis and poplar syntenic blocks, indicating that these genes arose prior to G. raimondii, Arabidopsis and poplar speciation. Finally, based on gene expression analysis of microarray and RNA-seq, we can speculate that some G. raimondii ALDH genes might respond to drought or waterlogging stresses. CONCLUSION: Genome-wide identification and analysis of the evolution and expression of ALDH genes in G. raimondii laid a foundation for studying this gene superfamily and offers new insights into the evolution history and speculated roles in Gossypium. These data can be used to inform functional genomic studies and molecular breeding in cotton.

[Advance in quantitative traits from phenotypic variation to gene discovery].

With the rapid development of molecular biotechnology, many improved approaches have been provided for studying on the genetic bases of quantitative traits. Many studies paid attention to dissect the genetic bases using different emendatory approaches from phenotypic variation to gene discovery, mainly containing QTL mapping, chemical mutagenesis etc. Especially, QTL analysis was executed for a lot of important quantitative traits. This paper summarizes the development of quantitative genetics, the recent progress on QTL mapping populations and methods, the status of QTL locating and QTG identification, positional cloning, and QTL application in breeding.

Localization of QTLs for yield and fiber quality traits of tetraploid cotton cultivar.

Using SSR and RAPD as molecular markers, and the 69 F2 families from a cross between Handan208 (Gossypium hirsutum) and Pima90 (Gossypium barbadense) as a mapping population, a linkage map comprising 126 markers was constructed. With an average spacing 13.7 cM between markers, the linkage map spanned 1,717.0 cM, which covers approximately 34.34% of the total recombinational length of cotton genome. Based on the linkage map and the F2:3 phenotypic data, overall genome QTL screening was conducted by composite interval mapping method, then 29 QTLs related to cotton yield traits and fiber quality traits were found. There are 16 loci related to yield traits and 13 loci related to fiber quality. Particularly, two QTLs relevant to fiber strength were detected, which could explain phenotypic variance of 34.15% and 13.86%, respectively.