VarFish - Collaborative and comprehensive variant analysis for diagnosis and research

VarFish is a user-friendly web application for the quality control, filtering, prioritization, analysis, and user-based annotation of panel and exome variant data for rare disease genetics. It is capable of processing variant call files with single or multiple samples. The variants are automatically annotated with population frequencies, molecular impact, and presence in databases such as ClinVar. Further, it provides support for pathogenicity scores including CADD, MutationTaster, and phenotypic similarity scores. Users can filter variants based on these annotations and presumed inheritance pattern and sort the results by these scores. Filtered variants are listed with their annotations and many useful link-outs to genome browsers, other gene/variant data portals, and external tools for variant assessment. VarFish allows user to create their own annotations including support for variant assessment following ACMG-AMP guidelines. In close collaboration with medical practitioners, VarFish was designed for variant analysis and prioritization in diagnostic and research settings as described in the software’s extensive manual. The user interface has been optimized for supporting these protocols. Users can install VarFish on their own in-house servers where it provides additional lab notebook features for collaborative analysis and allows re-analysis of cases, e.g., after update of genotype or phenotype databases

Privacy in crowdsourcing:a systematic review

The advent of crowdsourcing has brought with it multiple privacy challenges. For example, essential monitoring activities, while necessary and unavoidable, also potentially compromise contributor privacy. We conducted an extensive literature review of the research related to the privacy aspects of crowdsourcing. Our investigation revealed interesting gender differences and also differences in terms of individual perceptions. We conclude by suggesting a number of future research directions.</p

Complete lung agenesis caused by complex genomic rearrangements with neo-TAD formation at the SHH locus

During human organogenesis, lung development is a timely and tightly regulated developmental process under the control of a large number of signaling molecules. Understanding how genetic variants can disturb normal lung development causing different lung malformations is a major goal for dissecting molecular mechanisms during embryogenesis. Here, through exome sequencing (ES), array CGH, genome sequencing (GS) and Hi-C, we aimed at elucidating the molecular basis of bilateral isolated lung agenesis in three fetuses born to a non-consanguineous family. We detected a complex genomic rearrangement containing duplicated, triplicated and deleted fragments involving the SHH locus in fetuses presenting complete agenesis of both lungs and near-complete agenesis of the trachea, diagnosed by ultrasound screening and confirmed at autopsy following termination. The rearrangement did not include SHH itself, but several regulatory elements for lung development, such as MACS1, a major SHH lung enhancer, and the neighboring genes MNX1 and NOM1. The rearrangement incorporated parts of two topologically associating domains (TADs) including their boundaries. Hi-C of cells from one of the affected fetuses showed the formation of two novel TADs each containing SHH enhancers and the MNX1 and NOM1 genes. Hi-C together with GS indicate that the new 3D conformation is likely causative for this condition by an inappropriate activation of MNX1 included in the neo-TADs by MACS1 enhancer, further highlighting the importance of the 3D chromatin conformation in human disease

Variable pulmonary manifestations in Chitayat syndrome: Six additional affected individuals

Hand hyperphalangism leading to shortened index fingers with ulnar deviation, hallux valgus, mild facial dysmorphism and respiratory compromise requiring assisted ventilation are the key features of Chitayat syndrome. This condition results from the recurrent heterozygous missense variant NM_006494.2:c.266A>G; p.(Tyr89Cys) in ERF on chromosome 19q13.2, encoding the ETS2 repressor factor (ERF) protein. The pathomechanism of Chitayat syndrome is unknown. To date, seven individuals with Chitayat syndrome and the recurrent pathogenic ERF variant have been reported in the literature. Here, we describe six additional individuals, among them only one presenting with a history of assisted ventilation, and the remaining presenting with variable pulmonary phenotypes, including one individual without any obvious pulmonary manifestations. Our findings widen the phenotype spectrum caused by the recurrent pathogenic variant in ERF, underline Chitayat syndrome as a cause of isolated skeletal malformations and therefore contribute to the improvement of diagnostic strategies in individuals with hand hyperphalangism

Integration of Hi-C with short and long-read genome sequencing reveals the structure of germline rearranged genomes

Structural variants are a common cause of disease and contribute to a large extent to inter-individual variability, but their detection and interpretation remain a challenge. Here, we investigate 11 individuals with complex genomic rearrangements including germline chromothripsis by combining short- and long-read genome sequencing (GS) with Hi-C. Large-scale genomic rearrangements are identified in Hi-C interaction maps, allowing for an independent assessment of breakpoint calls derived from the GS methods, resulting in >300 genomic junctions. Based on a comprehensive breakpoint detection and Hi-C, we achieve a reconstruction of whole rearranged chromosomes. Integrating information on the three-dimensional organization of chromatin, we observe that breakpoints occur more frequently than expected in lamina-associated domains (LADs) and that a majority reshuffle topologically associating domains (TADs). By applying phased RNA-seq, we observe an enrichment of genes showing allelic imbalanced expression (AIG) within 100 kb around the breakpoints. Interestingly, the AIGs hit by a breakpoint (19/22) display both up- and downregulation, thereby suggesting different mechanisms at play, such as gene disruption and rearrangements of regulatory information. However, the majority of interpretable genes located 200 kb around a breakpoint do not show significant expression changes. Thus, there is an overall robustness in the genome towards large-scale chromosome rearrangements

