Genetic variants in mRNA untranslated regions

Machine learning-based classification reveals distinct clusters of non-coding genomic allelic variations associated with Erm-mediated antibiotic resistance

The erythromycin resistance RNA methyltransferase (erm) confers cross-resistance to all therapeutically important macrolides, lincosamides, and streptogramins (MLS phenotype). The expression of erm is often induced by the macrolide-mediated ribosome stalling in the upstream co-transcribed leader sequence, thereby triggering a conformational switch of the intergenic RNA hairpins to allow the translational initiation of erm. We investigated the evolutionary emergence of the upstream erm regulatory elements and the impact of allelic variation on erm expression and the MLS phenotype. Through systematic profiling of the upstream regulatory sequences across all known erm operons, we observed that specific erm subfamilies, such as ermB and ermC, have independently evolved distinct configurations of small upstream ORFs and palindromic repeats. A population-wide genomic analysis of the upstream ermB regions revealed substantial non-random allelic variation at numerous positions. Utilizing machine learning-based classification coupled with RNA structure modeling, we found that many alleles cooperatively influence the stability of alternative RNA hairpin structures formed by the palindromic repeats, which, in turn, affects the inducibility of ermB expression and MLS phenotypes. Subsequent experimental validation of 11 randomly selected variants demonstrated an impressive 91% accuracy in predicting MLS phenotypes. Furthermore, we uncovered a mixed distribution of MLS-sensitive and MLS-resistant ermB loci within the evolutionary tree, indicating repeated and independent evolution of MLS resistance. Taken together, this study not only elucidates the evolutionary processes driving the emergence and development of MLS resistance but also highlights the potential of using non-coding genomic allele data to predict antibiotic resistance phenotypes.

IMPORTANCE

Antibiotic resistance (AR) poses a global health threat as the efficacy of available antibiotics has rapidly eroded due to the widespread transmission of AR genes. Using Erm-dependent MLS resistance as a model, this study highlights the significance of non-coding genomic allelic variations. Through a comprehensive analysis of upstream regulatory elements within the erm family, we elucidated the evolutionary emergence and development of AR mechanisms. Leveraging population-wide machine learning (ML)-based genomic analysis, we transformed substantial non-random allelic variations into discernible clusters of elements, enabling precise prediction of MLS phenotypes from non-coding regions. These findings offer deeper insight into AR evolution and demonstrate the potential of harnessing non-coding genomic allele data for accurately predicting AR phenotypes.

CGG repeats in the human FMR1 gene regulate mRNA localization and cellular stress in developing neurons

The human genome has many short tandem repeats, yet the normal functions of these repeats are unclear. The 5' untranslated region (UTR) of the fragile X messenger ribonucleoprotein 1 (FMR1) gene contains polymorphic CGG repeats, the length of which has differing effects on FMR1 expression and human health, including the neurodevelopmental disorder fragile X syndrome. We deleted the CGG repeats in the FMR1 gene (0CGG) in human stem cells and examined the effects on differentiated neurons. 0CGG neurons have altered subcellular localization of FMR1 mRNA and protein, and differential expression of cellular stress proteins compared with neurons with normal repeats (31CGG). In addition, 0CGG neurons have altered responses to glucocorticoid receptor (GR) activation, including FMR1 mRNA localization, GR chaperone HSP90α expression, GR localization, and cellular stress protein levels. Therefore, the CGG repeats in the FMR1 gene are important for the homeostatic responses of neurons to stress signals.

A DIO2 missense mutation and its impact on fetal response to PRRSV infection

BACKGROUND

Porcine reproductive and respiratory syndrome virus 2 (PRRSV-2) infection during late gestation substantially lowers fetal viability and survival. In a previous genome-wide association study, a single nucleotide polymorphism on chromosome 7 was significantly associated with probability of fetuses being viable in response to maternal PRRSV-2 infection at 21 days post maternal inoculation. The iodothyronine deiodinase 2 (DIO2) gene, located ~ 14 Kilobase downstream of this SNP, was selected as a priority candidate related to fetal susceptibility following maternal PRRSV-2 infection. Our objectives were to identify mutation(s) within the porcine DIO2 gene and to determine if they were associated with fetal outcomes after PRRSV-2 challenge. Sequencing of the DIO2, genotyping identified variants, and association of DIO2 genotypes with fetal phenotypes including DIO2 mRNA levels, viability, survival, viral loads, cortisol and thyroid hormone levels, and growth measurements were conducted.

RESULTS

A missense variant (p.Asn91Ser) was identified in the parental populations from two independent PRRSV-2 challenge trials. This variant was further genotyped to determine association with fetal PRRS outcomes. DIO2 mRNA levels in fetal heart and kidney differed by the genotypes of Asn91Ser substitution with significantly greater DIO2 mRNA expression in heterozygotes compared with wild-type homozygotes (P < 0.001 for heart, P = 0.002 for kidney). While Asn91Ser did not significantly alter fetal viability and growth measurements, interaction effects of the variant with fetal sex or trial were identified for fetal viability or crown rump length, respectively. However, this mutation was not related to dysregulation of the hypothalamic-pituitary-adrenal and thyroid axis, indicated by no differences in circulating cortisol, T4, and T3 levels in fetuses of the opposing genotypes following PRRSV-2 infection.

CONCLUSIONS

The present study suggests that a complex relationship among DIO2 genotype, DIO2 expression, fetal sex, and fetal viability may exist during the course of fetal PRRSV infection. Our study also proposes the increase in cortisol levels, indicative of fetal stress response, may lead to fetal complications, such as fetal compromise, fetal death, or premature farrowing, during PRRSV infection.

Haemostatic gene variations in cervical cancer-associated venous thrombosis: considerations for clinical strategies

Venous thromboembolism (VTE) is a life-threatening haemostatic disease frequently diagnosed among the cancer population. The Khorana Score is currently the primal risk assessment model to stratify oncological patients according to their susceptibility to VTE, however, it displays a limited performance. Meanwhile, intensive research on VTE pathophysiology in the general population has uncovered a range of single-nucleotide polymorphisms (SNPs) associated with the condition. Nonetheless, their predictive ability concerning cancer-associated thrombosis (CAT) is controversial. Cervical cancer (CC) patients undergoing chemoradiotherapy often experience VTE, which negatively affects their survival. Thus, aiming for an improvement in thromboprophylaxis, new thrombotic biomarkers, including SNPs, are currently under investigation. In this study, the predictive capability of haemostatic gene SNPs on CC-related VTE and their prognostic value regardless of VTE were explored. Six SNPs in haemostatic genes were evaluated. A total of 401 CC patients undergoing chemoradiotherapy were enrolled in a retrospective cohort study. The implications for the time to VTE occurrence and overall survival (OS) were assessed. CAT considerably impacted the CC patients' OS (log-rank test, P < 0.001). SERPINE1 rs2070682 (T > C) showed a significant association with the risk of CC-related VTE (CC/CT vs. TT, log-rank test, P = 0.002; C allele, Cox model, hazard ratio (HR) = 6.99 and P = 0.009), while F2 rs1799963 (G > A) demonstrated an important prognostic value regardless of VTE (AA/AG vs. GG, log-rank test, P = 0.020; A allele, Cox model, HR = 2.76 and P = 0.026). For the remaining SNPs, no significant associations were detected. The polymorphisms SERPINE1 rs2070682 and F2 rs1799963 could be valuable tools in clinical decision-making, aiding in thromboprophylaxis and CC management, respectively.

Biological relevance of alternative splicing in hematologic malignancies

Alternative splicing (AS) is a strictly regulated process that generates multiple mRNA variants from a single gene, thus contributing to proteome diversity. Transcriptome-wide sequencing studies revealed networks of functionally coordinated splicing events, which produce isoforms with distinct or even opposing functions. To date, several mechanisms of AS are deregulated in leukemic cells, mainly due to mutations in splicing and/or epigenetic regulators and altered expression of splicing factors (SFs). In this review, we discuss aberrant splicing events induced by mutations affecting SFs (SF3B1, U2AF1, SRSR2, and ZRSR2), spliceosome components (PRPF8, LUC7L2, DDX41, and HNRNPH1), and epigenetic modulators (IDH1 and IDH2). Finally, we provide an extensive overview of the biological relevance of aberrant isoforms of genes involved in the regulation of apoptosis (e. g. BCL-X, MCL-1, FAS, and c-FLIP), activation of key cellular signaling pathways (CASP8, MAP3K7, and NOTCH2), and cell metabolism (PKM).

Generation of Rhesus Macaque Embryos with Expanded CAG Trinucleotide Repeats in the Huntingtin Gene

