Multi-omics approach dissects cis-regulatory mechanisms underlying North Carolina macular dystrophy, a retinal enhanceropathy.

North Carolina macular dystrophy (NCMD) is a rare autosomal-dominant disease affecting macular development. The disease is caused by non-coding single-nucleotide variants (SNVs) in two hotspot regions near PRDM13 and by duplications in two distinct chromosomal loci, overlapping DNase I hypersensitive sites near either PRDM13 or IRX1. To unravel the mechanisms by which these variants cause disease, we first established a genome-wide multi-omics retinal database, RegRet. Integration of UMI-4C profiles we generated on adult human retina then allowed fine-mapping of the interactions of the PRDM13 and IRX1 promoters and the identification of eighteen candidate cis-regulatory elements (cCREs), the activity of which was investigated by luciferase and Xenopus enhancer assays. Next, luciferase assays showed that the non-coding SNVs located in the two hotspot regions of PRDM13 affect cCRE activity, including two NCMD-associated non-coding SNVs that we identified herein. Interestingly, the cCRE containing one of these SNVs was shown to interact with the PRDM13 promoter, demonstrated in vivo activity in Xenopus, and is active at the developmental stage when progenitor cells of the central retina exit mitosis, suggesting that this region is a PRDM13 enhancer. Finally, mining of single-cell transcriptional data of embryonic and adult retina revealed the highest expression of PRDM13 and IRX1 when amacrine cells start to synapse with retinal ganglion cells, supporting the hypothesis that altered PRDM13 or IRX1 expression impairs interactions between these cells during retinogenesis. Overall, this study provides insight into the cis-regulatory mechanisms of NCMD and supports that this condition is a retinal enhanceropathy.

Combining full-length gene assay and SpliceAI to interpret the splicing impact of all possible SPINK1 coding variants.

BACKGROUND: Single-nucleotide variants (SNVs) within gene coding sequences can significantly impact pre-mRNA splicing, bearing profound implications for pathogenic mechanisms and precision medicine. In this study, we aim to harness the well-established full-length gene splicing assay (FLGSA) in conjunction with SpliceAI to prospectively interpret the splicing effects of all potential coding SNVs within the four-exon SPINK1 gene, a gene associated with chronic pancreatitis. RESULTS: Our study began with a retrospective analysis of 27 SPINK1 coding SNVs previously assessed using FLGSA, proceeded with a prospective analysis of 35 new FLGSA-tested SPINK1 coding SNVs, followed by data extrapolation, and ended with further validation. In total, we analyzed 67 SPINK1 coding SNVs, which account for 9.3% of the 720 possible coding SNVs. Among these 67 FLGSA-analyzed SNVs, 12 were found to impact splicing. Through detailed comparison of FLGSA results and SpliceAI predictions, we inferred that the remaining 653 untested coding SNVs in the SPINK1 gene are unlikely to significantly affect splicing. Of the 12 splice-altering events, nine produced both normally spliced and aberrantly spliced transcripts, while the remaining three only generated aberrantly spliced transcripts. These splice-impacting SNVs were found solely in exons 1 and 2, notably at the first and/or last coding nucleotides of these exons. Among the 12 splice-altering events, 11 were missense variants (2.17% of 506 potential missense variants), and one was synonymous (0.61% of 164 potential synonymous variants). Notably, adjusting the SpliceAI cut-off to 0.30 instead of the conventional 0.20 would improve specificity without reducing sensitivity. CONCLUSIONS: By integrating FLGSA with SpliceAI, we have determined that less than 2% (1.67%) of all possible coding SNVs in SPINK1 significantly influence splicing outcomes. Our findings emphasize the critical importance of conducting splicing analysis within the broader genomic sequence context of the study gene and highlight the inherent uncertainties associated with intermediate SpliceAI scores (0.20 to 0.80). This study contributes to the field by being the first to prospectively interpret all potential coding SNVs in a disease-associated gene with a high degree of accuracy, representing a meaningful attempt at shifting from retrospective to prospective variant analysis in the era of exome and genome sequencing.

MRAS in coronary artery disease-Unchartered territory.

