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

Analysis of Extreme Phenotype Bulk Copy Number Variation (XP-CNV) Identified the Association of rp1 with Resistance to Goss's Wilt of Maize

Goss's wilt (GW) of maize is caused by the Gram-positive bacterium Clavibacter michiganensis subsp. nebraskensis (Cmn) and has spread in recent years throughout the Great Plains, posing a threat to production. The genetic basis of plant resistance is unknown. Here, a simple method for quantifying disease symptoms was developed and used to select cohorts of highly resistant and highly susceptible lines known as extreme phenotypes (XP). Copy number variation (CNV) analyses using whole genome sequences of bulked XP revealed 141 genes containing CNV between the two XP groups. The CNV genes include the previously identified common rust resistant locus rp1. Multiple Rp1 accessions with distinct rp1 haplotypes in an otherwise susceptible accession exhibited hypersensitive responses upon inoculation. GW provides an excellent system for the genetic dissection of diseases caused by closely related subspecies of C. michiganesis. Further work will facilitate breeding strategies to control GW and provide needed insight into the resistance mechanism of important related diseases such as bacterial canker of tomato and bacterial ring rot of potato.

Sequence analysis of European maize inbred line F2 provides new insights into molecular and chromosomal characteristics of presence/absence variants

BACKGROUND

Maize is well known for its exceptional structural diversity, including copy number variants (CNVs) and presence/absence variants (PAVs), and there is growing evidence for the role of structural variation in maize adaptation. While PAVs have been described in this important crop species, they have been only scarcely characterized at the sequence level and the extent of presence/absence variation and relative chromosomal landscape of inbred-specific regions remain to be elucidated.

RESULTS

De novo genome sequencing of the French F2 maize inbred line revealed 10,044 novel genomic regions larger than 1 kb, making up 88 Mb of DNA, that are present in F2 but not in B73 (PAV). This set of maize PAV sequences allowed us to annotate PAV content and to analyze sequence breakpoints. Using PAV genotyping on a collection of 25 temperate lines, we also analyzed Linkage Disequilibrium in PAVs and flanking regions, and PAV frequencies within maize genetic groups.

CONCLUSIONS

We highlight the possible role of MMEJ-type double strand break repair in maize PAV formation and discover 395 new genes with transcriptional support. Pattern of linkage disequilibrium within PAVs strikingly differs from this of flanking regions and is in accordance with the intuition that PAVs may recombine less than other genomic regions. We show that most PAVs are ancient, while some are found only in European Flint material, thus pinpointing structural features that may be at the origin of adaptive traits involved in the success of this material. Characterization of such PAVs will provide useful material for further association genetic studies in European and temperate maize.

Copy number variation and disease resistance in plants

Plant genome diversity varies from single nucleotide polymorphisms to large-scale deletions, insertions, duplications, or re-arrangements. These re-arrangements of sequences resulting from duplication, gains or losses of DNA segments are termed copy number variations (CNVs). During the last decade, numerous studies have emphasized the importance of CNVs as a factor affecting human phenotype; in particular, CNVs have been associated with risks for several severe diseases. In plants, the exploration of the extent and role of CNVs in resistance against pathogens and pests is just beginning. Since CNVs are likely to be associated with disease resistance in plants, an understanding of the distribution of CNVs could assist in the identification of novel plant disease-resistance genes. In this paper, we review existing information about CNVs; their importance, role and function, as well as their association with disease resistance in plants.

Current Status and Challenges in Identifying Disease Resistance Genes in Brassica napus

Brassica napus is an economically important crop across different continents including temperate and subtropical regions in Europe, Canada, South Asia, China and Australia. Its widespread cultivation also brings setbacks as it plays host to fungal, oomycete and chytrid pathogens that can lead to serious yield loss. For sustainable crop production, identification of resistance (R) genes in B. napus has become of critical importance. In this review, we discuss four key pathogens affecting Brassica crops: Clubroot (Plasmodiophora brassicae), Blackleg (Leptosphaeria maculans and L. biglobosa), Sclerotinia Stem Rot (Sclerotinia sclerotiorum), and Downy Mildew (Hyaloperonospora parasitica). We first review current studies covering prevalence of these pathogens on Brassica crops and highlight the R genes and QTL that have been identified from Brassica species against these pathogens. Insights into the relationships between the pathogen and its Brassica host, the unique host resistance mechanisms and how these affect resistance outcomes is also presented. We discuss challenges in identification and deployment of R genes in B. napus in relation to highly specific genetic interactions between host subpopulations and pathogen pathotypes and emphasize the need for common or shared techniques and research materials or tighter collaboration between researchers to reconcile the inconsistencies in the research outcomes. Using current genomics tools, we provide examples of how characterization and cloning of R genes in B. napus can be carried out more effectively. Lastly, we put forward strategies to breed resistant cultivars through introgressions supported by genomic approaches and suggest prospects that can be implemented in the future for a better, pathogen-resistant B. napus.

Transposable Elements Direct The Coevolution between Plants and Microbes

Transposable elements are powerful drivers of genome evolution in many eukaryotes. Although they are mostly considered as 'selfish' genetic elements, increasing evidence suggests that they contribute to genetic variability; particularly under stress conditions. Over the past few years, the role of transposable elements during host-microbe interactions has been recognised. It has been proposed that many pathogenic microbes have evolved a 'two-speed' genome with regions that show increased variability and that are enriched in transposable elements and pathogenicity-related genes. Plants similarly display structured genomes with transposable-element-rich regions that mediate accelerated evolution. Immune receptor genes typically reside in such regions. Various mechanisms have recently been identified through which transposable elements contribute to the coevolution between plants and their associated microbes.

Crop Breeding Chips and Genotyping Platforms: Progress, Challenges, and Perspectives

There is a rapidly rising trend in the development and application of molecular marker assays for gene mapping and discovery in field crops and trees. Thus far, more than 50 SNP arrays and 15 different types of genotyping-by-sequencing (GBS) platforms have been developed in over 25 crop species and perennial trees. However, much less effort has been made on developing ultra-high-throughput and cost-effective genotyping platforms for applied breeding programs. In this review, we discuss the scientific bottlenecks in existing SNP arrays and GBS technologies and the strategies to develop targeted platforms for crop molecular breeding. We propose that future practical breeding platforms should adopt automated genotyping technologies, either array or sequencing based, target functional polymorphisms underpinning economic traits, and provide desirable prediction accuracy for quantitative traits, with universal applications under wide genetic backgrounds in crops. The development of such platforms faces serious challenges at both the technological level due to cost ineffectiveness, and the knowledge level due to large genotype-phenotype gaps in crop plants. It is expected that such genotyping platforms will be achieved in the next ten years in major crops in consideration of (a) rapid development in gene discovery of important traits, (b) deepened understanding of quantitative traits through new analytical models and population designs, (c) integration of multi-layer -omics data leading to identification of genes and pathways responsible for important breeding traits, and (d) improvement in cost effectiveness of large-scale genotyping. Crop breeding chips and genotyping platforms will provide unprecedented opportunities to accelerate the development of cultivars with desired yield potential, quality, and enhanced adaptation to mitigate the effects of climate change.

Investigation of orthologous pathogen recognition gene-rich regions in solanaceous species

