Development of 1,030 genomic SSR markers in switchgrass

Genome-Wide Comparative Analysis of Five Amaranthaceae Species Reveals a Large Amount of Repeat Content

Amaranthus is a genus of C4 dicotyledonous herbaceous plant species that are widely distributed in Asia, Africa, Australia, and Europe and are used as grain, vegetables, forages, and ornamental plants. Amaranth species have gained significant attention nowadays as potential sources of nutritious food and industrial products. In this study, we performed a comparative genome analysis of five amaranth species, namely, Amaranthus hypochondriacus, Amaranthus tuberculatus, Amaranthus hybridus, Amaranthus palmeri, and Amaranthus cruentus. The estimated repeat content ranged from 54.49% to 63.26% and was not correlated with the genome sizes. Out of the predicted repeat classes, the majority of repetitive sequences were Long Terminal Repeat (LTR) elements, which account for about 13.91% to 24.89% of all amaranth genomes. Phylogenetic analysis based on 406 single-copy orthologous genes revealed that A. hypochondriacus is most closely linked to A. hybridus and distantly related to A. cruentus. However, dioecious amaranth species, such as A. tuberculatus and A. palmeri, which belong to the subgenera Amaranthus Acnida, have formed their distinct clade. The comparative analysis of genomic data of amaranth species will be useful to identify and characterize agronomically important genes and their mechanisms of action. This will facilitate genomics-based, evolutionary studies, and breeding strategies to design faster, more precise, and predictable crop improvement programs.

Development of EST-SSR markers based on transcriptome sequencing for germplasm evaluation of 65 lilies (Lilium)

BACKGROUND

Lilium genus consists of approximately 100 species and numerous varieties, many of which are interspecific hybrids, which result in a complicated genetic background. The germplasm identification, genetic relationship analysis, and breeding of Lilium rely on exploiting genetic information among different accessions. Hence, an attempt was made to develop new EST-SSR markers and study the molecular divergence among 65 genotypes of Lilium.

METHODS AND RESULTS

A total of 5509 EST-SSRs were identified from the high-throughput sequencing database of L. 'Elodie'. After primer screening, six primer pairs with the most abundant polymorphic bands were selected from 100 primer pairs. Combined with the other 10 reported SSR primers, a total of 16 pairs detected genetic information with an average PIC value of 0.7583. The number of alleles per locus varied from four to 33, the expected heterozygosity varied from 0.3289 to 0.9231, and the observed heterozygosity varied from 0.2857 to 0.5122. Based on the phylogenic results, 22 Asiatic hybrids (A), seven Longiflorum × Asiatic hybrids (LA), as well as two native species were grouped. Eighteen Oriental hybrids (O) and nine Oriental × Trumpet (OT) hybrids, four native species, one LO, and one LL (L. pardalinum × L. longiflorum) variety were grouped.

CONCLUSIONS

Two major clusters were reported and a large number of genotypes were grouped based on UPGMA and STRUCTURE analysis methods. The PIC value as well as other parameters revealed that the EST-SSR markers selected were informative. In addition, the clustering pattern displayed better agreement with the cultivar's pedigree. The newly identified SSRs in this study provides molecular markers for germplasm characterization and genetic diversity for Lilium.

Mapping QTLs for spring green-up, plant vigor, and plant biomass in two lowland switchgrass populations

Switchgrass (Panicum virgatum L.) is an important perennial C4 species due to its large potential for cellulosic bioenergy feedstock production. Identification of quantitative trait loci (QTL) controlling important developmental traits is valuable to understanding the genetic basis and using marker-assisted selection (MAS) in switchgrass breeding. One F1 hybrid population derived from NL94 (♀) × SL93 (♂) and one S1 (first-generation selfed) population from NL94 were used in this study. Both the populations showed significant variations for genotype and genotype by environment interactions for three traits studied: plant vigor, spring green-up, and plant biomass. Plant vigor had strong and positive correlations with plant biomass in both populations. Broad-sense heritability estimates for plant vigor ranged from 0.46 to 0.74 and 0.45 to 0.74 in the hybrid and selfed population, respectively. Spring green-up had similar heritability estimates, 0.42-0.78 in the hybrid population, and 0.47-0.82 in the selfed population. Heritability of plant biomass was 0.54-0.64 in the hybrid population and 0.64-0.74 in the selfed population. Fifteen QTLs for spring green-up, 6 QTLs for plant vigor, and 3 QTLs for biomass yield were detected in the hybrid population, whereas 4 QTLs for spring green-up, 4 QTLs for plant vigor, and 1 QTL for biomass yield were detected in the selfed population. Markers associated with these QTLs can be used in MAS to accelerate switchgrass breeding program. This study provided new information in understanding the genetic control of biomass components and demonstrated substantial heterotic vigor that could be explored for breeding hybrid cultivars in switchgrass.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01296-7.

Full-Length SMRT Transcriptome Sequencing and SSR Analysis of Bactrocera dorsalis (Hendel)

Simple Summary

In this study, a full-length transcriptome was analyzed with single-molecule real-time (SMRT) sequencing, which was first used to discover simple sequence repeat (SSR) genetic markers from B. dorsalis. Moreover, SSR markers from isoforms were screened for the identification of species diversity. These results could provide molecular biology methods for further population research.

Abstract

Bactrocera dorsalis (Hendel), as one of the most notorious and destructive invasive agricultural pests in the world, causes damage to over 250 different types of fruits and vegetables throughout tropical and subtropical areas. PacBio single-molecule real-time (SMRT) sequencing was used to generate the full-length transcriptome data of B. dorsalis. A total of 40,319,890 subreads (76.6 Gb, clean reads) were generated, including 535,241 circular consensus sequences (CCSs) and 386,916 full-length non-concatemer reads (FLNCs). Transcript cluster analysis of the FLNC reads revealed 22,780 high-quality reads (HQs). In total, 12,274 transcripts were functionally annotated based on four different databases. A total of 1978 SSR loci were distributed throughout 1714 HQ transcripts, of which 1926 were complete SSRs and 52 were complex SSRs. Among the total SSR loci, 2–3 nucleotide repeats were dominant, occupying 83.62%, of which di- and tri- nucleotide repeats were 39.38% and 44.24%, respectively. We detected 105 repeat motifs, of which AT/AT (50.19%), AC/GT (39.15%), CAA/TTG (32.46%), and ACA/TGT (10.86%) were the most common in di- and tri-nucleotide repeats. The repeat SSR motifs were 12–190 bp in length, and 1638 (88.02%) were shorter than 20 bp. According to the randomly selected microsatellite sequence, 80 pairs of primers were designed, and 174 individuals were randomly amplified by PCR using primers. The number of primers that had amplification products with clear bands and showed good polymorphism came to 41, indicating that this was a feasible way to explore SSR markers from the transcriptomic data of B. dorsalis. These results lay a foundation for developing highly polymorphic microsatellites for researching the functional genomics, population genetic structure, and genetic diversity of B. dorsalis.