Genome-wide association studies (GWAS) have identified coronary artery disease (CAD) susceptibility locus on chromosome 3q22.3. This locus contains a cluster of several genes that includes muscle rat sarcoma virus (MRAS). Common MRAS variants are also associated with CAD causing risk factors such as hypertension, dyslipidemia, obesity, and type II diabetes. The MRAS gene is an oncogene that encodes a membrane-bound small GTPase. It is involved in a variety of signaling pathways, regulating cell differentiation and cell survival (mitogen-activated protein kinase [MAPK]/extracellular signal-regulated kinase and phosphatidylinositol 3-kinase) as well as acute phase response signaling (tumor necrosis factor [TNF] and interleukin 6 [IL6] signaling). In this review, we will summarize the role of genetic MRAS variants in the etiology of CAD and its comorbidities with the focus on tissue distribution of MRAS isoforms, cell type/tissue specificity, and mode of action of single nucleotide variants in MRAS associated complex traits. Finally, we postulate that CAD risk variants in the MRAS locus are specific to smooth muscle cells and lead to higher levels of MRAS, particularly in arterial and cardiac tissue, resulting in MAPK-dependent tissue hypertrophy or hyperplasia.

Selected cachaça yeast strains share a genomic profile related to traits relevant to industrial fermentation processes.

The isolation and selection of yeast strains to improve the quality of the cachaça-Brazilian Spirit-have been studied in our research group. Our strategy considers Saccharomyces cerevisiae as the predominant species involved in sugarcane juice fermentation and the presence of different stressors (osmolarity, temperature, ethanol content, and competition with other microorganisms). It also considers producing balanced concentrations of volatile compounds (higher alcohols and acetate and/or ethyl esters), flocculation capacity, and ethanol production. Since the genetic bases behind these traits of interest are not fully established, the whole genome sequencing of 11 different Saccharomyces cerevisiae strains isolated and selected from different places was analyzed to identify the presence of a specific genetic variation common to cachaça yeast strains. We have identified 20,128 single-nucleotide variants shared by all genomes. Of these shared variants, 37 were new variants (being six missenses), and 4,451 were identified as missenses. We performed a detailed functional annotation (using enrichment analysis, protein-protein interaction network analysis, and database and in-depth literature searches) of these new and missense variants. Many genes carrying these variations were involved in the phenotypes of flocculation, tolerance to fermentative stresses, and production of volatile compounds and ethanol. These results demonstrate the existence of a genetic profile shared by the 11 strains under study that could be associated with the applied selective strategy. Thus, this study points out genes and variants that may be used as molecular markers for selecting strains well suited to the fermentation process, including genetic improvement by genome editing, ultimately producing high-quality beverages and adding value.IMPORTANCEThis work demonstrates the existence of new genetic markers related to different phenotypes used to select yeast strains and mutations in genes directly involved in producing flavoring compounds and ethanol, and others related to flocculation and stress resistance.

Genetic characteristics involved in COVID-19 severity. The CARGENCORS case-control study and meta-analysis.

Pre-existing coronary artery disease (CAD), and thrombotic, inflammatory, or virus infectivity response phenomena have been associated with COVID-19 disease severity. However, the association of candidate single nucleotide variants (SNVs) related to mechanisms of COVID-19 complications has been seldom analysed. Our aim was to test and validate the effect of candidate SNVs on COVID-19 severity. CARGENCORS (CARdiovascular GENetic risk score for Risk Stratification of patients positive for SARS-CoV-2 [COVID-19] virus) is an age- and sex-matched case-control study with 818 COVID-19 cases hospitalized with hypoxemia, and 1636 controls with COVID-19 treated at home. The association between severity and SNVs related to CAD (n = 32), inflammation (n = 19), thrombosis (n = 14), virus infectivity (n = 11), and two published to be related to COVID-19 severity was tested with adjusted logistic regression models. Two external independent cohorts were used for meta-analysis (SCOURGE and UK Biobank). After adjustment for potential confounders, 14 new SNVs were associated with COVID-19 severity in the CARGENCORS Study. These SNVs were related to CAD (n = 10), thrombosis (n = 2), and inflammation (n = 2). We also confirmed eight SNVs previously related to severe COVID-19 and virus infectivity. The meta-analysis showed five SNVs associated with severe COVID-19 in adjusted analyses (rs11385942, rs1561198, rs6632704, rs6629110, and rs12329760). We identified 14 novel SNVs and confirmed eight previously related to COVID-19 severity in the CARGENCORS data. In the meta-analysis, five SNVs were significantly associated to COVID-19 severity, one of them previously related to CAD.

