Variation in shoot architecture traits and their relationship to canopy coverage and light interception in soybean (Glycine max).

BACKGROUND: In soybeans, faster canopy coverage (CC) is a highly desirable trait but a fully covered canopy is unfavorable to light interception at lower levels in the canopy with most of the incident radiation intercepted at the top of the canopy. Shoot architecture that influences CC is well studied in crops such as maize and wheat, and altering architectural traits has resulted in enhanced yield. However, in soybeans the study of shoot architecture has not been as extensive. RESULTS: This study revealed significant differences in CC among the selected soybean accessions. The rate of CC was found to decrease at the beginning of the reproductive stage (R1) followed by an increase during the R2-R3 stages. Most of the accessions in the study achieved maximum rate of CC between R2-R3 stages. We measured Light interception (LI), defined here as the ratio of Photosynthetically Active Radiation (PAR) transmitted through the canopy to the incoming PAR or the radiation above the canopy. LI was found to be significantly correlated with CC parameters, highlighting the relationship between canopy structure and light interception. The study also explored the impact of plant shape on LI and CO2 assimilation. Plant shape was characterized into distinct quantifiable parameters and by modeling the impact of plant shape on LI and CO2 assimilation, we found that plants with broad and flat shapes at the top maybe more photosynthetically efficient at low light levels, while conical shapes were likely more advantageous when light was abundant. Shoot architecture of plants in this study was described in terms of whole plant, branching and leaf-related traits. There was significant variation for the shoot architecture traits between different accessions, displaying high reliability. We found that that several shoot architecture traits such as plant height, and leaf and internode-related traits strongly influenced CC and LI. CONCLUSION: In conclusion, this study provides insight into the relationship between soybean shoot architecture, canopy coverage, and light interception. It demonstrates that novel shoot architecture traits we have defined here are genetically variable, impact CC and LI and contribute to our understanding of soybean morphology. Correlations between different architecture traits, CC and LI suggest that it is possible to optimize soybean growth without compromising on light transmission within the soybean canopy. In addition, the study underscores the utility of integrating low-cost 2D phenotyping as a practical and cost-effective alternative to more time-intensive 3D or high-tech low-throughput methods. This approach offers a feasible means of studying basic shoot architecture traits at the field level, facilitating a broader and efficient assessment of plant morphology.

Comparative transcriptomics of soybean genotypes with partial resistance towards Phytophthora sojae, Conrad, and M92-220, to moderately susceptible fast neutron mutant soybeans and Sloan.

The breeding of disease-resistant soybeans cultivars to manage Phytophthora root and stem rot caused by the pathogen Phytophthora sojae involves combining quantitative disease resistance (QDR) and Rps gene-mediated resistance. To identify and confirm potential mechanisms of QDR towards P. sojae, we conducted a time course study comparing changes in gene expression among Conrad and M92-220 with high QDR to susceptible genotypes, Sloan and 3 mutants derived from fast neutron (FN) irradiation of M92-220. Differentially expressed genes from Conrad and M92-220 indicated several shared defense-related pathways at the transcriptomic level, but also defense pathways unique to each cultivar such as stilbenoid, diarylheptanoid and gingerol biosynthesis, and monobactam biosynthesis. Gene Ontology pathway analysis showed that the susceptible FN mutants lacked enrichment of three terpenoid related-pathways and two cell wall-related pathways at either one or both timepoints, in contrast to M92-220. The susceptible mutants also lacked enrichment of potentially important KEGG pathways at either one or both timepoints, including sesquiterpenoid and triterpenoid biosynthesis, thiamine metabolism, arachidonic acid, stilbenoid, diarylheptanoid and gingerol biosynthesis, and monobactam biosynthesis. Additionally, thirty-one genes which were differentially expressed in M92-220 following P. sojae infection were not expressed in the mutants. These 31 genes have annotations related to unknown proteins, valine, leucine, and isoleucine biosynthesis and protein and lipid metabolic processes. The results of this study confirm previously proposed mechanisms of QDR, provide evidence for potential novel QDR pathways in M92-220, and furthers our understanding of the complex network associated with QDR mechanisms in soybean towards P. sojae.

Bowman-Birk Inhibitor Mutants of Soybean Generated by CRISPR-Cas9 Reveal Drastic Reductions in Trypsin and Chymotrypsin Inhibitor Activities.

