A preliminary analysis of genome structure and composition in Gossypium hirsutum

Transcriptomics for Drought Stress Mediated by Biological Processes in-relation to Key Regulated Pathways in Gossypium darwinii

BACKGROUND

Wild cotton Gossypium darwinii, an allotetraploid harbours important traits useful for tolerating abiotic stress, i.e., drought, salt and good genetic stability, hence these characteristics can be transferred to cultivated cotton for genetic improvement.

MATERIALS AND METHODS

In this study, we analyzed the RNA-seq transcriptomes from leaves of G. darwinii seedlings with and without drought stress. A total of 86.7 million valid reads with an average length of 95.79 bp were generated from the two samples and 58,960 transcripts with a length of more than 500 bp were assembled. We searched the known proteins on the strength of sequence similarity; these transcripts were annotated with COG, KEGG and GO functional categories. According to gene expression abundance RPKM value, we carried out RT-qPCR analysis to determine the expression pattern of the obtained transcription factors.

RESULTS

A total of 58,960 genes was differentially expressed (DEG), with 32,693 and 25,919 genes found to be upregulated and downregulated, respectively. Through gene ontology and KEGG pathways, the upregulated genes were found to associate with all the GO terms, molecular functions (MF), biological process (BP) and cellular components (CC), which are highly linked to enhancing drought stress tolerance.

CONCLUSION

The study provides an in-depth knowledge of regulation of pathways and genes involved in photosynthesis during drought stress in G. darwinii. These pathways and genes were found to be significantly downregulated and this information could be further utilized by cotton breeders in developing a more drought tolerant cotton germplasm.

Map-Based Functional Analysis of the GhNLP Genes Reveals Their Roles in Enhancing Tolerance to N-Deficiency in Cotton

Nitrogen is a key macronutrient needed by plants to boost their production, but the development of cotton genotypes through conventional approaches has hit a bottleneck due to the narrow genetic base of the elite cotton cultivars, due to intensive selection and inbreeding. Based on our previous research, in which the BC₂F₂ generations developed from two upland cotton genotypes, an abiotic stress-tolerant genotype, G. tomentosum (donor parent) and a highly-susceptible, and a highly-susceptible, but very productive, G. hirsutum (recurrent parent), were profiled under drought stress conditions. The phenotypic and the genotypic data generated through genotyping by sequencing (GBS) were integrated to map drought-tolerant quantitative trait loci (QTLs). Within the stable QTLs region for the various drought tolerance traits, a nodule-inception-like protein (NLP) gene was identified. We performed a phylogenetic analysis of the NLP proteins, mapped their chromosomal positions, intron-exon structures and conducted ds/dn analysis, which showed that most NLP genes underwent negative or purifying selection. Moreover, the functions of one of the highly upregulated genes, Gh_A05G3286 (Gh NLP5), were evaluated using the virus gene silencing (VIGS) mechanism. A total of 226 proteins encoded by the NLP genes were identified, with 105, 61, and 60 in Gossypium hirsutum, G. raimondii, and G. arboreum, respectively. Comprehensive Insilico analysis revealed that the proteins encoded by the NLP genes had varying molecular weights, protein lengths, isoelectric points (pI), and grand hydropathy values (GRAVY). The GRAVY values ranged from a negative one to zero, showing that proteins were hydrophilic. Moreover, various cis-regulatory elements that are the binding sites for stress-associated transcription factors were found in the promoters of various NLP genes. In addition, many miRNAs were predicted to target NLP genes, notably miR167a, miR167b, miR160, and miR167 that were previously shown to target five NAC genes, including NAC1 and CUC1, under N-limited conditions. The real-time quantitative polymerase chain reaction (RT-qPCR) analysis, revealed that five genes, Gh_D02G2018, Gh_A12G0439, Gh_A03G0493, Gh_A03G1178, and Gh_A05G3286 were significantly upregulated and perhaps could be the key NLP genes regulating plant response under N-limited conditions. Furthermore, the knockdown of the Gh_A05G3286 (GhNLP5) gene by virus-induced silencing (VIGS) significantly reduced the ability of these plants to the knockdown of the Gh_A05G3286 (GhNLP5) gene by virus-induced gene silencing (VIGS) significantly reduced the ability of the VIGS-plants to tolerate N-limited conditions compared to the wild types (WT). The VIGS-plants registered lower chlorophyll content, fresh shoot biomass, and fresh root biomass, addition to higher levels of malondialdehyde (MDA) and significantly reduced levels of proline, and superoxide dismutase (SOD) compared to the WT under N-limited conditions. Subsequently, the expression levels of the Nitrogen-stress responsive genes, GhTap46, GhRPL18A, and GhKLU were shown to be significantly downregulated in VIGS-plants compared to their WT under N-limited conditions. The downregulation of the nitrogen-stress responsive genes provided evidence that the silenced gene had an integral role in enhancing cotton plant tolerance to N-limited conditions.

Gossypium barbadense and Gossypium hirsutum genomes provide insights into the origin and evolution of allotetraploid cotton

Allotetraploid cotton is an economically important natural-fiber-producing crop worldwide. After polyploidization, Gossypium hirsutum L. evolved to produce a higher fiber yield and to better survive harsh environments than Gossypium barbadense, which produces superior-quality fibers. The global genetic and molecular bases for these interspecies divergences were unknown. Here we report high-quality de novo-assembled genomes for these two cultivated allotetraploid species with pronounced improvement in repetitive-DNA-enriched centromeric regions. Whole-genome comparative analyses revealed that species-specific alterations in gene expression, structural variations and expanded gene families were responsible for speciation and the evolutionary history of these species. These findings help to elucidate the evolution of cotton genomes and their domestication history. The information generated not only should enable breeders to improve fiber quality and resilience to ever-changing environmental conditions but also can be translated to other crops for better understanding of their domestication history and use in improvement.

Gossypium barbadense and Gossypium hirsutum genomes provide insights into the origin and evolution of allotetraploid cotton

Allotetraploid cotton is an economically important natural-fiber-producing crop worldwide. After polyploidization, Gossypium hirsutum L. evolved to produce a higher fiber yield and to better survive harsh environments than Gossypium barbadense, which produces superior-quality fibers. The global genetic and molecular bases for these interspecies divergences were unknown. Here we report high-quality de novo-assembled genomes for these two cultivated allotetraploid species with pronounced improvement in repetitive-DNA-enriched centromeric regions. Whole-genome comparative analyses revealed that species-specific alterations in gene expression, structural variations and expanded gene families were responsible for speciation and the evolutionary history of these species. These findings help to elucidate the evolution of cotton genomes and their domestication history. The information generated not only should enable breeders to improve fiber quality and resilience to ever-changing environmental conditions but also can be translated to other crops for better understanding of their domestication history and use in improvement.

Recent insights into cotton functional genomics: progress and future perspectives

Functional genomics has transformed from futuristic concept to well-established scientific discipline during the last decade. Cotton functional genomics promise to enhance the understanding of fundamental plant biology to systematically exploit genetic resources for the improvement of cotton fibre quality and yield, as well as utilization of genetic information for germplasm improvement. However, determining the cotton gene functions is a much more challenging task, which has not progressed at a rapid pace. This article presents a comprehensive overview of the recent tools and resources available with the major advances in cotton functional genomics to develop elite cotton genotypes. This effort ultimately helps to filter a subset of genes that can be used to assemble a final list of candidate genes that could be employed in future novel cotton breeding programme. We argue that next stage of cotton functional genomics requires the draft genomes refinement, re-sequencing broad diversity panels with the development of high-throughput functional genomics tools and integrating multidisciplinary approaches in upcoming cotton improvement programmes.

Fine mapping and candidate gene analysis of qFL-chr1, a fiber length QTL in cotton

