Genotyping-by-sequencing approaches to characterize crop genomes: choosing the right tool for the right application

The efficient physiological strategy of a novel tomato genotype to adapt to chronic combined water and heat stress

Climate change is increasing the frequency of high temperature shocks and water shortages, pointing to the need to develop novel tolerant varieties and to understand the mechanisms employed to withstand combined abiotic stresses. Two tomato genotypes, a heat-tolerant Solanum lycopersicum accession (LA3120) and a novel genotype (E42), previously selected as a stable yielding genotype under high temperatures, were exposed to single and combined water and heat stress. Plant functional traits, pollen viability and physiological (leaf gas exchange and chlorophyll a fluorescence emission measurements) and biochemical (antioxidant content and antioxidant enzyme activity) measurements were carried out. A Reduced Representation Sequencing approach allowed exploration of the genetic variability of both genotypes to identify candidate genes that could regulate stress responses. Both abiotic stresses had a severe impact on plant growth parameters and on the reproductive phase of development. Growth parameters and leaf gas exchange measurements revealed that the two genotypes used different physiological strategies to overcome individual and combined stresses, with E42 having a more efficient capacity to utilize the limiting water resources. Activation of antioxidant defence mechanisms seemed to be critical for both genotypes to counteract combined abiotic stresses. Candidate genes were identified that could explain the different physiological responses to stress observed in E42 compared with LA3120. Results here obtained have shown how new tomato genetic resources can be a valuable source of traits for adaptation to combined abiotic stresses and should be used in breeding programmes to improve stress tolerance in commercial varieties.

Genetic Diversity and Population Structure Analysis of the USDA Olive Germplasm Using Genotyping-By-Sequencing (GBS)

Olives are one of the most important fruit and woody oil trees cultivated in many parts of the world. Olive oil is a critical component of the Mediterranean diet due to its importance in heart health. Olives are believed to have been brought to the United States from the Mediterranean countries in the 18th century. Despite the increase in demand and production areas, only a few selected olive varieties are grown in most traditional or new growing regions in the US. By understanding the genetic background, new sources of genetic diversity can be incorporated into the olive breeding programs to develop regionally adapted varieties for the US market. This study aimed to explore the genetic diversity and population structure of 90 olive accessions from the USDA repository along with six popular varieties using genotyping-by-sequencing (GBS)-generated SNP markers. After quality filtering, 54,075 SNP markers were retained for the genetic diversity analysis. The average gene diversity (GD) and polymorphic information content (PIC) values of the SNPs were 0.244 and 0.206, respectively, indicating a moderate genetic diversity for the US olive germplasm evaluated in this study. The structure analysis showed that the USDA collection was distributed across seven subpopulations; 63% of the accessions were grouped into an identifiable subpopulation. The phylogenetic and principal coordinate analysis (PCoA) showed that the subpopulations did not align with the geographical origins or climatic zones. An analysis of the molecular variance revealed that the major genetic variation sources were within populations. These findings provide critical information for future olive breeding programs to select genetically distant parents and facilitate future gene identification using genome-wide association studies (GWAS) or a marker-assisted selection (MAS) to develop varieties suited to production in the US.

Population dynamics linked to glacial cycles in Cercis chuniana F. P. Metcalf (Fabaceae) endemic to the montane regions of subtropical China

The mountains of subtropical China are an excellent system for investigating the processes driving the geographical distribution of biodiversity and radiation of plant populations in response to Pleistocene climate fluctuations. How the major mountain ranges in subtropical China have affected the evolution of plant species in the subtropical evergreen broadleaved forest is an issue with long-term concern. Here, we focused on Cercis chuniana, a woody species endemic to the southern mountain ranges in subtropical China, to elucidate its population dynamics. We used genotyping by sequencing (GBS) to investigate the spatial pattern of genetic variation among 11 populations. Geographical isolation was detected between the populations located in adjacent mountain ranges, thought to function as geographical barriers due to their complex physiography. Bayesian time estimation revealed that population divergence occurred in the middle Pleistocene, when populations in the Nanling Mts. separated from those to the east. The orientation and physiography of the mountain ranges of subtropical China appear to have contributed to the geographical pattern of genetic variation between the eastern and western populations of C. chuniana. Complex physiography plus long-term stable ecological conditions across glacial cycles facilitated the demographic expansion in the Nanling Mts., from which contemporary migration began. The Nanling Mts. are thus considered as a suitable area for preserving population diversity and large population sizes of C. chuniana compared with other regions. As inferred by ecological niche modeling and coalescent simulations, secondary contact occurred during the warm Lushan-Tali Interglacial period, with intensified East Asia summer monsoon and continuous habitat available for occupation. Our data support the strong influence of both climatic history and topographic characteristics on the high regional phytodiversity of the subtropical evergreen broadleaved forest in subtropical China.

Breeding Canola (Brassica napus L.) for Protein in Feed and Food

Interest in canola (Brassica napus L.). In response to this interest, scientists have been tasked with altering and optimizing the protein production chain to ensure canola proteins are safe for consumption and economical to produce. Specifically, the role of plant breeders in developing suitable varieties with the necessary protein profiles is crucial to this interdisciplinary endeavour. In this article, we aim to provide an overarching review of the canola protein chain from the perspective of a plant breeder, spanning from the genetic regulation of seed storage proteins in the crop to advancements of novel breeding technologies and their application in improving protein quality in canola. A review on the current uses of canola meal in animal husbandry is presented to underscore potential limitations for the consumption of canola meal in mammals. General discussions on the allergenic potential of canola proteins and the regulation of novel food products are provided to highlight some of the challenges that will be encountered on the road to commercialization and general acceptance of canola protein as a dietary protein source.

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.

Robust and efficient software for reference-free genomic diversity analysis of genotyping-by-sequencing data on diploid and polyploid species

Genotyping-by-sequencing (GBS) is a widely used and cost-effective technique for obtaining large numbers of genetic markers from populations by sequencing regions adjacent to restriction cut sites. Although a standard reference-based pipeline can be followed to analyse GBS reads, a reference genome is still not available for a large number of species. Hence, reference-free approaches are required to generate the genetic variability information that can be obtained from a GBS experiment. Unfortunately, available tools to perform de novo analysis of GBS reads face issues of usability, accuracy and performance. Furthermore, few available tools are suitable for analysing data sets from polyploid species. In this manuscript, we describe a novel algorithm to perform reference-free variant detection and genotyping from GBS reads. Nonexact searches on a dynamic hash table of consensus sequences allow for efficient read clustering and sorting. This algorithm was integrated in the Next Generation Sequencing Experience Platform (NGSEP) to integrate the state-of-the-art variant detector already implemented in this tool. We performed benchmark experiments with three different empirical data sets of plants and animals with different population structures and ploidies, and sequenced with different GBS protocols at different read depths. These experiments show that NGSEP has comparable and in some cases better accuracy and always better computational efficiency compared to existing solutions. We expect that this new development will be useful for many research groups conducting population genetic studies in a wide variety of species.

Genome-wide association mapping of lentil (Lens culinaris Medikus) prebiotic carbohydrates toward improved human health and crop stress tolerance

