Agile Genetics: Single gene resolution without the fuss.

Gene discovery reveals new biology, expands the utility of marker-assisted selection, and enables targeted mutagenesis. Still, such discoveries can take over a decade. We present a general strategy, "Agile Genetics," that uses nested, structured populations to overcome common limits on gene resolution. Extensive simulation work on realistic genetic architectures shows that, at population sizes of >5000 samples, single gene-resolution can be achieved using bulk segregant pools. At this scale, read depth and technical replication become major drivers of resolution. Emerging enrichment methods to address coverage are on the horizon; we describe one possibility - iterative depth sequencing (ID-seq). In addition, graph-based pangenomics in experimental populations will continue to maximize accuracy and improve interpretation. Based on this merger of agronomic scale with molecular and bioinformatic innovation, we predict a new age of rapid gene discovery.

Aspergillus flavus pangenome (AflaPan) uncovers novel aflatoxin and secondary metabolite associated gene clusters.

BACKGROUND: Aspergillus flavus is an important agricultural and food safety threat due to its production of carcinogenic aflatoxins. It has high level of genetic diversity that is adapted to various environments. Recently, we reported two reference genomes of A. flavus isolates, AF13 (MAT1-2 and highly aflatoxigenic isolate) and NRRL3357 (MAT1-1 and moderate aflatoxin producer). Where, an insertion of 310 kb in AF13 included an aflatoxin producing gene bZIP transcription factor, named atfC. Observations of significant genomic variants between these isolates of contrasting phenotypes prompted an investigation into variation among other agricultural isolates of A. flavus with the goal of discovering novel genes potentially associated with aflatoxin production regulation. Present study was designed with three main objectives: (1) collection of large number of A. flavus isolates from diverse sources including maize plants and field soils; (2) whole genome sequencing of collected isolates and development of a pangenome; and (3) pangenome-wide association study (Pan-GWAS) to identify novel secondary metabolite cluster genes. RESULTS: Pangenome analysis of 346 A. flavus isolates identified a total of 17,855 unique orthologous gene clusters, with mere 41% (7,315) core genes and 59% (10,540) accessory genes indicating accumulation of high genomic diversity during domestication. 5,994 orthologous gene clusters in accessory genome not annotated in either the A. flavus AF13 or NRRL3357 reference genomes. Pan-genome wide association analysis of the genomic variations identified 391 significant associated pan-genes associated with aflatoxin production. Interestingly, most of the significantly associated pan-genes (94%; 369 associations) belonged to accessory genome indicating that genome expansion has resulted in the incorporation of new genes associated with aflatoxin and other secondary metabolites. CONCLUSION: In summary, this study provides complete pangenome framework for the species of Aspergillus flavus along with associated genes for pathogen survival and aflatoxin production. The large accessory genome indicated large genome diversity in the species A. flavus, however AflaPan is a closed pangenome represents optimum diversity of species A. flavus. Most importantly, the newly identified aflatoxin producing gene clusters will be a new source for seeking aflatoxin mitigation strategies and needs new attention in research.

A haplotype-resolved chromosome-level assembly and annotation of European hazelnut (C. avellana cv. Jefferson) provides insight into mechanisms of eastern filbert blight resistance.

European hazelnut (Corylus avellana L.) is an important tree nut crop. Hazelnut production in North America is currently limited in scalability due to Anisogramma anomala, a fungal pathogen that causes Eastern Filbert Blight (EFB) disease in hazelnut. Successful deployment of EFB resistant cultivars has been limited to the state of Oregon, where the breeding program at Oregon State University (OSU) has released cultivars with a dominant allele at a single resistance locus identified by classical breeding, linkage mapping, and molecular markers. C. avellana cultivar "Jefferson" is resistant to the predominant EFB biotype in Oregon and has been selected by the OSU breeding program as a model for hazelnut genetic and genomic research. Here, we present a near complete, haplotype-resolved chromosome-level hazelnut genome assembly for "Jefferson". This new assembly is a significant improvement over a previously published genome draft. Analysis of genomic regions linked to EFB resistance and self-incompatibility confirmed haplotype splitting and identified new gene candidates that are essential for downstream molecular marker development, thereby facilitating breeding efforts.

