Next generation breeding

Genome-wide association analysis of the resistance to infectious hematopoietic necrosis virus in two rainbow trout aquaculture lines confirms oligogenic architecture with several moderate effect quantitative trait loci

Infectious hematopoietic necrosis (IHN) is a disease of salmonid fish that is caused by the IHN virus (IHNV), which can cause substantial mortality and economic losses in rainbow trout aquaculture and fisheries enhancement hatchery programs. In a previous study on a commercial rainbow trout breeding line that has undergone selection, we found that genetic resistance to IHNV is controlled by the oligogenic inheritance of several moderate and many small effect quantitative trait loci (QTL). Here we used genome wide association analyses in two different commercial aquaculture lines that were naïve to previous exposure to IHNV to determine whether QTL were shared across lines, and to investigate whether there were major effect loci that were still segregating in the naïve lines. A total of 1,859 and 1,768 offspring from two commercial aquaculture strains were phenotyped for resistance to IHNV and genotyped with the rainbow trout Axiom 57K SNP array. Moderate heritability values (0.15-0.25) were estimated. Two statistical methods were used for genome wide association analyses in the two populations. No major QTL were detected despite the naïve status of the two lines. Further, our analyses confirmed an oligogenic architecture for genetic resistance to IHNV in rainbow trout. Overall, 17 QTL with notable effect (≥1.9% of the additive genetic variance) were detected in at least one of the two rainbow trout lines with at least one of the two statistical methods. Five of those QTL were mapped to overlapping or adjacent chromosomal regions in both lines, suggesting that some loci may be shared across commercial lines. Although some of the loci detected in this GWAS merit further investigation to better understand the biological basis of IHNV disease resistance across populations, the overall genetic architecture of IHNV resistance in the two rainbow trout lines suggests that genomic selection may be a more effective strategy for genetic improvement in this trait.

Population genomics analysis reveals footprints of selective breeding in a rapid-growth variety of Paulownia fortunei with apical dominance

Paulownia fortunei is an ecologically and economically valuable tree cultivated for its rapid growth and high-quality timber. To enhance Paulownia germplasm, we have developed the elite variety QingT with patented advantages in growth rate and apical dominance. To illuminate the genetic basis of QingT's superior traits, here we harness comparative population genomics to analyze genomic variation patterns between QingT and common Paulownia. We performed whole-genome re-sequencing of 30 QingT and 30 common samples, detecting 15.6 million SNPs and 2.6 million indels. Phylogeny and population structure analyses robustly partitioned common and QingT into distinct groups which indicate robust genome stabilization. QingT exhibited reduced heterozygosity and linkage disequilibrium decay compared to common Paulownia, reflecting high recombination, indicating hybridizing effects with common white-flowered string is the source of its patented advantages. Genome selection scans uncovered 25 regions of 169 genes with elevated nucleotide diversity, indicating selection sweeps among groups. Functional analysis of sweep genes revealed upregulation of ribosomal, biosynthesis, and growth pathways in QingT, implicating enhanced protein production and developmental processes in its rapid growth phenotype. This study's insights comprehensively chart genomic variation during Paulownia breeding, localizing candidate loci governing agronomic traits, and underpinnings of future molecular breeding efforts to boost productivity.

Guidelines for Performing CRISPR/Cas9 Genome Editing for Gene Validation and Trait Improvement in Crops

With the rapid advances in plant genome editing techniques over the past 10 years, more efficient and powerful crop genome editing applications are now possible. Candidate genes for key traits can be validated using CRISPR/Cas9-based knockouts and through the up- and down-regulation of gene expression. Likewise, new trait improvement approaches can take advantage of targeted editing to improve stress tolerance, disease resistance, and nutritional traits. However, several key steps in the process can prove tricky for researchers who might be new to plant genome editing. Here, we present step-by-step guidelines and best practices for a crop genome editing pipeline that should help to improve the rate of success. Important factors in the process include proper target sequence analysis and single guide RNA (sgRNA) design, sequencing of the target site in the genotypes of interest, performing an in vitro CRISPR/Cas9 ribonucleoprotein (RNP) assay to validate the designed sgRNAs, preparing the transformation constructs, considering a protoplast editing step as further validation, and, finally, stable plant transformation and mutation detection by Sanger and/or next-generation sequencing. With these detailed guidelines, a new user should be able to quickly set up a genome editing pipeline in their crop of interest and start making progress with the different CRISPR/Cas-based editing variants for gene validation and trait improvement purposes.

Diversity analyses in two ornamental and large-genome Ranunculaceae species based on a low-cost Klenow NGS-based protocol

Persian buttercup (Ranunculus asiaticus L.) and poppy anemone (Anemone coronaria L.) are ornamental, outcrossing, perennial species belonging to the Ranunculaceae family, characterized by large and highly repetitive genomes. We applied K-seq protocol in both species to generate high-throughput sequencing data and produce a large number of genetic polymorphisms. The technique entails the application of Klenow polymerase-based PCR using short primers designed by analyzing k-mer sets in the genome sequence. To date the genome sequence of both species has not been released, thus we designed primer sets based on the reference the genome sequence of the related species Aquilegia oxysepala var. kansuensis (Brühl). A whole of 11,542 SNPs were selected for assessing genetic diversity of eighteen commercial varieties of R. asiaticus, while 1,752 SNPs for assessing genetic diversity in six cultivars of A. coronaria. UPGMA dendrograms were constructed and in R. asiaticus integrated in with PCA analysis. This study reports the first molecular fingerprinting within Persian buttercup, while the results obtained in poppy anemone were compared with a previously published SSR-based fingerprinting, proving K-seq to be an efficient protocol for the genotyping of complex genetic backgrounds.

Comparative Metabolomic Analysis Reveals the Role of OsHPL1 in the Cold-Induced Metabolic Changes in Rice

Cytochrome P450 (CYP74) family members participate in the generation of oxylipins and play essential roles in plant adaptation. However, the metabolic reprogramming mediated by CYP74s under cold stress remains largely unexplored. Herein, we report how cold-triggered OsHPL1, a member of the CYP74 family, modulates rice metabolism. Cold stress significantly induced the expression of OsHPL1 and the accumulation of OPDA (12-oxo-phytodienoic acid) and jasmonates in the wild-type (WT) plants. The absence of OsHPL1 attenuates OPDA accumulation to a low temperature. Then, we performed a widely targeted metabolomics study covering 597 structurally annotated compounds. In the WT and hpl1 plants, cold stress remodeled the metabolism of lipids and amino acids. Although the WT and hpl1 mutants shared over one hundred cold-affected differentially accumulated metabolites (DAMs), some displayed distinct cold-responding patterns. Furthermore, we identified 114 and 56 cold-responding DAMs, specifically in the WT and hpl1 mutants. In conclusion, our work characterized cold-triggered metabolic rewiring and the metabolic role of OsHPL1 in rice.

Diversity of Phytosterols in Leaves of Wild Brassicaceae Species as Compared to Brassica napus Cultivars: Potential Traits for Insect Resistance and Abiotic Stress Tolerance

Phytosterols are natural compounds found in all higher plants that have a wide variety of roles in plant growth regulation and stress tolerance. The phytosterol composition can also influence the development and reproductive rate of strict herbivorous insects and other important agronomic traits such as temperature and drought tolerance in plants. In this study, we analysed the phytosterol composition in 18 Brassica napus (Rapeseed/canola) cultivars and 20 accessions belonging to 10 related wild Brassicaceae species to explore diverse and novel phytosterol profiles. Plants were grown in a controlled phytotron environment and their phytosterols were analysed using a saponification extraction method followed by GC-MS from the leaf samples. The B. napus cultivars showed slight diversity in eight phytosterols (>0.02%) due to the genotypic effect, whereas the wild accessions showed significant variability in their phytosterol profiles. Of interest, a number of wild accessions were found with high levels of campesterol (HIN20, HIN23, HUN27, HIN30, SARS2, and UPM6563), stigmasterol (UPM6813, UPM6563, ALBA17, and ALBA2), and isofucosterol (SARS12, SAR6, and DMU2). These changes in individual phytosterols, or ratios of phytosterols, can have a significant implication in plant tolerance to abiotic stress and plant insect resistance properties, which can be used in breeding for crop improvement.

A high-density integrated map for grapevine based on three mapping populations genotyped by the Vitis18K SNP chip

