Copy number variation and disease resistance in plants

The mutagenic effect of cold plasma on Medicago sativa L

To explore the feasibility of using cold plasma as a mutagenesis breeding technology for forage crops, in this study we used the Medicago sativa L. cultivar, Zhongcao No. 3, as the experimental material. The effects of plasma treatments on Medicago sativa L. were analyzed through the use of plasma and activated water. Treatments with plasma and activated water inhibited plant height but promoted root growth. By creating a closed environment, adding a dielectric barrier plate, and combining these two treatment methods, the greatest impact can be had on the growth of Medicago sativa L. seeds. After treatment, the plant heights were approximately half those of the control group, and the root lengths were approximately 1.6 times those of the control group. Through emission spectroscopy, it was found that active particles such as O, NO2, and N2* were present and could be considered to have produced plasma-activated water through contact with the water surface, thus affecting the survival and growth of the seeds. Whole-genome resequencing (WGRS) was performed on the wild-type and selected mutants after treatment, with an average sequencing depth of 115.93×, an average genome alignment rate of 91.72 %, and an average genome coverage rate of 91.85 %. Various types of mutations were detected and annotated. After filtering, 7,822,324 SNP (single nucleotide polymorphisms) sites, 2,161,917 indel sites, 200,544 SV sites, 238 CNV (copy number variation) sites. The SNPs, indels (insertions/deletions), and SVs (structural variations) were mainly heterozygous, with heterozygosity rates of 87.13 %, 92.16 %, and 83.49 %, respectively. The CNVs were dominated by low copy numbers, accounting for 53.77 %. These results indicate that plasma treatment has significant effects on plant growth and genome of Medicago sativa L.

The genome of Lupinus mutabilis: Evolution and genetics of an emerging bio-based crop

Lupinus mutabilis is an under-domesticated legume species from the Andean region of South America. It belongs to the New World lupins clade, which groups several lupin species displaying large genetic variation and adaptability to highly different environments. L. mutabilis is attracting interest as a potential multipurpose crop to diversify the European supply of plant proteins, increase agricultural biodiversity, and fulfill bio-based applications. This study reports the first high-quality L. mutabilis genome assembly, which is also the first sequenced assembly of a New World lupin species. Through comparative genomics and phylogenetics, the evolution of L. mutabilis within legumes and lupins is described, highlighting both genomic similarities and patterns specific to L. mutabilis, potentially linked to environmental adaptations. Furthermore, the assembly was used to study the genetics underlying important traits for the establishment of L. mutabilis as a novel crop, including protein and quinolizidine alkaloids contents in seeds, genomic patterns of classic resistance genes, and genomic properties of L. mutabilis mycorrhiza-related genes. These analyses pointed out copy number variation, differential genomic gene contexts, and gene family expansion through tandem duplications as likely important drivers of the genomic diversity observed for these traits between L. mutabilis and other lupins and legumes. Overall, the L. mutabilis genome assembly will be a valuable resource to conduct genetic research and enable genomic-based breeding approaches to turn L. mutabilis into a multipurpose legume crop.

Genetic diversity, disease resistance, and environmental adaptation of Arachis duranensis L.: New insights from landscape genomics

The genetic diversity that exists in natural populations of Arachis duranensis, the wild diploid donor of the A subgenome of cultivated tetraploid peanut, has the potential to improve crop adaptability, resilience to major pests and diseases, and drought tolerance. Despite its potential value for peanut improvement, limited research has been focused on the association between allelic variation, environmental factors, and response to early (ELS) and late leaf spot (LLS) diseases. The present study implemented a landscape genomics approach to gain a better understanding of the genetic variability of A. duranensis represented in the ex-situ peanut germplasm collection maintained at the U.S. Department of Agriculture, which spans the entire geographic range of the species in its center of origin in South America. A set of 2810 single nucleotide polymorphism (SNP) markers allowed a high-resolution genome-wide characterization of natural populations. The analysis of population structure showed a complex pattern of genetic diversity with five putative groups. The incorporation of bioclimatic variables for genotype-environment associations, using the latent factor mixed model (LFMM2) method, provided insights into the genomic signatures of environmental adaptation, and led to the identification of SNP loci whose allele frequencies were correlated with elevation, temperature, and precipitation-related variables (q < 0.05). The LFMM2 analysis for ELS and LLS detected candidate SNPs and genomic regions on chromosomes A02, A03, A04, A06, and A08. These findings highlight the importance of the application of landscape genomics in ex situ collections of peanut and other crop wild relatives to effectively identify favorable alleles and germplasm for incorporation into breeding programs. We report new sources of A. duranensis germplasm harboring adaptive allelic variation, which have the potential to be utilized in introgression breeding for a single or multiple environmental factors, as well as for resistance to leaf spot diseases.

Evolutionary and synteny analysis of HIS1, BADH2, GBSS1, and GBSS2 in rice: insights for effective introgression breeding strategies

The key genes BADH2, GBSS1, GBSS2, and HIS1 regulate the fragrance, starch synthesis, and herbicide resistance in rice. Although the molecular functions of four genes have been investigated in the Oryza sativa species, little is known regarding their evolutionary history in the Oryza genus. Here, we studied the evolution of four focal genes in 10 Oryza species using phylogenetic and syntenic approaches. The HIS1 family underwent several times of tandem duplication events in the Oryza species, resulting in copy number variation ranging from 2 to 7. At most one copy of BADH2, GBSS1, and GBSS2 orthologs were identified in each Oryza species, and gene loss events of BADH2 and GBSS2 were identified in three Oryza species. Gene transfer analysis proposed that the functional roles of GBSS1 and GBSS2 were developed in the Asian and African regions, respectively, and most allelic variations of BADH2 in japonica rice emerged after the divergence between the Asian and African rice groups. These results provide clues to determine the origin and evolution of the key genes in rice breeding as well as valuable information for molecular breeders and scientists to develop efficient strategies to simultaneously improve grain quality and yield potential in rice.

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

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