Mitochondrial genome variants associated with amyotrophic lateral sclerosis and their haplogroup distribution.

INTRODUCTION/AIMS: Amyotrophic lateral sclerosis (ALS) may be familial or sporadic, and twin studies have revealed that even sporadic forms have a significant genetic component. Variants in 55 nuclear genes have been associated with ALS and although mitochondrial dysfunction is observed in ALS, variants in mitochondrial genomes (mitogenomes) have not yet been tested for association with ALS. The aim of this study was to determine whether mitogenome variants are associated with ALS. METHODS: We conducted a genome-wide association study (GWAS) in mitogenomes of 1965 ALS patients and 2547 controls. RESULTS: We identified 51 mitogenome variants with p values <10-7, of which 13 had odds ratios (ORs) >1, in genes RNR1, ND1, CO1, CO3, ND5, ND6, and CYB, while 38 variants had OR <1 in genes RNR1, RNA2, ND1, ND2, CO2, ATP8, ATP6, CO3, ND3, ND4, ND5, ND6, and CYB. The frequencies of haplogroups H, U, and L, the most frequent in our ALS data set, were the same in different onset sites (bulbar, limb, spinal, and axial). Also, intra-haplogroup GWAS revealed unique ALS-associated variants in haplogroups L and U. DISCUSSION: Our study shows that mitogenome single nucleotide variants (SNVs) are associated with ALS and suggests that these SNVs could be included in routine genetic testing for ALS and that mitochondrial replacement therapy has the potential to serve as a basis for ALS treatment.

Obtaining the best igRNAs for bystander-less correction of all ABE-reversible pathogenic SNVs using high-throughput screening.

Imperfect -gRNA (igRNA) provides a simple strategy for single-base editing of a base editor. However, a significant number of igRNAs need to be generated and tested for each target locus to achieve efficient single-base reversion of pathogenic single nucleotide variations (SNVs), which hinders the direct application of this technology. To provide ready-to-use igRNAs for single-base and bystander-less correction of all the adenine base editor (ABE)-reversible pathogenic SNVs, we employed a high-throughput method to edit all 5,253 known ABE-reversible pathogenic SNVs, each with multiple systematically designed igRNAs, and two libraries of 96,000 igRNAs were tested. A total of 1,988 SNV loci could be single-base reversed by igRNA with a >30% efficiency. Among these 1,988 loci, 378 SNV loci exhibited an efficiency of more than 90%. At the same time, the bystander editing efficiency of 76.62% of the SNV loci was reduced to 0%, while remaining below 1% for another 18.93% of the loci. These ready-to-use igRNAs provided the best solutions for a substantial portion of the 4,657 pathogenic/likely pathogenic SNVs. In this work, we overcame one of the most significant obstacles of base editors and provide a ready-to-use platform for the genetic treatment of diseases caused by ABE-reversible SNVs.

Assembly of In-Situ Gel Containing Nano-Spanlastics of an Angiotensin II Inhibitor as a Novel Epitome for Hypertension Management: Factorial Design Optimization, In-vitro Gauging, Pharmacokinetics, and Pharmacodynamics Appraisal.

