The UCSC Genome Browser database: 2024 update.

The UCSC Genome Browser (https://genome.ucsc.edu) is a web-based genomic visualization and analysis tool that serves data to over 7,000 distinct users per day worldwide. It provides annotation data on thousands of genome assemblies, ranging from human to SARS-CoV2. This year, we have introduced new data from the Human Pangenome Reference Consortium and on viral genomes including SARS-CoV2. We have added 1,200 new genomes to our GenArk genome system, increasing the overall diversity of our genomic representation. We have added support for nine new user-contributed track hubs to our public hub system. Additionally, we have released 29 new tracks on the human genome and 11 new tracks on the mouse genome. Collectively, these new features expand both the breadth and depth of the genomic knowledge that we share publicly with users worldwide.

The UCSC Genome Browser database: 2023 update.

The UCSC Genome Browser (https://genome.ucsc.edu) is an omics data consolidator, graphical viewer, and general bioinformatics resource that continues to serve the community as it enters its 23rd year. This year has seen an emphasis in clinical data, with new tracks and an expanded Recommended Track Sets feature on hg38 as well as the addition of a single cell track group. SARS-CoV-2 continues to remain a focus, with regular annotation updates to the browser and continued curation of our phylogenetic sequence placing tool, hgPhyloPlace, whose tree has now reached over 12M sequences. Our GenArk resource has also grown, offering over 2500 hubs and a system for users to request any absent assemblies. We have expanded our bigBarChart display type and created new ways to visualize data via bigRmsk and dynseq display. Displaying custom annotations is now easier due to our chromAlias system which eliminates the requirement for renaming sequence names to the UCSC standard. Users involved in data generation may also be interested in our new tools and trackDb settings which facilitate the creation and display of their custom annotations.

Methylation of the 4q35 D4Z4 repeat defines disease status in facioscapulohumeral muscular dystrophy.

Genetic diagnosis of facioscapulohumeral muscular dystrophy (FSHD) remains a challenge in clinical practice as it cannot be detected by standard sequencing methods despite being the third most common muscular dystrophy. The conventional diagnostic strategy addresses the known genetic parameters of FSHD: the required presence of a permissive haplotype, a size reduction of the D4Z4 repeat of chromosome 4q35 (defining FSHD1) or a pathogenic variant in an epigenetic suppressor gene (consistent with FSHD2). Incomplete penetrance and epistatic effects of the underlying genetic parameters as well as epigenetic parameters (D4Z4 methylation) pose challenges to diagnostic accuracy and hinder prediction of clinical severity. In order to circumvent the known limitations of conventional diagnostics and to complement genetic parameters with epigenetic ones, we developed and validated a multistage diagnostic workflow that consists of a haplotype analysis and a high-throughput methylation profile analysis (FSHD-MPA). FSHD-MPA determines the average global methylation level of the D4Z4 repeat array as well as the regional methylation of the most distal repeat unit by combining bisulphite conversion with next-generation sequencing and a bioinformatics pipeline and uses these as diagnostic parameters. We applied the diagnostic workflow to a cohort of 148 patients and compared the epigenetic parameters based on FSHD-MPA to genetic parameters of conventional genetic testing. In addition, we studied the correlation of repeat length and methylation level within the most distal repeat unit with age-corrected clinical severity and age at disease onset in FSHD patients. The results of our study show that FSHD-MPA is a powerful tool to accurately determine the epigenetic parameters of FSHD, allowing discrimination between FSHD patients and healthy individuals, while simultaneously distinguishing FSHD1 and FSHD2. The strong correlation between methylation level and clinical severity indicates that the methylation level determined by FSHD-MPA accounts for differences in disease severity among individuals with similar genetic parameters. Thus, our findings further confirm that epigenetic parameters rather than genetic parameters represent FSHD disease status and may serve as a valuable biomarker for disease status.

The UCSC Genome Browser database: 2022 update.

The UCSC Genome Browser, https://genome.ucsc.edu, is a graphical viewer for exploring genome annotations. The website provides integrated tools for visualizing, comparing, analyzing, and sharing both publicly available and user-generated genomic datasets. Data highlights this year include a collection of easily accessible public hub assemblies on new organisms, now featuring BLAT alignment and PCR capabilities, and new and updated clinical tracks (gnomAD, DECIPHER, CADD, REVEL). We introduced a new Track Sets feature and enhanced variant displays to aid in the interpretation of clinical data. We also added a tool to rapidly place new SARS-CoV-2 genomes in a global phylogenetic tree enabling researchers to view the context of emerging mutations in our SARS-CoV-2 Genome Browser. Other new software focuses on usability features, including more informative mouseover displays and new fonts.

The UCSC Genome Browser database: 2021 update.

For more than two decades, the UCSC Genome Browser database (https://genome.ucsc.edu) has provided high-quality genomics data visualization and genome annotations to the research community. As the field of genomics grows and more data become available, new modes of display are required to accommodate new technologies. New features released this past year include a Hi-C heatmap display, a phased family trio display for VCF files, and various track visualization improvements. Striving to keep data up-to-date, new updates to gene annotations include GENCODE Genes, NCBI RefSeq Genes, and Ensembl Genes. New data tracks added for human and mouse genomes include the ENCODE registry of candidate cis-regulatory elements, promoters from the Eukaryotic Promoter Database, and NCBI RefSeq Select and Matched Annotation from NCBI and EMBL-EBI (MANE). Within weeks of learning about the outbreak of coronavirus, UCSC released a genome browser, with detailed annotation tracks, for the SARS-CoV-2 RNA reference assembly.

Variant interpretation: UCSC Genome Browser Recommended Track Sets.

The UCSC Genome Browser has been an important tool for genomics and clinical genetics since the sequence of the human genome was first released in 2000. As it has grown in scope to display more types of data it has also grown more complicated. The data, which are dispersed at many locations worldwide, are collected into one view on the Browser, where the graphical interface presents the data in one location. This supports the expertise of the researcher to interpret variants in the genome. Because the analysis of single nucleotide variants and copy number variants require interpretation of data at very different genomic scales, different data resources are required. We present here several Recommended Track Sets designed to facilitate the interpretation of variants in the clinic, offering quick access to datasets relevant to the appropriate scale.

Liquid Biopsy Hotspot Variant Assays: Analytical Validation for Application in Residual Disease Detection and Treatment Monitoring.

BACKGROUND: Analysis of circulating tumor DNA (ctDNA) in plasma is a powerful approach to guide decisions in personalized cancer treatment. Given the low concentration of ctDNA in plasma, highly sensitive methods are required to reliably identify clinically relevant variants. METHODS: We evaluated the suitability of 5 droplet digital PCR (ddPCR) assays targeting KRAS, BRAF, and EGFR variants for ctDNA analysis in clinical use. RESULTS: We investigated assay performance characteristics for very low amounts of variants, showing that the assays had very low limits of blank (0% to 0.11% variant allele frequency, VAF) and limits of quantification (0.41% to 0.7% VAF). Nevertheless, striking differences in detection and quantification of low mutant VAFs between the 5 tested assays were observed, highlighting the need for assay-specific analytical validation. Besides in-depth evaluation, a guide for clinical interpretation of obtained VAFs in plasma was developed, depending on the limits of blank and limits of quantification values. CONCLUSION: It is possible to provide comprehensive clinical reports on actionable variants, allowing minimal residual disease detection and treatment monitoring in liquid biopsy.

Critical assessment of secondary findings in genes linked to primary arrhythmia syndromes.

As comprehensive sequencing technologies gain widespread use, questions about so-called secondary findings (SF) require urgent consideration. The American College of Medical Genetics and Genomics has recommended to report SF in 59 genes (ACMG SF v2.0) including four actionable genes associated with inherited primary arrhythmia syndromes (IPAS) such as catecholaminergic polymorphic ventricular tachycardia, long QT syndrome, and Brugada syndrome. Databases provide conflicting results for the purpose of identifying pathogenic variants in SF associated with IPAS at a level of sufficient evidence for clinical return. As IPAS account for a significant proportion of sudden cardiac deaths (SCD) in young and apparently healthy individuals, variant interpretation has a great impact on diagnosis and prevention of disease. Of 6381 individuals, 0.4% carry pathogenic variants in one of the four actionable genes related to IPAS: RYR2, KCNQ1, KCNH2, and SCN5A. Comparison of the databases ClinVar, Leiden Open-source Variant Database, and Human Gene Mutation Database showed impactful differences (0.2% to 1.3%) in variant interpretation improvable by expert-curation depending on database and classification system used. These data further highlight the need for international consensus regarding the variant interpretation, and subsequently management of SF in particular with regard to treatable arrhythmic disorders with increased risk of SCD.

Comprehensive analysis of the MLH1 promoter region in 480 patients with colorectal cancer and 1150 controls reveals new variants including one with a heritable constitutional MLH1 epimutation.

BACKGROUND: Germline defects in MLH1, MSH2, MSH6 and PMS2 predisposing for Lynch syndrome (LS) are mainly based on sequence changes, whereas a constitutional epimutation of MLH1(CEM) is exceptionally rare. This abnormal MLH1 promoter methylation is not hereditary when arising de novo, whereas a stably heritable and variant-induced CEM was described for one single allele. We searched for MLH1 promoter variants causing a germline or somatic methylation induction or transcriptional repression. METHODS: We analysed the MLH1 promoter sequence in five different patient groups with colorectal cancer (CRC) (n=480) composed of patients with i) CEM (n=16), ii) unsolved loss of MLH1 expression in CRC (n=37), iii) CpG-island methylator-phenotype CRC (n=102), iv) patients with LS (n=83) and v) MLH1-proficient CRC (n=242) as controls. 1150 patients with non-LS tumours also served as controls to correctly judge the results. RESULTS: We detected 10 rare MLH1 promoter variants. One novel, complex MLH1 variant c.-63_-58delins18 is present in a patient with CRC with CEM and his sister, both showing a complete allele-specific promoter methylation and transcriptional silencing. The other nine promoter variants detected in 17 individuals were not associated with methylation. For four of these, a normal, biallelic MLH1 expression was found in the patients' cDNA. CONCLUSION: We report the second promoter variant stably inducing a hereditary CEM. Concerning the classification of promoter variants, we discuss contradictory results from the literature for two variants, describe classification discrepancies between existing rules for five variants, suggest the (re-)classification of five promoter variants to (likely) benign and regard four variants as functionally unclear.