The Triticeae CBF Gene Cluster-To Frost Resistance and Beyond

The pivotal role of CBF/DREB1 transcriptional factors in Triticeae crops involved in the abiotic stress response has been highlighted. The CBFs represent an important hub in the ICE-CBF-COR pathway, which is one of the most relevant mechanisms capable of activating the adaptive response to cold and drought in wheat, barley, and rye. Understanding the intricate mechanisms and regulation of the cluster of CBF genes harbored by the homoeologous chromosome group 5 entails significant potential for the genetic improvement of small grain cereals. Triticeae crops seem to share common mechanisms characterized, however, by some peculiar aspects of the response to stress, highlighting a combined landscape of single-nucleotide variants and copy number variation involving CBF members of subgroup IV. Moreover, while chromosome 5 ploidy appears to confer species-specific levels of resistance, an important involvement of the ICE factor might explain the greater tolerance of rye. By unraveling the genetic basis of abiotic stress tolerance, researchers can develop resilient varieties better equipped to withstand extreme environmental conditions. Hence, advancing our knowledge of CBFs and their interactions represents a promising avenue for improving crop resilience and food security.

Subtelomeric 5-enolpyruvylshikimate-3-phosphate synthase copy number variation confers glyphosate resistance in Eleusine indica

Genomic structural variation (SV) has profound effects on organismal evolution; often serving as a source of novel genetic variation. Gene copy number variation (CNV), one type of SV, has repeatedly been associated with adaptive evolution in eukaryotes, especially with environmental stress. Resistance to the widely used herbicide, glyphosate, has evolved through target-site CNV in many weedy plant species, including the economically important grass, Eleusine indica (goosegrass); however, the origin and mechanism of these CNVs remain elusive in many weed species due to limited genetic and genomic resources. To study this CNV in goosegrass, we present high-quality reference genomes for glyphosate-susceptible and -resistant goosegrass lines and fine-assembles of the duplication of glyphosate's target site gene 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS). We reveal a unique rearrangement of EPSPS involving chromosome subtelomeres. This discovery adds to the limited knowledge of the importance of subtelomeres as genetic variation generators and provides another unique example for herbicide resistance evolution.

M. Drummond E, Van Tassel DL, Warschefsky EJ. The next era of crop domestication starts now

Current food systems are challenged by relying on a few input-intensive, staple crops. The prioritization of yield and the loss of diversity during the recent history of domestication has created contemporary crops and cropping systems that are ecologically unsustainable, vulnerable to climate change, nutrient poor, and socially inequitable. For decades, scientists have proposed diversity as a solution to address these challenges to global food security. Here, we outline the possibilities for a new era of crop domestication, focused on broadening the palette of crop diversity, that engages and benefits the three elements of domestication: crops, ecosystems, and humans. We explore how the suite of tools and technologies at hand can be applied to renew diversity in existing crops, improve underutilized crops, and domesticate new crops to bolster genetic, agroecosystem, and food system diversity. Implementing the new era of domestication requires that researchers, funders, and policymakers boldly invest in basic and translational research. Humans need more diverse food systems in the Anthropocene-the process of domestication can help build them.

Occurrence and mechanism of target-site resistance to bensulfuron-methyl in Monochoria korsakowii from China

Monochoria korsakowii is an increasingly significant threat to rice production across China, particularly in Liaoning province. Few studies have reported herbicide resistance in M. korsakowii, and resistance status and mechanisms are poorly understood. Here, thirty field populations of M. korsakowii were collected from 11 rice-growing regions of Liaoning, and 97% of populations had evolved resistance to bensulfuron-methyl (BM), with majority (24 of 28) showing high resistance levels (RI > 10). The first in-depth analysis of molecular features of AHAS1 and AHAS2 in BM-resistant populations showed that four Pro197 mutations (Pro197 to His, Ala, Leu or Ser) in AHAS1 and one mutation (Pro197Ser) in AHAS2 were identified. Notably, novel double Pro197Ser mutations co-occurred in both AHAS1 and AHAS2 in the most resistant line LN-20. Furthermore, resistant mutants were used to investigate the effect of Pro197 mutations on AHAS functionality, binding modes, gene expression and cross-resistance in M. korsakowii. All the detected Pro197 mutations considerably reduced in vitro AHAS sensitivity to BM by weakening hydrogen bonds and hydrophobic interactions in the predicted BM-AHAS complexes, especially the double Pro197Ser mutations. This novel resistance mutation combination slightly impacted the extractable AHAS activity, and increased the affinity and catalytic rate of pyruvate. Also, the AHAS expression level was significantly up-regulated. Moreover, all mutations provided resistance only to other sulfonylureas herbicides but not triazolopyrimidine or pyrimidinyl-benzoates herbicides. In conclusion, bensulfuron-methyl resistance in M. korsakowii was grim in Liaoning, China, and amino acid mutations on AHAS isozymes were the primary resistance mechanism. Double Pro197Ser mutations in both AHAS1 and AHAS2 confer higher herbicide resistance than single mutations in AHAS1. Thus, this work deepens our understanding of resistance status and mechanisms of M. korsakowii.

Genomic and Molecular Characterization of Wheat Streak Mosaic Virus Resistance Locus 2 (Wsm2) in Common Wheat (Triticum aestivum L.)