Development of EST-SSR markers for genetic diversity analysis in coconut (Cocos nucifera L.)

Genetic improvement in coconut relies on exploiting the vast existing diversity among coconut accessions. Robust molecular markers are a pre-requisite for efficient characterization of genetic diversity. Microsatellites or simple sequence repeats (SSRs), mined from expressed sequence tags (ESTs), constitute an important resource for analysis of genetic diversity as they are abundant, polymorphic and represent function regions of the genome. We have identified a total of 318,528 putative EST-SSRs from 130,942 unigenes utilizing a leaf transcriptome dataset of coconut. Among the EST-SSRs, dinucleotide repeats were abundant (219,912; 69.04%) followed by trinucleotide (70,722; 22.2%) and tetra-nucleotide repeats (6281; 1.9%). Among the dinucleotide repeat motifs, the dominant repeat was AG/CT (35.87%), followed by AT/AT (18.59%), while the dominant trinucleotide repeat was AAG/CTT (4.59%). One hundred and twenty EST-SSR primer pairs were designed and utilized to amplify six DNA samples of coconut accessions. Fifty primers (41.7%) produced reproducible polymorphic fragments of expected sizes, from which a total of 10 primers were selected for the diversity assessment in 186 palms of 50 coconut accessions, comprising of 25 each of tall and dwarf accessions. A total of 137 alleles were detected with an average of 13.7 alleles per SSR locus. The number of alleles observed at each locus in the data set ranged from 7 to 22. All the loci showed 100% polymorphism with respect to the samples screened. The average observed heterozygosity was 0.46. The PIC values ranged from 0.79 (CnKGDEST129 and CnKGDEST100) to 0.91 (CnKGDEST117 and CnKGDEST122) with a mean value of 0.85, indicating the capacity of the EST-SSR markers to detect high levels of polymorphism. The cluster analysis revealed that accessions were generally clustered based on their relative similarity and irrespective of their geographic origins. The present study demonstrates the usefulness of transcriptome sequencing as a rapid and cost-effective methodology for the development of molecular markers. The EST-SSR markers generated through this study constitute useful and reliable tools for assessment of genetic diversity and marker-assisted selection in coconut.

J, A. S, Bhat KV, Gaikwad AB. Genome wide identification and characterization of microsatellite markers in black pepper (Piper nigrum): A valuable resource for boosting genomics applications

Black pepper is one of the most valued and widely used spices in the world and dominates multi-billion dollar global spices trade. India is amongst the major producers, consumers and exporters of black pepper. In spite of its commercial and cultural importance, black pepper has received meagre attention in terms of generation of genomic resources. Availability of markers distributed throughout the genome would facilitate and accelerate genetic studies, QTL identification, genetic enhancement and crop improvement in black pepper. In this perspective, the sequence information from the recently sequenced black pepper (Piper nigrum) genome has been used for identification and characterisation of Simple Sequence Repeats (SSRs). Total 69,126 SSRs were identified from assembled genomic sequence of P. nigrum. The SSR frequency was 158 per MB making it, one SSR for every 6.3 kb in the assembled genome. Among the different types of microsatellite repeat motifs, dinucleotides were the most abundant (48.6%), followed by trinucleotide (23.7%) and compound repeats (20.62%). A set of 85 SSRs were used for validation, of which 74 produced amplification products of expected size. Genetic diversity of 30 black pepper accessions using 50 SSRs revealed four distinct clusters. Further, the cross species transferability of the SSRs was checked in nine other Piper species. Out of 50 SSRs used, 19 and 31 SSRs were amplified in nine and seven species, respectively. Thus the identified SSRs may have application in other species of the genus Piper where genome sequence is not available yet. Present study reports the first NGS based genomic SSRs in black pepper and thus constitute a valuable resource for a whole fleet of applications in genetics and plant breeding studies such as genetic map construction, QTL identification, map-based gene cloning, marker-assisted selection and evolutionary studies in Piper nigrum and related species.

A Novel Software and Method for the Efficient Development of Polymorphic SSR Loci Based on Transcriptome Data

Traditional methods for developing polymorphic microsatellite loci without reference sequences are time-consuming and labor-intensive, and the polymorphisms of simple sequence repeat (SSR) loci developed from expressed sequence tag (EST) databases are generally poor. To address this issue, in this study, we developed a new software (PSSRdt) and established an effective method for directly obtaining polymorphism details of SSR loci by analyzing diverse transcriptome data. The new method includes three steps, raw data processing, PSSRdt application, and loci extraction and verification. To test the practicality of the method, we successfully obtained 1940 potential polymorphic SSRs from the transcript dataset combined with 44 pea aphid transcriptomes. Fifty-two SSR loci obtained by the new method were selected for validating the polymorphic characteristics by genotyping in pea aphid individuals. The results showed that over 92% of SSR loci were polymorphic and 73.1% of loci were highly polymorphic. Our new software and method provide an innovative approach to microsatellite development based on RNA-seq data, and open a new path for the rapid mining of numerous loci with polymorphism to add to the body of research on microsatellites.

Genome Survey Sequencing of Luffa Cylindrica L. and Microsatellite High Resolution Melting (SSR-HRM) Analysis for Genetic Relationship of Luffa Genotypes

Luffa cylindrica (L.) Roem. is an economically important vegetable crop in China. However, the genomic information on this species is currently unknown. In this study, for the first time, a genome survey of L. cylindrica was carried out using next-generation sequencing (NGS) technology. In total, 43.40 Gb sequence data of L. cylindrica, about 54.94× coverage of the estimated genome size of 789.97 Mb, were obtained from HiSeq 2500 sequencing, in which the guanine plus cytosine (GC) content was calculated to be 37.90%. The heterozygosity of genome sequences was only 0.24%. In total, 1,913,731 contigs (>200 bp) with 525 bp N₅₀ length and 1,410,117 scaffolds (>200 bp) with 885.01 Mb total length were obtained. From the initial assembled L. cylindrica genome, 431,234 microsatellites (SSRs) (≥5 repeats) were identified. The motif types of SSR repeats included 62.88% di-nucleotide, 31.03% tri-nucleotide, 4.59% tetra-nucleotide, 0.96% penta-nucleotide and 0.54% hexa-nucleotide. Eighty genomic SSR markers were developed, and 51/80 primers could be used in both “Zheda 23” and “Zheda 83”. Nineteen SSRs were used to investigate the genetic diversity among 32 accessions through SSR-HRM analysis. The unweighted pair group method analysis (UPGMA) dendrogram tree was built by calculating the SSR-HRM raw data. SSR-HRM could be effectively used for genotype relationship analysis of Luffa species.