The groundnut improvement network for Africa (GINA) germplasm collection: a unique genetic resource for breeding and gene discovery.

Cultivated peanut or groundnut (Arachis hypogaea L.) is a grain legume grown in many developing countries by smallholder farmers for food, feed, and/or income. The speciation of the cultivated species, that involved polyploidization followed by domestication, greatly reduced its variability at the DNA level. Mobilizing peanut diversity is a prerequisite for any breeding program for overcoming the main constraints that plague production and for increasing yield in farmer fields. In this study, the Groundnut Improvement Network for Africa assembled a collection of 1,049 peanut breeding lines, varieties, and landraces from 9 countries in Africa. The collection was genotyped with the Axiom_Arachis2 48K SNP array and 8,229 polymorphic single nucleotide polymorphism (SNP) markers were used to analyze the genetic structure of this collection and quantify the level of genetic diversity in each breeding program. A supervised model was developed using dapc to unambiguously assign 542, 35, and 172 genotypes to the Spanish, Valencia, and Virginia market types, respectively. Distance-based clustering of the collection showed a clear grouping structure according to subspecies and market types, with 73% of the genotypes classified as fastigiata and 27% as hypogaea subspecies. Using STRUCTURE, the global structuration was confirmed and showed that, at a minimum membership of 0.8, 76% of the varieties that were not assigned by dapc were actually admixed. This was particularly the case of most of the genotype of the Valencia subgroup that exhibited admixed genetic heritage. The results also showed that the geographic origin (i.e. East, Southern, and West Africa) did not strongly explain the genetic structure. The gene diversity managed by each breeding program, measured by the expected heterozygosity, ranged from 0.25 to 0.39, with the Niger breeding program having the lowest diversity mainly because only lines that belong to the fastigiata subspecies are used in this program. Finally, we developed a core collection composed of 300 accessions based on breeding traits and genetic diversity. This collection, which is composed of 205 genotypes of fastigiata subspecies (158 Spanish and 47 Valencia) and 95 genotypes of hypogaea subspecies (all Virginia), improves the genetic diversity of each individual breeding program and is, therefore, a unique resource for allele mining and breeding.

A derived weedy rice × ancestral cultivar cross identifies evolutionarily relevant weediness QTLs.

Weedy rice (Oryza spp.) is a weedy relative of the cultivated rice that competes with the crop and causes significant production loss. The BHA (blackhull awned) US weedy rice group has evolved from aus cultivated rice and differs from its ancestors in several important weediness traits, including flowering time, plant height and seed shattering. Prior attempts to determine the genetic basis of weediness traits in plants using linkage mapping approaches have not often considered weed origins. However, the timing of divergence between crossed parents can affect the detection of quantitative trait loci (QTL) relevant to the evolution of weediness. Here, we used a QTL-seq approach that combines bulked segregant analysis and high-throughput whole genome resequencing to map the three important weediness traits in an F2 population derived from a cross between BHA weedy rice with an ancestral aus cultivar. We compared these QTLs with those previously detected in a cross of BHA with a more distantly related crop, indica. We identified multiple QTLs that overlapped with regions under selection during the evolution of weedy BHA rice and some candidate genes possibly underlying the evolution weediness traits in BHA. We showed that QTLs detected with ancestor-descendant crosses are more likely to be involved in the evolution of weediness traits than those detected from crosses of more diverged taxa.

Integrating de novo QTL-seq and linkage mapping to identify quantitative trait loci conditioning physiological resistance and avoidance to white mold disease in dry bean.

