Next-generation transcriptome sequencing, SNP discovery and validation in four market classes of peanut, Arachis hypogaea L.

Single-nucleotide polymorphisms, which can be identified in the thousands or millions from comparisons of transcriptome or genome sequences, are ideally suited for making high-resolution genetic maps, investigating population evolutionary history, and discovering marker-trait linkages. Despite significant results from their use in human genetics, progress in identification and use in plants, and particularly polyploid plants, has lagged. As part of a long-term project to identify and use SNPs suitable for these purposes in cultivated peanut, which is tetraploid, we generated transcriptome sequences of four peanut cultivars, namely OLin, New Mexico Valencia C, Tamrun OL07 and Jupiter, which represent the four major market classes of peanut grown in the world, and which are important economically to the US southwest peanut growing region. CopyDNA libraries of each genotype were used to generate 2 × 54 paired-end reads using an Illumina GAIIx sequencer. Raw reads were mapped to a custom reference consisting of Tifrunner 454 sequences plus peanut ESTs in GenBank, compromising 43,108 contigs; 263,840 SNP and indel variants were identified among four genotypes compared to the reference. A subset of 6 variants was assayed across 24 genotypes representing four market types using KASP chemistry to assess the criteria for SNP selection. Results demonstrated that transcriptome sequencing can identify SNPs usable as selectable DNA-based markers in complex polyploid species such as peanut. Criteria for effective use of SNPs as markers are discussed in this context.

Sampling nucleotide diversity in cotton.

BACKGROUND: Cultivated cotton is an annual fiber crop derived mainly from two perennial species, Gossypium hirsutum L. or upland cotton, and G. barbadense L., extra long-staple fiber Pima or Egyptian cotton. These two cultivated species are among five allotetraploid species presumably derived monophyletically between G. arboreum and G. raimondii. Genomic-based approaches have been hindered by the limited variation within species. Yet, population-based methods are being used for genome-wide introgression of novel alleles from G. mustelinum and G. tomentosum into G. hirsutum using combinations of backcrossing, selfing, and inter-mating. Recombinant inbred line populations between genetics standards TM-1, (G. hirsutum) x 3-79 (G. barbadense) have been developed to allow high-density genetic mapping of traits. RESULTS: This paper describes a strategy to efficiently characterize genomic variation (SNPs and indels) within and among cotton species. Over 1000 SNPs from 270 loci and 279 indels from 92 loci segregating in G. hirsutum and G. barbadense were genotyped across a standard panel of 24 lines, 16 of which are elite cotton breeding lines and 8 mapping parents of populations from six cotton species. Over 200 loci were genetically mapped in a core mapping population derived from TM-1 and 3-79 and in G. hirsutum breeding germplasm. CONCLUSION: In this research, SNP and indel diversity is characterized for 270 single-copy polymorphic loci in cotton. A strategy for SNP discovery is defined to pre-screen loci for copy number and polymorphism. Our data indicate that the A and D genomes in both diploid and tetraploid cotton remain distinct from each such that paralogs can be distinguished. This research provides mapped DNA markers for intra-specific crosses and introgression of exotic germplasm in cotton.

Aberrant Expression of Critical Genes during Secondary Cell Wall Biogenesis in a Cotton Mutant, Ligon Lintless-1 (Li-1).

Over ninety percent of the value of cotton comes from its fiber; however, the genetic mechanisms governing fiber development are poorly understood. Due to their biochemical and morphological diversity in fiber cells cotton fiber mutants have been useful in examining fiber development; therefore, using the Ligon Lintless (Li-1) mutant, a monogenic dominant cotton mutant with very short fibers, we employed the high throughput approaches of microarray technology and real time PCR to gain insights into what genes were critical during the secondary cell wall synthesis stage. Comparative transcriptome analysis of the normal TM-1 genotype and the near isogenic Li-1 revealed that over 100 transcripts were differentially expressed at least 2-fold during secondary wall biogenesis, although the genetic profile of the expansion phase showed no significant differences in the isolines. Of particular note, we identified three candidate gene families-expansin, sucrose synthase, and tubulin-whose expression in Li-1 deviates from normal expression patterns of its parent, TM-1. These genes may contribute to retarded growth of fibers in Li-1 since they are fiber-expressed structural and metabolic genes. This work provides more details into the mechanisms of fiber development, and suggests the Li gene is active during the later stages of fiber development.

