Enrichment of genomic DNA for polymorphism detection in a non-model highly polyploid crop plant

A Large-Scale Candidate-Gene Association Mapping for Drought Tolerance and Agronomic Traits in Sugarcane

Dissection of the genetic loci controlling drought tolerance traits with a complex genetic inheritance is important for drought-tolerant sugarcane improvement. In this study, we conducted a large-scale candidate gene association study of 649 candidate genes in a sugarcane diversity panel to identify genetic variants underlying agronomic traits and drought tolerance indices evaluated in plant cane and ratoon cane under water-stressed (WS) and non-stressed (NS) environments. We identified 197 significant marker-trait associations (MTAs) in 141 candidate genes associated with 18 evaluated traits with the Bonferroni correction threshold (α = 0.05). Out of the total, 95 MTAs in 78 candidate genes and 62 MTAs in 58 candidate genes were detected under NS and WS conditions, respectively. Most MTAs were found only in specific water regimes and crop seasons. These MTAs explained 7.93-30.52% of phenotypic variation. Association mapping results revealed that 34, 59, and 104 MTAs involved physiological and molecular adaptation, phytohormone metabolism, and drought-inducible genes. They identified 19 pleiotropic genes associated with more than one trait and many genes related to drought tolerance indices. The genetic and genomic resources identified in this study will enable the combining of yield-related traits and sugar-related traits with agronomic value to optimize the yield of sugarcane cultivars grown under drought-stressed and non-stressed environments.

Assembly of Long-Adapter Single-Strand Oligonucleotide (LASSO) Probes for Massively Parallel Capture of Kilobase Size DNA Targets

Genome DNA sequencing has become an affordable means to resolve questions about the genetic background of life. However, the biological functions of many DNA-encoded sequences are still relatively unknown. A highly scalable and cost-effective cloning method to select natural DNA targets from genomic templates is therefore urgently needed to enable rapid understanding of the biological products of genomes. One such method involves LASSO probes, which are long single-stranded DNA oligonucleotides designed with a universal adapter that is used to link two sequences that are complementary to a genomic target of interest. Through a pooled assembly method, LASSOs can be made for multiplex DNA capture. Herein, we describe a robust, efficient method to assemble LASSO probe libraries using a Cre-recombinase-mediated reaction and a protocol for multiplex genome target capture. The starting components are a pre-LASSO probe library comprising short DNA oligo pools designed in silico and an Escherichia coli plasmid (pLASSO) that incorporates the pre-LASSO library. Through internal recombination of pLASSO with its inserts, a mature LASSO library in final configuration can be made with high purity. Assembly of a LASSO probe library takes 4 days, and target capture can be performed in a single day. With an exponentially growing list of new genomes available for investigation, this method can enable the rapid production of ORFeome libraries for high-throughput screening to identify biological functions as a complementary approach to understand genome functional biology. © 2021 Wiley Periodicals LLC. Basic Protocol 1: Assembly of LASSO probes Support Protocol 1: Generation of pLASSO vectors Support Protocol 2: Preparation of pre-LASSOs Basic Protocol 2: Massively parallel capture of large DNAs using LASSO probes.

An overview of the transcriptional responses of two tolerant and susceptible sugarcane cultivars to borer (Diatraea saccharalis) infestation

Diatraea saccharalis constitutes a threat to the sugarcane productivity, and obtaining borer tolerant cultivars is an alternative method of control. Although there are studies about the relationship between the interaction of D. saccharalis with sugarcane, little is known about the molecular and genomic basis of defense mechanisms that confer tolerance to sugarcane cultivars. Here, we analyzed the transcriptional profile of two sugarcane cultivars in response to borer attack, RB867515 and SP80-3280, which are considered tolerant and sensitive to the borer attack, respectively. A sugarcane genome and transcriptome were used for read mapping. Differentially expressed transcripts and genes were identified and termed to as DETs and DEGs, according to the sugarcane database adopted. A total of 745 DETs and 416 DEGs were identified (log₂|ratio| > 0.81; FDR corrected P value ≤ 0.01) after borer infestation. Following annotation of up- and down-regulated DETs and DEGs by similarity searches, the sugarcane cultivars demonstrated an up-regulation of jasmonic acid (JA), ethylene (ET), and defense protein genes, as well as a down-regulation of pathways involved in photosynthesis and energy metabolism. The expression analysis also highlighted that RB867515 cultivar is possibly more transcriptionally activated after 12 h from infestation than SP80-3280, which could imply in quicker responses by probably triggering more defense-related genes and mediating metabolic pathways to cope with borer attack.

