In the course of evolution, genomes continue to retain most genes with occasional duplications, while losing some genes (Blomme et al., 2006; Fernández and Gabaldón, 2020; Shen et al., 2018). This retention and loss can be interpreted as gene fate; genes are stably retained in the genome, but some factors may cause them to transition to a state where deletion occurs. Accordingly, identification of the factors allowing gene loss may facilitate our understanding of gene fate. Gene retention or loss has generally been considered to depend largely on the functional importance of the particular gene from the perspective of molecular evolutionary biology (Albalat and Cañestro, 2016; Bartha et al., 2018; Blanc et al., 2012; Liu et al., 2015; Olson, 1999; Sharma et al., 2018; Shen et al., 2018). Genes with indispensable functions have usually been retained with highly conserved sequences in genomes through rapid elimination of alleles that impair gene functions (Hirsh and Fraser, 2001; Krylov et al., 2003; Miyata et al., 1980; Pál et al., 2006). On the contrary, genes with less important functions are likely to accept more mutations and structural variations, which can degrade the original functions, leading to gene loss through pseudogenization or genomic deletion (Jordan et al., 2002; Yang et al., 2003). To date, gene loss has been imputed to the relaxation of functional constraints of individual genes. Gene loss has further been revealed to drive phenotypic adaptation in various organisms (Albalat and Cañestro, 2016; Olson, 1999), as well as in a gene knockout collection of yeasts in culture (Giaever and Nislow, 2014; Maclean et al., 2017).

To uncover the association between fates and functional importance of the genes, molecular evolutionary analyses have been conducted at various scales, from gene-by-gene to genome-wide. A number of studies have revealed that the genes with reduced non-synonymous substitution rates (or K_A values) and ratios of non-synonymous to synonymous substitution rates (K_A/K_S ratios) are less likely to be lost (Jordan et al., 2002; Yang et al., 2003). A genome-wide comparison of duplicated genes in yeast revealed larger K_A values for those lost in multiple lineages than those retained by all the species investigated (Byrne and Wolfe, 2007). Other comprehensive studies of gene loss across metazoans and teleosts revealed that the genes expressed in the central nervous system are less prone to loss (Fernández and Gabaldón, 2020; Roux et al., 2017). These observations again suggest that gene fate depends on the functional constraints of a particular gene.

Besides functional constraints, several studies have identified the genes lost independently in multiple lineages, revealing that the genomic regions containing these genes ‘prefer’ particular characteristics associated with structural instability (Cortez et al., 2014; Hughes et al., 2012; Lewin et al., 2021; Maeso et al., 2016). In mammals, tandemly arrayed homeobox genes derived from the Crx gene family were lost in multiple species (Lewin et al., 2021; Maeso et al., 2016). The findings suggest that genomic features containing tandem duplications facilitate unequal crossing over, leading to frequent gene loss. Mammalian chromosome Y, which contains abundant repetitive elements and continues to reduce in size, has lost a considerable number of genes (Cortez et al., 2014; Hughes et al., 2012). In the stickleback genome, a Pitx1 enhancer was independently lost in multiple lineages inhabiting freshwater due to its genomic location in a structurally fragile site, leading to recurrent loss of pelvic fins (Xie et al., 2019). Genes and genomic elements in such particular regions may be prone to loss in a more neutral manner than the relaxation of functional importance or via functional adaptations. Accordingly, these studies focusing on the particular genomic regions led us to search for the common features in genomes that potentially facilitate gene loss. Genome-wide scans have revealed heterogeneous distributions of a variety of sequence and structural features so far, for example, base composition (Bernardi and Bernardi, 1986; Cohen et al., 2005; Katzman et al., 2011), the frequency of repetitive elements (Korenberg and Rykowski, 1988; Medstrand et al., 2002), and DNA-damage sensitivity induced by replication inhibitors (Debatisse et al., 2012; Helmrich et al., 2006). However, the extent to which these characteristics are associated with gene fates has not been understood well at a genome-wide level.

The accumulation of near-complete genome assemblies for various organisms facilitates comprehensive taxon-wide analysis of gene loss (Fernández and Gabaldón, 2020; Guijarro-Clarke et al., 2020; Rice and McLysaght, 2017). Along with this motivation, we recently performed a comprehensive analysis on the fate of paralogs generated via the two-round whole-genome duplications in early vertebrates (Hara et al., 2018a). The results revealed that the genes retained by reptiles but lost in mammals and Aves rapidly accumulated not only non-synonymous but also synonymous substitutions in comparison with the counterparts retained by almost all the vertebrates examined, indicating that those genes prone to loss show increasing mutation rates. Furthermore, these loss-prone genes were located in genomic regions with high GC contents, high gene densities, and high repetitive element frequencies. These findings suggest that the fates of those genes are influenced not only by functional constraints but also by intrinsic genomic characteristics. Because the findings were restricted to a set of particular genes, they prompted us to examine whether this trend is associated with gene fates on a genome-wide scale.

In this study, we inferred molecular phylogenies of vertebrate orthologs to systematically search for the genes harboring different fates in the human genome. We previously referred to the nature of genes prone to loss as ‘elusive’ (Hara et al., 2018a; Hara et al., 2018b). In this study, we define the elusive genes as those that are retained by modern humans but have been lost independently in multiple mammalian lineages. As a comparison of the elusive genes, we retrieved the genes that were retained by almost all of the mammalian species examined and defined them as ‘non-elusive,’ representing those persistent in the genomes. We conducted a careful search for gene loss to reduce the false discovery rate (FDR), which is usually caused by incomplete sequence information (Botero-Castro et al., 2017; Deutekom et al., 2019). By comparing the genomic regions containing these genes, we uncovered genomic characteristics relevant to gene loss. We associated the elusive genes with a variety of findings from deep sequencing analyses of the human genome, including transcriptomics, epigenomics, and genetic variations. These data assisted us to understand how intrinsic genomic features may affect gene fate, leading to gene loss by decreasing the expression level and eventually relaxing the functional importance of ‘elusive’ genes.

Here we identified elusive genes that were lost in multiple lineages during mammalian evolution using a comprehensive scan of gene phylogenies. To identify gene loss events, absence of evidence (i.e. missing genes caused by incomplete genome assemblies and gene annotations) should be reviewed meticulously (Deutekom et al., 2019). Additionally, gene loss might be detected erroneously because of failure in similarity searches for orthologs of rapidly evolving genes (Moyers and Zhang, 2015). In this study, we aimed to reduce these false discoveries through our multifaceted approaches (Figure 1). We selected those species with highly complete gene annotations through integration of multiple gene annotations. Using these improved gene annotations, we created orthologous groups by employing a highly sensitive homology search with MMSeqs2 (Steinegger and Söding, 2017) and merged them into those identified in the Ensembl database. Furthermore, we restricted the loss events that were observed as gene absence in all species examined within all hierarchical levels of the selected taxonomic groups (Figure 1b). This absence is likely to have occurred as a gene loss in the common ancestor of the particular taxon rather than as a false discovery of gene loss in the individual species independently. Genuine continuous (e.g. telomere-to-telomere) genome assemblies are now available using modern sequencing technologies (Nurk et al., 2022). These genomic assemblies may help relieve the labor of examining for information losses, thereby facilitating the identification of genuine gene loss in any given species.

In the human genome, the elusive genes and their flanking regions harbor particular characteristics, including high GC content and high gene density, that may have originated long before the emergence of mammals (Figure 3). Frequent synonymous variations across modern humans in the elusive genes, consistent with higher synonymous substitution rates between the vertebrate orthologs, suggest that the genomic regions including elusive genes have been subject to rapid evolution for approximately 500 million years (Figures 2 and 4). Our findings indicate that heterogeneous genomic characteristics potentially affect the fate of genes at the latest period of vertebrate evolution. Analyses with large numbers of germline mutations in the human genome have illustrated the heterogeneity of mutation rates (Campbell and Eichler, 2013; Seplyarskiy and Sunyaev, 2021; Terekhanova et al., 2017). High GC content in the elusive genes may have facilitated an elevation of the mutation rate, as observed in the enrichment of rare variants in high-GC regions in the human genome (Schaibley et al., 2013). In addition, some of the elusive genes appear to have retained particular epigenetic marks including sparse ATAC-seq peaks, late replication timings, and location within LADs (Figure 6—figure supplements 1–4); these epigenetic marks are relevant to an increase in the mutation rate. Genomic regions with late replication timing exhibit increased mutation rates because of their unstable structure during the S-phase of the cell cycle (Koren et al., 2012; Stamatoyannopoulos et al., 2009). LADs retain more G-to-T mutations because of their susceptibility to oxidative damage in the nuclear periphery resulting in high levels of 8-oxoguanine (Yoshihara et al., 2014). Close coordination of the studies on gene evolution with germline mutation repertoires and spectra, which can be approximated from the collection of de novo mutations obtained by trio sequencing, may further facilitate our understanding of gene fates driven by heterogeneous genomic features—this would be viewed as ‘mutation-driven’ evolution (Nei, 2013).

The epigenetic marks of elusive genes are relevant to the suppression of gene expression (van Steensel and Belmont, 2017), and indeed, these genes harbor weakened and spatially restricted expression profiles (Figures 5 and 6 and Figure 6—figure supplements 1–4). However, the genomic features associated with these epigenetic marks usually exhibit lower GC contents and reduced gene density (Gilbert et al., 2004; Rao et al., 2014; van Steensel and Belmont, 2017). This discrepancy may be caused in part by a gain of local heterochromatin accompanied by suppression of the expression of transposable elements, as observed in various eukaryotic genomes (Choi and Lee, 2020; Fiston-Lavier et al., 2007; Grewal and Jia, 2007; Rangasamy, 2013; Slotkin and Martienssen, 2007; Underwood et al., 2017). Previous analyses showed frequent heterochromatinization of the human genomic regions where KRAB zinc finger genes colocalize with L1 retrotransposons (Imbeault et al., 2017; O’Geen et al., 2007; Vogel et al., 2006). One of the genomic regions found in human chromosome region 19p12 also contains many elusive genes (Vogel et al., 2006; Figure 7). Closer attention to the local gene and repeat contents including repetitive elements and tandem gene clusters might facilitate our understanding of heterochromatinization in restricted genomic regions, although we excluded such gene clusters in our search for elusive genes (Figure 1a).

A chromosomal-scale view of the distribution of elusive genes illuminated their significant correlation with the genes whose chicken orthologs are located on microchromosomes (Figure 7a). More importantly, genomic regions rich in elusive genes were traced back to the microchromosomes of the ancestral gnathostomes by reconstructing chromosomes of the ancestral genomes (Figure 7b). This inference of ancestral karyotypes augments our observations that some elusive natures of genomic sequences have been retained for hundreds of millions of years (Figure 3). In other words, the result suggests that the disparity of genomic regions that allows the ‘elusiveness’ for the genes has been retained during vertebrate evolution. On the other hand, comparisons of the expression profiles between the orthologs of the elusive and non-elusive genes for non-mammalian vertebrates suggest that the orthologs of the elusive genes have been associated with a reduction in pleiotropy of gene expression since vertebrate ancestors but acquired the diverse expressions in chicken and frog (Figure 5—figure supplement 1). Additionally, in the chicken genome, the diverse expressions of the chicken orthologs of the human elusive genes may be correlated with the abundance of ATAC-seq peaks (Figure 6—figure supplement 5). These findings again suggest that the chicken orthologs of the human elusive genes have increased pleiotropy of gene expression, which may lead to a lineage-specific acquisition of functional indispensability. It should be noted that the choices of tissues used in these analyses were largely different between the human and non-mammalian vertebrates (Tables S3 and S4; Supplementary file 1b and c). The chicken ATAC-seq data could be obtained only from developing embryos, while the human ATAC-seq in ENCODE were performed with cell lines. Therefore, the aforementioned interpretation should be treated carefully.

Finally, we note the potential evolutionary courses that facilitate the transition of gene fate from retention to loss. One possible course is a decrease in essential functions because of rapid sequence evolution in local genomic regions. The elusive genes located in those genomic regions with rapidly evolving characteristics are likely to accumulate neutral or even moderately harmful mutations in coding regions frequently, resulting in impaired essential functions. Another factor is the spatiotemporal suppression of gene expression via epigenetic constraints. Previous studies showed that lowly expressed genes are associated with low functional essentiality (Cherry, 2010; Gout et al., 2010), as shown for elusive genes in our study. Elusive genes with reduced pleiotropy may have limited opportunities to function, potentially leading to loss of their important roles. The extent of these evolutionary forces may have varied with time and lineages, resulting in a patchy loss of elusive genes phylogenetically. Interestingly, a recent large-scale scan of de novo mutations in Arabidopsis indicates the association of mutation rates with epigenetic features and functional essentiality of genes (Monroe et al., 2022). Further investigation of the association of genes with the surrounding genomic regions in various taxa may provide a common understanding of genomic and epigenomic features that potentially alter the fate of genes. Although epigenetic features are plastic, our findings indicate that the disparities of genomic regions are reflected in the heterogeneity of evolutionary forces and have been retained for hundreds of millions of years. This idea prompts us to explore evolutionary constraints on more global genomic regions that are potentially associated with structural characteristics including chromosomal composition and locations within the nucleus.

