Coding Region Research Articles

Knockdown of PHOX2B in the retrotrapezoid nucleus reduces the central CO2 chemoreflex in rats

PHOX2B is a transcription factor essential for the development of different classes of neurons in the central and peripheral nervous system. Heterozygous mutations in the PHOX2B coding region are responsible for the occurrence of Congenital Central Hypoventilation Syndrome (CCHS), a rare neurological disorder characterised by inadequate chemosensitivity and life-threatening sleep-related hypoventilation. Animal studies suggest that chemoreflex defects are caused in part by the improper development or function of PHOX2B expressing neurons in the retrotrapezoid nucleus (RTN), a central hub for CO2 chemosensitivity. Although the function of PHOX2B in rodents during development is well established, its role in the adult respiratory network remains unknown. In this study, we investigated whether reduction in PHOX2B expression in chemosensitive neuromedin-B (NMB) expressing neurons in the RTN altered respiratory function. Four weeks following local RTN injection of a lentiviral vector expressing the short hairpin RNA (shRNA) targeting Phox2b mRNA, a reduction of PHOX2B expression was observed in Nmb neurons compared to both naive rats and rats injected with the non-target shRNA. PHOX2B knockdown did not affect breathing in room air or under hypoxia, but ventilation was significantly impaired during hypercapnia. PHOX2B knockdown did not alter Nmb expression but it was associated with reduced expression of both Task2 and Gpr4, two CO2/pH sensors in the RTN. We conclude that PHOX2B in the adult brain has an important role in CO2 chemoreception and reduced PHOX2B expression in CCHS beyond the developmental period may contribute to the impaired central chemoreflex function.

Acquisition of different transcriptional shear mRNA and biological function of porcine interleukin 18 binding protein in PRRSV infection.

Interleukin-18 binding protein (IL-18BP), a natural regulator molecule of the pro-inflammatory cytokine interleukin-18 (IL-18), plays an important role in regulating the expression of the cellular immunity factor interferon-γ (IFN-γ). In a previous RNA-seq analysis of porcine alveolar macrophages (PAM) infected with the TIM and TJ strains of porcine reproductive and respiratory syndrome virus (PRRSV), we unexpectedly found that the mRNA expression of porcine interleukin 18-binding protein (pIL-18BP) in PAM cells infected with the TJM strain was significantly higher than that infected with the TJ strain. Studies have shown that human interleukin-18 binding protein (hIL-18bp) plays an important role in regulating cellular immunity in the course of the disease. However, there is a research gap on pIL-18BP. At the same time, PRRSV infection in pigs triggers weak cellular immune response problems. To explore the expression and the role of pIL-18BP in the cellular immune response induced by PRRSV, we strived to acquire the pIL-18BP gene from PAM or peripheral blood mononuclear cell (PBMC) with RT-PCR and sequencing. Furthermore, pIL-18BP and pIL-18 were both expressed prokaryotically and eukaryotically. The colocalization and interaction based on recombinant pIL-18BP and pIL-18 on cells were confirmed in vitro. Finally, the expression of pIL-18BP, pIL-18, and pIFN-γ was explored in pigs with different PRRSV infection states to interpret the biological function of pIL-18BP in vivo. The results showed there were five shear mutants of pIL-18BP. The mutant with the longest coding region was selected for subsequent functional validation. First, it was demonstrated that TJM-induced pIL-18BP mRNA expression was higher than that of TJ. A direct interaction between pIL-18BP and pIL-18 was confirmed through fluorescence colocalization, bimolecular fluorescent complimentary (BIFC), and co-immunoprecipitation (CO-IP). pIL-18BP also can regulate pIFN-γ mRNA expression. Finally, the expression of pIL-18BP, pIL-18, and pIFN-γ was explored in different PRRSV infection states. Surprisingly, both mRNA and protein expression of pIL-18 were suppressed. These findings fill the gap in understanding the roles played by pIL-18BP in PRRSV infection and provide a foundation for further research.IMPORTANCEPRRSV-infected pigs elicit a weak cellular immune response and the mechanisms of cellular immune regulation induced by PRRSV have not yet been fully elucidated. In this study, we investigated the role of pIL-18BP in PRRSV-induced immune response referring to the regulation of human IL-18BP to human interferon-gamma (hIFN-γ). This is expected to be used as a method to enhance the cellular immune response induced by the PRRSV vaccine. Here, we mined five transcripts of the pIL-18BP gene and demonstrated that it interacts with pIL-18 and regulates pIFN-γ mRNA expression. Surprisingly, we also found that both mRNA and protein expression of pIL-18 were suppressed under different PRRSV strains of infection status. These results have led to a renewed understanding of the roles of pIL-18BP and pIL-18 in cellular immunity induced by PRRSV infection, which has important implications for the prevention and control of PRRS.