Target enrichment sequencing of 307 germplasm accessions identified ancestry of ancient and modern hybrids and signatures of adaptation and selection in sugarcane (Saccharum spp.), a 'sweet' crop with 'bitter' genomes

Sugarcane (Saccharum spp.) is a highly energy-efficient crop primarily for sugar and bio-ethanol production. Sugarcane genetics and cultivar improvement have been extremely challenging largely due to its complex genomes with high polyploidy levels. In this study, we deeply sequenced the coding regions of 307 sugarcane germplasm accessions. Nearly five million sequence variations were catalogued. The average of 98× sequence depth enabled different allele dosages of sequence variation to be differentiated in this polyploid collection. With selected high-quality genome-wide SNPs, we performed population genomic studies and environmental association analysis. Results illustrated that the ancient sugarcane hybrids, S. barberi and S. sinense, and modern sugarcane hybrids are significantly different in terms of genomic compositions, hybridization processes and their potential ancestry contributors. Linkage disequilibrium (LD) analysis showed a large extent of LD in sugarcane, with 962.4 Kbp, 2739.2 Kbp and 3573.6 Kbp for S. spontaneum, S. officinarum and modern S. hybrids respectively. Candidate selective sweep regions and genes were identified during domestication and historical selection processes of sugarcane in addition to genes associated with environmental variables at the original locations of the collection. This research provided an extensive amount of genomic resources for sugarcane community and the in-depth population genomic analyses shed light on the breeding and evolution history of sugarcane, a highly polyploid species.

Advances in Clinical Sample Preparation for Identification and Characterization of Bacterial Pathogens Using Metagenomics

Whole genome sequencing (WGS) plays an increasing role in communicable disease control through high-resolution outbreak tracing, laboratory surveillance and diagnostics. However, WGS has traditionally relied on microbial culture in order to obtain pathogen specific DNA for sequencing. This has severely limited the application of whole genome sequencing on pathogens with fastidious culturing requirements. In addition, the widespread adoption of culture-independent diagnostic tests has reduced availability of cultured isolates for confirmatory testing and surveillance. These recent developments have created demand for the implementation of techniques enabling direct sequencing of microbial genomes in clinical samples without having to culture an isolate. However, sequencing of specific organisms from clinical samples can be affected by high levels of contaminating DNA from the host and other commensal microorganisms. Several methods have been introduced for selective lysis of host cells and/or separate specific organisms from a clinical sample. This review examines the different approaches for sample preparation that have been used in diagnostic and public health laboratories for metagenomic sequencing.

The Challenge of Analyzing the Sugarcane Genome

Reference genome sequences have become key platforms for genetics and breeding of the major crop species. Sugarcane is probably the largest crop produced in the world (in weight of crop harvested) but lacks a reference genome sequence. Sugarcane has one of the most complex genomes in crop plants due to the extreme level of polyploidy. The genome of modern sugarcane hybrids includes sub-genomes from two progenitors Saccharum officinarum and S. spontaneum with some chromosomes resulting from recombination between these sub-genomes. Advancing DNA sequencing technologies and strategies for genome assembly are making the sugarcane genome more tractable. Advances in long read sequencing have allowed the generation of a more complete set of sugarcane gene transcripts. This is supporting transcript profiling in genetic research. The progenitor genomes are being sequenced. A monoploid coverage of the hybrid genome has been obtained by sequencing BAC clones that cover the gene space of the closely related sorghum genome. The complete polyploid genome is now being sequenced and assembled. The emerging genome will allow comparison of related genomes and increase understanding of the functioning of this polyploidy system. Sugarcane breeding for traditional sugar and new energy and biomaterial uses will be enhanced by the availability of these genomic resources.

Mining sequence variations in representative polyploid sugarcane germplasm accessions

Background

Sugarcane (Saccharum spp.) is one of the most important economic crops because of its high sugar production and biofuel potential. Due to the high polyploid level and complex genome of sugarcane, it has been a huge challenge to investigate genomic sequence variations, which are critical for identifying alleles contributing to important agronomic traits. In order to mine the genetic variations in sugarcane, genotyping by sequencing (GBS), was used to genotype 14 representative Saccharum complex accessions. GBS is a method to generate a large number of markers, enabled by next generation sequencing (NGS) and the genome complexity reduction using restriction enzymes.

Results