We present a high-density integrated map for grapevine, allowing refinement and improved understanding of the grapevine genome, while demonstrating the applicability of the Vitis18K SNP chip for linkage mapping. The improvement of grapevine through biotechnology requires identification of the molecular bases of target traits by studying marker-trait associations. The Vitis18K SNP chip provides a useful genotyping tool for genome-wide marker analysis. Most linkage maps are based on single mapping populations, but an integrated map can increase marker density and show order conservation. Here we present an integrated map based on three mapping populations. The parents consist of the well-known wine cultivars 'Cabernet Sauvignon', 'Corvina' and 'Rhine Riesling', the lesser-known wine variety 'Deckrot', and a table grape selection, G1-7720. Three high-density population maps with an average inter-locus gap ranging from 0.74 to 0.99 cM were developed. These maps show high correlations (0.9965-0.9971) with the reference assembly, containing only 93 markers with large order discrepancies compared to expected physical positions, of which a third is consistent across multiple populations. Moreover, the genetic data aid the further refinement of the grapevine genome assembly, by anchoring 104 yet unanchored scaffolds. From these population maps, an integrated map was constructed which includes 6697 molecular markers and reduces the inter-locus gap distance to 0.60 cM, resulting in the densest integrated map for grapevine thus far. A small number of discrepancies, mainly of short distance, involve 88 markers that remain conflictual across maps. The integrated map shows similar collinearity to the reference assembly (0.9974) as the single maps. This high-density map increases our understanding of the grapevine genome and provides a useful tool for its further characterization and the dissection of complex traits.

Genetic Determinants of Fiber-Associated Traits in Flax Identified by Omics Data Integration

In this paper, we explore potential genetic factors in control of flax phenotypes associated with fiber by mining a collection of 306 flax accessions from the Federal Research Centre of the Bast Fiber Crops, Torzhok, Russia. In total, 11 traits were assessed in the course of 3 successive years. A genome-wide association study was performed for each phenotype independently using six different single-locus models implemented in the GAPIT3 R package. Moreover, we applied a multivariate linear mixed model implemented in the GEMMA package to account for trait correlations and potential pleiotropic effects of polymorphisms. The analyses revealed a number of genomic variants associated with different fiber traits, implying the complex and polygenic control. All stable variants demonstrate a statistically significant allelic effect across all 3 years of the experiment. We tested the validity of the predicted variants using gene expression data available for the flax fiber studies. The results shed new light on the processes and pathways associated with the complex fiber traits, while the pinpointed candidate genes may be further used for marker-assisted selection.

AirMeasurer: open-source software to quantify static and dynamic traits derived from multiseason aerial phenotyping to empower genetic mapping studies in rice

Low-altitude aerial imaging, an approach that can collect large-scale plant imagery, has grown in popularity recently. Amongst many phenotyping approaches, unmanned aerial vehicles (UAVs) possess unique advantages as a consequence of their mobility, flexibility and affordability. Nevertheless, how to extract biologically relevant information effectively has remained challenging. Here, we present AirMeasurer, an open-source and expandable platform that combines automated image analysis, machine learning and original algorithms to perform trait analysis using 2D/3D aerial imagery acquired by low-cost UAVs in rice (Oryza sativa) trials. We applied the platform to study hundreds of rice landraces and recombinant inbred lines at two sites, from 2019 to 2021. A range of static and dynamic traits were quantified, including crop height, canopy coverage, vegetative indices and their growth rates. After verifying the reliability of AirMeasurer-derived traits, we identified genetic variants associated with selected growth-related traits using genome-wide association study and quantitative trait loci mapping. We found that the AirMeasurer-derived traits had led to reliable loci, some matched with published work, and others helped us to explore new candidate genes. Hence, we believe that our work demonstrates valuable advances in aerial phenotyping and automated 2D/3D trait analysis, providing high-quality phenotypic information to empower genetic mapping for crop improvement.

Integrating speed breeding with artificial intelligence for developing climate-smart crops

INTRODUCTION

In climate change, breeding crop plants with improved productivity, sustainability, and adaptability has become a daunting challenge to ensure global food security for the ever-growing global population. Correspondingly, climate-smart crops are also the need to regulate biomass production, which is imperative for the maintenance of ecosystem services worldwide. Since conventional breeding technologies for crop improvement are limited, time-consuming, and involve laborious selection processes to foster new and improved crop varieties. An urgent need is to accelerate the plant breeding cycle using artificial intelligence (AI) to depict plant responses to environmental perturbations in real-time.

MATERIALS AND METHODS

The review is a collection of authorized information from various sources such as journals, books, book chapters, technical bulletins, conference papers, and verified online contents.

CONCLUSIONS

Speed breeding has emerged as an essential strategy for accelerating the breeding cycles of crop plants by growing them under artificial light and temperature conditions. Furthermore, speed breeding can also integrate marker-assisted selection and cutting-edged gene-editing tools for early selection and manipulation of essential crops with superior agronomic traits. Scientists have recently applied next-generation AI to delve deeper into the complex biological and molecular mechanisms that govern plant functions under environmental cues. In addition, AIs can integrate, assimilate, and analyze complex OMICS data sets, an essential prerequisite for successful speed breeding protocol implementation to breed crop plants with superior yield and adaptability.

Developing and deploying an efficient genotyping workflow for accelerating maize improvement in developing countries

Background: Molecular breeding is an essential tool for accelerating genetic gain in crop improvement towards meeting the need to feed an ever-growing world population. Establishing low-cost, flexible genotyping platforms in small, public and regional laboratories can stimulate the application of molecular breeding in developing countries. These laboratories can serve plant breeding projects requiring low- to medium-density markers for marker-assisted selection (MAS) and quality control (QC) activities. Methods: We performed two QC and MAS experiments consisting of 637 maize lines, using an optimised genotyping workflow involving an in-house competitive allele-specific PCR (KASP) genotyping system with an optimised sample collection, preparation, and DNA extraction and quantitation process. A smaller volume of leaf-disc size plant samples was collected directly in 96-well plates for DNA extraction, using a slightly modified CTAB-based DArT DNA extraction protocol. DNA quality and quantity analyses were performed using a microplate reader, and the KASP genotyping and data analysis was performed in our laboratory. Results: Applying the optimized genotyping workflow expedited the QC and MAS experiments from over five weeks (when outsourcing) to two weeks and eliminated the shipping cost. Using a set of 28 KASP single nucleotide polymorphisms (SNPs) validated for maize, the QC experiment revealed the genetic identity of four maize varieties taken from five seed sources. Another set of 10 KASP SNPs was sufficient in verifying the parentage of 390 F1 lines. The KASP-based MAS was successfully applied to a maize pro-vitamin A (PVA) breeding program and for introgressing the aflatoxin resistance gene into elite tropical maize lines. Conclusion: This improved workflow has helped accelerate maize improvement activities of IITA's Maize Improvement Program and facilitated DNA fingerprinting for tracking improved crop varieties. National Agricultural Research Systems (NARS) in developing countries can adopt this workflow to fast-track molecular marker-based genotyping for crop improvement.

Developing and deploying an efficient genotyping workflow for accelerating maize improvement in developing countries

Background: Molecular breeding is an essential tool for accelerating genetic gain in crop improvement, towards meeting the need to feed an ever-growing world population. Establishing low-cost, flexible genotyping platforms in small, public and regional laboratories can stimulate the application of molecular breeding in developing countries. These laboratories can serve plant breeding projects requiring low- to medium-density markers for marker-assisted selection (MAS) and quality control (QC) activities. Methods: We performed two QC and MAS experiments consisting of 637 maize lines, using an optimised genotyping workflow involving an in-house competitive allele-specific PCR (KASP) genotyping system with an optimised sample collection, preparation, and DNA extraction and quantitation process. A smaller volume of leaf-disc size plant samples was collected directly in 96-well plates for DNA extraction, using a slightly modified CTAB-based DArT DNA extraction protocol. DNA quality and quantity analyses were performed using a microplate reader, and the KASP genotyping and data analysis was performed in our laboratory. Results: Applying the optimized genotyping workflow expedited the QC and MAS experiments from over five weeks (when outsourcing) to two weeks and eliminated the shipping cost. Using a set of 28 KASP single nucleotide polymorphisms (SNPs) validated for maize, the QC experiment revealed the genetic identity of four maize varieties taken from five seed sources. Another set of 10 KASP SNPs was sufficient in verifying the parentage of 390 F1 lines. The KASP-based MAS was successfully applied to a maize pro-vitamin A (PVA) breeding program and for introgressing the aflatoxin resistance gene into elite tropical maize lines. Conclusion: This improved workflow has helped accelerate maize improvement activities of IITA's Maize Improvement Program and facilitated DNA fingerprinting for tracking improved crop varieties. National Agricultural Research Systems (NARS) in developing countries can adopt this workflow to fast-track molecular marker-based genotyping for crop improvement.

