Construction and characterization of the first bacterial artificial chromosome library for the cotton species Gossypium barbadense L

Transcriptome, Ectopic Expression and Genetic Population Analysis Identify Candidate Genes for Fiber Quality Improvement in Cotton

Comparative transcriptome analysis of fiber tissues between Gossypium barbadense and Gossypium hirsutum could reveal the molecular mechanisms underlying high-quality fiber formation and identify candidate genes for fiber quality improvement. In this study, 759 genes were found to be strongly upregulated at the elongation stage in G. barbadense, which showed four distinct expression patterns (I-IV). Among them, the 346 genes of group IV stood out in terms of the potential to promote fiber elongation, in which we finally identified 42 elongation-related candidate genes by comparative transcriptome analysis between G. barbadense and G. hirsutum. Subsequently, we overexpressed GbAAR3 and GbTWS1, two of the 42 candidate genes, in Arabidopsis plants and validated their roles in promoting cell elongation. At the secondary cell wall (SCW) biosynthesis stage, 2275 genes were upregulated and exhibited five different expression profiles (I-V) in G. barbadense. We highlighted the critical roles of the 647 genes of group IV in SCW biosynthesis and further picked out 48 SCW biosynthesis-related candidate genes by comparative transcriptome analysis. SNP molecular markers were then successfully developed to distinguish the SCW biosynthesis-related candidate genes from their G. hirsutum orthologs, and the genotyping and phenotyping of a BC3F5 population proved their potential in improving fiber strength and micronaire. Our results contribute to the better understanding of the fiber quality differences between G. barbadense and G. hirsutum and provide novel alternative genes for fiber quality improvement.

Evolution, expression and functional analysis of cultivated allotetraploid cotton DIR genes

BACKGROUND

Dirigent (DIR) proteins mediate regioselectivity and stereoselectivity during lignan biosynthesis and are also involved in lignin, gossypol and pterocarpan biosynthesis. This gene family plays a vital role in enhancing stress resistance and in secondary cell-wall development, but systematical understanding is lacking in cotton.

RESULTS

In this study, 107 GbDIRs and 107 GhDIRs were identified in Gossypium barbadense and Gossypium hirsutum, respectively. Most of these genes have a classical gene structure without intron and encode proteins containing a signal peptide. Phylogenetic analysis showed that cotton DIR genes were classified into four distinct subfamilies (a, b/d, e, and f). Of these groups, DIR-a and DIR-e were evolutionarily conserved, and segmental and tandem duplications contributed equally to their formation. In contrast, DIR-b/d mainly expanded by recent tandem duplications, accompanying with a number of gene clusters. With the rapid evolution, DIR-b/d-III was a Gossypium-specific clade involved in atropselective synthesis of gossypol. RNA-seq data highlighted GhDIRs in response to Verticillium dahliae infection and suggested that DIR gene family could confer Verticillium wilt resistance. We also identified candidate DIR genes related to fiber development in G. barbadense and G. hirsutum and revealed their differential expression. To further determine the involvement of DIR genes in fiber development, we overexpressed a fiber length-related gene GbDIR78 in Arabidopsis and validated its function in trichomes and hypocotyls.

CONCLUSIONS

These findings contribute novel insights towards the evolution of DIR gene family and provide valuable information for further understanding the roles of DIR genes in cotton fiber development as well as in stress responses.

Construction and characterization of a bacterial artificial chromosome library for Gossypium mustelinum

A bacterial artificial chromosome (BAC) library for G. mustelinum Miers ex G. Watt (AD₄) was constructed. Intact nuclei from G. mustelinum (AD₄) were used to isolate high molecular weight DNA, which was partially cleaved with Hind III and cloned into pSMART BAC (Hind III) vectors. The BAC library consisted of 208,182 clones arrayed in 542 384-microtiter plates, with an average insert size of 121.72 kb ranging from 100 to 150 kb. About 2% of the clones did not contain inserts. Based on an estimated genome size of 2372 Mb for G. mustelinum, the BAC library was estimated to have a total coverage of 10.50 × genome equivalents. The high capacity library of G. mustelinum will serve as a giant gene resource for map-based cloning of quantitative trait loci or genes associated with important agronomic traits or resistance to Verticillium wilt, physical mapping and comparative genome analysis.

Construction of BAC contig maps of homoeologous chromosomes A12 and D12 of Gossypium hirsutum L. acc. TM-1