Genome wide characterization of simple sequence repeats in watermelon genome and their application in comparative mapping and genetic diversity analysis

BACKGROUND

Microsatellite markers are one of the most informative and versatile DNA-based markers used in plant genetic research, but their development has traditionally been difficult and costly. The whole genome sequencing with next-generation sequencing (NGS) technologies provides large amounts of sequence data to develop numerous microsatellite markers at whole genome scale. SSR markers have great advantage in cross-species comparisons and allow investigation of karyotype and genome evolution through highly efficient computation approaches such as in silico PCR. Here we described genome wide development and characterization of SSR markers in the watermelon (Citrullus lanatus) genome, which were then use in comparative analysis with two other important crop species in the Cucurbitaceae family: cucumber (Cucumis sativus L.) and melon (Cucumis melo L.). We further applied these markers in evaluating the genetic diversity and population structure in watermelon germplasm collections.

RESULTS

A total of 39,523 microsatellite loci were identified from the watermelon draft genome with an overall density of 111 SSRs/Mbp, and 32,869 SSR primers were designed with suitable flanking sequences. The dinucleotide SSRs were the most common type representing 34.09 % of the total SSR loci and the AT-rich motifs were the most abundant in all nucleotide repeat types. In silico PCR analysis identified 832 and 925 SSR markers with each having a single amplicon in the cucumber and melon draft genome, respectively. Comparative analysis with these cross-species SSR markers revealed complicated mosaic patterns of syntenic blocks among the genomes of three species. In addition, genetic diversity analysis of 134 watermelon accessions with 32 highly informative SSR loci placed these lines into two groups with all accessions of C.lanatus var. citorides and three accessions of C. colocynthis clustered in one group and all accessions of C. lanatus var. lanatus and the remaining accessions of C. colocynthis clustered in another group. Furthermore, structure analysis was consistent with the dendrogram indicating the 134 watermelon accessions were classified into two populations.

CONCLUSION

The large number of genome wide SSR markers developed herein from the watermelon genome provides a valuable resource for genetic map construction, QTL exploration, map-based gene cloning and marker-assisted selection in watermelon which has a very narrow genetic base and extremely low polymorphism among cultivated lines. Furthermore, the cross-species transferable SSR markers identified herein should also have practical uses in many applications in species of Cucurbitaceae family whose whole genome sequences are not yet available.

Assessment of genetic diversity, population structure and relationships in Indian and non-Indian genotypes of finger millet (Eleusine coracana (L.) Gaertn) using genomic SSR markers

We evaluated the genetic variation and population structure in Indian and non-Indian genotypes of finger millet using 87 genomic SSR primers. The 128 finger millet genotypes were collected and genomic DNA was isolated. Eighty-seven genomic SSR primers with 60–70 % GC contents were used for PCR analysis of 128 finger millet genotypes. The PCR products were separated and visualized on a 6 % polyacrylamide gel followed by silver staining. The data were used to estimate major allele frequency using Power Marker v3.0. Dendrograms were constructed based on the Jaccard’s similarity coefficient. Statistical fitness and population structure analyses were performed to find the genetic diversity. The mean major allele frequency was 0.92; the means of polymorphic alleles were 2.13 per primer and 1.45 per genotype; the average polymorphism was 59.94 % per primer and average PIC value was 0.44 per primer. Indian genotypes produced an additional 0.21 allele than non-Indian genotypes. Gene diversity was in the range from 0.02 to 0.35. The average heterozygosity was 0.11, close to 100 % homozygosity. The highest inbreeding coefficient was observed with SSR marker UGEP67. The Jaccard’s similarity coefficient value ranged from 0.011 to 0.836. The highest similarity value was 0.836 between genotypes DPI009-04 and GPU-45. Indian genotypes were placed in Eleusine coracana major cluster (EcMC) 1 along with 6 non-Indian genotypes. AMOVA showed that molecular variance in genotypes from various geographical regions was 4 %; among populations it was 3 % and within populations it was 93 %. PCA scatter plot analysis showed that GPU-28, GPU-45 and DPI009-04 were closely dispersed in first component axis. In structural analysis, the genotypes were divided into three subpopulations (SP1, SP2 and SP3). All the three subpopulations had an admixture of alleles and no pure line was observed. These analyses confirmed that all the genotypes were genetically diverse and had been grouped based on their geographic regions.

Construction of a Genetic Linkage Map and Identification of QTLs for Seed Weight and Seed Size Traits in Lentil (Lens culinaris Medik.)

Seed weight and seed size both are quantitative traits and have been considered as important components of grain yield, thus identification of quantitative trait loci (QTL) for seed traits in lentil (Lens culinaris) would be beneficial for the improvement of grain yield. Hence the main objective of this study was to identify QTLs for seed traits using an intraspecific mapping population derived from a cross between L. culinaris cv. Precoz (seed weight-5.1g, seed size-5.7mm) and L. culinaris cv. L830 (seed weight-2.2g, seed size-4mm) comprising 126 F8-RILs. For this, two microsatellite genomic libraries enriched for (GA/CT) and (GAA/CTT) motif were constructed which resulted in the development of 501 new genomic SSR markers. Six hundred forty seven SSR markers (including 146 previously published) were screened for parental polymorphism and 219 (33.8%) were found to be polymorphic among the parents. Of these 216 were mapped on seven linkage groups at LOD4.0 spanning 1183.7cM with an average marker density of 5.48cM. Phenotypic data from the RILs was used to identify QTLs for the seed weight and seed size traits by single marker analysis (SMA) followed by composite interval mapping (CIM) which resulted in one QTL each for the 2 traits (qSW and qSS) that were co-localized on LG4 and explained 48.4% and 27.5% of phenotypic variance respectively. The current study would serve as a strong foundation for further validation and fine mapping for utilization in lentil breeding programs.

Development of genic SSR markers from transcriptome sequencing of pear buds

