Screening Mutations of the Monogenic Syndromic High Myopia by Whole Exome Sequencing From MAGIC Project.

Purpose: This observational study aimed to identify mutations in monogenic syndromic high myopia (msHM) using data from reported samples (n = 9370) of the Myopia Associated Genetics and Intervention Consortium (MAGIC) project. Methods: The targeted panel containing 298 msHM-related genes was constructed and screening of clinically actionable variants was performed based on whole exome sequencing. Capillary sequencing was used to verify the identified gene mutations in the probands and perform segregation analysis with their relatives. Results: A total of 381 candidate variants in 84 genes and 85 eye diseases were found to contribute to msHM in 3.6% (335/9370) of patients with HM. Among them, the 22 genes with the most variations accounted for 62.7% of the diagnostic cases. In the genotype-phenotype association analysis, 60% (201/335) of suspected msHM cases were recalled and 25 patients (12.4%) received a definitive genetic diagnosis. Pathogenic variants were distributed in 18 msHM-related diseases, mainly involving retinal dystrophy genes (e.g. TRPM1, CACNA1F, and FZD4), connective tissue disease genes (e.g. FBN1 and COL2A1), corneal or lens development genes (HSF4, GJA8, and MIP), and other genes (TEK). The msHM gene mutation types were allocated to four categories: nonsense mutations (36%), missense mutations (36%), frameshift mutations (20%), and splice site mutations (8%). Conclusions: This study highlights the importance of thorough molecular subtyping of msHM to provide appropriate genetic counselling and multispecialty care for children and adolescents with HM.

In-depth mining of single-cell transcriptome reveals the key immune-regulated loops in age-related macular degeneration.

Introduction: Age-related macular degeneration (AMD), an ever-increasing ocular disease, has become one of the leading causes of irreversible blindness. Recent advances in single-cell genomics are improving our understanding of the molecular mechanisms of AMD. However, the pathophysiology of this multifactorial disease is complicated and still an ongoing challenge. To better understand disease pathogenesis and identify effective targets, we conducted an in-depth analysis of the single-cell transcriptome of AMD. Methods: The cell expression specificity of the gene (CESG) was selected as an index to identify the novel cell markers. A computational framework was designed to explore the cell-specific TF regulatory loops, containing the interaction of gene pattern signatures, transcription factors regulons, and differentially expressed genes. Results: Three potential novel cell markers were DNASE1L3 for endothelial cells, ABCB5 for melanocytes, and SLC39A12 for RPE cells detected. We observed a notable change in the cell abundance and crosstalk of fibroblasts cells, melanocytes, schwann cells, and T/NK cells between AMD and controls, representing a complex cellular ecosystem in disease status. Finally, we identified six cell type related and three disease-associated ternary loops and elaborated on the robust association between key immune-pathway and AMD. Discussion: In conclusion, this study facilitates the optimization of screening for AMD-related receptor ligand pathways and proposes to further improve the interpretability of disease associations from single-cell data. It illuminated that immune-related regulation paths could be used as potential diagnostic markers for AMD, and in the future, also as therapeutic targets, providing insights into AMD diagnosis and potential interventions.

Sequencing of 19,219 exomes identifies a low-frequency variant in FKBP5 promoter predisposing to high myopia in a Han Chinese population.

High myopia (HM) is one of the leading causes of visual impairment and blindness worldwide. Here, we report a whole-exome sequencing (WES) study in 9,613 HM cases and 9,606 controls of Han Chinese ancestry to pinpoint HM-associated risk variants. Single-variant association analysis identified three newly identified -genetic loci associated with HM, including an East Asian ancestry-specific low-frequency variant (rs533280354) in FKBP5. Multi-ancestry meta-analysis with WES data of 2,696 HM cases and 7,186 controls of European ancestry from the UK Biobank discerned a newly identified European ancestry-specific rare variant in FOLH1. Functional experiments revealed a mechanism whereby a single G-to-A transition at rs533280354 disrupted the binding of transcription activator KLF15 to the promoter of FKBP5, resulting in decreased transcription of FKBP5. Furthermore, burden tests showed a significant excess of rare protein-truncating variants among HM cases involved in retinal blood vessel morphogenesis and neurotransmitter transport.

Impact of the RNA allosteric effect triggered by single nucleotide polymorphisms on the splicing process.