Wheat streak mosaic virus (WSMV) is an economically important viral pathogen that threatens global wheat production, particularly in the Great Plains of the United States. The Wsm2 locus confers resistance to WSMV and has been widely deployed in common wheat varieties adapted to this region. Characterizing the underlying causative genetic variant would contribute to our understanding of viral resistance mechanisms in wheat and aid the development of perfect markers for breeding. In this study, linkage mapping in a doubled-haploid (DH) mapping population confirmed Wsm2 as a major locus conferring WSMV resistance in wheat. The Wsm2 flanking markers were mapped to a 4.0 Mbp region at the distal end of chromosome 3BS containing 142 candidate genes. Eight haplotypes were identified from seventeen wheat genotypes collected from different agroecological zones, indicating that Wsm2 lies in a dynamic region of the genome with extensive structural variation and that it is likely a rare allele in most available genome assemblies of common wheat varieties. Exome sequencing of the variety "Snowmass", which carries Wsm2, revealed several loss-of-function mutations and copy number variants in the 142 candidate genes within the Wsm2 interval. Six of these genes are differentially expressed in "Snowmass" compared to "Antero," a variety lacking Wsm2, including a gene that encodes a nucleotide-binding site leucine-rich repeat (NBS-LRR) type protein with homology to RPM1. A de novo assembly of unmapped RNA-seq reads identified nine transcripts expressed only in "Snowmass," three of which are also induced in response to WSMV inoculation. This study sheds light on the variation underlying Wsm2 and provides a list of candidate genes for subsequent validation.

Novel Fusarium wilt resistance genes uncovered in natural and cultivated strawberry populations are found on three non-homoeologous chromosomes

Several Fusarium wilt resistance genes were discovered, genetically and physically mapped, and rapidly deployed via marker-assisted selection to develop cultivars resistant to Fusarium oxysporum f. sp. fragariae, a devastating soil-borne pathogen of strawberry. Fusarium wilt, a soilborne disease caused by Fusarium oxysporum f. sp. fragariae, poses a significant threat to strawberry (Fragaria [Formula: see text] ananassa) production in many parts of the world. This pathogen causes wilting, collapse, and death in susceptible genotypes. We previously identified a dominant gene (FW1) on chromosome 2B that confers resistance to race 1 of the pathogen, and hypothesized that gene-for-gene resistance to Fusarium wilt was widespread in strawberry. To explore this, a genetically diverse collection of heirloom and modern cultivars and octoploid ecotypes were screened for resistance to Fusarium wilt races 1 and 2. Here, we show that resistance to both races is widespread in natural and domesticated populations and that resistance to race 1 is conferred by partially to completely dominant alleles among loci (FW1, FW2, FW3, FW4, and FW5) found on three non-homoeologous chromosomes (1A, 2B, and 6B). The underlying genes have not yet been cloned and functionally characterized; however, plausible candidates were identified that encode pattern recognition receptors or other proteins known to confer gene-for-gene resistance in plants. High-throughput genotyping assays for SNPs in linkage disequilibrium with FW1-FW5 were developed to facilitate marker-assisted selection and accelerate the development of race 1 resistant cultivars. This study laid the foundation for identifying the genes encoded by FW1-FW5, in addition to exploring the genetics of resistance to race 2 and other races of the pathogen, as a precaution to averting a Fusarium wilt pandemic.

Identification and expression analysis of arabinogalactan protein genes in cotton reveal the function of GhAGP15 in Verticillium dahliae resistance

Verticillium wilt, primarily caused by the fungal pathogen Verticillium dahliae, is a serious disease in cotton. Arabinogalactan proteins (AGPs), a class of hydroxyproline-rich glycoproteins, have been widely implicated in plant growth and environmental adaptation. The purpose of this study is to identify and characterize AGP members in cotton plants and explore their roles in responding to environmental stressors. In total, 65 GhAGP members were identified in upland cotton (Gossypium hirsutum), along with 43, 35, and 37 AGP members that were also identified in G. barbadense, G. arboreum, and G. raimondii, respectively. According to gene structure and protein domains analysis, GhAGP genes in upland cotton are highly conserved. Meanwhile, tandem duplication events have occurred frequently throughout cotton's evolutionary history. Expression analysis showed that GhAGP genes were widely expressed during growth and development and in response to abiotic stressors. Many cis-elements related to hormonal responses and environmental stressors were detected in GhAGP promoter regions. GhAGP genes participate in responding to cold, drought, and salt stress, and were sensitive to ET signaling. Furthermore, the expression level of GhAGP15 was elevated during V. dahliae infection and resistance against V. dahliae in upland cotton was significantly weakened by silencing GhAGP15 using a virus-induced gene silencing (VIGS) approach. Our results further suggest that the function of GhAGP15 in V. dahliae resistance might be involved in regulation of the JA, SA, and reactive oxygen species (ROS) pathways. The comprehensive analysis of AGP genes in cotton performed in this study provides a basic framework for further functional research of these genes.

Oxford Nanopore and Bionano Genomics technologies evaluation for plant structural variation detection

BACKGROUND

Structural Variations (SVs) are genomic rearrangements derived from duplication, deletion, insertion, inversion, and translocation events. In the past, SVs detection was limited to cytological approaches, then to Next-Generation Sequencing (NGS) short reads and partitioned assemblies. Nowadays, technologies such as DNA long read sequencing and optical mapping have revolutionized the understanding of SVs in genomes, due to the enhancement of the power of SVs detection. This study aims to investigate performance of two techniques, 1) long-read sequencing obtained with the MinION device (Oxford Nanopore Technologies) and 2) optical mapping obtained with Saphyr device (Bionano Genomics) to detect and characterize SVs in the genomes of the two ecotypes of Arabidopsis thaliana, Columbia-0 (Col-0) and Landsberg erecta 1 (Ler-1).

