Gene Ontology Enrichment Analysis of Renal Agenesis: Improving Prenatal Molecular Diagnosis

Uni- or bilateral renal agenesis (RA) is a commonly occurring major congenital anomaly impacting fetal and neonatal outcomes. Since the etiology is highly heterogeneous, our aim was to provide a logically structured approach by highlighting the genes in which variants have been identified to be associated with RA and to define the pathways involved in this type of abnormal kidney development. We used Phenolyzer to collect a list of all the genes known as causative for RA. Using ClueGO gene enrichment analysis, we classified the relationship between these genes and the biological processes defined by gene ontology. We identified 287 genes and 69 groups of enriched biological processes. About 50% included pathways directly related to the development of urogenital organ tissues. Several ciliary, axis specification, hindgut development, and endocrine pathways were enriched, which may relate to different clinical presentations of RA. Our gene ontology enrichment analysis shows that genes representing distinct biological pathways are significantly enriched. This knowledge will lead to an improved molecular diagnosis in clinical care when applying genome-wide sequencing approaches. The findings will also allow to further study the biological pathways involved in RA and to identify novel candidate genes and pathways.

Identification of DNA methylation markers for early detection of CRC indicates a role for nervous system-related genes in CRC

Purpose Colonoscopy and the fecal immunochemical test (FIT) are currently the most widely used screening modalities for colorectal cancer (CRC), however, both with their own limitations. Here we aim to identify and validate stool-based DNA methylation markers for the early detection of CRC and investigate the biological pathways prone to DNA methylation. Methods DNA methylation marker discovery was performed using The Cancer Genome Atlas (TCGA) colon adenocarcinoma data set consisting of normal and primary colon adenocarcinoma tissue. The performance of the five best candidate markers and a previously identified marker, NDRG4, was evaluated on tissues and whole stool samples of healthy subjects and CRC patients using quantitative MSP assays. The results were compared and combined with FIT data. Finally, pathway and gene ontology enrichment analyses were performed using ToppFun, GOrilla and clusterProfiler. Results GDNF, HAND2, SLC35F3, SNAP91 and SORCS1 were ranked as the best performing markers. Gene combinations of all five markers, NDRG4 and FIT were evaluated to establish the biomarker panel with the highest diagnostic potential, resulting in the identification of GDNF/SNAP91/NDRG4/FIT as the best performing marker panel. Pathway and gene ontology enrichment analyses revealed that genes associated with the nervous system were enriched in the set of best performing CRC-specific biomarkers. Conclusion In silico discovery analysis using TCGA-derived data yielded a novel DNA-methylation-based assay for the early detection of CRC, potentially improving current screening modalities. Additionally, nervous system-related pathways were enriched in the identified genes, indicating an epigenetic regulation of neuronal genes in CRC.

Essential Amino Acid Enrichment and Positive Selection Highlight Endosymbiont's Role in a Global Virus-Vectoring Pest

ABSTRACT Host-associated microbes display remarkable convergence in genome repertoire resulting from selection to supplement missing host functions. Nutritional supplementation has been proposed in the verrucomicrobial endosymbiont Xiphinematobacter sp., which lives within a globally widespread group of plant-parasitic nematodes that vector damaging nepoviruses to plants. Only one genome sequence has been published from this symbiont, leaving unanswered questions about its diversity, host range, role, and selective pressures within its hosts. Because its hosts are exceptionally resistant to culturing, this symbiont is best studied through advanced genomic approaches. To analyze the role of Xiphinematobacter sp. in its host, sequencing was performed on nematode communities, and then genomes were extracted for comparative genomics, gene ontology enrichment tests, polymorphism analysis, de Bruijn-based genome-wide association studies, and tests of pathway- and site-specific selection on genes predicted play a role in the symbiosis. Results showed a closely clustered set of Xiphinematobacter isolates with reduced genomes of ∼917 kbp, for which a new species was proposed. Symbionts shared only 2.3% of genes with outgroup Verrucomicrobia, but comparative analyses showed high conservation of all 10 essential amino acid (EAA) biosynthesis pathways plus several vitamin pathways. These findings were supported by gene ontology enrichment tests and high polymorphisms in these pathways compared with background. Genome-wide association analysis confirmed high between-species fixation of alleles with significant functional enrichment for EAA and thiamine synthesis. Strong positive selection was detected on sites within these pathways, despite several being under increased purifying selection. Together, these results suggest that supplementation of EAAs missing in the host diet may drive this widespread symbiosis. IMPORTANCE Xiphinematobacter spp. are distinctly evolved intracellular symbionts in the phylum Verrucomicrobia, which includes the important human gut-associated microbe Akkermansia muciniphila and many highly abundant free-living soil microbes. Like Akkermansia sp., Xiphinematobacter sp. is obligately associated with the gut of its hosts, which in this case consists of a group of plant-parasitic nematodes that are among the top 10 most destructive species to global agriculture, by vectoring plant viruses. This study examined the hypothesis that the key to this symbiont’s stable evolutionary association with its host is through provisioning nutrients that its host cannot make that may be lacking in the nematode’s plant phloem diet, such as essential amino acids and several vitamins. The significance of our research is in demonstrating, using population genomics, the signatures of selective pressure on these hypothesized roles to ultimately learn how this independently evolved symbiont functionally mirrors symbionts of phloem-feeding insects.

