Insights into Land Plant Evolution Garnered from the Marchantia polymorpha Genome.

The evolution of land flora transformed the terrestrial environment. Land plants evolved from an ancestral charophycean alga from which they inherited developmental, biochemical, and cell biological attributes. Additional biochemical and physiological adaptations to land, and a life cycle with an alternation between multicellular haploid and diploid generations that facilitated efficient dispersal of desiccation tolerant spores, evolved in the ancestral land plant. We analyzed the genome of the liverwort Marchantia polymorpha, a member of a basal land plant lineage. Relative to charophycean algae, land plant genomes are characterized by genes encoding novel biochemical pathways, new phytohormone signaling pathways (notably auxin), expanded repertoires of signaling pathways, and increased diversity in some transcription factor families. Compared with other sequenced land plants, M. polymorpha exhibits low genetic redundancy in most regulatory pathways, with this portion of its genome resembling that predicted for the ancestral land plant. PAPERCLIP.

Evolving New Skeletal Traits by cis-Regulatory Changes in Bone Morphogenetic Proteins.

Changes in bone size and shape are defining features of many vertebrates. Here we use genetic crosses and comparative genomics to identify specific regulatory DNA alterations controlling skeletal evolution. Armor bone-size differences in sticklebacks map to a major effect locus overlapping BMP family member GDF6. Freshwater fish express more GDF6 due in part to a transposon insertion, and transgenic overexpression of GDF6 phenocopies evolutionary changes in armor-plate size. The human GDF6 locus also has undergone distinctive regulatory evolution, including complete loss of an enhancer that is otherwise highly conserved between chimps and other mammals. Functional tests show that the ancestral enhancer drives expression in hindlimbs but not forelimbs, in locations that have been specifically modified during the human transition to bipedalism. Both gain and loss of regulatory elements can localize BMP changes to specific anatomical locations, providing a flexible regulatory basis for evolving species-specific changes in skeletal form.

Long-read genome sequencing and variant reanalysis increase diagnostic yield in neurodevelopmental disorders.

Variant detection from long-read genome sequencing (lrGS) has proven to be more accurate and comprehensive than variant detection from short-read genome sequencing (srGS). However, the rate at which lrGS can increase molecular diagnostic yield for rare disease is not yet precisely characterized. We performed lrGS using Pacific Biosciences HiFi technology on 96 short-read-negative probands with rare diseases that were suspected to be genetic. We generated hg38-aligned variants and de novo phased genome assemblies, and subsequently annotated, filtered, and curated variants using clinical standards. New disease-relevant or potentially relevant genetic findings were identified in 16/96 (16.7%) probands, nine of which (8/96, ~9.4%) harbored pathogenic or likely pathogenic variants. Nine probands (~9.4%) had variants that were accurately called in both srGS and lrGS and represent changes to clinical interpretation, mostly from recently published gene-disease associations. Seven cases included variants that were only correctly interpreted in lrGS, including copy-number variants, an inversion, a mobile element insertion, two low-complexity repeat expansions, and a 1 bp deletion. While evidence for each of these variants is, in retrospect, visible in srGS, they were either not called within srGS data, were represented by calls with incorrect sizes or structures, or failed quality-control and filtration. Thus, while reanalysis of older srGS data clearly increases diagnostic yield, we find that lrGS allows for substantial additional yield (7/96, 7.3%) beyond srGS. We anticipate that as lrGS analysis improves, and as lrGS datasets grow allowing for better variant frequency annotation, the additional lrGS-only rare disease yield will grow over time.

Extraordinary preservation of gene collinearity over three hundred million years revealed in homosporous lycophytes.

