Structural variants in the soybean genome localize to clusters of biotic stress-response genes

Universal signatures of transposable element compartmentalization across eukaryotic genomes

The evolutionary mechanisms that drive the emergence of genome architecture remain poorly understood but can now be assessed with unprecedented power due to the massive accumulation of genome assemblies spanning phylogenetic diversity1,2. Transposable elements (TEs) are a rich source of large-effect mutations since they directly and indirectly drive genomic structural variation and changes in gene expression3. Here, we demonstrate universal patterns of TE compartmentalization across eukaryotic genomes spanning ~1.7 billion years of evolution, in which TEs colocalize with gene families under strong predicted selective pressure for dynamic evolution and involved in specific functions. For non-pathogenic species these genes represent families involved in defense, sensory perception and environmental interaction, whereas for pathogenic species, TE-compartmentalized genes are highly enriched for pathogenic functions. Many TE-compartmentalized gene families display signatures of positive selection at the molecular level. Furthermore, TE-compartmentalized genes exhibit an excess of high-frequency alleles for polymorphic TE insertions in fruit fly populations. We postulate that these patterns reflect selection for adaptive TE insertions as well as TE-associated structural variants. This process may drive the emergence of a shared TE-compartmentalized genome architecture across diverse eukaryotic lineages.

In Silico Comparison of WRKY Transcription Factors in Wild and Cultivated Soybean and Their Co-expression Network Arbitrating Disease Resistance

WRKY Transcription factors (TFs) play critical roles in plant defence mechanisms that are activated in response to biotic and abiotic stresses. However, information on the Glycine soja WRKYs (GsoWRKYs) is scarce. Owing to its importance in soybean breeding, here we identified putative WRKY TFs in wild soybean, and compared the results with Glycine max WRKYs (GmaWRKYs) by phylogenetic, conserved motif, and duplication analyses. Moreover, we explored the expression trends of WRKYs in G. max (oomycete, fungi, virus, bacteria, and soybean cyst nematode) and G. soja (soybean cyst nematode), and identified commonly expressed WRKYs and their co-expressed genes. We identified, 181 and 180 putative WRKYs in G. max and G. soja, respectively. Though the number of WRKYs in both studied species is almost the same, they differ in many ways, i.e., the number of WRKYs on corresponding chromosomes, conserved domain structures, WRKYGQK motif variants, and zinc-finger motifs. WRKYs in both species grouped in three major clads, i.e., I-III, where group-II had sub-clads IIa-IIe. We found that GsoWRKYs expanded mostly through segmental duplication. A large number of WRKYs were expressed in response to biotic stresses, i.e., Phakospora pachyrhizi, Phytoplasma, Heterodera glycines, Macrophomina phaseolina, and Soybean mosaic virus; 56 GmaWRKYs were commonly expressed in soybean plants infected with these diseases. Finally, 30 and 63 GmaWRKYs and GsoWRKYs co-expressed with 205 and 123 non-WRKY genes, respectively, indicating that WRKYs play essential roles in biotic stress tolerance in Glycine species.

Distinct structural variants and repeat landscape shape the genomes of the ancient grapes Aglianico and Falanghina

Mounting evidence recognizes structural variations (SVs) and repetitive DNA sequences as crucial players in shaping the existing grape phenotypic diversity at intra- and inter-species levels. To deepen our understanding on the abundance, diversity, and distribution of SVs and repetitive DNAs, including transposable elements (TEs) and tandemly repeated satellite DNA (satDNAs), we re-sequenced the genomes of the ancient grapes Aglianico and Falanghina. The analysis of large copy number variants (CNVs) detected candidate polymorphic genes that are involved in the enological features of these varieties. In a comparative analysis of Aglianico and Falanghina sequences with 21 publicly available genomes of cultivated grapes, we provided a genome-wide annotation of grape TEs at the lineage level. We disclosed that at least two main clusters of grape cultivars could be identified based on the TEs content. Multiple TEs families appeared either significantly enriched or depleted. In addition, in silico and cytological analyses provided evidence for a diverse chromosomal distribution of several satellite repeats between Aglianico, Falanghina, and other grapes. Overall, our data further improved our understanding of the intricate grape diversity held by two Italian traditional varieties, unveiling a pool of unique candidate genes never so far exploited in breeding for improved fruit quality.

Accurate prediction of quantitative traits with failed SNP calls in canola and maize

In modern plant breeding, genomic selection is becoming the gold standard to select superior genotypes in large breeding populations that are only partially phenotyped. Many breeding programs commonly rely on single-nucleotide polymorphism (SNP) markers to capture genome-wide data for selection candidates. For this purpose, SNP arrays with moderate to high marker density represent a robust and cost-effective tool to generate reproducible, easy-to-handle, high-throughput genotype data from large-scale breeding populations. However, SNP arrays are prone to technical errors that lead to failed allele calls. To overcome this problem, failed calls are often imputed, based on the assumption that failed SNP calls are purely technical. However, this ignores the biological causes for failed calls-for example: deletions-and there is increasing evidence that gene presence-absence and other kinds of genome structural variants can play a role in phenotypic expression. Because deletions are frequently not in linkage disequilibrium with their flanking SNPs, permutation of missing SNP calls can potentially obscure valuable marker-trait associations. In this study, we analyze published datasets for canola and maize using four parametric and two machine learning models and demonstrate that failed allele calls in genomic prediction are highly predictive for important agronomic traits. We present two statistical pipelines, based on population structure and linkage disequilibrium, that enable the filtering of failed SNP calls that are likely caused by biological reasons. For the population and trait examined, prediction accuracy based on these filtered failed allele calls was competitive to standard SNP-based prediction, underlying the potential value of missing data in genomic prediction approaches. The combination of SNPs with all failed allele calls or the filtered allele calls did not outperform predictions with only SNP-based prediction due to redundancy in genomic relationship estimates.

The haplotype-resolved T2T reference genome highlights structural variation underlying agronomic traits of melon

Melon (Cucumis melo L.) is an important vegetable crop that has an extensive history of cultivation. However, the genome of wild and semi-wild melon types that can be used for the analysis of agronomic traits is not yet available. Here we report a chromosome-level T2T genome assembly for 821 (C. melo ssp. agrestis var. acidulus), a semi-wild melon with two haplotypes of ~373 Mb and ~364 Mb, respectively. Comparative genome analysis discovered a significant number of structural variants (SVs) between melo (C. melo ssp. melo) and agrestis (C. melo ssp. agrestis) genomes, including a copy number variation located in the ToLCNDV resistance locus on chromosome 11. Genome-wide association studies detected a significant signal associated with climacteric ripening and identified one candidate gene CM_ac12g14720.1 (CmABA2), encoding a cytoplasmic short chain dehydrogenase/reductase, which controls the biosynthesis of abscisic acid. This study provides valuable genetic resources for future research on melon breeding.

Structural variation of GL1 gene determines the trichome formation in Brassica juncea

A 448 kb region on chromosome B02 was delimited to be associated with trichome trait in Brassica juncea, in which the BjuVB02G54610 gene with a structural variation of 3 kb structure variation (SV) encoding a MYB transcription factor was predicted as the possible candidate gene. Mustards (Brassica juncea) are allopolyploid crops in the worldwide, and trichomes are essential quality attributes that significantly influence its taste and palpability in vegetable-use cultivars. As important accessory tissues from specialized epidermal cells, trichomes also play an important role in mitigating biotic and abiotic stresses. In this study, we constructed a F₂ segregating population using YJ27 with intensive trichome leaves and 03B0307 with glabrous leaves as parents. By bulked segregant analysis (BSA-seq), we obtained a 2.1 Mb candidate region on B02 chromosome associated with the trichome or glabrous trait formation. Then, we used 13 Kompetitive Allele Specific PCR (KASP) markers for fine mapping and finally narrowed down the candidate region to about 448 kb in length. Interestingly, among the region, there was a 3 kb sequence deletion that located on the BjuVB02G54610 gene in the F₂ individuals with trichome leaves. Genotyping results of F₂ populations confirmed this deletion (R² = 81.44%) as a major QTL. Natural population re-sequencing analysis and genotyping results further validated the key role of the 3 kb structure variation (SV) of insertion/deletion type in trichome development in B. juncea. Our findings provide important information on the formation of trichomes and potential target gene for breeding vegetable mustards.

