Multiple omics levels of chronic lymphocytic leukemia.

Chronic lymphocytic leukemia (CLL) is a lymphoproliferative malignancy characterized by the proliferation of functionally mature but incompetent B cells. It is the most prevalent type of leukemia in Western populations, accounting for approximately 25% of new leukemia cases. While recent advances, such as ibrutinib and venetoclax treatment have improved patient outlook, aggressive forms of CLL such as Richter transformation still pose a significant challenge. This discrepancy may be due to the heterogeneity of factors contributing to CLL development at multiple -omics levels. However, information on the omics of CLL is fragmented, hindering multi-omics-based research into potential treatment options. To address this, we aggregated and presented a selection of important aspects of various omics levels of the disease in this review. The purpose of the present literature analysis is to portray examples of CLL studies from different omics levels, including genomics, epigenomics, transcriptomics, epitranscriptomics, proteomics, epiproteomics, metabolomics, glycomics and lipidomics, as well as those identified by multi-omics approaches. The review includes the list of 102 CLL-associated genes with relevant genomics information. While single-omics studies yield substantial and useful data, they omit a significant level of complex biological interplay present in the disease. As multi-omics studies integrate several different layers of data, they may be better suited for complex diseases such as CLL and have thus far yielded promising results. Future multi-omics studies may assist clinicians in improved treatment choices based on CLL subtypes as well as allow the identification of novel biomarkers and targets for treatments.

Transcriptomic signatures for human male infertility.

Introduction: Male infertility is a common, complex disorder. A better understanding of pathogenesis and etiology is needed for timely diagnosis and treatment. The aim of this study, therefore, was to identify genes involved in the pathogenesis of idiopathic male infertility based on data from transcriptomic level supported with data from genomic level. Materials and methods: First, we performed whole gene expression analysis in 20 testis biopsy samples of patients with severely impaired (10) and normal spermatogenesis (10). Further, we have performed systematic review of comparable male infertility studies and overlapped the most significantly expressed genes identified in our study with the most differentially expressed genes from selected studies. Gene Ontology analysis and KEGG functional enrichment have been performed with Enrichr analysis tool. Additionally, we have overlapped these genes with the genes where rare variants have been identified previously. Results: In 10 patients with severely impaired spermatogenesis and 10 controls, we identified more than 1,800 differentially expressed genes (p < 0.001). With the systematic review of three previously performed microarray studies that have met inclusion criteria we identified 257 overlapped differentialy expressed genes (144 downregulated and 113 upregulated). Intersection of genes from transcriptomic studies with genes with identified rare variants revealed a total of 7 genes linked with male infertility phenotype (CYP11A1, CYP17A1, RSPH3, TSGA10, AKAP4, CCIN, NDNF). Conclusion: Our comprehensive study highlighted the role of four genes in pathogenesis of male infertility and provided supporting evidence for three promising candidate genes which dysfunction may result in a male infertility disorder.

Shared Genetic Risk Factors Between Cancer and Cardiovascular Diseases.

Cancer and cardiovascular diseases (CVD) account for approximately 27.5 million deaths every year. While they share some common environmental risk factors, their shared genetic risk factors are not yet fully understood. The aim of the present study was to aggregate genetic risk factors associated with the comorbidity of cancer and CVDs. For this purpose, we: (1) created a catalog of genes associated with cancer and CVDs, (2) visualized retrieved data as a gene-disease network, and (3) performed a pathway enrichment analysis. We performed screening of PubMed database for literature reporting genetic risk factors in patients with both cancer and CVD. The gene-disease network was visualized using Cytoscape and the enrichment analysis was conducted using Enrichr software. We manually reviewed the 181 articles fitting the search criteria and included 13 articles in the study. Data visualization revealed a highly interconnected network containing a single subnetwork with 56 nodes and 146 edges. Genes in the network with the highest number of disease interactions were JAK2, TTN, TET2, and ATM. The pathway enrichment analysis revealed that genes included in the study were significantly enriched in DNA damage repair (DDR) pathways, such as homologous recombination. The role of DDR mechanisms in the development of CVDs has been studied in previously published research; however, additional functional studies are required to elucidate their contribution to the pathophysiology to CVDs.

MicroRNAs in Leukemias: A Clinically Annotated Compendium.