More than 1 billion people worldwide suffer from hypertension; therefore, hypertension management has been categorized as a global health priority. Losartan potassium (LP) is an antihypertensive drug with a limited oral bioavailability of about 33% since it undergoes the initial metabolic cycle. Thus, nasal administration is a unique route to overcome first-pass metabolism. The investigation focused on the potential effects of LP-loaded spanlastic vesicles (SNVs) on LP pharmacodynamics and pharmacokinetic parameters, utilizing a thin-film hydration methodology established on a 3122 full factorial design. Entrapment efficiency (EE%) ranged from 39.8 ± 3.87.8 to 83.8 ± 2.92% for LP-SNVs. Vesicle size (VS) varied from 205.5 ± 6.5.10 to 445.1 ± 13.52 nm, and the percentage of LP released after 8 h (Q8h) ranged from 30.8 ± 3.10 to 68.8 ± 1.45%. LP permeated through the nasal mucosa during 24 h and flocculated from 194.1 ± 4.90 to 435.3 ± 13.53 µg/cm2. After twenty-four hours, the optimal LP-SNVs in-situ gel showed 2.35 times more permeation through the nasal mucosa than the LP solution. It also lowered systolic blood pressure, so it is thought to be better than the reference formulation in terms of pharmacodynamics. The pharmacokinetics studies demonstrated that the intranasal LP-SNVs gel boosted its bioavailability approximately 6.36 times compared to the oral LP solution. Our research showed that intranasal LP-SNVs could be a good nanoplatform because they are well-tolerated and have possible pharmacokinetics and pharmacodynamics.

CVID-Associated B Cell Activating Factor Receptor Variants Change Receptor Oligomerization, Ligand Binding, and Signaling Responses.

PURPOSE: Binding of the B cell activating factor (BAFF) to its receptor (BAFFR) activates in mature B cells many essential pro-survival functions. Null mutations in the BAFFR gene result in complete BAFFR deficiency and cause a block in B cell development at the transition from immature to mature B cells leading therefore to B lymphopenia and hypogammaglobulinemia. In addition to complete BAFFR deficiency, single nucleotide variants encoding BAFFR missense mutations were found in patients suffering from common variable immunodeficiency (CVID), autoimmunity, or B cell lymphomas. As it remained unclear to which extent such variants disturb the activity of BAFFR, we performed genetic association studies and developed a cellular system that allows the unbiased analysis of BAFFR variants regarding oligomerization, signaling, and ectodomain shedding. METHODS: In addition to genetic association studies, the BAFFR variants P21R, A52T, G64V, DUP92-95, P146S, and H159Y were expressed by lentiviral gene transfer in DG-75 Burkitt's lymphoma cells and analyzed for their impacts on BAFFR function. RESULTS: Binding of BAFF to BAFFR was affected by P21R and A52T. Spontaneous oligomerization of BAFFR was disturbed by P21R, A52T, G64V, and P146S. BAFF-dependent activation of NF-κB2 was reduced by P21R and P146S, while interactions between BAFFR and the B cell antigen receptor component CD79B and AKT phosphorylation were impaired by P21R, A52T, G64V, and DUP92-95. P21R, G64V, and DUP92-95 interfered with phosphorylation of ERK1/2, while BAFF-induced shedding of the BAFFR ectodomain was only impaired by P21R. CONCLUSION: Although all variants change BAFFR function and have the potential to contribute as modifiers to the development of primary antibody deficiencies, autoimmunity, and lymphoma, P21R is the only variant that was found to correlate positively with CVID.

Two single nucleotide variants in the miR-23a promoter affect granulosa cell apoptosis.

microRNAs (miRNAs) are well known to be important in mammalian female fertility. However, the genetic regulation of miRNAs associated with female fertility remains largely unknown. Here, we report that two single-nucleotide variants (SNVs) in the miR-23a promoter strongly influence miR-23a transcription and function in granulosa cell (GC) apoptosis. Two novel SNVs, g.-283G>C and g.-271C>T, were detected in the porcine miR-23a promoter by pooled-DNA sequencing. Furthermore, SNVs in the promoter region influenced miR-23a transcription in porcine GCs by altering its promoter activity. Functionally, SNVs in the promoter strongly influenced miR-23a regulation of early apoptosis in porcine GCs cultured in vitro. In addition, a preliminary association analysis showed that the combined genotypes of the two SNVs, rather than a single SNV, were tentatively associated with sow fertility traits in a Large White population. Overall, our findings suggest that the SNVs g.-283G>C and g.-271C>T in the miR-23a promoter are causal variants affecting GC apoptosis and miR-23a may be a potential small-molecule nonhormonal drug for regulating female fertility.

