Phylogenomic discovery of deleterious mutations facilitates hybrid potato breeding.

Hybrid potato breeding will transform the crop from a clonally propagated tetraploid to a seed-reproducing diploid. Historical accumulation of deleterious mutations in potato genomes has hindered the development of elite inbred lines and hybrids. Utilizing a whole-genome phylogeny of 92 Solanaceae and its sister clade species, we employ an evolutionary strategy to identify deleterious mutations. The deep phylogeny reveals the genome-wide landscape of highly constrained sites, comprising ∼2.4% of the genome. Based on a diploid potato diversity panel, we infer 367,499 deleterious variants, of which 50% occur at non-coding and 15% at synonymous sites. Counterintuitively, diploid lines with relatively high homozygous deleterious burden can be better starting material for inbred-line development, despite showing less vigorous growth. Inclusion of inferred deleterious mutations increases genomic-prediction accuracy for yield by 24.7%. Our study generates insights into the genome-wide incidence and properties of deleterious mutations and their far-reaching consequences for breeding.

Genome design of hybrid potato.

Reinventing potato from a clonally propagated tetraploid into a seed-propagated diploid, hybrid potato, is an important innovation in agriculture. Due to deleterious mutations, it has remained a challenge to develop highly homozygous inbred lines, a prerequisite to breed hybrid potato. Here, we employed genome design to develop a generation of pure and fertile potato lines and thereby the uniform, vigorous F1s. The metrics we applied in genome design included the percentage of genome homozygosity and the number of deleterious mutations in the starting material, the number of segregation distortions in the S1 population, the haplotype information to infer the break of tight linkage between beneficial and deleterious alleles, and the genome complementarity of the parental lines. This study transforms potato breeding from a slow, non-accumulative mode into a fast-iterative one, thereby potentiating a broad spectrum of benefits to farmers and consumers.

Clonal Decomposition and DNA Replication States Defined by Scaled Single-Cell Genome Sequencing.

Accurate measurement of clonal genotypes, mutational processes, and replication states from individual tumor-cell genomes will facilitate improved understanding of tumor evolution. We have developed DLP+, a scalable single-cell whole-genome sequencing platform implemented using commodity instruments, image-based object recognition, and open source computational methods. Using DLP+, we have generated a resource of 51,926 single-cell genomes and matched cell images from diverse cell types including cell lines, xenografts, and diagnostic samples with limited material. From this resource we have defined variation in mitotic mis-segregation rates across tissue types and genotypes. Analysis of matched genomic and image measurements revealed correlations between cellular morphology and genome ploidy states. Aggregation of cells sharing copy number profiles allowed for calculation of single-nucleotide resolution clonal genotypes and inference of clonal phylogenies and avoided the limitations of bulk deconvolution. Finally, joint analysis over the above features defined clone-specific chromosomal aneuploidy in polyclonal populations.

Development of a Comprehensive Genotype-to-Fitness Map of Adaptation-Driving Mutations in Yeast.

Adaptive evolution plays a large role in generating the phenotypic diversity observed in nature, yet current methods are impractical for characterizing the molecular basis and fitness effects of large numbers of individual adaptive mutations. Here, we used a DNA barcoding approach to generate the genotype-to-fitness map for adaptation-driving mutations from a Saccharomyces cerevisiae population experimentally evolved by serial transfer under limiting glucose. We isolated and measured the fitness of thousands of independent adaptive clones and sequenced the genomes of hundreds of clones. We found only two major classes of adaptive mutations: self-diploidization and mutations in the nutrient-responsive Ras/PKA and TOR/Sch9 pathways. Our large sample size and precision of measurement allowed us to determine that there are significant differences in fitness between mutations in different genes, between different paralogs, and even between different classes of mutations within the same gene.

Stable Chromosome Condensation Revealed by Chromosome Conformation Capture.