White mold (WM), caused by the ubiquitous fungus Sclerotinia sclerotiorum, is a devastating disease that limits production and quality of dry bean globally. In the present study, classic linkage mapping combined with QTL-seq were employed in two recombinant inbred line (RIL) populations, "Montrose"/I9365-25 (M25) and "Raven"/I9365-31 (R31), with the initial goal of fine-mapping QTL WM5.4 and WM7.5 that condition WM resistance. The RILs were phenotyped for WM reactions under greenhouse (straw test) and field environments. The general region of WM5.4 and WM7.5 were reconfirmed with both mapping strategies within each population. Combining the results from both mapping strategies, WM5.4 was delimited to a 22.60-36.25 Mb interval in the heterochromatic regions on Pv05, while WM7.5 was narrowed to a 0.83 Mb (3.99-4.82 Mb) region on the Pv07 chromosome. Furthermore, additional QTL WM2.2a (3.81-7.24 Mb), WM2.2b (11.18-17.37 Mb, heterochromatic region), and WM2.2c (23.33-25.94 Mb) were mapped to a narrowed genomic interval on Pv02 and WM4.2 in a 0.89 Mb physical interval at the distal end of Pv04 chromosome. Gene models encoding gibberellin 2-oxidase proteins regulating plant architecture are likely candidate genes associated with WM2.2a resistance. Nine gene models encoding a disease resistance protein (quinone reductase family protein and ATWRKY69) found within the WM5.4 QTL interval are putative candidate genes. Clusters of 13 and 5 copies of gene models encoding cysteine-rich receptor-like kinase and receptor-like protein kinase-related family proteins, respectively, are potential candidate genes associated with WM7.5 resistance and most likely trigger physiological resistance to WM. Acquired knowledge of the narrowed major QTL intervals, flanking markers, and candidate genes provides promising opportunities to develop functional molecular markers to implement marker-assisted selection for WM resistant dry bean cultivars.

Linked candidate genes of different functions for white mold resistance in common bean (Phaseolus vulgaris L) are identified by multiple QTL mapping approaches.

White mold (WM) is a major disease in common bean (Phaseolus vulgaris L.), and its complex quantitative genetic control limits the development of WM resistant cultivars. WM2.2, one of the nine meta-QTL with a major effect on WM tolerance, explains up to 35% of the phenotypic variation and was previously mapped to a large genomic interval on Pv02. Our objective was to narrow the interval of this QTL using combined approach of classic QTL mapping and QTL-based bulk segregant analysis (BSA), and confirming those results with Khufu de novo QTL-seq. The phenotypic and genotypic data from two RIL populations, 'Raven'/I9365-31 (R31) and 'AN-37'/PS02-029C-20 (Z0726-9), were used to select resistant and susceptible lines to generate subpopulations for bulk DNA sequencing. The QTL physical interval was determined by considering overlapping interval of the identified QTL or peak region in both populations by three independent QTL mapping analyses. Our findings revealed that meta-QTL WM2.2 consists of three regions, WM2.2a (4.27-5.76 Mb; euchromatic), WM 2.2b (12.19 to 17.61 Mb; heterochromatic), and WM2.2c (23.01-25.74 Mb; heterochromatic) found in both populations. Gene models encoding for gibberellin 2-oxidase 8, pentatricopeptide repeat, and heat-shock proteins are the likely candidate genes associated with WM2.2a resistance. A TIR-NBS-LRR class of disease resistance protein (Phvul.002G09200) and LRR domain containing family proteins are potential candidate genes associated with WM2.2b resistance. Nine gene models encoding disease resistance protein [pathogenesis-related thaumatin superfamily protein and disease resistance-responsive (dirigent-like protein) family protein etc] found within the WM2.2c QTL interval are putative candidate genes. WM2.2a region is most likely associated with avoidance mechanisms while WM2.2b and WM2.2c regions trigger physiological resistance based on putative candidate genes.

Genome-wide association studies reveal novel loci for resistance to groundnut rosette disease in the African core groundnut collection.

KEY MESSAGE: We identified markers associated with GRD resistance after screening an Africa-wide core collection across three seasons in Uganda Groundnut is cultivated in several African countries where it is a major source of food, feed and income. One of the major constraints to groundnut production in Africa is groundnut rosette disease (GRD), which is caused by a complex of three agents: groundnut rosette assistor luteovirus, groundnut rosette umbravirus and its satellite RNA. Despite several years of breeding for GRD resistance, the genetics of the disease is not fully understood. The objective of the current study was to use the African core collection to establish the level of genetic variation in their response to GRD, and to map genomic regions responsible for the observed resistance. The African groundnut core genotypes were screened across two GRD hotspot locations in Uganda (Nakabango and Serere) for 3 seasons. The Area Under Disease Progress Curve combined with 7523 high quality SNPs were analyzed to establish marker-trait associations (MTAs). Genome-Wide Association Studies based on Enriched Compressed Mixed Linear Model detected 32 MTAs at Nakabango: 21 on chromosome A04, 10 on B04 and 1 on B08. Two of the significant markers were localised on the exons of a putative TIR-NBS-LRR disease resistance gene on chromosome A04. Our results suggest the likely involvement of major genes in the resistance to GRD but will need to be further validated with more comprehensive phenotypic and genotypic datasets. The markers identified in the current study will be developed into routine assays and validated for future genomics-assisted selection for GRD resistance in groundnut.