Lentil, a cool-season food legume, is rich in protein and micronutrients with a range of prebiotic carbohydrates, such as raffinose-family oligosaccharides (RFOs), fructooligosaccharides (FOSs), sugar alcohols (SAs), and resistant starch (RS), which contribute to lentil's health benefits. Beneficial microorganisms ferment prebiotic carbohydrates in the colon, which impart health benefits to the consumer. In addition, these carbohydrates are vital to lentil plant health associated with carbon transport, storage, and abiotic stress tolerance. Thus, lentil prebiotic carbohydrates are a potential nutritional breeding target for increasing crop resilience to climate change with increased global nutritional security. This study phenotyped a total of 143 accessions for prebiotic carbohydrates. A genome-wide association study (GWAS) was then performed to identify associated variants and neighboring candidate genes. All carbohydrates analyzed had broad-sense heritability estimates (H2) ranging from 0.22 to 0.44, comparable to those reported in the literature. Concentration ranges corresponded to percent recommended daily allowances of 2-9% SAs, 7-31% RFOs, 51-111% RS, and 57-116% total prebiotic carbohydrates. Significant SNPs and associated genes were identified for numerous traits, including a galactosyltransferase (Lcu.2RBY.1g019390) known to aid in RFO synthesis. Further studies in multiple field locations are necessary. Yet, these findings suggest the potential for molecular-assisted breeding for prebiotic carbohydrates in lentil to support human health and crop resilience to increase global food security.

Candidate Rlm6 resistance genes against Leptosphaeria. maculans identified through a genome-wide association study in Brassica juncea (L.) Czern

One hundred and sixty-seven B. juncea varieties were genotyped on the 90K Brassica assay (42,914 SNPs), which led to the identification of sixteen candidate genes for Rlm6. Brassica species are at high risk of severe crop loss due to pathogens, especially Leptosphaeria maculans (the causal agent of blackleg). Brassica juncea (L.) Czern is an important germplasm resource for canola improvement, due to its good agronomic traits, such as heat and drought tolerance and high blackleg resistance. The present study is the first using genome-wide association studies to identify candidate genes for blackleg resistance in B. juncea based on genome-wide SNPs obtained from the Illumina Infinium 90 K Brassica SNP array. The verification of Rlm6 in B. juncea was performed through a cotyledon infection test. Genotyping 42,914 single nucleotide polymorphisms (SNPs) in a panel of 167 B. juncea lines revealed a total of seven SNPs significantly associated with Rlm6 on chromosomes A07 and B04 in B. juncea. Furthermore, 16 candidate Rlm6 genes were found in these regions, defined as nucleotide binding site leucine-rich-repeat (NLR), leucine-rich repeat RLK (LRR-RLK) and LRR-RLP genes. This study will give insights into the blackleg resistance in B. juncea and facilitate identification of functional blackleg resistance genes which can be used in Brassica breeding.

Genotyping-by-Sequencing of the regional Pacific abalone (Haliotis discus) genomes reveals population structures and patterns of gene flow

Continuous monitoring of the present genetic status is essential to preserve the genetic resource of wild populations. In this study, we sequenced regional Pacific abalone Haliotis discus samples from three different locations around the Korean peninsula to assess population structure, utilizing Genotyping-by-Sequencing (GBS) method. Using PstI enzyme for genome reduction, we demonstrated the resultant library represented the whole genome region with even spacing, and as a result 16,603 single nucleotide variants (SNVs) were produced. Genetic diversity and population structure were investigated using several methods, and a strong genetic heterogeneity was observed in the Korean abalone populations. Additionally, by comparison of the variant sets among population groups, we were able to discover 26 Korean abalone population-specific SNVs, potentially associated with phenotype differences. This is the first study demonstrating the feasibility of GBS for population genetic study on H. discus. Our results will provide valuable data for the genetic conservation and management of wild abalone populations in Korea and help future GBS studies on the marine mollusks.

Genome resequencing reveals DNA polymorphisms associated with seed size/weight determination in chickpea

Crop productivity in legumes is determined by number and size/weight of seeds. To understand the genetic basis of seed size/weight in chickpea, we performed genome resequencing of 13 small- and 5 large-seeded genotypes using Illumina platform. Single nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) differentiating small- and large-seeded genotypes were identified. A total of 17,902 SNPs and 2594 InDels located in promoter and/or coding regions that may contribute to seed size/weight were detected. Of these, 266 SNPs showed significant association with seed size/weight trait. Twenty-three genes including those involved in cell growth/division, encoding transcription factors and located within QTLs associated with seed size/weight harbored SNPs within transcription factor binding motif(s) and/or coding region. The non-synonymous SNPs were found to affect the mutational sensitivity and stability of the encoded proteins. Overall, we provided a high-quality SNP map for large-scale genotyping applications and identified candidate genes that determine seed size/weight in chickpea.

Mapping Quantitative Trait Loci Associated With Graft (In)Compatibility in Apricot (Prunus armeniaca L.)

Graft incompatibility (GI) between the most popular Prunus rootstocks and apricot cultivars is one of the major problems for rootstock usage and improvement. Failure in producing long-leaving healthy grafts greatly affects the range of available Prunus rootstocks for apricot cultivation. Despite recent advances related to the molecular mechanisms of a graft-union formation between rootstock and scion, information on genetic control of this trait in woody plants is essentially missing because of a lack of hybrid crosses, segregating for the trait. In this study, we have employed the next-generation sequencing technology to generate the single-nucleotide polymorphism (SNP) markers and construct parental linkage maps for an apricot F₁ population "Moniqui (Mo)" × "Paviot (Pa)" segregating for ability to form successful grafts with universal Prunus rootstock "Marianna 2624". To localize genomic regions associated with this trait, we genotyped 138 individuals from the "Mo × Pa" cross and constructed medium-saturated genetic maps. The female "Mo" and male "Pa" maps were composed of 557 and 501 SNPs and organized in eight linkage groups that covered 780.2 and 690.4 cM of genetic distance, respectively. Parental maps were aligned to the Prunus persica v2.0 genome and revealed a high colinearity with the Prunus reference map. Two-year phenotypic data for characters associated with unsuccessful grafting such as necrotic line (NL), bark and wood discontinuities (BD and WD), and an overall estimate of graft (in)compatibility (GI) were collected for mapping quantitative trait loci (QTLs) on both parental maps. On the map of the graft-compatible parent "Pa", two genomic regions on LG5 (44.9-60.8 cM) and LG8 (33.2-39.2 cM) were associated with graft (in)compatibility characters at different significance level, depending on phenotypic dataset. Of these, the LG8 QTL interval was most consistent between the years and supported by two significant and two putative QTLs. To our best knowledge, this is the first report on QTLs for graft (in)compatibility in woody plants. Results of this work will provide a valuable genomic resource for apricot breeding programs and facilitate future efforts focused on candidate genes discovery for graft (in)compatibility in apricot and other Prunus species.

Whole-genome sequence diversity and association analysis of 198 soybean accessions in mini-core collections

We performed whole-genome Illumina resequencing of 198 accessions to examine the genetic diversity and facilitate the use of soybean genetic resources and identified 10 million single nucleotide polymorphisms and 2.8 million small indels. Furthermore, PacBio resequencing of 10 accessions was performed, and a total of 2,033 structure variants were identified. Genetic diversity and structure analysis congregated the 198 accessions into three subgroups (Primitive, World, and Japan) and showed the possibility of a long and relatively isolated history of cultivated soybean in Japan. Additionally, the skewed regional distribution of variants in the genome, such as higher structural variations on the R gene clusters in the Japan group, suggested the possibility of selective sweeps during domestication or breeding. A genome-wide association study identified both known and novel causal variants on the genes controlling the flowering period. Novel candidate causal variants were also found on genes related to the seed coat colour by aligning together with Illumina and PacBio reads. The genomic sequences and variants obtained in this study have immense potential to provide information for soybean breeding and genetic studies that may uncover novel alleles or genes involved in agronomically important traits.