Leukemias are a group of malignancies of the blood and bone marrow. Multiple types of leukemia are known, however reliable treatments have not been developed for most leukemia types. Furthermore, even relatively reliable treatments can result in relapses. MicroRNAs (miRNAs) are a class of short, noncoding RNAs responsible for epigenetic regulation of gene expression and have been proposed as a source of potential novel therapeutic targets for leukemias. In order to identify central miRNAs for leukemia, we conducted data synthesis using two databases: miRTarBase and DISNOR. A total of 137 unique miRNAs associated with 16 types of leukemia were retrieved from miRTarBase and 86 protein-coding genes associated with leukemia were retrieved from the DISNOR database. Based on these data, we formed a visual network of 248 miRNA-target interactions (MTI) between leukemia-associated genes and miRNAs associated with ≥4 leukemia types. We then manually reviewed the literature describing these 248 MTIs for interactions identified in leukemia studies. This manually curated data was then used to visualize a network of 64 MTIs identified in leukemia patients, cell lines and animal models. We also formed a visual network of miRNA-leukemia associations. Finally, we compiled leukemia clinical trials from the ClinicalTrials database. miRNAs with the highest number of MTIs were miR-125b-5p, miR-155-5p, miR-181a-5p and miR-19a-3p, while target genes with the highest number of MTIs were TP53, BCL2, KIT, ATM, RUNX1 and ABL1. The analysis of 248 MTIs revealed a large, highly interconnected network. Additionally, a large MTI subnetwork was present in the network visualized from manually reviewed data. The interconnectedness of the MTI subnetwork suggests that certain miRNAs represent central disease molecules for multiple leukemia types. Additional studies on miRNAs, their target genes and associated biological pathways are required to elucidate the therapeutic potential of miRNAs in leukemia.

Transcriptomics of receptive endometrium in women with sonographic features of adenomyosis.

BACKGROUND: Women with uterine adenomyosis seeking assisted reproduction have been associated with compromised endometrial receptivity to embryo implantation. To understand the mechanisms involved in this process, we aimed to compare endometrial transcriptome profiles during the window of implantation (WOI) between women with and without adenomyosis. METHODS: We obtained endometrial biopsies LH-timed to the WOI from women with sonographic features of adenomyosis (n=10) and controls (n=10). Isolated RNA samples were subjected to RNA sequencing (RNA-seq) by the Illumina NovaSeq 6000 platform and endometrial receptivity classification with a molecular tool for menstrual cycle phase dating (beREADY®, CCHT). The program language R and Bioconductor packages were applied to analyse RNA-seq data in the setting of the result of accurate endometrial dating. To suggest robust candidate pathways, the identified differentially expressed genes (DEGs) associated with the adenomyosis group in the receptive phase were further integrated with 151, 173 and 42 extracted genes from published studies that were related to endometrial receptivity in healthy uterus, endometriosis and adenomyosis, respectively. Enrichment analyses were performed using Cytoscape ClueGO and CluePedia apps. RESULTS: Out of 20 endometrial samples, 2 were dated to the early receptive phase, 13 to the receptive phase and 5 to the late receptive phase. Comparison of the transcriptomics data from all 20 samples provided 909 DEGs (p<0.05; nonsignificant after adjusted p value) in the adenomyosis group but only 4 enriched pathways (Bonferroni p value < 0.05). The analysis of 13 samples only dated to the receptive phase provided suggestive 382 DEGs (p<0.05; nonsignificant after adjusted p value) in the adenomyosis group, leading to 33 enriched pathways (Bonferroni p value < 0.05). These included pathways were already associated with endometrial biology, such as "Expression of interferon (IFN)-induced genes" and "Response to IFN-alpha". Data integration revealed pathways indicating a unique effect of adenomyosis on endometrial molecular organization (e.g., "Expression of IFN-induced genes") and its interference with endometrial receptivity establishment (e.g., "Extracellular matrix organization" and "Tumour necrosis factor production"). CONCLUSIONS: Accurate endometrial dating and RNA-seq analysis resulted in the identification of altered response to IFN signalling as the most promising candidate of impaired uterine receptivity in adenomyosis.

Integration and Visualization of Regulatory Elements and Variations of the EPAS1 Gene in Human.