The current manuscript is a computational study, so no data have been generated for this manuscript. Data from ENCODE Project was used and is available at https://www.encodeproject.org, IDs are shown in Table S3 in Supplementary File 1.

1. 1. van Schaik T
   2. Vos M
   3. Peric-Hupkes D
   4. Hn Celie P
   5. van Steensel B

   (2020) 4D Nucleome Data Portal

   ID f1218a92-1f37-4519-85d6-ccedd5f7ad39. Dara from: Cell cycle dynamics of lamina-associated DNA.

   https://data.4dnucleome.org/publications/f1218a92-1f37-4519-85d6-ccedd5f7ad39

1. 1. Albalat R
   2. Cañestro C

   (2016) Evolution by gene loss

   Nature Reviews. Genetics 17:379–391.

   https://doi.org/10.1038/nrg.2016.39

   • PubMed
   • Google Scholar
2. 1. Altschul SF
   2. Madden TL
   3. Schäffer AA
   4. Zhang J
   5. Zhang Z
   6. Miller W
   7. Lipman DJ

   (1997) Gapped blast and PSI-BLAST: a new generation of protein database search programs

   Nucleic Acids Research 25:3389–3402.

   https://doi.org/10.1093/nar/25.17.3389

   • PubMed
   • Google Scholar
3. 1. Bartha I
   2. di Iulio J
   3. Venter JC
   4. Telenti A

   (2018) Human gene essentiality

   Nature Reviews. Genetics 19:51–62.

   https://doi.org/10.1038/nrg.2017.75

   • PubMed
   • Google Scholar
4. 1. Bastian FB
   2. Roux J
   3. Niknejad A
   4. Comte A
   5. Fonseca Costa SS
   6. de Farias TM
   7. Moretti S
   8. Parmentier G
   9. de Laval VR
   10. Rosikiewicz M
   11. Wollbrett J
   12. Echchiki A
   13. Escoriza A
   14. Gharib WH
   15. Gonzales-Porta M
   16. Jarosz Y
   17. Laurenczy B
   18. Moret P
   19. Person E
   20. Roelli P
   21. Sanjeev K
   22. Seppey M
   23. Robinson-Rechavi M

   (2021) The bgee suite: integrated curated expression atlas and comparative transcriptomics in animals

   Nucleic Acids Research 49:D831–D847.

   https://doi.org/10.1093/nar/gkaa793

   • PubMed
   • Google Scholar
5. 1. Bernardi G
   2. Bernardi G

   (1986) Compositional constraints and genome evolution

   Journal of Molecular Evolution 24:1–11.

   https://doi.org/10.1007/BF02099946

   • PubMed
   • Google Scholar
6. 1. Blanc G
   2. Agarkova I
   3. Grimwood J
   4. Kuo A
   5. Brueggeman A
   6. Dunigan DD
   7. Gurnon J
   8. Ladunga I
   9. Lindquist E
   10. Lucas S
   11. Pangilinan J
   12. Pröschold T
   13. Salamov A
   14. Schmutz J
   15. Weeks D
   16. Yamada T
   17. Lomsadze A
   18. Borodovsky M
   19. Claverie J-M
   20. Grigoriev IV
   21. Van Etten JL

   (2012) The genome of the polar eukaryotic microalga coccomyxa subellipsoidea reveals traits of cold adaptation

   Genome Biology 13:R39.

   https://doi.org/10.1186/gb-2012-13-5-r39

   • Google Scholar
7. 1. Blomme T
   2. Vandepoele K
   3. De Bodt S
   4. Simillion C
   5. Maere S
   6. Van de Peer Y

   (2006) The gain and loss of genes during 600 million years of vertebrate evolution

   Genome Biology 7:R43.

   https://doi.org/10.1186/gb-2006-7-5-r43

   • PubMed
   • Google Scholar
8. 1. Botero-Castro F
   2. Figuet E
   3. Tilak MK
   4. Nabholz B
   5. Galtier N

   (2017) Avian genomes revisited: hidden genes uncovered and the rates versus traits paradox in birds

   Molecular Biology and Evolution 34:3123–3131.

   https://doi.org/10.1093/molbev/msx236

   • PubMed
   • Google Scholar
9. 1. Byrne KP
   2. Wolfe KH

   (2007) Consistent patterns of rate asymmetry and gene loss indicate widespread neofunctionalization of yeast genes after whole-genome duplication

   Genetics 175:1341–1350.

   https://doi.org/10.1534/genetics.106.066951

   • PubMed
   • Google Scholar
10. 1. Campbell CD
    2. Eichler EE

    (2013) Properties and rates of germline mutations in humans

    Trends in Genetics 29:575–584.

    https://doi.org/10.1016/j.tig.2013.04.005

    • PubMed
    • Google Scholar
11. 1. Cao J
    2. O’Day DR
    3. Pliner HA
    4. Kingsley PD
    5. Deng M
    6. Daza RM
    7. Zager MA
    8. Aldinger KA
    9. Blecher-Gonen R
    10. Zhang F
    11. Spielmann M
    12. Palis J
    13. Doherty D
    14. Steemers FJ
    15. Glass IA
    16. Trapnell C
    17. Shendure J

    (2020) A human cell atlas of fetal gene expression

    Science 370:eaba7721.

    https://doi.org/10.1126/science.aba7721

    • PubMed
    • Google Scholar
12. 1. Capella-Gutiérrez S
    2. Silla-Martínez JM
    3. Gabaldón T

    (2009) TrimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses

    Bioinformatics 25:1972–1973.

    https://doi.org/10.1093/bioinformatics/btp348

    • PubMed
    • Google Scholar
13. 1. Cherry JL

    (2010) Expression level, evolutionary rate, and the cost of expression

    Genome Biology and Evolution 2:757–769.

    https://doi.org/10.1093/gbe/evq059

    • PubMed
    • Google Scholar
14. 1. Choi JY
    2. Lee YCG

    (2020) Double-Edged sword: the evolutionary consequences of the epigenetic silencing of transposable elements

    PLOS Genetics 16:e1008872.

    https://doi.org/10.1371/journal.pgen.1008872

    • PubMed
    • Google Scholar
15. 1. Cohen N
    2. Dagan T
    3. Stone L
    4. Graur D

    (2005) Gc composition of the human genome: in search of isochores

    Molecular Biology and Evolution 22:1260–1272.

    https://doi.org/10.1093/molbev/msi115

    • PubMed
    • Google Scholar
16. 1. Cortez D
    2. Marin R
    3. Toledo-Flores D
    4. Froidevaux L
    5. Liechti A
    6. Waters PD
    7. Grützner F
    8. Kaessmann H

    (2014) Origins and functional evolution of Y chromosomes across mammals

    Nature 508:488–493.

    https://doi.org/10.1038/nature13151

    • PubMed
    • Google Scholar
17. 1. Debatisse M
    2. Le Tallec B
    3. Letessier A
    4. Dutrillaux B
    5. Brison O

    (2012) Common fragile sites: mechanisms of instability revisited

    Trends in Genetics 28:22–32.

    https://doi.org/10.1016/j.tig.2011.10.003

    • PubMed
    • Google Scholar
18. 1. Deutekom ES
    2. Vosseberg J
    3. van Dam TJP
    4. Snel B

    (2019) Measuring the impact of gene prediction on gene loss estimates in eukaryotes by quantifying falsely inferred absences

    PLOS Computational Biology 15:e1007301.

    https://doi.org/10.1371/journal.pcbi.1007301

    • PubMed
    • Google Scholar
19. 1. Emms DM
    2. Kelly S

    (2019) OrthoFinder: phylogenetic orthology inference for comparative genomics

    Genome Biology 20:238.

    https://doi.org/10.1186/s13059-019-1832-y

    • PubMed
    • Google Scholar
20. 1. Fernández R
    2. Gabaldón T

    (2020) Gene gain and loss across the metazoan tree of life

    Nature Ecology & Evolution 4:524–533.

    https://doi.org/10.1038/s41559-019-1069-x

    • PubMed
    • Google Scholar
21. 1. Fiston-Lavier AS
    2. Anxolabehere D
    3. Quesneville H

    (2007) A model of segmental duplication formation in Drosophila melanogaster

    Genome Research 17:1458–1470.

    https://doi.org/10.1101/gr.6208307

    • PubMed
    • Google Scholar
22. 1. Fu L
    2. Niu B
    3. Zhu Z
    4. Wu S
    5. Li W

    (2012) CD-HIT: accelerated for clustering the next-generation sequencing data

    Bioinformatics 28:3150–3152.

    https://doi.org/10.1093/bioinformatics/bts565

    • Google Scholar
23. 1. Giaever G
    2. Nislow C

    (2014) The yeast deletion collection: a decade of functional genomics

    Genetics 197:451–465.

    https://doi.org/10.1534/genetics.114.161620

    • PubMed
    • Google Scholar
24. 1. Gilbert N
    2. Boyle S
    3. Fiegler H
    4. Woodfine K
    5. Carter NP
    6. Bickmore WA

    (2004) Chromatin architecture of the human genome: gene-rich domains are enriched in open chromatin fibers

    Cell 118:555–566.

    https://doi.org/10.1016/j.cell.2004.08.011

    • PubMed
    • Google Scholar
25. 1. Gout JF
    2. Kahn D
    3. Duret L
    4. Pritchard JK
    5. Paramecium Post-Genomics Consortium

    (2010) The relationship among gene expression, the evolution of gene dosage, and the rate of protein evolution

    PLOS Genetics 6:e1000944.

    https://doi.org/10.1371/journal.pgen.1000944

    • PubMed
    • Google Scholar
26. 1. Grewal SIS
    2. Jia S

    (2007) Heterochromatin revisited

    Nature Reviews. Genetics 8:35–46.

    https://doi.org/10.1038/nrg2008

    • PubMed
    • Google Scholar
27. 1. Groenen MAM
    2. Wahlberg P
    3. Foglio M
    4. Cheng HH
    5. Megens HJ
    6. Crooijmans R
    7. Besnier F
    8. Lathrop M
    9. Muir WM
    10. Wong GKS
    11. Gut I
    12. Andersson L

    (2009) A high-density SNP-based linkage map of the chicken genome reveals sequence features correlated with recombination rate

    Genome Research 19:510–519.

    https://doi.org/10.1101/gr.086538.108

    • PubMed
    • Google Scholar
28. 1. Guijarro-Clarke C
    2. Holland PWH
    3. Paps J

    (2020) Widespread patterns of gene loss in the evolution of the animal Kingdom

    Nature Ecology & Evolution 4:519–523.

    https://doi.org/10.1038/s41559-020-1129-2

    • Google Scholar
29. 1. Hansen RS
    2. Thomas S
    3. Sandstrom R
    4. Canfield TK
    5. Thurman RE
    6. Weaver M
    7. Dorschner MO
    8. Gartler SM
    9. Stamatoyannopoulos JA

    (2010) Sequencing newly replicated DNA reveals widespread plasticity in human replication timing

    PNAS 107:139–144.

    https://doi.org/10.1073/pnas.0912402107

    • PubMed
    • Google Scholar
30. 1. Hao Z
    2. Lv D
    3. Ge Y
    4. Shi J
    5. Weijers D
    6. Yu G
    7. Chen J

    (2020) rideogram: drawing SVG graphics to visualize and MAP genome-wide data on the idiograms

    PeerJ. Computer Science 6:e251.

    https://doi.org/10.7717/peerj-cs.251

    • PubMed
    • Google Scholar
31. 1. Hara Y
    2. Tatsumi K
    3. Yoshida M
    4. Kajikawa E
    5. Kiyonari H
    6. Kuraku S

    (2015) Optimizing and benchmarking de novo transcriptome sequencing: from library preparation to assembly evaluation

    BMC Genomics 16:977.

    https://doi.org/10.1186/s12864-015-2007-1

    • Google Scholar
32. 1. Hara Y
    2. Takeuchi M
    3. Kageyama Y
    4. Tatsumi K
    5. Hibi M
    6. Kiyonari H
    7. Kuraku S

    (2018a) Madagascar ground gecko genome analysis characterizes asymmetric fates of duplicated genes

    BMC Biology 16:40.

    https://doi.org/10.1186/s12915-018-0509-4

    • PubMed
    • Google Scholar
33. 1. Hara Y
    2. Yamaguchi K
    3. Onimaru K
    4. Kadota M
    5. Koyanagi M
    6. Keeley SD
    7. Tatsumi K
    8. Tanaka K
    9. Motone F
    10. Kageyama Y
    11. Nozu R
    12. Adachi N
    13. Nishimura O
    14. Nakagawa R
    15. Tanegashima C
    16. Kiyatake I
    17. Matsumoto R
    18. Murakumo K
    19. Nishida K
    20. Terakita A
    21. Kuratani S
    22. Sato K
    23. Hyodo S
    24. Kuraku S

    (2018b) Shark genomes provide insights into elasmobranch evolution and the origin of vertebrates

    Nature Ecology & Evolution 2:1761–1771.

    https://doi.org/10.1038/s41559-018-0673-5

    • PubMed
    • Google Scholar
