Search

Exploring non-coding regions is increasingly gaining importance in the diagnosis of inherited retinal dystrophies. Deep-intronic variants causing aberrant splicing have been identified, prompting the development of antisense oligonucleotides (ASOs) to modulate splicing. We performed a screening of five previously described USH2A deep-intronic variants among USH2A monoallelic patients with Usher syndrome (USH) or isolated retinitis pigmentosa. Sequencing of entire USH2A or USH genes was then conducted in unresolved or newly monoallelic cases. The splicing impact of identified variants was assessed using minigene assays, and ASOs were designed to correct splicing. The screening allowed to diagnose 30.95% of the studied patients. The sequencing of USH genes revealed 16 new variants predicted to affect splicing, with four confirmed to affect splicing through minigene assays. Two of them were unreported deep-intronic variants and predicted to include a pseudoexon in the pre-mRNA, and the other two could alter a regulatory cis-element. ASOs designed for three USH2A deep-intronic variants successfully redirected splicing in vitro. Our study demonstrates the improvement in genetic characterization of IRDs when analyzing non-coding regions, highlighting that deep-intronic variants significantly contribute to USH2A pathogenicity. Furthermore, successful splicing modulation through ASOs highlights their therapeutic potential for patients carrying deep-intronic variants.

IntroductionLong-read sequencing (LRS) enables accurate structural variant detection and variant phasing. When a molecular diagnosis is suspected, target enrichment can reduce the cost and duration of sequencing.MethodsLRS was conducted in five inherited retinal dystrophy (IRD) patients harboring a monoallelic variant in RPE65 that remained uncharacterized after clinical exome sequencing (CES). CRISPR-Cas9 guide RNA probes were designed to target a 31 kb region, including the entire RPE65 locus. The DNA was sequenced on a MinION platform. Short-read ×30 whole-genome sequencing (WGS) was performed for five patients to validate nanopore results.ResultsThe nanopore sequencing process yielded a median of 271 reads within the targeted region, with a mean depth of 109 and a median read size of 8 kb. All variants identified by CES have been detected using this approach, and no additional RPE65 gene causative variants were found. Nanopore variant detection demonstrated performance akin to short-read WGS at similar coverage levels, although exhibiting increased false positive calls at lower coverage.DiscussionIn this study, we explore the advantages of using a targeted approach together with long-read sequencing to identify variants associated with IRD. The results underscore the utility of targeted long reads for characterizing patients affected by rare diseases when first-tier diagnostic tests are non-conclusive.

The genetic basis of severe COVID-19 has been thoroughly studied, and many genetic risk factors shared between populations have been identified. However, reduced sample sizes from non-European groups have limited the discovery of population-specific common risk loci. In this second study nested in the SCOURGE consortium, we conducted a genome-wide association study (GWAS) for COVID-19 hospitalization in admixed Americans, comprising a total of 4702 hospitalized cases recruited by SCOURGE and seven other participating studies in the COVID-19 Host Genetic Initiative. We identified four genome-wide significant associations, two of which constitute novel loci and were first discovered in Latin American populations (BAZ2B and DDIAS). A trans-ethnic meta-analysis revealed another novel cross-population risk locus in CREBBP. Finally, we assessed the performance of a cross-ancestry polygenic risk score in the SCOURGE admixed American cohort. This study constitutes the largest GWAS for COVID-19 hospitalization in admixed Latin Americans conducted to date. This allowed to reveal novel risk loci and emphasize the need of considering the diversity of populations in genomic research.

Duodeno-duodenostomy (DD) has been proposed as a more physiological alternative to conventional duodeno-jejunostomy (DJ) for pancreas transplantation. Accessibility of percutaneous biopsies in these grafts has not yet been assessed. We conducted a retrospective study including all pancreatic percutaneous graft biopsies requested between November 2009 and July 2021. Whenever possible, biopsies were performed under ultrasound (US) guidance or computed tomography (CT) guidance when the US approach failed. Patients were classified into two groups according to surgical technique (DJ and DD). Accessibility, success for histological diagnosis and complications were compared. Biopsy was performed in 93/136 (68.4%) patients in the DJ group and 116/132 (87.9%) of the DD group (p = 0.0001). The graft was not accessible for biopsy mainly due to intestinal loop interposition (n = 29 DJ, n = 10 DD). Adequate sample for histological diagnosis was obtained in 86/93 (92.5%) of the DJ group and 102/116 (87.9%) of the DD group (p = 0.2777). One minor complication was noted in the DD group. The retrocolic position of the DD pancreatic graft does not limit access to percutaneous biopsy. This is a safe technique with a high histological diagnostic success rate.