A multi-omics study on cutaneous and uveal melanoma

AIM: To present the multi-omics landscape of cutaneous melanoma (CM) and uveal melanoma (UM) from The Cancer Genome Atlas (TCGA). METHODS: The differentially expressed genes (DEGs) between CM and UM were found and integrated into a gene ontology enrichment analysis. Besides, the differentially expressed miRNAs were also identified. We also compared the methylation level of CM with UM and identified the differentially methylated regions to integrate with the DEGs to display the relationship between the gene expression and DNA methylation. The differentially expressed transcription factors (TFs) were identified. RESULTS: Though CM had more mutational burden than UM, they shared several similarities such as the same rankings in diverse variant types. Except GNAQ and GNA11, the other top 18 mutated genes of the combined group were mostly detected in CM instead of UM. On the transcriptomic level, 4610 DEGs were found and integrated into a gene ontology enrichment analysis. We also identified 485 differentially expressed miRNAs. The methylation analysis showed that UM had a significantly higher methylation level than CM. The integration of differentially methylated regions and DEGs demonstrated that most DEGs were downregulated in UM and the hypo- and hypermethylation presented no obvious difference within these DEGs. Finally, 116 hypermethylated TFs and 114 hypomethylated TFs were identified as differentially expressed TFs in CM when compared with UM. CONCLUSION: This multi-omics study on comparing CM with UM confirms that they differ in all analyzed levels. Of notice, the results also offer new insights with implications for elucidating certain unclear problems such as the distinct role of epithelial mesenchymal transition in two melanomas, the different metastatic routes of CM and UM and the liver tropism of metastatic UM.

Molecular signatures of the rediae, cercariae and adult stages in the complex life cycles of parasitic flatworms (Digenea: Psilostomatidae)

Background Parasitic flatworms (Trematoda: Digenea) represent one of the most remarkable examples of drastic morphological diversity among the stages within a life cycle. Which genes are responsible for extreme differences in anatomy, physiology, behavior, and ecology among the stages? Here we report a comparative transcriptomic analysis of parthenogenetic and amphimictic generations in two evolutionary informative species of Digenea belonging to the family Psilostomatidae. Methods In this study the transcriptomes of rediae, cercariae and adult worm stages of Psilotrema simillimum and Sphaeridiotrema pseudoglobulus, were sequenced and analyzed. High-quality transcriptomes were generated, and the reference sets of protein-coding genes were used for differential expression analysis in order to identify stage-specific genes. Comparative analysis of gene sets, their expression dynamics and Gene Ontology enrichment analysis were performed for three life stages within each species and between the two species. Results Reference transcriptomes for P. simillimum and S. pseudoglobulus include 21,433 and 46,424 sequences, respectively. Among 14,051 orthologous groups (OGs), 1354 are common and specific for two analyzed psilostomatid species, whereas 13 and 43 OGs were unique for P. simillimum and S. pseudoglobulus, respectively. In contrast to P. simillimum, where more than 60% of analyzed genes were active in the redia, cercaria and adult worm stages, in S. pseudoglobulus less than 40% of genes had such a ubiquitous expression pattern. In general, 7805 (36.41%) and 30,622 (65.96%) of genes were preferentially expressed in one of the analyzed stages of P. simillimum and S. pseudoglobulus, respectively. In both species 12 clusters of co-expressed genes were identified, and more than a half of the genes belonging to the reference sets were included into these clusters. Functional specialization of the life cycle stages was clearly supported by Gene Ontology enrichment analysis. Conclusions During the life cycles of the two species studied, most of the genes change their expression levels considerably, consequently the molecular signature of a stage is not only a unique set of expressed genes, but also the specific levels of their expression. Our results indicate unexpectedly high level of plasticity in gene regulation between closely related species. Transcriptomes of P. simillimum and S. pseudoglobulus provide high quality reference resource for future evolutionary studies and comparative analyses.