Endothelial PAS domain-containing protein 1 (EPAS1), also HIF2α, is an alpha subunit of hypoxia-inducible transcription factor (HIF), which mediates cellular and systemic response to hypoxia. EPAS1 has an important role in the transcription of many hypoxia-responsive genes, however, it has been less researched than HIF1α. The aim of this study was to integrate an increasing number of data on EPAS1 into a map of diverse OMICs elements. Publications, databases, and bioinformatics tools were examined, including Ensembl, MethPrimer, STRING, miRTarBase, COSMIC, and LOVD. The EPAS1 expression, stability, and activity are tightly regulated on several OMICs levels to maintain complex oxygen homeostasis. In the integrative EPAS1 map we included: 31 promoter-binding proteins, 13 interacting miRNAs and one lncRNA, and 16 post-translational modifications regulating EPAS1 protein abundance. EPAS1 has been associated with various cancer types and other diseases. The development of neuroendocrine tumors and erythrocytosis was shown to be associated with 11 somatic and 20 germline variants. The integrative map also includes 12 EPAS1 target genes and 27 interacting proteins. The study introduced the first integrative map of diverse genomics, transcriptomics, proteomics, regulomics, and interactomics data associated with EPAS1, to enable a better understanding of EPAS1 activity and regulation and support future research.

Integrative Map of HIF1A Regulatory Elements and Variations.

Hypoxia-inducible factor (HIF) family of transcription factors (HIF1A, EPAS1, and HIF3A) are regulators of the cellular response to hypoxia. They have been shown to be involved in development of various diseases such as cancer, diabetes, and erythrocytosis. A complete map of connections between HIF family of genes with various omics types has not yet been developed. The main aim of the present analysis was to construct the integrative map of genomic elements associated with HIF1A gene and prioritize potentially deleterious variants. Various genomic databases and bioinformatics tools were used, including Ensembl, MirTarBase, STRING, Cytoscape, MethPrimer, CADD, SIFT, and UALCAN. Integrative HIF1A gene map was visualized and includes transcriptional and post-transcriptional regulators, downstream targets, and genetic variants. One CpG island overlaps transcription start site of the HIF1A gene. Out of over 450 missense variants, four have predicted deleterious effect on protein function by at least five bioinformatics tools. Currently there are 85 miRNAs reported to target HIF1A. HIF1A downstream targets include protein-coding genes, long noncoding RNAs, and microRNAs (hypoxamiRs). The study presents the first integration of heterogeneous molecular interactions associated with HIF1A gene enabling a holistic view of the gene and lays the groundwork for supplementing the data in the future.

Top Trends in Multiomics Research: Evaluation of 52 Published Studies and New Ways of Thinking Terminology and Visual Displays.

Multiomics study designs have significantly increased understanding of complex biological systems. The multiomics literature is rapidly expanding and so is their heterogeneity. However, the intricacy and fragmentation of omics data are impeding further research. To examine current trends in multiomics field, we reviewed 52 articles from PubMed and Web of Science, which used an integrated omics approach, published between March 2006 and January 2021. From studies, data regarding investigated loci, species, omics type, and phenotype were extracted, curated, and streamlined according to standardized terminology, and summarized in a previously developed graphical summary. Evaluated studies included 21 omics types or applications of omics technology such as genomics, transcriptomics, metabolomics, epigenomics, environmental omics, and pharmacogenomics, species of various phyla including human, mouse, Arabidopsis thaliana, Saccharomyces cerevisiae, and various phenotypes, including cancer and COVID-19. In the analyzed studies, diverse methods, protocols, results, and terminology were used and accordingly, assessment of the studies was challenging. Adoption of standardized multiomics data presentation in the future will further buttress standardization of terminology and reporting of results in systems science. This shall catalyze, we suggest, innovation in both science communication and laboratory medicine by making available scientific knowledge that is easier to grasp, share, and harness toward medical breakthroughs.

Identification of Variants Associated With Rare Hematological Disorder Erythrocytosis Using Targeted Next-Generation Sequencing Analysis.