Abstract Background 5’ untranslated regions (5’UTRs) are essential modulators of protein translation. Predicting the impact of 5’UTR variants is challenging and rarely performed in routine diagnostics. Here, we present a combined approach of a comprehensive prioritization strategy and functional assays to evaluate 5’UTR variation in two large cohorts of patients with inherited retinal diseases (IRDs). Methods We performed an isoform-level re-analysis of retinal RNA-seq data to identify the protein-coding transcripts of 378 IRD genes with highest expression in retina. We evaluated the coverage of their 5’UTRs by different whole exome sequencing (WES) kits. The selected 5’UTRs were analyzed in whole genome sequencing (WGS) and WES data from IRD sub-cohorts from the 100,000 Genomes Project (n = 2397 WGS) and an in-house database (n = 1682 WES), respectively. Identified variants were annotated for 5’UTR-relevant features and classified into seven categories based on their predicted functional consequence. We developed a variant prioritization strategy by integrating population frequency, specific criteria for each category, and family and phenotypic data. A selection of candidate variants underwent functional validation using diverse approaches. Results Isoform-level re-quantification of retinal gene expression revealed 76 IRD genes with a non-canonical retina-enriched isoform, of which 20 display a fully distinct 5’UTR compared to that of their canonical isoform. Depending on the probe design, 3–20% of IRD genes have 5’UTRs fully captured by WES. After analyzing these regions in both cohorts, we prioritized 11 (likely) pathogenic variants in 10 genes (ARL3, MERTK, NDP, NMNAT1, NPHP4, PAX6, PRPF31, PRPF4, RDH12, RD3), of which 7 were novel. Functional analyses further supported the pathogenicity of three variants. Mis-splicing was demonstrated for the PRPF31:c.-9+1G>T variant. The MERTK:c.-125G>A variant, overlapping a transcriptional start site, was shown to significantly reduce both luciferase mRNA levels and activity. The RDH12:c.-123C>T variant was found in cis with the hypomorphic RDH12:c.701G>A (p.Arg234His) variant in 11 patients. This 5’UTR variant, predicted to introduce an upstream open reading frame, was shown to result in reduced RDH12 protein but unaltered mRNA levels. Conclusions This study demonstrates the importance of 5’UTR variants implicated in IRDs and provides a systematic approach for 5’UTR annotation and validation that is applicable to other inherited diseases.

Abstract Background Haploinsufficiency of the transcription factor PAX6 is the main cause of congenital aniridia, a genetic disorder characterized by iris and foveal hypoplasia. 11p13 microdeletions altering PAX6 or its downstream regulatory region (DRR) are present in about 25% of patients; however, only a few complex rearrangements have been described to date. Here, we performed nanopore-based whole-genome sequencing to assess the presence of cryptic structural variants (SVs) on the only two unsolved “PAX6-negative” cases from a cohort of 110 patients with congenital aniridia after unsuccessfully short-read sequencing approaches. Results Long-read sequencing (LRS) unveiled balanced chromosomal rearrangements affecting the PAX6 locus at 11p13 in these two patients and allowed nucleotide-level breakpoint analysis. First, we identified a cryptic 4.9 Mb de novo inversion disrupting intron 7 of PAX6, further verified by targeted polymerase chain reaction amplification and sequencing and FISH-based cytogenetic analysis. Furthermore, LRS was decisive in correctly mapping a t(6;11) balanced translocation cytogenetically detected in a second proband with congenital aniridia and considered non-causal 15 years ago. LRS resolved that the breakpoint on chromosome 11 was indeed located at 11p13, disrupting the DNase I hypersensitive site 2 enhancer within the DRR of PAX6, 161 Kb from the causal gene. Patient-derived RNA expression analysis demonstrated PAX6 haploinsufficiency, thus supporting that the 11p13 breakpoint led to a positional effect by cleaving crucial enhancers for PAX6 transactivation. LRS analysis was also critical for mapping the exact breakpoint on chromosome 6 to the highly repetitive centromeric region at 6p11.1. Conclusions In both cases, the LRS-based identified SVs have been deemed the hidden pathogenic cause of congenital aniridia. Our study underscores the limitations of traditional short-read sequencing in uncovering pathogenic SVs affecting low-complexity regions of the genome and the value of LRS in providing insight into hidden sources of variation in rare genetic diseases.

Retinitis pigmentosa (RP) is the most common inherited retinal disease (IRD) and is characterized by photoreceptor degeneration and progressive vision loss. We report 4 patients presenting with RP from 3 unrelated families with variants in TBC1D32, which to date has never been associated with an IRD. To validate TBC1D32 as a putative RP causative gene, we combined Xenopus in vivo approaches and human induced pluripotent stem cell–derived (iPSC-derived) retinal models. Our data showed that TBC1D32 was expressed during retinal development and that it played an important role in retinal pigment epithelium (RPE) differentiation. Furthermore, we identified a role for TBC1D32 in ciliogenesis of the RPE. We demonstrated elongated ciliary defects that resulted in disrupted apical tight junctions, loss of functionality (delayed retinoid cycling and altered secretion balance), and the onset of an epithelial-mesenchymal transition–like phenotype. Last, our results suggested photoreceptor differentiation defects, including connecting cilium anomalies, that resulted in impaired trafficking to the outer segment in cones and rods in TBC1D32 iPSC-derived retinal organoids. Overall, our data highlight a critical role for TBC1D32 in the retina and demonstrate that TBC1D32 mutations lead to RP. We thus identify TBC1D32 as an IRD-causative gene.

Introduction: Inherited retinal dystrophies (IRDs) can be caused by variants in more than 280 genes. The ATP-binding cassette transporter type A4 (ABCA4) gene is one of these genes and has been linked to Stargardt disease type 1 (STGD1), fundus flavimaculatus, cone–rod dystrophy (CRD), and pan-retinal CRD. Approximately 25% of the reported ABCA4 variants affect RNA splicing. In most cases, it is necessary to perform a functional assay to determine the effect of these variants.Methods: Whole genome sequencing (WGS) was performed in one Spanish proband with Stargardt disease. The putative pathogenicity of c.6480-35A>G on splicing was investigated both in silico and in vitro. The in silico approach was based on the deep-learning tool SpliceAI. For the in vitro approach we used a midigene splice assay in HEK293T cells, based on a previously established wild-type midigene (BA29) containing ABCA4 exons 46 to 48.Results: Through the analysis of WGS data, we identified two candidate variants in ABCA4 in one proband: a previously described deletion, c.699_768+342del (p.(Gln234Phefs*5)), and a novel branchpoint variant, c.6480-35A>G. Segregation analysis confirmed that the variants were in trans. For the branchpoint variant, SpliceAI predicted an acceptor gain with a high score (0.47) at position c.6480-47. A midigene splice assay in HEK293T cells revealed the inclusion of the last 47 nucleotides of intron 47 creating a premature stop codon and allowed to categorize the variant as moderately severe. Subsequent analysis revealed the presence of this variant as a second allele besides c.1958G>A p.(Arg653His) in an additional Spanish proband in a large cohort of IRD cases.Conclusion: A splice-altering effect of the branchpoint variant, confirmed by the midigene splice assay, along with the identification of this variant in a second unrelated individual affected with STGD, provides sufficient evidence to classify the variant as likely pathogenic. In addition, this research highlights the importance of studying non-coding regions and performing functional assays to provide a conclusive molecular diagnosis.