High-Throughput Genotyping Technologies in Plant Taxonomy

Molecular markers provide researchers with a powerful tool for variation analysis between plant genomes. They are heritable and widely distributed across the genome and for this reason have many applications in plant taxonomy and genotyping. Over the last decade, molecular marker technology has developed rapidly and is now a crucial component for genetic linkage analysis, trait mapping, diversity analysis, and association studies. This chapter focuses on molecular marker discovery, its application, and future perspectives for plant genotyping through pangenome assemblies. Included are descriptions of automated methods for genome and sequence distance estimation, genome contaminant analysis in sequence reads, genome structural variation, and SNP discovery methods.

Analysis of the genetic diversity of Russian common oat varieties using alleles of avenin-coding loci

Alleles of avenin-coding loci were identified in 24 varieties of common oat origin from 6 different breeding centers of the Russian Federation. It was found that 33% of the studied varieties are homogeneous. Heterogeneous varieties contain from 2 to 10 biotypes. The analysis revealed 67 different genotypes, ten of which are found in the spectra of several varieties at once. Despite the presence of matching genotypes, all the studied varieties differ in the genetic formulas of avenin, the number and frequency of occurrence of biotypes. This makes it possible to effectively identify and distinguish oat varieties sown by Russian breeding using the method of prolamin electrophoresis. The analyzed groups of varieties are characterized by high values of genetic (0.54-0.79) and intra-population diversity, which indicates the stability of populations. The exception is varieties originating from the Irkutsk region, low values of genetic and intra-population diversity in the population of which may indicate the process of genetic erosion. Analysis of the frequency distribution of alleles of avenin-coding loci does not allow to determine whether oat varieties belong to certain breeding centers, which is associated with the introduction of the same genotypes into the breeding programs of different regions.

Genome-Wide Association Studies Identifying Multiple Loci Associated With Alfalfa Forage Quality

Autotetraploid alfalfa is a major hay crop planted all over the world due to its adaptation in different environments and high quality for animal feed. However, the genetic basis of alfalfa quality is not fully understood. In this study, a diverse panel of 200 alfalfa accessions were planted in field trials using augmented experimental design at three locations in 2018 and 2019. Thirty-four quality traits were evaluated by Near Infrared Reflectance Spectroscopy (NIRS). The plants were genotyped using a genotyping by sequencing (GBS) approach and over 46,000 single nucleotide polymorphisms (SNPs) were obtained after variant calling and filtering. Genome-wide association studies (GWAS) identified 28 SNP markers associated with 16 quality traits. Among them, most of the markers were associated with fiber digestibility and protein content. Phenotypic variations were analyzed from three locations and different sets of markers were identified by GWAS when using phenotypic data from different locations, indicating that alfalfa quality traits were also affected by environmental factors. Among different sets of markers identified by location, two markers were associated with nine traits of fiber digestibility. One marker associated with lignin content was identified consistently in multiple environments. Putative candidate genes underlying fiber-related loci were identified and they are involved in the lignin and cell wall biosynthesis. The DNA markers and associated genes identified in this study will be useful for the genetic improvement of forage quality in alfalfa after the validation of the markers.

Genotyping-by-Sequencing Technology in Plant Taxonomy and Phylogeny

Genotyping-by-sequencing (GBS) is a method to discover and genotype simultaneous genome-wide high-throughput single nucleotide polymorphisms (SNPs). GBS is based on reducing genome complexity with restriction enzymes. Here we describe a method developed by Elshire et al. for constructing simplified GBS libraries and recent bioinformatic approaches developed to analyze the large volume of polymorphism data generated by this method. GBS approach is suitable for population studies, taxonomic and phylogenic studies, germplasm characterization, and breeding and trait mapping for a wide range of organisms, including plants with complex genomes.

A High-Density SNP Genetic Map Construction Using ddRAD-Seq and Mapping of Capsule Shattering Trait in Sesame

The seed-bearing capsule of sesame shatters at harvest. This wildish trait makes the crop unsuitable for mechanized harvesting and also restricts its commercial potential by limiting the cultivation for countries that have no access to low-cost labor. Therefore, the underlying genetic basis of the capsule shattering trait is highly important in order to develop mechanization-ready varieties for sustainable sesame farming. In the present study, we generated a sesame F₂ population derived from a cross between a capsule shattering cultivar (Muganli-57) and a non-shattering mutant (PI 599446), which was used to construct a genetic map based on double-digest restriction-site-associated DNA sequencing. The resulting high-density genetic map contained 782 single-nucleotide polymorphisms (SNPs) and spanned a length of 697.3 cM, with an average marker interval of 0.89 cM. Based on the reference genome, the capsule shattering trait was mapped onto SNP marker S8_5062843 (78.9 cM) near the distal end of LG8 (chromosome 8). In order to reveal genes potentially controlling the shattering trait, the marker region (S8_5062843) was examined, and a candidate gene including six CDSs was identified. Annotation showed that the gene encodes a protein with 440 amino acids, sharing ∼99% homology with transcription repressor KAN1. Compared with the capsule shattering allele, the SNP change and altered splicing in the flanking region of S8_5062843 caused a frameshift mutation in the mRNA, resulting in the loss of function of this gene in the mutant parent and thus in non-shattering capsules and leaf curling. With the use of genomic data, InDel and CAPS markers were developed to differentiate shattering and non-shattering capsule genotypes in marker-assisted selection studies. The obtained results in the study can be beneficial in breeding programs to improve the shattering trait and enhance sesame productivity.

Whole-genome resequencing reveals aberrant autosomal SNPs affect chicken feathering rate

Previous studies have shown that the feather growth rate of chicks is determined by two alleles located on the sex chromosome Z; however, in chicken production, feathering is usually not consistently controlled by the sex chromosome. To identify whether the feathering rate is related to autosomal inheritance, whole-genome resequencing was performed in eight chickens with slow- and fast-feathering rate. A total of 54,984 autosomal single nucleotide polymorphisms (SNPs) were identified, including 393 and 376 exonic SNPs in slow-feathering and fast-feathering chickens, respectively. Mutated genes were mainly involved in response to stimuli and growth and reproduction processes. Mutated genes related to slow-feathering rate were mainly involved in wingless-type MMTV integration site signaling pathway and mitogen-activated protein kinase signaling pathway, whereas mutated genes associated with fast-feathering rate were primarily enriched in autophagy, calcium signaling pathway, extracellular matrix-receptor interaction, and Focal adhesion processes. Importantly, two SNPs, involved in feather development, were found in the exonic regions of Wnt signaling genes. These results shed new light on the relationship between genetic mutation and feather growth rate from the perspective of autosomal inheritance and may have economic significance in chicken breeding.

Quick and efficient approach to develop genomic resources in orphan species: Application in Lavandula angustifolia

