Genotyping-by-Sequencing SNP Identification for Crops without a Reference Genome: Using Transcriptome Based Mapping as an Alternative Strategy

Genome-wide genetic architecture for plant maturity and drought tolerance in diploid potatoes

Cultivated potato (Solanum tuberosum) is known to be highly susceptible to drought. With climate change and its frequent episodes of drought, potato growers will face increased challenges to achieving their yield goals. Currently, a high proportion of untapped potato germplasm remains within the diploid potato relatives, and the genetic architecture of the drought tolerance and maturity traits of diploid potatoes is still unknown. As such, a panel of 384 ethyl methanesulfonate-mutagenized diploid potato clones were evaluated for drought tolerance and plant maturity under field conditions. Genome-wide association studies (GWAS) were conducted to dissect the genetic architecture of the traits. The results obtained from the genetic structure analysis of the panel showed five main groups and seven subgroups. Using the Genome Association and Prediction Integrated Tool-mixed linear model GWAS statistical model, 34 and 17 significant quantitative trait nucleotides (QTNs) were found associated with maturity and drought traits, respectively. Chromosome 5 carried most of the QTNs, some of which were also detected by using the restricted two-stage multi-locus multi-allele-GWAS haploblock-based model, and two QTNs were found to be pleiotropic for both maturity and drought traits. Using the non-parametric U-test, one and three QTNs, with 5.13%-7.4% phenotypic variations explained, showed favorable allelic effects that increase the maturity and drought trait values. The quantitaive trait loci (QTLs)/QTNs associated with maturity and drought trait were found co-located in narrow (0.5-1 kb) genomic regions with 56 candidate genes playing roles in plant development and senescence and in abiotic stress responses. A total of 127 potato clones were found to be late maturing and tolerant to drought, while nine were early to moderate-late maturing and tolerant to drought. Taken together, the data show that the studied germplasm panel and the identified candidate genes are prime genetic resources for breeders and biologists in conventional breeding and targeted gene editing as climate adaptation tools.

Genotyping by Sequencing (GBS) for Genome-Wide SNP Identification in Plants

Marker-assisted selection has played a pivotal role in developing several elite varieties in the past two decades. Molecular markers employed in plant breeding programs have recently shifted from microsatellites or simple sequence repeats (SSRs) to single nucleotide polymorphisms (SNPs) due to the ubiquity of SNP markers in the genome and the availability of various high-throughput SNP genotyping platforms. Rapid advances in sequencing technologies and the reduction in sequencing cost have facilitated SNP discovery in several plant species including non-model organisms with little or no genomic resources. Despite the lower cost of sequencing, genome complexity reduction approaches are still useful for SNP identification because many applications do not require every base of the genome to be sequenced. Genotyping-by-sequencing (GBS) is a quick and affordable reduced representation method that can simultaneously identify and genotype a large number of SNPs that has been successfully applied to a wide range of plant species. This chapter describes a robust two-enzyme GBS method for SNP discovery and genotyping that has been verified in non-model plant species.

A genome-wide association study of the occurrence of genetic variations in Edwardsiella piscicida, Vibrio harveyi, and Streptococcus parauberis under stressed environments

Bacterial mutation and genetic diversity in aquaculture have led to increasing phenotypic variances, which can weaken or invalidate strategies for controlling diseases. However, few studies have monitored the degree of mutation in fish bacterial pathogens caused by environmental pressure within a short period. In this study, transcriptomic sequences from Edwardsiella piscicida, Vibrio harveyi and Streptococcus parauberis under stressed environments were used for investigating the emergence of variants. In detail, a sub-inhibitory concentration of formalin and phenol for E. piscicida, sea water at 30°C for V. harveyi and flounder serum for S. parauberis were used as stressed environments, and significant single-nucleotide polymorphisms (SNPs) and/or mutation sites were investigated after culture in the ordinary liquid media (control) and the stressed environment through a genome-wide association study. As results, several SNPs or mutations during incubation were observed under different environments in E. piscicida and/or V. harveyi in the genes relevant to flagella, fimbria type 3 secretion systems, and outer and inner membranes that have been directly exposed to external environments. In particular, given that flagella and fimbriae are considered important factors in differentiating the serotypes in some bacterial pathogens, it can be speculated that different environmental pressures are the source of phenotypic or serotypic differentiation from the same origin. On the other hands, S. parauberis did not exhibit notable changes for 4 h when inoculated in the serum from olive flounder. The results presented in this study provide examples of possible molecular evolution in pathogens relevant to the aquaculture industry as a response to different environmental pressure.