An erythrocytosis is present when the red blood cell mass is increased, demonstrated as elevated hemoglobin and hematocrit in the laboratory evaluation. Congenital predispositions for erythrocytosis are rare, with germline variants in several genes involved in oxygen sensing (VHL, EGLN1, and EPAS1), signaling for hematopoietic cell maturation (EPOR and EPO), and oxygen transfer (HBB, HBA1, HBA2, and BPGM) that were already associated with the eight congenital types (ECYT1-8). Screening for variants in known congenital erythrocytosis genes with classical sequencing approach gives a correct diagnosis for only up to one-third of the patients. The genetic background of erythrocytosis is more heterogeneous, and additional genes involved in erythropoiesis and iron metabolism could have a putative effect on the development of erythrocytosis. This study aimed to detect variants in patients with yet unexplained erythrocytosis using the next-generation sequencing (NGS) approach, targeting genes associated with erythrocytosis and increased iron uptake and implementing the diagnostics of congenital erythrocytosis in Slovenia. Selected 25 patients with high hemoglobin, high hematocrit, and no acquired causes were screened for variants in the 39 candidate genes. We identified one pathogenic variant in EPAS1 gene and three novel variants with yet unknown significance in genes EPAS1, JAK2, and SH2B3. Interestingly, a high proportion of patients were heterozygous carriers for two variants in HFE gene, otherwise pathogenic for the condition of iron overload. The association between the HFE variants and the development of erythrocytosis is not clearly understood. With a targeted NGS approach, we determined an actual genetic cause for the erythrocytosis in one patient and contributed to better management of the disease for the patient and his family. The effect of variants of unknown significance on the enhanced production of red blood cells needs to be further explored with functional analysis. This study is of great significance for the improvement of diagnosis of Slovenian patients with unexplained erythrocytosis and future research on the etiology of this rare hematological disorder.

Identification of Novel RNA Binding Proteins Influencing Circular RNA Expression in Hepatocellular Carcinoma.

Circular RNAs (circRNAs) are increasingly recognized as having a role in cancer development. Their expression is modified in numerous cancers, including hepatocellular carcinoma (HCC); however, little is known about the mechanisms of their regulation. The aim of this study was to identify regulators of circRNAome expression in HCC. Using publicly available datasets, we identified RNA binding proteins (RBPs) with enriched motifs around the splice sites of differentially expressed circRNAs in HCC. We confirmed the binding of some of the candidate RBPs using ChIP-seq and eCLIP datasets in the ENCODE database. Several of the identified RBPs were found to be differentially expressed in HCC and/or correlated with the overall survival of HCC patients. According to our bioinformatics analyses and published evidence, we propose that NONO, PCPB2, PCPB1, ESRP2, and HNRNPK are candidate regulators of circRNA expression in HCC. We confirmed that the knocking down the epithelial splicing regulatory protein 2 (ESRP2), known to be involved in the maintenance of the adult liver phenotype, significantly changed the expression of candidate circRNAs in a model HCC cell line. By understanding the systemic changes in transcriptome splicing, we can identify new proteins involved in the molecular pathways leading to HCC development and progression.

Erythrocytosis: genes and pathways involved in disease development.

Erythrocytosis is a blood disorder characterised by an increased red blood cell mass. The most common causes of erythrocytosis are acquired and caused by diseases and conditions that are accompanied by hypoxaemia or overproduction of erythropoietin. More rarely, erythrocytosis has a known genetic background, such as for polycythaemia vera and familial erythrocytosis. The majority of cases of polycythaemia vera are associated with acquired variants in JAK2, while familial erythrocytosis is a group of congenital disorders. Familial erythrocytosis type 1 is associated with hypersensitivity to erythropoietin (variants in EPOR), types 2-5 with defects in oxygen-sensing pathways (variants in VHL, EGLN1, EPAS1, EPO), and types 6-8 with an increased affinity of haemoglobin for oxygen (variants in HBB, HBA1, HBA2, BPGM). Due to a heterogenic genetic background, the causes of disease are not fully discovered and in more than 70% of patients the condition remains labelled idiopathic.The transfer of next-generation sequencing into clinical practice is becoming a reality enabling detection of various variants in a single rapid test. In this review, we describe the current research on erythrocytosis gene variants and the mechanisms associated with disease development, along with the currently used diagnostic tests.

Determining the Molecular Background of Endometrial Receptivity in Adenomyosis.

BACKGROUND: Adenomyosis is a gynaecological condition with limited evidence of negative impact to endometrial receptivity. It is commonly associated with endometriosis, which has been shown to alter endometrial expression patterns. Therefore, the candidate genes identified in endometriosis could serve as a source to study endometrial function in adenomyosis. METHODS: Transcripts/proteins associated with endometrial receptivity in women with adenomyosis or endometriosis and healthy women were obtained from publications and their nomenclature was adopted according to the HUGO Gene Nomenclature Committee (HGNC). Retrieved genes were analysed for enriched pathways using Cytoscape/Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) and Reactome tools to prioritise candidates for endometrial receptivity. These were used for validation on women with (n = 9) and without (n = 13) adenomyosis. RESULTS: Functional enrichment analysis of 173, 42 and 151 genes associated with endometriosis, adenomyosis and healthy women, respectively, revealed signalling by interleukins and interleukin-4 and interleukin-13 signalling pathways, from which annotated LIF, JUNB, IL6, FOS, IL10 and SOCS3 were prioritised. Selected genes showed downregulated expression levels in adenomyosis compared to the control group, but without statistical significance. CONCLUSION: This is the first integrative study providing putative candidate genes and pathways characterising endometrial receptivity in women with adenomyosis in comparison to healthy women and women with endometriosis.

