Genetic variants in mRNA untranslated regions

Heritable Genetic Variability in Ovarian Tumours: Exploring Venous Thromboembolism Susceptibility and Cancer Prognosis in a Hospital-Based Study

Venous thromboembolism (VTE) is a frequently encountered paraneoplastic syndrome in patients with ovarian cancer (OC), an inflamm-aging entity. VTE is known to exacerbate their already poor prognosis, which is partially attributed to the contribution of the haemostatic system to ovarian tumourigenesis. In the past decade, numerous single-nucleotide polymorphisms (SNPs) implicated in VTE pathways have been proposed to influence tumour susceptibility and progression. These SNPs represent potential tools to improve the prognosis accuracy of OC patients. Hence, this study explored the influence of 12 haemostasis-associated SNPs on the risk for VTE, risk of OC progression and related death among 98 OC patients. The findings revealed a 20.5 % incidence of VTE, which was associated with more rapid disease progression and shorter survival times (log-rank test, p < 0.05). PROCR rs10747514 (AA/AG vs. GG; odds ratio (OR) = 3.67, p = 0.037) and SERPINE1 rs2070682 (CC/CT vs. TT; OR = 9.28, p = 0.040) were predictors of OC-related VTE development. Regarding patients' prognosis regardless of venous thrombogenesis, RGS7 rs2502448, F3 rs1361600, FGG rs2066865, and SERPINE1 rs2070682 were the most relevant biomarkers in different patient groups. These genetic variants might constitute attractive prognostic indicators among OC patients, offering insights to refine disease management strategies. However, due to the small cohort size and the study's retrospective nature, external validation is necessary to assess the generalisation of the findings.

Multiple, Single Trait GWAS and Supervised Machine Learning Reveal the Genetic Architecture of Fraxinus excelsior Tolerance to Ash Dieback in Europe

Common ash (Fraxinus excelsior) is under intensive attack from the invasive alien pathogenic fungus Hymenoscyphus fraxineus, causing ash dieback at epidemic levels throughout Europe. Previous studies have found significant genetic variation among genotypes in ash dieback susceptibility and that host phenology, such as autumn yellowing, is correlated with susceptibility of ash trees to H. fraxineus; however, the genomic basis of ash dieback tolerance in F. excelsior requires further investigation. Here, we integrate quantitative genetics based on multiple replicates and genome-wide association analyses with machine learning to reveal the genetic architecture of ash dieback tolerance and of phenological traits in F. excelsior populations in six European countries (Austria, Denmark, Germany, Ireland, Lithuania, Sweden). Based on phenotypic data of 486 F. excelsior replicated genotypes we observed negative genotypic correlations between crown damage caused by ash dieback and intensity of autumn leaf yellowing within multiple sampling sites. Our results suggest that the examined traits are polygenic and using genomic prediction models, with ranked single nucleotide polymorphisms (SNPs) based on GWAS associations as input, a large proportion of the variation was predicted by unlinked SNPs. Based on 100 unlinked SNPs, we can predict 55% of the variation in disease tolerance among genotypes (as phenotyped in genetic trials), increasing to a maximum of 63% when predicted from 9155 SNPs. In autumn leaf yellowing, 52% of variation is predicted by 100 unlinked SNPs, reaching a peak of 72% using 3740 SNPs. Based on feature permutations within genomic prediction models, a total of eight nonsynonymous SNPs linked to ash dieback crown damage and autumn leaf yellowing (three and five SNPs, respectively) were identified, these were located within genes related to plant defence (pattern triggered immunity, pathogen detection) and phenology (regulation of flowering and seed maturation, auxin transport). We did not find an overlap between genes associated with crown damage level and autumn leaf yellowing. Hence, our results shed light on the difference in the genomic basis of ADB tolerance and autumn leaf yellowing despite these two traits being correlated in quantitative genetic analysis. Overall, our methods show the applicability of genomic prediction models when combined with GWAS to reveal the genomic architecture of polygenic disease tolerance enabling the identification of ash dieback tolerant trees for breeding or conservation purposes.