Identification of structural variations related to drought tolerance in wheat (Triticum aestivum L.)

Structural variations are common in plant genomes, affecting meiotic recombination and distorted segregation in wheat. And presence/absence variations can significantly affect drought tolerance in wheat. Drought is a major abiotic stress limiting wheat production. Common wheat has a complex genome with three sub-genomes, which host large numbers of structural variations (SVs). SVs play critical roles in understanding the genetic contributions of plant domestication and phenotypic plasticity, but little is known about their genomic characteristics and their effects on drought tolerance. In the present study, high-resolution karyotypes of 180 doubled haploids (DHs) were developed. Signal polymorphisms between the parents involved with 8 presence-absence variations (PAVs) of tandem repeats (TR) distributed on the 7 (2A, 4A, 5A, 7A, 3B, 7B, and 2D) of 21 chromosomes. Among them, PAV on chromosome 2D showed distorted segregation, others transmit normal conforming to a 1:1 segregation ration in the population; and a PAVs recombination occurred on chromosome 2A. Association analysis of PAV and phenotypic traits under different water regimes, we found PAVs on chromosomes 4A, 5A, and 7B showed negative effect on grain length (GL) and grain width (GW); PAV.7A had opposite effect on grain thickness (GT) and spike length (SL), with the effect on traits differing under different water regimes. PAVs on linkage group 2A, 4A, 7A, 2D, and 7B associated with the drought tolerance coefficients (DTCs), and significant negative effect on drought resistance values (D values) were detected in PAV.7B. Additionally, quantitative trait loci (QTL) associated with phenotypic traits using the 90 K SNP array showed QTL for DTCs and grain-related traits in chromosomes 4A, and 5A, 3B were co-localized in differential regions of PAVs. These PAVs can cause the differentiation of the target region of SNP and could be used for genetic improvement of agronomic traits under drought stress via marker-assisted selection (MAS) breeding.

An NBS-LRR protein in the Rpp1 locus negates the dominance of Rpp1-mediated resistance against Phakopsora pachyrhizi in soybean

The soybean Rpp1 locus confers resistance to Phakopsora pachyrhizi, causal agent of rust, and resistance is usually dominant over susceptibility. However, dominance of Rpp1-mediated resistance is lost when a resistant genotype (Rpp1 or Rpp1b) is crossed with susceptible line TMG06_0011, and the mechanism of this dominant susceptibility (DS) is unknown. Sequencing the Rpp1 region reveals that the TMG06_0011 Rpp1 locus has a single nucleotide-binding site leucine-rich repeat (NBS-LRR) gene (DS-R), whereas resistant PI 594760B (Rpp1b) is similar to PI 200492 (Rpp1) and has three NBS-LRR resistance gene candidates. Evidence that DS-R is the cause of DS was reflected in virus-induced gene silencing of DS-R in Rpp1b/DS-R or Rpp1/DS-R heterozygous plants with resistance partially restored. In heterozygous Rpp1b/DS-R plants, expression of Rpp1b candidate genes was not significantly altered, indicating no effect of DS-R on transcription. Physical interaction of the DS-R protein with candidate Rpp1b resistance proteins was supported by yeast two-hybrid studies and in silico modeling. Thus, we conclude that suppression of resistance most likely does not occur at the transcript level, but instead probably at the protein level, possibly with Rpp1 function inhibited by binding to the DS-R protein. The DS-R gene was found in other soybean lines, with an estimated allele frequency of 6% in a diverse population, and also found in wild soybean (Glycine soja). The identification of a dominant susceptible NBS-LRR gene provides insight into the behavior of NBS-LRR proteins and serves as a reminder to breeders that the dominance of an R gene can be influenced by a susceptibility allele.

The soybean ubiquitin-proteasome system: Current knowledge and future perspective

Increasing soybean [Glycine max (L.) Merr.] yield has become a worldwide scientific problem in the world. Many studies have shown that ubiquitination plays a key role in stress response and yield formation. In the UniProtKB database, 2,429 ubiquitin-related proteins were predicted in soybean, however, <20 were studied. One key way to address this lack of progress in increasing soybean yield will be a deeper understanding of the ubiquitin-proteasome system (UPS) in soybean. In this review, we summarized the current knowledge about soybean ubiquitin-related proteins and discussed the method of combining phenotype, mutant library, transgenic system, genomics, and proteomics approaches to facilitate the exploration of the soybean UPS. We also proposed the strategy of applying the UPS in soybean improvement based on related studies in model plants. Our review will be helpful for soybean scientists to learn current research progress of the soybean UPS and further lay a theoretical reference for the molecular improvement of soybean in future research by use of this knowledge.

N. Lopes I, Abdelnoor RV, Marcelino-Guimarães FC. Mapping of a soybean rust resistance in PI 594756 at the Rpp1 locus

Asian soybean rust (ASR), caused by the fungus Phakopsora pachyrhizi, is the main disease affecting soybean in Brazil. This study aimed at investigating and mapping the resistance of the PI 594756 to P. pachyrhizi, by using Bulked Segregant Analysis (BSA). The PI 594756 and the susceptible PI 594891 were crossed and the resulting F2 and F2:3 populations (208 and 1770 plants, respectively) were tested against ASR. Also, these PIs and differential varieties were tested against a panel of monosporic isolates. Plants presenting tan lesions were classified as susceptible (S) while plants presenting reddish-brown (RB) lesions were classified as resistant. DNA bulks were genotyped with Infinium BeadChips and the genomic region identified was further analyzed in the F2 individuals with target GBS (tGBS). PI 594,56 presented a unique resistance profile compared to the differential varieties. The resistance was monogenic dominant; however, it was classified as incompletely dominant when quantitatively studied. Genetic and QTL mapping placed the PI 594756 gene between the genomic region located at 55,863,741 and 56,123,516 bp of chromosome 18. This position is slightly upstream mapping positions of Rpp1 (PI 200492) and Rpp1-b (PI 594538A). Finally, we performed a haplotype analysis in a whole genomic sequencing-SNP database composed of Brazilian historical germplasm and sources of Rpp genes. We found SNPs that successfully differentiated the new PI 594756 allele from Rpp1 and Rpp1-b sources. The haplotype identified can be used as a tool for marker-assisted selection (MAS).

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01358-4.

A high-quality pseudo-phased genome for Melaleuca quinquenervia shows allelic diversity of NLR-type resistance genes

BACKGROUND

Melaleuca quinquenervia (broad-leaved paperbark) is a coastal wetland tree species that serves as a foundation species in eastern Australia, Indonesia, Papua New Guinea, and New Caledonia. While extensively cultivated for its ornamental value, it has also become invasive in regions like Florida, USA. Long-lived trees face diverse pest and pathogen pressures, and plant stress responses rely on immune receptors encoded by the nucleotide-binding leucine-rich repeat (NLR) gene family. However, the comprehensive annotation of NLR encoding genes has been challenging due to their clustering arrangement on chromosomes and highly repetitive domain structure; expansion of the NLR gene family is driven largely by tandem duplication. Additionally, the allelic diversity of the NLR gene family remains largely unexplored in outcrossing tree species, as many genomes are presented in their haploid, collapsed state.

RESULTS

We assembled a chromosome-level pseudo-phased genome for M. quinquenervia and described the allelic diversity of plant NLRs using the novel FindPlantNLRs pipeline. Analysis reveals variation in the number of NLR genes on each haplotype, distinct clustering patterns, and differences in the types and numbers of novel integrated domains.

