Cited by 392Genomics and Phylogenetic StudiesChromosomal and Genetic VariationsGenome Rearrangement Algorithms

The DOI matches its Crossref record; 53 of the 54 references it cites were found in a registry or an index; 16 of the 34 sources its claims cite were checked for support.

GENESPACE tracks regions of interest and gene copy number variation across multiple genomes

Open access

Published version of a preprint

Abstract

The development of multiple chromosome-scale reference genome sequences in many taxonomic groups has yielded a high-resolution view of the patterns and processes of molecular evolution. Nonetheless, leveraging information across multiple genomes remains a significant challenge in nearly all eukaryotic systems. These challenges range from studying the evolution of chromosome structure, to finding candidate genes for quantitative trait loci, to testing hypotheses about speciation and adaptation. Here, we present GENESPACE, which addresses these challenges by integrating conserved gene order and orthology to define the expected physical position of all genes across multiple genomes. We demonstrate this utility by dissecting presence–absence, copy-number, and structural variation at three levels of biological organization: spanning 300 million years of vertebrate sex chromosome evolution, across the diversity of the Poaceae (grass) plant family, and among 26 maize cultivars. The methods to build and visualize syntenic orthology in the GENESPACE R package offer a significant addition to existing gene family and synteny programs, especially in polyploid, outbred, and other complex genomes.

Overview

AI snapshot

Written from the paper’s full text.

TL;DR

GENESPACE is an R package/pipeline that integrates conserved gene order (synteny) and orthology to build pan-genome annotations and track presence–absence, copy-number, and structural variation across multiple genomes. It outperformed OrthoFinder and MCScanX in single-copy orthogroup and syntenic block recovery across diverse plant and vertebrate comparisons, including polyploid and recently diverged genomes.

Key findings

  • Across seven genome pairs, GENESPACE increased the percentage of genes in single-copy orthogroups and single-copy syntenic blocks versus OrthoFinder/MCScanX; for cotton, single-copy orthogroups rose from 35.6% (OrthoFinder) to 85.7% (GENESPACE) and single-copy syntenic blocks from 14.1% (orthogroup-constrained MCScanX) to 96.2% (GENESPACE).
  • In vertebrate sex chromosome comparisons, a 91.7 M bp region is shared between human X and sand lizard Z, representing 59.0% of human X genic sequence; 107.2 M bp (68.8%) of human X is syntenic with 77.8 M bp (93.9%) of Tasmanian devil X; chicken Z is retained across five avian genomes except for a budgie partial fusion with a 19.5 M bp segment of chicken chromosome 11.
  • For 86 C3/C4 photosynthesis genes, absence was lower than background (9.7% vs. 38.2%, odds ratio = 5.7, p < 1 × 10−16), but C3 species had more absences than C4 species (15.3% vs. 5.5%, odds ratio = 3.1, p = 6.25 × 10−8); among 190 maize NAM QTLs, three had single-parent outlier effects, and the two Mo18W QTL regions fell within an 11.7 M bp Mo18W-specific inversion.
  • Rho WGD paralogs were less conserved than primary syntenic orthologs (median peptide identity 90.6% vs. 71.4%, 26.9% more similar; Wilcoxon p < 1 × 10−16), except one over-retained S. viridis chromosome 8 region where paralogs were not significantly less conserved (91.6% vs. 91.9%, p = 0.13).

Methods

  • GENESPACE takes gff3 annotations and peptide FASTAs, computes DIAMOND2 BLAST-like hits, runs OrthoFinder for orthogroups, extracts syntenic regions with MCScanX-based and graph/cluster approaches, condenses tandem arrays to representative genes, and collapses syntenic orthogroups into a reference-anchored pan-genome annotation.
  • Analyses used seven haploid genome pairs spanning <0.01 Mya maize to 310 Mya human–chicken, 15 avian/mammalian genomes plus two reptile outgroups, eight grass genomes, and 26 maize NAM founders; public annotations were downloaded on or before October 8, 2021.
  • Comparisons calculated percent genes in single-copy orthogroups and single-copy syntenic blocks; statistics included Wilcoxon signed rank tests, ungapped percent peptide identity from Needleman–Wunsch alignments, and Grubbs tests for outliers, using R 4.1.2 and GENESPACE v0.9.3.

Limitations

  • Stated: only a single C3/C4 phylogenetic split prevents testing evolutionary hypotheses about causes of PAV; the phylogenetic scale precludes inference about the origin of shared sex chromosome sequences; and it is not clear that within-block OrthoFinder improves much over global orthogroups for purely haploid genomes.
  • Stated/evident: GENESPACE relies on syntenic regions, so errors in syntenic block coordinates can have major downstream effects; true orthologs can exist outside syntenic regions, with 1096 non-syntenic human–chimpanzee orthologs (6.3% of all orthologs) and 9002 non-syntenic Sorghum–Brachypodium orthologs.
  • Evident: results depend on parameterization, outgroup choice, and ploidy treatment; polyploids had only a primary haplotype assembled into chromosomes, maize was treated as diploid in the multispecies run but haploid in the NAM run, and some maize alternative haplotype scaffolds were dropped.