Knockdown of PHOX2B in the retrotrapezoid nucleus reduces the central CO2 chemoreflex in rats.

PHOX2B is a transcription factor essential for the development of different classes of neurons in the central and peripheral nervous system. Heterozygous mutations in the PHOX2B coding region are responsible for the occurrence of Congenital Central Hypoventilation Syndrome (CCHS), a rare neurological disorder characterised by inadequate chemosensitivity and life-threatening sleep-related hypoventilation. Animal studies suggest that chemoreflex defects are caused in part by the improper development or function of PHOX2B expressing neurons in the retrotrapezoid nucleus (RTN), a central hub for CO2 chemosensitivity. Although the function of PHOX2B in rodents during development is well established, its role in the adult respiratory network remains unknown. In this study, we investigated whether reduction in PHOX2B expression in chemosensitive neuromedin-B (NMB) expressing neurons in the RTN altered respiratory function. Four weeks following local RTN injection of a lentiviral vector expressing the short hairpin RNA (shRNA) targeting Phox2b mRNA, a reduction of PHOX2B expression was observed in Nmb neurons compared to both naive rats and rats injected with the non-target shRNA. PHOX2B knockdown did not affect breathing in room air or under hypoxia, but ventilation was significantly impaired during hypercapnia. PHOX2B knockdown did not alter Nmb expression but it was associated with reduced expression of both Task2 and Gpr4, two CO2/pH sensors in the RTN. We conclude that PHOX2B in the adult brain has an important role in CO2 chemoreception and reduced PHOX2B expression in CCHS beyond the developmental period may contribute to the impaired central chemoreflex function.

Ribosome decision graphs for the representation of eukaryotic RNA translation complexity.

The application of ribosome profiling has revealed an unexpected abundance of translation in addition to that responsible for the synthesis of previously annotated protein-coding regions. Multiple short sequences have been found to be translated within single RNA molecules, within both annotated protein-coding and noncoding regions. The biological significance of this translation is a matter of intensive investigation. However, current schematic or annotation-based representations of mRNA translation generally do not account for the apparent multitude of translated regions within the same molecules. They also do not take into account the stochasticity of the process that allows alternative translations of the same RNA molecules by different ribosomes. There is a need for formal representations of mRNA complexity that would enable the analysis of quantitative information on translation and more accurate models for predicting the phenotypic effects of genetic variants affecting translation. To address this, we developed a conceptually novel abstraction that we term ribosome decision graphs (RDGs). RDGs represent translation as multiple ribosome paths through untranslated and translated mRNA segments. We termed the latter "translons." Nondeterministic events, such as initiation, reinitiation, selenocysteine insertion, or ribosomal frameshifting, are then represented as branching points. This representation allows for an adequate representation of eukaryotic translation complexity and focuses on locations critical for translation regulation. We show how RDGs can be used for depicting translated regions and for analyzing genetic variation and quantitative genome-wide data on translation for characterization of regulatory modulators of translation.

Antisense RNA regulates glutamine synthetase in a heterocyst-forming cyanobacterium.

Glutamine synthetase (GS) is a key enzyme involved in nitrogen assimilation and the maintenance of C/N balance, and it is strictly regulated in all bacteria. In cyanobacteria, glutamine synthetase expression is controlled by nitrogen control A (NtcA) transcription factor, which operates global nitrogen regulation in these photosynthetic organisms. Furthermore, post-translational regulation of GS is operated by protein-protein interaction with GS inactivating factors (IFs). Here, we describe an additional regulatory mechanism involving an antisense RNA. In Nostoc sp. PCC 7120, the GS inactivating factor A (gifA) gene (encoding GS inactivating factor IF7) is transcribed downstream of the glutamine synthetase (glnA) gene, from the opposite strand, and the gifA mRNA extends into the glnA coding sequence in antisense orientation. Therefore, the dual RNA transcript that encodes gifA constitutes two functional regions: a 5' protein-coding region, encoding IF7, and a 3' untranslated region that acts as an antisense to glnA. By increasing the levels of such antisense RNA either in cis or in trans, we demonstrate that the amount of GS activity can be modulated by the presence of the antisense RNA. The tail-to-tail disposition of the glnA and gifA genes observed in many cyanobacterial strains from the Nostocales clade suggests the prevalence of such antisense RNA-mediated regulation of GS in this group of cyanobacteria.