Next-Generation Sequencing (NGS) technologies, by reducing the cost and increasing the throughput of sequencing, have opened doors to generate genomic data in a range of previously poorly studied species. In this study, we propose a method for the rapid development of a large-scale molecular resources for orphan species. We studied as an example the true lavender (Lavandula angustifolia Mill.), a perennial sub-shrub plant native from the Mediterranean region and whose essential oil have numerous applications in cosmetics, pharmaceuticals, and alternative medicines. The heterozygous clone “Maillette” was used as a reference for DNA and RNA sequencing. We first built a reference Unigene, compound of coding sequences, thanks to de novo RNA-seq assembly. Then, we reconstructed the complete genes sequences (with introns and exons) using an Unigene-guided DNA-seq assembly approach. This aimed to maximize the possibilities of finding polymorphism between genetically close individuals despite the lack of a reference genome. Finally, we used these resources for SNP mining within a collection of 16 commercial lavender clones and tested the SNP within the scope of a genetic distance analysis. We obtained a cleaned reference of 8, 030 functionally in silico annotated genes. We found 359K polymorphic sites and observed a high SNP frequency (mean of 1 SNP per 90 bp) and a high level of heterozygosity (more than 60% of heterozygous SNP per genotype). On overall, we found similar genetic distances between pairs of clones, which is probably related to the out-crossing nature of the species and the restricted area of cultivation. The proposed method is transferable to other orphan species, requires little bioinformatics resources and can be realized within a year. This is also the first reported large-scale SNP development on Lavandula angustifolia. All the genomics resources developed herein are publicly available and provide a rich pool of molecular resources to explore and exploit lavender genetic diversity in breeding programs.

GRAS-Di system facilitates high-density genetic map construction and QTL identification in recombinant inbred lines of the wheat progenitor Aegilops tauschii

Due to large and complex genomes of Triticeae species, skim sequencing approaches have cost and analytical advantages for detecting genetic markers and building linkage maps. Here, we develop a high-density linkage map and identify quantitative trait loci (QTLs) for recombinant inbred lines of Aegilops tauschii, a D-genome donor of bread wheat, using the recently developed genotyping by Random Amplicon Sequencing-Direct (GRAS-Di) system, which facilitates skimming of the large and complicated genome and generates a large number of genetic markers. The deduced linkage groups based on the GRAS-Di genetic markers corresponded to the chromosome number of Ae. tauschii. We successfully identified stable QTLs for flowering time and spikelet shape-related traits. Genotype differences of RILs at the QTL-linked markers were significantly associated with the trait variations. In particular, one of the QTL-linked markers for flowering time was mapped close to VRN3 (also known as FLOWERING LOCUS T), which controls flowering. The GRAS-Di system is, therefore, an efficient and useful application for genotyping and linkage mapping in species with large and complex genomes, such as Triticeae species.

Linkage mapping and QTL analysis of flowering time using ddRAD sequencing with genotype error correction in Brassica napus

BACKGROUND

Brassica napus is an important oilseed crop cultivated worldwide. During domestication and breeding of B. napus, flowering time has been a target of selection because of its substantial impact on yield. Here we use double digest restriction-site associated DNA sequencing (ddRAD) to investigate the genetic basis of flowering in B. napus. An F2 mapping population was derived from a cross between an early-flowering spring type and a late-flowering winter type.

RESULTS

Flowering time in the mapping population differed by up to 25 days between individuals. High genotype error rates persisted after initial quality controls, as suggested by a genotype discordance of ~ 12% between biological sequencing replicates. After genotype error correction, a linkage map spanning 3981.31 cM and compromising 14,630 single nucleotide polymorphisms (SNPs) was constructed. A quantitative trait locus (QTL) on chromosome C2 was detected, covering eight flowering time genes including FLC.

CONCLUSIONS

These findings demonstrate the effectiveness of the ddRAD approach to sample the B. napus genome. Our results also suggest that ddRAD genotype error rates can be higher than expected in F2 populations. Quality filtering and genotype correction and imputation can substantially reduce these error rates and allow effective linkage mapping and QTL analysis.

Genome-wide SNP discovery from Darjeeling tea cultivars - their functional impacts and application toward population structure and trait associations

Genotyping by sequencing and identification of functionally relevant nucleotide variations in crop accessions are the key steps to unravel genetic control of desirable traits. Elite cultivars of Darjeeling tea were undergone SNP genotyping by double-digest restriction-site associated DNA sequencing method. This study reports a set of 54,206 high-quality SNP markers discovered from ~10.4 GB sequence data, encompassing 15 chromosomes of the reference tea genome. Genetic relatedness among the accessions conforms to the analyses of Bayesian clustering, UPGMA, and PCoA methods. Genomic positions of the discovered SNPs and their putative effect on annotated genes designated a thoughtful understanding of their functional aspects in tea system biology. A group of 95 genes was identified to be affected by high impact variants. Genome-wide association analyses of 21 agronomic and biochemical phenotypes resulted in trait-linked polymorphic loci with strong confidence (p < 0.05 and 0.001).

Innovative transcriptome-based genotyping highlights environmentally responsive genes for phenology, growth and yield in a non-model grain legume

The narrow-leafed lupin, Lupinus angustifolius L., is a grain legume crop, cultivated both as a green manure and as a source of protein for animal feed and human food production. During its domestication process, numerous agronomic traits were improved, however, only two trait-related genes were identified hitherto, both by linkage mapping. Genome-wide association studies (GWAS), exploiting genomic sequencing, did not select any novel candidate gene. In the present study, an innovative method of 3'-end reduced representation transcriptomic profiling, a massive analysis of cDNA ends, has been used for genotyping of 126 L. angustifolius lines surveyed by field phenotyping. Significant genotype × environment interactions were identified for all phenology and yield traits analysed. Principal component analysis of population structure evidenced European domestication bottlenecks, visualized by clustering of breeding materials and cultivars. GWAS provided contribution towards deciphering vernalization pathway in legumes, and, apart from highlighting known domestication loci (Ku/Julius and mol), designated novel candidate genes for L. angustifolius traits. Early phenology was associated with genes from vernalization, cold-responsiveness and phosphatidylinositol signalling pathways whereas high yield with genes controlling photosynthesis performance and abiotic stress (drought or heat) tolerance. PCR-based toolbox was developed and validated to enable tracking desired alleles in marker-assisted selection. Narrow-leafed lupin was genotyped with an innovative method of transcriptome profiling and phenotyped for phenology, growth and yield traits in field. Early phenology was found associated with genes from cold-response, vernalization and phosphatidylinositol signalling pathways, whereas high yield with genes running photosystem II and drought or heat stress response. Key loci were supplied with PCR-based toolbox for marker-assisted selection.

Understanding Host-Pathogen Interactions in Brassica napus in the Omics Era

Brassica napus (canola/oilseed rape/rapeseed) is an economically important crop, mostly found in temperate and sub-tropical regions, that is cultivated widely for its edible oil. Major diseases of Brassica crops such as Blackleg, Clubroot, Sclerotinia Stem Rot, Downy Mildew, Alternaria Leaf Spot and White Rust have caused significant yield and economic losses in rapeseed-producing countries worldwide, exacerbated by global climate change, and, if not remedied effectively, will threaten global food security. To gain further insights into the host-pathogen interactions in relation to Brassica diseases, it is critical that we review current knowledge in this area and discuss how omics technologies can offer promising results and help to push boundaries in our understanding of the resistance mechanisms. Omics technologies, such as genomics, proteomics, transcriptomics and metabolomics approaches, allow us to understand the host and pathogen, as well as the interaction between the two species at a deeper level. With these integrated data in multi-omics and systems biology, we are able to breed high-quality disease-resistant Brassica crops in a more holistic, targeted and accurate way.

An extended KASP-SNP resource for molecular breeding in Chinese cabbage(Brassica rapa L. ssp. pekinensis)

