Skim sequencing: an advanced NGS technology for crop improvement

Using genotype imputation to integrate Canola populations for genome-wide association and genomic prediction of blackleg resistance

BACKGROUND

Integrating germplasm populations genotyped by different genotyping platforms via genotype imputation is a way to utilize accumulated genetic resources. In this study, we used 278 canola samples genotyped via whole-genome sequencing (WGS) at 10× coverage to evaluate the imputation accuracy of three imputation approaches. The optimal imputation methods were used to impute and integrate two Canola genotype datasets: a diverse canola collection genotyped by genotyping-by-sequencing via transcriptome (GBS-t) and a double haploid (DH) line collection genotyped with low-coverage WGS (skim-WGS). The genomic predictive ability (GP) and detection power of marker‒trait association (GWAS) of the combined population for blackleg resistance were evaluated.

RESULTS

The empirical imputation accuracy (r2) measured as the squared correlation between observed and imputed genotypes was moderate for Minimac3 when imputing from the GBS-t density to the WGS. The accuracy dramatically improved from 0.64 to 0.82 by removing SNPs with poor Minimac3-reported Rsq (Rsq < 0.2) quality statistics. The r2 for GLIMPSE was higher than that for Beagle when imputing from different low-coverage to full-coverage WGS. We imputed and integrated the diverse canola collection and the DH lines, and the combined population showed similar or slightly greater predictive ability (PA) for blackleg resistance traits than did each of the single populations with ~ 921 K SNPs. Higher marker-trait association (MTA) detection powers were indicated with the combined population; however, similar numbers of MTAs were discovered when each single population was combined in a meta-GWAS.

CONCLUSION

It is feasible to impute and integrate germplasms from different sequencing platforms for downstream analyses. However, genetic heterogeneity across populations could add complexity to the analysis. Increasing the sample size by combining datasets showed slightly greater predictive ability and greater detection power in GWASs in the present study.

Low-pass nanopore sequencing for measurement of global methylation levels in plants

Nanopore sequencing enables detection of DNA methylation at the same time as identification of canonical sequence. A recent study validated low-pass nanopore sequencing to accurately estimate global methylation levels in vertebrates with sequencing coverage as low as 0.01x. We investigated the applicability of this approach to plants by testing three plant species and analysed the effect of technical and biological parameters on estimate precision and accuracy. Our results indicate that higher coverage (0.1x) is required to achieve accuracy in assessing plant global methylation comparable to that in vertebrates. Shorter read length and a closer sequence match between sample and reference genome improved measurement accuracy. Application of this method in Vitis vinifera showed consistent global methylation levels across different leaf sizes, and different sample preservation and DNA extraction methods, whereas different varieties and tissue types did exhibit methylation differences. Similarly, distinct methylation patterns were observed in different genomic features. Our findings suggest the suitability of this method as a low-cost screening tool for validation of experimental parameters, developmental time courses, and to assess methylation status for different modification types and sequence contexts at the level of whole genome or for abundant genomic features such as transposable elements.

Genome-wide association analyses using multilocus models on bananas (Musa spp.) reveal candidate genes related to morphology, fruit quality, and yield

Bananas (Musa spp.) are an essential fruit worldwide and rank as the fourth most significant food crop for addressing malnutrition due to their rich nutrients and starch content. The potential of their genetic diversity remains untapped due to limited molecular breeding tools. Our study examined a phenotypically diverse group of 124 accessions from the Colombian Musaceae Collection conserved in AGROSAVIA. We assessed 12 traits categorized into morphology, fruit quality, and yield, alongside sequence data. Our sequencing efforts provided valuable insights, with an average depth of about 7× per accession, resulting in 187,133 single-nucleotide polymorphisms (SNPs) against Musa acuminata (A genome) and 220,451 against Musa balbisiana (B genome). Population structure analysis grouped samples into four and five clusters based on the reference genome. By using different association models, we identified marker-trait associations (MTAs). The mixed linear model revealed four MTAs, while the Bayesian-information and linkage-disequilibrium iteratively nested keyway and fixed and random model for circulating probability unification models identified 82 and 70 MTAs, respectively. We identified 38 and 40 candidate genes in linkage proximity to significant MTAs for the A genome and B genome, respectively. Our findings provide insights into the genetic underpinnings of morphology, fruit quality, and yield. Once validated, the SNP markers and candidate genes can potentially drive advancements in genomic-guided breeding strategies to enhance banana crop improvement.

Polly: An R package for genotyping microsatellites and detecting highly polymorphic DNA markers from short-read data

Highly polymorphic markers, such as microsatellites, are invaluable for the study of natural populations. However, contemporary methods for genotyping highly polymorphic variants have serious drawbacks that impede their efficiency. We created Polly, an R package with C++ source code that uses Illumina short-read data to genotype microsatellites, detect highly polymorphic variants and identify clusters of highly polymorphic SNPs, indels and microsatellites. We tested Polly on short-read data from Xiphophorus birchmanni (Teleostei: Poeciliidae) and Arabidopsis thaliana, finding it to be efficient and accurate both for microsatellite genotyping and polymorphic marker detection. This program can be applied to any diploid population for which there exists short-read data and at least one scaffolded reference genome.

