Mining for single nucleotide polymorphisms and insertions/deletions in maize expressed sequence tag data

Epistasis and pleiotropy-induced variation for plant breeding

Epistasis refers to nonallelic interaction between genes that cause bias in estimates of genetic parameters for a phenotype with interactions of two or more genes affecting the same trait. Partitioning of epistatic effects allows true estimation of the genetic parameters affecting phenotypes. Multigenic variation plays a central role in the evolution of complex characteristics, among which pleiotropy, where a single gene affects several phenotypic characters, has a large influence. While pleiotropic interactions provide functional specificity, they increase the challenge of gene discovery and functional analysis. Overcoming pleiotropy-based phenotypic trade-offs offers potential for assisting breeding for complex traits. Modelling higher order nonallelic epistatic interaction, pleiotropy and non-pleiotropy-induced variation, and genotype × environment interaction in genomic selection may provide new paths to increase the productivity and stress tolerance for next generation of crop cultivars. Advances in statistical models, software and algorithm developments, and genomic research have facilitated dissecting the nature and extent of pleiotropy and epistasis. We overview emerging approaches to exploit positive (and avoid negative) epistatic and pleiotropic interactions in a plant breeding context, including developing avenues of artificial intelligence, novel exploitation of large-scale genomics and phenomics data, and involvement of genes with minor effects to analyse epistatic interactions and pleiotropic quantitative trait loci, including missing heritability.

Understanding genetic diversity in drought-adaptive hybrid parental lines in pearl millet

Information on genetic diversity and population structure is helpful to strategize enhancing the genetic base of hybrid parental lines in breeding programs. The present study determined the population structure and genetic diversity of 109 pearl millet hybrid parental lines, known for their better adaptation and performance in drought-prone environments, using 16,472 single nucleotide polymorphic (SNP) markers generated from GBS (genotyping-by-sequencing) platforms. The SNPs were distributed uniformly across the pearl millet genome and showed considerable genetic diversity (0.337), expected heterozygosity (0.334), and observed heterozygosity (0.031). Most of the pairs of lines (78.36%) had Identity-by-State (IBS) based genetic distances of more than 0.3, indicating a significant amount of genetic diversity among the parental lines. Bayesian model-based population stratification, neighbor-joining phylogenetic analysis, and principal coordinate analysis (PCoA) differentiated all hybrid parental lines into two clear-cut major groups, one each for seed parents (B-lines) and pollinators (R-lines). Majority of parental lines sharing common parentages were found grouped in the same cluster. Analysis of molecular variance (AMOVA) revealed 7% of the variation among subpopulations, and 93% of the variation was attributable to within sub-populations. Chromosome 3 had the highest number of LD regions. Genomic LD decay distance was 0.69 Mb and varied across the different chromosomes. Genetic diversity based on 11 agro-morphological and grain quality traits also suggested that the majority of the B- and R-lines were grouped into two major clusters with few overlaps. In addition, the combined analysis of phenotypic and genotypic data showed similarities in the population grouping patterns. The present study revealed the uniqueness of most of the inbred lines, which can be a valuable source of new alleles and help breeders to utilize these inbred lines for the development of hybrids in drought-prone environments.

Mapping of quantitative trait loci for the nutritional value of fresh market tomato

The incidence of many diseases, such as cancer, cardiovascular diseases, and diabetes, is associated with malnutrition and an unbalanced daily diet. Vegetables are an important source of vitamins and essential compounds for human health. As a result, such metabolites have increasingly become the focus of breeding programs. Tomato is one of the most popular components of our daily diet. Therefore, the improvement of tomato's nutritional quality is an important goal. In the present study, we performed targeted metabolic profiling of an interspecific Solanum pimpinellifolium × S. lycopersicum inbred backcross line (IBL) population and identified quantitative trait loci responsible for the nutritional value of tomato. Transgressive segregation was apparent for many of the nutritional compounds such that some IBLs had extremely high levels of various amino acids and vitamins compared to their parents. A total of 117 QTLs for nutritional traits including 62 QTLs for amino acids, 18 QTLs for fatty acids, 12 QTLs for water-soluble vitamins, and 25 QTLs for fat-soluble vitamins were identified. Moreover, almost 24% of identified QTLs were confirmed in previous studies, and 40 possible gene candidates were found for 18 identified QTLs. These findings can help breeders to improve the nutritional value of tomato.

Bioinformatics approaches and big data analytics opportunities in improving fisheries and aquaculture

Aquaculture has witnessed an excellent growth rate during the last two decades and offers huge potential to provide nutritional as well as livelihood security. Genomic research has contributed significantly toward the development of beneficial technologies for aquaculture. The existing high throughput technologies like next-generation technologies generate oceanic data which requires extensive analysis using appropriate tools. Bioinformatics is a rapidly evolving science that involves integrating gene based information and computational technology to produce new knowledge for the benefit of aquaculture. Bioinformatics provides new opportunities as well as challenges for information and data processing in new generation aquaculture. Rapid technical advancements have opened up a world of possibilities for using current genomics to improve aquaculture performance. Understanding the genes that govern economically relevant characteristics, necessitates a significant amount of additional research. The various dimensions of data sources includes next-generation DNA sequencing, protein sequencing, RNA sequencing gene expression profiles, metabolic pathways, molecular markers, and so on. Appropriate bioinformatics tools are developed to mine the biologically relevant and commercially useful results. The purpose of this scoping review is to present various arms of diverse bioinformatics tools with special emphasis on practical translation to the aquaculture industry.

Resequencing of 410 Sesame Accessions Identifies SINST1 as the Major Underlying Gene for Lignans Variation

Sesame is a promising oilseed crop that produces specific lignans of clinical importance. Hence, a molecular description of the regulatory mechanisms of lignan biosynthesis is essential for crop improvement. Here, we resequence 410 sesame accessions and identify 5.38 and 1.16 million SNPs (single nucleotide polymorphisms) and InDels, respectively. Population genomic analyses reveal that sesame has evolved a geographic pattern categorized into northern (NC), middle (MC), and southern (SC) groups, with potential origin in the southern region and subsequent introduction to the other regions. Selective sweeps analysis uncovers 120 and 75 significant selected genomic regions in MC and NC groups, respectively. By screening these genomic regions, we unveiled 184 common genes positively selected in these subpopulations for exploitation in sesame improvement. Genome-wide association study identifies 17 and 72 SNP loci for sesamin and sesamolin variation, respectively, and 11 candidate causative genes. The major pleiotropic SNPC/A locus for lignans variation is located in the exon of the gene SiNST1. Further analyses revealed that this locus was positively selected in higher lignan content sesame accessions, and the "C" allele is favorable for a higher accumulation of lignans. Overexpression of SiNST1C in sesame hairy roots significantly up-regulated the expression of SiMYB58, SiMYB209, SiMYB134, SiMYB276, and most of the monolignol biosynthetic genes. Consequently, the lignans content was significantly increased, and the lignin content was slightly increased. Our findings provide insights into lignans and lignin regulation in sesame and will facilitate molecular breeding of elite varieties and marker-traits association studies.

Low-Density Reference Fingerprinting SNP Dataset of CIMMYT Maize Lines for Quality Control and Genetic Diversity Analyses

CIMMYT maize lines (CMLs), which represent the tropical maize germplasm, are freely available worldwide. All currently released 615 CMLs and fourteen temperate maize inbred lines were genotyped with 180 kompetitive allele-specific PCR single nucleotide polymorphisms to develop a reference fingerprinting SNP dataset that can be used to perform quality control (QC) and genetic diversity analyses. The QC analysis identified 25 CMLs with purity, identity, or mislabeling issues. Further field observation, purification, and re-genotyping of these CMLs are required. The reference fingerprinting SNP dataset was developed for all of the currently released CMLs with 152 high-quality SNPs. The results of principal component analysis and average genetic distances between subgroups showed a clear genetic divergence between temperate and tropical maize, whereas the three tropical subgroups partially overlapped with one another. More than 99% of the pairs of CMLs had genetic distances greater than 0.30, showing their high genetic diversity, and most CMLs are distantly related. The heterotic patterns, estimated with the molecular markers, are consistent with those estimated using pedigree information in two major maize breeding programs at CIMMYT. These research findings are helpful for ensuring the regeneration and distribution of the true CMLs, via QC analysis, and for facilitating the effective utilization of the CMLs, globally.

Somatic mutations during rapid clonal domestication of Populus alba var. pyramidalis