Chemical cross-linking and DNA sequencing have revealed regions of intra-chromosomal interaction, referred to as topologically associating domains (TADs), interspersed with regions of little or no interaction, in interphase nuclei. We find that TADs and the regions between them correspond with the bands and interbands of polytene chromosomes of Drosophila. We further establish the conservation of TADs between polytene and diploid cells of Drosophila. From direct measurements on light micrographs of polytene chromosomes, we then deduce the states of chromatin folding in the diploid cell nucleus. Two states of folding, fully extended fibers containing regulatory regions and promoters, and fibers condensed up to 10-fold containing coding regions of active genes, constitute the euchromatin of the nuclear interior. Chromatin fibers condensed up to 30-fold, containing coding regions of inactive genes, represent the heterochromatin of the nuclear periphery. A convergence of molecular analysis with direct observation thus reveals the architecture of interphase chromosomes.

Population genomics of post-glacial western Eurasia.

Western Eurasia witnessed several large-scale human migrations during the Holocene1-5. Here, to investigate the cross-continental effects of these migrations, we shotgun-sequenced 317 genomes-mainly from the Mesolithic and Neolithic periods-from across northern and western Eurasia. These were imputed alongside published data to obtain diploid genotypes from more than 1,600 ancient humans. Our analyses revealed a 'great divide' genomic boundary extending from the Black Sea to the Baltic. Mesolithic hunter-gatherers were highly genetically differentiated east and west of this zone, and the effect of the neolithization was equally disparate. Large-scale ancestry shifts occurred in the west as farming was introduced, including near-total replacement of hunter-gatherers in many areas, whereas no substantial ancestry shifts happened east of the zone during the same period. Similarly, relatedness decreased in the west from the Neolithic transition onwards, whereas, east of the Urals, relatedness remained high until around 4,000 BP, consistent with the persistence of localized groups of hunter-gatherers. The boundary dissolved when Yamnaya-related ancestry spread across western Eurasia around 5,000 BP, resulting in a second major turnover that reached most parts of Europe within a 1,000-year span. The genetic origin and fate of the Yamnaya have remained elusive, but we show that hunter-gatherers from the Middle Don region contributed ancestry to them. Yamnaya groups later admixed with individuals associated with the Globular Amphora culture before expanding into Europe. Similar turnovers occurred in western Siberia, where we report new genomic data from a 'Neolithic steppe' cline spanning the Siberian forest steppe to Lake Baikal. These prehistoric migrations had profound and lasting effects on the genetic diversity of Eurasian populations.

A draft human pangenome reference.

Here the Human Pangenome Reference Consortium presents a first draft of the human pangenome reference. The pangenome contains 47 phased, diploid assemblies from a cohort of genetically diverse individuals1. These assemblies cover more than 99% of the expected sequence in each genome and are more than 99% accurate at the structural and base pair levels. Based on alignments of the assemblies, we generate a draft pangenome that captures known variants and haplotypes and reveals new alleles at structurally complex loci. We also add 119 million base pairs of euchromatic polymorphic sequences and 1,115 gene duplications relative to the existing reference GRCh38. Roughly 90 million of the additional base pairs are derived from structural variation. Using our draft pangenome to analyse short-read data reduced small variant discovery errors by 34% and increased the number of structural variants detected per haplotype by 104% compared with GRCh38-based workflows, which enabled the typing of the vast majority of structural variant alleles per sample.

The giant diploid faba genome unlocks variation in a global protein crop.

Increasing the proportion of locally produced plant protein in currently meat-rich diets could substantially reduce greenhouse gas emissions and loss of biodiversity1. However, plant protein production is hampered by the lack of a cool-season legume equivalent to soybean in agronomic value2. Faba bean (Vicia faba L.) has a high yield potential and is well suited for cultivation in temperate regions, but genomic resources are scarce. Here, we report a high-quality chromosome-scale assembly of the faba bean genome and show that it has expanded to a massive 13 Gb in size through an imbalance between the rates of amplification and elimination of retrotransposons and satellite repeats. Genes and recombination events are evenly dispersed across chromosomes and the gene space is remarkably compact considering the genome size, although with substantial copy number variation driven by tandem duplication. Demonstrating practical application of the genome sequence, we develop a targeted genotyping assay and use high-resolution genome-wide association analysis to dissect the genetic basis of seed size and hilum colour. The resources presented constitute a genomics-based breeding platform for faba bean, enabling breeders and geneticists to accelerate the improvement of sustainable protein production across the Mediterranean, subtropical and northern temperate agroecological zones.