Kompetitive allele-specific PCR (KASP) is a cost-effective single-step SNP genotyping technology, With an objective to enhance the marker repertoire and develop high efficient KASP-SNP markers in Chinese cabbage, we re-sequenced four Chinese cabbage doubled haploid (DH) lines, Y177-47, Y635-10, Y510-1 and Y510-9, and generated a total of more than 38.5 billion clean base pairs. A total of 827,720 SNP loci were identified with an estimated density of 3,217 SNPs/Mb. Further, a total of 387,354 SNPs with at least 30 bp to the next most adjacent SNPs on either side were selected as resource for KASP markers. From this resource, 258 (96.27%) of 268 SNP loci were successfully transformed into KASP-SNP markers using a Roche LightCycler 480-II instrument. Among these markers, 221 (85.66%) were co-dominant markers, 220 (85.27%) were non-synonymous SNPs, and 257 (99.6%) were newly developed markers. In addition, 53 markers were applied for genotyping of 34 Brassica rapa accessions. Cluster analysis separated these 34 accessions into three clusters based on heading types. The millions of SNP loci, a large set of resource for KASP markers, as well as the newly developed KASP markers in this study may facilitate further genetic and molecular breeding studies in Brassica rapa.

Applications of genotyping-by-sequencing (GBS) in maize genetics and breeding

Genotyping-by-Sequencing (GBS) is a low-cost, high-throughput genotyping method that relies on restriction enzymes to reduce genome complexity. GBS is being widely used for various genetic and breeding applications. In the present study, 2240 individuals from eight maize populations, including two association populations (AM), backcross first generation (BC1), BC1F2, F2, double haploid (DH), intermated B73 × Mo17 (IBM), and a recombinant inbred line (RIL) population, were genotyped using GBS. A total of 955,120 of raw data for SNPs was obtained for each individual, with an average genotyping error of 0.70%. The rate of missing genotypic data for these SNPs was related to the level of multiplex sequencing: ~ 25% missing data for 96-plex and ~ 55% for 384-plex. Imputation can greatly reduce the rate of missing genotypes to 12.65% and 3.72% for AM populations and bi-parental populations, respectively, although it increases total genotyping error. For analysis of genetic diversity and linkage mapping, unimputed data with a low rate of genotyping error is beneficial, whereas, for association mapping, imputed data would result in higher marker density and would improve map resolution. Because imputation does not influence the prediction accuracy, both unimputed and imputed data can be used for genomic prediction. In summary, GBS is a versatile and efficient SNP discovery approach for homozygous materials and can be effectively applied for various purposes in maize genetics and breeding.

The Use of Genetic and Gene Technologies in Shaping Modern Rapeseed Cultivars (Brassica napus L.)

Since their domestication, Brassica oilseed species have undergone progressive transformation allied with the development of breeding and molecular technologies. The canola (Brassica napus) crop has rapidly expanded globally in the last 30 years with intensive innovations in canola varieties, providing for a wider range of markets apart from the food industry. The breeding efforts of B. napus, the main source of canola oil and canola meal, have been mainly focused on improving seed yield, oil quality, and meal quality along with disease resistance, abiotic stress tolerance, and herbicide resistance. The revolution in genetics and gene technologies, including genetic mapping, molecular markers, genomic tools, and gene technology, especially gene editing tools, has allowed an understanding of the complex genetic makeup and gene functions in the major bioprocesses of the Brassicales, especially Brassica oil crops. Here, we provide an overview on the contributions of these technologies in improving the major traits of B. napus and discuss their potential use to accomplish new improvement targets.

An enrichment method for mapping ambiguous reads to the reference genome for NGS analysis

Mapping short reads to a reference genome is an essential step in many next-generation sequencing (NGS) analyses. In plants with large genomes, a large fraction of the reads can align to multiple locations of the genome with equally good alignment scores. How to map these ambiguous reads to the genome is a challenging problem with big impacts on the downstream analysis. Traditionally, the default method is to assign an ambiguous read randomly to one of the many potential locations. In this study, we explore two alternative methods that are based on the hypothesis that the possibility of an ambiguous read being generated by a location is proportional to the total number of reads produced by that location: (1) the enrichment method that assigns an ambiguous read to the location that has produced the most reads among all the potential locations, (2) the probability method that assigns an ambiguous read to a location based on a probability proportional to the number of reads the location produces. We systematically compared the performance of the proposed methods with that of the default random method. Our results showed that the enrichment method produced better results than the default random method and the probability method in the discovery of single nucleotide polymorphisms (SNPs). Not only did it produce more SNP markers, but it also produced SNP markers with better quality, which was demonstrated using multiple mainstay genomic analyses, including genome-wide association studies (GWAS), minor allele distribution, population structure, and genomic prediction.

Skim-Sequencing Based Genotyping Reveals Genetic Divergence of the Wild and Domesticated Population of Black Tiger Shrimp (Penaeus monodon) in the Indo-Pacific Region

The domestication of a wild-caught aquatic animal is an evolutionary process, which results in genetic discrimination at the genomic level in response to strong artificial selection. Although black tiger shrimp (Penaeus monodon) is one of the most commercially important aquaculture species, a systematic assessment of genetic divergence and structure of wild-caught and domesticated broodstock populations of the species is yet to be documented. Therefore, we used skim sequencing (SkimSeq) based genotyping approach to investigate the genetic structure of 50 broodstock individuals of P. monodon species, collected from five sampling sites (n = 10 in each site) across their distribution in Indo-Pacific regions. The wild-caught P. monodon broodstock population were collected from Malaysia (MS) and Japan (MJ), while domesticated broodstock populations were collected from Madagascar (MMD), Hawaii, HI, USA (MMO), and Thailand (MT). After various filtering process, a total of 194,259 single nucleotide polymorphism (SNP) loci were identified, in which 4983 SNP loci were identified as putatively adaptive by the pcadapt approach. In both datasets, pairwise F_ST estimates high genetic divergence between wild and domesticated broodstock populations. Consistently, different spatial clustering analyses in both datasets categorized divergent genetic structure into two clusters: (1) wild-caught populations (MS and MJ), and (2) domesticated populations (MMD, MMO and MT). Among 4983 putatively adaptive SNP loci, only 50 loci were observed to be in the coding region. The gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses suggested that non-synonymous mutated genes might be associated with the energy production, metabolic functions, respiration regulation and developmental rates, which likely act to promote adaptation to the strong artificial selection during the domestication process. This study has demonstrated the applicability of SkimSeq in a highly duplicated genome of P. monodon specifically, across a range of genetic backgrounds and geographical distributions, and would be useful for future genetic improvement program of this species in aquaculture.

A high-density SNP-based linkage map using genotyping-by-sequencing and its utilization for improved genome assembly of chickpea (Cicer arietinum L.)

Genotyping-by-sequencing (GBS) allows rapid identification of markers for use in development of linkage maps, which expedite efficient breeding programs. In the present study, we have utilized GBS approach to identify and genotype single-nucleotide polymorphism (SNP) markers in an inter-specific RIL population of Cicer arietinum L. X C. reticulatum. A total of 141,639 raw SNPs were identified using the TASSEL-GBS pipeline. After stringent filtering, 8208 candidate SNPs were identified of which ~ 37% were localized in the intragenic regions followed by genic regions (~ 30%) and intergenic regions (~ 27%). We then utilized 6920 stringent selected SNPs from present study and 6714 SNPs and microsatellite markers available from previous studies for construction of linkage map. The resulting high-density linkage map comprising of eight linkage groups contained 13,590 markers which spanned 1299.14 cM of map length with an average marker density of 0.095 cM. Further, the derived linkage map was used to improve the available assembly of desi chickpea genome by anchoring 443 previously unplaced scaffolds onto eight linkage groups. The present efforts have refined anchoring of the desi chickpea genome assembly to 55.57% of the ~ 520 Mb of assembled desi genome. To the best of our knowledge, the linkage map generated in the present study represents one of the most dense linkage map developed for the crop till date. It will serve as a valuable resource for fine mapping and positional cloning of important quantitative trait loci (QTLs) associated with agronomical traits and also for anchoring and ordering of future genome sequence assemblies.

