Functional regulatory variants implicate distinct transcriptional networks in dementia.

Predicting the function of noncoding variation is a major challenge in modern genetics. In this study, we used massively parallel reporter assays to screen 5706 variants identified from genome-wide association studies for both Alzheimer's disease (AD) and progressive supranuclear palsy (PSP), identifying 320 functional regulatory variants (frVars) across 27 loci, including the complex 17q21.31 region. We identified and validated multiple risk loci using CRISPR interference or excision, including complement 4 (C4A) and APOC1 in AD and PLEKHM1 and KANSL1 in PSP. Functional variants disrupt transcription factor binding sites converging on enhancers with cell type-specific activity in PSP and AD, implicating a neuronal SP1-driven regulatory network in PSP pathogenesis. These analyses suggest that noncoding genetic risk is driven by common genetic variants through their aggregate activity on specific transcriptional programs.

An Association between EMX2 Variations and Mayer-Rokitansky-Küster-Hauser Syndrome: A Case-Control Study of Chinese Women.

Mayer-Rokitansky-Küster-Hauser (MRKH) syndrome is characterized by congenital malformations of Müllerian structures, including the uterus and upper two-thirds of the vagina in women. Until now, the etiology of this disease has remained unknown. We hypothesized that EMX2 (the human homologue of Drosophila empty spiracles gene (2) might be a candidate gene for MRKH syndrome because it plays an important role in the development of the urogenital system. Through sequence analysis of EMX2 in forty patients with MRKH syndrome and one hundred and forty healthy women controls, we identified eleven variations in total. Four novel variations were only found in MRKH patients, and seven single nucleotide polymorphisms were identified in both patients and controls. In silico analyses suggested that the novel variations in the 5'UTR (untranslated region) and 3'UTR might affect transcriptional activity of the EMX2 promoter or posttranscriptional processing. In conclusion, our study suggests an association between noncoding variations in the EMX2 gene and MRKH syndrome in a Chinese Han population.

Extended ensemble simulations of a SARS-CoV-2 nsp1-5'-UTR complex.

Nonstructural protein 1 (nsp1) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a 180-residue protein that blocks translation of host mRNAs in SARS-CoV-2-infected cells. Although it is known that SARS-CoV-2's own RNA evades nsp1's host translation shutoff, the molecular mechanism underlying the evasion was poorly understood. We performed an extended ensemble molecular dynamics simulation to investigate the mechanism of the viral RNA evasion. Simulation results suggested that the stem loop structure of the SARS-CoV-2 RNA 5'-untranslated region (SL1) binds to both nsp1's N-terminal globular region and intrinsically disordered region. The consistency of the results was assessed by modeling nsp1-40S ribosome structure based on reported nsp1 experiments, including the X-ray crystallographic structure analysis, the cryo-EM electron density map, and cross-linking experiments. The SL1 binding region predicted from the simulation was open to the solvent, yet the ribosome could interact with SL1. Cluster analysis of the binding mode and detailed analysis of the binding poses suggest residues Arg124, Lys47, Arg43, and Asn126 may be involved in the SL1 recognition mechanism, consistent with the existing mutational analysis.

The 3'UTR-derived sRNA RsaG coordinates redox homeostasis and metabolism adaptation in response to glucose-6-phosphate uptake in Staphylococcus aureus.

Staphylococcus aureus RsaG is a 3'-untranslated region (3'UTR) derived sRNA from the conserved uhpT gene encoding a glucose-6-phosphate (G6P) transporter expressed in response to extracellular G6P. The transcript uhpT-RsaG undergoes degradation from 5'- to 3'-end by the action of the exoribonucleases J1/J2, which are blocked by a stable hairpin structure at the 5'-end of RsaG, leading to its accumulation. RsaG together with uhpT is induced when bacteria are internalized into host cells or in the presence of mucus-secreting cells. Using MS2-affinity purification coupled with RNA sequencing, several RNAs were identified as targets including mRNAs encoding the transcriptional factors Rex, CcpA, SarA, and the sRNA RsaI. Our data suggested that RsaG contributes to the control of redox homeostasis and adjusts metabolism to changing environmental conditions. RsaG uses different molecular mechanisms to stabilize, degrade, or repress the translation of its mRNA targets. Although RsaG is conserved only in closely related species, the uhpT 3'UTR of the ape pathogen S. simiae harbors an sRNA, whose sequence is highly different, and which does not respond to G6P levels. Our results hypothesized that the 3'UTRs from UhpT transporter encoding mRNAs could have rapidly evolved to enable adaptation to host niches.