For many clonally propagated species, the accumulation of somatic mutations is the principal driver of declines in yield and quality. However, somatic mutations may also promote genetic diversification. Thus, elucidating somatic mutation rates and patterns is important to understand the genetic basis undergirding the emergence of commercially valuable traits and developmental processes. In this study, we studied the effect of short-time clonal domestication of Populus alba var. pyramidalis, a species that has been propagated by cutting for the last 67 years. We found that: (1) the somatic mutation rate for P. alba var. pyramidalis is 9.24 × 10-9, which is higher than rates observed in related species; (2) there were more mutations near heterozygous regions, and a larger proportion of CpG and CHG sites were associated with somatic mutations, which may be related to the blocking of DNA repair by methylation; and (3) deleterious mutations were not shared by multiple individuals, and all occurred in heterozygous states, demonstrating the strong selective pressures that act against deleterious mutations. Taken together, the results of our study provide a global view of somatic mutation that will aid efforts to understand the genetic basis of commercially valuable traits and to improve clonally breeding species.

Whole genome re-sequencing of indian wheat genotypes for identification of genomic variants for grain iron and zinc content

BACKGROUND

Whole-genome sequencing information which is of abundant significance for genetic evolution, and breeding of crops. Wheat (Triticum spp) is most widely grown and consumed crops globally. Micronutrients are very essential for healthy development of human being and their sufficient consumption in diet is essential for various metabolic functions. Biofortification of wheat grains with iron (Fe) and zinc (Zn) has proved the most reliable and effective way to combat micronutrient associated deficiency. Genetic variability for grain micronutrient could provide insight to dissect the traits.

METHODS AND RESULTS

In the current study, 1300 wheat lines were screened for grain Fe and Zn content, out of which only five important Indian wheat genotypes were selected on the basis of Fe and Zn contents. These lines were multiplied during at the National Agri-Food Biotechnology Institute (NABI) and re-sequenced to identify genomic variants in candidate genes for Fe and Zn between the genotypes. Whole genome sequencing generated ̴ 12 Gb clean data. Comparative genome analysis identified 254 genomic variants in the candidate genes associated with deleterious effect on protein function.

CONCLUSIONS

The present study demonstrated the fundamental in understanding the genomic variations for Fe and Zn enrichment to generate healthier wheat grains.

Exploring Genomic Variations in Nematode-Resistant Mutant Rice Lines

Rice (Oryza sativa) production is seriously affected by the root-knot nematode Meloidogyne graminicola, which has emerged as a menace in upland and irrigated rice cultivation systems. Previously, activation tagging in rice was utilized to identify candidate gene(s) conferring resistance against M. graminicola. T-DNA insertional mutants were developed in a rice landrace (acc. JBT 36/14), and four mutant lines showed nematode resistance. Whole-genome sequencing of JBT 36/14 was done along with the four nematode resistance mutant lines to identify the structural genetic variations that might be contributing to M. graminicola resistance. Sequencing on Illumina NovaSeq 6000 platform identified 482,234 genetic variations in JBT 36/14 including 448,989 SNPs and 33,245 InDels compared to reference indica genome. In addition, 293,238-553,648 unique SNPs and 32,395-65,572 unique InDels were found in the four mutant lines compared to their JBT 36/14 background, of which 93,224 SNPs and 8,170 InDels were common between all the mutant lines. Functional annotation of genes containing these structural variations showed that the majority of them were involved in metabolism and growth. Trait analysis revealed that most of these genes were involved in morphological traits, physiological traits and stress resistance. Additionally, several families of transcription factors, such as FAR1, bHLH, and NAC, and putative susceptibility (S) genes, showed the presence of SNPs and InDels. Our results indicate that subject to further genetic validations, these structural genetic variations may be involved in conferring nematode resistance to the rice mutant lines.

Genome-wide SNP and InDel analysis of three Philippine mango species inferred from whole-genome sequencing

Background

The Philippines is among the top 10 major exporters of mango worldwide. However, genomic studies of Philippine mangoes remain largely unexplored and lacking. Here, we sequenced the whole genome of the three Philippine mango species, namely, Mangifera odorata (Huani), Mangifera altissima (Paho), and Mangifera indica “Carabao” variety using Illumina HiSeq 2500, to identify and analyze their genome-wide variants (SNPs and InDels).

Results

The high confidence variants were identified by successfully mapping 93–95% of the quality-filtered reads to the Alphonso and Tommy Atkins mango reference genomes. Using these two currently available mango genomes, most variants were observed in M. odorata (4,353,063 and 4,277,287), followed by M. altissima (3,392,763 and 3,449,917), and lastly, M. indica Carabao (2,755,267 and 2,852,480). Approximately 50, 46, and 38% of the variants were unique in the three Philippine mango genomes. The analysis of variant effects and functional annotation across the three mango species revealed 56,982 variants with high-impact effects mapped onto 37,746 genes, of which 25% were found to be novel. The affected mango genes include those with potential economic importance such as 6945 genes for defense/resistance/immune response, 323 genes for fruit development, and 338 genes for anthocyanin production.

Conclusions

To date, this is the first sequencing effort to comprehensively analyze genome-wide variants essential for the development of genome-wide markers specific to these mango species native to the Philippines. This study provides an important genomic resource that can be used for the genetic improvement of mangoes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s43141-022-00326-3.

Assessment of Genetic Diversity and Relatedness in an Andean Potato Collection from Argentina by High-Density Genotyping

Native potatoes are the most diverse among cultivated potato species and thus constitute a valuable source for identifying genes for potato improvement. Nevertheless, high-density mapping, needed to reveal allelic diversity, has not been performed for native Argentinian potatoes. We present a study of the genetic variability and population structure of 96 Andigena potatoes from Northwestern Argentina performed using a subset of 5035 SNPs with no missing data and full reproducibility. These high-density markers are distributed across the genome and present a good coverage of genomic regions. A Bayesian approach revealed the presence of: (I) a major group comprised of most of the Andean accessions; (II) a smaller group containing the out-group cv. Spunta and the sequenced genotype DM; and (III) a third group containing colored flesh potatoes. This grouping was also consistent when maximum likelihood trees were constructed and further confirmed by a principal coordinate analysis. A group of 19 accessions stored as Andean varieties clustered consistently with group Tuberosum accessions. This was in agreement with previous studies and we hypothesize that they may be reintroductions of European-bred long day-adapted potatoes. The present study constitutes a valuable source for allele mining of genes of interest and thus provides a tool for association mapping studies.

Comparative transcriptome analyses for metribuzin tolerance provide insights into key genes and mechanisms restoring photosynthetic efficiency in bread wheat (Triticum aestivum L.)

Weeds are the biggest threat to cropping system sustainability in wheat. Metribuzin is a versatile herbicide for broad-spectrum weed management. Understanding key genes, mechanisms and functional markers are essential to develop higher metribuzin tolerant wheats. We identified Chuan Mai 25 (tolerant) and Ritchie (susceptible) as contrasting genotypes to metribuzin stress through dose-response analyses. Transcriptome sequencing using NovaSeq 6000 RNA-Seq platform identified a total of 77,443 genes; 59,915 known genes and 17,528 novel genes. The functional enrichment analysis at 0 h, 24 h and 60 h herbicide exposure revealed that endogenous increase of metabolic enzymes, light-harvesting chlorophyll proteins, PSII stability factor HCF136 and glucose metabolism conferred metribuzin tolerance. The validation of DEGs using RT-qPCR and QTL mapping confirmed their responsiveness to metribuzin. Transcription factors MYB, AP2-EREBP, ABI3VP1, bHLH, NAC are significantly expressed during metribuzin stress. Transcripts with significant enrichments revealed 114 SSRs for genomic selection. The master regulators provide promising avenues for enhancing metribuzin tolerance.

Discovery of new genetic determinants of morphological plasticity in rice roots and shoots under phosphate starvation using GWAS