Rethinking GWAS: how lessons from genetic screens and artificial intelligence could reveal biological mechanisms

MOTIVATION

Modern single-cell omics data are key to unraveling the complex mechanisms underlying risk for complex diseases revealed by genome-wide association studies (GWAS). Phenotypic screens in model organisms have several important parallels to GWAS which the author explores in this essay.

RESULTS

The author provides the historical context of such screens, comparing and contrasting similarities to association studies, and how these screens in model organisms can teach us what to look for. Then the author considers how the results of GWAS might be exhaustively interrogated to interpret the biological mechanisms underpinning disease processes. Finally, the author proposes a general framework for tackling this problem computationally, and explore the data, mechanisms, and technology (both existing and yet to be invented) that are necessary to complete the task.

AVAILABILITY AND IMPLEMENTATION

There are no data or code associated with this article.

Small-molecule RNA therapeutics to target prostate cancer

Tuning protein expression by targeting RNA structure using small molecules is an unexplored avenue for cancer treatment. To understand whether this vulnerability could be therapeutically targeted in the most lethal form of prostate cancer, castration-resistant prostate cancer (CRPC), we use a clinical small molecule, zotatifin, that targets the RNA helicase and translation factor eukaryotic initiation factor 4A (eIF4A). Zotatifin represses tumorigenesis in patient-derived and xenograft models and prolonged survival in vivo alongside hormone therapy. Genome-wide transcriptome, translatome, and proteomic analysis reveals two important translational targets: androgen receptor (AR), a key oncogene in CRPC, and hypoxia-inducible factor 1A (HIF1A), an essential cancer modulator in hypoxia. We solve the structure of the 5' UTRs of these oncogenic mRNAs and strikingly observe complex structural remodeling of these select mRNAs by this small molecule. Remarkably, tumors treated with zotatifin become more sensitive to anti-androgen therapy and radiotherapy. Therefore, "translatome therapy" provides additional strategies to treat the deadliest cancers.

Massively parallel reporter assays identify functional enhancer variants at QT interval GWAS loci

Genome-wide association studies (GWAS) have identified >30 loci with multiple common noncoding variants explaining interindividual electrocardiographic QT interval (QTi) variation. Of the many types of noncoding functional elements, here we sought to identify transcriptional enhancers with sequence variation and their cognate transcription factors (TFs) that alter the expression of proximal cardiac genes to affect QTi variation. We used massively parallel reporter assays (MPRA) in mouse cardiomyocyte HL-1 cells to screen for functional enhancer variants among 1,018 QTi-associated GWAS variants that overlap candidate cardiac enhancers across 31 loci. We identified 445 GWAS variant-containing enhancers of which 79 showed significant allelic difference in enhancer activity across 21 GWAS loci, with multiple enhancer variants per locus. Of these, we predicted differential binding by cardiac TFs, including AP-1, ATF-1, GATA2, MEF2, NKX2.5, SRF and TBX5 which are known to play key roles in development and homeostasis, at 49 enhancer variants. Finally, we used expression quantitative trait locus mapping and predicted promoter-enhancer contacts to identify 14 candidate target genes through analyses of 36 enhancer variants at 16 loci. This study provides strong evidence for 14 cardiac genes, 10 of them novel, impacting on QTi variation, beyond explaining observed genetic associations.

Enhancing the annotation of small ORF-altering variants using MORFEE: introducing MORFEEdb, a comprehensive catalog of SNVs affecting upstream ORFs in human 5'UTRs

Non-canonical small open reading frames (sORFs) are among the main regulators of gene expression. The most studied of these are upstream ORFs (upORFs) located in the 5'-untranslated region (UTR) of coding genes. Internal ORFs (intORFs) in the coding sequence and downstream ORFs (dORFs) in the 3'UTR have received less attention. Different bioinformatics tools permit the prediction of single nucleotide variants (SNVs) altering upORFs, mainly those creating AUGs or deleting stop codons, but no tool predicts variants altering non-canonical translation initiation sites and those altering intORFs or dORFs. We propose an upgrade of our MORFEE bioinformatics tool to identify SNVs that may alter all types of sORFs in coding transcripts from a VCF file. Moreover, we generate an exhaustive catalog, named MORFEEdb, reporting all possible SNVs altering existing upORFs or creating new ones in human transcripts, and provide an R script for visualizing the results. MORFEEdb has been implemented in the public platform Mobidetails. Finally, the annotation of ClinVar variants with MORFEE reveals that > 45% of UTR-SNVs can alter upORFs or dORFs. In conclusion, MORFEE and MORFEEdb have the potential to improve the molecular diagnosis of rare human diseases and to facilitate the identification of functional variants from genome-wide association studies of complex traits.

