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

A Pipeline and Recommendations for Population and Individual Diagnostic SNP Selection in Non-Model Species

Despite substantial reductions in the cost of sequencing over the last decade, genetic panels remain relevant due to their cost-effectiveness and flexibility across a variety of sample types. In particular, single nucleotide polymorphism (SNP) panels are increasingly favoured for conservation applications. SNP panels are often used because of their adaptability, effectiveness with low-quality samples, and cost-efficiency for population monitoring and forensics. However, the selection of diagnostic SNPs for population assignment and individual identification can be challenging. The consequences of poor SNP selection are under-powered panels, inaccurate results, and monetary loss. Here, we develop a novel and user-friendly SNP selection pipeline (mPCRselect) that can be used to select SNPs for population assignment and/or individual identification. mPCRselect allows any researcher, who has sufficient SNP-level data, to design a successful and cost-effective SNP panel for a diploid species of conservation concern.

Morpho-Phenological, Chemical, and Genetic Characterization of Italian Maize Landraces from the Lazio Region

In the framework of a Collaboration Agreement between CREA and ARSIAL, a morpho-phenological, chemical, and genetic characterization of maize populations native to the Lazio region was carried out. During 2022 and 2023, a set of 50 accessions, belonging both to ARSIAL and CREA maize collections, were multiplied in Bergamo. Morpho-phenological descriptors were recorded in the field: plant height, ear height, and male and female flowering time. The grain chemical composition in terms of protein, lipid, starch, ash and fiber was evaluated by near-infrared spectroscopy (NIRS). A double-digest restriction-site-associated DNA sequencing (ddRADseq) strategy was used to genotype the landraces. The two collections were not significantly different in terms of grain chemical composition. On the other hand, the ARSIAL and CREA germplasm showed a different distribution in the three cluster-based population structure obtained by ddRADseq, which largely corresponded to the distribution map of their collection sites. The materials from the Lazio region maintained by ARSIAL and CREA were revealed to be different. The comparison between the two groups of landraces showed the importance of characterizing germplasm collections to promote the recovery and valorization of local biodiversity.

Fine-scale adaptive divergence and population genetic structure of Aedes aegypti in Metropolitan Manila, Philippines

BACKGROUND

The adaptive divergence of Aedes aegypti populations to heterogeneous environments can be a driving force behind the recent expansion of their habitat distribution and outbreaks of dengue disease in urbanized areas. In this study, we investigated the population genomics of Ae. aegypti at a regional scale in Metropolitan Manila, Philippines.

METHODS

We used the Pool-Seq double digestion restriction-site association DNA sequencing (ddRAD-Seq) approach to generate a high number of single nucleotide polymorphisms (SNPs), with the aim to determine local adaptation and compare the population structure with 11 microsatellite markers. A total of 217 Ae. aegypti individuals from seven female and seven male populations collected from Metropolitan Manila were used in the assays.

RESULTS

We detected 65,473 SNPs across the populations, of which 76 were non-neutral SNPs. Of these non-neutral SNPs, the multivariate regression test associated 50 with eight landscape variables (e.g. open space, forest, etc.) and 29 with five climate variables (e.g. air temperature, humidity, etc.) (P-value range 0.005-0.045) in female and male populations separately. Male and female populations exhibited contrasting spatial divergence, with males exhibiting greater divergence than females, most likely reflecting the different dispersal abilities of male and female mosquitoes. In the comparative analysis of the same Ae. aegypti individuals, the pairwise FST values of 11 microsatellite markers were lower than those of the neutral SNPs, indicating that the neutral SNPs generated via pool ddRAD-Seq were more sensitive in terms of detecting genetic differences between populations at fine-spatial scales.

CONCLUSIONS

Overall, our study demonstrates the utility of pool ddRAD-Seq for examining genetic differences in Ae. aegypti populations in areas at fine-spatial scales that could inform vector control programs such as Wolbachia-infected mosquito mass-release programs. This in turn would provide information on mosquito population dispersal patterns and the potential barriers to mosquito movement within and around the release area. In addition, the potential of environmental adaptability observed in Ae. aegypti could help population control efforts.