A fiber length QTL, qFL-chr1, was fine mapped to a 0.9 cM interval of cotton chromosome 1. Two positional candidate genes showed positive correlation between gene expression level and fiber length. Prior analysis of a backcross-self mapping population derived from a cross between Gossypium hirsutum L. and G. barbadense L. revealed a QTL on chromosome 1 associated with increased fiber length (qFL-chr1), which was confirmed in three independent populations of near-isogenic introgression lines (NIILs). Here, a single NIIL, R01-40-08, was used to develop a large population segregating for the target region. Twenty-two PCR-based polymorphic markers used to genotype 1672 BC₄F₂ plants identified 432 recombinants containing breakpoints in the target region. Substitution mapping using 141 informative recombinants narrowed the position of qFL-chr1 to a 1.0-cM interval between SSR markers MUSS084 and CIR018. To exclude possible effects of non-target introgressions on fiber length, different heterozygous BC₄F₃ plants introgressed between SSR markers NAU3384 and CGR5144 were selected to develop sub-NILs. The qFL-chr1 was further mapped at 0.9-cM interval between MUSS422 and CIR018 by comparisons of sub-NIL phenotype, and increased fiber length by ~1 mm. The 2.38-Mb region between MUSS422 and CIR018 in G. barbadense contained 19 annotated genes. Expression levels of two of these genes, GOBAR07705 (encoding 1-aminocyclopropane-1-carboxylate synthase) and GOBAR25992 (encoding amino acid permease), were positively correlated with fiber length in a small F₂ population, supporting these genes as candidates for qFL-chr1.

High-density linkage map construction and QTL analysis for earliness-related traits in Gossypium hirsutum L

Background

Gossypium hirsutum L., or upland cotton, is an important renewable resource for textile fiber. To enhance understanding of the genetic basis of cotton earliness, we constructed an intra-specific recombinant inbred line population (RIL) containing 137 lines, and performed linkage map construction and quantitative trait locus (QTL) mapping.

Results

Using restriction-site associated DNA sequencing, a genetic map composed of 6,434 loci, including 6,295 single nucleotide polymorphisms and 139 simple sequence repeat loci, was developed from RIL population. This map spanned 4,071.98 cM, with an average distance of 0.63 cM between adjacent markers. A total of 247 QTLs for six earliness-related traits were detected in 6 consecutive years. In addition, 55 QTL coincidence regions representing more than 60 % of total QTLs were found on 22 chromosomes, which indicated that several earliness-related traits might be simultaneously improved. Fine-mapping of a 2-Mb region on chromosome D3 associated with five stable QTLs between Marker25958 and Marker25963 revealed that lines containing alleles derived from CCRI36 in this region exhibited smaller phenotypes and earlier maturity. One candidate gene (EMF2) was predicted and validated by quantitative real-time PCR in early-, medium- and late-maturing cultivars from 3- to 6-leaf stages, with highest expression level in early-maturing cultivar, CCRI74, lowest expression level in late-maturing cultivar, Bomian1.

Conclusions

We developed an SNP-based genetic map, and this map is the first high-density genetic map for short-season cotton and has the potential to provide deeper insights into earliness. Cotton earliness-related QTLs and QTL coincidence regions will provide useful materials for QTL fine mapping, gene positional cloning and MAS. And the gene, EMF2, is promising for further study.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3269-y) contains supplementary material, which is available to authorized users.

A Cotton Annexin Affects Fiber Elongation and Secondary Cell Wall Biosynthesis Associated with Ca2+ Influx, ROS Homeostasis, and Actin Filament Reorganization

Annexins play pivotal roles in a variety of cellular processes as well as in fiber development; however, the functional mechanisms of their activities are unclear. Here, an annexin gene that is preferentially expressed in fibers, GhFAnnxA, was found to be significantly associated with various cotton (Gossypium hirsutum) fiber traits. Transgenic analysis demonstrated that GhFAnnxA affected cotton fiber elongation and was involved in secondary cell wall (SCW) biosynthesis. Functional studies demonstrated that GhFAnnxA may act as a Ca(2+) conductance regulator and that reactive oxygen species (ROS) produced by Rbohs in a Ca(2+)-dependent manner may determine fiber elongation caused by elevated intracellular turgor and cell wall loosening. However, excessive hydrogen peroxide (H2O2) inhibited cotton fiber elongation in vitro. We speculate that a positive feedback loop involving ROS and Ca(2+) is regulated by GhCDPK1 and regulates fiber cell elongation. Furthermore, the convergence of actin filaments is altered by their interaction with GhFAnnxA, and this also may contribute to fiber elongation. Moreover, GhFAnnxA may affect SCW biosynthesis through changes in cell wall components caused by an increase in H2O2 levels. These results not only provide new insights into the signaling pathways of GhFAnnxA in fiber development but also clarify the role of ROS in fiber development.

A Genetic Map Between Gossypium hirsutum and the Brazilian Endemic G. mustelinum and Its Application to QTL Mapping

Among the seven tetraploid cotton species, little is known about transmission genetics and genome organization in Gossypium mustelinum, the species most distant from the source of most cultivated cotton, G. hirsutum. In this research, an F₂ population was developed from an interspecific cross between G. hirsutum and G. mustelinum (HM). A genetic linkage map was constructed mainly using simple sequence repeat (SSRs) and restriction fragment length polymorphism (RFLP) DNA markers. The arrangements of most genetic loci along the HM chromosomes were identical to those of other tetraploid cotton species. However, both major and minor structural rearrangements were also observed, for which we propose a parsimony-based model for structural divergence of tetraploid cottons from common ancestors. Sequences of mapped markers were used for alignment with the 26 scaffolds of the G. hirsutum draft genome, and showed high consistency. Quantitative trait locus (QTL) mapping of fiber elongation in advanced backcross populations derived from the same parents demonstrated the value of the HM map. The HM map will serve as a valuable resource for QTL mapping and introgression of G. mustelinum alleles into G. hirsutum, and help clarify evolutionary relationships between the tetraploid cotton genomes.

Detection and validation of one stable fiber strength QTL on c9 in tetraploid cotton

Fiber strength is an essential trait of fiber property in cotton, and it is quantitatively inherited. Identification of stable quantitative trait loci (QTL) contributing to fiber strength would provide the key basis for marker-assisted selection (MAS) in cotton breeding. In this study, four interspecific hybridization populations were established with a common G. barbadense parent Pima 90-53 and two G. hirsutum parents (CCRI 8 and Handan 208), each of which had fiber strength characteristic. Based on the phenotypic data of fiber strength from seven environments, a stable QTL, qFS-c9-1, was detected and validated on c9 in a marker interval between SSR markers NAU2395 and NAU1092. The QTL explaining 14.4-17.9 % of the phenotypic variation was firstly detected in two populations (CCRI 8 × Pima 90-53, BC1F1 and BC1F2) and its derived lines in four environments. And it accounting for 12.1-14.8 % of the phenotypic variation was further confirmed in two populations (Handan 208 × Pima 90-53, BC1F1, and F2) under one environment. In silico mapping using three sequenced cotton genomes indicated that homologous genes, anchored by NAU2395 and NAU1092, were aligned to the G. arboreum genome within a physical distance between 81.10 Mbps and 87.07 Mbps. In that interval, several genes were confirmed in literatures to associate with fiber development. Among these genes, seven genes were further selected for an expression analysis through fiber development transcriptome database, revealing unique expression patterns across different stages of fiber development between CCRI 8 and Pima 90-53. The genes underlying qFS-c9-1 were favorable to fine mapping and cloning. The current study results provided valuable evidence for mapping stable QTL of fiber strength utilizing multiple populations and environments, as well as map-based cloning the candidate gene underlying the QTL for future prospective research directions.

Structure, evolution, and comparative genomics of tetraploid cotton based on a high-density genetic linkage map

A high-density linkage map was constructed using 1,885 newly obtained loci and 3,747 previously published loci, which included 5,152 loci with 4696.03 cM in total length and 0.91 cM in mean distance. Homology analysis in the cotton genome further confirmed the 13 expected homologous chromosome pairs and revealed an obvious inversion on Chr10 or Chr20 and repeated inversions on Chr07 or Chr16. In addition, two reciprocal translocations between Chr02 and Chr03 and between Chr04 and Chr05 were confirmed. Comparative genomics between the tetraploid cotton and the diploid cottons showed that no major structural changes exist between D_T and D chromosomes but rather between A_T and A chromosomes. Blast analysis between the tetraploid cotton genome and the mixed genome of two diploid cottons showed that most AD chromosomes, regardless of whether it is from the A_T or D_T genome, preferentially matched with the corresponding homologous chromosome in the diploid A genome, and then the corresponding homologous chromosome in the diploid D genome, indicating that the diploid D genome underwent converted evolution by the diploid A genome to form the D_T genome during polyploidization. In addition, the results reflected that a series of chromosomal translocations occurred among Chr01/Chr15, Chr02/Chr14, Chr03/Chr17, Chr04/Chr22, and Chr05/Chr19.