Cross-cohort analysis of expression and splicing quantitative trait loci in TOPMed

Most genetic variants associated with complex traits and diseases occur in non-coding genomic regions and are hypothesized to regulate gene expression. To understand the genetics underlying gene expression variability, we characterize 14,324 ancestrally diverse RNA-sequencing samples from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program and integrate whole genome sequencing data to perform cis and trans expression and splicing quantitative trait locus (cis-/trans-e/sQTL) analyses in six tissues and cell types, most notably whole blood (N=6,454) and lung (N=1,291). We show this dataset enables greater detection of secondary cis-e/sQTL signals than was achieved in previous studies, and that secondary cis-eQTL and primary trans-eQTL signal discovery is not saturated even though eGene discovery is. Most TOPMed trans-eQTL signals colocalize with cis-e/sQTL signals, suggesting many trans signals are mediated by cis signals. We fine-map European UK BioBank GWAS signals from 164 traits and colocalize the resulting 34,107 fine-mapped GWAS signals with TOPMed e/sQTL signals, finding that of 10,611 GWAS signals with a colocalization, 7,096 GWAS signals colocalize with at least one secondary e/sQTL signal. These results demonstrate that larger e/sQTL analyses will continue to uncover secondary e/sQTL signals, and that these new signals will benefit GWAS interpretation.

Multiplexed assays of human disease-relevant mutations reveal UTR dinucleotide composition as a major determinant of RNA stability

Untranslated regions (UTRs) contain crucial regulatory elements for RNA stability, translation and localization, so their integrity is indispensable for gene expression. Approximately 3.7% of genetic variants associated with diseases occur in UTRs, yet a comprehensive understanding of UTR variant functions remains limited due to inefficient experimental and computational assessment methods. To systematically evaluate the effects of UTR variants on RNA stability, we established a massively parallel reporter assay on 6555 UTR variants reported in human disease databases. We examined the RNA degradation patterns mediated by the UTR library in two cell lines, and then applied LASSO regression to model the influential regulators of RNA stability. We found that UA dinucleotides and UA-rich motifs are the most prominent destabilizing element. Gain of UA dinucleotide outlined mutant UTRs with reduced stability. Studies on endogenous transcripts indicate that high UA-dinucleotide ratios in UTRs promote RNA degradation. Conversely, elevated GC content and protein binding on UA dinucleotides protect high-UA RNA from degradation. Further analysis reveals polarized roles of UA-dinucleotide-binding proteins in RNA protection and degradation. Furthermore, the UA-dinucleotide ratio of both UTRs is a common characteristic of genes in innate immune response pathways, implying a coordinated stability regulation through UTRs at the transcriptomic level. We also demonstrate that stability-altering UTRs are associated with changes in biobank-based health indices, underscoring the importance of precise UTR regulation for wellness. Our study highlights the importance of RNA stability regulation through UTR primary sequences, paving the way for further exploration of their implications in gene networks and precision medicine.

DNA methylation as a contributor to dysregulation of STX6 and other frontotemporal lobar degeneration genetic risk-associated loci

