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

A Strategy of Assessing Gene Copy Number Differentiation Between Populations Using Ultra-Fast De Novo Assembly of Next-Generation Sequencing Data

Gene duplication and loss play pivotal roles in the evolutionary dynamics of genomes, contributing to species phenotypic diversity and adaptation. However, detecting copy number variations (CNVs) in homoploid populations and newly-diverged species using short reads from next-generation sequencing (NGS) with traditional methods can often be challenging due to uneven read coverage caused by variations in GC content and the presence of repetitive sequences. To address these challenges, we developed a novel pipeline, ST4gCNV, which leverages ultra-fast de novo assemblies of NGS data to detect gene-specific CNVs between populations. The pipeline effectively reduces the variance of read coverage due to technical factors such as GC bias, providing a reliable CNV detection with a minimum sequencing depth of 10. We successfully apply ST4gCNV to the resequencing analysis of homoploid species Nelumbo nucifera and Nelumbo lutea (lotus). We reveal significant CNV-driven differentiation between these species, particularly in genes related to petal colour diversity such as those involved in the anthocyanin pathway. By highlighting the extensive gene duplication and loss events in Nelumbo, our study demonstrates the utility of ST4gCNV in population genomics and underscores its potential of integrating genomic CNV analysis with traditional SNP-based resequencing analysis.

A comprehensive map of DNA-segment copy number variation in 491 genomes of common wheat uncovers genes associated with multiple agronomic traits

DNA-segment copy number variations (DSCNVs), such as deletions and duplications, are important sources of genomic structural variation. However, the types and sizes of DSCNVs, as well as their genome-wide distribution and potential functions, are poorly understood in wheat. Here, we identified 198 985 DSCNVs by investigating 491 genomes of common wheat and found that they account for 20% of the entire genome. Interestingly, approximately 38% of genes are linked to DSCNVs. The number of DSCNVs within each accession ranges from 47 366 to 96 342, and their total sizes vary from 421.3 to 1267.9 Mb. We found that 957 and 1304 DSCNVs have been favored by breeders in China and the United States, respectively. By conducting DSCNV-based genome-wide association studies for the principal components of plant developmental and yield-component traits, we identified 34 loci as directly or indirectly involved in controlling the formation of multiple traits. Notably, a newly discovered DSCNV covering TaFT-D1 is significantly associated with flowering time and other agronomic traits. Overall, our findings highlight the potential of DSCNVs to drive fundamental discoveries in plant science. The comprehensive DSCNV map and the DSCNV-associated genes will also facilitate future research efforts to improve wheat yield, quality, and adaptation.

Induced and natural variation affect traits independently in hybrid Populus

The genetic control of many plant traits can be highly complex. Both allelic variation (sequence change) and dosage variation (copy number change) contribute to a plant's phenotype. While numerous studies have investigated the effect of allelic or dosage variation, very few have documented both within the same system, leaving their relative contribution to phenotypic effects unclear. The Populus genome is highly polymorphic, and poplars are fairly tolerant of gene dosage variation. Here, using a previously established Populus hybrid F1 population, we assessed and compared the effect of natural allelic variation and induced dosage variation on biomass, phenology, and leaf morphology traits. We identified QTLs for many of these traits, but our results indicate limited overlap between the QTLs associated with natural allelic variation and induced dosage variation. Additionally, the integration of data from both allelic and dosage variation identifies a larger set of QTLs that together explain a larger percentage of the phenotypic variance. Finally, our results suggest that the effect of the large indels might mask that of allelic QTLs. Our study helps clarify the relationship between allelic and dosage variation and their effects on quantitative traits.

De novo genome assembly of white clover (Trifolium repens L.) reveals the role of copy number variation in rapid environmental adaptation

BACKGROUND

White clover (Trifolium repens) is a globally important perennial forage legume. This species also serves as an eco-evolutionary model system for studying within-species chemical defense variation; it features a well-studied polymorphism for cyanogenesis (HCN release following tissue damage), with higher frequencies of cyanogenic plants favored in warmer locations worldwide. Using a newly generated haplotype-resolved genome and two other long-read assemblies, we tested the hypothesis that copy number variants (CNVs) at cyanogenesis genes play a role in the ability of white clover to rapidly adapt to local environments. We also examined questions on subgenome evolution in this recently evolved allotetraploid species and on chromosomal rearrangements in the broader IRLC legume clade.