Huntington's disease (HD) arises from expanded CAG repeats in exon 1 of the Huntingtin (HTT) gene. The resultant misfolded HTT protein accumulates within neuronal cells, negatively impacting their function and survival. Ultimately, HTT accumulation results in cell death, causing the development of HD. A nonhuman primate (NHP) HD model would provide important insight into disease development and the generation of novel therapies due to their genetic and physiological similarity to humans. For this purpose, we tested CRISPR/Cas9 and a single-stranded DNA (ssDNA) containing expanded CAG repeats in introducing an expanded CAG repeat into the HTT gene in rhesus macaque embryos. Analyses were conducted on arrested embryos and trophectoderm (TE) cells biopsied from blastocysts to assess the insertion of the ssDNA into the HTT gene. Genotyping results demonstrated that 15% of the embryos carried an expanded CAG repeat. The integration of an expanded CAG repeat region was successfully identified in five blastocysts, which were cryopreserved for NHP HD animal production. Some off-target events were observed in biopsies from the cryopreserved blastocysts. NHP embryos were successfully produced, which will help to establish an NHP HD model and, ultimately, may serve as a vital tool for better understanding HD's pathology and developing novel treatments.

The 5' untranslated region variant rs3811050 C/T of the interleukin-38 encoding gene is associated with susceptibility to rheumatoid arthritis in Iraqi women

BACKGROUND

Interleukin (IL)-38, the latest member of the IL-1 cytokine family, is proposed to have a pathogenic role in rheumatoid arthritis (RA). It is encoded by the IL1F10 gene, which harbors single nucleotide polymorphisms (SNPs) that may predict the risk of autoimmune diseases. Among them are 5' untranslated region (UTR) SNPs, which play a key role in post-transcriptional control, but have not been studied in Iraqi RA patients.

METHODS

Two novel IL1F10 5'UTR SNPs (rs3811050 C/T and rs3811051 T/G) were explored in RA and control women (n = 120 and 110, respectively). SNPs were genotyped using TaqMan assay. An ELISA kit was used to measure serum IL-38 concentrations.

RESULTS

A reduced risk of RA was associated with rs3811050 T allele and CT genotype (corrected probability [pc] = 0.01 and < 0.001, respectively), while there was no significant association with rs3811051. Haplotype analysis demonstrated that C-T haplotype was associated with a 1.65-fold greater risk of RA, whereas a reduced risk was linked to T-G haplotype. IL-38 concentrations were higher in patients than in controls (p < 0.001). In addition, IL-38 showed acceptable performance in distinguishing between RA and control women (p < 0.001). When IL-38 concentrations were stratified according to SNP genotypes, no significant differences were found.

CONCLUSIONS

The rs3811050 variant was more likely to affect RA susceptibility in Iraqi women, and the T allele may play a role in reducing disease risk. IL-38 concentrations were elevated in RA patients, but were not affected by the rs3811050 and rs3811051 genotypes.

Massively parallel screen uncovers many rare 3' UTR variants regulating mRNA abundance of cancer driver genes

Understanding the function of rare non-coding variants represents a significant challenge. Using MapUTR, a screening method, we studied the function of rare 3' UTR variants affecting mRNA abundance post-transcriptionally. Among 17,301 rare gnomAD variants, an average of 24.5% were functional, with 70% in cancer-related genes, many in critical cancer pathways. This observation motivated an interrogation of 11,929 somatic mutations, uncovering 3928 (33%) functional mutations in 155 cancer driver genes. Functional MapUTR variants were enriched in microRNA- or protein-binding sites and may underlie outlier gene expression in tumors. Further, we introduce untranslated tumor mutational burden (uTMB), a metric reflecting the amount of somatic functional MapUTR variants of a tumor and show its potential in predicting patient survival. Through prime editing, we characterized three variants in cancer-relevant genes (MFN2, FOSL2, and IRAK1), demonstrating their cancer-driving potential. Our study elucidates the function of tens of thousands of non-coding variants, nominates non-coding cancer driver mutations, and demonstrates their potential contributions to cancer.

Bioinformatics pipeline for the systematic mining genomic and proteomic variation linked to rare diseases: The example of monogenic diabetes

Monogenic diabetes is characterized as a group of diseases caused by rare variants in single genes. Like for other rare diseases, multiple genes have been linked to monogenic diabetes with different measures of pathogenicity, but the information on the genes and variants is not unified among different resources, making it challenging to process them informatically. We have developed an automated pipeline for collecting and harmonizing data on genetic variants linked to monogenic diabetes. Furthermore, we have translated variant genetic sequences into protein sequences accounting for all protein isoforms and their variants. This allows researchers to consolidate information on variant genes and proteins linked to monogenic diabetes and facilitates their study using proteomics or structural biology. Our open and flexible implementation using Jupyter notebooks enables tailoring and modifying the pipeline and its application to other rare diseases.

A Comprehensive Allele Specific Expression Resource for the Equine Transcriptome

Background

Allele-specific expression (ASE) analysis provides a nuanced view of cis-regulatory mechanisms affecting gene expression.

Results

An equine ASE analysis was performed, using integrated Iso-seq and short-read RNA sequencing data from four healthy Thoroughbreds (2 mares and 2 stallions) across 9 tissues from the Functional Annotation of Animal Genomes (FAANG) project. Allele expression was quantified by haplotypes from long-read data, with 42,900 allele expression events compared. Within these events, 635 (1.48%) demonstrated ASE, with liver tissue containing the highest proportion. Genetic variants within ASE events were in histone modified regions 64.2% of the time. Validation of allele-specific variants, using a set of 66 equine liver samples from multiple breeds, confirmed that 97% of variants demonstrated ASE.

Conclusions

This valuable publicly accessible resource is poised to facilitate investigations into regulatory variation in equine tissues. Our results highlight the tissue-specific nature of allelic imbalance in the equine genome.

Lung Cancer in Ever- and Never-Smokers: Findings from Multi-Population GWAS Studies

BACKGROUND

Clinical, molecular, and genetic epidemiology studies displayed remarkable differences between ever- and never-smoking lung cancer.

METHODS

We conducted a stratified multi-population (European, East Asian, and African descent) association study on 44,823 ever-smokers and 20,074 never-smokers to identify novel variants that were missed in the non-stratified analysis. Functional analysis including expression quantitative trait loci (eQTL) colocalization and DNA damage assays, and annotation studies were conducted to evaluate the functional roles of the variants. We further evaluated the impact of smoking quantity on lung cancer risk for the variants associated with ever-smoking lung cancer.

RESULTS

Five novel independent loci, GABRA4, intergenic region 12q24.33, LRRC4C, LINC01088, and LCNL1 were identified with the association at two or three populations (P < 5 × 10-8). Further functional analysis provided multiple lines of evidence suggesting the variants affect lung cancer risk through excessive DNA damage (GABRA4) or cis-regulation of gene expression (LCNL1). The risk of variants from 12 independent regions, including the well-known CHRNA5, associated with ever-smoking lung cancer was evaluated for never-smokers, light-smokers (packyear ≤ 20), and moderate-to-heavy-smokers (packyear > 20). Different risk patterns were observed for the variants among the different groups by smoking behavior.

CONCLUSIONS

We identified novel variants associated with lung cancer in only ever- or never-smoking groups that were missed by prior main-effect association studies.

IMPACT

Our study highlights the genetic heterogeneity between ever- and never-smoking lung cancer and provides etiologic insights into the complicated genetic architecture of this deadly cancer.

The Effects of De Novo Mutation on Gene Expression and the Consequences for Fitness in Chlamydomonas reinhardtii

Mutation is the ultimate source of genetic variation, the bedrock of evolution. Yet, predicting the consequences of new mutations remains a challenge in biology. Gene expression provides a potential link between a genotype and its phenotype. But the variation in gene expression created by de novo mutation and the fitness consequences of mutational changes to expression remain relatively unexplored. Here, we investigate the effects of >2,600 de novo mutations on gene expression across the transcriptome of 28 mutation accumulation lines derived from 2 independent wild-type genotypes of the green algae Chlamydomonas reinhardtii. We observed that the amount of genetic variance in gene expression created by mutation (Vm) was similar to the variance that mutation generates in typical polygenic phenotypic traits and approximately 15-fold the variance seen in the limited species where Vm in gene expression has been estimated. Despite the clear effect of mutation on expression, we did not observe a simple additive effect of mutation on expression change, with no linear correlation between the total expression change and mutation count of individual MA lines. We therefore inferred the distribution of expression effects of new mutations to connect the number of mutations to the number of differentially expressed genes (DEGs). Our inferred DEE is highly L-shaped with 95% of mutations causing 0-1 DEG while the remaining 5% are spread over a long tail of large effect mutations that cause multiple genes to change expression. The distribution is consistent with many cis-acting mutation targets that affect the expression of only 1 gene and a large target of trans-acting targets that have the potential to affect tens or hundreds of genes. Further evidence for cis-acting mutations can be seen in the overabundance of mutations in or near differentially expressed genes. Supporting evidence for trans-acting mutations comes from a 15:1 ratio of DEGs to mutations and the clusters of DEGs in the co-expression network, indicative of shared regulatory architecture. Lastly, we show that there is a negative correlation with the extent of expression divergence from the ancestor and fitness, providing direct evidence of the deleterious effects of perturbing gene expression.

Dissection of the E8 locus in two early maturing Canadian soybean populations