Phosphorus is an essential nutrient for plants that is often in short supply. In rice (Oryza sativa L.), inorganic phosphate (P_i) deficiency leads to various physiological disorders that consequently affect plant productivity. In this study, a large-scale phenotyping experiment using 160 Vietnamese rice landraces was performed under greenhouse conditions, by employing an alpha lattice design with three replicates, to identify quantitative trait loci (QTLs) associated with plant growth inhibition caused by P_i deficiency. Rice plantlets were grown for six weeks in the PVC sand column (16 cm diameter × 80 cm height) supplied with P_i-deficient medium (10 μM P) or full-P_i Yoshida medium (320 μM P). The effects of P_i deficiency on the number of crown roots, root length, shoot length, root weight, shoot weight and total weight were studied. From 36 significant markers identified using a genome-wide association study, 21 QTLs associated with plant growth inhibition under P_i starvation were defined. In total, 158 candidate genes co-located with the defined QTLs were identified. Interestingly, one QTL (qRST9.14) was associated with all three weight-traits. The co-located gene GLYCEROPHOSPHODIESTER PHOSPHODIESTERASE 13 was found to be potentially involved in P_i transport. Understanding the molecular mechanisms of Pi-starvation responses, and identifying the potential QTLs responsible for low-P_i stress tolerance, will provide valuable information for developing new varieties tolerant of low-P_i conditions.

Resequencing of global Tartary buckwheat accessions reveals multiple domestication events and key loci associated with agronomic traits

BACKGROUND

Tartary buckwheat (Fagopyrum tataricum) is a nutritionally balanced and flavonoid-rich crop plant that has been in cultivation for 4000 years and is now grown globally. Despite its nutraceutical and agricultural value, the characterization of its genetics and its domestication history is limited.

RESULTS

Here, we report a comprehensive database of Tartary buckwheat genomic variation based on whole-genome resequencing of 510 germplasms. Our analysis suggests that two independent domestication events occurred in southwestern and northern China, resulting in diverse characteristics of modern Tartary buckwheat varieties. Genome-wide association studies for important agricultural traits identify several candidate genes, including FtUFGT3 and FtAP2YT1 that significantly correlate with flavonoid accumulation and grain weight, respectively.

CONCLUSIONS

We describe the domestication history of Tartary buckwheat and provide a detailed resource of genomic variation to allow for genomic-assisted breeding in the improvement of elite cultivars.

Field Screen and Genotyping of Phaseolus vulgaris against Two Begomoviruses in Georgia, USA

Simple Summary