A total of 8375 genic simple sequence repeat (SSR) loci were discovered from a unigene set assembled from 116282 transcriptomic unigenes in this study. Dinucleotide repeat motifs were the most common with a frequency of 65.11%, followed by trinucleotide (32.81%). A total of 4100 primer pairs were designed from the SSR loci. Of these, 343 primer pairs (repeat length ≥15 bp) were synthesized with an M13 tail and tested for stable amplification and polymorphism in four Pyrus accessions. After the preliminary test, 104 polymorphic genic SSR markers were developed; dinucleotide and trinucleotide repeats represented 97.11% (101) of these. Twenty-eight polymorphic genic SSR markers were selected randomly to further validate genetic diversity among 28 Pyrus accessions. These markers displayed a high level of polymorphism. The number of alleles at these SSR loci ranged from 2 to 17, with a mean of 9.43 alleles per locus, and the polymorphism information content (PIC) values ranged from 0.26 to 0.91. The UPGMA (unweighted pair-group method with arithmetic average) cluster analysis grouped the 28 Pyrus accessions into two groups: Oriental pears and Occidental pears, which are congruent to the traditional taxonomy, demonstrating their effectiveness in analyzing Pyrus phylogenetic relationships, enriching rare Pyrus EST-SSR resources, and confirming the potential value of a pear transcriptome database for the development of new SSR markers.

Isolation and characterization of microsatellite markers and analysis of genetic diversity in Chinese jujube (Ziziphus jujuba Mill.)

Chinese jujube (Ziziphus jujuba Mill, 2n = 2× = 24, Rhamnaceae) is an economically important Chinese native species. It has high nutritional value, and its medicinal properties have led to extensive use in traditional oriental medicine. The characterization of genotypes using molecular markers is important for genetic studies and plant breeding. However, few simple sequence repeat (SSR) markers are available for this species. In this study, 1,488 unique SSR clones were isolated from Z. jujuba 'Dongzao' using enriched genomic libraries coupled with a three-primer colony PCR screening strategy, yielding a high enrichment rate of 73.3%. Finally, 1,188 (80.87%) primer pairs were amplified successfully in the size expected for 'Dongzao'. A total of 350 primer pairs were further selected and evaluated for their ability to detect polymorphisms across a panel of six diverse cultivars; among these, 301 primer pairs detected polymorphisms, and the polymorphism information content (PIC) value across all loci ranged from 0.15 to 0.82, with an average of 0.52. An analysis of 76 major cultivars employed in Chinese jujube production using 31 primer pairs revealed comparatively high genetic diversity among these cultivars. Within-population differences among individuals accounted for 98.2% of the observed genetic variation. Neighbor-joining clustering divided the cultivars into three main groups, none of which correspond to major geographic regions, suggesting that the genetics and geographical origin of modern Chinese jujube cultivars might not be linked. The current work firstly reports the large-scale development of Chinese jujube SSR markers. The development of these markers and their polymorphic information represent a significant improvement in the available Chinese jujube genomic resources and will facilitate both genetic and breeding applications, further accelerating the development of new cultivars.

Expanding the repertoire of microsatellite markers for polymorphism studies in Indian accessions of mung bean (Vigna radiata L. Wilczek)

Limited availability of validated, polymorphic microsatellite markers in mung bean (Vigna radiata), an important food legume of India, has been a major hurdle towards its improvement and higher yield. The present study was undertaken in order to develop a new set of microsatellite markers and utilize them for the analysis of genetic diversity within mung bean accessions from India. A GA/CT enriched library was constructed from V. radiata which resulted in 1,250 putative recombinant clones of which 850 were sequenced. SSR motifs were identified and their flanking sequences were utilized to design 328 SSR primer pairs. Of these, 48 SSR markers were employed for assessing genetic diversity among 76 mung bean accessions from various geographical locations in India. Two hundred and thirty four alleles with an average of 4.85 alleles per locus were detected at 48 loci. The polymorphic information content (PIC) per locus varied from 0.1 to 0.88 (average: 0.49 per locus). The observed and expected heterozygosities ranged from 0.40 to 0.95 and 0.40 to 0.81 respectively. Based on Jaccard's similarity matrix, a dendrogram was constructed using the unweighted pair-group method with arithmetic averages (UPGMA) analysis which revealed that one accession from Bundi, Rajasthan was clustered out separately while remaining accessions were grouped into two major clusters. The markers generated in this study will help in expanding the repertoire of the available SSR markers thereby facilitating analysis of genetic diversity, molecular mapping and ultimately broadening the scope for genetic improvement of this legume.

Exploring genetic variability within lentil (Lens culinaris Medik.) and across related legumes using a newly developed set of microsatellite markers

Lentil (Lens culinaris Medik.) is an economically important grain legume, yet the genetic and genomic resources remain largely uncharacterized and unexploited in this crop. Microsatellites have become markers of choice for crop improvement applications. Hence, simple sequence repeat (SSR) markers were developed for lentil through the construction of genomic library enriched for GA/CT motifs. As a result 122 functional SSR primer pairs were developed from 151 microsatellite loci and validated in L. culinaris cv. Precoz. Thirty three SSR markers were utilized for the analysis of genetic relationships between cultivated and wild species of Lens and related legumes. A total of 123 alleles were amplified at 33 loci ranging from 2-5 alleles with an average of 3.73 alleles per locus. Polymorphic information content (PIC) for all the loci ranged from 0.13 to 0.99 with an average of 0.66 per locus. Varied levels of cross genera transferability were obtained ranging from 69.70 % across Pisum sativum to 12.12 % across Vigna radiata. The UPGMA based dendrogram was able to establish the uniqueness of each genotype and grouped them into two major clusters clearly resolving the genetic relationships within lentil and related species. The new set of SSR markers reported here were efficient and highly polymorphic and would add to the existing repertoire of lentil SSR markers to be utilized in molecular breeding. Moreover, the improved knowledge about intra- and inter-specific genetic relationships would facilitate germplasm utilization for lentil improvement.

Development of simple sequence repeat (SSR) markers of sesame (Sesamum indicum) from a genome survey

Sesame (Sesamum indicum), an important oil crop, is widely grown in tropical and subtropical regions. It provides part of the daily edible oil allowance for almost half of the world's population. A limited number of co-dominant markers has been developed and applied in sesame genetic diversity and germplasm identity studies. Here we report for the first time a whole genome survey used to develop simple sequence repeat (SSR) markers and to detect the genetic diversity of sesame germplasm. From the initial assembled sesame genome, 23,438 SSRs (≥5 repeats) were identified. The most common repeat motif was dinucleotide with a frequency of 84.24%, followed by 13.53% trinucleotide, 1.65% tetranucleotide, 0.3% pentanucleotide and 0.28% hexanucleotide motifs. From 1500 designed and synthesised primer pairs, 218 polymorphic SSRs were developed and used to screen 31 sesame accessions that from 12 countries. STRUCTURE and phylogenetic analyses indicated that all sesame accessions could be divided into two groups: one mainly from China and another from other countries. Cluster analysis classified Chinese major sesame varieties into three groups. These novel SSR markers are a useful tool for genetic linkage map construction, genetic diversity detection, and marker-assisted selective sesame breeding.