Functional differentiation and genetic diversity of rice cation exchanger (CAX) genes and their potential use in rice improvement

Cation exchanger (CAX) genes play an important role in plant growth/development and response to biotic and abiotic stresses. Here, we tried to obtain important information on the functionalities and phenotypic effects of CAX gene family by systematic analyses of their expression patterns, genetic diversity (gene CDS haplotypes, structural variations, gene presence/absence variations) in 3010 rice genomes and nine parents of 496 Huanghuazhan introgression lines, the frequency shifts of the predominant gcHaps at these loci to artificial selection during modern breeding, and their association with tolerances to several abiotic stresses. Significant amounts of variation also exist in the cis-regulatory elements (CREs) of the OsCAX gene promoters in 50 high-quality rice genomes. The functional differentiation of OsCAX gene family were reflected primarily by their tissue and development specific expression patterns and in varied responses to different treatments, by unique sets of CREs in their promoters and their associations with specific agronomic traits/abiotic stress tolerances. Our results indicated that OsCAX1a and OsCAX2 as general signal transporters were in many processes of rice growth/development and responses to diverse environments, but they might be of less value in rice improvement. OsCAX1b, OsCAX1c, OsCAX3 and OsCAX4 was expected to be of potential value in rice improvement because of their associations with specific traits, responsiveness to specific abiotic stresses or phytohormones, and relatively high gcHap and CRE diversity. Our strategy was demonstrated to be highly efficient to obtain important genetic information on genes/alleles of specific gene family and can be used to systematically characterize the other rice gene families.

Evolutionary history of an Alpine Archaeognath (Machilis pallida): Insights from different variant

Reconstruction of species histories is a central aspect of evolutionary biology. Patterns of genetic variation within and among populations can be leveraged to elucidate evolutionary processes and demographic histories. However, interpreting genetic signatures and unraveling the contributing processes can be challenging, in particular for non-model organisms with complex reproductive modes and genome organization. One way forward is the combined consideration of patterns revealed by different molecular markers (nuclear vs. mitochondrial) and types of variants (common vs. rare) that differ in their age, mode, and rate of evolution. Here, we applied this approach to RNAseq data generated for Machilis pallida (Archaeognatha), an Alpine jumping bristletail considered parthenogenetic and triploid. We generated de novo transcriptome and mitochondrial assemblies to obtain high-density data to investigate patterns of mitochondrial and common and rare nuclear variation in 17 M. pallida individuals sampled from all known populations. We find that the different variant types capture distinct aspects of the evolutionary history and discuss the observed patterns in the context of parthenogenesis, polyploidy, and survival during glaciation. This study highlights the potential of different variant types to gain insights into evolutionary scenarios even from challenging but often available data and the suitability of M. pallida and the genus Machilis as a study system for the evolution of sexual strategies and polyploidization during environmental change. We also emphasize the need for further research which will be stimulated and facilitated by these newly generated resources and insights.

Genetic dissection of thousand-seed weight in linseed (Linum usitatissimum L.) using multi-locus genome-wide association study

Flaxseed/linseed is an important oilseed crop having applications in the food, nutraceutical, and paint industry. Seed weight is one of the most crucial determinants of seed yield in linseed. Here, quantitative trait nucleotides (QTNs) associated with thousand-seed weight (TSW) have been identified using multi-locus genome-wide association study (ML-GWAS). Field evaluation was carried out in five environments in multi-year-location trials. SNP genotyping information of the AM panel of 131 accessions comprising 68,925 SNPs was employed for ML-GWAS. From the six ML-GWAS methods employed, five methods helped identify a total of 84 unique significant QTNs for TSW. QTNs identified in ≥ 2 methods/environments were designated as stable QTNs. Accordingly, 30 stable QTNs have been identified for TSW accounting up to 38.65% trait variation. Alleles with positive effect on trait were analyzed for 12 strong QTNs with r ² ≥ 10.00%, which showed significant association of specific alleles with higher trait value in three or more environments. A total of 23 candidate genes have been identified for TSW, which included B3 domain-containing transcription factor, SUMO-activating enzyme, protein SCARECROW, shaggy-related protein kinase/BIN2, ANTIAUXIN-RESISTANT 3, RING-type E3 ubiquitin transferase E4, auxin response factors, WRKY transcription factor, and CBS domain-containing protein. In silico expression analysis of candidate genes was performed to validate their possible role in different stages of seed development process. The results from this study provide significant insight and elevate our understanding on genetic architecture of TSW trait in linseed.