Structural and functional analysis of somatic coding and UTR indels in breast and lung cancer genomes.

Insertions and deletions (Indels) represent one of the major variation types in the human genome and have been implicated in diseases including cancer. To study the features of somatic indels in different cancer genomes, we investigated the indels from two large samples of cancer types: invasive breast carcinoma (BRCA) and lung adenocarcinoma (LUAD). Besides mapping somatic indels in both coding and untranslated regions (UTRs) from the cancer whole exome sequences, we investigated the overlap between these indels and transcription factor binding sites (TFBSs), the key elements for regulation of gene expression that have been found in both coding and non-coding sequences. Compared to the germline indels in healthy genomes, somatic indels contain more coding indels with higher than expected frame-shift (FS) indels in cancer genomes. LUAD has a higher ratio of deletions and higher coding and FS indel rates than BRCA. More importantly, these somatic indels in cancer genomes tend to locate in sequences with important functions, which can affect the core secondary structures of proteins and have a bigger overlap with predicted TFBSs in coding regions than the germline indels. The somatic CDS indels are also enriched in highly conserved nucleotides when compared with germline CDS indels.

In Silico Analysis of Dengue Virus Serotype 2 Mutations Detected at the Intrahost Level in Patients with Different Clinical Outcomes.

Intrahost genetic diversity is thought to facilitate arbovirus adaptation to changing environments and hosts, and it may also be linked to viral pathogenesis. Intending to shed light on the viral determinants for severe dengue pathogenesis, we previously analyzed the DENV-2 intrahost genetic diversity in 68 patients clinically classified as dengue fever (n = 31), dengue with warning signs (n = 19), and severe dengue (n = 18), performing viral whole-genome deep sequencing from clinical samples with an amplicon-free approach. From it, we identified a set of 141 relevant mutations distributed throughout the viral genome that deserved further attention. Therefore, we employed molecular modeling to recreate three-dimensional models of the viral proteins and secondary RNA structures to map the mutations and assess their potential effects. Results showed that, in general lines, disruptive variants were identified primarily among dengue fever cases. In contrast, potential immune-escape variants were associated mainly with warning signs and severe cases, in line with the latter's longer intrahost evolution times. Furthermore, several mutations were located on protein-surface regions, with no associated function. They could represent sites of further investigation, as the interaction of viral and host proteins is critical for both host immunomodulation and virus hijacking of the cellular machinery. The present analysis provides new information about the implications of the intrahost genetic diversity of DENV-2, contributing to the knowledge about the viral factors possibly involved in its pathogenesis within the human host. Strengthening our results with functional studies could allow many of these variants to be considered in the design of therapeutic or prophylactic compounds and the improvement of diagnostic assays. IMPORTANCE Previous evidence showed that intrahost genetic diversity in arboviruses may be linked to viral pathogenesis and that one or a few amino acid replacements within a single protein are enough to modify a biological feature of an RNA virus. To assess dengue virus serotype 2 determinants potentially involved in pathogenesis, we previously analyzed the intrahost genetic diversity of the virus in patients with different clinical outcomes and identified a set of 141 mutations that deserved further study. Thus, through a molecular modeling approach, we showed that disruptive variants were identified primarily among cases with mild dengue fever, while potential immune-escape variants were mainly associated with cases of greater severity. We believe that some of the variants pointed out in this study were attractive enough to be potentially considered in future intelligent designs of therapeutic or prophylactic compounds or the improvement of diagnostic tools. The present analysis provides new information about DENV-2 viral factors possibly involved in its pathogenesis within the human host.

Detecting and Profiling Endogenous RNA G-Quadruplexes in the Human Transcriptome.

