Non-synonymous Variants Research Articles

A Novel High-Throughput Method for Identifying T-Cell Receptor: Minor Histocompatibility Antigen Interactions in Mouse Models of Graft-Vs-Host Disease

Graft-vs-host disease is a complication of allogeneic stem cell transplant (alloSCT) in which donor T cells target alloantigens expressed by healthy recipient tissues. In major-histocompatibility-matched alloSCT, CD8-mediated GVHD is caused by T cell receptor (TCR) allorecognition of class-I presented minor histocompatibility antigens (miHAs). The vast majority of miHAs remain unknown between any given donor-recipient pairing. However, determining the identities of these miHAs and their cognate TCRs would prove useful towards the development of allografts that minimize GVHD yet preserve the therapeutic benefits of alloSCT. Here, we implement a high-throughput TCR cloning method, TCXpress, to screen alloreactive TCRs from a mouse model of MHC-matched alloSCT against putative miHA targets determined informatically. 129 (H-2 b) mice, which express the immunodominant miHA H60, were lethally irradiated and reconstituted with bone marrow and CD4 + and CD8 + T cells from B6 (H-2 b) donors. At day +10, recipient spleens were harvested and CD44 +PD-1 + donor CD8 + cells were single-cell sorted according to H-2Kb:H60 tetramer (Tet H60) staining into unknown miHA-reactive (Tet H60-) and control H60-reactive (Tet H60+) cohorts for TCR cloning using the TCXpress platform (BlueSphere Bio). Clonotyping was performed by sequencing the CDR3 region of TCRs. Putative miHAs were identified informatically by filtering non-synonymous variants between B6 and 129 reference exomes on predicted class I binding affinities. To interrogate TCR target specificity, Jurkat reporter cells, which express GFP and CD69 upon TCR activation, were transduced with selected alloreactive TCRs and reacted in a multiplexed fashion against a panel of B6WT3 antigen-presenting cells (H-2 b) transduced with tandem minigenes (TMG) of predicted miHAs. To test TCR reactivity against hematopoietically-expressed miHAs independent of our calling strategy, Jurkat reporters were also reacted against bone-marrow derived dendritic cells (BMDCs) from control B6 and 129 mice. From 3 recipient mice, 966 alloreactive TCRs were cloned and sequenced (735 Tet H60- and 231 Tet H60+) corresponding to 553 unique clonotypes (437 Tet H60- and 116 Tet H60+). Sixty-seven percent (N=668) of TCRs sequenced had clonotype occurrences of ≥ 2, indicating clonal expansion post-alloSCT, whereas there were no repeated clonotypes among donor B6 CD8 + cells pre-transplant (N=173 TCRs). A total of 47 Tet H60- TCR clones and 3 Tet H60+ control TCR clones were selected for transduction into Jurkat reporters and screened against a limited TMG library of 158 putative B6→129 miHAs, as well as H60. Expectedly, all TCR clonotypes isolated from the Tet H60+ sort gate reacted against H60-expressing B6WT3 cells and 129 BMDCs. From the Tet H60- TCRs, one clone was validated against another documented 129 miHA, H4, while an additional 4 clones were found to react against novel B6→129 miHAs (Serpinb8 and Cdt1). Interestingly, 2 highly over-represented clones were autoreactive against B6 tissues (B6WT3 cells and BMDCs), but not 129 BMDCs. In summary, we demonstrate the capacity of a novel TCR cloning platform to validate informatically-predicted miHAs and their cognate TCRs in a mouse model of GVHD. This approach may be readily adapted to develop new TCR:miHA therapeutics and gain a deeper understanding of GVH responses targeting miHAs.

Primary lipoblastic nerve sheath tumor in an inguinal lymph node mimicking metastatic tumor: a case report and literature review

Lipoblastic nerve sheath tumors of soft tissue are characterized as schwannoma tumors that exhibit adipose tissue and lipoblast-like cells with signet-ring morphology. They have been documented to arise in various anatomic locations, including the thigh, groin, shoulder, and retroperitoneum. However, to our knowledge, this tumor has not been previously reported as a lymph node primary. We present herein the first case of a benign primary lipoblastic nerve sheath tumor arising in an inguinal lymph node in a 69-year-old man. Microscopic examination revealed a multinodular tumor comprising fascicles of spindle cells, as well as adipocytic and lipoblast-like signet-ring cell component in the context of schwannoma. Despite the presence of some bizarre cells with nuclear atypia, no obvious mitotic activity or necrosis was observed. Immunohistochemical analysis showed strong and diffuse expression of S-100, SOX10, CD56, and NSE in the spindle cells as well as in the signet-ring lipoblast-like cells and the mature adipocytes. Sequencing analysis of the neoplasm identified six non-synonymous single nucleotide variant genes, specifically NF1, BRAF, ECE1, AMPD3, CRYAB, and NPHS1, as well as four nonsense mutation genes including MRE11A, CEP290, OTOA, and ALOXE3. The patient remained alive and well with no evidence of recurrence over a period of ten-year follow-up.