To use GBS for high throughput genotyping highly polyploid sugarcane, the GBS analysis pipelines in 14 Saccharum complex accessions were established by evaluating different alignment methods, sequence variants callers, and sequence depth for single nucleotide polymorphism (SNP) filtering. By using the established pipeline, a total of 76,251 non-redundant SNPs, 5642 InDels, 6380 presence/absence variants (PAVs), and 826 copy number variations (CNVs) were detected among the 14 accessions. In addition, non-reference based universal network enabled analysis kit and Stacks de novo called 34,353 and 109,043 SNPs, respectively. In the 14 accessions, the percentages of single dose SNPs ranged from 38.3% to 62.3% with an average of 49.6%, much more than the portions of multiple dosage SNPs. Concordantly called SNPs were used to evaluate the phylogenetic relationship among the 14 accessions. The results showed that the divergence time between the Erianthus genus and the Saccharum genus was more than 10 million years ago (MYA). The Saccharum species separated from their common ancestors ranging from 0.19 to 1.65 MYA.

Conclusions

The GBS pipelines including the reference sequences, alignment methods, sequence variant callers, and sequence depth were recommended and discussed for the Saccharum complex and other related species. A large number of sequence variations were discovered in the Saccharum complex, including SNPs, InDels, PAVs, and CNVs. Genome-wide SNPs were further used to illustrate sequence features of polyploid species and demonstrated the divergence of different species in the Saccharum complex. The results of this study showed that GBS was an effective NGS-based method to discover genomic sequence variations in highly polyploid and heterozygous species.

Target enrichment sequencing in cultivated peanut (Arachis hypogaea L.) using probes designed from transcript sequences

Enabled by the next generation sequencing, target enrichment sequencing (TES) is a powerful method to enrich genomic regions of interest and to identify sequence variations. The objective of this study was to explore the feasibility of probe design from transcript sequences for TES application in calling sequence variants in peanut, an important allotetraploid crop with a large genome size. In this study, we applied an in-solution hybridization method to enrich DNA sequences of seven peanut genotypes. Our results showed that it is feasible to apply TES with probes designed from transcript sequences in polyploid peanut. Using a set of 31,123 probes, a total of 5131 and 7521 genes were targeted in peanut A and B genomes, respectively. For each genotype used in this study, the probe target capture regions were efficiently covered with high depth. The average on-target rate of sequencing reads was 42.47%, with a significant amount of off-target reads coming from genomic regions homologous to target regions. In this study, when given predefined genomic regions of interest and the same amount of sequencing data, TES provided the highest coverage of target regions when compared to whole genome sequencing, RNA sequencing, and genotyping by sequencing. Single nucleotide polymorphism (SNP) calling and subsequent validation revealed a high validation rate (85.71%) of homozygous SNPs, providing valuable markers for peanut genotyping. This study demonstrated the success of applying TES for SNP identification in peanut, which shall provide valuable suggestions for TES application in other non-model species without a genome reference available.

Targeted Capture of Complete Coding Regions across Divergent Species

Despite continued advances in sequencing technologies, there is a need for methods that can efficiently sequence large numbers of genes from diverse species. One approach to accomplish this is targeted capture (hybrid enrichment). While these methods are well established for genome resequencing projects, cross-species capture strategies are still being developed and generally focus on the capture of conserved regions, rather than complete coding regions from specific genes of interest. The resulting data is thus useful for phylogenetic studies, but the wealth of comparative data that could be used for evolutionary and functional studies is lost. Here, we design and implement a targeted capture method that enables recovery of complete coding regions across broad taxonomic scales. Capture probes were designed from multiple reference species and extensively tiled in order to facilitate cross-species capture. Using novel bioinformatics pipelines we were able to recover nearly all of the targeted genes with high completeness from species that were up to 200 myr divergent. Increased probe diversity and tiling for a subset of genes had a large positive effect on both recovery and completeness. The resulting data produced an accurate species tree, but importantly this same data can also be applied to studies of molecular evolution and function that will allow researchers to ask larger questions in broader phylogenetic contexts. Our method demonstrates the utility of cross-species approaches for the capture of full length coding sequences, and will substantially improve the ability for researchers to conduct large-scale comparative studies of molecular evolution and function.

From Genomes to GENE-omes: Exome Sequencing Concept and Applications in Crop Improvement

