Major impacts of widespread structural variation on sorghum.

Genetic diversity is critical to crop breeding and improvement, and dissection of the genomic variation underlying agronomic traits can both assist breeding and give insight into basic biological mechanisms. Although recent genome analyses in plants reveal many structural variants (SVs), most current studies of crop genetic variation are dominated by single-nucleotide polymorphisms (SNPs). The extent of the impact of SVs on global trait variation, as well as their utility in genome-wide selection, is not yet understood. In this study, we built an SV data set based on whole-genome resequencing of diverse sorghum lines (n = 363), validated the correlation of photoperiod sensitivity and variety type, and identified SV hotspots underlying the divergent evolution of cellulosic and sweet sorghum. In addition, we showed the complementary contribution of SVs for heritability of traits related to sorghum adaptation. Importantly, inclusion of SV polymorphisms in association studies revealed genotype-phenotype associations not observed with SNPs alone. Three-way genome-wide association studies (GWAS) based on whole-genome SNP, SV, and integrated SNP + SV data sets showed substantial associations between SVs and sorghum traits. The addition of SVs to GWAS substantially increased heritability estimates for some traits, indicating their important contribution to functional allelic variation at the genome level. Our discovery of the widespread impacts of SVs on heritable gene expression variation could render a plausible mechanism for their disproportionate impact on phenotypic variation. This study expands our knowledge of SVs and emphasizes the extensive impacts of SVs on sorghum.

Haplotype-resolved genome assembly of Populus tremula × P. alba reveals aspen-specific megabase satellite DNA.

Populus species play a foundational role in diverse ecosystems and are important renewable feedstocks for bioenergy and bioproducts. Hybrid aspen Populus tremula × P. alba INRA 717-1B4 is a widely used transformation model in tree functional genomics and biotechnology research. As an outcrossing interspecific hybrid, its genome is riddled with sequence polymorphisms which present a challenge for sequence-sensitive analyses. Here we report a telomere-to-telomere genome for this hybrid aspen with two chromosome-scale, haplotype-resolved assemblies. We performed a comprehensive analysis of the repetitive landscape and identified both tandem repeat array-based and array-less centromeres. Unexpectedly, the most abundant satellite repeats in both haplotypes lie outside of the centromeres, consist of a 147 bp monomer PtaM147, frequently span >1 megabases, and form heterochromatic knobs. PtaM147 repeats are detected exclusively in aspens (section Populus) but PtaM147-like sequences occur in LTR-retrotransposons of closely related species, suggesting their origin from the retrotransposons. The genomic resource generated for this transformation model genotype has greatly improved the design and analysis of genome editing experiments that are highly sensitive to sequence polymorphisms. The work should motivate future hypothesis-driven research to probe into the function of the abundant and aspen-specific PtaM147 satellite DNA.

DNA hypomethylation of the host tree impairs interaction with mutualistic ectomycorrhizal fungus.

Ectomycorrhizas are an intrinsic component of tree nutrition and responses to environmental variations. How epigenetic mechanisms might regulate these mutualistic interactions is unknown. By manipulating the level of expression of the chromatin remodeler DECREASE IN DNA METHYLATION 1 (DDM1) and two demethylases DEMETER-LIKE (DML) in Populus tremula × Populus alba lines, we examined how host DNA methylation modulates multiple parameters of the responses to root colonization with the mutualistic fungus Laccaria bicolor. We compared the ectomycorrhizas formed between transgenic and wild-type (WT) trees and analyzed their methylomes and transcriptomes. The poplar lines displaying lower mycorrhiza formation rate corresponded to hypomethylated overexpressing DML or RNAi-ddm1 lines. We found 86 genes and 288 transposable elements (TEs) differentially methylated between WT and hypomethylated lines (common to both OX-dml and RNAi-ddm1) and 120 genes/1441 TEs in the fungal genome suggesting a host-induced remodeling of the fungal methylome. Hypomethylated poplar lines displayed 205 differentially expressed genes (cis and trans effects) in common with 17 being differentially methylated (cis). Our findings suggest a central role of host and fungal DNA methylation in the ability to form ectomycorrhizas including not only poplar genes involved in root initiation, ethylene and jasmonate-mediated pathways, and immune response but also terpenoid metabolism.

High-resolution mapping reveals hotspots and sex-biased recombination in Populus trichocarpa.