34. Software

    1. Hara Y

    (2022) Elusivegenes, version swh:1:rev:4c4d279c77c838ec6acb86032bdb51c514c3cc60

    Software Heritage.

    https://archive.softwareheritage.org/swh:1:dir:67291d403a94399fcb9479be45d35fee4585becc;origin=https://github.com/yuichiroharajpn/ElusiveGenes;visit=swh:1:snp:d53bb1fab81dbd98069b50871fd11e474fe72156;anchor=swh:1:rev:4c4d279c77c838ec6acb86032bdb51c514c3cc60
35. 1. Helmrich A
    2. Stout-Weider K
    3. Hermann K
    4. Schrock E
    5. Heiden T

    (2006) Common fragile sites are conserved features of human and mouse chromosomes and relate to large active genes

    Genome Research 16:1222–1230.

    https://doi.org/10.1101/gr.5335506

    • PubMed
    • Google Scholar
36. 1. Hinrichs AS
    2. Karolchik D
    3. Baertsch R
    4. Barber GP
    5. Bejerano G
    6. Clawson H
    7. Diekhans M
    8. Furey TS
    9. Harte RA
    10. Hsu F
    11. Hillman-Jackson J
    12. Kuhn RM
    13. Pedersen JS
    14. Pohl A
    15. Raney BJ
    16. Rosenbloom KR
    17. Siepel A
    18. Smith KE
    19. Sugnet CW
    20. Sultan-Qurraie A
    21. Thomas DJ
    22. Trumbower H
    23. Weber RJ
    24. Weirauch M
    25. Zweig AS
    26. Haussler D
    27. Kent WJ

    (2006) The UCSC genome browser database: update 2006

    Nucleic Acids Research 34:D590–D598.

    https://doi.org/10.1093/nar/gkj144

    • PubMed
    • Google Scholar
37. 1. Hirsh AE
    2. Fraser HB

    (2001) Protein dispensability and rate of evolution

    Nature 411:1046–1049.

    https://doi.org/10.1038/35082561

    • PubMed
    • Google Scholar
38. 1. Huerta-Cepas J
    2. Dopazo H
    3. Dopazo J
    4. Gabaldón T

    (2007) The human phylome

    Genome Biology 8:R109.

    https://doi.org/10.1186/gb-2007-8-6-r109

    • PubMed
    • Google Scholar
39. 1. Huerta-Cepas J
    2. Serra F
    3. Bork P

    (2016) Ete 3: reconstruction, analysis, and visualization of phylogenomic data

    Molecular Biology and Evolution 33:1635–1638.

    https://doi.org/10.1093/molbev/msw046

    • PubMed
    • Google Scholar
40. 1. Hughes JF
    2. Skaletsky H
    3. Brown LG
    4. Pyntikova T
    5. Graves T
    6. Fulton RS
    7. Dugan S
    8. Ding Y
    9. Buhay CJ
    10. Kremitzki C
    11. Wang Q
    12. Shen H
    13. Holder M
    14. Villasana D
    15. Nazareth LV
    16. Cree A
    17. Courtney L
    18. Veizer J
    19. Kotkiewicz H
    20. Cho TJ
    21. Koutseva N
    22. Rozen S
    23. Muzny DM
    24. Warren WC
    25. Gibbs RA
    26. Wilson RK
    27. Page DC

    (2012) Strict evolutionary conservation followed rapid gene loss on human and rhesus Y chromosomes

    Nature 483:82–86.

    https://doi.org/10.1038/nature10843

    • PubMed
    • Google Scholar
41. 1. Imbeault M
    2. Helleboid PY
    3. Trono D

    (2017) Krab zinc-finger proteins contribute to the evolution of gene regulatory networks

    Nature 543:550–554.

    https://doi.org/10.1038/nature21683

    • PubMed
    • Google Scholar
42. 1. International Chicken Genome Sequencing Consortium

    (2004) Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution

    Nature 432:695–716.

    https://doi.org/10.1038/nature03154

    • PubMed
    • Google Scholar
43. 1. Jordan IK
    2. Rogozin IB
    3. Wolf YI
    4. Koonin EV

    (2002) Essential genes are more evolutionarily conserved than are nonessential genes in bacteria

    Genome Research 12:962–968.

    https://doi.org/10.1101/gr.87702

    • PubMed
    • Google Scholar
44. 1. Kalyaanamoorthy S
    2. Minh BQ
    3. Wong TKF
    4. von Haeseler A
    5. Jermiin LS

    (2017) ModelFinder: fast model selection for accurate phylogenetic estimates

    Nature Methods 14:587–589.

    https://doi.org/10.1038/nmeth.4285

    • PubMed
    • Google Scholar
45. 1. Karczewski KJ
    2. Francioli LC
    3. Tiao G
    4. Cummings BB
    5. Alföldi J
    6. Wang Q
    7. Collins RL
    8. Laricchia KM
    9. Ganna A
    10. Birnbaum DP
    11. Gauthier LD
    12. Brand H
    13. Solomonson M
    14. Watts NA
    15. Rhodes D
    16. Singer-Berk M
    17. England EM
    18. Seaby EG
    19. Kosmicki JA
    20. Walters RK
    21. Tashman K
    22. Farjoun Y
    23. Banks E
    24. Poterba T
    25. Wang A
    26. Seed C
    27. Whiffin N
    28. Chong JX
    29. Samocha KE
    30. Pierce-Hoffman E
    31. Zappala Z
    32. O’Donnell-Luria AH
    33. Minikel EV
    34. Weisburd B
    35. Lek M
    36. Ware JS
    37. Vittal C
    38. Armean IM
    39. Bergelson L
    40. Cibulskis K
    41. Connolly KM
    42. Covarrubias M
    43. Donnelly S
    44. Ferriera S
    45. Gabriel S
    46. Gentry J
    47. Gupta N
    48. Jeandet T
    49. Kaplan D
    50. Llanwarne C
    51. Munshi R
    52. Novod S
    53. Petrillo N
    54. Roazen D
    55. Ruano-Rubio V
    56. Saltzman A
    57. Schleicher M
    58. Soto J
    59. Tibbetts K
    60. Tolonen C
    61. Wade G
    62. Talkowski ME
    63. Genome Aggregation Database Consortium
    64. Neale BM
    65. Daly MJ
    66. MacArthur DG

    (2021) Author correction: the mutational constraint spectrum quantified from variation in 141,456 humans

    Nature 590:E53.

    https://doi.org/10.1038/s41586-020-03174-8

    • PubMed
    • Google Scholar
46. 1. Katoh K
    2. Standley DM

    (2013) MAFFT multiple sequence alignment software version 7: improvements in performance and usability

    Molecular Biology and Evolution 30:772–780.

    https://doi.org/10.1093/molbev/mst010

    • PubMed
    • Google Scholar
47. 1. Katzman S
    2. Capra JA
    3. Haussler D
    4. Pollard KS

    (2011) Ongoing GC-biased evolution is widespread in the human genome and enriched near recombination hot spots

    Genome Biology and Evolution 3:614–626.

    https://doi.org/10.1093/gbe/evr058

    • PubMed
    • Google Scholar
48. 1. Koren A
    2. Polak P
    3. Nemesh J
    4. Michaelson JJ
    5. Sebat J
    6. Sunyaev SR
    7. McCarroll SA

    (2012) Differential relationship of DNA replication timing to different forms of human mutation and variation

    The American Journal of Human Genetics 91:1033–1040.

    https://doi.org/10.1016/j.ajhg.2012.10.018

    • Google Scholar
49. 1. Korenberg JR
    2. Rykowski MC

    (1988) Human genome organization: Alu, lines, and the molecular structure of metaphase chromosome bands

    Cell 53:391–400.

    https://doi.org/10.1016/0092-8674(88)90159-6

    • Google Scholar
50. 1. Krylov DM
    2. Wolf YI
    3. Rogozin IB
    4. Koonin EV

    (2003) Gene loss, protein sequence divergence, gene dispensability, expression level, and interactivity are correlated in eukaryotic evolution

    Genome Research 13:2229–2235.

    https://doi.org/10.1101/gr.1589103

    • PubMed
    • Google Scholar
51. 1. Lewin TD
    2. Royall AH
    3. Holland PWH

    (2021) Dynamic molecular evolution of mammalian homeobox genes: duplication, loss, divergence and gene conversion sculpt PRD class repertoires

    Journal of Molecular Evolution 89:396–414.

    https://doi.org/10.1007/s00239-021-10012-6

    • PubMed
    • Google Scholar
52. 1. Liu G
    2. Yong MYJ
    3. Yurieva M
    4. Srinivasan KG
    5. Liu J
    6. Lim JSY
    7. Poidinger M
    8. Wright GD
    9. Zolezzi F
    10. Choi H
    11. Pavelka N
    12. Rancati G

    (2015) Gene essentiality is a quantitative property linked to cellular evolvability

    Cell 163:1388–1399.

    https://doi.org/10.1016/j.cell.2015.10.069

    • Google Scholar
53. 1. MacDonald JR
    2. Ziman R
    3. Yuen RKC
    4. Feuk L
    5. Scherer SW

    (2014) The database of genomic variants: a curated collection of structural variation in the human genome

    Nucleic Acids Research 42:D986–D992.

    https://doi.org/10.1093/nar/gkt958

    • Google Scholar
54. 1. Maclean CJ
    2. Metzger BPH
    3. Yang JR
    4. Ho WC
    5. Moyers B
    6. Zhang J

    (2017) Deciphering the genic basis of yeast fitness variation by simultaneous forward and reverse genetics

    Molecular Biology and Evolution 34:2486–2502.

    https://doi.org/10.1093/molbev/msx151

    • Google Scholar
55. 1. Maeso I
    2. Dunwell TL
    3. Wyatt CDR
    4. Marlétaz F
    5. Vető B
    6. Bernal JA
    7. Quah S
    8. Irimia M
    9. Holland PWH

    (2016) Evolutionary origin and functional divergence of totipotent cell homeobox genes in eutherian mammals

    BMC Biology 14:45.

    https://doi.org/10.1186/s12915-016-0267-0

    • PubMed
    • Google Scholar
56. 1. McQueen HA
    2. Siriaco G
    3. Bird AP

    (1998) Chicken microchromosomes are hyperacetylated, early replicating, and gene rich

    Genome Research 8:621–630.

    https://doi.org/10.1101/gr.8.6.621

    • PubMed
    • Google Scholar
57. 1. Medstrand P
    2. van de Lagemaat LN
    3. Mager DL

    (2002) Retroelement distributions in the human genome: variations associated with age and proximity to genes

    Genome Research 12:1483–1495.

    https://doi.org/10.1101/gr.388902

    • PubMed
    • Google Scholar
58. 1. Miyata T
    2. Yasunaga T
    3. Nishida T

    (1980) Nucleotide sequence divergence and functional constraint in mRNA evolution

    PNAS 77:7328–7332.

    https://doi.org/10.1073/pnas.77.12.7328

    • Google Scholar
59. 1. Monroe JG
    2. Srikant T
    3. Carbonell-Bejerano P
    4. Becker C
    5. Lensink M
    6. Exposito-Alonso M
    7. Klein M
    8. Hildebrandt J
    9. Neumann M
    10. Kliebenstein D
    11. Weng ML
    12. Imbert E
    13. Ågren J
    14. Rutter MT
    15. Fenster CB
    16. Weigel D

    (2022) Mutation bias reflects natural selection in Arabidopsis thaliana

    Nature 602:101–105.

    https://doi.org/10.1038/s41586-021-04269-6

    • PubMed
    • Google Scholar
60. 1. Moyers BA
    2. Zhang J

    (2015) Phylostratigraphic bias creates spurious patterns of genome evolution

    Molecular Biology and Evolution 32:258–267.

    https://doi.org/10.1093/molbev/msu286

    • PubMed
    • Google Scholar
61. 1. Nakatani Y
    2. Shingate P
    3. Ravi V
    4. Pillai NE
    5. Prasad A
    6. McLysaght A
    7. Venkatesh B

    (2021) Reconstruction of proto-vertebrate, proto-cyclostome and proto-gnathostome genomes provides new insights into early vertebrate evolution

    Nature Communications 12:4489.

    https://doi.org/10.1038/s41467-021-24573-z

    • PubMed
    • Google Scholar
62. Book

    1. Nei M.

    (2013)

    Mutation-Driven Evolution

    Oxford University Press.

    • Google Scholar
63. 1. Nguyen L-T
    2. Schmidt HA
    3. von Haeseler A
    4. Minh BQ

    (2015) IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies

    Molecular Biology and Evolution 32:268–274.

    https://doi.org/10.1093/molbev/msu300

    • PubMed
    • Google Scholar