Genetic Variations of Ultraconserved Elements in the Human Genome.

Ultraconserved elements (UCEs) are among the most popular DNA markers for phylogenomic analysis. In at least three of five placental mammalian genomes (human, dog, cow, mouse, and rat), 2189 UCEs of at least 200 bp in length that are identical have been identified. Most of these regions have not yet been functionally annotated, and their associations with diseases remain largely unknown. This is an important knowledge gap in human genomics with regard to UCE roles in physiologically critical functions, and by extension, their relevance for shared susceptibilities to common complex diseases across several mammalian organisms in the event of their polymorphic variations. In the present study, we remapped the genomic locations of these UCEs to the latest human genome assembly, and examined them for documented polymorphisms in sequenced human genomes. We identified 29,983 polymorphisms within analyzed UCEs, but revealed that a vast majority exhibits very low minor allele frequencies. Notably, only 112 of the identified polymorphisms are associated with a phenotype in the Ensembl genome browser. Through literature analyses, we confirmed associations of 37 (i.e., out of the 112) polymorphisms within 23 UCEs with 25 diseases and phenotypic traits, including, muscular dystrophies, eye diseases, and cancers (e.g., familial adenomatous polyposis). Most reports of UCE polymorphism-disease associations appeared to be not cognizant that their candidate polymorphisms were actually within UCEs. The present study offers strategic directions and knowledge gaps for future computational and experimental work so as to better understand the thus far intriguing and puzzling role(s) of UCEs in mammalian genomes.

A Map of the microRNA Regulatory Networks Identified by Experimentally Validated microRNA-Target Interactions in Five Domestic Animals: Cattle, Pig, Sheep, Dog, and Chicken.

Domestic animals are members of the broader ecological context, in which humans are situated. Yet, genomics and systems science research have lagged behind and been relatively underappreciated in domestic animals compared to human genetics/genomics. Harnessing big data calls for omics data mapping studies in a broad range of mammals. To this end, microRNAs (miRNAs) regulate posttranscriptional expression of target genes, hence, governing different biological pathways and physiological processes. The knowledge of miRNA regulatory networks and maps is important for understanding regulation of gene expression and functions in both humans and domestic animals. However, complete miRNA regulatory networks have not yet been described in all species, particularly in domestic animals. We report here an original analysis so as to map the miRNA regulatory networks in domestic animals based on miRNA-target interactions (MTIs). Validated MTIs for five species; cattle, pig, sheep, dog, and chicken were extracted from the miRTarBase. miRNA regulomes were visualized using the Cytoscape software. The data in cattle, chicken, and pig were sufficient to visualize networks, identify central molecules, and subnetworks associated with the same phenotype; however, the MTI data in dog and sheep are still limited. We found several hub genes with large number of interactions, for example, 1 miRNA (bta-miR-17-5p) interacting with 27 genes and 7 miRNAs interacting with the same gene (tumor necrosis factor [TNF]) in cattle. In addition, two single-nucleotide polymorphisms were identified within the seed region of a previously demonstrated MTI, namely, between HMGB3 (high mobility group box 3) gene and bta-miR-17-5p. In summary, this miRNA regulome mapping study will enable and guide further studies of genome function in mammals with a view to applications in human as well as veterinary medicine. Furthermore, these miRNA regulomes can help to clarify fundamental pathways in cell biology and reveal molecular insights on phenotypic trait variability in common complex diseases and response phenotypes of drugs or other health interventions for precision medicine in the future.

Vascular endothelial growth factor (VEGF)-related polymorphisms rs10738760 and rs6921438 are not risk factors for proliferative diabetic retinopathy (PDR) in patients with type 2 diabetes mellitus (T2DM).