RESULTS

We described the SVs detected from the alignment of the best ONT assembly and DLE-1 optical maps of A. thaliana Ler-1 against the public reference genome Col-0 TAIR10.1. After filtering (SV > 1 kb), 1184 and 591 Ler-1 SVs were retained from ONT and Bionano technologies respectively. A total of 948 Ler-1 ONT SVs (80.1%) corresponded to 563 Bionano SVs (95.3%) leading to 563 common locations. The specific locations were scrutinized to assess improvement in SV detection by either technology. The ONT SVs were mostly detected near TE and gene features, and resistance genes seemed particularly impacted.

CONCLUSIONS

Structural variations linked to ONT sequencing error were removed and false positives limited, with high quality Bionano SVs being conserved. When compared with the Col-0 TAIR10.1 reference genome, most of the detected SVs discovered by both technologies were found in the same locations. ONT assembly sequence leads to more specific SVs than Bionano one, the latter being more efficient to characterize large SVs. Even if both technologies are complementary approaches, ONT data appears to be more adapted to large scale populations studies, while Bionano performs better in improving assembly and describing specificity of a genome compared to a reference.

Construction of a Quantitative Genomic Map, Identification and Expression Analysis of Candidate Genes for Agronomic and Disease-Related Traits in Brassica napus

Rapeseed is the second most important oil crop in the world. Improving seed yield and seed oil content are the two main highlights of the research. Unfortunately, rapeseed development is frequently affected by different diseases. Extensive research has been made through many years to develop elite cultivars with high oil, high yield, and/or disease resistance. Quantitative trait locus (QTL) analysis has been one of the most important strategies in the genetic deciphering of agronomic characteristics. To comprehend the distribution of these QTLs and to uncover the key regions that could simultaneously control multiple traits, 4,555 QTLs that have been identified during the last 25 years were aligned in one unique map, and a quantitative genomic map which involved 128 traits from 79 populations developed in 12 countries was constructed. The present study revealed 517 regions of overlapping QTLs which harbored 2,744 candidate genes and might affect multiple traits, simultaneously. They could be selected to customize super-rapeseed cultivars. The gene ontology and the interaction network of those candidates revealed genes that highly interacted with the other genes and might have a strong influence on them. The expression and structure of these candidate genes were compared in eight rapeseed accessions and revealed genes of similar structures which were expressed differently. The present study enriches our knowledge of rapeseed genome characteristics and diversity, and it also provided indications for rapeseed molecular breeding improvement in the future.

Genomic Variations and Mutational Events Associated with Plant-Pathogen Interactions

Phytopathologists are actively researching the molecular basis of plant-pathogen interactions. The mechanisms of responses to pathogens have been studied extensively in model crop plant species and natural populations. Today, with the rapid expansion of genomic technologies such as DNA sequencing, transcriptomics, proteomics, and metabolomics, as well as the development of new methods and protocols, data analysis, and bioinformatics, it is now possible to assess the role of genetic variation in plant-microbe interactions and to understand the underlying molecular mechanisms of plant defense and microbe pathogenicity with ever-greater resolution and accuracy. Genetic variation is an important force in evolution that enables organisms to survive in stressful environments. Moreover, understanding the role of genetic variation and mutational events is essential for crop breeders to produce improved cultivars. This review focuses on genetic variations and mutational events associated with plant-pathogen interactions and discusses how these genome compartments enhance plants' and pathogens' evolutionary processes.

Dual activity of Meloidogyne incognita-regulated Musa acuminata Pathogenesis-related-10 (MaPR-10) gene

Plant immunity to pathogen infections is a dynamic response that involves multiple organelles and defence signalling systems such as hypersensitive response (HR) and systemic acquired resistance (SAR). The latter requires the function of Pathogenesis-related (PR) proteins, a common plant protein family with diverse roles in plant innate immunity. Our previous proteomics study showed that a PR gene (ITC1587_Bchr9_P26466_MUSBA) was differentially regulated during a compatible banana-M. incognita interaction, substantiating the isolation of this gene in the current study. Here, we successfully isolated and characterised Pathogenesis-related-10 (PR10) gene with β-1,3-glucanase and ribonuclease (RNase) activities from two Musa acuminata cultivars (denoted as MaPR10) namely Berangan and Grand Naine (ITC1256). We found that MaPR10 cloned sequences possess glycine-rich loop domain and shared conserved motifs specific to PR10 gene group, confirming its identity as a member of this group. Interestingly, we also found a catalytic domain sequence for glycoside hydrolase family 16 (EXDXXE), unique only to MaPR10 cloned sequences. Two peptide variants closely related to the reference sequence ITC1587_Bchr9_P26466_MUSBA namely MaPR10-BeB5 and MaPR10-GNA5 were overexpressed and purified to test for their functionality. Here, we confirmed that both protein variants possess β-1,3-glucanase and ribonuclease (RNase) activities, and inhibit the growth of Aspergillus fumigatus, a human opportunistic pathogen. To our knowledge, this is the first PR10 plant proteins with such properties to be reported thus far.

Genome-wide analysis of deletions in maize population reveals abundant genetic diversity and functional impact

