Massively parallel assessment of human variants with base editor screens.

Understanding the functional consequences of single-nucleotide variants is critical to uncovering the genetic underpinnings of diseases, but technologies to characterize variants are limiting. Here, we leverage CRISPR-Cas9 cytosine base editors in pooled screens to scalably assay variants at endogenous loci in mammalian cells. We benchmark the performance of base editors in positive and negative selection screens, identifying known loss-of-function mutations in BRCA1 and BRCA2 with high precision. To demonstrate the utility of base editor screens to probe small molecule-protein interactions, we screen against BH3 mimetics and PARP inhibitors, identifying point mutations that confer drug sensitivity or resistance. We also create a library of single guide RNAs (sgRNAs) predicted to generate 52,034 ClinVar variants in 3,584 genes and conduct screens in the presence of cellular stressors, identifying loss-of-function variants in numerous DNA damage repair genes. We anticipate that this screening approach will be broadly useful to readily and scalably functionalize genetic variants.

Machine learning-based reclassification of germline variants of unknown significance: The RENOVO algorithm.

The increasing scope of genetic testing allowed by next-generation sequencing (NGS) dramatically increased the number of genetic variants to be interpreted as pathogenic or benign for adequate patient management. Still, the interpretation process often fails to deliver a clear classification, resulting in either variants of unknown significance (VUSs) or variants with conflicting interpretation of pathogenicity (CIP); these represent a major clinical problem because they do not provide useful information for decision-making, causing a large fraction of genetically determined disease to remain undertreated. We developed a machine learning (random forest)-based tool, RENOVO, that classifies variants as pathogenic or benign on the basis of publicly available information and provides a pathogenicity likelihood score (PLS). Using the same feature classes recommended by guidelines, we trained RENOVO on established pathogenic/benign variants in ClinVar (training set accuracy = 99%) and tested its performance on variants whose interpretation has changed over time (test set accuracy = 95%). We further validated the algorithm on additional datasets including unreported variants validated either through expert consensus (ENIGMA) or laboratory-based functional techniques (on BRCA1/2 and SCN5A). On all datasets, RENOVO outperformed existing automated interpretation tools. On the basis of the above validation metrics, we assigned a defined PLS to all existing ClinVar VUSs, proposing a reclassification for 67% with >90% estimated precision. RENOVO provides a validated tool to reduce the fraction of uninterpreted or misinterpreted variants, tackling an area of unmet need in modern clinical genetics.

Sub-genic intolerance, ClinVar, and the epilepsies: A whole-exome sequencing study of 29,165 individuals.

Both mild and severe epilepsies are influenced by variants in the same genes, yet an explanation for the resulting phenotypic variation is unknown. As part of the ongoing Epi25 Collaboration, we performed a whole-exome sequencing analysis of 13,487 epilepsy-affected individuals and 15,678 control individuals. While prior Epi25 studies focused on gene-based collapsing analyses, we asked how the pattern of variation within genes differs by epilepsy type. Specifically, we compared the genetic architectures of severe developmental and epileptic encephalopathies (DEEs) and two generally less severe epilepsies, genetic generalized epilepsy and non-acquired focal epilepsy (NAFE). Our gene-based rare variant collapsing analysis used geographic ancestry-based clustering that included broader ancestries than previously possible and revealed novel associations. Using the missense intolerance ratio (MTR), we found that variants in DEE-affected individuals are in significantly more intolerant genic sub-regions than those in NAFE-affected individuals. Only previously reported pathogenic variants absent in available genomic datasets showed a significant burden in epilepsy-affected individuals compared with control individuals, and the ultra-rare pathogenic variants associated with DEE were located in more intolerant genic sub-regions than variants associated with non-DEE epilepsies. MTR filtering improved the yield of ultra-rare pathogenic variants in affected individuals compared with control individuals. Finally, analysis of variants in genes without a disease association revealed a significant burden of loss-of-function variants in the genes most intolerant to such variation, indicating additional epilepsy-risk genes yet to be discovered. Taken together, our study suggests that genic and sub-genic intolerance are critical characteristics for interpreting the effects of variation in genes that influence epilepsy.

Specifications of the ACMG/AMP Variant Classification Guidelines for Germline DICER1 Variant Curation.