Identification and Functional Testing of ERCC2 Mutations in a Multi-national Cohort of Patients with Familial Breast- and Ovarian Cancer.

The increasing application of gene panels for familial cancer susceptibility disorders will probably lead to an increased proposal of susceptibility gene candidates. Using ERCC2 DNA repair gene as an example, we show that proof of a possible role in cancer susceptibility requires a detailed dissection and characterization of the underlying mutations for genes with diverse cellular functions (in this case mainly DNA repair and basic cellular transcription). In case of ERCC2, panel sequencing of 1345 index cases from 587 German, 405 Lithuanian and 353 Czech families with breast and ovarian cancer (BC/OC) predisposition revealed 25 mutations (3 frameshift, 2 splice-affecting, 20 missense), all absent or very rare in the ExAC database. While 16 mutations were unique, 9 mutations showed up repeatedly with population-specific appearance. Ten out of eleven mutations that were tested exemplarily in cell-based functional assays exert diminished excision repair efficiency and/or decreased transcriptional activation capability. In order to provide evidence for BC/OC predisposition, we performed familial segregation analyses and screened ethnically matching controls. However, unlike the recently published RECQL example, none of our recurrent ERCC2 mutations showed convincing co-segregation with BC/OC or significant overrepresentation in the BC/OC cohort. Interestingly, we detected that some deleterious founder mutations had an unexpectedly high frequency of > 1% in the corresponding populations, suggesting that either homozygous carriers are not clinically recognized or homozygosity for these mutations is embryonically lethal. In conclusion, we provide a useful resource on the mutational landscape of ERCC2 mutations in hereditary BC/OC patients and, as our key finding, we demonstrate the complexity of correct interpretation for the discovery of "bonafide" breast cancer susceptibility genes.