Microscopic and Transcriptomic Analyses of Dalbergoid Legume Peanut Reveal a Divergent Evolution Leading to Nod-Factor-Dependent Epidermal Crack-Entry and Terminal Bacteroid Differentiation.

Root nodule symbiosis (RNS) is the pillar behind sustainable agriculture and plays a pivotal role in the environmental nitrogen cycle. Most of the genetic, molecular, and cell-biological knowledge on RNS comes from model legumes that exhibit a root-hair mode of bacterial infection, in contrast to the Dalbergoid legumes exhibiting crack-entry of rhizobia. As a step toward understanding this important group of legumes, we have combined microscopic analysis and temporal transcriptome to obtain a dynamic view of plant gene expression during Arachis hypogaea (peanut) nodule development. We generated comprehensive transcriptome data by mapping the reads to A. hypogaea, and two diploid progenitor genomes. Additionally, we performed BLAST searches to identify nodule-induced yet-to-be annotated peanut genes. Comparison between peanut, Medicago truncatula, Lotus japonicus, and Glycine max showed upregulation of 61 peanut orthologs among 111 tested known RNS-related genes, indicating conservation in mechanisms of nodule development among members of the Papilionoid family. Unlike model legumes, recruitment of class 1 phytoglobin-derived symbiotic hemoglobin (SymH) in peanut indicates diversification of oxygen-scavenging mechanisms in the Papilionoid family. Finally, the absence of cysteine-rich motif-1-containing nodule-specific cysteine-rich peptide (NCR) genes but the recruitment of defensin-like NCRs suggest a diverse molecular mechanism of terminal bacteroid differentiation. In summary, our work describes genetic conservation and diversification in legume-rhizobia symbiosis in the Papilionoid family, as well as among members of the Dalbergoid legumes.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Legacy genetics of Arachis cardenasii in the peanut crop shows the profound benefits of international seed exchange.

The narrow genetics of most crops is a fundamental vulnerability to food security. This makes wild crop relatives a strategic resource of genetic diversity that can be used for crop improvement and adaptation to new agricultural challenges. Here, we uncover the contribution of one wild species accession, Arachis cardenasii GKP 10017, to the peanut crop (Arachis hypogaea) that was initiated by complex hybridizations in the 1960s and propagated by international seed exchange. However, until this study, the global scale of the dispersal of genetic contributions from this wild accession had been obscured by the multiple germplasm transfers, breeding cycles, and unrecorded genetic mixing between lineages that had occurred over the years. By genetic analysis and pedigree research, we identified A. cardenasii-enhanced, disease-resistant cultivars in Africa, Asia, Oceania, and the Americas. These cultivars provide widespread improved food security and environmental and economic benefits. This study emphasizes the importance of wild species and collaborative networks of international expertise for crop improvement. However, it also highlights the consequences of the implementation of a patchwork of restrictive national laws and sea changes in attitudes regarding germplasm that followed in the wake of the Convention on Biological Diversity. Today, the botanical collections and multiple seed exchanges which enable benefits such as those revealed by this study are drastically reduced. The research reported here underscores the vital importance of ready access to germplasm in ensuring long-term world food security.

Transcriptome Profile Reveals Drought-Induced Genes Preferentially Expressed in Response to Water Deficit in Cultivated Peanut (Arachis hypogaea L.).