Homosporous lycophytes (Lycopodiaceae) are a deeply diverged lineage in the plant tree of life, having split from heterosporous lycophytes (Selaginella and Isoetes) ~400 Mya. Compared to the heterosporous lineage, Lycopodiaceae has markedly larger genome sizes and remains the last major plant clade for which no chromosome-level assembly has been available. Here, we present chromosomal genome assemblies for two homosporous lycophyte species, the allotetraploid Huperzia asiatica and the diploid Diphasiastrum complanatum. Remarkably, despite that the two species diverged ~350 Mya, around 30% of the genes are still in syntenic blocks. Furthermore, both genomes had undergone independent whole genome duplications, and the resulting intragenomic syntenies have likewise been preserved relatively well. Such slow genome evolution over deep time is in stark contrast to heterosporous lycophytes and is correlated with a decelerated rate of nucleotide substitution. Together, the genomes of H. asiatica and D. complanatum not only fill a crucial gap in the plant genomic landscape but also highlight a potentially meaningful genomic contrast between homosporous and heterosporous species.

The Chlamydomonas Genome Project, version 6: Reference assemblies for mating-type plus and minus strains reveal extensive structural mutation in the laboratory.

Five versions of the Chlamydomonas reinhardtii reference genome have been produced over the last two decades. Here we present version 6, bringing significant advances in assembly quality and structural annotations. PacBio-based chromosome-level assemblies for two laboratory strains, CC-503 and CC-4532, provide resources for the plus and minus mating-type alleles. We corrected major misassemblies in previous versions and validated our assemblies via linkage analyses. Contiguity increased over ten-fold and >80% of filled gaps are within genes. We used Iso-Seq and deep RNA-seq datasets to improve structural annotations, and updated gene symbols and textual annotation of functionally characterized genes via extensive manual curation. We discovered that the cell wall-less classical reference strain CC-503 exhibits genomic instability potentially caused by deletion of the helicase RECQ3, with major structural mutations identified that affect >100 genes. We therefore present the CC-4532 assembly as the primary reference, although this strain also carries unique structural mutations and is experiencing rapid proliferation of a Gypsy retrotransposon. We expect all laboratory strains to harbor gene-disrupting mutations, which should be considered when interpreting and comparing experimental results. Collectively, the resources presented here herald a new era of Chlamydomonas genomics and will provide the foundation for continued research in this important reference organism.

Unlocking saponin biosynthesis in soapwort.

Soapwort (Saponaria officinalis) is a flowering plant from the Caryophyllaceae family with a long history of human use as a traditional source of soap. Its detergent properties are because of the production of polar compounds (saponins), of which the oleanane-based triterpenoid saponins, saponariosides A and B, are the major components. Soapwort saponins have anticancer properties and are also of interest as endosomal escape enhancers for targeted tumor therapies. Intriguingly, these saponins share common structural features with the vaccine adjuvant QS-21 and, thus, represent a potential alternative supply of saponin adjuvant precursors. Here, we sequence the S. officinalis genome and, through genome mining and combinatorial expression, identify 14 enzymes that complete the biosynthetic pathway to saponarioside B. These enzymes include a noncanonical cytosolic GH1 (glycoside hydrolase family 1) transglycosidase required for the addition of D-quinovose. Our results open avenues for accessing and engineering natural and new-to-nature pharmaceuticals, drug delivery agents and potential immunostimulants.

The barley pan-genome reveals the hidden legacy of mutation breeding.

Genetic diversity is key to crop improvement. Owing to pervasive genomic structural variation, a single reference genome assembly cannot capture the full complement of sequence diversity of a crop species (known as the 'pan-genome'1). Multiple high-quality sequence assemblies are an indispensable component of a pan-genome infrastructure. Barley (Hordeum vulgare L.) is an important cereal crop with a long history of cultivation that is adapted to a wide range of agro-climatic conditions2. Here we report the construction of chromosome-scale sequence assemblies for the genotypes of 20 varieties of barley-comprising landraces, cultivars and a wild barley-that were selected as representatives of global barley diversity. We catalogued genomic presence/absence variants and explored the use of structural variants for quantitative genetic analysis through whole-genome shotgun sequencing of 300 gene bank accessions. We discovered abundant large inversion polymorphisms and analysed in detail two inversions that are frequently found in current elite barley germplasm; one is probably the product of mutation breeding and the other is tightly linked to a locus that is involved in the expansion of geographical range. This first-generation barley pan-genome makes previously hidden genetic variation accessible to genetic studies and breeding.

Assembly, comparative analysis, and utilization of a single haplotype reference genome for soybean.