Background

The Gossypium hirsutum homoeologous chromosome 12 encodes important genes that contribute to fiber fuzz, lethality, gland development and male sterility. In this study a physical map of the cotton TM-1 chromosome 12 was constructed. A number of large-insert cotton genome libraries are available, and genome-wide physical mapping using large insert segments combined with bacterial cloning is a thriving area of genome research. However, sequencing of the cotton genome is difficult due to sequence repeats and homoeologous regions. In order to effectively distinguish the homologous segments, a new method for adjusting the parameters of the FPC software was applied for contig map construction.

Results

All available markers on chromosomes A12 and D12 were used to screen the TM-1 BAC library by PCR. A total of 775 clones (387 for A12, 388 for D12) were obtained using Hind III fingerprinting and used for construction of the contig map. Seven pairs of SSR markers located on A12 and D12 were chosen for contig analysis. Following optimization of the tolerance (10) and cutoff (1e-12) parameters, combining all clones from A12 and D12 produced two separate contigs.

Conclusions

The BAC contig map of chromosomes A12 and D12 was constructed and FPC software parameters were optimized for analysis. The resulting approach is a powerful platform for genome-wide and evolutionary research on cotton.

Transcriptome profiling of Gossypium barbadense inoculated with Verticillium dahliae provides a resource for cotton improvement

BACKGROUND

Verticillium wilt, caused by the fungal pathogen Verticillium dahliae, is the most severe disease in cotton (Gossypium spp.), causing great lint losses worldwide. Disease management could be achieved in the field if genetically improved, resistant plants were used. However, the interaction between V. dahliae and cotton is a complicated process, and its molecular mechanism remains obscure. To understand better the defense response to this pathogen as a means for obtaining more tolerant cultivars, we monitored the transcriptome profiles of roots from resistant plants of G. barbadense cv. Pima90-53 that were challenged with V. dahliae.

RESULTS

In all, 46,192 high-quality expressed sequence tags (ESTs) were generated from a full-length cDNA library of G. barbadense. They were clustered and assembled into 23126 unigenes that comprised 2661 contigs and 20465 singletons. Those unigenes were assigned Gene Ontology terms and mapped to 289 KEGG pathways. A total of 3027 unigenes were found to be homologous to known defense-related genes in other plants. They were assigned to the functional classification of plant-pathogen interactions, including disease defenses and signal transduction. The branch of "SA→NPR1→TGA→PR-1→Disease resistance" was first discovered in the interaction of cotton-V. dahliae, indicating that this wilt process includes both biotrophic and necrotrophic stages. In all, 4936 genes coding for putative transcription factors (TF) were identified in our library. The most abundant TF family was the NAC group (527), followed by G2-like (440), MYB (372), BHLH (331), bZIP (271) ERF, C3H, and WRKY. We also analyzed the expression of genes involved in pathogen-associated molecular pattern (PAMP) recognition, the activation of effector-triggered immunity, TFs, and hormone biosynthesis, as well as genes that are pathogenesis-related, or have roles in signaling/regulatory functions and cell wall modification. Their differential expression patterns were compared among mock-/inoculated- and resistant/susceptible cotton. Our results suggest that the cotton defense response has significant transcriptional complexity and that large accumulations of defense-related transcripts may contribute to V. dahliae resistance in cotton. Therefore, these data provide a resource for cotton improvement through molecular breeding approaches.

CONCLUSIONS

This study generated a substantial amount of cotton transcript sequences that are related to defense responses against V. dahliae. These genomics resources and knowledge of important related genes contribute to our understanding of host-pathogen interactions and the defense mechanisms utilized by G. barbadense, a non-model plant system. These tools can be applied in establishing a modern breeding program that uses marker-assisted selections and oligonucleotide arrays to identify candidate genes that can be linked to valuable agronomic traits in cotton, including disease resistance.

Construction of a plant-transformation-competent BIBAC library and genome sequence analysis of polyploid Upland cotton (Gossypium hirsutum L.)

Background

Cotton, one of the world’s leading crops, is important to the world’s textile and energy industries, and is a model species for studies of plant polyploidization, cellulose biosynthesis and cell wall biogenesis. Here, we report the construction of a plant-transformation-competent binary bacterial artificial chromosome (BIBAC) library and comparative genome sequence analysis of polyploid Upland cotton (Gossypium hirsutum L.) with one of its diploid putative progenitor species, G. raimondii Ulbr.