Two read depth methods were jointly used in next-generation sequencing data to identify deletions in maize population. GWAS by deletions were analyzed for gene expression pattern and classical traits, respectively. Many studies have confirmed that structural variation (SV) is pervasive throughout the maize genome. Deletion is one type of SV that may impact gene expression and cause phenotypic changes in quantitative traits. In this study, two read count approaches were used to analyze the deletions in the whole-genome sequencing data of 270 maize inbred lines. A total of 19,754 deletion windows overlapped 12,751 genes, which were unevenly distributed across the genome. The deletions explained population structure well and correlated with genomic features. The deletion proportion of genes was determined to be negatively correlated with its expression. The detection of gene expression quantitative trait loci (eQTL) indicated that local eQTL were fewer but had larger effects than distant ones. The common associated genes were related to basic metabolic processes, whereas unique associated genes with eQTL played a role in the stress or stimulus responses in multiple tissues. Compared with the eQTL detected by SNPs derived from the same sequencing data, 89.4% of the associated genes could be detected by both markers. The effect of top eQTL detected by SNPs was usually larger than that detected by deletions for the same gene. A genome-wide association study (GWAS) on flowering time and plant height illustrated that only a few loci could be consistently captured by SNPs, suggesting that combining deletion and SNP for GWAS was an excellent strategy to dissect trait architecture. Our findings will provide insights into characteristic and biological function of genome-wide deletions in maize.

High-quality reference genome and annotation aids understanding of berry development for evergreen blueberry (Vaccinium darrowii)

Vaccinium darrowii Camp (2n = 2x = 24) is a native North American blueberry species and an important source of traits such as low chill requirement in commercial southern highbush blueberry breeding (Vaccinium corymbosum, 2n = 4x = 48). We present a chromosomal-scale genome of V. darrowii generated by the combination of PacBio sequencing and high throughput chromatin conformation capture (Hi-C) scaffolding technologies, yielding a total length of 1.06 Gigabases (Gb). Over 97.8% of the genome sequences are scaffolded into 24 chromosomes representing the two haplotypes. The primary haplotype assembly of V. darrowii contains 34,809 protein-coding genes. Comparison to a V. corymbosum haplotype assembly reveals high collinearity between the two genomes with small intrachromosomal rearrangements in eight chromosome pairs. With small RNA sequencing, the annotation was further expanded to include more than 200,000 small RNA loci and 638 microRNAs expressed in berry tissues. Transcriptome analysis across fruit development stages indicates that genes involved in photosynthesis are downregulated, while genes involved in flavonoid and anthocyanin biosynthesis are significantly increased at the late stage of berry ripening. A high-quality reference genome and accompanying annotation of V. darrowii is a significant new resource for assessing the evergreen blueberry contribution to the breeding of southern highbush blueberries.

Current status of structural variation studies in plants

Structural variations (SVs) including gene presence/absence variations and copy number variations are a common feature of genomes in plants and, together with single nucleotide polymorphisms and epigenetic differences, are responsible for the heritable phenotypic diversity observed within and between species. Understanding the contribution of SVs to plant phenotypic variation is important for plant breeders to assist in producing improved varieties. The low resolution of early genetic technologies and inefficient methods have previously limited our understanding of SVs in plants. However, with the rapid expansion in genomic technologies, it is possible to assess SVs with an ever-greater resolution and accuracy. Here, we review the current status of SV studies in plants, examine the roles that SVs play in phenotypic traits, compare current technologies and assess future challenges for SV studies.

Defence by duplication: The relation between phenotypic glyphosate resistance and EPSPS gene copy number variation in Amaranthus palmeri

Gene copy number variation (CNV) has been increasingly associated with organismal responses to environmental stress, but we know little about the quantitative relation between CNV and phenotypic variation. In this study we quantify the relation between variation in EPSPS (5-enolpyruvylshikimate-3-phosphate synthase) copy number using digital drop PCR and variation in phenotypic glyphosate resistance in 22 populations of Amaranthus palmeri (Palmer Amaranth), a range-expanding agricultural weed. Overall, we detected a significant positive relation between population mean copy number and resistance. The majority of populations exhibited high glyphosate resistance yet maintained low-resistance individuals, resulting in bimodality in many populations. We also investigated threshold models for the relation between copy number and resistance, and found evidence for a threshold of ~15 EPSPS copies: there was a steep increase in resistance below the threshold, followed by a much shallower increase. Across 924 individuals, as copy number increased the range of variation in resistance decreased, yielding an increasing frequency of high phenotypic resistance individuals. Among populations we detected a decline in variation (s.d.) as mean phenotypic resistance increased from moderate to high, consistent with the prediction that as phenotypic resistance increases in populations, stabilizing selection decreases variation in the trait. Our study demonstrates that populations of A. palmeri can harbour wide variation in EPSPS copy number and phenotypic glyphosate resistance, reflecting the history of, and template for future, resistance evolution.

High resolution estimates of relative gene abundance with quantitative ratiometric regression PCR (qRR-PCR)

Quantification of the relative abundance of genetic traits has broad applications for biomarker discovery, diagnostics, and assessing gene expression in biological research. Relative quantification of genes is traditionally done with the 2-ΔΔCT method using quantitative real-time polymerase chain reaction (qPCR) data, which is often limited in resolution beyond orders of magnitude difference. The latest techniques for quantification of nucleic acids employ digital PCR or microarrays which involve lengthy sample preparation and complex instrumentation. In this work, we describe a quantitative ratiometric regression PCR (qRR-PCR) method for computing relative abundance of genetic traits in a sample with high resolution from a single duplexed real-time quantitative PCR assay. Instead of comparing the individual cycle threshold (Ct) values as is done for the 2-ΔΔCT method, our qRR-PCR algorithm leverages the innate relationship of co-amplified PCR targets to measure their relative quantities using characteristic curves derived from the normalized ratios of qPCR fluorescence curves. We demonstrate the utility of this technique for discriminating the fractional abundance of mixed alleles with resolution below 5%.