Cultivar Williams 82 has served as the reference genome for the soybean research community since 2008, but is known to have areas of genomic heterogeneity among different sub-lines. This work provides an updated assembly (version Wm82.a6) derived from a specific sub-line known as Wm82-ISU-01 (seeds available under USDA accession PI 704477). The genome was assembled using Pacific BioSciences HiFi reads and integrated into chromosomes using HiC. The 20 soybean chromosomes assembled into a genome of 1.01Gb, consisting of 36 contigs. The genome annotation identified 48 387 gene models, named in accordance with previous assembly versions Wm82.a2 and Wm82.a4. Comparisons of Wm82.a6 with other near-gapless assemblies of Williams 82 reveal large regions of genomic heterogeneity, including regions of differential introgression from the cultivar Kingwa within approximately 30 Mb and 25 Mb segments on chromosomes 03 and 07, respectively. Additionally, our analysis revealed a previously unknown large (>20 Mb) heterogeneous region in the pericentromeric region of chromosome 12, where Wm82.a6 matches the 'Williams' haplotype while the other two near-gapless assemblies do not match the haplotype of either parent of Williams 82. In addition to the Wm82.a6 assembly, we also assembled the genome of 'Fiskeby III,' a rich resource for abiotic stress resistance genes. A genome comparison of Wm82.a6 with Fiskeby III revealed the nucleotide and structural polymorphisms between the two genomes within a QTL region for iron deficiency chlorosis resistance. The Wm82.a6 and Fiskeby III genomes described here will enhance comparative and functional genomics capacities and applications in the soybean community.

The western redcedar genome reveals low genetic diversity in a self-compatible conifer.

We assembled the 9.8-Gbp genome of western redcedar (WRC; Thuja plicata), an ecologically and economically important conifer species of the Cupressaceae. The genome assembly, derived from a uniquely inbred tree produced through five generations of self-fertilization (selfing), was determined to be 86% complete by BUSCO analysis, one of the most complete genome assemblies for a conifer. Population genomic analysis revealed WRC to be one of the most genetically depauperate wild plant species, with an effective population size of approximately 300 and no significant genetic differentiation across its geographic range. Nucleotide diversity, π, is low for a continuous tree species, with many loci showing zero diversity, and the ratio of π at zero- to fourfold degenerate sites is relatively high (approximately 0.33), suggestive of weak purifying selection. Using an array of genetic lines derived from up to five generations of selfing, we explored the relationship between genetic diversity and mating system. Although overall heterozygosity was found to decline faster than expected during selfing, heterozygosity persisted at many loci, and nearly 100 loci were found to deviate from expectations of genetic drift, suggestive of associative overdominance. Nonreference alleles at such loci often harbor deleterious mutations and are rare in natural populations, implying that balanced polymorphisms are maintained by linkage to dominant beneficial alleles. This may account for how WRC remains responsive to natural and artificial selection, despite low genetic diversity.

JGI Plant Gene Atlas: an updateable transcriptome resource to improve functional gene descriptions across the plant kingdom.