Results

We constructed the cotton BIBAC library in a vector competent for high-molecular-weight DNA transformation in different plant species through either Agrobacterium or particle bombardment. The library contains 76,800 clones with an average insert size of 135 kb, providing an approximate 99% probability of obtaining at least one positive clone from the library using a single-copy probe. The quality and utility of the library were verified by identifying BIBACs containing genes important for fiber development, fiber cellulose biosynthesis, seed fatty acid metabolism, cotton-nematode interaction, and bacterial blight resistance. In order to gain an insight into the Upland cotton genome and its relationship with G. raimondii, we sequenced nearly 10,000 BIBAC ends (BESs) randomly selected from the library, generating approximately one BES for every 250 kb along the Upland cotton genome. The retroelement Gypsy/DIRS1 family predominates in the Upland cotton genome, accounting for over 77% of all transposable elements. From the BESs, we identified 1,269 simple sequence repeats (SSRs), of which 1,006 were new, thus providing additional markers for cotton genome research. Surprisingly, comparative sequence analysis showed that Upland cotton is much more diverged from G. raimondii at the genomic sequence level than expected. There seems to be no significant difference between the relationships of the Upland cotton D- and A-subgenomes with the G. raimondii genome, even though G. raimondii contains a D genome (D₅).

Conclusions

The library represents the first BIBAC library in cotton and related species, thus providing tools useful for integrative physical mapping, large-scale genome sequencing and large-scale functional analysis of the Upland cotton genome. Comparative sequence analysis provides insights into the Upland cotton genome, and a possible mechanism underlying the divergence and evolution of polyploid Upland cotton from its diploid putative progenitor species, G. raimondii.

A high density consensus genetic map of tetraploid cotton that integrates multiple component maps through molecular marker redundancy check

A consensus genetic map of tetraploid cotton was constructed using six high-density maps and after the integration of a sequence-based marker redundancy check. Public cotton SSR libraries (17,343 markers) were curated for sequence redundancy using 90% as a similarity cutoff. As a result, 20% of the markers (3,410) could be considered as redundant with some other markers. The marker redundancy information had been a crucial part of the map integration process, in which the six most informative interspecific Gossypium hirsutum×G. barbadense genetic maps were used for assembling a high density consensus (HDC) map for tetraploid cotton. With redundant markers being removed, the HDC map could be constructed thanks to the sufficient number of collinear non-redundant markers in common between the component maps. The HDC map consists of 8,254 loci, originating from 6,669 markers, and spans 4,070 cM, with an average of 2 loci per cM. The HDC map presents a high rate of locus duplications, as 1,292 markers among the 6,669 were mapped in more than one locus. Two thirds of the duplications are bridging homoeologous A(T) and D(T) chromosomes constitutive of allopolyploid cotton genome, with an average of 64 duplications per A(T)/D(T) chromosome pair. Sequences of 4,744 mapped markers were used for a mutual blast alignment (BBMH) with the 13 major scaffolds of the recently released Gossypium raimondii genome indicating high level of homology between the diploid D genome and the tetraploid cotton genetic map, with only a few minor possible structural rearrangements. Overall, the HDC map will serve as a valuable resource for trait QTL comparative mapping, map-based cloning of important genes, and better understanding of the genome structure and evolution of tetraploid cotton.

Preparations of meiotic pachytene chromosomes and extended DNA fibers from cotton suitable for fluorescence in situ hybridization

Fluorescence in situ hybridization (FISH) has become one of the most important techniques applied in plant molecular cytogenetics. However, the application of this technique in cotton has lagged behind because of difficulties in chromosome preparation. The focus of this article was FISH performed not only on cotton pachytene chromosomes, but also on cotton extended DNA fibers. The cotton pollen mother cells (PMCs) instead of buds or anthers were directly digested in enzyme to completely breakdown the cell wall. Before the routine acetic acid treatment, PMCs were incubated in acetic acid and enzyme mixture to remove the cytoplasm and clear the background. The method of ice-cold Carnoy's solution spreading chromosome was adopted instead of nitrogen removed method to avoid chromosomes losing and fully stretch chromosome. With the above-improved steps, the high-quality well-differentiated pachytene chromosomes with clear background were obtained. FISH results demonstrated that a mature protocol of cotton pachytene chromosomes preparation was presented. Intact and no debris cotton nuclei were obtained by chopping from etiolation cotyledons instead of the conventional liquid nitrogen grinding method. After incubating the nuclei with nucleus lysis buffer on slide, the parallel and clear background DNA fibers were acquired along the slide. This method overcomes the twist, accumulation and fracture of DNA fibers compared with other methods. The entire process of DNA fibers preparation requires only 30 min, in contrast, it takes 3 h with routine nitrogen grinding method. The poisonous mercaptoethanol in nucleus lysis buffer is replaced by nonpoisonous dithiothreitol. PVP40 in nucleus isolation buffer is used to prevent oxidation. The probability of success in isolating nuclei for DNA fiber preparation is almost 100% tested with this method in cotton. So a rapid, safe, and efficient method for the preparation of cotton extended DNA fibers suitable for FISH was established.