Soybean [Glycine max (L.) Merr.] is a short-day crop for which breeders want to expand the cultivation range to more northern agro-environments by introgressing alleles involved in early reproductive traits. To do so, we investigated quantitative trait loci (QTL) and expression quantitative trait loci (eQTL) regions comprised within the E8 locus, a large undeciphered region (~7.0 Mbp to 44.5 Mbp) associated with early maturity located on chromosome GM04. We used a combination of two mapping algorithms, (i) inclusive composite interval mapping (ICIM) and (ii) genome-wide composite interval mapping (GCIM), to identify major and minor regions in two soybean populations (QS15524F2:F3 and QS15544RIL) having fixed E1, E2, E3, and E4 alleles. Using this approach, we identified three main QTL regions with high logarithm of the odds (LODs), phenotypic variation explained (PVE), and additive effects for maturity and pod-filling within the E8 region: GM04:16,974,874-17,152,230 (E8-r1); GM04:35,168,111-37,664,017 (E8-r2); and GM04:41,808,599-42,376,237 (E8-r3). Using a five-step variant analysis pipeline, we identified Protein far-red elongated hypocotyl 3 (Glyma.04G124300; E8-r1), E1-like-a (Glyma.04G156400; E8-r2), Light-harvesting chlorophyll-protein complex I subunit A4 (Glyma.04G167900; E8-r3), and Cycling dof factor 3 (Glyma.04G168300; E8-r3) as the most promising candidate genes for these regions. A combinatorial eQTL mapping approach identified significant regulatory interactions for 13 expression traits (e-traits), including Glyma.04G050200 (Early flowering 3/E6 locus), with the E8-r3 region. Four other important QTL regions close to or encompassing major flowering genes were also detected on chromosomes GM07, GM08, and GM16. In GM07:5,256,305-5,404,971, a missense polymorphism was detected in the candidate gene Glyma.07G058200 (Protein suppressor of PHYA-105). These findings demonstrate that the locus known as E8 is regulated by at least three distinct genomic regions, all of which comprise major flowering genes.

Variant Classification Discordance: Contributing Factors and Predictive Models

An ever-growing catalog of human variants is hosted in the ClinVar database. In this database, submissions on a variant are combined into a multisubmitter record; and in the case of discordance in variant classification between submitters, the record is labeled as conflicting. The current study used ClinVar data to identify characteristics that would make variants more likely to be associated with the conflict class of variants. Furthermore, the Extreme Gradient Boosting algorithm was used to train classifier models to provide prediction of classification discordance for single submission variants in ClinVar database. Population allele frequency, the gene harboring the variant, variant type, consequence on protein, variant deleteriousness score, first submitter identity, and submission count were associated with conflict in variant classification. Using such features, the optimized classifier showed accuracy on the test set of 88% with the weighted average of precision, recall, and f1-score of 0.84, 0.88, and 0.85, respectively. There were pronounced associations between variant classification discordance and allele frequency, gene type, and the identity of the first submitter. The study provides the predicted discordance status for single-submitter variants deposited in ClinVar. This approach can be used to assess whether single-submitter variants are likely to be supported, or in conflict with, future entries; this knowledge may help laboratories with clinical variant assessment.

Exploring noncoding variants in genetic diseases: from detection to functional insights

Previous studies on genetic diseases predominantly focused on protein-coding variations, overlooking the vast noncoding regions in the human genome. The development of high-throughput sequencing technologies and functional genomics tools has enabled the systematic identification of functional noncoding variants. These variants can impact gene expression, regulation, and chromatin conformation, thereby contributing to disease pathogenesis. Understanding the mechanisms that underlie the impact of noncoding variants on genetic diseases is indispensable for the development of precisely targeted therapies and the implementation of personalized medicine strategies. The intricacies of noncoding regions introduce a multitude of challenges and research opportunities. In this review, we introduce a spectrum of noncoding variants involved in genetic diseases, along with research strategies and advanced technologies for their precise identification and in-depth understanding of the complexity of the noncoding genome. We will delve into the research challenges and propose potential solutions for unraveling the genetic basis of rare and complex diseases.

Combining a prioritization strategy and functional studies nominates 5'UTR variants underlying inherited retinal disease

BACKGROUND

5' untranslated regions (5'UTRs) are essential modulators of protein translation. Predicting the impact of 5'UTR variants is challenging and rarely performed in routine diagnostics. Here, we present a combined approach of a comprehensive prioritization strategy and functional assays to evaluate 5'UTR variation in two large cohorts of patients with inherited retinal diseases (IRDs).

METHODS

We performed an isoform-level re-analysis of retinal RNA-seq data to identify the protein-coding transcripts of 378 IRD genes with highest expression in retina. We evaluated the coverage of their 5'UTRs by different whole exome sequencing (WES) kits. The selected 5'UTRs were analyzed in whole genome sequencing (WGS) and WES data from IRD sub-cohorts from the 100,000 Genomes Project (n = 2397 WGS) and an in-house database (n = 1682 WES), respectively. Identified variants were annotated for 5'UTR-relevant features and classified into seven categories based on their predicted functional consequence. We developed a variant prioritization strategy by integrating population frequency, specific criteria for each category, and family and phenotypic data. A selection of candidate variants underwent functional validation using diverse approaches.

RESULTS

Isoform-level re-quantification of retinal gene expression revealed 76 IRD genes with a non-canonical retina-enriched isoform, of which 20 display a fully distinct 5'UTR compared to that of their canonical isoform. Depending on the probe design, 3-20% of IRD genes have 5'UTRs fully captured by WES. After analyzing these regions in both cohorts, we prioritized 11 (likely) pathogenic variants in 10 genes (ARL3, MERTK, NDP, NMNAT1, NPHP4, PAX6, PRPF31, PRPF4, RDH12, RD3), of which 7 were novel. Functional analyses further supported the pathogenicity of three variants. Mis-splicing was demonstrated for the PRPF31:c.-9+1G>T variant. The MERTK:c.-125G>A variant, overlapping a transcriptional start site, was shown to significantly reduce both luciferase mRNA levels and activity. The RDH12:c.-123C>T variant was found in cis with the hypomorphic RDH12:c.701G>A (p.Arg234His) variant in 11 patients. This 5'UTR variant, predicted to introduce an upstream open reading frame, was shown to result in reduced RDH12 protein but unaltered mRNA levels.

CONCLUSIONS

This study demonstrates the importance of 5'UTR variants implicated in IRDs and provides a systematic approach for 5'UTR annotation and validation that is applicable to other inherited diseases.

A Comprehensive Allele Specific Expression Resource for the Equine Transcriptome

Background

Allele-specific expression (ASE) analysis provides a nuanced view of cis-regulatory mechanisms affecting gene expression.

Results

In this work, we introduce and highlight the significance of an equine ASE analysis, containing integrated long- and short-read RNA sequencing data, along with insight from histone modification data, from four healthy Thoroughbreds (2 mares and 2 stallions) across 9 tissues.

Conclusions

This valuable publicly accessible resource is poised to facilitate investigations into regulatory variation in equine tissues and foster a deeper understanding of the impact of allelic imbalance in equine health and disease at the molecular level.

GCLiPP: global crosslinking and protein purification method for constructing high-resolution occupancy maps for RNA binding proteins

GCLiPP is a global RNA interactome capture method that detects RNA-binding protein (RBP) occupancy transcriptome-wide. GCLiPP maps RBP-occupied sites at a higher resolution than phase separation-based techniques. GCLiPP sequence tags correspond with known RBP binding sites and are enriched for sites detected by RBP-specific crosslinking immunoprecipitation (CLIP) for abundant cytosolic RBPs. Comparison of human Jurkat T cells and mouse primary T cells uncovers shared peaks of GCLiPP signal across homologous regions of human and mouse 3' UTRs, including a conserved mRNA-destabilizing cis-regulatory element. GCLiPP signal overlapping with immune-related SNPs uncovers stabilizing cis-regulatory regions in CD5, STAT6, and IKZF1.

Gonadal transcriptome analysis of paradise fish Macropodus opercularis to reveal sex-related genes

Macropodus opercularis is an ornamental fish species endemic to China, with obvious sexual dimorphism in phenotype. To obtain the gene expression profile of the gonads of M. opercularis and explore its sex-related genes, six cDNA libraries were constructed from the sexually mature M. opercularis, and RNA-seq analysis was performed. The sequenced clean data were assembled by de novo splicing to generate 171,415 unigenes, and differentially expressed genes (DEGs) screening revealed that there were 41,638 DEGs in the gonads of M. opercularis. By comparing those DEGS in the ovary with the testis, we found 29,870 DEGs were upregulated and 11,768 DEGs were downregulated. Kyoto Encyclopedia of Genes and Genomes (KEGG) and gene ontology (GO) enrichment analysis showed that GO terms related to cell cycle and gamete formation were enriched, and pathway signals related to sex differences, such as FoxO signalling pathway and PI3K-Akt signalling pathway, were also detected. Reverse transcript fluorescence quantitative PCR (RT-qPCR) validation of 14 DEGs associated with sex differences showed that the RT-qPCR results were consistent with RNA-Seq analysis, and five genes, foxl2, sox3, foxo, zar1, cyp19a1, were significantly expressed in the ovaries. dmrt1, cyp11b, amh, sf1, sox9, gdf6, dmrt3, fstl1 and hsd11b2, a total of nine genes were significantly expressed in the testis. The results of this study provide a basis for the study of gonadal differentiation, developmental mechanisms and related functional genes in M. opercularis.