RESULTS

Integration of PacBio HiFi, Omni-C, Illumina, and linkage map data yielded a completely de novo genome assembly for white clover (created without a priori sequence assignment to subgenomes). We find that white clover has undergone extensive transposon diversification since its origin but otherwise shows highly conserved genome organization and composition with its diploid progenitors. Unlike some other clover species, its chromosomal structure is conserved with other IRLC legumes. We further find extensive evidence of CNVs at the major cyanogenesis loci; these contribute to quantitative variation in the cyanogenic phenotype and to local adaptation across wild North American populations.

CONCLUSIONS

This work provides a case study documenting the role of CNVs in local adaptation in a plant species, and it highlights the value of pan-genome data for identifying contributions of structural variants to adaptation in nature.

Population genomics and epigenomics of Spirodela polyrhiza provide insights into the evolution of facultative asexuality

Many plants are facultatively asexual, balancing short-term benefits with long-term costs of asexuality. During range expansion, natural selection likely influences the genetic controls of asexuality in these organisms. However, evidence of natural selection driving asexuality is limited, and the evolutionary consequences of asexuality on the genomic and epigenomic diversity remain controversial. We analyzed population genomes and epigenomes of Spirodela polyrhiza, (L.) Schleid., a facultatively asexual plant that flowers rarely, revealing remarkably low genomic diversity and DNA methylation levels. Within species, demographic history and the frequency of asexual reproduction jointly determined intra-specific variations of genomic diversity and DNA methylation levels. Genome-wide scans revealed that genes associated with stress adaptations, flowering and embryogenesis were under positive selection. These data are consistent with the hypothesize that natural selection can shape the evolution of asexuality during habitat expansions, which alters genomic and epigenomic diversity levels.

Integration of eQTL and GWAS analysis uncovers a genetic regulation of natural ionomic variation in Arabidopsis

KEY MESSAGE

This study provided important insights into the genetic architecture of variations in A. thaliana leaf ionome in a cell-type-specific manner. The functional interpretation of traits associated variants by expression quantitative trait loci (eQTL) analysis is usually performed in bulk tissue samples. While the regulation of gene expression is context-dependent, such as cell-type-specific manner. In this study, we estimated cell-type abundances from 728 bulk tissue samples using single-cell RNA-sequencing dataset, and performed cis-eQTL mapping to identify cell-type-interaction eQTL (cis-eQTLs(ci)) in A. thaliana. Also, we performed Genome-wide association studies (GWAS) analyses for 999 accessions to identify the genetic basis of variations in A. thaliana leaf ionome. As a result, a total of 5,664 unique eQTL genes and 15,038 unique cis-eQTLs(ci) were significant. The majority (62.83%) of cis-eQTLs(ci) were cell-type-specific eQTLs. Using colocalization, we uncovered one interested gene AT2G25590 in Phloem cell, encoding a kind of plant Tudor-like protein with possible chromatin-associated functions, which colocalized with the most significant cis-eQTL(ci) of a Mo-related locus (Chr2:10,908,806:A:C; P = 3.27 × 10-27). Furthermore, we prioritized eight target genes associated with AT2G25590, which were previously reported in regulating the concentration of Mo element in A. thaliana. This study revealed the genetic regulation of ionomic variations and provided a foundation for further studies on molecular mechanisms of genetic variants controlling the A. thaliana ionome.

Repurposed inhibitor of bacterial dihydrodipicolinate reductase exhibits effective herbicidal activity