Model-written summary. Check it against the paper before relying on it.

Integrity signals

Registry

Crossref · 10.7554/elife.78526

Cites

53 of 53 references verified0 citation errors · 0 not checkable

Statistics

Nothing to recompute.

Figures

Figure 1. GENESPACE synteny and pan-genome annotation methods. (A, grey panel) GENESPACE runs and parses OrthoFinder results into a synteny-constrained pan-genome annotation. (B, purple panel) Chromosome, gene rank order, and orthogroup membership are added to BLAST hits, which allows direct integration between estimates of orthology and synteny. The three dotplots present the efficacy of GENESPACE syntenic blocks by exploring a particularly challenging region on human (x-axis) and chimpanzee (y-axis) chr. 6. Each point is a BLAST hit rank-order position, colored by syntenic block; colors are recycled if there are more than eight blocks. (C, green panel) Synteny-constrained orthogroups and optionally non-syntenic orthologs are decomposed into a pan-genome annotation where each orthogroup is placed at its inferred syntenic position.
Figure 1. GENESPACE synteny and pan-genome annotation methods. (A, grey panel) GENESPACE runs and parses OrthoFinder results into a synteny-constrained pan-genome annotation. (B, purple panel) Chromosome, gene rank order, and orthogroup membership are added to BLAST hits, which allows direct integration between estimates of orthology and synteny. The three dotplots present the efficacy of GENESPACE syntenic blocks by exploring a particularly challenging region on human (x-axis) and chimpanzee (y-axis) chr. 6. Each point is a BLAST hit rank-order position, colored by syntenic block; colors are recycled if there are more than eight blocks. (C, green panel) Synteny-constrained orthogroups and optionally non-syntenic orthologs are decomposed into a pan-genome annotation where each orthogroup is placed at its inferred syntenic position.CC BY · PubMed Central
Figure 2. Sex chromosome syntenic network across 17 representative vertebrate genomes. The plot was generated by the plot_riparian GENESPACE function. Genomes are ordered vertically to minimize the number of translocations between each pairwise combination. Chromosomes are ordered horizontally to maximize synteny with the human chromosomes [X, 1–22]. Regions containing syntenic orthogroup members to the mammalian X (gold) or avian Z (blue) chromosomes are highlighted. All sex chromosomes are represented by red segments while autosomes are white. Chromosome segment sizes are scaled by the total number of genes in syntenic networks and positions of the braids are the gene order along the chromosome sequence. See Figure 2—figure supplement 1 for the full synteny graph including autosomes and chromosome labels.
Figure 2. Sex chromosome syntenic network across 17 representative vertebrate genomes. The plot was generated by the plot_riparian GENESPACE function. Genomes are ordered vertically to minimize the number of translocations between each pairwise combination. Chromosomes are ordered horizontally to maximize synteny with the human chromosomes [X, 1–22]. Regions containing syntenic orthogroup members to the mammalian X (gold) or avian Z (blue) chromosomes are highlighted. All sex chromosomes are represented by red segments while autosomes are white. Chromosome segment sizes are scaled by the total number of genes in syntenic networks and positions of the braids are the gene order along the chromosome sequence. See Figure 2—figure supplement 1 for the full synteny graph including autosomes and chromosome labels.CC BY · PubMed Central
Figure 2—figure supplement 1. Complete map of synteny, color-coded by synteny with human chromosomes X, 1-22.
Figure 2—figure supplement 1. Complete map of synteny, color-coded by synteny with human chromosomes X, 1-22.CC BY · PubMed Central
Figure 3. Comparative–quantitative genomics in the grasses. (A) The GENESPACE syntenic map (‘riparian plot’) of orthologous regions among eight grass genomes. Chromosomes are ordered horizontally to maximize synteny with rice and ribbons are color coded by synteny to rice chromosomes. Genomes are ordered vertically by general phylogenetic positions. (B) The upper bars display the proportion of maize gene models without syntenic orthologs (‘absent’) in each genome, split by the full background (dark colors) and 86 C3/C4 genes (light colors). (C) The proportion of absent genes is higher in the C3 genomes (green bars), even when controlling for more global gene absences (lower odds ratios). (D) Syntenic orthologs, excluding homeologs among the 26 maize nested association mapping (NAM) founder genomes, with two quantitative trait loci (QTL) intervals highlighted on chromosome 3 (‘Chr3’) and chromosome 6 (’Chr6’). (E) Focal QTL regions that affect productivity in drought where only the genome that drives the QTL effect (middle), the top (B73) and bottom (Tzi8) genomes are presented and the region plotted is restricted to the physical B73 QTL interval and a 25 M bp buffer on either side. Note that the Chr3 QTL disarticulates into two intervals. Due to a larger number of potential candidate genes, the larger Chr3 region, flagged with **, is explored separately in Figure 3—figure supplement 1. (F) Presence–absence and copy number variation are presented for two of the three intervals as heatmaps where each row is a genome (order following panel D), each column is a pan-genome entry (see Figure 1), and the color of each tile indicates absence (gray), single copy (light blue), and multicopy (dark blue). PAV/CNV of the focal genome is outlined. For each interval, the estimated QTL allelic effect relative to B73 of each genome is plotted as bars to the right of the heatmap.
Figure 3. Comparative–quantitative genomics in the grasses. (A) The GENESPACE syntenic map (‘riparian plot’) of orthologous regions among eight grass genomes. Chromosomes are ordered horizontally to maximize synteny with rice and ribbons are color coded by synteny to rice chromosomes. Genomes are ordered vertically by general phylogenetic positions. (B) The upper bars display the proportion of maize gene models without syntenic orthologs (‘absent’) in each genome, split by the full background (dark colors) and 86 C3/C4 genes (light colors). (C) The proportion of absent genes is higher in the C3 genomes (green bars), even when controlling for more global gene absences (lower odds ratios). (D) Syntenic orthologs, excluding homeologs among the 26 maize nested association mapping (NAM) founder genomes, with two quantitative trait loci (QTL) intervals highlighted on chromosome 3 (‘Chr3’) and chromosome 6 (’Chr6’). (E) Focal QTL regions that affect productivity in drought where only the genome that drives the QTL effect (middle), the top (B73) and bottom (Tzi8) genomes are presented and the region plotted is restricted to the physical B73 QTL interval and a 25 M bp buffer on either side. Note that the Chr3 QTL disarticulates into two intervals. Due to a larger number of potential candidate genes, the larger Chr3 region, flagged with **, is explored separately in Figure 3—figure supplement 1. (F) Presence–absence and copy number variation are presented for two of the three intervals as heatmaps where each row is a genome (order following panel D), each column is a pan-genome entry (see Figure 1), and the color of each tile indicates absence (gray), single copy (light blue), and multicopy (dark blue). PAV/CNV of the focal genome is outlined. For each interval, the estimated QTL allelic effect relative to B73 of each genome is plotted as bars to the right of the heatmap.CC BY · PubMed Central
Figure 3—figure supplement 1. Map of PAV in the larger MO18W chromosome 3 QTL.
Figure 3—figure supplement 1. Map of PAV in the larger MO18W chromosome 3 QTL.CC BY · PubMed Central
Figure 4. Analysis of the grass Rho WGD. (A) BLAST hits between P. hallii and S. viridis where the target and query genes were in the same orthogroup are plotted and color coded by sequence similarity. Two over-retained regions are highlighted in the red and yellow boxes. (B) The protein identity of S. viridis chromosome 8 primary orthologous (blue line) hits against P. hallii chromosome 8 and the secondary hits (orange line) against P. hallii chromosome 3 demonstrate sequence conservation heterogeneity. The region between the two red vertical lines corresponds to the red-boxed over-retained primary block in panel A. (C) The two boxed regions in panel A were tracked from P. hallii chromosomes 3 (red) and 8 (yellow); 50% transparency of the braids means that overlapping regions appear orange.
Figure 4. Analysis of the grass Rho WGD. (A) BLAST hits between P. hallii and S. viridis where the target and query genes were in the same orthogroup are plotted and color coded by sequence similarity. Two over-retained regions are highlighted in the red and yellow boxes. (B) The protein identity of S. viridis chromosome 8 primary orthologous (blue line) hits against P. hallii chromosome 8 and the secondary hits (orange line) against P. hallii chromosome 3 demonstrate sequence conservation heterogeneity. The region between the two red vertical lines corresponds to the red-boxed over-retained primary block in panel A. (C) The two boxed regions in panel A were tracked from P. hallii chromosomes 3 (red) and 8 (yellow); 50% transparency of the braids means that overlapping regions appear orange.CC BY · PubMed Central