Influence of Genetic Variations in miRNA and Genes Encoding Proteins in the miRNA Synthesis Complex on Toxicity of the Treatment of Pediatric B-Cell ALL in the Brazilian Amazon.

Acute lymphoblastic leukemia (ALL) is the most common childhood cancer in the world. Single nucleotide variants (SNVs) in miRNA and genes encoding proteins of the miRNA synthesis complex (SC) may affect the processing of drugs used in the treatment of ALL, resulting in treatment-related toxicities (TRTs). We investigated the role of 25 SNVs in microRNA genes and genes encoding proteins of the miRNA SC, in 77 patients treated for ALL-B from the Brazilian Amazon. The 25 SNVs were investigated using the TaqMan® OpenArray™ Genotyping System. SNVs rs2292832 (MIR149), rs2043556 (MIR605), and rs10505168 (MIR2053) were associated with an increased risk of developing Neurological Toxicity, while rs2505901 (MIR938) was associated with protection from this toxicity. MIR2053 (rs10505168) and MIR323B (rs56103835) were associated with protection from gastrointestinal toxicity, while DROSHA (rs639174) increased the risk of development. The rs2043556 (MIR605) variant was related to protection from infectious toxicity. SNVs rs12904 (MIR200C), rs3746444 (MIR499A), and rs10739971 (MIRLET7A1) were associated with a lower risk for severe hematologic toxicity during ALL treatment. These findings reveal the potential for the use of these genetic variants to understand the development of toxicities related to the treatment of ALL in patients from the Brazilian Amazon region.

Molecular Landscape of Pelvic Organ Prolapse Provides Insights into Disease Etiology.

Pelvic organ prolapse (POP) represents a major health care burden in women, but its underlying pathophysiological mechanisms have not been elucidated. We first used a case-control design to perform an exome chip study in 526 women with POP and 960 control women to identify single nucleotide variants (SNVs) associated with the disease. We then integrated the functional interactions between the POP candidate proteins derived from the exome chip study and other POP candidate molecules into a molecular landscape. We found significant associations between POP and SNVs in 54 genes. The proteins encoded by 26 of these genes fit into the molecular landscape, together with 43 other POP candidate molecules. The POP landscape is located in and around epithelial cells and fibroblasts of the urogenital tract and harbors four interacting biological processes-epithelial-mesenchymal transition, immune response, modulation of the extracellular matrix, and fibroblast function-that are regulated by sex hormones and TGFB1. Our findings were corroborated by enrichment analyses of differential gene expression data from an independent POP cohort. Lastly, based on the landscape and using vaginal fibroblasts from women with POP, we predicted and showed that metformin alters gene expression in these fibroblasts in a beneficial direction. In conclusion, our integrated molecular landscape of POP provides insights into the biological processes underlying the disease and clues towards novel treatments.

Detection of Nonsynonymous Single Variants in Human HLA-DRB1 Exon 2 Associated with Renal Transplant Rejection.

Background: HLA-DRB1 is the most polymorphic gene in the human leukocyte antigen (HLA) class II, and exon 2 is critical because it encodes antigen-binding sites. This study aimed to detect functional or marker genetic variants of HLA-DRB1 exon 2 in renal transplant recipients (acceptance and rejection) using Sanger sequencing. Methods: This hospital-based case-control study collected samples from two hospitals over seven months. The 60 participants were equally divided into three groups: rejection, acceptance, and control. The target regions were amplified and sequenced by PCR and Sanger sequencing. Several bioinformatics tools have been used to assess the impact of non-synonymous single-nucleotide variants (nsSNVs) on protein function and structure. The sequences data that support the findings of this study with accession numbers (OQ747803-OQ747862) are available in National Center for Biotechnology Information (GenBank database). Results: Seven SNVs were identified, two of which were novel (chr6(GRCh38.p12): 32584356C>A (K41N) and 32584113C>A (R122R)). Three of the seven SNVs were non-synonymous and found in the rejection group (chr6(GRCh38.p12): 32584356C>A (K41N), 32584304A>G (Y59H), and 32584152T>A (R109S)). The nsSNVs had varying effects on protein function, structure, and physicochemical parameters and could play a role in renal transplant rejection. The chr6(GRCh38.p12):32584152T>A variant showed the greatest impact. This is because of its conserved nature, main domain location, and pathogenic effects on protein structure, function, and stability. Finally, no significant markers were identified in the acceptance samples. Conclusion: Pathogenic variants can affect intramolecular/intermolecular interactions of amino acid residues, protein function/structure, and disease risk. HLA typing based on functional SNVs could be a comprehensive, accurate, and low-cost method for covering all HLA genes while shedding light on previously unknown causes in many graft rejection cases.