CONCLUSIONS

The high-quality M. quinquenervia genome assembly establishes a new framework for functional and evolutionary studies of this significant tree species. Our findings suggest that maintaining allelic diversity within the NLR gene family is crucial for enabling responses to environmental stress, particularly in long-lived plants.

A mutant allele of the flowering promoting factor 1 gene at the tomato BRACHYTIC locus reduces plant height with high quality fruit

Reduced plant height due to shortened stems is beneficial for improving crop yield potential, better resilience to biotic/abiotic stresses, and rapid crop producer adoption of the agronomic and management practices. Breeding tomato plants with a reduced height, however, poses a particular challenge because this trait is often associated with a significant fruit size (weight) reduction. The tomato BRACHYTIC (BR) locus controls plant height. Genetic mapping and genome assembly revealed three flowering promoting factor 1 (FPF1) genes located within the BR mapping interval, and a complete coding sequence deletion of the telomere proximal FPF1 (Solyc01g066980) was found in the br allele but not in BR. The knock-out of Solyc01g066980 in BR large-fruited fresh-market tomato reduced the height and fruit yield, but the ability to produce large size fruits was retained. However, concurrent yield evaluation of a pair of sister lines with or without the br allele revealed that artificial selection contributes to commercially acceptable yield potential in br tomatoes. A network analysis of gene-expression patterns across genotypes, tissues, and the gibberellic acid (GA) treatment revealed that member(s) of the FPF1 family may play a role in the suppression of the GA biosynthesis in roots and provided a framework for identifying the responsible molecular signaling pathways in br-mediated phenotypic changes. Lastly, mutations of br homologs also resulted in reduced height. These results shed light on the genetic and physiological mechanisms by which the br allele alters tomato architecture.

Identification of Candidate Genes for a Major Quantitative Disease Resistance Locus From Soybean PI 427105B for Resistance to Phytophthora sojae

Phytophthora root and stem rot is a yield-limiting soybean disease caused by the soil-borne oomycete Phytophthora sojae. Although multiple quantitative disease resistance loci (QDRL) have been identified, most explain <10% of the phenotypic variation (PV). The major QDRL explaining up to 45% of the PV were previously identified on chromosome 18 and represent a valuable source of resistance for soybean breeding programs. Resistance alleles from plant introductions 427105B and 427106 significantly increase yield in disease-prone fields and result in no significant yield difference in fields with less to no disease pressure. In this study, high-resolution mapping reduced the QDRL interval to 3.1 cm, and RNA-seq analysis of near-isogenic lines (NILs) varying at QDRL-18 pinpointed a single gene of interest which was downregulated in inoculated NILs carrying the resistant allele compared to inoculated NILs with the susceptible allele. This gene of interest putatively encodes a serine-threonine kinase (STK) related to the AtCR4 family and may be acting as a susceptibility factor, based on the specific increase of jasmonic acid concentration in inoculated NILs. This work facilitates further functional analyses and marker-assisted breeding efforts by prioritizing candidate genes and narrowing the targeted region for introgression.

Oil plant genomes: current state of the science

Vegetable oils are an indispensable nutritional component of the human diet as well as important raw materials for a variety of industrial applications such as pharmaceuticals, cosmetics, oleochemicals, and biofuels. Oil plant genomes are highly diverse, and their genetic variation leads to a diversity in oil biosynthesis and accumulation along with agronomic traits. This review discusses plant oil biosynthetic pathways, current state of genome assembly, polyploidy and asymmetric evolution of genomes of oil plants and their wild relatives, and research progress of pan-genomics in oil plants. The availability of complete high-resolution genomes and pan-genomes has enabled the identification of structural variations in the genomes that are associated with the diversity of agronomic and environment fitness traits. These and future genomes also provide powerful tools to understand crop evolution and to harvest the rich natural variations to improve oil crops for enhanced productivity, oil quality, and adaptability to changing environments.

Genomic structural variation in tomato and its role in plant immunity

It is well known that large genomic variations can greatly impact the phenotype of an organism. Structural Variants (SVs) encompass any genomic variation larger than 30 base pairs, and include changes caused by deletions, inversions, duplications, transversions, and other genome modifications. Due to their size and complex nature, until recently, it has been difficult to truly capture these variations. Recent advances in sequencing technology and computational analyses now permit more extensive studies of SVs in plant genomes. In tomato, advances in sequencing technology have allowed researchers to sequence hundreds of genomes from tomatoes, and tomato relatives. These studies have identified SVs related to fruit size and flavor, as well as plant disease response, resistance/susceptibility, and the ability of plants to detect pathogens (immunity). In this review, we discuss the implications for genomic structural variation in plants with a focus on its role in tomato immunity. We also discuss how advances in sequencing technology have led to new discoveries of SVs in more complex genomes, the current evidence for the role of SVs in biotic and abiotic stress responses, and the outlook for genetic modification of SVs to advance plant breeding objectives.

Combined use of Oxford Nanopore and Illumina sequencing yields insights into soybean structural variation biology

BACKGROUND

Structural variants (SVs), including deletions, insertions, duplications, and inversions, are relatively long genomic variations implicated in a diverse range of processes from human disease to ecology and evolution. Given their complex signatures, tendency to occur in repeated regions, and large size, discovering SVs based on short reads is challenging compared to single-nucleotide variants. The increasing availability of long-read technologies has greatly facilitated SV discovery; however, these technologies remain too costly to apply routinely to population-level studies. Here, we combined short-read and long-read sequencing technologies to provide a comprehensive population-scale assessment of structural variation in a panel of Canadian soybean cultivars.

RESULTS

We used Oxford Nanopore long-read sequencing data (~12× mean coverage) for 17 samples to both benchmark SV calls made from Illumina short-read data and predict SVs that were subsequently genotyped in a population of 102 samples using Illumina data. Benchmarking results show that variants discovered using Oxford Nanopore can be accurately genotyped from the Illumina data. We first use the genotyped deletions and insertions for population genetics analyses and show that results are comparable to those based on single-nucleotide variants. We observe that the population frequency and distribution within the genome of deletions and insertions are constrained by the location of genes. Gene Ontology and PFAM domain enrichment analyses also confirm previous reports that genes harboring high-frequency deletions and insertions are enriched for functions in defense response. Finally, we discover polymorphic transposable elements from the deletions and insertions and report evidence of the recent activity of a Stowaway MITE.

CONCLUSIONS

We show that structural variants discovered using Oxford Nanopore data can be genotyped with high accuracy from Illumina data. Our results demonstrate that long-read and short-read sequencing technologies can be efficiently combined to enhance SV analysis in large populations, providing a reusable framework for their study in a wider range of samples and non-model species.

Progress in soybean functional genomics over the past decade

Soybean is one of the most important oilseed and fodder crops. Benefiting from the efforts of soybean breeders and the development of breeding technology, large number of germplasm has been generated over the last 100 years. Nevertheless, soybean breeding needs to be accelerated to meet the needs of a growing world population, to promote sustainable agriculture and to address future environmental changes. The acceleration is highly reliant on the discoveries in gene functional studies. The release of the reference soybean genome in 2010 has significantly facilitated the advance in soybean functional genomics. Here, we review the research progress in soybean omics (genomics, transcriptomics, epigenomics and proteomics), germplasm development (germplasm resources and databases), gene discovery (genes that are responsible for important soybean traits including yield, flowering and maturity, seed quality, stress resistance, nodulation and domestication) and transformation technology during the past decade. At the end, we also briefly discuss current challenges and future directions.

Entailing the Next-Generation Sequencing and Metabolome for Sustainable Agriculture by Improving Plant Tolerance