Frontotemporal Lobar Degeneration (FTLD) represents a spectrum of clinically, genetically, and pathologically heterogeneous neurodegenerative disorders characterised by progressive atrophy of the frontal and temporal lobes of the brain. The two major FTLD pathological subgroups are FTLD-TDP and FTLD-tau. While the majority of FTLD cases are sporadic, heterogeneity also exists within the familial cases, typically involving mutations in MAPT, GRN or C9orf72, which is not fully explained by known genetic mechanisms. We sought to address this gap by investigating the effect of epigenetic modifications, specifically DNA methylation variation, on genes associated with FTLD genetic risk in different FTLD subtypes. We compiled a list of genes associated with genetic risk of FTLD using text-mining databases and literature searches. Frontal cortex DNA methylation profiles were derived from three FTLD datasets containing different subgroups of FTLD-TDP and FTLD-tau: FTLD1m (N = 23) containing FTLD-TDP type A C9orf72 mutation carriers and TDP Type C sporadic cases, FTLD2m (N = 48) containing FTLD-Tau MAPT mutation carriers, FTLD-TDP Type A GRN mutation carriers, and FTLD-TDP Type B C9orf72 mutation carriers and FTLD3m (N = 163) progressive supranuclear palsy (PSP) cases, and corresponding controls. To investigate the downstream effects of DNA methylation further, we then leveraged transcriptomic and proteomic datasets for FTLD cases and controls to examine gene and protein expression levels. Our analysis revealed shared promoter region hypomethylation in STX6 across FTLD-TDP and FTLD-tau subtypes, though the largest effect size was observed in the PSP cases compared to controls (delta-beta = -32%, adjusted-p value=0.002). We also observed dysregulation of the STX6 gene and protein expression across FTLD subtypes. Additionally, we performed a detailed examination of MAPT, GRN and C9orf72 in subtypes with and without the presence of the genetic mutations and observed nominally significant differentially methylated CpGs in variable positions across the genes, often with unique patterns and downstream consequences in gene/protein expression in mutation carriers. We highlight the contribution of DNA methylation at different gene regions in regulating the expression of genes previously associated with genetic risk of FTLD, including STX6. We analysed the relationship of subtypes and presence of mutations with this epigenetic mechanism to increase our understanding of how these mechanisms interact in FTLD.

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 located 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.

mRNA Transcript Variants Expressed in Mammalian Cells

Gene expression or gene regulation studies often assume one gene expresses one mRNA. However, contrary to the conventional idea, a single gene in mammalian cells can express multiple transcript variants translated into several different proteins. The transcript variants are generated through transcription from alternative start sites and alternative post-transcriptional processing of the precursor mRNA (pre-mRNA). In addition, gene mutations and RNA editing further enhance the diversity of the transcript variants. The transcript variants can encode proteins with various domains, expanding the functional repertoire of a single gene. Some transcript variants may not encode proteins but function as non-coding RNAs and regulate gene expression. The expression level of the transcript variants may vary between cell types or within the same cells under different biological conditions. Transcript variants are characteristic of cell differentiation in a particular tissue, and the variants may play a key role in normal development and aging. Studies also reported that some transcript variants may have roles in disease pathogenesis. The biological significances urge studying the complexity of gene expression at the transcript level. This article updates the molecular basis of transcript variants in mammalian cells, including the formation mechanisms and potential roles in host biology. Gaining insight into the transcript variants will not only identify novel mechanisms of gene regulation but also unravel the role of the variants in health and disease.

Association of selected gene variants with nonsyndromic orofacial clefts in Kuwait

INTRODUCTION AND OBJECTIVES

Non-syndromic orofacial clefts (NSOFCs) are complex congenital abnormalities involving both environmental and genetic factors involved in orofacial development. This study aimed to investigate the genetic association of specific genetic variants at different CYRIA gene loci with the development of NSOFCs in Kuwait.

METHODS

Four genetic variants (rs7552, rs3758249, rs3821949, and rs3917201) at four selected gene loci (CYRIA, FOXE1, MSX1, and TGFB3) were genotyped in a total of 240 DNA samples (patients (n = 114) and random controls (n = 126)) employing TaqMan® allele discrimination assay. For each variant and its genotype, the frequencies were determined and tested for Hardy-Weinberg Equilibrium. Genotype frequencies was compared between patients and controls using Pearson's test. Logistic regression analyses were employed to test for the associations of the four selected variants with the occurrence of NSOFCSs.

RESULTS

