Kidney ‘pathway-orphan’ genes as a potential overlooked source of novel insights into kidney health and disease

Dianne Acoba, Sophie Gunnarsson, Anna Reznichenko
doi: https://doi.org/10.1101/2024.07.09.24309717

ABSTRACT

Hundreds of protein-coding genes still have unknown functions and no biological pathway annotations. Mechanistic studies typically investigate well-known genes, despite growing evidence of associations between disease and some understudied genes. We hypothesized that examining these overlooked genes lacking pathway annotation could lead to new insights on chronic kidney disease (CKD) pathogenesis. Human protein-coding genes (N=19,258) from HGNC were annotated with pathway membership using a comprehensive collection of pathway databases (including but not limited to GO, KEGG, MSigDB, Reactome, WikiPathways) to reveal ‘pathway-orphan genes’— genes that are not members of any biological pathways. Expression and enrichment of pathway-orphan genes in healthy kidneys were established using GTEx data. Kidney-expressed pathway-orphan genes were tested for differential case-control expression using publicly available CKD datasets with kidney tissue RNA-seq transcriptomics profiling (GSE98422, GSE142025, GSE175759, GSE197307, Levin et al., 2020). A total of 286 genes in the human genome currently lack any biological pathway membership and are identified as pathway-orphan genes. We have determined 97 of these pathway-orphan genes are expressed in healthy kidneys, with several showing kidney-specific enrichment. Furthermore, 34 pathway-orphan genes show significant modulation of expression (FDR < 0.05) in CKD kidney, out of which 10 robustly demonstrate concordant directionality of change in more than one RNA-seq study. Through interrogating multiple lines of evidence, we showed how possible physiological functions of the pathway-orphan genes in kidney health and disease can be inferred. A substantial number of kidney-expressed genes remain ‘pathway-orphan’ while displaying clear signals of kidney relevance, such as enrichment of expression in the kidney, significant modulation in CKD, and genetic associations with kidney function. Directing mechanistic studies into this overlooked gene group might broaden our biological understanding of kidney physiology and highlight novel disease drivers.

INTRODUCTION

Giving biological pathway context is paramount to biomedical research but pathway analysis is limited by the pre-existing pathway entries in annotation databases, such as Gene Ontology (GO), Reactome, and Kyoto Encyclopedia of Genes and Genomes (KEGG) [1]. Despite the human genome being published more than 20 years ago, there are still protein-coding genes with no pathway membership in annotation databases as functional studies typically focus on well-studied genes [2]. Such inequality in research attention gets amplified with time as the understudied genes remain overlooked and the ‘rich get richer’ with pathway and functional annotations [1].

This attention bias is commonplace in biological research and is also exhibited by publication trends—95% of all life science publications focus on only 5 000 human proteins [1–4]. Factors perpetuating this inequality are limited prior knowledge, biased and unequal annotation databases, reagent availability, funding granted to projects with preceding evidence, risk-averse researchers, preference for experimentally well-accessible genes, and a shortage of large-scale perturbation studies [3, 5].

Disease understanding and drug discovery for chronic kidney disease (CKD) has been historically difficult due to its heterogeneous etiologies, which may involve more than 100 pathways in different crosstalk architectures across diagnoses [6, 7]. This molecular complexity, however, while a challenge, also allows the use of high-throughput or omics data from large patient cohorts to extensively identify potential disease drivers and candidate drug targets [6, 8]. These omics studies allow researchers to form hypotheses for follow-up mechanistic studies, by triaging and prioritizing genes with available evidence of their possible pathophysiological role [1].

In this study, we hence set out to systematically identify and investigate protein-coding genes without any known biological pathway annotations in the kidneys of both healthy individuals and CKD patients. We in silico characterize and compare these ‘pathway-orphan’ genes to ‘pathway-annotated’ and housekeeping genes. Our hypothesis is that examining this previously overlooked group of genes may lead to new insights into the disease pathogenesis of CKD.

MATERIALS AND METHODS

All databases were accessed on 26 January 2024, and the versions utilized are provided in Supplementary Table 1.

A list of protein-coding genes (N=19,258) was downloaded from the HUGO Gene Nomenclature Committee (HGNC) and annotated using a comprehensive collection of pathway databases: Gene Ontology (GO) Biological Process (BP) and Molecular Function (MF) terms, Reactome, Kyoto Encyclopedia of Genes and Genomes (KEGG), Hallmark and Curated gene sets from Molecular Signatures Database (MSigDB), WikiPathways, Ingenuity Pathway Analysis (IPA), and 28 other pathway repositories in ConsensusPathDB (summarized in Supplementary Table 2).

