Comparative genomics of Gossypium and Arabidopsis: unraveling the consequences of both ancient and recent polyploidy

GhIMP10D, an inositol monophosphates family gene, enhances ascorbic acid and antioxidant enzyme activities to confer alkaline tolerance in Gossypium hirsutum L

BACKGROUND

Inositol monophosphates (IMP) are key enzymes in the ascorbic acid (AsA) synthesis pathways, which play vital roles in regulating plant growth and development and stresses tolerance. To date, no comprehensive analysis of the expression profile of IMP genes and their functions under abiotic stress in cotton has been reported.

RESULTS

In this study, the genetic characteristics, phylogenetic evolution, cis-acting elements and expression patterns of IMP gene family in cotton were systematically analyzed. A total of 28, 27, 13 and 13 IMP genes were identified in Gossypium hirsutum (G. hirsutum), Gossypium barbadense (G. barbadense), Gossypium arboreum (G. arboreum), and Gossypium raimondii (G. raimondii), respectively. Phylogenetic analysis showed that IMP family genes could cluster into 3 clades. Structure analysis of genes showed that GhIMP genes from the same subgroup had similar genetic structure and exon number. And most GhIMP family members contained hormone-related elements (abscisic acid response element, MeJA response element, gibberellin response element) and stress-related elements (low temperature response element, defense and stress response element, wound response element). After exogenous application of abscisic acid (ABA), some GhIMP genes containing ABA response elements positively responded to alkaline stress, indicating that ABA response elements played an important role in response to alkaline stress. qRT-PCR showed that most of GhIMP genes responded positively to alkaline stress, and GhIMP10D significantly upregulated under alkaline stress, with the highest up-regulated expression level. Virus-induced gene silencing (VIGS) experiment showed that compared with 156 plants, MDA content of pYL156:GhIMP10D plants increased significantly, while POD, SOD, chlorophyII and AsA content decreased significantly.

CONCLUSIONS

This study provides a thorough overview of the IMP gene family and presents a new perspective on the evolution of this gene family. In particular, some IMP genes may be involved in alkaline stress tolerance regulation, and GhIMP10D showed high expression levels in leaves, stems and roots under alkaline stress, and preliminary functional verification of GhIMP10D gene suggested that it may regulate tolerance to alkaline stress by regulating the activity of antioxidant enzymes and the content of AsA. This study contributes to the subsequent broader discussion of the structure and alkaline resistance of IMP genes in cotton.

Transcriptome Meta-Analysis Associated Targeting Hub Genes and Pathways of Drought and Salt Stress Responses in Cotton (Gossypium hirsutum): A Network Biology Approach

Abiotic stress tolerance is an intricate feature controlled through several genes and networks in the plant system. In abiotic stress, salt, and drought are well known to limit cotton productivity. Transcriptomics meta-analysis has arisen as a robust method to unravel the stress-responsive molecular network in crops. In order to understand drought and salt stress tolerance mechanisms, a meta-analysis of transcriptome studies is crucial. To confront these issues, here, we have given details of genes and networks associated with significant differential expression in response to salt and drought stress. The key regulatory hub genes of drought and salt stress conditions have notable associations with functional drought and salt stress-responsive (DSSR) genes. In the network study, nodulation signaling pathways 2 (NSP2), Dehydration-responsive element1 D (DRE1D), ethylene response factor (ERF61), cycling DOF factor 1 (CDF1), and tubby like protein 3 (TLP3) genes in drought and tubby like protein 1 (TLP1), thaumatin-like proteins (TLP), ethylene-responsive transcription factor ERF109 (EF109), ETS-Related transcription Factor (ELF4), and Arabidopsis thaliana homeodomain leucine-zipper gene (ATHB7) genes in salt showed the significant putative functions and pathways related to providing tolerance against drought and salt stress conditions along with the significant expression values. These outcomes provide potential candidate genes for further in-depth functional studies in cotton, which could be useful for the selection of an improved genotype of Gossypium hirsutum against drought and salt stress conditions.

A genome-wide comparative analysis of bZIP transcription factors in G. arboreum and G. raimondii (Diploid ancestors of present-day cotton)

Basic leucine zipper motif (bZIP) transcription factors (TFs) are involved in plant growth regulation, development, and environmental stress responses. These genes have been well characterized in model plants. In current study, a genome-wide analysis of bZIP genes was performed in Gossypium raimondii and Gossypium arboreum taking Arabidopsis thaliana as a reference genome. In total, 85 members of G. raimondii and 87 members of G. arboreum were identified and designated as GrbZIPs and GabZIPs respectively. Phylogenetic analysis clustered bZIP genes into 11 subgroups (A, B, C, D, F, G, H, I, S and X). Gene structure analysis to find the intro-exon structures revealed 1-14 exons in both species. The maximum number of introns were present in subgroup G and D while genes in subgroup S were intron-less except GrbZIP78, which is a unique characteristic as compared to other groups. Results of motif analysis predicted that all three species share a common bZIP motif. A detailed comparison of bZIPs gene distribution on chromosomes has shown a diverse arrangement of genes in both cotton species. Moreover, the functional similarity with orthologs was also predicted. The findings of this study revealed close similarity in gene structure of both cotton species and diversity in gene distribution on chromosomes. This study supports the divergence of both species from the common ancestor and later diversity in gene distribution on chromosomes due to evolutionary changes. Additionally, this work will facilitate the functional characterization of bZIP genes in cotton. Outcomes of this study represent foundation research on the bZIP TFs family in cotton and as a reference for other crops.

An EMS-induced mutation in a tetratricopeptide repeat-like superfamily protein gene (Ghir_A12G008870) on chromosome A12 is responsible for the liy short fiber phenotype in cotton