Genetic susceptibility to post-endoscopic retrograde cholangiopancreatography pancreatitis identified in propensity score-matched analysis.

A previous history of post-endoscopic retrograde cholangiopancreatography (ERCP) pancreatitis (PEP) is a risk factor for PEP, suggesting that there may be a genetic predisposition to PEP. However, nothing is known about this yet. The aim of this study was to identify genetic variations associated with PEP. A cohort of high-risk PEP patients was queried from December 2016 to January 2019. For each PEP case, two propensity score-matched controls were selected. Whole exome sequencing was performed using blood samples. Genetic variants reported to be related to pancreatitis were identified. To discover genetic variants that predispose to PEP, a logistic regression analysis with clinical adjustment was performed. Gene-wise analyses were also conducted. Totals of 25 PEP patients and 50 matched controls were enrolled. Among the genetic variants reported to be associated with pancreatitis, only CASR rs1042636 was identified, and it showed no significant difference between the case and control groups. A total of 54,269 non-synonymous variants from 14,313 genes was identified. Logistic regression analysis of these variants showed that the IRF2BP1 rs60158447 GC genotype was significantly associated with the occurrence of PEP (odds ratio 2.248, FDR q value = 0.005). Gene-wise analyses did not show any significant results. This study found that the IRF2BP1 gene variant was significantly associated with PEP. This genetic variant is a highly targeted PEP risk factor candidate and can be used for screening high-risk PEP groups before ERCP through future validation. (ClinicalTrials.gov no. NCT02928718).

Clinical relevance of genetic variants in the von Willebrand factor according to in-silico methods.

Clinical interpretation of genetic variants in the context of the patient's phenotype is a time-consuming and costly process. In-silico analysis using in-silico prediction tools, and molecular modeling have been developed to predict the influence of genetic variants on the quality and/or quantity of the resulting translated protein, and in this way, to alert clinicians of disease likelihood in the absence of previous evidence. Our objectives were to evaluate the success rate of the in-silico analysis in predicting the disease-causing variants as pathogenic and the single-nucleotide variants as neutral, and to establish the reliability of in-silico analysis for determining pathogenicity or neutrality of von Willebrand factor gene-associated genetic variants. Using in-silico analysis, we studied pathogenicity in 31 disease-causing variants, and neutrality in 61 single-nucleotide variants from patients previously diagnosed as type 2 von Willebrand disease. Disease-causing variants and non-synonymous single-nucleotide variants were explored by in-silico tools that analyze the amino acidic sequence. Intronic and synonymous single-nucleotide variants were analyzed by in-silico methods that evaluate the nucleotidic sequence. We found a consistent agreement between predictions achieved by in-silico prediction tools and molecular modeling, both for defining the pathogenicity of disease-causing variants and the neutrality of single-nucleotide variants. Based on our results, the in-silico analysis would help to define the pathogenicity or neutrality in novel genetic variants observed in patients with clinical and laboratory phenotypes suggestive of von Willebrand disease.

Mutations in Coding and Non-Coding Regions in Varicella-Zoster Virus Causing Fatal Hemorrhagic Fever Without Rash in an Immunocompetent Patient: Case Report.

We report the case of a fatal hemorrhagic varicella primary infection in an immunocompetent man and whole-genome characterization of the virus for the investigation of biomarkers of virulence. A 38-year-old patient born in Nigeria presented to the emergency department with abdominal pain and subsequently developed fatal hemorrhagic disease without skin rash. Extensive laboratory tests including serology and PCR for arenaviruses, bunyaviruses and ebolaviruses were negative. Varicella-zoster virus (VZV) PCR of sera, liver and spleen tissue samples from autopsy revealed the presence of VZV DNA. Primary infection by varicella-zoster virus with hemorrhagic manifestations was diagnosed after virological testing. The VZV genome was sequenced using a mWGS approach. Bioinformatic analysis showed 53 mutations across the genome, 33 of them producing non-synonymous variants affecting up to 14 genes. Some of them, such as ORF11 and ORF 62, encoded for essential functions related to skin or neurotropism. To our knowledge, the mutations reported here have never been described in a VZV causing such a devastating outcome. In immunocompetent patients, viral factors should be considered in patients with uncommon symptoms or severe diseases. Some relevant mutations revealed by using whole genome sequencing (WGS) directly from clinical samples may be involved in this case and deserves further investigation. Differential diagnosis of varicella-zoster virus in immunocompetent adults should be considered among patients with suspected VHF, even if the expected vesicular rash is not present at admission and does not arise thereafter. Whole genome sequencing of strains causing uncommon symptoms and/or mortality is needed for epidemiological surveillance and further characterization of putative markers of virulence. Additionally, this report highlights the recommendation for a VZV vaccination policy in non-immunized migrants from developing countries.