3D-GBS: a universal genotyping-by-sequencing approach for genomic selection and other high-throughput low-cost applications in species with small to medium-sized genomes

Despite the increased efficiency of sequencing technologies and the development of reduced-representation sequencing (RRS) approaches allowing high-throughput sequencing (HTS) of multiplexed samples, the per-sample genotyping cost remains the most limiting factor in the context of large-scale studies. For example, in the context of genomic selection (GS), breeders need genome-wide markers to predict the breeding value of large cohorts of progenies, requiring the genotyping of thousands candidates. Here, we introduce 3D-GBS, an optimized GBS procedure, to provide an ultra-high-throughput and ultra-low-cost genotyping solution for species with small to medium-sized genome and illustrate its use in soybean. Using a combination of three restriction enzymes (PstI/NsiI/MspI), the portion of the genome that is captured was reduced fourfold (compared to a "standard" ApeKI-based protocol) while reducing the number of markers by only 40%. By better focusing the sequencing effort on limited set of restriction fragments, fourfold more samples can be genotyped at the same minimal depth of coverage. This GBS protocol also resulted in a lower proportion of missing data and provided a more uniform distribution of SNPs across the genome. Moreover, we investigated the optimal number of reads per sample needed to obtain an adequate number of markers for GS and QTL mapping (500-1000 markers per biparental cross). This optimization allows sequencing costs to be decreased by ~ 92% and ~ 86% for GS and QTL mapping studies, respectively, compared to previously published work. Overall, 3D-GBS represents a unique and affordable solution for applications requiring extremely high-throughput genotyping where cost remains the most limiting factor.

Identification of QTNs Associated With Flowering Time, Maturity, and Plant Height Traits in Linum usitatissimum L. Using Genome-Wide Association Study

Early flowering, maturity, and plant height are important traits for linseed to fit in rice fallows, for rainfed agriculture, and for economically viable cultivation. Here, Multi-Locus Genome-Wide Association Study (ML-GWAS) was undertaken in an association mapping panel of 131 accessions, genotyped using 68,925 SNPs identified by genotyping by sequencing approach. Phenotypic evaluation data of five environments comprising 3 years and two locations were used. GWAS was performed for three flowering time traits including days to 5%, 50%, and 95% flowering, days to maturity, and plant height by employing five ML-GWAS methods: FASTmrEMMA, FASTmrMLM, ISIS EM-BLASSO, mrMLM, and pLARmEB. A total of 335 unique QTNs have been identified for five traits across five environments. 109 QTNs were stable as observed in ≥2 methods and/or environments, explaining up to 36.6% phenotypic variance. For three flowering time traits, days to maturity, and plant height, 53, 30, and 27 stable QTNs, respectively, were identified. Candidate genes having roles in flower, pollen, embryo, seed and fruit development, and xylem/phloem histogenesis have been identified. Gene expression of candidate genes for flowering and plant height were studied using transcriptome of an early maturing variety Sharda (IC0523807). The present study unravels QTNs/candidate genes underlying complex flowering, days to maturity, and plant height traits in linseed.

Lens orientalis Contributes Quantitative Trait Loci and Candidate Genes Associated With Ascochyta Blight Resistance in Lentil