Fine-scale meiotic recombination is fundamental to the outcome of natural and artificial selection. Here, dense genetic mapping and haplotype reconstruction were used to estimate recombination for a full factorial Populus trichocarpa cross of 7 males and 7 females. Genomes of the resulting 49 full-sib families (N = 829 offspring) were resequenced, and high-fidelity biallelic SNP/INDELs and pedigree information were used to ascertain allelic phase and impute progeny genotypes to recover gametic haplotypes. The 14 parental genetic maps contained 1,820 SNP/INDELs on average that covered 376.7 Mb of physical length across 19 chromosomes. Comparison of parental and progeny haplotypes allowed fine-scale demarcation of cross-over regions, where 38,846 cross-over events in 1,658 gametes were observed. Cross-over events were positively associated with gene density and negatively associated with GC content and long-terminal repeats. One of the most striking findings was higher rates of cross-overs in males in 8 out of 19 chromosomes. Regions with elevated male cross-over rates had lower gene density and GC content than windows showing no sex bias. High-resolution analysis identified 67 candidate cross-over hotspots spread throughout the genome. DNA sequence motifs enriched in these regions showed striking similarity to those of maize, Arabidopsis, and wheat. These findings, and recombination estimates, will be useful for ongoing efforts to accelerate domestication of this and other biomass feedstocks, as well as future studies investigating broader questions related to evolutionary history, perennial development, phenology, wood formation, vegetative propagation, and dioecy that cannot be studied using annual plant model systems.

Salinity tolerance mechanisms of an Arctic Pelagophyte using comparative transcriptomic and gene expression analysis.

Little is known at the transcriptional level about microbial eukaryotic adaptations to short-term salinity change. Arctic microalgae are exposed to low salinity due to sea-ice melt and higher salinity with brine channel formation during freeze-up. Here, we investigate the transcriptional response of an ice-associated microalgae over salinities from 45 to 8. Our results show a bracketed response of differential gene expression when the cultures were exposed to progressively decreasing salinity. Key genes associated with salinity changes were involved in specific metabolic pathways, transcription factors and regulators, protein kinases, carbohydrate active enzymes, and inorganic ion transporters. The pelagophyte seemed to use a strategy involving overexpression of Na+-H+ antiporters and Na+ -Pi symporters as salinity decreases, but the K+ channel complex at higher salinities. Specific adaptation to cold saline arctic conditions was seen with differential expression of several antifreeze proteins, an ice-binding protein and an acyl-esterase involved in cold adaptation.

Multiplex knockout of trichome-regulating MYB duplicates in hybrid poplar using a single gRNA.

As the focus for CRISPR/Cas-edited plants moves from proof-of-concept to real-world applications, precise gene manipulation will increasingly require concurrent multiplex editing for polygenic traits. A common approach for editing across multiple sites is to design one guide RNA (gRNA) per target; however, this complicates construct assembly and increases the possibility of off-target mutations. In this study, we utilized one gRNA to target MYB186, a known positive trichome regulator, as well as its paralogs MYB138 and MYB38 at a consensus site for mutagenesis in hybrid poplar (Populus tremula × P. alba INRA 717-1B4). Unexpected duplications of MYB186 and MYB138 resulted in eight alleles for the three targeted genes in the hybrid poplar. Deep sequencing and polymerase chain reaction analyses confirmed editing across all eight targets in nearly all of the resultant glabrous mutants, ranging from small indels to large genomic dropouts, with no off-target activity detected at four potential sites. This highlights the effectiveness of a single gRNA targeting conserved exonic regions for multiplex editing. Additionally, cuticular wax and whole-leaf analyses showed a complete absence of triterpenes in the trichomeless mutants, hinting at a previously undescribed role for the nonglandular trichomes of poplar.

Evolutionary transition to the ectomycorrhizal habit in the genomes of a hyperdiverse lineage of mushroom-forming fungi.

The ectomycorrhizal (ECM) symbiosis has independently evolved from diverse types of saprotrophic ancestors. In this study, we seek to identify genomic signatures of the transition to the ECM habit within the hyperdiverse Russulaceae. We present comparative analyses of the genomic architecture and the total and secreted gene repertoires of 18 species across the order Russulales, of which 13 are newly sequenced, including a representative of a saprotrophic member of Russulaceae, Gloeopeniophorella convolvens. The genomes of ECM Russulaceae are characterized by a loss of genes for plant cell wall-degrading enzymes (PCWDEs), an expansion of genome size through increased transposable element (TE) content, a reduction in secondary metabolism clusters, and an association of small secreted proteins (SSPs) with TE 'nests', or dense aggregations of TEs. Some PCWDEs have been retained or even expanded, mostly in a species-specific manner. The genome of G. convolvens possesses some characteristics of ECM genomes (e.g. loss of some PCWDEs, TE expansion, reduction in secondary metabolism clusters). Functional specialization in ECM decomposition may drive diversification. Accelerated gene evolution predates the evolution of the ECM habit, indicating that changes in genome architecture and gene content may be necessary to prime the evolutionary switch.