Dissection of the practical soybean breeding pipeline by developing ZDX1, a high-throughput functional array

KEY MESSAGE

We developed the ZDX1 high-throughput functional soybean array for high accuracy evaluation and selection of both parents and progeny, which can greatly accelerate soybean breeding. Microarray technology facilitates rapid, accurate, and economical genotyping. Here, using resequencing data from 2214 representative soybean accessions, we developed the high-throughput functional array ZDX1, containing 158,959 SNPs, covering 90.92% of soybean genes and sites related to important traits. By application of the array, a total of 817 accessions were genotyped, including three subpopulations of candidate parental lines, parental lines and their progeny from practical breeding. The fixed SNPs were identified in progeny, indicating artificial selection during the breeding process. By identifying functional sites of target traits, novel soybean cyst nematode-resistant progeny and maturity-related novel sources were identified by allele combinations, demonstrating that functional sites provide an efficient method for the rapid screening of desirable traits or gene sources. Notably, we found that the breeding index (BI) was a good indicator for progeny selection. Superior progeny were derived from the combination of distantly related parents, with at least one parent having a higher BI. Furthermore, new combinations based on good performance were proposed for further breeding after excluding redundant and closely related parents. Genomic best linear unbiased prediction (GBLUP) analysis was the best analysis method and achieved the highest accuracy in predicting four traits when comparing SNPs in genic regions rather than whole genomic or intergenic SNPs. The prediction accuracy was improved by 32.1% by using progeny to expand the training population. Collectively, a versatile assay demonstrated that the functional ZDX1 array provided efficient information for the design and optimization of a breeding pipeline for accelerated soybean breeding.

In rice splice variants that restore the reading frame after frameshifting indel introduction are common, often induced by the indels and sometimes lead to organism-level rescue

The introduction of frameshifting non-3n indels enables the identification of gene-trait associations. However, it has been hypothesised that recovery of the original reading frame owing to usage of non-canonical splice forms could cause rescue. To date there is very little evidence for organism-level rescue by such a mechanism and it is unknown how commonly indels induce, or are otherwise associated with, frame-restoring splice forms. We perform CRISPR/Cas9 editing of randomly selected loci in rice to investigate these issues. We find that the majority of loci have a frame-restoring isoform. Importantly, three quarters of these isoforms are not seen in the absence of the indels, consistent with indels commonly inducing novel isoforms. This is supported by analysis in the context of NMD knockdowns. We consider in detail the two top rescue candidates, in wax deficient anther 1 (wda1) and brittle culm (bc10), finding that organismal-level rescue in both cases is strong but owing to different splice modification routes. More generally, however, as frame-restoring isoforms are low abundance and possibly too disruptive, such rescue we suggest to be the rare exception, not the rule. Nonetheless, assuming that indels commonly induce frame-restoring isoforms, these results emphasize the need to examine RNA level effects of non-3n indels and suggest that multiple non-3n indels in any given gene are advisable to probe a gene’s trait associations.

As protein coding genes are read in units of three (codons), insertions or deletions (indels) that are not a multiple of three long (non 3n) are expected to be especially harmful. Whether they are is important both for interpreting the results of non-3n indel experiments to probe a gene’s functional importance and for diagnostics. Particularly enigmatic are incidences where some non-3n changes in a gene compromise phenotypes while other seemingly comparable ones do not. One explanation for the latter is that a non-3n indel might be rescued via a frame-restoring splice form. Here we examine this hypothesis by inducing non-3n indels in many genes in rice and find that many non-3n indels are associated with a splice form that restores the reading frame. In the majority of these cases the indel appears to induce the potential rescuing splice form. We examine two top hit cases in detail and show functional rescue by splice modification. More generally, the frame-restoring forms are, however, low abundance and probably result in compromised proteins. We conclude then that splice mediated rescue is possible, but probably uncommon. Nonetheless it should not be overlooked in experimental design and interpretation.

Cross-species transcriptomic analyses reveals common and opposite responses in Arabidopsis, rice and barley following oxidative stress and hormone treatment

BACKGROUND

For translational genomics, a roadmap is needed to know the molecular similarities or differences between species, such as model species and crop species. This knowledge is invaluable for the selection of target genes and pathways to alter downstream in response to the same stimuli. Here, the transcriptomic responses to six treatments including hormones (abscisic acid - ABA and salicylic acid - SA); treatments that cause oxidative stress (3-amino-1,2,4-triazole - 3AT, methyl viologen - MV); inhibit respiration (antimycin A - AA) or induce genetic damage (ultraviolet radiation -UV) were analysed and compared between Arabidopsis (Arabidopsis thaliana), barley (Hordeum vulgare) and rice (Oryza sativa).

RESULTS

Common and opposite responses were identified between species, with the number of differentially expressed genes (DEGs) varying greatly between treatments and species. At least 70% of DEGs overlapped with at least one other treatment within a species, indicating overlapping response networks. Remarkably, 15 to 34% of orthologous DEGs showed opposite responses between species, indicating diversity in responses, despite orthology. Orthologous DEGs with common responses to multiple treatments across the three species were correlated with experimental data showing the functional importance of these genes in biotic/abiotic stress responses. The mitochondrial dysfunction response was revealed to be highly conserved in all three species in terms of responsive genes and regulation via the mitochondrial dysfunction element.

CONCLUSIONS

The orthologous DEGs that showed a common response between species indicate conserved transcriptomic responses of these pathways between species. However, many genes, including prominent salt-stress responsive genes, were oppositely responsive in multiple-stresses, highlighting fundamental differences in the responses and regulation of these genes between species. This work provides a resource for translation of knowledge or functions between species.

Developing and deploying an efficient genotyping workflow for accelerating maize improvement in developing countries

Background: Molecular breeding is an essential tool for accelerating genetic gain in crop improvement, towards meeting the need to feed an ever-growing world population. Establishing low-cost, flexible genotyping platforms in small, public and regional laboratories can stimulate the application of molecular breeding in developing countries. These laboratories can serve plant breeding projects requiring low- to medium-density markers for marker-assisted selection (MAS) and quality control (QC) activities. Methods: We performed two QC and MAS experiments consisting of 637 maize lines, using an optimised genotyping workflow involving an in-house competitive allele-specific PCR (KASP) genotyping system with an optimised sample collection, preparation, and DNA extraction and quantitation process. A smaller volume of leaf-disc size plant samples was collected directly in 96-well plates for DNA extraction, using a slightly modified CTAB-based DArT DNA extraction protocol. DNA quality and quantity analyses were performed using a microplate reader, and the KASP genotyping and data analysis was performed in our laboratory. Results: Applying the optimized genotyping workflow expedited the QC and MAS experiments from over five weeks (when outsourcing) to two weeks and eliminated the shipping cost. Using a set of 28 KASP single nucleotide polymorphisms (SNPs) validated for maize, the QC experiment revealed the genetic identity of four maize varieties taken from five seed sources. Another set of 10 KASP SNPs was sufficient in verifying the parentage of 390 F1 lines. The KASP-based MAS was successfully applied to a maize pro-vitamin A (PVA) breeding program and for introgressing the aflatoxin resistance gene into elite tropical maize lines. Conclusion: This improved workflow has helped accelerate maize improvement activities of IITA's Maize Improvement Program and facilitated DNA fingerprinting for tracking improved crop varieties. National Agricultural Research Systems (NARS) in developing countries can adopt this workflow to fast-track molecular marker-based genotyping for crop improvement.

Preparation and Curation of Phenotypic Datasets

