Targeted enrichment by solution-based hybrid capture to identify genetic sequence variants in barley

Rare but diverse off-target and somatic mutations found in field and greenhouse grown trees expressing CRISPR/Cas9

Introduction: CRISPR gene editing, while highly efficient in creating desired mutations, also has the potential to cause off-target mutations. This risk is especially high in clonally propagated plants, where editing reagents may remain in the genome for long periods of time or in perpetuity. We studied a diverse population of Populus and Eucalyptus trees that had CRISPR/Cas9-containing transgenes that targeted one or two types of floral development genes, homologs of LEAFY and AGAMOUS. Methods: Using a targeted sequence approach, we studied approximately 20,000 genomic sites with degenerate sequence homology of up to five base pairs relative to guide RNA (gRNA) target sites. We analyzed those sites in 96 individual tree samples that represented 37 independent insertion events containing one or multiples of six unique gRNAs. Results: We found low rates of off-target mutations, with rates of 1.2 × 10-9 in poplar and 3.1 × 10-10 in eucalypts, respectively, comparable to that expected due to sexual reproduction. The rates of mutation were highly idiosyncratic among sites and not predicted by sequence similarity to the target sites; a subset of two gRNAs showed off-target editing of four unique genomic sites with up to five mismatches relative to the true target sites, reaching fixation in some gene insertion events and clonal ramets. The location of off-target mutations relative to the PAM site were essentially identical to that seen with on-target CRISPR mutations. Discussion: The low rates observed support many other studies in plants that suggest that the rates of off-target mutagenesis from CRISPR/Cas9 transgenes are negligible; our study extends this conclusion to trees and other long-lived plants where CRISPR/Cas9 transgenes were present in the genome for approximately four years.

Single-molecule long-read methylation profiling reveals regional DNA methylation regulated by Elongator Complex Subunit 2 in Arabidopsis roots experiencing spaceflight

BACKGROUND

The Advanced Plant Experiment-04 - Epigenetic Expression (APEX-04-EpEx) experiment onboard the International Space Station examined the spaceflight-altered cytosine methylation in two genetic lines of Arabidopsis thaliana, wild-type Col-0 and the mutant elp2-5, which is deficient in an epigenetic regulator Elongator Complex Subunit 2 (ELP2). Whole-genome bisulfite sequencing (WGBS) revealed distinct spaceflight associated methylation differences, presenting the need to explore specific space-altered methylation at single-molecule resolution to associate specific changes over large regions of spaceflight related genes. To date, tools of multiplexed targeted DNA methylation sequencing remain limited for plant genomes.

RESULTS

To provide methylation data at single-molecule resolution, Flap-enabled next-generation capture (FENGC), a novel targeted multiplexed DNA capture and enrichment technique allowing cleavage at any specified sites, was applied to survey spaceflight-altered DNA methylation in genic regions of interest. The FENGC capture panel contained 108 targets ranging from 509 to 704 nt within the promoter or gene body regions of gene targets derived from spaceflight whole-genome data sets. In addition to genes with significant changes in expression and average methylation levels between spaceflight and ground control, targets with space-altered distributions of the proportion of methylated cytosines per molecule were identified. Moreover, trends of co-methylation of different cytosine contexts were exhibited in the same DNA molecules. We further identified significant DNA methylation changes in three previously biological process-unknown genes, and loss-of-function mutants of two of these genes (named as EMO1 and EMO2 for ELP2-regulated Methylation in Orbit 1 and 2) showed enhanced root growth rate.

CONCLUSIONS

FENGC simplifies and reduces the cost of multiplexed, targeted, single-molecule profiling of methylation in plants, providing additional resolution along each DNA molecule that is not seen in population-based short-read data such as WGBS. This case study has revealed spaceflight-altered regional modification of cytosine methylation occurring within single DNA molecules of cell subpopulations, which were not identified by WGBS. The single-molecule survey by FENGC can lead to identification of novel functional genes. The newly identified EMO1 and EMO2 are root growth regulators which may be epigenetically involved in plant adaptation to spaceflight.