Exome sequencing represents targeted capture and sequencing of 1-2% of 'high-value genomic regions' (subset of the genome) which are enriched for functional variants and harbors low level of repetitive regions. We discuss here an overview of exome sequencing, ways to approach plant exomes, and advantages and applicability of this powerful approach in deciphering functional regions of genomes. Though initially this approach was developed as an alternative to whole genome sequencing (WGS), but the multitude of benefits conferred by sequence capture via hybridization approaches created a niche for itself to solve many of biological riddles, particularly for resolving phylogenetic distances. The technique has also proved to be successful in understanding the basis of natural and induced molecular variation, marker development and developing genomic resources for complex, wild and non-model species, which are still intractable for WGS efforts. Thus, with profound applications of this powerful sequencing strategy, near future is expected to witness a collective expansion of both techniques, i.e., sequence capture via hybridization for evolutionary and ecological research and WGS approaches for its universal accessibility.

The two chromosomes of the mitochondrial genome of a sugarcane cultivar: assembly and recombination analysis using long PacBio reads

Sugarcane accounts for a large portion of the worlds sugar production. Modern commercial cultivars are complex hybrids of S. officinarum and several other Saccharum species. Historical records identify New Guinea as the origin of S. officinarum and that a small number of plants originating from there were used to generate all modern commercial cultivars. The mitochondrial genome can be a useful way to identify the maternal origin of commercial cultivars. We have used the PacBio RSII to sequence and assemble the mitochondrial genome of a South East Asian commercial cultivar, known as Khon Kaen 3. The long read length of this sequencing technology allowed for the mitochondrial genome to be assembled into two distinct circular chromosomes with all repeat sequences spanned by individual reads. Comparison of five commercial hybrids, two S. officinarum and one S. spontaneum to our assembly reveals no structural rearrangements between our assembly, the commercial hybrids and an S. officinarum from New Guinea. The S. spontaneum, from India, and one sample of S. officinarum (unknown origin) are substantially rearranged and have a large number of homozygous variants. This supports the record that S. officinarum plants from New Guinea are the maternal source of all modern commercial hybrids.

Sequence capture by hybridization to explore modern and ancient genomic diversity in model and nonmodel organisms

The recent expansion of next-generation sequencing has significantly improved biological research. Nevertheless, deep exploration of genomes or metagenomic samples remains difficult because of the sequencing depth and the associated costs required. Therefore, different partitioning strategies have been developed to sequence informative subsets of studied genomes. Among these strategies, hybridization capture has proven to be an innovative and efficient tool for targeting and enriching specific biomarkers in complex DNA mixtures. It has been successfully applied in numerous areas of biology, such as exome resequencing for the identification of mutations underlying Mendelian or complex diseases and cancers, and its usefulness has been demonstrated in the agronomic field through the linking of genetic variants to agricultural phenotypic traits of interest. Moreover, hybridization capture has provided access to underexplored, but relevant fractions of genomes through its ability to enrich defined targets and their flanking regions. Finally, on the basis of restricted genomic information, this method has also allowed the expansion of knowledge of nonreference species and ancient genomes and provided a better understanding of metagenomic samples. In this review, we present the major advances and discoveries permitted by hybridization capture and highlight the potency of this approach in all areas of biology.

Co-expression network analysis reveals transcription factors associated to cell wall biosynthesis in sugarcane

Sugarcane is a hybrid of Saccharum officinarum and Saccharum spontaneum, with minor contributions from other species in Saccharum and other genera. Understanding the molecular basis of cell wall metabolism in sugarcane may allow for rational changes in fiber quality and content when designing new energy crops. This work describes a comparative expression profiling of sugarcane ancestral genotypes: S. officinarum, S. spontaneum and S. robustum and a commercial hybrid: RB867515, linking gene expression to phenotypes to identify genes for sugarcane improvement. Oligoarray experiments of leaves, immature and intermediate internodes, detected 12,621 sense and 995 antisense transcripts. Amino acid metabolism was particularly evident among pathways showing natural antisense transcripts expression. For all tissues sampled, expression analysis revealed 831, 674 and 648 differentially expressed genes in S. officinarum, S. robustum and S. spontaneum, respectively, using RB867515 as reference. Expression of sugar transporters might explain sucrose differences among genotypes, but an unexpected differential expression of histones were also identified between high and low Brix° genotypes. Lignin biosynthetic genes and bioenergetics-related genes were up-regulated in the high lignin genotype, suggesting that these genes are important for S. spontaneum to allocate carbon to lignin, while S. officinarum allocates it to sucrose storage. Co-expression network analysis identified 18 transcription factors possibly related to cell wall biosynthesis while in silico analysis detected cis-elements involved in cell wall biosynthesis in their promoters. Our results provide information to elucidate regulatory networks underlying traits of interest that will allow the improvement of sugarcane for biofuel and chemicals production.