Australian lentil production is affected by several major biotic constraints including Ascochyta blight (AB), caused by Ascochyta lentis, a devastating fungal disease. Cultivation of AB resistant cultivars, alongside agronomic management including fungicide application, is the current most economically viable control strategy. However, the breakdown of AB resistance in cultivars, such as Northfield and Nipper, suggests the need for introgression of new and diverse resistance genes. Successful introgression entails an understanding of the genetic basis of resistance. In this context, a biparental mapping population derived from a cross between a recently identified AB resistant accession ILWL 180 (Lens orientalis) and a susceptible cultivar ILL 6002 was produced. A genetic linkage map was constructed from single-nucleotide polymorphism markers generated using a genotyping-by-sequencing transcript approach. Genetic dissection of the mapping population revealed a major quantitative trait loci (QTL) region nested with three QTLs on linkage group 5 and explained 9.5-11.5 percent (%) of phenotypic variance for AB resistance. Another QTL was identified on LG2 with phenotypic variance of 9.6%. The identified QTL regions harbored putative candidate genes potentially associated with defense responses to A. lentis infection. The QTL analysis and the candidate gene information are expected to contribute to the development of diagnostic markers and enable marker-assisted resistance selection in lentil breeding programmes.

Genetic control of tracheid properties in Norway spruce wood

Through the use of genome-wide association studies (GWAS) mapping it is possible to establish the genetic basis of phenotypic trait variation. Our GWAS study presents the first such effort in Norway spruce (Picea abies (L). Karst.) for the traits related to wood tracheid characteristics. The study employed an exome capture genotyping approach that generated 178 101 Single Nucleotide Polymorphisms (SNPs) from 40 018 probes within a population of 517 Norway spruce mother trees. We applied a least absolute shrinkage and selection operator (LASSO) based association mapping method using a functional multi-locus mapping approach, with a stability selection probability method as the hypothesis testing approach to determine significant Quantitative Trait Loci (QTLs). The analysis has provided 30 significant associations, the majority of which show specific expression in wood-forming tissues or high ubiquitous expression, potentially controlling tracheids dimensions, their cell wall thickness and microfibril angle. Among the most promising candidates based on our results and prior information for other species are: Picea abies BIG GRAIN 2 (PabBG2) with a predicted function in auxin transport and sensitivity, and MA_373300g0010 encoding a protein similar to wall-associated receptor kinases, which were both associated with cell wall thickness. The results demonstrate feasibility of GWAS to identify novel candidate genes controlling industrially-relevant tracheid traits in Norway spruce.

Recommendations for Choosing the Genotyping Method and Best Practices for Quality Control in Crop Genome-Wide Association Studies

High-throughput genotyping boosts genome-wide association studies (GWAS) in crop species, leading to the identification of single-nucleotide polymorphisms (SNPs) associated with economically important traits. Choosing a cost-effective genotyping method for crop GWAS requires careful examination of several aspects, namely, the purpose and the scale of the study, crop-specific genomic features, and technical and economic matters associated with each genotyping option. Once genotypic data have been obtained, quality control (QC) procedures must be applied to avoid bias and false signals in genotype–phenotype association tests. QC for human GWAS has been extensively reviewed; however, QC for crop GWAS may require different actions, depending on the GWAS population type. Here, we review most popular genotyping methods based on next-generation sequencing (NGS) and array hybridization, and report observations that should guide the investigator in the choice of the genotyping method for crop GWAS. We provide recommendations to perform QC in crop species, and deliver an overview of bioinformatics tools that can be used to accomplish all needed tasks. Overall, this work aims to provide guidelines to harmonize those procedures leading to SNP datasets ready for crop GWAS.

SECNVs: A Simulator of Copy Number Variants and Whole-Exome Sequences From Reference Genomes

Copy number variants are duplications and deletions of the genome that play an important role in phenotypic changes and human disease. Many software applications have been developed to detect copy number variants using either whole-genome sequencing or whole-exome sequencing data. However, there is poor agreement in the results from these applications. Simulated datasets containing copy number variants allow comprehensive comparisons of the operating characteristics of existing and novel copy number variant detection methods. Several software applications have been developed to simulate copy number variants and other structural variants in whole-genome sequencing data. However, none of the applications reliably simulate copy number variants in whole-exome sequencing data. We have developed and tested Simulator of Exome Copy Number Variants (SECNVs), a fast, robust and customizable software application for simulating copy number variants and whole-exome sequences from a reference genome. SECNVs is easy to install, implements a wide range of commands to customize simulations, can output multiple samples at once, and incorporates a pipeline to output rearranged genomes, short reads and BAM files in a single command. Variants generated by SECNVs are detected with high sensitivity and precision by tools commonly used to detect copy number variants. SECNVs is publicly available at https://github.com/YJulyXing/SECNVs.