Precise correction of a spectrum of β-thalassemia mutations in coding and non-coding regions by base editors

β-thalassemia/HbE results from mutations in the β-globin locus that impede the production of functional adult hemoglobin. Base editors (BEs) could facilitate the correction of the point mutations with minimal or no indel creation, but its efficiency and bystander editing for the correction of β-thalassemia mutations in coding and non-coding regions remains unexplored. Here, we screened BE variants in HUDEP-2 cells for their ability to correct a spectrum of β-thalassemia mutations that were integrated into the genome as fragments ofHBB. The identified targets were introduced into their endogenous genomic location using BEs and Cas9/homology-directed repair (HDR) to create cellular models with β-thalassemia/HbE. These β-thalassemia/HbE models were then used to assess the efficiency of correction in the native locus and functional β-globin restoration. Most bystander edits produced near target sites did not interfere with adult hemoglobin expression and are not predicted to be pathogenic. Further, the effectiveness of BE was validated for the correction of the pathogenic HbE variant in severe β0/βE-thalassaemia patient cells. Overall, our study establishes a novel platform to screen and select optimal BE tools for therapeutic genome editing by demonstrating the precise, efficient, and scarless correction of pathogenic point mutations spanning multiple regions of HBB including the promoter, intron, and exons.

Complete genome sequence of Paraburkholderia terrae strain KU-46, a 2,4-dinitrophenol-degrading bacterium.

Paraburkholderia terrae strain KU-46 has been studied for its capability to degrade 2,4-dinitrophenol. Here, we present the complete 10,833,180bp genome of this microorganism, comprising five circular chromosomes housing 9,797 protein-coding sequences. The genes responsible for 2,4-dinitrophenol and 4-nitrophenol degradation are located on chromosome 2.

Molecular Evolution of RAMOSA1 (RA1) in Land Plants.

RAMOSA1 (RA1) is a Cys2-His2-type (C2H2) zinc finger transcription factor that controls plant meristem fate and identity and has played an important role in maize domestication. Despite its importance, the origin of RA1 is unknown, and the evolution in plants is only partially understood. In this paper, we present a well-resolved phylogeny based on 73 amino acid sequences from 48 embryophyte species. The recovered tree topology indicates that, during grass evolution, RA1 arose from two consecutive SUPERMAN duplications, resulting in three distinct grass sequence lineages: RA1-like A, RA1-like B, and RA1; however, most of these copies have unknown functions. Our findings indicate that RA1 and RA1-like play roles in the nucleus despite lacking a traditional nuclear localization signal. Here, we report that copies diversified their coding region and, with it, their protein structure, suggesting different patterns of DNA binding and protein-protein interaction. In addition, each of the retained copies diversified regulatory elements along their promoter regions, indicating differences in their upstream regulation. Taken together, the evidence indicates that the RA1 and RA1-like gene families in grasses underwent subfunctionalization and neofunctionalization enabled by gene duplication.

The complete chloroplast genome of Durio zibethinus L. cultivar Ri6 (Helicteroideae, Malvaceae)

Durian, a member of the Malvaceae family, is famous for its delicious fruits, which have strong scents and are rich in nutrients. In this study, we sequenced and characterized the complete chloroplast genome of Durio zibethinus L. 1774 cultivar Ri6, a popular durian cultivar in Vietnam, using the Illumina Hiseq platform. The results showed a circular chloroplast genome composed of a large single copy of 96,115 bp, a small single copy of 20,819 bp, and two inverted repeat regions of 24,185 bp. This genome consisted of 79 protein-coding genes, 30 transfer RNA genes, and four ribosomal RNA genes. The overall GC content of this genome was 35.7%. Phylogenetic analysis inferred from 78 protein-coding regions revealed monophyly of Durio species and a close relationship between D. zibethinus cultivar Ri6 and cultivar Mongthong. This study provides essential information for further studies examining genetic population, breedings, and species identification among Durio taxa and cultivars.