Diversity of Gibberellin 2-oxidase genes in the barley genome offers opportunities for genetic improvement

INTRODUCTION

Gibberellin (GA) is a vital phytohormone in regulating plant growth and development. During the "Green Revolution", modification of GA-related genes created semi-dwarfing phenotype in cereal crops but adversely affected grain weight. Gibberellin 2-oxidases (GA2oxs) in barley act as key catabolic enzymes in deactivating GA, but their functions are still less known.

OBJECTIVES

This study investigates the physiological function of two HvGA2ox genes in barley and identifies novel semi-dwarf alleles with minimum impacts on other agronomic traits.

METHODS

Virus-induced gene silencing and CRISPR/Cas9 technology were used to manipulate gene expression of HvGA2ox9 and HvGA2ox8a in barley and RNA-seq was conducted to compare the transcriptome between wild type and mutants. Also, field trials in multiple environments were performed to detect the functional haplotypes.

RESULTS

There were ten GA2oxs that distinctly expressed in shoot, tiller, inflorescence, grain, embryo and root. Knockdown of HvGA2ox9 did not affect plant height, while ga2ox8a mutants generated by CRISPR/Cas9 increased plant height and significantly altered seed width and weight due to the increased bioactive GA4 level. RNA-seq analysis revealed that genes involved in starch and sucrose metabolism were significantly decreased in the inflorescence of ga2ox8a mutants. Furthermore, haplotype analysis revealed one naturally occurring HvGA2ox8a haplotype was associated with decreased plant height, early flowering and wider and heavier seed.

CONCLUSION

Our results demonstrate the potential of manipulating GA2ox genes to fine tune GA signalling and biofunctions in desired plant tissues and open a promising avenue for minimising the trade-off effects of Green Revolution semi-dwarfing genes on grain size and weight. The knowledge will promote the development of next generation barley cultivars with better adaptation to a changing climate.

Next-generation fungal identification using target enrichment and Nanopore sequencing

BACKGROUND

Rapid and accurate pathogen identification is required for disease management. Compared to sequencing entire genomes, targeted sequencing may be used to direct sequencing resources to genes of interest for microbe identification and mitigate the low resolution that single-locus molecular identification provides. This work describes a broad-spectrum fungal identification tool developed to focus high-throughput Nanopore sequencing on genes commonly employed for disease diagnostics and phylogenetic inference.

RESULTS

Orthologs of targeted genes were extracted from 386 reference genomes of fungal species spanning six phyla to identify homologous regions that were used to design the baits used for enrichment. To reduce the cost of producing probes without diminishing the phylogenetic power, DNA sequences were first clustered, and then consensus sequences within each cluster were identified to produce 26,000 probes that targeted 114 genes. To test the efficacy of our probes, we applied the technique to three species representing Ascomycota and Basidiomycota fungi. The efficiency of enrichment, quantified as mean target coverage over the mean genome-wide coverage, ranged from 200 to 300. Furthermore, enrichment of long reads increased the depth of coverage across the targeted genes and into non-coding flanking sequence. The assemblies generated from enriched samples provided well-resolved phylogenetic trees for taxonomic assignment and molecular identification.

CONCLUSIONS

Our work provides data to support the utility of targeted Nanopore sequencing for fungal identification and provides a platform that may be extended for use with other phytopathogens.

New semi-dwarfing alleles with increased coleoptile length by gene editing of gibberellin 3-oxidase 1 using CRISPR-Cas9 in barley (Hordeum vulgare L.)