Despite the high quality of soybean protein, raw soybeans and soybean meal cannot be directly included in animal feed mixtures due to the presence of Kunitz (KTi) and Bowman-Birk protease inhibitors (BBis), which reduces animal productivity. Heat treatment can substantially inactivate trypsin and chymotrypsin inhibitors (BBis), but such treatment is energy-intensive, adds expense, and negatively impacts the quality of seed proteins. As an alternative approach, we have employed CRISPR/Cas9 gene editing to create mutations in BBi genes to drastically lower the protease inhibitor content in soybean seed. Agrobacterium-mediated transformation was used to generate several stable transgenic soybean events. These independent CRISPR/Cas9 events were examined in comparison to wild-type plants using Sanger sequencing, proteomic analysis, trypsin/chymotrypsin inhibitor activity assays, and qRT-PCR. Collectively, our results demonstrate the creation of an allelic series of loss-of-function mutations affecting the major BBi gene in soybean. Mutations in two of the highly expressed seed-specific BBi genes lead to substantial reductions in both trypsin and chymotrypsin inhibitor activities.

A chromosome 16 deletion conferring a high sucrose phenotype in soybean.

KEY MESSAGE: Sucrose in soybean seeds is desirable for many end-uses. Increased sucrose contents were discovered to associate with a chromosome 16 deletion resulting from fast neutron irradiation. Soybean is one of the most economically important crops in the United States. A primary end-use of soybean is for livestock feed. Therefore, genetic improvement of seed composition is one of the most important goals in soybean breeding programs. Sucrose is desired in animal feed due to its role as an easily digestible energy source. An elite soybean line was irradiated with fast neutrons and the seed from plants were screened for altered seed composition with near-infrared spectroscopy (NIR). One mutant line, G15FN-54, was found to have higher sucrose content (8-9%) than the parental line (5-6%). Comparative genomic hybridization (CGH) revealed three large deletions on chromosomes (Chrs) 10, 13, and 16 in the mutant, which were confirmed through whole genome sequencing (WGS). A bi-parental population derived from the mutant G15FN-54 and the cultivar Benning was developed to conduct a bulked segregant analysis (BSA) with SoySNP50K BeadChips, revealing that the deletion on Chr 16 might be responsible for the altered phenotype. The mapping result using the bi-parental population confirmed that the deletion on Chr 16 conferred elevated sucrose content and a total of 21 genes are located within this Chr 16 deletion. NIR and high-pressure liquid chromatography (HPLC) were used to confirm the stability of the phenotype across generations in the bi-parental population. The mutation will be useful to understand the genetic control of soybean seed sucrose content.

Plant Genome Editing and the Relevance of Off-Target Changes.

Site-directed nucleases (SDNs) used for targeted genome editing are powerful new tools to introduce precise genetic changes into plants. Like traditional approaches, such as conventional crossing and induced mutagenesis, genome editing aims to improve crop yield and nutrition. Next-generation sequencing studies demonstrate that across their genomes, populations of crop species typically carry millions of single nucleotide polymorphisms and many copy number and structural variants. Spontaneous mutations occur at rates of ∼10-8 to 10-9 per site per generation, while variation induced by chemical treatment or ionizing radiation results in higher mutation rates. In the context of SDNs, an off-target change or edit is an unintended, nonspecific mutation occurring at a site with sequence similarity to the targeted edit region. SDN-mediated off-target changes can contribute to a small number of additional genetic variants compared to those that occur naturally in breeding populations or are introduced by induced-mutagenesis methods. Recent studies show that using computational algorithms to design genome editing reagents can mitigate off-target edits in plants. Finally, crops are subject to strong selection to eliminate off-type plants through well-established multigenerational breeding, selection, and commercial variety development practices. Within this context, off-target edits in crops present no new safety concerns compared to other breeding practices. The current generation of genome editing technologies is already proving useful to develop new plant varieties with consumer and farmer benefits. Genome editing will likely undergo improved editing specificity along with new developments in SDN delivery and increasing genomic characterization, further improving reagent design and application.

Mining Fiskeby III and Mandarin (Ottawa) Expression Profiles to Understand Iron Stress Tolerant Responses in Soybean.