Gene-rich islands for fiber development in the cotton genome.

Cotton fiber is an economically important seed trichome and the world's leading natural fiber used in the manufacture of textiles. As a step toward elucidating the genomic organization and distribution of gene networks responsible for cotton fiber development, we investigated the distribution of fiber genes in the cotton genome. Results revealed the presence of gene-rich islands for fiber genes with a biased distribution in the tetraploid cotton (Gossypium hirsutum L.) genome that was also linked to discrete fiber developmental stages based on expression profiles. There were 3 fiber gene-rich islands associated with fiber initiation on chromosome 5, 3 islands for the early to middle elongation stage on chromosome 10, 3 islands for the middle to late elongation stage on chromosome 14, and 1 island on chromosome 15 for secondary cell wall deposition, for a total of 10 fiber gene-rich islands. Clustering of functionally related gene clusters in the cotton genome displaying similar transcriptional regulation indicates an organizational hierarchy with significant implications for the genetic enhancement of particular fiber quality traits. The relationship between gene-island distribution and functional expression profiling suggests for the first time the existence of functional coupling gene clusters in the cotton genome.

Double feature selection and cluster analyses in mining of microarray data from cotton.

BACKGROUND: Cotton fiber is a single-celled seed trichome of major biological and economic importance. In recent years, genomic approaches such as microarray-based expression profiling were used to study fiber growth and development to understand the developmental mechanisms of fiber at the molecular level. The vast volume of microarray expression data generated requires a sophisticated means of data mining in order to extract novel information that addresses fundamental questions of biological interest. One of the ways to approach microarray data mining is to increase the number of dimensions/levels to the analysis, such as comparing independent studies from different genotypes. However, adding dimensions also creates a challenge in finding novel ways for analyzing multi-dimensional microarray data. RESULTS: Mining of independent microarray studies from Pima and Upland (TM1) cotton using double feature selection and cluster analyses identified species-specific and stage-specific gene transcripts that argue in favor of discrete genetic mechanisms that govern developmental programming of cotton fiber morphogenesis in these two cultivated species. Double feature selection analysis identified the highest number of differentially expressed genes that distinguish the fiber transcriptomes of developing Pima and TM1 fibers. These results were based on the finding that differences in fibers harvested between 17 and 24 day post-anthesis (dpa) represent the greatest expressional distance between the two species. This powerful selection method identified a subset of genes expressed during primary (PCW) and secondary (SCW) cell wall biogenesis in Pima fibers that exhibits an expression pattern that is generally reversed in TM1 at the same developmental stage. Cluster and functional analyses revealed that this subset of genes are primarily regulated during the transition stage that overlaps the termination of PCW and onset of SCW biogenesis, suggesting that these particular genes play a major role in the genetic mechanism that underlies the phenotypic differences in fiber traits between Pima and TM1. CONCLUSION: The novel application of double feature selection analysis led to the discovery of species- and stage-specific genetic expression patterns, which are biologically relevant to the genetic programs that underlie the differences in the fiber phenotypes in Pima and TM1. These results promise to have profound impacts on the ongoing efforts to improve cotton fiber traits.

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.

A 3347-locus genetic recombination map of sequence-tagged sites reveals features of genome organization, transmission and evolution of cotton (Gossypium).