Herbicide resistance represents one of the biggest threats to our natural environment and agricultural sector. Thus, new herbicides are urgently needed to tackle the rise in herbicide-resistant weeds. Here, we employed a novel strategy to repurpose a 'failed' antibiotic into a new and target-specific herbicidal compound. Specifically, we identified an inhibitor of bacterial dihydrodipicolinate reductase (DHDPR), an enzyme involved in lysine biosynthesis in plants and bacteria, that exhibited no antibacterial activity but severely attenuated germination of the plant Arabidopsis thaliana. We confirmed that the inhibitor targets plant DHDPR orthologues in vitro, and exhibits no toxic effects against human cell lines. A series of analogues were then synthesised with improved efficacy in germination assays and against soil-grown A. thaliana. We also showed that our lead compound is the first lysine biosynthesis inhibitor with activity against both monocotyledonous and dicotyledonous weed species, by demonstrating its effectiveness at reducing the germination and growth of Lolium rigidum (rigid ryegrass) and Raphanus raphanistrum (wild radish). These results provide proof-of-concept that DHDPR inhibition may represent a much-needed new herbicide mode of action. Furthermore, this study exemplifies the untapped potential of repurposing 'failed' antibiotic scaffolds to fast-track the development of herbicide candidates targeting the respective plant enzymes.

A comprehensive analysis of copy number variations in diverse apple populations

BACKGROUND

As an important source of genetic variation, copy number variation (CNV) can alter the dosage of DNA segments, which in turn may affect gene expression level and phenotype. However, our knowledge of CNV in apple is still limited. Here, we obtained high-confidence CNVs and investigated their functional impact based on genome resequencing data of two apple populations, cultivars and wild relatives.

RESULTS

In this study, we identified 914,610 CNVs comprising 14,839 CNV regions (CNVRs) from 346 apple accessions, including 289 cultivars and 57 wild relatives. CNVRs summed to 71.19 Mb, accounting for 10.03% of the apple genome. Under the low linkage disequilibrium (LD) with nearby SNPs, they could also accurately reflect the population structure of apple independent of SNPs. Furthermore, A total of 3,621 genes were covered by CNVRs and functionally involved in biological processes such as defense response, reproduction and metabolic processes. In addition, the population differentiation index ([Formula: see text]) analysis between cultivars and wild relatives revealed 127 CN-differentiated genes, which may contribute to trait differences in these two populations.

CONCLUSIONS

This study was based on identification of CNVs from 346 diverse apple accessions, which to our knowledge was the largest dataset for CNV analysis in apple. Our work presented the first comprehensive CNV map and provided valuable resources for understanding genomic variations in apple.

Differences in the intraspecies copy number variation of Arabidopsis thaliana conserved and nonconserved miRNA genes

MicroRNAs (miRNAs) regulate gene expression by RNA interference mechanism. In plants, miRNA genes (MIRs) which are grouped into conserved families, i.e. they are present among the different plant taxa, are involved in the regulation of many developmental and physiological processes. The roles of the nonconserved MIRs-which are MIRs restricted to one plant family, genus, or even species-are less recognized; however, many of them participate in the responses to biotic and abiotic stresses. Both over- and underproduction of miRNAs may influence various biological processes. Consequently, maintaining intracellular miRNA homeostasis seems to be crucial for the organism. Deletions and duplications in the genomic sequence may alter gene dosage and/or activity. We evaluated the extent of copy number variations (CNVs) among Arabidopsis thaliana (Arabidopsis) MIRs in over 1000 natural accessions, using population-based analysis of the short-read sequencing data. We showed that the conserved MIRs were unlikely to display CNVs and their deletions were extremely rare, whereas nonconserved MIRs presented moderate variation. Transposon-derived MIRs displayed exceptionally high diversity. Conversely, MIRs involved in the epigenetic control of transposons reactivated during development were mostly invariable. MIR overlap with the protein-coding genes also limited their variability. At the expression level, a higher rate of nonvariable, nonconserved miRNAs was detectable in Col-0 leaves, inflorescence, and siliques compared to nonconserved variable miRNAs, although the expression of both groups was much lower than that of the conserved MIRs. Our data indicate that CNV rate of Arabidopsis MIRs is related with their age, function, and genomic localization.

Extensive sequence duplication in Arabidopsis revealed by pseudo-heterozygosity

BACKGROUND

It is apparent that genomes harbor much structural variation that is largely undetected for technical reasons. Such variation can cause artifacts when short-read sequencing data are mapped to a reference genome. Spurious SNPs may result from mapping of reads to unrecognized duplicated regions. Calling SNP using the raw reads of the 1001 Arabidopsis Genomes Project we identified 3.3 million (44%) heterozygous SNPs. Given that Arabidopsis thaliana (A. thaliana) is highly selfing, and that extensively heterozygous individuals have been removed, we hypothesize that these SNPs reflected cryptic copy number variation.