Intra- and interspecific diversity analyses in the genus Eremurus in Iran using genotyping-by-sequencing reveal geographic population structure

Eremurus species, better known as 'Foxtail Lily' or 'Desert Candle', are important worldwide in landscaping and the cut-flower industry. One of the centers of highest diversity of the genus Eremurus is Iran, which has seven species. However, little is known about the genetic diversity within the genus Eremurus. With the advent of genotyping-by-sequencing (GBS), it is possible to develop and employ single nucleotide polymorphism (SNP) markers in a cost-efficient manner in any species, regardless of its ploidy level, genome size or availability of a reference genome. Population structure and phylogeographic analyses of the genus Eremurus in Iran using a minimum of 3002 SNP markers identified either at the genus level or at the species level from GBS data showed longitudinal geographic structuring at the country scale for the genus and for the species E. spectabilis and E. luteus, and at the regional scale for E. olgae. Our analyses furthermore showed a close genetic relatedness between E. olgae and E. stenophyllus to the extent that they should be considered subspecies within an E. olgae/stenophyllus species complex. Their close genetic relatedness may explain why crosses between these two (sub)species have been found in the wild and are exploited extensively as ornamentals. Last, current species identification, while robust, relies on flower morphology. A subset of seven SNPs with species-specific (private) alleles were selected that differentiate the seven Eremurus species. The markers will be especially useful for cultivar protection and in hybrid production, where true hybrids could be identified at the seedling stage.

Development of species-specific molecular markers in Vanilla for seedling selection of hybrids

Vanilla planifolia is the primary botanical source of vanilla extract used globally in various foods and beverages. V. planifolia has a global distribution based on a few foundational clones and therefore has limited genetic diversity. Many Vanilla species easily hybridize with V. planifolia and could be a source of valuable genetic traits like increased vanillin content, disease resistance, or early flowering. While breeding Vanilla hybrids may improve plant performance, basic molecular tools for this species are lacking. DNA-based molecular markers are the most efficient method to validate hybrid progeny, detect hybrids in commercial plantings, and identify unknown accessions. This study used publicly available sequence data to develop species-specific, qRT-PCR-based molecular markers for Vanilla. Over 580,000 assembled sequence fragments were filtered for species specificity and twenty-two targets were selected for qRT-PCR screening. Ten targets differentially amplified among V. planifolia, V. pompona, V. phaeantha, and V. palmarum with ΔCT values as high as 17.58 between species. The ten targets were used to validate the parentage of hybrid progeny from controlled crosses with most hybrid progeny showing amplification patterns similar to both parents. The ten targets were also used to screen sixteen Vanilla species for specificity, and supported species assignments for unknown accessions including the detection of putative hybrids. This is the first report using species-specific, qRT-PCR-based molecular markers in Vanilla. These markers are inexpensive, simple to develop, and can rapidly screen large populations. These methods will enable the further development of species-specific molecular markers when creating Vanilla interspecific hybrid populations.

Genome-Wide SNP discovery and genomic characterization in avocado (Persea americana Mill.)