Based on case studies, in this chapter we discuss the extent to which the number and identity of quantitative trait loci (QTL) identified from genome-wide association studies (GWAS) are affected by curation and analysis of phenotypic data. The chapter demonstrates through examples the impact of (1) cleaning of outliers, and of (2) the choice of statistical method for estimating genotypic mean values of phenotypic inputs in GWAS. No cleaning of outliers resulted in the highest number of dubious QTL, especially at loci with highly unbalanced allelic frequencies. A trade-off was identified between the risk of false positives and the risk of missing interesting, yet rare alleles. The choice of the statistical method to estimate genotypic mean values also affected the output of GWAS analysis, with reduced QTL overlap between methods. Using mixed models that capture spatial trends, among other features, increased the narrow-sense heritability of traits, the number of identified QTL and the overall power of GWAS analysis. Cleaning and choosing robust statistical models for estimating genotypic mean values should be included in GWAS pipelines to decrease both false positive and false negative rates of QTL detection.

A novel insight into transcriptional and epigenetic regulation underlying sex expression and flower development in melon (Cucumis melo L.)

Melon (Cucumis melo L.) is an important cucurbit and has been considered as a model plant for studying sex determination. The four most common sexual morphotypes in melon are monoecious (A-G-M), gynoecious (--ggM-), andromonoecious (A-G-mm), and hermaphrodite (--ggmm). Sex expression in melons is complex, as the genes and associated networks that govern the sex expression are not fully explored. Recently, RNA-seq transcriptomic profiling, ChIP-qPCR analysis integrated with gene ontology annotation and Kyoto Encyclopedia of Genes and Genomes pathways predicted the differentially expressed genes including sex-specific ACS and ACO genes, in regulating the sex-expression, phytohormonal cross-talk, signal transduction, and secondary metabolism in melons. Integration of transcriptional control through genetic interaction in between the ACS7, ACS11, and WIP1 in epistatic or hypostatic manner, along with the recruitment of H3K9ac and H3K27me3, epigenetically, overall determine sex expression. Alignment of protein sequences for establishing phylogenetic evolution, motif comparison, and protein-protein interaction supported the structural conservation while presence of the conserved hydrophilic and charged residues across the diverged evolutionary group predicted the functional conservation of the ACS protein. Presence of the putative cis-binding elements or DNA motifs, and its further comparison with DAP-seq-based cistrome and epicistrome of Arabidopsis, unraveled strong ancestry of melons with Arabidopsis. Motif comparison analysis also characterized putative genes and transcription factors involved in ethylene biosynthesis, signal transduction, and hormonal cross-talk related to sex expression. Overall, we have comprehensively reviewed research findings for a deeper insight into transcriptional and epigenetic regulation of sex expression and flower development in melons.

Weed competitive ability in wheat: a peek through in its functional significance, present status and future prospects

Weed competitive ability of a crop is one of the most widely explored aspects in the current scenario of aftermaths of synthetic herbicides such as herbicide resistant weeds emergence, residue accumulation in trophic levels; increased demands of organic produce, global climatic shifts, and other environmental issues. Further weed infestations are known to cause much more economic losses relative to crop attacks by pests. To understand the basic characteristics and underlying processes governing the competitive ability of a crop is therefore prudent, particularly in staples such as wheat. We discuss here an overview of the existing attributes of wheat-weed environment, the significance of crop competitiveness and various associated above-ground and below-ground traits (pertaining to early seed vigor and early seedling germination) discerned through biological, classical genetics and high throughput omics toolbox to provide numerous resources in terms of genome and transcriptome sequences, potential QTLs, genetic variation, molecular markers, association mapping studies, and others. Competitiveness is a cumulative response manifested as morphological, physiological, biochemical or allelochemical response ultimately driven through genetic architecture of a crop and its interaction with environment. Development of wheat competitive cultivar thus requires interdisciplinary approaches and germplasm screening to identify potential donors for competitiveness is an attractive and feasible alternative. For which utilization of landraces and other wild species, already proven to house sufficient genetic heterogeneity, thus poses a competitive advantage. Further, the availability of novel breeding techniques such as rapid generation advance could speed up the development of competitive wheat ideotype.

Combining UAV-RGB high-throughput field phenotyping and genome-wide association study to reveal genetic variation of rice germplasms in dynamic response to drought stress

Accurate and high-throughput phenotyping of the dynamic response of a large rice population to drought stress in the field is a bottleneck for genetic dissection and breeding of drought resistance. Here, high-efficiency and high-frequent image acquisition by an unmanned aerial vehicle (UAV) was utilized to quantify the dynamic drought response of a rice population under field conditions. Deep convolutional neural networks (DCNNs) and canopy height models were applied to extract highly correlated phenotypic traits including UAV-based leaf-rolling score (LRS_uav), plant water content (PWC_uav) and a new composite trait, drought resistance index by UAV (DRI_uav). The DCNNs achieved high accuracy (correlation coefficient R = 0.84 for modeling set and R = 0.86 for test set) to replace manual leaf-rolling rating. PWC_uav values were precisely estimated (correlation coefficient R = 0.88) and DRI_uav was modeled to monitor the drought resistance of rice accessions dynamically and comprehensively. A total of 111 significantly associated loci were detected by genome-wide association study for the three dynamic traits, and 30.6% of them were not detected in previous mapping studies using nondynamic drought response traits. Unmanned aerial vehicle and deep learning are confirmed effective phenotyping techniques for more complete genetic dissection of rice dynamic responses to drought and exploration of valuable alleles for drought resistance improvement.

Loss of Function of OsFBX267 and OsGA20ox2 in Rice Promotes Early Maturing and Semi-Dwarfism in γ-Irradiated IWP and Genome-Edited Pusa Basmati-1

Targeted mutagenesis is now becoming the most favored methodology to improve traits in popular rice cultivars selectively. Understanding the genetic basis of already available mutants could be the first step in designing such experiment. Improved White Ponni (IWP), a popularly grown South Indian rice variety, was subjected to γ irradiation to develop WP-22-2, an M₆ line superior in semi-dwarfism, early flowering, and high yield, and it has grain qualities similar to those of IWP. The exogenous application of gibberellic acid (GA₃) on WP-22-2 resulted in the elongation of shorter internodes to a level similar to IWP. The expression profiling of six genes regulating plant height showed their differential expression pattern at different time points post GA₃ treatment. Furthermore, the sequencing of WP-22-2 and IWP genomes revealed several single nucleotide polymorphisms (SNPs) and large-scale deletions in WP-22-2. The conversion of functional codons to stop codons was observed in OsGA20ox2 and OsFBX267, which have been reported to have roles in regulating semi-dwarfism and early flowering, respectively. The loss of function of OsGA20ox2 and OsFBX267 in WP-22-2 resulted in reduced plant height as well as early flowering, and the same has been confirmed by editing OsGA20ox2 in the rice variety Pusa Basmati1 (PB1) using the CRISPR-Cas9 approach. The targeted editing of OsGA20ox2 in PB1 conferred shorter plant height to the edited lines compared with the wild type. Altogether, the study provides evidence on mutating OsGA20ox2 and OsFBX267 genes to develop early maturing and semi-dwarf varieties that can be released to farmers after functional characterization and field trials.

Plant Variety Protection: Current Practices and Insights

Breeders persistently supply farmers with the best varieties in order to exceed consumer demand through plant-breeding processes that are resource-intensive. In order to motivate continuous innovation in variety development, a system needs to provide incentives for plant breeders to develop superior varieties, for example, exclusive ownership to produce and market those varieties. The most common system is the acquisition of intellectual property protection through plant variety protection, also known as the breeder’s right. Most countries have adopted the system established by the International Union for the Protection of New Varieties of Plants (UPOV). To be granted plant variety protection, the variety should prove to be unique by meeting three requirements: distinctness, uniformity, and stability. This review summarizes (1) the plant variety protection via UPOV convention, (2) technical methods for distinctness, uniformity, and stability testing via phenotype, molecular markers, and sequencing as well as their challenges and potentiality, and (3) additional discussions in essentially derived variety, value for cultivation and use testing, and open source seed initiative.

Chromosome-level reference genome of the soursop (Annona muricata): A new resource for Magnoliid research and tropical pomology