RESULTS

The heterozygosity we observe consists of particular SNPs being heterozygous across individuals in a manner that strongly suggests it reflects shared segregating duplications rather than random tracts of residual heterozygosity due to occasional outcrossing. Focusing on such pseudo-heterozygosity in annotated genes, we use genome-wide association to map the position of the duplicates. We identify 2500 putatively duplicated genes and validate them using de novo genome assemblies from six lines. Specific examples included an annotated gene and nearby transposon that transpose together. We also demonstrate that cryptic structural variation produces highly inaccurate estimates of DNA methylation polymorphism.

CONCLUSIONS

Our study confirms that most heterozygous SNP calls in A. thaliana are artifacts and suggest that great caution is needed when analyzing SNP data from short-read sequencing. The finding that 10% of annotated genes exhibit copy-number variation, and the realization that neither gene- nor transposon-annotation necessarily tells us what is actually mobile in the genome suggests that future analyses based on independently assembled genomes will be very informative.

Rab GTPases, tethers, and SNAREs work together to regulate Arabidopsis cell plate formation

Cell plates are transient structures formed by the fusion of vesicles at the center of the dividing plane; furthermore, these are precursors to new cell walls and are essential for cytokinesis. Cell plate formation requires a highly coordinated process of cytoskeletal rearrangement, vesicle accumulation and fusion, and membrane maturation. Tethering factors have been shown to interact with the Ras superfamily of small GTP binding proteins (Rab GTPases) and soluble N-ethylmaleimide-sensitive factor attachment protein receptors (SNAREs), which are essential for cell plate formation during cytokinesis and are fundamental for maintaining normal plant growth and development. In Arabidopsis thaliana, members of the Rab GTPases, tethers, and SNAREs are localized in cell plates, and mutations in the genes encoding these proteins result in typical cytokinesis-defective phenotypes, such as the formation of abnormal cell plates, multinucleated cells, and incomplete cell walls. This review highlights recent findings on vesicle trafficking during cell plate formation mediated by Rab GTPases, tethers, and SNAREs.

Analysis of Arabidopsis non-reference accessions reveals high diversity of metabolic gene clusters and discovers new candidate cluster members

Metabolic gene clusters (MGCs) are groups of genes involved in a common biosynthetic pathway. They are frequently formed in dynamic chromosomal regions, which may lead to intraspecies variation and cause phenotypic diversity. We examined copy number variations (CNVs) in four Arabidopsis thaliana MGCs in over one thousand accessions with experimental and bioinformatic approaches. Tirucalladienol and marneral gene clusters showed little variation, and the latter was fixed in the population. Thalianol and especially arabidiol/baruol gene clusters displayed substantial diversity. The compact version of the thalianol gene cluster was predominant and more conserved than the noncontiguous version. In the arabidiol/baruol cluster, we found a large genomic insertion containing divergent duplicates of the CYP705A2 and BARS1 genes. The BARS1 paralog, which we named BARS2, encoded a novel oxidosqualene synthase. The expression of the entire arabidiol/baruol gene cluster was altered in the accessions with the duplication. Moreover, they presented different root growth dynamics and were associated with warmer climates compared to the reference-like accessions. In the entire genome, paired genes encoding terpene synthases and cytochrome P450 oxidases were more variable than their nonpaired counterparts. Our study highlights the role of dynamically evolving MGCs in plant adaptation and phenotypic diversity.

Copy Number Variation among Resistance Genes Analogues in Brassica napus

Copy number variations (CNVs) are defined as deletions, duplications and insertions among individuals of a species. There is growing evidence that CNV is a major factor underlining various autoimmune disorders and diseases in humans; however, in plants, especially oilseed crops, the role of CNVs in disease resistance is not well studied. Here, we investigate the genome-wide diversity and genetic properties of CNVs in resistance gene analogues (RGAs) across eight Brassica napus lines. A total of 1137 CNV events (704 deletions and 433 duplications) were detected across 563 RGAs. The results show CNVs are more likely to occur across clustered RGAs compared to singletons. In addition, 112 RGAs were linked to a blackleg resistance QTL, of which 25 were affected by CNV. Overall, we show that the presence and abundance of CNVs differ between lines, suggesting that in B. napus, the distribution of CNVs depends on genetic background. Our findings advance the understanding of CNV as an important type of genomic structural variation in B. napus and provide a resource to support breeding of advanced canola lines.