Introduction The incidence of hepatocellular carcinoma (HCC) in Budd-Chiari syndrome (BCS) is unknown and there is no validated diagnostic work-up to define the liver nodules with arterial phase hyperenhancement (APHE) suggesting malignancy. This prospective study evaluates HCC incidence in a Western cohort of patients with BCS and assesses the performance of MRI with hepatobiliary contrast (HB-MRI) for nodule characterization. Methods Patients with BCS followed in our hospital were prospectively evaluated by MRI with extracellular contrast (EC-MRI). Nodules with APHE categorized as non-conclusively benign by 2 radiologists, were studied by HB-MRI and reviewed by 2 radiologists blinded to the EC-MRI results. A new EC-MRI one year later and clinical, analytical and sonographic follow-up every six months for a median of 10 years were performed. Results A total of 55 non-conclusively benign nodules with APHE were detected at EC-MRI in 41 patients. While 32 of them were suggestive of HCC by EC-MRI, all the 55 nodules showed increased uptake of hepatobiliary contrast. An unequivocal central scar was seen in 12/55 nodules at HB-MRI regardless of it was not detected on the EC-MRI. None of the nodules was hypointense in the hepatobiliary phase. HCC was not detected during a median of 10 years of follow-up. Conclusions Detection of nodules with APHE is frequent in patients with BCS, but HCC is rare in Western patients with BCS. While EC-MRI may detect nodules suggesting malignancy, the identification of contrast uptake in the hepatobiliary phase at HB-MRI may help to categorize them as benign.

Inherited macular dystrophies (iMDs) are a group of genetic disorders, which affect the central region of the retina. To investigate the genetic basis of iMDs, we used single-molecule Molecular Inversion Probes to sequence 105 maculopathy-associated genes in 1352 patients diagnosed with iMDs. Within this cohort, 39.8% of patients were considered genetically explained by 460 different variants in 49 distinct genes of which 73 were novel variants, with some affecting splicing. The top five most frequent causative genes were ABCA4 (37.2%), PRPH2 (6.7%), CDHR1 (6.1%), PROM1 (4.3%) and RP1L1 (3.1%). Interestingly, variants with incomplete penetrance were revealed in almost one-third of patients considered solved (28.1%), and therefore, a proportion of patients may not be explained solely by the variants reported. This includes eight previously reported variants with incomplete penetrance in addition to CDHR1:c.783G>A and CNGB3:c.1208G>A. Notably, segregation analysis was not routinely performed for variant phasing—a limitation, which may also impact the overall diagnostic yield. The relatively high proportion of probands without any putative causal variant (60.2%) highlights the need to explore variants with incomplete penetrance, the potential modifiers of disease and the genetic overlap between iMDs and age-related macular degeneration. Our results provide valuable insights into the genetic landscape of iMDs and warrant future exploration to determine the involvement of other maculopathy genes.

Abstract Bardet–Biedl syndrome (BBS) is an autosomal recessive ciliopathy characterized by extensive inter- and intra-familial variability, in which oligogenic interactions have been also reported. Our main goal is to elucidate the role of mutational load in the clinical variability of BBS. A cohort of 99 patients from 77 different families with biallelic pathogenic variants in a BBS-associated gene was retrospectively recruited. Human Phenotype Ontology terms were used in the annotation of clinical symptoms. The mutational load in 39 BBS-related genes was studied in index cases using different molecular and next-generation sequencing (NGS) approaches. Candidate allele combinations were analysed using the in silico tools ORVAL and DiGePred. After clinical annotation, 76 out of the 99 cases a priori fulfilled established criteria for diagnosis of BBS or BBS-like. BBS1 alleles, found in 42% of families, were the most represented in our cohort. An increased mutational load was excluded in 41% of the index cases (22/54). Oligogenic inheritance was suspected in 52% of the screened families (23/45), being 40 tested by means of NGS data and 5 only by traditional methods. Together, ORVAL and DiGePred platforms predicted an oligogenic effect in 44% of the triallelic families (10/23). Intrafamilial variable severity could be clinically confirmed in six of the families. Our findings show that the presence of more than two alleles in BBS-associated genes correlated in six families with a more severe phenotype and associated with specific findings, highlighting the role of the mutational load in the management of BBS cases.

Abstract Rare variants affecting host defense against pathogens could be involved in COVID-19 severity and may help explain fatal outcomes in young and middle-aged patients. Our aim was to report the presence of rare genetic variants in certain genes, by using whole exome sequencing, in a selected group of COVID-19 patients under 65 years who required intubation or resulting in death (n = 44). To this end, different etiopathogenic mechanisms were explored using gene prioritization-based analysis in which genes involved in immune response, immunodeficiencies or blood coagulation were studied. We detected 44 different variants of interest, in 29 different patients (66%). Some of these variants were previously described as pathogenic and were located in genes mainly involved in immune response. A network analysis, including the 42 genes with candidate variants, showed three main components, consisting of 25 highly interconnected genes related to immune response and two additional networks composed by genes enriched in carbohydrate metabolism and in DNA metabolism and repair processes. In conclusion, we have detected candidate variants that may potentially influence COVID-19 outcome in our cohort of patients. Further studies are needed to confirm the ultimate role of the genetic variants described in the present study on COVID-19 severity.