The flowering plant family Annonaceae includes important commercially grown tropical crops, but development of promising species is hindered by a lack of genomic resources to build breeding programs. Annonaceae are part of the magnoliids, an ancient lineage of angiosperms for which evolutionary relationships with other major clades remain unclear. To provide resources to breeders and evolutionary researchers, we report a chromosome-level genome assembly of the soursop (Annona muricata). We assembled the genome using 444.32 Gb of DNA sequences (676× sequencing depth) from PacBio and Illumina short-reads, in combination with 10× Genomics and Bionano data (v1). A total of 949 scaffolds were assembled to a final size of 656.77 Mb, with a scaffold N50 of 3.43 Mb (v1), and then further improved to seven pseudo-chromosomes using Hi-C sequencing data (v2; scaffold N50: 93.2 Mb, total size in chromosomes: 639.6 Mb). Heterozygosity was very low (0.06%), while repeat sequences accounted for 54.87% of the genome, and 23,375 protein-coding genes with an average of 4.79 exons per gene were annotated using de novo, RNA-seq and homology-based approaches. Reconstruction of the historical population size showed a slow continuous contraction, probably related to Cenozoic climate changes. The soursop is the first genome assembled in Annonaceae, supporting further studies of floral evolution in magnoliids, providing an essential resource for delineating relationships of ancient angiosperm lineages. Both genome-assisted improvement and conservation efforts will be strengthened by the availability of the soursop genome. As a community resource, this assembly will further strengthen the role of Annonaceae as model species for research on the ecology, evolution and domestication potential of tropical species in pomology and agroforestry.

Advancing knowledge of the plant nuclear periphery and its application for crop science

In this review, we explore recent advances in knowledge of the structure and dynamics of the plant nuclear envelope. As a paradigm, we focused our attention on the Linker of Nucleoskeleton and Cytoskeleton (LINC) complex, a structurally conserved bridging complex comprising SUN domain proteins in the inner nuclear membrane and KASH domain proteins in the outer nuclear membrane. Studies have revealed that this bridging complex has multiple functions with structural roles in positioning the nucleus within the cell, conveying signals across the membrane and organizing chromatin in the 3D nuclear space with impact on gene transcription. We also provide an up-to-date survey in nuclear dynamics research achieved so far in the model plant Arabidopsis thaliana that highlights its potential impact on several key plant functions such as growth, seed maturation and germination, reproduction and response to biotic and abiotic stress. Finally, we bring evidences that most of the constituents of the LINC Complex and associated components are, with some specificities, conserved in monocot and dicot crop species and are displaying very similar functions to those described for Arabidopsis. This leads us to suggest that a better knowledge of this system and a better account of its potential applications will in the future enhance the resilience and productivity of crop plants.

Genomic Selection in Tropical Forage Grasses: Current Status and Future Applications

The world population is expected to be larger and wealthier over the next few decades and will require more animal products, such as milk and beef. Tropical regions have great potential to meet this growing global demand, where pasturelands play a major role in supporting increased animal production. Better forage is required in consonance with improved sustainability as the planted area should not increase and larger areas cultivated with one or a few forage species should be avoided. Although, conventional tropical forage breeding has successfully released well-adapted and high-yielding cultivars over the last few decades, genetic gains from these programs have been low in view of the growing food demand worldwide. To guarantee their future impact on livestock production, breeding programs should leverage genotyping, phenotyping, and envirotyping strategies to increase genetic gains. Genomic selection (GS) and genome-wide association studies play a primary role in this process, with the advantage of increasing genetic gain due to greater selection accuracy, reduced cycle time, and increased number of individuals that can be evaluated. This strategy provides solutions to bottlenecks faced by conventional breeding methods, including long breeding cycles and difficulties to evaluate complex traits. Initial results from implementing GS in tropical forage grasses (TFGs) are promising with notable improvements over phenotypic selection alone. However, the practical impact of GS in TFG breeding programs remains unclear. The development of appropriately sized training populations is essential for the evaluation and validation of selection markers based on estimated breeding values. Large panels of single-nucleotide polymorphism markers in different tropical forage species are required for multiple application targets at a reduced cost. In this context, this review highlights the current challenges, achievements, availability, and development of genomic resources and statistical methods for the implementation of GS in TFGs. Additionally, the prediction accuracies from recent experiments and the potential to harness diversity from genebanks are discussed. Although, GS in TFGs is still incipient, the advances in genomic tools and statistical models will speed up its implementation in the foreseeable future. All TFG breeding programs should be prepared for these changes.

Exploring genetic architecture for pod-related traits in soybean using image-based phenotyping

Mature pod color (PC) and pod size (PS) served as important characteristics are used in the soybean breeding programs. However, manual phenotyping of such complex traits is time-consuming, laborious, and expensive for breeders. Here, we collected pod images from two different populations, namely, a soybean association panel (SAP) consisting of 187 accessions and an inter-specific recombinant inbred line (RIL) population containing 284 RILs. An image-based phenotyping method was developed and used to extract the pod color- and size-related parameters from images. Genome-wide association study (GWAS) and linkage mapping were performed to decipher the genetic control of pod color- and size-related traits across 2 successive years. Both populations exhibited wide phenotypic variations and continuous distribution in pod color- and size-related traits, indicating quantitative polygenic inheritance of these traits. GWAS and linkage mapping approaches identified the two major quantitative trait loci (QTL) underlying the pod color parameters, i.e., qPC3 and qPC19, located to chromosomes 3 and 19, respectively, and 12 stable QTLs for pod size-related traits across nine chromosomes. Several genes residing within the genomic region of stable QTL were identified as potential candidates underlying these pod-related traits based on the gene annotation and expression profiling data. Our results provide the useful information for fine-mapping/map-based cloning of QTL and marker-assisted selection of elite varieties with desirable pod traits.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01223-2.

Engineering crops of the future: CRISPR approaches to develop climate-resilient and disease-resistant plants

To meet increasing global food demand, breeders and scientists aim to improve the yield and quality of major food crops. Plant diseases threaten food security and are expected to increase because of climate change. CRISPR genome-editing technology opens new opportunities to engineer disease resistance traits. With precise genome engineering and transgene-free applications, CRISPR is expected to resolve the major challenges to crop improvement. Here, we discuss the latest developments in CRISPR technologies for engineering resistance to viruses, bacteria, fungi, and pests. We conclude by highlighting current concerns and gaps in technology, as well as outstanding questions for future research.

Photosynthesis in a Changing Global Climate: Scaling Up and Scaling Down in Crops

Photosynthesis is the major process leading to primary production in the Biosphere. There is a total of 7000bn tons of CO2 in the atmosphere and photosynthesis fixes more than 100bn tons annually. The CO2 assimilated by the photosynthetic apparatus is the basis of crop production and, therefore, of animal and human food. This has led to a renewed interest in photosynthesis as a target to increase plant production and there is now increasing evidence showing that the strategy of improving photosynthetic traits can increase plant yield. However, photosynthesis and the photosynthetic apparatus are both conditioned by environmental variables such as water availability, temperature, [CO2], salinity, and ozone. The "omics" revolution has allowed a better understanding of the genetic mechanisms regulating stress responses including the identification of genes and proteins involved in the regulation, acclimation, and adaptation of processes that impact photosynthesis. The development of novel non-destructive high-throughput phenotyping techniques has been important to monitor crop photosynthetic responses to changing environmental conditions. This wealth of data is being incorporated into new modeling algorithms to predict plant growth and development under specific environmental constraints. This review gives a multi-perspective description of the impact of changing environmental conditions on photosynthetic performance and consequently plant growth by briefly highlighting how major technological advances including omics, high-throughput photosynthetic measurements, metabolic engineering, and whole plant photosynthetic modeling have helped to improve our understanding of how the photosynthetic machinery can be modified by different abiotic stresses and thus impact crop production.

Recovery, Assessment, and Molecular Characterization of Minor Olive Genotypes in Tunisia

Olive is one of the oldest cultivated species in the Mediterranean Basin, including Tunisia, where it has a wide diversity, with more than 200 cultivars, of both wild and feral forms. Many minor cultivars are still present in marginal areas of Tunisia, where they are maintained by farmers in small local groves, but they are poorly characterized and evaluated. In order to recover this neglected germplasm, surveys were conducted in different areas, and 31 genotypes were collected, molecularly characterized with 12 nuclear microsatellite (simple sequence repeat (SSR)) markers, and compared with 26 reference cultivars present in the Tunisian National Olive collection. The analysis revealed an overall high genetic diversity of this olive's germplasm, but also discovered the presence of synonymies and homonymies among the commercialized varieties. The structure analysis showed the presence of different gene pools in the analyzed germplasm. In particular, the marginal germplasm from Ras Jbal and Azmour is characterized by gene pools not present in commercial (Nurseries) varieties, pointing out the very narrow genetic base of the commercialized olive material in Tunisia, and the need to broaden it to avoid the risk of genetic erosion of this species in this country.

Reduced fire blight susceptibility in apple cultivars using a high-efficiency CRISPR/Cas9-FLP/FRT-based gene editing system