Snap bean (Phaseolus vulgaris) production and quality have been negatively impacted by two whitefly-transmitted begomoviruses: cucurbit leaf crumple virus (CuLCrV) and sida golden mosaic Florida virus (SiGMFV), which often appear as a mixed infection in Georgia. However, there is no information available in terms of resistance to these two viruses in commercial cultivars/genotypes. Hence, commercially available snap bean varieties/genotypes (n = 84 in 2018; n = 80 in 2019; most of the genotypes were common in both years (with a few exceptions) were screened in two field seasons of 2018 and 2019. We also included two commonly grown Lima bean (Phaseolus lunatus) varieties in our field screening. As a result of this screening, we identified twenty Phaseolus genotypes with high-to-moderate levels of resistance and twenty-one genotypes with high levels of susceptibility. While there were differences among the Phaseolus spp. in severity of viral symptoms, suggesting differential susceptibility to viruses (CuLCrV and SiGMFV) and potential field resistance, the resistance mechanism is yet to be characterized. However, based on the greenhouse evaluation with two genotypes-each (susceptible vs. resistant) exposed to viruliferous whiteflies infected with CuLCrV and SiGMFV, we observed that the susceptible genotypes accumulated higher copy numbers of both viruses and displayed severe crumple severity compared to the resistant genotypes, indicating that resistant might potentially be against the virus complex than against the whiteflies. Adult whitefly counts differed among the Phaseolus spp. in both the years, indicating variability in host preference. We further sequenced 82 genotypes (80 snap bean and two Lima bean) to unravel the variations within the genomes. Genome sequencing followed by bioinformatic analyses revealed a considerable number of sequence variants, single nucleotide polymorphisms (SNPs), and insertions and deletions (InDels) in the genomes. Considering the variations in disease response and the underlying variations in the sequenced genomes, it can be speculated that some of the phenotypic variations (against CuLCrV and SiGMFV) could be due to a high level of genomic variation in the host. Future genome-wide association studies with the identified genomic variants may shed some light on this.

Abstract

The production and quality of Phaseolusvulgaris (snap bean) have been negatively impacted by leaf crumple disease caused by two whitefly-transmitted begomoviruses: cucurbit leaf crumple virus (CuLCrV) and sida golden mosaic Florida virus (SiGMFV), which often appear as a mixed infection in Georgia. Host resistance is the most economical management strategy against whitefly-transmitted viruses. Currently, information is not available with respect to resistance to these two viruses in commercial cultivars. In two field seasons (2018 and 2019), we screened Phaseolus spp. genotypes (n = 84 in 2018; n = 80 in 2019; most of the genotypes were common in both years with a few exceptions) for resistance against CuLCrV and/or SiGMFV. We also included two commonly grown Lima bean (Phaseolus lunatus) varieties in our field screening. Twenty Phaseolus spp. genotypes with high to moderate-levels of resistance (disease severity ranging from 5%–50%) to CuLCrV and/or SiGMFV were identified. Twenty-one Phaseolus spp. genotypes were found to be highly susceptible with a disease severity of ≥66%. Furthermore, based on the greenhouse evaluation with two genotypes-each (two susceptible and two resistant; identified in field screen) exposed to viruliferous whiteflies infected with CuLCrV and SiGMFV, we observed that the susceptible genotypes accumulated higher copy numbers of both viruses and displayed severe crumple severity compared to the resistant genotypes, indicating that resistance might potentially be against the virus complex rather than against the whiteflies. Adult whitefly counts differed significantly among Phaseolus genotypes in both years. The whole genome of these Phaseolus spp. [snap bean (n = 82); Lima bean (n = 2)] genotypes was sequenced and genetic variability among them was identified. Over 900 giga-base (Gb) of filtered data were generated and >88% of the resulting data were mapped to the reference genome, and SNP and Indel variants in Phaseolus spp. genotypes were obtained. A total of 645,729 SNPs and 68,713 Indels, including 30,169 insertions and 38,543 deletions, were identified, which were distributed in 11 chromosomes with chromosome 02 harboring the maximum number of variants. This phenotypic and genotypic information will be helpful in genome-wide association studies that will aid in identifying the genetic basis of resistance to these begomoviruses in Phaseolus spp.

Prospects of next generation sequencing in lentil breeding

Lentil is an important food legume crop that has large and complex genome. During past years, considerable attention has been given on the use of next generation sequencing for enriching the genomic resources including identification of SSR and SNP markers, development of unigenes, transcripts, and identification of candidate genes for biotic and abiotic stresses, analysis of genetic diversity and identification of genes/ QTLs for agronomically important traits. However, in other crops including pulses, next generation sequencing has revolutionized the genomic research and helped in genomic assisted breeding rapidly and cost effectively. The present review discuss current status and future prospects of the use NGS based breeding in lentil.

Genome-Wide Discovery of InDel Markers in Sesame (Sesamum indicum L.) Using ddRADSeq

The development and validation of different types of molecular markers is crucial to conducting marker-assisted sesame breeding. Insertion-deletion (InDel) markers are highly polymorphic and suitable for low-cost gel-based genotyping. From this perspective, this study aimed to discover and develop InDel markers through bioinformatic analysis of double digest restriction site-associated DNA sequencing (ddRADSeq) data from 95 accessions belonging to the Mediterranean sesame core collection. Bioinformatic analysis indicated the presence of 7477 InDel positions genome wide. Deletions accounted for 61% of the InDels and short deletions (1-2 bp) were the most abundant type (94.9%). On average, InDels of at least 2 bp in length had a frequency of 2.99 InDels/Mb. The 86 InDel sites having length ≥8 bp were detected in genome-wide analysis. These regions can be used for the development of InDel markers considering low-cost genotyping with agarose gels. In order to validate these InDels, a total of 38 InDel regions were selected and primers were successfully amplified. About 13% of these InDels were in the coding sequences (CDSs) and in the 3'- and 5'- untranslated regions (UTRs). Furthermore, the efficiencies of these 16 InDel markers were assessed on 32 sesame accessions. The polymorphic information content (PIC) of these 16 markers ranged from 0.06 to 0.62 (average: 0.33). These results demonstrated the success of InDel identification and marker development for sesame with the use of ddRADSeq data. These agarose-resolvable InDel markers are expected to be useful for sesame breeders.

Transcriptional dynamics of Zn-accumulation in developing kernels of maize reveals important Zn-uptake mechanisms

In the present study, transcriptomic analysis of 10-days old baby kernels of two contrasting maize genotypes, namely VQL-2 (high kernel Zn accumulator) and CM-145 (low kernel Zn accumulator), under low- and optimum- soil Zn conditions generated 1948 differentially expressed transcripts. Among these, 666 and 437 transcripts were up-regulated and down-regulated respectively in VQL-2; whereas, 437 and 408 transcripts were up-regulated and down-regulated respectively in CM-145. Remarkably, 135 transcription factors and 77 known Zn transporters expressed differentially. By comparing the transcripts differentially expressed between the optimum-Zn and low-Zn libraries of the contrasting genotypes, we identified 21,986 and 26,871 SNPs, respectively. Similarly, 6810 and 8192 InDels were found between optimum- and low-Zn growing conditions, respectively. Further, 21 differentially expressed genes were co-localized with already known QTLs associated with Zn uptake, such as qZn10, CQZnK9-1 and YNZnK6. These findings will be useful to develop high Zn-accumulator maize through marker-assisted breeding in future.

In Silico Identification of QTL-Based Polymorphic Genes as Salt-Responsive Potential Candidates through Mapping with Two Reference Genomes in Rice

Recent advances in next generation sequencing have created opportunities to directly identify genetic loci and candidate genes for abiotic stress responses in plants. With the objective of identifying candidate genes within the previously identified QTL-hotspots, the whole genomes of two divergent cultivars for salt responses, namely At 354 and Bg 352, were re-sequenced using Illumina Hiseq 2500 100PE platform and mapped to Nipponbare and R498 genomes. The sequencing results revealed approximately 2.4 million SNPs and 0.2 million InDels with reference to Nipponbare while 1.3 million and 0.07 million with reference to R498 in two parents. In total, 32,914 genes were reported across all rice chromosomes of this study. Gene mining within QTL hotspots revealed 1236 genes, out of which 106 genes were related to abiotic stress. In addition, 27 abiotic stress-related genes were identified in non-QTL regions. Altogether, 32 genes were identified as potential genes containing polymorphic non-synonymous SNPs or InDels between two parents. Out of 10 genes detected with InDels, tolerant haplotypes of Os01g0581400, Os10g0107000, Os11g0655900, Os12g0622500, and Os12g0624200 were found in the known salinity tolerant donor varieties. Our findings on different haplotypes would be useful in developing resilient rice varieties for abiotic stress by haplotype-based breeding studies.

Genome-wide discovery of DNA variants in cucumber somaclonal lines

The emergence of somaclonal variability in in vitro cultures is undesirable during micropropagation, but this phenomenon may be a source of genetic variability sought by breeders. The main factors that affect the appearance of variability are known, but the exact mechanism has not yet been determined. In this paper, we used next-generation sequencing and comparative genomics to study changes in the genomes of cucumber lines resulting from in vitro regeneration and somaclonal mutation in comparison to a reference, the highly inbred B10 line. The total number of obtained polymorphisms differed between the three somaclonal lines S1, S2 and S3, with 8369, 7591 and 44510, respectively. Polymorphisms occurred most frequently in non-coding regions and were mainly SNPs. High-impact changes accounted for 1%-3% of all polymorphisms and most often caused an open reading frame shift. Functional analysis of genes affected by high impact variants showed that they were related to transport, biosynthetic processes, nucleotide-containing compounds and cellular protein modification processes. The obtained results indicated significant factors affecting somaclonal variability and the appearance of changes in the genome, and demonstrated a lack of dependence between phenotype and the number of genomic polymorphisms.

GWAS Discovery Of Candidate Genes for Yield-Related Traits in Peanut and Support from Earlier QTL Mapping Studies

Peanut (Arachis hypogaea L.) is one of the most important oil crops worldwide, and its yet increasing market demand may be met by genetic improvement of yield related traits, which may be facilitated by a good understanding of the underlying genetic base of these traits. Here, we have carried out a genome-wide association study (GWAS) with the aim to identify genomic regions and the candidate genes within these regions that may be involved in determining the phenotypic variation at seven yield-related traits in peanut. For the GWAS analyses, 195 peanut accessions were phenotyped and/or genotyped; the latter was done using a genotyping-by-sequencing approach, which produced a total of 13,435 high-quality single nucleotide polymorphisms (SNPs). Analyses of these SNPs show that the analyzed peanut accessions can be approximately grouped into two big groups that, to some extent, agree with the botanical classification of peanut at the subspecies level. By taking this genetic structure as well as the relationships between the analyzed accessions into consideration, our GWAS analyses have identified 93 non-overlapping peak SNPs that are significantly associated with four of the studied traits. Gene annotation of the genome regions surrounding these peak SNPs have found a total of 311 unique candidate genes. Among the 93 yield-related-trait-associated SNP peaks, 12 are found to be co-localized with the quantitative trait loci (QTLs) that were identified by earlier related QTL mapping studies, and these 12 SNP peaks are only related to three traits and are almost all located on chromosomes Arahy.05 and Arahy.16. Gene annotation of these 12 co-localized SNP peaks have found 36 candidates genes, and a close examination of these candidate genes found one very interesting gene (arahy.RI9HIF), the rice homolog of which produces a protein that has been shown to improve rice yield when over-expressed. Further tests of the arahy.RI9HIF gene, as well as other candidate genes especially those within the more confident co-localized genomic regions, may hold the potential for significantly improving peanut yield.

Whole Genome Sequencing and Comparative Genomic Analysis Reveal Allelic Variations Unique to a Purple Colored Rice Landrace (Oryza sativa ssp. indica cv. Purpleputtu)

Purpleputtu (Oryza sativa ssp. indica cv. Purpleputtu) is a unique rice landrace from southern India that exhibits predominantly purple color. This study reports the underlying genetic complexity of the trait, associated domestication and de-domestication processes during its coevolution with present day cultivars. Along-with genome level allelic variations in the entire gene repertoire associated with the purple, red coloration of grain and other plant parts. Comparative genomic analysis using 'a panel of 108 rice lines' revealed a total of 3,200,951 variants including 67,774 unique variations in Purpleputtu (PP) genome. Multiple sequence alignment uncovered a 14 bp deletion in Rc (Red colored, a transcription factor of bHLH class) locus of PP, a key regulatory gene of anthocyanin biosynthetic pathway. Interestingly, this deletion in Rc gene is a characteristic feature of the present-day white pericarped rice cultivars. Phylogenetic analysis of Rc locus revealed a distinct clade showing proximity to the progenitor species Oryza rufipogon and O. nivara. In addition, PP genome exhibits a well conserved 4.5 Mbp region on chromosome 5 that harbors several loci associated with domestication of rice. Further, PP showed 1,387 unique when SNPs compared to 3,023 lines of rice (SNP-Seek database). The results indicate that PP genome is rich in allelic diversity and can serve as an excellent resource for rice breeding for a variety of agronomically important traits such as disease resistance, enhanced nutritional values, stress tolerance, and protection from harmful UV-B rays.

Construction of a highly saturated linkage map in Japanese plum (Prunus salicina L.) using GBS for SNP marker calling

This study reports the construction of high density linkage maps of Japanese plum (Prunus salicina Lindl.) using single nucleotide polymorphism markers (SNPs), obtained with a GBS strategy. The mapping population (An x Au) was obtained by crossing cv. "Angeleno" (An) as maternal line and cv. "Aurora" (Au) as the pollen donor. A total of 49,826 SNPs were identified using the peach genome V2.1 as a reference. Then a stringent filtering was carried out, which revealed 1,441 high quality SNPs in 137 An x Au offspring, which were mapped in eight linkage groups. Finally, the consensus map was built using 732 SNPs which spanned 617 cM with an average of 0.96 cM between adjacent markers. The majority of the SNPs were distributed in the intragenic region in all the linkage groups. Considering all linkage groups together, 85.6% of the SNPs were located in intragenic regions and only 14.4% were located in intergenic regions. The genetic linkage analysis was able to co-localize two to three SNPs over 37 putative orthologous genes in eight linkage groups in the Japanese plum map. These results indicate a high level of synteny and collinearity between Japanese plum and peach genomes.

Finding invisible quantitative trait loci with missing data

Evolutionary processes during plant polyploidization and speciation have led to extensive presence-absence variation (PAV) in crop genomes, and there is increasing evidence that PAV associates with important traits. Today, high-resolution genetic analysis in major crops frequently implements simple, cost-effective, high-throughput genotyping from single nucleotide polymorphism (SNP) hybridization arrays; however, these are normally not designed to distinguish PAV from failed SNP calls caused by hybridization artefacts. Here, we describe a strategy to recover valuable information from single nucleotide absence polymorphisms (SNaPs) by population-based quality filtering of SNP hybridization data to distinguish patterns associated with genuine deletions from those caused by technical failures. We reveal that including SNaPs in genetic analyses elucidate segregation of small to large-scale structural variants in nested association mapping populations of oilseed rape (Brassica napus), a recent polyploid crop with widespread structural variation. Including SNaP markers in genomewide association studies identified numerous quantitative trait loci, invisible using SNP markers alone, for resistance to two major fungal diseases of oilseed rape, Sclerotinia stem rot and blackleg disease. Our results indicate that PAV has a strong influence on quantitative disease resistance in B. napus and that SNaP analysis using cost-effective SNP array data can provide extensive added value from 'missing data'. This strategy might also be applicable for improving the precision of genetic mapping in many important crop species.

Genome-wide discovery of DNA polymorphisms by whole genome sequencing differentiates weedy and cultivated rice

Analyzing the genome level DNA polymorphisms between weedy and cultivated rice is crucial to elucidate the molecular basis of weedy and agronomic traits, which in turn can enhance our ability to control weedy rice and its utilization for rice improvement. Here, we presented the genome-wide genetic variations between a weedy rice accession PSRR-1 and two cultivated rice accessions, Bengal and Nona Bokra, belonging to japonica and indica subspecies, respectively. The total number of SNPs and InDels in PSRR/Bengal was similar to that of Nona Bokra/Bengal, but was three times greater than that of PSRR/Nona Bokra. There were 11546 large-effect SNPs/InDels affecting 5673 genes, which most likely differentiated weedy rice from cultivated rice. These large effect DNA polymorphisms were mostly resulted in stop codon gain and least by start codon loss. Analysis of the molecular functions and biological processes of weedy rice specific SNPs/InDels indicated that most of these genes were involved in protein modification/phosphorylation, protein kinase activity, and protein/nucleotide binding. By integrating previous QTL mapping results with the DNA polymorphisms data, the candidate genes for seed dormancy and seed shattering were narrowed down. The genomic resource generated in this study will facilitate discovery of functional variants for weedy and agronomic traits.

Genomics-Assisted Identification and Characterization of the Genetic Variants Underlying Differential Nitrogen Use Efficiencies in Allotetraploid Rapeseed Genotypes

Nitrogen (N) is a non-mineral macronutrient essential for plant growth and development. Oilseed rape (A_nA_nC_nC_n, 2n = 4x = 38) has a high requirement for N nutrients whereas showing the lowest N use efficiency (NUE) among crops. The mechanisms underlying NUE regulation in Brassica napus remain unclear because of genome complexity. In this study, we performed high-depth and -coverage whole-genome re-sequencing (WGS) of an N-efficient (higher NUE) genotype “XY15” and an N-inefficient (lower NUE) genotype “814” of rapeseed. More than 687 million 150-bp paired-end reads were generated, which provided about 93% coverage and 50× depth of the rapeseed genome. Applying stringent parameters, we identified a total of 1,449,157 single-nucleotide polymorphisms (SNPs), 335,228 InDels, 175,602 structure variations (SVs) and 86,280 copy number variations (CNVs) between the N-efficient and -inefficient genotypes. The largest proportion of various DNA polymorphisms occurred in the inter-genic regions. Unlike CNVs, the SNP/InDel and SV polymorphisms showed variation bias of the A_n and C_n subgenomes, respectively. Gene ontology analysis showed the genetic variants were mapped onto the genes involving N compound transport and ATPase complex metabolism, but not including N assimilation-related genes. On basis of identification of N-starvation responsive genes through high-throughput expression profiling, we also mapped these variants onto some key NUE-regulating genes, and validated their significantly differential expression between the N-efficient and -inefficient genotypes through qRT-PCR assays. Our data provide genome-wide high resolution DNA variants underlying NUE divergence in allotetraploid rapeseed genotypes, which would expedite the effective identification and functional validation of key NUE-regulating genes through genomics-assisted improvement of crop nutrient efficiency.

Assessment of Genetic Diversity and Population Structure in Iranian Cannabis Germplasm

Cannabis sativa has a complex history reflected in both selection on naturally occurring compounds and historical trade routes among humans. Iran is a rich resource of natural populationswhich hold the promise to characterize historical patterns of population structure and genetic diversity within Cannabis. Recent advances in high-throughput DNA sequencing technologies have dramatically increased our ability to produce information to the point that it is now feasible to inexpensively obtain population level genotype information at a large scale. In the present investigation, we have explored the use of Genotyping-By-Sequencing (GBS) in Iranian cannabis. We genotyped 98 cannabis samples 36 from Iranian locations and 26 accessions from two germplasm collections. In total, 24,710 high-quality Single Nucleotide Polymorphisms (SNP) were identified. Clustering analysis by Principal Component Analysis (PCA) identified two genetic clusters among Iranian populations and fineSTRUCTURE analysis identified 19 populations with some geographic partitioning. We defined Iranian cannabis in two main groups using the results of the PCA and discovered some strong signal to define some locations as population according to fineSTRUCTURE analyses. However, single nucleotide variant analysis uncovered a relatively moderate level of variation among Iranian cannabis.

Uncovering SNP and indel variations of tetraploid cottons by SLAF-seq

Background

Cotton (Gossypium spp.), as the world’s most utilized textile fibre source, is an important, economically valuable crop worldwide. Understanding the genomic variation of tetraploid cotton species is important for exploitation of the excellent characteristics of wild cotton and for improving the diversity of cotton in breeding. However, the discovery of DNA polymorphisms in tetraploid cotton genomes has lagged behind other important crops.

Results

A total of 111,795,823 reads, 467,735 specific length amplified fragment (SLAF) tags and 139,176 high-quality DNA polymorphisms were identified using specific length amplified fragment sequencing (SLAF-seq), including 132,880 SNPs and 6,296 InDels between the reference genome (TM-1) and the five tetraploid cotton species. Intriguingly, gene ontology (GO) enrichment analysis revealed that a number of significant terms were related to reproduction in G. barbadense acc. 3–79. Based on the new data sets, we reconstructed phylogenetic trees that showed a high concordance to the phylogeny of diploid and polyploid cottons. A large amount of interspecific genetic variations were identified, and some of them were validated by the single-strand conformation polymorphism (SSCP) method, which will be applied in introgression genetics and breeding with G. hirsutum cv. Emian22 as the receptor and the other species as donors.

Conclusions

Using SLAF-seq, a large number of DNA polymorphisms were identified. The comprehensive analysis of DNA polymorphisms provided invaluable insights into the different tetraploid cotton species. More importantly, the identification of numerous interspecific genetic variations provides the basis and is very practical for future introgression breeding. The results presented herein provide a valuable genomic resource for new insights into the genetics and breeding of cotton.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3643-4) contains supplementary material, which is available to authorized users.

