Genetic variants in mRNA untranslated regions

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.

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.