The susceptibility single nucleotide polymorphisms (SNPs) obtained by genome-wide association studies leave some thorny questions, such as prioritization, false positives and unknown pathogenesis. Previous studies suggested that genetic variation may perturb the RNA secondary structure, influence protein recruitment and binding and ultimately affect splicing processes. Therefore, exploring the perturbation of SNPs to structure-function correlations may provide an effective bridge toward understanding the genetic contribution to diseases. Here, aiming to decipher the regulatory mechanism of myopia susceptibility variants, we systematically evaluated the roles of SNP-induced structural changes during splicing. In addition, 7.53% of myopia-related SNPs exhibited significant global structural changes, 19.53% presented noteworthy local structural disturbance and there were wide-ranging structural perturbations in the splice-related motifs. We established a comprehensive evaluation system for structural disturbance in the splicing-related motifs and gave the priority ranking for the SNPs at RNA structural level. These high-priority SNPs were revealed to widely disturb the molecular interaction properties between splicing-related proteins and pre-mRNAs by HDOCK. Moreover, mini-gene assays confirmed that structural perturbation could influence splicing efficiency through structural remodelling. This study deepens our understanding of the potential molecular regulatory mechanisms of susceptible SNPs in myopia and contributes to personalized diagnosis, personalized medicine, disease-risk prediction and functional verification study by guiding the prioritization of the susceptibility SNPs.

RNAsmc: A integrated tool for comparing RNA secondary structure and evaluating allosteric effects.

RNA structure plays a crucial role in gene regulation, in RNA stability and the essential biological processes. RNA secondary structure (RSS) motifs are the basic building blocks for investigating the biological mechanisms of structure. Here, we present a strategy for structural motif-based dynamic alignment, namely, RNA secondary-structural motif-comparing (RNAsmc), to identify structural motifs and quantitatively evaluate their underlying molecular functions. RNAsmc also has strong robustness to sequence length, folding protocol and RNA structural profile by chemical probing. Notably, it is also applicable to quantify structural variation in special RNA editing events (SNVs or SNPs, fragment insertion or deletion, etc.). The findings indicate that RNAsmc can uncover the heterogeneity of RNA secondary structure and score for similarities among components, which provides an impetus to cluster RNA families and evaluate allosteric effects. We find that RNAsmc exhibits remarkable detection efficiency for experimentally-derived RiboSNitches. Finally, the pipeline was assembled into an R software package to serve as an automated toolkit to explore, align, and cluster RSS. It is freely available for download at https://CRAN.R-project.org/package=RNAsmc.

Simulations to Assess the Performance of Multifactor Risk Scores for Predicting Myopia Prevalence in Children and Adolescents in China.

Background: Myopia is the most common visual impairment among Chinese children and adolescents. The purpose of this study is to explore key interventions for myopia prevalence, especially for early-onset myopia and high myopia. Methods: Univariate and multivariate analyses were conducted to evaluate potential associations between risk factor exposure and myopia. LASSO was performed to prioritize the risk features, and the selected leading factors were used to establish the assembled simulation model. Finally, two forecasting models were constructed to predict the risk of myopia and high myopia. Results: Children and adolescents with persistently incorrect posture had a high risk of myopia (OR 7.205, 95% CI 5.999-8.652), which was 2.8 times higher than that in students who always maintained correct posture. In the cohort with high myopia, sleep time of less than 7 h per day (OR 9.789, 95% CI 6.865-13.958), incorrect sitting posture (OR 8.975, 95% CI 5.339-15.086), and siblings with spherical equivalent <-6.00 D (OR 8.439, 95% CI 5.420-13.142) were the top three risk factors. The AUCs of integrated simulation models for myopia and high myopia were 0.8716 and 0.8191, respectively. Conclusion: The findings illustrate that keeping incorrect posture is the leading risk factor for myopia onset, while the onset age of myopia is the primary factor affecting high myopia progression. The age between 8 and 12 years is the crucial stage for clinical intervention, especially for children with parental myopia.

Design, methodology, and baseline of whole city-million scale children and adolescents myopia survey (CAMS) in Wenzhou, China.

BACKGROUND: Myopia is the most common visual impairment in children and adolescents worldwide. This study described an economical and effective population-based screening pipeline and performed the project of a million scale children and adolescents myopia survey (CAMS), which will shed light on the further study of myopia from the level of epidemiology and precision medicine. METHODS: We developed a novel population-based screening pattern, an intelligent screening process and internet-based information transmission and analysis system to carry out the survey consisting of school children in Wenzhou, China. The examination items include unaided distance visual acuity, presenting distance visual acuity, and non-cycloplegic autorefraction. Myopia and high myopia were defined as spherical equivalent (SE) ≤ - 1.00 diopters (D) and SE ≤ - 6.00 D, respectively. Next, the reports of the vision checking were automatically sent to parents and the related departments. The CAMS project will be done two to four times annually with the support of the government. An online eyesight status information management system (OESIMS) was developed to construct comprehensive and efficient electronic vision health records (EVHRs) for myopia information inquiry, risk pre-warning, and further study. RESULTS: The CAMS completed the first-round of screening within 30 days for 99.41% of Wenzhou students from districts and counties, in June 2019. A total of 1,060,925 participants were eligible for CAMS and 1,054,251 (99.37% participation rate) were selected through data quality control, which comprised 1305 schools, and 580,609, 251,050 and 170,967 elementary, middle, and high school students. The mean age of participants was 12.21 ± 3.32 years (6-20 years), the female-to-male ratio was 0.82. The prevalence of myopia in elementary, middle, and high school students was 38.16%, 77.52%, and 84.00%, respectively, and the high myopia incidence was 0.95%, 6.90%, and 12.98%. CONCLUSIONS: The CAMS standardized myopia screening model involves automating large-scale information collection, data transmission, data analysis and early warning, thereby supporting myopia prevention and control. The entire survey reduced 90% of staff, cost, and time consumption compared with previous surveys. This will provide new insights for decision support for public health intervention.