The green revolution was based on genetic modification of the gibberellin (GA) hormone system with "dwarfing" gene mutations that reduces GA signals, conferring shorter stature, thus enabling plant adaptation to modern farming conditions. Strong GA-related mutants with shorter stature often have reduced coleoptile length, discounting yield gain due to their unsatisfactory seedling emergence under drought conditions. Here we present gibberellin (GA) 3-oxidase1 (GA3ox1) as an alternative semi-dwarfing gene in barley that combines an optimal reduction in plant height without restricting coleoptile and seedling growth. Using large-scale field trials with an extensive collection of barley accessions, we showed that a natural GA3ox1 haplotype moderately reduced plant height by 5-10 cm. We used CRISPR/Cas9 technology, generated several novel GA3ox1 mutants and validated the function of GA3ox1. We showed that altered GA3ox1 activities changed the level of active GA isoforms and consequently increased coleoptile length by an average of 8.2 mm, which could provide essential adaptation to maintain yield under climate change. We revealed that CRISPR/Cas9-induced GA3ox1 mutations increased seed dormancy to an ideal level that could benefit the malting industry. We conclude that selecting HvGA3ox1 alleles offers a new opportunity for developing barley varieties with optimal stature, longer coleoptile and additional agronomic traits.

Open Government Data: Usage trends and metadata quality

Open Government Data (OGD) have the potential to support social and economic progress. However, this potential can be frustrated if these data remain unused. Although the literature suggests that OGD data sets’ metadata quality is one of the main factors affecting their use, to the best of our knowledge, no quantitative study provided evidence of this relationship. Considering about 400,000 data sets of 28 national, municipal and international OGD portals, we have programmatically analysed their usage, their metadata quality and the relationship between the two. Our analysis has highlighted three main findings. First, regardless of their size, the software platform adopted, and their administrative and territorial coverage, most OGD data sets are underutilised. Second, OGD portals pay varying attention to the quality of their data sets’ metadata. Third, we did not find clear evidence that data sets’ usage is positively correlated to better metadata publishing practices. Finally, we have considered other factors, such as data sets’ category, and some demographic characteristics of the OGD portals, and analysed their relationship with data sets’ usage, obtaining partially affirmative answers.

Development and application of the Faba_bean_130K targeted next-generation sequencing SNP genotyping platform based on transcriptome sequencing

KEY MESSAGE

Large-scale faba bean transcriptome data are available, and the first genotyping platform based on liquid-phase probe targeted capture technology was developed for genetic and molecular breeding studies. Faba bean (Vicia faba L., 2n = 12) is an important food legume crop that is widely grown for multiple uses worldwide. However, no reference genome is currently available due to its very large genome size (approximately 13 Gb) and limited single nucleotide polymorphism (SNP) markers as well as highly efficient genotyping tools have been reported for faba bean. In this study, 16.7 billion clean reads were obtained from transcriptome libraries of flowers and leaves of 102 global faba bean accessions. A total of 243,120 unigenes were de novo assembled and functionally annotated. Moreover, a total of 1,579,411 SNPs were identified and further filtered according to a selection pipeline to develop a high-throughput, flexible, low-cost Faba_bean_130K targeted next-generation sequencing (TNGS) genotyping platform. A set of 69 Chinese faba bean accessions were genotyped with the TNGS genotyping platform, and the average mapping rate of captured reads to reference transcripts was 93.14%, of which 53.23% were located in the targeted regions. The TNGS genotyping results were validated by Sanger sequencing and the average consistency rate reached 93.6%. Comprehensive population genetic analysis was performed on the 69 Chinese faba bean accessions and identified four genetic subgroups correlated with the geographic distribution. This study provides valuable genomic resources and a reliable genotyping tool that could be implemented in genetic and molecular breeding studies to accelerate new cultivar development and improvement in faba bean.

Genomic structural equation modelling provides a whole-system approach for the future crop breeding