Pathogen receptor proteins such as receptor-like protein (RLP), receptor-like kinase (RLK), and nucleotide-binding leucine-rich repeat (NLR) play a leading role in plant immunity activation. The genome architecture of such genes has been extensively investigated in several plant species. However, we still know little about their elaborate reorganization that arose during the plant speciation process. Using recently released pepper and eggplant genome sequences, we were able to identify 1097 pathogen recognition genes (PRGs) in the cultivated pepper Zunla-1 and 775 in the eggplant line Nakate-Shinkuro. The retrieved genes were analysed for their tendency to cluster, using different methods to infer the means of grouping. Orthologous relationships among clustering loci were found, and interesting reshuffling within given loci was observed for each analysed species. The information obtained was integrated into a comparative map to highlight the evolutionary dynamics in which the PRG loci were involved. Diversification of 14 selected PRG-rich regions was also explored using a DNA target-enrichment approach. A large number of gene variants were found as well as rearrangements of sequences encoding single protein domain and changes in chromosome gene order among species. Gene duplication and transposition activity have clearly influenced plant genome R-gene architecture and diversification. Our findings contribute to addressing several biological questions concerning the parallel evolution that occurred between genomes of the family Solanaceae. Moreover, the integration of different methods proved a powerful approach to reconstruct the evolutionary history in plant families and to transfer important biology findings among plant genomes.

Fitness effects of mutation: testing genetic redundancy in Arabidopsis thaliana

Screens of organisms with disruptive mutations in a single gene often fail to detect phenotypic consequences for the majority of mutants. One explanation for this phenomenon is that the presence of paralogous loci provides genetic redundancy. However, it is also possible that the assayed traits are affected by few loci, that effects could be subtle or that phenotypic effects are restricted to certain environments. We assayed a set of T-DNA insertion mutant lines of Arabidopsis thaliana to determine the frequency with which mutation affected fitness-related phenotypes. We found that between 8% and 42% of the assayed lines had altered fitness from the wild type. Furthermore, many of these lines exhibited fitness greater than the wild type. In a second experiment, we grew a subset of the lines in multiple environments and found whether a T-DNA insert increased or decreased fitness traits depended on the assay environment. Overall, our evidence contradicts the hypothesis that genetic redundancy is a common phenomenon in A. thaliana for fitness traits. Evidence for redundancy from prior screens of knockout mutants may often be an artefact of the design of the phenotypic assays which have focused on less complex phenotypes than fitness and have used single environments. Finally, our study adds to evidence that beneficial mutations may represent a significant component of the mutational spectrum of A. thaliana.

An Induced Chromosomal Translocation in Soybean Disrupts a KASI Ortholog and Is Associated with a High-Sucrose and Low-Oil Seed Phenotype

Mutagenesis is a useful tool in many crop species to induce heritable genetic variability for trait improvement and gene discovery. In this study, forward screening of a soybean fast neutron (FN) mutant population identified an individual that produced seed with nearly twice the amount of sucrose (8.1% on dry matter basis) and less than half the amount of oil (8.5% on dry matter basis) as compared to wild type. Bulked segregant analysis (BSA), comparative genomic hybridization, and genome resequencing were used to associate the seed composition phenotype with a reciprocal translocation between chromosomes 8 and 13. In a backcross population, the translocation perfectly cosegregated with the seed composition phenotype and exhibited non-Mendelian segregation patterns. We hypothesize that the translocation is responsible for the altered seed composition by disrupting a β-ketoacyl-[acyl carrier protein] synthase 1 (KASI) ortholog. KASI is a core fatty acid synthesis enzyme that is involved in the conversion of sucrose into oil in developing seeds. This finding may lead to new research directions for developing soybean cultivars with modified carbohydrate and oil seed composition.

Soybean domestication: the origin, genetic architecture and molecular bases

Domestication provides an important model for the study of evolution, and information learned from domestication research aids in the continued improvement of crop species. Recent progress in de novo assembly and whole-genome resequencing of wild and cultivated soybean genomes, in addition to new archeological discoveries, sheds light on the origin of this important crop and provides a clearer view on the modes of artificial selection that drove soybean domestication and diversification. This novel genomic information enables the search for polymorphisms that underlie variation in agronomic traits and highlights genes that exhibit a signature of selection, leading to the identification of a number of candidate genes that may have played important roles in soybean domestication, diversification and improvement. These discoveries provide a novel point of comparison on the evolutionary bases of important agronomic traits among different crop species.

MLPA-Based Analysis of Copy Number Variation in Plant Populations

Copy number variants (CNVs) are intraspecies duplications/deletions of large DNA segments (>1 kb). A growing number of reports highlight the functional and evolutionary impact of CNV in plants, increasing the need for appropriate tools that enable locus-specific CNV genotyping on a population scale. Multiplex ligation-dependent probe amplification (MLPA) is considered a gold standard in genotyping CNV in humans. Consequently, numerous commercial MLPA assays for CNV-related human diseases have been created. We routinely genotype complex multiallelic CNVs in human and plant genomes using the modified MLPA procedure based on fully synthesized oligonucleotide probes (90-200 nt), which greatly simplifies the design process and allows for the development of custom assays. Here, we present a step-by-step protocol for gene-specific MLPA probe design, multiplexed assay setup and data analysis in a copy number genotyping experiment in plants. As a case study, we present the results of a custom assay designed to genotype the copy number status of 12 protein coding genes in a population of 80 Arabidopsis accessions. The genes were pre-selected based on whole genome sequencing data and are localized in the genomic regions that display different levels of population-scale variation (non-variable, biallelic, or multiallelic, as well as CNVs overlapping whole genes or their fragments). The presented approach is suitable for population-scale validation of the CNV regions inferred from whole genome sequencing data analysis and for focused analysis of selected genes of interest. It can also be very easily adopted for any plant species, following optimization of the template amount and design of the appropriate control probes, according to the general guidelines presented in this paper.

CNVs into the wild: screening the genomes of conifer trees (Picea spp.) reveals fewer gene copy number variations in hybrids and links to adaptation

Background

Copy number variations (CNVs) have been linked to different phenotypes in human, including many diseases. A genome-scale understanding of CNVs is available in a few plants but none are wild species, leaving a knowledge gap regarding their genome biology and evolutionary role. We developed a reliable CNV detection method for species lacking contiguous reference genome. We selected multiple probes within 14,078 gene sequences and developed comparative genome hybridization on arrays. Gene CNVs were assessed in three full-sib families from species with 20 Gb genomes, i.e., white and black spruce, and interior spruce - a natural hybrid.

Results

We discovered hundreds of gene CNVs in each species, 3612 in total, which were enriched in functions related to stress and defense responses and narrow expression profiles, indicating a potential role in adaptation. The number of shared CNVs was in accordance with the degree of relatedness between individuals and species. The genetically mapped subset of these genes showed a wide distribution across the genome, implying numerous structural variations. The hybrid family presented significantly fewer CNVs, suggesting that the admixture of two species within one genome reduces the occurrence of CNVs.

Conclusions

The approach we developed is of particular interest in non-model species lacking a reference genome. Our findings point to a role for CNVs in adaptation. Their reduced abundance in the hybrid may limit genetic variability and evolvability of hybrids. We propose that CNVs make a qualitatively distinct contribution to adaptation which could be important for short term change.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3458-8) contains supplementary material, which is available to authorized users.

Analysis of Plant Pan-Genomes and Transcriptomes with GET_HOMOLOGUES-EST, a Clustering Solution for Sequences of the Same Species

The pan-genome of a species is defined as the union of all the genes and non-coding sequences found in all its individuals. However, constructing a pan-genome for plants with large genomes is daunting both in sequencing cost and the scale of the required computational analysis. A more affordable alternative is to focus on the genic repertoire by using transcriptomic data. Here, the software GET_HOMOLOGUES-EST was benchmarked with genomic and RNA-seq data of 19 Arabidopsis thaliana ecotypes and then applied to the analysis of transcripts from 16 Hordeum vulgare genotypes. The goal was to sample their pan-genomes and classify sequences as core, if detected in all accessions, or accessory, when absent in some of them. The resulting sequence clusters were used to simulate pan-genome growth, and to compile Average Nucleotide Identity matrices that summarize intra-species variation. Although transcripts were found to under-estimate pan-genome size by at least 10%, we concluded that clusters of expressed sequences can recapitulate phylogeny and reproduce two properties observed in A. thaliana gene models: accessory loci show lower expression and higher non-synonymous substitution rates than core genes. Finally, accessory sequences were observed to preferentially encode transposon components in both species, plus disease resistance genes in cultivated barleys, and a variety of protein domains from other families that appear frequently associated with presence/absence variation in the literature. These results demonstrate that pan-genome analyses are useful to explore germplasm diversity.