Mutations in EXOSC2 are associated with a novel syndrome characterised by retinitis pigmentosa, progressive hearing loss, premature ageing, short stature, mild intellectual disability and distinctive gestalt.

BACKGROUND: Retinitis pigmentosa in combination with hearing loss can be a feature of different Mendelian disorders. We describe a novel syndrome caused by biallelic mutations in the 'exosome component 2' (EXOSC2) gene. METHODS: Clinical ascertainment of three similar affected patients followed by whole exome sequencing. RESULTS: Three individuals from two unrelated German families presented with a novel Mendelian disorder encompassing childhood myopia, early onset retinitis pigmentosa, progressive sensorineural hearing loss, hypothyroidism, short stature, brachydactyly, recognisable facial gestalt, premature ageing and mild intellectual disability. Whole exome sequencing revealed homozygous or compound heterozygous missense variants in the EXOSC2 gene in all three patients. EXOSC2 encodes the 'ribosomal RNA-processing protein 4' (RRP4)-one of the core components of the RNA exosome. The RNA exosome is a multiprotein complex that plays key roles in RNA processing and degradation. Intriguingly, the EXOSC2-associated phenotype shows only minimal overlap with the previously reported diseases associated with mutations in the RNA exosome core component genes EXOSC3 and EXOSC8. CONCLUSION: We report a novel condition that is probably caused by altered RNA exosome function and expands the spectrum of clinical consequences of impaired RNA metabolism.

