(a) Pipeline of ortholog group clustering and gene loss detection. (b) Definition of an elusive gene schematized with ortholog presence/absence pattern referring to a taxonomic hierarchy. Red and …
Distributions of non-synonymous, synonymous, and intronic nucleotide substitution rates, namely KA (a), KS (b), and KI (c) values, respectively, between the human–chimpanzee orthologs of the elusive …
Distributions of KA and KS values between the orthologs of human elusive and non-elusive genes of closely related vertebrates. Correction for multiple testing was performed for comparison in each …
Retention of the genomic and evolutionary characteristics of the human elusive genes across vertebrates. The individual round squares with arrows indicate significant increases or decreases of the …
(a) Number of retained orthologs of the human elusive and non-elusive genes that originated in the common ancestors of the gnathostomes or older age. (b) Intersections of the retained orthologs …
A 2×2 contingency table in Figure 3—figure supplement 1.
Distribution of (a) gene length and (b) GC content of the orthologs of the human elusive and non-elusive genes and (c) distribution of the gene density of the genomic regions where the orthologs of …
Comparison of the density of rare single-nucleotide variants (SNVs) (a), deletion copy number variants (CNVs) (b), duplication CNVs (c), and Z-scores of synonymous (d), missense (e), and …
The figure shows density scatter plots of the expression quantity and divergence of elusive and non-elusive genes. The numbers of the elusive/non-elusive genes and those for which the expression …
Density scatter plots of the expression quantity and divergence of elusive and non-elusive genes. The total numbers of the elusive/non-elusive genes and the number of them for which expression data …
Comparison of the distribution of ATAC-seq peak density (a), length of the topologically associating domains (TADs) including the elusive or non-elusive genes (b), the replication timing indicator …
Comparison of the distribution of ATAC-seq peak density between the elusive and non-elusive genes across multiple cell types. In the elusive gene panels (orange bar), purple bar indicates the …
Comparison of the distribution of length of TADs including the elusive or non-elusive genes across multiple cell types. In the elusive gene panels (orange bar), purple bar indicates the elusive …
Comparison of the distribution of replication timing indicator based on Repli-seq between the elusive and non-elusive genes across multiple cell types. In the elusive gene panels (orange bar), …
Comparison of frequency of overlap with LADs computed from Lamin B1 ChIP-seq data between the elusive and non-elusive genes across multiple data. In the elusive gene panels (orange bar), purple bar …
Comparison of the distribution of ATAC-seq peak density between the orthologs of elusive (210) and non-elusive (7218) genes in the chicken genome. Correction for multiple testing was performed for …
Red and dark blue horizontal bars beside the chromosome ideogram represent the location of elusive genes with restricted expression (Shannon’s H′ < 1) and the other elusive genes, respectively. (a) …
Features | Non-significant in theinitial gene set(1,081 elusive genes) | Non-significant in thecurrent gene set(813) | Non-significant in the current gene setexcluding chr19(669) |
---|---|---|---|
Gene density in the turkey genome | ✓ | ✓ | ✓ |
Gene density in the green anole genome | ✓ | ||
Gene density in the bamboo shark genome | ✓ | ||
Gene density in the whale shark genome | ✓ | ✓ | ✓ |
KS in avians | ✓ | ✓ | |
KA in avians | ✓ | ✓ | ✓ |
KA in sharks | ✓ | ✓ | ✓ |
ATAC-seq peak density for GM23338 | ✓ | ||
Lamin B1 ChIP-seq peak density for K562 | ✓ |
The other features showed significantly different trends between the elusive and non-elusive genes for all of the elusive gene sets and thus are not included in this table.
Elusive | Non-elusive | |
---|---|---|
Macrochromosome (chr1-10, Z) | 93 | 4211 |
Microchromosome (chr11-, W) | 68 | 2078 |
Genes in non-chromosome scaffolds and mitochondrial genome were excluded Fisher's exact test p=0.0175, odds ratio=1.48
Elusive | Non-elusive | |
---|---|---|
Ancestral macrochromosome | 395 | 5950 |
Ancestral microchromosome | 296 | 1929 |
Genes in the genomic regions that did not correspond to the ancestral macro/microchromsomes were excluded. Fisher's exact test p=9.50x10-24, odds ratio=2.31.
Supplementary Tables S1, S3, S4.
(a) Supplementary Table S1. Vertebrate species used for creating gene phylogenies. (b) Supplementary Table S3. ENCODE accession ID list used for epigenomic analyses. (c) Supplementary Table S4. RNA-seq and ATAC-seq samples of non-mammalian vertebrates.
Supplementary Table S2.
Characteristics of the elusive and non-elusive genes in the human genome.