Development and evaluation of polygenic risk scores for prediction of endometrial cancer risk in European women.

PURPOSE: Single-nucleotide variations (SNVs) (formerly single-nucleotide polymorphism [SNV]) influence genetic predisposition to endometrial cancer. We hypothesized that a polygenic risk score (PRS) comprising multiple SNVs may improve endometrial cancer risk prediction for targeted screening and prevention. METHODS: We developed PRSs from SNVs identified from a systematic review of published studies and suggestive SNVs from the Endometrial Cancer Association Consortium. These were tested in an independent study of 555 surgically-confirmed endometrial cancer cases and 1202 geographically-matched controls from Manchester, United Kingdom and validated in 1676 cases and 116,960 controls from the UK Biobank (UKBB). RESULTS: Age and body mass index predicted endometrial cancer in both data sets (Manchester: area under the receiver operator curve [AUC] = 0.77, 95% CI = 0.74-0.80; UKBB: AUC = 0.74, 95% CI = 0.73-0.75). The AUC for PRS19, PRS24, and PRS72 were 0.58, 0.55, and 0.57 in the Manchester study and 0.56, 0.54, and 0.54 in UKBB, respectively. For PRS19, women in the third tertile had a 2.1-fold increased risk of endometrial cancer compared with those in the first tertile of the Manchester study (odds ratio = 2.08, 95% CI = 1.61-2.68, Ptrend = 5.75E-9). Combining PRS19 with age and body mass index improved discriminatory power (Manchester study: AUC = 0.79, 95% CI = 0.76-0.82; UKBB: AUC =0.75, 95% CI = 0.73-0.76). CONCLUSION: An endometrial cancer risk prediction model incorporating a PRS derived from multiple SNVs may help stratify women for screening and prevention strategies.

Single-Nucleotide Variants in microRNAs Sequences or in their Target Genes Might Influence the Risk of Epilepsy: A Review.

Single-nucleotide variant (SNV) is a single base mutation at a specific location in the genome and may play an import role in epilepsy pathophysiology. The aim of this study was to review case-control studies that have investigated the relationship between SNVs within microRNAs (miRs) sequences or in their target genes and epilepsy susceptibility from January 1, 2010 to October 31, 2020. Nine case-control studies were included in the present review. The mainly observed SNVs associated with drug-resistant epilepsy (DRE) risk were SNVs n.60G > C (rs2910164) and n.-411A > G (rs57095329), both located at miR-146a mature sequence and promoter region, respectively. In addition, the CC haplotype (rs987195-rs969885) and the AA genotype at rs4817027 in the MIR155HG/miR-155 tagSNV were also genetic susceptibility markers for early-onset epilepsy. MiR-146a has been observed as upregulated in human astrocytes in epileptogenesis and it regulates inflammatory process through NF-κB signaling by targeting tumor necrosis factor-associated factor 6 (TRAF6) gene. The SNVs rs2910164 and rs57095329 may modify the expression level of mature miR-146a and the risk for epilepsy and SNVs located at rs987195-rs969885 haplotype and at rs4817027 in the MIR155HG/miR-155 tagSNV could interfere in the miR-155 expression modulating inflammatory pathway genes involved in the development of early-onset epilepsy. In addition, SNVs rs662702, rs3208684, and rs35163679 at 3'untranslated region impairs the ability of miR-328, let-7b, and miR-200c binding affinity with paired box protein PAX-6 (PAX6), BCL2 like 1 (BCL2L1), and DNA methyltransferase 3 alpha (DNMT3A) target genes. The SNV rs57095329 might be correlated with DRE when a larger number of patients are evaluated. Thus, we concluded that the main drawback of most of studies is the small number of individuals enrolled, which lacks sample power.