64. 1. Nurk S
    2. Koren S
    3. Rhie A
    4. Rautiainen M
    5. Bzikadze AV
    6. Mikheenko A
    7. Vollger MR
    8. Altemose N
    9. Uralsky L
    10. Gershman A
    11. Aganezov S
    12. Hoyt SJ
    13. Diekhans M
    14. Logsdon GA
    15. Alonge M
    16. Antonarakis SE
    17. Borchers M
    18. Bouffard GG
    19. Brooks SY
    20. Caldas GV
    21. Chen N-C
    22. Cheng H
    23. Chin C-S
    24. Chow W
    25. de Lima LG
    26. Dishuck PC
    27. Durbin R
    28. Dvorkina T
    29. Fiddes IT
    30. Formenti G
    31. Fulton RS
    32. Fungtammasan A
    33. Garrison E
    34. Grady PGS
    35. Graves-Lindsay TA
    36. Hall IM
    37. Hansen NF
    38. Hartley GA
    39. Haukness M
    40. Howe K
    41. Hunkapiller MW
    42. Jain C
    43. Jain M
    44. Jarvis ED
    45. Kerpedjiev P
    46. Kirsche M
    47. Kolmogorov M
    48. Korlach J
    49. Kremitzki M
    50. Li H
    51. Maduro VV
    52. Marschall T
    53. McCartney AM
    54. McDaniel J
    55. Miller DE
    56. Mullikin JC
    57. Myers EW
    58. Olson ND
    59. Paten B
    60. Peluso P
    61. Pevzner PA
    62. Porubsky D
    63. Potapova T
    64. Rogaev EI
    65. Rosenfeld JA
    66. Salzberg SL
    67. Schneider VA
    68. Sedlazeck FJ
    69. Shafin K
    70. Shew CJ
    71. Shumate A
    72. Sims Y
    73. Smit AFA
    74. Soto DC
    75. Sović I
    76. Storer JM
    77. Streets A
    78. Sullivan BA
    79. Thibaud-Nissen F
    80. Torrance J
    81. Wagner J
    82. Walenz BP
    83. Wenger A
    84. Wood JMD
    85. Xiao C
    86. Yan SM
    87. Young AC
    88. Zarate S
    89. Surti U
    90. McCoy RC
    91. Dennis MY
    92. Alexandrov IA
    93. Gerton JL
    94. O’Neill RJ
    95. Timp W
    96. Zook JM
    97. Schatz MC
    98. Eichler EE
    99. Miga KH
    100. Phillippy AM

    (2022) The complete sequence of a human genome

    Science 376:44–53.

    https://doi.org/10.1126/science.abj6987

    • PubMed
    • Google Scholar
65. 1. O’Geen H
    2. Squazzo SL
    3. Iyengar S
    4. Blahnik K
    5. Rinn JL
    6. Chang HY
    7. Green R
    8. Farnham PJ

    (2007) Genome-Wide analysis of KAP1 binding suggests autoregulation of KRAB-znfs

    PLOS Genetics 3:e89.

    https://doi.org/10.1371/journal.pgen.0030089

    • PubMed
    • Google Scholar
66. 1. Olson MV

    (1999) When less is more: gene loss as an engine of evolutionary change

    American Journal of Human Genetics 64:18–23.

    https://doi.org/10.1086/302219

    • PubMed
    • Google Scholar
67. 1. Pál C
    2. Papp B
    3. Lercher MJ

    (2006) An integrated view of protein evolution

    Nature Reviews. Genetics 7:337–348.

    https://doi.org/10.1038/nrg1838

    • PubMed
    • Google Scholar
68. 1. Paradis E
    2. Schliep K

    (2019) Ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R

    Bioinformatics 35:526–528.

    https://doi.org/10.1093/bioinformatics/bty633

    • PubMed
    • Google Scholar
69. 1. Perry BW
    2. Schield DR
    3. Adams RH
    4. Castoe TA
    5. Arkhipova I

    (2020) Microchromosomes exhibit distinct features of vertebrate chromosome structure and function with underappreciated ramifications for genome evolution

    Molecular Biology and Evolution 38:904–910.

    https://doi.org/10.1093/molbev/msaa253

    • Google Scholar
70. 1. Rangasamy D

    (2013) Distinctive patterns of epigenetic marks are associated with promoter regions of mouse LINE-1 and LTR retrotransposons

    Mobile DNA 4:27.

    https://doi.org/10.1186/1759-8753-4-27

    • PubMed
    • Google Scholar
71. 1. Rao SSP
    2. Huntley MH
    3. Durand NC
    4. Stamenova EK
    5. Bochkov ID
    6. Robinson JT
    7. Sanborn AL
    8. Machol I
    9. Omer AD
    10. Lander ES
    11. Aiden EL

    (2014) A 3D map of the human genome at kilobase resolution reveals principles of chromatin looping

    Cell 159:1665–1680.

    https://doi.org/10.1016/j.cell.2014.11.021

    • PubMed
    • Google Scholar
72. 1. Rice AM
    2. McLysaght A

    (2017) Dosage sensitivity is a major determinant of human copy number variant pathogenicity

    Nature Communications 8:14366.

    https://doi.org/10.1038/ncomms14366

    • PubMed
    • Google Scholar
73. 1. Roux J
    2. Liu J
    3. Robinson-Rechavi M

    (2017) Selective constraints on coding sequences of nervous system genes are a major determinant of duplicate gene retention in vertebrates

    Molecular Biology and Evolution 34:2773–2791.

    https://doi.org/10.1093/molbev/msx199

    • PubMed
    • Google Scholar
74. 1. Schaibley VM
    2. Zawistowski M
    3. Wegmann D
    4. Ehm MG
    5. Nelson MR
    6. St Jean PL
    7. Abecasis GR
    8. Novembre J
    9. Zöllner S
    10. Li JZ

    (2013) The influence of genomic context on mutation patterns in the human genome inferred from rare variants

    Genome Research 23:1974–1984.

    https://doi.org/10.1101/gr.154971.113

    • PubMed
    • Google Scholar
75. 1. Schield DR
    2. Card DC
    3. Hales NR
    4. Perry BW
    5. Pasquesi GM
    6. Blackmon H
    7. Adams RH
    8. Corbin AB
    9. Smith CF
    10. Ramesh B
    11. Demuth JP
    12. Betrán E
    13. Tollis M
    14. Meik JM
    15. Mackessy SP
    16. Castoe TA

    (2019) The origins and evolution of chromosomes, dosage compensation, and mechanisms underlying venom regulation in snakes

    Genome Research 29:590–601.

    https://doi.org/10.1101/gr.240952.118

    • PubMed
    • Google Scholar
76. 1. Seplyarskiy VB
    2. Sunyaev S

    (2021) The origin of human mutation in light of genomic data

    Nature Reviews Genetics 22:672–686.

    https://doi.org/10.1038/s41576-021-00376-2

    • Google Scholar
77. 1. Sharma V
    2. Hecker N
    3. Roscito JG
    4. Foerster L
    5. Langer BE
    6. Hiller M

    (2018) A genomics approach reveals insights into the importance of gene losses for mammalian adaptations

    Nature Communications 9:1215.

    https://doi.org/10.1038/s41467-018-03667-1

    • Google Scholar
78. 1. Shen XX
    2. Opulente DA
    3. Kominek J
    4. Zhou X
    5. Steenwyk JL
    6. Buh KV
    7. Haase MAB
    8. Wisecaver JH
    9. Wang M
    10. Doering DT
    11. Boudouris JT
    12. Schneider RM
    13. Langdon QK
    14. Ohkuma M
    15. Endoh R
    16. Takashima M
    17. Manabe RI
    18. Čadež N
    19. Libkind D
    20. Rosa CA
    21. DeVirgilio J
    22. Hulfachor AB
    23. Groenewald M
    24. Kurtzman CP
    25. Hittinger CT
    26. Rokas A

    (2018) Tempo and mode of genome evolution in the budding yeast subphylum

    Cell 175:1533–1545.

    https://doi.org/10.1016/j.cell.2018.10.023

    • PubMed
    • Google Scholar
79. 1. Sherry ST
    2. Ward MH
    3. Kholodov M
    4. Baker J
    5. Phan L
    6. Smigielski EM
    7. Sirotkin K

    (2001) DbSNP: the NCBI database of genetic variation

    Nucleic Acids Research 29:308–311.

    https://doi.org/10.1093/nar/29.1.308

    • PubMed
    • Google Scholar
80. 1. Simão FA
    2. Waterhouse RM
    3. Ioannidis P
    4. Kriventseva EV
    5. Zdobnov EM

    (2015) BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs

    Bioinformatics 31:3210–3212.

    https://doi.org/10.1093/bioinformatics/btv351

    • PubMed
    • Google Scholar
81. 1. Slotkin RK
    2. Martienssen R

    (2007) Transposable elements and the epigenetic regulation of the genome

    Nature Reviews Genetics 8:272–285.

    https://doi.org/10.1038/nrg2072

    • Google Scholar
82. 1. Stamatoyannopoulos JA
    2. Adzhubei I
    3. Thurman RE
    4. Kryukov GV
    5. Mirkin SM
    6. Sunyaev SR

    (2009) Human mutation rate associated with DNA replication timing

    Nature Genetics 41:393–395.

    https://doi.org/10.1038/ng.363

    • PubMed
    • Google Scholar
83. 1. Steinegger M
    2. Söding J

    (2017) MMseqs2 enables sensitive protein sequence searching for the analysis of massive data sets

    Nature Biotechnology 35:1026–1028.

    https://doi.org/10.1038/nbt.3988

    • Google Scholar
84. 1. Terekhanova NV
    2. Seplyarskiy VB
    3. Soldatov RA
    4. Bazykin GA

    (2017) Evolution of local mutation rate and its determinants

    Molecular Biology and Evolution 34:msx060.

    https://doi.org/10.1093/molbev/msx060

    • Google Scholar
85. 1. The ENCODE Project Consortium

    (2012) An integrated encyclopedia of DNA elements in the human genome

    Nature 489:57–74.

    https://doi.org/10.1038/nature11247

    • Google Scholar
86. 1. The GTEx Consortium
    2. Aguet F
    3. Anand S
    4. Ardlie KG
    5. Gabriel S
    6. Getz GA
    7. Graubert A
    8. Hadley K
    9. Handsaker RE
    10. Huang KH
    11. Kashin S
    12. Li X
    13. MacArthur DG
    14. Meier SR
    15. Nedzel JL
    16. Nguyen DT
    17. Segrè AV
    18. Todres E
    19. Balliu B
    20. Barbeira AN
    21. Battle A
    22. Bonazzola R
    23. Brown A
    24. Brown CD
    25. Castel SE
    26. Conrad DF
    27. Cotter DJ
    28. Cox N
    29. Das S
    30. de Goede OM
    31. Dermitzakis ET
    32. Einson J
    33. Engelhardt BE
    34. Eskin E
    35. Eulalio TY
    36. Ferraro NM
    37. Flynn ED
    38. Fresard L
    39. Gamazon ER
    40. Garrido-Martín D
    41. Gay NR
    42. Gloudemans MJ
    43. Guigó R
    44. Hame AR
    45. He Y
    46. Hoffman PJ
    47. Hormozdiari F
    48. Hou L
    49. Im HK
    50. Jo B
    51. Kasela S
    52. Kellis M
    53. Kim-Hellmuth S
    54. Kwong A
    55. Lappalainen T
    56. Li X
    57. Liang Y
    58. Mangul S
    59. Mohammadi P
    60. Montgomery SB
    61. Muñoz-Aguirre M
    62. Nachun DC
    63. Nobel AB
    64. Oliva M
    65. Park Y
    66. Park Y
    67. Parsana P
    68. Rao AS
    69. Reverter F
    70. Rouhana JM
    71. Sabatti C
    72. Saha A
    73. Stephens M
    74. Stranger BE
    75. Strober BJ
    76. Teran NA
    77. Viñuela A
    78. Wang G
    79. Wen X
    80. Wright F
    81. Wucher V
    82. Zou Y
    83. Ferreira PG
    84. Li G
    85. Melé M
    86. Yeger-Lotem E
    87. Barcus ME
    88. Bradbury D
    89. Krubit T
    90. McLean JA
    91. Qi L
    92. Robinson K
    93. Roche NV
    94. Smith AM
    95. Sobin L
    96. Tabor DE
    97. Undale A
    98. Bridge J
    99. Brigham LE
    100. Foster BA
    101. Gillard BM
    102. Hasz R
    103. Hunter M
    104. Johns C
    105. Johnson M
    106. Karasik E
    107. Kopen G
    108. Leinweber WF
    109. McDonald A
    110. Moser MT
    111. Myer K
    112. Ramsey KD
    113. Roe B
    114. Shad S
    115. Thomas JA
    116. Walters G
    117. Washington M
    118. Wheeler J
    119. Jewell SD
    120. Rohrer DC
    121. Valley DR
    122. Davis DA
    123. Mash DC
    124. Branton PA
    125. Barker LK
    126. Gardiner HM
    127. Mosavel M
    128. Siminoff LA
    129. Flicek P
    130. Haeussler M
    131. Juettemann T
    132. Kent WJ
    133. Lee CM
    134. Powell CC
    135. Rosenbloom KR
    136. Ruffier M
    137. Sheppard D
    138. Taylor K
    139. Trevanion SJ
    140. Zerbino DR
    141. Abell NS
    142. Akey J
    143. Chen L
    144. Demanelis K
    145. Doherty JA
    146. Feinberg AP
    147. Hansen KD
    148. Hickey PF
    149. Jasmine F
    150. Jiang L
    151. Kaul R
    152. Kibriya MG
    153. Li JB
    154. Li Q
    155. Lin S
    156. Linder SE
    157. Pierce BL
    158. Rizzardi LF
    159. Skol AD
    160. Smith KS
    161. Snyder M
    162. Stamatoyannopoulos J
    163. Tang H
    164. Wang M
    165. Carithers LJ
    166. Guan P
    167. Koester SE
    168. Little AR
    169. Moore HM
    170. Nierras CR
    171. Rao AK
    172. Vaught JB
    173. Volpi S

    (2020) The gtex Consortium atlas of genetic regulatory effects across human tissues

    Science 369:1318–1330.

    https://doi.org/10.1126/science.aaz1776

    • Google Scholar