The bacterium Erwinia amylovora, the causal agent of fire blight disease in apple, triggers its infection through the DspA/E effector which interacts with the apple susceptibility protein MdDIPM4. In this work, MdDIPM4 knockout has been produced in two Malus × domestica susceptible cultivars using the CRISPR/Cas9 system delivered via Agrobacterium tumefaciens. Fifty-seven transgenic lines were screened to identify CRISPR/Cas9-induced mutations. An editing efficiency of 75% was obtained. Seven edited lines with a loss-of-function mutation were inoculated with the pathogen. Highly significant reduction in susceptibility was observed compared to control plants. Sequencing of five potential off-target sites revealed no mutation event. Moreover, our construct contained a heat-shock inducible FLP/FRT recombination system designed specifically to remove the T-DNA harbouring the expression cassettes for CRISPR/Cas9, the marker gene and the FLP itself. Six plant lines with reduced susceptibility to the pathogen were heat-treated and screened by real-time PCR to quantify the exogenous DNA elimination. The T-DNA removal was further validated by sequencing in one plant line. To our knowledge, this work demonstrates for the first time the development and application of a CRISPR/Cas9-FLP/FRT gene editing system for the production of edited apple plants carrying a minimal trace of exogenous DNA.

Photosynthetic Metabolism under Stressful Growth Conditions as a Bases for Crop Breeding and Yield Improvement

Increased periods of water shortage and higher temperatures, together with a reduction in nutrient availability, have been proposed as major factors that negatively impact plant development. Photosynthetic CO2 assimilation is the basis of crop production for animal and human food, and for this reason, it has been selected as a primary target for crop phenotyping/breeding studies. Within this context, knowledge of the mechanisms involved in the response and acclimation of photosynthetic CO2 assimilation to multiple changing environmental conditions (including nutrients, water availability, and rising temperature) is a matter of great concern for the understanding of plant behavior under stress conditions, and for the development of new strategies and tools for enhancing plant growth in the future. The current review aims to analyze, from a multi-perspective approach (ranging across breeding, gas exchange, genomics, etc.) the impact of changing environmental conditions on the performance of the photosynthetic apparatus and, consequently, plant growth.

Strong temporal dynamics of QTL action on plant growth progression revealed through high-throughput phenotyping in canola

A major challenge of plant biology is to unravel the genetic basis of complex traits. We took advantage of recent technical advances in high-throughput phenotyping in conjunction with genome-wide association studies to elucidate genotype-phenotype relationships at high temporal resolution. A diverse Brassica napus population from a commercial breeding programme was analysed by automated non-invasive phenotyping. Time-resolved data for early growth-related traits, including estimated biovolume, projected leaf area, early plant height and colour uniformity, were established and complemented by fresh and dry weight biomass. Genome-wide SNP array data provided the framework for genome-wide association analyses. Using time point data and relative growth rates, multiple robust main effect marker-trait associations for biomass and related traits were detected. Candidate genes involved in meristem development, cell wall modification and transcriptional regulation were detected. Our results demonstrate that early plant growth is a highly complex trait governed by several medium and many small effect loci, most of which act only during short phases. These observations highlight the importance of taking the temporal patterns of QTL/allele actions into account and emphasize the need for detailed time-resolved analyses to effectively unravel the complex and stage-specific contributions of genes affecting growth processes that operate at different developmental phases.

A Systematic Review of Durum Wheat: Enhancing Production Systems by Exploring Genotype, Environment, and Management (G × E × M) Synergies

According to the UN-FAO, agricultural production must increase by 50% by 2050 to meet global demand for food. This goal can be accomplished, in part, by the development of improved cultivars coupled with modern best management practices. Overall, wheat production on farms will have to increase significantly to meet future demand, and in the face of a changing climate that poses risk to even current rates of production. Durum wheat [Triticum turgidum L. ssp. durum (Desf.)] is used largely for pasta, couscous and bulgur production. Durum producers face a range of factors spanning abiotic (frost damage, drought, and sprouting) and biotic (weed, disease, and insect pests) stresses that impact yields and quality specifications desired by export market end-users. Serious biotic threats include Fusarium head blight (FHB) and weed pest pressures, which have increased as a result of herbicide resistance. While genetic progress for yield and quality is on pace with common wheat (Triticum aestivum L.), development of resistant durum cultivars to FHB is still lagging. Thus, successful biotic and abiotic threat mitigation are ideal case studies in Genotype (G) × Environment (E) × Management (M) interactions where superior cultivars (G) are grown in at-risk regions (E) and require unique approaches to management (M) for sustainable durum production. Transformational approaches to research are needed in order for agronomists, breeders and durum producers to overcome production constraints. Designing robust agronomic systems for durum demands scientific creativity and foresight based on a deep understanding of constitutive components and their innumerable interactions with each other and the environment. This encompasses development of durum production systems that suit specific agro-ecozones and close the yield gap between genetic potential and on-farm achieved yield. Advances in individual technologies (e.g., genetic improvements, new pesticides, seeding technologies) are of little benefit until they are melded into resilient G × E × M systems that will flourish in the field under unpredictable conditions of prairie farmlands. We explore how recent genetic progress and selected management innovations can lead to a resilient and transformative durum production system.

Insights into deployment of DNA markers in plant variety protection and registration

The efficiency of phenotype-based assessments of plant variety protection and registration could be improved by the integration of DNA-based testing. We review the current and proposed models in the era of next-generation breeding. The current plant variety protection system relies on morphological description of plant varieties. Distinctness, uniformity, and stability (DUS) assessments determine whether a new variety is distinguishable from common knowledge varieties and exhibits sufficient phenotypic uniformity and stability during two independent growing cycles. However, DUS assessment can be costly, time-consuming and often restricted to a relatively small number of traits that can be influenced by environmental conditions. This calls for the adoption of a DNA-based system which is endorsed by the International Union for the Protection of New Varieties of Plants (UPOV). This could enable examiners to deploy trait-specific DNA markers in DUS testing as well as using such genetic markers to manage reference collections. Within UPOV's system, breeders can freely use protected varieties in breeding programs. However, breeders of protected varieties may seek sharing in ownership of essentially derived varieties once it is proven that they, with the exception of a few distinctive DUS trait(s), conform to parental varieties in essential characteristics. As well as their complementary role in DUS testing, DNA markers have been known as a good replacement of morphological traits in defining boundaries between independently and essentially derived varieties. With the advent of new breeding technologies that allow minor modification in varieties with outcomes of specific merit or utility, detecting distinctness between varieties may become increasingly challenging. This, together with the ever-increasing number of varieties with which to compare new candidate varieties, supports the potential utility of using DNA-based approaches in variety description.

Genome-wide mapping of adult plant stripe rust resistance in wheat cultivar Toni

Two adult plant stripe rust resistance QTL, QYrto.swust-3AS and QYrto.swust-3BS, were identified and mapped in common wheat cultivar Toni. The two QTL were located to corresponding positions in the wheat physical map position based on flanking SNP markers. Stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), is one of the most important foliar diseases of wheat. Characterization and utilization of resistance genes are the most effective, economic and environmental-friendly way to control the disease. The wheat cultivar Toni resistant at the adult plant stage to predominant Chinese Pst races was crossed with the susceptible genotype Mingxian 169. A recombinant inbred line population comprising 171 lines was tested in the field at three locations in the 2016 and 2017 crop seasons. The Affymetrix Axiom^® 35 K single-nucleotide polymorphism (SNP) Wheat Breeder's Genotyping Array was used to map quantitative trait loci (QTL) for adult plant resistance to stripe rust. Inclusive composite interval mapping identified stable QTL QYrto.swust-3AS and QYrto.swust-3BS that explained 31.6-48.2% and 21.9-56.3% of the variation in stripe rust severity and infection type, respectively. The two QTL regions were anchored to the wheat IWGSC Ref Seq v1.0 sequence. QYrto.swust-3AS was localized to a 2.22-Mb interval flanked by SNP markers AX-95240191 and AX-94828890. Among 65 HC (high confidence) annotated genes in this region, 11 (16.9%) contained NB-ARC domains and 9 (13.8%) contained protein kinase domains and thus could contribute to disease resistance. QYrto.swust-3BS was localized to a 4.77-Mb interval flanked by SNP markers AX-94509749 and AX-94998050. One hundred and thirty three HC genes are annotated in this region. Among them, 14 (10.5%) protein kinase domain genes may contribute to disease resistance. The linked markers should be useful for marker-assisted selection in breeding for resistance.