Increased ANKRD1 Levels in Early Senescence Mediated by RBMS1-Elicited ANKRD1 mRNA Stabilization.

Cellular senescence is a dynamic biological process triggered by sublethal cell damage and driven by specific changes in gene expression programs. We recently identified ANKRD1 (ankyrin repeat domain 1) as a protein strongly elevated after triggering senescence in fibroblasts. Here, we set out to investigate the mechanisms driving the elevated production of ANKRD1 in the early stages of senescence. Our results indicated that the rise in ANKRD1 levels after triggering senescence using etoposide (Eto) was the result of moderate increases in transcription and translation, and robust mRNA stabilization. Antisense oligomer (ASO) pulldown followed by mass spectrometry revealed a specific interaction of the RNA-binding protein RBMS1 with ANKRD1 mRNA that was confirmed by ribonucleoprotein immunoprecipitation analysis. RBMS1 abundance decreased in the nucleus and increased in the cytoplasm during Eto-induced senescence; in agreement with the hypothesis that RBMS1 may participate in post-transcriptional stabilization of ANKRD1 mRNA, silencing RBMS1 reduced, while overexpressing RBMS1 enhanced ANKRD1 mRNA half-life after Eto treatment. A segment proximal to the ANKRD1 coding region was identified as binding RBMS1 and conferring RBMS1-dependent increased expression of a heterologous reporter. We propose that RBMS1 increases expression of ANKRD1 during the early stages of senescence by stabilizing ANKRD1 mRNA.

Novel polymorphisms in the octopamine receptor gene of amitraz resistant population of Rhipicephalus sanguineus sensu lato, in south India

The brown dog tick or Rhipicephalus sanguineus sensu lato is an ixodid tick, responsible for the dissemination of pathogens that cause canine infectious diseases besides inflicting the direct effects of tick bite. The hot humid climate of Kerala, a south Indian state, is favorable for propagation of tick vectors and acaricides are the main stay of tick control. Though the resistance against synthetic pyrethroids is reported among these species, the status of amitraz resistance in R. sanguineus s. l. in the country is uncertain due to the lack of molecular characterisation data and scarce literature reports. Hence the present study was focused on the phenotypic detection and preliminary genotypic characterisation of amitraz resistance in the R. sanguineus s. l. A modified larval packet test (LPT) on a susceptible isolate was performed to determine the discriminating dose (DD). Further LPT-DD on 35 tick isolates was carried out to detect amitraz resistance robustly, along with that full dose response bioassays on the resistant isolates were performed. The results indicated that amitraz resistance is prevalent with 49 per cent of the samples being resistant. Amplification of exon 3 of octopamine receptor gene from both the susceptible and resistant larval isolates was carried out. Amplicons of ten pooled amitraz susceptible and ten pooled amitraz resistant representative samples were sequenced and analysed, unveiling a total of three novel non-synonymous mutations in the partial coding region at positions V32A, N41D and V58I in phenotypically resistant larval DNA samples. In silico analysis by homology modelling and molecular docking of the mutated and unmutated receptors showed that these mutations had reduced the binding affinity to amitraz. However, lack of mutations in the octopamine receptor gene in three of the pooled low order resistant R. sanguineus s. l. larval samples could be suggestive of other mechanisms associated with amitraz resistance in the region. Hence, further association studies should be carried out to confirm the association of these mutations with target insensitivity in R. sanguineus s. l. ticks, along with exploring the status of metabolic resistance and other mechanisms of resistance.

Evolutionary Rate Shifts in Coding and Regulatory Regions Underpin Repeated Adaptation to Sulfidic Streams in Poeciliid Fishes.