Candidate Genes in Testing Strategies for Linkage Analysis and Bioinformatic Sorting of Whole Genome Sequencing Data in Three Small Japanese Families with Idiopathic Superior Oblique Muscle Palsy.

Idiopathic superior oblique muscle palsy is a major type of paralytic, non-comitant strabismus and presents vertical and cyclo-torsional deviation of one eye against the other eye, with a large vertical fusion range and abnormal head posture such as head tilt. Genetic background is considered to play a role in its development, as patients with idiopathic superior oblique muscle palsy have varying degrees of muscle hypoplasia and, rarely, the complete absence of the muscle, that is, aplasia. In this study, whole genome sequencing was performed, and single nucleotide variations and short insertions/deletions (SNVs/InDels) were annotated in two patients each in three small families (six patients in total) with idiopathic superior oblique muscle palsy, in addition to three normal individuals in one family. At first, linkage analysis was carried out in the three families and SNVs/InDels in chromosomal loci with negative LOD scores were excluded. Next, SNVs/InDels shared by the six patients, but not by the three normal individuals, were chosen. SNVs/InDels were further narrowed down by choosing low-frequency (<1%) or non-registered SNVs/InDels in four databases for the Japanese population, and then by choosing SNVs/InDels with functional influence, leading to one candidate gene, SSTR5-AS1 in chromosome 16. The six patients were heterozygous for 13-nucleotide deletion in SSTR5-AS1, except for one homozygous patient, while the three normal individuals were wild type. Targeted polymerase chain reaction (PCR) and direct sequencing of PCR products confirmed the 13-nucleotide deletion in SSTR5-AS1. In the face of newly-registered SSTR5-AS1 13-nucleotide deletion at a higher frequency in a latest released database for the Japanese population, the skipping of low-frequency and non-registration sorting still resulted in only 13 candidate genes including SSTR5-AS1 as common variants. The skipping of linkage analysis also led to the same set of 13 candidate genes. Different testing strategies that consisted of linkage analysis and simple unintentional bioinformatics could reach candidate genes in three small families with idiopathic superior oblique muscle palsy.

FiNGS: high quality somatic mutations using filters for next generation sequencing.

BACKGROUND: Somatic variant callers are used to find mutations in sequencing data from cancer samples. They are very sensitive and have high recall, but also may produce low precision data with a large proportion of false positives. Further ad hoc filtering is commonly performed after variant calling and before further analysis. Improving the filtering of somatic variants in a reproducible way represents an unmet need. We have developed Filters for Next Generation Sequencing (FiNGS), software written specifically to address these filtering issues. RESULTS: Developed and tested using publicly available sequencing data sets, we demonstrate that FiNGS reliably improves upon the precision of default variant caller outputs and performs better than other tools designed for the same task. CONCLUSIONS: FiNGS provides researchers with a tool to reproducibly filter somatic variants that is simple to both deploy and use, with filters and thresholds that are fully configurable by the user. It ingests and emits standard variant call format (VCF) files and will slot into existing sequencing pipelines. It allows users to develop and implement their own filtering strategies and simple sharing of these with others.

Evolution of SARS-CoV-2 genome from December 2019 to late March 2020: Emerged haplotypes and informative Tag nucleotide variations.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes serious disease in humans. First identified in November/December 2019 in China, it has rapidly spread worldwide. We analyzed 2790 SARS-CoV-2 genome sequences from 56 countries that were available on April 2, 2020, to assess the evolution of the virus during this early phase of its expansion. We aimed to assess sequence variations that had evolved in virus genomes, giving the greatest attention to the S gene. We also aimed to identify haplotypes that the variations may define and consider their geographic and chronologic distribution. Variations at 1930 positions that together cause 1203 amino acid changes were identified. The frequencies of changes normalized to the lengths of genes and encoded proteins were relatively high in ORF3a and relatively low in M. A variation that causes an Asp614Gly near the receptor-binding domain of S were found at a high frequency, and it was considered that this may contribute to the rapid spread of viruses with this variation. Our most important findings relate to haplotypes. Sixty-six haplotypes that constitute thirteen haplotype groups (H1-H13) were identified, and 84.6% of the 2790 sequences analyzed were associated with these haplotypes. The majority of the sequences (75.1%) were associated with haplotype groups H1-H3. The distribution pattern of the haplotype groups differed in various geographic regions. A few were country/territory specific. The location and time of emergence of some haplotypes are discussed. Importantly, nucleotide variations that define the various haplotypes and Tag/signature variations for most of the haplotypes are reported. The practical applications of these variations are discussed.