We report genetic maps for diploid (D) and tetraploid (AtDt) Gossypium genomes composed of sequence-tagged sites (STS) that foster structural, functional, and evolutionary genomic studies. The maps include, respectively, 2584 loci at 1.72-cM ( approximately 600 kb) intervals based on 2007 probes (AtDt) and 763 loci at 1.96-cM ( approximately 500 kb) intervals detected by 662 probes (D). Both diploid and tetraploid cottons exhibit negative crossover interference; i.e., double recombinants are unexpectedly abundant. We found no major structural changes between Dt and D chromosomes, but confirmed two reciprocal translocations between At chromosomes and several inversions. Concentrations of probes in corresponding regions of the various genomes may represent centromeres, while genome-specific concentrations may represent heterochromatin. Locus duplication patterns reveal all 13 expected homeologous chromosome sets and lend new support to the possibility that a more ancient polyploidization event may have predated the A-D divergence of 6-11 million years ago. Identification of SSRs within 312 RFLP sequences plus direct mapping of 124 SSRs and exploration for CAPS and SNPs illustrate the "portability" of these STS loci across populations and detection systems useful for marker-assisted improvement of the world's leading fiber crop. These data provide new insights into polyploid evolution and represent a foundation for assembly of a finished sequence of the cotton genome.

Identification of salt responsive genes using comparative microarray analysis in Upland cotton (Gossypium hirsutum L.).

Salinity negatively impacts plant growth and productivity, and little is known about salt responsive genes in cotton. In this study, an intra-specific backcross population of cotton (Gossypium hirsutum L.) was treated with 200 mM NaCl after which differentially expressed genes were identified by comparison between salt tolerant and susceptible segregant bulks using comparative microarray analysis. Microarray analysis identified 720 salt-responsive genes, of which 695 were down-regulated and only 25 were up-regulated in the salt tolerant bulk. Gene ontology of annotated genes revealed that at least some of the identified salt responsive transcripts belong to pathways known to be associated with salt stress including osmolyte and lipid metabolism, cell wall structure, and membrane synthesis. About 48% of all salt-responsive genes were functionally unknown. Quantitative RT-PCR was used to validate 17 selected salt responsive genes. This work represents the first study in employing microarray to investigate the possible mechanisms of the salt response in cotton. Further analysis of salt-responsive genes associated with salt tolerance in cotton will assist in laying a foundation for molecular manipulation in development of new cultivars with improved salt tolerance.

Transcriptome profiling, sequence characterization, and SNP-based chromosomal assignment of the EXPANSIN genes in cotton.

The knowledge of biological significance associated with DNA markers is very limited in cotton. SNPs are potential functional marker to tag genes of biological importance. Plant expansins are a group of extracellular proteins that directly modify the mechanical properties of cell walls, enable turgor-driven cell extension, and likely affect length and quality of cotton fibers. Here, we report the expression profiles of EXPANSIN transcripts during fiber elongation and the discovery of SNP markers, assess the SNP characteristics, and localize six EXPANSIN A genes to chromosomes. Transcriptome profiling of cotton fiber oligonucleotide microarrays revealed that seven EXPANSIN transcripts were differentially expressed when there was parallel polar elongation during morphogenesis at early stage of fiber development, suggesting that major and minor isoforms perform discrete functions during polar elongation and lateral expansion. Ancestral and homoeologous relationships of the six EXPANSIN A genes were revealed by phylogenetic grouping and comparison to extant A- and D-genome relatives of contemporary AD-genome cottons. The average rate of SNP per nucleotide was 2.35% (one SNP per 43 bp), with 1.74 and 3.99% occurring in coding and noncoding regions, respectively, in the selected genotypes. An unequal evolutionary rate of the EXPANSIN A genes at the subgenome level of tetraploid cotton was recorded. Chromosomal locations for each of the six EXPANSIN A genes were established by gene-specific SNP markers. Results revealed a strategy for discovering SNP markers in a polyploidy species like cotton. These markers could be useful to associate candidate genes with the complex fiber traits in MAS.

A global assembly of cotton ESTs.