Adaptation to extreme environments often involves the evolution of dramatic physiological changes. To better understand how organisms evolve these complex phenotypic changes, the repeatability and predictability of evolution, and possible constraints on adapting to an extreme environment, it is important to understand how adaptive variation has evolved. Poeciliid fishes represent a particularly fruitful study system for investigations of adaptation to extreme environments due to their repeated colonization of toxic hydrogen sulfide-rich springs across multiple species within the clade. Previous investigations have highlighted changes in the physiology and gene expression in specific species that are thought to facilitate adaptation to hydrogen sulfide-rich springs. However, the presence of adaptive nucleotide variation in coding and regulatory regions and the degree to which convergent evolution has shaped the genomic regions underpinning sulfide tolerance across taxa are unknown. By sampling across seven independent lineages in which nonsulfidic lineages have colonized and adapted to sulfide springs, we reveal signatures of shared evolutionary rate shifts across the genome. We found evidence of genes, promoters, and putative enhancer regions associated with both increased and decreased convergent evolutionary rate shifts in hydrogen sulfide-adapted lineages. Our analysis highlights convergent evolutionary rate shifts in sulfidic lineages associated with the modulation of endogenous hydrogen sulfide production and hydrogen sulfide detoxification. We also found that regions with shifted evolutionary rates in sulfide spring fishes more often exhibited convergent shifts in either the coding region or the regulatory sequence of a given gene, rather than both.

Genome-wide patterns of noncoding and protein-coding sequence variation in the major fungal pathogen Aspergillus fumigatus.

Aspergillus fumigatus is a deadly fungal pathogen, responsible for >400,000 infections/year and high mortality rates. A. fumigatus strains exhibit variation in infection-relevant traits, including in their virulence. However, most A. fumigatus protein-coding genes, including those that modulate its virulence, are shared between A. fumigatus strains and closely related nonpathogenic relatives. We hypothesized that A. fumigatus genes exhibit substantial genetic variation in the noncoding regions immediately upstream to the start codons of genes, which could reflect differences in gene regulation between strains. To begin testing this hypothesis, we identified 5,812 single-copy orthologs across the genomes of 263 A. fumigatus strains. In general, A. fumigatus noncoding regions showed higher levels of sequence variation compared with their corresponding protein-coding regions. Focusing on 2,482 genes whose protein-coding sequence identity scores ranged between 75 and 99%, we identified 478 total genes with signatures of positive selection only in their noncoding regions and 65 total genes with signatures only in their protein-coding regions. Twenty-eight of the 478 noncoding regions and 5 of the 65 protein-coding regions under selection are associated with genes known to modulate A. fumigatus virulence. Noncoding region variation between A. fumigatus strains included single-nucleotide polymorphisms and insertions or deletions of at least a few nucleotides. These results show that noncoding regions of A. fumigatus genes harbor greater sequence variation than protein-coding regions, raising the hypothesis that this variation may contribute to A. fumigatus phenotypic heterogeneity.

Insertion-Deletion Events Are Depleted in Protein Regions with Predicted Secondary Structure.

A fundamental goal in evolutionary biology and population genetics is to understand how selection shapes the fate of new mutations. Here, we test the null hypothesis that insertion-deletion (indel) events in protein-coding regions occur randomly with respect to secondary structures. We identified indels across 11,444 sequence alignments in mouse, rat, human, chimp, and dog genomes and then quantified their overlap with four different types of secondary structure-alpha helices, beta strands, protein bends, and protein turns-predicted by deep-learning methods of AlphaFold2. Indels overlapped secondary structures 54% as much as expected and were especially underrepresented over beta strands, which tend to form internal, stable regions of proteins. In contrast, indels were enriched by 155% over regions without any predicted secondary structures. These skews were stronger in the rodent lineages compared to the primate lineages, consistent with population genetic theory predicting that natural selection will be more efficient in species with larger effective population sizes. Nonsynonymous substitutions were also less common in regions of protein secondary structure, although not as strongly reduced as in indels. In a complementary analysis of thousands of human genomes, we showed that indels overlapping secondary structure segregated at significantly lower frequency than indels outside of secondary structure. Taken together, our study shows that indels are selected against if they overlap secondary structure, presumably because they disrupt the tertiary structure and function of a protein.

Abstract PO4-14-11: Findings of clinically significant variants (Tier IA) with OmniSeq INSIGHT ® in a breast cancer cohort of 987 patients