87. 1. Trapnell C
    2. Roberts A
    3. Goff L
    4. Pertea G
    5. Kim D
    6. Kelley DR
    7. Pimentel H
    8. Salzberg SL
    9. Rinn JL
    10. Pachter L

    (2012) Differential gene and transcript expression analysis of RNA-Seq experiments with tophat and cufflinks

    Nature Protocols 7:562–578.

    https://doi.org/10.1038/nprot.2012.016

    • Google Scholar
88. 1. Underwood CJ
    2. Henderson IR
    3. Martienssen RA

    (2017) Genetic and epigenetic variation of transposable elements in Arabidopsis

    Current Opinion in Plant Biology 36:135–141.

    https://doi.org/10.1016/j.pbi.2017.03.002

    • Google Scholar
89. 1. van Schaik T
    2. Vos M
    3. Peric-Hupkes D
    4. Hn Celie P
    5. van Steensel B

    (2020) Cell cycle dynamics of lamina-associated DNA

    EMBO Reports 21:e50636.

    https://doi.org/10.15252/embr.202050636

    • PubMed
    • Google Scholar
90. 1. van Steensel B
    2. Belmont AS

    (2017) Lamina-Associated domains: links with chromosome architecture, heterochromatin, and gene repression

    Cell 169:780–791.

    https://doi.org/10.1016/j.cell.2017.04.022

    • Google Scholar
91. 1. Vogel MJ
    2. Guelen L
    3. de Wit E
    4. Peric-Hupkes D
    5. Lodén M
    6. Talhout W
    7. Feenstra M
    8. Abbas B
    9. Classen A-K
    10. van Steensel B

    (2006) Human heterochromatin proteins form large domains containing KRAB-ZNF genes

    Genome Research 16:1493–1504.

    https://doi.org/10.1101/gr.5391806

    • PubMed
    • Google Scholar
92. 1. Waters PD
    2. Patel HR
    3. Ruiz-Herrera A
    4. Álvarez-González L
    5. Lister NC
    6. Simakov O
    7. Ezaz T
    8. Kaur P
    9. Frere C
    10. Grützner F
    11. Georges A
    12. Graves JAM

    (2021) Microchromosomes are building blocks of bird, reptile, and mammal chromosomes

    PNAS 118:e2112494118.

    https://doi.org/10.1073/pnas.2112494118

    • PubMed
    • Google Scholar
93. 1. Xie KT
    2. Wang G
    3. Thompson AC
    4. Wucherpfennig JI
    5. Reimchen TE
    6. MacColl ADC
    7. Schluter D
    8. Bell MA
    9. Vasquez KM
    10. Kingsley DM

    (2019) Dna fragility in the parallel evolution of pelvic reduction in stickleback fish

    Science 363:81–84.

    https://doi.org/10.1126/science.aan1425

    • PubMed
    • Google Scholar
94. 1. Yang J
    2. Gu Z
    3. Li WH

    (2003) Rate of protein evolution versus fitness effect of gene deletion

    Molecular Biology and Evolution 20:772–774.

    https://doi.org/10.1093/molbev/msg078

    • PubMed
    • Google Scholar
95. 1. Yang Z

    (2007) PAML 4: phylogenetic analysis by maximum likelihood

    Molecular Biology and Evolution 24:1586–1591.

    https://doi.org/10.1093/molbev/msm088

    • PubMed
    • Google Scholar
96. 1. Yoshihara M
    2. Jiang L
    3. Akatsuka S
    4. Suyama M
    5. Toyokuni S

    (2014) Genome-Wide profiling of 8-oxoguanine reveals its association with spatial positioning in nucleus

    DNA Research 21:603–612.

    https://doi.org/10.1093/dnares/dsu023

    • PubMed
    • Google Scholar
97. 1. Zhang J
    2. Yang JR

    (2015) Determinants of the rate of protein sequence evolution

    Nature Reviews. Genetics 16:409–420.

    https://doi.org/10.1038/nrg3950

    • PubMed
    • Google Scholar
98. 1. Zheng X
    2. Hu J
    3. Yue S
    4. Kristiani L
    5. Kim M
    6. Sauria M
    7. Taylor J
    8. Kim Y
    9. Zheng Y

    (2018) Lamins organize the global three-dimensional genome from the nuclear periphery

    Molecular Cell 71:802–815.

    https://doi.org/10.1016/j.molcel.2018.05.017

    • PubMed
    • Google Scholar

Decision letter after peer review:

Thank you for submitting your article "Gene fate spectrum as a reflection of local genomic properties" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, including Wenfeng Qian as the Reviewing Editor and Reviewer #1, and the evaluation has been overseen by George Perry as the Senior Editor.

The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.

Essential revisions:

1) The identification of "elusive genes" in the current manuscript requires additional scrutinization, given it is the foundation of the whole study. Please check published pipelines in the identification of gene losses (e.g., TOGA – https://github.com/hillerlab/TOGA) and use additional tools such as BLASTX search to test for known technical artifacts when calling genes (homology detection failure, refer to https://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.3000862). Please give a reason for the parameters used in the analysis (e.g., CD-Hit clustering) or examine if the conclusions remain supported by various parameters in the computational pipeline. Also, please take a look at the enrichment of elusive genes in human chr19, and use the synteny-based age estimation of the elusive genes (Shao et al. 2019).

2) Please also present the features of other genes (the genomic background other than elusive and non-elusive genes). Are these genes show intermediate patterns between those of elusive and non-elusive genes? Please edit the manuscript accordingly if the definition of non-elusiveness is actually equivalent to the genomic background.

3) It would be informative to test the links between recombination rate / LD and the genomic locations of elusive genes (compared against randomly sampled genes).

4) Please control confounding factors such as gene expression level and confirm whether the proxy of mutation rate (i.e., Ks) is actually confounded by gene importance.

5) Please consider extending the analyses on fish and birds to other genomic features.

6) The authors should reconcile the findings in this study with previous reports about microchromosomes.

7) Please consider improving the clarity/presentation of Figures 5 and 7 and examine whether the pattern remains robust using various parameter sets.

8) Please think of a better term than "elusive gene" to describe the genes that were lost independently in different lineages. Please also clearly define other terms in the manuscript (e.g., functionally indispensable vs. importance, are they the same concept?)

9) Please consider presenting the current study in the framework of mutation-selection balance, and better explain the novelty of the study over previous tremendous studies about gene losses.

Reviewer #1 (Recommendations for the authors):

Line 18. "neutral factor" is better replaced by "factors independent of gene dispensability".

Line 47. "However" should be "on the contrary"?

Lines 113-114. As indicated in the weaknesses part, I am not fully convinced these genomic, epigenomic, and transcriptomic features are completely independent of gene function.

Figure 1b. Define the red and orange crosses in the legend.

Figure 7 appeared first in the Discussion section. Can this part be moved to the Results section?

Reviewer #2 (Recommendations for the authors):

Overall, I believe that this is an interesting study. However, this version of the manuscript could be significantly improved in terms of logical depth and methodological stringency.

1. Authors actually support the concept of mutation-driven evolution, i.e., the high mutation rate in genomic regions harboring the elusive genes would predispose their fate toward death. To increase the significance of their work, I suggest authors cite (Nei 2013; Xie et al. 2019) and put their work in a bigger context.

2. Authors mentioned that elusive genes are less important and thus more prone to loss. In my view, pleiotropy is a better term compared to importance. That is, elusive genes are less pleiotropic [e.g. narrowly expressed, Figure 5] and thus their loss are more tolerable or easily compensated by other genes. Actually, narrow expression breadth has been observed to be correlated with gene loss in both humans and flies (MacArthur et al. 2012; Yang et al. 2015).

3. I am generally convinced that authors reliably identified elusive genes by identifying gene loss events in the common ancestor of multiple descendant species (to control for errors induced by assembly or annotation, Figure 1). However, Figure 7 shows the enrichment of elusive genes in human chr19. This chromosome is well known to be enriched with tandemly duplicated Krueppel-associated box C2H2 zinc-finger protein family (KZNF), many of which are primate-specific (Shao et al. 2019). I suspect that tree-based strategy implemented in Figure 1 could not be able to dissect the evolution of this super complex gene family. I am proposing two specific analyses: how many elusive genes encoded by chr19 are KZNFs? how many of them have Ensembl one-to-one orthologs across mammals?

4. With the patterns in Figure 3 and 7, authors argued that features of elusive genes are deeply ancient and could be inherited from the microchromosomes of early vertebrates. This statement has multiple problems.

a) Figure 3 only show genomic level features (e.g., high GC content) conserved in multiple vertebrates including shark and chicken. Epigenetic features analyzed in Figure 5 to 6 were only based on human data. I suggest authors to extend these analyses to shark or chicken. Although some epigenetic data could not be available for these species, transcriptome data analyzed in Figure 5 should be available for at least some species.

b) In Figure 7, authors propose the concept related with microchromosomes. These chromosomes have been extensively studied, especially in birds. Some features of microchromosomes are consistent with that of elusive genes [e.g., high GC, (Bravo et al. 2021)]. However, microchromosomes are conserved in terms of gene order and their genes generally show high protein-level constraints as shown by low Ka/Ks (Waters et al. 2021; Li et al. 2022). Authors need to reconcile their discovery with the previous rich literature.

c) Line (L) 204, among 982 human elusive genes, only 540~390 are shared by other species (e.g., shark). I suggest taking advantage of genome-level synteny based age data generated in Shao et al. 2019 to examine the age distribution of human elusive genes. If a high proportion of them are dated as being old (e.g., shared by jawed vertebrates), the statement that these genes have an ancient origin could be better supported.

References

Bravo GA, Schmitt CJ, Edwards SV. 2021. What have we learned from the first 500 avian genomes. Annu Rev Ecol Evol Syst 52: 611-639.

Li M, Sun C, Xu N, Bian P, Tian X, Wang X, Wang Y, Jia X, Heller R, Wang M et al. 2022. de novo Assembly of 20 Chicken Genomes Reveals the Undetectable Phenomenon for Thousands of Core Genes on Microchromosomes and Subtelomeric Regions. Mol Biol Evol 39.

MacArthur DG, Balasubramanian S, Frankish A, Huang N, Morris J, Walter K, Jostins L, Habegger L, Pickrell JK, Montgomery SB et al. 2012. A systematic survey of loss-of-function variants in human protein-coding genes. Science 335: 823-828.

Nei M. 2013. Mutation-driven evolution. OUP Oxford.

Shao Y, Chen C, Shen H, He BZ, Yu D, Jiang S, Zhao S, Gao Z, Zhu Z, Chen X et al. 2019. GenTree, an integrated resource for analyzing the evolution and function of primate-specific coding genes. Genome Res 29: 682-696.

Waters PD, Patel HR, Ruiz-Herrera A, Alvarez-Gonzalez L, Lister NC, Simakov O, Ezaz T, Kaur P, Frere C, Grutzner F et al. 2021. Microchromosomes are building blocks of bird, reptile, and mammal chromosomes. Proc Natl Acad Sci U S A 118.

Xie KT, Wang G, Thompson AC, Wucherpfennig JI, Reimchen TE, MacColl AD, Schluter D, Bell MA, Vasquez KM, Kingsley DM. 2019. DNA fragility in the parallel evolution of pelvic reduction in stickleback fish. Science 363: 81-84.

Yang H, He BZ, Ma H, Tsaur SC, Ma C, Wu Y, Ting CT, Zhang YE. 2015. Expression profile and gene age jointly shaped the genome-wide distribution of premature termination codons in a Drosophila melanogaster population. Mol Biol Evol 32: 216-228.