Genetic variants associated with hantavirus infection in a reservoir host are related to regulation of inflammation and immune surveillance

The immunogenetics of wildlife populations influence the epidemiology and evolutionary dynamic of the host-pathogen system. Profiling immune gene diversity present in wildlife may be especially important for those species that, while not at risk of disease or extinction themselves, are host to diseases that are a threat to humans, other wildlife, or livestock. Hantaviruses (genus: Orthohantavirus) are globally distributed zoonotic RNA viruses with pathogenic strains carried by a diverse group of rodent hosts. The marsh rice rat (Oryzomys palustris) is the reservoir host of Orthohantavirus bayoui, a hantavirus that causes fatal cases of hantavirus cardiopulmonary syndrome in humans. We performed a genome wide association study (GWAS) using the rice rat "immunome" (i.e., all exons related to the immune response) to identify genetic variants associated with infection status in wild-caught rice rats naturally infected with their endemic strain of hantavirus. First, we created an annotated reference genome using 10× Chromium Linked Reads sequencing technology. This reference genome was used to create custom baits which were then used to target enrich prepared rice rat libraries (n = 128) and isolate their immunomes prior to sequencing. Top SNPs in the association test were present in four genes (Socs5, Eprs, Mrc1, and Il1f8) which have not been previously implicated in hantavirus infections. However, these genes correspond with other loci or pathways with established importance in hantavirus susceptibility or infection tolerance in reservoir hosts: the JAK/STAT, MHC, and NFκB. These results serve as informative markers for future exploration and highlight the importance of immune pathways that repeatedly emerge across hantavirus systems. Our work aids in creating cross-species comparisons for better understanding mechanisms of genetic susceptibility and host-pathogen coevolution in hantavirus systems.

Functional characterization of the disease-associated CCL2 rs1024611G-rs13900T haplotype: The role of the RNA-binding protein HuR

CC-chemokine ligand 2 (CCL2) is involved in the pathogenesis of several diseases associated with monocyte/macrophage recruitment, such as HIV-associated neurocognitive disorder (HAND), tuberculosis, and atherosclerosis. The rs1024611 (alleles:A>G; G is the risk allele) polymorphism in the CCL2 cis-regulatory region is associated with increased CCL2 expression in vitro and ex vivo, leukocyte mobilization in vivo, and deleterious disease outcomes. However, the molecular basis for the rs1024611-associated differential CCL2 expression remains poorly characterized. It is conceivable that genetic variant(s) in linkage disequilibrium (LD) with rs1024611 could mediate such effects. Previously, we used rs13900 (alleles:_C>T) in the CCL2 3' untranslated region (3' UTR) that is in perfect LD with rs1024611 to demonstrate allelic expression imbalance (AEI) of CCL2 in heterozygous individuals. Here we tested the hypothesis that the rs13900 could modulate CCL2 expression by altering mRNA turnover and/or translatability. The rs13900 T allele conferred greater stability to the CCL2 transcript when compared to the rs13900 C allele. The rs13900 T allele also had increased binding to Human Antigen R (HuR), an RNA-binding protein, in vitro and ex vivo. The rs13900 alleles imparted differential activity to reporter vectors and influenced the translatability of the reporter transcript. We further demonstrated a role for HuR in mediating allele-specific effects on CCL2 expression in overexpression and silencing studies. The presence of the rs1024611G-rs13900T conferred a distinct transcriptomic signature related to inflammation and immunity. Our studies suggest that the differential interactions of HuR with rs13900 could modulate CCL2 expression and explain the interindividual differences in CCL2-mediated disease susceptibility.

LDLR and PCSK9 3´UTR variants and their putative effects on microRNA molecular interactions in familial hypercholesterolemia: a computational approach

BACKGROUND

Familial hypercholesterolemia (FH) is caused by pathogenic variants in low-density lipoprotein (LDL) receptor (LDLR) or its associated genes, including apolipoprotein B (APOB), proprotein convertase subtilisin/kexin type 9 (PCSK9), and LDLR adaptor protein 1 (LDLRAP1). However, approximately 40% of the FH patients clinically diagnosed (based on FH phenotypes) may not carry a causal variant in a FH-related gene. Variants located at 3' untranslated region (UTR) of FH-related genes could elucidate mechanisms involved in FH pathogenesis. This study used a computational approach to assess the effects of 3'UTR variants in FH-related genes on miRNAs molecular interactions and to explore the association of these variants with molecular diagnosis of FH.

METHODS AND RESULTS

Exons and regulatory regions of FH-related genes were sequenced in 83 FH patients using an exon-target gene sequencing strategy. In silico prediction tools were used to study the effects of 3´UTR variants on interactions between miRNAs and target mRNAs. Pathogenic variants in FH-related genes (molecular diagnosis) were detected in 44.6% FH patients. Among 59 3'UTR variants identified, LDLR rs5742911 and PCSK9 rs17111557 were associated with molecular diagnosis of FH, whereas LDLR rs7258146 and rs7254521 and LDLRAP1 rs397860393 had an opposite effect (p < 0.05). 3´UTR variants in LDLR (rs5742911, rs7258146, rs7254521) and PCSK9 (rs17111557) disrupt interactions with several miRNAs, and more stable bindings were found with LDLR (miR-4435, miR-509-3 and miR-502) and PCSK9 (miR-4796).

CONCLUSION

LDLR and PCSK9 3´UTR variants disturb miRNA:mRNA interactions that could affect gene expression and are potentially associated with molecular diagnosis of FH.

Polymorphisms in corticotrophin-releasing hormone-proopiomalanocortin (CRH-POMC) system genes: Neuroimmune contributions to psoriasis disease

BACKGROUND

Skin is a target organ and source of the corticotropin-releasing hormone-proopiomelanocortin (CRH-POMC) system, operating as a coordinator and executor of responses to stress. Environmental stress exacerbates and triggers inflammatory skin diseases through modifying the cellular components of the immune system supporting the importance of CRH-POMC system in the pathogenesis of psoriasis. The aim of this study was to analyse the association of CRH-POMC polymorphisms with psoriasis and evaluate transcript expression of lesional psoriatic and normal skin in RNA-seq data.

METHODS

Samples of 104 patients with psoriasis and 174 healthy controls were genotyped for 42 single nucleotide polymorphisms (SNPs) of CRH-POMC using Applied Biosystems SNPlex™ method. The transcript quantification was performed using Salmon software v1.3.0.

RESULTS

This study demonstrated the associations between melanocortin 1 receptor (MC1R) polymorphisms rs2228479, rs3212369, dopachrome tautomerase (DCT) polymorphisms rs7987802, rs2031526, rs9524501 and psoriasis in the Tatar population. Very strong association was evident for the SNP rs7987802 in the DCT gene (pc  = 5.95е-006) in psoriasis patients. Additionally, the haplotype analysis provided AT DCT (rs7992630 and rs7987802) and AGA MC1R (rs3212358, 2228479 and 885479) haplotypes significantly associated (pc  ˂ 0.05) with psoriasis in the Tatar population, supporting the involvement of DCT and MC1R to the psoriasis susceptibility. Moreover, MC1R-203 and DCT-201 expression levels were decreased in psoriasis lesional skin compared with healthy control skin.

CONCLUSIONS

This study is the first to identify genetic variants of the MC1R and DCT genes significantly associated with psoriasis in Tatar population. Our results support potential roles of CRH-POMC system genes and DCT in the pathogenesis of psoriasis.

A curated census of pathogenic and likely pathogenic UTR variants and evaluation of deep learning models for variant effect prediction

Introduction: Variants in 5' and 3' untranslated regions (UTR) contribute to rare disease. While predictive algorithms to assist in classifying pathogenicity can potentially be highly valuable, the utility of these tools is often unclear, as it depends on carefully selected training and validation conditions. To address this, we developed a high confidence set of pathogenic (P) and likely pathogenic (LP) variants and assessed deep learning (DL) models for predicting their molecular effects. Methods: 3' and 5' UTR variants documented as P or LP (P/LP) were obtained from ClinVar and refined by reviewing the annotated variant effect and reassessing evidence of pathogenicity following published guidelines. Prediction scores from sequence-based DL models were compared between three groups: P/LP variants acting though the mechanism for which the model was designed (model-matched), those operating through other mechanisms (model-mismatched), and putative benign variants. PhyloP was used to compare conservation scores between P/LP and putative benign variants. Results: 295 3' and 188 5' UTR variants were obtained from ClinVar, of which 26 3' and 68 5' UTR variants were classified as P/LP. Predictions by DL models achieved statistically significant differences when comparing modelmatched P/LP variants to both putative benign variants and modelmismatched P/LP variants, as well as when comparing all P/LP variants to putative benign variants. PhyloP conservation scores were significantly higher among P/LP compared to putative benign variants for both the 3' and 5' UTR. Discussion: In conclusion, we present a high-confidence set of P/LP 3' and 5' UTR variants spanning a range of mechanisms and supported by detailed pathogenicity and molecular mechanism evidence curation. Predictions from DL models further substantiate these classifications. These datasets will support further development and validation of DL algorithms designed to predict the functional impact of variants that may be implicated in rare disease.