Vascular endothelial growth factor (VEGF) is an important regulator of angiogenesis and has been investigated as a candidate gene in a number of conditions, including diabetes and its microvascular complications (e.g., retinopathy and nephropathy). Several VEGF-related polymorphisms have been shown to contribute to nearly half of the variability in circulating VEGF levels in healthy individuals. Our aim was to assess the association between VEGF-related rs10738760 and rs6921438 polymorphisms and proliferative diabetic retinopathy (PDR) in Slovenian patients with type 2 diabetes mellitus (T2DM). We also investigated the effect of these polymorphisms on VEGF receptor 2 (VEGFR-2) expression in fibrovascular membranes (FVMs) from patients with PDR. This case-control study enrolled 505 unrelated patients with T2DM: 143 diabetic patients with PDR as a study group, and 362 patients with T2DM of >10 years duration and with no clinical signs of PDR as a control group. Patient clinical and laboratory data were obtained from their medical records. rs10738760 and rs6921438 polymorphisms were genotyped using TaqMan SNP Genotyping assay. VEGFR-2 expression was assessed by immunohistochemistry in 20 FVMs from patients with PDR, and numerical areal density of VEGFR-2-positive cells was calculated. The occurrence of PDR was 1.7 times higher in diabetic patients carrying GA genotype of rs6921438 compared to patients with GG genotype, with a borderline statistical significance (OR = 1.7, 95% CI = 1.00 - 2.86, p = 0.05). In addition, A allele of rs6921438 was associated with increased VEGFR-2 expression in FVMs from PDR patients. However, we observed no association between AA genotype of rs6921438 nor between rs10738760 variants and PDR, indicating that the two polymorphisms are not genetic risk factors for PDR.

Molecular Mechanisms of Syndromic Cryptorchidism: Data Synthesis of 50 Studies and Visualization of Gene-Disease Network.

Background: Cryptorchidism is one of the most frequent congenital birth defects in male children and is present in 2-4% of full-term male births. It has several possible health effects including reduced fertility, increased risk for testicular neoplasia, testicular torsion, and psychological consequences. Cryptorchidism is often diagnosed as comorbid; copresent with other diseases. It is also present in clinical picture of several syndromes. However, this field has not been systematically studied. The aim of the present study was to catalog published cases of syndromes which include cryptorchidism in the clinical picture and associated genomic information. Methods: The literature was extracted from Public/Publisher MEDLINE and Web of Science databases, using the keywords including: syndrome, cryptorchidism, undescended testes, loci, and gene. The obtained data was organized in a table according to the previously proposed standardized data format. The results of the study were visually represented using Gephi and karyotype view. Results: Fifty publications had sufficient data for analysis. Literature analysis resulted in 60 genomic loci, associated with 44 syndromes that have cryptorchidism in clinical picture. Genomic loci included 38 protein-coding genes and 22 structural variations containing microdeletions and microduplications. Loci, associated with syndromic cryptorchidism are located on 16 chromosomes. Visualization of retrieved data is presented in a gene-disease network. Conclusions: The study is ongoing and further studies will be needed to develop a complete catalog with the data from upcoming publications. Additional studies will also be needed for revealing of molecular mechanisms associated with syndromic cryptorchidism and revealing complete diseasome network.

Harnessing Omics Big Data in Nine Vertebrate Species by Genome-Wide Prioritization of Sequence Variants with the Highest Predicted Deleterious Effect on Protein Function.

Harnessing the genomics big data requires innovation in how we extract and interpret biologically relevant variants. Currently, there is no established catalog of prioritized missense variants associated with deleterious protein function phenotypes. We report in this study, to the best of our knowledge, the first genome-wide prioritization of sequence variants with the most deleterious effect on protein function (potentially deleterious variants [pDelVars]) in nine vertebrate species: human, cattle, horse, sheep, pig, dog, rat, mouse, and zebrafish. The analysis was conducted using the Ensembl/BioMart tool. Genes comprising pDelVars in the highest number of examined species were identified using a Python script. Multiple genomic alignments of the selected genes were built to identify interspecies orthologous potentially deleterious variants, which we defined as the "ortho-pDelVars." Genome-wide prioritization revealed that in humans, 0.12% of the known variants are predicted to be deleterious. In seven out of nine examined vertebrate species, the genes encoding the multiple PDZ domain crumbs cell polarity complex component (MPDZ) and the transforming acidic coiled-coil containing protein 2 (TACC2) comprise pDelVars. Five interspecies ortho-pDelVars were identified in three genes. These findings offer new ways to harness genomics big data by facilitating the identification of functional polymorphisms in humans and animal models and thus provide a future basis for optimization of protocols for whole genome prioritization of pDelVars and screening of orthologous sequence variants. The approach presented here can inform various postgenomic applications such as personalized medicine and multiomics study of health interventions (iatromics).