Reviewer #3 (Recommendations for the authors):

– In recent years several gene loss pipelines were already published (e.g. TOGA – https://github.com/hillerlab/TOGA) and it would be highly beneficial for this study to compare their gene loss reports with output obtained from existing pipelines (which also address the false discovery rate issue)

– We did not appreciate Figure 5. We strongly recommend finding a more quantitative approach to visualise these results, since heat maps are misleading and show different x-axis and y-axis ranges (zoom-ins/ zoom-outs)

– All Figures and Supplementary Figures showing violin plots need to report the number of genes that underly these distributions.

[Editors' note: further revisions were suggested prior to acceptance, as described below.]

Thank you for resubmitting your work entitled "Gene fate spectrum as a reflection of local genomic properties" for further consideration by eLife. Your revised article has been evaluated by George Perry (Senior Editor) and a Reviewing Editor. A previous reviewer also read and commented on the revised manuscript.

Apparently, the manuscript has been significantly improved but there are some remaining issues that need to be addressed, as outlined below. In view of these comments, we kindly request that you consider revising the manuscript once more. Our hope is that, through this additional revision, the manuscript will be written more clearly and rigorously. Thank you for your understanding and continued efforts in improving your submission.

Reviewer #2 made the following comment. Please consider it during revision.

"The authors attempt to argue that the elusive status is ancient ("Thus, the heterogeneous genomic features driving gene fates toward loss have been in place since the ancestral vertebrates"). However, in response to my previous suggestion regarding chicken microchromosomes, the authors present mixed results. They observed high GC content, high gene density, and short gene length in chicken, similar to the findings in humans (Figure 3). Yet, the critical functional data between the two species are conflicting: human elusive genes exhibit low expression and fewer ATAC-seq peaks, while their chicken counterparts display the opposite pattern. In other words, chicken elusive genes exhibit higher pleiotropy, which may decrease the likelihood of their loss. Thus, these genes are not elusive, and high GC content, high gene density, and short gene length do not necessarily predict elusiveness. Given that the authors only analyzed the functional data of human and chicken genomes, it is not possible to determine whether the "elusive" status is ancient or derived from the human or mammalian lineage. I suggest that the authors analyze the transcriptome data in shark or spotted gar to provide further phylogenetic context. Otherwise, the authors should significantly tone down their statement."

The reviewing editor also has some comments on the title and abstract.

1. Please consider revising the title to "The Impact of Local Genomic Properties on the Evolutionary Fate of Genes Across Vertebrates".

2. Please consider adding a sentence (or something similar) at the end of the abstract. "This study sheds light on the complex interplay between gene function and local genomic properties in shaping gene evolution across vertebrate lineages."

Essential revisions:

1) The identification of "elusive genes" in the current manuscript requires additional scrutinization, given it is the foundation of the whole study. Please check published pipelines in the identification of gene losses (e.g., TOGA – https://github.com/hillerlab/TOGA) and use additional tools such as BLASTX search to test for known technical artifacts when calling genes (homology detection failure, refer to https://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.3000862). Please give a reason for the parameters used in the analysis (e.g., CD-Hit clustering) or examine if the conclusions remain supported by various parameters in the computational pipeline. Also, please take a look at the enrichment of elusive genes in human chr19, and use the synteny-based age estimation of the elusive genes (Shao et al. 2019).

Based on this suggestion, we have repeated the analysis, and to more rigorously exclude possible false positives, and as a result, we have obtained a refined set of 813 elusive genes. The validation of gene loss has employed the ortholog annotations implemented in the RefSeq gene prediction pipeline, a synteny-based ortholog clustering platform other than TOGA, as well as similarity sequence search by MMSeqs2 that performs superior to NCBI BLAST.

For the revised manuscript, we have refined the elusive gene set as the reviewer suggested. In the genome assemblies, we have searched for the orthologs of the elusive genes for the species in which they were missing. The search has been conducted by querying amino acid sequences of the elusive genes with tblastn as well as MMSeqs2 that performed superior to tblastn in sensitivity and computational speed. In addition, regarding a comment by Reviewer 3 we have searched for the orthologs by referring to existing ortholog annotations. We used the ortholog annotations implemented in RefSeq instead of those from the TOGA pipeline: both employ synteny conservation. We have coordinated the identified orthologs with our gene loss criteria–absence from all the species used in a particular taxon–and excluded 268 genes from the original elusive gene set. These genes contain those missing in the previous gene annotations used in the original manuscript but present in the latest ones, as well as those falsely missing due to incorrect inference of gene trees. Finally, the refined set of 813 elusive genes were subject to comparisons with the non-elusive genes. Importantly, these comparisons retained the significantly different trends of the particular genomic, transcriptomic, and epigenomic features between them except for very few cases (Author response table 1). This indicates that both initial and revised sets of the elusive genes reflect the nature of the ‘elusiveness,’ though the initial set contained some noises. We have modified the numbers of elusive genes in the corresponding parts of the manuscript including figures and tables. Additionally, we have added the validation procedures in Methods.

cd-hit clustering with 100% sequence identity only clusters those with identical (and sometimes truncated) sequences, and, in the cluster, the sequences other than the representative are discarded. This means that the sequences remain if they are not identical to the other ones. If the similarity threshold is lowered, both identical and highly similar sequences are clustered with each other, and more sequences are discarded. Therefore, our approach that employs clustering with 100% similarity may minimize false positive gene loss.

Please refer to the reply to [Rev2 point 3] for the abundance of the elusive genes on chromosome 19. To examine this point, we have performed the comparisons between the elusive and non-elusive genes excluding the genes on chromosome 19, and the characteristics remained unchanged even when chromosome 19 was excluded.

2) Please also present the features of other genes (the genomic background other than elusive and non-elusive genes). Are these genes show intermediate patterns between those of elusive and non-elusive genes? Please edit the manuscript accordingly if the definition of non-elusiveness is actually equivalent to the genomic background.

Our aim in this study is to extract the characteristics of the genes that differentiate their fates from retention to loss. To achieve this, we compared the genes with clearly different phylogenetic signatures for gene loss, namely elusive and non-elusive genes.

The remainders excluding the elusive and non-elusive genes do not necessarily exhibit intermediate features. Because our definitions of the elusive and non-elusive genes were stringent, the remainders may contain considerable numbers of genes with loss in more restricted taxa than our criterion. In addition, the reminders contain those with other particular phylogenetic signals (e.g., frequently duplicated). These do not necessarily exhibit intermediate features, and at least the former is rather closer to elusive genes.

3) It would be informative to test the links between recombination rate / LD and the genomic locations of elusive genes (compared against randomly sampled genes).

We have retrieved fine-scale recombination rate data of males and females from https://www.decode.com/addendum/ (Suppl. Data of Kong, A et al., Nature, 467:1099–1103, 2010) and have compared them between the gene regions of the elusive and non-elusive genes. Both comparisons show no significant differences: average 0.829 and 0.900 recombinations/kb for the elusive and non-elusive genes, respectively, p=0.898, for males; average 0.836 and 0.846 recombinations/kb for the elusive and non-elusive genes, respectively, p=0.256, for females.

4) Please control confounding factors such as gene expression level and confirm whether the proxy of mutation rate (i.e., Ks) is actually confounded by gene importance.

We thank the reviewer for this important comment. We totally agree that transcriptomic and epigenomic features cannot be easily distinguished from gene dispensability and do not think that these features of the elusive genes can be explained solely by intrinsic properties of the genomes. Our motivation for investigating the expression profiles of the elusive gene is to understand how they lost their functional indispensability (original manuscript L285-286 in Results). We also discussed the possibility that sequence composition and genomic location of elusive genes may be associated with epigenetic features for expression depression, which may result in a decrease of functional constraints (original manuscript L470-474 in Discussion). Nevertheless, we think that the original manuscript may have contained misleading wordings, and thus we have edited them to better convey our view that gene expression and epigenomic features are related to gene function.

(P.2, Introduction) “This evolutionary fate of a gene can also be affected by factors independent of gene dispensability, including the mutability of genomic positions, but such features have not been examined well.”

(P6, Introduction) “These data assisted us to understand how intrinsic genomic features may affect gene fate, leading to gene loss by decreasing the expression level and eventually relaxing the functional importance of ʻelusiveʼ genes.”

(P33, Discussion) “Another factor is the spatiotemporal suppression of gene expression via epigenetic constraints. Previous studies showed that lowly expressed genes reduce their functional dispensability (Cherry, 2010; Gout et al., 2010), and so do the elusive genes.”

Additionally, responding to the advices from Reviewers 1 and 2, we have added a new section Elusive gene orthologs in the chicken microchromosomes in which we describe the relationship between the elusive genes and chicken microchromosomes. In this section, we also argue for the relationship between the genomic feature of the elusive genes and their transcriptomic and epigenomic characteristics. In the chicken genome, elusive genes did not show reduced pleiotropy of gene expression nor the epigenetic features relevant with the reduction, consistently with the moderation of nucleotide substitution rates. This also suggests that the relaxation of the ‘elusiveness’ is associated with the increase of functional indispensability.

(P27, Elusive gene orthologs in the chicken microchromosomes in Results) “Our analyses indicates that the genomic features of the elusive genes such as high GC and high nucleotide substitutions do not always correlate with a reduction in pleiotropy of gene expression that potentially leads to an increase in functional dispensability, although these features have been well conserved across vertebrates. In addition, the avian orthologs of the elusive genes did not show higher K_A and K_S values than those of the non-elusive genes (Figure 3; Figure 3—figure supplement 1), likely consistent with similar expression levels between them (Figure 5—figure supplement 1) (Cherry, 2010; Zhang and Yang, 2015). With respect to the chicken genome, the sequence features of the elusive genes themselves might have been relaxed during evolution.”

Also, please refer to [Rev1- point2] for the mutability of the elusive genes. To examine this point, we have computed nucleotide sequence differences in introns, namely K_I, between the human and chimpanzee genomes. This analysis revealed higher K_I values in the elusive genes than in the non-elusive genes, which is in line with our original hypothesis. The results have been added in the revised manuscript.

5) Please consider extending the analyses on fish and birds to other genomic features.

Please refer to the reply to [Rev2 point 4a] for the comparison between the elusive and non-elusive gene orthologs of nonmammalian vertebrates. We analyzed expression profiles and ATAC-seq peak densities of the elusive and non-elusive gene of these animals. We have created a new section Elusive gene orthologs in the chicken microchromosomes in Results and described the results in this section and Discussion.

We analyzed expression profiles and ATAC-seq peak densities of the elusive and non-elusive gene of these animals. We have created a new section Elusive gene orthologs in the chicken microchromosomes in Results and described the results in this section and Discussion.

We appreciate the reviewer for this meaningful suggestion. As a response, we have computed the differences in intron sequences between the human and chimpanzee genomes and compared them between the elusive and non-elusive genes. As expected, we found larger sequence differences in introns for the elusive genes than for the non-elusive genes. In Figure 2c of the revised manuscript, we have included the distribution of K_I, sequence differences in introns between the human and chimpanzee genomes for the elusive and non-elusive genes. Additionally, we have added the corresponding texts to Results and the procedure to Methods as shown below.

(P11, Identification of human ‘elusive’ genes in Results) “In addition, we computed nucleotide substitution rates for introns (K_I) between human and chimpanzee (Pan troglodytes) orthologs and compared them between the elusive and non-elusive genes.”

(P11, Identification of human ‘elusive’ genes in Results) “Our analysis further illuminated larger K_S and K_I values for the elusive genes than in the non-elusive genes (Figure 2b, c; Figure 2—figure supplement 1). Importantly, the higher rate of synonymous and intronic nucleotide substitutions, which may not affect changes in amino acid residues, indicates that the elusive genes are also susceptible to genomic characteristics independent of selective constraints on gene functions.”

(P39, Methods) “To compute nucleotide sequence differences of the individual introns, we extracted 473 elusive and 4,626 non-elusive genes that harbored introns aligned with the chimpanzee genome assembly. The nucleotide differences were calculated via the whole genome alignments of hg38 and panTro6 retrieved from the UCSC genome browser.”

6) The authors should reconcile the findings in this study with previous reports about microchromosomes.

Please refer to the reply to [Rev2-point 4b] for the reconciliation of our findings of the elusive genes with the features of microchromosomes. In the human genome, we found a significant overlap between the elusive genes and the genes whose chicken orthologs are located on microchromosomes. Although the chicken microchromosomes shared some sequence features such as high GC content and high K_S values in common with the elusive genes in our sense but exhibited opposite transcriptomic and epigenetic trends. Although the result does not change the basis of our study in the human genome, it indicates that, in the course of evolution, genomic features of the elusive genes are not always associated with a reduction of pleiotropy of gene expression. The results have been described in the newly created section Elusive gene orthologs in the chicken microchromosomes in Results.