Modern crop breeding is based on the use of genetically and phenotypically diverse plant material and, consequently, a proper understanding of population structure and genetic diversity is essential for the effective development of breeding programs. An example is avocado, a woody perennial fruit crop native to Mesoamerica with an increasing popularity worldwide. Despite its commercial success, there are important gaps in the molecular tools available to support on-going avocado breeding programs. In order to fill this gap, in this study, an avocado 'Hass' draft assembly was developed and used as reference to study 71 avocado accessions which represent the three traditionally recognized avocado horticultural races or subspecies (Mexican, Guatemalan and West Indian). An average of 5.72 M reads per individual and a total of 7,108 single nucleotide polymorphism (SNP) markers were produced for the 71 accessions analyzed. These molecular markers were used in a study of genetic diversity and population structure. The results broadly separate the accessions studied according to their botanical race in four main groups: Mexican, Guatemalan, West Indian and an additional group of Guatemalan × Mexican hybrids. The high number of SNP markers developed in this study will be a useful genomic resource for the avocado community.

Transcriptome analysis and annotation: SNPs identified from single copy annotated unigenes of three polyploid blueberry crops

Blueberry is a kind of new rising popular perennial fruit with high healthful quality. It is of utmost importance to develop new blueberry varieties for different climatic zones to satisfy the demand of people in the world. Molecular marker assisted breeding is believed to be an ideal method for the development of new blueberry varieties for its shorter breeding cycle than the conventional breeding. Simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs) markers are widely used molecular tools for marker assisted breeding, which could be detected at large scale by the transcriptome sequencing. Here, we sequenced the leaves transcriptome of 19 rabbiteye (Vaccinium ashei Reade), 13 southern highbush (Vaccinium. corymbosum L × native southern Vaccinium Spp) and 22 cultivars of northern highbush blueberry (Vaccinium corymbosum L) by using next generation sequencing technologies. A total of 80.825 Gb clean data with an average of about 12.525 million reads per cultivar were obtained. We assembled 58,968, 55,973 and 53,887 unigenes by using the clean data from rabbiteye, southern highbush and northern highbush blueberry cultivars, respectively. Among these unigenes, 3599, 3495 and 3513 unigenes were detected as candidate resistance genes in three blueberry crops. Moreover, we identified more than 8756, 9020, and 9198 SSR markers from these unigenes, and 7665, 4861, 13,063 SNPs from the annotated single copy unigenes, respectively. The results will be helpful for the molecular genetics and association analysis of blueberry and the basic molecular information of pest and disease resistance of blueberry, and would also offer huge number of molecular tools for the marker assisted breeding to produce blueberry cultivars with different adaptive characteristics.

New insights into adaptation and population structure of cork oak using genotyping by sequencing

Species respond to global climatic changes in a local context. Understanding this process, including its speed and intensity, is paramount due to the pace at which such changes are currently occurring. Tree species are particularly interesting to study in this regard due to their long generation times, sedentarism, and ecological and economic importance. Quercus suber L. is an evergreen forest tree species of the Fagaceae family with an essentially Western Mediterranean distribution. Despite frequent assessments of the species' evolutionary history, large-scale genetic studies have mostly relied on plastidial markers, whereas nuclear markers have been used on studies with locally focused sampling strategies. In this work, "Genotyping by sequencing" is used to derive 1,996 single nucleotide polymorphism markers to assess the species' evolutionary history from a nuclear DNA perspective, gain insights into how local adaptation is shaping the species' genetic background, and to forecast how Q. suber may respond to global climatic changes from a genetic perspective. Results reveal (a) an essentially unstructured species, where (b) a balance between gene flow and local adaptation keeps the species' gene pool somewhat homogeneous across its distribution, but still allowing (c) variation clines for the individuals to cope with local conditions. "Risk of Non-Adaptedness" (RONA) analyses suggest that for the considered variables and most sampled locations, (d) the cork oak should not require large shifts in allele frequencies to survive the predicted climatic changes. Future directions include integrating these results with ecological niche modeling perspectives, improving the RONA methodology, and expanding its use to other species. With the implementation presented in this work, the RONA can now also be easily assessed for other organisms.

SNP discovery of Korean short day onion inbred lines using double digest restriction site-associated DNA sequencing