SNP Discovery Using a Pangenome: Has the Single Reference Approach Become Obsolete?

Increasing evidence suggests that a single individual is insufficient to capture the genetic diversity within a species due to gene presence absence variation. In order to understand the extent to which genomic variation occurs in a species, the construction of its pangenome is necessary. The pangenome represents the complete set of genes of a species; it is composed of core genes, which are present in all individuals, and variable genes, which are present only in some individuals. Aside from variations at the gene level, single nucleotide polymorphisms (SNPs) are also an important form of genetic variation. The advent of next-generation sequencing (NGS) coupled with the heritability of SNPs make them ideal markers for genetic analysis of human, animal, and microbial data. SNPs have also been extensively used in crop genetics for association mapping, quantitative trait loci (QTL) analysis, analysis of genetic diversity, and phylogenetic analysis. This review focuses on the use of pangenomes for SNP discovery. It highlights the advantages of using a pangenome rather than a single reference for this purpose. This review also demonstrates how extra information not captured in a single reference alone can be used to provide additional support for linking genotypic data to phenotypic data.

Insight into MAS: A Molecular Tool for Development of Stress Resistant and Quality of Rice through Gene Stacking

Rice yield is subjected to severe losses due to adverse effect of a number of stress factors. The most effective method of controlling reduced crop production is utilization of host resistance. Recent technological advances have led to the improvement of DNA based molecular markers closely linked to genes or QTLs in rice chromosome that bestow tolerance to various types of abiotic stresses and resistance to biotic stress factors. Transfer of several genes with potential characteristics into a single genotype is possible through the process of marker assisted selection (MAS), which can quicken the advancement of tolerant/resistant cultivars in the lowest number of generations with the utmost precision through the process of gene pyramiding. Overall, this review presented various types of molecular tools including MAS that can be reasonable and environmental friendly approach for the improvement of abiotic and biotic stress resistant rice with enhanced quality.