Crop production is a serious challenge to provide food for the 10 billion individuals forecasted to live across the globe in 2050. The scientists' emphasize establishing an equilibrium among diversity and quality of crops by enhancing yield to fulfill the increasing demand for food supply sustainably. The exploitation of genetic resources using genomics and metabolomics strategies can help generate resilient plants against stressors in the future. The innovation of the next-generation sequencing (NGS) strategies laid the foundation to unveil various plants' genetic potential and help us to understand the domestication process to unmask the genetic potential among wild-type plants to utilize for crop improvement. Nowadays, NGS is generating massive genomic resources using wild-type and domesticated plants grown under normal and harsh environments to explore the stress regulatory factors and determine the key metabolites. Improved food nutritional value is also the key to eradicating malnutrition problems around the globe, which could be attained by employing the knowledge gained through NGS and metabolomics to achieve suitability in crop yield. Advanced technologies can further enhance our understanding in defining the strategy to obtain a specific phenotype of a crop. Integration among bioinformatic tools and molecular techniques, such as marker-assisted, QTLs mapping, creation of reference genome, de novo genome assembly, pan- and/or super-pan-genomes, etc., will boost breeding programs. The current article provides sequential progress in NGS technologies, a broad application of NGS, enhancement of genetic manipulation resources, and understanding the crop response to stress by producing plant metabolites. The NGS and metabolomics utilization in generating stress-tolerant plants/crops without deteriorating a natural ecosystem is considered a sustainable way to improve agriculture production. This highlighted knowledge also provides useful research that explores the suitable resources for agriculture sustainability.

In-Depth Sequence Analysis of Bread Wheat VRN1 Genes

The VERNALIZATION1 (VRN1) gene encodes a MADS-box transcription factor and plays an important role in the cold-induced transition from the vegetative to reproductive stage. Allelic variability of VRN1 homoeologs has been associated with large differences in flowering time. The aim of this study was to investigate the genetic variability of VRN1 homoeologs (VRN-A1, VRN-B1 and VRN-D1). We performed an in-depth sequence analysis of VRN1 homoeologs in a panel of 105 winter and spring varieties of hexaploid wheat. We describe the novel allele Vrn-B1f with an 836 bp insertion within intron 1 and show its specific expression pattern associated with reduced heading time. We further provide the complete sequence of the Vrn-A1b allele, revealing a 177 bp insertion in intron 1, which is transcribed into an alternative splice variant. Copy number variation (CNV) analysis of VRN1 homoeologs showed that VRN-B1 and VRN-D1 are present in only one copy. The copy number of recessive vrn-A1 ranged from one to four, while that of dominant Vrn-A1 was one or two. Different numbers of Vrn-A1a copies in the spring cultivars Branisovicka IX/49 and Bastion did not significantly affect heading time. We also report on the deletion of secondary structures (G-quadruplex) in promoter sequences of cultivars with more vrn-A1 copies.

High quality haplotype-resolved genome assemblies of Populus tomentosa Carr., a stabilized interspecific hybrid species widespread in Asia

Populus has a wide ecogeographical range spanning the Northern Hemisphere, and interspecific hybrids are common. Populus tomentosa Carr. is widely distributed and cultivated in the eastern region of Asia, where it plays multiple important roles in forestry, agriculture, conservation, and urban horticulture. Reference genomes are available for several Populus species, however, our goals were to produce a very high quality de novo chromosome-level genome assembly in P. tomentosa genome that could serve as a reference for evolutionary and ecological studies of hybrid speciation throughout the genus. Here, combining long-read sequencing and Hi-C scaffolding, we present a high-quality, haplotype-resolved genome assembly. The genome size was 740.2 Mb, with a contig N50 size of 5.47 Mb and a scaffold N50 size of 46.68 Mb, consisting of 38 chromosomes, as expected with the known diploid chromosome number (2n = 2x = 38). A total of 59,124 protein-coding genes were identified. Phylogenomic analyses revealed that P. tomentosa is comprised of two distinct subgenomes, which we deomonstrate is likely to have resulted from hybridization between Populus adenopoda as the female parent and Populus alba var. pyramidalis as the male parent, with an origin of approximately 3.93 Ma. Although highly colinear, significant structural variation was found between the two subgenomes. Our study provides a valuable resource for ecological genetics and forest biotechnology.

Expanding the conservation genomics toolbox: Incorporating structural variants to enhance genomic studies for species of conservation concern

Structural variants (SVs) are large rearrangements (>50 bp) within the genome that impact gene function and the content and structure of chromosomes. As a result, SVs are a significant source of functional genomic variation, that is, variation at genomic regions underpinning phenotype differences, that can have large effects on individual and population fitness. While there are increasing opportunities to investigate functional genomic variation in threatened species via single nucleotide polymorphism (SNP) data sets, SVs remain understudied despite their potential influence on fitness traits of conservation interest. In this future-focused Opinion, we contend that characterizing SVs offers the conservation genomics community an exciting opportunity to complement SNP-based approaches to enhance species recovery. We also leverage the existing literature-predominantly in human health, agriculture and ecoevolutionary biology-to identify approaches for readily characterizing SVs and consider how integrating these into the conservation genomics toolbox may transform the way we manage some of the world's most threatened species.

Next-Generation Breeding Strategies for Climate-Ready Crops

Climate change is a threat to global food security due to the reduction of crop productivity around the globe. Food security is a matter of concern for stakeholders and policymakers as the global population is predicted to bypass 10 billion in the coming years. Crop improvement via modern breeding techniques along with efficient agronomic practices innovations in microbiome applications, and exploiting the natural variations in underutilized crops is an excellent way forward to fulfill future food requirements. In this review, we describe the next-generation breeding tools that can be used to increase crop production by developing climate-resilient superior genotypes to cope with the future challenges of global food security. Recent innovations in genomic-assisted breeding (GAB) strategies allow the construction of highly annotated crop pan-genomes to give a snapshot of the full landscape of genetic diversity (GD) and recapture the lost gene repertoire of a species. Pan-genomes provide new platforms to exploit these unique genes or genetic variation for optimizing breeding programs. The advent of next-generation clustered regularly interspaced short palindromic repeat/CRISPR-associated (CRISPR/Cas) systems, such as prime editing, base editing, and de nova domestication, has institutionalized the idea that genome editing is revamped for crop improvement. Also, the availability of versatile Cas orthologs, including Cas9, Cas12, Cas13, and Cas14, improved the editing efficiency. Now, the CRISPR/Cas systems have numerous applications in crop research and successfully edit the major crop to develop resistance against abiotic and biotic stress. By adopting high-throughput phenotyping approaches and big data analytics tools like artificial intelligence (AI) and machine learning (ML), agriculture is heading toward automation or digitalization. The integration of speed breeding with genomic and phenomic tools can allow rapid gene identifications and ultimately accelerate crop improvement programs. In addition, the integration of next-generation multidisciplinary breeding platforms can open exciting avenues to develop climate-ready crops toward global food security.

Dissecting nematode resistance regions in soybean revealed pleiotropic effect of soybean cyst and reniform nematode resistance genes

Reniform nematode (RN, Rotylenchulus reniformis Linford & Oliveira) has emerged as one of the most important plant parasitic nematodes of soybean [Glycine max (L.) Merr.]. Planting resistant varieties is the most effective strategy for nematode management. The objective of this study was to identify quantitative trait loci (QTL) for RN resistance in an exotic soybean line, PI 438489B, using two linkage maps constructed from the Universal Soybean Linkage Panel (USLP 1.0) and next-generation whole-genome resequencing (WGRS) technology. Two QTL controlling RN resistance were identified-the soybean cyst nematode (SCN, Heterodera glycines) resistance gene GmSNAP18 at the rhg1 locus and its paralog GmSNAP11. Strong association between resistant phenotype and haplotypes of the GmSNAP11 and GmSNAP18 was observed. The results indicated that GmSNAP11 possibly could have epistatic effect on GmSNAP18, or vice versa, with the presence of a significant correlation in RN resistance of rhg1-a GmSNAP18 vs. rhg1-b GmSNAP18. Most importantly, our preliminary data suggested that GmSNAP18 and GmSNAP11 proteins physically interact in planta, suggesting that they belong to the same pathway for resistance. Unlike GmSNAP18, no indication of GmSNAP11 copy number variation was found. Moreover, gene-based single nucleotide polymorphism (SNP) markers were developed for rapid detection of RN or SCN resistance at these loci. Our analysis substantiates synergic interaction between GmSNAP11 and GmSNAP18 genes and confirms their roles in RN as well as SCN resistance. These results could contribute to a better understanding of evolution and subfunctionalization of genes conferring resistance to multiple nematode species and provide a framework for further investigations.