Genome plasticity in Paramecium bursaria revealed by population genomics

Background

Ciliates are an ancient and diverse eukaryotic group found in various environments. A unique feature of ciliates is their nuclear dimorphism, by which two types of nuclei, the diploid germline micronucleus (MIC) and polyploidy somatic macronucleus (MAC), are present in the same cytoplasm and serve different functions. During each sexual cycle, ciliates develop a new macronucleus in which newly fused genomes are extensively rearranged to generate functional minichromosomes. Interestingly, each ciliate species seems to have its way of processing genomes, providing a diversity of resources for studying genome plasticity and its regulation. Here, we sequenced and analyzed the macronuclear genome of different strains of Paramecium bursaria, a highly divergent species of the genus Paramecium which can stably establish endosymbioses with green algae.

Results

We assembled a high-quality macronuclear genome of P. bursaria and further refined genome annotation by comparing population genomic data. We identified several species-specific expansions in protein families and gene lineages that are potentially associated with endosymbiosis. Moreover, we observed an intensive chromosome breakage pattern that occurred during or shortly after sexual reproduction and contributed to highly variable gene dosage throughout the genome. However, patterns of copy number variation were highly correlated among genetically divergent strains, suggesting that copy number is adjusted by some regulatory mechanisms or natural selection. Further analysis showed that genes with low copy number variation among populations tended to function in basic cellular pathways, whereas highly variable genes were enriched in environmental response pathways.

Conclusions

We report programmed DNA rearrangements in the P. bursaria macronuclear genome that allow cells to adjust gene copy number globally according to individual gene functions. Our results suggest that large-scale gene copy number variation may represent an ancient mechanism for cells to adapt to different environments.

Supplementary information

The online version contains supplementary material available at 10.1186/s12915-020-00912-2.

Whole-genome sequencing reveals uniqueness of black-hulled and straw-hulled weedy rice genomes

Both SH and BHA weedy rice genotypes evolved independently and have distinct genomic composition. Different genetic mechanisms may be responsible for their competitiveness and adaptation to diverse environmental conditions. Two major types of weedy rice are recognized in the USA based on morphology: straw-hull (SH) and black-hull awned (BHA) weedy rice. We performed whole-genome resequencing of a SH weedy rice 'PSRR-1', a BHA weedy rice 'BHA1115', and a japonica cultivar 'Cypress' to delineate genome-wide differences and their relevance to genetics and evolution of weedy attributes. The high-quality reads were uniformly distributed with 82-88% genome coverage. The number of genotype-specific SNPs and InDels was highest in Cypress, followed by BHA1115 and PSRR-1. However, more genes were affected in BHA1115 compared with other two genotypes which is evident from the number of high-impact SNPs and InDels. Haplotype analysis of selected genes involved in domestication, adaptation, and agronomic performance not only differentiated SH from BHA weedy rice and supported evolution of weedy rice through de-domestication, but also validated the function of several genes such as qAn-1, qAn-2, Bh4, Rc, SD1, OsLG1, and OsC1. Several candidate genes were identified for previously reported seed dormancy and seed shattering QTLs. The SH and BHA weedy rice have distinct genomic composition, and the BHA weedy rice likely diverged earlier than SH weedy rice. The accumulation of plant development, reproduction, and defense-related genes in weedy rice possibly helped them to compete, survive, and spread under a wide range of environmental conditions by employing novel and diverse mechanisms. The genomic resources will be useful for both weed management and rice improvement by exploring the molecular basis of key agronomic, adaptive, and domestication attributes.

Influence of CNV on transcript levels of HvCBF genes at Fr-H2 locus revealed by resequencing in resistant barley cv. 'Nure' and expression analysis

Resequencing in resistant cultivar 'Nure' and structural comparison with the same region of susceptible 'Morex' was performed in order to gain a better insight into barley Frost-resistance-H2 locus. Accurate annotation showed copy number variation (CNV) in the proximal part of the locus. In 'Nure', two exact copies of the HvCBF4-HvCBF2A region and one of the HvCBF4-HvCBF2B segment were observed, while in 'Morex' the corresponding region harboured a single HvCBF4-HvCBF2A (22 kb) segment. Abundance and diversity of repetitive element classes, gene function gain/losses, regulatory motifs and SNPs in gene sequences were identified. An expression study of key HvCBFs with/without CNV on selected genotypes contrasting for frost resistance and estimated HvCBF4-HvCBF2B copy number (2-10 copies) was also performed. Under light stimulus at warm temperature (23 °C), CNV of HvCBF2A and HvCBF4 correlated with their expression levels and reported frost resistance of genotypes; moreover, expression levels of HvCBF2A and HvCBF14 were strongly correlated (r = 0.908, p < 0.01). On the other hand, frost resistance correlated to HvCBF14 expression (r = 0.871, p < 0.01) only after cold induction (6°C) in the dark. A complex interplay of HvCBFs expression levels under different light/temperature stimuli is discussed in light of CNV and presence/number of regulatory elements that integrate different signal transduction pathways.

Pangenomics Comes of Age: From Bacteria to Plant and Animal Applications

The pangenome refers to a collection of genomic sequence found in the entire species or population rather than in a single individual; the sequence can be core, present in all individuals, or accessory (variable or dispensable), found in a subset of individuals only. While pangenomic studies were first undertaken in bacterial species, developments in genome sequencing and assembly approaches have allowed construction of pangenomes for eukaryotic organisms, fungi, plants, and animals, including two large-scale human pangenome projects. Analysis of the these pangenomes revealed key differences, most likely stemming from divergent evolutionary histories, but also surprising similarities.