Cultivated peanut (Arachis hypogaea) is one of the most widely grown food legumes in the world, being valued for its high protein and unsaturated oil contents. Drought stress is one of the major constraints that limit peanut production. This study's objective was to identify the drought-responsive genes preferentially expressed under drought stress in different peanut genotypes. To accomplish this, four genotypes (drought tolerant: C76-16 and 587; drought susceptible: Tifrunner and 506) subjected to drought stress in a rainout shelter experiment were examined. Transcriptome sequencing analysis identified that all four genotypes shared a total of 2,457 differentially expressed genes (DEGs). A total of 139 enriched gene ontology terms consisting of 86 biological processes and 53 molecular functions, with defense response, reproductive process, and signaling pathways, were significantly enriched in the common DEGs. In addition, 3,576 DEGs were identified only in drought-tolerant lines in which a total of 74 gene ontology terms were identified, including 55 biological processes and 19 molecular functions, mainly related to protein modification process, pollination, and metabolic process. These terms were also found in shared genes in four genotypes, indicating that tolerant lines adjusted more related genes to respond to drought. Forty-three significantly enriched Kyoto Encyclopedia of Genes and Genomes pathways were also identified, and the most enriched pathways were those processes involved in metabolic pathways, biosynthesis of secondary metabolites, plant circadian rhythm, phenylpropanoid biosynthesis, and starch and sucrose metabolism. This research expands our current understanding of the mechanisms that facilitate peanut drought tolerance and shed light on breeding advanced peanut lines to combat drought stress.

Draft Genome Sequences of One Aspergillus parasiticus Isolate and Nine Aspergillus flavus Isolates with Varying Stress Tolerance and Aflatoxin Production.

Aspergillus flavus and Aspergillus parasiticus produce carcinogenic aflatoxins during crop infection, with extensive variations in production among isolates, ranging from atoxigenic to highly toxigenic. Here, we report draft genome sequences of one A. parasiticus isolate and nine A. flavus isolates from field environments for use in comparative, functional, and phylogenetic studies.

Two New Aspergillus flavus Reference Genomes Reveal a Large Insertion Potentially Contributing to Isolate Stress Tolerance and Aflatoxin Production.

Efforts in genome sequencing in the Aspergillus genus have led to the development of quality reference genomes for several important species including A. nidulans, A. fumigatus, and A. oryzae However, less progress has been made for A. flavus As part of the effort of the USDA-ARS Annual Aflatoxin Workshop Fungal Genome Project, the isolate NRRL3357 was sequenced and resulted in a scaffold-level genome released in 2005. Our goal has been biologically driven, focusing on two areas: isolate variation in aflatoxin production and drought stress exacerbating aflatoxin production by A. flavus Therefore, we developed two reference pseudomolecule genome assemblies derived from chromosome arms for two isolates: AF13, a MAT1-2, highly stress tolerant, and highly aflatoxigenic isolate; and NRRL3357, a MAT1-1, less stress tolerant, and moderate aflatoxin producer in comparison to AF13. Here, we report these two reference-grade assemblies for these isolates through a combination of PacBio long-read sequencing and optical mapping, and coupled them with comparative, functional, and phylogenetic analyses. This analysis resulted in the identification of 153 and 45 unique genes in AF13 and NRRL3357, respectively. We also confirmed the presence of a unique 310 Kb insertion in AF13 containing 60 genes. Analysis of this insertion revealed the presence of a bZIP transcription factor, named atfC, which may contribute to isolate pathogenicity and stress tolerance. Phylogenomic analyses comparing these and other available assemblies also suggest that the species complex of A. flavus is polyphyletic.

Comparison of SNP Calling Pipelines and NGS Platforms to Predict the Genomic Regions Harboring Candidate Genes for Nodulation in Cultivated Peanut.