Abstract Clinical exome (CE) sequencing has become a first-tier diagnostic test for hereditary diseases; however, its diagnostic rate is around 30–50%. In this study, we aimed to increase the diagnostic yield of CE using a custom reanalysis algorithm. Sequencing data were available for three cohorts using two commercial protocols applied as part of the diagnostic process. Using these cohorts, we compared the performance of general and clinically relevant variant calling and the efficacy of an in-house bioinformatic protocol (FJD-pipeline) in detecting causal variants as compared to commercial protocols. On the whole, the FJD-pipeline detected 99.74% of the causal variants identified by the commercial protocol in previously solved cases. In the unsolved cases, FJD-pipeline detects more INDELs and non-exonic variants, and is able to increase the diagnostic yield in 2.5% and 3.2% in the re-analysis of 78 cancer and 62 cardiovascular cases. These results were considered to design a reanalysis, filtering and prioritization algorithm that was tested by reassessing 68 inconclusive cases of monoallelic autosomal recessive retinal dystrophies increasing the diagnosis by 4.4%. In conclusion, a guided NGS reanalysis of unsolved cases increases the diagnostic yield in genetic disorders, making it a useful diagnostic tool in medical genetics.

Abstract Most consensus recommendations for the genetic diagnosis of neurodevelopmental disorders (NDDs) do not include the use of next generation sequencing (NGS) and are still based on chromosomal microarrays, such as comparative genomic hybridization array (aCGH). This study compares the diagnostic yield obtained by aCGH and clinical exome sequencing in NDD globally and its spectrum of disorders. To that end, 1412 patients clinically diagnosed with NDDs and studied with aCGH were classified into phenotype categories: global developmental delay/intellectual disability (GDD/ID); autism spectrum disorder (ASD); and other NDDs. These categories were further subclassified based on the most frequent accompanying signs and symptoms into isolated forms, forms with epilepsy; forms with micro/macrocephaly and syndromic forms. Two hundred and forty-five patients of the 1412 were subjected to clinical exome sequencing. Diagnostic yield of aCGH and clinical exome sequencing, expressed as the number of solved cases, was compared for each phenotype category and subcategory. Clinical exome sequencing was superior than aCGH for all cases except for isolated ASD, with no additional cases solved by NGS. Globally, clinical exome sequencing solved 20% of cases (versus 5.7% by aCGH) and the diagnostic yield was highest for all forms of GDD/ID and lowest for Other NDDs (7.1% versus 1.4% by aCGH) and ASD (6.1% versus 3% by aCGH). In the majority of cases, diagnostic yield was higher in the phenotype subcategories than in the mother category. These results suggest that NGS could be used as a first-tier test in the diagnostic algorithm of all NDDs followed by aCGH when necessary.

Abstract Background Mesenchymal stromal cells (MSCs) constitute one of the cell types most frequently used in cell therapy. Although several studies have shown the efficacy of these cells to modulate inflammation in different animal models, the results obtained in human clinical trials have been more modest. Here, we aimed at improving the therapeutic properties of MSCs by inducing a transient expression of two molecules that could enhance two different properties of these cells. With the purpose of improving MSC migration towards inflamed sites, we induced a transient expression of the C-X-C chemokine receptor type 4 (CXCR4). Additionally, to augment the anti-inflammatory properties of MSCs, a transient expression of the anti-inflammatory cytokine, interleukin 10 (IL10), was also induced. Methods Human adipose tissue-derived MSCs were transfected with messenger RNAs carrying the codon-optimized versions of CXCR4 and/or IL10. mRNA-transfected MSCs were then studied, first to evaluate whether the characteristic phenotype of MSCs was modified. Additionally, in vitro and also in vivo studies in an LPS-induced inflamed pad model were conducted to evaluate the impact associated to the transient expression of CXCR4 and/or IL10 in MSCs. Results Transfection of MSCs with CXCR4 and/or IL10 mRNAs induced a transient expression of these molecules without modifying the characteristic phenotype of MSCs. In vitro studies then revealed that the ectopic expression of CXCR4 significantly enhanced the migration of MSCs towards SDF-1, while an increased immunosuppression was associated with the ectopic expression of IL10. Finally, in vivo experiments showed that the co-expression of CXCR4 and IL10 increased the homing of MSCs into inflamed pads and induced an enhanced anti-inflammatory effect, compared to wild-type MSCs. Conclusions Our results demonstrate that the transient co-expression of CXCR4 and IL10 enhances the therapeutic potential of MSCs in a local inflammation mouse model, suggesting that these mRNA-modified cells may constitute a new step in the development of more efficient cell therapies for the treatment of inflammatory diseases.

Abstract Inherited retinal dystrophies (IRD) are a highly heterogeneous group of rare diseases with a molecular diagnostic rate of >50%. Reclassification of variants of uncertain significance (VUS) poses a challenge for IRD diagnosis. We collected 668 IRD cases analyzed by our geneticists using two different clinical exome-sequencing tests. We identified 114 unsolved cases pending reclassification of 125 VUS and studied their genomic, functional, and laboratory-specific features, comparing them to pathogenic and likely pathogenic variants from the same cohort (N = 390). While the clinical exome used did not show differences in diagnostic rate, the more IRD-experienced geneticist reported more VUS (p = 4.07e-04). Significantly fewer VUS were reported in recessive cases (p = 2.14e-04) compared to other inheritance patterns, and of all the genes analyzed, ABCA4 and IMPG2 had the lowest and highest VUS frequencies, respectively (p = 3.89e-04, p = 6.93e-03). Moreover, few frameshift and stop-gain variants were found to be informed VUS (p = 6.73e-08 and p = 2.93e-06). Last, we applied five pathogenicity predictors and found there is a significant proof of deleteriousness when all score for pathogenicity in missense variants. Altogether, these results provided input for a set of rules that correctly reclassified ~70% of VUS as pathogenic in validation datasets. Disease- and setting-specific features influence VUS reporting. Comparison with pathogenic and likely pathogenic variants can prioritize VUS more likely to be reclassified as causal.