Introgressing the Aegilops tauschii genome into wheat as a basis for cereal improvement

Increasing crop production is necessary to feed the world's expanding population, and crop breeders often utilize genetic variations to improve crop yield and quality. However, the narrow diversity of the wheat D genome seriously restricts its selective breeding. A practical solution is to exploit the genomic variations of Aegilops tauschii via introgression. Here, we established a rapid introgression platform for transferring the overall genetic variations of A. tauschii to elite wheats, thereby enriching the wheat germplasm pool. To accelerate the process, we assembled four new reference genomes, resequenced 278 accessions of A. tauschii and constructed the variation landscape of this wheat progenitor species. Genome comparisons highlighted diverse functional genes or novel haplotypes with potential applications in wheat improvement. We constructed the core germplasm of A. tauschii, including 85 accessions covering more than 99% of the species' overall genetic variations. This was crossed with elite wheat cultivars to generate an A. tauschii-wheat synthetic octoploid wheat (A-WSOW) pool. Laboratory and field analysis with two examples of the introgression lines confirmed its great potential for wheat breeding. Our high-quality reference genomes, genomic variation landscape of A. tauschii and the A-WSOW pool provide valuable resources to facilitate gene discovery and breeding in wheat.

A novel deletion in FLOWERING LOCUS T modulates flowering time in winter oilseed rape

A novel structural variant was discovered in the FLOWERING LOCUS T orthologue BnaFT.A02 by long-read sequencing. Nested association mapping in an elite winter oilseed rape population revealed that this 288 bp deletion associates with early flowering, putatively by modification of binding-sites for important flowering regulation genes. Perfect timing of flowering is crucial for optimal pollination and high seed yield. Extensive previous studies of flowering behavior in Brassica napus (canola, rapeseed) identified mutations in key flowering regulators which differentiate winter, semi-winter and spring ecotypes. However, because these are generally fixed in locally adapted genotypes, they have only limited relevance for fine adjustment of flowering time in elite cultivar gene pools. In crosses between ecotypes, the ecotype-specific major-effect mutations mask minor-effect loci of interest for breeding. Here, we investigated flowering time in a multiparental mapping population derived from seven elite winter oilseed rape cultivars which are fixed for major-effect mutations separating winter-type rapeseed from other ecotypes. Association mapping revealed eight genomic regions on chromosomes A02, C02 and C03 associating with fine modulation of flowering time. Long-read genomic resequencing of the seven parental lines identified seven structural variants coinciding with candidate genes for flowering time within chromosome regions associated with flowering time. Segregation patterns for these variants in the elite multiparental population and a diversity set of winter types using locus-specific assays revealed significant associations with flowering time for three deletions on chromosome A02. One of these was a previously undescribed 288 bp deletion within the second intron of FLOWERING LOCUS T on chromosome A02, emphasizing the advantage of long-read sequencing for detection of structural variants in this size range. Detailed analysis revealed the impact of this specific deletion on flowering-time modulation under extreme environments and varying day lengths in elite, winter-type oilseed rape.

New insights into homoeologous copy number variations in the hexaploid wheat genome

Bread wheat is an allohexaploid species originating from two successive and recent rounds of hybridization between three diploid species that were very similar in terms of chromosome number, genome size, TE content, gene content and synteny. As a result, it has long been considered that most of the genes were in three pairs of homoeologous copies. However, these so-called triads represent only one half of wheat genes, while the remaining half belong to homoeologous groups with various number of copies across subgenomes. In this study, we examined and compared the distribution, conservation, function, expression and epigenetic profiles of triads with homoeologous groups having undergone a deletion (dyads) or a duplication (tetrads) in one subgenome. We show that dyads and tetrads are mostly located in distal regions and have lower expression level and breadth than triads. Moreover, they are enriched in functions related to adaptation and more associated with the repressive H3K27me3 modification. Altogether, these results suggest that triads mainly correspond to housekeeping genes and are part of the core genome, while dyads and tetrads belong to the Triticeae dispensable genome. In addition, by comparing the different categories of dyads and tetrads, we hypothesize that, unlike most of the allopolyploid species, subgenome dominance and biased fractionation are absent in hexaploid wheat. Differences observed between the three subgenomes are more likely related to two successive and ongoing waves of post-polyploid diploidization, that had impacted A and B more significantly than D, as a result of the evolutionary history of hexaploid wheat.

The Genetic Basis of Tomato Aroma

Tomato (Solanum lycopersicum L.) aroma is determined by the interaction of volatile compounds (VOCs) released by the tomato fruits with receptors in the nose, leading to a sensorial impression, such as "sweet", "smoky", or "fruity" aroma. Of the more than 400 VOCs released by tomato fruits, 21 have been reported as main contributors to the perceived tomato aroma. These VOCs can be grouped in five clusters, according to their biosynthetic origins. In the last decades, a vast array of scientific studies has investigated the genetic component of tomato aroma in modern tomato cultivars and their relatives. In this paper we aim to collect, compare, integrate and summarize the available literature on flavour-related QTLs in tomato. Three hundred and 5ifty nine (359) QTLs associated with tomato fruit VOCs were physically mapped on the genome and investigated for the presence of potential candidate genes. This review makes it possible to (i) pinpoint potential donors described in literature for specific traits, (ii) highlight important QTL regions by combining information from different populations, and (iii) pinpoint potential candidate genes. This overview aims to be a valuable resource for researchers aiming to elucidate the genetics underlying tomato flavour and for breeders who aim to improve tomato aroma.

High-Throughput Genotyping Technologies in Plant Taxonomy

Molecular markers provide researchers with a powerful tool for variation analysis between plant genomes. They are heritable and widely distributed across the genome and for this reason have many applications in plant taxonomy and genotyping. Over the last decade, molecular marker technology has developed rapidly and is now a crucial component for genetic linkage analysis, trait mapping, diversity analysis, and association studies. This chapter focuses on molecular marker discovery, its application, and future perspectives for plant genotyping through pangenome assemblies. Included are descriptions of automated methods for genome and sequence distance estimation, genome contaminant analysis in sequence reads, genome structural variation, and SNP discovery methods.

Prospects of next generation sequencing in lentil breeding

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

Identifying barley pan-genome sequence anchors using genetic mapping and machine learning

We identified 1.844 million barley pan-genome sequence anchors from 12,306 genotypes using genetic mapping and machine learning. There is increasing evidence that genes from a given crop genotype are far to cover all genes in that species; thus, building more comprehensive pan-genomes is of great importance in genetic research and breeding. Obtaining a thousand-genotype scale pan-genome using deep-sequencing data is currently impractical for species like barley which has a huge and highly repetitive genome. To this end, we attempted to identify barley pan-genome sequence anchors from a large quantity of genotype-by-sequencing (GBS) datasets by combining genetic mapping and machine learning algorithms. Based on the GBS sequences from 11,166 domesticated and 1140 wild barley genotypes, we identified 1.844 million pan-genome sequence anchors. Of them, 532,253 were identified as presence/absence variation (PAV) tags. Through aligning these PAV tags to the genome of hulless barley genotype Zangqing320, our analysis resulted in a validation of 83.6% of them from the domesticated genotypes and 88.6% from the wild barley genotypes. Association analyses against flowering time, plant height and kernel size showed that the relative importance of the PAV and non-PAV tags varied for different traits. The pan-genome sequence anchors based on GBS tags can facilitate the construction of a comprehensive pan-genome and greatly assist various genetic studies including identification of structural variation, genetic mapping and breeding in barley.