7) Please consider improving the clarity/presentation of Figures 5 and 7 and examine whether the pattern remains robust using various parameter sets.

Please refer to the replies to [Rev3- point 2] for this point. First, following this suggestion, we have conducted a statistical test to see whether the elusive genes contain more genes with restricted expression profiles (H’ < 1) than the non-elusive genes, and this trend was statistically supported. We have added the gene numbers of the individual categories and the result of this statistical tests to Figure 5. Therefore, we think the classification of the elusive genes based on the threshold (H’ = 1) is reasonable in Figure 7.

In addition, we have conducted statistical tests in a similar way with different thresholds for H’ (2, 3, and 0.5) and found that the pattern remains robust (p = 3.85x10^-67, 1.29x10^-50, and 9.40x10^-46 setting thresholds H’ at 2, 3, and 0.5, respectively, for the GTEx dataset and p = 8.35x10^-61, 2.01x10⁶¹, and 1.25x10^-24 setting thresholds H’ at 2, 3, and 0.5, respectively, for the Descartes dataset).

To use Figure 7 in a new section in Results, we have added an ideogram showing the distribution of the genes that retain the chicken orthologs in microchromosomes. In response to the comment by Reviewer 2 [Rev2- point 4b], we have performed statistical tests and found that the elusive genes were significantly more abundant in orthologs in microchromosomes than the non-elusive genes. Furthermore, the observation that the elusive genes prefer to be located in gene-rich regions was already statistically supported (Figure 2f).

As shown in Figure 5, Shannon’s H' ranged from zero to approximately 4 (exact maximum value is 3.97) and 5 (5.11) for the GTEx and Descartes gene expression datasets, respectively. Although the threshold H'=1 was an arbitrarily set, we think that it is reasonable to classify the genes with high pleiotropy from those with low pleiotropy.

8) Please think of a better term than "elusive gene" to describe the genes that were lost independently in different lineages. Please also clearly define other terms in the manuscript (e.g., functionally indispensable vs. importance, are they the same concept?)

The phrase ‘elusive gene’ was already used in our previous paper to follow the advice of peer-reviewers for our past manuscripts, and for consistency, we would like use this term, while we have modified the sentence introducing this term to more carefully explain it.

9) Please consider presenting the current study in the framework of mutation-selection balance, and better explain the novelty of the study over previous tremendous studies about gene losses.

Please refer to the replies below to convey our novelty. We have cited existing literature in the revised manuscript.

Reviewer #2 (Recommendations for the authors): Overall, I believe that this is an interesting study. However, this version of the manuscript could be significantly improved in terms of logical depth and methodological stringency.

1. Authors actually support the concept of mutation-driven evolution, i.e., the high mutation rate in genomic regions harboring the elusive genes would predispose their fate toward death. To increase the significance of their work, I suggest authors cite (Nei 2013; Xie et al. 2019) and put their work in a bigger context.

We appreciate the reviewer for considering a way to enhance the significance of our study. We now recognize the importance of the study by Xie et al. (2019) reporting the recurrent loss of an enhancer that modulates fin morphogenesis in stickleback. As suggested, in the revised manuscript, we have cited this paper in Introduction.

(P4, Introduction) “In the stickleback genome, a Pitx1 enhancer was independently lost in multiple lineages inhabiting freshwater due to its genomic location in a structurally fragile site, leading to recurrent loss of pelvic fins (Xie et al., 2019). Genes and genomic elements in such particular regions may be prone to loss in a more neutral manner than the relaxation of functional importance or via functional adaptations.”

Additionally, to enhance the broad interests of our study, we have cited the Dr. Nei’s book in Discussion as shown below.

(P31, Discussion) “Close coordination of the studies on gene evolution with germline mutation repertoires and spectra, which can be approximated from the collection of de novo mutations obtained by trio sequencing, may further facilitate our understanding of gene fates driven by heterogeneous genomic features—this would be viewed as ‘mutation-driven’ evolution (Nei, 2013).”

2. Authors mentioned that elusive genes are less important and thus more prone to loss. In my view, pleiotropy is a better term compared to importance. That is, elusive genes are less pleiotropic [e.g. narrowly expressed, Figure 5] and thus their loss are more tolerable or easily compensated by other genes. Actually, narrow expression breadth has been observed to be correlated with gene loss in both humans and flies (MacArthur et al. 2012; Yang et al. 2015).

We thank the reviewer for the thoughtful suggestion. We agree that the word pleiotropy is more suitable in our manuscript. We have modified the manuscript as shown below.

(P20, Transcriptomic natures of elusive genes in Results) “Our findings demonstrate that some elusive genes harbor low-level and spatially-restricted expression profiles, i.e., less pleiotropic states, which are rarely observed in the non-elusive genes.”

(P23, Epigenetic nature of elusive genes in Results) “… we classified the elusive genes into two groups based on the pleiotropy in terms of gene expressions: that is, 589 elusive genes with Shannon’s diversity index H¢ ≥ 1 were ubiquitously expressed, i.e, more pleiotropic, and 150 of those with H¢ < 1 were expressed in only a few or none of the tissues examined, i.e., less pleiotropic (Figure 5).”

(P34, Discussion) “Elusive genes with reduced pleiotropy may have limited opportunities to function, potentially leading to loss of their important roles.”

3. I am generally convinced that authors reliably identified elusive genes by identifying gene loss events in the common ancestor of multiple descendant species (to control for errors induced by assembly or annotation, Figure 1). However, Figure 7 shows the enrichment of elusive genes in human chr19. This chromosome is well known to be enriched with tandemly duplicated Krueppel-associated box C2H2 zinc-finger protein family (KZNF), many of which are primate-specific (Shao et al. 2019). I suspect that tree-based strategy implemented in Figure 1 could not be able to dissect the evolution of this super complex gene family. I am proposing two specific analyses: how many elusive genes encoded by chr19 are KZNFs? how many of them have Ensembl one-to-one orthologs across mammals?

Responding to this comment, we investigated the KZNF genes on chromosome 19. We identified 206 those genes in chromosome 19, of which 75 were found to be elusive. Of these, 30 genes retained one-to-one orthologs of the mouse or dog. Although we excluded the nearly identical paralogs in our pipeline in the original manuscript and have refined the elusive gene set with synteny-based ortholog annotation in the current manuscript, the elusive KZNF genes have remained in this refined gene set. The elusive KZNF genes on chromosome 19 were lost in some mammalian lineages and duplicated in early primates.

Motivated by this comment, we conceived a possibility whether the enrichment of the paralogs of the elusive genes on chromosome 19 overrepresents the features relevant to the ‘elusiveness’. We thus have created another set of elusive genes, those excluding the genes on chromosome 19, and performed comparisons on the genomic, transcriptomic, and epigenomic features of the elusive and non-elusive genes. The results showed that the significant/non-significant differences have been maintained except for a few cases (see Author response table 1), indicating that the characteristics remained unchanged even when chromosome 19 is excluded.

4. With the patterns in Figure 3 and 7, authors argued that features of elusive genes are deeply ancient and could be inherited from the microchromosomes of early vertebrates. This statement has multiple problems.

a) Figure 3 only show genomic level features (e.g., high GC content) conserved in multiple vertebrates including shark and chicken. Epigenetic features analyzed in Figure 5 to 6 were only based on human data. I suggest authors to extend these analyses to shark or chicken. Although some epigenetic data could not be available for these species, transcriptome data analyzed in Figure 5 should be available for at least some species.

In response to this reviewer’s comment, we have investigated the transcriptomic and epigenetic characteristics of the orthologs of the elusive genes in non-mammalian vertebrates as we had done for human in the original manuscript. We have retrieved gene expression profiles of normal tissues and early embryos from bulk RNA-seq data of chicken, tropical clawed frog, coelacanth, and spotted gar, respectively, from the Bgee database (https://bgee.org/). We have then compared expression profiles between the orthologs of the elusive and non-elusive genes for the individual species using the same procedure as those in the manuscript. Only the anole, coelacanth, and gar orthologs of the elusive genes show the enrichment of low H’ values (newly created figure in Figure 5—figure supplement 1). This indicates that the low pleiotropy of expression of the elusive genes is not always observed in the non-mammalian species.

Furthermore, we have compared epigenomic properties between the orthologs of the elusive and non-elusive genes in the chicken genome. We have retrieved ATAC-seq narrow peak data for tissues of chicken embryos from NCBI GEO and compared the density of ATAC peaks between the orthologs of the elusive and non-elusive genes. The result indicates that, in five samples out of the eight, the orthologs of the elusive genes retained more ATAC peaks than those of the non-elusive genes, and that the reminders did not show this difference (newly created figure in Figure 6—figure supplement 5). This observation may remind us of a link between the reduction of the ‘elusiveness’ and the decrease of functional dispensability in gene evolution. However, it should be interpreted carefully, as different sets of tissues were used for the transcriptomic and epigenomic analyses between human and chicken. As described above, the chicken ATAC-seq experiments were mainly performed with developing embryos, while the human ATAC-seq used in our study were performed with cell lines. Nevertheless, the cross-species comparison of the epigenetic features suggests that the sequence features relevant to the elusive genes do not always induce the epigenetic conditions for gene expression depletion.

We have newly created Figure 5—figure supplement 1 and Figure 6—figure supplement 5 as shown below, described the results in the newly created section Elusive gene orthologs in the chicken microchromosomes in Results, and discussed them in Discussion. Also, we described the procedures in Methods.

(P-26-27, Elusive gene orthologs in the chicken microchromosomes in Results)

“Elusive gene orthologs in the chicken microchromosomes

The heterogeneous locations of the elusive genes can also be examined from a chromosome-scale viewpoint (Figure 7; Figure 7—figure supplement 1). The visualization via chromosome ideograms indicated an overlap of the elusive genes with the genomic regions enriched for the genes whose chicken orthologs are on the microchromosomes (chromosomes 11-38 and W), providing a statistical support for this trend (p=0.0175; Figure 7a). Indeed, microchromosomes of the chicken and other vertebrate exhibit genomic features including high GC-content, high gene density, and rapid nucleotide substitutions in comparison with their macrochromosomes (Groenen et al., 2009; International Chicken Genome Sequencing Consortium, 2004; Schield et al., 2019; Waters et al., 2021), which also characterize genomic regions containing elusive genes. On the contrary, previous studies revealed that the chicken microchromosomes are preferentially located in the A compartments of the nucleus (Perry et al., 2020) and are early replicating (McQueen et al., 1998). These characteristics associated with the microchromosomes were opposite characteristics to the human genomic regions preferentially containing the elusive genes.

We further analyzed the ATAC-seq peaks in the chicken genome and found more peaks in the genomic regions including the elusive gene orthologs than in those containing non-elusive gene orthologs in four samples out of eight and no significant differences in the peak density in the four remaining samples (Figure 6; Figure 6—figure supplement 5). These observations indicate that, in an epigenetic manner, the chicken orthologs of the elusive genes are not regulated to reduce their expression level. This idea was further supported by a comparison of the expression profiles between the chicken orthologs of the elusive and non-elusive genes, showing no significant differences between them (Figure 5; Figure 5—figure supplement 1). Our analyses indicate that the genomic features of the elusive genes such as high GC and high nucleotide substitutions do not always correlate with a reduction in pleiotropy of gene expression that potentially leads to an increase in functional dispensability in the course of vertebrate evolution. In addition, avian orthologs of the elusive genes did not show higher K_A and K_S values than those of the non-elusive genes (Figure 3; Figure 3—figure supplement 1), likely consistent with not significant difference in gene expression levels between them in the species (Figure 5—figure supplement 1) (Cherry, 2010; Zhang and Yang, 2015). With respect to the chicken genome, the sequence features of the elusive genes might have been relaxed during evolution.”

(P32-33, Discussion) “A chromosomal-scale view of the distribution of elusive genes illuminated their significant correlation with the genes whose chicken orthologs are located on microchromosomes (Figure 7a). More importantly, genomic regions rich in elusive genes were traced back to the microchromosomes of the ancestral gnathostomes by reconstructing chromosomes of the ancestral genomes (Figure 7b). This inference of ancestral karyotypes augments our observations that some elusive natures of genomic sequences have been retained for hundreds of millions of years (Figure 3). In other words, the result suggests that the disparity of genomic regions which allows the ‘elusiveness’ for the genes has been retained during vertebrate evolution. On the other hand, comparisons of the expression profiles between the orthologs of the elusive and non-elusive genes for non-mammalian vertebrates indicate that the elusive genes are not always associated with the restricted expression profiles (Figure 5; Figure 5—figure supplement 1). Additionally, in the chicken genome, this trend in gene expression may be correlated with the abundance of ATAC-seq peaks in the elusive genes (Figure 6—figure supplement 5). These findings again suggest that the elusive genes are not always associated with a reduction in pleiotropy of gene expression, which may lead to an increase of functional dispensability during evolution. It should be noted that the choices of tissues used in these analyses were largely different between the human and non-mammalian vertebrates (Tables S3, S4). The chicken ATAC-seq data could be obtained only from developing embryos, while the human ATAC-seq in ENCODE were performed with cell lines. Therefore, the aforementioned interpretation should be treated carefully.”

(P40, Methods) “We also compared expression profiles and ATAC-seq peak densities between the orthologs of the elusive and non-elusive genes in nonmammalian vertebrates in a similar way as we did with the human datasets.

Normalized gene expression profiles from RNA-seq data of normal adult tissues and early embryos for chicken, green anole, Western clawed frog, coelacanth, and spotted gar were obtained from the Bgee version 15 database (Bastian et al., 2021) (Table S4). ATAC-seq narrow peak signals of chicken tissues were retrieved from NCBI GEO (Table S4) followed by coordination of the genome assembly with galGal5 with the UCSC liftOver tool (Hinrichs et al., 2006) as needed.”

b) In Figure 7, authors propose the concept related with microchromosomes. These chromosomes have been extensively studied, especially in birds. Some features of microchromosomes are consistent with that of elusive genes [e.g., high GC, (Bravo et al. 2021)]. However, microchromosomes are conserved in terms of gene order and their genes generally show high protein-level constraints as shown by low Ka/Ks (Waters et al. 2021; Li et al. 2022). Authors need to reconcile their discovery with the previous rich literature.