Abstract Background: There were an estimated 2.3 million new breast cancer cases worldwide in 2020, with breast cancer now the most diagnosed type of cancer representing 25% of cancer cases and 17% of cancer deaths (1). The GLOBOCAN Cancer Tomorrow prediction tool estimates that incidence will increase by > 45% by 2040 (2). Optimal treatment for breast cancer is increasingly dependent upon knowing a patient’s somatic and germline genomic alteration status. To highlight the importance of these genomic alterations, we provide an overview of the genomic findings in 987 consecutive breast cancer patients tested in the course of routine clinical care. Methods: Comprehensive Genomic and Immune Profiling (CGIP) was performed on 987 qualified breast cancer samples at a CAP/CLIA and NYS CLEP certified reference laboratory with the OmniSeq INSIGHT ® test (3). OmniSeq INSIGHT ® is a next generation sequencing-based laboratory developed test for the detection of genomic variants, signatures, HLA Class I genotypes, and immune gene expression in formalin-fixed paraffin-embedded (FFPE) tumor tissue. DNA is sequenced to detect small variants in the full exonic coding region of 523 genes, copy number alterations in 59 genes (gains and losses), as well as analysis of microsatellite instability (MSI) and tumor mutational burden (TMB). RNA is sequenced to detect fusions and splice variants in 55 genes, in addition to mRNA expression in 64 immune genes. The resultant information, along with PD-L1 protein expression by immunohistochemistry (IHC), is intended for use by qualified health care professionals in accordance with professional guidelines in oncology for management of patients with solid neoplasms. Tier IA variants include variants with strong clinical significance as per AMP–ASCO–CAP recommendations (4). Results: There were 384 Tier IA variants detected in 370 of 987 breast cancer patients (~37.5%) with 251 in PIK3CA (65.3%), 60 in ERBB2 (15.6%), 46 in BRCA2 (11.9%), 19 in BRCA1 (4.9%), 4 in NTRK3 (1.0%), 1 in NTRK1 (0.26%) and 1 in PALB2 (0.26%) including copy number variations (CNV), gene fusion and single nucleotide variants (SNV). At the individual gene level for Tier IA variants, PIK3CA had 251 SNVs detected, ERBB2 had 55 CNVs, 1 gene fusion and 4 SNVs detected, BRCA2 had 39 SNVs and 8 CNVs detected, BRCA1 had 19 SNVs, NTRK3 had 4 gene fusions, NTRK1 had 1 gene fusion and PALB2 had 1 SNV detected. At the variant class level, there were 63 CNVs, 6 gene fusions and 315 SNVs observed. Conclusions: CGIP for breast cancer patients identified one or more clinically significant Tier IA genomic alterations that directs targeted therapy in ~ 37.5% of patients in a cohort of real world patients tested during the standard course of clinical care. This highlights the need for comprehensive genomic testing in breast cancer patients to drive therapeutic decision making.

A unified classification system for HIV-1 5' long terminal repeats.

The HIV-1 provirus mainly consists of internal coding region flanked by 1 long terminal repeats (LTRs) at each terminus. The LTRs play important roles in HIV-1 reverse transcription, integration, and transcription. However, despite of the significant study advances of the internal coding regions of HIV-1 by using definite reference classification, there are no systematic and phylogenetic classifications for HIV-1 5' LTRs, which hinders our elaboration on 5' LTR and a better understanding of the viral origin, spread and therapy. Here, by analyzing all available resources of 5' LTR sequences in public databases following 4 recognized principles for the reference classification, 83 representatives and 14 consensus sequences were identified as representatives of 2 groups, 6 subtypes, 6 sub-subtypes, and 9 CRFs. To test the reliability of the supplemented classification system, the constructed references were applied to identify the 5' LTR assignment of the 22 clinical isolates in China. The results revealed that 16 out of 22 tested strains showed a consistent subtype classification with the previous LTR-independent classification system. However, 6 strains, for which recombination events within 5' LTR were demonstrated, unexpectedly showed a different subtype classification, leading a significant change of binding sites for important transcription factors including SP1, p53, and NF-κB. The binding change of these transcriptional factors would probably affect the transcriptional activity of 5' LTR. This study supplemented a unified classification system for HIV-1 5' LTRs, which will facilitate HIV-1 characterization and be helpful for both basic and clinical research fields.

Evolutionary trends of respiratory syncytial viruses: Insights from large-scale surveillance and molecular dynamics of G glycoprotein