A genome-wide survey of copy number variations reveals an asymmetric evolution of duplicated genes in rice

Background

Copy number variations (CNVs) are an important type of structural variations in the genome that usually affect gene expression levels by gene dosage effect. Understanding CNVs as part of genome evolution may provide insights into the genetic basis of important agricultural traits and contribute to the crop breeding in the future. While available methods to detect CNVs utilizing next-generation sequencing technology have helped shed light on prevalence and effects of CNVs, the complexity of crop genomes poses a major challenge and requires development of additional tools.

Results

Here, we generated genomic and transcriptomic data of 93 rice (Oryza sativa L.) accessions and developed a comprehensive pipeline to call CNVs in this large-scale dataset. We analyzed the correlation between CNVs and gene expression levels and found that approximately 13% of the identified genes showed a significant correlation between their expression levels and copy numbers. Further analysis showed that about 36% of duplicate pairs were involved in pseudogenetic events while only 5% of them showed functional differentiation. Moreover, the offspring copy mainly contributed to the expression levels and seemed more likely to become a pseudogene, whereas the parent copy tended to maintain the function of ancestral gene.

Conclusion

We provide a high-accuracy CNV dataset that will contribute to functional genomics studies and molecular breeding in rice. We also showed that gene dosage effect of CNVs in rice is not exponential or linear. Our work demonstrates that the evolution of duplicated genes is asymmetric in both expression levels and gene fates, shedding a new insight into the evolution of duplicated genes.

AthCNV: A Map of DNA Copy Number Variations in the Arabidopsis Genome

Copy number variations (CNVs) greatly contribute to intraspecies genetic polymorphism and phenotypic diversity. Recent analyses of sequencing data for >1000 Arabidopsis (Arabidopsis thaliana) accessions focused on small variations and did not include CNVs. Here, we performed genome-wide analysis and identified large indels (50 to 499 bp) and CNVs (500 bp and larger) in these accessions. The CNVs fully overlap with 18.3% of protein-coding genes, with enrichment for evolutionarily young genes and genes involved in stress and defense. By combining analysis of both genes and transposable elements (TEs) affected by CNVs, we revealed that the variation statuses of genes and TEs are tightly linked and jointly contribute to the unequal distribution of these elements in the genome. We also determined the gene copy numbers in a set of 1060 accessions and experimentally validated the accuracy of our predictions by multiplex ligation-dependent probe amplification assays. We then successfully used the CNVs as markers to analyze population structure and migration patterns. Finally, we examined the impact of gene dosage variation triggered by a CNV spanning the SEC10 gene on SEC10 expression at both the transcript and protein levels. The catalog of CNVs, CNV-overlapping genes, and their genotypes in a top model dicot will stimulate the exploration of the genetic basis of phenotypic variation.

Early Drought-Responsive Genes Are Variable and Relevant to Drought Tolerance

Drought stress is an important crop yield limiting factor worldwide. Plant physiological responses to drought stress are driven by changes in gene expression. While drought-responsive genes (DRGs) have been identified in maize, regulation patterns of gene expression during progressive water deficits remain to be elucidated. In this study, we generated time-series transcriptomic data from the maize inbred line B73 under well-watered and drought conditions. Comparisons between the two conditions identified 8,626 DRGs and the stages (early, middle, and late drought) at which DRGs occurred. Different functional groups of genes were regulated at the three stages. Specifically, early and middle DRGs display higher copy number variation among diverse Zea mays lines, and they exhibited stronger associations with drought tolerance as compared to late DRGs. In addition, correlation of expression between small RNAs (sRNAs) and DRGs from the same samples identified 201 negatively sRNA/DRG correlated pairs, including genes showing high levels of association with drought tolerance, such as two glutamine synthetase genes, gln2 and gln6 The characterization of dynamic gene responses to progressive drought stresses indicates important adaptive roles of early and middle DRGs, as well as roles played by sRNAs in gene expression regulation upon drought stress.

SMRT sequencing of the Oryza rufipogon genome reveals the genomic basis of rice adaptation

Asian cultivated rice is believed to have been domesticated from a wild progenitor, Oryza rufipogon, offering promising sources of alleles for world rice improvement. Here we first present a high-quality chromosome-scale genome of the typical O. rufipogon. Comparative genomic analyses of O. sativa and its two wild progenitors, O. nivara and O. rufipogon, identified many dispensable genes functionally enriched in the reproductive process. We detected millions of genomic variants, of which large-effect mutations could affect agronomically relevant traits. We demonstrate how lineage-specific expansion of gene families may have contributed to the formation of reproduction isolation. We document thousands of genes with signatures of positive selection that are mainly involved in the reproduction and response to biotic- and abiotic stresses. We show that selection pressures may serve as forces to govern substantial genomic alterations that form the genetic basis of rapid evolution of mating and reproductive systems under diverse habitats.

Li et al. provide a chromosome-based genome of the ancestral rice relative Oryza rufipogon using SMRT, 10x and Hi-C. Using comparative genomics analysis, they identify genes important for agronomical traits and provide insights for the evolution of reproductive systems which can facilitate rice breeding programs.

Structural genome analysis in cultivated potato taxa

KEY MESSAGE

Twelve potato accessions were selected to represent two principal views on potato taxonomy. The genomes were sequenced and analyzed for structural variation (copy number variation) against three published potato genomes. The common potato (Solanum tuberosum L.) is an important staple crop with a highly heterozygous and complex tetraploid genome. The other taxa of cultivated potato contain varying ploidy levels (2X-5X), and structural variations are common in the genomes of these species, likely contributing to the diversification or agronomic traits during domestication. Increased understanding of the genomes and genomic variation will aid in the exploration of novel agronomic traits. Thus, sequencing data from twelve potato landraces, representing the four ploidy levels, were used to identify structural genomic variation compared to the two currently available reference genomes, a double monoploid potato genome and a diploid inbred clone of S. chacoense. The results of a copy number variation analysis showed that in the majority of the genomes, while the number of deletions is greater than the number of duplications, the number of duplicated genes is greater than the number of deleted ones. Specific regions in the twelve potato genomes have a high density of CNV events. Further, the auxin-induced SAUR genes (involved in abiotic stress), disease resistance genes and the 2-oxoglutarate/Fe(II)-dependent oxygenase superfamily proteins, among others, had increased copy numbers in these sequenced genomes relative to the references.

Functional analysis and development of a CRISPR/Cas9 allelic series for a CPR5 ortholog necessary for proper growth of soybean trichomes

Developments in genomic and genome editing technologies have facilitated the mapping, cloning, and validation of genetic variants underlying trait variation. This study combined bulked-segregant analysis, array comparative genomic hybridization, and CRISPR/Cas9 methodologies to identify a CPR5 ortholog essential for proper trichome growth in soybean (Glycine max). A fast neutron mutant line exhibited short trichomes with smaller trichome nuclei compared to its parent line. A fast neutron-induced deletion was identified within an interval on chromosome 6 that co-segregated with the trichome phenotype. The deletion encompassed six gene models including an ortholog of Arabidopsis thaliana CPR5. CRISPR/Cas9 was used to mutate the CPR5 ortholog, resulting in five plants harboring a total of four different putative knockout alleles and two in-frame alleles. Phenotypic analysis of the mutants validated the candidate gene, and included intermediate phenotypes that co-segregated with the in-frame alleles. These findings demonstrate that the CPR5 ortholog is essential for proper growth and development of soybean trichomes, similar to observations in A. thaliana. Furthermore, this work demonstrates the value of using CRISPR/Cas9 to generate an allelic series and intermediate phenotypes for functional analysis of candidate genes and/or the development of novel traits.