The pangenome of an agronomically important crop plant Brassica oleracea

There is an increasing awareness that as a result of structural variation, a reference sequence representing a genome of a single individual is unable to capture all of the gene repertoire found in the species. A large number of genes affected by presence/absence and copy number variation suggest that it may contribute to phenotypic and agronomic trait diversity. Here we show by analysis of the Brassica oleracea pangenome that nearly 20% of genes are affected by presence/absence variation. Several genes displaying presence/absence variation are annotated with functions related to major agronomic traits, including disease resistance, flowering time, glucosinolate metabolism and vitamin biosynthesis.

Brassica oleracea is a single species that includes diverse crops such as cabbage, broccoli and Brussels sprouts. Here, the authors identify genes not captured in existing B. oleracea reference genomes by the assembly of a pangenome and show variations in gene content that may be related to important agronomic traits

Arabidopsis thaliana population analysis reveals high plasticity of the genomic region spanning MSH2, AT3G18530 and AT3G18535 genes and provides evidence for NAHR-driven recurrent CNV events occurring in this location

Background

Intraspecies copy number variations (CNVs), defined as unbalanced structural variations of specific genomic loci, ≥1 kb in size, are present in the genomes of animals and plants. A growing number of examples indicate that CNVs may have functional significance and contribute to phenotypic diversity. In the model plant Arabidopsis thaliana at least several hundred protein-coding genes might display CNV; however, locus-specific genotyping studies in this plant have not been conducted.

Results

We analyzed the natural CNVs in the region overlapping MSH2 gene that encodes the DNA mismatch repair protein, and AT3G18530 and AT3G18535 genes that encode poorly characterized proteins. By applying multiplex ligation-dependent probe amplification and droplet digital PCR we genotyped those genes in 189 A. thaliana accessions. We found that AT3G18530 and AT3G18535 were duplicated (2–14 times) in 20 and deleted in 101 accessions. MSH2 was duplicated in 12 accessions (up to 12-14 copies) but never deleted. In all but one case, the MSH2 duplications were associated with those of AT3G18530 and AT3G18535. Considering the structure of the CNVs, we distinguished 5 genotypes for this region, determined their frequency and geographical distribution. We defined the CNV breakpoints in 35 accessions with AT3G18530 and AT3G18535 deletions and tandem duplications and showed that they were reciprocal events, resulting from non-allelic homologous recombination between 99 %-identical sequences flanking these genes. The widespread geographical distribution of the deletions supported by the SNP and linkage disequilibrium analyses of the genomic sequence confirmed the recurrent nature of this CNV.

Conclusions

We characterized in detail for the first time the complex multiallelic CNV in Arabidopsis genome. The region encoding MSH2, AT3G18530 and AT3G18535 genes shows enormous variation of copy numbers among natural ecotypes, being a remarkable example of high Arabidopsis genome plasticity. We provided the molecular insight into the mechanism underlying the recurrent nature of AT3G18530-AT3G18535 duplications/deletions. We also performed the first direct comparison of the two leading experimental methods, suitable for assessing the DNA copy number status. Our comprehensive case study provides foundation information for further analyses of CNV evolution in Arabidopsis and other plants, and their possible use in plant breeding.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3221-1) contains supplementary material, which is available to authorized users.

Inter-varietal structural variation in grapevine genomes

Grapevine (Vitis vinifera L.) is one of the world's most important crop plants, which is of large economic value for fruit and wine production. There is much interest in identifying genomic variations and their functional effects on inter-varietal, phenotypic differences. Using an approach developed for the analysis of human and mammalian genomes, which combines high-throughput sequencing, array comparative genomic hybridization, fluorescent in situ hybridization and quantitative PCR, we created an inter-varietal atlas of structural variations and single nucleotide variants (SNVs) for the grapevine genome analyzing four economically and genetically relevant table grapevine varieties. We found 4.8 million SNVs and detected 8% of the grapevine genome to be affected by genomic variations. We identified more than 700 copy number variation (CNV) regions and more than 2000 genes subjected to CNV as potential candidates for phenotypic differences between varieties.

Copy number variation at the HvCBF4-HvCBF2 genomic segment is a major component of frost resistance in barley

A family of CBF transcription factors plays a major role in reconfiguring the plant transcriptome in response to low-freezing temperature in temperate cereals. In barley, more than 13 HvCBF genes map coincident with the major QTL FR-H2 suggesting them as candidates to explain the function of the locus. Variation in copy number (CNV) of specific HvCBFs was assayed in a panel of 41 barley genotypes using RT-qPCR. Taking advantage of an accurate phenotyping that combined Fv/Fm and field survival, resistance-associated variants within FR-H2 were identified. Genotypes with an increased copy number of HvCBF4 and HvCBF2 (at least ten and eight copies, respectively) showed greater frost resistance. A CAPS marker able to distinguish the CBF2A, CBF2B and CBF2A/B forms was developed and showed that all the higher-ranking genotypes in term of resistance harbour only CBF2A, while other resistant winter genotypes harbour also CBF2B, although at a lower CNV. In addition to the major involvement of the HvCBF4-HvCBF2 genomic segment in the proximal cluster of CBF elements, a negative role of HvCBF3 in the distal cluster was identified. Multiple linear regression models taking into account allelic variation at FR-H1/VRN-H1 explained 0.434 and 0.550 (both at p < 0.001) of the phenotypic variation for Fv/Fm and field survival respectively, while no interaction effect between CNV at the HvCBFs and FR-H1/VRN-H1 was found. Altogether our data suggest a major involvement of the CBF genes located in the proximal cluster, with no apparent involvement of the central cluster contrary to what was reported for wheat.

Fast neutron-induced structural rearrangements at a soybean NAP1 locus result in gnarled trichomes

KEY MESSAGE

Three adjacent and distinct sequence rearrangements were identified at a NAP1 locus in a soybean mutant. Genetic dissection and validation revealed the function of this gene in soybean trichome development. A soybean (Glycine max (L.) Merr.) gnarled trichome mutant, exhibiting stunted trichomes compared to wild-type, was identified in a fast neutron mutant population. Genetic mapping using whole genome sequencing-based bulked segregant analysis identified a 26.6 megabase interval on chromosome 20 that co-segregated with the phenotype. Comparative genomic hybridization analysis of the mutant indicated that the chromosome 20 interval included a small structural variant within the coding region of a soybean ortholog (Glyma.20G019300) of Arabidopsis Nck-Associated Protein 1 (NAP1), a regulator of actin nucleation during trichome morphogenesis. Sequence analysis of the candidate allele revealed multiple rearrangements within the coding region, including two deletions (approximately 1-2 kb each), a translocation, and an inversion. Further analyses revealed that the mutant allele perfectly co-segregated with the phenotype, and a wild-type soybean NAP1 transgene functionally complemented an Arabidopsis nap1 mutant. In addition, mapping and exon sequencing of NAP1 in a spontaneous soybean gnarled trichome mutant (T31) identified a frame shift mutation resulting in a truncation of the coding region. These data indicate that the soybean NAP1 gene is essential for proper trichome development and show the utility of the soybean fast neutron population for forward genetic approaches for identifying genes.

Characterization of the Poplar Pan-Genome by Genome-Wide Identification of Structural Variation

Many recent studies have emphasized the important role of structural variation (SV) in determining human genetic and phenotypic variation. In plants, studies aimed at elucidating the extent of SV are still in their infancy. Evidence has indicated a high presence and an active role of SV in driving plant genome evolution in different plant species.