Genome-wide SNP discovery and population structure analysis in pepper (Capsicum annuum) using genotyping by sequencing

Background

Knowledge on population structure and genetic diversity in vegetable crops is essential for association mapping studies and genomic selection. Genotyping by sequencing (GBS) represents an innovative method for large scale SNP detection and genotyping of genetic resources. Herein we used the GBS approach for the genome-wide identification of SNPs in a collection of Capsicum spp. accessions and for the assessment of the level of genetic diversity in a subset of 222 cultivated pepper (Capsicum annum) genotypes.

Results

GBS analysis generated a total of 7,568,894 master tags, of which 43.4% uniquely aligned to the reference genome CM334. A total of 108,591 SNP markers were identified, of which 105,184 were in C. annuum accessions. In order to explore the genetic diversity of C. annuum and to select a minimal core set representing most of the total genetic variation with minimum redundancy, a subset of 222 C. annuum accessions were analysed using 32,950 high quality SNPs. Based on Bayesian and Hierarchical clustering it was possible to divide the collection into three clusters. Cluster I had the majority of varieties and landraces mainly from Southern and Northern Italy, and from Eastern Europe, whereas clusters II and III comprised accessions of different geographical origins. Considering the genome-wide genetic variation among the accessions included in cluster I, a second round of Bayesian (K = 3) and Hierarchical (K = 2) clustering was performed. These analysis showed that genotypes were grouped not only based on geographical origin, but also on fruit-related features.

Conclusions

GBS data has proven useful to assess the genetic diversity in a collection of C. annuum accessions. The high number of SNP markers, uniformly distributed on the 12 chromosomes, allowed the accessions to be distinguished according to geographical origin and fruit-related features. SNP markers and information on population structure developed in this study will undoubtedly support genome-wide association mapping studies and marker-assisted selection programs.

Electronic supplementary material

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

Towards a molecular taxonomic key of the Aurantioideae subfamily using chloroplastic SNP diagnostic markers of the main clades genotyped by competitive allele-specific PCR

Background

Chloroplast DNA is a primary source of molecular variations for phylogenetic analysis of photosynthetic eukaryotes. However, the sequencing and analysis of multiple chloroplastic regions is difficult to apply to large collections or large samples of natural populations. The objective of our work was to demonstrate that a molecular taxonomic key based on easy, scalable and low-cost genotyping method should be developed from a set of Single Nucleotide Polymorphisms (SNPs) diagnostic of well-established clades. It was applied to the Aurantioideae subfamily, the largest group of the Rutaceae family that includes the cultivated citrus species.

Results

The publicly available nucleotide sequences of eight plastid genomic regions were compared for 79 accessions of the Aurantioideae subfamily to search for SNPs revealing taxonomic differentiation at the inter-tribe, inter-subtribe, inter-genus and interspecific levels. Diagnostic SNPs (DSNPs) were found for 46 of the 54 clade levels analysed. Forty DSNPs were selected to develop KASPar markers and their taxonomic value was tested by genotyping 108 accessions of the Aurantioideae subfamily. Twenty-seven markers diagnostic of 24 clades were validated and they displayed a very high rate of transferability in the Aurantioideae subfamily (only 1.2 % of missing data on average). The UPGMA from the validated markers produced a cladistic organisation that was highly coherent with the previous phylogenetic analysis based on the sequence data of the eight plasmid regions. In particular, the monophyletic origin of the “true citrus” genera plus Oxanthera was validated. However, some clarification remains necessary regarding the organisation of the other wild species of the Citreae tribe.

Conclusions

We validated the concept that with well-established clades, DSNPs can be selected and efficiently transformed into competitive allele-specific PCR markers (KASPar method) allowing cost-effective highly efficient cladistic analysis in large collections at subfamily level. The robustness of this genotyping method is an additional decisive advantage for network collaborative research. The availability of WGS data for the main “true citrus” species should soon make it possible to develop a set of DSNP markers allowing very fine resolution of this very important horticultural group.

Electronic supplementary material

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

Single nucleotide polymorphism in sugar pathway and disease resistance genes in sugarcane

Single nucleotide polymorphism in sugar pathway and disease resistance genes showing genetic association with sugar content and red rot resistance would be useful in marker-assisted genetic improvement of sugarcane. Validation and genotyping of potential sequence variants in candidate genes are necessary to understand their functional significance and trait association potential. We discovered, characterized, validated and genotyped SNPs and InDels in sugar pathway and disease resistance genes of Saccharum complex and sugarcane varieties using amplicon sequencing and CAPS assays. The SNPs were abundant in the non-coding 3'UTRs than 5'UTRs and coding sequences depicting a strong bias toward C to T transition substitutions than transversions. Sequencing of cloned amplicons validated 61.6 and 45.2 % SNPs detected in silico in 21 sugar pathway and 16 disease resistance genes, respectively. Sixteen SNPs in four sugar pathway genes and 10 SNPs in nine disease resistance genes were validated through cost-effective CAPS assay. Functional and adaptive significance of SNP and protein haplotypes identified in sugar pathway and disease resistance genes was assessed by correlating their allelic variation with missense amino acid substitutions in the functional domains, alteration in protein structure models and possible modulation of catalytic enzyme activity in contrasting high and low sugar and moderately red rot resistant and highly susceptible sugarcane genotypes. A strong genetic association of five SNPs in the sugar pathway and disease resistance genes, and an InDel marker in the promoter sequence of sucrose synthase-2 gene, with sugar content and red rot resistance, was evident. The functionally relevant SNPs and InDels, detected and validated in sugar pathway and disease resistance genes, and genic CAPS markers designed, would be of immense use in marker-assisted genetic improvement of sugarcane for sugar content and disease resistance.

WheatGenome.info: A Resource for Wheat Genomics Resource

An integrated database with a variety of Web-based systems named WheatGenome.info hosting wheat genome and genomic data has been developed to support wheat research and crop improvement. The resource includes multiple Web-based applications, which are implemented as a variety of Web-based systems. These include a GBrowse2-based wheat genome viewer with BLAST search portal, TAGdb for searching wheat second generation genome sequence data, wheat autoSNPdb, links to wheat genetic maps using CMap and CMap3D, and a wheat genome Wiki to allow interaction between diverse wheat genome sequencing activities. This portal provides links to a variety of wheat genome resources hosted at other research organizations. This integrated database aims to accelerate wheat genome research and is freely accessible via the web interface at http://www.wheatgenome.info/ .

Genetic Divergence between Camellia sinensis and Its Wild Relatives Revealed via Genome-Wide SNPs from RAD Sequencing

Tea is one of the most popular beverages across the world and is made exclusively from cultivars of Camellia sinensis. Many wild relatives of the genus Camellia that are closely related to C. sinensis are native to Southwest China. In this study, we first identified the distinct genetic divergence between C. sinensis and its wild relatives and provided a glimpse into the artificial selection of tea plants at a genome-wide level by analyzing 15,444 genomic SNPs that were identified from 18 cultivated and wild tea accessions using a high-throughput genome-wide restriction site-associated DNA sequencing (RAD-Seq) approach. Six distinct clusters were detected by phylogeny inferrence and principal component and genetic structural analyses, and these clusters corresponded to six Camellia species/varieties. Genetic divergence apparently indicated that C. taliensis var. bangwei is a semi-wild or transient landrace occupying a phylogenetic position between those wild and cultivated tea plants. Cultivated accessions exhibited greater heterozygosity than wild accessions, with the exception of C. taliensis var. bangwei. Thirteen genes with non-synonymous SNPs exhibited strong selective signals that were suggestive of putative artificial selective footprints for tea plants during domestication. The genome-wide SNPs provide a fundamental data resource for assessing genetic relationships, characterizing complex traits, comparing heterozygosity and analyzing putatitve artificial selection in tea plants.

Molecular marker development from transcript sequences and germplasm evaluation for cultivated peanut (Arachis hypogaea L.)