Machine Learning Application Identifies Germline Markers of Hypertension in Patients With Ovarian Cancer Treated With Carboplatin, Taxane, and Bevacizumab

Pharmacogenomics studies how genes influence a person's response to treatment. When complex phenotypes are influenced by multiple genetic variations with little effect, a single piece of genetic information is often insufficient to explain this variability. The application of machine learning (ML) in pharmacogenomics holds great potential - namely, it can be used to unravel complicated genetic relationships that could explain response to therapy. In this study, ML techniques were used to investigate the relationship between genetic variations affecting more than 60 candidate genes and carboplatin-induced, taxane-induced, and bevacizumab-induced toxicities in 171 patients with ovarian cancer enrolled in the MITO-16A/MaNGO-OV2A trial. Single-nucleotide variation (SNV, formerly SNP) profiles were examined using ML to find and prioritize those associated with drug-induced toxicities, specifically hypertension, hematological toxicity, nonhematological toxicity, and proteinuria. The Boruta algorithm was used in cross-validation to determine the significance of SNVs in predicting toxicities. Important SNVs were then used to train eXtreme gradient boosting models. During cross-validation, the models achieved reliable performance with a Matthews correlation coefficient ranging from 0.375 to 0.410. A total of 43 SNVs critical for predicting toxicity were identified. For each toxicity, key SNVs were used to create a polygenic toxicity risk score that effectively divided individuals into high-risk and low-risk categories. In particular, compared with low-risk individuals, high-risk patients were 28-fold more likely to develop hypertension. The proposed method provided insightful data to improve precision medicine for patients with ovarian cancer, which may be useful for reducing toxicities and improving toxicity management.

SIGIRR gene variants in term newborns with congenital heart defects and necrotizing enterocolitis

Background

Necrotizing enterocolitis (NEC) is a common gastrointestinal emergency among neonates which is characterized by acute intestinal inflammation and necrosis. The main risk factors for NEC are prematurity, low birth weight, and some preexisting health conditions such as congenital heart defects (CHDs). Investigation of the potential genetic predisposition to NEC is a promising approach that might provide new insights into its pathogenesis. One of the most important proteins that play a significant role in the pathogenesis of NEC is Toll-like receptor 4 (TLR4) which recognizes lipopolysaccharide found in Gram-negative bacteria. In intestinal epithelial cells, a protein encoded by the SIGIRR gene is a major inhibitor of TLR4 signaling. A few SIGIRR variants, including rare p.Y168X and p.S80Y, have already been identified in preterm infants with NEC, but their pathogenic significance remains unclear. This study aimed to investigate the spectrum of SIGIRR genetic variants in term newborns with CHD and to assess their potential association with NEC.

Methods and Results

A total of 93 term newborns with critical CHD were enrolled in this study, 33 of them developed NEC. SIGIRR genetic variants were determined by Sanger sequencing of all exons. In total, eight SIGIRR genetic variants were identified, two of which were found only in newborns with NEC (P = 0.12). The rare missense p.S80Y (rs117739035) variant in exon 4 was found in two infants with NEC stage IIA. Two infants with NEC stage III and stage IB carried a novel duplication c. 102_121dup (rs552367848) variant in exon 10 that has not been previously associated with any clinical phenotype.

Conclusions

The presence of both variants only in neonates who developed NEC, together with earlier published data, may suggest their potential contribution to the risk of developing NEC in term infants with CHD and allow planning larger cohort studies to clarify their relevance.

Analysis of subcellular RNA fractions demonstrates significant genetic regulation of gene expression in human brain post-transcriptionally

Gaining insight into the genetic regulation of gene expression in human brain is key to the interpretation of genome-wide association studies for major neurological and neuropsychiatric diseases. Expression quantitative trait loci (eQTL) analyses have largely been used to achieve this, providing valuable insights into the genetic regulation of steady-state RNA in human brain, but not distinguishing between molecular processes regulating transcription and stability. RNA quantification within cellular fractions can disentangle these processes in cell types and tissues which are challenging to model in vitro. We investigated the underlying molecular processes driving the genetic regulation of gene expression specific to a cellular fraction using allele-specific expression (ASE). Applying ASE analysis to genomic and transcriptomic data from paired nuclear and cytoplasmic fractions of anterior prefrontal cortex, cerebellar cortex and putamen tissues from 4 post-mortem neuropathologically-confirmed control human brains, we demonstrate that a significant proportion of genetic regulation of gene expression occurs post-transcriptionally in the cytoplasm, with genes undergoing this form of regulation more likely to be synaptic. These findings have implications for understanding the structure of gene expression regulation in human brain, and importantly the interpretation of rapidly growing single-nucleus brain RNA-sequencing and eQTL datasets, where cytoplasm-specific regulatory events could be missed.

Time trees and clock genes: a systematic review and comparative analysis of contemporary avian migration genetics

Timing is a crucial aspect for survival and reproduction in seasonal environments leading to carefully scheduled annual programs of migration in many species. But what are the exact mechanisms through which birds (class: Aves) can keep track of time, anticipate seasonal changes, and adapt their behaviour? One proposed mechanism regulating annual behaviour is the circadian clock, controlled by a highly conserved set of genes, collectively called 'clock genes' which are well established in controlling the daily rhythmicity of physiology and behaviour. Due to diverse migration patterns observed within and among species, in a seemingly endogenously programmed manner, the field of migration genetics has sought and tested several candidate genes within the clock circuitry that may underlie the observed differences in breeding and migration behaviour. Among others, length polymorphisms within genes such as Clock and Adcyap1 have been hypothesised to play a putative role, although association and fitness studies in various species have yielded mixed results. To contextualise the existing body of data, here we conducted a systematic review of all published studies relating polymorphisms in clock genes to seasonality in a phylogenetically and taxonomically informed manner. This was complemented by a standardised comparative re-analysis of candidate gene polymorphisms of 76 bird species, of which 58 are migrants and 18 are residents, along with population genetics analyses for 40 species with available allele data. We tested genetic diversity estimates, used Mantel tests for spatial genetic analyses, and evaluated relationships between candidate gene allele length and population averages for geographic range (breeding- and non-breeding latitude), migration distance, timing of migration, taxonomic relationships, and divergence times. Our combined analysis provided evidence (i) of a putative association between Clock gene variation and autumn migration as well as a putative association between Adcyap1 gene variation and spring migration in migratory species; (ii) that these candidate genes are not diagnostic markers to distinguish migratory from sedentary birds; and (iii) of correlated variability in both genes with divergence time, potentially reflecting ancestrally inherited genotypes rather than contemporary changes driven by selection. These findings highlight a tentative association between these candidate genes and migration attributes as well as genetic constraints on evolutionary adaptation.

A Pilot Genome-Wide Association Study Identifies Novel Markers of Metabolic Syndrome in Patients with Psoriasis

BACKGROUND

Recent studies have reported that psoriasis is associated with the development of metabolic syndrome. Genome-wide association studies have been used to discover gene variant markers that occur frequently in case group in relation to specific diseases.

OBJECTIVE

The aim of the present study was to investigate the variants of specific genes involved in metabolic syndrome associated with psoriasis.

METHODS

A total of 95 psoriasis patients were recruited and divided into two groups: one with metabolic syndrome (38 patients) and the other without (57 patients). After genotyping, imputation, and quality checking, the association between the several single nucleotide polymorphisms and metabolic syndrome in psoriasis was tested, followed by gene set enrichment analysis.

RESULTS

We found 76 gene polymorphisms that conferred an increased risk for metabolic syndrome in patients with psoriasis. Four single nucleotide polymorphisms (rs17154774 of FRMD4A, rs77498336 of GPR116, rs75949580 and rs187682251 of MAPK4) showed the strongest association between metabolic syndrome and psoriasis. The epidermal growth factor receptor protein was located at the center of the protein interactions for the gene polymorphisms.

CONCLUSION

This study identified several previously unknown polymorphisms associated with metabolic syndrome in psoriasis. These results highlight the potential for future genetic studies to elucidate the development, and ultimately prevent the onset, of metabolic syndrome in patients with psoriasis.

Investigating the NRAS 5' UTR as a target for small molecules