With the aim of characterizing the size and the composition of the poplar pan-genome, we performed a genome-wide analysis of structural variation in three intercrossable poplar species: Populus nigra, Populus deltoides, and Populus trichocarpa. We detected a total of 7,889 deletions and 10,586 insertions relative to the P. trichocarpa reference genome, covering respectively 33.2 Mb and 62.9 Mb of genomic sequence, and 3,230 genes affected by copy number variation (CNV). The majority of the detected variants are inter-specific in agreement with a recent origin following separation of species.

Insertions and deletions (INDELs) were preferentially located in low-gene density regions of the poplar genome and were, for the majority, associated with the activity of transposable elements. Genes affected by SV showed lower-than-average expression levels and higher levels of dN/dS, suggesting that they are subject to relaxed selective pressure or correspond to pseudogenes.

Functional annotation of genes affected by INDELs showed over-representation of categories associated with transposable elements activity, while genes affected by genic CNVs showed enrichment in categories related to resistance to stress and pathogens. This study provides a genome-wide catalogue of SV and the first insight on functional and structural properties of the poplar pan-genome.

Using Genotyping by Sequencing to Map Two Novel Anthracnose Resistance Loci in Sorghum bicolor

Colletotrichum sublineola is an aggressive fungal pathogen that causes anthracnose in sorghum [Sorghum bicolor (L.) Moench]. The obvious symptoms of anthracnose are leaf blight and stem rot. Sorghum, the fifth most widely grown cereal crop in the world, can be highly susceptible to the disease, most notably in hot and humid environments. In the southeastern United States the acreage of sorghum has been increasing steadily in recent years, spurred by growing interest in producing biofuels, bio-based products, and animal feed. Resistance to anthracnose is, therefore, of paramount importance for successful sorghum production in this region. To identify anthracnose resistance loci present in the highly resistant cultivar ‘Bk7’, a biparental mapping population of F_3:4 and F_4:5 sorghum lines was generated by crossing ‘Bk7’ with the susceptible inbred ‘Early Hegari-Sart’. Lines were phenotyped in three environments and in two different years following natural infection. The population was genotyped by sequencing. Following a stringent custom filtering protocol, totals of 5186 and 2759 informative SNP markers were identified in the two populations. Segregation data and association analysis identified resistance loci on chromosomes 7 and 9, with the resistance alleles derived from ‘Bk7’. Both loci contain multiple classes of defense-related genes based on sequence similarity and gene ontologies. Genetic analysis following an independent selection experiment of lines derived from a cross between ‘Bk7’ and sweet sorghum ‘Mer81-4’ narrowed the resistance locus on chromosome 9 substantially, validating this QTL. As observed in other species, sorghum appears to have regions of clustered resistance genes. Further characterization of these regions will facilitate the development of novel germplasm with resistance to anthracnose and other diseases.

Candidate Gene Sequence Analyses toward Identifying -Type Resistance to

is one of three genetic loci conferring strain-specific resistance to (SMV). The locus has been mapped to a 154-kb region on chromosome 14, containing a cluster of five nucleotide-binding leucine-rich repeat (NB-LRR) resistance genes. High sequence similarity between the candidate genes challenges fine mapping of the locus. Among the five, Glyma14g38533 showed the highest transcript abundance in 1 to 3 h of SMV-G7 inoculation. Comparative sequence analyses were conducted with the five candidate NB-LRR genes from susceptible (-type) soybean [ (L.) Merr.] cultivar Williams 82, resistant (-type) cultivar Hwangkeum, and resistant lines L29 and RRR. Sequence comparisons revealed that Glyma14g38533 had far more polymorphisms than the other candidate genes. Interestingly, Glyma14g38533 gene from -type lines exhibited 150 single-nucleotide polymorphism (SNP and six insertion-deletion (InDel) markers relative to -type line, Furthermore, the polymorphisms identified in three -type lines were highly conserved. Several polymorphisms were validated in 18 -type resistant and six -type susceptible lines and were found associated with their disease response. The majority of the polymorphisms were located in LRR domain encoding region, which is involved in pathogen recognition via protein-protein interactions. These findings associating Glyma14g38533 with -type resistance to SMV suggest it is the most likely candidate gene for .

Genomic variation and DNA repair associated with soybean transgenesis: a comparison to cultivars and mutagenized plants

BACKGROUND

The safety of mutagenized and genetically transformed plants remains a subject of scrutiny. Data gathered and communicated on the phenotypic and molecular variation induced by gene transfer technologies will provide a scientific-based means to rationally address such concerns. In this study, genomic structural variation (e.g. large deletions and duplications) and single nucleotide polymorphism rates were assessed among a sample of soybean cultivars, fast neutron-derived mutants, and five genetically transformed plants developed through Agrobacterium based transformation methods.

RESULTS

On average, the number of genes affected by structural variations in transgenic plants was one order of magnitude less than that of fast neutron mutants and two orders of magnitude less than the rates observed between cultivars. Structural variants in transgenic plants, while rare, occurred adjacent to the transgenes, and at unlinked loci on different chromosomes. DNA repair junctions at both transgenic and unlinked sites were consistent with sequence microhomology across breakpoints. The single nucleotide substitution rates were modest in both fast neutron and transformed plants, exhibiting fewer than 100 substitutions genome-wide, while inter-cultivar comparisons identified over one-million single nucleotide polymorphisms.

CONCLUSIONS

Overall, these patterns provide a fresh perspective on the genomic variation associated with high-energy induced mutagenesis and genetically transformed plants. The genetic transformation process infrequently results in novel genetic variation and these rare events are analogous to genetic variants occurring spontaneously, already present in the existing germplasm, or induced through other types of mutagenesis. It remains unclear how broadly these results can be applied to other crops or transformation methods.

Genes and Small RNA Transcripts Exhibit Dosage-Dependent Expression Pattern in Maize Copy-Number Alterations

Copy-number alterations are widespread in animal and plant genomes, but their immediate impact on gene expression is still unclear. In animals, copy-number alterations usually exhibit dosage effects, except for sex chromosomes which tend to be dosage compensated. In plants, genes within small duplications (<100 kb) often exhibit dosage-dependent expression, whereas large duplications (>50 Mb) are more often dosage compensated. However, little or nothing is known about expression in moderately-sized (1-50 Mb) segmental duplications, and about the response of small RNAs to dosage change. Here, we compared maize (Zea mays) plants with two, three, and four doses of a 14.6-Mb segment of chromosome 1 that contains ∼300 genes. Plants containing the duplicated segment exhibit dosage-dependent effects on ear length and flowering time. Transcriptome analyses using GeneChip and RNA-sequencing methods indicate that most expressed genes and unique small RNAs within the duplicated segments exhibit dosage-dependent transcript levels. We conclude that dosage effect is the predominant regulatory response for both genes and unique small RNA transcripts in the segmental dosage series we tested. To our knowledge this is the first analysis of small RNA expression in plant gene dosage variants. Because segmental duplications comprise a significant proportion of eukaryotic genomes, these findings provide important new insight into the regulation of genes and small RNAs in response to dosage changes.

Genome-wide analysis of Dongxiang wild rice (Oryza rufipogon Griff.) to investigate lost/acquired genes during rice domestication

BACKGROUND

It is widely accepted that cultivated rice (Oryza sativa L.) was domesticated from common wild rice (Oryza rufipogon Griff.). Compared to other studies which concentrate on rice origin, this study is to genetically elucidate the substantially phenotypic and physiological changes from wild rice to cultivated rice at the whole genome level.

RESULTS

Instead of comparing two assembled genomes, this study directly compared the Dongxiang wild rice (DXWR) Illumina sequencing reads with the Nipponbare (O. sativa) complete genome without assembly of the DXWR genome. Based on the results from the comparative genomics analysis, structural variations (SVs) between DXWR and Nipponbare were determined to locate deleted genes which could have been acquired by Nipponbare during rice domestication. To overcome the limit of the SV detection, the DXWR transcriptome was also sequenced and compared with the Nipponbare transcriptome to discover the genes which could have been lost in DXWR during domestication. Both 1591 Nipponbare-acquired genes and 206 DXWR-lost transcripts were further analyzed using annotations from multiple sources. The NGS data are available in the NCBI SRA database with ID SRP070627.

CONCLUSIONS

These results help better understanding the domestication from wild rice to cultivated rice at the whole genome level and provide a genomic data resource for rice genetic research or breeding. One finding confirmed transposable elements contribute greatly to the genome evolution from wild rice to cultivated rice. Another finding suggested the photophosphorylation and oxidative phosphorylation system in cultivated rice could have adapted to environmental changes simultaneously during domestication.

Towards plant pangenomics

As an increasing number of genome sequences become available for a wide range of species, there is a growing understanding that the genome of a single individual is insufficient to represent the gene diversity within a whole species. Many studies examine the sequence diversity within genes, and this allelic variation is an important source of phenotypic variation which can be selected for by man or nature. However, the significant gene presence/absence variation that has been observed within species and the impact of this variation on traits is only now being studied in detail. The sum of the genes for a species is termed the pangenome, and the determination and characterization of the pangenome is a requirement to understand variation within a species. In this review, we explore the current progress in pangenomics as well as methods and approaches for the characterization of pangenomes for a wide range of plant species.

Exploring structural variants in environmentally sensitive gene families

Environmentally sensitive plant gene families like NBS-LRRs, receptor kinases, defensins and others, are known to be highly variable. However, most existing strategies for discovering and describing structural variation in complex gene families provide incomplete and imperfect results. The move to de novo genome assemblies for multiple accessions or individuals within a species is enabling more comprehensive and accurate insights about gene family variation. Earlier array-based genome hybridization and sequence-based read mapping methods were limited by their reliance on a reference genome and by misplacement of paralogous sequences. Variant discovery based on de novo genome assemblies overcome the problems arising from a reference genome and reduce sequence misplacement. As de novo genome sequencing moves to the use of longer reads, artifacts will be minimized, intact tandem gene clusters will be constructed accurately, and insights into rapid evolution will become feasible.

The impact and origin of copy number variations in the Oryza species

Background

Copy number variation (CNV), a complex genomic rearrangement, has been extensively studied in humans and other organisms. In plants, CNVs of several genes were found to be responsible for various important traits; however, the cause and consequence of CNVs remains largely unknown. Recently released next-generation sequencing (NGS) data provide an opportunity for a genome-wide study of CNVs in rice.

Results

Here, by an NGS-based approach, we generated a CNV map comprising 9,196 deletions compared to the reference genome ‘Nipponbare’. Using Oryza glaberrima as the outgroup, 80 % of the CNV events turned out to be insertions in Nipponbare. There were 2,806 annotated genes affected by these CNV events. We experimentally validated 28 functional CNV genes including OsMADS56, BPH14, OsDCL2b and OsMADS30, implying that CNVs might have contributed to phenotypic variations in rice. Most CNV genes were found to be located in non-co-linear positions by comparison to O. glaberrima. One of the origins of these non-co-linear genes was genomic duplications caused by transposon activity or double-strand break repair. Comprehensive analysis of mutation mechanisms suggested an abundance of CNVs formed by non-homologous end-joining and mobile element insertion.

Conclusions

This study showed the impact and origin of copy number variations in rice on a genomic scale.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2589-2) contains supplementary material, which is available to authorized users.