An efficient CRISPR-Cas9 enrichment sequencing strategy for characterizing complex and highly duplicated genomic regions. A case study in the Prunus salicina LG3-MYB10 genes cluster

BACKGROUND

Genome complexity is largely linked to diversification and crop innovation. Examples of regions with duplicated genes with relevant roles in agricultural traits are found in many crops. In both duplicated and non-duplicated genes, much of the variability in agronomic traits is caused by large as well as small and middle scale structural variants (SVs), which highlights the relevance of the identification and characterization of complex variability between genomes for plant breeding.

RESULTS

Here we improve and demonstrate the use of CRISPR-Cas9 enrichment combined with long-read sequencing technology to resolve the MYB10 region in the linkage group 3 (LG3) of Japanese plum (Prunus salicina). This region, which has a length from 90 to 271 kb according to the P. salicina genomes available, is associated with fruit color variability in Prunus species. We demonstrate the high complexity of this region, with homology levels between Japanese plum varieties comparable to those between Prunus species. We cleaved MYB10 genes in five plum varieties using the Cas9 enzyme guided by a pool of crRNAs. The barcoded fragments were then pooled and sequenced in a single MinION Oxford Nanopore Technologies (ONT) run, yielding 194 Mb of sequence. The enrichment was confirmed by aligning the long reads to the plum reference genomes, with a mean read on-target value of 4.5% and a depth per sample of 11.9x. From the alignment, 3261 SNPs and 287 SVs were called and phased. A de novo assembly was constructed for each variety, which also allowed detection, at the haplotype level, of the variability in this region.

CONCLUSIONS

CRISPR-Cas9 enrichment is a versatile and powerful tool for long-read targeted sequencing even on highly duplicated and/or polymorphic genomic regions, being especially useful when a reference genome is not available. Potential uses of this methodology as well as its limitations are further discussed.

A composite strategy of genome-wide association study and copy number variation analysis for carcass traits in a Duroc pig population

BACKGROUND

Carcass traits are important in pig breeding programs for improving pork production. Understanding the genetic variants underlies complex phenotypes can help explain trait variation in pigs. In this study, we integrated a weighted single-step genome-wide association study (wssGWAS) and copy number variation (CNV) analyses to map genetic variations and genes associated with loin muscle area (LMA), loin muscle depth (LMD) and lean meat percentage (LMP) in Duroc pigs.

RESULTS

Firstly, we performed a genome-wide analysis for CNV detection using GeneSeek Porcine SNP50 Bead chip data of 3770 pigs. A total of 11,100 CNVs were detected, which were aggregated by overlapping 695 CNV regions (CNVRs). Next, we investigated CNVs of pigs from the same population by whole-genome resequencing. A genome-wide analysis of 21 pigs revealed 23,856 CNVRs that were further divided into three categories (851 gain, 22,279 loss, and 726 mixed), which covered 190.8 Mb (~ 8.42%) of the pig autosomal genome. Further, the identified CNVRs were used to determine an overall validation rate of 68.5% for the CNV detection accuracy of chip data. CNVR association analyses identified one CNVR associated with LMA, one with LMD and eight with LMP after applying stringent Bonferroni correction. The wssGWAS identified eight, six and five regions explaining more than 1% of the additive genetic variance for LMA, LMD and LMP, respectively. The CNVR analyses and wssGWAS identified five common regions, of which three regions were associated with LMA and two with LMP. Four genes (DOK7, ARAP1, ELMO2 and SLC13A3) were highlighted as promising candidates according to their function.

CONCLUSIONS

We determined an overall validation rate for the CNV detection accuracy of low-density chip data and constructed a genomic CNV map for Duroc pigs using resequencing, thereby proving a value genetic variation resource for pig genome research. Furthermore, our study utilized a composite genetic strategy for complex traits in pigs, which will contribute to the study for elucidating the genetic architecture that may be influenced and regulated by multiple forms of variations.

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.