Verification

Claims the paper makes

32 Own interpretation · 30 Own result · 24 No basis given · 9 Partly supported · 6 Sources not found · 5 Not checked yet · 4 Own method · 3 No full text of the sources · 2 Supported · accuracy not yet measured

Verify all 4 references — free account

Show all 58 claims

  • Probably not a substantial claim. 35
  • It may state two separate claims as one. 30
57 of 57
  • The authors state this as an assumption.

    Homologs should be exactly single copy within any syntenic region between a pair of genomes.

    “GENESPACE resolves this circular problem by operating on a foundational assumption: homologs should be exactly single copy within any syntenic region between a pair of genomes.”
  • The paper relies on this without support.

    The three vertebrate genome comparisons have no recent history of WGDs, while WGDs predated evolutionary divergence of the four plant contrasts.

    “The three comparisons between vertebrate genomes have no recent history of WGDs, while WGDs predated evolutionary divergence of the four plant contrasts: maize is a 12 M ya paleo-tetraploid, cotton is a 1.6 M ya meso-tetraploid, and the ~70 M ya Rho WGD predated grass diversification.”
  • The paper relies on this without support.

    The genomes in each comparison share a common ancestor that arose after the WGD.

    “As such, we treated all genomes as haploid and did not include an outgroup; therefore, orthogroups should not include homeologs since the genomes share a common ancestor that arose after the WGD (see ‘additional considerations’ methods section: Outgroups and the phylogenetic context of orthology inference).”
  • The paper relies on this without support.

    The heteromorphic chromosomes (Y and W) are often unassembled, or, where assemblies exist, lack sufficient synteny to provide a useful metric for comparative genomics.

    “The heteromorphic chromosomes (Y and W) are often unassembled, or, where assemblies exist, lack sufficient synteny to provide a useful metric for comparative genomics.”
  • The paper relies on this without support.

    Many genetic mapping populations now have reference genome sequences for all parents.

    “However, many genetic mapping populations now have reference genome sequences for all parents; this offers an opportunity to explore variation among functional alleles and PAV, which would be impossible with a single reference genome.”
  • The paper relies on this without support.

    Inversions reduce recombination.

    “Since inversions reduce recombination, it is possible that multiple Mo18W causal variants have been fixed in linkage disequilibrium in this NAM population.”
  • The paper relies on this without support.

    In most pairwise comparisons between genomes, some syntenic orthogroups will be missing in the reference.

    “In most pairwise comparisons between genomes, some syntenic orthogroups will be missing in the reference.”
  • The paper relies on this without support.

    Most plant lineages have undergone one or more whole-genome duplications.

    “most plant lineages have undergone one or more WGDs”
  • The paper relies on this without support.

    The coalescence of any two maize genotypes occurred more recently than the ~12M ya whole-genome duplication.

    “the coalescence of any two maize genotypes occurred more recently than the ~12M ya WGD”