Optimization of the genotyping-by-sequencing SNP calling for diversity analysis in cape gooseberry (Physalis peruviana L.) and related taxa

A robust Genotyping-By-Sequencing (GBS) pipeline platform was examined to provide accurate discovery of Single Nucleotide Polymorphisms (SNPs) in a cape gooseberry (Physalis peruviana L.) and related taxa germplasm collection. A total of 176 accessions representing, wild, weedy, and commercial cultivars as well as related taxa from the Colombian germplasm bank and other world repositories were screened using GBS. The pipeline parameters mnLCov of 0.5 and a mnScov of 0.7, tomato and potato genomes, and cape gooseberry transcriptome for read alignments, were selected to better assess diversity and population structure in cape gooseberry and related taxa. A total of 7,425 SNPs, derived from P. peruviana common tags (unique 64 bp sequences shared between selected species), were used. Within P. peruviana, five subpopulations with a high genetic diversity and allele fixation (H_E: 0.35 to 0.36 and F_IS: -0.11 to -0.01, respectively) were detected. Conversely, low genetic differentiation (F_ST: 0.01 to 0.05) was also observed, indicating a high gene flow among subpopulations. These results contribute to the establishment of adequate conservation and breeding strategies for Cape gooseberry and closely related Physalis species.

SILEX: a fast and inexpensive high-quality DNA extraction method suitable for multiple sequencing platforms and recalcitrant plant species

BACKGROUND

The use of sequencing and genotyping platforms has undergone dramatic improvements, enabling the generation of a wealth of genomic information. Despite this progress, the availability of high-quality genomic DNA (gDNA) in sufficient concentrations is often a main limitation, especially for third-generation sequencing platforms. A variety of DNA extraction methods and commercial kits are available. However, many of these are costly and frequently give either low yield or low-quality DNA, inappropriate for next generation sequencing (NGS) platforms. Here, we describe a fast and inexpensive DNA extraction method (SILEX) applicable to a wide range of plant species and tissues.

RESULTS

SILEX is a high-throughput DNA extraction protocol, based on the standard CTAB method with a DNA silica matrix recovery, which allows obtaining NGS-quality high molecular weight genomic plant DNA free of inhibitory compounds. SILEX was compared with a standard CTAB extraction protocol and a common commercial extraction kit in a variety of species, including recalcitrant ones, from different families. In comparison with the other methods, SILEX yielded DNA in higher concentrations and of higher quality. Manual extraction of 48 samples can be done in 96 min by one person at a cost of 0.12 €/sample of reagents and consumables. Hundreds of tomato gDNA samples obtained with either SILEX or the commercial kit were successfully genotyped with Single Primer Enrichment Technology (SPET) with the Illumina HiSeq 2500 platform. Furthermore, DNA extracted from Solanum elaeagnifolium using this protocol was assessed by Pulsed-field gel electrophoresis (PFGE), obtaining a suitable size ranges for most sequencing platforms that required high-molecular-weight DNA such as Nanopore or PacBio.

CONCLUSIONS

A high-throughput, fast and inexpensive DNA extraction protocol was developed and validated for a wide variety of plants and tissues. SILEX offers an easy, scalable, efficient and inexpensive way to extract DNA for various next-generation sequencing applications including SPET and Nanopore among others.

Two Clubroot-Resistance Genes, Rcr3 and Rcr9wa, Mapped in Brassica rapa Using Bulk Segregant RNA Sequencing

Genetic resistance is widely used to manage clubroot (Plasmodiophora brassicae) in brassica crops, but new pathotypes have recently been identified on canola (Brassica napus) on the Canadian prairies. Resistance effective against both the most prevalent pathotype (3H, based on the Canadian Clubroot Differential system) and the new pathotypes is needed. BC₁ plants of Brassica rapa from a cross of line 96-6990-2 (clubroot resistance originating from turnip cultivar ‘Waaslander’) and a susceptible doubled-haploid line, ACDC, exhibited a 1:1 segregation for resistance against pathotypes 3H and 5X. A resistance gene designated as Rcr3 was mapped initially based on the percentage of polymorphic variants using bulked segregant RNA sequencing (BSR-Seq) and further mapped using Kompetitive Allele Specific PCR. DNA variants were identified by assembling short reads against a reference genome of B. rapa. Rcr3 was mapped into chromosome A08. It was flanked by single nucleotide polymorphisms (SNP) markers (A90_A08_SNP_M12 and M16) between 10.00 and 10.23 Mb, in an interval of 231.6 Kb. There were 32 genes in the Rcr3 interval. Three genes (Bra020951, Bra020974, and Bra020979) were annotated with disease resistance mechanisms, which are potential candidates for Rcr3. Another resistance gene, designated as Rcr9^wa, for resistance to pathotype 5X was mapped, with the flanking markers (A90_A08_SNP_M28 and M79) between 10.85 and 11.17 Mb using the SNP sites identified through BSR-Seq for Rcr3. There were 44 genes in the Rcr9^wa interval, three of which (Bra020827, Bra020828, Bra020814) were annotated as immune-system-process related genes, which are potential candidates for Rcr9^wa.

Fostering Conservation via an Integrated Use of Conventional Approaches and High-Throughput SPET Genotyping: A Case Study Using the Endangered Canarian Endemics Solanum lidii and S. vespertilio (Solanaceae)

Islands provide unique opportunities to integrated research approaches to study evolution and conservation because boundaries are circumscribed, geological ages are often precise, and many taxa are greatly imperiled. We combined morphological and hybridization studies with high-throughput genotyping platforms to streamline relationships in the endangered monophyletic and highly diverse lineage of Solanum in the Canarian archipelago, where three endemic taxa are currently recognized. Inter-taxa hybridizations were performed, and morphological expression was assessed with a common-garden approach. Using the eggplant Single Primer Enrichment Technology (SPET) platform with 5,093 probes, 74 individuals of three endemic taxa (Solanum lidii, S. vespertilio subsp. vespertilio, and S. vespertilio subsp. doramae) were sampled for SNPs. While morphological and breeding studies showed clear distinctions and some continuous variation, inter-taxon hybrids were fertile and heterotic for vigor traits. SPET genotyping revealed 1,421 high-quality SNPs and supported four, not three, distinct taxonomic entities associated with post-emergence geological, ecological and geographic factors of the islands. Given the lack of barriers to hybridization among all the taxa and their molecular differences, great care must be taken in population management. Conservation strategies must take account of the sexual and breeding systems and genotypic distribution among populations to successfully conserve and restore threatened/endangered island taxa, as exemplified by Solanum on the Canary Islands.

Adapting Genotyping-by-Sequencing and Variant Calling for Heterogeneous Stock Rats

The heterogeneous stock (HS) is an outbred rat population derived from eight inbred rat strains. HS rats are ideally suited for genome wide association studies; however, only a few genotyping microarrays have ever been designed for rats and none of them are currently in production. To address the need for an efficient and cost effective method of genotyping HS rats, we have adapted genotype-by-sequencing (GBS) to obtain genotype information at large numbers of single nucleotide polymorphisms (SNPs). In this paper, we have outlined the laboratory and computational steps we took to optimize double digest genotype-by-sequencing (ddGBS) for use in rats. We evaluated multiple existing computational tools and explain the workflow we have used to call and impute over 3.7 million SNPs. We have also compared various rat genetic maps, which are necessary for imputation, including a recently developed map specific to the HS. Using our approach, we obtained concordance rates of 99% with data obtained using data from a genotyping array. The principles and computational pipeline that we describe could easily be adapted for use in other species for which reliable reference genome sets are available.