Neuroblastoma RAS (NRAS) is an oncogene that is deregulated and highly mutated in cancers including melanomas and acute myeloid leukemias. The 5' untranslated region (UTR) (5' UTR) of the NRAS mRNA contains a G-quadruplex (G4) that regulates translation. Here we report a novel class of small molecule that binds to the G4 structure located in the 5' UTR of the NRAS mRNA. We used a small molecule microarray screen to identify molecules that selectively bind to the NRAS-G4 with submicromolar affinity. One compound inhibits the translation of NRAS in vitro but showed only moderate effects on the NRAS levels in cellulo. Rapid Amplification of cDNA Ends and RT-PCR analysis revealed that the predominant NRAS transcript does not possess the G4 structure. Thus, although NRAS transcripts lack a G4 in many cell lines the concept of targeting folded regions within 5' UTRs to control translation remains a highly attractive strategy.

Role of non-coding variants in cardiovascular disease

Cardiovascular diseases (CVDs) constitute one of the significant causes of death worldwide. Different pathological states are linked to CVDs, which despite interventions and treatments, still have poor prognoses. The genetic component, as a beneficial tool in the risk stratification of CVD development, plays a role in the pathogenesis of this group of diseases. The emergence of genome-wide association studies (GWAS) have led to the identification of non-coding parts associated with cardiovascular traits and disorders. Variants located in functional non-coding regions, including promoters/enhancers, introns, miRNAs and 5'/3' UTRs, account for 90% of all identified single-nucleotide polymorphisms associated with CVDs. Here, for the first time, we conducted a comprehensive review on the reported non-coding variants for different CVDs, including hypercholesterolemia, cardiomyopathies, congenital heart diseases, thoracic aortic aneurysms/dissections and coronary artery diseases. Additionally, we present the most commonly reported genes involved in each CVD. In total, 1469 non-coding variants constitute most reports on familial hypercholesterolemia, hypertrophic cardiomyopathy and dilated cardiomyopathy. The application and identification of non-coding variants are beneficial for the genetic diagnosis and better therapeutic management of CVDs.

Genetic bases of C7 deficiency: systematic review and report of a novel deletion determining functional hemizygosity

Primary complement system (C) deficiencies are rare but notably associated with an increased risk of infections, autoimmunity, or immune disorders. Patients with terminal pathway C-deficiency have a 1,000- to 10,000-fold-higher risk of Neisseria meningitidis infections and should be therefore promptly identified to minimize the likelihood of further infections and to favor vaccination. In this paper, we performed a systematic review about clinical and genetic patterns of C7 deficiency starting from the case of a ten-year old boy infected by Neisseria meningitidis B and with clinical presentation suggestive of reduced C activity. Functional assay via Wieslab ELISA Kit confirmed a reduction in total C activity of the classical (0.6% activity), lectin (0.2% activity) and alternative (0.1% activity) pathways. Western blot analysis revealed the absence of C7 in patient serum. Sanger sequencing of genomic DNA extracted from peripheral blood of the patient allowed the identification of two pathogenetic variants in the C7 gene: the already well-characterized missense mutation G379R and a novel heterozygous deletion of three nucleotides located at the 3'UTR (c.*99_*101delTCT). This mutation resulted in an instability of the mRNA; thus, only the allele containing the missense mutation was expressed, making the proband a functional hemizygote for the expression of the mutated C7 allele.

The role of noncoding genetic variants in cardiomyopathy

Cardiomyopathies remain one of the leading causes of morbidity and mortality worldwide. Environmental risk factors and genetic predisposition account for most cardiomyopathy cases. As with all complex diseases, there are significant challenges in the interpretation of the molecular mechanisms underlying cardiomyopathy-associated genetic variants. Given the technical improvements and reduced costs of DNA sequence technologies, an increasing number of patients are now undergoing genetic testing, resulting in a continuously expanding list of novel mutations. However, many patients carry noncoding genetic variants, and although emerging evidence supports their contribution to cardiac disease, their role in cardiomyopathies remains largely understudied. In this review, we summarize published studies reporting on the association of different types of noncoding variants with various types of cardiomyopathies. We focus on variants within transcriptional enhancers, promoters, intronic sites, and untranslated regions that are likely associated with cardiac disease. Given the broad nature of this topic, we provide an overview of studies that are relatively recent and have sufficient evidence to support a significant degree of causality. We believe that more research with additional validation of noncoding genetic variants will provide further mechanistic insights on the development of cardiac disease, and noncoding variants will be increasingly incorporated in future genetic screening tests.

The non-coding genome in Autism Spectrum Disorders

Autism Spectrum Disorders (ASD) are a group of neurodevelopmental disorders (NDDs) characterized by difficulties in social interaction and communication, repetitive behavior, and restricted interests. While ASD have been proven to have a strong genetic component, current research largely focuses on coding regions of the genome. However, non-coding DNA, which makes up for ∼99% of the human genome, has recently been recognized as an important contributor to the high heritability of ASD, and novel sequencing technologies have been a milestone in opening up new directions for the study of the gene regulatory networks embedded within the non-coding regions. Here, we summarize current progress on the contribution of non-coding alterations to the pathogenesis of ASD and provide an overview of existing methods allowing for the study of their functional relevance, discussing potential ways of unraveling ASD's "missing heritability".

Processing genome-wide association studies within a repository of heterogeneous genomic datasets

BACKGROUND

Genome Wide Association Studies (GWAS) are based on the observation of genome-wide sets of genetic variants - typically single-nucleotide polymorphisms (SNPs) - in different individuals that are associated with phenotypic traits. Research efforts have so far been directed to improving GWAS techniques rather than on making the results of GWAS interoperable with other genomic signals; this is currently hindered by the use of heterogeneous formats and uncoordinated experiment descriptions.

RESULTS

To practically facilitate integrative use, we propose to include GWAS datasets within the META-BASE repository, exploiting an integration pipeline previously studied for other genomic datasets that includes several heterogeneous data types in the same format, queryable from the same systems. We represent GWAS SNPs and metadata by means of the Genomic Data Model and include metadata within a relational representation by extending the Genomic Conceptual Model with a dedicated view. To further reduce the gap with the descriptions of other signals in the repository of genomic datasets, we perform a semantic annotation of phenotypic traits. Our pipeline is demonstrated using two important data sources, initially organized according to different data models: the NHGRI-EBI GWAS Catalog and FinnGen (University of Helsinki). The integration effort finally allows us to use these datasets within multi-sample processing queries that respond to important biological questions. These are then made usable for multi-omic studies together with, e.g., somatic and reference mutation data, genomic annotations, epigenetic signals.

CONCLUSIONS

As a result of the our work on GWAS datasets, we enable 1) their interoperable use with several other homogenized and processed genomic datasets in the context of the META-BASE repository; 2) their big data processing by means of the GenoMetric Query Language and associated system. Future large-scale tertiary data analysis may extensively benefit from the addition of GWAS results to inform several different downstream analysis workflows.

3' Untranslated Region Structural Elements in CYP24A1 Are Associated With Infantile Hypercalcemia Type 1

Loss-of-function mutations in the CYP24A1 protein-coding region causing reduced 25 hydroxyvitamin D (25OHD) and 1,25 dihydroxyvitamin D (1,25(OH)₂ D) catabolism have been observed in some cases of infantile hypercalcemia type 1 (HCINF1), which can manifest as nephrocalcinosis, hypercalcemia and adult-onset hypercalciuria, and renal stone formation. Some cases present with apparent CYP24A1 phenotypes but do not exhibit pathogenic mutations. Here, we assessed the molecular mechanisms driving apparent HCINF1 where there was a lack of CYP24A1 mutation. We obtained blood samples from 47 patients with either a single abnormality of no obvious cause or a combination of hypercalcemia, hypercalciuria, and nephrolithiasis as part of our metabolic and stone clinics. We used liquid chromatography tandem mass spectrometry (LC-MS/MS) to determine serum vitamin D metabolites and direct sequencing to confirm CYP24A1 genotype. Six patients presented with profiles characteristic of altered CYP24A1 function but lacked protein-coding mutations in CYP24A1. Analysis upstream and downstream of the coding sequence showed single nucleotide variants (SNVs) in the CYP24A1 3' untranslated region (UTR). Bioinformatics approaches revealed that these 3' UTR abnormalities did not result in microRNA silencing but altered the CYP24A1 messenger RNA (mRNA) secondary structure, which negatively impacted translation. Our experiments showed that mRNA misfolding driven by these 3' UTR sequence-dependent structural elements was associated with normal 25OHD but abnormal 1,25(OH)₂ D catabolism. Using CRISPR-Cas9 gene editing, we developed an in vitro mutant model for future CYP24A1 studies. Our results form a basis for future studies investigating structure-function relationships and novel CYP24A1 mutations producing a semifunctional protein. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Recent Trends in Antisense Therapies for Duchenne Muscular Dystrophy

Duchenne muscular dystrophy (DMD) is a debilitating and fatal genetic disease affecting 1/5000 boys globally, characterized by progressive muscle breakdown and eventual death, with an average lifespan in the mid–late twenties. While no cure yet exists for DMD, gene and antisense therapies have been heavily explored in recent years to better treat this disease. Four antisense therapies have received conditional FDA approval, and many more exist in varying stages of clinical trials. These upcoming therapies often utilize novel drug chemistries to address limitations of existing therapies, and their development could herald the next generation of antisense therapy. This review article aims to summarize the current state of development for antisense-based therapies for the treatment of Duchenne muscular dystrophy, exploring candidates designed for both exon skipping and gene knockdown.