Using genomic structural equation modelling, this research demonstrates an efficient way to identify genetically correlating traits and provides an effective proxy for multi-trait selection to consider the joint genetic architecture of multiple interacting traits in crop breeding. Breeding crop cultivars with optimal value across multiple traits has been a challenge, as traits may negatively correlate due to pleiotropy or genetic linkage. For example, grain yield and grain protein content correlate negatively with each other in cereal crops. Future crop breeding needs to be based on practical yet accurate evaluation and effective selection of beneficial trait to retain genes with the best agronomic score for multiple traits. Here, we test the framework of whole-system-based approach using structural equation modelling (SEM) to investigate how one trait affects others to guide the optimal selection of a combination of agronomically important traits. Using ten traits and genome-wide SNP profiles from a worldwide barley panel and SEM analysis, we revealed a network of interacting traits, in which tiller number contributes positively to both grain yield and protein content; we further identified common genetic factors affecting multiple traits in the network of interaction. Our method demonstrates an efficient way to identify genetically correlating traits and underlying pleiotropic genetic factors and provides an effective proxy for multi-trait selection within a whole-system framework that considers the joint genetic architecture of multiple interacting traits in crop breeding. Our findings suggest the promise of a whole-system approach to overcome challenges such as the negative correlation of grain yield and protein content to facilitating quantitative and objective breeding decisions in future crop breeding.

A global barley panel revealing genomic signatures of breeding in modern Australian cultivars

The future of plant cultivar improvement lies in the evaluation of genetic resources from currently available germplasm. Today's gene pool of crop genetic diversity has been shaped during domestication and more recently by breeding. Recent efforts in plant breeding have been aimed at developing new and improved varieties from poorly adapted crops to suit local environments. However, the impact of these breeding efforts is poorly understood. Here, we assess the contributions of both historical and recent breeding efforts to local adaptation and crop improvement in a global barley panel by analysing the distribution of genetic variants with respect to geographic region or historical breeding category. By tracing the impact that breeding had on the genetic diversity of Hordeum vulgare (barley) released in Australia, where the history of barley production is relatively young, we identify 69 candidate regions within 922 genes that were under selection pressure. We also show that modern Australian barley varieties exhibit 12% higher genetic diversity than historical cultivars. Finally, field-trialling and phenotyping for agriculturally relevant traits across a diverse range of Australian environments suggests that genomic regions under strong breeding selection and their candidate genes are closely associated with key agronomic traits. In conclusion, our combined data set and germplasm collection provide a rich source of genetic diversity that can be applied to understanding and improving environmental adaptation and enhanced yields.

Genome-Wide Association Study and Identification of Candidate Genes for Nitrogen Use Efficiency in Barley (Hordeum vulgare L.)

Nitrogen (N) fertilizer is largely responsible for barley grain yield potential and quality, yet excessive application leads to environmental pollution and high production costs. Therefore, efficient use of N is fundamental for sustainable agriculture. In the present study, we investigated the performance of 282 barley accessions through hydroponic screening using optimal and low NH4NO3 treatments. Low-N treatment led to an average shoot dry weight reduction of 50%, but there were significant genotypic differences among the accessions. Approximately 20% of the genotypes showed high (>75%) relative shoot dry weight under low-N treatment and were classified as low-N tolerant, whereas 20% were low-N sensitive (≤55%). Low-N tolerant accessions exhibited well-developed root systems with an average increase of 60% in relative root dry weight to facilitate more N absorption. A genome-wide association study (GWAS) identified 66 significant marker trait associations (MTAs) conferring high nitrogen use efficiency, four of which were stable across experiments. These four MTAs were located on chromosomes 1H(1), 3H(1), and 7H(2) and were associated with relative shoot length, relative shoot and root dry weight. Genes corresponding to the significant MTAs were retrieved as candidate genes, including members of the asparagine synthetase gene family, several transcription factor families, protein kinases, and nitrate transporters. Most importantly, the high-affinity nitrate transporter 2.7 (HvNRT2.7) was identified as a promising candidate on 7H for root and shoot dry weight. The identified candidate genes provide new insights into our understanding of the molecular mechanisms driving nitrogen use efficiency in barley and represent potential targets for genetic improvement.

CRISPR-Cas9 enrichment and long read sequencing for fine mapping in plants

BACKGROUND

Genomic methods for identifying causative variants for trait loci applicable to a wide range of germplasm are required for plant biologists and breeders to understand the genetic control of trait variation.

RESULTS