ddRAD Sequencing-Based Identification of Genomic Boundaries and Permeability in Quercus ilex and Q. suber Hybrids

Hybridization and its relevance is a hot topic in ecology and evolutionary biology. Interspecific gene flow may play a key role in species adaptation to environmental change, as well as in the survival of endangered populations. Despite the fact that hybridization is quite common in plants, many hybridizing species, such as Quercus spp., maintain their integrity, while precise determination of genomic boundaries between species remains elusive. Novel high throughput sequencing techniques have opened up new perspectives in the comparative analysis of genomes and in the study of historical and current interspecific gene flow. In this work, we applied ddRADseq technique and developed an ad hoc bioinformatics pipeline for the study of ongoing hybridization between two relevant Mediterranean oaks, Q. ilex and Q. suber. We adopted a local scale approach, analyzing adult hybrids (sensu lato) identified in a mixed stand and their open-pollinated progenies. We have identified up to 9,435 markers across the genome and have estimated individual introgression levels in adults and seedlings. Estimated contribution of Q. suber to the genome is higher, on average, in hybrid progenies than in hybrid adults, suggesting preferential backcrossing with this parental species, maybe followed by selection during juvenile stages against individuals with higher Q. suber genomic contribution. Most discriminating markers seem to be scattered throughout the genome, suggesting that a large number of small genomic regions underlie boundaries between these species. A noticeable proportion of the markers (26%) showed allelic frequencies in adult hybrids very similar to one of the parental species, and very different from the other; a finding that seems relevant for understanding the hybridization process and the occurrence of adaptive introgression. Candidate marker databases developed in this study constitute a valuable resource to design large scale re-sequencing experiments in Mediterranean sclerophyllous oak species and could provide insight in species boundaries and on adaptive introgression between Q. suber and Q. ilex.

Genome-wide association study identified novel candidate loci affecting wood formation in Norway spruce

Norway spruce is a boreal forest tree species of significant ecological and economic importance. Hence there is a strong imperative to dissect the genetics underlying important wood quality traits in the species. We performed a functional genome-wide association study (GWAS) of 17 wood traits in Norway spruce using 178 101 single nucleotide polymorphisms (SNPs) generated from exome genotyping of 517 mother trees. The wood traits were defined using functional modelling of wood properties across annual growth rings. We applied a Least Absolute Shrinkage and Selection Operator (LASSO-based) association mapping method using a functional multilocus mapping approach that utilizes latent traits, with a stability selection probability method as the hypothesis testing approach to determine a significant quantitative trait locus. The analysis provided 52 significant SNPs from 39 candidate genes, including genes previously implicated in wood formation and tree growth in spruce and other species. Our study represents a multilocus GWAS for complex wood traits in Norway spruce. The results advance our understanding of the genetics influencing wood traits and identifies candidate genes for future functional studies.

Fragmentation of Pooled PCR Products for Highly Multiplexed TILLING

Improvements to massively parallel sequencing have allowed the routine recovery of natural and induced sequence variants. A broad range of biological disciplines have benefited from this, ranging from plant breeding to cancer research. The need for high sequence coverage to accurately recover single nucleotide variants and small insertions and deletions limits the applicability of whole genome approaches. This is especially true in organisms with a large genome size or for applications requiring the screening of thousands of individuals, such as the reverse-genetic technique known as TILLING. Using PCR to target and sequence chosen genomic regions provides an attractive alternative as the vast reduction in interrogated bases means that sample size can be dramatically increased through amplicon multiplexing and multi-dimensional sample pooling while maintaining suitable coverage for recovery of small mutations. Direct sequencing of PCR products is limited, however, due to limitations in read lengths of many next generation sequencers. In the present study we show the optimization and use of ultrasonication for the simultaneous fragmentation of multiplexed PCR amplicons for TILLING highly pooled samples. Sequencing performance was evaluated in a total of 32 pooled PCR products produced from 4096 chemically mutagenized Hordeum vulgare DNAs pooled in three dimensions. Evaluation of read coverage and base quality across amplicons suggests this approach is suitable for high-throughput TILLING and other applications employing highly pooled complex sampling schemes. Induced mutations previously identified in a traditional TILLING screen were recovered in this dataset further supporting the efficacy of the approach.