Genotype and phenotype data standardization, utilization and integration in the big data era for agricultural sciences

Large-scale genotype and phenotype data have been increasingly generated to identify genetic markers, understand gene function and evolution and facilitate genomic selection. These datasets hold immense value for both current and future studies, as they are vital for crop breeding, yield improvement and overall agricultural sustainability. However, integrating these datasets from heterogeneous sources presents significant challenges and hinders their effective utilization. We established the Genotype-Phenotype Working Group in November 2021 as a part of the AgBioData Consortium (https://www.agbiodata.org) to review current data types and resources that support archiving, analysis and visualization of genotype and phenotype data to understand the needs and challenges of the plant genomic research community. For 2021-22, we identified different types of datasets and examined metadata annotations related to experimental design/methods/sample collection, etc. Furthermore, we thoroughly reviewed publicly funded repositories for raw and processed data as well as secondary databases and knowledgebases that enable the integration of heterogeneous data in the context of the genome browser, pathway networks and tissue-specific gene expression. Based on our survey, we recommend a need for (i) additional infrastructural support for archiving many new data types, (ii) development of community standards for data annotation and formatting, (iii) resources for biocuration and (iv) analysis and visualization tools to connect genotype data with phenotype data to enhance knowledge synthesis and to foster translational research. Although this paper only covers the data and resources relevant to the plant research community, we expect that similar issues and needs are shared by researchers working on animals. Database URL: https://www.agbiodata.org.

A multitiered haplotype strategy to enhance phased assembly and fine mapping of a disease resistance locus

Fine mapping of quantitative trait loci (QTL) to dissect the genetic basis of traits of interest is essential to modern breeding practice. Here, we employed a multitiered haplotypic marker system to increase fine mapping accuracy by constructing a chromosome-level, haplotype-resolved parental genome, accurate detection of recombination sites, and allele-specific characterization of the transcriptome. In the first tier of this system, we applied the preexisting panel of 2,000 rhAmpSeq core genome markers that is transferable across the entire Vitis genus and provides a genomic resolution of 200 kb to 1 Mb. The second tier consisted of high-density haplotypic markers generated from Illumina skim sequencing data for samples enriched for relevant recombinations, increasing the potential resolution to hundreds of base pairs. We used this approach to dissect a novel Resistance to Plasmopara viticola-33 (RPV33) locus conferring resistance to grapevine downy mildew, narrowing the candidate region to only 0.46 Mb. In the third tier, we used allele-specific RNA-seq analysis to identify a cluster of 3 putative disease resistance RPP13-like protein 2 genes located tandemly in a nonsyntenic insertion as candidates for the disease resistance trait. In addition, combining the rhAmpSeq core genome haplotype markers and skim sequencing-derived high-density haplotype markers enabled chromosomal-level scaffolding and phasing of the grape Vitis × doaniana 'PI 588149' assembly, initially built solely from Pacific Biosciences (PacBio) high-fidelity (HiFi) reads, leading to the correction of 16 large-scale phasing errors. Our mapping strategy integrates high-density, phased genetic information with individual reference genomes to pinpoint the genetic basis of QTLs and will likely be widely adopted in highly heterozygous species.

DNA storage-from natural biology to synthetic biology

Natural DNA storage allows cellular differentiation, evolution, the growth of our children and controls all our ecosystems. Here, we discuss the fundamental aspects of DNA storage and recent advances in this field, with special emphasis on natural processes and solutions that can be exploited. We point out new ways of efficient DNA and nucleotide storage that are inspired by nature. Within a few years DNA-based information storage may become an attractive and natural complementation to current electronic data storage systems. We discuss rapid and directed access (e.g. DNA elements such as promotors, enhancers), regulatory signals and modulation (e.g. lncRNA) as well as integrated high-density storage and processing modules (e.g. chromosomal territories). There is pragmatic DNA storage for use in biotechnology and human genetics. We examine DNA storage as an approach for synthetic biology (e.g. light-controlled nucleotide processing enzymes). The natural polymers of DNA and RNA offer much for direct storage operations (read-in, read-out, access control). The inbuilt parallelism (many molecules at many places working at the same time) is important for fast processing of information. Using biology concepts from chromosomal storage, nucleic acid processing as well as polymer material sciences such as electronical effects in enzymes, graphene, nanocellulose up to DNA macramé , DNA wires and DNA-based aptamer field effect transistors will open up new applications gradually replacing classical information storage methods in ever more areas over time (decades).

High-resolution quantitative trait locus mapping for rice grain quality traits using genotyping by sequencing

Rice is a major food crop that sustains approximately half of the world population. Recent worldwide improvements in the standard of living have increased the demand for high-quality rice. Accurate identification of quantitative trait loci (QTLs) for rice grain quality traits will facilitate rice quality breeding and improvement. In the present study, we performed high-resolution QTL mapping for rice grain quality traits using a genotyping-by-sequencing approach. An F₂ population derived from a cross between an elite japonica variety, Koshihikari, and an indica variety, Nona Bokra, was used to construct a high-density genetic map. A total of 3,830 single nucleotide polymorphism markers were mapped to 12 linkage groups spanning a total length of 2,456.4 cM, with an average genetic distance of 0.82 cM. Seven grain quality traits-the percentage of whole grain, percentage of head rice, percentage of area of head rice, transparency, percentage of chalky rice, percentage of chalkiness area, and degree of chalkiness-of the F₂ population were investigated. In total, 15 QTLs with logarithm of the odds (LOD) scores >4 were identified, which mapped to chromosomes 6, 7, and 9. These loci include four QTLs for transparency, four for percentage of chalky rice, four for percentage of chalkiness area, and three for degree of chalkiness, accounting for 0.01%-61.64% of the total phenotypic variation. Of these QTLs, only one overlapped with previously reported QTLs, and the others were novel. By comparing the major QTL regions in the rice genome, several key candidate genes reported to play crucial roles in grain quality traits were identified. These findings will expedite the fine mapping of these QTLs and QTL pyramiding, which will facilitate the genetic improvement of rice grain quality.

High-Density Mapping of Quantitative Trait Loci Controlling Agronomically Important Traits in Quinoa (Chenopodium quinoa Willd.)

Quinoa is a pseudocereal originating from the Andean regions. Despite quinoa's long cultivation history, genetic analysis of this crop is still in its infancy. We aimed to localize quantitative trait loci (QTL) contributing to the phenotypic variation of agronomically important traits. We crossed the Chilean accession PI-614889 and the Peruvian accession CHEN-109, which depicted significant differences in days to flowering, days to maturity, plant height, panicle length, and thousand kernel weight (TKW), saponin content, and mildew susceptibility. We observed sizeable phenotypic variation across F2 plants and F3 families grown in the greenhouse and the field, respectively. We used Skim-seq to genotype the F2 population and constructed a high-density genetic map with 133,923 single nucleotide polymorphism (SNPs). Fifteen QTL were found for ten traits. Two significant QTL, common in F2 and F3 generations, depicted pleiotropy for days to flowering, plant height, and TKW. The pleiotropic QTL harbored several putative candidate genes involved in photoperiod response and flowering time regulation. This study presents the first high-density genetic map of quinoa that incorporates QTL for several important agronomical traits. The pleiotropic loci can facilitate marker-assisted selection in quinoa breeding programs.

Meta-QTL analysis and candidate genes identification for various abiotic stresses in maize (Zea mays L.) and their implications in breeding programs

Global climate change leads to the concurrence of a number of abiotic stresses including moisture stress (drought, waterlogging), temperature stress (heat, cold), and salinity stress, which are the major factors affecting maize production. To develop abiotic stress tolerance in maize, many quantitative trait loci (QTL) have been identified, but very few of them have been utilized successfully in breeding programs. In this context, the meta-QTL analysis of the reported QTL will enable the identification of stable/real QTL which will pave a reliable way to introgress these QTL into elite cultivars through marker-assisted selection. In this study, a total of 542 QTL were summarized from 33 published studies for tolerance to different abiotic stresses in maize to conduct meta-QTL analysis using BiomercatorV4.2.3. Among those, only 244 major QTL with more than 10% phenotypic variance were preferably utilised to carry out meta-QTL analysis. In total, 32 meta-QTL possessing 1907 candidate genes were detected for different abiotic stresses over diverse genetic and environmental backgrounds. The MQTL2.1, 5.1, 5.2, 5.6, 7.1, 9.1, and 9.2 control different stress-related traits for combined abiotic stress tolerance. The candidate genes for important transcription factor families such as ERF, MYB, bZIP, bHLH, NAC, LRR, ZF, MAPK, HSP, peroxidase, and WRKY have been detected for different stress tolerances. The identified meta-QTL are valuable for future climate-resilient maize breeding programs and functional validation of candidate genes studies, which will help to deepen our understanding of the complexity of these abiotic stresses.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01294-9.

Novel Insights Into Refugia at the Southern Margin of the Distribution Range of the Endangered Species Ulmus laevis

Riparian ecosystems, in long-time developed regions, are among the most heavily impacted by human activities; therefore, the distribution of tree riparian species, such as Ulmus laevis, is highly affected. This phenomenon is particularly relevant at the margins of the natural habitat of the species, where populations are small and rare. In these cases, it is difficult to distinguish between relics or introductions, but it is relevant for the restoration of natural habitats and conservation strategies. The aim of this study was to study the phylogeography of the southern distribution of the species. We sequenced the entire chloroplast (cp) genomes of 54 individuals from five sampled populations across different European regions to highlight polymorphisms and analyze their distribution. Thirty-two haplotypes were identified. All the sampled populations showed private haplotypes that can be considered an indicator of long-term residency, given the low mutation rate of organellar DNA. The network of all haplotypes showed a star-like topology, and Serbian haplotypes were present in all branches. The Balkan population showed the highest level of nucleotide and genetic diversity. Low genetic differentiation between populations was observed but we found a significant differentiation among Serbia vs. other provenances. Our estimates of divergent time of U. laevis samples highlight the early split of above all Serbian individuals from other populations, emphasizing the reservoir role of white elm genetic diversity of Serbian population.