The EMS-induced threonine/isoleucine substitution in a tetratricopeptide repeat-like superfamily protein encoded by gene Ghir_A12G008870 is responsible for the Ligon-lintless-y (li_y) short fiber phenotype in cotton. A short fiber mutant Ligon-lintless-y was created through treating the seeds of the cotton line MD15 with ethyl methanesulfonate. Genetic analysis indicated that the short fiber phenotype is controlled by a single recessive locus designated li_y. From F₂ populations derived from crosses between the mutant and its wild type (WT), we selected 132 short fiber progeny (li_y/li_y) and made two DNA bulks. We sequenced these DNA bulks along with the two parents of the population. The li_y locus was located on chromosome A12. Using multiple F₂ populations and F₃ progeny plants, we mapped the li_y locus within a genomic region of 1.18 Mb. In this region, there is only one gene, i.e., Ghir_A12G008870 encoding a tetratricopeptide repeat-like superfamily protein that has a non-synonymous mutation between the li_y mutant and its WT. Analysis of a SNP marker representing this gene in the F₂ and F₃ progeny plants demonstrated its complete linkage with the li_y short fiber phenotype. We further analyzed this SNP marker in a panel of 384 cotton varieties. The mutant allele is absent in all varieties analyzed. RNAseq and RT-qPCR analysis of the gene Ghir_A12G008870 during fiber development showed a significant expression difference between the li_y mutant and its WT in developing fiber cells beginning at 12 days post-anthesis. Virus-induced gene silencing of the gene Ghir_A12G008870 significantly reduced the fiber length of the WT cotton line MD15. Taken together, our results suggest that the gene Ghir_A12G008870 is involved in the cotton fiber cell elongation process and is a promising candidate gene responsible for the li_y short fiber phenotype.

Genome-wide comparative and evolutionary analysis of Calmodulin-binding Transcription Activator (CAMTA) family in Gossypium species

The CAMTA gene family is crucial in managing both biotic and abiotic stresses in plants. Our comprehensive analysis of this gene family in cotton resulted in the identification of 6, 7 and 9 CAMTAs in three sequenced cotton species, i.e., Gossypium arboreum, Gossypium raimondii, and Gossypium hirsutum, respectively. All cotton CAMTAs were localized in the nucleus and possessed calmodulin-binding domain (CaMBD) as identified computationally. Phylogenetically four significant groups of cotton CAMTAs were identified out of which, Group II CAMTAs experienced higher evolutionary pressure, leading to a faster evolution in diploid cotton. The expansion of cotton CAMTAs in the genome was mainly due to segmental duplication. Purifying selection played a significant role in the evolution of cotton CAMTAs. Expression profiles of GhCAMTAs revealed that GhCAMTA2A.2 and GhCAMTA7A express profoundly in different stages of cotton fiber development. Positive correlation between expression of these two CAMTAs and fiber strength confirmed their functional relevance in fiber development. The promoter region of co-expressing genes network of GhCAMTA2A.2 and GhCAMTA7A showed a higher frequency of occurrence of CAMTA binding motifs. Our present study thus contributes to broad probing into the structure and probable function of CAMTA genes in Gossypium species.

Microdissection of the Ah01 chromosome in upland cotton and microcloning of resistance gene anologs from the single chromosome

BACKGROUND

Chromosome microdissection is one of the most important techniques in molecular cytogenetic research. Cotton (Gossypium Linnaeus, 1753) is the main natural fiber crop in the world. The resistance gene analog (RGA) cloning after its single chromosome microdissection can greatly promote cotton genome research and breeding.

RESULTS

Using the linker adaptor PCR (LA-PCR) with the primers of rice disease-resistance homologues, three nucleotide sequences PS016 (KU051681), PS054 (KU051682), and PS157 (KU051680) were obtained from the chromosome A_h01 of upland cotton (cv. TM-1). The Blast results showed that the three sequences are the nucleotide binding site-leucine rich repeat (NBS-LRR) type RGAs. Clustering results indicated that they are homologous to these published RGAs. Thus, the three RGAs can definitely be confirmed as NBS-LRR class of RGAs in upland cotton.

CONCLUSIONS

Using single chromosome microdissection technique, DNA libraries containing cotton RGAs were obtained. This technique can promote cotton gene cloning, marker development and even the improvement of cotton genome research and breeding.

Toward Coalescing Gene Expression and Function with QTLs of Water-Deficit Stress in Cotton

Cotton exhibits moderately high vegetative tolerance to water-deficit stress but lint production is restricted by the available rainfed and irrigation capacity. We have described the impact of water-deficit stress on the genetic and metabolic control of fiber quality and production. Here we examine the association of tentative consensus sequences (TCs) derived from various cotton tissues under irrigated and water-limited conditions with stress-responsive QTLs. Three thousand sixteen mapped sequence-tagged-sites were used as anchored targets to examine sequence homology with 15,784 TCs to test the hypothesis that putative stress-responsive genes will map within QTLs associated with stress-related phenotypic variation more frequently than with other genomic regions not associated with these QTLs. Approximately 1,906 of 15,784 TCs were mapped to the consensus map. About 35% of the annotated TCs that mapped within QTL regions were genes involved in an abiotic stress response. By comparison, only 14.5% of the annotated TCs mapped outside these QTLs were classified as abiotic stress genes. A simple binomial probability calculation of this degree of bias being observed if QTL and non-QTL regions are equally likely to contain stress genes was P (x ≥ 85) = 7.99  × 10(-15). These results suggest that the QTL regions have a higher propensity to contain stress genes.

Exploiting a Reference Genome in Terms of Duplications: The Network of Paralogs and Single Copy Genes in Arabidopsis thaliana

Arabidopsis thaliana became the model organism for plant studies because of its small diploid genome, rapid lifecycle and short adult size. Its genome was the first among plants to be sequenced, becoming the reference in plant genomics. However, the Arabidopsis genome is characterized by an inherently complex organization, since it has undergone ancient whole genome duplications, followed by gene reduction, diploidization events and extended rearrangements, which relocated and split up the retained portions. These events, together with probable chromosome reductions, dramatically increased the genome complexity, limiting its role as a reference. The identification of paralogs and single copy genes within a highly duplicated genome is a prerequisite to understand its organization and evolution and to improve its exploitation in comparative genomics. This is still controversial, even in the widely studied Arabidopsis genome. This is also due to the lack of a reference bioinformatics pipeline that could exhaustively identify paralogs and singleton genes. We describe here a complete computational strategy to detect both duplicated and single copy genes in a genome, discussing all the methodological issues that may strongly affect the results, their quality and their reliability. This approach was used to analyze the organization of Arabidopsis nuclear protein coding genes, and besides classifying computationally defined paralogs into networks and single copy genes into different classes, it unraveled further intriguing aspects concerning the genome annotation and the gene relationships in this reference plant species. Since our results may be useful for comparative genomics and genome functional analyses, we organized a dedicated web interface to make them accessible to the scientific community