Early-onset leukoencephalopathy due to a homozygous missense mutation in the DARS2 gene.

Mutations in the DARS2 gene are known to cause leukoencephalopathy with brainstem and spinal cord involvement and lactate elevation (LBSL), a rare autosomal recessive neurological disorder. It was originally described as juvenile-onset slowly progressive ataxia and spasticity, but recent reports suggest a broader clinical spectrum. Most patients were found to carry compound heterozygous DARS2 mutations, and only very few patients with homozygous mutations have been described so far. We present here an 8-month-old boy carrying a homozygous missense mutation in DARS2 who clinically showed severe neurological deterioration after a respiratory tract infection, followed by an almost complete remission of symptoms. This report further extends the knowledge about the clinical and molecular genetic spectrum of LBSL.

DMP1 mutations in autosomal recessive hypophosphatemia implicate a bone matrix protein in the regulation of phosphate homeostasis.

Hypophosphatemia is a genetically heterogeneous disease. Here, we mapped an autosomal recessive form (designated ARHP) to chromosome 4q21 and identified homozygous mutations in DMP1 (dentin matrix protein 1), which encodes a non-collagenous bone matrix protein expressed in osteoblasts and osteocytes. Intact plasma levels of the phosphaturic protein FGF23 were clearly elevated in two of four affected individuals, providing a possible explanation for the phosphaturia and inappropriately normal 1,25(OH)2D levels and suggesting that DMP1 may regulate FGF23 expression.

A mutation in VPS35, encoding a subunit of the retromer complex, causes late-onset Parkinson disease.

To identify rare causal variants in late-onset Parkinson disease (PD), we investigated an Austrian family with 16 affected individuals by exome sequencing. We found a missense mutation, c.1858G>A (p.Asp620Asn), in the VPS35 gene in all seven affected family members who are alive. By screening additional PD cases, we saw the same variant cosegregating with the disease in an autosomal-dominant mode with high but incomplete penetrance in two further families with five and ten affected members, respectively. The mean age of onset in the affected individuals was 53 years. Genotyping showed that the shared haplotype extends across 65 kilobases around VPS35. Screening the entire VPS35 coding sequence in an additional 860 cases and 1014 controls revealed six further nonsynonymous missense variants. Three were only present in cases, two were only present in controls, and one was present in cases and controls. The familial mutation p.Asp620Asn and a further variant, c.1570C>T (p.Arg524Trp), detected in a sporadic PD case were predicted to be damaging by sequence-based and molecular-dynamics analyses. VPS35 is a component of the retromer complex and mediates retrograde transport between endosomes and the trans-Golgi network, and it has recently been found to be involved in Alzheimer disease.

Autosomal recessive cortical myoclonic tremor and epilepsy: association with a mutation in the potassium channel associated gene CNTN2.