Significant differences in the distribution of genotypes between cases and controls, rs7552, rs3821949, and rs3917201 were found to have a positive association with NSOFCs. After adjusting for gender, the GG genotype of the rs7552 variant, the AG genotype of the rs3821949 variant, and the CC genotype of the rs3917201 variant showed nearly a two-fold increased risk of NSOFC (p < 0.05).

CONCLUSION

This study reports significant findings on the contribution and modest effect of CYRIA rs7552, MSX1 rs3821949, and TGFB3 rs3917201 in the development of NSOFCs. Our findings provide further evidence on the molecular mechanism and the role of the selected genes in NSOFCs.

Genetic Variants at PRKCG Splice and UTR Sites Promote Cancer Susceptibility by Disrupting Epigenetic and miRNA Regulatory Network

The changes in the protein kinase C gamma gene (PRKCG) expression are associated with both coding and non-coding variants. No studies have specifically established the association between PRKCG 3'UTR, 5'UTR, donor and acceptor splice variants with post-transcriptional changes through utilizing in-silico tools. The current study intends to uncover this linkage. In total, 419 3' and 5'UTR variants were retrieved. 325 of these variant IDs were annotated as functionally significant. 18 variants impacted the transcription factors binding and therefore influenced the post-transcriptional regulatory activity while 7 variants affected regulatory mechanisms through histone modifications. 2 rsIDs (rs373228, rs446795) potentially impacted the interactions with RNA binding proteins. In addition to that, PRKCG showed high expression in brain cells and had variable expression in TCGA tumors, respectively. Furthermore, 5 3' UTR variants were identified to be targeted by miRNAs. In total, 5 of these miRNAs (hsa-miR-663a, hsa-miR-324-5p, hsa-miR-646, hsa-miR-1205 and hsa-miR-4270) that targeted 3'UTRs (rs57483118, rs181418157 and rs60891969) showed differential expressions in distinct cancer types. The presence of 3'UTR variants likely altered the secondary structure of mRNA. The 7 rsIDs at 3' UTR site caused the loss of function of authentic splice site at 10 positions was noted; at 1 position, gain of function was observed while at 2 positions no effect was identified. Moreover, the loss of donor and acceptor splice site was evident. Our results highlight the importance of non-coding regions that might boost our research capacity to predict and construct targeted therapeutic approaches.

Alu insertion-mediated dsRNA structure formation with pre-existing Alu elements as a disease-causing mechanism

We previously identified a homozygous Alu insertion variant (Alu_Ins) in the 3'-untranslated region (3'-UTR) of SPINK1 as the cause of severe infantile isolated exocrine pancreatic insufficiency. Although we established that Alu_Ins leads to the complete loss of SPINK1 mRNA expression, the precise mechanisms remained elusive. Here, we aimed to elucidate these mechanisms through a hypothesis-driven approach. Initially, we speculated that, owing to its particular location, Alu_Ins could independently disrupt mRNA 3' end formation and/or affect other post-transcriptional processes such as nuclear export and translation. However, employing a 3'-UTR luciferase reporter assay, Alu_Ins was found to result in only an ∼50% reduction in luciferase activity compared to wild type, which is insufficient to account for the severe pancreatic deficiency in the Alu_Ins homozygote. We then postulated that double-stranded RNA (dsRNA) structures formed between Alu elements, an upstream mechanism regulating gene expression, might be responsible. Using RepeatMasker, we identified two Alu elements within SPINK1's third intron, both oriented oppositely to Alu_Ins. Through RNAfold predictions and full-length gene expression assays, we investigated orientation-dependent interactions between these Alu repeats. We provide compelling evidence to link the detrimental effect of Alu_Ins to extensive dsRNA structures formed between Alu_Ins and pre-existing intronic Alu sequences, including the restoration of SPINK1 mRNA expression by aligning all three Alu elements in the same orientation. Given the widespread presence of Alu elements in the human genome and the potential for new Alu insertions at almost any locus, our findings have important implications for detecting and interpreting Alu insertions in disease genes.