We implemented Cas9-targeted sequencing for fine-mapping in apple, a method combining CRISPR-Cas9 targeted cleavage of a region of interest, followed by enrichment and long-read sequencing using the Oxford Nanopore Technology (ONT). We demonstrated the capability of this methodology to specifically cleave and enrich a plant genomic locus spanning 8 kb. The repeated mini-satellite motif located upstream of the Malus × domestica (apple) MYB10 transcription factor gene, causing red fruit colouration when present in a heterozygous state, was our exemplar to demonstrate the efficiency of this method: it contains a genomic region with a long structural variant normally ignored by short-read sequencing technologiesCleavage specificity of the guide RNAs was demonstrated using polymerase chain reaction products, before using them to specify cleavage of high molecular weight apple DNA. An enriched library was subsequently prepared and sequenced using an ONT MinION flow cell (R.9.4.1). Of the 7,056 ONT reads base-called using both Albacore2 (v2.3.4) and Guppy (v3.2.4), with a median length of 9.78 and 9.89 kb, respectively, 85.35 and 91.38%, aligned to the reference apple genome. Of the aligned reads, 2.98 and 3.04% were on-target with read depths of 180 × and 196 × for Albacore2 and Guppy, respectively, and only five genomic loci were off-target with read depth greater than 25 × , which demonstrated the efficiency of the enrichment method and specificity of the CRISPR-Cas9 cleavage.

CONCLUSIONS

We demonstrated that this method can isolate and resolve single-nucleotide and structural variants at the haplotype level in plant genomic regions. The combination of CRISPR-Cas9 target enrichment and ONT sequencing provides a more efficient technology for fine-mapping loci than genome-walking approaches.

The Potential of HTS Approaches for Accurate Genotyping in Grapevine (Vitis vinifera L.)

The main challenge associated with genotyping based on conventional length polymorphisms is the cross-laboratory standardization of allele sizes. This step requires the inclusion of standards and manual sizing to avoid false results. Capillary electrophoresis (CE) approaches limit the information to the length polymorphism and do not allow the determination of a complete marker sequence. As an alternative, high-throughput sequencing (HTS) offers complete information regarding marker sequences and their flanking regions. In this work, we investigated the suitability of a semi-quantitative sequencing approach for microsatellite genotyping using Illumina paired-end technology. Twelve microsatellite loci that are well established for grapevine CE typing were analysed on 96 grapevine samples from six different countries. We redesigned primers to the length of the amplicon for short sequencing (~100 bp). The primer pair was flanked with a 10 bp overhang for the introduction of barcodes on both sides of the amplicon to enable high multiplexing. The highest data peaks were determined as simple sequence repeat (SSR) alleles and compared with the CE dataset based on 12 reference samples. The comparison showed that HTS SSR genotyping can successfully replace the CE system in further experiments. We believe that, with next-generation sequencing, genotyping can be improved in terms of its speed, accuracy, and price.

Genome-wide association mapping reveals novel genes associated with coleoptile length in a worldwide collection of barley

BACKGROUND

Drought is projected to become more frequent and severe in a changing climate, which requires deep sowing of crop seeds to reach soil moisture. Coleoptile length is a key agronomic trait in cereal crops such as barley, as long coleoptiles are linked to drought tolerance and improved seedling establishment under early water-limited growing conditions.

RESULTS

In this study, we detected large genetic variation in a panel of 328 diverse barley (Hordeum vulgare L.) accessions. To understand the overall genetic basis of barley coleoptile length, all accessions were germinated in the dark and phenotyped for coleoptile length after 2 weeks. The investigated barleys had significant variation for coleoptile length. We then conducted genome-wide association studies (GWASs) with more than 30,000 molecular markers and identified 8 genes and 12 intergenic loci significantly associated with coleoptile length in our barley panel. The Squamosa promoter-binding-like protein 3 gene (SPL3) on chromosome 6H was identified as a major candidate gene. The missense variant on the second exon changed serine to alanine in the conserved SBP domain, which likely impacted its DNA-binding activity.