We characterize a consanguineous Egyptian family with an autosomal recessively inherited familial cortical myoclonic tremor and epilepsy. We used multipoint linkage analysis to map the causative mutation to a 12.7 megabase interval within 1q31.3-q32.2 with a log of odds score of 3.6. For further investigation of the linked region in an efficient and unbiased manner, we performed exome sequencing. Within the suspected region we identified a homozygous single base pair deletion (c.503_503delG) leading to a frameshift in the coding region of the sixth exon of CNTN2 alias TAG-1 (p.Trp168fs), which segregated in the respective family. Many studies point towards an important role of the CNTN2 product contactin 2 in neuronal excitability. Contactin 2, a glycosylphosphatidylinositol-anchored neuronal membrane protein, and another transmembrane protein called contactin associated protein-like 2 (CNTNAP2 alias CASPR2) are together necessary to maintain voltage-gated potassium channels at the juxtaparanodal region. CNTN2 knockout mice were previously reported to suffer from spontaneous seizures and mutations in the CNTNAP2 gene have been described to cause epilepsy in humans. To further delineate the role of CNTN2 in patients with epilepsy, we sequenced the coding exons in 189 Caucasian patients with epilepsy. No recessive mutation was detected and heterozygote carriers of rare CNTN2 variants do not seem to be predisposed to epilepsy. Given the severity of the mutation and the proposed function of the gene, we consider this mutation as the most likely cause for cortical myoclonic tremor and epilepsy in this family.

Optical Genome Mapping as a Potential Routine Clinical Diagnostic Method.

Chromosome analysis (CA) and chromosomal microarray analysis (CMA) have been successfully used to diagnose genetic disorders. However, many conditions remain undiagnosed due to limitations in resolution (CA) and detection of only unbalanced events (CMA). Optical genome mapping (OGM) has the potential to address these limitations by capturing both structural variants (SVs) resulting in copy number changes and balanced rearrangements with high resolution. In this study, we investigated OGM's concordance using 87 SVs previously identified by CA, CMA, or Southern blot. Overall, OGM was 98% concordant with only three discordant cases: (1) uncalled translocation with one breakpoint in a centromere; (2) uncalled duplication with breakpoints in the pseudoautosomal region 1; and (3) uncalled mosaic triplication originating from a marker chromosome. OGM provided diagnosis for three previously unsolved cases: (1) disruption of the SON gene due to a balanced reciprocal translocation; (2) disruption of the NBEA gene due to an inverted insertion; (3) disruption of the TSC2 gene due to a mosaic deletion. We show that OGM is a valid method for the detection of many types of SVs in a single assay and is highly concordant with legacy cytogenomic methods; however, it has limited SV detection capabilities in centromeric and pseudoautosomal regions.

Assessment of the genomic variation in a cattle population by re-sequencing of key animals at low to medium coverage.

BACKGROUND: Genome- and population-wide re-sequencing would allow for most efficient detection of causal trait variants. However, despite a strong decrease of costs for next-generation sequencing in the last few years, re-sequencing of large numbers of individuals is not yet affordable. We therefore resorted to re-sequencing of a limited number of bovine animals selected to explain a major proportion of the population's genomic variation, so called key animals, in order to provide a catalogue of functional variants and a substrate for population- and genome-wide imputation of variable sites. RESULTS: Forty-three animals accounting for about 69 percent of the genetic diversity of the Fleckvieh population, a cattle breed of Southern Germany and Austria, were sequenced with coverages ranging from 4.17 to 24.98 and averaging 7.46. After alignment to the reference genome (UMD3.1) and multi-sample variant calling, more than 17 million variant positions were identified, about 90 percent biallelic single nucleotide variants (SNVs) and 10 percent short insertions and deletions (InDels). The comparison with high-density chip data revealed a sensitivity of at least 92 percent and a specificity of 81 percent for sequencing based genotyping, and 97 percent and 93 percent when a imputation step was included. There are 91,733 variants in coding regions of 18,444 genes, 46 percent being non-synonymous exchanges, of which 575 variants are predicted to cause premature stop codons. Three variants are listed in the OMIA database as causal for specific phenotypes. CONCLUSIONS: Low- to medium-coverage re-sequencing of individuals explaining a major fraction of a population's genomic variation allows for the efficient and reliable detection of most variants. Imputation strongly improves genotype quality of lowly covered samples and thus enables maximum density genotyping by sequencing. The functional annotation of variants provides the basis for exhaustive genotype imputation in the population, e.g., for highest-resolution genome-wide association studies.

Quality assurance within the context of genome diagnostics (a german perspective).