Germline pathogenic variants in DICER1 predispose individuals to develop a variety of benign and malignant tumors. Accurate variant curation and classification is essential for reliable diagnosis of DICER1-related tumor predisposition and identification of individuals who may benefit from surveillance. Since 2015, most labs have followed the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG/AMP) sequence variant classification guidelines for DICER1 germline variant curation. However, these general guidelines lack gene-specific nuances and leave room for subjectivity. Consequently, a group of DICER1 experts joined ClinGen to form the DICER1 and miRNA-Processing Genes Variant Curation Expert Panel (VCEP), to create DICER1- specific ACMG/AMP guidelines for germline variant curation. The VCEP followed the FDA-approved ClinGen protocol for adapting and piloting these guidelines. A diverse set of 40 DICER1 variants were selected for piloting, including 14 known Pathogenic/Likely Pathogenic (P/LP) variants, 12 known Benign/Likely Benign (B/LB) variants, and 14 variants classified as variants of uncertain significance (VUS) or with conflicting interpretations in ClinVar. Clinically meaningful classifications (i.e., P, LP, LB, or B) were achieved for 82.5% (33/40) of the pilot variants, with 100% concordance among the known P/LP and known B/LB variants. Half of the VUS or conflicting variants were resolved with four variants classified as LB and three as LP. These results demonstrate that the DICER1-specific guidelines for germline variant curation effectively classify known pathogenic and benign variants while reducing the frequency of uncertain classifications. Individuals and labs curating DICER1 variants should consider adopting this classification framework to encourage consistency and improve objectivity.

Heme Interactions as Regulators of the Alternative Pathway Complement Responses and Implications for Heme-Associated Pathologies.

Heme (Fe2+-protoporphyrin IX) is a pigment of life, and as a prosthetic group in several hemoproteins, it contributes to diverse critical cellular processes. While its intracellular levels are tightly regulated by networks of heme-binding proteins (HeBPs), labile heme can be hazardous through oxidative processes. In blood plasma, heme is scavenged by hemopexin (HPX), albumin and several other proteins, while it also interacts directly with complement components C1q, C3 and factor I. These direct interactions block the classical pathway (CP) and distort the alternative pathway (AP). Errors or flaws in heme metabolism, causing uncontrolled intracellular oxidative stress, can lead to several severe hematological disorders. Direct interactions of extracellular heme with alternative pathway complement components (APCCs) may be implicated molecularly in diverse conditions at sites of abnormal cell damage and vascular injury. In such disorders, a deregulated AP could be associated with the heme-mediated disruption of the physiological heparan sulphate-CFH coat of stressed cells and the induction of local hemostatic responses. Within this conceptual frame, a computational evaluation of HBMs (heme-binding motifs) aimed to determine how heme interacts with APCCs and whether these interactions are affected by genetic variation within putative HBMs. Combined computational analysis and database mining identified putative HBMs in all of the 16 APCCs examined, with 10 exhibiting disease-associated genetic (SNPs) and/or epigenetic variation (PTMs). Overall, this article indicates that among the pleiotropic roles of heme reviewed, the interactions of heme with APCCs could induce differential AP-mediated hemostasis-driven pathologies in certain individuals.

A Solve-RD ClinVar-based reanalysis of 1522 index cases from ERN-ITHACA reveals common pitfalls and misinterpretations in exome sequencing.