Interleukin-37: associations of plasma levels and genetic variants in gout

ObjectivesIL-37 is an anti-inflammatory cytokine involved in inflammatory and autoimmune diseases. We aimed to investigate the association between IL-37 genetic variants, IL-37 plasma levels, and various clinical phases of gout.MethodsThe study included a control group with no history of primary hyperuricemia/gout, (n = 50), asymptomatic hyperuricemia (n = 74), intercritical gout (n = 200), acute gouty flare (n = 18), and chronic tophaceous gout (n = 30). Plasma IL-37 was analysed using enzyme-linked immunosorbent assay. All coding regions and intron–exon boundaries of IL-37 and exons 1–5 were amplified and sequenced.ResultsPlasma levels of IL-37 were significantly higher in asymptomatic hyperuricemic (p = 0.045), intercritical gout (p = 0.001), and chronic tophaceous gout (p = 0.021) cohorts when compared to control group. The levels of IL-37 in patients with acute gouty flare were comparable to control group (p = 0.061). We identified 15 genetic variants of IL-37: eight intron (rs2708959, rs2723170, rs2708958, rs2723169 rs2466448, rs3811045, rs3811048, rs2708944) and seven non-synonymous allelic variants (rs3811046, rs3811047, rs2708943, rs2723183, rs2723187, rs2708947, rs27231927), of which rs2708959 showed an over-presentation in gouty and acute flare cohorts (p = 0.003 and 0.033, respectively) compared to European population (minor allelic frequency MAF = 0.05) but not in control and hyperuricemic cohorts (p/MAF = 0.17/0.08 and 0.71/0.05, respectively).. On the contrary, rs3811045, rs3811046, rs3811047, and rs3811048 were underrepresented among individuals with tophaceous gout (MAF = 0.57) compared to European MAF 0.70–0.71, but not compared to the control cohort (MAF = 0.67).ConclusionsWe demonstrated the up-regulation of IL-37 levels across the clinical phases of gout: asymptomatic hyperuricemia, intercritical, and chronic tophaceous gout compared to control. Moreover, 15 genetic variants of IL-37 were identified and their associations with the clinical variants of gout were evaluated.

A whole exome sequencing study to identify rare variants in multiplex families with alcohol use disorder.

Alcohol use disorder (AUD) runs in families and is accompanied by genetic variation. Some families exhibit an extreme susceptibility in which multiple cases are found and often with an early onset of the disorder. Large scale genome-wide association studies have identified several genes with impressive statistical probabilities. Most of these genes are common variants. Our goal was to perform exome sequencing in families characterized by multiple cases (multiplex families) to determine if rare variants might be segregating with disease status. A case-control approach was used to leverage the power of a large control sample of unrelated individuals (N = 8,983) with exome sequencing [Institute for Genomic Medicine (IGM)], for comparison with probands with AUD (N = 53) from families selected for AUD multiplex status. The probands were sequenced at IGM using similar protocols to those used for the archival controls. Specifically, the presence of a same-sex pair of adult siblings with AUD was the minimal criteria for inclusion. Using a gene-based collapsing analysis strategy, a search for qualifying variants within the sequence data was undertaken to identify ultra-rare non-synonymous variants. We searched 18,666 protein coding genes to identify an excess of rare deleterious genetic variation using whole exome sequence data in the 53 AUD individuals from a total of 282 family members. To complete a case/control analysis of unrelated individuals, probands were compared to unrelated controls. Case enrichment for 16 genes with significance at 10-4 and one at 10-5 are plausible candidates for follow-up studies. Six genes were ultra rare [minor allele frequency (MAF) of 0.0005]: CDSN, CHRNA9, IFT43, TLR6, SELENBP1, and GMPPB. Eight genes with MAF of 0.001: ZNF514, OXGR1, DIEXF, TMX4, MTBP, PON2, CRHBP, and ANKRD46 were identified along with three protein-truncating variants associated with loss-of-function: AGTRAP, ANKRD46, and PPA1. Using an ancestry filtered control group (N = 2,814), nine genes were found; three were also significant in the comparison to the larger control group including CHRNA9 previously implicated in alcohol and nicotine dependence. This study implicates ultra-rare loss-of-function genes in AUD cases. Among the genes identified include those previously reported for nicotine and alcohol dependence (CHRNA9 and CRHBP).

Novel Melanocortin-3 and -4 Receptor Functional Variants in Asian Children With Severe Obesity.