G-quadruplexes are the non-canonical nucleic acid structures that are preferentially formed in G-rich regions. This structure has been shown to be associated with many biological functions. Regardless of the broad efforts on DNA G-quadruplexes, we still have limited knowledge on RNA G-quadruplexes, especially in a transcriptome-wide manner. Herein, by integrating the DMS-seq and the bioinformatics pipeline, we profiled and depicted the RNA G-quadruplexes in the human transcriptome. The genes that contain RNA G-quadruplexes in their specific regions are significantly related to immune pathways and the COVID-19-related gene sets. Bioinformatics analysis reveals the potential regulatory functions of G-quadruplexes on miRNA targeting at the scale of the whole transcriptome. In addition, the G-quadruplexes are depleted in the putative, not the real, PAS-strong poly(A) sites, which may weaken the possibility of such sites being the real cleaved sites. In brief, our study provides insight into the potential function of RNA G-quadruplexes in post-transcription.

Characterization of a novel botoulivirus isolated from the phytopathogenic fungus Sclerotinia sclerotiorum.

Sclerotinia sclerotiorum ourmiavirus 17 (SsOV17) was isolated from the hypovirulent strain GF3 of Sclerotinia sclerotiorum. The genome of SsOV17 is 2,802 nt in length and contains a single long open reading frame (ORF) flanked by a short structured 5'-untranslated region (5'-UTR) (28 nt) and a long 3'-UTR (788 nt), respectively. The ORF encodes a protein with 663 amino acids and a predicted molecular mass of 75.0 kDa. A BLASTp search indicated that the protein encoded by SsOV17 is closely related to the putative RNA-dependent RNA polymerase (RdRp) of Sclerotinia sclerotiorum ourmiavirus 13 (71% identity). A multiple sequence alignment indicated that eight conserved amino acid motifs were present in the RdRp conserved region of SsOV17. Phylogenetic analysis demonstrated that SsOV17 clustered with members of the genus Botoulivirus.

RAD51-UTR haplotype genetic polymorphisms and susceptibility to breast cancer in women from Jordanian population.

BACKGROUND: Genetic predisposition to breast cancer (BC) has been extensively explored to achieve an enhanced understanding of the biology of BC. Targeting candidate genes to screen different genetic variants such as RAD51 gene that plays a critical role in DNA repair pathways including the double-strand break repair system is an important task. AIM: To study several single nucleotide polymorphisms (SNPs) within RAD51-UTR gene and to find their relationship with BC risk and prognosis among Jordanian females. MATERIALS AND METHODS: In this case-control study, DNA sequencing technique was used to screen SNPs within the untranslated region (UTR) of RAD51 in 206 cases and 185 controls and the resulting data were statistically analyzed using different types of genetic analyses. Patients' clinical and pathological features were obtained from their medical records to perform genotype-phenotype association analysis. RESULTS: Our findings show a significant association between both SNPs rs528590644, rs1801320 and BC risk (p = 0.016). We estimated the correlation between many of BC prognostic factors and BC risk, and we found an association between rs1801321 and age at first menstruation (p = 0.032) in addition to a strong correlation between age at BC diagnosis and rs1801320 (p = 0.008). CONCLUSION: RAD51-UTR polymorphisms may be involved in BC development and progression.

Synergistic Effect between 3'-Terminal Noncoding and Adjacent Coding Regions of the Influenza A Virus Hemagglutinin Segment on Template Preference.

The influenza A virus genome is comprised of eight single-stranded negative-sense viral RNA (vRNA) segments. Each of the eight vRNA segments contains segment-specific nonconserved noncoding regions (NCRs) of similar sequence and length in different influenza A virus strains. However, in the subtype-determinant segments, encoding hemagglutinin (HA) and neuraminidase (NA), the segment-specific noncoding regions are subtype specific, varying significantly in sequence and length at both the 3' and 5' termini among different subtypes. The significance of these subtype-specific noncoding regions (ssNCR) in the influenza virus replication cycle is not fully understood. In this study, we show that truncations of the 3'-end H1-subtype-specific noncoding region (H1-ssNCR) resulted in recombinant viruses with decreased HA vRNA replication and attenuated growth phenotype, although the vRNA replication was not affected in single-template RNP reconstitution assays. The attenuated viruses were unstable, and point mutations at nucleotide position 76 or 56 in the adjacent coding region of HA vRNA were found after serial passage. The mutations restored the HA vRNA replication and reversed the attenuated virus growth phenotype. We propose that the terminal noncoding and adjacent coding regions act synergistically to ensure optimal levels of HA vRNA replication in a multisegment environment. These results provide novel insights into the role of the 3'-end nonconserved noncoding regions and adjacent coding regions on template preference in multiple-segmented negative-strand RNA viruses. IMPORTANCE While most influenza A virus vRNA segments contain segment-specific nonconserved noncoding regions of similar length and sequence, these regions vary considerably both in length and sequence in the segments encoding HA and NA, the two major antigenic determinants of influenza A viruses. In this study, we investigated the function of the 3'-end H1-ssNCR and observed a synergistic effect between the 3'-end H1-ssNCR nucleotides and adjacent coding nucleotide(s) of the HA segment on template preference in a multisegment environment. The results unravel an additional level of complexity in the regulation of RNA replication in multiple-segmented negative-strand RNA viruses.