Identification and Characterization of EI (Elongated Internode) Gene in Tomato (Solanum lycopersicum)

Internode length is an important agronomic trait affecting plant architecture and crop yield. However, few genes for internode elongation have been identified in tomato. In this study, we characterized an elongated internode inbred line P502, which is a natural mutant of the tomato cultivar 05T606. The mutant P502 exhibits longer internode and higher bioactive GA concentration compared with wild-type 05T606. Genetic analysis suggested that the elongated internode trait is controlled by quantitative trait loci (QTL). Then, we identified a major QTL on chromosome 2 based on molecular markers and bulked segregant analysis (BSA). The locus was designated as EI (Elongated Internode), which explained 73.6% genetic variance. The EI was further mapped to a 75.8-kb region containing 10 genes in the reference Heinz 1706 genome. One single nucleotide polymorphism (SNP) in the coding region of solyc02g080120.1 was identified, which encodes gibberellin 2-beta-dioxygenase 7 (SlGA2ox7). SlGA2ox7, orthologous to AtGA2ox7 and AtGA2ox8, is involved in the regulation of GA degradation. Overexpression of the wild EI gene in mutant P502 caused a dwarf phenotype with a shortened internode. The difference of EI expression levels was not significant in the P502 and wild-type, but the expression levels of GA biosynthetic genes including CPS, KO, KAO, GA20ox1, GA20ox2, GA20ox4, GA3ox1, GA2ox1, GA2ox2, GA2ox4, and GA2ox5, were upregulated in mutant P502. Our results may provide a better understanding of the genetics underlying the internode elongation and valuable information to improve plant architecture of the tomato.

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

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

The evolving EU regulatory framework for precision breeding

Plant breeding has always relied on progress in various scientific disciplines to generate and enable access to genetic variation. Until the 1970s, available techniques generated mostly random genetic alterations that were subject to a selection procedure in the plant material. Recombinant nucleic acid technology, however, started a new era of targeted genetic alterations, or precision breeding, enabling a much more targeted approach to trait management. More recently, developments in genome editing are now providing yet more control by enabling alterations at exact locations in the genome. The potential of recombinant nucleic acid technology fueled discussions about potentially new associated risks and, starting in the late 1980s, biosafety legislation for genetically modified organisms (GMOs) has developed in the European Union. However, the last decade has witnessed a lot of discussions as to whether or not genome editing and other precision breeding techniques should be encompassed by the EU GMO legislation. A recent ruling from the Court of Justice of the European Union indicated that directed mutagenesis techniques should be subject to the provisions of the GMO Directive, essentially putting many precision breeding techniques in the same regulatory basket. This review outlines the evolving EU regulatory framework for GMOs and discusses some potential routes that the EU may take for the regulation of precision breeding.

Enhancing the rate of genetic gain in public-sector plant breeding programs: lessons from the breeder's equation

The integration of new technologies into public plant breeding programs can make a powerful step change in agricultural productivity when aligned with principles of quantitative and Mendelian genetics. The breeder's equation is the foundational application of quantitative genetics to crop improvement. Guided by the variables that describe response to selection, emerging breeding technologies can make a powerful step change in the effectiveness of public breeding programs. The most promising innovations for increasing the rate of genetic gain without greatly increasing program size appear to be related to reducing breeding cycle time, which is likely to require the implementation of parent selection on non-inbred progeny, rapid generation advance, and genomic selection. These are complex processes and will require breeding organizations to adopt a culture of continuous optimization and improvement. To enable this, research managers will need to consider and proactively manage the, accountability, strategy, and resource allocations of breeding teams. This must be combined with thoughtful management of elite genetic variation and a clear separation between the parental selection process and product development and advancement process. With an abundance of new technologies available, breeding teams need to evaluate carefully the impact of any new technology on selection intensity, selection accuracy, and breeding cycle length relative to its cost of deployment. Finally breeding data management systems need to be well designed to support selection decisions and novel approaches to accelerate breeding cycles need to be routinely evaluated and deployed.

Progress in plant genome sequencing: research directions

Since the first plant genome of Arabidopsis thaliana has been sequenced and published, genome sequencing technologies have undergone significant changes. New algorithms, sequencing technologies and bioinformatic approaches were adopted to obtain genome, transcriptome and exome sequences for model and crop species, which have permitted deep inferences into plant biology. As a result of an improved genome assembly and analysis methods, genome sequencing costs plummeted and the number of high-quality plant genome sequences is constantly growing. Consequently, more than 300 plant genome sequences have been published over the past twenty years. Although many of the published genomes are considered incomplete, they proved to be a valuable tool for identifying genes involved in the formation of economically valuable plant traits, for marker-assisted and genomic selection and for comparative analysis of plant genomes in order to determine the basic patterns of origin of various plant species. Since a high coverage and resolution of a genome sequence is not enough to detect all changes in complex samples, targeted sequencing, which consists in the isolation and sequencing of a specific region of the genome, has begun to develop. Targeted sequencing has a higher detection power (the ability to identify new differences/variants) and resolution (up to one basis). In addition, exome sequencing (the method of sequencing only protein-coding genes regions) is actively developed, which allows for the sequencing of non-expressed alleles and genes that cannot be found with RNA-seq. In this review, an analysis of sequencing technologies development and the construction of “reference” genomes of plants is performed. A comparison of the methods of targeted sequencing based on the use of the reference DNA sequence is accomplished.

Converging phenomics and genomics to study natural variation in plant photosynthetic efficiency

In recent years developments in plant phenomic approaches and facilities have gradually caught up with genomic approaches. An opportunity lies ahead to dissect complex, quantitative traits when both genotype and phenotype can be assessed at a high level of detail. This is especially true for the study of natural variation in photosynthetic efficiency, for which forward genetics studies have yielded only a little progress in our understanding of the genetic layout of the trait. High-throughput phenotyping, primarily from chlorophyll fluorescence imaging, should help to dissect the genetics of photosynthesis at the different levels of both plant physiology and development. Specific emphasis should be directed towards understanding the acclimation of the photosynthetic machinery in fluctuating environments, which may be crucial for the identification of genetic variation for relevant traits in food crops. Facilities should preferably be designed to accommodate phenotyping of photosynthesis-related traits in such environments. The use of forward genetics to study the genetic architecture of photosynthesis is likely to lead to the discovery of novel traits and/or genes that may be targeted in breeding or bio-engineering approaches to improve crop photosynthetic efficiency. In the near future, big data approaches will play a pivotal role in data processing and streamlining the phenotype-to-gene identification pipeline.

Sunflower pan-genome analysis shows that hybridization altered gene content and disease resistance

Domesticated plants and animals often display dramatic responses to selection, but the origins of the genetic diversity underlying these responses remain poorly understood. Despite domestication and improvement bottlenecks, the cultivated sunflower remains highly variable genetically, possibly due to hybridization with wild relatives. To characterize genetic diversity in the sunflower and to quantify contributions from wild relatives, we sequenced 287 cultivated lines, 17 Native American landraces and 189 wild accessions representing 11 compatible wild species. Cultivar sequences failing to map to the sunflower reference were assembled de novo for each genotype to determine the gene repertoire, or 'pan-genome', of the cultivated sunflower. Assembled genes were then compared to the wild species to estimate origins. Results indicate that the cultivated sunflower pan-genome comprises 61,205 genes, of which 27% vary across genotypes. Approximately 10% of the cultivated sunflower pan-genome is derived through introgression from wild sunflower species, and 1.5% of genes originated solely through introgression. Gene ontology functional analyses further indicate that genes associated with biotic resistance are over-represented among introgressed regions, an observation consistent with breeding records. Analyses of allelic variation associated with downy mildew resistance provide an example in which such introgressions have contributed to resistance to a globally challenging disease.

Regulation and Evolution of NLR Genes: A Close Interconnection for Plant Immunity

NLR (NOD-like receptor) genes belong to one of the largest gene families in plants. Their role in plants' resistance to pathogens has been clearly described for many members of this gene family, and dysregulation or overexpression of some of these genes has been shown to induce an autoimmunity state that strongly affects plant growth and yield. For this reason, these genes have to be tightly regulated in their expression and activity, and several regulatory mechanisms are described here that tune their gene expression and protein levels. This gene family is subjected to rapid evolution, and to maintain diversity at NLRs, a plethora of genetic mechanisms have been identified as sources of variation. Interestingly, regulation of gene expression and evolution of this gene family are two strictly interconnected aspects. Indeed, some examples have been reported in which mechanisms of gene expression regulation have roles in promotion of the evolution of this gene family. Moreover, co-evolution of the NLR gene family and other gene families devoted to their control has been recently demonstrated, as in the case of miRNAs.