The mosaic oat genome gives insights into a uniquely healthy cereal crop.

Cultivated oat (Avena sativa L.) is an allohexaploid (AACCDD, 2n = 6x = 42) thought to have been domesticated more than 3,000 years ago while growing as a weed in wheat, emmer and barley fields in Anatolia1,2. Oat has a low carbon footprint, substantial health benefits and the potential to replace animal-based food products. However, the lack of a fully annotated reference genome has hampered efforts to deconvolute its complex evolutionary history and functional gene dynamics. Here we present a high-quality reference genome of A. sativa and close relatives of its diploid (Avena longiglumis, AA, 2n = 14) and tetraploid (Avena insularis, CCDD, 2n = 4x = 28) progenitors. We reveal the mosaic structure of the oat genome, trace large-scale genomic reorganizations in the polyploidization history of oat and illustrate a breeding barrier associated with the genome architecture of oat. We showcase detailed analyses of gene families implicated in human health and nutrition, which adds to the evidence supporting oat safety in gluten-free diets, and we perform mapping-by-sequencing of an agronomic trait related to water-use efficiency. This resource for the Avena genus will help to leverage knowledge from other cereal genomes, improve understanding of basic oat biology and accelerate genomics-assisted breeding and reanalysis of quantitative trait studies.

Semi-automated assembly of high-quality diploid human reference genomes.

The current human reference genome, GRCh38, represents over 20 years of effort to generate a high-quality assembly, which has benefitted society1,2. However, it still has many gaps and errors, and does not represent a biological genome as it is a blend of multiple individuals3,4. Recently, a high-quality telomere-to-telomere reference, CHM13, was generated with the latest long-read technologies, but it was derived from a hydatidiform mole cell line with a nearly homozygous genome5. To address these limitations, the Human Pangenome Reference Consortium formed with the goal of creating high-quality, cost-effective, diploid genome assemblies for a pangenome reference that represents human genetic diversity6. Here, in our first scientific report, we determined which combination of current genome sequencing and assembly approaches yield the most complete and accurate diploid genome assembly with minimal manual curation. Approaches that used highly accurate long reads and parent-child data with graph-based haplotype phasing during assembly outperformed those that did not. Developing a combination of the top-performing methods, we generated our first high-quality diploid reference assembly, containing only approximately four gaps per chromosome on average, with most chromosomes within ±1% of the length of CHM13. Nearly 48% of protein-coding genes have non-synonymous amino acid changes between haplotypes, and centromeric regions showed the highest diversity. Our findings serve as a foundation for assembling near-complete diploid human genomes at scale for a pangenome reference to capture global genetic variation from single nucleotides to structural rearrangements.

Constructing telomere-to-telomere diploid genome by polishing haploid nanopore-based assembly.

Draft genomes generated from Oxford Nanopore Technologies (ONT) long reads are known to have a higher error rate. Although existing genome polishers can enhance their quality, the error rate (including mismatches, indels and switching errors between paternal and maternal haplotypes) can be significant. Here, we develop two polishers, hypo-short and hypo-hybrid to address this issue. Hypo-short utilizes Illumina short reads to polish an ONT-based draft assembly, resulting in a high-quality assembly with low error rates and switching errors. Expanding on this, hypo-hybrid incorporates ONT long reads to further refine the assembly into a diploid representation. Leveraging on hypo-hybrid, we have created a diploid genome assembly pipeline called hypo-assembler. Hypo-assembler automates the generation of highly accurate, contiguous and nearly complete diploid assemblies using ONT long reads, Illumina short reads and optionally Hi-C reads. Notably, our solution even allows for the production of telomere-to-telomere diploid genomes with additional manual steps. As a proof of concept, we successfully assembled a fully phased telomere-to-telomere diploid genome of HG00733, achieving a quality value exceeding 50.

Scalable telomere-to-telomere assembly for diploid and polyploid genomes with double graph.

Despite advances in long-read sequencing technologies, constructing a near telomere-to-telomere assembly is still computationally demanding. Here we present hifiasm (UL), an efficient de novo assembly algorithm combining multiple sequencing technologies to scale up population-wide near telomere-to-telomere assemblies. Applied to 22 human and two plant genomes, our algorithm produces better diploid assemblies at a cost of an order of magnitude lower than existing methods, and it also works with polyploid genomes.