Contribution of structural accessibility to the cooperative relationship of TF-lncRNA in myopia.

Transcriptional regulation is associated with complicated mechanisms including multiple molecular interactions and collaborative drive. Long noncoding RNAs (lncRNAs) have highly structured characteristics and play vital roles in the regulation of transcription in organisms. However, the specific contributions of conformation feature and underlying molecular mechanisms are still unclear. In the present paper, a hypothesis regarding molecular structure effect is presented, which proposes that lncRNAs fold into a complex spatial architecture and act as a skeleton to recruit transcription factors (TF) targeted binding, and which is involved in cooperative regulation. A candidate set of TF-lncRNA coregulation was constructed, and it was found that structural accessibility affected molecular binding force. In addition, transcription factor binding site (TFBS) regions of myopia-related lncRNA transcripts were disturbed, and it was discovered that base mutations affected the occurrence of significant molecular allosteric changes in important elements and variable splicing regions, mediating the onset and development of myopia. The results originated from structureomics and interactionomics and created conditions for systematic research on the mechanisms of structure-mediated TF-lncRNA coregulation in transcriptional regulation. Finally, these findings will help further the understanding of key regulatory roles of molecular allostery in cell physiological and pathological processes.

Identification of the RNase-binding site of SARS-CoV-2 RNA for anchor primer-PCR detection of viral loading in 306 COVID-19 patients.

The pandemic of coronavirus disease 2019 (COVID-19) urgently calls for more sensitive molecular diagnosis to improve sensitivity of current viral nuclear acid detection. We have developed an anchor primer (AP)-based assay to improve viral RNA stability by bioinformatics identification of RNase-binding site of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA and implementing AP dually targeting the N gene of SARS-CoV-2 RNA and RNase 1, 3, 6. The arbitrarily primed polymerase chain reaction (AP-PCR) improvement of viral RNA integrity was supported by (a) the AP increased resistance of the targeted gene (N gene) of SARS-CoV-2 RNA to RNase treatment; (b) the detection of SARS-CoV-2 RNA by AP-PCR with lower cycle threshold values (-2.7 cycles) compared to two commercially available assays; (c) improvement of the viral RNA stability of the ORF gene upon targeting of the N gene and RNase. Furthermore, the improved sensitivity by AP-PCR was demonstrated by detection of SARS-CoV-2 RNA in 70-80% of sputum, nasal, pharyngeal swabs and feces and 36% (4/11) of urine of the confirmed cases (n = 252), 7% convalescent cases (n = 54) and none of 300 negative cases. Lastly, AP-PCR analysis of 306 confirmed and convalescent cases revealed prolonged presence of viral loading for >20 days after the first positive diagnosis. Thus, the AP dually targeting SARS-CoV-2 RNA and RNase improves molecular detection by preserving SARS-CoV-2 RNA integrity and reveals the prolonged viral loading associated with older age and male gender in COVID-19 patients.

Detection of Allosteric Effects of lncRNA Secondary Structures Altered by SNPs in Human Diseases.

Recent studies have shown that structuralized long non-coding RNAs (lncRNAs) play important roles in genetic and epigenetic processes. The spatial structures of most lncRNAs can be altered by distinct in vivo and in vitro cellular environments, as well as by DNA structural variations, such as single-nucleotide polymorphisms (SNPs) and variants (SNVs). In the present study, we extended candidate SNPs that had linkage disequilibria with those significantly associated with lung diseases in genome-wide association studies in order to investigate potential disease mechanisms originating from SNP structural changes of host lncRNAs. Following accurate alignments, we recognized 115 ternary-relationship pairs among 41 SNPs, 10 lncRNA transcripts, and 1 type of lung disease (adenocarcinoma of the lung). Then, we evaluated the structural heterogeneity induced by SNP alleles by developing a local-RNA-structure alignment algorithm and employing randomized strategies to determine the significance of structural variation. We identified four ternary-relationship pairs that were significantly associated with SNP-induced lncRNA allosteric effects. Moreover, these conformational changes disrupted the interactive regions and binding affinities of lncRNA-HCG23 and TF-E2F6, suggesting that these may represent regulatory mechanisms in lung diseases. Taken together, our findings support that SNP-induced changes in lncRNA conformations regulate many biological processes, providing novel insight into the role of the lncRNA "structurome" in human diseases.