Genomic Signatures of a Major Adaptive Event in the Pathogenic Fungus Melampsora larici-populina.

The recent availability of genome-wide sequencing techniques has allowed systematic screening for molecular signatures of adaptation, including in nonmodel organisms. Host-pathogen interactions constitute good models due to the strong selective pressures that they entail. We focused on an adaptive event which affected the poplar rust fungus Melampsora larici-populina when it overcame a resistance gene borne by its host, cultivated poplar. Based on 76 virulent and avirulent isolates framing narrowly the estimated date of the adaptive event, we examined the molecular signatures of selection. Using an array of genome scan methods based on different features of nucleotide diversity, we detected a single locus exhibiting a consistent pattern suggestive of a selective sweep in virulent individuals (excess of differentiation between virulent and avirulent samples, linkage disequilibrium, genotype-phenotype statistical association, and long-range haplotypes). Our study pinpoints a single gene and further a single amino acid replacement which may have allowed the adaptive event. Although our samples are nearly contemporary to the selective sweep, it does not seem to have affected genome diversity further than the immediate vicinity of the causal locus, which can be explained by a soft selective sweep (where selection acts on standing variation) and by the impact of recombination in mitigating the impact of selection. Therefore, it seems that properties of the life cycle of M. larici-populina, which entails both high genetic diversity and outbreeding, has facilitated its adaptation.

Genetic determinants of endophytism in the Arabidopsis root mycobiome.

The roots of Arabidopsis thaliana host diverse fungal communities that affect plant health and disease states. Here, we sequence the genomes of 41 fungal isolates representative of the A. thaliana root mycobiota for comparative analysis with other 79 plant-associated fungi. Our analyses indicate that root mycobiota members evolved from ancestors with diverse lifestyles and retain large repertoires of plant cell wall-degrading enzymes (PCWDEs) and effector-like small secreted proteins. We identify a set of 84 gene families associated with endophytism, including genes encoding PCWDEs acting on xylan (family GH10) and cellulose (family AA9). Transcripts encoding these enzymes are also part of a conserved transcriptional program activated by phylogenetically-distant mycobiota members upon host contact. Recolonization experiments with individual fungi indicate that strains with detrimental effects in mono-association with the host colonize roots more aggressively than those with beneficial activities, and dominate in natural root samples. Furthermore, we show that the pectin-degrading enzyme family PL1_7 links aggressiveness of endophytic colonization to plant health.

Pests, diseases, and aridity have shaped the genome of Corymbia citriodora.

Corymbia citriodora is a member of the predominantly Southern Hemisphere Myrtaceae family, which includes the eucalypts (Eucalyptus, Corymbia and Angophora; ~800 species). Corymbia is grown for timber, pulp and paper, and essential oils in Australia, South Africa, Asia, and Brazil, maintaining a high-growth rate under marginal conditions due to drought, poor-quality soil, and biotic stresses. To dissect the genetic basis of these desirable traits, we sequenced and assembled the 408 Mb genome of Corymbia citriodora, anchored into eleven chromosomes. Comparative analysis with Eucalyptus grandis reveals high synteny, although the two diverged approximately 60 million years ago and have different genome sizes (408 vs 641 Mb), with few large intra-chromosomal rearrangements. C. citriodora shares an ancient whole-genome duplication event with E. grandis but has undergone tandem gene family expansions related to terpene biosynthesis, innate pathogen resistance, and leaf wax formation, enabling their successful adaptation to biotic/abiotic stresses and arid conditions of the Australian continent.

Experimentally Validated Reconstruction and Analysis of a Genome-Scale Metabolic Model of an Anaerobic Neocallimastigomycota Fungus.