Genome-wide detection of intervals of genetic heterogeneity associated with complex traits

Motivation: Genetic heterogeneity, the fact that several sequence variants give rise to the same phenotype, is a phenomenon that is of the utmost interest in the analysis of complex phenotypes. Current approaches for finding regions in the genome that exhibit genetic heterogeneity suffer from at least one of two shortcomings: (i) they require the definition of an exact interval in the genome that is to be tested for genetic heterogeneity, potentially missing intervals of high relevance, or (ii) they suffer from an enormous multiple hypothesis testing problem due to the large number of potential candidate intervals being tested, which results in either many false positives or a lack of power to detect true intervals.

Results: Here, we present an approach that overcomes both problems: it allows one to automatically find all contiguous sequences of single nucleotide polymorphisms in the genome that are jointly associated with the phenotype. It also solves both the inherent computational efficiency problem and the statistical problem of multiple hypothesis testing, which are both caused by the huge number of candidate intervals. We demonstrate on Arabidopsis thaliana genome-wide association study data that our approach can discover regions that exhibit genetic heterogeneity and would be missed by single-locus mapping.

Conclusions: Our novel approach can contribute to the genome-wide discovery of intervals that are involved in the genetic heterogeneity underlying complex phenotypes.

Availability and implementation: The code can be obtained at: http://www.bsse.ethz.ch/mlcb/research/bioinformatics-and-computational-biology/sis.html.

Contact:felipe.llinares@bsse.ethz.ch

Supplementary information:Supplementary data are available at Bioinformatics online.

Diversity in boron toxicity tolerance of Australian barley (Hordeum vulgare L.) genotypes

BACKGROUND

Boron (B) is an important micronutrient for plant growth, but is toxic when levels are too high. This commonly occurs in environments with alkaline soils and relatively low rainfall, including many of the cereal growing regions of southern Australia. Four major genetic loci controlling tolerance to high soil B have been identified in the landrace barley, Sahara 3771. Genes underlying two of the loci encode the B transporters HvBot1 and HvNIP2;1.

RESULTS

We investigated sequence and expression level diversity in HvBot1 and HvNIP2;1 across barley germplasm, and identified five novel coding sequence alleles for HvBot1. Lines were identified containing either single or multiple copies of the Sahara HvBot1 allele. We established that only the tandemly duplicated Sahara allele conferred B tolerance, and this duplicated allele was found only in a set of nine lines accessioned in Australian collections as Sahara 3763-3771. HvNIP2;1 coding sequences were highly conserved across barley germplasm. We identified the likely causative SNP in the 5'UTR of Sahara HvNIP2;1, and propose that the creation of a small upstream open reading frame interferes with HvNIP2;1 translation in Sahara 3771. Similar to HvBot1, the tolerant HvNIP2;1 allele was unique to the Sahara barley accessions. We identified a new source of the 2H B tolerance allele controlling leaf symptom development, in the landrace Ethiopia 756.

CONCLUSIONS

Ethiopia 756, as well as the cultivar Sloop Vic which carries both the 2H and HvBot1 B tolerance alleles derived from Sahara 3771, may be valuable as alternative parents in breeding programs targeted to high soil B environments. There is significant diversity in B toxicity tolerance among contemporary Australian barley varieties but this is not related to variation at any of the four known B tolerance loci, indicating that novel, as yet undiscovered, sources of tolerance exist.

Transposon Insertions, Structural Variations, and SNPs Contribute to the Evolution of the Melon Genome

The availability of extensive databases of crop genome sequences should allow analysis of crop variability at an unprecedented scale, which should have an important impact in plant breeding. However, up to now the analysis of genetic variability at the whole-genome scale has been mainly restricted to single nucleotide polymorphisms (SNPs). This is a strong limitation as structural variation (SV) and transposon insertion polymorphisms are frequent in plant species and have had an important mutational role in crop domestication and breeding. Here, we present the first comprehensive analysis of melon genetic diversity, which includes a detailed analysis of SNPs, SV, and transposon insertion polymorphisms. The variability found among seven melon varieties representing the species diversity and including wild accessions and highly breed lines, is relatively high due in part to the marked divergence of some lineages. The diversity is distributed nonuniformly across the genome, being lower at the extremes of the chromosomes and higher in the pericentromeric regions, which is compatible with the effect of purifying selection and recombination forces over functional regions. Additionally, this variability is greatly reduced among elite varieties, probably due to selection during breeding. We have found some chromosomal regions showing a high differentiation of the elite varieties versus the rest, which could be considered as strongly selected candidate regions. Our data also suggest that transposons and SV may be at the origin of an important fraction of the variability in melon, which highlights the importance of analyzing all types of genetic variability to understand crop genome evolution.