Abstract Inherited retinal diseases (IRDs), defined by dysfunction or progressive loss of photoreceptors, are disorders characterized by elevated heterogeneity, both at the clinical and genetic levels. Our main goal was to address the genetic landscape of IRD in the largest cohort of Spanish patients reported to date. A retrospective hospital-based cross-sectional study was carried out on 6089 IRD affected individuals (from 4403 unrelated families), referred for genetic testing from all the Spanish autonomous communities. Clinical, demographic and familiar data were collected from each patient, including family pedigree, age of appearance of visual symptoms, presence of any systemic findings and geographical origin. Genetic studies were performed to the 3951 families with available DNA using different molecular techniques. Overall, 53.2% (2100/3951) of the studied families were genetically characterized, and 1549 different likely causative variants in 142 genes were identified. The most common phenotype encountered is retinitis pigmentosa (RP) (55.6% of families, 2447/4403). The most recurrently mutated genes were PRPH2, ABCA4 and RS1 in autosomal dominant (AD), autosomal recessive (AR) and X-linked (XL) NON-RP cases, respectively; RHO, USH2A and RPGR in AD, AR and XL for non-syndromic RP; and USH2A and MYO7A in syndromic IRD. Pathogenic variants c.3386G > T (p.Arg1129Leu) in ABCA4 and c.2276G > T (p.Cys759Phe) in USH2A were the most frequent variants identified. Our study provides the general landscape for IRD in Spain, reporting the largest cohort ever presented. Our results have important implications for genetic diagnosis, counselling and new therapeutic strategies to both the Spanish population and other related populations.

Replication of SARS-CoV-2 in the human population is defined by distributions of mutants that are present at different frequencies within the infected host and can be detected by ultra-deep sequencing techniques. In this study, we examined the SARS-CoV-2 mutant spectra of amplicons from the spike-coding (S-coding) region of 5 nasopharyngeal isolates derived from patients with vaccine breakthrough. Interestingly, all patients became infected with the Alpha variant, but amino acid substitutions that correspond to the Delta Plus, Iota, and Omicron variants were present in the mutant spectra of the resident virus. Deep sequencing analysis of SARS-CoV-2 from patients with vaccine breakthrough revealed a rich reservoir of mutant types and may also identify tolerated substitutions that can be represented in epidemiologically dominant variants.

ABSTRACT Mutant spectra of RNA viruses are important to understand viral pathogenesis and response to selective pressures. There is a need to characterize the complexity of mutant spectra in coronaviruses sampled from infected patients. In particular, the possible relationship between SARS-CoV-2 mutant spectrum complexity and disease associations has not been established. In the present study, we report an ultradeep sequencing (UDS) analysis of the mutant spectrum of amplicons from the nsp12 (polymerase)- and spike (S)-coding regions of 30 nasopharyngeal isolates (diagnostic samples) of SARS-CoV-2 of the first COVID-19 pandemic wave (Madrid, Spain, April 2020) classified according to the severity of ensuing COVID-19. Low-frequency mutations and deletions, counted relative to the consensus sequence of the corresponding isolate, were overwhelmingly abundant. We show that the average number of different point mutations, mutations per haplotype, and several diversity indices was significantly higher in SARS-CoV-2 isolated from patients who developed mild disease than in those associated with moderate or severe disease (exitus). No such bias was observed with RNA deletions. Location of amino acid substitutions in the three-dimensional structures of nsp12 (polymerase) and S suggest significant structural or functional effects. Thus, patients who develop mild symptoms may be a richer source of genetic variants of SARS-CoV-2 than patients with moderate or severe COVID-19. IMPORTANCE The study shows that mutant spectra of SARS-CoV-2 from diagnostic samples differ in point mutation abundance and complexity and that significantly larger values were observed in virus from patients who developed mild COVID-19 symptoms. Mutant spectrum complexity is not a uniform trait among isolates. The nature and location of low-frequency amino acid substitutions present in mutant spectra anticipate great potential for phenotypic diversification of SARS-CoV-2.

Abstract Background Mutations in CRYAA, which encodes the α-crystallin protein, are associated with a spectrum of congenital cataract–microcornea syndromes. Results In this study, we performed clinical examination and subsequent genetic analysis in two unrelated sporadic cases of different geographical origins presenting with a complex phenotype of ocular malformation. Both cases manifested bilateral microphthalmia and severe anterior segment dysgenesis, primarily characterized by congenital aphakia, microcornea, and iris hypoplasia/aniridia. NGS-based analysis revealed two novel single nucleotide variants occurring de novo and affecting the translation termination codon of the CRYAA gene, c.520T > C and c.521A > C. Both variants are predicted to elongate the C-terminal protein domain by one-third of the original length. Conclusions Our report not only expands the mutational spectrum of CRYAA but also identifies the genetic cause of the unusual ocular phenotype described in this report.