Approximately 185,000 Gossypium EST sequences comprising >94,800,000 nucleotides were amassed from 30 cDNA libraries constructed from a variety of tissues and organs under a range of conditions, including drought stress and pathogen challenges. These libraries were derived from allopolyploid cotton (Gossypium hirsutum; A(T) and D(T) genomes) as well as its two diploid progenitors, Gossypium arboreum (A genome) and Gossypium raimondii (D genome). ESTs were assembled using the Program for Assembling and Viewing ESTs (PAVE), resulting in 22,030 contigs and 29,077 singletons (51,107 unigenes). Further comparisons among the singletons and contigs led to recognition of 33,665 exemplar sequences that represent a nonredundant set of putative Gossypium genes containing partial or full-length coding regions and usually one or two UTRs. The assembly, along with their UniProt BLASTX hits, GO annotation, and Pfam analysis results, are freely accessible as a public resource for cotton genomics. Because ESTs from diploid and allotetraploid Gossypium were combined in a single assembly, we were in many cases able to bioinformatically distinguish duplicated genes in allotetraploid cotton and assign them to either the A or D genome. The assembly and associated information provide a framework for future investigation of cotton functional and evolutionary genomics.

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.

Genetic mapping of new cotton fiber loci using EST-derived microsatellites in an interspecific recombinant inbred line cotton population.

There is an immediate need for a high-density genetic map of cotton anchored with fiber genes to facilitate marker-assisted selection (MAS) for improved fiber traits. With this goal in mind, genetic mapping with a new set of microsatellite markers [comprising both simple (SSR) and complex (CSR) sequence repeat markers] was performed on 183 recombinant inbred lines (RILs) developed from the progeny of the interspecific cross Gossypium hirsutum L. cv. TM1 x Gossypium barbadense L. Pima 3-79. Microsatellite markers were developed using 1557 ESTs-containing SSRs (> or = 10 bp) and 5794 EST-containing CSRs (> or = 12 bp) obtained from approximately 14,000 consensus sequences derived from fiber ESTs generated from the cultivated diploid species Gossypium arboreum L. cv AKA8401. From a total of 1232 EST-derived SSR (MUSS) and CSR (MUCS) primer-pairs, 1019 (83%) successfully amplified PCR products from a survey panel of six Gossypium species; 202 (19.8%) were polymorphic between the G. hirsutum L. and G. barbadense L. parents of the interspecific mapping population. Among these polymorphic markers, only 86 (42.6%) showed significant sequence homology to annotated genes with known function. The chromosomal locations of 36 microsatellites were associated with 14 chromosomes and/or 13 chromosome arms of the cotton genome by hypoaneuploid deficiency analysis, enabling us to assign genetic linkage groups (LG) to specific chromosomes. The resulting genetic map consists of 193 loci, including 121 new fiber loci not previously mapped. These fiber loci were mapped to 19 chromosomes and 11 LG spanning 1277 cM, providing approximately 27% genome coverage. Preliminary quantitative trait loci analysis suggested that chromosomes 2, 3, 15, and 18 may harbor genes for traits related to fiber quality. These new PCR-based microsatellite markers derived from cotton fiber ESTs will facilitate the development of a high-resolution integrated genetic map of cotton for structural and functional study of fiber genes and MAS of genes that enhance fiber quality.

Characterization of the V-type H((+))-ATPase in the resurrection plant Tortula ruralis: accumulation and polysomal recruitment of the proteolipid c subunit in response to salt-stress.