Copy number variation at the GL7 locus contributes to grain size diversity in rice

Copy number variants (CNVs) are associated with changes in gene expression levels and contribute to various adaptive traits. Here we show that a CNV at the Grain Length on Chromosome 7 (GL7) locus contributes to grain size diversity in rice (Oryza sativa L.). GL7 encodes a protein homologous to Arabidopsis thaliana LONGIFOLIA proteins, which regulate longitudinal cell elongation. Tandem duplication of a 17.1-kb segment at the GL7 locus leads to upregulation of GL7 and downregulation of its nearby negative regulator, resulting in an increase in grain length and improvement of grain appearance quality. Sequence analysis indicates that allelic variants of GL7 and its negative regulator are associated with grain size diversity and that the CNV at the GL7 locus was selected for and used in breeding. Our work suggests that pyramiding beneficial alleles of GL7 and other yield- and quality-related genes may improve the breeding of elite rice varieties.

Gaining insight into soybean defense responses using functional genomics approaches

Soybean pathogens significantly impact yield, resulting in over $4 billion dollars in lost revenue annually in the United States. Despite the deployment of improved soybean cultivars, pathogens continue to evolve to evade plant defense responses. Thus, there is an urgent need to identify and characterize gene networks controlling defense responses to harmful pathogens. In this review, we focus on major advances that have been made in identifying the genes and gene networks regulating defense responses with an emphasis on soybean-pathogen interactions that have been amenable to gene function analyses using gene silencing technologies. Further we describe new research striving to identify genes involved in durable broad-spectrum resistance. Finally, we consider future prospects for functional genomic studies in soybean and demonstrate that understanding soybean disease and stress tolerance will be expedited at an unprecedented pace.

Genome-Wide Mapping of Structural Variations Reveals a Copy Number Variant That Determines Reproductive Morphology in Cucumber

Structural variations (SVs) represent a major source of genetic diversity. However, the functional impact and formation mechanisms of SVs in plant genomes remain largely unexplored. Here, we report a nucleotide-resolution SV map of cucumber (Cucumis sativas) that comprises 26,788 SVs based on deep resequencing of 115 diverse accessions. The largest proportion of cucumber SVs was formed through nonhomologous end-joining rearrangements, and the occurrence of SVs is closely associated with regions of high nucleotide diversity. These SVs affect the coding regions of 1676 genes, some of which are associated with cucumber domestication. Based on the map, we discovered a copy number variation (CNV) involving four genes that defines the Female (F) locus and gives rise to gynoecious cucumber plants, which bear only female flowers and set fruit at almost every node. The CNV arose from a recent 30.2-kb duplication at a meiotically unstable region, likely via microhomology-mediated break-induced replication. The SV set provides a snapshot of structural variations in plants and will serve as an important resource for exploring genes underlying key traits and for facilitating practical breeding in cucumber.

Homologues of potato chromosome 5 show variable collinearity in the euchromatin, but dramatic absence of sequence similarity in the pericentromeric heterochromatin

BACKGROUND

In flowering plants it has been shown that de novo genome assemblies of different species and genera show a significant drop in the proportion of alignable sequence. Within a plant species, however, it is assumed that different haplotypes of the same chromosome align well. In this paper we have compared three de novo assemblies of potato chromosome 5 and report on the sequence variation and the proportion of sequence that can be aligned.

RESULTS

For the diploid potato clone RH89-039-16 (RH) we produced two linkage phase controlled and haplotype-specific assemblies of chromosome 5 based on BAC-by-BAC sequencing, which were aligned to each other and compared to the 52 Mb chromosome 5 reference sequence of the doubled monoploid clone DM 1-3 516 R44 (DM). We identified 17.0 Mb of non-redundant sequence scaffolds derived from euchromatic regions of RH and 38.4 Mb from the pericentromeric heterochromatin. For 32.7 Mb of the RH sequences the correct position and order on chromosome 5 was determined, using genetic markers, fluorescence in situ hybridisation and alignment to the DM reference genome. This ordered fraction of the RH sequences is situated in the euchromatic arms and in the heterochromatin borders. In the euchromatic regions, the sequence collinearity between the three chromosomal homologs is good, but interruption of collinearity occurs at nine gene clusters. Towards and into the heterochromatin borders, absence of collinearity due to structural variation was more extensive and was caused by hemizygous and poorly aligning regions of up to 450 kb in length. In the most central heterochromatin, a total of 22.7 Mb sequence from both RH haplotypes remained unordered. These RH sequences have very few syntenic regions and represent a non-alignable region between the RH and DM heterochromatin haplotypes of chromosome 5.

CONCLUSIONS

Our results show that among homologous potato chromosomes large regions are present with dramatic loss of sequence collinearity. This stresses the need for more de novo reference assemblies in order to capture genome diversity in this crop. The discovery of three highly diverged pericentric heterochromatin haplotypes within one species is a novelty in plant genome analysis. The possible origin and cytogenetic implication of this heterochromatin haplotype diversity are discussed.

High-resolution genetic mapping of maize pan-genome sequence anchors

In addition to single-nucleotide polymorphisms, structural variation is abundant in many plant genomes. The structural variation across a species can be represented by a 'pan-genome', which is essential to fully understand the genetic control of phenotypes. However, the pan-genome's complexity hinders its accurate assembly via sequence alignment. Here we demonstrate an approach to facilitate pan-genome construction in maize. By performing 18 trillion association tests we map 26 million tags generated by reduced representation sequencing of 14,129 maize inbred lines. Using machine-learning models we select 4.4 million accurately mapped tags as sequence anchors, 1.1 million of which are presence/absence variations. Structural variations exhibit enriched association with phenotypic traits, indicating that it is a significant source of adaptive variation in maize. The ability to efficiently map ultrahigh-density pan-genome sequence anchors enables fine characterization of structural variation and will advance both genetic research and breeding in many crops.

PAV markers in Sorghum bicolour: genome pattern, affected genes and pathways, and genetic linkage map construction

KEY MESSAGE

5,511 genic small-size PAVs in sorghum were identified and examined, including the pattern and the function enrichment of PAV genes. 325 PAV markers were developed to construct a genetic map. Presence/absence variants (PAVs) correlate closely to the phenotypic variation, by impacting plant genome sizes and the adaption to the environment. To shed more light on their genome-wide patterns, functions and the possibility of using them as molecular markers, we generated next generation genome sequencing data for four sorghum inbred lines and used associated bioinformatic pipelines to identify small-size PAVs (40-10 kb). Five thousand five hundreds and eleven genic PAVs (40-10 kb) were identified and found to affect 3,238 genes. These PAVs were mainly distributed on the sub-telomeric regions, but the highest proportions occurred in the vicinity of the centromeric regions. One of the prominent features of the PAVs is the high occurrence of long terminal repeats retrotransposons and DNA transposons. PAVs caused various alterations to gene structure, primarily including the coding sequence variants, intron variants, transcript ablation, and initiator codon changes. The genes affected by PAVs were significantly enriched in those involved in stress responses and protein modification. We used 325 PAVs polymorphic between two sorghum inbred lines Ji2731 and E-Tian, together with 49 SSR markers, and constructed a genetic map, which consisted of 10 linkage groups corresponding to the 10 chromosomes of sorghum and spanned 1,430.3 cM in length covering 97% of the physical genome. The resources reported here should be useful for genetic study and breeding of sorghum and related species.

Evolution and selection of Rhg1, a copy-number variant nematode-resistance locus