BAC-pool sequencing and analysis of large segments of A12 and D12 homoeologous chromosomes in upland cotton

Although new and emerging next-generation sequencing (NGS) technologies have reduced sequencing costs significantly, much work remains to implement them for de novo sequencing of complex and highly repetitive genomes such as the tetraploid genome of Upland cotton (Gossypium hirsutum L.). Herein we report the results from implementing a novel, hybrid Sanger/454-based BAC-pool sequencing strategy using minimum tiling path (MTP) BACs from Ctg-3301 and Ctg-465, two large genomic segments in A12 and D12 homoeologous chromosomes (Ctg). To enable generation of longer contig sequences in assembly, we implemented a hybrid assembly method to process ~35x data from 454 technology and 2.8-3x data from Sanger method. Hybrid assemblies offered higher sequence coverage and better sequence assemblies. Homology studies revealed the presence of retrotransposon regions like Copia and Gypsy elements in these contigs and also helped in identifying new genomic SSRs. Unigenes were anchored to the sequences in Ctg-3301 and Ctg-465 to support the physical map. Gene density, gene structure and protein sequence information derived from protein prediction programs were used to obtain the functional annotation of these genes. Comparative analysis of both contigs with Arabidopsis genome exhibited synteny and microcollinearity with a conserved gene order in both genomes. This study provides insight about use of MTP-based BAC-pool sequencing approach for sequencing complex polyploid genomes with limited constraints in generating better sequence assemblies to build reference scaffold sequences. Combining the utilities of MTP-based BAC-pool sequencing with current longer and short read NGS technologies in multiplexed format would provide a new direction to cost-effectively and precisely sequence complex plant genomes.

A whole-genome DNA marker map for cotton based on the D-genome sequence of Gossypium raimondii L

We constructed a very-high-density, whole-genome marker map (WGMM) for cotton by using 18,597 DNA markers corresponding to 48,958 loci that were aligned to both a consensus genetic map and a reference genome sequence. The WGMM has a density of one locus per 15.6 kb, or an average of 1.3 loci per gene. The WGMM was anchored by the use of colinear markers to a detailed genetic map, providing recombinational information. Mapped markers occurred at relatively greater physical densities in distal chromosomal regions and lower physical densities in the central regions, with all 1 Mb bins having at least nine markers. Hotspots for quantitative trait loci and resistance gene analog clusters were aligned to the map and DNA markers identified for targeting of these regions of high practical importance. Based on the cotton D genome reference sequence, the locations of chromosome structural rearrangements plotted on the map facilitate its translation to other Gossypium genome types. The WGMM is a versatile genetic map for marker assisted breeding, fine mapping and cloning of genes and quantitative trait loci, developing new genetic markers and maps, genome-wide association mapping, and genome evolution studies.

Repeated polyploidization of Gossypium genomes and the evolution of spinnable cotton fibres

Polyploidy often confers emergent properties, such as the higher fibre productivity and quality of tetraploid cottons than diploid cottons bred for the same environments. Here we show that an abrupt five- to sixfold ploidy increase approximately 60 million years (Myr) ago, and allopolyploidy reuniting divergent Gossypium genomes approximately 1-2 Myr ago, conferred about 30-36-fold duplication of ancestral angiosperm (flowering plant) genes in elite cottons (Gossypium hirsutum and Gossypium barbadense), genetic complexity equalled only by Brassica among sequenced angiosperms. Nascent fibre evolution, before allopolyploidy, is elucidated by comparison of spinnable-fibred Gossypium herbaceum A and non-spinnable Gossypium longicalyx F genomes to one another and the outgroup D genome of non-spinnable Gossypium raimondii. The sequence of a G. hirsutum A(t)D(t) (in which 't' indicates tetraploid) cultivar reveals many non-reciprocal DNA exchanges between subgenomes that may have contributed to phenotypic innovation and/or other emergent properties such as ecological adaptation by polyploids. Most DNA-level novelty in G. hirsutum recombines alleles from the D-genome progenitor native to its New World habitat and the Old World A-genome progenitor in which spinnable fibre evolved. Coordinated expression changes in proximal groups of functionally distinct genes, including a nuclear mitochondrial DNA block, may account for clusters of cotton-fibre quantitative trait loci affecting diverse traits. Opportunities abound for dissecting emergent properties of other polyploids, particularly angiosperms, by comparison to diploid progenitors and outgroups.

Analyses of the sucrose synthase gene family in cotton: structure, phylogeny and expression patterns

BACKGROUND

In plants, sucrose synthase (Sus) is widely considered as a key enzyme involved in sucrose metabolism. Several paralogous genes encoding different isozymes of Sus have been identified and characterized in multiple plant genomes, while limited information of Sus genes is available to date for cotton.

RESULTS

Here, we report the molecular cloning, structural organization, phylogenetic evolution and expression profiles of seven Sus genes (GaSus1 to 7) identified from diploid fiber cotton (Gossypium arboreum). Comparisons between cDNA and genomic sequences revealed that the cotton GaSus genes were interrupted by multiple introns. Comparative screening of introns in homologous genes demonstrated that the number and position of Sus introns are highly conserved among Sus genes in cotton and other more distantly related plant species. Phylogenetic analysis showed that GaSus1, GaSus2, GaSus3, GaSus4 and GaSus5 could be clustered together into a dicot Sus group, while GaSus6 and GaSus7 were separated evenly into other two groups, with members from both dicot and monocot species. Expression profiles analyses of the seven Sus genes indicated that except GaSus2, of which the transcripts was undetectable in all tissues examined, and GaSus7, which was only expressed in stem and petal, the other five paralogues were differentially expressed in a wide ranges of tissues, and showed development-dependent expression profiles in cotton fiber cells.