"Unveiling the genetic symphony: Diversity and expression of chicken IFITM genes in Aseel and Kadaknath breeds"

In this investigation, single nucleotide variants (SNVs) within the chicken interferon-inducible transmembrane protein (chIFITM) genes were explored in Aseel and Kadaknath breeds. Comparative analysis with the GRCg6a reference genome revealed 9 and 16 SNVs in the chIFITM locus for Aseel and Kadaknath breeds, respectively. When referencing the Genome Reference Consortium GRCg7b, Kadaknath exhibited 10 variants, contrasting with none in Aseel. Notably, 17, 8, 2, and 5 SNVs were identified in chIFITM1, chIFITM2, chIFITM3, and chIFITM5 genes, with chIFITM1 showing the highest polymorphism in Kadaknath, featuring 10 intronic variants, including three SNVs (rs16457112, rs16457111, and rs313341707) common to both breeds. Two synonymous exonic variants (g.1817767C > A and g.1819102C > T) were also noted in chIFITM1. Although chIFITM protein sequences were generally conserved, genetic variations clustered predominantly in UTR and intronic regions. Examination of immune response dynamics in live embryos uncovered notable variations in chIFITM gene expression across diverse organs and chicken breeds. Specifically, chIFITM1 mRNA was abundant in cecal tonsils for both breeds and bursa of Aseel (7.61 folds), but it was absent in the heart and lung tissues of both breeds. Conversely, chIFITM3 consistently exhibited heightened expression, particularly in bursa of Aseel (10.23 folds). Whereas mRNA of the chIFITM2 gene was found to be abundant in the heart of Kadaknath (11.03 folds) and lung of both breeds. Furthermore, the expression pattern of chIFITM5 diverged between the two breeds, the heart of Kadaknath chickens showed highest (10.45 folds). The study discovered that breed-specific genetic variants within these genes present a potential pathway for selection and breeding to improve disease resistance in chicken. The observed genetic variation among chicken populations highlights the critical importance of these variants in reinforcing virus resistance, exhibiting applicability across a wide range of breeds.

An exploratory analysis of associations of genetic variation with the efficacy of tocilizumab in severe COVID-19 patients. A pharmacogenetic study based on next-generation sequencing

Background

In the context of the cytokine storm the takes place in severe COVID-19 patients, the Interleukin 6 (IL6) pathway emerges as one of the key pathways involved in the pathogenesis of this hyperinflammatory state. The strategy of blocking the inflammatory storm by targeting the IL6 is a promising therapy to mitigate mortality. The use of Tocilizumab was recommended by the World Health Organization (WHO) to treat severe COVID-19 patients. However, the efficacy of Tocilizumab is variable. We hypothesize that the genetic background could be behind the efficacy of Tocilizumab in terms of mortality.

Methods

We performed a targeted-next generation sequencing of 287 genes, of which 264 belong to a community panel of ThermoFisher for the study of genetic causes of primary immunodeficiency disorders, and 23 additional genes mostly related to inflammation, not included in the original community panel. This panel was sequenced in an initial cohort of 425 COVID-19 patients, of which 232 were treated with Tocilizumab and standard therapy, and 193 with standard therapy only. Selected genetic variants were genotyped by single base extension in additional 245 patients (95 treated with Tocilizumab and 150 non-treated with Tocilizumab). Appropriate statistical analyses and internal validation, including logistic regression models, with the interaction between Tocilizumab and genetic variants, were applied to assess the impact of these genetic variants in the efficacy of Tocilizumab in terms of mortality.

Results

Age (p < 0.001) and cardiovascular disease (p < 0.001) are risk factors for mortality in COVID-19 patients. The presence of GG and TT genotypes at IL10Rβ (rs2834167) and IL1β (rs1143633) genes significantly associates with a reduced risk of mortality in patients treated with Tocilizumab (OR = 0.111; 95%CI = 0.015-0.829; p = 0.010 and OR = 0.378; 95%CI = 0.154-0.924; p = 0.028 respectively). The presence of CC genotype at IL1RN (rs2234679) significantly associates with an increased risk of mortality, but only in patients not treated with Tocilizumab (OR = 3.200; 95%CI = 1.512-6.771; p = 0.002). Exhaustive internal validation using a bootstrap method (B = 500 replicates) validated the accuracy of the predictive models.