Gene length information was extracted from GENCODE, while protein information (length, mass, domains and families) were downloaded from UniProt and InterPro. DeepLoc 2.0 was additionally used to predict subcellular location. The counts of gene-associated publications indexed on PubMed were extracted using the gene2pubmed file obtained from the NCBI website (https://ftp.ncbi.nlm.nih.gov/gene/DATA/gene2pubmed.gz).

Expression levels of the protein-coding genome in healthy kidneys were established using RNA-seq data retrieved from Genotype-Tissue Expression (GTEx) v8 and Human Protein Atlas (HPA) v23. The presence of detectable gene expression in the kidney was defined as TPM > 0 in at least 75% of GTEx samples (63 out of 85 cortex samples and 3 out of 4 medulla samples). To compute a measure of kidney-specific enrichment of expression as compared to other tissues, the per-tissue mean expression for each gene was subjected to Z-transformation across tissues and then to a second Z-transformation across genes to bring all Z-scores to the same scale.

A list of 407 housekeeping genes (genes ubiquitously expressed in all tissue and cell types) was extracted from literature [9]. Pathway-annotated genes were defined as protein-coding genes that are neither pathway-orphan nor housekeeping genes. A list of genes with protein products enriched in the kidney was also extracted from HPA. Modulation of the expression in disease was tested using publicly available CKD kidney tissue RNA-seq transcriptomics datasets (GSE142025, GSE175759, GSE197307, GSE98422, Levin et al [10]). DESeq2 was used to perform differential expression analysis contrasting CKD vs control, and genes with false discovery rate (FDR) < 0.05 were considered differentially expressed [11]. Pathway-orphan genes with case-control differential expression in at least two CKD studies were further in silico characterized.

Expression during fetal kidney development was determined using the DESCARTES and Bgee databases and a study by Lindström et al [12]. GWAS Catalog was used to query for genetic associations, HumanBase for predicted functions, BioGrid for known interactors, Alliance of Genome Resources (AGR) for orthologs in other species, and Pharos for other predicted annotations. Nephroseq v5 was used to check for pathway-orphan differential expression and gene expression correlation with clinical variables, such as glomerular filtration rate (GFR), proteinuria, and blood urea nitrogen (BUN). GeneHancer was used to determine if pathway-orphan genes have regulatory enhancers with single-nucleotide polymorphisms (SNPs) that affect kidney function phenotypes. NephQTL2 and GTEx were queried to check for cis-eQTLs present within 1 megabase pairs of the pathway-orphan genes in the kidney.

All analyses were performed in R v4.1.3 [13].

RESULTS

Pathway-orphan genes identification and general characteristics

Following our analysis strategy (Figure 1), we reveal that almost two percent of the protein-coding human genome—286 genes—do not belong to any of the currently known pathways or gene sets. We now refer to these unannotated genes as ‘pathway-orphan genes’.

Figure 1.
Figure 1. Analysis scheme.

The flowchart shows the stepwise process of pathway-orphan genes identification and annotation in the chronic kidney disease context.

Pathway-orphan genes are distributed across all 24 human chromosomes, with the X chromosome harboring the highest number comprising over 6% of its total gene count (Figure 2A). As expected, pathway-orphan genes have substantially fewer associated publications compared to both pathway-annotated and housekeeping genes, with 17% having zero publications on PubMed (covering all articles indexed on PubMed) (Figure 2B).

Figure 2.
Figure 2. Pathway-orphan genes characteristics and trends.

(A) Chromosomal location of pathway-orphan genes. The barplot shows the pathway-orphan genes distribution across all 24 chromosomes, as reported by HGNC. The x-axis corresponds to the chromosome number and the y-axis to the pathway-orphan gene count. (B) Boxplots show the values and distribution of the number of publications indexed on PubMed of the three protein-coding gene types. The y-axis corresponds to the log10-transformed publication count. (C) Boxplots show the values and distribution of the gene lengths of the three protein-coding gene types. The y-axis corresponds to the log10-transformed gene length in base pair units. (D) Boxplots show the values and distribution of the protein lengths of the three protein-coding gene types. The y-axis corresponds to the log10-transformed protein length in number of amino acids. (E) Boxplots show the values and distribution of the protein mass of the three protein-coding gene types. The y-axis corresponds to the log10-transformed protein mass in Daltons. (F) The barplot shows how the different protein-coding gene types are localized in each subcellular compartment percentage-wise over total gene count of each type. (G) Donut chart shows the most common domains in pathway-orphan proteins. DUF: Domain of Unknown Functions.

With respect to gene length, pathway-orphan genes are generally significantly shorter, with mean length of 23.82 kilobase pairs (kbp), compared to housekeeping genes (mean: 39.10 kbp, P < 0.001) and pathway-annotated genes (mean: 69.92 kbp, P < 0.001). Consequently, the same tendency is observed in protein length where the mean length of pathway-orphan proteins is 320 amino acids long, compared to 514 (P < 0.001) and 582 (P < 0.001) for housekeeping and pathway-annotated proteins, respectively. Pathway-orphan proteins are also lighter in molecular weight, with mean mass of 35.61 kilodaltons (kDa), whereas housekeeping proteins have a mean of 57.50 kDa (P < 0.001) and 64.77 kDa (P < 0.001) for pathway-annotated proteins. Figures 2C-E show the value distribution of these attributes.

For subcellular localization of the expressed proteins, pathway-orphan genes are predicted to be predominantly cytoplasmic (N = 110) and nuclear localization comes in close second (N = 101). Distributions of protein localization for the three classes of protein-coding genes are similar and are shown in Figure 2F. Some genes are assumed to be present in more than one subcellular compartment and the localization breakdown is in Supplementary Table 3. Due to the limitations of DeepLoc, 120 protein-coding genes have no available subcellular localization prediction (26 pathway-orphan genes and 94 pathway-annotated genes).

Functional domain prediction could be accomplished for 151 out of the 286 pathway-orphan proteins. Out of these, 123 have at least one identified domain, while 28 have domains of unknown functions (Figure 2G). A full summary of the domain predictions of the pathway-orphan proteins is presented in Supplementary Table 4.

Pathway-orphan genes in the healthy kidney

Analysis of the reference kidney expression data from GTEx reveals that 97 pathway-orphan genes have detectable expression levels in the kidney tissue, 80 in the cortex and 92 in the medulla, with an overlap of 75. Supplementary Table 5 provides transcript expression and protein information on the 97 kidney-expressed pathway-orphan genes.

Kidney-expressed pathway-orphan genes retain the same characteristics of pathway-orphan genes in general, being significantly shorter in gene and protein length, with lighter protein molecular weight, predominantly cytoplasmic localization, and similar trends for predicted protein domains and families. Interestingly, pathway-orphan genes in the kidney have significantly lower mRNA expression levels as compared to pathway-annotated or housekeeping genes both in the cortex (Figure 3A) and medulla (Figure 3B). Of note, five pathway-orphan genes show specific enrichment of expression (specificity z-score > 2) [14] in the kidney cortex (LY6L, C10orf106, MYOCOS, FAM240A, and SEC14L6) and six in the medulla (ERVV-2, FRG2C, MYOCOS, SEC14L6, FAM240A, and MED14OS) as compared to other tissues and organs. Figures 3C-D show the distribution of tissue-specific enrichment, with known kidney-enriched genes as positive control.

Figure 3.
Figure 3. Kidney-expressed pathway-orphan genes characteristics and trends.

(A, B) Boxplots show the values and distribution of the kidney gene expression values from GTEx of the three protein-coding gene types in the (A) cortex and (B) medulla. The y-axis corresponds to the log10-transformed transcript expression in TPM or transcripts-per-million. (C,D) Density plots show the distribution of the (C) cortex and (D) medulla enrichment z-scores calculated from GTEx of the three protein-coding gene types with kidney-enriched transcripts as positive control. The Z-scores of kidney-enriched genes extracted from the Human Protein Atlas were also plotted as positive control.

Pathway-orphan genes modulation in CKD

Analysis of RNA-seq data from five independent studies identifies 34 pathway-orphan genes differentially expressed in the kidney tissues of CKD patients compared to healthy controls (Figure 4). Out of these, 10 pathway-orphan genes show concordant modulation directionality in two or more studies, 6 of which are downregulated (C17orf107, C22orf31, C2CD4D, FAM229A, MED14OS, and STPG2) and 4 that are upregulated (C10orf105, IQANK1, NBPF26, and SEC14L6) in disease vs control. In a dataset with micro-dissected kidney tissues, C22orf31, FAM229A, and STPG2 are downregulated and SEC14L6 is upregulated in both glomeruli and tubulointerstitium of CKD patients. Focusing on the 10 pathway-orphan genes that have concordant modulation in multiple studies, Table 1 summarizes the evidence on the potential role of these pathway in CKD pathogenesis. Table 2 shows differential expression in CKD and kidney measures correlation evidence extracted from Nephroseq.

Figure 4.
Figure 4. Dysregulation of pathway-orphan genes in chronic kidney disease.

Dotplot shows the differentially expressed pathway-orphan genes when comparing diseased and healthy control kidney tissues. Only statistically significant values are shown. Size of the dots reflects the significance level, color indicates directionality of change. Genes in bold are the pathway-orphan genes with case-control differential expression in two or more studies.

View this table:
Table 1. Summary of kidney-relevant evidence for the pathway-orphan genes robustly modulated in CKD.
View this table:
Table 2. Compiled in silico evidence of pathway-orphan genes in CKD from Nephroseq.

Genetic association evidence for pathway-orphan genes in CKD

The GWAS Catalog was systematically queried for associations with kidney-relevant phenotypes and the results are summarized in Table 3. These genetic associations are based on predominantly European ancestry cohorts, such as the UK Biobank. Lists of cis-eQTLs within 1mbp of the pathway-orphan genes are provided in Supplementary Table 7. Using the significance threshold of P=5E-8, we have found four cis-eQTLs for SEC14L6 in the glomerular tissue and 27 SNPs in the tubulointerstitial compartment based on NephQTL2 data.

View this table:
Table 3. Variants mapping to pathway-orphan genes and their associated kidney-relevant traits.

DISCUSSION

Pathway or gene set enrichment analysis has long become a routine step in the ‘-omics’ data bioinformatics analysis pipeline to infer the biological significance of genes or proteins found to be dysregulated in disease [1]. However, pathway analysis relies on predefined pathway annotations found in existing databases, thus making the method dependent on previously established pathway memberships. Moreover, some genes have been allocated to multiple pathways (e.g., complement factor C2 is included in 1 466 pathways and gene sets from ConsensusPathDB) resulting in their overrepresentation and inflated statistics, while some others are completely lacking pathway membership (‘pathway-orphan’ genes) and thus are routinely filtered out from the subsequent analysis. This study aimed to shed light on these previously overlooked genes in the kidney and their potential as an untapped source of novel disease biological insights.

We performed a data-driven exploration of the functionally enigmatic protein-coding genome focusing on those expressed in healthy kidneys and dysregulated in CKD, and observed that despite being uncharacterized, pathway-orphan genes have potential physiological roles in the kidney and pathophysiological function in CKD. Signals of kidney relevance include enriched expression in the kidney, significant modulation in CKD, and genetic associations with kidney function. These support our hypothesis that pathway-orphan genes warrant further functional characterization and mechanistic studies to determine how they can be involved in disease pathogenesis.

We report that many pathway-orphan genes are located in the X chromosome. Due to its high-sequence identity regions and transmission pattern that causes technical artifacts, the X chromosome is typically excluded from genomic analyses [43]. However, the X chromosome comprises 5% of the human genome and it also bears several genes involved in the immune response [43]. Moreover, the X chromosome is an important factor in disease, often contributing in a sex-specific manner. Pathway-orphan genes in the X chromosome thus deserve further studies, such as MED14OS as it is also dysregulated in CKD.

While pathway-orphan proteins are generally shorter and lighter than housekeeping and pathway-annotated proteins, they should not be overlooked as recent studies have shown how microproteins—proteins typically less than 100 amino acids in length—are involved in diverse biological processes. These microproteins encoded by short open reading frames can function as allosteric regulators of other proteins, signaling molecules, or effector proteins [44]. Eleven kidney-expressed pathway-orphan proteins are less than 100 amino acids long, which includes SMIM35. Other small integral membrane proteins (e.g. SMIM22, SMIM43) have been found to have roles as regulators of cytoskeletal organization and glucose transport [44].

Out of the 10 pathway-orphan genes dysregulated in CKD, C10orf105 and SEC14L6 have multiple lines of evidence for potential involvement in CKD. Both are upregulated at the mRNA level in diseased kidneys and are genetically linked to glomerular filtration rate, either through its regulatory enhancers (C10orf105) or its own genetic variants (SEC14L6). While genome-wide significant associations of gene variants to CKD or other kidney function traits do not confirm causation of pathology, they shed light on genes that might be involved in pathogenesis despite lack of mechanistic understanding. Single-pass membrane proteins are transmembrane proteins that play key roles in signal transduction, cell communication, immunity, transport, and energy conversion [45]. C10orf105 has a transmembrane helix domain and is a single-pass membrane protein, which suggests its potential function. For SEC14L6, the presence of the CRAL-TRIO lipid binding domain could indicate a role in lipid traffic [26]. SEC14L6 also has significant cis-eQTLs from glomeruli and tubulointerstitium derived from kidney biopsies of nephrotic syndrome patients.

Looking at which protein family these pathway-orphan genes proteins belong to and which protein domains they contain also gives us clues on their possible roles. The C2 domain in C2CD4D, in general, senses the cellular lipid microenvironment and can regulate lipid signal transduction and membrane trafficking, which could be its role in the kidney [46]. IQANK1 has ankyrin repeats, commonly found in Notch receptors, which is involved in CKD pathophysiology [47]. There are seven pathway-orphan genes with experimentally validated interactors and one (SEC14L6) with a known ubiquitination site, suggesting that these proteins are part of protein networks and could serve a function.

Majority of the pathway-orphan genes are also expressed in the human fetal kidney (varying expression levels from Week 9 to 21) [12, 16], hinting about their involvement in kidney development. Defective kidney development has been linked to kidney disease and multiple evidence suggest that nephron deficiencies from possibly developmental defects can lead to CKD [48].

The remaining 24 pathway-orphan genes dysregulated in CKD could still have potential pathophysiological roles and are also worth investigating. Increasing the number of CKD RNA-seq studies analyzed could reveal additional evidence for these pathway-orphan genes.

We acknowledge several limitations of this study. Databases and annotation repositories are dynamic and are updated periodically, therefore, the validity of the results is restricted to the accession dates. Database identifiers can also be inaccurate, redundant, or obsolete. Some missingness accounts for a small proportion of protein-coding genes without expression level and functional data resulting in them having been missed in our analyses. Lastly, RNA-seq-based gene expression detection is dependent upon the sequencing depth, increasing this might potentially reveal additional low-abundant transcripts.

Despite these limitations, our data-driven investigation of this overlooked group of genes uses an analytical approach to characterize them. Previous studies have reiterated that data-driven hypotheses and analyses can help alleviate the annotation inequality we observe in all these databases [1, 3]. Our design should also apply to other tissues and organs, as well as other diseases.

In conclusion, a number of routinely overlooked ‘pathway-orphan’ genes are likely to have a plausible link to kidney physiology or play a role in disease. We hope that our results will help rectify a bias in functional gene annotation and spark subsequent hypothesis-testing investigations, for example, through mechanistic experimental or Mendelian Randomization type of studies, that are warranted to verify their potential roles in CKD pathogenesis.

Data availability statement

All data analyzed in this study are publicly available and cited within the manuscript.

Funding

This research was supported by the European Union Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No. 860977 titled TrainCKDis.

Authors‘ contributions

Dianne Acoba: Data curation, Formal analysis, Investigation, Methodology, Resources, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

Sophie Gunnarsson: Data curation, Formal analysis, Investigation, Methodology, Resources, Software, Visualization, Writing – review & editing

Anna Reznichenko: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing – review & editing

Conflict of interest statement

DA, SG, and AR are AstraZeneca employees. DA is an industrial PhD student at AstraZeneca.

Acknowledgements

Initial preliminary results from this work were previously shared as an abstract at the ASN Kidney Week 2022 (November 3 – 6, 2022) and as an abstract and poster at the World Congress of Nephrology 2023 (March 30 – April 2, 2023).

Footnotes

  • Abstract was updated for easier reading.

REFERENCES

  1. 1.↵
    Haynes WA, Tomczak A, Khatri P. Gene annotation bias impedes biomedical research. Sci Rep. 2018;8(1):1362. doi:10.1038/s41598-018-19333-x
    OpenUrlCrossRef
  2. 2.↵
    Sinha S, Eisenhaber B, Jensen LJ, Kalbuaji B, Eisenhaber F. Darkness in the Human Gene and Protein Function Space: Widely Modest or Absent Illumination by the Life Science Literature and the Trend for Fewer Protein Function Discoveries Since 2000. Proteomics. 2018;18(21-22):e1800093. doi:10.1002/pmic.201800093
    OpenUrlCrossRef
  3. 3.↵
    Stoeger T, Gerlach M, Morimoto RI, Nunes Amaral LA. Large-scale investigation of the reasons why potentially important genes are ignored. PLoS Biol. 2018;16(9):e2006643. doi:10.1371/journal.pbio.2006643
    OpenUrlCrossRefPubMed
  4. 4.↵
    Wood V, Lock A, Harris MA, Rutherford K, Bahler J, Oliver SG. Hidden in plain sight: what remains to be discovered in the eukaryotic proteome? Open Biol. 2019;9(2):180241. doi:10.1098/rsob.180241
    OpenUrlCrossRefPubMed
  5. 5.↵
    Rocha JJ, Jayaram SA, Stevens TJ, et al. Functional unknomics: Systematic screening of conserved genes of unknown function. PLoS Biol. 2023;21(8):e3002222. doi:10.1371/journal.pbio.3002222
    OpenUrlCrossRef
  6. 6.↵
    Buvall L, Menzies RI, Williams J, et al. Selecting the right therapeutic target for kidney disease. Front Pharmacol. 2022;13. doi:10.3389/fphar.2022.971065
    OpenUrlCrossRef
  7. 7.↵
    Martini S, Nair V, Keller BJ, et al. Integrative Biology Identifies Shared Transcriptional Networks in CKD. J Am Soc Nephrol. 2014;25(11). doi: 10.1681/ASN.2013080906
    OpenUrlAbstract/FREE Full Text
  8. 8.↵
    Cisek K, Krochmal M, Klein J, Mischak H. The application of multi-omics and systems biology to identify therapeutic targets in chronic kidney disease. Nephrol Dial Transplant. 2016;31(12):2003–11. doi:10.1093/ndt/gfv364
    OpenUrlCrossRefPubMed
  9. 9.↵
    Zhu J, He F, Song S, Wang J, Yu J. How many human genes can be defined as housekeeping with current expression data? BMC Genomics. 2008;9(1):172. doi:10.1186/1471-2164-9-172
    OpenUrlCrossRefPubMed
  10. 10.↵
    Levin A, Reznichenko A, Witasp A, et al. Novel insights into the disease transcriptome of human diabetic glomeruli and tubulointerstitium. Nephrol Dial Transpl. 2020;35(12):2059–72. doi:10.1093/ndt/gfaa121
    OpenUrlCrossRef
  11. 11.↵
    Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550. doi:10.1186/s13059-014-0550-8
    OpenUrlCrossRefPubMed
  12. 12.↵
    Lindstrom NO, McMahon JA, Guo J, et al. Conserved and Divergent Features of Human and Mouse Kidney Organogenesis. J Am Soc Nephrol. 2018;29(3):785–805. doi:10.1681/ASN.2017080887
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    R Core Team. R: The R Project for Statistical Computing. R: The R Project for Statistical Computing. http://r-project.org
  14. 14.↵
    Ryaboshapkina M, Hammar M. Tissue-specific genes as an underutilized resource in drug discovery. Sci Rep. 2019;9(1). doi:10.1038/s41598-019-43829-9
    OpenUrlCrossRefPubMed
  15. 15.
    UniProt C. UniProt: the Universal Protein Knowledgebase in 2023. Nucleic Acids Res. 2023;51(D1):D523–31. doi:10.1093/nar/gkac1052
    OpenUrlCrossRefPubMed
  16. 16.↵
    Cao J, O’Day DR, Pliner HA, et al. A human cell atlas of fetal gene expression. Science. 2020;370(6518):eaba7721. doi:10.1126/science.aba7721
    OpenUrlAbstract/FREE Full Text
  17. 17.
    Greene CS, Krishnan A, Wong AK, et al. Understanding multicellular function and disease with human tissue-specific networks. Nat Genet. 2015;47(6):569–76. doi:10.1038/ng.3259
    OpenUrlCrossRefPubMed
  18. 18.
    Fishilevich S, Nudel R, Rappaport N, et al. GeneHancer: genome-wide integration of enhancers and target genes in GeneCards. Database Oxf. 2017;2017. doi:10.1093/database/bax028
    OpenUrlCrossRefPubMed
  19. 19.
    Oughtred R, Rust J, Chang C, et al. The BioGRID database: A comprehensive biomedical resource of curated protein, genetic, and chemical interactions. Protein Sci. 2021;30(1):187–200. doi:10.1002/pro.3978
    OpenUrlCrossRefPubMed
  20. 20.
    Bastian FB, Roux J, Niknejad A, et al. The Bgee suite: integrated curated expression atlas and comparative transcriptomics in animals. Nucleic Acids Res. 2021;49(D1):D831–47. doi:10.1093/nar/gkaa793
    OpenUrlCrossRefPubMed
  21. 21.
    Sollis E, Mosaku A, Abid A, et al. The NHGRI-EBI GWAS Catalog: knowledgebase and deposition resource. Nucleic Acids Res. 2023;51(D1):D977–85. doi:10.1093/nar/gkac1010
    OpenUrlCrossRef
  22. 22.
    Kelleher KJ, Sheils TK, Mathias SL, et al. Pharos 2023: an integrated resource for the understudied human proteome. Nucleic Acids Res. 2023;51(D1):D1405–16. doi:10.1093/nar/gkac1033
    OpenUrlCrossRef
  23. 23.
    He D, Pan C, Zhao Y, et al. Exome-wide screening identifies novel rare risk variants for bone mineral density. Osteoporos Int. 2023;34(5):965–75. doi:10.1007/s00198-023-06710-0
    OpenUrlCrossRef
  24. 24.
    Alliance of Genome Resources C. Harmonizing model organism data in the Alliance of Genome Resources. Genetics. 2022;220(4). doi:10.1093/genetics/iyac022
    OpenUrlCrossRefPubMed
  25. 25.
    Sinnott-Armstrong N, Tanigawa Y, Amar D, et al. Genetics of 35 blood and urine biomarkers in the UK Biobank. Nat Genet. 2021;53(2):185–94. doi:10.1038/s41588-020-00757-z
    OpenUrlCrossRefPubMed
  26. 26.↵
    Paysan-Lafosse T, Blum M, Chuguransky S, et al. InterPro in 2022. Nucleic Acids Res. 2023;51(D1):D418–27. doi:10.1093/nar/gkac993
    OpenUrlCrossRefPubMed
  27. 27.
    Lonsdale J, Thomas J, Salvatore M, et al. The Genotype-Tissue Expression (GTEx) project. Nat Genet. 2013;45(6):580–5. doi:10.1038/ng.2653
    OpenUrlCrossRefPubMed
  28. 28.
    Liu H, Doke T, Guo D, et al. Epigenomic and transcriptomic analyses define core cell types, genes and targetable mechanisms for kidney disease. Nat Genet. 2022;54(7):950–62. doi:10.1038/s41588-022-01097-w
    OpenUrlCrossRef
  29. 29.
    Han SK, McNulty MT, Benway CJ, et al. Mapping genomic regulation of kidney disease and traits through high-resolution and interpretable eQTLs. Nat Commun. 2023;14(1):2229. doi:10.1038/s41467-023-37691-7
    OpenUrlCrossRef
  30. 30.
    Hood RL, Schenkel LC, Nikkel SM, et al. The defining DNA methylation signature of Floating-Harbor Syndrome. Sci Rep. 2016;6(1). doi:10.1038/srep38803
    OpenUrlCrossRef
  31. 31.
    Mariani LH, Eddy S, AlAkwaa FM, et al. Precision nephrology identified tumor necrosis factor activation variability in minimal change disease and focal segmental glomerulosclerosis. Kidney Int. 2023;103(3):565–79. doi:10.1016/j.kint.2022.10.023
    OpenUrlCrossRef
  32. 32.
    Sampson MG, Robertson CC, Martini S, et al. Integrative Genomics Identifies Novel Associations with APOL1 Risk Genotypes in Black NEPTUNE Subjects. J Am Soc Nephrol. 2016;27(3). doi:10.1681/ASN.2014111131
    OpenUrlAbstract/FREE Full Text
  33. 33.
    Nakagawa S, Nishihara K, Miyata H, et al. Molecular Markers of Tubulointerstitial Fibrosis and Tubular Cell Damage in Patients with Chronic Kidney Disease. PLOS ONE. 2015;10(8):e0136994. doi:10.1371/journal.pone.0136994
    OpenUrlCrossRefPubMed
  34. 34.
    Woroniecka KI, Park ASD, Mohtat D, Thomas DB, Pullman JM, Susztak K. Transcriptome Analysis of Human Diabetic Kidney Disease. Diabetes. 2011;60(9):2354–69. doi:10.2337/db10-1181
    OpenUrlAbstract/FREE Full Text
  35. 35.
    Ju W, Greene CS, Eichinger F, et al. Defining cell-type specificity at the transcriptional level in human disease. Genome Res. 2013;23(11):1862–73. doi:10.1101/gr.155697.113
    OpenUrlAbstract/FREE Full Text
  36. 36.
    Ju W, Nair V, Smith S, et al. Tissue transcriptome-driven identification of epidermal growth factor as a chronic kidney disease biomarker. Sci Transl Med. 2015;7(316):316ra193. doi:10.1126/scitranslmed.aac7071
    OpenUrlAbstract/FREE Full Text
  37. 37.
    Neusser MA, Lindenmeyer MT, Moll AG, et al. Human Nephrosclerosis Triggers a Hypoxia-Related Glomerulopathy. Am J Pathol. 2010;176(2):594–607. doi:10.2353/ajpath.2010.090268
    OpenUrlCrossRefPubMedWeb of Science
  38. 38.
    Schmid H, Boucherot A, Yasuda Y, et al. Modular Activation of Nuclear Factor-κB Transcriptional Programs in Human Diabetic Nephropathy. Diabetes. 2006;55(11):2993–3003. doi:10.2337/db06-0477
    OpenUrlAbstract/FREE Full Text
  39. 39.
    Berthier CC, Bethunaickan R, Gonzalez-Rivera T, et al. Cross-Species Transcriptional Network Analysis Defines Shared Inflammatory Responses in Murine and Human Lupus Nephritis. J Immunol. 2012;189(2):988–1001. doi:10.4049/jimmunol.1103031
    OpenUrlAbstract/FREE Full Text
  40. 40.
    Lin BM, Nadkarni GN, Tao R, et al. Genetics of Chronic Kidney Disease Stages Across Ancestries: The PAGE Study. Front Genet. 2019;10. doi: 10.3389/fgene.2019.00494
    OpenUrlCrossRef
  41. 41.
    He D, Liu H, Wei W, et al. A longitudinal genome-wide association study of bone mineral density mean and variability in the UK Biobank. Osteoporos Int. 2023;34(11):1907–16. doi:10.1007/s00198-023-06852-1
    OpenUrlCrossRef
  42. 42.
    Al-Khelaifi F, Diboun I, Donati F, et al. Metabolic GWAS of elite athletes reveals novel genetically-influenced metabolites associated with athletic performance. Sci Rep. 2019;9(1):19889. doi:10.1038/s41598-019-56496-7
    OpenUrlCrossRef
  43. 43.↵
    Pinto BJ, O’Connor B, Schatz MC, Zarate S, Wilson MA, Hahn M. Concerning the eXclusion in human genomics: the choice of sex chromosome representation in the human genome drastically affects the number of identified variants. G3 Genes Genomes Genet. 2023;13(10). doi:10.1093/g3journal/jkad169
    OpenUrlCrossRef
  44. 44.↵
    Hassel KR, Brito-Estrada O, Makarewich CA. Microproteins: Overlooked regulators of physiology and disease. iScience. 2023;26(6):106781. doi:10.1016/j.isci.2023.106781
    OpenUrlCrossRef
  45. 45.↵
    Pogozheva ID, Lomize AL. Evolution and adaptation of single-pass transmembrane proteins. Biochim Biophys Acta Biomembr. 2018;1860(2):364–77. doi:10.1016/j.bbamem.2017.11.002
    OpenUrlCrossRef
  46. 46.↵
    Leonard TA. C2 Domain Proteins. In: Kretsinger RH, Uversky VN, Permyakov EA, eds. Encyclopedia of Metalloproteins. Springer; 2013:309–18.
  47. 47.↵
    Edeling M, Ragi G, Huang S, Pavenstadt H, Susztak K. Developmental signalling pathways in renal fibrosis: the roles of Notch, Wnt and Hedgehog. Nat Rev Nephrol. 2016;12(7):426–39. doi:10.1038/nrneph.2016.54
    OpenUrlCrossRefPubMed
  48. 48.↵
    Schedl A. Renal abnormalities and their developmental origin. Nat Rev Genet. 2007;8(10):791–802. doi:10.1038/nrg2205
    OpenUrlCrossRefPubMedWeb of Science
Back to top
Posted July 11, 2024.
Download PDF

Supplementary Material

Data/Code
Email
Kidney ‘pathway-orphan’ genes as a potential overlooked source of novel insights into kidney health and disease
Dianne Acoba, Sophie Gunnarsson, Anna Reznichenko
medRxiv 2024.07.09.24309717; doi: https://doi.org/10.1101/2024.07.09.24309717
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area