The rhg1-a (Rhg1 low-copy) nematode resistance source harbors a copia-family retrotransposon within the Rhg1-encoded α-SNAP gene

Soybean growers widely use the Resistance to Heterodera glycines 1 (Rhg1) locus to reduce yield losses caused by soybean cyst nematode (SCN). Rhg1 is a tandemly repeated four gene block. Two classes of SCN resistance-conferring Rhg1 haplotypes are recognized: rhg1-a ("Peking-type," low-copy number, three or fewer Rhg1 repeats) and rhg1-b ("PI 88788-type," high-copy number, four or more Rhg1 repeats). The rhg1-a and rhg1-b haplotypes encode α-SNAP (alpha-Soluble NSF Attachment Protein) variants α-SNAP Rhg1 LC and α-SNAP Rhg1 HC, respectively, with differing atypical C-terminal domains, that contribute to SCN resistance. Here we report that rhg1-a soybean accessions harbor a copia retrotransposon within their Rhg1 Glyma.18G022500 (α-SNAP-encoding) gene. We termed this retrotransposon "RAC," for Rhg1 alpha-SNAP copia. Soybean carries multiple RAC-like retrotransposon sequences. The Rhg1 RAC insertion is in the Glyma.18G022500 genes of all true rhg1-a haplotypes we tested and was not detected in any examined rhg1-b or Rhg1WT (single-copy) soybeans. RAC is an intact element residing within intron 1, anti-sense to the rhg1-a α-SNAP open reading frame. RAC has intrinsic promoter activities, but overt impacts of RAC on transgenic α-SNAP Rhg1 LC mRNA and protein abundance were not detected. From the native rhg1-a RAC+ genomic context, elevated α-SNAP Rhg1 LC protein abundance was observed in syncytium cells, as was previously observed for α-SNAP Rhg1 HC (whose rhg1-b does not carry RAC). Using a SoySNP50K SNP corresponding with RAC presence, just ~42% of USDA accessions bearing previously identified rhg1-a SoySNP50K SNP signatures harbor the RAC insertion. Subsequent analysis of several of these putative rhg1-a accessions lacking RAC revealed that none encoded α-SNAPRhg1LC, and thus, they are not rhg1-a. rhg1-a haplotypes are of rising interest, with Rhg4, for combating SCN populations that exhibit increased virulence against the widely used rhg1-b resistance. The present study reveals another unexpected structural feature of many Rhg1 loci, and a selectable feature that is predictive of rhg1-a haplotypes.

Whole-genome re-sequencing reveals the impact of the interaction of copy number variants of the rhg1 and Rhg4 genes on broad-based resistance to soybean cyst nematode

Soybean cyst nematode (SCN) is the most devastating plant-parasitic nematode. Most commercial soybean varieties with SCN resistance are derived from PI88788. Resistance derived from PI88788 is breaking down due to narrow genetic background and SCN population shift. PI88788 requires mainly the rhg1-b locus, while 'Peking' requires rhg1-a and Rhg4 for SCN resistance. In the present study, whole genome re-sequencing of 106 soybean lines was used to define the Rhg haplotypes and investigate their responses to the SCN HG-Types. The analysis showed a comprehensive profile of SNPs and copy number variations (CNV) at these loci. CNV of rhg1 (GmSNAP18) only contributed towards resistance in lines derived from PI88788 and 'Cloud'. At least 5.6 copies of the PI88788-type rhg1 were required to confer SCN resistance, regardless of the Rhg4 (GmSHMT08) haplotype. However, when the GmSNAP18 copies dropped below 5.6, a 'Peking'-type GmSHMT08 haplotype was required to ensure SCN resistance. This points to a novel mechanism of epistasis between GmSNAP18 and GmSHMT08 involving minimum requirements for copy number. The presence of more Rhg4 copies confers resistance to multiple SCN races. Moreover, transcript abundance of the GmSHMT08 in root tissue correlates with more copies of the Rhg4 locus, reinforcing SCN resistance. Finally, haplotype analysis of the GmSHMT08 and GmSNAP18 promoters inferred additional levels of the resistance mechanism. This is the first report revealing the genetic basis of broad-based resistance to SCN and providing new insight into epistasis, haplotype-compatibility, CNV, promoter variation and its impact on broad-based disease resistance in plants.

The poplar pangenome provides insights into the evolutionary history of the genus

The genus Populus comprises a complex amalgam of ancient and modern species that has become a prime model for evolutionary and taxonomic studies. Here we sequenced the genomes of 10 species from five sections of the genus Populus, identified 71 million genomic variations, and observed new correlations between the single-nucleotide polymorphism-structural variation (SNP-SV) density and indel-SV density to complement the SNP-indel density correlation reported in mammals. Disease resistance genes (R genes) with heterozygous loss-of-function (LOF) were significantly enriched in the 10 species, which increased the diversity of poplar R genes during evolution. Heterozygous LOF mutations in the self-incompatibility genes were closely related to the self-fertilization of poplar, suggestive of genomic control of self-fertilization in dioecious plants. The phylogenetic genome-wide SNPs tree also showed possible ancient hybridization among species in sections Tacamahaca, Aigeiros, and Leucoides. The pangenome resource also provided information for poplar genetics and breeding.

Structural variants in 3000 rice genomes

Investigation of large structural variants (SVs) is a challenging yet important task in understanding trait differences in highly repetitive genomes. Combining different bioinformatic approaches for SV detection, we analyzed whole-genome sequencing data from 3000 rice genomes and identified 63 million individual SV calls that grouped into 1.5 million allelic variants. We found enrichment of long SVs in promoters and an excess of shorter variants in 5′ UTRs. Across the rice genomes, we identified regions of high SV frequency enriched in stress response genes. We demonstrated how SVs may help in finding causative variants in genome-wide association analysis. These new insights into rice genome biology are valuable for understanding the effects SVs have on gene function, with the prospect of identifying novel agronomically important alleles that can be utilized to improve cultivated rice.

Copy Number Variation in Domestication

Domesticated plants have long served as excellent models for studying evolution. Many genes and mutations underlying important domestication traits have been identified, and most causal mutations appear to be SNPs. Copy number variation (CNV) is an important source of genetic variation that has been largely neglected in studies of domestication. Ongoing work demonstrates the importance of CNVs as a source of genetic variation during domestication, and during the diversification of domesticated taxa. Here, we review how CNVs contribute to evolutionary processes underlying domestication, and review examples of domestication traits caused by CNVs. We draw from examples in plant species, but also highlight cases in animal systems that could illuminate the roles of CNVs in the domestication process.

Connecting genome structural variation with complex traits in crop plants

Structural genome variation is a major determinant of useful trait diversity. We describe how genome analysis methods are enabling discovery of trait-associated structural variants and their potential impact on breeding. As our understanding of complex crop genomes continues to grow, there is growing evidence that structural genome variation plays a major role in determining traits important for breeding and agriculture. Identifying the extent and impact of structural variants in crop genomes is becoming increasingly feasible with ongoing advances in the sophistication of genome sequencing technologies, particularly as it becomes easier to generate accurate long sequence reads on a genome-wide scale. In this article, we discuss the origins of structural genome variation in crops from ancient and recent genome duplication and polyploidization events and review high-throughput methods to assay such variants in crop populations in order to find associations with phenotypic traits. There is increasing evidence from such studies that gene presence-absence and copy number variation resulting from segmental chromosome exchanges may be at the heart of adaptive variation of crops to counter abiotic and biotic stress factors. We present examples from major crops that demonstrate the potential of pangenomic diversity as a key resource for future plant breeding for resilience and sustainability.

Exploring and Exploiting Pan-genomics for Crop Improvement