The soybean (Glycine max L. merr) genotype Fiskeby III is highly resistant to a multitude of abiotic stresses, including iron deficiency, incurring only mild yield loss during stress conditions. Conversely, Mandarin (Ottawa) is highly susceptible to disease and suffers severe phenotypic damage and yield loss when exposed to abiotic stresses such as iron deficiency, a major challenge to soybean production in the northern Midwestern United States. Using RNA-seq, we characterize the transcriptional response to iron deficiency in both Fiskeby III and Mandarin (Ottawa) to better understand abiotic stress tolerance. Previous work by our group identified a quantitative trait locus (QTL) on chromosome 5 associated with Fiskeby III iron efficiency, indicating Fiskeby III utilizes iron deficiency stress mechanisms not previously characterized in soybean. We targeted 10 of the potential candidate genes in the Williams 82 genome sequence associated with the QTL using virus-induced gene silencing. Coupling virus-induced gene silencing with RNA-seq, we identified a single high priority candidate gene with a significant impact on iron deficiency response pathways. Characterization of the Fiskeby III responses to iron stress and the genes underlying the chromosome 5 QTL provides novel targets for improved abiotic stress tolerance in soybean.

Construction and comparison of three reference-quality genome assemblies for soybean.

We report reference-quality genome assemblies and annotations for two accessions of soybean (Glycine max) and for one accession of Glycine soja, the closest wild relative of G. max. The G. max assemblies provided are for widely used US cultivars: the northern line Williams 82 (Wm82) and the southern line Lee. The Wm82 assembly improves the prior published assembly, and the Lee and G. soja assemblies are new for these accessions. Comparisons among the three accessions show generally high structural conservation, but nucleotide difference of 1.7 single-nucleotide polymorphisms (snps) per kb between Wm82 and Lee, and 4.7 snps per kb between these lines and G. soja. snp distributions and comparisons with genotypes of the Lee and Wm82 parents highlight patterns of introgression and haplotype structure. Comparisons against the US germplasm collection show placement of the sequenced accessions relative to global soybean diversity. Analysis of a pan-gene collection shows generally high conservation, with variation occurring primarily in genomically clustered gene families. We found approximately 40-42 inversions per chromosome between either Lee or Wm82v4 and G. soja, and approximately 32 inversions per chromosome between Wm82 and Lee. We also investigated five domestication loci. For each locus, we found two different alleles with functional differences between G. soja and the two domesticated accessions. The genome assemblies for multiple cultivated accessions and for the closest wild ancestor of soybean provides a valuable set of resources for identifying causal variants that underlie traits for the domestication and improvement of soybean, serving as a basis for future research and crop improvement efforts for this important crop species.

Integration, abundance, and transmission of mutations and transgenes in a series of CRISPR/Cas9 soybean lines.

BACKGROUND: As with many plant species, current genome editing strategies in soybean are initiated by stably transforming a gene that encodes an engineered nuclease into the genome. Expression of the transgene results in a double-stranded break and repair at the targeted locus, oftentimes resulting in mutation(s) at the intended site. As soybean is a self-pollinating species with 20 chromosome pairs, the transgene(s) in the T0 plant are generally expected to be unlinked to the targeted mutation(s), and the transgene(s)/mutation(s) should independently assort into the T1 generation, resulting in Mendellian combinations of transgene presence/absence and allelic states within the segregating family. This prediction, however, is not always consistent with observed results. RESULTS: In this study, we investigated inheritance patterns among three different CRISPR/Cas9 transgenes and their respective induced mutations in segregating soybean families. Next-generation resequencing of four T0 plants and four T1 progeny plants, followed by broader assessments of the segregating families, revealed both expected and unexpected patterns of inheritance among the different lineages. These unexpected patterns included: (1) A family in which T0 transgenes and mutations were not transmitted to progeny; (2) A family with four unlinked transgene insertions, including two respectively located at paralogous CRISPR target break sites; (3) A family in which mutations were observed and transmitted, but without evidence of transgene integration nor transmission. CONCLUSIONS: Genome resequencing provides high-resolution of transgene integration structures and gene editing events. Segregation patterns of these events can be complicated by several potential mechanisms. This includes, but is not limited to, plant chimeras, multiple unlinked transgene integrations, editing of intended and paralogous targets, linkage between the transgene integration and target site, and transient expression of the editing reagents without transgene integration into the host genome.

Screening populations for copy number variation using genotyping-by-sequencing: a proof of concept using soybean fast neutron mutants.