Evolutionary and biomedical insights from a marmoset diploid genome assembly.

The accurate and complete assembly of both haplotype sequences of a diploid organism is essential to understanding the role of variation in genome functions, phenotypes and diseases1. Here, using a trio-binning approach, we present a high-quality, diploid reference genome, with both haplotypes assembled independently at the chromosome level, for the common marmoset (Callithrix jacchus), an primate model system that is widely used in biomedical research2,3. The full spectrum of heterozygosity between the two haplotypes involves 1.36% of the genome-much higher than the 0.13% indicated by the standard estimation based on single-nucleotide heterozygosity alone. The de novo mutation rate is 0.43 × 10-8 per site per generation, and the paternal inherited genome acquired twice as many mutations as the maternal. Our diploid assembly enabled us to discover a recent expansion of the sex-differentiation region and unique evolutionary changes in the marmoset Y chromosome. In addition, we identified many genes with signatures of positive selection that might have contributed to the evolution of Callithrix biological features. Brain-related genes were highly conserved between marmosets and humans, although several genes experienced lineage-specific copy number variations or diversifying selection, with implications for the use of marmosets as a model system.

Aneuploidy renders cancer cells vulnerable to mitotic checkpoint inhibition.

Selective targeting of aneuploid cells is an attractive strategy for cancer treatment1. However, it is unclear whether aneuploidy generates any clinically relevant vulnerabilities in cancer cells. Here we mapped the aneuploidy landscapes of about 1,000 human cancer cell lines, and analysed genetic and chemical perturbation screens2-9 to identify cellular vulnerabilities associated with aneuploidy. We found that aneuploid cancer cells show increased sensitivity to genetic perturbation of core components of the spindle assembly checkpoint (SAC), which ensures the proper segregation of chromosomes during mitosis10. Unexpectedly, we also found that aneuploid cancer cells were less sensitive than diploid cells to short-term exposure to multiple SAC inhibitors. Indeed, aneuploid cancer cells became increasingly sensitive to inhibition of SAC over time. Aneuploid cells exhibited aberrant spindle geometry and dynamics, and kept dividing when the SAC was inhibited, resulting in the accumulation of mitotic defects, and in unstable and less-fit karyotypes. Therefore, although aneuploid cancer cells could overcome inhibition of SAC more readily than diploid cells, their long-term proliferation was jeopardized. We identified a specific mitotic kinesin, KIF18A, whose activity was perturbed in aneuploid cancer cells. Aneuploid cancer cells were particularly vulnerable to depletion of KIF18A, and KIF18A overexpression restored their response to SAC inhibition. Our results identify a therapeutically relevant, synthetic lethal interaction between aneuploidy and the SAC.

Interspecific transfer of genetic information through polyploid bridges.

Hybridization blurs species boundaries and leads to intertwined lineages resulting in reticulate evolution. Polyploidy, the outcome of whole genome duplication (WGD), has more recently been implicated in promoting and facilitating hybridization between polyploid species, potentially leading to adaptive introgression. However, because polyploid lineages are usually ephemeral states in the evolutionary history of life it is unclear whether WGD-potentiated hybridization has any appreciable effect on their diploid counterparts. Here, we develop a model of cytotype dynamics within mixed-ploidy populations to demonstrate that polyploidy can in fact serve as a bridge for gene flow between diploid lineages, where introgression is fully or partially hampered by the species barrier. Polyploid bridges emerge in the presence of triploid organisms, which despite critically low levels of fitness, can still allow the transfer of alleles between diploid states of independently evolving mixed-ploidy species. Notably, while marked genetic divergence prevents polyploid-mediated interspecific gene flow, we show that increased recombination rates can offset these evolutionary constraints, allowing a more efficient sorting of alleles at higher-ploidy levels before introgression into diploid gene pools. Additionally, we derive an analytical approximation for the rate of gene flow at the tetraploid level necessary to supersede introgression between diploids with nonzero introgression rates, which is especially relevant for plant species complexes, where interspecific gene flow is ubiquitous. Altogether, our results illustrate the potential impact of polyploid bridges on the (re)distribution of genetic material across ecological communities during evolution, representing a potential force behind reticulation.