Cultivated peanut (Arachis hypogaea L.) forms root nodules to enable a symbiotic relationship with rhizobia for biological nitrogen fixation. To understand the genetic factors of peanut nodulation, it is fundamental to genetically map and clone the genes involved in nodulation. For genetic mapping, high throughput genotyping with a large number of polymorphic markers is critical. In this study, two sets of sister recombinant inbred lines (RILs), each containing a nodulating (Nod+) and non-nodulating (Nod-) line, and their Nod+ parental lines were extensively genotyped. Several next generation sequencing (NGS) methods including target enrichment sequencing (TES), RNA-sequencing (RNA-seq), genotyping by sequencing (GBS), and the 48K Axiom Arachis2 SNP array, and various analysis pipelines were applied to identify single nucleotide polymorphisms (SNP) among the two sets of RILs and their parents. TES revealed the largest number of homozygous SNPs (15,947) between the original parental lines, followed by the Axiom Arachis2 SNP array (1,887), RNA-seq (1,633), and GBS (312). Among the five SNP analysis pipelines applied, the alignment to A/B genome followed by HAPLOSWEEP revealed the largest number of homozygous SNPs and highest concordance rate (79%) with the array. A total of 222 and 1,200 homozygous SNPs were polymorphic between the Nod+ and Nod- sister RILs and between their parents, respectively. A graphical genotype map of the sister RILs was constructed with these SNPs, which demonstrated the candidate genomic regions harboring genes controlling nodulation across the whole genome. Results of this study mainly provide the pros and cons of NGS and SNP genotyping platforms for genetic mapping in peanut, and also provide potential genetic resources to narrow down the genomic regions controlling peanut nodulation, which would lay the foundation for gene cloning and improvement of nitrogen fixation in peanut.

Identification of QTLs for resistance to leaf spots in cultivated peanut (Arachis hypogaea L.) through GWAS analysis.

KEY MESSAGE: Two QTLs on ChrB09 significantly associated with both early and late leaf spots were identified by genome-wide association study in the US peanut mini-core collection. Early leaf spot (ELS) and late leaf spot (LLS) are two serious peanut diseases in the USA, causing tens of millions of dollars of annual economic losses. However, the genetic factors underlying resistance to those diseases in peanuts have not been well-studied. We conducted a genome-wide association study for the two peanut diseases using Affymetrix version 2.0 SNP array with 120 genotypes mainly coming from the US peanut mini-core collection. A total of 46 quantitative trait loci (QTLs) were identified with phenotypic variation explained (PVE) from 10.19 to 24.11%, in which eighteen QTLs are for resistance to ELS and 28 QTLs for LLS. Among the 46 QTLs, there were four and two major QTLs with PVE higher than 16.99% for resistance ELS and LLS, respectively. Of the six major QTLs, five were located on the B sub-genome and only one was on the A sub-genome, which suggested that the B sub-genome has more potential resistance genomic regions than the A sub-genome. In addition, two genomic regions on chromosome B09 were found to provide significant resistance to both ELS and LLS. A total of 74 non-redundant genes were identified as resistance genes, among which, twelve candidate genes were in significant genomic regions including two candidate genes for both ELS and LLS, and other ten candidate genes for ELS. The QTLs and candidate genes obtained from this study will be useful to breed peanuts for resistances to the diseases.

Publisher Correction: A recombination bin-map identified a major QTL for resistance to Tomato Spotted Wilt Virus in peanut (Arachis hypogaea).

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

A recombination bin-map identified a major QTL for resistance to Tomato Spotted Wilt Virus in peanut (Arachis hypogaea).

Tomato spotted wilt virus (TSWV) is a devastating disease to peanut growers in the South-eastern region of the United States. Newly released peanut cultivars in recent years are crucial as they have some levels of resistance to TSWV. One mapping population of recombinant inbred line (RIL) used in this study was derived from peanut lines of SunOleic 97R and NC94022. A whole genome re-sequencing approach was used to sequence these two parents and 140 RILs. A recombination bin-based genetic map was constructed, with 5,816 bins and 20 linkage groups covering a total length of 2004 cM. Using this map, we identified three QTLs which were colocalized on chromosome A01. One QTL had the largest effect of 36.51% to the phenotypic variation and encompassed 89.5 Kb genomic region. This genome region had a cluster of genes, which code for chitinases, strictosidine synthase-like, and NBS-LRR proteins. SNPs linked to this QTL were used to develop Kompetitive allele specific PCR (KASP) markers, and the validated KASP markers showed expected segregation of alleles coming from resistant and susceptible parents within the population. Therefore, this bin-map and QTL associated with TSWV resistance made it possible for functional gene mapping, map-based cloning, and marker-assisted breeding. This study identified the highest number of SNP makers and demonstrated recombination bin-based map for QTL identification in peanut. The chitinase gene clusters and NBS-LRR disease resistance genes in this region suggest the possible involvement in peanut resistance to TSWV.