Cloning and characterization of a Verticillium wilt resistance gene from Gossypium barbadense and functional analysis in Arabidopsis thaliana

Verticillium wilt causes enormous loss to yield or quality in many crops. In an effort to help controlling this disease through genetic engineering, we first cloned and characterized a Verticillium wilt resistance gene (GbVe) from cotton (Gossypium barbadense) and analyzed its function in Arabidopsis thaliana. Its nucleotide sequence is 3,819 bp long, with an open reading frame of 3,387 bp, and encoding an 1,128-aa protein precursor. Sequence analysis shows that GbVe produces a leucine-rich repeat receptor-like protein. It shares identities of 55.9% and 57.4% with tomato Ve1 and Ve2, respectively. Quantitative real-time PCR indicated that the Ve gene expression pattern was different between the resistant and susceptible cultivars. In the resistant Pima90-53, GbVe was quickly induced and reached to a peak at 2 h after inoculation, two-fold higher than that of control. We localized the GbVe-GFP fusion protein to the cytomembrane in onion epidermal cells. By inserting GbVe into Arabidopsis via Agrobacterium-mediated transformation, T(3) transgenic lines were obtained. Compared with the wild-type control, GbVe-overexpressing plants had greater levels of resistance to V. dahliae. This suggests that GbVe is a useful gene for improving the plant resistance against fungal diseases.

Construction and characterization of a bacterial artificial chromosome library for the A-genome of cotton (G. arboreum L.)

A bacterial artificial chromosome (BAC) library for the A-genome of cotton has been constructed from the leaves of G. arboreum L cv. Jianglinzhongmian. It is used as elite A-genome germplasm resources in the present cotton breeding program and has been used to build a genetic reference map of cotton. The BAC library consists of 123,648 clones stored in 322 384-well plates. Statistical analysis of a set of 103 randomly selected BAC clones indicated that each clone has an average insert length of 100.2 kb per plasmid, with a range of 30 to 190 kb. Theoretically, this represents 7.2 haploid genome equivalents based on an A-genome size of 1697 Mb. The BAC library has been arranged in column pools and superpools allowing screening with various PCR-based markers. In the future, the A-genome cotton BAC library will serve as both a giant gene resource and a valuable tool for map-based gene isolation, physical mapping and comparative genome analysis.

Construction of a bacterial artificial chromosome library of TM-1, a standard line for genetics and genomics in Upland cotton

A bacterial artificial chromosome (BAC) library was constructed for Gossypium hirsutum acc. TM-1, a genetic and genomic standard line for Upland cotton. The library consists of 147 456 clones with an average insert size of 122.8 kb ranging from 97 to 240 kb. About 96.0% of the clones have inserts over 100 kb. Therefore, this library represents theoretically 7.4 haploid genome equivalents based on an AD genome size of 2 425 Mb. Clones were stored in 384 384(-) well plates and arrayed into multiplex pools for rapid and reliable library screening. BAC screening was carried out by four-round polymerase chain reactions using 23 simple sequence repeats (SSR) markers, three sequence-related amplified polymorphism markers and one pair of primers for a gene associated with fiber development to test the quality of the library. Correspondingly, in total 92 positive BAC clones were identified with an average four positive clones per SSR marker, ranging from one to eight hits. Additionally, since these SSR markers have been localized to chromosome 12 (A12) and 26 (D12) according to the genetic map, these BAC clones are expected to serve as seeds for the physical mapping of these two homologous chromosomes, sequentially map-based cloning of quantitative trait loci or genes associated with important agronomic traits.