Targeted sequencing of genes associated with the mismatch repair pathway in patients with endometrial cancer.

Germline variants inactivating the mismatch repair (MMR) genes MLH1, MSH2, MSH6 and PMS2 cause Lynch syndrome that implies an increased cancer risk, where colon and endometrial cancer are the most frequent. Identification of these pathogenic variants is important to identify endometrial cancer patients with inherited increased risk of new cancers, in order to offer them lifesaving surveillance. However, several other genes are also part of the MMR pathway. It is therefore relevant to search for variants in additional genes that may be associated with cancer risk by including all known genes involved in the MMR pathway. Next-generation sequencing was used to screen 22 genes involved in the MMR pathway in constitutional DNA extracted from full blood from 199 unselected endometrial cancer patients. Bioinformatic pipelines were developed for identification and functional annotation of variants, using several different software tools and custom programs. This facilitated identification of 22 exonic, 4 UTR and 9 intronic variants that could be classified according to pathogenicity. This study has identified several germline variants in genes of the MMR pathway that potentially may be associated with an increased risk for cancer, in particular endometrial cancer, and therefore are relevant for further investigation. We have also developed bioinformatics strategies to analyse targeted sequencing data, including low quality data and genomic regions outside of the protein coding exons of the relevant genes.

Pretransplant Genetic Susceptibility: Clinical Relevance in Transplant-Associated Thrombotic Microangiopathy.

Transplant-associated thrombotic microangiopathy (TA-TMA) is a life-threatening complication of allogeneic hematopoietic cell transplantation (HCT). We hypothesized that pretransplant genetic susceptibility is evident in adult TA-TMA and further investigated the association of TMA-associated variants with clinical outcomes. We studied 40 patients with TA-TMA, donors of 18 patients and 40 control non-TMA HCT recipients, without significant differences in transplant characteristics. Genomic DNA from pretransplant peripheral blood was sequenced for TMA-associated genes. Donors presented significantly lower frequency of rare variants and variants in exonic/splicing/untranslated region (UTR) regions, compared with TA-TMA patients. Controls also showed a significantly lower frequency of rare variants in ADAMTS13, CD46, CFH, and CFI. The majority of TA-TMA patients (31/40) presented with pathogenic or likely pathogenic variants. Patients refractory to conventional treatment (62%) and patients that succumbed to transplant-related mortality (65%) were significantly enriched for variants in exonic/splicing/UTR regions. In conclusion, increased incidence of pathogenic, rare and variants in exonic/splicing/UTR regions of TA-TMA patients suggests genetic susceptibility not evident in controls or donors. Notably, variants in exonic/splicing/UTR regions were associated with poor response and survival. Therefore, pretransplant genomic screening may be useful to intensify monitoring and early intervention in patients at high risk for TA-TMA.

Molecular and clinical profile of type 2 von Willebrand disease in Iran: a thirteen-year experience.