CRISPR/Cas9-mediated efficient targeted mutagenesis in grape in the first generation

The clustered regularly interspaced short palindromic repeats-associated protein 9 (CRISPR/Cas9) system is a powerful tool for editing plant genomes. Efficient genome editing of grape (Vitis vinifera) suspension cells using the type II CRISPR/Cas9 system has been demonstrated; however, it has not been established whether this system can be applied to get biallelic mutations in the first generation of grape. In this current study, we designed four guide RNAs for the VvWRKY52 transcription factor gene for using with the CRISPR/Cas9 system, and obtained transgenic plants via Agrobacterium-mediated transformation, using somatic embryos of the Thompson Seedless cultivar. Analysis of the first-generation transgenic plants verified 22 mutant plants of the 72 T-DNA-inserted plants. Of these, 15 lines carried biallelic mutations and seven were heterozygous. A range of RNA-guided editing events, including large deletions, were found in the mutant plants, while smaller deletions comprised the majority of the detected mutations. Sequencing of potential off-target sites for all four targets revealed no off-target events. In addition, knockout of VvWRKY52 in grape increased the resistance to Botrytis cinerea. We conclude that the CRISPR/Cas9 system allows precise genome editing in the first generation of grape and represents a useful tool for gene functional analysis and grape molecular breeding.

A practical introduction to Random Forest for genetic association studies in ecology and evolution

Large genomic studies are becoming increasingly common with advances in sequencing technology, and our ability to understand how genomic variation influences phenotypic variation between individuals has never been greater. The exploration of such relationships first requires the identification of associations between molecular markers and phenotypes. Here, we explore the use of Random Forest (RF), a powerful machine-learning algorithm, in genomic studies to discern loci underlying both discrete and quantitative traits, particularly when studying wild or nonmodel organisms. RF is becoming increasingly used in ecological and population genetics because, unlike traditional methods, it can efficiently analyse thousands of loci simultaneously and account for nonadditive interactions. However, understanding both the power and limitations of Random Forest is important for its proper implementation and the interpretation of results. We therefore provide a practical introduction to the algorithm and its use for identifying associations between molecular markers and phenotypes, discussing such topics as data limitations, algorithm initiation and optimization, as well as interpretation. We also provide short R tutorials as examples, with the aim of providing a guide to the implementation of the algorithm. Topics discussed here are intended to serve as an entry point for molecular ecologists interested in employing Random Forest to identify trait associations in genomic data sets.

Decision Tree Algorithm-Generated Single-Nucleotide Polymorphism Barcodes of rbcL Genes for 38 Brassicaceae Species Tagging

DNA barcode sequences are accumulating in large data sets. A barcode is generally a sequence larger than 1000 base pairs and generates a computational burden. Although the DNA barcode was originally envisioned as straightforward species tags, the identification usage of barcode sequences is rarely emphasized currently. Single-nucleotide polymorphism (SNP) association studies provide us an idea that the SNPs may be the ideal target of feature selection to discriminate between different species. We hypothesize that SNP-based barcodes may be more effective than the full length of DNA barcode sequences for species discrimination. To address this issue, we tested a ribulose diphosphate carboxylase (rbcL) SNP barcoding (RSB) strategy using a decision tree algorithm. After alignment and trimming, 31 SNPs were discovered in the rbcL sequences from 38 Brassicaceae plant species. In the decision tree construction, these SNPs were computed to set up the decision rule to assign the sequences into 2 groups level by level. After algorithm processing, 37 nodes and 31 loci were required for discriminating 38 species. Finally, the sequence tags consisting of 31 rbcL SNP barcodes were identified for discriminating 38 Brassicaceae species based on the decision tree-selected SNP pattern using RSB method. Taken together, this study provides the rational that the SNP aspect of DNA barcode for rbcL gene is a useful and effective sequence for tagging 38 Brassicaceae species.

Trans-species synthetic gene design allows resistance pyramiding and broad-spectrum engineering of virus resistance in plants

To infect plants, viruses rely heavily on their host's machinery. Plant genetic resistances based on host factor modifications can be found among existing natural variability and are widely used for some but not all crops. While biotechnology can supply for the lack of natural resistance alleles, new strategies need to be developed to increase resistance spectra and durability without impairing plant development. Here, we assess how the targeted allele modification of the Arabidopsis thaliana translation initiation factor eIF4E1 can lead to broad and efficient resistance to the major group of potyviruses. A synthetic Arabidopsis thaliana eIF4E1 allele was designed by introducing multiple amino acid changes associated with resistance to potyvirus in naturally occurring Pisum sativum alleles. This new allele encodes a functional protein while maintaining plant resistance to a potyvirus isolate that usually hijacks eIF4E1. Due to its biological functionality, this synthetic allele allows, at no developmental cost, the pyramiding of resistances to potyviruses that selectively use the two major translation initiation factors, eIF4E1 or its isoform eIFiso4E. Moreover, this combination extends the resistance spectrum to potyvirus isolates for which no efficient resistance has so far been found, including resistance-breaking isolates and an unrelated virus belonging to the Luteoviridae family. This study is a proof-of-concept for the efficiency of gene engineering combined with knowledge of natural variation to generate trans-species virus resistance at no developmental cost to the plant. This has implications for breeding of crops with broad-spectrum and high durability resistance using recent genome editing techniques.

Comparison of Marker-Based Genomic Estimated Breeding Values and Phenotypic Evaluation for Selection of Bacterial Spot Resistance in Tomato

Bacterial spot affects tomato crops (Solanum lycopersicum) grown under humid conditions. Major genes and quantitative trait loci (QTL) for resistance have been described, and multiple loci from diverse sources need to be combined to improve disease control. We investigated genomic selection (GS) prediction models for resistance to Xanthomonas euvesicatoria and experimentally evaluated the accuracy of these models. The training population consisted of 109 families combining resistance from four sources and directionally selected from a population of 1,100 individuals. The families were evaluated on a plot basis in replicated inoculated trials and genotyped with single nucleotide polymorphisms (SNP). We compared the prediction ability of models developed with 14 to 387 SNP. Genomic estimated breeding values (GEBV) were derived using Bayesian least absolute shrinkage and selection operator regression (BL) and ridge regression (RR). Evaluations were based on leave-one-out cross validation and on empirical observations in replicated field trials using the next generation of inbred progeny and a hybrid population resulting from selections in the training population. Prediction ability was evaluated based on correlations between GEBV and phenotypes (r_g), percentage of coselection between genomic and phenotypic selection, and relative efficiency of selection (r_g/r_p). Results were similar with BL and RR models. Models using only markers previously identified as significantly associated with resistance but weighted based on GEBV and mixed models with markers associated with resistance treated as fixed effects and markers distributed in the genome treated as random effects offered greater accuracy and a high percentage of coselection. The accuracy of these models to predict the performance of progeny and hybrids exceeded the accuracy of phenotypic selection.

The first identification of genomic loci in plants associated with resistance to galling insects: a case study in Eucalyptus L'Hér. (Myrtaceae)

Genomic loci related with resistance to gall-inducing insects have not been identified in any plants. Here, association mapping was used to identify molecular markers for resistance to the gall wasp Leptocybe invasa in two Eucalyptus species. A total of 86 simple sequence repeats (SSR) markers were screened out from 839 SSRs and used for association mapping in E. grandis. By applying the mixed linear model, seven markers were identified to be associated significantly (P ≤ 0.05) with the gall wasp resistance in E. grandis, including two validated with a correction of permutation test (P ≤ 0.008). The proportion of the variance in resistance explained by a significant marker ranged from 3.3% to 37.8%. Four out of the seven significant associations in E. grandis were verified and also validated (P ≤ 0.073 in a permutation test) in E. tereticornis, with the variation explained ranging from 24.3% to 48.5%. Favourable alleles with positive effect were also mined from the significant markers in both species. These results provide insight into the genetic control of gall wasp resistance in plants and have great potential for marker-assisted selection for resistance to L. invasa in the important tree genus Eucalyptus.