Alcohol reverses the effects of KCNJ6 (GIRK2) noncoding variants on excitability of human glutamatergic neurons

Synonymous and noncoding single nucleotide polymorphisms (SNPs) in the KCNJ6 gene, encoding G protein-gated inwardly rectifying potassium channel subunit 2 (GIRK2), have been linked with increased electroencephalographic frontal theta event-related oscillations (ERO) in subjects diagnosed with alcohol use disorder (AUD). To identify molecular and cellular mechanisms while retaining the appropriate genetic background, we generated induced excitatory glutamatergic neurons (iN) from iPSCs derived from four AUD-diagnosed subjects with KCNJ6 variants ("Affected: AF") and four control subjects without variants ("Unaffected: UN"). Neurons were analyzed for changes in gene expression, morphology, excitability and physiological properties. Single-cell RNA sequencing suggests that KCNJ6 AF variant neurons have altered patterns of synaptic transmission and cell projection morphogenesis. Results confirm that AF neurons express lower levels of GIRK2, have greater neurite area, and elevated excitability. Interestingly, exposure to intoxicating concentrations of ethanol induces GIRK2 expression and reverses functional effects in AF neurons. Ectopic overexpression of GIRK2 alone mimics the effect of ethanol to normalize induced excitability. We conclude that KCNJ6 variants decrease GIRK2 expression and increase excitability and that this effect can be minimized or reduced with ethanol.

UTRdb 2.0: a comprehensive, expert curated catalog of eukaryotic mRNAs untranslated regions

The 5' and 3' untranslated regions of eukaryotic mRNAs (UTRs) play crucial roles in the post-transcriptional regulation of gene expression through the modulation of nucleo-cytoplasmic mRNA transport, translation efficiency, subcellular localization, and message stability. Since 1996, we have developed and maintained UTRdb, a specialized database of UTR sequences. Here we present UTRdb 2.0, a major update of UTRdb featuring an extensive collection of eukaryotic 5' and 3' UTR sequences, including over 26 million entries from over 6 million genes and 573 species, enriched with a curated set of functional annotations. Annotations include CAGE tags and polyA signals to label the completeness of 5' and 3'UTRs, respectively. In addition, uORFs and IRES are annotated in 5'UTRs as well as experimentally validated miRNA targets in 3'UTRs. Further annotations include evolutionarily conserved blocks, Rfam motifs, ADAR-mediated RNA editing events, and m6A modifications. A web interface allowing a flexible selection and retrieval of specific subsets of UTRs, selected according to a combination of criteria, has been implemented which also provides comprehensive download facilities. UTRdb 2.0 is accessible at http://utrdb.cloud.ba.infn.it/utrdb/.

A haplotype in the dipeptidyl peptidase 4 gene impacts glycemic-related traits of Brazilian older adults

Dipeptidyl peptidase 4 (DPP4) regulates various physiological pathways and has a pivotal role in glucose homeostasis. The objective of this study was to verify the association of a haplotype constituted by two single nucleotide polymorphisms (rs2268894 and rs6741949) in the DPP4 gene with type 2 diabetes mellitus (T2DM) and fasting glycemia-related variables in a sample of Brazilian older adults, taking serum levels and enzymatic activity of DPP4 into account. Clinical, biochemical, and anthropometric characteristics as well as DPP4 serum levels and enzymatic activity were determined in 800 elderly (≥60 years old) individuals. Assessment of polymorphic sites was performed by real-time PCR whereas haplotypes were inferred from genotypic frequencies. Statistical analyses compared measures and proportions according to T2DM diagnosis and DPP4 haplotypic groups. The most common haplotype consisted of the T-rs2268894/G-rs6741949 string, which was 20% more frequent among non-diabetics. Considering non-diabetic patients alone, carriers of the T/G haplotype had significantly lower levels of blood glucose, insulin, HOMA-IR index, and DPP4 activity. Among diabetic patients, the T/G haplotype was associated with lower DPP4 levels whereas glycemic scores were not affected by allelic variants. Our results suggested that the genetic architecture of DPP4 affects the glycemic profile and DPP4 serum levels and activity among elderly individuals according to the presence or absence of T2DM, with a possible implication of the T/G haplotype to the risk of T2DM onset.

PLK2 Single Nucleotide Variant in Gastric Cancer Patients Affects miR-23b-5p Binding

PURPOSE

Chromosomal instability is a hallmark of gastric cancer (GC). It can be driven by single nucleotide variants (SNVs) in cell cycle genes. We investigated the associations between SNVs in candidate genes, PLK2, PLK3, and ATM, and GC risk and clinicopathological features.

MATERIALS AND METHODS

The genotyping study included 542 patients with GC and healthy controls. Generalized linear models were used for the risk and clinicopathological association analyses. Survival analysis was performed using the Kaplan-Meier method. The binding of candidate miRs was analyzed using a luciferase reporter assay.

RESULTS

The PLK2 Crs15009-Crs963615 haplotype was under-represented in the GC group compared to that in the control group (Pcorr=0.050). Male patients with the PLK2 rs963615 CT genotype had a lower risk of GC, whereas female patients had a higher risk (P=0.023; P=0.026). The PLK2 rs963615 CT genotype was associated with the absence of vascular invasion (P=0.012). The PLK3 rs12404160 AA genotype was associated with a higher risk of GC in the male population (P=0.015). The ATM Trs228589-Ars189037-Grs4585 haplotype was associated with a higher risk of GC (P<0.001). The ATM rs228589, rs189037, and rs4585 genotypes TA+AA, AG+GG, and TG+GG were associated with the absence of perineural invasion (P=0.034). In vitro analysis showed that the cancer-associated miR-23b-5p mimic specifically bound to the PLK2 rs15009 G allele (P=0.0097). Moreover, low miR-23b expression predicted longer 10-year survival (P=0.0066) in patients with GC.

CONCLUSIONS

PLK2, PLK3, and ATM SNVs could potentially be helpful for the prediction of GC risk and clinicopathological features. PLK2 rs15009 affects the binding of miR-23b-5p. MiR-23b-5p expression status could serve as a prognostic marker for survival in patients with GC.

Lack of association between SOCS3 and SOCS7 polymorphisms and psoriasis

Background

Psoriasis is a common, chronic, inflammatory skin disease that involves changes taking place as a result of activation of the immune system. Suppressor of cytokine signaling proteins (SOCS) are intracellular proteins that act as endogenous inhibitors of proinflammatory pathways triggered by various cytokines. In this study, the relationship between psoriasis disease and SOCS gene polymorphisms is investigated in relation to the pathogenesis of psoriasis to clarify the psoriasis susceptibility profile.

Methods

The SOCS3 rs4969169 and SOCS7 rs3748726 polymorphisms were detected using the polymerase chain reaction–restriction fragment length polymorphism (PCR‐RFLP) method. The study was approved by the Clinical Research Ethics Committee of Bulent Ecevit University and performed in accordance with the ethical standards established in the 1964 Declaration of Helsinki and later amendments. All participants were informed of the parameters of the study, and they signed consent forms before being included. Statistical analysis was performed using the SPSS 18.0 (SPSS Inc.) package program.

Results

For the SOCS3 rs4969169 genotype frequency, the CC/CT genotypes represented 67%/33% in the patient group and 73%/27% in the control group. For the SOCS7 rs3748726 genotype frequency, the TT/TC/CC genotypes made up 89%/9%/1% in the patient group and 91%/8%/1% in the control group.

Conclusion

The polymorphisms of SOCS3 rs4969169 and SOCS7 rs3748726 were found to have no effective role in the pathogenesis of psoriasis. This is the first study to investigate this topic, and further studies with larger, more ethnically diverse samples are encouraged.

Disease-associated human genetic variation through the lens of precursor and mature RNA structure

Disease-associated variants (DAVs) are commonly considered either through a genomic lens that describes variant function at the DNA level, or at the protein function level if the variant is translated. Although the genomic and proteomic effects of variation are well-characterized, genetic variants disrupting post-transcriptional regulation is another mechanism of disease that remains understudied. Specific RNA sequence motifs mediate post-transcriptional regulation both in the nucleus and cytoplasm of eukaryotic cells, often by binding to RNA-binding proteins or other RNAs. However, many DAVs map far from these motifs, which suggests deeper layers of post-transcriptional mechanistic control. Here, we consider a transcriptomic framework to outline the importance of post-transcriptional regulation as a mechanism of disease-causing single-nucleotide variation in the human genome. We first describe the composition of the human transcriptome and the importance of abundant yet overlooked components such as introns and untranslated regions (UTRs) of messenger RNAs (mRNAs). We present an analysis of Human Gene Mutation Database variants mapping to mRNAs and examine the distribution of causative disease-associated variation across the transcriptome. Although our analysis confirms the importance of post-transcriptional regulatory motifs, a majority of DAVs do not directly map to known regulatory motifs. Therefore, we review evidence that regions outside these well-characterized motifs can regulate function by RNA structure-mediated mechanisms in all four elements of an mRNA: exons, introns, 5' and 3' UTRs. To this end, we review published examples of riboSNitches, which are single-nucleotide variants that result in a change in RNA structure that is causative of the disease phenotype. In this review, we present the current state of knowledge of how DAVs act at the transcriptome level, both through altering post-transcriptional regulatory motifs and by the effects of RNA structure.