Natural Allelic Variations in Highly Polyploidy Saccharum Complex

Sugarcane (Saccharum spp.) is an important sugar and biofuel crop with high polyploid and complex genomes. The Saccharum complex, comprised of Saccharum genus and a few related genera, are important genetic resources for sugarcane breeding. A large amount of natural variation exists within the Saccharum complex. Though understanding their allelic variation has been challenging, it is critical to dissect allelic structure and to identify the alleles controlling important traits in sugarcane. To characterize natural variations in Saccharum complex, a target enrichment sequencing approach was used to assay 12 representative germplasm accessions. In total, 55,946 highly efficient probes were designed based on the sorghum genome and sugarcane unigene set targeting a total of 6 Mb of the sugarcane genome. A pipeline specifically tailored for polyploid sequence variants and genotype calling was established. BWA-mem and sorghum genome approved to be an acceptable aligner and reference for sugarcane target enrichment sequence analysis, respectively. Genetic variations including 1,166,066 non-redundant SNPs, 150,421 InDels, 919 gene copy number variations, and 1,257 gene presence/absence variations were detected. SNPs from three different callers (Samtools, Freebayes, and GATK) were compared and the validation rates were nearly 90%. Based on the SNP loci of each accession and their ploidy levels, 999,258 single dosage SNPs were identified and most loci were estimated as largely homozygotes. An average of 34,397 haplotype blocks for each accession was inferred. The highest divergence time among the Saccharum spp. was estimated as 1.2 million years ago (MYA). Saccharum spp. diverged from Erianthus and Sorghum approximately 5 and 6 MYA, respectively. The target enrichment sequencing approach provided an effective way to discover and catalog natural allelic variation in highly polyploid or heterozygous genomes.

Potential for Genetic Improvement of Sugarcane as a Source of Biomass for Biofuels

Sugarcane (Saccharum spp. hybrids) has great potential as a major feedstock for biofuel production worldwide. It is considered among the best options for producing biofuels today due to an exceptional biomass production capacity, high carbohydrate (sugar + fiber) content, and a favorable energy input/output ratio. To maximize the conversion of sugarcane biomass into biofuels, it is imperative to generate improved sugarcane varieties with better biomass degradability. However, unlike many diploid plants, where genetic tools are well developed, biotechnological improvement is hindered in sugarcane by our current limited understanding of the large and complex genome. Therefore, understanding the genetics of the key biofuel traits in sugarcane and optimization of sugarcane biomass composition will advance efficient conversion of sugarcane biomass into fermentable sugars for biofuel production. The large existing phenotypic variation in Saccharum germplasm and the availability of the current genomics technologies will allow biofuel traits to be characterized, the genetic basis of critical differences in biomass composition to be determined, and targets for improvement of sugarcane for biofuels to be established. Emerging options for genetic improvement of sugarcane for the use as a bioenergy crop are reviewed. This will better define the targets for potential genetic manipulation of sugarcane biomass composition for biofuels.

Captured metagenomics: large-scale targeting of genes based on 'sequence capture' reveals functional diversity in soils

Microbial enzyme diversity is a key to understand many ecosystem processes. Whole metagenome sequencing (WMG) obtains information on functional genes, but it is costly and inefficient due to large amount of sequencing that is required. In this study, we have applied a captured metagenomics technique for functional genes in soil microorganisms, as an alternative to WMG. Large-scale targeting of functional genes, coding for enzymes related to organic matter degradation, was applied to two agricultural soil communities through captured metagenomics. Captured metagenomics uses custom-designed, hybridization-based oligonucleotide probes that enrich functional genes of interest in metagenomic libraries where only probe-bound DNA fragments are sequenced. The captured metagenomes were highly enriched with targeted genes while maintaining their target diversity and their taxonomic distribution correlated well with the traditional ribosomal sequencing. The captured metagenomes were highly enriched with genes related to organic matter degradation; at least five times more than similar, publicly available soil WMG projects. This target enrichment technique also preserves the functional representation of the soils, thereby facilitating comparative metagenomics projects. Here, we present the first study that applies the captured metagenomics approach in large scale, and this novel method allows deep investigations of central ecosystem processes by studying functional gene abundances.