Extensions of BLUP Models for Genomic Prediction in Heterogeneous Populations: Application in a Diverse Switchgrass Sample

Genomic prediction is a useful tool to accelerate genetic gain in selection using DNA marker information. However, this technology typically relies on standard prediction procedures, such as genomic BLUP, that are not designed to accommodate population heterogeneity resulting from differences in marker effects across populations. In this study, we assayed different prediction procedures to capture marker-by-population interactions in genomic prediction models. Prediction procedures included genomic BLUP and two kernel-based extensions of genomic BLUP which explicitly accounted for population heterogeneity. To model population heterogeneity, dissemblance between populations was either depicted by a unique coefficient (as previously reported), or a more flexible function of genetic distance between populations (proposed herein). Models under investigation were applied in a diverse switchgrass sample under two validation schemes: whole-sample calibration, where all individuals except selection candidates are included in the calibration set, and cross-population calibration, where the target population is entirely excluded from the calibration set. First, we showed that using fixed effects, from principal components or putative population groups, appeared detrimental to prediction accuracy, especially in cross-population calibration. Then we showed that modeling population heterogeneity by our proposed procedure resulted in highly significant improvements in model fit. In such cases, gains in accuracy were often positive. These results suggest that population heterogeneity may be parsimoniously captured by kernel methods. However, in cases where improvement in model fit by our proposed procedure is null-to-moderate, ignoring heterogeneity should probably be preferred due to the robustness and simplicity of the standard genomic BLUP model.

Validation of Genotyping by Sequencing Using Transcriptomics for Diversity and Application of Genomic Selection in Tetraploid Potato

Potato is an important food crop due to its increasing consumption, and as a result, there is demand for varieties with improved production. However, the current status of breeding for improved varieties is a long process which relies heavily on phenotypic evaluation and dated molecular techniques and has little emphasis on modern genotyping approaches. Evaluation and selection before a cultivar is commercialized typically takes 10-15 years. Molecular markers have been developed for disease and pest resistance, resulting in initial marker-assisted selection in breeding. This study has evaluated and implemented a high-throughput transcriptome sequencing method for dense marker discovery in potato for the application of genomic selection. An Australian relevant collection of commercial cultivars was selected, and identification and distribution of high quality SNPs were examined using standard bioinformatic pipelines and a custom approach for the prediction of allelic dosage. As a result, a large number of SNP markers were identified and filtered to generate a high-quality subset that was then combined with historic phenotypic data to assess the approach for genomic selection. Genomic selection potential was predicted for highly heritable traits and the approach demonstrated advantages over the previously used technologies in terms of markers identified as well as costs incurred. The high-quality SNP list also provided acceptable genome coverage which demonstrates its applicability for much larger future studies. This SNP list was also annotated to provide an indication of function and will serve as a resource for the community in future studies. Genome wide marker tools will provide significant benefits for potato breeding efforts and the application of genomic selection will greatly enhance genetic progress.

Evaluation and Recommendations for Routine Genotyping Using Skim Whole Genome Re-sequencing in Canola