CONCLUSIONS

This is a comprehensive study of the Sus gene family in cotton plant. The results presented in this work provide new insights into the evolutionary conservation and sub-functional divergence of the cotton Sus gene family in response to cotton fiber growth and development.

The cytonuclear dimension of allopolyploid evolution: an example from cotton using rubisco

During allopolyploid speciation, two divergent nuclear genomes merge, yet only one (usually the maternal) of the two sets of progenitor organellar genomes is maintained. Rubisco (1,5-bisphosphate carboxylase/oxygenase) is composed of nuclear-encoded small subunits (SSUs) and plastome-encoded large subunits (LSUs), providing an ideal system to explore the evolutionary process of cytonuclear accommodation. Here, we take initial steps in this direction, using Gossypium allopolyploids as our model. SSU copies from divergent (5-10 My) progenitor diploids ("A" and "D" genomes) were combined at the time of polyploid formation 1-2 Ma, with the LSU encoded by the maternal A-genome parent. LSU genes from A- and D-genome diploids and AD-genome allopolyploids were sequenced, revealing several nonsynonymous substitutions and suggesting the possibility of differential selection on the nuclear-encoded rbcS partner following allopolyploid formation. Sequence data for the rbcS gene family revealed nonreciprocal homoeologous recombination between A- and D-rbcS homoeologs in all polyploid species but not in a synthetic intergenomic F1 hybrid, demonstrating "gene conversion" during allopolyploid evolution. All progenitor rbcS genes are retained and expressed in the five extant allopolyploid species, but analysis of the leaf transcriptome showed that A-homoeologs are preferentially expressed in both the allopolyploid and hybrid, consistent with the maternal origin of rbcL. Although rbcS genes from both progenitor genomes are expressed, some appear to have experienced mutations that may represent cytonuclear coevolution.

Ancestral synteny shared between distantly-related plant species from the asterid (Coffea canephora and Solanum Sp.) and rosid (Vitis vinifera) clades

BACKGROUND

Coffee trees (Rubiaceae) and tomato (Solanaceae) belong to the Asterid clade, while grapevine (Vitaceae) belongs to the Rosid clade. Coffee and tomato separated from grapevine 125 million years ago, while coffee and tomato diverged 83-89 million years ago. These long periods of divergent evolution should have permitted the genomes to reorganize significantly. So far, very few comparative mappings have been performed between very distantly related species belonging to different clades. We report the first multiple comparison between species from Asterid and Rosid clades, to examine both macro-and microsynteny relationships.

RESULTS

Thanks to a set of 867 COSII markers, macrosynteny was detected between coffee, tomato and grapevine. While coffee and tomato genomes share 318 orthologous markers and 27 conserved syntenic segments (CSSs), coffee and grapevine also share a similar number of syntenic markers and CSSs: 299 and 29 respectively. Despite large genome macrostructure reorganization, several large chromosome segments showed outstanding macrosynteny shedding new insights into chromosome evolution between Asterids and Rosids. We also analyzed a sequence of 174 kb containing the ovate gene, conserved in a syntenic block between coffee, tomato and grapevine that showed a high-level of microstructure conservation. A higher level of conservation was observed between coffee and grapevine, both woody and long life-cycle plants, than between coffee and tomato. Out of 16 coffee genes of this syntenic segment, 7 and 14 showed complete synteny between coffee and tomato or grapevine, respectively.

CONCLUSIONS

These results show that significant conservation is found between distantly related species from the Asterid (Coffea canephora and Solanum sp.) and Rosid (Vitis vinifera) clades, at the genome macrostructure and microstructure levels. At the ovate locus, conservation did not decline in relation to increasing phylogenetic distance, suggesting that the time factor alone does not explain divergences. Our results are considerably useful for syntenic studies between supposedly remote species for the isolation of important genes for agronomy.

Use of diversity arrays technology markers for integration into a cotton reference map and anchoring to a recombinant inbred line map

A diversity array technology (DArT) marker platform was developed for the cotton genome, to evaluate the use of DArT markers compared with AFLP markers in mapping and transferability across the mapping populations. We used a reference genetic map of tetraploid Gossypium L. that already contained ~5000 loci, which coalesced into 26 chromosomes, to anchor newly developed DArT and AFLP markers with the aim of further improving utility and map resolution. Our results indicated that the percentage of polymorphic DArT markers that could be genetically mapped (78.15%) was much higher than that of AFLP markers (22.28%). Sequence analysis of DArT markers indicated that a majority matched known expressed sequence tag (EST) sequences from tetraploid and diploid Gossypium species. A total of 794 Arabidopsis genes were homologous with various DArT marker sequences. Chromosomes 5(A), 7(A), 19(D), 23(D), and 24(D) had more Arabidopsis syntenic DArT markers than the other chromosomes. Anchoring DArT markers from the reference map to a recombinant inbred line (RIL) map indicated that DArT markers will speed the building of maps in de novo RIL populations.

Silencing GhNDR1 and GhMKK2 compromises cotton resistance to Verticillium wilt

Cotton is an important cash crop worldwide, and is a significant source of fiber, feed, foodstuff, oil and biofuel products. Considerable effort has been expended to increase sustainable yield and quality through molecular breeding and genetic engineering of new cotton cultivars. Given the recent availability of the whole-genome sequence of cotton, it is necessary to develop molecular tools and resources for large-scale analysis of gene functions at the genome-wide level. We have successfully developed an Agrobacterium-mediated virus-induced gene silencing (VIGS) assay in several cotton cultivars with various genetic backgrounds. The genes of interest were potently and readily silenced within 2 weeks after inoculation at the seedling stage. Importantly, we showed that silencing GhNDR1 and GhMKK2 compromised cotton resistance to the infection by Verticillium dahliae, a fungal pathogen causing Verticillium wilt. Furthermore, we developed a cotton protoplast system for transient gene expression to study gene functions by a gain-of-function approach. The viable protoplasts were isolated from green cotyledons, etiolated cotyledons and true leaves, and responded to a wide range of pathogen elicitors and phytohormones. Remarkably, cotton plants possess conserved, but also distinct, MAP kinase activation with Arabidopsis upon bacterial elicitor flagellin perception. Thus, using gene silencing assays, we have shown that GhNDR1 and GhMKK2 are required for Verticillium resistance in cotton, and have developed high throughput loss-of-function and gain-of-function assays for functional genomic studies in cotton.