Anaerobic gut fungi in the phylum Neocallimastigomycota typically inhabit the digestive tracts of large mammalian herbivores, where they play an integral role in the decomposition of raw lignocellulose into its constitutive sugar monomers. However, quantitative tools to study their physiology are lacking, partially due to their complex and unresolved metabolism that includes the largely uncharacterized fungal hydrogenosome. Modern omics approaches combined with metabolic modeling can be used to establish an understanding of gut fungal metabolism and develop targeted engineering strategies to harness their degradation capabilities for lignocellulosic bioprocessing. Here, we introduce a high-quality genome of the anaerobic fungus Neocallimastix lanati from which we constructed the first genome-scale metabolic model of an anaerobic fungus. Relative to its size (200 Mbp, sequenced at 62× depth), it is the least fragmented publicly available gut fungal genome to date. Of the 1,788 lignocellulolytic enzymes annotated in the genome, 585 are associated with the fungal cellulosome, underscoring the powerful lignocellulolytic potential of N. lanati The genome-scale metabolic model captures the primary metabolism of N. lanati and accurately predicts experimentally validated substrate utilization requirements. Additionally, metabolic flux predictions are verified by 13C metabolic flux analysis, demonstrating that the model faithfully describes the underlying fungal metabolism. Furthermore, the model clarifies key aspects of the hydrogenosomal metabolism and can be used as a platform to quantitatively study these biotechnologically important yet poorly understood early-branching fungi.IMPORTANCE Recent genomic analyses have revealed that anaerobic gut fungi possess both the largest number and highest diversity of lignocellulolytic enzymes of all sequenced fungi, explaining their ability to decompose lignocellulosic substrates, e.g., agricultural waste, into fermentable sugars. Despite their potential, the development of engineering methods for these organisms has been slow due to their complex life cycle, understudied metabolism, and challenging anaerobic culture requirements. Currently, there is no framework that can be used to combine multi-omic data sets to understand their physiology. Here, we introduce a high-quality PacBio-sequenced genome of the anaerobic gut fungus Neocallimastix lanati Beyond identifying a trove of lignocellulolytic enzymes, we use this genome to construct the first genome-scale metabolic model of an anaerobic gut fungus. The model is experimentally validated and sheds light on unresolved metabolic features common to gut fungi. Model-guided analysis will pave the way for deepening our understanding of anaerobic gut fungi and provides a systematic framework to guide strain engineering efforts of these organisms for biotechnological use.

Desert truffle genomes reveal their reproductive modes and new insights into plant-fungal interaction and ectendomycorrhizal lifestyle.

Desert truffles are edible hypogeous fungi forming ectendomycorrhizal symbiosis with plants of Cistaceae family. Knowledge about the reproductive modes of these fungi and the molecular mechanisms driving the ectendomycorrhizal interaction is lacking. Genomes of the highly appreciated edible desert truffles Terfezia claveryi Chatin and Tirmania nivea Trappe have been sequenced and compared with other Pezizomycetes. Transcriptomes of T. claveryi × Helianthemum almeriense mycorrhiza from well-watered and drought-stressed plants, when intracellular colonizations is promoted, were investigated. We have identified the fungal genes related to sexual reproduction in desert truffles and desert-truffles-specific genomic and secretomic features with respect to other Pezizomycetes, such as the expansion of a large set of gene families with unknown Pfam domains and a number of species or desert-truffle-specific small secreted proteins differentially regulated in symbiosis. A core set of plant genes, including carbohydrate, lipid-metabolism, and defence-related genes, differentially expressed in mycorrhiza under both conditions was found. Our results highlight the singularities of desert truffles with respect to other mycorrhizal fungi while providing a first glimpse on plant and fungal determinants involved in ecto to endo symbiotic switch that occurs in desert truffle under dry conditions.

Large-scale genome sequencing of mycorrhizal fungi provides insights into the early evolution of symbiotic traits.

Mycorrhizal fungi are mutualists that play crucial roles in nutrient acquisition in terrestrial ecosystems. Mycorrhizal symbioses arose repeatedly across multiple lineages of Mucoromycotina, Ascomycota, and Basidiomycota. Considerable variation exists in the capacity of mycorrhizal fungi to acquire carbon from soil organic matter. Here, we present a combined analysis of 135 fungal genomes from 73 saprotrophic, endophytic and pathogenic species, and 62 mycorrhizal species, including 29 new mycorrhizal genomes. This study samples ecologically dominant fungal guilds for which there were previously no symbiotic genomes available, including ectomycorrhizal Russulales, Thelephorales and Cantharellales. Our analyses show that transitions from saprotrophy to symbiosis involve (1) widespread losses of degrading enzymes acting on lignin and cellulose, (2) co-option of genes present in saprotrophic ancestors to fulfill new symbiotic functions, (3) diversification of novel, lineage-specific symbiosis-induced genes, (4) proliferation of transposable elements and (5) divergent genetic innovations underlying the convergent origins of the ectomycorrhizal guild.