Type 2 von Willebrand disease (VWD) is the most common congenital bleeding disorder, with variable bleeding tendency and a complex laboratory phenotype. In the current study, we report the clinical and molecular profile of a large number of Iranian patients with type 2 VWD. All exons, intron-exon boundaries, and untranslated regions were sequenced by Sanger sequencing for direct mutation detection. All identified mutations were confirmed in family members and by relevant bioinformatics studies. A total of 136 patients with type 2 VWD were diagnosed, including 42 (30.9%), 32 (23.6%), 38 (27.9%), and 24 (17.6%) patients with type 2A, type 2B, type 2M, and type 2N, respectively. Epistaxis (49%), gum bleeding (30.2%), ecchymosis (23.2%), and menorrhagia (16.3%) were the most common clinical presentations, while miscarriage (2.3%) and umbilical cord bleeding (0.8%) were the rarest. Thirty mutations were identified within the VWF gene, nine (30%) being novel, with p.Arg1379Cys (n = 20), p.Val1316Met (n = 13), p.Arg1597Trp (n = 13), p.Arg1374Cys (n = 10), p.Ser1506Leu (n = 10), and p.Arg1308Cys (n = 9) the most common. Type 2 VWD is a hemorrhagic disorder with variable bleeding tendency and a heterogeneous molecular basis in patients in Iran.

Response to IL-17A inhibitors secukinumab and ixekizumab cannot be explained by genetic variation in the protein-coding and untranslated regions of the IL-17A gene: results from a multicentre study of four European psoriasis cohorts.

BACKGROUND: Genetic predictors for treatment response could optimize allocation of biological treatment in patients with psoriasis. There is minimal knowledge about pharmacogenetics of anti-IL-17 agents. OBJECTIVES: To assess whether genetic variants in the protein-coding region or untranslated regions of the IL-17A gene are associated with response to IL-17A inhibitors in patients with psoriasis. METHODS: This was a multicenter European cohort study investigating pharmacogenetics of IL-17A inhibitors in patients with psoriasis. Patients with plaque psoriasis treated with secukinumab or ixekizumab in daily practice were included. For all participants, the protein-coding region and untranslated regions of the IL-17A gene were analysed using Sanger sequencing. Identified genetic variants were tested for association with response to secukinumab/ixekizumab, measured as ∆PASI, after 12 weeks (primary outcome) and after 24 weeks (secondary outcome). Association was tested using a linear regression model with correction for baseline PASI as a fixed covariate and for biological naivety and body mass index as additional covariates. RESULTS: In total, 134 patients treated with secukinumab or ixekizumab were included. Genotyping of the cohort identified genetic variants present in untranslated regions and intronic DNA, but not in the protein-coding region of the IL-17A gene. Five genetic variants in non-coding DNA with a known or suspected functional effect on IL-17A expression were selected for association analyses: rs2275913, rs8193037, rs3819025, rs7747909 and rs3748067. After 12 weeks, 62% of patients achieved PASI75 and 39% achieved PASI90. At week 24, PASI75 and PASI90 response rates were 72% and 62%, respectively. No associations were found between the five genetic variants and ∆PASI, PASI75 or PASI90 after 12 and 24 weeks of anti-IL-17A treatment. CONCLUSIONS: Response to IL-17A inhibitors secukinumab and ixekizumab cannot be explained by genetic variation in the protein-coding and untranslated regions of the IL-17A gene. Pharmacogenetics of IL-17A inhibitors in the treatment of psoriasis requires further exploration.

HLA-G Ins/Del polymorphism and +3142C/G SNP are not related to neuromyelitis optica spectrum disorder (NMOSD) development, disability status or anti-aquaporin 4 presence in Brazilian patients.

We analyzed the association of polymorphisms from the 3' untranslated region of the HLA-G gene in 70 neuromyelitis optica spectrum disorder (NMOSD) patients and 162 healthy controls. No associations were found between the polymorphisms in NMOSD when compared to healthy controls, serology of the anti-AQP4 NMOSD biomarker and Expanded Disability Status Scale (EDSS). In conclusion, the 3' untranslated region 14 bp Ins/Del and +3142C/G polymorphisms seem not to be associated with NMOSD susceptibility, autoantibody production, nor a neurological deficit in patients.

The role of noncoding mutations in blood cancers.