Genome Wide SSR High Density Genetic Map Construction from an Interspecific Cross of Gossypium hirsutum × Gossypium tomentosum

A high density genetic map was constructed using F2 population derived from an interspecific cross of G. hirsutum × G. tomentosum. The map consisted of 3093 marker loci distributed across all the 26 chromosomes and covered 4365.3 cM of cotton genome with an average inter-marker distance of 1.48 cM. The maximum length of chromosome was 218.38 cM and the minimum was 122.09 cM with an average length of 167.90 cM. A sub-genome covers more genetic distance (2189.01 cM) with an average inter loci distance of 1.53 cM than D sub-genome which covers a length of 2176.29 cM with an average distance of 1.43 cM. There were 716 distorted loci in the map accounting for 23.14% and most distorted loci were distributed on D sub-genome (25.06%), which were more than on A sub-genome (21.23%). In our map 49 segregation hotspots (SDR) were distributed across the genome with more on D sub-genome as compared to A genome. Two post-polyploidization reciprocal translocations of "A2/A3 and A4/A5" were suggested by seven pairs of duplicate loci. The map constructed through these studies is one of the three densest genetic maps in cotton however; this is the first dense genome wide SSR interspecific genetic map between G. hirsutum and G. tomentosum.

Functional divergence of GhCFE5 homoeologs revealed in cotton fiber and Arabidopsis root cell development

In GhCFE5 homoeologs, GhCFE5D interacted with more actin homologs and stronger interaction activity than GhCFE5A. GhCFE5D - but not GhCFE5A -overexpression severely disrupted actin cytoskeleton organization and significantly suppressed cell elongation. Homoeologous genes are common in polyploid plants; however, their functional divergence is poorly elucidated. Allotetraploid Upland cotton (Gossypium hirsutum, AADD) is the most widely cultivated cotton; accounting for more than 90 % of the world's cotton production. Here, we characterized GhCFE5A and GhCFE5D homoeologs from G. hirsutum acc TM-1. GhCFE5 homoeologs are expressed preferentially in fiber cells; and a significantly greater accumulation of GhCFE5A mRNA than GhCFE5D mRNA was found in all tested tissues. Overexpression of GhCFE5D but not GhCFE5A seriously inhibits the Arabidopsis hypocotyl and root cell elongation. Yeast two-hybrid assay and bimolecular fluorescence complementation (BiFC) analysis showed that compared with GhCFE5A, GhCFE5D interacts with more actin homologs and has a stronger interaction activity both from Arabidopsis and Upland cotton. Interestingly, subcellular localization showed that GhCFE5 resides on the cortical endoplasmic reticulum (ER) network and is colocalized with actin cables. The interaction activities between GhCFE5 homoeologs and actin differ in their effects on F-actin structure in transgenic Arabidopsis root cells. The F-actin changed direction from vertical to lateral, and the actin cytoskeleton organization was severely disrupted in GhCFE5D-overexpressing root cells. These data support the functional divergence of GhCFE5 homoeologs in the actin cytoskeleton structure and cell elongation, implying an important role for GhCFE5 in the evolution and selection of cotton fiber.

Evaluation and Exploration of Favorable QTL Alleles for Salt Stress Related Traits in Cotton Cultivars (G. hirsutum L.)

Soil salinization is one of the major problems in global agricultural production. Cotton is a pioneer crop with regard to salt stress tolerance, and can be used for saline-alkali land improvement. The large-scale detection of salt tolerance traits in cotton accessions, and the identification of elite quantitative trait loci (QTLs)/genes for salt-tolerance have been very important in salt tolerance breeding. Here, 43 advanced salt-tolerant and 31 highly salt-sensitive cultivars were detected by analyzing ten salt tolerance related traits in 304 upland cotton cultivars. Among them, 11 advanced salt-tolerance and eight highly salt-sensitive cultivars were consistent with previously reported results. Association analysis of ten salt-tolerance related traits and 145 SSRs was performed, and a total of 95 significant associations were detected; 17, 41, and 37 of which were associated with germinative index, seedling stage physiological index, and four seedling stage biochemical indexes, respectively. Of these associations, 20 SSR loci were simultaneously associated with two or more traits. Furthermore, we detected 117 elite alleles associated with salt-tolerance traits, 4 of which were reported previously. Among these loci, 44 (37.60%) were rare alleles with a frequency of less than 5%, 6 only existed in advanced salt-tolerant cultivars, and 2 only in highly salt-sensitive cultivars. As a result, 13 advanced salt-tolerant cultivars were selected to assemble the optimal cross combinations by computer simulation for the development of salt-tolerant accessions. This study lays solid foundations for further improvements in cotton salt-tolerance by referencing elite germplasms, alleles associated with salt-tolerance traits, and optimal crosses.

Partial Dominance, Overdominance and Epistasis as the Genetic Basis of Heterosis in Upland Cotton (Gossypium hirsutum L.)

Determination of genetic basis of heterosis may promote hybrid production in Upland cotton (Gossypium hirsutum L.). This study was designed to explore the genetic mechanism of heterosis for yield and yield components in F2: 3 and F2: 4 populations derived from a hybrid 'Xinza No. 1'. Replicated yield field trials of the progenies were conducted in 2008 and 2009. Phenotypic data analyses indicated overdominance in F1 for yield and yield components. Additive and dominance effects at single-locus level and digenic epistatic interactions at two-locus level were analyzed by 421 marker loci spanning 3814 cM of the genome. A total of 38 and 49 QTLs controlling yield and yield components were identified in F2: 3 and F2: 4 populations, respectively. Analyses of these QTLs indicated that the effects of partial dominance and overdominance contributed to heterosis in Upland cotton simultaneously. Most of the QTLs showed partial dominance whereas 13 QTLs showing overdominance in F2:3 population, and 19 QTLs showed overdominance in F2:4. Among them, 21 QTLs were common in both F2: 3 and F2: 4 populations. A large number of two-locus interactions for yield and yield components were detected in both generations. AA (additive × additive) epistasis accounted for majority portion of epistatic effects. Thirty three complementary two-locus homozygotes (11/22 and 22/11) were the best genotypes for AA interactions in terms of bolls per plant. Genotypes of double homozygotes, 11/22, 22/11 and 22/22, performed best for AD/DA interactions, while genotype of 11/12 performed best for DD interactions. These results indicated that (1) partial dominance and overdominance effects at single-locus level and (2) epistasis at two-locus level elucidated the genetic basis of heterosis in Upland cotton.

Sequence composition of BAC clones and SSR markers mapped to Upland cotton chromosomes 11 and 21 targeting resistance to soil-borne pathogens