Genetic variants in melanocortin 3 receptor (MC3R) and melanocortin 4 receptor (MC4R) genes are strongly associated with childhood obesity. This study aims to identify and functionally characterize MC3R and MC4R variants in an Asian cohort of children with severe early-onset obesity. Whole-exome sequencing was performed to screen for MC3R and MC4R coding variants in 488 Asian children with severe early-onset obesity (body mass index for age ≥97th percentile). Functionality of the identified variants were determined via measurement of intracellular cyclic adenosine monophosphate (cAMP) concentrations and luciferase activity. Four MC3R and 2 MC4R heterozygous nonsynonymous rare variants were detected. There were 3 novel variants: MC3R c.151G > C (p.Val51Leu), MC4R c.127C > A (p.Gln43Lys), and MC4R c.272T > G (p.Met91Arg), and 3 previously reported variants: MC3R c.127G > A (p.Glu43Lys), MC3R c.97G > A (p.Ala33Thr), and MC3R c.437T > A (p.Ile146Asn). Both MC3R c.127G > A (p.Glu43Lys) and MC4R c.272T > G (p.Met91Arg) variants demonstrated defective downstream cAMP signaling activity. The MC4R c.127C > A (p.Gln43Lys) variant showed reduced cAMP signaling activity at low substrate concentration but the signaling activity was restored at high substrate concentration. The MC3R c.151G > C (p.Val51Leu) variant did not show a significant reduction in cAMP signaling activity compared to wild-type (WT) MC3R. Coexpression studies of the WT and variant MC3R/MC4R showed that the heterozygous variants did not exhibit dominant negative effect. Our functional assays demonstrated that MC3R c.127G > A (p.Glu43Lys) and MC4R c.272T > G (p.Met91Arg) variants might predispose individuals to early-onset obesity, and further studies are needed to establish the causative effect of these variants in the pathogenesis of obesity.

THU209 Small INDELs In The Androgen Receptor Gene: Phenotype Implications And Mechanisms Of Mutagenesis In Androgen Insensitivity Syndrome

Abstract Disclosure: R.L. Batista: None. R.M. Ramos: None. A.F. Afonso: None. R.J. Petroli: None. M.Y. Nishi: None. M.M. Ferrari: None. S. Domenice: None. B.B. Mendonca: None. Small indels are insertions or deletions of one or more nucleotides into the DNA sequence. Indels are highly abundant in human genomes and have contributed massively to genome evolution. However, the role of indels and their underlined mechanisms of mutagenesis in Mendelian disorders still needs to be appreciated. To explore the contribution of indels in a monogenic disorder, we analyzed all indels in the androgen receptor gene, including three novel indels, using Androgen Insensitivity Syndrome (AIS) as a model. The goal was to profile indels in AIS. To build this profile, we analyzed the indel distribution through the AR coding region, compared with genomic and exonic populational data, and analyzed their impact on the AIS phenotype and their likely mechanism of mutagenesis. We present 82 different indels in the AR gene reported in individuals with AIS. In most cases (n=78; 95.1%), the phenotype is complete AIS (CAIS). There are four AR indels reported in mild AIS. Most short indels in the AR gene among AIS are frameshift indels (n=64; 79%). Indels are distributed along the entire AR gene, with NTD domain predominance (n=45; 54.9%), followed by LTD domain (n=28; 34.1%). Most are small deletions (n=49; 59.7%), followed by small insertions (n=29; 35.3%) and complex frameshift indels (n=4; 4.8%). Most indels are short (<10 bp; n=81 = 98.8%), with an average of 2.4 nucleotides by indel. Most indels are at non-repetitive sequences with no changes in copy counting (n = 68; 82.9%). Thirteen indels were identified within tandem repeats (15.8%), whereas seven indels were at homopolymers run (8.5%). Trinucleotide regions (Poly-Q, Poly-G) are prone to indels, but indels also happen in other trinucleotide stretches (Poli-L, Poly-P). Regarding the change in copy counting, thirty-nine (47.5%) indels are CCC indels. Indel reported in MAIS is all located at tandem repeats, while all at non-repetitive sequences caused CAIS. 3.65% of AR indels are at a 6bp consensus sequence (TG A/G A/G G/T A/C). In summary, indels represent a significant portion of the variants identified in the AR gene in patients with AIS. Unlike the non-synonymous variants related to this gene, there is a genotype-phenotype correlation (indel - CAIS congruency). However, small NFS indels at small triplet runs can slightly compromise AR function leading to the mild phenotype. Homopolymer runs, and trinucleotide repeats are exonic regions at risk for indels. Slippage explains half of the cases of AR indels. These results increase the molecular understanding of the androgen receptor disruption and the underlined mechanisms of indel mutagenesis in exonic regions and their implication on phenotype. Presentation: Thursday, June 15, 2023

In-silico prediction of TGF-β1 non-synonymous variants and their impact on binding affinity to Fresolimumab