Succession of physiological stages hallmarks the transcriptomic response of the fungus Aspergillus niger to lignocellulose.

BACKGROUND: Understanding how fungi degrade lignocellulose is a cornerstone of improving renewables-based biotechnology, in particular for the production of hydrolytic enzymes. Considerable progress has been made in investigating fungal degradation during time-points where CAZyme expression peaks. However, a robust understanding of the fungal survival strategies over its life time on lignocellulose is thereby missed. Here we aimed to uncover the physiological responses of the biotechnological workhorse and enzyme producer Aspergillus niger over its life time to six substrates important for biofuel production. RESULTS: We analysed the response of A. niger to the feedstock Miscanthus and compared it with our previous study on wheat straw, alone or in combination with hydrothermal or ionic liquid feedstock pretreatments. Conserved (substrate-independent) metabolic responses as well as those affected by pretreatment and feedstock were identified via multivariate analysis of genome-wide transcriptomics combined with targeted transcript and protein analyses and mapping to a metabolic model. Initial exposure to all substrates increased fatty acid beta-oxidation and lipid metabolism transcripts. In a strain carrying a deletion of the ortholog of the Aspergillus nidulans fatty acid beta-oxidation transcriptional regulator farA, there was a reduction in expression of selected lignocellulose degradative CAZyme-encoding genes suggesting that beta-oxidation contributes to adaptation to lignocellulose. Mannan degradation expression was wheat straw feedstock-dependent and pectin degradation was higher on the untreated substrates. In the later life stages, known and novel secondary metabolite gene clusters were activated, which are of high interest due to their potential to synthesize bioactive compounds. CONCLUSION: In this study, which includes the first transcriptional response of Aspergilli to Miscanthus, we highlighted that life time as well as substrate composition and structure (via variations in pretreatment and feedstock) influence the fungal responses to lignocellulose. We also demonstrated that the fungal response contains physiological stages that are conserved across substrates and are typically found outside of the conditions with high CAZyme expression, as exemplified by the stages that are dominated by lipid and secondary metabolism.

Marker-free carotenoid-enriched rice generated through targeted gene insertion using CRISPR-Cas9.

Targeted insertion of transgenes at pre-determined plant genomic safe harbors provides a desirable alternative to insertions at random sites achieved through conventional methods. Most existing cases of targeted gene insertion in plants have either relied on the presence of a selectable marker gene in the insertion cassette or occurred at low frequency with relatively small DNA fragments (<1.8 kb). Here, we report the use of an optimized CRISPR-Cas9-based method to achieve the targeted insertion of a 5.2 kb carotenoid biosynthesis cassette at two genomic safe harbors in rice. We obtain marker-free rice plants with high carotenoid content in the seeds and no detectable penalty in morphology or yield. Whole-genome sequencing reveals the absence of off-target mutations by Cas9 in the engineered plants. These results demonstrate targeted gene insertion of marker-free DNA in rice using CRISPR-Cas9 genome editing, and offer a promising strategy for genetic improvement of rice and other crops.

Modulation of Plant and Fungal Gene Expression Upon Cd Exposure and Symbiosis in Ericoid Mycorrhizal Vaccinium myrtillus.

The success of Ericaceae in stressful habitats enriched in heavy metals has been ascribed to the distinctive abilities of their mycorrhizal fungal partners to withstand heavy metal stress and to enhance metal tolerance in the host plant. Whereas heavy metal tolerance has been extensively investigated in some ericoid mycorrhizal (ERM) fungi, the molecular and cellular mechanisms that extend tolerance to the host plant are currently unknown. Here, we show a reduced Cd content in Cd-exposed mycorrhizal roots of Vaccinium myrtillus colonized by a metal tolerant isolate of the fungus Oidiodendron maius as compared to non-mycorrhizal roots. To better understand this phenotype, we applied Next Generation Sequencing technologies to analyze gene expression in V. myrtillus and O. maius Zn grown under normal and Cd-stressed conditions, in the free living and in the mycorrhizal status. The results clearly showed that Cd had a stronger impact on plant gene expression than symbiosis, whereas fungal gene expression was mainly regulated by symbiosis. The higher abundance of transcripts coding for stress related proteins in non-mycorrhizal roots may be related to the higher Cd content. Regulated plant metal transporters have been identified that may play a role in reducing Cd content in mycorrhizal roots exposed to this metal.