Genetic and physical framework mapping in cotton (Gossypium spp.) were used to discover putative gene sequences involved in resistance to common soil-borne pathogens. Chromosome (Chr) 11 and its homoeologous Chr 21 of Upland cotton (G. hirsutum) are foci for discovery of resistance (R) or pathogen-induced R (PR) genes underlying QTLs involved in response to root-knot nematode (Meloidogyne incognita), reniform nematode (Rotylenchulus reniformis), Fusarium wilt (Fusarium oxysporum f.sp. vasinfectum), Verticillium wilt (Verticillium dahliae), and black root rot (Thielaviopsis basicola). Simple sequence repeat (SSR) markers and bacterial artificial chromosome (BAC) clones from a BAC library developed from the Upland cotton Acala Maxxa were mapped on Chr 11 and Chr 21. DNA sequence through Gene Ontology (GO) of 99 of 256 Chr 11 and 109 of 239 Chr 21 previously mapped SSRs revealed response elements to internal and external stimulus, stress, signaling process, and cell death. The reconciliation between genetic and physical mapping of gene annotations from new DNA sequences of 20 BAC clones revealed 467 (Chr 11) and 285 (Chr 21) G. hirsutum putative coding sequences, plus 146 (Chr 11) and 98 (Chr 21) predicted genes. GO functional profiling of Unigenes uncovered genes involved in different metabolic functions and stress response elements (SRE). Our results revealed that Chrs 11 and 21 harbor resistance gene rich genomic regions. Sequence comparisons with the ancestral diploid D5 (G. raimondii), A2 (G. arboreum) and domesticated tetraploid TM-1 AD1 (G. hirsutum) genomes revealed abundance of transposable elements and confirmed the richness of resistance gene motifs in these chromosomes. The sequence information of SSR markers and BAC clones and the genetic mapping of BAC clones provide enhanced genetic and physical frameworks of resistance gene-rich regions of the cotton genome, thereby aiding discovery of R and PR genes and breeding for resistance to cotton diseases.

A high-density SSR genetic map constructed from a F2 population of Gossypium hirsutum and Gossypium darwinii

The cultivated allotetraploid species Gossypium hirsutum, accounts for 90% of the world cotton production, has narrow genetic basis that's why its yield, quality or stress resistance breeding is stagnant. It is therefore, essential to explore desirable genes from Gossypium darwinii which has enviable traits such as high fiber fineness, drought tolerance, fusarium and verticillium resistance. We used G. darwinii as primary plant materials in this study not only to enrich the genetic diversity of exiting germplasm but also to better understand its genome structure. An interspecific high density linkage map of allotetraploid cotton was constructed using F2 population (G. hirsutum×G. darwinii). The map was based entirely on genome-wide simple sequence repeat (SSR) markers. A total of 2763 markers were mapped in 26 linkage groups (chromosomes) covering a genome length of 4176.7cM with an average inter-locus distance of 1.5cM. The length of the chromosomes ranged from 84.7 to 238.5cM with an average length of 160.6cM. At subgenome length was 2160.7cM with an average distance of 1.6cM, where as Dt genome length was 2016cM with an average distance of 1.4cM. There were 601 distorted SSR loci. Less number of segregation distortion loci were located in At subgenome than in Dt subgenome. Two post-polyploidization reciprocal translocations of "A2/A3 and A4/A5" were suggested by 44 pairs of duplicate loci.

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.

The Basic/Helix-Loop-Helix Protein Family in Gossypium: Reference Genes and Their Evolution during Tetraploidization

Basic/helix-loop-helix (bHLH) proteins comprise one of the largest transcription factor families and play important roles in diverse cellular and molecular processes. Comprehensive analyses of the composition and evolution of the bHLH family in cotton are essential to elucidate their functions and the molecular basis of cotton development. By searching bHLH homologous genes in sequenced diploid cotton genomes (Gossypium raimondii and G. arboreum), a set of cotton bHLH reference genes containing 289 paralogs were identified and named as GobHLH001-289. Based on their phylogenetic relationships, these cotton bHLH proteins were clustered into 27 subfamilies. Compared to those in Arabidopsis and cacao, cotton bHLH proteins generally increased in number, but unevenly in different subfamilies. To further uncover evolutionary changes of bHLH genes during tetraploidization of cotton, all genes of S5a and S5b subfamilies in upland cotton and its diploid progenitors were cloned and compared, and their transcript profiles were determined in upland cotton. A total of 10 genes of S5a and S5b subfamilies (doubled from A- and D-genome progenitors) maintained in tetraploid cottons. The major sequence changes in upland cotton included a 15-bp in-frame deletion in GhbHLH130D and a long terminal repeat retrotransposon inserted in GhbHLH062A, which eliminated GhbHLH062A expression in various tissues. The S5a and S5b bHLH genes of A and D genomes (except GobHLH062) showed similar transcription patterns in various tissues including roots, stems, leaves, petals, ovules, and fibers, while the A- and D-genome genes of GobHLH110 and GobHLH130 displayed clearly different transcript profiles during fiber development. In total, this study represented a genome-wide analysis of cotton bHLH family, and revealed significant changes in sequence and expression of these genes in tetraploid cottons, which paved the way for further functional analyses of bHLH genes in the cotton genus.

Sequencing of allotetraploid cotton (Gossypium hirsutum L. acc. TM-1) provides a resource for fiber improvement

Upland cotton is a model for polyploid crop domestication and transgenic improvement. Here we sequenced the allotetraploid Gossypium hirsutum L. acc. TM-1 genome by integrating whole-genome shotgun reads, bacterial artificial chromosome (BAC)-end sequences and genotype-by-sequencing genetic maps. We assembled and annotated 32,032 A-subgenome genes and 34,402 D-subgenome genes. Structural rearrangements, gene loss, disrupted genes and sequence divergence were more common in the A subgenome than in the D subgenome, suggesting asymmetric evolution. However, no genome-wide expression dominance was found between the subgenomes. Genomic signatures of selection and domestication are associated with positively selected genes (PSGs) for fiber improvement in the A subgenome and for stress tolerance in the D subgenome. This draft genome sequence provides a resource for engineering superior cotton lines.

Molecular Mapping of Restriction-Site Associated DNA Markers In Allotetraploid Upland Cotton

Upland cotton (Gossypium hirsutum L., 2n = 52, AADD) is an allotetraploid, therefore the discovery of single nucleotide polymorphism (SNP) markers is difficult. The recent emergence of genome complexity reduction technologies based on the next-generation sequencing (NGS) platform has greatly expedited SNP discovery in crops with highly repetitive and complex genomes. Here we applied restriction-site associated DNA (RAD) sequencing technology for de novo SNP discovery in allotetraploid cotton. We identified 21,109 SNPs between the two parents and used these for genotyping of 161 recombinant inbred lines (RILs). Finally, a high dense linkage map comprising 4,153 loci over 3500-cM was developed based on the previous result. Using this map quantitative trait locus (QTLs) conferring fiber strength and Verticillium Wilt (VW) resistance were mapped to a more accurate region in comparison to the 1576-cM interval determined using the simple sequence repeat (SSR) genetic map. This suggests that the newly constructed map has more power and resolution than the previous SSR map. It will pave the way for the rapid identification of the marker-assisted selection in cotton breeding and cloning of QTL of interest traits.

GhCFE1A, a dynamic linker between the ER network and actin cytoskeleton, plays an important role in cotton fibre cell initiation and elongation

Fibre cell initiation and elongation is critical for cotton fibre development. However, little is known about the regulation of initiation and elongation during fibre cell development. Here, the regulatory role of a novel protein GhCFE1A was uncovered. GhCFE1A is preferentially expressed at initiation and rapid elongation stages during fibre development; in addition, much higher expression of GhCFE1A was detected at the fibre initiation stage in fibreless cotton mutants than in the fibre-bearing TM-1 wild-type. Importantly, overexpression of GhCFE1A in cotton not only delayed fibre cell elongation but also significantly reduced the density of lint and fuzz fibre initials and stem trichomes. Yeast two-hybrid assay showed that GhCFE1A interacted with several actin proteins, and the interaction was further confirmed by co-sedimentation assay. Interestingly, a subcellular localization assay showed that GhCFE1A resided on the cortical endoplasmic reticulum (ER) network and co-localized with actin cables. Moreover, the density of F-actin filaments was shown to be reduced in GhCFE1A-overexpressing fibres at the rapid elongation stage compared with the wild-type control. Taken together, the results demonstrate that GhCFE1A probably functions as a dynamic linker between the actin cytoskeleton and the ER network, and plays an important role in fibre cell initiation and elongation during cotton fibre development.

Development of chromosome-specific markers with high polymorphism for allotetraploid cotton based on genome-wide characterization of simple sequence repeats in diploid cottons (Gossypium arboreum L. and Gossypium raimondii Ulbrich)

Background