Selection and validation of reference genes for gene expression analysis in switchgrass (Panicum virgatum) using quantitative real-time RT-PCR

Switchgrass (Panicum virgatum) has received a lot of attention as a forage and bioenergy crop during the past few years. Gene expression studies are in progress to improve new traits and develop new cultivars. Quantitative real time PCR (qRT-PCR) has emerged as an important technique to study gene expression analysis. For accurate and reliable results, normalization of data with reference genes is essential. In this work, we evaluate the stability of expression of genes to use as reference for qRT-PCR in the grass P. virgatum. Eleven candidate reference genes, including eEF-1α, UBQ6, ACT12, TUB6, eIF-4a, GAPDH, SAMDC, TUA6, CYP5, U2AF, and FTSH4, were validated for qRT-PCR normalization in different plant tissues and under different stress conditions. The expression stability of these genes was verified by the use of two distinct algorithms, geNorm and NormFinder. Differences were observed after comparison of the ranking of the candidate reference genes identified by both programs but eEF-1α, eIF-4a, CYP5 and U2AF are ranked as the most stable genes in the samples sets under study. Both programs discard the use of SAMDC and TUA6 for normalization. Validation of the reference genes proposed by geNorm and NormFinder were performed by normalization of transcript abundance of a group of target genes in different samples. Results show similar expression patterns when the best reference genes selected by both programs were used but differences were detected in the transcript abundance of the target genes. Based on the above research, we recommend the use of different statistical algorithms to identify the best reference genes for expression data normalization. The best genes selected in this study will help to improve the quality of gene expression data in a wide variety of samples in switchgrass.

Development and characterization of genomic SSR markers in Cynodon transvaalensis Burtt-Davy

Simple sequence repeat (SSR) markers are a major molecular tool for genetic and genomic research that have been extensively developed and used in major crops. However, few are available in African bermudagrass (Cynodon transvaalensis Burtt-Davy), an economically important warm-season turfgrass species. African bermudagrass is mainly used for hybridizations with common bermudagrass [C. dactylon var. dactylon (L.) Pers.] in the development of superior interspecific hybrid turfgrass cultivars. Accordingly, the major objective of this study was to develop and characterize a large set of SSR markers. Genomic DNA of C. transvaalensis '4200TN 24-2' from an Oklahoma State University (OSU) turf nursery was extracted for construction of four SSR genomic libraries enriched with [CA](n), [GA](n), [AAG](n), and [AAT](n) as core repeat motifs. A total of 3,064 clones were sequenced at the OSU core facility. The sequences were categorized into singletons and contiguous sequences to exclude redundancy. From the two sequence categories, 1,795 SSR loci were identified. After excluding duplicate SSRs by comparison with previously developed SSR markers using a nucleotide basic local alignment tool, 1,426 unique primer pairs (PPs) were designed. Out of the 1,426 designed PPs, 981 (68.8 %) amplified alleles of the expected size in the donor DNA. Polymorphisms of the SSR PPs tested in eight C. transvaalensis plants were 93 % polymorphic with 544 markers effective in all genotypes. Inheritance of the SSRs was examined in six F(1) progeny of African parents 'T577' × 'Uganda', indicating 917 markers amplified heritable alleles. The SSR markers developed in the study are the first large set of co-dominant markers in African bermudagrass and should be highly valuable for molecular and traditional breeding research.

Advances in biotechnology and genomics of switchgrass

Switchgrass (Panicum virgatum L.) is a C4 perennial warm season grass indigenous to the North American tallgrass prairie. A number of its natural and agronomic traits, including adaptation to a wide geographical distribution, low nutrient requirements and production costs, high water use efficiency, high biomass potential, ease of harvesting, and potential for carbon storage, make it an attractive dedicated biomass crop for biofuel production. We believe that genetic improvements using biotechnology will be important to realize the potential of the biomass and biofuel-related uses of switchgrass. Tissue culture techniques aimed at rapid propagation of switchgrass and genetic transformation protocols have been developed. Rapid progress in genome sequencing and bioinformatics has provided efficient strategies to identify, tag, clone and manipulate many economically-important genes, including those related to higher biomass, saccharification efficiency, and lignin biosynthesis. Application of the best genetic tools should render improved switchgrass that will be more economically and environmentally sustainable as a lignocellulosic bioenergy feedstock.

Development of an integrated transcript sequence database and a gene expression atlas for gene discovery and analysis in switchgrass (Panicum virgatum L.)

Switchgrass (Panicum virgatum L.) is a perennial C4 grass with the potential to become a major bioenergy crop. To help realize this potential, a set of RNA-based resources were developed. Expressed sequence tags (ESTs) were generated from two tetraploid switchgrass genotypes, Alamo AP13 and Summer VS16. Over 11.5 million high-quality ESTs were generated with 454 sequencing technology, and an additional 169 079 Sanger sequences were obtained from the 5' and 3' ends of 93 312 clones from normalized, full-length-enriched cDNA libraries. AP13 and VS16 ESTs were assembled into 77 854 and 30 524 unique transcripts (unitranscripts), respectively, using the Newbler and pave programs. Published Sanger-ESTs (544 225) from Alamo, Kanlow, and 15 other cultivars were integrated with the AP13 and VS16 assemblies to create a universal switchgrass gene index (PviUT1.2) with 128 058 unitranscripts, which were annotated for function. An Affymetrix cDNA microarray chip (Pvi_cDNAa520831) containing 122 973 probe sets was designed from PviUT1.2 sequences, and used to develop a Gene Expression Atlas for switchgrass (PviGEA). The PviGEA contains quantitative transcript data for all major organ systems of switchgrass throughout development. We developed a web server that enables flexible, multifaceted analyses of PviGEA transcript data. The PviGEA was used to identify representatives of all known genes in the phenylpropanoid-monolignol biosynthesis pathway.

Development of a genome-wide multiple duplex-SSR protocol and its applications for the identification of selfed progeny in switchgrass

Background

Switchgrass (Panicum virgatum) is a herbaceous crop for the cellulosic biofuel feedstock development in the USA and Europe. As switchgrass is a naturally outcrossing species, accurate identification of selfed progeny is important to producing inbreds, which can be used in the production of heterotic hybrids. Development of a technically reliable, time-saving and easily used marker system is needed to quantify and characterize breeding origin of progeny plants of targeted parents.