Scattered differentiation of unlinked loci across the genome underlines ecological divergence of the selfing grass Brachypodium stacei.

Ecological divergence without geographic isolation, as an early speciation process that may lead finally to reproductive isolation through natural selection, remains a captivating topic in evolutionary biology. However, the pattern of genetic divergence underlying this process across the genome may vary between species and mating systems. Here, we present evidence that Brachypodium stacei, an annual and highly selfing grass model species, has undergone sympatric ecological divergence without geographic isolation. Genomic, transcriptomic, and metabolomic analyses together with lab experiments mimicking the two opposite environmental conditions suggest that diploid B. stacei populations have diverged sympatrically in two slopes characterized by distinct biomes at Evolution Canyon I (ECI), Mount Carmel, Israel. Despite ongoing gene flow, primarily facilitated by seed dispersal, the level of gene flow has progressively decreased over time. This local adaptation involves the scattered divergence of many unlinked loci across the total genome that include both coding genes and noncoding regions. Additionally, we have identified significant differential expressions of genes related to the ABA signaling pathway and contrasting metabolome composition between the arid- vs. forest-adapted B. stacei populations in ECI. These results suggest that multiple small loci involved in environmental responses act additively to account for ecological adaptations by this selfing species in contrasting environments.

Redundant microtubule crosslinkers prevent meiotic spindle bending to ensure diploid offspring in C. elegans.

Oocyte meiotic spindles mediate the expulsion of ¾ of the genome into polar bodies to generate diploid zygotes in nearly all animal species. Failures in this process result in aneuploid or polyploid offspring that are typically inviable. Accurate meiotic chromosome segregation and polar body extrusion require the spindle to elongate while maintaining its structural integrity. Previous studies have implicated three hypothetical activities during this process, including microtubule crosslinking, microtubule sliding and microtubule polymerization. However, how these activities regulate spindle rigidity and elongation as well as the exact proteins involved in the activities remain unclear. We discovered that C. elegans meiotic anaphase spindle integrity is maintained through redundant microtubule crosslinking activities of the Kinesin-5 family motor BMK-1, the microtubule bundling protein SPD-1/PRC1, and the Kinesin-4 family motor, KLP-19. Using time-lapse imaging, we found that single depletion of KLP-19KIF4A, SPD-1PRC1 or BMK-1Eg5 had minimal effects on anaphase B spindle elongation velocity. In contrast, double depletion of SPD-1PRC1 and BMK-1Eg5 or double depletion of KLP-19KIF4A and BMK-1Eg5 resulted in spindles that elongated faster, bent in a myosin-dependent manner, and had a high rate of polar body extrusion errors. Bending spindles frequently extruded both sets of segregating chromosomes into two separate polar bodies. Normal anaphase B velocity was observed after double depletion of KLP-19KIF4A and SPD-1PRC1. These results suggest that KLP-19KIF4A and SPD-1PRC1 act in different pathways, each redundant with a separate BMK-1Eg5 pathway in regulating meiotic spindle elongation. Depletion of ZYG-8, a doublecortin-related microtubule binding protein, led to slower anaphase B spindle elongation. We found that ZYG-8DCLK1 acts by excluding SPD-1PRC1 from the spindle. Thus, three mechanistically distinct microtubule regulation modules, two based on crosslinking, and one based on exclusion of crosslinkers, power the mechanism that drives spindle elongation and structural integrity during anaphase B of C.elegans female meiosis.

Genome architecture and genetic diversity of allopolyploid okra (Abelmoschus esculentus).