ATP-binding cassette transporters in nonmodel plants

Knowledge about plant ATP-binding cassette (ABC) proteins is of great value for sustainable agriculture, economic yield, and the generation of high-quality products, especially under unfavorable growth conditions. We have learned much about ABC proteins in model organisms, notably Arabidopsis thaliana; however, the importance of research dedicated to these transporters extends far beyond Arabidopsis biology. Recent progress in genomic and transcriptomic approaches for nonmodel and noncanonical model plants allows us to look at ABC transporters from a wider perspective and consider chemodiversity and functionally driven adaptation as distinctive mechanisms during their evolution. Here, by considering several representatives from agriculturally important families and recent progress in functional characterization of nonArabidopsis ABC proteins, we aim to bring attention to understanding the evolutionary background, distribution among lineages and possible mechanisms underlying the adaptation of this versatile transport system for plant needs. Increasing the knowledge of ABC proteins in nonmodel plants will facilitate breeding and development of new varieties based on, for example, genetic variations of endogenous genes and/or genome editing, representing an alternative to transgenic approaches.

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.

Diverse phosphate and auxin transport loci distinguish phosphate tolerant from sensitive Arabidopsis accessions

Phosphorus (P) is an essential element for plant growth often limiting agroecosystems. To identify genetic determinants of performance under variable phosphate (Pi) supply, we conducted genome-wide association studies on five highly predictive Pi starvation response traits in 200 Arabidopsis (Arabidopsis thaliana) accessions. Pi concentration in Pi-limited organs had the strongest, and primary root length had the weakest genetic component. Of 70 trait-associated candidate genes, 17 responded to Pi withdrawal. The PHOSPHATE TRANSPORTER1 gene cluster on chromosome 5 comprises PHT1;1, PHT1;2, and PHT1;3 with known impact on P status. A second locus featured uncharacterized endomembrane-associated auxin efflux carrier encoding PIN-LIKES7 (PILS7) which was more strongly suppressed in Pi-limited roots of Pi-starvation sensitive accessions. In the Col-0 background, Pi uptake and organ growth were impaired in both Pi-limited pht1;1 and two pils7 T-DNA insertion mutants, while Pi -limited pht1;2 had higher biomass and pht1;3 was indistinguishable from wild-type. Copy number variation at the PHT1 locus with loss of the PHT1;3 gene and smaller scale deletions in PHT1;1 and PHT1;2 predicted to alter both protein structure and function suggest diversification of PHT1 is a key driver for adaptation to P limitation. Haplogroup analysis revealed a phosphorylation site in the protein encoded by the PILS7 allele from stress-sensitive accessions as well as additional auxin-responsive elements in the promoter of the "stress tolerant" allele. The former allele's inability to complement the pils7-1 mutant in the Col-0 background implies the presence of a kinase signaling loop controlling PILS7 activity in accessions from P-rich environments, while survival in P-poor environments requires fine-tuning of stress-responsive root auxin signaling.

Calling large indels in 1047 Arabidopsis with IndelEnsembler

Large indels greatly impact the observable phenotypes in different organisms including plants and human. Hence, extracting large indels with high precision and sensitivity is important. Here, we developed IndelEnsembler to detect large indels in 1047 Arabidopsis whole-genome sequencing data. IndelEnsembler identified 34 093 deletions, 12 913 tandem duplications and 9773 insertions. Our large indel dataset was more comprehensive and accurate compared with the previous dataset of AthCNV (1). We captured nearly twice of the ground truth deletions and on average 27% more ground truth duplications compared with AthCNV, though our dataset has less number of large indels compared with AthCNV. Our large indels were positively correlated with transposon elements across the Arabidopsis genome. The non-homologous recombination events were the major formation mechanism of deletions in Arabidopsis genome. The Neighbor joining (NJ) tree constructed based on IndelEnsembler's deletions clearly divided the geographic subgroups of 1047 Arabidopsis. More importantly, our large indels represent a previously unassessed source of genetic variation. Approximately 49% of the deletions have low linkage disequilibrium (LD) with surrounding single nucleotide polymorphisms. Some of them could affect trait performance. For instance, using deletion-based genome-wide association study (DEL-GWAS), the accessions containing a 182-bp deletion in AT1G11520 had delayed flowering time and all accessions in north Sweden had the 182-bp deletion. We also found the accessions with 65-bp deletion in the first exon of AT4G00650 (FRI) flowered earlier than those without it. These two deletions cannot be detected in AthCNV and, interestingly, they do not co-occur in any Arabidopsis thaliana accession. By SNP-GWAS, surrounding SNPs of these two deletions do not correlate with flowering time. This example demonstrated that existing large indel datasets miss phenotypic variations and our large indel dataset filled in the gap.

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.