Results

Genome-wide screening of 915 mapped microsatellite (simple sequence repeat, SSR) markers was conducted, and 842 (92.0%) produced clear and scorable bands on a pooled DNA sample of eight switchgrass varieties. A total of 166 primer pairs were selected on the basis of their relatively even distribution in switchgrass genome and PCR amplification quality on 16 tetraploid genotypes. Mean polymorphic information content value for the 166 markers was 0.810 ranging from 0.116 to 0.959. From them, a core set of 48 loci, which had been mapped on 17 linkage groups, was further tested and optimized to develop 24 sets of duplex markers. Most of (up to 87.5%) targeted, but non-allelic amplicons within each duplex were separated by more than 10-bp. Using the established duplex PCR protocol, selfing ratio (i.e., selfed/all progeny x100%) was identified as 0% for a randomly selected open-pollinated ‘Kanlow’ genotype grown in the field, 15.4% for 22 field-grown plants of bagged inflorescences, and 77.3% for a selected plant grown in a growth chamber.

Conclusions

The study developed a duplex SSR-based PCR protocol consisting of 48 markers, providing ample choices of non-tightly-linked loci in switchgrass whole genome, and representing a powerful, time-saving and easily used method for the identification of selfed progeny in switchgrass. The protocol should be a valuable tool in switchgrass breeding efforts.

Development of SSR markers and construction of a linkage map in jute

Jute is an important natural fibre crop, which is only second to cotton in its importance at the global level. It is mostly grown in Indian subcontinent and has been recently used for the development of genomics resources.We recently initiated a programme to develop simple sequence repeat markers and reported a set of 2469 SSR that were developed using four SSR-enriched libraries (Mir et al. 2009). In this communication, we report an additional set of 607 novel SSR in 393 SSR containing sequences. However, primers could be designed for only 417 potentially useful SSR. Polymorphism survey was carried out for 374 primer pairs using two parental genotypes (JRO 524 and PPO4) of a mapping population developed for fibre fineness; only 66 SSR were polymorphic. Owing to a low level of polymorphism between the parental genotypes and a high degree of segregation distortion in recombinant inbred lines, genotypic data of only 53 polymorphic SSR on the mapping population consisting of 120 RIL could be used for the construction of a linkage map; 36 SSR loci were mapped on six linkage groups that covered a total genetic distance of 784.3 cM. Hopefully, this map will be enriched with more SSR loci in future and will prove useful for identification of quantitative trait loci/genes for molecular breeding involving improvement of fibre fineness and other related traits in jute.

Development of simple sequence repeat (SSR) markers from a genome survey of Chinese bayberry (Myrica rubra)

BACKGROUND

Chinese bayberry (Myrica rubra Sieb. and Zucc.) is a subtropical evergreen tree originating in China. It has been cultivated in southern China for several thousand years, and annual production has reached 1.1 million tons. The taste and high level of health promoting characters identified in the fruit in recent years has stimulated its extension in China and introduction to Australia. A limited number of co-dominant markers have been developed and applied in genetic diversity and identity studies. Here we report, for the first time, a survey of whole genome shotgun data to develop a large number of simple sequence repeat (SSR) markers to analyse the genetic diversity of the common cultivated Chinese bayberry and the relationship with three other Myrica species.

RESULTS

The whole genome shotgun survey of Chinese bayberry produced 9.01Gb of sequence data, about 26x coverage of the estimated genome size of 323 Mb. The genome sequences were highly heterozygous, but with little duplication. From the initial assembled scaffold covering 255 Mb sequence data, 28,602 SSRs (≥5 repeats) were identified. Dinucleotide was the most common repeat motif with a frequency of 84.73%, followed by 13.78% trinucleotide, 1.34% tetranucleotide, 0.12% pentanucleotide and 0.04% hexanucleotide. From 600 primer pairs, 186 polymorphic SSRs were developed. Of these, 158 were used to screen 29 Chinese bayberry accessions and three other Myrica species: 91.14%, 89.87% and 46.84% SSRs could be used in Myrica adenophora, Myrica nana and Myrica cerifera, respectively. The UPGMA dendrogram tree showed that cultivated Myrica rubra is closely related to Myrica adenophora and Myrica nana, originating in southwest China, and very distantly related to Myrica cerifera, originating in America. These markers can be used in the construction of a linkage map and for genetic diversity studies in Myrica species.

CONCLUSION

Myrica rubra has a small genome of about 323 Mb with a high level of heterozygosity. A large number of SSRs were identified, and 158 polymorphic SSR markers developed, 91% of which can be transferred to other Myrica species.

Exploring the switchgrass transcriptome using second-generation sequencing technology

Background

Switchgrass (Panicum virgatum L.) is a C4 perennial grass and widely popular as an important bioenergy crop. To accelerate the pace of developing high yielding switchgrass cultivars adapted to diverse environmental niches, the generation of genomic resources for this plant is necessary. The large genome size and polyploid nature of switchgrass makes whole genome sequencing a daunting task even with current technologies. Exploring the transcriptional landscape using next generation sequencing technologies provides a viable alternative to whole genome sequencing in switchgrass.

Principal Findings

Switchgrass cDNA libraries from germinating seedlings, emerging tillers, flowers, and dormant seeds were sequenced using Roche 454 GS-FLX Titanium technology, generating 980,000 reads with an average read length of 367 bp. De novo assembly generated 243,600 contigs with an average length of 535 bp. Using the foxtail millet genome as a reference greatly improved the assembly and annotation of switchgrass ESTs. Comparative analysis of the 454-derived switchgrass EST reads with other sequenced monocots including Brachypodium, sorghum, rice and maize indicated a 70–80% overlap. RPKM analysis demonstrated unique transcriptional signatures of the four tissues analyzed in this study. More than 24,000 ESTs were identified in the dormant seed library. In silico analysis indicated that there are more than 2000 EST-SSRs in this collection. Expression of several orphan ESTs was confirmed by RT-PCR.

Significance

We estimate that about 90% of the switchgrass gene space has been covered in this analysis. This study nearly doubles the amount of EST information for switchgrass currently in the public domain. The celerity and economical nature of second-generation sequencing technologies provide an in-depth view of the gene space of complex genomes like switchgrass. Sequence analysis of closely related members of the NAD⁺-malic enzyme type C4 grasses such as the model system Setaria viridis can serve as a viable proxy for the switchgrass genome.