Development of highly reliable in silico SNP resource and genotyping assay from exome capture and sequencing: an example from black spruce (Picea mariana)

Picea mariana is a widely distributed boreal conifer across Canada and the subject of advanced breeding programmes for which population genomics and genomic selection approaches are being developed. Targeted sequencing was achieved after capturing P. mariana exome with probes designed from the sequenced transcriptome of Picea glauca, a distant relative. A high capture efficiency of 75.9% was reached although spruce has a complex and large genome including gene sequences interspersed by some long introns. The results confirmed the relevance of using probes from congeneric species to perform successfully interspecific exome capture in the genus Picea. A bioinformatics pipeline was developed including stringent criteria that helped detect a set of 97,075 highly reliable in silico SNPs. These SNPs were distributed across 14,909 genes. Part of an Infinium iSelect array was used to estimate the rate of true positives by validating 4267 of the predicted in silico SNPs by genotyping trees from P. mariana populations. The true positive rate was 96.2% for in silico SNPs, compared to a genotyping success rate of 96.7% for a set 1115 P. mariana control SNPs recycled from previous genotyping arrays. These results indicate the high success rate of the genotyping array and the relevance of the selection criteria used to delineate the new P. mariana in silico SNP resource. Furthermore, in silico SNPs were generally of medium to high frequency in natural populations, thus providing high informative value for future population genomics applications.

Culture-independent genome sequencing of clinical samples reveals an unexpected heterogeneity of infections by Chlamydia pecorum

Chlamydia pecorum is an important global pathogen of livestock, and it is also a significant threat to the long-term survival of Australia's koala populations. This study employed a culture-independent DNA capture approach to sequence C. pecorum genomes directly from clinical swab samples collected from koalas with chlamydial disease as well as from sheep with arthritis and conjunctivitis. Investigations into single-nucleotide polymorphisms within each of the swab samples revealed that a portion of the reads in each sample belonged to separate C. pecorum strains, suggesting that all of the clinical samples analyzed contained mixed populations of genetically distinct C. pecorum isolates. This observation was independent of the anatomical site sampled and the host species. Using the genomes of strains identified in each of these samples, whole-genome phylogenetic analysis revealed that a clade containing a bovine and a koala isolate is distinct from other clades comprised of livestock or koala C. pecorum strains. Providing additional evidence to support exposure of koalas to Australian livestock strains, two minor strains assembled from the koala swab samples clustered with livestock strains rather than koala strains. Culture-independent probe-based genome capture and sequencing of clinical samples provides the strongest evidence yet to suggest that naturally occurring chlamydial infections are comprised of multiple genetically distinct strains.

Future of portable devices for plant pathogen diagnosis

The demand for rapid and accurate diagnosis of plant diseases has risen in the last decade. On-site diagnosis of single or multiple pathogens using portable devices is the first step in this endeavour. Despite extensive attempts to develop portable devices for pathogen detection, current technologies are still restricted to detecting known pathogens with limited detection accuracy. Developing new detection techniques for rapid and accurate detection of multiple plant pathogens and their associated variants is essential. Recent single DNA sequencing technologies are a promising new avenue for developing future portable devices for plant pathogen detection. In this review, we detail the current progress in portable devices and technologies used for detecting plant pathogens, the current position of emerging sequencing technologies for analysis of plant genomics, and the future of portable devices for rapid pathogen diagnosis.

Exploring natural selection to guide breeding for agriculture

Climate change threatens reduced crop production and poses major challenges to food security. The breeding of climate-resilient crop varieties is increasingly urgent. Wild plant populations evolve to cope with changes in their environment due to the forces of natural selection. This adaptation may be followed over time in populations at the same site or explored by examining differences between populations growing in different environments or across an environmental gradient. Survival in the wild has important differences to the objective of agriculture to maximize crop yields. However, understanding the nature of adaptation in wild populations at the whole genome level may suggest strategies for crop breeding to deliver agricultural production with more resilience to climate variability.

Sugarcane genome sequencing by methylation filtration provides tools for genomic research in the genus Saccharum