Tetraploid cotton contains two sets of homologous chromosomes, the At- and Dt-subgenomes. Consequently, many markers in cotton were mapped to multiple positions during linkage genetic map construction, posing a challenge to anchoring linkage groups and mapping economically-important genes to particular chromosomes. Chromosome-specific markers could solve this problem. Recently, the genomes of two diploid species were sequenced whose progenitors were putative contributors of the At- and Dt-subgenomes to tetraploid cotton. These sequences provide a powerful tool for developing chromosome-specific markers given the high level of synteny among tetraploid and diploid cotton genomes. In this study, simple sequence repeats (SSRs) on each chromosome in the two diploid genomes were characterized. Chromosome-specific SSRs were developed by comparative analysis and proved to distinguish chromosomes.

Results

A total of 200,744 and 142,409 SSRs were detected on the 13 chromosomes of Gossypium arboreum L. and Gossypium raimondii Ulbrich, respectively. Chromosome-specific SSRs were obtained by comparing SSR flanking sequences from each chromosome with those from the other 25 chromosomes. The average was 7,996 per chromosome. To confirm their chromosome specificity, these SSRs were used to distinguish two homologous chromosomes in tetraploid cotton through linkage group construction. The chromosome-specific SSRs and previously-reported chromosome markers were grouped together, and no marker mapped to another homologous chromosome, proving that the chromosome-specific SSRs were unique and could distinguish homologous chromosomes in tetraploid cotton. Because longer dinucleotide AT-rich repeats were the most polymorphic in previous reports, the SSRs on each chromosome were sorted by motif type and repeat length for convenient selection. The primer sequences of all chromosome-specific SSRs were also made publicly available.

Conclusion

Chromosome-specific SSRs are efficient tools for chromosome identification by anchoring linkage groups to particular chromosomes during genetic mapping and are especially useful in mapping of qualitative-trait genes or quantitative trait loci with just a few markers. The SSRs reported here will facilitate a number of genetic and genomic studies in cotton, including construction of high-density genetic maps, positional gene cloning, fingerprinting, and genetic diversity and comparative evolutionary analyses among Gossypium species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1265-2) contains supplementary material, which is available to authorized users.

Enrichment of an intraspecific genetic map of upland cotton by developing markers using parental RAD sequencing

RAD sequencing was performed using DH962 and Jimian5 as upland cotton mapping parents. Sequencing data for DH962 and Jimian5 were assembled into the genome sequences of ≈55.27 and ≈57.06 Mb, respectively. Analysing genome sequences of the two parents, 1,323 SSR, 3,838 insertion/deletion (InDel), and 9,366 single-nucleotide polymorphism (SNP) primer pairs were developed. All of the SSRs, 121 InDels, 441 SNPs, and other 6,747 primer pairs were screened in the two parents, and a total of 535 new polymorphic loci were identified. A genetic map including 1,013 loci was constructed using these results and 506 loci previously published for this population. Twenty-seven new QTLs for yield and fibre quality were identified, indicating that the efficiency of QTL detection was greatly improved by the increase in map density. Comparative genomics showed there to be considerable homology and collinearity between the A_T and A₂ genomes and between the D_T and D₅ genomes, although there were a few exchanges and introgressions among the chromosomes of the A₂ genome. Here, the development of markers using parental RAD sequencing was effective, and a high-density intraspecific genetic map was constructed. This map can be used for molecular marker-assisted selection in cotton.

Molecular characterization of the Gossypium Diversity Reference Set of the US National Cotton Germplasm Collection

A core marker set containing markers developed to be informative within a single commercial cotton species can elucidate diversity structure within a multi-species subset of the Gossypium germplasm collection. An understanding of the genetic diversity of cotton (Gossypium spp.) as represented in the US National Cotton Germplasm Collection is essential to develop strategies for collecting, conserving, and utilizing these germplasm resources. The US collection is one of the largest world collections and includes not only accessions with improved yield and fiber quality within cultivated species, but also accessions possessing sources of abiotic and biotic stress resistance often found in wild species. We evaluated the genetic diversity of a subset of 272 diploid and 1,984 tetraploid accessions in the collection (designated the Gossypium Diversity Reference Set) using a core set of 105 microsatellite markers. Utility of the core set of markers in differentiating intra-genome variation was much greater in commercial tetraploid genomes (99.7 % polymorphic bands) than in wild diploid genomes (72.7 % polymorphic bands), and may have been influenced by pre-selection of markers for effectiveness in the commercial species. Principal coordinate analyses revealed that the marker set differentiated interspecific variation among tetraploid species, but was only capable of partially differentiating among species and genomes of the wild diploids. Putative species-specific marker bands in G. hirsutum (73) and G. barbadense (81) were identified that could be used for qualitative identification of misclassifications, redundancies, and introgression within commercial tetraploid species. The results of this broad-scale molecular characterization are essential to the management and conservation of the collection and provide insight and guidance in the use of the collection by the cotton research community in their cotton improvement efforts.

Identification of centromeric regions on the linkage map of cotton using centromere-related repeats

Centromere usually contains high-copy-number retrotransposons and satellite repeats, which are difficult to map, clone and sequence. Currently, very little is known about the centromere in cotton. Here, we sequenced a bacterial artificial chromosome (BAC) mapping to the centromeric region and predicted four long-terminal-repeat (LTR) retrotransposons. They were located in the heterochromatic centromeric regions of all 52 pachytene chromosomes in Gossypium hirsutum. Fiber-FISH mapping revealed that these retrotransposons span an area of at least 1.8Mb in the centromeric region. Comparative analysis showed that these retrotransposons generated similar, strong fluorescent signals in the D progenitor Gossypium raimondii but not in the A progenitor Gossypium herbaceum, suggesting that the centromere sequence of tetraploid cotton might be derived from the D progenitor. Centromeric regions were anchored on 13 chromosomes of D-genome sequence. Characterization of these centromere-related repeats and regions will enhance cotton centromere mapping, sequencing and evolutionary studies.

Characterization of the global transcriptome for cotton (Gossypium hirsutum L.) anther and development of SSR marker

Cotton is an important fiber plant, and it's attractive to elucidate the molecular mechanism of anther development due to the close relationship between the anther fertility and boll-setting, and also fiber yield. In the present paper, 47.2 million paired-end reads with average length of 82.87 bp from the anthers of TM-1 (Gossypium hirsutum L.), a genetic standard line, were generated through transcriptome sequencing, and 210,965 unigenes of more than 100 bp were obtained. BLAST, KEGG, COG, and GO analyses showed that the genes were enriched in the processes of transcription, translation, and post-translation as well as hormone signal transduction, the transcription factor families, and cell wall-related genes mainly participating in cell expansion and carbohydrate metabolism. Further analysis identified 11,153 potential SSRs. A suit of 5122 primer pair sequences were designed, and 82 of 300 randomly selected primer pairs produced reproducible amplicons that were polymorphic among 22 cotton accessions from G. hirsutum, Gossypium barbadense and Gossypium arboreum. The UPGMA clustering analysis further confirmed high quality and effectiveness of these novel SSR markers. The present study provided insights into the transcriptome profile of the cotton and established a public information platform for functional genomics and molecular breeding.

Study on the activity of the signaling pathways regulating hepatocytes from G0 phase into G1 phase during rat liver regeneration

Under normal physiological conditions, the majority of hepatocytes are in the functional state (G0 phase). After injury or liver partial hepatectomy (PH), hepatocytes are rapidly activated to divide. To understand the mechanism underlying hepatocyte G0/G1 transition during rat liver regeneration, we used the Rat Genome 230 2.0 Array to determine the expression changes of genes, then searched the GO and NCBI databases for genes associated with the G0/G1 transition, and QIAGEN and KEGG databases for the G0/G1 transition signaling pathways. We used expression profile function (E t ) to calculate the activity level of the known G0/G1 transition signal pathways, and Ingenuity Pathway Analysis 9.0 (IPA) to determine the interactions among these signaling pathways. The results of our study show that the activity of the signaling pathways of HGF, IL-10 mediated by p38MAPK, IL-6 mediated by STAT3, and JAK/STAT mediated by Ras/ERK and STAT3 are significantly increased during the priming phase (2-6 h after PH) of rat liver regeneration. This leads us to conclude that during rat liver regeneration, the HGF, IL-10, IL-6 and JAK/STAT signaling pathways play a major role in promoting hepatocyte G0/G1 transition in the regenerating liver.