Surviving in a Hostile World: Plant Strategies to Resist Pests and Diseases

As primary producers, plants are under constant pressure to defend themselves against potentially deadly pathogens and herbivores. In this review, we describe short- and long-term strategies that enable plants to cope with these stresses. Apart from internal immunological strategies that involve physiological and (epi)genetic modifications at the cellular level, plants also employ external strategies that rely on recruitment of beneficial organisms. We discuss these strategies along a gradient of increasing timescales, ranging from rapid immune responses that are initiated within seconds to (epi)genetic adaptations that occur over multiple plant generations. We cover the latest insights into the mechanistic and evolutionary underpinnings of these strategies and present explanatory models. Finally, we discuss how knowledge from short-lived model species can be translated to economically and ecologically important perennials to exploit adaptive plant strategies and mitigate future impacts of pests and diseases in an increasingly interconnected and changing world.

"Doubled-haploid" allohexaploid Brassica lines lose fertility and viability and accumulate genetic variation due to genomic instability

Microspore culture stimulates immature pollen grains to develop into plants via tissue culture and is used routinely in many crop species to produce "doubled haploids": homozygous, true-breeding lines. However, microspore culture is also often used on material that does not have stable meiosis, such as interspecific hybrids. In this case, the resulting progeny may lose their "doubled haploid" homozygous status as a result of chromosome missegregation and homoeologous exchanges. However, little is known about the frequency of these effects. We assessed fertility, meiosis and genetic variability in self-pollinated progeny sets (the MDL2 population) resulting from first-generation plants (the MDL1 population) derived from microspores of a near-allohexaploid interspecific hybrid from the cross (Brassica napus × B. carinata) × B. juncea. Allelic inheritance and copy number variation were predicted using single nucleotide polymorphism marker data from the Illumina Infinium 60K Brassica array. Seed fertility and viability decreased substantially from the MDL1 to the MDL2 generation. In the MDL2 population, 87% of individuals differed genetically from their MDL1 parent. These genetic differences resulted from novel homoeologous exchanges between chromosomes, chromosome loss and gain, and segregation and instability of pre-existing karyotype abnormalities. Novel karyotype change was extremely common, with 2.2 new variants observed per MDL2 individual. Significant differences between progeny sets in the number of novel genetic variants were also observed. Meiotic instability clearly has the potential to dramatically change karyotypes (often without detectable effects on the presence or absence of alleles) in putatively homozygous, microspore-derived lines, resulting in loss of fertility and viability.

Copy number variation of Ppd-B1 is the major determinant of heading time in durum wheat

Background

Heading time is an important adaptive trait in durum wheat. In hexaploid wheat, Photoperiod-1 (Ppd) loci are essential regulators of heading time, with Ppd-B1 conferring photoperiod insensitivity through copy number variations (CNV). In tetraploid wheat, the D-genome Ppd-D1 locus is absent and generally, our knowledge on the genetic architecture underlying heading time lacks behind that of bread wheat.

Results

In this study, we employed a panel of 328 diverse European durum genotypes that were evaluated for heading time at five environments. Genome-wide association mapping identified six putative QTL, with a major QTL on chromosome 2B explaining 26.2% of the genotypic variance. This QTL was shown to correspond to copy number variation at Ppd-B1, for which two copy number variants appear to be present. The higher copy number confers earlier heading and was more frequent in the heat and drought prone countries of lower latitude. In addition, two other QTL, corresponding to Vrn-B3 (TaFT) and Ppd-A1, were found to explain 9.5 and 5.3% of the genotypic variance, respectively.

Conclusions

Our results revealed the yet unknown role of copy number variation of Ppd-B1 as the major source underlying the variation in heading time in European durum wheat. The observed geographic patterns underline the adaptive value of this polymorphism and suggest that it is already used in durum breeding to tailor cultivars to specific target environments. In a broader context our findings provide further support for a more widespread role of copy number variation in mediating abiotic and biotic stress tolerance in plants.

Electronic supplementary material

The online version of this article (10.1186/s12863-019-0768-2) contains supplementary material, which is available to authorized users.

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.

A reference-grade wild soybean genome

Efficient crop improvement depends on the application of accurate genetic information contained in diverse germplasm resources. Here we report a reference-grade genome of wild soybean accession W05, with a final assembled genome size of 1013.2 Mb and a contig N50 of 3.3 Mb. The analytical power of the W05 genome is demonstrated by several examples. First, we identify an inversion at the locus determining seed coat color during domestication. Second, a translocation event between chromosomes 11 and 13 of some genotypes is shown to interfere with the assignment of QTLs. Third, we find a region containing copy number variations of the Kunitz trypsin inhibitor (KTI) genes. Such findings illustrate the power of this assembly in the analysis of large structural variations in soybean germplasm collections. The wild soybean genome assembly has wide applications in comparative genomic and evolutionary studies, as well as in crop breeding and improvement programs.

Modeling copy number variation in the genomic prediction of maize hybrids