Many economically important crops have large and complex genomes that hamper their sequencing by standard methods such as whole genome shotgun (WGS). Large tracts of methylated repeats occur in plant genomes that are interspersed by hypomethylated gene-rich regions. Gene-enrichment strategies based on methylation profiles offer an alternative to sequencing repetitive genomes. Here, we have applied methyl filtration with McrBC endonuclease digestion to enrich for euchromatic regions in the sugarcane genome. To verify the efficiency of methylation filtration and the assembly quality of sequences submitted to gene-enrichment strategy, we have compared assemblies using methyl-filtered (MF) and unfiltered (UF) libraries. The use of methy filtration allowed a better assembly by filtering out 35% of the sugarcane genome and by producing 1.5× more scaffolds and 1.7× more assembled Mb in length compared with unfiltered dataset. The coverage of sorghum coding sequences (CDS) by MF scaffolds was at least 36% higher than by the use of UF scaffolds. Using MF technology, we increased by 134× the coverage of gene regions of the monoploid sugarcane genome. The MF reads assembled into scaffolds that covered all genes of the sugarcane bacterial artificial chromosomes (BACs), 97.2% of sugarcane expressed sequence tags (ESTs), 92.7% of sugarcane RNA-seq reads and 98.4% of sorghum protein sequences. Analysis of MF scaffolds from encoded enzymes of the sucrose/starch pathway discovered 291 single-nucleotide polymorphisms (SNPs) in the wild sugarcane species, S. spontaneum and S. officinarum. A large number of microRNA genes was also identified in the MF scaffolds. The information achieved by the MF dataset provides a valuable tool for genomic research in the genus Saccharum and for improvement of sugarcane as a biofuel crop.

A high-density gene map of loblolly pine (Pinus taeda L.) based on exome sequence capture genotyping

Loblolly pine (Pinus taeda L.) is an economically and ecologically important conifer for which a suite of genomic resources is being generated. Despite recent attempts to sequence the large genome of conifers, their assembly and the positioning of genes remains largely incomplete. The interspecific synteny in pines suggests that a gene-based map would be useful to support genome assemblies and analysis of conifers. To establish a reference gene-based genetic map, we performed exome sequencing of 14729 genes on a mapping population of 72 haploid samples, generating a resource of 7434 sequence variants segregating for 3787 genes. Most markers are single-nucleotide polymorphisms, although short insertions/deletions and multiple nucleotide polymorphisms also were used. Marker segregation in the population was used to generate a high-density, gene-based genetic map. A total of 2841 genes were mapped to pine’s 12 linkage groups with an average of one marker every 0.58 cM. Capture data were used to detect gene presence/absence variations and position 65 genes on the map. We compared the marker order of genes previously mapped in loblolly pine and found high agreement. We estimated that 4123 genes had enough sequencing depth for reliable detection of markers, suggesting a high marker conversation rate of 92% (3787/4123). This is possible because a significant portion of the gene is captured and sequenced, increasing the chances of identifying a polymorphic site for characterization and mapping. This sub-centiMorgan genetic map provides a valuable resource for gene positioning on chromosomes and guide for the assembly of a reference pine genome.

Reduced representation approaches to interrogate genome diversity in large repetitive plant genomes

Technology and software improvements in the last decade now provide methodologies to access the genome sequence of not only a single accession, but also multiple accessions of plant species. This provides a means to interrogate species diversity at the genome level. Ample diversity among accessions in a collection of species can be found, including single-nucleotide polymorphisms, insertions and deletions, copy number variation and presence/absence variation. For species with small, non-repetitive rich genomes, re-sequencing of query accessions is robust, highly informative, and economically feasible. However, for species with moderate to large sized repetitive-rich genomes, technical and economic barriers prevent en masse genome re-sequencing of accessions. Multiple approaches to access a focused subset of loci in species with larger genomes have been developed, including reduced representation sequencing, exome capture and transcriptome sequencing. Collectively, these approaches have enabled interrogation of diversity on a genome scale for large plant genomes, including crop species important to worldwide food security.

Capturing sequence variation among flowering-time regulatory gene homologs in the allopolyploid crop species Brassica napus

Flowering, the transition from the vegetative to the generative phase, is a decisive time point in the lifecycle of a plant. Flowering is controlled by a complex network of transcription factors, photoreceptors, enzymes and miRNAs. In recent years, several studies gave rise to the hypothesis that this network is also strongly involved in the regulation of other important lifecycle processes ranging from germination and seed development through to fundamental developmental and yield-related traits. In the allopolyploid crop species Brassica napus, (genome AACC), homoeologous copies of flowering time regulatory genes are implicated in major phenological variation within the species, however the extent and control of intraspecific and intergenomic variation among flowering-time regulators is still unclear. To investigate differences among B. napus morphotypes in relation to flowering-time gene variation, we performed targeted deep sequencing of 29 regulatory flowering-time genes in four genetically and phenologically diverse B. napus accessions. The genotype panel included a winter-type oilseed rape, a winter fodder rape, a spring-type oilseed rape (all B. napus ssp. napus) and a swede (B. napus ssp. napobrassica), which show extreme differences in winter-hardiness, vernalization requirement and flowering behavior. A broad range of genetic variation was detected in the targeted genes for the different morphotypes, including non-synonymous SNPs, copy number variation and presence-absence variation. The results suggest that this broad variation in vernalization, clock and signaling genes could be a key driver of morphological differentiation for flowering-related traits in this recent allopolyploid crop species.