Human respiratory syncytial virus (RSV) is an underlying cause of lower respiratory illnesses in children, elderly and immunocompromised adults. RSV contains multiple structural and non-structural proteins with two major glycoproteins that control the initial phase of infection, fusion glycoprotein and the attachment (G) glycoprotein. G protein attaches to the ciliated cells of airways initiating the infection. The hypervariable G protein plays a vital role in evolution of RSV strains. We employed multiple bioinformatics tools on systematically accessed large-scale data to evaluate mutations, evolutionary history, and phylodynamics of RSV. Mutational analysis of central conserved region (CCR) on G protein-coding sequences between 163 and 189 positions revealed frequent mutations at site 178 in human RSV (hRSV) A while arginine to glutamine substitutions at site 180 positions in hRSV B, remained prevalent from 2009 to 2014. Phylogenetic analysis indicates multiple signature mutations within G protein responsible for diversification of clades. The USA and China have highest number of surveillance records, followed by Kenya. Markov Chain Monte Carlo Bayesian skyline plot revealed that RSV A evolved steadily from 1990 to 2000, and rapidly between 2003 and 2005. Evolution of RSV B continued from 2003 to 2022, with a high evolution stage from 2016 to 2020. Throughout evolution, cysteine residues maintained their strict conserved states while CCR has an entropy value of 0.0039(±0.0005). This study concludes the notion that RSV G glycoprotein is continuously evolving while the CCR region of G protein maintains its conserved state providing an opportunity for CCR-specific monoclonal antibodys (mAbs) and inhibitors as potential candidates for immunoprophylaxis.

Regional-specific calibration enables application of computational evidence for clinical classification of 5′ cis-regulatory variants in Mendelian disease

To date, clinical genetic testing for Mendelian disease variants has focused heavily on exonic coding and intronic gene regions. This multi-step study was undertaken to provide an evidence base for selecting and applying computational approaches for use in clinical classification of 5' cis-regulatory region variants. Curated datasets of clinically reported disease-causing 5' cis-regulatory region variants and variants from matched genomic regions in population controls were used to calibrate six bioinformatic tools as predictors of variant pathogenicity. Likelihood ratio estimates were aligned to code weights following ClinGen recommendations for application of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG/AMP) classification scheme. Considering code assignment across all reference dataset variants, performance was best for CADD (81.2%) and REMM (81.5%). Optimized thresholds provided moderate evidence toward pathogenicity (CADD, REMM) and moderate (CADD) or supporting (REMM) evidence against pathogenicity. Both sensitivity and specificity of prediction were improved when further categorizing variants based on location in an EPDnew-defined promoter region. Combining predictions (CADD, REMM, and location in a promoter region) increased specificity at the expense of sensitivity. Importantly, the optimal CADD thresholds for assigning ACMG/AMP codes PP3 (≥10) and BP4 (≤8) were vastly different from recommendations for protein-coding variants (PP3≥25.3; BP4≤22.7); CADD <22.7 would incorrectly assign BP4 for >90% of reported disease-causing cis-regulatory region variants. Our results demonstrate the need to consider a tiered approach and tailored score thresholds to optimize bioinformatic impact prediction for clinical classification of 5' cis-regulatory region variants.

Signaling of GPR17 in Hepatic and Renal Tissues