The rapid and dynamic implementation of Next-Generation Sequencing (NGS)-based assays has revolutionized genetic testing, and in the near future, nearly all molecular alterations of the human genome will be diagnosable via massive parallel sequencing. While this progress will further corroborate the central role of human genetics in the multidisciplinary management of patients with genetic disorders, it must be accompanied by quality assurance measures in order to allow the safe and optimal use of knowledge ascertained from genome diagnostics. To achieve this, several valuable tools and guidelines have been developed to support the quality of genome diagnostics. In this paper, authors with experience in diverse aspects of genomic analysis summarize the current status of quality assurance in genome diagnostics, with the aim of facilitating further standardization and quality improvement in one of the core competencies of the field.

Closing the Gap - Detection of 5q-Spinal Muscular Atrophy by Short-Read Next-Generation Sequencing and Unexpected Results in a Diagnostic Patient Cohort.

BACKGROUND: The importance of early diagnosis of 5q-Spinal muscular atrophy (5q-SMA) has heightened as early intervention can significantly improve clinical outcomes. In 96% of cases, 5q-SMA is caused by a homozygous deletion of SMN1. Around 4 % of patients carry a SMN1 deletion and a single-nucleotide variant (SNV) on the other allele. Traditionally, diagnosis is based on multiplex ligation probe amplification (MLPA) to detect homozygous or heterozygous exon 7 deletions in SMN1. Due to high homologies within the SMN1/SMN2 locus, sequence analysis to identify SNVs of the SMN1 gene is unreliable by standard Sanger or short-read next-generation sequencing (srNGS) methods. OBJECTIVE: The objective was to overcome the limitations in high-throughput srNGS with the aim of providing SMA patients with a fast and reliable diagnosis to enable their timely therapy. METHODS: A bioinformatics workflow to detect homozygous SMN1 deletions and SMN1 SNVs on srNGS analysis was applied to diagnostic whole exome and panel testing for suggested neuromuscular disorders (1684 patients) and to fetal samples in prenatal diagnostics (260 patients). SNVs were detected by aligning sequencing reads from SMN1 and SMN2 to an SMN1 reference sequence. Homozygous SMN1 deletions were identified by filtering sequence reads for the ,, gene-determining variant" (GDV). RESULTS: 10 patients were diagnosed with 5q-SMA based on (i) SMN1 deletion and hemizygous SNV (2 patients), (ii) homozygous SMN1 deletion (6 patients), and (iii) compound heterozygous SNVs in SMN1 (2 patients). CONCLUSIONS: Applying our workflow in srNGS-based panel and whole exome sequencing (WES) is crucial in a clinical laboratory, as otherwise patients with an atypical clinical presentation initially not suspected to suffer from SMA remain undiagnosed.

CpG-methylation regulates a class of Epstein-Barr virus promoters.

DNA methylation is the major modification of eukaryotic genomes and plays an essential role in mammalian gene regulation. In general, cytosine-phosphatidyl-guanosine (CpG)-methylated promoters are transcriptionally repressed and nuclear proteins such as MECP2, MBD1, MBD2, and MBD4 bind CpG-methylated DNA and contribute to epigenetic silencing. Methylation of viral DNA also regulates gene expression of Epstein-Barr virus (EBV), which is a model of herpes virus latency. In latently infected human B cells, the viral DNA is CpG-methylated, the majority of viral genes is repressed and virus synthesis is therefore abrogated. EBV's BZLF1 encodes a transcription factor of the AP-1 family (Zta) and is the master gene to overcome viral gene repression. In a genome-wide screen, we now identify and characterize those viral genes, which Zta regulates. Among them are genes essential for EBV's lytic phase, which paradoxically depend on strictly CpG-methylated promoters for their Zta-induced expression. We identified novel DNA recognition motifs, termed meZRE (methyl-Zta-responsive element), which Zta selectively binds in order to 'read' DNA in a methylation- and sequence-dependent manner unlike any other known protein. Zta is a homodimer but its binding characteristics to meZREs suggest a sequential, non-palindromic and bipartite DNA recognition element, which confers superior DNA binding compared to CpG-free ZREs. Our findings indicate that Zta has evolved to transactivate cytosine-methylated, hence repressed, silent promoters as a rule to overcome epigenetic silencing.