Gene functional descriptions offer a crucial line of evidence for candidate genes underlying trait variation. Conversely, plant responses to environmental cues represent important resources to decipher gene function and subsequently provide molecular targets for plant improvement through gene editing. However, biological roles of large proportions of genes across the plant phylogeny are poorly annotated. Here we describe the Joint Genome Institute (JGI) Plant Gene Atlas, an updateable data resource consisting of transcript abundance assays spanning 18 diverse species. To integrate across these diverse genotypes, we analyzed expression profiles, built gene clusters that exhibited tissue/condition specific expression, and tested for transcriptional response to environmental queues. We discovered extensive phylogenetically constrained and condition-specific expression profiles for genes without any previously documented functional annotation. Such conserved expression patterns and tightly co-expressed gene clusters let us assign expression derived additional biological information to 64 495 genes with otherwise unknown functions. The ever-expanding Gene Atlas resource is available at JGI Plant Gene Atlas (https://plantgeneatlas.jgi.doe.gov) and Phytozome (https://phytozome.jgi.doe.gov/), providing bulk access to data and user-specified queries of gene sets. Combined, these web interfaces let users access differentially expressed genes, track orthologs across the Gene Atlas plants, graphically represent co-expressed genes, and visualize gene ontology and pathway enrichments.

A generalist-specialist trade-off between switchgrass cytotypes impacts climate adaptation and geographic range.

Polyploidy results from whole-genome duplication and is a unique form of heritable variation with pronounced evolutionary implications. Different ploidy levels, or cytotypes, can exist within a single species, and such systems provide an opportunity to assess how ploidy variation alters phenotypic novelty, adaptability, and fitness, which can, in turn, drive the development of unique ecological niches that promote the coexistence of multiple cytotypes. Switchgrass, Panicum virgatum, is a widespread, perennial C4 grass in North America with multiple naturally occurring cytotypes, primarily tetraploids (4×) and octoploids (8×). Using a combination of genomic, quantitative genetic, landscape, and niche modeling approaches, we detect divergent levels of genetic admixture, evidence of niche differentiation, and differential environmental sensitivity between switchgrass cytotypes. Taken together, these findings support a generalist (8×)­specialist (4×) trade-off. Our results indicate that the 8× represent a unique combination of genetic variation that has allowed the expansion of switchgrass' ecological niche and thus putatively represents a valuable breeding resource.

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.

Transcriptome and DNA methylome divergence of inflorescence development between 2 ecotypes in Panicum hallii.

The morphological diversity of the inflorescence determines flower and seed production, which is critical for plant adaptation. Hall's panicgrass (Panicum hallii, P. hallii) is a wild perennial grass that has been developed as a model to study perennial grass biology and adaptive evolution. Highly divergent inflorescences have evolved between the 2 major ecotypes in P. hallii, the upland ecotype (P. hallii var hallii, HAL2 genotype) with compact inflorescence and large seed and the lowland ecotype (P. hallii var filipes, FIL2 genotype) with an open inflorescence and small seed. Here we conducted a comparative analysis of the transcriptome and DNA methylome, an epigenetic mark that influences gene expression regulation, across different stages of inflorescence development using genomic references for each ecotype. Global transcriptome analysis of differentially expressed genes (DEGs) and co-expression modules underlying the inflorescence divergence revealed the potential role of cytokinin signaling in heterochronic changes. Comparing DNA methylome profiles revealed a remarkable level of differential DNA methylation associated with the evolution of P. hallii inflorescence. We found that a large proportion of differentially methylated regions (DMRs) were located in the flanking regulatory regions of genes. Intriguingly, we observed a substantial bias of CHH hypermethylation in the promoters of FIL2 genes. The integration of DEGs, DMRs, and Ka/Ks ratio results characterized the evolutionary features of DMR-associated DEGs that contribute to the divergence of the P. hallii inflorescence. This study provides insights into the transcriptome and epigenetic landscape of inflorescence divergence in P. hallii and a genomic resource for perennial grass biology.

Long-read sequence assembly: a technical evaluation in barley.

Sequence assembly of large and repeat-rich plant genomes has been challenging, requiring substantial computational resources and often several complementary sequence assembly and genome mapping approaches. The recent development of fast and accurate long-read sequencing by circular consensus sequencing (CCS) on the PacBio platform may greatly increase the scope of plant pan-genome projects. Here, we compare current long-read sequencing platforms regarding their ability to rapidly generate contiguous sequence assemblies in pan-genome studies of barley (Hordeum vulgare). Most long-read assemblies are clearly superior to the current barley reference sequence based on short-reads. Assemblies derived from accurate long reads excel in most metrics, but the CCS approach was the most cost-effective strategy for assembling tens of barley genomes. A downsampling analysis indicated that 20-fold CCS coverage can yield very good sequence assemblies, while even five-fold CCS data may capture the complete sequence of most genes. We present an updated reference genome assembly for barley with near-complete representation of the repeat-rich intergenic space. Long-read assembly can underpin the construction of accurate and complete sequences of multiple genomes of a species to build pan-genome infrastructures in Triticeae crops and their wild relatives.

Gene disruption by structural mutations drives selection in US rice breeding over the last century.

The genetic basis of general plant vigor is of major interest to food producers, yet the trait is recalcitrant to genetic mapping because of the number of loci involved, their small effects, and linkage. Observations of heterosis in many crops suggests that recessive, malfunctioning versions of genes are a major cause of poor performance, yet we have little information on the mutational spectrum underlying these disruptions. To address this question, we generated a long-read assembly of a tropical japonica rice (Oryza sativa) variety, Carolina Gold, which allowed us to identify structural mutations (>50 bp) and orient them with respect to their ancestral state using the outgroup, Oryza glaberrima. Supporting prior work, we find substantial genome expansion in the sativa branch. While transposable elements (TEs) account for the largest share of size variation, the majority of events are not directly TE-mediated. Tandem duplications are the most common source of insertions and are highly enriched among 50-200bp mutations. To explore the relative impact of various mutational classes on crop fitness, we then track these structural events over the last century of US rice improvement using 101 resequenced varieties. Within this material, a pattern of temporary hybridization between medium and long-grain varieties was followed by recent divergence. During this long-term selection, structural mutations that impact gene exons have been removed at a greater rate than intronic indels and single-nucleotide mutations. These results support the use of ab initio estimates of mutational burden, based on structural data, as an orthogonal predictor in genomic selection.

Widespread polycistronic gene expression in green algae.

Polycistronic gene expression, common in prokaryotes, was thought to be extremely rare in eukaryotes. The development of long-read sequencing of full-length transcript isomers (Iso-Seq) has facilitated a reexamination of that dogma. Using Iso-Seq, we discovered hundreds of examples of polycistronic expression of nuclear genes in two divergent species of green algae: Chlamydomonas reinhardtii and Chromochloris zofingiensis Here, we employ a range of independent approaches to validate that multiple proteins are translated from a common transcript for hundreds of loci. A chromatin immunoprecipitation analysis using trimethylation of lysine 4 on histone H3 marks confirmed that transcription begins exclusively at the upstream gene. Quantification of polyadenylated [poly(A)] tails and poly(A) signal sequences confirmed that transcription ends exclusively after the downstream gene. Coexpression analysis found nearly perfect correlation for open reading frames (ORFs) within polycistronic loci, consistent with expression in a shared transcript. For many polycistronic loci, terminal peptides from both ORFs were identified from proteomics datasets, consistent with independent translation. Synthetic polycistronic gene pairs were transcribed and translated in vitro to recapitulate the production of two distinct proteins from a common transcript. The relative abundance of these two proteins can be modified by altering the Kozak-like sequence of the upstream gene. Replacement of the ORFs with selectable markers or reporters allows production of such heterologous proteins, speaking to utility in synthetic biology approaches. Conservation of a significant number of polycistronic gene pairs between C. reinhardtii, C. zofingiensis, and five other species suggests that this mechanism may be evolutionarily ancient and biologically important in the green algal lineage.

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.

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.

Phylogenomic structure and speciation in an emerging model: the Sphagnum magellanicum complex (Bryophyta).

Sphagnum magellanicum is one of two Sphagnum species for which a reference-quality genome exists to facilitate research in ecological genomics. Phylogenetic and comparative genomic analyses were conducted based on resequencing data from 48 samples and RADseq analyses based on 187 samples. We report herein that there are four clades/species within the S. magellanicum complex in eastern North America and that the reference genome belongs to Sphagnum divinum. The species exhibit tens of thousands (RADseq) to millions (resequencing) of fixed nucleotide differences. Two species, however, referred to informally as S. diabolicum and S. magni because they have not been formally described, are differentiated by only 100 (RADseq) to 1000 (resequencing) of differences. Introgression among species in the complex is demonstrated using D-statistics and f4 ratios. One ecologically important functional trait, tissue decomposability, which underlies peat (carbon) accumulation, does not differ between segregates in the S. magellanicum complex, although previous research showed that many closely related Sphagnum species have evolved differences in decomposability/carbon sequestration. Phylogenetic resolution and more accurate species delimitation in the S. magellanicum complex substantially increase the value of this group for studying the early evolutionary stages of climate adaptation and ecological evolution more broadly.