The soybean cyst nematode (SCN) resistance locus Rhg1 is a tandem repeat of a 31.2 kb unit of the soybean genome. Each 31.2-kb unit contains four genes. One allele of Rhg1, Rhg1-b, is responsible for protecting most US soybean production from SCN. Whole-genome sequencing was performed, and PCR assays were developed to investigate allelic variation in sequence and copy number of the Rhg1 locus across a population of soybean germplasm accessions. Four distinct sequences of the 31.2-kb repeat unit were identified, and some Rhg1 alleles carry up to three different types of repeat unit. The total number of copies of the repeat varies from 1 to 10 per haploid genome. Both copy number and sequence of the repeat correlate with the resistance phenotype, and the Rhg1 locus shows strong signatures of selection. Significant linkage disequilibrium in the genome outside the boundaries of the repeat allowed the Rhg1 genotype to be inferred using high-density single nucleotide polymorphism genotyping of 15 996 accessions. Over 860 germplasm accessions were found likely to possess Rhg1 alleles. The regions surrounding the repeat show indications of non-neutral evolution and high genetic variability in populations from different geographic locations, but without evidence of fixation of the resistant genotype. A compelling explanation of these results is that balancing selection is in operation at Rhg1.

Genetic basis and detection of unintended effects in genetically modified crop plants

In January 2014, an international meeting sponsored by the International Life Sciences Institute/Health and Environmental Sciences Institute and the Canadian Food Inspection Agency titled “Genetic Basis of Unintended Effects in Modified Plants” was held in Ottawa, Canada, bringing together over 75 scientists from academia, government, and the agro-biotech industry. The objectives of the meeting were to explore current knowledge and identify areas requiring further study on unintended effects in plants and to discuss how this information can inform and improve genetically modified (GM) crop risk assessments. The meeting featured presentations on the molecular basis of plant genome variability in general, unintended changes at the molecular and phenotypic levels, and the development and use of hypothesis-driven evaluations of unintended effects in assessing conventional and GM crops. The development and role of emerging “omics” technologies in the assessment of unintended effects was also discussed. Several themes recurred in a number of talks; for example, a common observation was that no system for genetic modification, including conventional methods of plant breeding, is without unintended effects. Another common observation was that “unintended” does not necessarily mean “harmful”. This paper summarizes key points from the information presented at the meeting to provide readers with current viewpoints on these topics.

Interspecific and intraspecific gene variability in a 1-Mb region containing the highest density of NBS-LRR genes found in the melon genome

BACKGROUND

Plant NBS-LRR -resistance genes tend to be found in clusters, which have been shown to be hot spots of genome variability. In melon, half of the 81 predicted NBS-LRR genes group in nine clusters, and a 1 Mb region on linkage group V contains the highest density of R-genes and presence/absence gene polymorphisms found in the melon genome. This region is known to contain the locus of Vat, an agronomically important gene that confers resistance to aphids. However, the presence of duplications makes the sequencing and annotation of R-gene clusters difficult, usually resulting in multi-gapped sequences with higher than average errors.

RESULTS

A 1-Mb sequence that contains the largest NBS-LRR gene cluster found in melon was improved using a strategy that combines Illumina paired-end mapping and PCR-based gap closing. Unknown sequence was decreased by 70% while about 3,000 SNPs and small indels were corrected. As a result, the annotations of 18 of a total of 23 NBS-LRR genes found in this region were modified, including additional coding sequences, amino acid changes, correction of splicing boundaries, or fussion of ORFs in common transcription units. A phylogeny analysis of the R-genes and their comparison with syntenic sequences in other cucurbits point to a pattern of local gene amplifications since the diversification of cucurbits from other families, and through speciation within the family. A candidate Vat gene is proposed based on the sequence similarity between a reported Vat gene from a Korean melon cultivar and a sequence fragment previously absent in the unrefined sequence.

CONCLUSIONS

A sequence refinement strategy allowed substantial improvement of a 1 Mb fragment of the melon genome and the re-annotation of the largest cluster of NBS-LRR gene homologues found in melon. Analysis of the cluster revealed that resistance genes have been produced by sequence duplication in adjacent genome locations since the divergence of cucurbits from other close families, and through the process of speciation within the family a candidate Vat gene was also identified using sequence previously unavailable, which demonstrates the advantages of genome assembly refinements when analyzing complex regions such as those containing clusters of highly similar genes.

Identical substitutions in magnesium chelatase paralogs result in chlorophyll-deficient soybean mutants

The soybean [Glycine max (L.) Merr.] chlorophyll-deficient line MinnGold is a spontaneous mutant characterized by yellow foliage. Map-based cloning and transgenic complementation revealed that the mutant phenotype is caused by a nonsynonymous nucleotide substitution in the third exon of a Mg-chelatase subunit gene (ChlI1a) on chromosome 13. This gene was selected as a candidate for a different yellow foliage mutant, T219H (Y11y11), that had been previously mapped to chromosome 13. Although the phenotypes of MinnGold and T219H are clearly distinct, sequencing of ChlI1a in T219H identified a different nonsynonymous mutation in the third exon, only six base pairs from the MinnGold mutation. This information, along with previously published allelic tests, were used to identify and clone a third yellow foliage mutation, CD-5, which was previously mapped to chromosome 15. This mutation was identified in the ChlI1b gene, a paralog of ChlI1a. Sequencing of the ChlI1b allele in CD-5 identified a nonsynonymous substitution in the third exon that confers an identical amino acid change as the T219H substitution at ChlI1a. Protein sequence alignments of the two Mg-chelatase subunits indicated that the sites of amino acid modification in MinnGold, T219H, and CD-5 are highly conserved among photosynthetic species. These results suggest that amino acid alterations in this critical domain may create competitive inhibitory interactions between the mutant and wild-type ChlI1a and ChlI1b proteins.

A comparative analysis of insertional effects in genetically engineered plants: considerations for pre-market assessments

During genetic engineering, DNA is inserted into a plant’s genome, and such insertions are often accompanied by the insertion of additional DNA, deletions and/or rearrangements. These genetic changes are collectively known as insertional effects, and they have the potential to give rise to unintended traits in plants. In addition, there are many other genetic changes that occur in plants both spontaneously and as a result of conventional breeding practices. Genetic changes similar to insertional effects occur in plants, namely as a result of the movement of transposable elements, the repair of double-strand breaks by non-homologous end-joining, and the intracellular transfer of organelle DNA. Based on this similarity, insertional effects should present a similar level of risk as these other genetic changes in plants, and it is within the context of these genetic changes that insertional effects must be considered. Increased familiarity with genetic engineering techniques and advances in molecular analysis techniques have provided us with a greater understanding of the nature and impact of genetic changes in plants, and this can be used to refine pre-market assessments of genetically engineered plants and food and feeds derived from genetically engineered plants.

De novo assembly of soybean wild relatives for pan-genome analysis of diversity and agronomic traits

Wild relatives of crops are an important source of genetic diversity for agriculture, but their gene repertoire remains largely unexplored. We report the establishment and analysis of a pan-genome of Glycine soja, the wild relative of cultivated soybean Glycine max, by sequencing and de novo assembly of seven phylogenetically and geographically representative accessions. Intergenomic comparisons identified lineage-specific genes and genes with copy number variation or large-effect mutations, some of which show evidence of positive selection and may contribute to variation of agronomic traits such as biotic resistance, seed composition, flowering and maturity time, organ size and final biomass. Approximately 80% of the pan-genome was present in all seven accessions (core), whereas the rest was dispensable and exhibited greater variation than the core genome, perhaps reflecting a role in adaptation to diverse environments. This work will facilitate the harnessing of untapped genetic diversity from wild soybean for enhancement of elite cultivars.

Exploration of presence/absence variation and corresponding polymorphic markers in soybean genome