Molecular markers are important tools for genotyping in genetic studies and molecular breeding. The SSR and SNP are two commonly used marker systems developed from genomic or transcript sequences. The objectives of this study were to: (1) assemble and annotate the publicly available ESTs in Arachis and the in-house short reads, (2) develop and validate SSR and SNP markers, and (3) investigate the genetic diversity and population structure of the peanut breeding lines and the U.S. peanut mini core collection using developed SSR markers. An NCBI EST dataset with 252,951 sequences and an in-house 454 RNAseq dataset with 288,701 sequences were assembled separately after trimming. Transcript sequence comparison and phylogenetic analysis suggested that peanut is closer to cowpea and scarlet bean than to soybean, common bean and Medicago. From these two datasets, 6455 novel SSRs and 11,902 SNPs were identified. Of the discovered SSRs, 380 representing various SSR types were selected for PCR validation. The amplification rate was 89.2 %. Twenty-two (6.5 %) SSRs were polymorphic between at least one pair of four genotypes. Sanger sequencing of PCR products targeting 110 SNPs revealed 13 true SNPs between tetraploid genotypes and 193 homoeologous SNPs within genotypes. Eight out of the 22 polymorphic SSR markers were selected to evaluate the genetic diversity of Florida peanut breeding lines and the U.S. peanut mini core collection. This marker set demonstrated high discrimination power by displaying an average polymorphism information content value of 0.783, a combined probability of identity of 10(-11), and a combined power of exclusion of 0.99991. The structure analysis revealed four sub-populations among the peanut accessions and lines evaluated. The results of this study enriched the peanut genomic resources, provided over 6000 novel SSR markers and the credentials for true peanut SNP marker development, and demonstrated the power of newly developed SSR markers in genotyping peanut germplasm and breeding materials.

CicArVarDB: SNP and InDel database for advancing genetics research and breeding applications in chickpea

Molecular markers are valuable tools for breeders to help accelerate crop improvement. High throughput sequencing technologies facilitate the discovery of large-scale variations such as single nucleotide polymorphisms (SNPs) and simple sequence repeats (SSRs). Sequencing of chickpea genome along with re-sequencing of several chickpea lines has enabled the discovery of 4.4 million variations including SNPs and InDels. Here we report a repository of 1.9 million variations (SNPs and InDels) anchored on eight pseudomolecules in a custom database, referred as CicArVarDB that can be accessed at http://cicarvardb.icrisat.org/. It includes an easy interface for users to select variations around specific regions associated with quantitative trait loci, with embedded webBLAST search and JBrowse visualisation. We hope that this database will be immensely useful for the chickpea research community for both advancing genetics research as well as breeding applications for crop improvement.

Database URL:http://cicarvardb.icrisat.org.

Genome-wide DNA polymorphisms in low Phosphate tolerant and sensitive rice genotypes

Soil Phosphorus (P) deficiency is one of the major challenges to rice crop world-wide. Modern rice genotypes are highly P-responsive and rely on high input of P fertilizers. However, low P tolerant traditional cultivars and landraces have genetic potential to sustain well under low P. Identification of high resolution DNA polymorphisms (SNPs and InDels) in such contrasting genotypes is largely missing for low P response at gene levels. Here, we report high quality DNA polymorphisms in low P sensitive genotype, PB1 and tolerant traditional genotype, Dular. We performed whole genome resequencing using Illumina NGS platform and identified a total of 5,157,939 sequence variants in PB1 and Dular with reference to Nipponbare genome. We have identified approximately 2.3 million and 2.9 million high quality polymorphisms in PB1 and Dular, respectively, with an average read depth of ≥24X. We further mapped several DNA polymorphisms (non-synonymous and regulatory variants) having potential functional significance to key Phosphate Starvation Responsive (PSR) and root architecture genes in Dular and Kasalath using a compiled list of low P responsive genes. These identified variants can serve as a useful source of genetic variability for improving low P tolerance and root architecture of high yielding modern genotypes.

Whole genome sequencing and analysis of Swarna, a widely cultivated indica rice variety with low glycemic index

Swarna is a popular cultivated indica rice variety with low glycemic index (GI) but its genetic basis is not known. The whole genome of Swarna was sequenced using Illumina’s paired-end technology, and the reads were mapped to the Nipponbare reference genome. Overall, 65,984 non-synonymous SNPs were identified in 20,350 genes, and in silico analysis predicted that 4,847 of them in 2,214 genes may have deleterious effect on protein functions. Polymorphisms were found in all the starch biosynthesis genes, except the gene for branching enzyme IIa. It was found that T/G SNP at position 246, ‘A’ at position 2,386, and ‘C’ at position 3,378 in the granule bound starch synthase I gene, and C/T SNP at position 1,188 in the glucose-6-phosphate translocator gene may contribute to the low GI phenotype in Swarna. All these variants were also found in the genome of another low GI indica rice variety from Columbia, Fedearroz 50. The whole genome analysis of Swarna helped to understand the genetic basis of GI in rice, which is a complex trait involving multiple factors.

ConPADE: genome assembly ploidy estimation from next-generation sequencing data

As a result of improvements in genome assembly algorithms and the ever decreasing costs of high-throughput sequencing technologies, new high quality draft genome sequences are published at a striking pace. With well-established methodologies, larger and more complex genomes are being tackled, including polyploid plant genomes. Given the similarity between multiple copies of a basic genome in polyploid individuals, assembly of such data usually results in collapsed contigs that represent a variable number of homoeologous genomic regions. Unfortunately, such collapse is often not ideal, as keeping contigs separate can lead both to improved assembly and also insights about how haplotypes influence phenotype. Here, we describe a first step in avoiding inappropriate collapse during assembly. In particular, we describe ConPADE (Contig Ploidy and Allele Dosage Estimation), a probabilistic method that estimates the ploidy of any given contig/scaffold based on its allele proportions. In the process, we report findings regarding errors in sequencing. The method can be used for whole genome shotgun (WGS) sequencing data. We also show applicability of the method for variant calling and allele dosage estimation. Results for simulated and real datasets are discussed and provide evidence that ConPADE performs well as long as enough sequencing coverage is available, or the true contig ploidy is low. We show that ConPADE may also be used for related applications, such as the identification of duplicated genes in fragmented assemblies, although refinements are needed.

Diploid organisms, such as human beings, have two “copies” of each chromosome, whereas polyploid organisms have multiple “copies” (we use quotes to stress that the “copies” are not identical). A key difference between diploid and polyploid organisms is that the “copies” tend to be less similar in polyploid organisms. This difference leads to important differences in the process of de novo genome assembly from short fragments of DNA. In particular, when assembling polyploid organisms, contigs corresponding to different copies of the chromosomes can be quite different, and merging them leads to loss of information. Thus, it is important to maintain distinct contigs, even though they correspond to copies of the same chromosomal region. An important step in doing so is to determine how many truly distinct copies of a chromosomal region are found in a single contig. For example, if there are 12 copies of a particular chromosome, the possible number of distinct copies could be anywhere from 1 to 12. We call this task “contig ploidy estimation”, and present a method for accomplishing it. This set of methods is useful for the de novo assembly of complex, polyploid genomes such as sugarcane, switchgrass, and wheat.

High quality SNPs/Indels mining and characterization in ginger from ESTs data base

Ginger (Zingiber officinale Rosc.) is an important herb of the family Zingiberaceae. It is accepted as a universal cure for a multitude of diseases in Indian systems of medicine and its rhizomes are equally popular as a spice ingredient throughout Asia. SNPs, the definitive genetic markers, representing the finest resolution of a DNA sequence, are abundantly found in populations having a lower rate of mutation and are used for genomic analysis. The public ESTs sequences mostly lack quality files, making high quality SNPs detection more difficult since it is exclusively based on sequence comparisons. In the present study, current dbESTs of NCBI was mined and 38115 ginger ESTs sequences were obtained and assembled into contigs using CAP3 program. In this analysis, recent software tool QualitySNP was used to detect 11523 potential SNPs sites, 8810 high quality SNPs and 1008 indels polymorphisms with a frequency of 1.61 SNPs / 10 kbp. Of ESTs libraries generated from three ginger tissues together, rhizomes had a frequency of 0.32 SNPs and 0.03 indels per 10 kbp whereas the leaves had a frequency of 2.51 SNPs and 0.23 indels per 10 kbp and root is showing relative frequency of 0.76/10 kbp SNPs and 0.02/10 kbp indels. The present analysis provides additional information about the tissue wise presence of haplotypes (222), distribution of high quality exonic (2355) and intronic (6455) SNPs and information about singletons (7538) in addition to contigs transitions and transversions ratio (0.57). Among all tissue detected SNPs, transversions number is higher in comparison to the number of transitions. Quality SNPs detected in this work can be used as markers for further ginger genetic experiments.