Whole genome sequencing offers genome wide, unbiased markers, and inexpensive library preparation. With the cost of sequencing decreasing rapidly, many plant genomes of modest size are amenable to skim whole genome resequencing (skim WGR). The use of skim WGR in diverse sample sets without the use of imputation was evaluated in silico in 149 canola samples representative of global diversity. Fastq files with an average of 10x coverage of the reference genome were used to generate skim samples representing 0.25x, 0.5x, 1x, 2x, 3x, 4x, and 5x sequencing coverage. Applying a pre-defined list of SNPs versus de novo SNP discovery was evaluated. As skim WGR is expected to result in some degree of insufficient allele sampling, all skim coverage levels were filtered at a range of minimum read depths from a relaxed minimum read depth of 2 to a stringent read depth of 5, resulting in 28 list-based SNP sets. As a broad recommendation, genotyping pre-defined SNPs between 1x and 2x coverage with relatively stringent depth filtering is appropriate for a diverse sample set of canola due to a balance between marker number, sufficient accuracy, and sequencing cost, but depends on the intended application. This was experimentally examined in two sample sets with different genetic backgrounds: 1x coverage of 1,590 individuals from 84 Australian spring type four-parent crosses aimed at maximizing diversity as well as one commercial F1 hybrid, and 2x coverage of 379 doubled haploids (DHs) derived from a subset of the four-parent crosses. To determine optimal coverage in a simpler genetic background, the DH sample sequence coverage was further down sampled in silico. The flexible and cost-effective nature of the protocol makes it highly applicable across a range of species and purposes.

Genotyping-by-sequencing identifies date palm clone preference in agronomics of the State of Qatar

Understanding the genetic diversity in a crop population is key to its targeted breeding for desired traits, such as higher yields, better fruit quality and resistance to disease and changing climates. Date fruits represent a major crop in the Middle East and are key to achieving future food independence in arid countries like Qatar. We previously determined the genome of the date palm Phoenix dactylifera and showed that date palm trees world-wide divide into two distinct subpopulations of Eastern and Western origins. Here we applied a resource of SNPs from 179 commercially available date fruits to assess the genetic diversity of date palm trees grown in the State of Qatar. We found that palm trees in Qatar are mainly of Eastern origin, and that their genetic diversity doesn’t associate with regions of the State. Together with targeted genetic assays, our resource can be used in the future for date palm cultivar identification, to aid selecting suitable cultivars for targeted breeding, to improve a country’s date palm genetic diversity, and to certify the origin of date fruits and trees.

Genotyping-by-sequencing through transcriptomics: implementation in a range of crop species with varying reproductive habits and ploidy levels

The application of genomics in crops has the ability to significantly improve genetic gain for agriculture. Many marker-dense tools have been developed, but few have seen broad adoption in plant genomics due to issues of significant variations of genome size, levels of ploidy, single nucleotide polymorphism (SNP) frequency and reproductive habit. When combined with limited breeding activities, small research communities and scant sequence resources, the suitability of popular systems is often suboptimal and routinely fails to effectively balance cost-effectiveness and sample throughput. Genotyping-by-sequencing (GBS) encompasses a range of protocols including resequencing of the transcriptome. This study describes a skim GBS-transcriptomics (GBS-t) approach developed to be broadly applicable, cost-effective and high-throughput while still assaying a significant number of SNP loci. A range of crop species with differing levels of ploidy and degree of inbreeding/outbreeding were chosen, including perennial ryegrass, a diploid outbreeding forage grass; phalaris, a putative segmental allotetraploid outbreeding forage grass; lentil, a diploid inbreeding grain legume; and canola, an allotetraploid partially outbreeding oilseed. GBS-t was validated as a simple and largely automated, cost-effective method which generates sufficient SNPs (from 89 738 to 231 977) with acceptable levels of missing data and even genome coverage from c. 3 million sequence reads per sample. GBS-t is therefore a broadly applicable system suitable for many crops, offering advantages over other systems. The correct choice of subsequent sequence analysis software is important, and the bioinformatics process should be iterative and tailored to the specific challenges posed by ploidy variation and extent of heterozygosity.

Toward Genomics-Based Breeding in C3 Cool-Season Perennial Grasses