Sequence-based SNP genotyping in durum wheat

Marker development for marker-assisted selection in plant breeding is increasingly based on next-generation sequencing (NGS). However, marker development in crops with highly repetitive, complex genomes is still challenging. Here we applied sequence-based genotyping (SBG), which couples AFLP®-based complexity reduction to NGS, for de novo single nucleotide polymorphisms (SNP) marker discovery in and genotyping of a biparental durum wheat population. We identified 9983 putative SNPs in 6372 contigs between the two parents and used these SNPs for genotyping 91 recombinant inbred lines (RILs). Excluding redundant information from multiple SNPs per contig, 2606 (41%) markers were used for integration in a pre-existing framework map, resulting in the integration of 2365 markers over 2607 cM. Of the 2606 markers available for mapping, 91% were integrated in the pre-existing map, containing 708 SSRs, DArT markers, and SNPs from CRoPS technology, with a map-size increase of 492 cM (23%). These results demonstrate the high quality of the discovered SNP markers. With this methodology, it was possible to saturate the map at a final marker density of 0.8 cM/marker. Looking at the binned marker distribution (Figure 2), 63 of the 268 10-cM bins contained only SBG markers, showing that these markers are filling in gaps in the framework map. As to the markers that could not be used for mapping, the main reason was the low sequencing coverage used for genotyping. We conclude that SBG is a valuable tool for efficient, high-throughput and high-quality marker discovery and genotyping for complex genomes such as that of durum wheat.

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.

Application of large-scale sequencing to marker discovery in plants

Advances in DNA sequencing provide tools for efficient large-scale discovery of markers for use in plants. Discovery options include large-scale amplicon sequencing, transcriptome sequencing, gene-enriched genome sequencing and whole genome sequencing. Examples of each of these approaches and their potential to generate molecular markers for specific applications have been described. Sequencing the whole genome of parents identifies all the polymorphisms available for analysis in their progeny. Sequencing PCR amplicons of sets of candidate genes from DNA bulks can be used to define the available variation in these genes that might be exploited in a population or germplasm collection. Sequencing of the transcriptomes of genotypes varying for the trait of interest may identify genes with patterns of expression that could explain the phenotypic variation. Sequencing genomic DNA enriched for genes by hybridization with probes for all or some of the known genes simplifies sequencing and analysis of differences in gene sequences between large numbers of genotypes and genes especially when working with complex genomes. Examples of application of the above-mentioned techniques have been described.

Accessing complex crop genomes with next-generation sequencing

Many important crop species have genomes originating from ancestral or recent polyploidisation events. Multiple homoeologous gene copies, chromosomal rearrangements and amplification of repetitive DNA within large and complex crop genomes can considerably complicate genome analysis and gene discovery by conventional, forward genetics approaches. On the other hand, ongoing technological advances in molecular genetics and genomics today offer unprecedented opportunities to analyse and access even more recalcitrant genomes. In this review, we describe next-generation sequencing and data analysis techniques that vastly improve our ability to dissect and mine genomes for causal genes underlying key traits and allelic variation of interest to breeders. We focus primarily on wheat and oilseed rape, two leading examples of major polyploid crop genomes whose size or complexity present different, significant challenges. In both cases, the latest DNA sequencing technologies, applied using quite different approaches, have enabled considerable progress towards unravelling the respective genomes. Our ability to discover the extent and distribution of genetic diversity in crop gene pools, and its relationship to yield and quality-related traits, is swiftly gathering momentum as DNA sequencing and the bioinformatic tools to deal with growing quantities of genomic data continue to develop. In the coming decade, genomic and transcriptomic sequencing, discovery and high-throughput screening of single nucleotide polymorphisms, presence-absence variations and other structural chromosomal variants in diverse germplasm collections will give detailed insight into the origins, domestication and available trait-relevant variation of polyploid crops, in the process facilitating novel approaches and possibilities for genomics-assisted breeding.