Conclusion

We developed a series of predictive models based on three genotypes in genes with a strong implication in the etiopathogenesis of COVID-19 disease capable of predicting the risk of mortality in patients treated with Tocilizumab.

Bioinformatics analysis to identify the relationship between human papillomavirus-associated cervical cancer, toll-like receptors and exomes: A genetic epidemiology study

INTRODUCTION

Genetic variants may influence Toll-like receptor (TLR) signaling in the immune response to human papillomavirus (HPV) infection and lead to cervical cancer. In this study, we investigated the pattern of TLR expression in the transcriptome of HPV-positive and HPV-negative cervical cancer samples and looked for variants potentially related to TLR gene alterations in exomes from different populations.

MATERIALS AND METHODS

A cervical tissue sample from 28 women, which was obtained from the Gene Expression Omnibus database, was used to examine TLR gene expression. Subsequently, the transcripts related to the TLRs that showed significant gene expression were queried in the Genome Aggregation Database to search for variants in more than 5,728 exomes from different ethnicities.

RESULTS

Cancer and HPV were found to be associated (p<0.0001). TLR1(p = 0.001), TLR3(p = 0.004), TLR4(221060_s_at)(p = 0.001), TLR7(p = 0.001;p = 0.047), TLR8(p = 0.002) and TLR10(p = 0.008) were negatively regulated, while TLR4(1552798_at)(p<0.0001) and TLR6(p = 0.019) were positively regulated in HPV-positive patients (p<0.05). The clinical significance of the variants was statistically significant for TLR1, TLR3, TLR6 and TLR8 in association with ethnicity. Genetic variants in different TLRs have been found in various ethnic populations. Variants of the TLR gene were of the following types: TLR1(5_prime_UTR), TLR4(start_lost), TLR8(synonymous;missense) and TLR10(3_prime_UTR). The "missense" variant was found to have a risk of its clinical significance being pathogenic in South Asian populations (OR = 56,820[95%CI:40,206,80,299]).

CONCLUSION

The results of this study suggest that the variants found in the transcriptomes of different populations may lead to impairment of the functional aspect of TLRs that show significant gene expression in cervical cancer samples caused by HPV.

Phenotypic correlates between clock genes and phenology among populations of Diederik cuckoo, Chrysococcyx caprius

The Diederik cuckoo, Chrysococcyx caprius, is a small Afrotropical bird in the family Cuculidae. It is taxonomically related to 13 other species within the genus Chrysococcyx and is migratory in sub-Saharan Africa. It has a unique breeding behaviour of being a brood parasite: Breeding pairs lay their eggs in the nests of a host species and hatchlings expel the eggs of the host species. The aim of the present study was to investigate diversity in two circadian clock genes, Clock and Adcyap1, to probe for a relationship between genetic polymorphisms and their role in circannual timing and habitat selection (phenology) in intra-African migrants. DNA extracted from blood was used for the PCR amplification and sequencing of clock genes in 30 Diederik cuckoos. Three alleles were detected for Clock with similar genotypes between individuals from the Northern and Southern breeding ranges while 10 alleles were detected for Adcyap1, having shorter alleles in the North and longer alleles in the South. Population genetic analyses, including allele frequency and zygosity analysis, showed distinctly higher frequencies for the most abundant Clock allele, containing 10 polyglutamine repeats, as well as a high degree of homozygosity. In contrast, all individuals were heterozygous for Adcyap1 and alleles from both regions showed distinct differences in abundance. Comparisons between both clock genes and phenology found several phenotypic correlations. This included evidence of a relationship between the shorter alleles and habitat selection as well as a relationship between longer alleles and timing. In both instances, evidence is provided that these effects may be sex-specific. Given that these genes drive some of the synchronicity between environments and the life cycles of birds, they provide valuable insight into the fitness of species facing global challenges including climate change, urbanisation and expanding agricultural practices.

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.