Integration of Hi-C with short and long-read genome sequencing reveals the structure of germline rearranged genomes

Here the authors characterize structural variations (SVs) in a cohort of individuals with complex genomic rearrangements, identifying breakpoints by employing short- and long-read genome sequencing and investigate their impact on gene expression and the three-dimensional chromatin architecture. They find breakpoints are enriched in inactive regions and can result in chromatin domain fusions.Structural variants are a common cause of disease and contribute to a large extent to inter-individual variability, but their detection and interpretation remain a challenge. Here, we investigate 11 individuals with complex genomic rearrangements including germline chromothripsis by combining short- and long-read genome sequencing (GS) with Hi-C. Large-scale genomic rearrangements are identified in Hi-C interaction maps, allowing for an independent assessment of breakpoint calls derived from the GS methods, resulting in &gt;300 genomic junctions. Based on a comprehensive breakpoint detection and Hi-C, we achieve a reconstruction of whole rearranged chromosomes. Integrating information on the three-dimensional organization of chromatin, we observe that breakpoints occur more frequently than expected in lamina-associated domains (LADs) and that a majority reshuffle topologically associating domains (TADs). By applying phased RNA-seq, we observe an enrichment of genes showing allelic imbalanced expression (AIG) within 100 kb around the breakpoints. Interestingly, the AIGs hit by a breakpoint (19/22) display both up- and downregulation, thereby suggesting different mechanisms at play, such as gene disruption and rearrangements of regulatory information. However, the majority of interpretable genes located 200 kb around a breakpoint do not show significant expression changes. Thus, there is an overall robustness in the genome towards large-scale chromosome rearrangements

Highly multiplexed immune repertoire sequencing links multiple lymphocyte classes with severity of response to COVID-19

BACKGROUND: Disease progression of subjects with coronavirus disease 2019 (COVID-19) varies dramatically. Understanding the various types of immune response to SARS-CoV-2 is critical for better clinical management of coronavirus outbreaks and to potentially improve future therapies. Disease dynamics can be characterized by deciphering the adaptive immune response. METHODS: In this cross-sectional study we analyzed 117 peripheral blood immune repertoires from healthy controls and subjects with mild to severe COVID-19 disease to elucidate the interplay between B and T cells. We used an immune repertoire Primer Extension Target Enrichment method (immunoPETE) to sequence simultaneously human leukocyte antigen (HLA) restricted T cell receptor beta chain (TRB) and unrestricted T cell receptor delta chain (TRD) and immunoglobulin heavy chain (IgH) immune receptor repertoires. The distribution was analyzed of TRB, TRD and IgH clones between healthy and COVID-19 infected subjects. Using McFadden's Adjusted R2 variables were examined for a predictive model. The aim of this study is to analyze the influence of the adaptive immune repertoire on the severity of the disease (value on the World Health Organization Clinical Progression Scale) in COVID-19. FINDINGS: Combining clinical metadata with clonotypes of three immune receptor heavy chains (TRB, TRD, and IgH), we found significant associations between COVID-19 disease severity groups and immune receptor sequences of B and T cell compartments. Logistic regression showed an increase in shared IgH clonal types and decrease of TRD in subjects with severe COVID-19. The probability of finding shared clones of TRD clonal types was highest in healthy subjects (controls). Some specific TRB clones seems to be present in severe COVID-19 (Figure S7b). The most informative models (McFadden´s Adjusted R2=0.141) linked disease severity with immune repertoire measures across all three cell types, as well as receptor-specific cell counts, highlighting the importance of multiple lymphocyte classes in disease progression. INTERPRETATION: Adaptive immune receptor peripheral blood repertoire measures are associated with COVID-19 disease severity

Increased risk of severe clinical course of COVID-19 in carriers of HLA-C*04:01

BACKGROUND: Since the beginning of the coronavirus disease 2019 (COVID-19) pandemic, there has been increasing urgency to identify pathophysiological characteristics leading to severe clinical course in patients infected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Human leukocyte antigen alleles (HLA) have been suggested as potential genetic host factors that affect individual immune response to SARS-CoV-2. We sought to evaluate this hypothesis by conducting a multicenter study using HLA sequencing. METHODS: We analyzed the association between COVID-19 severity and HLAs in 435 individuals from Germany ((n) = 135), Spain ((n) = 133), Switzerland ((n) = 20) and the United States ((n) = 147), who had been enrolled from March 2020 to August 2020. This study included patients older than 18 years, diagnosed with COVID-19 and representing the full spectrum of the disease. Finally, we tested our results by meta-analysing data from prior genome-wide association studies (GWAS). FINDINGS: We describe a potential association of HLA-C*04:01 with severe clinical course of COVID-19. Carriers of HLA-C*04:01 had twice the risk of intubation when infected with SARS-CoV-2 (risk ratio 1.5 [95% CI 1.1-2.1], odds ratio 3.5 [95% CI 1.9-6.6], adjusted (p)-value = 0.0074). These findings are based on data from four countries and corroborated by independent results from GWAS. Our findings are biologically plausible, as HLA-C*04:01 has fewer predicted bindings sites for relevant SARS-CoV-2 peptides compared to other HLA alleles. INTERPRETATION: HLA-C*04:01 carrier state is associated with severe clinical course in SARS-CoV-2. Our findings suggest that HLA class I alleles have a relevant role in immune defense against SARS-CoV-2