GWAS for genetics of complex quantitative traits: Genome to pangenome and SNPs to SVs and k-mers

The development of improved methods for genome-wide association studies (GWAS) for genetics of quantitative traits has been an active area of research during the last 25 years. This activity initially started with the use of mixed linear model (MLM), which was variously modified. During the last decade, however, with the availability of high throughput next generation sequencing (NGS) technology, development and use of pangenomes and novel markers including structural variations (SVs) and k-mers for GWAS has taken over as a new thrust area of research. Pangenomes and SVs are now available in humans, livestock, and a number of plant species, so that these resources along with k-mers are being used in GWAS for exploring additional genetic variation that was hitherto not available for analysis. These developments have resulted in significant improvement in GWAS methodology for detection of marker-trait associations (MTAs) that are relevant to human healthcare and crop improvement.

Multi-omics approach highlights differences between RLP classes in Arabidopsis thaliana

BACKGROUND

The Leucine rich-repeat (LRR) receptor-like protein (RLP) family is a complex gene family with 57 members in Arabidopsis thaliana. Some members of the RLP family are known to be involved in basal developmental processes, whereas others are involved in defence responses. However, functional data is currently only available for a small subset of RLPs, leaving the remaining ones classified as RLPs of unknown function.

RESULTS

Using publicly available datasets, we annotated RLPs of unknown function as either likely defence-related or likely fulfilling a more basal function in plants. Then, using these categories, we can identify important characteristics that differ between the RLP subclasses. We found that the two classes differ in abundance on both transcriptome and proteome level, physical clustering in the genome and putative interaction partners. However, the classes do not differ in the genetic di versity of their individual members in accessible pan-genome data.

CONCLUSIONS

Our work has several implications for work related to functional studies on RLPs as well as for the understanding of RLP gene family evolution. Using our annotations, we can make suggestions on which RLPs can be identified as potential immune receptors using genetics tools and thereby complement disease studies. The lack of differences in nucleotide diversity between the two RLP subclasses further suggests that non-synonymous diversity of gene sequences alone cannot distinguish defence from developmental genes. By contrast, differences in transcript and protein abundance or clustering at genomic loci might also allow for functional annotations and characterisation in other plant species.

A systems genetics approach to deciphering the effect of dosage variation on leaf morphology in Populus

Gene copy number variation is frequent in plant genomes of various species, but the impact of such gene dosage variation on morphological traits is poorly understood. We used a large population of Populus carrying genomically characterized insertions and deletions across the genome to systematically assay the effect of gene dosage variation on a suite of leaf morphology traits. A systems genetics approach was used to integrate insertion and deletion locations, leaf morphology phenotypes, gene expression, and transcriptional network data, to provide an overview of how gene dosage influences morphology. Dosage-sensitive genomic regions were identified that influenced individual or pleiotropic morphological traits. We also identified cis-expression quantitative trait loci (QTL) within these dosage QTL regions, a subset of which modulated trans-expression QTL as well. Integration of data types within a gene co-expression framework identified co-expressed gene modules that are dosage sensitive, enriched for dosage expression QTL, and associated with morphological traits. Functional description of these modules linked dosage-sensitive morphological variation to specific cellular processes, as well as candidate regulatory genes. Together, these results show that gene dosage variation can influence morphological variation through complex changes in gene expression, and suggest that frequently occurring gene dosage variation has the potential to likewise influence quantitative traits in nature.

Construction of a chromosome-level genome and variation map for the Pacific oyster Crassostrea gigas