Genetic variation ranging from single-nucleotide polymorphisms to large structural variants (SVs) can cause variation of gene content among individuals within the same species. There is an increasing appreciation that a single reference genome is insufficient to capture the full landscape of genetic diversity of a species. Pan-genome analysis offers a platform to evaluate the genetic diversity of a species via investigation of its entire genome repertoire. Although a recent wave of pan-genomic studies has shed new light on crop diversity and improvement using advanced sequencing technology, the potential applications of crop pan-genomics in crop improvement are yet to be fully exploited. In this review, we highlight the progress achieved in understanding crop pan-genomics, discuss biological activities that cause SVs, review important agronomical traits affected by SVs, and present our perspective on the application of pan-genomics in crop improvement.

Finding invisible quantitative trait loci with missing data

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

Leveraging RNA-Seq to Characterize Resistance to Brown Stem Rot and the Rbs3 Locus in Soybean

Brown stem rot, caused by the fungus Phialophora gregata, reduces soybean yield by up to 38%. Although three dominant resistance loci have been identified (Rbs1 to Rbs3), the gene networks responsible for pathogen recognition and defense remain unknown. Further, identification and characterization of resistant and susceptible germplasm remains difficult. We conducted RNA-Seq of infected and mock-infected leaf, stem, and root tissues of a resistant (PI 437970, Rbs3) and susceptible (Corsoy 79) genotype. Combining historical mapping data with genotype expression differences allowed us to identify a cluster of receptor-like proteins that are candidates for the Rbs3 resistance gene. Reads mapping to the Rbs3 locus were used to identify potential novel single-nucleotide polymorphisms within candidate genes that could improve phenotyping and breeding efficiency. Comparing responses to infection revealed little overlap in differential gene expression between genotypes or tissues. Gene networks associated with defense, DNA replication, and iron homeostasis are hallmarks of resistance to P. gregata. This novel research demonstrates the utility of combining contrasting genotypes, gene expression, and classical genetic studies to characterize complex disease resistance loci.

Detection of large sequence insertions by a hybrid approach that combine de novo assembly and resequencing of medium-coverage genome sequences

Large sequence insertion (LSI) is one of the structural variations (SVs) that may cause phenotypic differences in plants. To identify the LSIs using medium-coverage sequencing data of four wild soybean (Glycine soja) genotypes, we designed a hybrid approach combining de novo assembly and read mapping. Total reads and reads with both ends unmapped were independently assembled into "ordinary contigs" and "orphan contigs", respectively, and subjected to pairwise alignment and stringent filtering. This approach predicted 24 LSIs averaging 2682 bp in size, with no overlap with SVs detected by Pindel, BreakDancer, or ScanIndel, and they were validated by PCR. Compared with the soybean (Glycine max) reference genome, 20 LSIs were located outside genic regions. One of the four LSIs within a genic region, LSI05, is located in the coding DNA sequence region of a protein kinase superfamily gene (Glyma.08G123500). It caused delayed translation initiation and loss of 24 amino acids in the wild soybean genotype CW12. LSI05 was more frequently observed in 29 G. soja accessions than in 34 G. max accessions. Identified LSIs would be genomic resources harboring novel gene contents for studying SVs and improving crops. Moreover, our cost-efficient approach may be applicable to other plant species.

Gene copy number variations involved in balsam poplar (Populus balsamifera L.) adaptive variations

Gene copy number variations (CNVs) involved in phenotypic variations have already been shown in plants, but genomewide testing of CNVs for adaptive variation was not doable until recent technological developments. Thus, reports of the genomic architecture of adaptation involving CNVs remain scarce to date. Here, we investigated F₁ progenies of an intraprovenance cross (north-north cross, 58th parallel) and an interprovenances cross (north-south cross, 58th/49th parallels) for CNVs using comparative genomic hybridization on arrays of probes targeting gene sequences in balsam poplar (Populus balsamifera L.), a widespread North American forest tree. A total of 1,721 genes were found in varying copy numbers over the set of 19,823 tested genes. These gene CNVs presented an estimated average size of 8.3 kb and were distributed over poplar's 19 chromosomes including 22 hotspot regions. Gene CNVs number was higher for the interprovenance progeny in accordance with an expected higher genetic diversity related to the composite origin of this family. Regression analyses between gene CNVs and seven adaptive trait variations resulted in 23 significant links; among these adaptive gene CNVs, 30% were located in hotspots. One-to-five gene CNVs were found related to each of the measured adaptive traits and annotated for both biotic and abiotic stress responses. These annotations can be related to the occurrence of a higher pathogenic pressure in the southern parts of balsam poplar's distribution, and higher photosynthetic assimilation rates and water-use efficiency at high latitudes. Overall, our findings suggest that gene CNVs typically having higher mutation rates than SNPs may in fact represent efficient adaptive variations against fast-evolving pathogens.

NLR Mutations Suppressing Immune Hybrid Incompatibility and Their Effects on Disease Resistance

Genetic divergence between populations can lead to reproductive isolation. Hybrid incompatibilities (HI) represent intermediate points along a continuum toward speciation. In plants, genetic variation in disease resistance (R) genes underlies several cases of HI. The progeny of a cross between Arabidopsis (Arabidopsis thaliana) accessions Landsberg erecta (Ler, Poland) and Kashmir2 (Kas2, central Asia) exhibits immune-related HI. This incompatibility is due to a genetic interaction between a cluster of eight TNL (TOLL/INTERLEUKIN1 RECEPTOR-NUCLEOTIDE BINDING-LEU RICH REPEAT) RPP1 (RECOGNITION OF PERONOSPORA PARASITICA1)-like genes (R1-R8) from Ler and central Asian alleles of a Strubbelig-family receptor-like kinase (SRF3) from Kas2. In characterizing mutants altered in Ler/Kas2 HI, we mapped multiple mutations to the RPP1-like Ler locus. Analysis of these suppressor of Ler/Kas2 incompatibility (sulki) mutants reveals complex, additive and epistatic interactions underlying RPP1-like Ler locus activity. The effects of these mutations were measured on basal defense, global gene expression, primary metabolism, and disease resistance to a local Hyaloperonospora arabidopsidis isolate (Hpa Gw) collected from Gorzów (Gw), where the Landsberg accession originated. Gene expression sectors and metabolic hallmarks identified for HI are both dependent and independent of RPP1-like Ler members. We establish that mutations suppressing immune-related Ler/Kas2 HI do not compromise resistance to Hpa Gw. QTL mapping analysis of Hpa Gw resistance point to RPP7 as the causal locus. This work provides insight into the complex genetic architecture of the RPP1-like Ler locus and immune-related HI in Arabidopsis and into the contributions of RPP1-like genes to HI and defense.

Mapping of a Major QTL for Ceratocystis Wilt Disease in an F1 Population of Theobroma cacao

Cacao is an important crop, its beans are key raw materials for the chocolate and cosmetic industries. Ceratocystis wilt of cacao (CWC) caused by Ceratocystis cacaofunesta is a lethal disease for the crop. Therefore, the selection of resistant cacao varieties is one of the viable ways to minimize losses in cacao production. In this paper, we described the identification of a major QTL associated with CWC in an F1 mapping population from a cross between a resistant, "TSH 1188," and a susceptible genotype, "CCN 51." A set of 266 trees were genotyped using 3,526 single nucleotide polymorphic markers and then multiple QTL mapping analyses were performed. Two QTLs were identified on chromosomes IV and VI. The major QTL was located at 20 cM from the top position of chromosome VI, accounting for more than 60% of the phenotypic variation. The favorable allele T1, with haplotype GTT, came from the "TSH 1188" parent. It was evident that the haplotype combination T1C2 on chromosome VI was the most significant for resistance, since 93% of resistant trees had this haplotype. The major QTL converged to a genomic region of 739.4 kb that harbored nine candidate genes, including two major classes of resistance genes, which would make them the primary candidates involved in the resistance to CWC. The haplotypes detected are now used to improve the efficiency and precision of the selection of resistant trees in cacao breeding.