Identification of genetic loci associated with forage quality in response to water deficit in autotetraploid alfalfa (Medicago sativa L.)

BACKGROUND

Alfalfa has been cultivated in many regions around the world as an important forage crop due to its nutritive value to livestock and ability to adapt to various environments. However, the genetic basis by which plasticity of quality-relevant traits influence alfalfa adaption to different water conditions remain largely unknown.

RESULTS

In the present study, 198 accessions of alfalfa of the core collection for drought tolerance were evaluated for 26 forage quality traits in a field trial under an imposed deficit irrigation gradient. Regression analysis between quality traits and water stress revealed that values of fiber-related traits were negatively correlated with values of energy-related traits as water deficit increased. More than one hundred significant markers associated with forage quality under different water treatments were identified using genome-wide association studies with genotyping by sequencing. Among them, 131 markers associated with multiple traits in all the water deficit treatments. Most of the associated markers were dependent to the levels of water deficit, suggesting genetic controls for forage quality traits were dependent to the stress treatment. Twenty-four loci associated with forage quality were annotated to functional genes that may play roles in cell development or in response to water stress.

CONCLUSIONS

This study addressed the genetic base of phenotypic variation of forage quality traits under water deficit. The SNP markers identified in this study will be useful in marker-assisted selection for the genetic improvement of alfalfa with enhanced drought tolerance while maintaining forage quality.

The development of genome-wide single nucleotide polymorphisms in blue wildebeest using the DArTseq platform

Blue wildebeest (Connochaetes taurinus taurinus) are economically important antelope that are widely utilised in the South African wildlife industry. However, very few genomic resources are available for blue wildebeest that can assist in breeding management and facilitate research. This study aimed to develop a set of genome-wide single nucleotide polymorphism (SNP) markers for blue wildebeest. The DArTseq genotyping platform, commonly used in polyploid plant species, was selected for SNP discovery. A limited number of published articles have described the use of the DArTseq platform in animals and, therefore, this study also provided a unique opportunity to assess the performance of the DArTseq platform in an animal species. A total of 20,563 SNPs, each located within a 69 bp sequence, were generated. The developed SNP markers had a high average scoring reproducibility (>99%) and a low percentage missing data (~9.21%) compared to other reduced representation sequencing approaches that have been used in animal studies. Furthermore, the number of candidate SNPs per nucleotide position decreased towards the 3' end of sequence reads, and the ratio of transitions (Ts) to transversions (Tv) remained similar for each read position. These observations indicate that there was no read position bias, such as the identification of false SNPs due to low sequencing quality, towards the tail-end of sequencing reads. The DArTseq platform was also successful in identifying a large number of informative SNPs with desirable polymorphism parameters such as a high minor allele frequency (MAF). The Bos taurus genome was used for the in silico mapping of the marker sequences and a total of 6020 (29.28%) sequences were successfully mapped against the bovine genome. The marker sequences mapped to all of the bovine chromosomes establishing the genome-wide distribution of the SNPs. Moreover, the high observed Ts:Tv ratio (2.84:1) indicate that the DArTseq platform targeted gene-rich regions of the blue wildebeest genome. Finally, functional annotation of the marker sequences revealed a wide range of different putative functions indicating that these SNP markers can be useful in functional gene studies. The DArTseq platform, therefore, represents a high-throughput, robust and cost-effective genotyping platform, which may find adoption in several other African antelope and animal species.

Construction of a high density linkage map in Oil Palm using SPET markers

A high-density genetic linkage map from a controlled cross of two oil palm (Elaeis guineensis) genotypes was constructed based on Single Primer Enrichment Technology (SPET) markers. A 5K panel of hybridization probes were used for this purpose which was derived from previously developed SNP primers in oil palm. Initially, 13,384 SNPs were detected which were reduced to 13,073 SNPs after filtering for only bi-allelic SNP. Around 75% of the markers were found to be monomorphic in the progeny, reducing the markers left for linkage mapping to 3,501. Using Lep-MAP3 software, a linkage map was constructed which contained initially 2,388 markers and had a total length of 1,370 cM. In many cases several adjacent SNP were located on the same locus, due to missing recombination events between them, leading to a total of 1,054 loci on the 16 LG. Nevertheless, the marker density of 1.74 markers per cM (0.57 cM/marker) should allow the detection of QTLs in the future. This study shows that cost efficient SPET markers are suitable for linkage map construction in oil palm and probably, also in other species.

Role of NGS and SNP genotyping methods in sugarcane improvement programs

Sugarcane (Saccharum spp.) is one of the most economically significant crops because of its high sucrose content and it is a promising biomass feedstock for biofuel production. Sugarcane genome sequencing and analysis is a difficult task due to its heterozygosity and polyploidy. Long sequence read technologies, PacBio Single-Molecule Real-Time (SMRT) sequencing, the Illumina TruSeq, and the Oxford Nanopore sequencing could solve the problem of genome assembly. On the applications side, next generation sequencing (NGS) technologies played a major role in the discovery of single nucleotide polymorphism (SNP) and the development of low to high throughput genotyping platforms. The two mainstream high throughput genotyping platforms are the SNP microarray and genotyping by sequencing (GBS). This paper reviews the NGS in sugarcane genomics, genotyping methodologies, and the choice of these methods. Array-based SNP genotyping is robust, provides consistent SNPs, and relatively easier downstream data analysis. The GBS method identifies large scale SNPs across the germplasm. A combination of targeted GBS and array-based genotyping methods should be used to increase the accuracy of genomic selection and marker-assisted breeding.

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.

The Wheat 660K SNP array demonstrates great potential for marker-assisted selection in polyploid wheat

The rapid development and application of molecular marker assays have facilitated genomic selection and genome-wide linkage and association studies in wheat breeding. Although PCR-based markers (e.g. simple sequence repeats and functional markers) and genotyping by sequencing have contributed greatly to gene discovery and marker-assisted selection, the release of a more accurate and complete bread wheat reference genome has resulted in the design of single-nucleotide polymorphism (SNP) arrays based on different densities or application targets. Here, we evaluated seven types of wheat SNP arrays in terms of their SNP number, distribution, density, associated genes, heterozygosity and application. The results suggested that the Wheat 660K SNP array contained the highest percentage (99.05%) of genome-specific SNPs with reliable physical positions. SNP density analysis indicated that the SNPs were almost evenly distributed across the whole genome. In addition, 229 266 SNPs in the Wheat 660K SNP array were located in 66 834 annotated gene or promoter intervals. The annotated genes revealed by the Wheat 660K SNP array almost covered all genes revealed by the Wheat 35K (97.44%), 55K (99.73%), 90K (86.9%) and 820K (85.3%) SNP arrays. Therefore, the Wheat 660K SNP array could act as a substitute for other 6 arrays and shows promise for a wide range of possible applications. In summary, the Wheat 660K SNP array is reliable and cost-effective and may be the best choice for targeted genotyping and marker-assisted selection in wheat genetic improvement.

MBKbase for rice: an integrated omics knowledgebase for molecular breeding in rice