BACKGROUND: The effective use of mutant populations for reverse genetic screens relies on the population-wide characterization of the induced mutations. Genome- and population-wide characterization of the mutations found in fast neutron populations has been hindered, however, by the wide range of mutations generated and the lack of affordable technologies to detect DNA sequence changes. In this study, we therefore aimed to test whether genotyping-by-sequencing (GBS) technology could be used to characterize copy number variation (CNV) induced by fast neutrons in a soybean mutant population. RESULTS: We called CNVs from GBS data in 79 soybean mutants and assessed the sensitivity and precision of this approach by validating our results against array comparative genomic hybridization (aCGH) data for 19 of these mutants as well as targeted PCR and ddPCR assays for a representative subset of the smallest events detected by GBS. Our GBS pipeline detected 55 of the 96 events found by aCGH, with approximate detection thresholds of 60 kb, 500 kb and 1 Mb for homozygous deletions, hemizygous deletions and duplications, respectively. Among the whole set of 79 mutants, the GBS data revealed 105 homozygous deletions, 32 hemizygous deletions and 19 duplications. This included several extremely large events, exhibiting maximum sizes of ~ 11.2 Mb for a homozygous deletion, ~ 11.6 Mb for a hemizygous deletion, and ~ 50 Mb for a duplication. CONCLUSIONS: This study provides a proof of concept that GBS can be used as an affordable high-throughput method for assessing CNVs in fast neutron mutants. The modularity of this GBS approach allows combining as many different libraries or sequencing runs as is necessary for reaching the goals of a particular study. This method should enable the low-cost genome-wide characterization of hundreds to thousands of individuals in fast neutron mutant populations or any population with large genomic deletions and duplications.

Genomic changes and biochemical alterations of seed protein and oil content in a subset of fast neutron induced soybean mutants.

BACKGROUND: Soybean is subjected to genetic manipulation by breeding, mutation, and transgenic approaches to produce value-added quality traits. Among those genetic approaches, mutagenesis through fast neutrons radiation is intriguing because it yields a variety of mutations, including single/multiple gene deletions and/or duplications. Characterizing the seed composition of the fast neutron mutants and its relationship with gene mutation is useful towards understanding oil and protein traits in soybean. RESULTS: From a large population of fast neutron mutagenized plants, we selected ten mutants based on a screening of total oil and protein content using near infra-red spectroscopy. These ten mutants were regrown, and the seeds were analyzed for oil by GC-MS, protein profiling by SDS-PAGE and gene mapping by comparative genomic hybridization. The mutant 2R29C14Cladecr233cMN15 (nicknamed in this study as L10) showed higher protein and lower oil content compared to the wild type, followed by three other lines (nicknamed in this study as L03, L05, and L06). We characterized the fatty acid methyl esters profile of the trans-esterified oil and found the presence of five major fatty acids (palmitic, stearic, oleic, linoleic, and linolenic acids) at varying proportions among the mutants. Protein profile using SDS-PAGE of the ten mutants did exhibit discernable variation between storage (glycinin and ß-conglycinin) and anti-nutritional factor (trypsin inhibitor) proteins. In addition, we physically mapped the position of the gene deletions or duplications in each mutant using comparative genomic hybridization. CONCLUSION: Characterization of oil and protein profile in soybean fast neutron mutants will assist scientist and breeders to develop new value-added soybeans with improved protein and oil quality traits.

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

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

Identification and characterization of a fast-neutron-induced mutant with elevated seed protein content in soybean.