The search for oncogenic mutations in haematological malignancies has largely focused on coding sequence variants. These variants have been critical in understanding these complex cancers in greater detail, ultimately leading to better disease monitoring, subtyping and prognostication. In contrast, the search for oncogenic variants in the noncoding genome has proven to be challenging given the vastness of the search space, the intrinsic difficulty in assessing the impact of variants that do not code for functional proteins, and our still primitive understanding of the function harboured by large parts of the noncoding genome. Recent studies have broken ground on this quest, identifying somatically acquired and recurrent mutations in the noncoding genome that activate the expression of proto-oncogenes. In this Review, we explore some of the best-characterised examples of noncoding mutations in haematological malignancies, and highlight how a significant majority of these variants impinge on gene regulation through the formation of aberrant enhancers and promoters. We delve into the challenges faced by those that embark on a search for noncoding driver mutations, and provide a framework distilled from studies that have successfully identified such variants to overcome some of the most salient hurdles. Finally, we discuss the current therapeutic strategies being explored to target the oncogenic mechanism supported by recurrent noncoding variants. We postulate that the continued discovery and functional characterisation of somatic variants in the noncoding genome will not only advance our understanding of haematological malignancies, but offer novel therapeutic avenues and provide important insights into transcriptional regulation on a broader scale.

Flavivirus RNA-Dependent RNA Polymerase Interacts with Genome UTRs and Viral Proteins to Facilitate Flavivirus RNA Replication.

Flaviviruses, most of which are emerging and re-emerging human pathogens and significant public health concerns worldwide, are positive-sense RNA viruses. Flavivirus replication occurs on the ER and is regulated by many mechanisms and factors. NS5, which consists of a C-terminal RdRp domain and an N-terminal methyltransferase domain, plays a pivotal role in genome replication and capping. The C-terminal RdRp domain acts as the polymerase for RNA synthesis and cooperates with diverse viral proteins to facilitate productive RNA proliferation within the replication complex. Here, we provide an overview of the current knowledge of the functions and characteristics of the RdRp, including the subcellular localization of NS5, as well as the network of interactions formed between the RdRp and genome UTRs, NS3, and the methyltransferase domain. We posit that a detailed understanding of RdRp functions may provide a target for antiviral drug discovery and therapeutics.

HLA alleles and haplotypes observed in 263 US families.

The 17th International HLA and Immunogenetics Workshop (IHIW) conducted a project entitled "The Study of Haplotypes in Families by NGS HLA". We investigated the HLA haplotypes of 1017 subjects in 263 nuclear families sourced from five US clinical immunogenetics laboratories, primarily as part of the evaluation of related donor candidates for hematopoietic stem cell and solid organ transplantation. The parents in these families belonged to five broad groups - African (72 parents), Asian (115), European (210), Hispanic (118) and "Other" (11). High-resolution HLA genotypes were generated for each subject using next-generation sequencing (NGS) HLA typing systems. We identified the HLA haplotypes in each family using HaplObserve, software that builds haplotypes in families by reviewing HLA allele segregation from parents to children. We calculated haplotype frequencies within each broad group, by treating the parents in each family as unrelated individuals. We also calculated standard measures of global linkage disequilibrium (LD) and conditional asymmetric LD for each ethnic group, and used untruncated and two-field allele names to investigate LD patterns. Finally we demonstrated the utility of consensus DNA sequences in identifying novel variants, confirming them using HLA allele segregation at the DNA sequence level.

Dynamic UTR Usage Regulates Alternative Translation to Modulate Gap Junction Formation during Stress and Aging.

Connexin43 (Cx43; gene name GJA1) is the most ubiquitously expressed gap junction protein, and understanding of its regulation largely falls under transcription and post-translational modification. In addition to Cx43, Gja1 mRNA encodes internally translated isoforms regulating gap junction formation, whose expression is modulated by TGF-ß. Here, using RLM-RACE, we identify distinct Gja1 transcripts differing only in 5' UTR length, of which two are upregulated during TGF-ß exposure and hypoxia. Introduction of these transcripts into Gja1-/- cells phenocopies the response of Gja1 to TGF-ß with reduced internal translation initiation. Inhibiting pathways downstream of TGF-ß selectively regulates levels of Gja1 transcript isoforms and translation products. Reporter assays reveal enhanced translation of full-length Cx43 from shorter Gja1 5' UTR isoforms. We also observe a correlation among UTR selection, translation, and reduced gap junction formation in aged heart tissue. These data elucidate a relationship between transcript isoform expression and translation initiation regulating intercellular communication.