Integration of Cot analysis, DNA cloning, and high-throughput sequencing facilitates genome characterization and gene discovery

Genotype by environment interaction, AMMI, GGE biplot, and mega environment analysis of elite Sorghum bicolor (L.) Moench genotypes in humid lowland areas of Ethiopia

This study aimed to evaluate high-yielding, stable sorghum genotypes and determine the ideal (representative and discriminating) testing environments for genotypes in the humid lowlands of Ethiopia. A total of forty-two sorghum genotypes were used for a field trial conducted in six different environments using a randomized complete block design. Yield stability, Additive main effect, multiplicative interaction (AMMI), and genotype and genotype by environment interaction (GGE) were computed. The AMMI analysis explained 62.85% of the G×E variance. The AMMI1 biplot revealed that (G4; Mok079 and (G16; Ba066) genotypes had higher grain yields. AMMI2 biplot suggested that genotypes (G18; Y0470),(G23;100620), (G29; PML981475), and (G11; ETSC300373-4) show higher sensitivity to environmental changes because of their strong genotype-by-environment interactions. The GGE captured 79.46% of the GGE variance, and the GGE biplot identified genotypes (G4; Mok079), (G10; Sl081) and (G16; Ba066) were the most stable genotypes whereas(G39; ETSC120051-3) was the least stable genotypes. The GGE biplot identified Assosa (AS20) as a suitable environment, whereas PW20 and JM20 were the most discriminating and non-representative environments. The GGE biplot was found to identify three main mega-environments for sorghum growing in the humid lowlands of Ethiopia., both the AMMI and GGE biplots revealed (G4; Mok079) had the highest level of adaptability to all tested environments and was approved by the National Variety Release Committee for release in 2022.

Construction of an experimental study and addition of adapter sequences using HiDi DNA polymerase for improving DNA normalization methods relevant to novel gene discovery

Microorganisms in the environment can be distinguished into dominant and rare microbial species based on their genes. It is difficult to obtain genetic information derived from rare microbial species (rare genes) because of the differences in relative abundance. DNA normalization is an approach that is used to obtain genetic information derived from rare microbial species from an environmental sample. This method involves the addition of adapter sequences for the amplification, denaturation, and reassociation of the DNA fragments and single-stranded DNA (ssDNA)/double-stranded DNA (dsDNA) separation. In this method, the amount of a high-copy-number of DNA fragments and a low-copy-number of DNA fragments can be equalized. Improvements in this technique are expected to provide novel genetic information or genes in rare microbial species. However, few model experimental systems have been reported to validate the DNA normalization techniques. This study is aimed to improve the DNA normalization technique used to obtain genetic information of rare genes from rare microbial species. An experimental study was constructed with two antibiotic resistance genes, whose copy numbers differed up to a million-fold. Both genes were mixed and the mixture of DNA fragments, of high- and low-copy-number, containing these genes was normalized by separating ssDNA/dsDNA fragments using hydroxyapatite. Normalized DNA fragments were introduced into Escherichia coli and DNA normalization was evaluated by counting colonies. Moreover, we improved the method to amplify a low-copy-number of DNA fragments by the addition of adapter sequences to DNA fragments using HiDi DNA polymerase to increase the efficiency of DNA normalization. This normalization method was achieved with a 100,000-fold difference. These methods allowed for quantitative evaluation of the DNA normalization efficiency. The experimental data and methods obtained in this study are expected to improve the DNA normalization efficiency to obtain novel genetic information or genes.

Ploidy Identification by Flow Cytometry and Application of the Method to Characterize Seasonal Ploidy Variation of Wild Populations of the Red Alga Gracilariopsis lemaneiformis

Gracilariopsis lemaneiformis (Gp. lemaneiformis) is an economically important alga. At present, there is no way to quickly and easily determine its ploidy and life cycle dynamics in wild populations, which affects the process of genetic breeding. In this study, we developed and verified a ploidy identification method using flow cytometry and then used it to explore the seasonal fluctuation of the ploidy ratio and the environmental factors that influence it in wild populations of this species. Of the three methods we tested for nucleus extraction, quick chopping was the best because of its high extraction efficiency, low debris background, obvious subcellular scatter plot, and clear typical histogram. Samples from the tip of the alga were more suitable for preparing the nuclear suspension than samples from the base. Generalized linear model analysis based on diagnosis of multicollinearity revealed a negative correlation between temperature and ploidy ratio. Among the environmental factors tested, temperature had the greatest influence on the ploidy ratio, whereas precipitation and sunshine duration had no effect on the ploidy ratio fluctuation. Our study will be useful for material collection and studies of utilization and life cycle dynamics. Moreover, understanding of ploidy dynamics may provide a theoretical basis for improving variety and breeding of Gp. lemaneiformis in the future.

Molecular biology applications of the red king crab duplex-specific nuclease

Duplex-specific nuclease (DSN) from hepatopancreas of the craboid Paralithodes camtschaticus (red king crab) has a unique combination of properties. Along with thermal stability and a high optimal temperature of catalysis, this enzyme exhibits high substrate selectivity, cleaving only DNA in duplexes (DNA-DNA or DNA-RNA). Accordingly, it digests neither single strands (nor single-stranded regions) of DNA, nor RNA strands with any secondary structure. Such properties make it possible to create unique protocols based on DSN, which is also an important object of fundamental research in the field of nuclease evolution. The review considers diverse applications of the red king crab DSN in modern methods of molecular biology.

Genotype by environment interaction, correlation, AMMI, GGE biplot and cluster analysis for grain yield and other agronomic traits in sorghum (Sorghum bicolor L. Moench)

Genotype by environment (G×E) interaction is a major factor limiting the success of germplasm selection and identification of superior genotypes for use in plant breeding programs. Similar to the case in other crops, G×E complicates the improvement of sorghum, and hence it should be determined and used in decision-making programs. The present study aimed at assessing the G×E interaction, and the correlation between traits for superior sorghum genotypes. Three hundred twenty sorghum landraces and four improved varieties were used in alpha lattice experimental design-based field trial across three environments (Melkassa, Mieso and Mehoni) in Ethiopia. Phenotypic data were collected for days to flowering (DTF), plant height (PH), panicle length (PALH), panicle width (PAWD), panicle weight (PAWT) and grain yield (GY). The results revealed that the variance due to genotype, environment and G×E interaction were highly significant (P < 0.001) for all traits. GY and PAWT were highly affected by environments and G×E whereas DTF, PALH, PAWD and PH were mainly affected by genotypic variation. Therefore, multi-environment testing is needed for taking care of G × E interaction to identify high yielding and stable sorghum landraces. GY and PAWT revealed highly significant positive correlations indicating the possibility of effective selection of the two traits simultaneously. Among the studied populations, South Wello, West Hararghe and Shewa zones had highly diverse genotypes that were distributed across all clusters. Hence, these areas can be considered as hotspots for identifying divergent sorghum landraces that could be used in breeding programs. Melkassa was the most representative environment whereas Mieso was the most discriminating. Five genotypes (G148, G123, G110, G203 and G73) were identified as superior across the test environments for grain yield with farmer-preferred trait, such as plant height. The identified stable and high yielding genotypes are valuable genetic resources that should be used in sorghum breeding programs.

Assessment of genetic diversity and variety identification based on developed retrotransposon-based insertion polymorphism (RBIP) markers in sweet potato (Ipomoea batatas (L.) Lam.)

Sweet potato, a dicotyledonous and perennial plant, is the third tuber/root crop species behind potato and cassava in terms of production. Long terminal repeat (LTR) retrotransposons are highly abundant in sweet potato, contributing to genetic diversity. Retrotransposon-based insertion polymorphism (RBIP) is a high-throughput marker system to study the genetic diversity of plant species. To date, there have been no transposon marker-based genetic diversity analyses of sweet potato. Here, we reported a structure-based analysis of the sweet potato genome, a total of 21555 LTR retrotransposons, which belonged to the main LTR-retrotransposon subfamilies Ty3-gypsy and Ty1-copia were identified. After searching and selecting using Hidden Markov Models (HMMs), 1616 LTR retrotransposon sequences containing at least two models were screened. A total of 48 RBIP primers were synthesized based on the high copy numbers of conserved LTR sequences. Fifty-six amplicons with an average polymorphism of 91.07% were generated in 105 sweet potato germplasm resources based on RBIP markers. A Unweighted Pair Group Method with Arithmatic Mean (UPGMA) dendrogram, a model-based genetic structure and principal component analysis divided the sweet potato germplasms into 3 groups containing 8, 53, and 44 germplasms. All the three analyses produced significant groupwise consensus. However, almost all the germplasms contained only one primary locus. The analysis of molecular variance (AMOVA) among the groups indicated higher intergroup genetic variation (53%) than intrapopulation genetic variation. In addition, long-term self-retention may cause some germplasm resources to exhibit variable segregation. These results suggest that these sweet potato germplasms are not well evolutionarily diversified, although geographic speciation could have occurred at a limited level. This study highlights the utility of RBIP markers for determining the intraspecies variability of sweet potato and have the potential to be used as core primer pairs for variety identification, genetic diversity assessment and linkage map construction. The results could provide a good theoretical reference and guidance for germplasm research and breeding.

Cytogenetic analysis of Bienertia sinuspersici Akhani as the first step in genome sequencing

BACKGROUND

C₄ plants are efficient in suppressing photorespiration and enhancing carbon gain as compared to C₃ plants. Bienertia sinuspersici Akhani is one of the few species in the family Amaranthaceae that can perform C₄ photosynthesis within individual chlorenchyma cells, without the conventional Kranz anatomy in its leaf. This plant is salt-tolerant and is well-adapted to thrive in hot and humid climates. To date, there have been no reported cytogenetic analyses yet on this species.

OBJECTIVE

This study aims to provide a cytogenetic analysis of B. sinuspersici as the first step in genome sequencing.

METHODS

Fluorescence in situ hybridization (FISH) karyotype analysis was conducted using the metaphase chromosomes of B. sinuspersici probed with 5S and 45S rDNA and Arabidopsis-type telomeric repeats.

RESULTS

Results of the cytogenetic analysis confirmed that B. sinuspersici has 2n = 2x = 18 consisting of nine pairs of metacentric chromosomes. Two loci of 45S rDNA were found on the distal regions of the short arm of chromosome 7. Nine loci of 5S rDNA were found in the pericentromeric regions of chromosomes 1, 3, 4, 6, and 8, which also colocalized with Arabidopsis-type telomeric repeats; while four loci in the interstitial regions of chromosome 5 and 8 can be observed. The single locus of 5S rDNA that was found in chromosome 8 appears to be hemizygous.

CONCLUSION

The FISH karyotype analysis, based on the combination of rDNAs, telomeric tandem repeat markers and C₀t DNA chromosome landmarks, allowed efficient chromosome identification and provided useful information in characterizing the genome of B. sinuspersici.

Transposable Elements: Classification, Identification, and Their Use As a Tool For Comparative Genomics

Most genomes are populated by hundreds of thousands of sequences originated from mobile elements. On the one hand, these sequences present a real challenge in the process of genome analysis and annotation. On the other hand, they are very interesting biological subjects involved in many cellular processes. Here we present an overview of transposable elements biodiversity, and we discuss different approaches to transposable elements detection and analyses.

Exploring the loblolly pine (Pinus taeda L.) genome by BAC sequencing and Cot analysis

Loblolly pine (LP; Pinus taeda L.) is an economically and ecologically important tree in the southeastern U.S. To advance understanding of the loblolly pine (LP; Pinus taeda L.) genome, we sequenced and analyzed 100 BAC clones and performed a Cot analysis. The Cot analysis indicates that the genome is composed of 57, 24, and 10% highly-repetitive, moderately-repetitive, and single/low-copy sequences, respectively (the remaining 9% of the genome is a combination of fold back and damaged DNA). Although single/low-copy DNA only accounts for 10% of the LP genome, the amount of single/low-copy DNA in LP is still 14 times the size of the Arabidopsis genome. Since gene numbers in LP are similar to those in Arabidopsis, much of the single/low-copy DNA of LP would appear to be composed of DNA that is both gene- and repeat-poor. Macroarrays prepared from a LP bacterial artificial chromosome (BAC) library were hybridized with probes designed from cell wall synthesis/wood development cDNAs, and 50 of the "targeted" clones were selected for further analysis. An additional 25 clones were selected because they contained few repeats, while 25 more clones were selected at random. The 100 BAC clones were Sanger sequenced and assembled. Of the targeted BACs, 80% contained all or part of the cDNA used to target them. One targeted BAC was found to contain fungal DNA and was eliminated from further analysis. Combinations of similarity-based and ab initio gene prediction approaches were utilized to identify and characterize potential coding regions in the 99 BACs containing LP DNA. From this analysis, we identified 154 gene models (GMs) representing both putative protein-coding genes and likely pseudogenes. Ten of the GMs (all of which were specifically targeted) had enough support to be classified as intact genes. Interestingly, the 154 GMs had statistically indistinguishable (α = 0.05) distributions in the targeted and random BAC clones (15.18 and 12.61 GM/Mb, respectively), whereas the low-repeat BACs contained significantly fewer GMs (7.08 GM/Mb). However, when GM length was considered, the targeted BACs had a significantly greater percentage of their length in GMs (3.26%) when compared to random (1.63%) and low-repeat (0.62%) BACs. The results of our study provide insight into LP evolution and inform ongoing efforts to produce a reference genome sequence for LP, while characterization of genes involved in cell wall production highlights carbon metabolism pathways that can be leveraged for increasing wood production.

Capturing the 'ome': the expanding molecular toolbox for RNA and DNA library construction

All sequencing experiments and most functional genomics screens rely on the generation of libraries to comprehensively capture pools of targeted sequences. In the past decade especially, driven by the progress in the field of massively parallel sequencing, numerous studies have comprehensively assessed the impact of particular manipulations on library complexity and quality, and characterized the activities and specificities of several key enzymes used in library construction. Fortunately, careful protocol design and reagent choice can substantially mitigate many of these biases, and enable reliable representation of sequences in libraries. This review aims to guide the reader through the vast expanse of literature on the subject to promote informed library generation, independent of the application.

Complete chloroplast genome sequence and comparative analysis of loblolly pine (Pinus taeda L.) with related species

Pinaceae, the largest family of conifers, has a diversified organization of chloroplast (cp) genomes with two typical highly reduced inverted repeats (IRs). In the current study, we determined the complete sequence of the cp genome of an economically and ecologically important conifer tree, the loblolly pine (Pinus taeda L.), using Illumina paired-end sequencing and compared the sequence with those of other pine species. The results revealed a genome size of 121,531 base pairs (bp) containing a pair of 830-bp IR regions, distinguished by a small single copy (42,258 bp) and large single copy (77,614 bp) region. The chloroplast genome of P. taeda encodes 120 genes, comprising 81 protein-coding genes, four ribosomal RNA genes, and 35 tRNA genes, with 151 randomly distributed microsatellites. Approximately 6 palindromic, 34 forward, and 22 tandem repeats were found in the P. taeda cp genome. Whole cp genome comparison with those of other Pinus species exhibited an overall high degree of sequence similarity, with some divergence in intergenic spacers. Higher and lower numbers of indels and single-nucleotide polymorphism substitutions were observed relative to P. contorta and P. monophylla, respectively. Phylogenomic analyses based on the complete genome sequence revealed that 60 shared genes generated trees with the same topologies, and P. taeda was closely related to P. contorta in the subgenus Pinus. Thus, the complete P. taeda genome provided valuable resources for population and evolutionary studies of gymnosperms and can be used to identify related species.

Whole-Genome Sequence Accuracy Is Improved by Replication in a Population of Mutagenized Sorghum

The accurate detection of induced mutations is critical for both forward and reverse genetics studies. Experimental chemical mutagenesis induces relatively few single base changes per individual. In a complex eukaryotic genome, false positive detection of mutations can occur at or above this mutagenesis rate. We demonstrate here, using a population of ethyl methanesulfonate (EMS)-treated Sorghum bicolor BTx623 individuals, that using replication to detect false positive-induced variants in next-generation sequencing (NGS) data permits higher throughput variant detection with greater accuracy. We used a lower sequence coverage depth (average of 7×) from 586 independently mutagenized individuals and detected 5,399,493 homozygous single nucleotide polymorphisms (SNPs). Of these, 76% originated from only 57,872 genomic positions prone to false positive variant calling. These positions are characterized by high copy number paralogs where the error-prone SNP positions are at copies containing a variant at the SNP position. The ability of short stretches of homology to generate these error-prone positions suggests that incompletely assembled or poorly mapped repeated sequences are one driver of these error-prone positions. Removal of these false positives left 1,275,872 homozygous and 477,531 heterozygous EMS-induced SNPs, which, congruent with the mutagenic mechanism of EMS, were >98% G:C to A:T transitions. Through this analysis, we generated a collection of sequence indexed mutants of sorghum. This collection contains 4035 high-impact homozygous mutations in 3637 genes and 56,514 homozygous missense mutations in 23,227 genes. Each line contains, on average, 2177 annotated homozygous SNPs per genome, including seven likely gene knockouts and 96 missense mutations. The number of mutations in a transcript was linearly correlated with the transcript length and also the G+C count, but not with the GC/AT ratio. Analysis of the detected mutagenized positions identified CG-rich patches, and flanking sequences strongly influenced EMS-induced mutation rates. This method for detecting false positive-induced mutations is generally applicable to any organism, is independent of the choice of in silico variant-calling algorithm, and is most valuable when the true mutation rate is likely to be low, such as in laboratory-induced mutations or somatic mutation detection in medicine.

Reproducible genomic DNA preparation from diverse crop species for molecular genetic applications

BACKGROUND

Several high-throughput molecular genetic analyses rely on high-quality genomic DNA. Copurification of other molecules can negatively impact the functionality of plant DNA preparations employed in these procedures. Isolating DNA from agronomically important crops, such as sugarcane, rice, citrus, potato and tomato is a challenge due to the presence of high fiber, polysaccharides, or secondary metabolites. We present a simplified, rapid and reproducible SDS-based method that provides high-quality and -quantity of DNA from small amounts of leaf tissue, as required by the emerging biotechnology and molecular genetic applications.

RESULTS

We developed the TENS-CO method as a simplified SDS-based isolation procedure with sequential steps of purification to remove polysaccharides and polyphenols using 2-mercaptoethanol and potassium acetate, chloroform partitioning, and sodium acetate/ethanol precipitation to yield high-quantity and -quality DNA consistently from small amounts of tissue (0.15 g) for different plant species. The method is simplified and rapid in terms of requiring minimal manipulation, smaller extraction volume, reduced homogenization time (20 s) and DNA precipitation (one precipitation for 1 h). The method has been demonstrated to accelerate screening of large amounts of plant tissues from species that are rich in polysaccharides and secondary metabolites for Southern blot analysis of reporter gene overexpressing lines, pathogen detection by quantitative PCR, and genotyping of disease-resistant plants using marker-assisted selection.

CONCLUSION

To facilitate molecular genetic studies in major agronomical crops, we have developed the TENS-CO method as a simple, rapid, reproducible and scalable protocol enabling efficient and robust isolation of high-quality and -quantity DNA from small amounts of tissue from sugarcane, rice, citrus, potato, and tomato, thereby reducing significantly the time and resources used for DNA isolation.

Satellite DNA: An Evolving Topic

Satellite DNA represents one of the most fascinating parts of the repetitive fraction of the eukaryotic genome. Since the discovery of highly repetitive tandem DNA in the 1960s, a lot of literature has extensively covered various topics related to the structure, organization, function, and evolution of such sequences. Today, with the advent of genomic tools, the study of satellite DNA has regained a great interest. Thus, Next-Generation Sequencing (NGS), together with high-throughput in silico analysis of the information contained in NGS reads, has revolutionized the analysis of the repetitive fraction of the eukaryotic genomes. The whole of the historical and current approaches to the topic gives us a broad view of the function and evolution of satellite DNA and its role in chromosomal evolution. Currently, we have extensive information on the molecular, chromosomal, biological, and population factors that affect the evolutionary fate of satellite DNA, knowledge that gives rise to a series of hypotheses that get on well with each other about the origin, spreading, and evolution of satellite DNA. In this paper, I review these hypotheses from a methodological, conceptual, and historical perspective and frame them in the context of chromosomal organization and evolution.

Rapid Gene Isolation Using MutChromSeq

MutChromSeq is an approach for isolation of genes and DNA sequences controlling gene expression in plants with complex and polyploid genomes. It involves a lossless complexity reduction by flow cytometric chromosome sorting and shotgun sequencing DNA from isolated chromosomes. Comparison of sequences from wild-type parental chromosome with chromosomes from multiple independently derived mutants identifies causative mutations in a single candidate gene or a noncoding sequence. MutChromSeq does not rely on recombination-based genetic mapping and does not exclude any DNA sequence from being targeted.

Generation and Characterisation of a Reference Transcriptome for Lentil (Lens culinaris Medik.)

RNA-Seq using second-generation sequencing technologies permits generation of a reference unigene set for a given species, in the absence of a well-annotated genome sequence, supporting functional genomics studies, gene characterisation and detailed expression analysis for specific morphophysiological or environmental stress response traits. A reference unigene set for lentil has been developed, consisting of 58,986 contigs and scaffolds with an N50 length of 1719 bp. Comparison to gene complements from related species, reference protein databases, previously published lentil transcriptomes and a draft genome sequence validated the current dataset in terms of degree of completeness and utility. A large proportion (98%) of unigenes were expressed in more than one tissue, at varying levels. Candidate genes associated with mechanisms of tolerance to both boron toxicity and time of flowering were identified, which can eventually be used for the development of gene-based markers. This study has provided a comprehensive, assembled and annotated reference gene set for lentil that can be used for multiple applications, permitting identification of genes for pathway-specific expression analysis, genetic modification approaches, development of resources for genotypic analysis, and assistance in the annotation of a future lentil genome sequence.

Methods for accurate quantification of LTR-retrotransposon copy number using short-read sequence data: a case study in Sorghum

Transposable elements (TEs) are ubiquitous in eukaryotic genomes and their mobility impacts genome structure and function in myriad ways. Because of their abundance, activity, and repetitive nature, the characterization and analysis of TEs remain challenging, particularly from short-read sequencing projects. To overcome this difficulty, we have developed a method that estimates TE copy number from short-read sequences. To test the accuracy of our method, we first performed an in silico analysis of the reference Sorghum bicolor genome, using both reference-based and de novo approaches. The resulting TE copy number estimates were strikingly similar to the annotated numbers. We then tested our method on real short-read data by estimating TE copy numbers in several accessions of S. bicolor and its close relative S. propinquum. Both methods effectively identify and rank similar TE families from highest to lowest abundance. We found that de novo characterization was effective at capturing qualitative variation, but underestimated the abundance of some TE families, specifically families of more ancient origin. Also, interspecific reference-based mapping of S. propinquum reads to the S. bicolor database failed to fully describe TE content in S. propinquum, indicative of recent TE activity leading to changes in the respective repetitive landscapes over very short evolutionary timescales. We conclude that reference-based analyses are best suited for within-species comparisons, while de novo approaches are more reliable for evolutionarily distant comparisons.

Ploidy Level and DNA Content of Erianthus arundinaceus as Determined by Flow Cytometry and the Association with Biological Characteristics

Erianthus arundinaceus is not only an important germplasm resource for sugarcane breeding but also a potential bioenergy plant. Making clear the distribution of the chromosome ploidy of wild E. arundinaceus in china is the premise of the research and utilization of this species. Therefore, the objectives of this study were to determine the ploidy level and DNA content of the 55 E. arundinaceus accessions using flow cytometry and to identify the correlation between ploidy and phenotypic traits. Among the 55 accessions, four tetraploids and 51 hexaploids were identified. The four tetraploids originated from Mengma Yunnan, Shuangjiang Yunnan, Gaozhou Guangdong and Chengle Sichuan. The mean DNA content was 4.82 pg/2C for the tetraploid and 7.30 pg/2C for the hexaploid plants. The ploidy was negatively correlated with cellulose content and positively correlated (P<0.05) with plant height, stem diameter, leaf width, dry weight per plant, fresh weight per plant and hemicellulose content. However, ploidy was not correlated with leaf length, tiller number and the ratio of dry weight and fresh weight. This study will be useful for revealing the distribution of the ploidy of wild E. arundinaceus in Chin, traits markers analysis, and utilization of this species, such as cultivar improvement and sugarcane breeding in the future.

Primitive Genepools of Asian Pears and Their Complex Hybrid Origins Inferred from Fluorescent Sequence-Specific Amplification Polymorphism (SSAP) Markers Based on LTR Retrotransposons

Recent evidence indicated that interspecific hybridization was the major mode of evolution in Pyrus. The genetic relationships and origins of the Asian pear are still unclear because of frequent hybrid events, fast radial evolution, and lack of informative data. Here, we developed fluorescent sequence-specific amplification polymorphism (SSAP) markers with lots of informative sites and high polymorphism to analyze the population structure among 93 pear accessions, including nearly all species native to Asia. Results of a population structure analysis indicated that nearly all Asian pear species experienced hybridization, and originated from five primitive genepools. Four genepools corresponded to four primary Asian species: P. betulaefolia, P. pashia, P. pyrifolia, and P. ussuriensis. However, cultivars of P. ussuriensis were not monophyletic and introgression occurred from P. pyrifolia. The specific genepool detected in putative hybrids between occidental and oriental pears might be from occidental pears. The remaining species, including P. calleryana, P. xerophila, P. sinkiangensis, P. phaeocarpa, P. hondoensis, and P. hopeiensis in Asia, were inferred to be of hybrid origins and their possible genepools were identified. This study will be of great help for understanding the origin and evolution of Asian pears.

Isolation and characterization of putative functional long terminal repeat retrotransposons in the Pyrus genome

Background

Long terminal repeat (LTR)-retrotransposons constitute 42.4 % of the genome of the ‘Suli’ pear (Pyrus pyrifolia white pear group), implying that retrotransposons have played important roles in Pyrus evolution. Therefore, further analysis of retrotransposons will enhance our understanding of the evolutionary history of Pyrus.

Results

We identified 1836 LTR-retrotransposons in the ‘Suli’ pear genome, of which 440 LTR-retrotransposons were predicted to contain at least two of three gene models (gag, integrase and reverse transcriptase). Because these were most likely to be functional transposons, we focused our analyses on this set of 440. Most of the LTR-retrotransposons were estimated to have inserted into the genome less than 2.5 million years ago. Sequence analysis showed that the reverse transcriptase component of the identified LTR-retrotransposons was highly heterogeneous. Analyses of transcripts assembled from RNA-Seq databases of two cultivars of Pyrus species showed that LTR-retrotransposons were expressed in the buds and fruit of Pyrus. A total of 734 coding sequences in the ‘Suli’ genome were disrupted by the identified LTR-retrotransposons. Five high-copy-number LTR-retrotransposon families were identified in Pyrus. These families were rarely found in the genomes of Malus and Prunus, but were distributed extensively in Pyrus and abundance varied between species.

Conclusions

We identified potentially functional, full-length LTR-retrotransposons with three gene models in the ‘Suli’ genome. The analysis of RNA-seq data demonstrated that these retrotransposons are expressed in the organs of pears. The differential copy number of LTR-retrotransposon families between Pyrus species suggests that the transposition of retrotransposons is an important evolutionary force driving the genetic divergence of species within the genus.

Electronic supplementary material

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

Epigenomics and bolting tolerance in sugar beet genotypes

In sugar beet (Beta vulgaris altissima), bolting tolerance is an essential agronomic trait reflecting the bolting response of genotypes after vernalization. Genes involved in induction of sugar beet bolting have now been identified, and evidence suggests that epigenetic factors are involved in their control. Indeed, the time course and amplitude of DNA methylation variations in the shoot apical meristem have been shown to be critical in inducing sugar beet bolting, and a few functional targets of DNA methylation during vernalization have been identified. However, molecular mechanisms controlling bolting tolerance levels among genotypes are still poorly understood. Here, gene expression and DNA methylation profiles were compared in shoot apical meristems of three bolting-resistant and three bolting-sensitive genotypes after vernalization. Using Cot fractionation followed by 454 sequencing of the isolated low-copy DNA, 6231 contigs were obtained that were used along with public sugar beet DNA sequences to design custom Agilent microarrays for expression (56k) and methylation (244k) analyses. A total of 169 differentially expressed genes and 111 differentially methylated regions were identified between resistant and sensitive vernalized genotypes. Fourteen sequences were both differentially expressed and differentially methylated, with a negative correlation between their methylation and expression levels. Genes involved in cold perception, phytohormone signalling, and flowering induction were over-represented and collectively represent an integrative gene network from environmental perception to bolting induction. Altogether, the data suggest that the genotype-dependent control of DNA methylation and expression of an integrative gene network participate in bolting tolerance in sugar beet, opening up perspectives for crop improvement.

Sequencing and assembly of the 22-gb loblolly pine genome

Conifers are the predominant gymnosperm. The size and complexity of their genomes has presented formidable technical challenges for whole-genome shotgun sequencing and assembly. We employed novel strategies that allowed us to determine the loblolly pine (Pinus taeda) reference genome sequence, the largest genome assembled to date. Most of the sequence data were derived from whole-genome shotgun sequencing of a single megagametophyte, the haploid tissue of a single pine seed. Although that constrained the quantity of available DNA, the resulting haploid sequence data were well-suited for assembly. The haploid sequence was augmented with multiple linking long-fragment mate pair libraries from the parental diploid DNA. For the longest fragments, we used novel fosmid DiTag libraries. Sequences from the linking libraries that did not match the megagametophyte were identified and removed. Assembly of the sequence data were aided by condensing the enormous number of paired-end reads into a much smaller set of longer “super-reads,” rendering subsequent assembly with an overlap-based assembly algorithm computationally feasible. To further improve the contiguity and biological utility of the genome sequence, additional scaffolding methods utilizing independent genome and transcriptome assemblies were implemented. The combination of these strategies resulted in a draft genome sequence of 20.15 billion bases, with an N50 scaffold size of 66.9 kbp.

Prediction of retrotransposons and assessment of genetic variability based on developed retrotransposon-based insertion polymorphism (RBIP) markers in Pyrus L

Interspecific hybridization has been considered the major mode of evolution in Pyrus (pear), and thus, the genetic relationships within this genus have not been well documented. Retrotransposons are ubiquitous components of plant genomes and 42.4 % of the pear genome was reported to be long terminal repeat (LTR) retrotransposons, implying that retrotransposons might be significant in the evolution of Pyrus. In this study, 1,836 putative full-length LTR retrotransposons were isolated and 196 retrotransposon-based insertion polymorphism (RBIP) primers were developed, of which 24 pairs to the Ppcr1 subfamily of copia retrotransposons were used to analyze genetic diversity among 110 Pyrus accessions from Eurasia. Our results showed that Ppcr1 replicated many times in the development of cultivated Asian pears. The genetic structure analysis and the unweighted pair group method with arithmetic mean (UPGMA) dendrogram indicated that all accessions could be divided into Oriental and Occidental groups. In Oriental pears, wild pea pears clustered separately into independent groups in accordance with their morphological classifications. Cultivars of P. ussuriensis Maxim, P. pyrifolia Nakai, and P. pyrifolia Chinese white pear were mingled together, which inferred that hybridization events occurred during the development of the cultivated Asian pears. In Occidental pears, two clades were obtained in the UPGMA dendrogram in accordance with their geographical distribution; one contained the European species and the other included species from North Africa and West Asia. New findings in this study will be important to further understand the phylogeny of Pyrus and origins of cultivated pears.

Chromosomal organizations of major repeat families on potato (Solanum tuberosum) and further exploring in its sequenced genome

One of the most powerful technologies in unraveling the organization of a eukaryotic plant genome is high-resolution Fluorescent in situ hybridization of repeats and single copy DNA sequences on pachytene chromosomes. This technology allows the integration of physical mapping information with chromosomal positions, including centromeres, telomeres, nucleolar-organizing region, and euchromatin and heterochromatin. In this report, we established chromosomal positions of different repeat fractions of the potato genomic DNA (Cot100, Cot500 and Cot1000) on the chromosomes. We also analysed various repeat elements that are unique to potato including the moderately repetitive P5 and REP2 elements, where the REP2 is part of a larger Gypsy-type LTR retrotransposon and cover most chromosome regions, with some brighter fluorescing spots in the heterochromatin. The most abundant tandem repeat is the potato genomic repeat 1 that covers subtelomeric regions of most chromosome arms. Extensive multiple alignments of these repetitive sequences in the assembled RH89-039-16 potato BACs and the draft assembly of the DM1-3 516 R44 genome shed light on the conservation of these repeats within the potato genome. The consensus sequences thus obtained revealed the native complete transposable elements from which they were derived.

Defining loci in restriction-based reduced representation genomic data from nonmodel species: sources of bias and diagnostics for optimal clustering

Next generation sequencing holds great promise for applications of phylogeography, landscape genetics, and population genomics in wild populations of nonmodel species, but the robustness of inferences hinges on careful experimental design and effective bioinformatic removal of predictable artifacts. Addressing this issue, we use published genomes from a tunicate, stickleback, and soybean to illustrate the potential for bioinformatic artifacts and introduce a protocol to minimize two sources of error expected from similarity-based de-novo clustering of stacked reads: the splitting of alleles into different clusters, which creates false homozygosity, and the grouping of paralogs into the same cluster, which creates false heterozygosity. We present an empirical application focused on Ciona savignyi, a tunicate with very high SNP heterozygosity (~0.05), because high diversity challenges the computational efficiency of most existing nonmodel pipelines while also potentially exacerbating paralog artifacts. The simulated and empirical data illustrate the advantages of using higher sequence difference clustering thresholds than is typical and demonstrate the utility of our protocol for efficiently identifying an optimum threshold from data without prior knowledge of heterozygosity. The empirical Ciona savignyi data also highlight null alleles as a potentially large source of false homozygosity in restriction-based reduced representation genomic data.

Fluorescent In Situ Hybridization (FISH) on pachytene chromosomes as a tool for genome characterization

A growing number of international genome consortia have initiated large-scale sequencing projects for most of the major crop species. This huge amount of information not only boosted genetic and physical mapping research, but it also enabled novel applications on the level of chromosome biology including molecular cytogenetics that supports plant genetics, genomics, and breeding programs. The simultaneous detection of a large number of BAC-based probes by multicolor fluorescent in situ hybridization (FISH) can provide a rapid overview of super-contig and gap distribution on euchromatin chromosome areas and will display directly and precisely the positions of chromosome rearrangements. Furthermore, hybridizations of BACs on the chromosomes of related species can confirm genomic colinearity, or the occurrence of inversions and translocations events. This cross-species FISH together with meiotic pairing studies is a powerful source of information that elucidates the nature of genome rearrangements, and the consequences of such rearrangements for introgressive hybridizations. In this chapter we describe a general-purpose protocol for FISH on pachytene chromosomes.

Whole-exome targeted sequencing of the uncharacterized pine genome

The large genome size of many species hinders the development and application of genomic tools to study them. For instance, loblolly pine (Pinus taeda L.), an ecologically and economically important conifer, has a large and yet uncharacterized genome of 21.7 Gbp. To characterize the pine genome, we performed exome capture and sequencing of 14 729 genes derived from an assembly of expressed sequence tags. Efficiency of sequence capture was evaluated and shown to be similar across samples with increasing levels of complexity, including haploid cDNA, haploid genomic DNA and diploid genomic DNA. However, this efficiency was severely reduced for probes that overlapped multiple exons, presumably because intron sequences hindered probe:exon hybridizations. Such regions could not be entirely avoided during probe design, because of the lack of a reference sequence. To improve the throughput and reduce the cost of sequence capture, a method to multiplex the analysis of up to eight samples was developed. Sequence data showed that multiplexed capture was reproducible among 24 haploid samples, and can be applied for high-throughput analysis of targeted genes in large populations. Captured sequences were de novo assembled, resulting in 11 396 expanded and annotated gene models, significantly improving the knowledge about the pine gene space. Interspecific capture was also evaluated with over 98% of all probes designed from P. taeda that were efficient in sequence capture, were also suitable for analysis of the related species Pinus elliottii Engelm.

cDNA normalization by hydroxyapatite chromatography to enrich transcriptome diversity in RNA-seq applications

Second-generation sequencing (SGS) has become the preferred method for RNA transcriptome profiling of organisms and single cells. However, SGS analysis of transcriptome diversity (including protein-coding transcripts and regulatory non-coding RNAs) is inefficient unless the sample of interest is first depleted of nucleic acids derived from ribosomal RNA (rRNA), which typically account for up to 95% of total intracellular RNA content. Here we describe a novel microscale hydroxyapatite chromatography (HAC) normalization method to remove eukaryotic and prokaryotic high abundant rRNA species, thereby increasing sequence coverage depth and transcript diversity across non-rRNA populations. RNA-seq analysis of Escherichia coli K-12 and human intracellular total RNA showed that HAC-based normalization enriched for all non-ribosomal RNA species regardless of RNA transcript abundance or length when compared with untreated controls. Microcolumn HAC normalization generated rRNA-depleted cDNA libraries comparable to the well-established duplex specific nuclease (DSN) normalization and Ribo-Zero rRNA-depletion methods, thus establishing microscale HAC as an effective, cost saving, and non-destructive alternative normalization technique.

Consequences of normalizing transcriptomic and genomic libraries of plant genomes using a duplex-specific nuclease and tetramethylammonium chloride

Several applications of high throughput genome and transcriptome sequencing would benefit from a reduction of the high-copy-number sequences in the libraries being sequenced and analyzed, particularly when applied to species with large genomes. We adapted and analyzed the consequences of a method that utilizes a thermostable duplex-specific nuclease for reducing the high-copy components in transcriptomic and genomic libraries prior to sequencing. This reduces the time, cost, and computational effort of obtaining informative transcriptomic and genomic sequence data for both fully sequenced and non-sequenced genomes. It also reduces contamination from organellar DNA in preparations of nuclear DNA. Hybridization in the presence of 3 M tetramethylammonium chloride (TMAC), which equalizes the rates of hybridization of GC and AT nucleotide pairs, reduced the bias against sequences with high GC content. Consequences of this method on the reduction of high-copy and enrichment of low-copy sequences are reported for Arabidopsis and lettuce.

The dynamics of LTR retrotransposon accumulation across 25 million years of panicoid grass evolution

Sample sequence analysis was employed to investigate the repetitive DNAs that were most responsible for the evolved variation in genome content across seven panicoid grasses with >5-fold variation in genome size and different histories of polyploidy. In all cases, the most abundant repeats were LTR retrotransposons, but the particular families that had become dominant were found to be different in the Pennisetum, Saccharum, Sorghum and Zea lineages. One element family, Huck, has been very active in all of the studied species over the last few million years. This suggests the transmittal of an active or quiescent autonomous set of Huck elements to this lineage at the founding of the panicoids. Similarly, independent recent activity of Ji and Opie elements in Zea and of Leviathan elements in Sorghum and Saccharum species suggests that members of these families with exceptional activation potential were present in the genome(s) of the founders of these lineages. In a detailed analysis of the Zea lineage, the combined action of several families of LTR retrotransposons were observed to have approximately doubled the genome size of Zea luxurians relative to Zea mays and Zea diploperennis in just the last few million years. One of the LTR retrotransposon amplification bursts in Zea may have been initiated by polyploidy, but the great majority of transposable element activations are not. Instead, the results suggest random activation of a few or many LTR retrotransposons families in particular lineages over evolutionary time, with some families especially prone to future activation and hyper-amplification.

Towards decoding the conifer giga-genome

Several new initiatives have been launched recently to sequence conifer genomes including pines, spruces and Douglas-fir. Owing to the very large genome sizes ranging from 18 to 35 gigabases, sequencing even a single conifer genome had been considered unattainable until the recent throughput increases and cost reductions afforded by next generation sequencers. The purpose of this review is to describe the context for these new initiatives. A knowledge foundation has been acquired in several conifers of commercial and ecological interest through large-scale cDNA analyses, construction of genetic maps and gene mapping studies aiming to link phenotype and genotype. Exploratory sequencing in pines and spruces have pointed out some of the unique properties of these giga-genomes and suggested strategies that may be needed to extract value from their sequencing. The hope is that recent and pending developments in sequencing technology will contribute to rapidly filling the knowledge vacuum surrounding their structure, contents and evolution. Researchers are also making plans to use comparative analyses that will help to turn the data into a valuable resource for enhancing and protecting the world's conifer forests.

Environmental heterogeneity and phenotypic divergence: can heritable epigenetic variation aid speciation?

The dualism of genetic predisposition and environmental influences, their interactions, and respective roles in shaping the phenotype have been a hot topic in biological sciences for more than two centuries. Heritable epigenetic variation mediates between relatively slowly accumulating mutations in the DNA sequence and ephemeral adaptive responses to stress, thereby providing mechanisms for achieving stable, but potentially rapidly evolving phenotypic diversity as a response to environmental stimuli. This suggests that heritable epigenetic signals can play an important role in evolutionary processes, but so far this hypothesis has not been rigorously tested. A promising new area of research focuses on the interaction between the different molecular levels that produce phenotypic variation in wild, closely-related taxa that lack genome-wide genetic differentiation. By pinpointing specific adaptive traits and investigating the mechanisms responsible for phenotypic differentiation, such study systems could allow profound insights into the role of epigenetics in the evolution and stabilization of phenotypic discontinuities, and could add to our understanding of adaptive strategies to diverse environmental conditions and their dynamics.

Rare target enrichment for ultrasensitive PCR detection using cot-rehybridization and duplex-specific nuclease

Nucleic acid detection by polymerase chain reaction (PCR) is invaluable for the detection of dilute and rare sequences, including pathogens and infrequent species in complex clinical and environmental backgrounds. The presence of excess complex background nucleic acid can reduce sensitivity and specificity. This is because mispriming can cause failure of the amplification reaction. Here we describe a new approach to ultrasensitive PCR detection, using enrichment of rare target nucleic acid from abundant background by combining the classic technique of cot-rehybridization to convert the abundant background to double-stranded form, with the use of a newly described, highly processive duplex-specific crab nuclease. We show that trace sequences in a vast excess of background DNA can be undetectable by PCR, independent of the amount of the mixture added to the PCR, and that these sequences can be made detectable by background suppression using this method.

Dental caries pathogenicity: a genomic and metagenomic perspective

In this review we address the subject of dental caries pathogenicity from a genomic and metagenomic perspective. The application of genomic technologies is certain to yield novel insights into the relationship between the bacterial flora, dental health and disease. Three primary attributes of bacterial species are thought to have direct impact on caries development, these include: adherence on tooth surfaces (biofilm formation), acid production and acid tolerance. Attempts to define the specific aetiological agents of dental caries have proven to be elusive, supporting the notion that caries aetiology is perhaps complex and multi-faceted. The recently introduced Human Microbiome Project (HMP) that endeavors to characterise the micro-organisms living in and on the human body is likely to shed new light on these questions and improve our understanding of polymicrobial disease, microbial ecology in the oral cavity and provide new avenues for therapeutic and molecular diagnostics developments.

Characterization of the genome of bald cypress

BACKGROUND

Bald cypress (Taxodium distichum var. distichum) is a coniferous tree of tremendous ecological and economic importance. It is a member of the family Cupressaceae which also includes cypresses, redwoods, sequoias, thujas, and junipers. While the bald cypress genome is more than three times the size of the human genome, its 1C DNA content is amongst the smallest of any conifer. To learn more about the genome of bald cypress and gain insight into the evolution of Cupressaceae genomes, we performed a Cot analysis and used Cot filtration to study Taxodium DNA. Additionally, we constructed a 6.7 genome-equivalent BAC library that we screened with known Taxodium genes and select repeats.

RESULTS

The bald cypress genome is composed of 90% repetitive DNA with most sequences being found in low to mid copy numbers. The most abundant repeats are found in fewer than 25,000 copies per genome. Approximately 7.4% of the genome is single/low-copy DNA (i.e., sequences found in 1 to 5 copies). Sequencing of highly repetitive Cot clones indicates that most Taxodium repeats are highly diverged from previously characterized plant repeat sequences. The bald cypress BAC library consists of 606,336 clones (average insert size of 113 kb) and collectively provides 6.7-fold genome equivalent coverage of the bald cypress genome. Macroarray screening with known genes produced, on average, about 1.5 positive clones per probe per genome-equivalent. Library screening with Cot-1 DNA revealed that approximately 83% of BAC clones contain repetitive sequences iterated 103 to 104 times per genome.

CONCLUSIONS

The BAC library for bald cypress is the first to be generated for a conifer species outside of the family Pinaceae. The Taxodium BAC library was shown to be useful in gene isolation and genome characterization and should be an important tool in gymnosperm comparative genomics, physical mapping, genome sequencing, and gene/polymorphism discovery. The single/low-copy (SL) component of bald cypress is 4.6 times the size of the Arabidopsis genome. As suggested for other gymnosperms, the large amount of SL DNA in Taxodium is likely the result of divergence among ancient repeat copies and gene/pseudogene duplication.

Sequence-based marker development in wheat: advances and applications to breeding

In the past two decades, the wheat community has made remarkable progress in developing molecular resources for breeding. A wide variety of molecular tools has been established to accelerate genetic and physical mapping for facilitating the efficient identification of molecular markers linked to genes and QTL of agronomic interest. Already, wheat breeders are benefiting from a wide range of techniques to follow the introgression of the most favorable alleles in elite material and develop improved varieties. Breeders soon will be able to take advantage of new technological developments based on Next Generation Sequencing. In this paper, we review the molecular toolbox available to wheat scientists and breeders for performing fundamental genomic studies and breeding. Special emphasis is given on the production and detection of single nucleotide polymorphisms (SNPs) that should enable a step change in saturating the wheat genome for more efficient genetic studies and for the development of new selection methods. The perspectives offered by the access to an ordered full genome sequence for further marker development and enhanced precision breeding is also discussed. Finally, we discuss the advantages and limitations of marker-assisted selection for supporting wheat improvement.

Inference of subgenomic origin of BACs in an interspecific hybrid sugarcane cultivar by overlapping oligonucleotide hybridizations

Sugarcane (Saccharum spp.) breeders in the early 20th century made remarkable progress in increasing yield and disease resistance by crossing Saccharum spontaneum L., a wild relative, to Saccharum officinarum L., a traditional cultivar. Modern sugarcane cultivars have approximately 71%-83% of their chromosomes originating from S. officinarum, approximately 10%-21% from S. spontaneum, and approximately 2%-13% recombinant or translocated chromosomes. In the present work, C(0)t-based cloning and sequencing (CBCS) was implemented to further explore highly repetitive DNA and to seek species-specific repeated DNA in both S. officinarum and S. spontaneum. For putatively species-specific sequences, overlappping oligonucleotide probes (overgos) were designed and hybridized to BAC filters from the interspecific hybrid sugarcane cultivar 'R570' to try to deduce parental origins of BAC clones. We inferred that 12 967 BACs putatively originated from S. officinarum and 5117 BACs from S. spontaneum. Another 1103 BACs were hybridized by both species-specific overgos, too many to account for by conventional recombination, thus suggesting ectopic recombination and (or) translocation of DNA elements. Constructing a low C(0)t library is useful to collect highly repeated DNA sequences and to search for potentially species-specific molecular markers, especially among recently diverged species. Even in the absence of repeat families that are species-specific in their entirety, the identification of localized variations within consensus sequences, coupled with the site specificity of short synthetic overgos, permits researchers to monitor species-specific or species-enriched variants.

Use of methylation filtration and C(0)t fractionation for analysis of genome composition and comparative genomics in bread wheat

We investigated the compositional and structural differences in sequences derived from different fractions of wheat genomic DNA obtained using methylation filtration and C(0)t fractionation. Comparative analysis of these sequences revealed large compositional and structural variations in terms of GC content, different structural elements including repeat sequences (e.g., transposable elements and simple sequence repeats), protein coding genes, and non-coding RNA genes. A correlation between methylation status [determined on the basis of selective inclusion/exclusion in methylation-filtered (MF) library] of different repeat elements and expression level was observed. The expression levels were determined by comparing MF sequences with expressed sequence tags (ESTs) available in the public domain. Only a limited overlap among MF, high C(0)t (HC), and ESTs was observed, suggesting that these sequences may largely either represent the low-copy non-transcribed sequences or include genes with low expression levels. Thus, these results indicated a need to study MF and HC sequences along with ESTs to fully appreciate complexity of wheat gene space.

Adventures in the enormous: a 1.8 million clone BAC library for the 21.7 Gb genome of loblolly pine

Loblolly pine (LP; Pinus taeda L.) is the most economically important tree in the U.S. and a cornerstone species in southeastern forests. However, genomics research on LP and other conifers has lagged behind studies on flowering plants due, in part, to the large size of conifer genomes. As a means to accelerate conifer genome research, we constructed a BAC library for the LP genotype 7-56. The LP BAC library consists of 1,824,768 individually-archived clones making it the largest single BAC library constructed to date, has a mean insert size of 96 kb, and affords 7.6X coverage of the 21.7 Gb LP genome. To demonstrate the efficacy of the library in gene isolation, we screened macroarrays with overgos designed from a pine EST anchored on LP chromosome 10. A positive BAC was sequenced and found to contain the expected full-length target gene, several gene-like regions, and both known and novel repeats. Macroarray analysis using the retrotransposon IFG-7 (the most abundant repeat in the sequenced BAC) as a probe indicates that IFG-7 is found in roughly 210,557 copies and constitutes about 5.8% or 1.26 Gb of LP nuclear DNA; this DNA quantity is eight times the Arabidopsis genome. In addition to its use in genome characterization and gene isolation as demonstrated herein, the BAC library should hasten whole genome sequencing of LP via next-generation sequencing strategies/technologies and facilitate improvement of trees through molecular breeding and genetic engineering. The library and associated products are distributed by the Clemson University Genomics Institute (www.genome.clemson.edu).

[Construction of a methylation filtration library in Hevea brasiliensis.]

In order to enrich gene encoding region of Hevea brasiliensis, a methylation filtration library was constructed using Escherichia coli McrBC restriction-modification system. The titers of the non-amplified library and the amplified library were 2.6×106 pfu/ml and 9.0×109, respectively. The rate of positive clones was 86.4%. The lengths of inserted DNA sequence ranged from 1 kb to 2.5 kb and the average size of inserts was 1.2 kb. One hundred clones were selected randomly for sequencing, resulting in splicing out of 81 non-redundant sequences, including 6 contigs and 75 singlets. The redundancy was 17.35%. Blast analysis showed that 39.5% of non-redundant sequences were homologous with the Nr database, 14.81% with the EST database, and 32.1% were unknown sequences. Some sequences were related genes for flowering, insect and disease resistance. Therefore, the rubber tree methylation library is helpful for discovery and cloning of functional genes.

Normalization of genomic DNA using duplex-specific nuclease

An application of duplex-specific nuclease (DSN) normalization technology to whole-genome shotgun sequencing of genomes with a large proportion of repetitive DNA is described. The method uses a thermostable DSN from the Kamchatka crab that specifically hydrolyzes dsDNA. In model experiments on human genomic DNA, we demonstrated that DSN normalization of double-stranded DNA formed during C0t analysis is effective against abundant repetitive sequences with high sequence identity, while retaining highly divergent repeats and coding regions at base-line levels. Thus, DSN normalization applied to C0t analysis can be used to eliminate evolutionarily young repetitive elements from genomic DNA before sequencing, and should prove invaluable in studies of large eukaryotic genomes, such as those of higher plants.

The Tomato Sequencing Project, the first cornerstone of the International Solanaceae Project (SOL)

The genome of tomato (Solanum lycopersicum) is being sequenced by an international consortium of 10 countries (Korea, China, the United Kingdom, India, The Netherlands, France, Japan, Spain, Italy and the United States) as part of a larger initiative called the ‘International Solanaceae Genome Project (SOL): Systems Approach to Diversity and Adaptation’. The goal of this grassroots initiative, launched in November 2003, is to establish a network of information, resources and scientists to ultimately tackle two of the most significant questions in plant biology and agriculture: (1) How can a common set of genes/proteins give rise to a wide range of morphologically and ecologically distinct organisms that occupy our planet? (2) How can a deeper understanding of the genetic basis of plant diversity be harnessed to better meet the needs of society in an environmentally friendly and sustainable manner? The Solanaceae and closely related species such as coffee, which are included in the scope of the SOL project, are ideally suited to address both of these questions. The first step of the SOL project is to use an ordered BAC approach to generate a high quality sequence for the euchromatic portions of the tomato as a reference for the Solanaceae. Due to the high level of macro and micro-synteny in the Solanaceae the BAC-by-BAC tomato sequence will form the framework for shotgun sequencing of other species. The starting point for sequencing the genome is BACs anchored to the genetic map by overgo hybridization and AFLP technology. The overgos are derived from approximately 1500 markers from the tomato high density F2-2000 genetic map (http://sgn.cornell.edu/). These seed BACs will be used as anchors from which to radiate the tiling path using BAC end sequence data. Annotation will be performed according to SOL project guidelines. All the information generated under the SOL umbrella will be made available in a comprehensive website. The information will be interlinked with the ultimate goal that the comparative biology of the Solanaceae—and beyond—achieves a context that will facilitate a systems biology approach.

Reassociation kinetics-based approach for partial genome sequencing of the cattle tick, Rhipicephalus (Boophilus) microplus

Background

The size and repetitive nature of the Rhipicephalus microplus genome makes obtaining a full genome sequence fiscally and technically problematic. To selectively obtain gene-enriched regions of this tick's genome, Cot filtration was performed, and Cot-filtered DNA was sequenced via 454 FLX pyrosequencing.

Results

The sequenced Cot-filtered genomic DNA was assembled with an EST-based gene index of 14,586 unique entries where each EST served as a potential "seed" for scaffold formation. The new sequence assembly extended the lengths of 3,913 of the 14,586 gene index entries. Over half of the extensions corresponded to extensions of over 30 amino acids. To survey the repetitive elements in the tick genome, the complete sequences of five BAC clones were determined. Both Class I and II transposable elements were found. Comparison of the BAC and Cot filtration data indicates that Cot filtration was highly successful in filtering repetitive DNA out of the genomic DNA used in 454 sequencing.

Conclusion

Cot filtration is a very useful strategy to incorporate into genome sequencing projects on organisms with large genome sizes and which contain high percentages of repetitive, difficult to assemble, genomic DNA. Combining the Cot selection approach with 454 sequencing and assembly with a pre-existing EST database as seeds resulted in extensions of 27% of the members of the EST database.

Analysis of a c0t-1 library enables the targeted identification of minisatellite and satellite families in Beta vulgaris

BACKGROUND

Repetitive DNA is a major fraction of eukaryotic genomes and occurs particularly often in plants. Currently, the sequencing of the sugar beet (Beta vulgaris) genome is under way and knowledge of repetitive DNA sequences is critical for the genome annotation. We generated a c0t-1 library, representing highly to moderately repetitive sequences, for the characterization of the major B. vulgaris repeat families. While highly abundant satellites are well-described, minisatellites are only poorly investigated in plants. Therefore, we focused on the identification and characterization of these tandemly repeated sequences.

RESULTS

Analysis of 1763 c0t-1 DNA fragments, providing 442 kb sequence data, shows that the satellites pBV and pEV are the most abundant repeat families in the B. vulgaris genome while other previously described repeats show lower copy numbers. We isolated 517 novel repetitive sequences and used this fraction for the identification of minisatellite and novel satellite families. Bioinformatic analysis and Southern hybridization revealed that minisatellites are moderately to highly amplified in B. vulgaris. FISH showed a dispersed localization along most chromosomes clustering in arrays of variable size and number with exclusion and depletion in distinct regions.

CONCLUSION

The c0t-1 library represents major repeat families of the B. vulgaris genome, and analysis of the c0t-1 DNA was proven to be an efficient method for identification of minisatellites. We established, so far, the broadest analysis of minisatellites in plants and observed their chromosomal localization providing a background for the annotation of the sugar beet genome and for the understanding of the evolution of minisatellites in plant genomes.

Genome-wide comparative analysis of the Brassica rapa gene space reveals genome shrinkage and differential loss of duplicated genes after whole genome triplication

Euchromatic regions of the Brassica rapa genome were sequenced and mapped onto the corresponding regions in the Arabidopsis thaliana genome.

Background

Brassica rapa is one of the most economically important vegetable crops worldwide. Owing to its agronomic importance and phylogenetic position, B. rapa provides a crucial reference to understand polyploidy-related crop genome evolution. The high degree of sequence identity and remarkably conserved genome structure between Arabidopsis and Brassica genomes enables comparative tiling sequencing using Arabidopsis sequences as references to select the counterpart regions in B. rapa, which is a strong challenge of structural and comparative crop genomics.

Results

We assembled 65.8 megabase-pairs of non-redundant euchromatic sequence of B. rapa and compared this sequence to the Arabidopsis genome to investigate chromosomal relationships, macrosynteny blocks, and microsynteny within blocks. The triplicated B. rapa genome contains only approximately twice the number of genes as in Arabidopsis because of genome shrinkage. Genome comparisons suggest that B. rapa has a distinct organization of ancestral genome blocks as a result of recent whole genome triplication followed by a unique diploidization process. A lack of the most recent whole genome duplication (3R) event in the B. rapa genome, atypical of other Brassica genomes, may account for the emergence of B. rapa from the Brassica progenitor around 8 million years ago.

Conclusions

This work demonstrates the potential of using comparative tiling sequencing for genome analysis of crop species. Based on a comparative analysis of the B. rapa sequences and the Arabidopsis genome, it appears that polyploidy and chromosomal diploidization are ongoing processes that collectively stabilize the B. rapa genome and facilitate its evolution.

Genome analysis and genetic enhancement of tomato

The Solanaceae is an important family of vegetable crops, ornamentals and medicinal plants. Tomato has served as a model member of this family largely because of its enriched cytogenetic, genetic, as well as physical, maps. Mapping has helped in cloning several genes of importance such as Pto, responsible for resistance against bacterial speck disease, Mi-1.2 for resistance against nematodes, and fw2.2 QTL for fruit weight. A high-throughput genome-sequencing program has been initiated by an international consortium of 10 countries. Since heterochromatin has been found to be concentrated near centromeres, the consortium is focusing on sequencing only the gene-rich euchromatic region. Genomes of the members of Solanaceae show a significant degree of synteny, suggesting that the tomato genome sequence would help in the cloning of genes for important traits from other Solanaceae members as well. ESTs from a large number of cDNA libraries have been sequenced, and microarray chips, in conjunction with wide array of ripening mutants, have contributed immensely to the understanding of the fruit-ripening phenomenon. Work on the analysis of the tomato proteome has also been initiated. Transgenic tomato plants with improved abiotic stress tolerance, disease resistance and insect resistance, have been developed. Attempts have also been made to develop tomato as a bioreactor for various pharmaceutical proteins. However, control of fruit quality and ripening remains an active and challenging area of research. Such efforts should pave the way to improve not only tomato, but also other solanaceous crops.

[Advances of gene enrichment in plant genome]

The genome size varies greatly in higher plants. Repetitive sequences account for most of the large plant genomes while low-copy or single copy genic sequences, referred to as gene space, take up only a small portion of the genomes. Considering the large amount of repetitive sequences, it is a great challenge to obtain genic sequences using high-throughout methods in non-model plants bearing large genomes. Currently, several approaches have been developed for gene enrichment on a genome-wide scale, such as cDNA library, methylation filtration library, high Cot library and transposon tagging. Here, we reviewed the technical principles, advantages and disadvantages of these methods, as well as the recent development of methylation filtration technology. An in-depth discussion was performed for selection of one method or combination of methods according to the research objectives and plant materials, especially for plants with large genomes.