KEY MESSAGE: Protein content of soybean is critical for utility of soybean meal. A fast-neutron-induced deletion on chromosome 12 was found to be associated with increased protein content. Soybean seed composition affects the utility of soybean, and improving seed composition is an essential breeding goal. Fast neutron radiation introduces genomic mutations resulting in novel variation for traits of interest. Two elite soybean lines were irradiated with fast neutrons and screened for altered seed composition. Twenty-three lines with altered protein, oil, or sucrose content were selected based on near-infrared spectroscopy data from five environments and yield tested at five locations. Mutants with significantly increased protein averaged 19.1-36.8 g kg-1 more protein than the parents across 10 environments. Comparative genomic hybridization (CGH) identified putative mutations in a mutant, G15FN-12, that has 36.8 g kg-1 higher protein than the parent genotype, and whole genome sequencing (WGS) of the mutant has confirmed these mutations. An F2:3 population was developed from G15FN-12 to determine association between genomic changes and increased protein content. Bulked segregant analysis of the population using the SoySNP50K BeadChip identified a CGH- and WGS-confirmed deletion on chromosome 12 to be responsible for elevated protein content. The population was genotyped using a KASP marker designed at the mutation region, and significant association (P < 0.0001) between the deletion on chromosome 12 and elevated protein content was observed and confirmed in the F3:4 generation. The F2 segregants homozygous for the deletion averaged 27 g kg-1 higher seed protein and 8 g kg-1 lower oil than homozygous wild-type segregants. Mutants with altered seed composition are a new resource for gene function studies and provide elite materials for genetic improvement of seed composition.

The importance of genotype identity, genetic heterogeneity, and bioinformatic handling for properly assessing genomic variation in transgenic plants.

BACKGROUND: The advent of -omics technologies has enabled the resolution of fine molecular differences among individuals within a species. DNA sequence variations, such as single nucleotide polymorphisms or small deletions, can be tabulated for many kinds of genotype comparisons. However, experimental designs and analytical approaches are replete with ways to overestimate the level of variation present within a given sample. Analytical pipelines that do not apply proper thresholds nor assess reproducibility among samples are susceptible to calling false-positive variants. Furthermore, issues with sample genotype identity or failing to account for heterogeneity in reference genotypes may lead to misinterpretations of standing variants as polymorphisms derived de novo. RESULTS: A recent publication that featured the analysis of RNA-sequencing data in three transgenic soybean event series appeared to overestimate the number of sequence variants identified in plants that were exposed to a tissue culture based transformation process. We reanalyzed these data with a stringent set of criteria and demonstrate three different factors that lead to variant overestimation, including issues related to the genetic identity of the background genotype, unaccounted genetic heterogeneity in the reference genome, and insufficient bioinformatics filtering. CONCLUSIONS: This study serves as a cautionary tale to users of genomic and transcriptomic data that wish to assess the molecular variation attributable to tissue culture and transformation processes. Moreover, accounting for the factors that lead to sequence variant overestimation is equally applicable to samples derived from other germplasm sources, including chemical or irradiation mutagenesis and genome engineering (e.g., CRISPR) processes.

CRISPR/Cas9 and TALENs generate heritable mutations for genes involved in small RNA processing of Glycine max and Medicago truncatula.

Processing of double-stranded RNA precursors into small RNAs is an essential regulator of gene expression in plant development and stress response. Small RNA processing requires the combined activity of a functionally diverse group of molecular components. However, in most of the plant species, there are insufficient mutant resources to functionally characterize each encoding gene. Here, mutations in loci encoding protein machinery involved in small RNA processing in soya bean and Medicago truncatula were generated using the CRISPR/Cas9 and TAL-effector nuclease (TALEN) mutagenesis platforms. An efficient CRISPR/Cas9 reagent was used to create a bi-allelic double mutant for the two soya bean paralogous Double-stranded RNA-binding2 (GmDrb2a and GmDrb2b) genes. These mutations, along with a CRISPR/Cas9-generated mutation of the M. truncatula Hua enhancer1 (MtHen1) gene, were determined to be germ-line transmissible. Furthermore, TALENs were used to generate a mutation within the soya bean Dicer-like2 gene. CRISPR/Cas9 mutagenesis of the soya bean Dicer-like3 gene and the GmHen1a gene was observed in the T0 generation, but these mutations failed to transmit to the T1 generation. The irregular transmission of induced mutations and the corresponding transgenes was investigated by whole-genome sequencing to reveal a spectrum of non-germ-line-targeted mutations and multiple transgene insertion events. Finally, a suite of combinatorial mutant plants were generated by combining the previously reported Gmdcl1a, Gmdcl1b and Gmdcl4b mutants with the Gmdrb2ab double mutant. Altogether, this study demonstrates the synergistic use of different genome engineering platforms to generate a collection of useful mutant plant lines for future study of small RNA processing in legume crops.

Validating Genome-Wide Association Candidates Controlling Quantitative Variation in Nodulation.