Our study indicates that copy variants may play an essential role in the phenotypic variation of complex traits in maize hybrids. Moreover, predicting hybrid phenotypes by combining additive-dominance effects with copy variants has the potential to be a viable predictive model. Non-additive effects resulting from the actions of multiple loci may influence trait variation in single-cross hybrids. In addition, complementation of allelic variation could be a valuable contributor to hybrid genetic variation, especially when crossing inbred lines with higher contents of copy gains. With this in mind, we aimed (1) to study the association between copy number variation (CNV) and hybrid phenotype, and (2) to compare the predictive ability (PA) of additive and additive-dominance genomic best linear unbiased prediction model when combined with the effects of CNV in two datasets of maize hybrids (USP and HELIX). In the USP dataset, we observed a significant negative phenotypic correlation of low magnitude between copy number loss and plant height, revealing a tendency that more copy losses lead to lower plants. In the same set, when CNV was combined with the additive plus dominance effects, the PA significantly increased only for plant height under low nitrogen. In this case, CNV effects explicitly capture relatedness between individuals and add extra information to the model. In the HELIX dataset, we observed a pronounced difference in PA between additive (0.50) and additive-dominance (0.71) models for predicting grain yield, suggesting a significant contribution of dominance. We conclude that copy variants may play an essential role in the phenotypic variation of complex traits in maize hybrids, although the inclusion of CNVs into datasets does not return significant gains concerning PA.

Extreme copy number variation at a tRNA ligase gene affecting phenology and fitness in yellow monkeyflowers

Copy number variation (CNV) is a major part of the genetic diversity segregating within populations, but remains poorly understood relative to single nucleotide variation. Here, we report on a tRNA ligase gene (Migut.N02091; RLG1a) exhibiting unprecedented, and fitness-relevant, CNV within an annual population of the yellow monkeyflower Mimulus guttatus. RLG1a variation was associated with multiple traits in pooled population sequencing (PoolSeq) scans of phenotypic and phenological cohorts. Resequencing of inbred lines revealed intermediate-frequency three-copy variants of RLG1a (trip+; 5/35 = 14%), and trip+ lines exhibited elevated RLG1a expression under multiple conditions. trip+ carriers, in addition to being over-represented in late-flowering and large-flowered PoolSeq populations, flowered later under stressful conditions in a greenhouse experiment (p < 0.05). In wild population samples, we discovered an additional rare RLG1a variant (high+) that carries 250-300 copies of RLG1a totalling ~5.7 Mb (20-40% of a chromosome). In the progeny of a high+ carrier, Mendelian segregation of diagnostic alleles and qPCR-based copy counts indicate that high+ is a single tandem array unlinked to the single-copy RLG1a locus. In the wild, high+ carriers had highest fitness in two particularly dry and/or hot years (2015 and 2017; both p < 0.01), while single-copy individuals were twice as fecund as either CNV type in a lush year (2016: p < 0.005). Our results demonstrate fluctuating selection on CNVs affecting phenological traits in a wild population, suggest that plant tRNA ligases mediate stress-responsive life-history traits, and introduce a novel system for investigating the molecular mechanisms of gene amplification.

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.

Single-cell copy number variant detection reveals the dynamics and diversity of adaptation

Copy number variants (CNVs) are a pervasive source of genetic variation and evolutionary potential, but the dynamics and diversity of CNVs within evolving populations remain unclear. Long-term evolution experiments in chemostats provide an ideal system for studying the molecular processes underlying CNV formation and the temporal dynamics with which they are generated, selected, and maintained. Here, we developed a fluorescent CNV reporter to detect de novo gene amplifications and deletions in individual cells. We used the CNV reporter in Saccharomyces cerevisiae to study CNV formation at the GAP1 locus, which encodes the general amino acid permease, in different nutrient-limited chemostat conditions. We find that under strong selection, GAP1 CNVs are repeatedly generated and selected during the early stages of adaptive evolution, resulting in predictable dynamics. Molecular characterization of CNV-containing lineages shows that the CNV reporter detects different classes of CNVs, including aneuploidies, nonreciprocal translocations, tandem duplications, and complex CNVs. Despite GAP1's proximity to repeat sequences that facilitate intrachromosomal recombination, breakpoint analysis revealed that short inverted repeat sequences mediate formation of at least 50% of GAP1 CNVs. Inverted repeat sequences are also found at breakpoints at the DUR3 locus, where CNVs are selected in urea-limited chemostats. Analysis of 28 CNV breakpoints indicates that inverted repeats are typically 8 nucleotides in length and separated by 40 bases. The features of these CNVs are consistent with origin-dependent inverted-repeat amplification (ODIRA), suggesting that replication-based mechanisms of CNV formation may be a common source of gene amplification. We combined the CNV reporter with barcode lineage tracking and found that 102-104 independent CNV-containing lineages initially compete within populations, resulting in extreme clonal interference. However, only a small number (18-21) of CNV lineages ever constitute more than 1% of the CNV subpopulation, and as selection progresses, the diversity of CNV lineages declines. Our study introduces a novel means of studying CNVs in heterogeneous cell populations and provides insight into their dynamics, diversity, and formation mechanisms in the context of adaptive evolution.

Attention sports fans! The far-reaching contributions of bud sport mutants to horticulture and plant biology

A bud sport is a lateral shoot, inflorescence or single flower/fruit with a visibly different phenotype from the rest of the plant. The new phenotype is often caused by a stable somatic mutation in a single cell that is passed on to its clonal descendants and eventually populates part or all of a meristem. In many cases, a bud sport can be vegetatively propagated, thereby preserving the novel phenotype without sexual reproduction. Bud sports provide new characteristics while retaining the desirable qualities of the parent plant, which is why many bud sports have been developed into popular cultivars. We present an overview of the history of bud sports, the causes and methods of detecting somaclonal variation, and the types of mutant phenotypes that have arisen spontaneously. We focus on examples where the molecular or cytological changes causing the phenotype have been identified. Analysis of these sports has provided valuable insight into developmental processes, gene function and regulation, and in some cases has revealed new information about layer-specific roles of some genes. Examination of the molecular changes causing a phenotype and in some cases reversion back to the original state has contributed to our understanding of the mechanisms that drive genomic evolution.