Tortula ruralis is an important experimental system for the study of plant vegetative desiccation tolerance. EST gene discovery efforts utilizing desiccated gametophytes have identified a cDNA Vac1 encoding a predicted polypeptide with significant similarity to the vacuolar H(+)-ATPase c subunit. VAC1, the 167 amino acid deduced polypeptide, has a predicted molecular mass of 16.9 kDa, and a predicted pI of 9.7. Phylogenetic analysis demonstrated that previously characterized proteolipid polypeptide sequences could be reproducibly grouped into two major clades and that VAC1 forms a discrete evolutionary group. RNA blot and Western blot hybridizations were used to analyse expression of Vac1 and accumulation of VAC1 in response to (1) desiccation and rehydration, (2) increased NaCl concentration, and (3) NaCl-shock. During a desiccation-rehydration cycle, Vac1 transcripts are expressed in both the total and polysomal RNA fractions in approximately equal amounts, and the steady-state transcript levels are unchanged. However, Vac1 transcript levels increased in response to both elevated NaCl concentration and NaCl-shock. There is a preferential accumulation of Vac1 transcripts within the polysomal RNA fraction in response to salt stress, and these data suggest that T. ruralis possesses a salinity-stress-dependent and desiccation-stress-independent mechanism for post-transcriptional gene control. Using a cotton anti-c subunit polyclonal antibody raised against the C-terminal domain, it was shown that the amount of Tortula 16 kDa proteolipid in the tonoplast protein fraction was unaffected by any stress treatment.

Functional genomics of cell elongation in developing cotton fibers.

Cotton fibers are single-celled seed trichomes of major economic importance. Factors that regulate the rate and duration of cell expansion control fiber morphology and important agronomic traits. For genetic characterization of rapid cell elongation in cotton fibers, approximately 14,000 unique genes were assembled from 46,603 expressed sequence tags (ESTs) from developmentally staged fiber cDNAs of a cultivated diploid species ( Gossypium arboreum L.). Conservatively, the fiber transcriptome represents 35-40% of the genes in the cotton genome. In silico expression analysis revealed that rapidly elongating fiber cells exhibit significant metabolic activity, with the bulk of gene transcripts, represented by three major functional groups - cell wall structure and biogenesis, the cytoskeleton and energy/carbohydrate metabolism. Oligonucleotide microarrays revealed dynamic changes in gene expression between primary and secondary cell wall biogenesis showing that fiber genes in the dbEST are highly stage-specific for cell expansion - a conclusion supported by the absence of known secondary cell wall-specific genes from our fiber dbEST. During the developmental switch from primary to secondary cell wall syntheses, 2553 "expansion-associated" fiber genes are significantly down regulated. Genes (81) significantly up-regulated during secondary cell wall synthesis are involved in cell wall biogenesis and energy/carbohydrate metabolism, which is consistent with the stage of cellulose synthesis during secondary cell wall modification in developing fibers. This work provides the first in-depth view of the genetic complexity of the transcriptome of an expanding cell, and lays the groundwork for studying fundamental biological processes in plant biology with applications in agricultural biotechnology.

Rate variation among nuclear genes and the age of polyploidy in Gossypium.

Molecular evolutionary rate variation in Gossypium (cotton) was characterized using sequence data for 48 nuclear genes from both genomes of allotetraploid cotton, models of its diploid progenitors, and an outgroup. Substitution rates varied widely among the 48 genes, with silent and replacement substitution levels varying from 0.018 to 0.162 and from 0.000 to 0.073, respectively, in comparisons between orthologous Gossypium and outgroup sequences. However, about 90% of the genes had silent substitution rates spanning a more narrow threefold range. Because there was no evidence of rate heterogeneity among lineages for any gene and because rates were highly correlated in independent tests, evolutionary rate is inferred to be a property of each gene or its genetic milieu rather than the clade to which it belongs. Evidence from approximately 200,000 nucleotides (40,000 per genome) suggests that polyploidy in Gossypium led to a modest enhancement in rates of nucleotide substitution. Phylogenetic analysis for each gene yielded the topology expected from organismal history, indicating an absence of gene conversion or recombination among homoeologs subsequent to allopolyploid formation. Using the mean synonymous substitution rate calculated across the 48 genes, allopolyploid cotton is estimated to have formed circa 1.5 million years ago (MYA), after divergence of the diploid progenitors about 6.7 MYA.