Genome-wide association (GWA) studies offer the opportunity to identify genes that contribute to naturally occurring variation in quantitative traits. However, GWA relies exclusively on statistical association, so functional validation is necessary to make strong claims about gene function. We used a combination of gene-disruption platforms (Tnt1 retrotransposons, hairpin RNA-interference constructs, and CRISPR/Cas9 nucleases) together with randomized, well-replicated experiments to evaluate the function of genes that an earlier GWA study in Medicago truncatula had identified as candidates contributing to variation in the symbiosis between legumes and rhizobia. We evaluated ten candidate genes found in six clusters of strongly associated single nucleotide polymorphisms, selected on the basis of their strength of statistical association, proximity to annotated gene models, and root or nodule expression. We found statistically significant effects on nodule production for three candidate genes, each validated in two independent mutants. Annotated functions of these three genes suggest their contributions to quantitative variation in nodule production occur through processes not previously connected to nodulation, including phosphorous supply and salicylic acid-related defense response. These results demonstrate the utility of GWA combined with reverse mutagenesis technologies to discover and validate genes contributing to naturally occurring variation in quantitative traits. The results highlight the potential for GWA to complement forward genetics in identifying the genetic basis of ecologically and economically important traits.

Strategies for optimizing BioNano and Dovetail explored through a second reference quality assembly for the legume model, Medicago truncatula.

BACKGROUND: Third generation sequencing technologies, with sequencing reads in the tens- of kilo-bases, facilitate genome assembly by spanning ambiguous regions and improving continuity. This has been critical for plant genomes, which are difficult to assemble due to high repeat content, gene family expansions, segmental and tandem duplications, and polyploidy. Recently, high-throughput mapping and scaffolding strategies have further improved continuity. Together, these long-range technologies enable quality draft assemblies of complex genomes in a cost-effective and timely manner. RESULTS: Here, we present high quality genome assemblies of the model legume plant, Medicago truncatula (R108) using PacBio, Dovetail Chicago (hereafter, Dovetail) and BioNano technologies. To test these technologies for plant genome assembly, we generated five assemblies using all possible combinations and ordering of these three technologies in the R108 assembly. While the BioNano and Dovetail joins overlapped, they also showed complementary gains in continuity and join numbers. Both technologies spanned repetitive regions that PacBio alone was unable to bridge. Combining technologies, particularly Dovetail followed by BioNano, resulted in notable improvements compared to Dovetail or BioNano alone. A combination of PacBio, Dovetail, and BioNano was used to generate a high quality draft assembly of R108, a M. truncatula accession widely used in studies of functional genomics. As a test for the usefulness of the resulting genome sequence, the new R108 assembly was used to pinpoint breakpoints and characterize flanking sequence of a previously identified translocation between chromosomes 4 and 8, identifying more than 22.7 Mb of novel sequence not present in the earlier A17 reference assembly. CONCLUSIONS: Adding Dovetail followed by BioNano data yielded complementary improvements in continuity over the original PacBio assembly. This strategy proved efficient and cost-effective for developing a quality draft assembly compared to traditional reference assemblies.

Exploring structural variation and gene family architecture with De Novo assemblies of 15 Medicago genomes.

BACKGROUND: Previous studies exploring sequence variation in the model legume, Medicago truncatula, relied on mapping short reads to a single reference. However, read-mapping approaches are inadequate to examine large, diverse gene families or to probe variation in repeat-rich or highly divergent genome regions. De novo sequencing and assembly of M. truncatula genomes enables near-comprehensive discovery of structural variants (SVs), analysis of rapidly evolving gene families, and ultimately, construction of a pan-genome. RESULTS: Genome-wide synteny based on 15 de novo M. truncatula assemblies effectively detected different types of SVs indicating that as much as 22% of the genome is involved in large structural changes, altogether affecting 28% of gene models. A total of 63 million base pairs (Mbp) of novel sequence was discovered, expanding the reference genome space for Medicago by 16%. Pan-genome analysis revealed that 42% (180 Mbp) of genomic sequences is missing in one or more accession, while examination of de novo annotated genes identified 67% (50,700) of all ortholog groups as dispensable - estimates comparable to recent studies in rice, maize and soybean. Rapidly evolving gene families typically associated with biotic interactions and stress response were found to be enriched in the accession-specific gene pool. The nucleotide-binding site leucine-rich repeat (NBS-LRR) family, in particular, harbors the highest level of nucleotide diversity, large effect single nucleotide change, protein diversity, and presence/absence variation. However, the leucine-rich repeat (LRR) and heat shock gene families are disproportionately affected by large effect single nucleotide changes and even higher levels of copy number variation. CONCLUSIONS: Analysis of multiple M. truncatula genomes illustrates the value of de novo assemblies to discover and describe structural variation, something that is often under-estimated when using read-mapping approaches. Comparisons among the de novo assemblies also indicate that different large gene families differ in the architecture of their structural variation.