A high-density simple sequence repeat-based genetic linkage map of switchgrass

Switchgrass (Panicum virgatum) has been identified as a promising cellulosic biofuel crop in the United States. Construction of a genetic linkage map is fundamental for switchgrass molecular breeding and the elucidation of its genetic mechanisms for economically important traits. In this study, a novel population consisting of 139 selfed progeny of a northern lowland genotype, NL 94 LYE 16X13, was used to construct a linkage map. A total of 2493 simple sequence repeat markers were screened for polymorphism. Of 506 polymorphic loci, 80.8% showed a goodness-of-fit of 1:2:1 segregation ratio. Among 469 linked loci on the framework map, 241 coupling vs. 228 repulsion phase linkages were detected that conformed to a 1:1 ratio, confirming disomic inheritance. A total of 499 loci were mapped to 18 linkage groups (LG), of which the cumulative length was 2085.2 cM, with an average marker interval of 4.2 cM. Nine homeologous LG pairs were identified based on multi-allele markers and comparative genomic analysis. Two clusters of segregation-distorted loci were identified on LG 5b and 9b, respectively. Comparative analysis indicated a one-to-one relationship between nine switchgrass homeologous groups and nine foxtail millet (Setaria italica) chromosomes, suggesting strong homology between the two species. The linkage map derived from selfing a heterozygous parent, instead of two separate maps usually constructed for a cross-fertilized species, provides a new genetic framework to facilitate genomics research, quantitative trait locus (QTL) mapping, and marker-assisted breeding.

Sequence-based novel genomic microsatellite markers for robust genotyping purposes in foxtail millet [Setaria italica (L.) P. Beauv]

The unavailability of microsatellite markers and saturated genetic linkage map has restricted the genetic improvement of foxtail millet [Setaria italica (L.) P. Beauv.], despite the fact that in recent times it has been documented as a new model species for biofuel grasses. With the objective to generate a good number of microsatellite markers in foxtail millet cultivar 'Prasad', 690 clones were sequenced which generated 112.95 kb high quality sequences obtained from three genomic libraries each enriched with different microsatellite repeat motifs. Microsatellites were identified in 512 (74.2%) of the 690 positive clones and 172 primer pairs (pp) were successfully designed from 249 (48.6%) unique SSR-containing clones. The efficacies of the microsatellite containing genomic sequences were established by superior primer designing ability (69%), PCR amplification efficiency (85.5%) and polymorphic potential (52%) in the parents of F(2) mapping population. Out of 172 pp, functional 147 markers showed high level of cross-species amplification (~74%) in six grass species. Higher polymorphism rate and broad range of genetic diversity (0.30-0.69 averaging 0.58) obtained in constructed phylogenetic tree using 52 microsatellite markers, demonstrated the utility of markers in germplasm characterizations. In silico comparative mapping of 147 foxtail millet microsatellite containing sequences against the mapping data of sorghum (~18%), maize (~16%) and rice (~5%) indicated the presence of orthologous sequences of the foxtail millet in the respective species. The result thus demonstrates the applicability of microsatellite markers in various genotyping applications, determining phylogenetic relationships and comparative mapping in several important grass species.

Genomic diversity in switchgrass (Panicum virgatum): from the continental scale to a dune landscape

Connecting broad-scale patterns of genetic variation and population structure to genetic diversity on a landscape is a key step towards understanding historical processes of migration and adaptation. New genomic approaches can be used to increase the resolution of phylogeographic studies while reducing locus sampling effects and circumventing ascertainment bias. Here, we use a novel approach based on high-throughput sequencing to characterize genetic diversity in complete chloroplast genomes and >10,000 nuclear loci in switchgrass, at continental and landscape scales. Switchgrass is a North American tallgrass species, which is widely used in conservation and perennial biomass production, and shows strong ecotypic adaptation and population structure across the continental range. We sequenced 40.9 billion base pairs from 24 individuals from across the species' range and 20 individuals from the Indiana Dunes. Analysis of plastome sequence revealed 203 variable SNP sites that define eight haplogroups, which are differentiated by 4-127 SNPs and confirmed by patterns of indel variation. These include three deeply divergent haplogroups, which correspond to the previously described lowland-upland ecotypic split and a novel upland haplogroup split that dates to the mid-Pleistocene. Most of the plastome haplogroup diversity present in the northern switchgrass range, including in the Indiana Dunes, originated in the mid- or upper Pleistocene prior to the most recent postglacial recolonization. Furthermore, a recently colonized landscape feature (approximately 150 ya) in the Indiana Dunes contains several deeply divergent upland haplogroups. Nuclear markers also support a deep lowland-upland split, followed by limited gene flow, and show extensive gene flow in the local population of the Indiana Dunes.

Discovery of EST-SSRs in lung cancer: tagged ESTs with SSRs lead to differential amino acid and protein expression patterns in cancerous tissues

Tandem repeats are found in both coding and non-coding sequences of higher organisms. These sequences can be used in cancer genetics and diagnosis to unravel the genetic basis of tumor formation and progression. In this study, a possible relationship between SSR distributions and lung cancer was studied by comparative analysis of EST-SSRs in normal and lung cancerous tissues. While the EST-SSR distribution was similar between tumorous tissues, this distribution was different between normal and tumorous tissues. Trinucleotides tandem repeats were highly different; the number of trinucleotides in ESTs of lung cancer was 3 times higher than normal tissue. Significant negative correlation between normal and cancerous tissue showed that cancerous tissue generates different types of trinucleotides. GGC and CGC were the more frequent expressed trinucleotides in cancerous tissue, but these SSRs were not expressed in normal tissue. Similar to the EST level, the expression pattern of EST-SSRs-derived amino acids was significantly different between normal and cancerous tissues. Arg, Pro, Ser, Gly, and Lys were the most abundant amino acids in cancerous tissues, and Leu, Cys, Phe, and His were significantly more abundant in normal tissues than in cancerous tissues. Next, the putative functions of triplet SSR-containing genes were analyzed. In cancerous tissue, EST-SSRs produce different types of proteins. Chromodomain helicase DNA binding proteins were one of the major protein products of EST-SSRs in the cancerous library, while these proteins were not produced from EST-SSRs in normal tissue. For the first time, the findings of this study confirmed that EST-SSRs in normal lung tissues are different than in unhealthy tissues, and tagged ESTs with SSRs cause remarkable differences in amino acid and protein expression patterns in cancerous tissue. We suggest that EST-SSRs and EST-SSRs differentially expressed in cancerous tissue may be suitable candidate markers for lung cancer diagnosis and prediction.