A rapid and effective method for silver staining of PCR products separated in polyacrylamide gels

With the development of molecular quantitative genetics, particularly, genetic linkage map construction, quantitative trait loci mapping or genes fine mapping and association analysis etc., more and more PCR products separated in polyacrylamide gels need to be silver-stained. However, conventional silver-staining procedures are complicated and time-consuming as they require a lot of preparation and handling of several solutions prior to use. In this study, a simple and rapid protocol for silver staining of PCR products was developed. The number of steps was reduced compared to conventional protocols, thus achieving detection of PCR products in 7 min, saving time and resources. Fixation and staining solution and developing solution in present staining procedure allowed a reutilization for 12 and 8 times, respectively, reducing the cost greatly. Meanwhile, the sensitivity was significantly improved with the improved method and the minimum of 0.097 ng/μL of DNA amount can be detected in denaturing polyacrylamide gel. The protocol developed in this study will facilitate the development of molecular quantitative genetics.

Molecular mapping and validation of a major QTL conferring resistance to a defoliating isolate of verticillium wilt in cotton (Gossypium hirsutum L.)

Verticillium wilt (VW) caused by Verticillium dahliae Kleb is one of the most destructive diseases of cotton. Development and use of a VW resistant variety is the most practical and effective way to manage this disease. Identification of highly resistant genes/QTL and the underlining genetic architecture is a prerequisite for developing a VW resistant variety. A major QTL qVW-c6-1 conferring resistance to the defoliating isolate V991 was identified on chromosome 6 in LHB22×JM11 F2∶3 population inoculated and grown in a greenhouse. This QTL was further validated in the LHB22×NNG F2∶3 population that was evaluated in an artificial disease nursery of V991 for two years and in its subsequent F4 population grown in a field severely infested by V991. The allele conferring resistance within the QTL qVW-c6-1 region originated from parent LHB22 and could explain 23.1-27.1% of phenotypic variation. Another resistance QTL qVW-c21-1 originated from the susceptible parent JM11 was mapped on chromosome 21, explaining 14.44% of phenotypic variation. The resistance QTL reported herein provides a useful tool for breeding a cotton variety with enhanced resistance to VW.

SSR marker-assisted improvement of fiber qualities in Gossypium hirsutum using G. barbadense introgression lines

This study demonstrates the first practical use of CSILs for the transfer of fiber quality QTLs into Upland cotton cultivars using SSR markers without detrimentally affecting desirable agronomic characteristics. Gossypium hirsutum is characterized by its high lint production and medium fiber quality compared to extra-long staple cotton G. barbadense. Transferring valuable traits or genes from G. barbadense into G. hirsutum is a promising but challenging approach through a traditional interspecific introgression strategy. We developed one set of chromosome segment introgression lines (CSILs), where TM-1, the genetic standard in G. hirsutum, was used as the recipient parent and the long staple cotton G. barbadense cv. Hai7124 was used as the donor parent by molecular marker-assisted selection (MAS). Among them, four CSILs, IL040-A4-1, IL080-D6-1, IL088-A7-3 and IL019-A2-6, found to be associated with superior fiber qualities including fiber length, strength and fineness QTL in Xinjiang were selected and backcrossed, and transferred these QTLs into three commercial Upland cotton cultivars such as Xinluzao (XLZ) 26, 41 and 42 grown in Xinjiang. By backcrossing and self-pollinating twice, five improved lines (3262-4, 3389-2, 3326-3, 3380-4 and 3426-5) were developed by MAS of background and introgressed segments. In diverse field trials, these QTLs consistently and significantly offered additive effects on the target phenotype. Furthermore, we also pyramided two segments from different CSILs (IL080-D6-1 and IL019-A2-6) into cultivar 0768 to accelerate breeding process purposefully with MAS. The improved lines pyramided by these two introgressed segments showed significant additive epistatic effects in four separate field trials. No significant alteration in yield components was observed in these modified lines. In summary, we first report that these CSILs have great potential to improve fiber qualities in Upland cotton MAS breeding programs.

Association analysis of fiber quality traits and exploration of elite alleles in Upland cotton cultivars/accessions (Gossypium hirsutum L.)

Exploring the elite alleles and germplasm accessions related to fiber quality traits will accelerate the breeding of cotton for fiber quality improvement. In this study, 99 Gossypium hirsutum L. accessions with diverse origins were used to perform association analysis of fiber quality traits using 97 polymorphic microsatellite marker primer pairs. A total of 107 significant marker-trait associations were detected for three fiber quality traits under three different environments, with 70 detected in two or three environments and 37 detected in only one environment. Among the 70 significant marker-trait associations, 52.86% were reported previously, implying that these are stable loci for target traits. Furthermore, we detected a large number of elite alleles associated simultaneously with two or three traits. These elite alleles were mainly from accessions collected in China, introduced to China from the United States, or rare alleles with a frequency of less than 5%. No one cultivar contained more than half of the elite alleles, but 10 accessions were collected from China and the two introduced from the United States did contain more than half of these alleles. Therefore, there is great potential for mining elite alleles from germplasm accessions for use in fiber quality improvement in modern cotton breeding.

First insights into the large genome of Epimedium sagittatum (Sieb. et Zucc) Maxim, a Chinese Ttaditional medicinal plant

Epimedium sagittatum (Sieb. et Zucc) Maxim is a member of the Berberidaceae family of basal eudicot plants, widely distributed and used as a traditional medicinal plant in China for therapeutic effects on many diseases with a long history. Recent data shows that E. sagittatum has a relatively large genome, with a haploid genome size of ~4496 Mbp, divided into a small number of only 12 diploid chromosomes (2n = 2x = 12). However, little is known about Epimedium genome structure and composition. Here we present the analysis of 691 kb of high-quality genomic sequence derived from 672 randomly selected plasmid clones of E. sagittatum genomic DNA, representing ~0.0154% of the genome. The sampled sequences comprised at least 78.41% repetitive DNA elements and 2.51% confirmed annotated gene sequences, with a total GC% content of 39%. Retrotransposons represented the major class of transposable element (TE) repeats identified (65.37% of all TE repeats), particularly LTR (Long Terminal Repeat) retrotransposons (52.27% of all TE repeats). Chromosome analysis and Fluorescence in situ Hybridization of Gypsy-Ty3 retrotransposons were performed to survey the E. sagittatum genome at the cytological level. Our data provide the first insights into the composition and structure of the E. sagittatum genome, and will facilitate the functional genomic analysis of this valuable medicinal plant.

Discovery and identification of a novel Ligon lintless-like mutant (Lix) similar to the Ligon lintless (Li1) in allotetraploid cotton

Mutants are a powerful resource for studying gene structure, function, and evolution. In this present study, a novel Ligon lintless-like mutant (Lix), that has short fibers and deformed leaves and stems, was isolated from the progeny of transgenic cottons. The Lix mutant is similar in morphology to the Ligon lintless (Li1) mutant. Genetic analysis and molecular mapping were performed for the Lix and Li1 mutants. These two mutants are monogenic dominant mutants with distorted growth of vegetative and reproductive structures. Seedlings of the dominant homozygote Li 1 Li 1 genotype are lethal, while LixLix plants are viable but show no reproductive growth. Molecular tagging showed that the Lix gene is located on Chr. 04 in a 30.9-cM region spanned by NAU8376 and NAU3469. In a previous study, the Li 1 gene was mapped to Chr. 22, and Chr. 04 and Chr. 22 are homoelogous chromosomes in tetraploid cotton. So, we propose that Lix and Li1 mutants have similar mutated morphology, and Lix is mapped to a homoelogous chromosome carrying Li 1 . The identification and genetic mapping of Lix/Li 1 genes using mutants provides a foundation for isolating these genes. In turn, this will permit studies to elucidate the functional and evolutionary roles for these genes in cotton growth and development.

Characterization of expressed sequence tags from developing fibers of Gossypium barbadense and evaluation of insertion-deletion variation in tetraploid cultivated cotton species

Background