Single nucleotide variants in immune-response genes and the tumor microenvironment composition predict progression of mantle cell lymphoma.

BACKGROUND: There is evidence to consider that the tumor microenvironment (TME) composition associates with antitumor immune response, and may predict the outcome of various non-Hodgkin lymphoma subtypes. However, in the case of mantle cell lymphoma (MCL), a rare and aggressive disease, there is lacking a detailed study of the TME components, as well as an integrative approach among them in patients' samples. Also, from the genetic point of view, it is known that single nucleotide variants (SNVs) in immune-response genes are among important regulators of immunity. At present, it is uncertain whether SNVs in candidate immune-response genes and the TME composition are able to alter the prognosis in MCL. METHODS: We assessed a detailed TME composition in 88 MCL biopsies using immunohistochemistry, which was automatically analyzed by pixel counting (Aperio system). We also genotyped SNVs located in candidate immune-response genes (IL12A, IL2, IL10, TGFB1, TGFBR1, TGFBR2, IL17A, IL17F) in 95 MCL patients. We tested whether the SNVs could modulate the respective protein expression and TME composition in the tumor compartment. Finally, we proposed survival models in rituximab-treated patients, considering immunohistochemical and SNV models. RESULTS: High FOXP3/CD3 ratios (p = 0.001), high IL17A levels (p = 0.003) and low IL2 levels (p = 0.03) were individual immunohistochemical predictors of poorer survival. A principal component, comprising high quantities of macrophages and high Ki-67 index, also worsened outcome (p = 0.02). In the SNV model, the CC haplotype of IL10 (p < 0.01), the GG genotype of IL2 rs2069762 (p = 0.02) and the AA+AG genotypes of TGFBR2 rs3087465 (p < 0.01) were independent predictors of outcome. Finally, the GG genotype of TGFB1 rs6957 associated with lower tumor TGFß levels (p = 0.03) and less CD163+ macrophages (p = 0.01), but did not modulate patients' survival. CONCLUSIONS: Our results indicate that the TME composition has relevant biological roles in MCL. In this setting, immunohistochemical detection of T-reg cells, IL17A and IL2, coupled with SNV genotyping in IL10, TGFBR2 and IL2, may represent novel prognostic factors in this disease, following future validations.

Novel MYT1 variants expose the complexity of oculo-auriculo-vertebral spectrum genetic mechanisms.

Oculo-auriculo-vertebral spectrum (OAVS) is a developmental disorder characterized by anomalies mainly involving the structures derived from the first and second pharyngeal arches. The spectrum presents with heterogeneous clinical features and complex etiology with genetic factors not yet completely understood. To date, MYT1 is the most important gene unambiguously associated with the spectrum and with functional data confirmation. In this work, we aimed to identify new single nucleotide variants (SNVs) affecting MYT1 in a cohort of 73 Brazilian patients diagnosed with OAVS. In addition, we investigated copy number variations (CNVs) encompassing this gene or its cis-regulatory elements and compared the frequency of these events in patients versus a cohort of 455 Brazilian control individuals. A new SNV, predicted as likely deleterious, was identified in five unrelated patients with OAVS. All five patients presented hearing impairment and orbital asymmetry suggesting an association with the variant. CNVs near MYT1, located in its neighboring topologically associating domain (TAD), were found to be enriched in patients when compared to controls, indicating a possible involvement of this region with OAVS pathogenicity. Our findings highlight the genetic complexity of the spectrum that seems to involve more than one variant type and inheritance patterns.