Metabolite and mineral analyses of cotton near-isogenic lines introgressed with QTLs for productivity and drought-related traits

Quantitative trait loci (QTLs) for yield and drought-related traits were exchanged via marker-assisted selection between elite cultivars of two cotton species, Gossypium barbadense (GB) cv. F-177 and Gossypium hirsutum (GH) cv. Siv'on. Three of the resultant near-isogenic lines (NILs), each introgressed with a different QTL region, expressed an advantage in osmotic adjustment (OA) and other drought-related traits relative to their recipient parents. These NILs and the parental genotypes were field-grown under well-watered and water-limited conditions, and characterized for their metabolic and mineral compositions. Comparisons were then made between (1) GB and GH genotypes, (2) the contrasting water regimes and (3) each NIL and its recipient parent. Hierarchical clustering analysis clearly distinguished between GB and GH genotypes based on either metabolite or mineral composition. Comparisons between well-watered and water-limited conditions in each of the genotypes showed differing trends in the various solutes. The greater concentrations of potassium, magnesium and calcium under water stress, when compared with well-watered conditions, may have enhanced OA or osmoprotection. All NILs exhibited significantly modified solute composition relative to their recipient parents. In particular, increased levels of alanine, aspartic acid, citric acid, malic acid, glycerol, myoinositol, threonic acid, potassium, magnesium and calcium were found under drought conditions in one or more of the NILs relative to their recipient parents. The increased values of these solutes could contribute to the superior capacity of these NILs to cope with drought.

Comparative analysis of Gossypium and Vitis genomes indicates genome duplication specific to the Gossypium lineage

Genetic mapping studies have suggested that diploid cotton (Gossypium) might be an ancient polyploid. However, further evidence is lacking due to the complexity of the genome and the lack of sequence resources. Here, we used the grape (Vitis vinifera) genome as an out-group in two different approaches to further explore evidence regarding ancient genome duplication (WGD) event(s) in the diploid Gossypium lineage and its (their) effects: a genome-level alignment analysis and a local-level sequence component analysis. Both studies suggest that at least one round of genome duplication occurred in the Gossypium lineage. Also, gene densities in corresponding regions from Gossypium raimondii, V. vinifera, Arabidopsis thaliana and Carica papaya genomes are similar, despite the huge difference in their genome sizes and the different number of WGDs each genome has experienced. These observations fit the model that differences in plant genome sizes are largely explained by transposon insertions into heterochromatic regions.

A draft physical map of a D-genome cotton species (Gossypium raimondii)

BACKGROUND

Genetically anchored physical maps of large eukaryotic genomes have proven useful both for their intrinsic merit and as an adjunct to genome sequencing. Cultivated tetraploid cottons, Gossypium hirsutum and G. barbadense, share a common ancestor formed by a merger of the A and D genomes about 1-2 million years ago. Toward the long-term goal of characterizing the spectrum of diversity among cotton genomes, the worldwide cotton community has prioritized the D genome progenitor Gossypium raimondii for complete sequencing.

RESULTS

A whole genome physical map of G. raimondii, the putative D genome ancestral species of tetraploid cottons was assembled, integrating genetically-anchored overgo hybridization probes, agarose based fingerprints and 'high information content fingerprinting' (HICF). A total of 13,662 BAC-end sequences and 2,828 DNA probes were used in genetically anchoring 1585 contigs to a cotton consensus genetic map, and 370 and 438 contigs, respectively to Arabidopsis thaliana (AT) and Vitis vinifera (VV) whole genome sequences.

CONCLUSION

Several lines of evidence suggest that the G. raimondii genome is comprised of two qualitatively different components. Much of the gene rich component is aligned to the Arabidopsis and Vitis vinifera genomes and shows promise for utilizing translational genomic approaches in understanding this important genome and its resident genes. The integrated genetic-physical map is of value both in assembling and validating a planned reference sequence.

The construction of a tetraploid cotton genome wide comprehensive reference map