CONCLUSION

This study provides genetic loci for seedling coleoptile length along with candidate genes for future potential incorporation in breeding programmes to enhance early vigour and yield potential in water-limited environments.

Genome-Wide Association Study of Salinity Tolerance During Germination in Barley (Hordeum vulgare L.)

Barley seeds need to be able to germinate and establish seedlings in saline soils in Mediterranean-type climates. Despite being a major cereal crop, barley has few reported quantitative trait loci (QTL) and candidate genes underlying salt tolerance at the germination stage. Breeding programs targeting salinity tolerance at germination require an understanding of genetic loci and alleles in the current germplasm. In this study, we investigated seed-germination-related traits under control and salt stress conditions in 350 diverse barley accessions. A genome-wide association study, using ~24,000 genetic markers, was undertaken to detect marker-trait associations (MTA) and the underlying candidate genes for salinity tolerance during germination. We detected 19 loci containing 52 significant salt-tolerance-associated markers across all chromosomes, and 4 genes belonging to 4 family functions underlying the predicted MTAs. Our results provide new genetic resources and information to improve salt tolerance at germination in future barley varieties via genomic and marker-assisted selection and to open up avenues for further functional characterization of the identified candidate genes.

Gene-set association and epistatic analyses reveal complex gene interaction networks affecting flowering time in a worldwide barley collection

Single-marker genome-wide association studies (GWAS) have successfully detected associations between single nucleotide polymorphisms (SNPs) and agronomic traits such as flowering time and grain yield in barley. However, the analysis of individual SNPs can only account for a small proportion of genetic variation, and can only provide limited knowledge on gene network interactions. Gene-based GWAS approaches provide enormous opportunity both to combine genetic information and to examine interactions among genetic variants. Here, we revisited a previously published phenotypic and genotypic data set of 895 barley varieties grown in two years at four different field locations in Australia. We employed statistical models to examine gene-phenotype associations, as well as two-way epistasis analyses to increase the capability to find novel genes that have significant roles in controlling flowering time in barley. Genetic associations were tested between flowering time and corresponding genotypes of 174 putative flowering time-related genes. Gene-phenotype association analysis detected 113 genes associated with flowering time in barley, demonstrating the unprecedented power of gene-based analysis. Subsequent two-way epistasis analysis revealed 19 pairs of gene×gene interactions involved in controlling flowering time. Our study demonstrates that gene-based association approaches can provide higher capacity for future crop improvement to increase crop performance and adaptation to different environments.

Targeted enrichment by solution-based hybrid capture to identify genetic sequence variants in barley

In barley and other cereal crops, phenological diversity drives adaptation to different cultivation areas. Improvement of barley yield and quality traits requires adaptation to specific production areas with introgression of favorable alleles dependent upon precise identification of the underlying genes. Combining targeted sequence capture systems with next-generation sequencing provides an efficient approach to explore target genetic regions at high resolution, and allows rapid discovery of thousands of genetic polymorphisms. Here, we apply a versatile target-capture method to detect genome-wide polymorphisms in 174 flowering time-related genes, chosen based on prior knowledge from barley, rice, and Arabidopsis thaliana. Sequences were generated across a phenologically diverse panel of 895 barley varieties, resulting a high mean depth coverage of ~25x allowing reliable discovery and calling of insertion-deletion (InDel) and single nucleotide polymorphisms (SNPs). Sequences of InDel and SNPs from the targeted enrichment were utilized to develop 67 Kompetitive Allele Specific PCR (KASP) markers for validation. This work provides researchers and breeders a comprehensive molecular toolkit for the selection of phenology-related traits in barley.

Design Type(s)	genetic polymorphism analysis objective • strain comparison design
Measurement Type(s)	genetic variation
Technology Type(s)	DNA sequencing
Factor Type(s)	cultivar
Sample Characteristic(s)	Hordeum vulgare subsp. vulgare • leaf

Machine-accessible metadata file describing the reported data (ISA-Tab format)