Bioinformatics: identification of markers from next-generation sequence data

With the advent of sequencing technology, next-generation sequencing (NGS) technology has dramatically revolutionized plant genomics. NGS technology combined with new software tools enables the discovery, validation, and assessment of genetic markers on a large scale. Among different markers systems, simple sequence repeats (SSRs) and Single nucleotide polymorphisms (SNPs) are the markers of choice for genetics and plant breeding. SSR markers have been a choice for large-scale characterization of germplasm collections, construction of genetic maps, and QTL identification. Similarly, SNPs are the most abundant genetic variations with higher frequencies throughout the genome of plant species. This chapter discusses various tools available for genome assembly and widely focuses on SSR and SNP marker discovery.

Molecular marker databases

The detection and analysis of genetic variation plays an important role in plant breeding and this role is increasing with the continued development of genome sequencing technologies. Molecular genetic markers are important tools to characterize genetic variation and assist with genomic breeding. Processing and storing the growing abundance of molecular marker data being produced requires the development of specific bioinformatics tools and advanced databases. Molecular marker databases range from species specific through to organism wide and often host a variety of additional related genetic, genomic, or phenotypic information. In this chapter, we will present some of the features of plant molecular genetic marker databases, highlight the various types of marker resources, and predict the potential future direction of crop marker databases.

Identification and characterization of more than 4 million intervarietal SNPs across the group 7 chromosomes of bread wheat

Despite being a major international crop, our understanding of the wheat genome is relatively poor due to its large size and complexity. To gain a greater understanding of wheat genome diversity, we have identified single nucleotide polymorphisms between 16 Australian bread wheat varieties. Whole-genome shotgun Illumina paired read sequence data were mapped to the draft assemblies of chromosomes 7A, 7B and 7D to identify more than 4 million intervarietal SNPs. SNP density varied between the three genomes, with much greater density observed on the A and B genomes than the D genome. This variation may be a result of substantial gene flow from the tetraploid Triticum turgidum, which possesses A and B genomes, during early co-cultivation of tetraploid and hexaploid wheat. In addition, we examined SNP density variation along the chromosome syntenic builds and identified genes in low-density regions which may have been selected during domestication and breeding. This study highlights the impact of evolution and breeding on the bread wheat genome and provides a substantial resource for trait association and crop improvement. All SNP data are publically available on a generic genome browser GBrowse at www.wheatgenome.info.

Development of EST-based SNP and InDel markers and their utilization in tetraploid cotton genetic mapping

Background

Availability of molecular markers has proven to be an efficient tool in facilitating progress in plant breeding, which is particularly important in the case of less researched crops such as cotton. Considering the obvious advantages of single nucleotide polymorphisms (SNPs) and insertion-deletion polymorphisms (InDels), expressed sequence tags (ESTs) were analyzed in silico to identify SNPs and InDels in this study, aiming to develop more molecular markers in cotton.

Results

A total of 1,349 EST-based SNP and InDel markers were developed by comparing ESTs between Gossypium hirsutum and G. barbadense, mining G. hirsutum unigenes, and analyzing 3′ untranslated region (3′UTR) sequences. The marker polymorphisms were investigated using the two parents of the mapping population based on the single-strand conformation polymorphism (SSCP) analysis. Of all the markers, 137 (10.16%) were polymorphic, and revealed 142 loci. Linkage analysis using a BC₁ population mapped 133 loci on the 26 chromosomes. Statistical analysis of base variations in SNPs showed that base transitions accounted for 55.78% of the total base variations and gene ontology indicated that cotton genes varied greatly in harboring SNPs ranging from 1.00 to 24.00 SNPs per gene. Sanger sequencing of three randomly selected SNP markers revealed discrepancy between the in silico predicted sequences and the actual sequencing results.

Conclusions

In silico analysis is a double-edged blade to develop EST-SNP/InDel markers. On the one hand, the designed markers can be well used in tetraploid cotton genetic mapping. And it plays a certain role in revealing transition preference and SNP frequency of cotton genes. On the other hand, the developmental efficiency of markers and polymorphism of designed primers are comparatively low.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1046) contains supplementary material, which is available to authorized users.

Genomic reduction assisted single nucleotide polymorphism discovery using 454-pyrosequencing

We report the development of a simple genomic reduction protocol based on 454-pyrosequencing technology that discovers large numbers of single nucleotide polymorphisms (SNP) from pooled DNA samples. The method is based on the conservation of restriction endonuclease sites across samples and biotin separation for genomic reduction and the addition of multiplex identifier (MID) barcodes to each of the pooled samples to allow for postsequencing deconvolution of the pooled DNA fragments and SNP discovery.

Exploration of presence/absence variation and corresponding polymorphic markers in soybean genome

This study was designed to reveal the genome-wide distribution of presence/absence variation (PAV) and to establish a database of polymorphic PAV markers in soybean. The 33 soybean whole-genome sequences were compared to each other with that of Williams 82 as a reference genome. A total of 33,127 PAVs were detected and 28,912 PAV markers with their primer sequences were designed as the database NJAUSoyPAV_1.0. The PAVs scattered on whole genome while only 518 (1.8%) overlapped with simple sequence repeats (SSRs) in BARCSOYSSR_1.0 database. In a random sample of 800 PAVs, 713 (89.13%) showed polymorphism among the 12 differential genotypes. Using 126 PAVs and 108 SSRs to test a Chinese soybean germplasm collection composed of 828 Glycine soja Sieb. et Zucc. and Glycine max (L.) Merr. accessions, the per locus allele number and its variation appeared less in PAVs than in SSRs. The distinctness among alleles/bands of PCR (polymerase chain reaction) products showed better in PAVs than in SSRs, potential in accurate marker-assisted allele selection. The association mapping results showed SSR + PAV was more powerful than any single marker systems. The NJAUSoyPAV_1.0 database has enriched the source of PCR markers, and may fit the materials with a range of per locus allele numbers, if jointly used with SSR markers.

Tissue culture-induced heritable genomic variation in rice, and their phenotypic implications

BACKGROUND

Somaclonal variation generally occurs in plants regenerated from tissue culture. However, fundamental issues regarding molecular characteristics, mutation rates and mutation spectra of plant somatic variation as well as their phenotypic relevance have been addressed only recently. Moreover, these studies have reported highly discrepant results in different plant species and even in the same plant genotype.

METHODOLOGY/PRINCIPAL FINDINGS

We investigated heritable genomic variation induced by tissue culture in rice by whole genome re-sequencing of an extensively selfed somaclonal line (TC-reg-2008) and its wild type (WT) donor (cv. Hitomebore). We computed the overall mutation rate, single nucleotide polymorphisms (SNPs), small scale insertions/deletions (Indels) and mobilization of transposable elements (TEs). We assessed chromosomal distribution of the various types of genomic variations, tested correlations between SNPs and Indels, and examined concomitancy between TE activity and its cytosine methylation states. We also performed gene ontology (GO) analysis of genes containing nonsynonymous mutations and large-effect mutations, and assayed effects of the genomic variations on phenotypes under both normal growing condition and several abiotic stresses. We found that heritable somaclonal genomic variation occurred extensively in rice. The genomic variations distributed non-randomly across each of the 12 rice chromosomes, and affected a large number of functional genes. The phenotypic penetrance of the genomic variations was condition-dependent.

CONCLUSIONS/SIGNIFICANCE

Tissue culture is a potent means to generate heritable genetic variations in rice, which bear distinct difference at least in space (chromosomal distribution) from those occurred under natural settings. Our findings have provided new information regarding the mutation rate and spectrum as well as chromosomal distribution pattern of somaclonal variation in rice. Our data also suggest that rice possesses a strong capacity to canalize genetic variations under normal growing conditions to maintain phenotypic robustness, which however can be released by certain abiotic stresses to generate variable phenotypes.

Molecular markers for genetics and plant breeding: the MFLP marker system and its application in narrow-leafed lupin (Lupinus angustifolius)

Since the development of molecular markers to tag genes of agronomic traits of interests, molecular markers have played an increasingly significant role in breeding programs. Molecular markers have been implemented for large-scale marker-assisted selection in the breeding program of many important crops including lupin. So far, more than a dozen molecular markers for disease resistance genes and for other agronomic traits of interest have been developed in lupin. The DNA fingerprinting method, "MFLP" has played a pivotal role in the success of lupin breeding program in Australia. Here, we describe the MFLP technique used in lupin breeding which could be easily transferable to other crop species.