Type 2 diabetes mellitus (T2DM) has been associated with an increased risk of cancer, including colon cancer (CC). However, we recently reported no influence of T2DM on CC prognosis, suggesting that any effect might be at the early stages of tumor development. We hypothesized that T2DM may create an environment in the healthy tissue, which acts as a carcinogenesis driver in agreement with the field of cancerization concept. Here, we focused on early carcinogenesis by analyzing paired tumor and normal colonic mucosa samples from the same patients. The proteome of CC and paired mucosa was quantitatively analyzed in 28 individuals (12 diabetics and 16 nondiabetics) by mass spectrometry with isobaric labeling. Out of 3076 identified proteins, 425 were differentially expressed at the tumor in diabetics compared with nondiabetics. In the adjacent mucosa, 143 proteins were differentially expressed in diabetics and nondiabetics. An enrichment analysis of this signature pointed to mitochondria, ribosome, and translation. Only six proteins were upregulated by diabetes both in tumor and mucosa, of which five were mitochondrial proteins. Differential expression in diabetic versus nondiabetic mucosa was confirmed for MRPL53, MRPL18, and TIMM8B. Higher levels of MRPL18, TIMM8B, and EIF1A were also found in normal colon epithelial cells exposed to high‐glucose conditions. We conclude that T2DM is associated with specific molecular changes in the normal mucosa of CC patients, consistent with field of cancerization in a diabetic environment. The mitochondrial protein signature identifies a potential therapeutic target that could underlie the higher risk of CC in diabetics.

The potential involvement of type 2 diabetes mellitus (T2DM) as a risk factor for colon cancer (CC) has been previously reported. While several clinical studies show a higher incidence of CC and a lower survival rate in diabetics, others report no association. Our own experience indicates that diabetes does not seem to worsen the prognosis once the tumor is present. Despite this controversy, there are no wide‐spectrum molecular studies that delve into the impact of T2DM‐related mechanisms in colon carcinogenesis. Here, we present a transcriptomic and proteomic profiling of paired tumor and normal colon mucosa samples in a cohort of 42 CC patients, 23 of which have T2DM. We used gene set enrichment and network approaches to extract relevant pathways in diabetics, referenced them to current knowledge, and tested them using in vitro techniques. Through our transcriptomics approach, we identified an unexpected overlap of pathways overrepresented in diabetics compared to nondiabetics, in both tumor and normal mucosa, including diabetes‐related metabolic and signaling processes. Proteomic approaches highlighted several cancer‐related signaling routes in diabetics found only in normal mucosa, not in tumors. An integration of the transcriptome and proteome analyses suggested the deregulation of key pathways related to colon carcinogenesis which converged on tumor initiation axis TEAD/YAP‐TAZ as a potential initiator of the process. In vitro studies confirmed upregulation of this pathway in nontumor colon cells under high‐glucose conditions. In conclusion, T2DM associates with deregulation of cancer‐related processes in normal colon mucosa adjacent to tissue which has undergone a malignant transformation. These data support that in diabetic patients, the local microenvironment in normal colon mucosa may be a factor driving field cancerization promoting carcinogenesis. Our results set a new framework to study links between diabetes and colon cancer, including a new role of the TEAD/YAP‐TAZ complex as a potential driver.

Populations of RNA viruses are composed of complex and dynamic mixtures of variant genomes that are termed mutant spectra or mutant clouds. This applies also to SARS-CoV-2, and mutations that are detected at low frequency in an infected individual can be dominant (represented in the consensus sequence) in subsequent variants of interest or variants of concern. Here we briefly review the main conclusions of our work on mutant spectrum characterization of hepatitis C virus (HCV) and SARS-CoV-2 at the nucleotide and amino acid levels and address the following two new questions derived from previous results: (i) how is the SARS-CoV-2 mutant and deletion spectrum composition in diagnostic samples, when examined at progressively lower cut-off mutant frequency values in ultra-deep sequencing; (ii) how the frequency distribution of minority amino acid substitutions in SARS-CoV-2 compares with that of HCV sampled also from infected patients. The main conclusions are the following: (i) the number of different mutations found at low frequency in SARS-CoV-2 mutant spectra increases dramatically (50- to 100-fold) as the cut-off frequency for mutation detection is lowered from 0.5% to 0.1%, and (ii) that, contrary to HCV, SARS-CoV-2 mutant spectra exhibit a deficit of intermediate frequency amino acid substitutions. The possible origin and implications of mutant spectrum differences among RNA viruses are discussed.

Resumen: El desarrollo de medicamentos para ciertas enfermedades raras puede verse truncado por la falta de financiación. En estos casos, en ocasiones, los propios pacientes —o sus familiares— financian el ensayo clínico en el que recibirán el medicamento experimental. Hay 3 modelos de autofinanciación: 2 de ellos, «pagar por probar» y «pagar por participar», ya se han puesto en práctica; el tercero, el modelo «plutocrático», que ha sido recientemente planteado, es, hasta la fecha, un modelo teórico. En este trabajo se repasan los beneficios y riesgos científicos, sociales y éticos de los 2 modelos de investigación clínica: «pagar por participar» y el modelo «plutocrático». Una manera frecuente de poder autofinanciar estos ensayos clínicos es la obtención de fondos por micromecenazgo. Los aspectos más controvertidos de esta modalidad de financiación también son abordados en este trabajo desde diversas perspectivas. Por último, se plantea un escenario futuro que permitiera en nuestro país la puesta en marcha de estos ensayos clínicos autofinanciados mediante el modelo plutocrático. Abstract: The development of medicines for certain rare diseases can be frustrated by lack of funding. In certain cases the patients themselves, or their relatives, occasionally fund the clinical trial in which they will be treated with the investigational medicine. There are 3 models of self-funded research: 2 of them, “pay to try” and “pay to participate”, have already been put into practice. The third, the “plutocratic” proposal, which has been recently put forward is still a theoretical model. In this work the scientific, social and ethical benefits and risks of the 2 clinical research models, “pay to participate” and the “plutocratic” proposal, are reviewed. Patient-funded clinical trials are frequently performed through crowdfunding. The most controversial aspects of this funding modality are also addressed in this article from several perspectives. Finally, a future scenario that would allow the launching of self-funded clinical trials in Spain by the “plutocratic” proposal is proposed.