This study was designed to reveal the genome-wide distribution of presence/absence variation (PAV) and to establish a database of polymorphic PAV markers in soybean. The 33 soybean whole-genome sequences were compared to each other with that of Williams 82 as a reference genome. A total of 33,127 PAVs were detected and 28,912 PAV markers with their primer sequences were designed as the database NJAUSoyPAV_1.0. The PAVs scattered on whole genome while only 518 (1.8%) overlapped with simple sequence repeats (SSRs) in BARCSOYSSR_1.0 database. In a random sample of 800 PAVs, 713 (89.13%) showed polymorphism among the 12 differential genotypes. Using 126 PAVs and 108 SSRs to test a Chinese soybean germplasm collection composed of 828 Glycine soja Sieb. et Zucc. and Glycine max (L.) Merr. accessions, the per locus allele number and its variation appeared less in PAVs than in SSRs. The distinctness among alleles/bands of PCR (polymerase chain reaction) products showed better in PAVs than in SSRs, potential in accurate marker-assisted allele selection. The association mapping results showed SSR + PAV was more powerful than any single marker systems. The NJAUSoyPAV_1.0 database has enriched the source of PCR markers, and may fit the materials with a range of per locus allele numbers, if jointly used with SSR markers.

Nucleotide polymorphism and copy number variant detection using exome capture and next-generation sequencing in the polyploid grass Panicum virgatum

Switchgrass (Panicum virgatum) is a polyploid, outcrossing grass species native to North America and has recently been recognized as a potential biofuel feedstock crop. Significant phenotypic variation including ploidy is present across the two primary ecotypes of switchgrass, referred to as upland and lowland switchgrass. The tetraploid switchgrass genome is approximately 1400 Mbp, split between two subgenomes, with significant repetitive sequence content limiting the efficiency of re-sequencing approaches for determining genome diversity. To characterize genetic diversity in upland and lowland switchgrass as a first step in linking genotype to phenotype, we designed an exome capture probe set based on transcript assemblies that represent approximately 50 Mb of annotated switchgrass exome sequences. We then evaluated and optimized the probe set using solid phase comparative genome hybridization and liquid phase exome capture followed by next-generation sequencing. Using the optimized probe set, we assessed variation in the exomes of eight switchgrass genotypes representing tetraploid lowland and octoploid upland cultivars to benchmark our exome capture probe set design. We identified ample variation in the switchgrass genome including 1,395,501 single nucleotide polymorphisms (SNPs), 8173 putative copy number variants and 3336 presence/absence variants. While the majority of the SNPs (84%) detected was bi-allelic, a substantial number was tri-allelic with limited occurrence of tetra-allelic polymorphisms consistent with the heterozygous and polyploid nature of the switchgrass genome. Collectively, these data demonstrate the efficacy of exome capture for discovery of genome variation in a polyploid species with a large, repetitive and heterozygous genome.

The standing pool of genomic structural variation in a natural population of Mimulus guttatus

Major unresolved questions in evolutionary genetics include determining the contributions of different mutational sources to the total pool of genetic variation in a species, and understanding how these different forms of genetic variation interact with natural selection. Recent work has shown that structural variants (SVs) (insertions, deletions, inversions, and transpositions) are a major source of genetic variation, often outnumbering single nucleotide variants in terms of total bases affected. Despite the near ubiquity of SVs, major questions about their interaction with natural selection remain. For example, how does the allele frequency spectrum of SVs differ when compared with single nucleotide variants? How often do SVs affect genes, and what are the consequences? To begin to address these questions, we have systematically identified and characterized a large set of submicroscopic insertion and deletion (indel) variants (between 1 and 200 kb in length) among ten inbred lines from a single natural population of the plant species Mimulus guttatus. After extensive computational filtering, we focused on a set of 4,142 high-confidence indels that showed an experimental validation rate of 73%. All but one of these indels were less than 200 kb. Although the largest were generally at lower frequencies in the population, a surprising number of large indels are at intermediate frequencies. Although indels overlapping with genes were much rarer than expected by chance, approximately 600 genes were affected by an indel. Nucleotide-binding site leucine-rich repeat (NBS–LRR) defense response genes were the most enriched among the gene families affected. Most indels associated with genes were rare and appeared to be under purifying selection, though we do find four high-frequency derived insertion alleles that show signatures of recent positive selection.

Genome diversity in Brachypodium distachyon: deep sequencing of highly diverse inbred lines

Brachypodium distachyon is small annual grass that has been adopted as a model for the grasses. Its small genome, high-quality reference genome, large germplasm collection, and selfing nature make it an excellent subject for studies of natural variation. We sequenced six divergent lines to identify a comprehensive set of polymorphisms and analyze their distribution and concordance with gene expression. Multiple methods and controls were utilized to identify polymorphisms and validate their quality. mRNA-Seq experiments under control and simulated drought-stress conditions, identified 300 genes with a genotype-dependent treatment response. We showed that large-scale sequence variants had extremely high concordance with altered expression of hundreds of genes, including many with genotype-dependent treatment responses. We generated a deep mRNA-Seq dataset for the most divergent line and created a de novo transcriptome assembly. This led to the discovery of >2400 previously unannotated transcripts and hundreds of genes not present in the reference genome. We built a public database for visualization and investigation of sequence variants among these widely used inbred lines.

Unraveling genomic complexity at a quantitative disease resistance locus in maize

Multiple disease resistance has important implications for plant fitness, given the selection pressure that many pathogens exert directly on natural plant populations and indirectly via crop improvement programs. Evidence of a locus conditioning resistance to multiple pathogens was found in bin 1.06 of the maize genome with the allele from inbred line "Tx303" conditioning quantitative resistance to northern leaf blight (NLB) and qualitative resistance to Stewart's wilt. To dissect the genetic basis of resistance in this region and to refine candidate gene hypotheses, we mapped resistance to the two diseases. Both resistance phenotypes were localized to overlapping regions, with the Stewart's wilt interval refined to a 95.9-kb segment containing three genes and the NLB interval to a 3.60-Mb segment containing 117 genes. Regions of the introgression showed little to no recombination, suggesting structural differences between the inbred lines Tx303 and "B73," the parents of the fine-mapping population. We examined copy number variation across the region using next-generation sequencing data, and found large variation in read depth in Tx303 across the region relative to the reference genome of B73. In the fine-mapping region, association mapping for NLB implicated candidate genes, including a putative zinc finger and pan1. We tested mutant alleles and found that pan1 is a susceptibility gene for NLB and Stewart's wilt. Our data strongly suggest that structural variation plays an important role in resistance conditioned by this region, and pan1, a gene conditioning susceptibility for NLB, may underlie the QTL.

Structural variations in plant genomes

Differences between plant genomes range from single nucleotide polymorphisms to large-scale duplications, deletions and rearrangements. The large polymorphisms are termed structural variants (SVs). SVs have received significant attention in human genetics and were found to be responsible for various chronic diseases. However, little effort has been directed towards understanding the role of SVs in plants. Many recent advances in plant genetics have resulted from improvements in high-resolution technologies for measuring SVs, including microarray-based techniques, and more recently, high-throughput DNA sequencing. In this review we describe recent reports of SV in plants and describe the genomic technologies currently used to measure these SVs.

A roadmap for functional structural variants in the soybean genome

Gene structural variation (SV) has recently emerged as a key genetic mechanism underlying several important phenotypic traits in crop species. We screened a panel of 41 soybean (Glycine max) accessions serving as parents in a soybean nested association mapping population for deletions and duplications in more than 53,000 gene models. Array hybridization and whole genome resequencing methods were used as complementary technologies to identify SV in 1528 genes, or approximately 2.8%, of the soybean gene models. Although SV occurs throughout the genome, SV enrichment was noted in families of biotic defense response genes. Among accessions, SV was nearly eightfold less frequent for gene models that have retained paralogs since the last whole genome duplication event, compared with genes that have not retained paralogs. Increases in gene copy number, similar to that described at the Rhg1 resistance locus, account for approximately one-fourth of the genic SV events. This assessment of soybean SV occurrence presents a target list of genes potentially responsible for rapidly evolving and/or adaptive traits.