9 of 9

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

Each contradiction between two passages listed here was confirmed by an independent check.

1 of 1

Citations

Cites

Show all 54 references

Cited by

Sorted by citations

Items 301–392 of 392

Metadata

Authors
Lovell, J. T., Sreedasyam, A., Schranz, M. E., Wilson, M., Carlson, J. W., Harkess, A., Emms, D., Goodstein, D. M., & Schmutz, J.
Affiliations
  • Lawrence Berkeley National Laboratory (US)
  • HudsonAlpha Institute for Biotechnology (US)
  • Wageningen University & Research (NL)
  • Arizona State University (US)
  • University of Oxford (GB)
Corresponding author
Lovell, John T
Venue
eLife
Published
2022-09-09
Volume
11
ISSN
2050-084X (ISSN-L)
OpenAlex
W4295078843
Cited by
391
References
54
Citations by year
2026: 144 · 2025: 130 · 2024: 78 · 2023: 37 · 2022: 2
FWCI
18.24
Type
journal article
Language
en
Open access
Locations
  • eLife (publishedVersion, cc-by)
  • bioRxiv (Cold Spring Harbor Laboratory) (submittedVersion, cc-by-nc-nd)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • Socio-Environmental Systems Modeling (acceptedVersion, cc-by)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion, cc-by-sa)
  • eScholarship (California Digital Library) (submittedVersion)
  • OSTI OAI (U.S. Department of Energy Office of Scientific and Technical Information) (submittedVersion)
  • Data Archiving and Networked Services (DANS) (submittedVersion, cc-by)
Funders
  • U.S. Department of Energy (DE-AC02-05CH1123)
  • National Institute of General Medical Sciences (R35GM124827)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Genomics and Phylogenetic Studies
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Genetics › Chromosomal and Genetic Variations
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Genetics › Genome Rearrangement Algorithms
Keywords
gene order, orthology, synteny, copy number variation, presence/absence variation, structural variants, sex chromosome evolution, Poaceae, maize varieties, polyploidy, comparative genomics, R package
MeSH terms
Gene Dosage, Genome, Plant, Evolution, Molecular, Synteny, Quantitative Trait Loci, DNA Copy Number Variations

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

A language model drew the claims from the paper on 2026-11-10, under version 3 of the claim definition (model deepseek-v4.1-flash, prompt extract-v3).

A second model labelled each claim on 2026-11-10 (model jev-1.13-20260917, definition revision gold-guideline-v0.2).

Cited sources are checked with a language model’s help (gpt-6-sol).

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review