The development of medicines for certain rare diseases can be cut short by lack of funding. In certain cases the patients themselves, or their relatives, occasionally fund the clinical trial in which they will be treated with the investigational medicine. There are three models of self-funded clinical research: two of them, ‘pay to try’ and ‘pay to participate’, have already been put into practice. The third, the ‘plutocratic’ proposal, which has been recently put forward is still a theoretical model. In this work the scientific, social and ethical benefits and risks of the two clinical research models, ‘pay to participate and the ‘plutocratic’ proposal, are reviewed. Patient-funded clinical trials are frequently performed through crowdfunding. The most controversial aspects of this funding modality are also addressed in this article from several perspectives. Finally, a future scenario that would allow the launching of self-funded clinical trials in Spain by the ‘plutocratic’ proposal is proposed. Resumen: El desarrollo de medicamentos para ciertas enfermedades raras puede verse truncado por la falta de financiación. En estos casos, en ocasiones, los propios pacientes o sus familiares financian el ensayo clínico en el que recibirán el medicamento experimental. Hay tres modelos de autofinanciación de ensayos clínicos: dos de ellos, ‘pagar por probar’ y ‘pagar por participar’, ya se han puesto en práctica; el tercero, el modelo ‘plutocrático’, que ha sido recientemente planteado es, hasta la fecha, un modelo teórico. En este trabajo se repasan los beneficios y riesgos científicos, sociales y éticos de los dos modelos de investigación clínica: ‘pagar por participar’ y el modelo ‘plutocrático’. Una manera frecuente de poder autofinanciar estos ensayos clínicos es la obtención de fondos por micromecenazgo. Los aspectos más controvertidos de esta modalidad de financiación también son abordados en este trabajo desde diversas perspectivas. Por último, se plantea un escenario futuro que permitiera en nuestro país la puesta en marcha de estos ensayos clínicos autofinanciados mediante el modelo plutocrático.

The human induced pluripotent stem cell (hiPSC) line, derived from dermal fibroblasts from Leber congenital amaurosis patient with homozygous mutation c.265 T > C, p.Cys89Arg in aryl hydrocarbon receptor-interacting protein-like 1 (AIPL1) was generated by Sendai virus reprogramming. The generated hiPSC line was free of Sendai virus genes, had stable karyotype, carried the homozygous mutation, was immunopositive to pluripotency markers and able to generate all three germ layers upon embryoid body formation.Resource table.Unlabelled TableUnique stem cell line identifierCIPFi001-AAlternative name(s) of stem cell lineLCA –FiPS4F1InstitutionResearch Center Principe Felipe, Eduardo Primo Yufera 3, Valencia, SpainContact information of distributorDunja Lukovic dlukovic@cipf.esType of cell lineiPSCOriginHumanAdditional origin infoAge:31Sex: FemaleEthnicity if known: CaucasianCell SourceDermal fibroblastsClonalityClonalMethod of reprogrammingSendai virusGenetic modificationN/AType of modificationN/AAssociated diseaseLeber congenital amaurosisGene/locusAIPL1 (NM_014336.3), Chr17: g.6337250A > G (hg19); Ex.2 c.[265 T > C];[265 T > C], p.[Cys89Arg];[Cys89Arg]Method of modificationN/AName of transgene or resistanceN/AInducible/constitutive systemN/ADate archived/stock dateApril 2017Cell line repository/bankhttp://www.isciii.es/ISCIII/es/contenidos/fd-el-instituto/fd-organizacion/fd-estructura-directiva/fd-subdireccion-general-investigacion-terapia-celular-medicina-regenerativa/fd-centros-unidades/fd-banco-nacional-lineas-celulares/fd-lineas-celulares-disponibles/lineas-de-celulas-iPS.shtmlEthical approvalEthics Review Board-competent authority approval obtained by the Valencian Authority for Stem Cell Research (Approval number: S:177–15)

A human iPSC line, IISHDOi004-A, from fibroblasts obtained from a patient with Usher syndrome, harboring a homozygous mutation in the USH2A gene (c.2276G>T; p.Cys759Phe) has been generated. Reprogramming factors Oct3/4, Sox2, Klf4, and c-Myc were delivered using Sendai virus.

The human iPSC cell line, GLC-FiPS4F1 (ESi047-A), derived from dermal fibroblast from the patient with congenital glaucoma caused by the mutation of the gene CYP1B1, was generated by non-integrative reprogramming technology using OCT3/4, SOX2, CMYC and KLF4 reprogramming factors.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

We have generated a human iPSC line IISHDOi003-A from fibroblasts of a patient with a dominant optic atrophy ‘plus’ phenotype, harbouring a heterozygous mutation, c.1635C>A; p.Ser545Arg, in the OPA1 gene. Reprogramming factors Oct3/4, Sox2, Klf4, and c-Myc were delivered using Sendai virus.