Most important food and feed crops in the world belong to the C3 grass family. The future of food security is highly reliant on achieving genetic gains of those grasses. Conventional breeding methods have already reached a plateau for improving major crops. Genomics tools and resources have opened an avenue to explore genome-wide variability and make use of the variation for enhancing genetic gains in breeding programs. Major C3 annual cereal breeding programs are well equipped with genomic tools; however, genomic research of C3 cool-season perennial grasses is lagging behind. In this review, we discuss the currently available genomics tools and approaches useful for C3 cool-season perennial grass breeding. Along with a general review, we emphasize the discussion focusing on forage grasses that were considered orphan and have little or no genetic information available. Transcriptome sequencing and genotype-by-sequencing technology for genome-wide marker detection using next-generation sequencing (NGS) are very promising as genomics tools. Most C3 cool-season perennial grass members have no prior genetic information; thus NGS technology will enhance collinear study with other C3 model grasses like Brachypodium and rice. Transcriptomics data can be used for identification of functional genes and molecular markers, i.e., polymorphism markers and simple sequence repeats (SSRs). Genome-wide association study with NGS-based markers will facilitate marker identification for marker-assisted selection. With limited genetic information, genomic selection holds great promise to breeders for attaining maximum genetic gain of the cool-season C3 perennial grasses. Application of all these tools can ensure better genetic gains, reduce length of selection cycles, and facilitate cultivar development to meet the future demand for food and fodder.

Genomics-based strategies for the use of natural variation in the improvement of crop metabolism

Next-generation genomics holds great potential in the study of plant phenotypic variation. With several crop reference genomes now available, the affordable costs of de novo genome assembly or target resequencing offer the opportunity to mine the enormous amount of genetic diversity hidden in crop wild relatives. Wide introgressions from these wild ancestors species or land races represent a possible strategy to improve cultivated varieties. In this review, we discuss the mechanisms underlying metabolic diversity within plant species and the possible strategies (and barriers) to introgress novel metabolic traits into cultivated varieties. We show how deep genomic surveys uncover various types of structural variants from extended gene pools of major crops and highlight how this variation may be used for the improvement of crop metabolism.

Targeted capture in evolutionary and ecological genomics

The rapid expansion of next-generation sequencing has yielded a powerful array of tools to address fundamental biological questions at a scale that was inconceivable just a few years ago. Various genome-partitioning strategies to sequence select subsets of the genome have emerged as powerful alternatives to whole-genome sequencing in ecological and evolutionary genomic studies. High-throughput targeted capture is one such strategy that involves the parallel enrichment of preselected genomic regions of interest. The growing use of targeted capture demonstrates its potential power to address a range of research questions, yet these approaches have yet to expand broadly across laboratories focused on evolutionary and ecological genomics. In part, the use of targeted capture has been hindered by the logistics of capture design and implementation in species without established reference genomes. Here we aim to (i) increase the accessibility of targeted capture to researchers working in nonmodel taxa by discussing capture methods that circumvent the need of a reference genome, (ii) highlight the evolutionary and ecological applications where this approach is emerging as a powerful sequencing strategy and (iii) discuss the future of targeted capture and other genome-partitioning approaches in the light of the increasing accessibility of whole-genome sequencing. Given the practical advantages and increasing feasibility of high-throughput targeted capture, we anticipate an ongoing expansion of capture-based approaches in evolutionary and ecological research, synergistic with an expansion of whole-genome sequencing.

Single Nucleotide Polymorphism Identification in Polyploids: A Review, Example, and Recommendations

Understanding the relationship between genotype and phenotype is a major biological question and being able to predict phenotypes based on molecular genotypes is integral to molecular breeding. Whole-genome duplications have shaped the history of all flowering plants and present challenges to elucidating the relationship between genotype and phenotype, especially in neopolyploid species. Although single nucleotide polymorphisms (SNPs) have become popular tools for genetic mapping, discovery and application of SNPs in polyploids has been difficult. Here, we summarize common experimental approaches to SNP calling, highlighting recent polyploid successes. To examine the impact of software choice on these analyses, we called SNPs among five peanut genotypes using different alignment programs (BWA-mem and Bowtie 2) and variant callers (SAMtools, GATK, and Freebayes). Alignments produced by Bowtie 2 and BWA-mem and analyzed in SAMtools shared 24.5% concordant SNPs, and SAMtools, GATK, and Freebayes shared 1.4% concordant SNPs. A subsequent analysis of simulated Brassica napus chromosome 1A and 1C genotypes demonstrated that, of the three software programs, SAMtools performed with the highest sensitivity and specificity on Bowtie 2 alignments. These results, however, are likely to vary among species, and we therefore propose a series of best practices for SNP calling in polyploids.