TGF-β1 is a potent immunoregulatory cytokine that plays diverse roles in development, bone healing, fibrosis, and cancer. However, characterizing TGF-β1 gene variants is challenging because the structural and functional consequences of these variants are still undetermined. In this study, we aimed to perform an in-silico analysis of TGF-β1 non-synonymous variants and their pathogenic effects on the TGF-β1 protein. A total of 10,252 TGF-β1 SNPs were collected from the NCBI dbSNP database and in-silico tools (SIFT, PROVEAN, Mutation Taster, ClinVar, PolyPhen-2, CScape, MutPred, and ConSurf) were used. The in-silico predicted potential variants were further investigated for their binding to the TGF-β1 targeting drug “Fresolimumab”. Molecular docking was performed using HADDOCK and confirmed by PRODIGY and PDBsum. The in-silico analysis predicted four potential TGF-β1 nsSNPs: E47G in the LAP domain of the propeptide and I22T, L28F, and E35D in the mature TGF-β1 peptide. HADDOCK and molecular dynamics simulations revealed that the I22T and E35D variants have higher binding affinity for Fresolimumab as compared to the wild type and L28F variants. Molecular dynamics simulations (100 ns) and principal component analysis showed that TGF-β1 variants influenced the protein structure and caused variations in the internal dynamics of protein complexes with the antibody. Among them, the E35D variant significantly destabilized the TGF-β1 protein structure, resulting in rearrangement in the binding site and affecting the interactions with the Fresolimumab. This study identified four variants that can affect the TGF-β1 protein structure and result in functional consequences such as impaired response to Fresolimumab. Communicated by Ramaswamy H. Sarma

Detection of genomic variations and selection signatures in Wagyu using whole-genome sequencing data.

Wagyu is recognized for producing marbled beef with high nutritional value and flavor. Reportedly, Wagyu has been widely used to improve the meat quality of local breeds around the world. However, studies on the genetic mechanism of meat quality in Wagyu at the whole-genome level are rarely reported. Here, whole-genome sequencing data of 11 Wagyu and 115 other individuals were used to explore the genomic variations and genes under selection pressure in Wagyu. A total of 31 349 non-synonymous variants and 53 102 synonymous variants were identified in Wagyu. The population structure analysis showed that Wagyu had the closest genetic relationship with Mishima-Ushi cattle and was apparently separated from other cattle breeds. Then, composite likelihood ratio (CLR), integrated haplotype score, fixation index and cross-population composite likelihood ratio (XP-CLR) tests were performed to identify the candidate genes under positive selection in Wagyu. In total, 770 regions containing 312 genes were identified by at least three methods. Among them, 97 regions containing 27 genes were detected by all four methods. We specifically illustrate a list of interesting genes, including LRP2BP, GAA, CACNG6, CXADR, GPCPD1, KLF2, KLF13, SOX5, MYBPC1, SLC25A10, ATP8A1 and MYH15, which are associated with lipid metabolism, fat deposition, muscle development, bone development, feed intake and growth traits in Wagyu. This is the first study to explore the genomic variations and selection signatures of Wagyu at the whole-genome level. These results will provide significant help to beef cattle improvement and breeding.

Genome-Wide Association Analysis of Protein-Coding Variants in IgA Nephropathy.

Genome-wide association studies have identified nearly 20 IgA nephropathy susceptibility loci. However, most nonsynonymous coding variants, particularly ones that occur rarely or at a low frequency, have not been well investigated. The authors performed a chip-based association study of IgA nephropathy in 8529 patients with the disorder and 23,224 controls. They identified a rare variant in the gene encoding vascular endothelial growth factor A (VEGFA) that was significantly associated with a two-fold increased risk of IgA nephropathy, which was further confirmed by sequencing analysis. They also identified a novel common variant in PKD1L3 that was significantly associated with lower haptoglobin protein levels. This study, which was well-powered to detect low-frequency variants with moderate to large effect sizes, helps expand our understanding of the genetic basis of IgA nephropathy susceptibility. Genome-wide association studies have identified nearly 20 susceptibility loci for IgA nephropathy. However, most nonsynonymous coding variants, particularly those occurring rarely or at a low frequency, have not been well investigated. We performed a three-stage exome chip-based association study of coding variants in 8529 patients with IgA nephropathy and 23,224 controls, all of Han Chinese ancestry. Sequencing analysis was conducted to investigate rare coding variants that were not covered by the exome chip. We used molecular dynamic simulation to characterize the effects of mutations of VEGFA on the protein's structure and function. We also explored the relationship between the identified variants and the risk of disease progression. We discovered a novel rare nonsynonymous risk variant in VEGFA (odds ratio, 1.97; 95% confidence interval [95% CI], 1.61 to 2.41; P = 3.61×10 -11 ). Further sequencing of VEGFA revealed twice as many carriers of other rare variants in 2148 cases compared with 2732 controls. We also identified a common nonsynonymous risk variant in PKD1L3 (odds ratio, 1.16; 95% CI, 1.11 to 1.21; P = 1.43×10 -11 ), which was associated with lower haptoglobin protein levels. The rare VEGFA mutation could cause a conformational change and increase the binding affinity of VEGFA to its receptors. Furthermore, this variant was associated with the increased risk of kidney disease progression in IgA nephropathy (hazard ratio, 2.99; 95% CI, 1.09 to 8.21; P = 0.03). Our study identified two novel risk variants for IgA nephropathy in VEGFA and PKD1L3 and helps expand our understanding of the genetic basis of IgA nephropathy susceptibility.