Integration of multiple genomic maps provides a higher density of markers and greater genome coverage, which not only facilitates the identification and positioning of QTLs and candidate genes, but it also provides a basic structure for the genome sequence assembly. However, the diversity in markers and populations used in individual mapping studies limits the ability to fully integrate the available data. By concentrating on marker orders rather than marker distances, published map data could be used to produce a comprehensive reference map (CRM) that includes a majority of known markers with optimally estimated order of those markers across the genome. In this study, a tetraploid cotton genome-wide CRM was constructed from 28 public cotton genetic maps. The initial CRM contained 7,424 markers and represented over 93% of the combined mapping information from the 28 individual maps. The current output is stored and displayed through CottonDB (http://www.cottondb.org), the public cotton genome database.

Gene copy number evolution during tetraploid cotton radiation

After polyploid formation, retention or loss of duplicated genes is not random. Genes with some functional domains are convergently restored to 'singleton' state after many independent genome duplications, and have been referred to as 'duplication-resistant' (DR) genes. To further explore the timeframe for their restoration to the singleton state, 27 cotton homologs of genes found to be 'DR' in Arabidopsis were selected based on diagnostic Pfam domains. Their copy numbers were studied using southern hybridization and sequence analysis in five tetraploid species and their ancestral A and D genome diploids. DR genes had significantly lower copy number than gene families hybridizing to randomly selected cotton ESTs. Three DR genes showed complete loss of D genome-derived homoeologs in some or all tetraploid species. Prior analysis has shown gene loss in polyploid cotton to be rare, and herein only one randomly selected gene showed loss of a homoeolog in only one of the five tetraploid species (Gossypium mustelinum). BAC sequencing confirmed two cases of gene loss in tetraploid cotton. Divergence among 5' sequences of DR genes amplified from G. arboreum, G. raimondii, and Gossypioides kirkii was correlated with gene copy number. These results show that genes containing Pfam domains associated with duplication resistance in Arabidopsis have also been preferentially restored to low copy number after a more recent polyploidization event in cotton. In tetraploid cotton, genes from the progenitor D genome seem to experience more gene copy number divergence than genes from the A genome. Together with D subgenome-biased alterations in gene expression, perhaps gene loss may contribute to the relatively larger portion of quantitative trait variation attributable to D than A subgenome chromosomes of tetraploid cotton.

Microsatellites in Brassica unigenes: relative abundance, marker design, and use in comparative physical mapping and genome analysis

Microsatellites present in the transcribed regions of the genome have the potential to reveal functional diversity. Unigene sequence databases are the sources of such genic microsatellites with unique flanking sequences and genomic locations even in complex polyploids. The present study was designed to assay the unigenes of Brassica napus and B. rapa for various microsatellite repeats, and to design markers and use them in comparative genome analysis and study of evolution. The average frequency of microsatellites in Brassica unigenes was one in every 7.25 kb of sequence, as compared with one in every 8.57 kb of sequence in Arabidopsis thaliana . Trinucleotide motifs coding for serine and the dinucleotide motif GA were most abundant. We designed 2374 and 347 unigene-based microsatellite (UGMS) markers including 541 and 58 class I types in B. napus and B. rapa, respectively, and evaluated their use across diverse species and genera. Most of these markers (93.3%) gave successful amplification of target microsatellite motifs, which was confirmed by sequencing. Interspecific polymorphism between B. napus and B. rapa detected in silico for the UGMS markers was 4.16 times higher in 5′ untranslated regions than in coding sequences. Physical anchoring of Brassica UGMS markers on the A. thaliana genome indicated their significance in studying the evolutionary history of A. thaliana genomic duplications in relation to speciation. Comparative physical mapping identified 85% of Brassica unigenes as single copy and gave clues for the presence of conserved primordial gene order. Complex chromosomal rearrangements such as inversions, tandem and segmental duplications, and insertions/deletions were evident between A. thaliana and B. rapa genomes. The results obtained have encouraging implications for the use of UGMS markers in comparative genome analysis and for understanding evolutionary complexities in the family Brassicaceae.

Resistance to Thielaviopsis basicola in the cultivated A genome cotton

Black root rot (BRR), incited by the soilborne pathogen Thielaviopsis basicola has the potential to cause significant economic loss in cotton (Gossypium spp.) production. Cultivated tetraploids of cotton (G. hirsutum and G. barbadense) are susceptible although resistant types have been identified in a possible tetraploid progenitor, G. herbaceum. Genetic mapping was used to detect the chromosomal locations of quantitative trait loci (QTL) that confer resistance to the BRR pathogen. A population of F(2) individuals (G. herbaceum x G. arboreum) and F(2:3) progeny families were examined. Phenotypic variation between resistant and susceptible reactions could be explained partly by three QTL. The BRR5.1, BRR9.1, and BRR13.1 QTL each explained 19.1, 10.3 and 8.5% of the total phenotypic variation, respectively. The combination of all three in a single genetic model explained 32.7% of the phenotypic variation. Comparative analysis was conducted on significant QTL regions to deduce the cotton-Arabidopsis synteny relationship and examine the correspondence between BRR QTL and Arabidopsis pathogen defense genes. Totally 20 Arabidopsis synteny segments corresponded within one of three BRR QTL regions. Each synteny segment contains many potential Arabidopsis candidate genes. A total of 624 Arabidopsis genes, including 22 pathogen defense and 36 stress response genes, could be placed within the syntenic regions corresponding to the BRR QTL. Fine mapping is needed to delineate each underlying BRR R-gene and possible Arabidopsis orthologs. Research and breeding activities to examine each QTL and underlying genes in Upland cotton (G. hirsutum) are ongoing.

Geminivirus-mediated gene silencing from Cotton leaf crumple virus is enhanced by low temperature in cotton

A silencing vector for cotton (Gossypium hirsutum) was developed from the geminivirus Cotton leaf crumple virus (CLCrV). The CLCrV coat protein gene was replaced by up to 500 bp of DNA homologous to one of two endogenous genes, the magnesium chelatase subunit I gene (ChlI) or the phytoene desaturase gene (PDS). Cotyledons of cotton cultivar 'Deltapine 5415' bombarded with the modified viral vectors manifested chlorosis due to silencing of either ChlI or PDS in approximately 70% of inoculated plants after 2 to 3 weeks. Use of the green fluorescence protein gene showed that replication of viral DNA was restricted to vascular tissue and that the viral vector could transmit to leaves, roots, and the ovule integument from which fibers originate. Temperature had profound effects on vector DNA accumulation and the spread of endogenous gene silencing. Consistent with reports that silencing against viruses increases at higher temperatures, plants grown at a 30 degrees C/26 degrees C day/night cycle had a greater than 10-fold reduction in viral DNA accumulation compared to plants grown at 22 degrees C/18 degrees C. However, endogenous gene silencing decreased at 30 degrees C/26 degrees C. There was an approximately 7 d delay in the onset of gene silencing at 22 degrees C/18 degrees C, but silencing was extensive and persisted throughout the life of the plant. The extent of silencing in new growth could be increased or decreased by changing temperature regimes at various times following the onset of silencing. Our experiments establish the use of the CLCrV silencing vector to study gene function in cotton and show that temperature can have a major impact on the extent of geminivirus-induced gene silencing.

Synteny conservation between two distantly-related Rosaceae genomes: Prunus (the stone fruits) and Fragaria (the strawberry)

BACKGROUND

The Rosaceae encompass a large number of economically-important diploid and polyploid fruit and ornamental species in many different genera. The basic chromosome numbers of these genera are x = 7, 8 and 9 and all have compact and relatively similar genome sizes. Comparative mapping between distantly-related genera has been performed to a limited extent in the Rosaceae including a comparison between Malus (subfamily Maloideae) and Prunus (subfamily Prunoideae); however no data has been published to date comparing Malus or Prunus to a member of the subfamily Rosoideae. In this paper we compare the genome of Fragaria, a member of the Rosoideae, to Prunus, a member of the Prunoideae.

RESULTS

The diploid genomes of Prunus (2n = 2x = 16) and Fragaria (2n = 2x = 14) were compared through the mapping of 71 anchor markers - 40 restriction fragment length polymorphisms (RFLPs), 29 indels or single nucleotide polymorphisms (SNPs) derived from expressed sequence tags (ESTs) and two simple-sequence repeats (SSRs) - on the reference maps of both genera. These markers provided good coverage of the Prunus (78%) and Fragaria (78%) genomes, with maximum gaps and average densities of 22 cM and 7.3 cM/marker in Prunus and 32 cM and 8.0 cM/marker in Fragaria.

CONCLUSION

Our results indicate a clear pattern of synteny, with most markers of each chromosome of one of these species mapping to one or two chromosomes of the other. A large number of rearrangements (36), most of which produced by inversions (27) and the rest (9) by translocations or fission/fusion events could also be inferred. We have provided the first framework for the comparison of the position of genes or DNA sequences of these two economically valuable and yet distantly-related genera of the Rosaceae.

Synteny conservation between two distantly-related Rosaceae genomes: Prunus (the stone fruits) and Fragaria (the strawberry)

BACKGROUND

The Rosaceae encompass a large number of economically-important diploid and polyploid fruit and ornamental species in many different genera. The basic chromosome numbers of these genera are x = 7, 8 and 9 and all have compact and relatively similar genome sizes. Comparative mapping between distantly-related genera has been performed to a limited extent in the Rosaceae including a comparison between Malus (subfamily Maloideae) and Prunus (subfamily Prunoideae); however no data has been published to date comparing Malus or Prunus to a member of the subfamily Rosoideae. In this paper we compare the genome of Fragaria, a member of the Rosoideae, to Prunus, a member of the Prunoideae.

RESULTS

The diploid genomes of Prunus (2n = 2x = 16) and Fragaria (2n = 2x = 14) were compared through the mapping of 71 anchor markers - 40 restriction fragment length polymorphisms (RFLPs), 29 indels or single nucleotide polymorphisms (SNPs) derived from expressed sequence tags (ESTs) and two simple-sequence repeats (SSRs) - on the reference maps of both genera. These markers provided good coverage of the Prunus (78%) and Fragaria (78%) genomes, with maximum gaps and average densities of 22 cM and 7.3 cM/marker in Prunus and 32 cM and 8.0 cM/marker in Fragaria.

CONCLUSION

Our results indicate a clear pattern of synteny, with most markers of each chromosome of one of these species mapping to one or two chromosomes of the other. A large number of rearrangements (36), most of which produced by inversions (27) and the rest (9) by translocations or fission/fusion events could also be inferred. We have provided the first framework for the comparison of the position of genes or DNA sequences of these two economically valuable and yet distantly-related genera of the Rosaceae.

Meta-analysis of polyploid cotton QTL shows unequal contributions of subgenomes to a complex network of genes and gene clusters implicated in lint fiber development

QTL mapping experiments yield heterogeneous results due to the use of different genotypes, environments, and sampling variation. Compilation of QTL mapping results yields a more complete picture of the genetic control of a trait and reveals patterns in organization of trait variation. A total of 432 QTL mapped in one diploid and 10 tetraploid interspecific cotton populations were aligned using a reference map and depicted in a CMap resource. Early demonstrations that genes from the non-fiber-producing diploid ancestor contribute to tetraploid lint fiber genetics gain further support from multiple populations and environments and advanced-generation studies detecting QTL of small phenotypic effect. Both tetraploid subgenomes contribute QTL at largely non-homeologous locations, suggesting divergent selection acting on many corresponding genes before and/or after polyploid formation. QTL correspondence across studies was only modest, suggesting that additional QTL for the target traits remain to be discovered. Crosses between closely-related genotypes differing by single-gene mutants yield profoundly different QTL landscapes, suggesting that fiber variation involves a complex network of interacting genes. Members of the lint fiber development network appear clustered, with cluster members showing heterogeneous phenotypic effects. Meta-analysis linked to synteny-based and expression-based information provides clues about specific genes and families involved in QTL networks.

Chromosome structural changes in diploid and tetraploid A genomes of Gossypium

The genus Gossypium, which comprises a divergent group of diploid species and several recently formed allotetraploids, offers an excellent opportunity to study polyploid genome evolution. In this study, chromosome structural variation among the A, At, and D genomes of Gossypium was evaluated by comparative genetic linkage mapping. We constructed a fully resolved RFLP linkage map for the diploid A genome consisting of 275 loci using an F2 interspecific Gossypium arboreum x Gossypium herbaceum family. The 13 chromosomes of the A genome are represented by 12 large linkage groups in our map, reflecting an expected interchromosomal translocation between G. arboreum and G. herbaceum. The A-genome chromosomes are largely collinear with the D genomes, save for a few small inversions. Although the 2 diploid mapping parents represent the closest living relatives of the allotetraploid At-genome progenitor, 2 translocations and 7 inversions were observed between the A and At genomes. The recombination rates are similar between the 2 diploid genomes; however, the At genome shows a 93% increase in recombination relative to its diploid progenitors. Elevated recombination in the Dt genome was reported previously. These data on the At genome thus indicate that elevated recombination was a general property of allotetraploidy in cotton.

Analysis of papaya BAC end sequences reveals first insights into the organization of a fruit tree genome

Papaya (Carica papaya L.) is a major tree fruit crop of tropical and subtropical regions with an estimated genome size of 372 Mbp. We present the analysis of 4.7% of the papaya genome based on BAC end sequences (BESs) representing 17 million high-quality bases. Microsatellites discovered in 5,452 BESs and flanking primer sequences are available to papaya breeding programs at http://www.genomics.hawaii.edu/papaya/BES . Sixteen percent of BESs contain plant repeat elements, the vast majority (83.3%) of which are class I retrotransposons. Several novel papaya-specific repeats were identified. Approximately 19.1% of the BESs have homology to Arabidopsis cDNA. Increasing numbers of completely sequenced plant genomes and BES projects enable novel approaches to comparative plant genomics. Paired BESs of Carica, Arabidopsis, Populus, Brassica and Lycopersicon were mapped onto the completed genomes of Arabidopsis and Populus. In general the level of microsynteny was highest between closely related organisms. However, papaya revealed a higher degree of apparent synteny with the more distantly related poplar than with the more closely related Arabidopsis. This, as well as significant colinearity observed between peach and poplar genome sequences, support recent observations of frequent genome rearrangements in the Arabidopsis lineage and suggest that the poplar genome sequence may be more useful for elucidating the papaya and other rosid genomes. These insights will play a critical role in selecting species and sequencing strategies that will optimally represent crop genomes in sequence databases.

Leafing through the genomes of our major crop plants: strategies for capturing unique information

Crop plants not only have economic significance, but also comprise important botanical models for evolution and development. This is reflected by the recent increase in the percentage of publicly available sequence data that are derived from angiosperms. Further genome sequencing of the major crop plants will offer new learning opportunities, but their large, repetitive, and often polyploid genomes present challenges. Reduced-representation approaches - such as EST sequencing, methyl filtration and Cot-based cloning and sequencing - provide increased efficiency in extracting key information from crop genomes without full-genome sequencing. Combining these methods with phylogenetically stratified sampling to allow comparative genomic approaches has the potential to further accelerate progress in angiosperm genomics.

Cotton genome mapping with new microsatellites from Acala ‘Maxxa’ BAC-ends

Fine mapping and positional cloning will eventually improve with the anchoring of additional markers derived from genomic clones such as BACs. From 2,603 new BAC-end genomic sequences from Gossypium hirsutum Acala 'Maxxa', 1,316 PCR primer pairs (designated as MUSB) were designed to flank microsatellite or simple sequence repeat motif sequences. Most (1164 or 88%) MUSB primer pairs successfully amplified DNA from three species of cotton with an average of three amplicons per marker and 365 markers (21%) were polymorphic between G. hirsutum and G. barbadense. An interspecific RIL population developed from the above two entries was used to map 433 marker loci and 46 linkage groups with a genetic distance of 2,126.3 cM covering approximately 45% of the cotton genome and an average distance between two loci of 4.9 cM. Based on genome-specific chromosomes identified in G. hirsutum tetraploid (A and D), 56.9% of the coverage was located on the A subgenome while 39.7% was assigned to the D subgenome in the genetic map, suggesting that the A subgenome may be more polymorphic and recombinationally active than originally thought. The linkage groups were assigned to 23 of the 26 chromosomes. This is the first genetic map in which the linkage groups A01 and A02/D03 have been assigned to specific chromosomes. In addition the MUSB-derived markers from BAC-end sequences markers allows fine genetic and QTL mapping of important traits and for the first time provides reconciliation of the genetic and physical maps. Limited QTL analyses suggested that loci on chromosomes 2, 3, 12, 15 and 18 may affect variation in fiber quality traits. The original BAC clones containing the newly mapped MUSB that tag the QTLs provide critical DNA regions for the discovery of gene sequences involved in biological processes such as fiber development and pest resistance in cotton.

A framework linkage map of bermudagrass (Cynodon dactylon x transvaalensis) based on single-dose restriction fragments

This study describes the first detailed linkage maps of two bermudagrass species, Cynodon dactylon (T89) and Cynodon transvaalensis (T574), based on single-dose restriction fragments (SDRFs). The mapping population consisted of 113 F1 progeny of a cross between the two parents. Loci were generated using 179 bermudagrass genomic clones and 50 heterologous cDNAs from Pennisetum and rice. The map of T89 is based on 155 SDRFs and 17 double-dose restriction fragments on 35 linkage groups, with an average marker spacing of 15.3 cM. The map of T574 is based on 77 SDRF loci on 18 linkage groups with an average marker spacing of 16.5 cM. About 16 T89 linkage groups were arranged into four complete and eight into four incomplete homologous sets, while 15 T574 linkage groups were arranged into seven complete homologous sets, all on the basis of multi-locus probes and repulsion linkages. Eleven T89 and three T574 linkage groups remain unassigned. In each parent consensus maps were built based on alignments of homologous linkage groups. Four ancestral chromosomes were inferred after aligning T89 and T574 parental consensus maps using multi-locus probes. The inferred ancestral marker orders were used in comparisons to a detailed Sorghum linkage map using 40 common probes, and to the rice genome sequence using 98 significant BLAST hits, to find regions of colinearity. Using these maps we have estimated the recombinational length of the T89 and T574 genomes at 3,012 and 1,569 cM, respectively, which are 61 and 62% covered by our maps.

Genetic mapping and comparative analysis of seven mutants related to seed fiber development in cotton

Mapping of genes that play major roles in cotton fiber development is an important step toward their cloning and manipulation, and provides a test of their relationships (if any) to agriculturally-important QTLs. Seven previously identified fiber mutants, four dominant (Li (1), Li (2), N (1) and Fbl) and three recessive (n (2), sma-4(h (a)), and sma-4(fz)), were genetically mapped in six F(2) populations comprising 124 or more plants each. For those mutants previously assigned to chromosomes by using aneuploids or by linkage to other morphological markers, all map locations were concordant except n (2), which mapped to the homoeolog of the chromosome previously reported. Three mutations with primary effects on fuzz fibers (N (1), Fbl, n (2)) mapped near the likelihood peaks for QTLs that affected lint fiber productivity in the same populations, perhaps suggesting pleiotropic effects on both fiber types. However, only Li (1) mapped within the likelihood interval for 191 previously detected lint fiber QTLs discovered in non-mutant crosses, suggesting that these mutations may occur in genes that played early roles in cotton fiber evolution, and for which new allelic variants are quickly eliminated from improved germplasm. A close positional association between sma-4(h ( a )), two leaf and stem-borne trichome mutants (t (1) , t (2)), and a gene previously implicated in fiber development, sucrose synthase, raises questions about the possibility that these genes may be functionally related. Increasing knowledge of the correspondence of the cotton and Arabidopsis genomes provides several avenues by which genetic dissection of cotton fiber development may be accelerated.