DNMT3A/miR-129-2-5p/Rac1 Is an Effector Pathway for SNHG1 to Drive Stem-Cell-like and Invasive Behaviors of Advanced Bladder Cancer Cells

The stem-cell-like behavior of cancer cells plays a central role in tumor heterogeneity and invasion and correlates closely with drug resistance and unfavorable clinical outcomes. However, the molecular underpinnings of cancer cell stemness remain incompletely defined. Here, we show that SNHG1, a long non-coding RNA that is over-expressed in ~95% of human muscle-invasive bladder cancers (MIBCs), induces stem-cell-like sphere formation and the invasion of cultured bladder cancer cells by upregulating Rho GTPase, Rac1. We further show that SNHG1 binds to DNA methylation transferase 3A protein (DNMT3A), and tethers DNMT3A to the promoter of miR-129-2, thus hyper-methylating and repressing miR-129-2-5p transcription. The reduced binding of miR-129-2 to the 3'-UTR of Rac1 mRNA leads to the stabilization of Rac1 mRNA and increased levels of Rac1 protein, which then stimulates MIBC cell sphere formation and invasion. Analysis of the Human Protein Atlas shows that a high expression of Rac1 is strongly associated with poor survival in patients with MIBC. Our data strongly suggest that the SNHG1/DNMT3A/miR-129-2-5p/Rac1 effector pathway drives stem-cell-like and invasive behaviors in MIBC, a deadly form of bladder cancer. Targeting this pathway, alone or in combination with platinum-based therapy, may reduce chemoresistance and improve longer-term outcomes in MIBC patients.

Neuromuscular disorders: finding the missing genetic diagnoses

Neuromuscular disorders (NMDs) are a wide-ranging group of diseases that seriously affect the quality of life of affected individuals. The development of next-generation sequencing revolutionized the diagnosis of NMD, enabling the discovery of hundreds of NMD genes and many more pathogenic variants. However, the diagnostic yield of genetic testing in NMD cohorts remains incomplete, indicating a large number of genetic diagnoses are not identified through current methods. Fortunately, recent advancements in sequencing technologies, analytical tools, and high-throughput functional screening provide an opportunity to circumvent current challenges. Here, we discuss reasons for missing genetic diagnoses in NMD, how emerging technologies and tools can overcome these hurdles, and examine future approaches to improving diagnostic yields in NMD.

The emerging roles of PHOSPHO1 and its regulated phospholipid homeostasis in metabolic disorders

Emerging evidence suggests that phosphoethanolamine/phosphocholine phosphatase 1 (PHOSPHO1), a specific phosphoethanolamine and phosphocholine phosphatase, is involved in energy metabolism. In this review, we describe the structure and regulation of PHOSPHO1, as well as current knowledge about the role of PHOSPHO1 and its related phospholipid metabolites in regulating energy metabolism. We also examine mechanistic evidence of PHOSPHO1- and phospholipid-mediated regulation of mitochondrial and lipid droplets functions in the context of metabolic homeostasis, which could be potentially targeted for treating metabolic disorders.

Mutational Analysis of the VPREB1 Gene of Pre-BCR Complex in a Cohort of Sporadic Pediatric Patients With B-Cell Acute Lymphoblastic Leukemia

During bone marrow B-cell development, the pre-B-cell receptor is formed by the association of the immunoglobulin heavy chain with a surrogate light chain, which is encoded by the VPREB1, and λ5 genes. It is known that pre-BCR signaling signifies a critical checkpoint at the pre-B-cell stage. Thus, failure pre-BCR signaling is proposed as a critical factor for the development of B-cell acute lymphoblastic leukemia (B-ALL). B‑ALL is the most common pediatric cancer and is one of the leading causes of death in children. Until now, several molecular analyses were performed for genomic alterations in B-ALL, but for genomic analysis of the VPREB1 gene and its rare variations, limited studies have been conducted. In this study, using polymerase chain reaction and direct sequencing of 88 pediatric patients with B-ALL, we investigated the genomic region of the VPREB1 gene to find sequence variations of this gene. Our study presented ten homozygous and heterozygous point mutations and heterozygous nucleotide deletions, in the VPREB1 gene in 36 boys and 32 girls' patients. Our Bioinformatics assay results presented that these variations may alter the RNA folding, protein structure, and therefore probable effect on the protein function. These results propose that nucleotide changes probably contribute to B-ALL pathogenesis.

Meta-analysis across Nellore cattle populations identifies common metabolic mechanisms that regulate feed efficiency-related traits

Background

Feed efficiency (FE) related traits play a key role in the economy and sustainability of beef cattle production systems. The accurate knowledge of the physiologic background for FE-related traits can help the development of more efficient selection strategies for them. Hence, multi-trait weighted GWAS (MTwGWAS) and meta-analyze were used to find genomic regions associated with average daily gain (ADG), dry matter intake (DMI), feed conversion ratio (FCR), feed efficiency (FE), and residual feed intake (RFI). The FE-related traits and genomic information belong to two breeding programs that perform the FE test at different ages: post-weaning (1,024 animals IZ population) and post-yearling (918 animals for the QLT population).

Results

The meta-analyze MTwGWAS identified 14 genomic regions (-log10(p -value) > 5) regions mapped on BTA 1, 2, 3, 4, 7, 8, 11, 14, 15, 18, 21, and 29. These regions explained a large proportion of the total genetic variance for FE-related traits across-population ranging from 20% (FCR) to 36% (DMI) in the IZ population and from 22% (RFI) to 28% (ADG) in the QLT population. Relevant candidate genes within these regions (LIPE, LPL, IGF1R, IGF1, IGFBP5, IGF2, INS, INSR, LEPR, LEPROT, POMC, NPY, AGRP, TGFB1, GHSR, JAK1, LYN, MOS, PLAG1, CHCD7, LCAT, and PLA2G15) highlighted that the physiological mechanisms related to neuropeptides and the metabolic signals controlling the body's energy balance are responsible for leading to greater feed efficiency. Integrated meta-analysis results and functional pathway enrichment analysis highlighted the major effect of biological functions linked to energy, lipid metabolism, and hormone signaling that mediates the effects of peptide signals in the hypothalamus and whole-body energy homeostasis affecting the genetic control of FE-related traits in Nellore cattle.

Conclusions

Genes and pathways associated with common signals for feed efficiency-related traits provide better knowledge about regions with biological relevance in physiological mechanisms associated with differences in energy metabolism and hypothalamus signaling. These pleiotropic regions would support the selection for feed efficiency-related traits, incorporating and pondering causal variations assigning prior weights in genomic selection approaches.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08671-w.

Transcriptome analysis of fasudil treatment in the APPswe/PSEN1dE9 transgenic (APP/PS1) mice model of Alzheimer's disease

Alzheimer's disease (AD) is the most common cause of progressive dementia. In the present study, we showed hippocampal tissue transcriptome analysis in APPswe/PSEN1dE9 (APP/PS1, AD model) mice treated with fasudil (ADF) and compared with AD mice treated with saline (ADNS) and wild type mice (WT). The competing endogenous RNA (ceRNA) network was constructed and validated the differential expression of mRNA, lncRNA, miRNA, and circRNA. Our study showed differentially expressed mRNAs (DEMs) between WT and ADNS, while enriched in cell growth and death and nervous system pathways. DEMs between ADNS-ADF were enriched in the nervous system, glycosaminoglycan biosynthesis-keratan sulfate (KS) and Quorum sensing pathways. We validated four genes with RT-PCR, whereas enrichment of Acyl-CoA Synthetase Long Chain Family Member 4 (Acsl4, ENSMUST00000112903) in Quorum sensing pathways, and BTG anti-proliferation factor 1 (Btg1, ENSMUST00000038377) in RNA degradation pathways were conducted. Expression of these two genes were higher in ADNS, but were significantly reduced in ADF. Histone H4 transcription factor (Hinfp, ENSMUST00000216508) orchestrate G1/S transition of mitotic cell cycle and co-expressed with mmu-miR-26a-2-3p-mediated ceRNA and mmu-miR-3065-5p-mediated ceRNA; Wnt family member 4 (Wnt4, ENSMUST00000045747) was enriched in mTOR, Hippo and Wnt signaling pathway. Expression of these two genes were significantly lower in ADNS, and fasudil treatment reverse it. The present studies demonstrated four genes: Acsl4, Btg1, Hinfp, Wnt4 could be potential biomarkers of AD and the targets of fasudil treatment. These results will pave a novel direction for future clinic studies for AD and fasudil treatment.