Genomic Characterization of Low and High Grade Cervical Intraepithelial Neoplasia in Comparison to Cervical Cancer

Screening reduces incidence of cervical cancer (CC) through identification and treatment of cervical intraepithelial neoplasia (CIN). Triage of CIN must balance between overtreatment versus progression as invasive therapies increase the risk of side effects including obstetric complications. The majority of CIN 1-2 lesions regress and a subset of CIN 3 lesions progress, however, due to inability to differentiate lesions likely to progress, the majority of CIN 2-3 lesions are treated with resection. Improved genomic and molecular identifiers of disease represent an urgent unmet clinical need. Using next generation sequencing of a targeted exome panel of 1109 genes (previously validated), we characterized somatic mutations in 36 CIN (14 CIN 1, 11 CIN 2, 11 CIN 3) and 13 CC samples. Sequencing of CIN samples was performed on exfoliated cervical cells. CIN diagnosis and grade was confirmed on biopsy. Mutation profiles between CIN grades and CC were compared to identify genomic patterns that distinguish these groups. Across the 49 samples sequenced, we identified a total of 5142 somatic mutations, including 2178 missense, 2522 synonymous, 171 nonsense, 62 splice site, 135 inframe indels, and 74 frameshift mutations. The mutation frequency was significantly higher in CC vs CIN1-3 (Table), p<0.01 (Wilcoxon signed-rank test). The difference in mutation frequency and type (Table) was not significant between CIN grades 1, 2 and 3 (p = 0.07, Kruskal-Wallis test). Cancer related pathway signatures were analyzed for the percentage of CC vs CIN samples with at least one altered variant as follows: RTK-RAS (92% vs 31%), PI3K (92% vs 17%), NOTCH (100% vs 39%), WNT (85% vs 14%) and cell cycle (53% vs 3%). Ninety percent of CC samples versus only 30% of CINs had nonsynonymous variants in the RTK-RAS and NOTCH pathways as well as the PI3K pathway which is implicated as a late event in cervical carcinogenesis. We observed recurrent missense variants in ABL1, IGF1R, TSC2 in the CIN2, CIN3 and CC samples, particularly in genes that belong to the RTK-RAS, PI3K, NOTCH pathway signatures. Potential driver mutations in EGFR, PIK3CA, ERBB4, MTOR, CSF1R genes were exclusive to CC samples. Our results show that there is a clear distinction of mutational burden and type between cervical cancer and CIN1-3. However, the pathologic grading of CIN is not consistently associated with dynamic changes in overall mutation burden, type, or specific genes with escalating CIN grade. This data supports the hypothesis that late genetic events in CIN accompany a transition to invasive cervical cancer.

Cysteine-Altering NOTCH3 Variants Are Associated with an Increased Risk of Autoimmune Diseases.

Autoimmune conditions have been reported among patients with cysteine-altering NOTCH3 variants and CADASIL. This study aimed to investigate the occurrence of autoimmune illnesses and markers of inflammation in such populations. Cases were identified who had a NOTCH3 cysteine-altering variant from the Geisinger MyCode® Community Health Initiative (MyCode®). We further performed external validation using the UK Biobank cohort. A cohort of 121 individuals with a NOTCH3 cysteine-altering variant from MyCode® was compared to a control group with no non-synonymous variation in NOTCH3 (n = 184). Medical records were evaluated for inflammatory markers and autoimmune conditions, which were grouped by the organ systems involved. A similar analysis was conducted using data from the UK Biobank (n~450,000). An overall increase in inflammatory markers among participants with a NOTCH3 cysteine-altering variant was observed when compared to an age- and sex-matched MyCode® control group (out of participants with laboratory testing: 50.9% versus 26.7%; p = 0.0047; out of total participants: 23.1% versus 10.9%; p = 0.004). Analysis of UK Biobank data indicated any autoimmune diagnosis (1.63 [1.14, 2.09], p= 2.665 × 10-3) and multiple sclerosis (3.42 [1.67, 6.02], p = 9.681 × 10-4) are associated with a NOTCH3 cysteine-altering variant in any domain. Our findings suggest a possible association between NOTCH3 cysteine-altering variants and autoimmune conditions.

Genic constraint against nonsynonymous variation across the mouse genome