Cotton is the leading fiber crop worldwide. Gossypium barbadense is an important species of cotton because of its extra-long staple fibers with superior luster and silkiness. However, a systematic analysis and utilization of cDNA sequences from G. barbadense fiber development remains understudied.

Results

A total of 21,079 high quality sequences were generated from two non-normalized cDNA libraries prepared by using a mixture of G. barbadense Hai7124 fibers and ovules. After assembly processing, a set of 8,653 unigenes were obtained. Of those, 7,786 were matched to known proteins and 7,316 were assigned to functional categories. The molecular functions of these unigenes were mostly related to binding and catalytic activity, and carbohydrate, amino acid, and energy metabolisms were major contributors among the subsets of metabolism. Sequences comparison between G. barbadense and G. hirsutum revealed that 8,245 unigenes from G. barbadense were detected the similarity with those released publicly in G. hirsutum, however, the remaining 408 sequences had no hits against G. hirsutum unigenes database. Furthermore, 13,275 putative ESTs InDels loci involved in the orthologous and/or homoeologous differences between/within G. barbadense and G. hirsutum were discovered by in silico analyses, and 2,160 InDel markers were developed by ESTs with more than five insertions or deletions. By gel electrophoresis combined with sequencing verification, 71.11% candidate InDel loci were reconfirmed orthologous and/or homoeologous loci polymorphisms using G. hirsutum acc TM-1 and G. barbadense cv Hai7124. Blastx result showed among 2,160 InDel loci, 81 with significant function similarity with known genes associated with secondary wall synthesis process, indicating the important roles in fiber quality in tetraploid cultivated cotton species.

Conclusion

Sequence comparisons and InDel markers development will lay the groundwork for promoting the identification of genes related to superior agronomic traits, genetic differentiation and comparative genomic studies between G. hirsutum and G. barbadense.

Molecular markers associated with the immature fiber (im) gene affecting the degree of fiber cell wall thickening in cotton (Gossypium hirsutum L.)

Cotton fiber fineness and maturity measured indirectly as micronaire (MIC) are important properties of determining fiber grades in the textile market. To understand the genetic control and molecular mechanisms of fiber fineness and maturity, we studied two near isogenic lines, Gossypium hirsutum, Texas Marker-1 wild type (TM-1) and immature fiber (im) mutant showing a significant difference in MIC values. The fibers from im mutant plants were finer and less mature with lower MIC values than those from the recurrent parent, TM-1. A comprehensive fiber property analysis of TM-1 and im mutant showed that the lower MIC of fibers in im mutant was due to the lower degree of fiber cell wall thickening as compared to the TM-1 fibers. Using an F(2) population comprising 366 progenies derived from a cross between TM-1 and im mutant, we confirmed that the immature fiber phenotype present in a mutant plant was controlled by one single recessive gene im. Furthermore, we identified 13 simple sequence repeat markers that were closely linked to the im gene located on chromosome 3. Molecular markers associated with the im gene will lay the foundation to further investigate genetic information required for improving cotton fiber fineness and maturity.

Toward allotetraploid cotton genome assembly: integration of a high-density molecular genetic linkage map with DNA sequence information

Background

Cotton is the world’s most important natural textile fiber and a significant oilseed crop. Decoding cotton genomes will provide the ultimate reference and resource for research and utilization of the species. Integration of high-density genetic maps with genomic sequence information will largely accelerate the process of whole-genome assembly in cotton.

Results

In this paper, we update a high-density interspecific genetic linkage map of allotetraploid cultivated cotton. An additional 1,167 marker loci have been added to our previously published map of 2,247 loci. Three new marker types, InDel (insertion-deletion) and SNP (single nucleotide polymorphism) developed from gene information, and REMAP (retrotransposon-microsatellite amplified polymorphism), were used to increase map density. The updated map consists of 3,414 loci in 26 linkage groups covering 3,667.62 cM with an average inter-locus distance of 1.08 cM. Furthermore, genome-wide sequence analysis was finished using 3,324 informative sequence-based markers and publicly-available Gossypium DNA sequence information. A total of 413,113 EST and 195 BAC sequences were physically anchored and clustered by 3,324 sequence-based markers. Of these, 14,243 ESTs and 188 BACs from different species of Gossypium were clustered and specifically anchored to the high-density genetic map. A total of 2,748 candidate unigenes from 2,111 ESTs clusters and 63 BACs were mined for functional annotation and classification. The 337 ESTs/genes related to fiber quality traits were integrated with 132 previously reported cotton fiber quality quantitative trait loci, which demonstrated the important roles in fiber quality of these genes. Higher-level sequence conservation between different cotton species and between the A- and D-subgenomes in tetraploid cotton was found, indicating a common evolutionary origin for orthologous and paralogous loci in Gossypium.

Conclusion

This study will serve as a valuable genomic resource for tetraploid cotton genome assembly, for cloning genes related to superior agronomic traits, and for further comparative genomic analyses in Gossypium.

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.

Genome-wide identification of differentially expressed genes under water deficit stress in upland cotton (Gossypium hirsutum L.)

BACKGROUND

Cotton is the world's primary fiber crop and is a major agricultural commodity in over 30 countries. Like many other global commodities, sustainable cotton production is challenged by restricted natural resources. In response to the anticipated increase of agricultural water demand, a major research direction involves developing crops that use less water or that use water more efficiently. In this study, our objective was to identify differentially expressed genes in response to water deficit stress in cotton. A global expression analysis using cDNA-Amplified Fragment Length Polymorphism was conducted to compare root and leaf gene expression profiles from a putative drought resistant cotton cultivar grown under water deficit stressed and well watered field conditions.

RESULTS

We identified a total of 519 differentially expressed transcript derived fragments. Of these, 147 transcript derived fragment sequences were functionally annotated according to their gene ontology. Nearly 70 percent of transcript derived fragments belonged to four major categories: 1) unclassified, 2) stress/defense, 3) metabolism, and 4) gene regulation. We found heat shock protein-related and reactive oxygen species-related transcript derived fragments to be among the major parts of functional pathways induced by water deficit stress. Also, twelve novel transcripts were identified as both water deficit responsive and cotton specific. A subset of differentially expressed transcript derived fragments was verified using reverse transcription-polymerase chain reaction. Differential expression analysis also identified five pairs of duplicated transcript derived fragments in which four pairs responded differentially between each of their two homologues under water deficit stress.

CONCLUSIONS

In this study, we detected differentially expressed transcript derived fragments from water deficit stressed root and leaf tissues in tetraploid cotton and provided their gene ontology, functional/biological distribution, and possible roles of gene duplication. This discovery demonstrates complex mechanisms involved with polyploid cotton's transcriptome response to naturally occurring field water deficit stress. The genes identified in this study will provide candidate targets to manipulate the water use characteristics of cotton at the molecular level.

The role of Kupffer cells in rat liver regeneration revealed by cell-specific microarray analysis

Liver regeneration after partial hepatectomy is a process with various types of cells involved. The role of Kupffer cells (KCs) in liver regeneration is still controversial. In this study we isolated KCs from regenerating liver and conducted cell-specific microarray analysis. The results demonstrated that the controversial role of KCs in liver regeneration could be explained with the expression patterns of TGF-α, IL-6, TNF, and possibly IL-18 during liver regeneration. IL-18 may play an important role in negative regulation of liver regeneration. The functional profiles of gene expression in KCs also indicated that KC signaling might play a negative role in cell proliferation: signaling genes were down regulated before cell division. Immune response genes in KCs were also down regulated during liver regeneration, demonstrating similar expression profiles to that of hepatocytes. The expression patterns of key genes in these functional categories were consistent with the temporal functional profiles.

Analysis of the role of the integrin signaling pathway in hepatocytes during rat liver regeneration