Usher syndrome (USH) is a rare autosomal recessive disease and the most common inherited form of combined visual and hearing impairment. Up to 13 genes are associated with this disorder, with USH2A being the most prevalent, due partially to the recurrence rate of the c.2299delG mutation. Excluding hearing aids or cochlear implants for hearing impairment, there are no medical solutions available to treat USH patients. The repair of specific mutations by gene editing is, therefore, an interesting strategy that can be explored using the CRISPR/Cas9 system. In this study, this method of gene editing is used to target the c.2299delG mutation on fibroblasts from an USH patient carrying the mutation in homozygosis. Successful in vitro mutation repair was demonstrated using locus-specific RNA-Cas9 ribonucleoproteins with subsequent homologous recombination repair induced by an engineered template supply. Effects on predicted off-target sites in the CRISPR-treated cells were discarded after a targeted deep-sequencing screen. The proven effectiveness and specificity of these correction tools, applied to the c.2299delG pathogenic variant of USH2A, indicates that the CRISPR system should be considered to further explore a potential treatment of USH. Keywords: Usher syndrome, USH2A, c.2299delG, CRISPR, gene editing, RNPs

Human iPSC line, IISHDOi006-A, was obtained from fibroblasts of a patient with Dominant Optic Atrophy (DOA) carrying a heterozygous mutation in the gene ACO2: c.1999G>A; p.Glu667Lys. Reprogramming factors Oct3/4, Sox2, Klf4, and c-Myc were delivered using a non-integrative methodology that involves the use of Sendai virus.

The human induced pluripotent stem cell (hiPSC) line RP1-FiPS4F1 generated from the patient with autosomal recessive retinitis pigmentosa (arRP) caused by homozygous Ser331Cysfs*5 mutation in Mer tyrosine kinase receptor (MERTK) was genetically corrected using CRISPR/Cas9 system. Two isogenic hiPSCs lines, with heterozygous and homozygous correction of c.992_993delCA mutation in the MERTK gene were generated. These cell lines demonstrate normal karyotype, maintain a pluripotent state, and can differentiate toward three germ layers in vitro. These genetically corrected hiPSCs represent accurate controls to study the contribution of the specific genetic change to the disease, and potentially therapeutic material for cell-replacement therapy.

Mutations in PAX6 are involved in several developmental eye disorders. These disorders have considerable phenotypic variability, ranging from panocular forms of congenital aniridia and microphthalmia to isolated anomalies of the anterior or posterior segment. Here, we describe 3 families with variable inter-generational ocular expression of aniridia, iris coloboma, or microphthalmia, and an unusual transmission of PAX6 mutations from an unaffected or mildly affected parent; all of which raised suspicion of gonosomal mosaicism. We first identified two previously known nonsense mutations and one novel likely pathogenic missense variant in PAX6 in probands by means of targeted NGS. The subsequent segregation analysis by Sanger sequencing evidenced the presence of highly probable mosaic events in paternal blood samples. Mosaicism was further confirmed by droplet digital PCR analysis in several somatic tissues of mosaic fathers. Quantification of the mutant allele fraction in parental samples showed a marked deviation from 50%, with a range between 12 and 29% depending on cell type. Gonosomal mosaicsm was definitively confirmed in one of the families thanks to the availability of a sperm sample from the mosaic father. Thus, the recurrence risk in this family was estimated to be about one-third. This is the first report confirming parental PAX6 mosaicism as a cause of disease recurrence in aniridia and other related phenotypes. In addition, we demonstrated that post-zygotic mosaicism is a frequent and underestimated pathogenic mechanism in aniridia, explaining intra-familial phenotypic variability in many cases. Our findings may have substantial implications for genetic counseling in congenital aniridia. Thus, we also highlight the importance of comprehensive genetic screening of parents for new sporadic cases with aniridia or related developmental eye disease to more accurately assess recurrence risk. In conclusion, somatic and/or gonosomal mosaicism should be taken into consideration as a genetic factor to explain not only families with unaffected parents despite multiple affected children but also variable expressivity, apparent de novo cases, and even uncharacterized cases of aniridia and related developmental eye disorders, apparently lacking PAX6 mutations.

INTRODUCTION:Mutations in USH2A cause both isolated Retinitis Pigmentosa (RP) and Usher syndrome (that implies RP and hearing impairment). One of the most frequent variants identified in this gene and among these patients is the p.(Cys759Phe) change. However, the pathogenic role of this allele has been questioned since it was found in homozygosity in two healthy siblings of a Spanish family. To assess the causative role of USH2A p.(Cys759Phe) in autosomal recessive RP (ARRP) and Usher syndrome type II (USH2) and to establish possible genotype-phenotype correlations associated with p.(Cys759Phe), we performed a comprehensive genetic and clinical study in patients suffering from any of the two above-mentioned diseases and carrying at least one p.(Cys759Phe) allele. MATERIALS AND METHODS:Diagnosis was set according to previously reported protocols. Genetic analyses were performed by using classical molecular and Next-Generation Sequencing approaches. Probands of 57 unrelated families were molecularly studied and 63 patients belonging to these families were phenotypically evaluated. RESULTS:Molecular analysis characterized 100% of the cases, identifying: 11 homozygous patients for USH2A p.(Cys759Phe), 42 compound heterozygous patients (12 of them with another missense USH2A pathogenic variant and 30 with a truncating USH2A variant), and 4 patients carrying the p.(Cys759Phe) allele and a pathogenic variant in another RP gene (PROM1, CNGB1 or RP1). No additional causative variants were identified in symptomatic homozygous patients. Statistical analysis of clinical differences between zygosity states yielded differences (p≤0.05) in age at diagnosis of RP and hypoacusis, and progression of visual field loss. Homozygosity of p.(Cys759Phe) and compound heterozygosity with another USH2A missense variant is associated with ARRP or ARRP plus late onset hypoacusis (OR = 20.62, CI = 95%, p = 0.041). CONCLUSIONS:The present study supports the role of USH2A p.(Cys759Phe) in ARRP and USH2 pathogenesis, and demonstrates the clinical differences between different zygosity states. Phenotype-genotype correlations may guide the genetic characterization based upon specific clinical signs and may advise on the clinical management and prognosis based upon a specific genotype.