BackgroundSelective constraint, the depletion of variation due to negative selection, provides insights into the functional impact of variants and disease mechanisms. However, its characterization in mice, the most commonly used mammalian model, remains limited. This study aims to quantify mouse gene constraint using a new metric called the nonsynonymous observed expected ratio (NOER) and investigate its relationship with gene function.ResultsNOER was calculated using whole-genome sequencing data from wild mouse populations (Mus musculus sp and Mus spretus). Positive correlations were observed between mouse gene constraint and the number of associated knockout phenotypes, indicating stronger constraint on pleiotropic genes. Furthermore, mouse gene constraint showed a positive correlation with the number of pathogenic variant sites in their human orthologues, supporting the relevance of mouse models in studying human disease variants.ConclusionsNOER provides a resource for assessing the fitness consequences of genetic variants in mouse genes and understanding the relationship between gene constraint and function. The study’s findings highlight the importance of pleiotropy in selective constraint and support the utility of mouse models in investigating human disease variants. Further research with larger sample sizes can refine constraint estimates in mice and enable more comprehensive comparisons of constraint between mouse and human orthologues.

Characterization of Common Genetic Variants in P2RX7 and Their Contribution to Chronic Pain Conditions

The adenosine triphosphate (ATP)-gated channel P2X7 is encoded by a gene enriched for common nonsynonymous variants. Many of these variants have functional cellular effects, and some have been implicated in chronic pain. In this study, we first systematically characterized all 17 common nonsynonymous variants using whole-cell patch clamp electrophysiology. Then, we analyzed these variants for statistical association with chronic pain phenotypes using both individual P2RX7 variants as predictors and cumulative allele counts of same-direction cellular effect in univariate models. Association and validation analyses were conducted in the Orofacial Pain: Prospective Evaluation and Risk Assessment (OPPERA) cohort (N = 3260) and in the Complex Persistent Pain Conditions (CPPC) cohort (N = 900), respectively. Our results showed an association between allele A of rs7958311 and an increased risk of chronic pelvic pain, with convergent evidence for contribution to fibromyalgia and irritable bowel syndrome, confirmed in a meta-analysis. This allelic variant produced a unique cellular phenotype: a gain-of-function in channel opening, and a loss-of-function in pore opening. A computational study using a 12-state Markov model of ATP binding to the P2X7 receptor suggested that this cellular phenotype arises from an increased ATP binding affinity and an increased open channel conductance combined with a loss of sensitization. Cumulative allele count analysis did not provide additional insights. In conclusion, our results go beyond reproducing association for rs7958311 with chronic pain and suggest that its unique combination of gain-of-function in channel and loss-of-function in pore activity may explain why it is likely the only common P2RX7 variant with contribution to chronic pain. PerspectiveThis study characterizes all common P2RX7 variants using cellular assays and statistical association analyses with chronic pain, with Markov state modeling of the most robustly associated variant.

Mitochondrial DNA heteroplasmy analysis in keratoconus patients from China.

Background: Mitochondrial DNA (mtDNA) variants have been implicated in keratoconus (KC). The present study aimed to characterize the mtDNA heteroplasmy profile in KC and explore the association of mitochondrial heteroplasmic levels with KC. Methods: Mitochondrial sequencing of peripheral blood samples and corneal tomography were conducted in 300 KC cases and 300 matched controls. The number of heteroplasmic and homoplasmic variants was calculated across the mitochondrial genome. Spearman's correlation was used to analyze the correlation between the number of heteroplasmic variants and age. The association of mtDNA heteroplasmic level with KC was analyzed by logistic regression analysis. Moreover, the relationship between mitochondrial heteroplasmic levels and clinical parameters was determined by linear regression analysis. Results: The distribution of mtDNA heteroplasmic variants showed the highest number of heteroplasmic variants in the non-coding region, while the COX3 gene exhibited the highest number in protein-coding genes. Comparisons of the number of heteroplasmic and homoplasmic non-synonymous variants in protein-coding genes revealed no significant differences between KC cases and controls (all p > 0.05). In addition, the number of heteroplasmic variants was positively associated with age in all subjects (r = 0.085, p = 0.037). The logistic regression analyses indicated that the heteroplasmic levels of m.16180_16181delAA was associated with KC (p < 0.005). Linear regression analyses demonstrated that the heteroplasmic levels of m.16180_16181delAA and m.302A>C were not correlated with thinnest corneal thickness (TCT), steep keratometry (Ks), and flat keratometry (Kf) (all p > 0.05) in KC cases and controls separately. Conclusion: The current study characterized the mtDNA heteroplasmy profile in KC, and revealed that the heteroplasmic levels of m.16180_16181delAA were associated with KC.

Genetic investigation of Nordic patients with complement-mediated kidney diseases.