The allopolyploid okra (Abelmoschus esculentus) unveiled telomeric repeats flanking distal gene-rich regions and short interstitial TTTAGGG telomeric repeats, possibly representing hallmarks of chromosomal speciation. Ribosomal RNA (rRNA) genes organize into 5S clusters, distinct from the 18S-5.8S-28S units, indicating an S-type rRNA gene arrangement. The assembly, in line with cytogenetic and cytometry observations, identifies 65 chromosomes and a 1.45 Gb genome size estimate in a haploid sibling. The lack of aberrant meiotic configurations implies limited to no recombination among sub-genomes. k-mer distribution analysis reveals 75% has a diploid nature and 15% heterozygosity. The configurations of Benchmarking Universal Single-Copy Ortholog (BUSCO), k-mer, and repeat clustering point to the presence of at least two sub-genomes one with 30 and the other with 35 chromosomes, indicating the allopolyploid nature of the okra genome. Over 130 000 putative genes, derived from mapped IsoSeq data and transcriptome data from public okra accessions, exhibit a low genetic diversity of one single nucleotide polymorphisms per 2.1 kbp. The genes are predominantly located at the distal chromosome ends, declining toward central scaffold domains. Long terminal repeat retrotransposons prevail in central domains, consistent with the observed pericentromeric heterochromatin and distal euchromatin. Disparities in paralogous gene counts suggest potential sub-genome differentiation implying possible sub-genome dominance. Amino acid query sequences of putative genes facilitated phenol biosynthesis pathway annotation. Comparison with manually curated reference KEGG pathways from related Malvaceae species reveals the genetic basis for putative enzyme coding genes that likely enable metabolic reactions involved in the biosynthesis of dietary and therapeutic compounds in okra.

Investigation of regulatory divergence between homoeologs in the recently formed allopolyploids, Tragopogon mirus and T. miscellus (Asteraceae).

Polyploidy is an important evolutionary process throughout eukaryotes, particularly in flowering plants. Duplicated gene pairs (homoeologs) in allopolyploids provide additional genetic resources for changes in molecular, biochemical, and physiological mechanisms that result in evolutionary novelty. Therefore, understanding how divergent genomes and their regulatory networks reconcile is vital for unraveling the role of polyploidy in plant evolution. Here, we compared the leaf transcriptomes of recently formed natural allotetraploids (Tragopogon mirus and T. miscellus) and their diploid parents (T. porrifolius X T. dubius and T. pratensis X T. dubius, respectively). Analysis of 35 400 expressed loci showed a significantly higher level of transcriptomic additivity compared to old polyploids; only 22% were non-additively expressed in the polyploids, with 5.9% exhibiting transgressive expression (lower or higher expression in the polyploids than in the diploid parents). Among approximately 7400 common orthologous regions (COREs), most loci in both allopolyploids exhibited expression patterns that were vertically inherited from their diploid parents. However, 18% and 20.3% of the loci showed novel expression bias patterns in T. mirus and T. miscellus, respectively. The expression changes of 1500 COREs were explained by cis-regulatory divergence (the condition in which the two parental subgenomes do not interact) between the diploid parents, whereas only about 423 and 461 of the gene expression changes represent trans-effects (the two parental subgenomes interact) in T. mirus and T. miscellus, respectively. The low degree of both non-additivity and trans-effects on gene expression may present the ongoing evolutionary processes of the newly formed Tragopogon polyploids (~80-90 years).

Aneuploidy Can Be an Evolutionary Diversion on the Path to Adaptation.

Aneuploidy is common in eukaryotes, often leading to decreased fitness. However, evidence from fungi and human tumur cells suggests that specific aneuploidies can be beneficial under stressful conditions and facilitate adaptation. In a previous evolutionary experiment with yeast, populations evolving under heat stress became aneuploid, only to later revert to euploidy after beneficial mutations accumulated. It was therefore suggested that aneuploidy is a "stepping stone" on the path to adaptation. Here, we test this hypothesis. We use Bayesian inference to fit an evolutionary model with both aneuploidy and mutation to the experimental results. We then predict the genotype frequency dynamics during the experiment, demonstrating that most of the evolved euploid population likely did not descend from aneuploid cells, but rather from the euploid wild-type population. Our model shows how the beneficial mutation supply-the product of population size and beneficial mutation rate-determines the evolutionary dynamics: with low supply, much of the evolved population descends from aneuploid cells; but with high supply, beneficial mutations are generated fast enough to outcompete aneuploidy due to its inherent fitness cost. Our results suggest that despite its potential fitness benefits under stress, aneuploidy can be an evolutionary "diversion" rather than a "stepping stone": it can delay, rather than facilitate, the adaptation of the population, and cells that become aneuploid may leave less descendants compared to cells that remain diploid.