Characterization and development of EST-derived SSR markers in cultivated sweetpotato (Ipomoea batatas)

BACKGROUND

Currently there exists a limited availability of genetic marker resources in sweetpotato (Ipomoea batatas), which is hindering genetic research in this species. It is necessary to develop more molecular markers for potential use in sweetpotato genetic research. With the newly developed next generation sequencing technology, large amount of transcribed sequences of sweetpotato have been generated and are available for identifying SSR markers by data mining.

RESULTS

In this study, we investigated 181,615 ESTs for the identification and development of SSR markers. In total, 8,294 SSRs were identified from 7,163 SSR-containing unique ESTs. On an average, one SSR was found per 7.1 kb of EST sequence with tri-nucleotide motifs (42.9%) being the most abundant followed by di- (41.2%), tetra- (9.2%), penta- (3.7%) and hexa-nucleotide (3.1%) repeat types. The top five motifs included AG/CT (26.9%), AAG/CTT (13.5%), AT/TA (10.6%), CCG/CGG (5.8%) and AAT/ATT (4.5%). After removing possible duplicate of published EST-SSRs of sweetpotato, a total of non-repeat 7,958 SSR motifs were identified. Based on these SSR-containing sequences, 1,060 pairs of high-quality SSR primers were designed and used for validation of the amplification and assessment of the polymorphism between two parents of one mapping population (E Shu 3 Hao and Guang 2k-30) and eight accessions of cultivated sweetpotatoes. The results showed that 816 primer pairs could yield reproducible and strong amplification products, of which 195 (23.9%) and 342 (41.9%) primer pairs exhibited polymorphism between E Shu 3 Hao and Guang 2k-30 and among the 8 cultivated sweetpotatoes, respectively.

CONCLUSION

This study gives an insight into the frequency, type and distribution of sweetpotato EST-SSRs and demonstrates successful development of EST-SSR markers in cultivated sweetpotato. These EST-SSR markers could enrich the current resource of molecular markers for the sweetpotato community and would be useful for qualitative and quantitative trait mapping, marker-assisted selection, evolution and genetic diversity studies in cultivated sweetpotato and related Ipomoea species.

New genomic resources for switchgrass: a BAC library and comparative analysis of homoeologous genomic regions harboring bioenergy traits

BACKGROUND

Switchgrass, a C4 species and a warm-season grass native to the prairies of North America, has been targeted for development into an herbaceous biomass fuel crop. Genetic improvement of switchgrass feedstock traits through marker-assisted breeding and biotechnology approaches calls for genomic tools development. Establishment of integrated physical and genetic maps for switchgrass will accelerate mapping of value added traits useful to breeding programs and to isolate important target genes using map based cloning. The reported polyploidy series in switchgrass ranges from diploid (2X = 18) to duodecaploid (12X = 108). Like in other large, repeat-rich plant genomes, this genomic complexity will hinder whole genome sequencing efforts. An extensive physical map providing enough information to resolve the homoeologous genomes would provide the necessary framework for accurate assembly of the switchgrass genome.

RESULTS

A switchgrass BAC library constructed by partial digestion of nuclear DNA with EcoRI contains 147,456 clones covering the effective genome approximately 10 times based on a genome size of 3.2 Gigabases (~1.6 Gb effective). Restriction digestion and PFGE analysis of 234 randomly chosen BACs indicated that 95% of the clones contained inserts, ranging from 60 to 180 kb with an average of 120 kb. Comparative sequence analysis of two homoeologous genomic regions harboring orthologs of the rice OsBRI1 locus, a low-copy gene encoding a putative protein kinase and associated with biomass, revealed that orthologous clones from homoeologous chromosomes can be unambiguously distinguished from each other and correctly assembled to respective fingerprint contigs. Thus, the data obtained not only provide genomic resources for further analysis of switchgrass genome, but also improve efforts for an accurate genome sequencing strategy.

CONCLUSIONS

The construction of the first switchgrass BAC library and comparative analysis of homoeologous harboring OsBRI1 orthologs present a glimpse into the switchgrass genome structure and complexity. Data obtained demonstrate the feasibility of using HICF fingerprinting to resolve the homoeologous chromosomes of the two distinct genomes in switchgrass, providing a robust and accurate BAC-based physical platform for this species. The genomic resources and sequence data generated will lay the foundation for deciphering the switchgrass genome and lead the way for an accurate genome sequencing strategy.

New genomic resources for switchgrass: a BAC library and comparative analysis of homoeologous genomic regions harboring bioenergy traits

BACKGROUND

Switchgrass, a C4 species and a warm-season grass native to the prairies of North America, has been targeted for development into an herbaceous biomass fuel crop. Genetic improvement of switchgrass feedstock traits through marker-assisted breeding and biotechnology approaches calls for genomic tools development. Establishment of integrated physical and genetic maps for switchgrass will accelerate mapping of value added traits useful to breeding programs and to isolate important target genes using map based cloning. The reported polyploidy series in switchgrass ranges from diploid (2X = 18) to duodecaploid (12X = 108). Like in other large, repeat-rich plant genomes, this genomic complexity will hinder whole genome sequencing efforts. An extensive physical map providing enough information to resolve the homoeologous genomes would provide the necessary framework for accurate assembly of the switchgrass genome.

RESULTS

A switchgrass BAC library constructed by partial digestion of nuclear DNA with EcoRI contains 147,456 clones covering the effective genome approximately 10 times based on a genome size of 3.2 Gigabases (~1.6 Gb effective). Restriction digestion and PFGE analysis of 234 randomly chosen BACs indicated that 95% of the clones contained inserts, ranging from 60 to 180 kb with an average of 120 kb. Comparative sequence analysis of two homoeologous genomic regions harboring orthologs of the rice OsBRI1 locus, a low-copy gene encoding a putative protein kinase and associated with biomass, revealed that orthologous clones from homoeologous chromosomes can be unambiguously distinguished from each other and correctly assembled to respective fingerprint contigs. Thus, the data obtained not only provide genomic resources for further analysis of switchgrass genome, but also improve efforts for an accurate genome sequencing strategy.

CONCLUSIONS

The construction of the first switchgrass BAC library and comparative analysis of homoeologous harboring OsBRI1 orthologs present a glimpse into the switchgrass genome structure and complexity. Data obtained demonstrate the feasibility of using HICF fingerprinting to resolve the homoeologous chromosomes of the two distinct genomes in switchgrass, providing a robust and accurate BAC-based physical platform for this species. The genomic resources and sequence data generated will lay the foundation for deciphering the switchgrass genome and lead the way for an accurate genome sequencing strategy.