Onion (Allium cepa L.) is an economically important vegetable crop around the world. Genetic and genomic research into various onion accessions will provide insights into the onion genome to enhance breeding strategies and improve crops. However, the onion’s large genome size means that studies of molecular markers are limited in onion. This study aimed to discover high quality single nucleotide polymorphisms (SNPs) from 192 onion inbred lines relating to short-day cultivation in Korea. Paired-end (PE) double digested restriction site-associated DNA sequencing (ddRAD-seq) was used to discover SNPs in onion. A total of 538,973,706 reads (25.9 GB), with an average of 2,658,491 high-quality reads, were generated using ddRAD-seq. With stringent filtering, 1904 SNPs were discovered based on onion reference scaffolds. Further, population structure and genetic relationship studies suggested that two well-differentiated sub-populations exist in onion lines. SNP-associated flanking sequences were also compared with a public non-redundant database for gene ontology and pathway analysis. To our knowledge, this is the first report to identify high-quality SNPs in onion based on reference sequences using the ddRAD-seq platform. The SNP markers identified will be useful for breeders and the research community to deepen their understanding, enhance breeding programs, and support the management of onion genomic resources.

Genome scan reveals selection acting on genes linked to stress response in wild pearl millet

Uncovering genomic regions involved in adaption is a major goal in evolutionary biology. High-throughput sequencing now makes it possible to tackle this challenge in nonmodel species. Yet, despite the increasing number of methods targeted to specifically detect genomic footprints of selection, the complex demography of natural populations often causes high rates of false positive in gene discoveries. The aim of this study was to identify climate adaptations in wild pearl millet populations, Cenchrus americanus ssp. monodii. We focused on two climate gradients, one in Mali and one in Niger. We used a two-step strategy to limit false-positive outliers. First, we considered gradients as biological replicates and performed RNA sequencing of four populations at the extremities. We combined four methods-three based on differentiation among populations and one based on diversity patterns within populations-to identify outlier SNPs from a set of 87 218 high-quality SNPs. Among 11 155 contigs of pearl millet reference transcriptome, 540 exhibited selection signals as evidenced by at least one of the four methods. In a second step, we genotyped 762 samples in 11 additional populations distributed along the gradients using SNPs from the detected contigs and random SNPs as control. We further assessed selection on this large data set using a differentiation-based method and a method based on correlations with environmental variables based. Four contigs displayed consistent signatures between the four extreme and 11 additional populations, two of which were linked to abiotic and biotic stress responses.

Comparative Transcriptome Analysis Revealed Genes Commonly Responsive to Varied Nitrate Stress in Leaves of Tibetan Hulless Barley

Nitrogen (N) deprivation or excess can lead to dramatic phenotype change, disrupt important biological processes, and ultimately limit plant productivity. To explore genes in Tibetan hulless barley responsive to varied N stress, we utilized a comparative transcriptomics method to investigate gene expression patterns under three nitrate treatments. The transcriptome of the control (optimal-nitrate, ON) sample was compared with that of free-nitrate (FN), low-nitrate (LN), and high-nitrate (HN) treatment samples, identifying 2428, 1274, and 1861 genes, respectively, that exhibited significant differences in transcript abundance. Among these, 9 genes encoding ribulose bisphosphate carboxylases exhibited up-regulated expression under varied N stress. We further compared FN vs. ON and LN vs. ON to investigate the impact of stress degree on gene expression. With the aggravation of stress, more genes were differentially expressed and thus possibly involved in the response to nitrogen deficiency. Cluster and functional enrichment analysis indicated that the differentially expressed genes (DEGs) in FN were highly enriched in response to stress, defense response, and gene expression regulation. Comprehensive comparison analysis further suggested that Tibetan hulless barley could respond to varied N stress by regulating multiple common biological processes and pathways such as nitrogen metabolism, carbon metabolism, and photosynthesis. A large number of specific DEGs involved in diverse biological processes were also detected, implying differences in the potential regulatory patterns of low- and high-N stress response. Notably, we also identified some NIN-like proteins and other transcription factors significantly modulated by these stresses, suggesting the involvement of these transcription factors in N stress response. To our knowledge, this study is the first investigation of the Tibetan hulless barley transcriptome under N stress. The identified N-stress-related genes may provide resources for genetic improvement and promote N use efficiency.