To explore the role of the integrin signaling pathway in hepatocytes during rat liver regeneration, the integrin signaling pathway-related gene expression profile in hepatocytes of regenerative liver was detected using Rat Genome 230 2.0 array. The chip data showed that 265 genes of the integrin signaling pathway were included by Rat Genome 230 2.0 array and 132 genes showed significant expression changes in hepatocytes of regenerative liver. The numbers of up-, down- and up/down-regulated genes were 110, 15 and 7 respectively. In addition, bioinformatics and systems biology methods were used to analyze the role of the integrin signaling pathway in hepatocytes. The analysis of gene synergy value indicated that paths 1, 8, 12, and 15 promoted hepatocyte proliferation at the priming phase of liver regeneration; paths 1, 3, 8, and 12-15 enhanced hepatocyte proliferation at the progressing phase; paths 11 and 14 promoted hepatocyte proliferation, while paths 12 and 13 reduced hepatocyte proliferation at the terminal phase. Additionally, the other 8 paths (2, 4, 5-7, 9-10, and 16) were not found to be related to liver regeneration. In conclusion, 132 genes and 8 cascades of the integrin signaling pathway participated in regulating hepatocyte proliferation during rat liver regeneration.

High-throughput analysis of tumor necrosis factor signaling pathways in eight cell types during rat hepatic regeneration

This study aims to clarify the relevance of tumor necrosis factor (TNFs) signaling pathways and liver regeneration (LR) at the cellular level. Eight liver cell types were isolated using Percoll density gradient centrifugation and immunomagnetic beads methods. Expressions of TNF signaling pathway-involved genes in each cell type after 2/3 hepatectomy (PH) were detected using gene chip. Results show the following: gene TNFα was upregulated in most cell types, especially in Kupffer cells (KC); TNFβ expression was insignificantly changed in eight liver cell types; the majority of genes involved in four TNFα signaling pathways showed increased expression during LR in hepatocytes (HC); TNFα-induced NFκB pathway-involved genes were upregulated preferentially between 2 and 24 h during LR in biliary epithelial cells (BECs); and TNFα-induced apoptotic pathway genes were downregulated preferentially at progressing phase of LR in dendritic cells (DCs). Referring to the above results, TNFα-mediated signaling pathways, in contrast to TNFβ, play the more proactive role in LR, and four TNFα-mediated signaling pathways seem helpful to regulate biological events in HC; BEC proliferation was partly controlled by TNFα-mediated NFκB pathway; and the impaired TNFα-mediated apoptotic pathway in DCs might contribute to the restoration of DC mass after PH. Briefly, the comparative analysis of genomewide expression profiles of TNF signaling pathways between different cell types is helpful in understanding the implication of TNF signaling in LR at the cellular level.

Localization of high level of sequence conservation and divergence regions in cotton

In a previous study, we observed that the variations in chromosome size are due to uneven expansion and contraction by comparing the structures and sizes of a pair of homoeologous high-resolution cytogenetic maps of chromosomes 12A and 12D in tetraploid cotton. To reveal the variation at the sequence level, in the present paper, we sequenced two pairs of homoeologous bacterial artificial chromosomes derived from high- to low-variable genomic regions. Comparisons of their sequence variations confirmed that the highly conserved and divergent sequences existed in the distal and pericentric regions, e.g., high- and low-variable genome size regions in these two pairs of cotton homoeologous chromosomes. Sequence analysis also confirmed that the differential accumulation of Gossypium retrotransposable gypsy-like element (Gorge3) accounted for the main contributions for the size difference between the pericentric regions. By fluorescence in situ hybridization analysis, we found that Gorge3 has a bias distribution in the A(T)/A proximal regions and is associated with the heterochromatin along the chromosomes in the entire Gossypium genome. These results indicate that, between A(T)/A and D(T)/D genomes, the distal and pericentric regions usually possess high level of sequence conservation and divergence, respectively, in cotton.

Inheritance of long staple fiber quality traits of Gossypium barbadense in G. hirsutum background using CSILs

Gossypium hirsutum is a high yield cotton species that exhibits only moderate performance in fiber qualities. A promising but challenging approach to improving its phenotypes is interspecific introgression, the transfer of valuable traits or genes from the germplasm of another species such as G. barbadense, an important cultivated extra long staple cotton species. One set of chromosome segment introgression lines (CSILs) was developed, where TM-1, the genetic standard in G. hirsutum, was used as the recipient parent and the long staple cotton G. barbadense Hai7124 was used as the donor parent by molecular marker-assisted selection (MAS) in BC(5)S(1–4) and BC(4)S(1–3) generations. After four rounds of MAS, the CSIL population was comprised of 174 lines containing 298 introgressed segments, of which 86 (49.4%) lines had single introgressed segments. The total introgressed segment length covered 2,948.7 cM with an average length of 16.7 cM and represented 83.3% of tetraploid cotton genome. The CSILs were highly varied in major fiber qualities. By integrated analysis of data collected in four environments, a total of 43 additive quantitative trait loci (QTL) and six epistatic QTL associated with fiber qualities were detected by QTL IciMapping 3.0 and multi-QTL joint analysis. Six stable QTL were detected in various environments. The CSILs developed and the analyses presented here will enhance the understanding of the genetics of fiber qualities in long staple G. barbadense and facilitate further molecular breeding to improve fiber quality in Upland cotton.

Overexpression of GhSusA1 increases plant biomass and improves cotton fiber yield and quality

Cotton (Gossypium spp.) is an important economic crop and the largest source of textile fiber in the world. However, to date, only a few genes have been identified that exhibit critical roles in fiber development, and few has shown positive effects on fiber yield and quality in transgenic cotton. Here, we report the characterization of a novel sucrose synthase (SusA1) gene from a superior quality fiber germplasm line 7235 in Gossypium hirsutum. By association analysis, GhSusA1 was highly correlated with fiber qualities in (7235× TM-1) recombinant inbred lines based on polymorphism of GhSusA1 between 7235 and TM-1. Subsequently, based on an interspecific population of 141 BC₁ individuals generated from the cross between TM-1 and Gossypium barbadense line, Hai7124, we further mapped GhSusA1 genes on homeologous chromosomes A8 (chro.8) and D8 (chro.24). Suppression of GhSusA1 in transgenic cotton reduced fiber quality and decreased the boll size and seed weight. Importantly, overexpression of this gene increased fiber length and strength, with the latter indicated by the enhanced thickening of cell wall during secondary wall formation stage. Moreover, increasing GhSusA1 transcript abundance in vegetative tissues led to elevated seedling biomass. Together, these findings identified GhSusA1 as a key regulator of sink strength in cotton, which is tightly associated with productivity, and hence a promising candidate gene that can be developed to increase cotton fiber yield and quality.

Evolutionary dynamics of the Pgk1 gene in the polyploid genus Kengyilia (Triticeae: Poaceae) and its diploid relatives

The level and pattern of nucleotide variation in duplicate gene provide important information on the evolutionary history of polyploids and divergent process between homoeologous loci within lineages. Kengyilia is a group of allohexaploid species with the StYP genomic constitutions in the wheat tribe. To investigate the evolutionary dynamics of the Pgk1 gene in Kengyilia and its diploid relatives, three copies of Pgk1 homoeologues were isolated from all sampled hexaploid Kengyilia species and analyzed with the Pgk1 sequences from 47 diploid taxa representing 18 basic genomes in Triticeae. Sequence diversity patterns and genealogical analysis suggested that (1) Kengyilia species from the Central Asia and the Qinghai-Tibetan plateau have independent origins with geographically differentiated P genome donors and diverged levels of nucleotide diversity at Pgk1 locus; (2) a relatively long-time sweep event has allowed the Pgk1 gene within Agropyron to adapt to cold climate triggered by the recent uplifts of the Qinghai-Tibetan Plateau; (3) sweep event and population expansion might result in the difference in the d_N/d_S value of the Pgk1 gene in allopatric Agropyron populations, and this difference may be genetically transmitted to Kengyilia lineages via independent polyploidization events; (4) an 83 bp MITE element insertion has shaped the Pgk1 loci in the P genome lineage with different geographical regions; (5) the St and P genomes in Kengyilia were donated by Pseudoroegneria and Agropyron, respectively, and the Y genome is closely related to the Xp genome of Peridictyon sanctum. The interplay of evolutionary forces involving diverged natural selection, population expansion, and transposable events in geographically differentiated P genome donors could attribute to geographical differentiation of Kengyilia species via independent origins.