Responding to this reviewer’s comment, we have first examined whether the chicken orthologs of the elusive genes tend to be located on microchromosomes. We have classified the chicken orthologs of the elusive and non-elusive genes into those on macro- and microchrosomes. Fisher’s exact test has indicated that elusive genes are enriched in the chicken microchromomes (p=0.0175, odds ratio=1.46; Author response table 2). In addition, the elusive genes are enriched in the genomic regions corresponding to the ancestral microchromosomes in early vertebrates (p=9.50x10^-24, odds ratio=2.31; Author response table 3).

As the reviewer noted, the genomic regions including the elusive genes share several genomic characteristics with microchromosomes: high GC content, high gene density, and high K_A and K_S values (although K_A/K_S values are lower for newly identified genes in the chicken microchromosomes in Li et al., 2022). On the contrary, previous analyses showed that the chicken microchromosomes tend to be early replicating and located in the A compartment in the nucleus, where the genes are actively transcribed. These observations are concordant with our finding that the chicken orthologs of the elusive genes are rich in ATAC peaks (see in the previous reply). We could not determine which epigenomic states are ancestral, but the genomic features of high GC-content and high nucleotide substitutions are not always correlated to reduced pleiotropy of gene expression, potentially leading to an increase in functional dispensability.

In the revised manuscript, we have created a new section Elusive gene orthologs in the chicken microchromosomes in Results to describe the enrichment of the elusive genes in the chicken microchromsomes and the association of their genomic characters with those of the chicken orthologs, while we cite the literature mentioned by the reviewer. In addition, we have mentioned the relevance of the elusive genes to the ancestral macrochromosomes in early vertebrates in Discussion. The corresponding sentences in the revised manuscript have already been included above in the previous.

c) Line (L) 204, among 982 human elusive genes, only 540~390 are shared by other species (e.g., shark). I suggest taking advantage of genome-level synteny based age data generated in Shao et al. 2019 to examine the age distribution of human elusive genes. If a high proportion of them are dated as being old (e.g., shared by jawed vertebrates), the statement that these genes have an ancient origin could be better supported.

We thank the reviewer for an opportunity to revisit the age of the elusive genes. We have examined this for the revised manuscript. Though the reviewer suggested using GenTree (http://gentree.ioz.ac.cn/) for this purpose, this database consists of euteleostome animals and does not include chondrichthyans and more distantly related animals. Instead, we have used the ortholog groups of jawed vertebrates built in this study and further added older branching dates of their orthologs and paralogs from the Ensembl Gene Tree. Of the 813 elusive genes in the revised dataset, 152 retain only mammalian orthologs, and this proportion is higher than that of the non-elusive genes (65 out of 8,050, p=2.50x10^-110). We have then extracted the 517 elusive and 7,900 non-elusive genes whose ancestors were dated at the early evolution of jawed vertebrates or older. This subset of old gene age allowed us to examine if the elusive genes also retain loss-prone nature in non-mammalian vertebrates. On average, 40% of the 517 elusive genes are found to be retained in the nonmammalian vertebrates. On the other hand, more than 90% of the 7,900 non-elusive genes are retained by these species. These results indicate that while relatively young genes are frequently found in the elusive genes, the loss-prone nature of the elusive genes in the non-mammalian vertebrates is recapitulated by using this gene set of old origin. We have replaced Figure 3—figure supplement 1 with a revised version.

(P14, Tracing elusiveness back along the vertebrate evolutionary tree in Results) “We found that 152 out of 813 elusive genes originated in mammalian lineages, and this proportion was larger than those of the elusive genes (65 out of 8050, p = 2.50 ´ 10^-110), indicating that the elusive genes are more abundant in recently born genes than non-elusive genes. We then selected 517 elusive and 7,900 non-elusive genes that originated in the common ancestors of jawed vertebrates or earlier. These subsets allowed us to examine the degree of retention of non-mammalian vertebrate orthologs in the elusive and non-elusive genes. On average, approximately 40% of these elusive genes were found to be retained by non-mammalian vertebrates, while this proportion increased up to 90% for the non-elusive genes. (Figure 3—figure supplement 1a). In the coelacanth, gar, and shark, the orthologs of the elusive genes were less frequently retained by all the species than those of the non-elusive ones (Figure 3—figure supplement 1b). The results suggest that the origins of the loss-prone propensity of the elusive genes potentially date back to the period long before the emergence of the Mammalia.”

Reviewer #3 (Recommendations for the authors):

– In recent years several gene loss pipelines were already published (e.g. TOGA – https://github.com/hillerlab/TOGA) and it would be highly beneficial for this study to compare their gene loss reports with output obtained from existing pipelines (which also address the false discovery rate issue)

Please refer to the reply to Essential revisions 1. We understand that TOGA is an ortholog search pipeline, whose output include gene loss information, employing synteny information. TOGA inputs genome alignments between the reference and all the target species, which requires tremendous computational time in our case involving a number of species. Instead, we have incorporated ortholog annotations implemented in RefSeq, another synteny-based ortholog detection platform, into our validation for gene loss.

– We did not appreciate Figure 5. We strongly recommend finding a more quantitative approach to visualise these results, since heat maps are misleading and show different x-axis and y-axis ranges (zoom-ins/ zoom-outs)

For the visibility of Figure 5, we have added the percentage of the genes of H'≥1 and H'<1 in the plots and performed the statistical test to see the fraction of the genes of H'≥1 and H'<1 is equal between the elusive and non-elusive genes. Additionally, we have coordinated x-axis and y-axis ranges between the elusive and non-elusive genes. We have replaced Figure 5 with a revised version containing these modified figure panels.

– All Figures and Supplementary Figures showing violin plots need to report the number of genes that underly these distributions.

We have added the sample numbers of the plots to all figures or figure legends.

[Editors' note: further revisions were suggested prior to acceptance, as described below.]

Apparently, the manuscript has been significantly improved but there are some remaining issues that need to be addressed, as outlined below. In view of these comments, we kindly request that you consider revising the manuscript once more. Our hope is that, through this additional revision, the manuscript will be written more clearly and rigorously. Thank you for your understanding and continued efforts in improving your submission.

Reviewer #2 made the following comment. Please consider it during revision.

"The authors attempt to argue that the elusive status is ancient ("Thus, the heterogeneous genomic features driving gene fates toward loss have been in place since the ancestral vertebrates"). However, in response to my previous suggestion regarding chicken microchromosomes, the authors present mixed results. They observed high GC content, high gene density, and short gene length in chicken, similar to the findings in humans (Figure 3). Yet, the critical functional data between the two species are conflicting: human elusive genes exhibit low expression and fewer ATAC-seq peaks, while their chicken counterparts display the opposite pattern. In other words, chicken elusive genes exhibit higher pleiotropy, which may decrease the likelihood of their loss. Thus, these genes are not elusive, and high GC content, high gene density, and short gene length do not necessarily predict elusiveness. Given that the authors only analyzed the functional data of human and chicken genomes, it is not possible to determine whether the "elusive" status is ancient or derived from the human or mammalian lineage. I suggest that the authors analyze the transcriptome data in shark or spotted gar to provide further phylogenetic context. Otherwise, the authors should significantly tone down their statement."

We are aware of your concern about the inconsistency of the results between human and chicken, and therefore already included comparisons of the gene expression profiles between the orthologs of the elusive and non-elusive genes with several non-mammalian vertebrates in the previous revision (Figure 5—figure supplement 1). This comparison indicated that the orthologs of the elusive genes were rich in those with lower pleiotropy in anole, coelacanth, and gar, while this tendency was not found in chicken and frog. The result suggests that the orthologs of the elusive genes are likely to have retained the “elusive” features in the common ancestor of bony vertebrates.

In addition, regarding the reviewer’s comment, we should have clearly stated that in non-mammalians, the genes we focused on are not necessarily ‘elusive genes’ but rather ‘the orthologs of the human elusive genes’ in the previous manuscript. These orthologs do not always retain the ‘elusiveness’ as the human genes do. Interestingly, the orthologs of the elusive genes in the chicken genome do not exhibit significant differences in K_A and K_S values from those of the non-elusive genes, which is potentially associated with increased pleiotropy of gene expression in these genes. These observations suggest that the orthologs of the human elusive genes we identified have increased functional importance in the lineage leading to chicken.

We have revised the following parts of Results and Discussion to describe transcriptomic characteristics of the orthologs of the human elusive genes in non-mammalians more clearly in the current manuscript as included below.

Elusive gene orthologs in the chicken microchromosomes in Results

“We further compared expression profiles between the orthologs of the human elusive and non-elusive genes in several non-mammalian vertebrates and found that the orthologs of the elusive genes tend to exhibit low pleiotropy in green anole, coelacanth, and gar but not in Western clawed frog. The result suggests that the low pleiotropy of the elusive genes has persisted at least since the bony vertebrate ancestors. With respect to the chicken genome, the ‘elusive’ features the genes orthologous to human elusive genes might have been relaxed —functional importance of the orthologs have increased—during evolution leading to chicken.”

Discussion

“On the other hand, comparisons of the expression profiles between the orthologs of the elusive and non-elusive genes for non-mammalian vertebrates suggest that the orthologs of the elusive genes have been associated with a reduction in pleiotropy of gene expression since vertebrate ancestors but acquired the diverse expressions in chicken and frog (Figure 5; Figure 5—figure supplement 1). Additionally, in the chicken genome, the diverse expressions of the chicken orthologs of the human elusive genes may be correlated with the abundance of ATAC-seq peaks (Figure 6—figure supplement 5). These findings again suggest that the chicken orthologs of the human elusive genes have increased pleiotropy of gene expression, which may lead to a lineage-specific acquisition of functional indispensability.”

The reviewing editor also has some comments on the title and abstract.

1. Please consider revising the title to "The Impact of Local Genomic Properties on the Evolutionary Fate of Genes Across Vertebrates".

Thank you for your suggestion. We would like to use the title that you suggested without the last two words (‘The Impact of Local Genomic Properties on the Evolutionary Fate of Genes’). This title is more likely to increase the accessibility of our paper to a broad readership beyond those who are primarily investigating vertebrates, potentially prompting them to investigate genomic properties associated with the fate of genes in diverse organisms such as invertebrates, fungi, and plants.

2. Please consider adding a sentence (or something similar) at the end of the abstract. "This study sheds light on the complex interplay between gene function and local genomic properties in shaping gene evolution across vertebrate lineages."

Thank you for your suggestion that strengthens our focus and conclusion. We have revised the abstract as follows.

Thus, the heterogeneous genomic features driving gene fates toward loss have been in place and may sometimes have relaxed the functional indispensability of such genes. This study sheds light on the complex interplay between gene function and local genomic properties in shaping gene evolution that has persisted since the vertebrate ancestor.

Author details

1. Yuichiro Hara

   Research Center for Genome & Medical Sciences, Tokyo Metropolitan Institute of Medical Science, Tokyo, Japan

   Contribution

   Conceptualization, Data curation, Funding acquisition, Writing – original draft, Writing – review and editing, Investigation, Methodology

   For correspondence

   hara-yi@igakuken.or.jp

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-7817-8963

2. Shigehiro Kuraku

   1. Molecular Life History Laboratory, Department of Genomics and Evolutionary Biology, National Institute of Genetics, Mishima, Japan
   2. Department of Genetics, Sokendai (Graduate University for Advanced Studies), Mishima, Japan
   3. RIKEN Center for Biosystems Dynamics Research, Kobe, Japan

   Contribution

   Conceptualization, Supervision, Writing – review and editing, Funding acquisition

   Competing interests

   Reviewing editor, eLife

   "This ORCID iD identifies the author of this article:" 0000-0003-1464-8388

• 1,231

  Page views

• 141

  Downloads

• 0

  Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.