G protein-coupled receptors (GPCRs) are the largest class of proteins in the human body and are considered to be excellent targets for drug discovery. However, the vast majority of GPCRs are understudied. Using a large-scale TaqMan array screen on both male and female C57BL6 murine livers, we identified GPR17 as one of the most highly expressed orphan receptors within the liver (ΔCt = 19.9, normalized to GAPDH). GPR17 expression in the liver was confirmed using reverse transcriptase polymerase chain reaction (RT-PCR) and was also identified in murine whole kidney from both males and females indicating that GPR17 RNA expression is not sex-dependent. While GPR17 has known roles in oligodendrocyte maturation and myelination in the central nervous system, its functions within the liver and kidney are undetermined. To investigate its physiological functions, we started by uncovering its localization within the liver and kidney and identifying its ligand profile. A series of RNAScope in-situ hybridization assays were performed to localize GPR17 within hepatic tissue. Gpr17 mRNA was found within hepatocytes indicating that its highest expression is in the parenchymal cells of the liver. To identify possible activators for GPR17, we cloned the coding region into both flag-tagged PME18s and the Presto-Tango vectors. Live-cell surface labeling of HEK293T expressing the GPR17 constructs confirmed robust surface expression on the plasma membrane facilitating ligand screening. GPR17 downstream signaling has yet to be studied. Thus, we used the PRESTO-TANGO assay, an exogenous chemical ligand screening assay that relies on a luminescence readout to indicate -arrestin recruitment and receptor activation. A large, non-biased chemical screen of over 500 physiological compounds was tested for activation of GPR17. To date, we identified seven exogenous activators of GPR17: 3-[methyl(1,2,3,4- tetrahydronaphthalen-2-yl)amino]pyrazine-2-carbonitrile, 4-(4-fluorophenyl)-5-methyl-N-(2-methyl-1-oxo-1,2,3,4-tetrahydroisoquinolin-7-yl)thiophene-2-carboxamide,({1-[(2-chlorophenyl)methyl]-3,5-dimethyl-1H-pyrazol-4-yl}methyl)(methyl)[(1-methyl-1H-imidazol-2-yl)methyl]amine,N,N-dimethyl-6-[(5-phenyl-1,3-oxazol-2-yl)methyl]-5H,6H,7H,8H-pyrido[4,3-d]pyrimidin-4-amine,4-phenyl-1-[6-(pyrrolidin-1-yl)-9H-purin-9-yl]butan-2-ol, 2-(3-fluorophenyl)-3-methyl-1-{[1,2,4]triazolo[4,3-a]pyridine-6-carbonyl}pyrrolidine, and 2,6,6-trimethyl-N-({2-methyl-4H,5H,6H,7H-pyrazolo[1,5-a]pyrimidin-6-yl}methyl)-4,5,6,7-tetrahydro-1-benzofuran-4-amine. Ongoing studies are aimed at narrowing down the in vivo activators of GPR17 to ultimately shed light on the physiological functions that this receptor plays in the kidney and liver. NIH/NCATS. This is the full abstract presented at the American Physiology Summit 2024 meeting and is only available in HTML format. There are no additional versions or additional content available for this abstract. Physiology was not involved in the peer review process.

Abstract 137: Identification Of A Noncoding Genetic Variant In The Tissue Factor Locus That Reverses Lethal Thrombosis In Mice

Background: Thrombosis is initiated by tissue factor (TF, gene name F3 ) binding to coagulationFVII, with tissue factor pathway inhibitor (TFPI) inhibiting this complex. Alterations in TF orTFPI expression significantly affect thrombosis. Reducing TFPI expression by 50% ( Tfpi +/-) inmice results in a perinatal lethal phenotype on the Factor V Leiden homozygous( F5 L/L )prothrombotic background. We used the F5 L/L Tfpi +/- lethal phenotype to conduct a dominantsensitized whole genome ENU mutagenesis screen to suppress the F5 L/L Tfpi +/- lethality. Weidentified a Modifier of Factor 5 Leiden 6 (MF5L6) line with 72% penetrance and 85 F5 L/L Tfpi +/- offspring. A significant linkage peak (LOD=4.35),explaining half the suppressing effect andcontaining F3 (Chromosome 3) was identified. Goals/Hypothesis: To identify the genomic variant controlling F3 expression in the MF5L6 line. Methods: To quantify F3 expression in the surviving mice from MF5L6, quantitative PCR onliver, lung, and heart tissues from MF5L6 was performed. We used Sanger DNA and highthroughput sequencing to identify candidate TF regulatory variants in the F3 locus. Theprothrombin time assay was used to test the effects of reduced TF expression on in vitro bloodcoagulation. Results: Two distinct expression profiles in the lung and liver of the MF5L6 mice wereobserved, those that had a 50% reduction in F3 mRNA and those that did not. Heart tissuesexhibited one expression profile, suggesting that the variant regulates F3 expressiontissue-specifically. Sanger sequencing of the F3 coding region revealed no coding mutations inMF5L6 mice. Whole genome sequencing identified two novel candidate variants (in unknownF3 regulatory elements) in the 200 kilobase upstream region of F3 . The 50% reduction in F3 resulted in significant changes in coagulation by the prothrombin assay (n=18,p&lt;0.0009). Conclusion: We identified novel candidate variants for regulating F3 gene expression and aredetermining their mechanism of action. Investigation of these variants will provide new insightsinto the regulation of F3 and enable us to identify the variant(s) responsible for the remainder ofthe thrombosis suppressing effect in MF5L6. Our findings provide new insights into the geneticregulation of thrombosis.