Complement activation in atypical hemolytic uremic syndrome (aHUS), C3 glomerulonephropathy (C3G) and immune complex-mediated membranoproliferative glomerulonephritis (IC-MPGN) may be associated with rare genetic variants. Here we describe gene variants in the Swedish and Norwegian populations. Patients with these diagnoses (N=141) were referred for genetic screening. Sanger or next-generation sequencing were performed to identify genetic variants in 16 genes associated with these conditions. Nonsynonymous genetic variants are described when they have a minor allele frequency of <1% or were previously reported as being disease-associated. In patients with aHUS (n=94, one also had IC-MPGN) 68 different genetic variants or deletions were identified in 60 patients, of which 18 were novel. Thirty-two patients had more than one genetic variant. In patients with C3G (n=40) 29 genetic variants, deletions or duplications were identified in 15 patients, of which 9 were novel. Eight patients had more than one variant. In patients with IC-MPGN (n=7) five genetic variants were identified in five patients. Factor H variants were the most frequent in aHUS and C3 variants in C3G. Seventeen variants occurred in more than one condition. Genetic screening of patients with aHUS, C3G and IC-MPGN is of paramount importance for diagnostics and treatment. In this study, we describe genetic assessment of Nordic patients in which 26 novel variants were found.

The identification and genetic characteristics of Quang Binh virus from field-captured Culex tritaeniorhynchus (Diptera: Culicidae) from Guizhou Province, China

BackgroundMosquitoes carry a variety of viruses that can cause disease in humans, animals and livestock. Surveys for viruses carried by wild mosquitoes can significantly contribute to surveillance efforts and early detection systems. In addition to mosquito-borne viruses, mosquitoes harbor many insect-specific viruses (ISVs). Quang Binh virus (QBV) is one such example, categorized as an ISV within the Flavivirus genus (family Flaviviridae). QBV has been specifically documented in Vietnam and China, with reports limited to several mosquito species.MethodsThe homogenate obtained from female mosquitoes was cultured on C6/36 (Aedes albopictus) and BHK-21 (baby hamster kidney) cell lines. Positive cultures were identified by reverse transcription-polymerase chain reaction (RT‒PCR) with taxon- or species-specific primers. Next-generation sequencing was employed to sequence the complete genomes of the identified positive samples. Subsequently, phylogenetic, gene homology, molecular evolutionary and genetic variation analyses were conducted.ResultIn 2021, a total of 32,177 adult female mosquitoes were collected from 15 counties in Guizhou Province, China. The predominant mosquito species identified were Culex tritaeniorhynchus, Armigeres subalbatus and Anopheles sinensis. Among the collected mosquitoes, three positive cultures were obtained from Cx. tritaeniorhynchus pools, revealing the presence of Quang Binh virus (QBV) RNA sequences. Phylogenetic analysis indicated that the three Guizhou isolates, along with the prototype isolate from Vietnam, formed distinct branches. These branches were primarily closely related to other QBV isolates reported in China. Comparative analysis revealed a high degree of nucleotide and amino acid homology between the Guizhou isolates and both Vietnamese and other indigenous Chinese isolates. Additionally, nonsynonymous single-nucleotide variants (SNVs) were observed in these strains compared to the QBV prototype strain.ConclusionThis study represents the first report of QBV presences in Cx. tritaeniorhynchus mosquitoes in Guizhou Province, China. Phylogenetic tree analysis showed that the three Guizhou isolates were most closely related to the QBV genes found in China. In addition, the study of the genetic characteristics and variation of this virus provided a deeper understanding of QBV and enriched the baseline data of these insect-specific flaviviruses (ISFVs).Graphical

An in-silico study to determine susceptibility to cancer by evaluating the coding and non-coding non-synonymous single nucleotide variants in the SOCS3 gene

Single Nucleotide Variant (SNVs) affect gene expression as well as protein structure and activity, leading to reduced signaling capabilities and ultimately, increasing cancer risk. SOCS3 (suppressor of cytokine signaling 3), a critical tumor suppressor providing a substantial part in the feedback loop of the JAK/STAT pathway, is abnormally suppressed in various cancer. This study aims to screen non-coding and potentially deleterious coding SNVs in the SOCS3 gene. We performed six programs: PredictSNP1.0 (predicting Deleterious nsSNVs), ConSurf (analyzing sequence conservation), ModPred (analyzing SNVS in PTMs sites), I-Mutant and MUpro (to analyze SNVs effecting protein stability), and molecular docking and molecular dynamics (MD) (to assess the consequences of SOCS3 genetic variations on JAK interactions) for coding regions and three programs (UTRSite, SNP2TFBS, miRNA SNP) (to analyze SNVs effecting the gene expression) in non-coding regions, respectively. After screening 2786 SOCS3 SNVs, we found 10 SNVs, as well as 49 SNPs that change the function of non-coding areas. Out of 10 selected nsSNVs, 3 SNVs (W48R, R71C, N198S) predicted to be the most damaging by all the software programs, as well as one nsSNV (R194W) could be highly deleterious from Molecular Docking analysis combined with MD Simulations. Our findings propose a procedure for studying the structure-related consequences of SNVs on protein function in the future. Communicated by Ramaswamy H. Sarma