Hybrid assembly with long and short reads improves discovery of gene family expansions.

BACKGROUND: Long-read and short-read sequencing technologies offer competing advantages for eukaryotic genome sequencing projects. Combinations of both may be appropriate for surveys of within-species genomic variation. METHODS: We developed a hybrid assembly pipeline called "Alpaca" that can operate on 20X long-read coverage plus about 50X short-insert and 50X long-insert short-read coverage. To preclude collapse of tandem repeats, Alpaca relies on base-call-corrected long reads for contig formation. RESULTS: Compared to two other assembly protocols, Alpaca demonstrated the most reference agreement and repeat capture on the rice genome. On three accessions of the model legume Medicago truncatula, Alpaca generated the most agreement to a conspecific reference and predicted tandemly repeated genes absent from the other assemblies. CONCLUSION: Our results suggest Alpaca is a useful tool for investigating structural and copy number variation within de novo assemblies of sampled populations.

The Role of Deleterious Substitutions in Crop Genomes.

Populations continually incur new mutations with fitness effects ranging from lethal to adaptive. While the distribution of fitness effects of new mutations is not directly observable, many mutations likely either have no effect on organismal fitness or are deleterious. Historically, it has been hypothesized that a population may carry many mildly deleterious variants as segregating variation, which reduces the mean absolute fitness of the population. Recent advances in sequencing technology and sequence conservation-based metrics for inferring the functional effect of a variant permit examination of the persistence of deleterious variants in populations. The issue of segregating deleterious variation is particularly important for crop improvement, because the demographic history of domestication and breeding allows deleterious variants to persist and reach moderate frequency, potentially reducing crop productivity. In this study, we use exome resequencing of 15 barley accessions and genome resequencing of 8 soybean accessions to investigate the prevalence of deleterious single nucleotide polymorphisms (SNPs) in the protein-coding regions of the genomes of two crops. We conclude that individual cultivars carry hundreds of deleterious SNPs on average, and that nonsense variants make up a minority of deleterious SNPs. Our approach identifies known phenotype-altering variants as deleterious more frequently than the genome-wide average, suggesting that putatively deleterious variants are likely to affect phenotypic variation. We also report the implementation of a SNP annotation tool BAD_Mutations that makes use of a likelihood ratio test based on alignment of all currently publicly available Angiosperm genomes.

Transcriptomic basis of genome by genome variation in a legume-rhizobia mutualism.

In the legume-rhizobia mutualism, the benefit each partner derives from the other depends on the genetic identity of both host and rhizobial symbiont. To gain insight into the extent of genome × genome interactions on hosts at the molecular level and to identify potential mechanisms responsible for the variation, we examined host gene expression within nodules (the plant organ where the symbiosis occurs) of four genotypes of Medicago truncatula grown with either Ensifer meliloti or E. medicae symbionts. These host × symbiont combinations show significant variation in nodule and biomass phenotypes. Likewise, combinations differ in their transcriptomes: host, symbiont and host × symbiont affected the expression of 70%, 27% and 21%, respectively, of the approximately 27,000 host genes expressed in nodules. Genes with the highest levels of expression often varied between hosts and/or symbiont strain and include leghemoglobins that modulate oxygen availability and hundreds of Nodule Cysteine-Rich (NCR) peptides involved in symbiont differentiation and viability in nodules. Genes with host × symbiont-dependent expression were enriched for functions related to resource exchange between partners (sulphate/iron/amino acid transport and dicarboxylate/amino acid synthesis). These enrichments suggest mechanisms for host control of the currencies of the mutualism. The transcriptome of M. truncatula accession HM101 (A17), the reference genome used for most molecular research, was less affected by symbiont identity than the other hosts. These findings underscore the importance of assessing the molecular basis of variation in ecologically important traits, particularly those involved in biotic interactions, in multiple genetic contexts.