Sarcoidosis Related Novel Candidate Genes Identified by Multi-Omics Integrative Analyses.

Sarcoidosis is a multifactorial systemic disease characterized by granulomatous inflammation and greatly impacting on global public health. The etiology and mechanisms of sarcoidosis are not fully understood. Recent high-throughput biological research has generated vast amounts of multi-omics big data on sarcoidosis, but their significance remains to be determined. We sought to identify novel candidate regions, and genes consistently altered in heterogeneous omics studies so as to reveal the underlying molecular mechanisms. We conducted a comprehensive integrative literature analysis on global data on sarcoidosis, including genomic, transcriptomic, proteomic, and phenomic studies. We performed positional integration analysis of 38 eligible datasets originating from 17 different biological layers. Using the integration interval length of 50 kb, we identified 54 regions reaching significance value p ≤ 0.0001 and 15 regions with significance value p ≤ 0.00001, when applying more stringent criteria. Secondary literature analysis of the top 20 regions, with the most significant accumulation of signals, revealed several novel candidate genes for which associations with sarcoidosis have not yet been established, but have considerable support for their involvement based on omic data. These new plausible candidate genes include NELFE, CFB, EGFL7, AGPAT2, FKBPL, NRC3, and NEU1. Furthermore, annotated data were prepared to enable custom visualization and browsing of these sarcoidosis related omics evidence in the University of California Santa Cruz (UCSC) Genome Browser. Further multi-omics approaches are called for sarcoidosis biomarkers and diagnostic and therapeutic innovation. Our approach for harnessing multi-omics data and the findings presented herein reflect important steps toward understanding the etiology and underlying pathological mechanisms of sarcoidosis.

Genome-wide screening for smallest regions of overlaps in cryptorchidism.

Cryptorchidism is a urogenital abnormality associated with increased rates of testicular neoplasia and impaired spermatogenesis. The field is facing expansion of genomics data; however, it lacks protocols for biomarker prioritization. Identification of smallest regions of overlap (SRO) presents an approach for candidate gene identification but has not yet been systematically conducted in cryptorchidism. The aim of this study was to conduct a genome-wide screening for SRO (GW-SRO) associated with cryptorchidism development. We updated the Cryptorchidism Gene Database to version 3, remapped genomic coordinates of loci from older assemblies to the GRCh38 and performed genome-wide screening for overlapping regions associated with cryptorchidism risk. A total of 73 chromosomal loci (68 involved in chromosomal mutations and five copy number variations) described in 37 studies associated with cryptorchidism risk in humans were used for SRO identification. Analysis resulted in 18 SRO, based on deletions, duplications, inversions, derivations and copy number variations. Screening for SRO was challenging owing to heterogeneous reporting of genomic locations. To our knowledge, this is the first GW-SRO study for cryptorchidism and it presents the basis for further narrowing of critical regions for cryptorchidism and planning functional experiments. The developed protocol could also be applied to other multifactorial diseases.

Insights from Ion Binding Site Network Analysis into Evolution and Functions of Proteins.

Many biological phenomena can be represented as complex networks. Using a protein binding site comparison approach, we generated a network of ion binding sites on the scale of all known protein structures from the Protein Data Bank. We found that this ion binding site similarity network is scale-free, indicating a network in which a few ion binding site scaffolds are the network hubs, and these are connected to hundreds of nodes, whereas the vast majority of nodes have only a few neighbors. Enrichment and statistical analysis of the network components and communities yielded insights into underlying processes from the functional and the structural perspective. Largest components and communities were observed to be closely related to basic metabolic processes and some of the most common structural folds, which, from the evolutionary point of view, indicates that they may be the oldest ones. Further, we derived the first comprehensive map of ion interchangeability, based on binding site similarity. Several highly interchangeable protein-ion binding site pairs emerged (e.g., Ca2+ and Mg2+ ), as well as structurally distinct ones. The constructed network of ion binding site similarities will aid in understanding the general principles of protein-ion binding sites structure, function and evolution. We demonstrate potential uses of the network on proteins involved in cancer development and immune response, where individual ions play prominent roles in disease development.