PURPOSE: Within the Solve-RD project (https://solve-rd.eu/), the European Reference Network for Intellectual disability, TeleHealth, Autism and Congenital Anomalies aimed to investigate whether a reanalysis of exomes from unsolved cases based on ClinVar annotations could establish additional diagnoses. We present the results of the "ClinVar low-hanging fruit" reanalysis, reasons for the failure of previous analyses, and lessons learned. METHODS: Data from the first 3576 exomes (1522 probands and 2054 relatives) collected from European Reference Network for Intellectual disability, TeleHealth, Autism and Congenital Anomalies was reanalyzed by the Solve-RD consortium by evaluating for the presence of single-nucleotide variant, and small insertions and deletions already reported as (likely) pathogenic in ClinVar. Variants were filtered according to frequency, genotype, and mode of inheritance and reinterpreted. RESULTS: We identified causal variants in 59 cases (3.9%), 50 of them also raised by other approaches and 9 leading to new diagnoses, highlighting interpretation challenges: variants in genes not known to be involved in human disease at the time of the first analysis, misleading genotypes, or variants undetected by local pipelines (variants in off-target regions, low quality filters, low allelic balance, or high frequency). CONCLUSION: The "ClinVar low-hanging fruit" analysis represents an effective, fast, and easy approach to recover causal variants from exome sequencing data, herewith contributing to the reduction of the diagnostic deadlock.

Clinical implications of conflicting variant interpretations in the cancer genetics clinic.

PURPOSE: The aim of this study was to describe the clinical impact of commercial laboratories issuing conflicting classifications of genetic variants. METHODS: Results from 2000 patients undergoing a multigene hereditary cancer panel by a single laboratory were analyzed. Clinically significant discrepancies between the laboratory-provided test reports and other major commercial laboratories were identified, including differences between pathogenic/likely pathogenic and variant of uncertain significance (VUS) classifications, via review of ClinVar archives. For patients carrying a VUS, clinical documentation was assessed for evidence of provider awareness of the conflict. RESULTS: Fifty of 975 (5.1%) patients with non-negative results carried a variant with a clinically significant conflict, 19 with a pathogenic/likely pathogenic variant reported in APC or MUTYH, and 31 with a VUS reported in CDKN2A, CHEK2, MLH1, MSH2, MUTYH, RAD51C, or TP53. Only 10 of 28 (36%) patients with a VUS with a clinically significant conflict had a documented discussion by a provider about the conflict. Discrepant counseling strategies were used for different patients with the same variant. Among patients with a CDKN2A variant or a monoallelic MUTYH variant, providers were significantly more likely to make recommendations based on the laboratory-reported classification. CONCLUSION: Our findings highlight the frequency of variant interpretation discrepancies and importance of clinician awareness. Guidance is needed on managing patients with discrepant variants to support accurate risk assessment.

Informing variant assessment using structured evidence from prior classifications (PS1, PM5, and PVS1 sequence variant interpretation criteria).

PURPOSE: This study aimed to explore whether evidence of pathogenicity from prior variant classifications in ClinVar could be used to inform variant interpretation using the American College of Medical Genetics and Genomics/Association for Molecular Pathology clinical guidelines. METHODS: We identified distinct single-nucleotide variants (SNVs) that are either similar in location or in functional consequence to pathogenic variants in ClinVar and analyzed evidence in support of pathogenicity using 3 interpretation criteria. RESULTS: Thousands of variants, including many in clinically actionable disease genes (American College of Medical Genetics and Genomics secondary findings v3.0), have evidence of pathogenicity from existing variant classifications, accounting for 2.5% of nonsynonymous SNVs within ClinVar. Notably, there are many variants with uncertain or conflicting classifications that cause the same amino acid substitution as other pathogenic variants (PS1, N = 323), variants that are predicted to cause different amino acid substitutions in the same codon as pathogenic variants (PM5, N = 7692), and loss-of-function variants that are present in genes in which many loss-of-function variants are classified as pathogenic (PVS1, N = 3635). Most of these variants have similar computational predictions of pathogenicity and splicing effect as their associated pathogenic variants. CONCLUSION: Broadly, for >1.4 million SNVs exome wide, information from previously classified variants could be used to provide evidence of pathogenicity. We have developed a pipeline to identify variants meeting these criteria that may inform interpretation efforts.

Variant Classification for Pompe disease; ACMG/AMP specifications from the ClinGen Lysosomal Diseases Variant Curation Expert Panel.

Accurate determination of the clinical significance of genetic variants is critical to the integration of genomics in medicine. To facilitate this process, the NIH-funded Clinical Genome Resource (ClinGen) has assembled Variant Curation Expert Panels (VCEPs), groups of experts and biocurators which provide gene- and disease- specifications to the American College of Medical Genetics & Genomics and Association for Molecular Pathology's (ACMG/AMP) variation classification guidelines. With the goal of classifying the clinical significance of GAA variants in Pompe disease (Glycogen storage disease, type II), the ClinGen Lysosomal Diseases (LD) VCEP has specified the ACMG/AMP criteria for GAA. Variant classification can play an important role in confirming the diagnosis of Pompe disease as well as in the identification of carriers. Furthermore, since the inclusion of Pompe disease on the Recommended Uniform Screening Panel (RUSP) for newborns in the USA in 2015, the addition of molecular genetic testing has become an important component in the interpretation of newborn screening results, particularly for asymptomatic individuals. To date, the LD VCEP has submitted classifications and supporting data on 243 GAA variants to public databases, specifically ClinVar and the ClinGen Evidence Repository. Here, we describe the ACMG/AMP criteria specification process for GAA, an update of the GAA-specific variant classification guidelines, and comparison of the ClinGen LD VCEP's GAA variant classifications with variant classifications submitted to ClinVar. The LD VCEP has added to the publicly available knowledge on the pathogenicity of variants in GAA by increasing the number of expert-curated GAA variants present in ClinVar, and aids in resolving conflicting classifications and variants of uncertain clinical significance.

Two novel RNA-binding proteins identification through computational prediction and experimental validation.

Since RBPs play important roles in the cell, it's particularly important to find new RBPs. We performed iRIP-seq and CLIP-seq to verify two proteins, CLIP1 and DMD, predicted by RBPPred whether are RBPs or not. The experimental results confirm that these two proteins have RNA-binding activity. We identified significantly enriched binding motifs UGGGGAGG, CUUCCG and CCCGU for CLIP1 (iRIP-seq), DMD (iRIP-seq) and DMD (CLIP-seq), respectively. The computational KEGG and GO analysis show that the CLIP1 and DMD share some biological processes and functions. Besides, we found that the SNPs between DMD and its RNA partners may be associated with Becker muscular dystrophy, Duchenne muscular dystrophy, Dilated cardiomyopathy 3B and Cardiovascular phenotype. Among the thirteen cancers data, CLIP1 and another 300 oncogenes always co-occur, and 123 of these 300 genes interact with CLIP1. These cancers may be associated with the mutations occurred in both CLIP1 and the genes it interacts with.

Analysis of hereditary cancer gene variant classifications from ClinVar indicates a need for regular reassessment of clinical assertions.

The clinical classification of variants may change with new information, however, there is limited guidance on how often significant changes in variant classification occur. We used ClinVar to examine how variant classification changes over time. We developed a custom parser and accessed variant data from ClinVar between January 2015 and July 2021. The ClinVar-assigned "aggregate" classification of variants in 121 hereditary cancer genes was harmonized across releases to align to the American College of Medical Genetics and Genomics and the Association for Molecular Pathology terms. Aggregate classification categories were grouped as: benign/likely benign (B/LB); likely pathogenic/pathogenic (LP/P); variant of uncertain significance (VUS); conflicting interpretations of pathogenicity (Conflicting); or Other. We profiled changes in aggregate variant classification between consecutive semi-annual ClinVar releases. The proportion of variants that changed aggregate classification between semi-annual ClinVar releases ranged from 0.6% to 6.4%. The most frequent changes were "VUS to conflicting," "other to LP/P," and "B/LB to Conflicting." A limited number of variants changed aggregate classification from "LP/P to B/LB," or vice versa. Our analysis indicates need for regular reassessment of clinical variant interpretations. The parser developed for this project will facilitate extraction of relevant interpretation data from ClinVar.

Genome Alert!: A standardized procedure for genomic variant reinterpretation and automated gene-phenotype reassessment in clinical routine.

PURPOSE: Retrospective interpretation of sequenced data in light of the current literature is a major concern of the field. Such reinterpretation is manual and both human resources and variable operating procedures are the main bottlenecks. METHODS: Genome Alert! method automatically reports changes with potential clinical significance in variant classification between releases of the ClinVar database. Using ClinVar submissions across time, this method assigns validity category to gene-disease associations. RESULTS: Between July 2017 and December 2019, the retrospective analysis of ClinVar submissions revealed a monthly median of 1247 changes in variant classification with potential clinical significance and 23 new gene-disease associations. Re-examination of 4929 targeted sequencing files highlighted 45 changes in variant classification, and of these classifications, 89% were expert validated, leading to 4 additional diagnoses. Genome Alert! gene-disease association catalog provided 75 high-confidence associations not available in the OMIM morbid list; of which, 20% became available in OMIM morbid list For more than 356 negative exome sequencing data that were reannotated for variants in these 75 genes, this elective approach led to a new diagnosis. CONCLUSION: Genome Alert! (https://genomealert.univ-grenoble-alpes.fr/) enables systematic and reproducible reinterpretation of acquired sequencing data in a clinical routine with limited human resource effect.

Carrier frequency of autosomal recessive genetic conditions in diverse populations: Lessons learned from the genome aggregation database.

An equitable approach by the American College of Medical Genetics and Genomics (ACMG) has recently recommended carrier screening for genes associated with moderate to severe autosomal recessive conditions with a carrier frequency of ≥1/200 in the Genome Aggregation Database exomes (gnomADv2.0.2). We analyzed carrier frequencies in gnomADv3.1.1 genomes representing diverse populations. ClinVar data on 35 996 pathogenic/likely pathogenic variants in 419 genes were used to estimate the gnomAD frequency of heterozygous carriers. We found that ninety-two genes had a carrier frequency of ≥1/200, of which 63 were shared between v3.1.1 and v2.0.2 and 29 were new in v3.1.1. Addition of new populations (Amish, Finnish and Middle Eastern) increased the number of new genes with a carrier frequency of ≥1/200 to 71. Changes in carrier frequencies were attributed to new gnomAD populations, different sample sizes, new ClinVar data, and technical differences between exomes and genomes. This study highlights the dynamic changes in carrier frequencies due to new datasets from diverse populations and provides updated carrier frequencies based on the combined data from 184 352 genomes and exomes in gnomAD. We recommend a periodic review for inclusion of new population data to update carrier screening panels in the future.

A Comprehensive Evaluation of the Performance of Prediction Algorithms on Clinically Relevant Missense Variants.

The rapid integration of genomic technologies in clinical diagnostics has resulted in the detection of a multitude of missense variants whose clinical significance is often unknown. As a result, a plethora of computational tools have been developed to facilitate variant interpretation. However, choosing an appropriate software from such a broad range of tools can be challenging; therefore, systematic benchmarking with high-quality, independent datasets is critical. Using three independent benchmarking datasets compiled from the ClinVar database, we evaluated the performance of ten widely used prediction algorithms with missense variants from 21 clinically relevant genes, including BRCA1 and BRCA2. A fourth dataset consisting of 1053 missense variants was also used to investigate the impact of type 1 circularity on their performance. The performance of the prediction algorithms varied widely across datasets. Based on Matthews Correlation Coefficient and Area Under the Curve, SNPs&GO and PMut consistently displayed an overall above-average performance across the datasets. Most of the tools demonstrated greater sensitivity and negative predictive values at the expense of lower specificity and positive predictive values. We also demonstrated that type 1 circularity significantly impacts the performance of these tools and, if not accounted for, may confound the selection of the best performing algorithms.

Verifying nomenclature of DNA variants in submitted manuscripts: Guidance for journals.

Documenting variation in our genomes is important for research and clinical care. Accuracy in the description of DNA variants is therefore essential. To address this issue, the Human Variome Project convened a committee to evaluate the feasibility of requiring authors to verify that all variants submitted for publication complied with a widely accepted standard for description. After a pilot study of two journals, the committee agreed that requiring authors to verify that variants complied with Human Genome Variation Society nomenclature is a reasonable step toward standardizing the worldwide inventory of human variation.

Identification of Misclassified ClinVar Variants via Disease Population Prevalence.

There is a significant interest in the standardized classification of human genetic variants. We used whole-genome sequence data from 10,495 unrelated individuals to contrast population frequency of pathogenic variants to the expected population prevalence of the disease. Analyses included the ACMG-recommended 59 gene-condition sets for incidental findings and 463 genes associated with 265 OrphaNet conditions. A total of 25,505 variants were used to identify patterns of inflation (i.e., excess genetic risk and misclassification). Inflation increases as the level of evidence supporting the pathogenic nature of the variant decreases. We observed up to 11.5% of genetic disorders with inflation in pathogenic variant sets and up to 92.3% for the variant set with conflicting interpretations. This improved to 7.7% and 57.7%, respectively, after filtering for disease-specific allele frequency. The patterns of inflation were replicated using public data from more than 138,000 genomes. The burden of rare variants was a main contributing factor of the observed inflation, indicating collective misclassified rare variants. We also analyzed the dynamics of re-classification of variant pathogenicity in ClinVar over time, which indicates progressive improvement in variant classification. The study shows that databases include a significant proportion of wrongly ascertained variants; however, it underscores the critical role of ClinVar to contrast claims and foster validation across submitters.

Comprehensive Analysis of Constraint on the Spatial Distribution of Missense Variants in Human Protein Structures.

The spatial distribution of genetic variation within proteins is shaped by evolutionary constraint and provides insight into the functional importance of protein regions and the potential pathogenicity of protein alterations. Here, we comprehensively evaluate the 3D spatial patterns of human germline and somatic variation in 6,604 experimentally derived protein structures and 33,144 computationally derived homology models covering 77% of all human proteins. Using a systematic approach, we quantify differences in the spatial distributions of neutral germline variants, disease-causing germline variants, and recurrent somatic variants. Neutral missense variants exhibit a general trend toward spatial dispersion, which is driven by constraint on core residues. In contrast, germline disease-causing variants are generally clustered in protein structures and form clusters more frequently than recurrent somatic variants identified from tumor sequencing. In total, we identify 215 proteins with significant spatial constraints on the distribution of disease-causing missense variants in experimentally derived protein structures, only 65 (30%) of which have been previously reported. This analysis identifies many clusters not detectable from sequence information alone; only 12% of proteins with significant clustering in 3D were identified from similar analyses of linear protein sequence. Furthermore, spatial analyses of mutations in homology-based structural models are highly correlated with those from experimentally derived structures, supporting the use of computationally derived models. Our approach highlights significant differences in the spatial constraints on different classes of mutations in protein structure and identifies regions of potential function within individual proteins.

Frequency spectrum of rare and clinically relevant markers in multiethnic Indian populations (ClinIndb): A resource for genomic medicine in India.

There have been concerted efforts toward cataloging rare and deleterious variants in different world populations using high-throughput genotyping and sequencing-based methods. The Indian population is underrepresented or its information with respect to clinically relevant variants is sparse in public data sets. The aim of this study was to estimate the burden of monogenic disease-causing variants in Indian populations. Toward this, we have assessed the frequency profile of monogenic phenotype-associated ClinVar variants. The study utilized a genotype data set (global screening array, Illumina) from 2795 individuals (multiple in-house genomics cohorts) representing diverse ethnic and geographically distinct Indian populations. Of the analyzed variants from Global Screening Array, ~9% were found to be informative and were either not known earlier or underrepresented in public databases in terms of their frequencies. These variants were linked to disorders, namely inborn errors of metabolism, monogenic diabetes, hereditary cancers, and various other hereditary conditions. We have also shown that our study cohort is genetically a better representative of the Indian population than its representation in the 1000 Genome Project (South Asians). We have created a database, ClinIndb, linked to the Leiden Open Variation Database, to help clinicians and researchers in diagnosis, counseling, and development of appropriate genetic screening tools relevant to the Indian populations and Indians living abroad.

InterVar: Clinical Interpretation of Genetic Variants by the 2015 ACMG-AMP Guidelines.

In 2015, the American College of Medical Genetics and Genomics (ACMG) and the Association for Molecular Pathology (AMP) published updated standards and guidelines for the clinical interpretation of sequence variants with respect to human diseases on the basis of 28 criteria. However, variability between individual interpreters can be extensive because of reasons such as the different understandings of these guidelines and the lack of standard algorithms for implementing them, yet computational tools for semi-automated variant interpretation are not available. To address these problems, we propose a suite of methods for implementing these criteria and have developed a tool called InterVar to help human reviewers interpret the clinical significance of variants. InterVar can take a pre-annotated or VCF file as input and generate automated interpretation on 18 criteria. Furthermore, we have developed a companion web server, wInterVar, to enable user-friendly variant interpretation with an automated interpretation step and a manual adjustment step. These tools are especially useful for addressing severe congenital or very early-onset developmental disorders with high penetrance. Using results from a few published sequencing studies, we demonstrate the utility of InterVar in significantly reducing the time to interpret the clinical significance of sequence variants.

MARRVEL: Integration of Human and Model Organism Genetic Resources to Facilitate Functional Annotation of the Human Genome.

One major challenge encountered with interpreting human genetic variants is the limited understanding of the functional impact of genetic alterations on biological processes. Furthermore, there remains an unmet demand for an efficient survey of the wealth of information on human homologs in model organisms across numerous databases. To efficiently assess the large volume of publically available information, it is important to provide a concise summary of the most relevant information in a rapid user-friendly format. To this end, we created MARRVEL (model organism aggregated resources for rare variant exploration). MARRVEL is a publicly available website that integrates information from six human genetic databases and seven model organism databases. For any given variant or gene, MARRVEL displays information from OMIM, ExAC, ClinVar, Geno2MP, DGV, and DECIPHER. Importantly, it curates model organism-specific databases to concurrently display a concise summary regarding the human gene homologs in budding and fission yeast, worm, fly, fish, mouse, and rat on a single webpage. Experiment-based information on tissue expression, protein subcellular localization, biological process, and molecular function for the human gene and homologs in the seven model organisms are arranged into a concise output. Hence, rather than visiting multiple separate databases for variant and gene analysis, users can obtain important information by searching once through MARRVEL. Altogether, MARRVEL dramatically improves efficiency and accessibility to data collection and facilitates analysis of human genes and variants by cross-disciplinary integration of 18 million records available in public databases to facilitate clinical diagnosis and basic research.