An ectomycorrhizal fungus alters sensitivity to jasmonate, salicylate, gibberellin, and ethylene in host roots.

The phytohormones jasmonate, gibberellin, salicylate, and ethylene regulate an interconnected reprogramming network integrating root development with plant responses against microbes. The establishment of mutualistic ectomycorrhizal symbiosis requires the suppression of plant defense responses against fungi as well as the modification of root architecture and cortical cell wall properties. Here, we investigated the contribution of phytohormones and their crosstalk to the ontogenesis of ectomycorrhizae (ECM) between grey poplar (Populus tremula x alba) roots and the fungus Laccaria bicolor. To obtain the hormonal blueprint of developing ECM, we quantified the concentrations of jasmonates, gibberellins, and salicylate via liquid chromatography-tandem mass spectrometry. Subsequently, we assessed root architecture, mycorrhizal morphology, and gene expression levels (RNA sequencing) in phytohormone-treated poplar lateral roots in the presence or absence of L. bicolor. Salicylic acid accumulated in mid-stage ECM. Exogenous phytohormone treatment affected the fungal colonization rate and/or frequency of Hartig net formation. Colonized lateral roots displayed diminished responsiveness to jasmonate but regulated some genes, implicated in defense and cell wall remodelling, that were specifically differentially expressed after jasmonate treatment. Responses to salicylate, gibberellin, and ethylene were enhanced in ECM. The dynamics of phytohormone accumulation and response suggest that jasmonate, gibberellin, salicylate, and ethylene signalling play multifaceted roles in poplar L. bicolor ectomycorrhizal development.

Genome sequence of the model rice variety KitaakeX.

BACKGROUND: The availability of thousands of complete rice genome sequences from diverse varieties and accessions has laid the foundation for in-depth exploration of the rice genome. One drawback to these collections is that most of these rice varieties have long life cycles, and/or low transformation efficiencies, which limits their usefulness as model organisms for functional genomics studies. In contrast, the rice variety Kitaake has a rapid life cycle (9 weeks seed to seed) and is easy to transform and propagate. For these reasons, Kitaake has emerged as a model for studies of diverse monocotyledonous species. RESULTS: Here, we report the de novo genome sequencing and analysis of Oryza sativa ssp. japonica variety KitaakeX, a Kitaake plant carrying the rice XA21 immune receptor. Our KitaakeX sequence assembly contains 377.6 Mb, consisting of 33 scaffolds (476 contigs) with a contig N50 of 1.4 Mb. Complementing the assembly are detailed gene annotations of 35,594 protein coding genes. We identified 331,335 genomic variations between KitaakeX and Nipponbare (ssp. japonica), and 2,785,991 variations between KitaakeX and Zhenshan97 (ssp. indica). We also compared Kitaake resequencing reads to the KitaakeX assembly and identified 219 small variations. The high-quality genome of the model rice plant KitaakeX will accelerate rice functional genomics. CONCLUSIONS: The high quality, de novo assembly of the KitaakeX genome will serve as a useful reference genome for rice and will accelerate functional genomics studies of rice and other species.

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.

Megaphylogeny resolves global patterns of mushroom evolution.

Mushroom-forming fungi (Agaricomycetes) have the greatest morphological diversity and complexity of any group of fungi. They have radiated into most niches and fulfil diverse roles in the ecosystem, including wood decomposers, pathogens or mycorrhizal mutualists. Despite the importance of mushroom-forming fungi, large-scale patterns of their evolutionary history are poorly known, in part due to the lack of a comprehensive and dated molecular phylogeny. Here, using multigene and genome-based data, we assemble a 5,284-species phylogenetic tree and infer ages and broad patterns of speciation/extinction and morphological innovation in mushroom-forming fungi. Agaricomycetes started a rapid class-wide radiation in the Jurassic, coinciding with the spread of (sub)tropical coniferous forests and a warming climate. A possible mass extinction, several clade-specific adaptive radiations and morphological diversification of fruiting bodies followed during the Cretaceous and the Paleogene, convergently giving rise to the classic toadstool morphology, with a cap, stalk and gills (pileate-stipitate morphology). This morphology is associated with increased rates of lineage diversification, suggesting it represents a key innovation in the evolution of mushroom-forming fungi. The increase in mushroom diversity started during the Mesozoic-Cenozoic radiation event, an era of humid climate when terrestrial communities dominated by gymnosperms and reptiles were also expanding.