Location deviations of DNA functional elements affected SNP mapping in the published databases and references.

The recent extensive application of next-generation sequencing has led to the rapid accumulation of multiple types of data for functional DNA elements. With the advent of precision medicine, the fine-mapping of risk loci based on these elements has become of paramount importance. In this study, we obtained the human reference genome (GRCh38) and the main DNA sequence elements, including protein-coding genes, miRNAs, lncRNAs and single nucleotide polymorphism flanking sequences, from different repositories. We then realigned these elements to identify their exact locations on the genome. Overall, 5%-20% of all sequence element locations deviated among databases, on the scale of kilobase-pair to megabase-pair. These deviations even affected the selection of genome-wide association study risk-associated genes. Our results implied that the location information for functional DNA elements may deviate among public databases. Researchers should take care when using cross-database sources and should perform pilot sequence alignments before element location-based studies.

Landscape of SNPs-mediated lncRNA structural variations and their implication in human complex diseases.

An increasing number of functional studies shows that long noncoding RNAs (lncRNAs) are involved in many aspects of cellular physiology and fulfills a wide variety of regulatory roles at almost every stage of gene expression. A major feature of lncRNAs is the highly folded modular domains in transcripts. With improved modeling and definition, it is now feasible to explore and gain novel insights into the structural-functional relationship of lncRNAs and their association with complex human diseases. In this study, we utilized an automatic computational pipeline to scan lncRNA architecture at the genome-wide scale and to obtain a landscape of functional domains. An accurate alignment algorithm was performed to identify 40 triple pairs between single-nucleotide polymorphisms (SNPs), lncRNAs and diseases. In order to detect the potential contribution of a lncRNA's modular character, we estimated and evaluated structural rearrangements, which were derived from disease-associated SNPs. In addition, we focused on annotating and comparing the global and local heterogeneity of the wild-type and mutant lncRNAs. Assessing lncRNA architecture has yielded how variations in structured regions impact the molecular mechanisms of lncRNAs and how SNPs disturb binding and recruiting ability. These observations are the first glimpse of the 'lncRNA structurome' and make it possible to robustly explore and assemble intricate space conformation and their stress variation. This result also successfully demonstrates that lncRNA transcripts contain a complex structural landscape and highlights the proposed contribution of lncRNA structure in controlling RNA functions and disease mechanisms.

Evaluation of drug efficacy based on the spatial position comparison of drug-target interaction centers.

The spatial position and interaction of drugs and their targets is the most important characteristics for understanding a drug's pharmacological effect, and it could help both in finding new and more precise treatment targets for diseases and in exploring the targeting effects of the new drugs. In this work, we develop a computational pipeline to confirm the spatial interaction relationship of the drugs and their targets and compare the drugs' efficacies based on the interaction centers. First, we produce a 100-sample set to reconstruct a stable docking model of the confirmed drug-target pairs. Second, we set 5.5 Å as the maximum distance threshold for the drug-amino acid residue atom interaction and construct 3-dimensional interaction surface models. Third, by calculating the spatial position of the 3-dimensional interaction surface center, we develop a comparison strategy for estimating the efficacy of different drug-target pairs. For the 1199 drug-target interactions of the 649 drugs and 355 targets, the drugs that have similar interaction center positions tend to have similar efficacies in disease treatment, especially in the analysis of the 37 targeted relationships between the 15 known anti-cancer drugs and 10 target molecules. Furthermore, the analysis of the unpaired anti-cancer drug and target molecules suggests that there is a potential application for discovering new drug actions using the sampling molecular docking and analyzing method. The comparison of the drug-target interaction center spatial position method better reflect the drug-target interaction situations and could support the discovery of new efficacies among the known anti-cancer drugs.

LncRNA Structural Characteristics in Epigenetic Regulation.

The rapid development of new generation sequencing technology has deepened the understanding of genomes and functional products. RNA-sequencing studies in mammals show that approximately 85% of the DNA sequences have RNA products, for which the length greater than 200 nucleotides (nt) is called long non-coding RNAs (lncRNA). LncRNAs now have been shown to play important epigenetic regulatory roles in key molecular processes, such as gene expression, genetic imprinting, histone modification, chromatin dynamics, and other activities by forming specific structures and interacting with all kinds of molecules. This paper mainly discusses the correlation between the structure and function of lncRNAs with the recent progress in epigenetic regulation, which is important to the understanding of the mechanism of lncRNAs in physiological and pathological processes.