To date, large amounts of genomic and phenotypic data have been accumulated in the fields of crop genetics and genomic research, and the data are increasing very quickly. However, the bottleneck to using big data in breeding is integrating the data and developing tools for revealing the relationship between genotypes and phenotypes. Here, we report a rice sub-database of an integrated omics knowledgebase (MBKbase-rice, www.mbkbase.org/rice), which integrates rice germplasm information, multiple reference genomes with a united set of gene loci, population sequencing data, phenotypic data, known alleles and gene expression data. In addition to basic data search functions, MBKbase provides advanced web tools for genotype searches at the population level and for visually displaying the relationship between genotypes and phenotypes. Furthermore, the database also provides online tools for comparing two samples by their genotypes and finding target germplasms by genotype or phenotype information, as well as for analyzing the user submitted SNP or sequence data to find important alleles in the germplasm. A soybean sub-database is planned for release in 3 months and wheat and maize will be added in 1–2 years. The data and tools integrated in MBKbase will facilitate research in crop functional genomics and molecular breeding.

Breeding Maize Maternal Haploid Inducers

Maize doubled haploid (DH) lines are usually created in vivo, through crosses with maternal haploid inducers. These inducers have the inherent ability of generating seeds with haploid embryos when used to pollinate other genotypes. The resulting haploid plants are treated with a doubling agent and self-pollinated, producing completely homozygous seeds. This rapid method of inbred line production reduces the length of breeding cycles and, consequently, increases genetic gain. Such advantages explain the wide adoption of this technique by large, well-established maize breeding programs. However, a slower rate of adoption was observed in medium to small-scale breeding programs. The high price and/or lack of environmental adaptation of inducers available for licensing, or the poor performance of those free of cost, might explain why smaller operations did not take full advantage of this technique. The lack of adapted inducers is especially felt in tropical countries, where inducer breeding efforts are more recent. Therefore, defining optimal breeding approaches for inducer development could benefit many breeding programs which are in the process of adopting the DH technique. In this manuscript, we review traits important to maize maternal haploid inducers, explain their genetic basis, listing known genes and quantitative trait loci (QTL), and discuss different breeding approaches for inducer development. The performance of haploid inducers has an important impact on the cost of DH line production.

Dynamics of the genetic diversity of oat varieties in the Tyumen region at avenin-coding loci

Molecular and biochemical markers are used to analyze the intraspecific genetic diversity of crops. Prolamin-coding loci are highly effective for assessing this indicator. On the basis of the Laboratory of Varietal Seed Identification of the State Agrarian University of the Northern Trans-Urals, 18 varieties of common oat included in the State Register of Selection Achievements in the Tyumen Region from the 1930s to 2019 were studied by electrophoresis in 2018–2019. The aim of the work was to study the dynamics of the genetic diversity of oat varieties at avenin-coding loci. For the analysis, 100 grains of each variety were used. Electrophoresis was carried out in vertical plates of 13.2 % polyacrylamide gel at a constant voltage of 500 V for 4.0–4.5 h. It was found that 44.4 % of the varieties are heterogeneous, each consisting of two biotypes. For three loci, 20 alleles were identified, 10 of which were detected for the first time. The allele frequency of avenin-coding loci varied with time. In the process of variety exchange, alleles that are characteristic of varieties of non-Russian origin were replaced by alleles present in domestic varieties and then in the varieties developed by local breeding institutions. The following alleles had the highest frequency in Tyumen varieties: Avn A4 (50.0 %), A2 (25.0 %), Avn B4 (50.0 %), Bnew6 (37.5 %), Avn C1 (37.5 %), C2 and C5 (25.0 %). These alleles are of great value as markers of agronomically and adaptively important characters for the region in question. The amount of genetic diversity of oats varied with time from 0.33 in 1929–1950 to up to 0.75 in 2019. The high value of genetic diversity in modern breeding varieties of the Scientific Research Institute of Agriculture of the Northern Trans-Urals and an increase in this indicator over the past 20 years are associated with the use of genetically heterogeneous source material in the breeding process. This allowed obtaining varieties with high adaptive potentials in the natural climatic conditions of the region.

Stacking up RADSeq assembly programs: From complete hit to completely abysmal

Decreasing sequencing costs have driven a rapid expansion of novel genotyping methods. One of these methods is the exploitation of restriction enzyme cut sites to generate genome-wide but reduced representation sequencing libraries (RRLs), alternatively termed genotyping by sequencing or restriction-site associated DNA sequencing. Without a reference genome, the resulting short sequence reads must be assembled de novo. There are many possible assembly programs, most not explicitly developed for RRL data, and we know little of their effectiveness. In this issue of Molecular Ecology Resources, LaCava et al. (2020) systematically evaluate six commonly used programs and two commonly varied parameters for complete and accurate assembly of RRLs, using simulated double digests of Homo sapiens and Arabidopsis thaliana genomes with varied mutation rates and types. The authors find substantial variation in performance across assembly programs. The most consistently high-performing assembler is infrequently used in their literature survey (CD-HIT; Li and Godzik, 2006), while several others fail to produce complete, accurate assemblies under many conditions. LaCava et al. additionally recommend best practices in parameter choice and evaluation of future assembly programs-advice that molecular ecologists working to assemble sequences of all kinds should take to heart.

Multiregional origins of the domesticated tetraploid wheats

We used genotyping-by-sequencing (GBS) to investigate the evolutionary history of domesticated tetraploid wheats. With a panel of 189 wild and domesticated wheats, we identified 1,172,469 single nucleotide polymorphisms (SNPs) with a read depth ≥3. Principal component analyses (PCAs) separated the Triticum turgidum and Triticum timopheevii accessions, as well as wild T. turgidum from the domesticated emmers and the naked wheats, showing that SNP typing by GBS is capable of providing robust information on the genetic relationships between wheat species and subspecies. The PCAs and a neighbour-joining analysis suggested that domesticated tetraploid wheats have closest affinity with wild emmers from the northern Fertile Crescent, consistent with the results of previous genetic studies on the origins of domesticated wheat. However, a more detailed examination of admixture and allele sharing between domesticates and different wild populations, along with genome-wide association studies (GWAS), showed that the domesticated tetraploid wheats have also received a substantial genetic input from wild emmers from the southern Levant. Taking account of archaeological evidence that tetraploid wheats were first cultivated in the southern Levant, we suggest that a pre-domesticated crop spread from this region to southeast Turkey and became mixed with a wild emmer population from the northern Fertile Crescent. Fixation of the domestication traits in this mixed population would account for the allele sharing and GWAS results that we report. We also propose that feralization of the component of the pre-domesticated population that did not acquire domestication traits has resulted in the modern wild population from southeast Turkey displaying features of both the domesticates and wild emmer from the southern Levant, and hence appearing to be the sole progenitor of domesticated tetraploids when the phylogenetic relationships are studied by methods that assume a treelike pattern of evolution.

Enhancing Genetic Gain through Genomic Selection: From Livestock to Plants

Although long-term genetic gain has been achieved through increasing use of modern breeding methods and technologies, the rate of genetic gain needs to be accelerated to meet humanity's demand for agricultural products. In this regard, genomic selection (GS) has been considered most promising for genetic improvement of the complex traits controlled by many genes each with minor effects. Livestock scientists pioneered GS application largely due to livestock's significantly higher individual values and the greater reduction in generation interval that can be achieved in GS. Large-scale application of GS in plants can be achieved by refining field management to improve heritability estimation and prediction accuracy and developing optimum GS models with the consideration of genotype-by-environment interaction and non-additive effects, along with significant cost reduction. Moreover, it would be more effective to integrate GS with other breeding tools and platforms for accelerating the breeding process and thereby further enhancing genetic gain. In addition, establishing an open-source breeding network and developing transdisciplinary approaches would be essential in enhancing breeding efficiency for small- and medium-sized enterprises and agricultural research systems in developing countries. New strategies centered on GS for enhancing genetic gain need to be developed.