The Pacific oyster (Crassostrea gigas) is a widely distributed marine bivalve of great ecological and economic importance. In this study, we provide a high-quality chromosome-level genome assembled using Pacific Bioscience long reads and Hi-C-based and linkage-map-based scaffolding technologies and a high-resolution variation map constructed using large-scale resequencing analysis. The 586.8 Mb genome consists of 10 pseudochromosome sequences ranging from 38.6 to 78.9 Mb, containing 301 contigs with an N50 size of 3.1 Mb. A total of 30,078 protein-coding genes were predicted, of which 22,757 (75.7%) were high-reliability annotations supported by a homologous match to a curated protein in the SWISS-PROT database or transcript expression. Although a medium level of repeat components (57.2%) was detected, the genomic content of the segmental duplications reached 26.2%, which is the highest among the reported genomes. By whole genome resequencing analysis of 495 Pacific oysters, a comprehensive variation map was built, comprised of 4.78 million single nucleotide polymorphisms, 0.60 million short insertions and deletions, and 49,333 copy number variation regions. The structural variations can lead to an average interindividual genomic divergence of 0.21, indicating their crucial role in shaping the Pacific oyster genome diversity. The large amount of mosaic distributed repeat elements, small variations, and copy number variations indicate that the Pacific oyster is a diploid organism with an extremely high genomic complexity at the intra- and interindividual level. The genome and variation maps can improve our understanding of oyster genome diversity and enrich the resources for oyster molecular evolution, comparative genomics, and genetic research.

A New Catalog of Structural Variants in 1,301 A. thaliana Lines from Africa, Eurasia, and North America Reveals a Signature of Balancing Selection at Defense Response Genes

Genomic variation in the model plant Arabidopsis thaliana has been extensively used to understand evolutionary processes in natural populations, mainly focusing on single-nucleotide polymorphisms. Conversely, structural variation has been largely ignored in spite of its potential to dramatically affect phenotype. Here, we identify 155,440 indels and structural variants ranging in size from 1 bp to 10 kb, including presence/absence variants (PAVs), inversions, and tandem duplications in 1,301 A. thaliana natural accessions from Morocco, Madeira, Europe, Asia, and North America. We show evidence for strong purifying selection on PAVs in genes, in particular for housekeeping genes and homeobox genes, and we find that PAVs are concentrated in defense-related genes (R-genes, secondary metabolites) and F-box genes. This implies the presence of a "core" genome underlying basic cellular processes and a "flexible" genome that includes genes that may be important in spatially or temporally varying selection. Further, we find an excess of intermediate frequency PAVs in defense response genes in nearly all populations studied, consistent with a history of balancing selection on this class of genes. Finally, we find that PAVs in genes involved in the cold requirement for flowering (vernalization) and drought response are strongly associated with temperature at the sites of origin.

AtMAD: Arabidopsis thaliana multi-omics association database

Integration analysis of multi-omics data provides a comprehensive landscape for understanding biological systems and mechanisms. The abundance of high-quality multi-omics data (genomics, transcriptomics, methylomics and phenomics) for the model organism Arabidopsis thaliana enables scientists to study the genetic mechanism of many biological processes. However, no resource is available to provide comprehensive and systematic multi-omics associations for Arabidopsis. Here, we developed an Arabidopsis thaliana Multi-omics Association Database (AtMAD, http://www.megabionet.org/atmad), a public repository for large-scale measurements of associations between genome, transcriptome, methylome, pathway and phenotype in Arabidopsis, designed for facilitating identification of eQTL, emQTL, Pathway-mQTL, Phenotype-pathway, GWAS, TWAS and EWAS. Candidate variants/methylations/genes were identified in AtMAD for specific phenotypes or biological processes, many of them are supported by experimental evidence. Based on the multi-omics association strategy, we have identified 11 796 cis-eQTLs and 10 119 trans-eQTLs. Among them, 68 837 environment-eQTL associations and 149 622 GWAS-eQTL associations were identified and stored in AtMAD. For expression–methylation quantitative trait loci (emQTL), we identified 265 776 emQTLs and 122 344 pathway-mQTLs. For TWAS and EWAS, we obtained 62 754 significant phenotype-gene associations and 3 993 379 significant phenotype-methylation associations, respectively. Overall, the multi-omics associated network in AtMAD will provide new insights into exploring biological mechanisms of plants at multi-omics levels.