Reliably Detecting Clinically Important Variants Requires Both Combined Variant Calls and Optimized Filtering Strategies

L-plastin associated syndrome of immune deficiency and hematologic cytopenia

BACKGROUND

LCP1 encodes L-plastin, an actin-bundling protein primarily expressed in hematopoietic cells. In mouse and fish models, LCP1 deficiency has been shown to result in hematologic and immune defects.

OBJECTIVE

This study aimed to determine the nature of a human inborn error of immunity resulting from a novel genetic variant of LCP1.

METHODS

We performed genetic, protein, and cellular analysis of PBMCs from a kindred with apparent autosomal dominant immune deficiency. We identified a candidate causal mutation in LCP1, which we evaluated by engineering the orthologous mutation in mice and Jurkat cells.

RESULTS

A splice-site variant in LCP1 segregated with lymphopenia, neutropenia, and thrombocytopenia. The splicing defect resulted in at least 2 aberrant transcripts, producing an in-frame deletion of 24 nucleotides, and a frameshift deletion of exon 8. Cellular analysis of the kindred revealed a proportionate reduction of T and B cells and a mild expansion of transitional B cells. Similarly, mice carrying the orthologous genetic variant exhibited the same in-frame aberrant transcript, reduced expression Lcp1 and gene dose-dependent leukopenia, mild thrombocytopenia, and lymphopenia, with a significant reduction of T-cell populations. Functional analysis revealed that LCP1c740-1G>A confers a defect in platelet development and function with aberrant spreading on collagen. Immunologic analysis revealed defective actin organization in T cells, reduced migration of PBMCs from patients, splenocytes from mutant mice, and a mutant Jurkat cell line in response to CXCL12; impaired germinal center B-cell expansion after immunization; and reduced cytokinesis during T cell proliferation.

CONCLUSIONS

We describe a unique human hematopoietic defect affecting neutrophils, lymphocytes, and platelets arising from partial LCP1 deficiency.

Identifying cancer driver genes in individual tumours

Cancer is a heterogeneous disease with a strong genetic component making it suitable for precision medicine approaches aimed at identifying the underlying molecular drivers within a tumour. Large scale population-level cancer sequencing consortia have identified many actionable mutations common across both cancer types and sub-types, resulting in an increasing number of successful precision medicine programs. Nonetheless, such approaches fail to consider the effects of mutations unique to an individual patient and may miss rare driver mutations, necessitating personalised approaches to driver-gene prioritisation. One approach is to quantify the functional importance of individual mutations in a single tumour based on how they affect the expression of genes in a gene interaction network (GIN). These GIN-based approaches can be broadly divided into those that utilise an existing reference GIN and those that construct de novo patient-specific GINs. These single-tumour approaches have several limitations that likely influence their results, such as use of reference cohort data, network choice, and approaches to mathematical approximation, and more research is required to evaluate the in vitro and in vivo applicability of their predictions. This review examines the current state of the art methods that identify driver genes in single tumours with a focus on GIN-based driver prioritisation.

A multimorphic mutation in IRF4 causes human autosomal dominant combined immunodeficiency

Interferon regulatory factor 4 (IRF4) is a transcription factor (TF) and key regulator of immune cell development and function. We report a recurrent heterozygous mutation in IRF4, p.T95R, causing an autosomal dominant combined immunodeficiency (CID) in seven patients from six unrelated families. The patients exhibited profound susceptibility to opportunistic infections, notably Pneumocystis jirovecii, and presented with agammaglobulinemia. Patients' B cells showed impaired maturation, decreased immunoglobulin isotype switching, and defective plasma cell differentiation, whereas their T cells contained reduced TH17 and TFH populations and exhibited decreased cytokine production. A knock-in mouse model of heterozygous T95R showed a severe defect in antibody production both at the steady state and after immunization with different types of antigens, consistent with the CID observed in these patients. The IRF4T95R variant maps to the TF's DNA binding domain, alters its canonical DNA binding specificities, and results in a simultaneous multimorphic combination of loss, gain, and new functions for IRF4. IRF4T95R behaved as a gain-of-function hypermorph by binding to DNA with higher affinity than IRF4WT. Despite this increased affinity for DNA, the transcriptional activity on IRF4 canonical genes was reduced, showcasing a hypomorphic activity of IRF4T95R. Simultaneously, IRF4T95R functions as a neomorph by binding to noncanonical DNA sites to alter the gene expression profile, including the transcription of genes exclusively induced by IRF4T95R but not by IRF4WT. This previously undescribed multimorphic IRF4 pathophysiology disrupts normal lymphocyte biology, causing human disease.

Bioinformatic Challenges Detecting Genetic Variation in Precision Medicine Programs

Precision medicine programs to identify clinically relevant genetic variation have been revolutionized by access to increasingly affordable high-throughput sequencing technologies. A decade of continual drops in per-base sequencing costs means it is now feasible to sequence an individual patient genome and interrogate all classes of genetic variation for < $1,000 USD. However, while advances in these technologies have greatly simplified the ability to obtain patient sequence information, the timely analysis and interpretation of variant information remains a challenge for the rollout of large-scale precision medicine programs. This review will examine the challenges and potential solutions that exist in identifying predictive genetic biomarkers and pharmacogenetic variants in a patient and discuss the larger bioinformatic challenges likely to emerge in the future. It will examine how both software and hardware development are aiming to overcome issues in short read mapping, variant detection and variant interpretation. It will discuss the current state of the art for genetic disease and the remaining challenges to overcome for complex disease. Success across all types of disease will require novel statistical models and software in order to ensure precision medicine programs realize their full potential now and into the future.

Compound heterozygous variants in OTULIN are associated with fulminant atypical late-onset ORAS

Autoinflammatory diseases are a heterogenous group of disorders defined by fever and systemic inflammation suggesting involvement of genes regulating innate immune responses. Patients with homozygous loss-of-function variants in the OTU-deubiquitinase OTULIN suffer from neonatal-onset OTULIN-related autoinflammatory syndrome (ORAS) characterized by fever, panniculitis, diarrhea, and arthritis. Here, we describe an atypical form of ORAS with distinct clinical manifestation of the disease caused by two new compound heterozygous variants (c.258G>A (p.M86I)/c.500G>C (p.W167S)) in the OTULIN gene in a 7-year-old affected by a life-threatening autoinflammatory episode with sterile abscess formation. On the molecular level, we find binding of OTULIN to linear ubiquitin to be compromised by both variants; however, protein stability and catalytic activity is most affected by OTULIN variant p.W167S. These molecular changes together lead to increased levels of linear ubiquitin linkages in patient-derived cells triggering the disease. Our data indicate that the spectrum of ORAS patients is more diverse than previously thought and, thus, supposedly asymptomatic individuals might also be affected. Based on our results, we propose to subdivide the ORAS into classical and atypical entities.

Deletions in VANGL1 are a risk factor for antibody-mediated kidney disease

We identify an intronic deletion in VANGL1 that predisposes to renal injury in high risk populations through a kidney-intrinsic process. Half of all SLE patients develop nephritis, yet the predisposing mechanisms to kidney damage remain poorly understood. There is limited evidence of genetic contribution to specific organ involvement in SLE.1,2 We identify a large deletion in intron 7 of Van Gogh Like 1 (VANGL1), which associates with nephritis in SLE patients. The same deletion occurs at increased frequency in an indigenous population (Tiwi Islanders) with 10-fold higher rates of kidney disease compared with non-indigenous populations. Vangl1 hemizygosity in mice results in spontaneous IgA and IgG deposition within the glomerular mesangium in the absence of autoimmune nephritis. Serum transfer into B cell-deficient Vangl1+/- mice results in mesangial IgG deposition indicating that Ig deposits occur in a kidney-intrinsic fashion in the absence of Vangl1. These results suggest that Vangl1 acts in the kidney to prevent Ig deposits and its deficiency may trigger nephritis in individuals with SLE.

Genetic Variation and the Distribution of Variant Types in the Horse

Genetic variation is a key contributor to health and disease. Understanding the link between an individual's genotype and the corresponding phenotype is a major goal of medical genetics. Whole genome sequencing (WGS) within and across populations enables highly efficient variant discovery and elucidation of the molecular nature of virtually all genetic variation. Here, we report the largest catalog of genetic variation for the horse, a species of importance as a model for human athletic and performance related traits, using WGS of 534 horses. We show the extent of agreement between two commonly used variant callers. In data from ten target breeds that represent major breed clusters in the domestic horse, we demonstrate the distribution of variants, their allele frequencies across breeds, and identify variants that are unique to a single breed. We investigate variants with no homozygotes that may be potential embryonic lethal variants, as well as variants present in all individuals that likely represent regions of the genome with errors, poor annotation or where the reference genome carries a variant. Finally, we show regions of the genome that have higher or lower levels of genetic variation compared to the genome average. This catalog can be used for variant prioritization for important equine diseases and traits, and to provide key information about regions of the genome where the assembly and/or annotation need to be improved.

Efficacy of computational predictions of the functional effect of idiosyncratic pharmacogenetic variants

Background

Pharmacogenetic variation is important to drug responses through diverse and complex mechanisms. Predictions of the functional impact of missense pharmacogenetic variants primarily rely on the degree of sequence conservation between species as a primary discriminator. However, idiosyncratic or off-target drug-variant interactions sometimes involve effects that are peripheral or accessory to the central systems in which a gene functions. Given the importance of sequence conservation to functional prediction tools-these idiosyncratic pharmacogenetic variants may violate the assumptions of predictive software commonly used to infer their effect.

Methods

Here we exhaustively assess the effectiveness of eleven missense mutation functional inference tools on all known pharmacogenetic missense variants contained in the Pharmacogenomics Knowledgebase (PharmGKB) repository. We categorize PharmGKB entries into sub-classes to catalog likely off-target interactions, such that we may compare predictions across different variant annotations.

Results

As previously demonstrated, functional inference tools perform variably across the complete set of PharmGKB variants, with large numbers of variants incorrectly classified as 'benign'. However, we find substantial differences amongst PharmGKB variant sub-classes, particularly in variants known to cause off-target, type B adverse drug reactions, that are largely unrelated to the main pharmacological action of the drug. Specifically, variants associated with off-target effects (hence referred to as off-target variants) were most often incorrectly classified as 'benign'. These results highlight the importance of understanding the underlying mechanism of pharmacogenetic variants and how variants associated with off-target effects will ultimately require new predictive algorithms.

Conclusion

In this work we demonstrate that functional inference tools perform poorly on pharmacogenetic variants, particularly on subsets enriched for variants causing off-target, type B adverse drug reactions. We describe how to identify variants associated with off-target effects within PharmGKB in order to generate a training set of variants that is needed to develop new algorithms specifically for this class of variant. Development of such tools will lead to more accurate functional predictions and pave the way for the increased wide-spread adoption of pharmacogenetics in clinical practice.

A rare variant in EZH2 is associated with prostate cancer risk

Prostate cancer (PrCa) is highly heritable, and although rare variants contribute significantly to PrCa risk, few have been identified to date. Herein, whole-genome sequencing was performed in a large PrCa family featuring multiple affected relatives spanning several generations. A rare, predicted splice site EZH2 variant, rs78589034 (G > A), was identified as segregating with disease in all but two individuals in the family, one of whom was affected with lymphoma and bowel cancer and a female relative. This variant was significantly associated with disease risk in combined familial and sporadic PrCa datasets (n = 1551; odds ratio [OR] = 3.55, P = 1.20 × 10^-5 ). Transcriptome analysis was performed on prostate tumour needle biopsies available for two rare variant carriers and two wild-type cases. Although no allele-dependent differences were detected in EZH2 transcripts, a distinct differential gene expression signature was observed when comparing prostate tissue from the rare variant carriers with the wild-type samples. The gene expression signature comprised known downstream targets of EZH2 and included the top-ranked genes, DUSP1, FOS, JUNB and EGR1, which were subsequently validated by qPCR. These data provide evidence that rs78589034 is associated with increased PrCa risk in Tasmanian men and further, that this variant may be associated with perturbed EZH2 function in prostate tissue. Disrupted EZH2 function is a driver of tumourigenesis in several cancers, including prostate, and is of significant interest as a therapeutic target.

Chromosome-length genome assembly and structural variations of the primal Basenji dog (Canis lupus familiaris) genome

BACKGROUND

Basenjis are considered an ancient dog breed of central African origins that still live and hunt with tribesmen in the African Congo. Nicknamed the barkless dog, Basenjis possess unique phylogeny, geographical origins and traits, making their genome structure of great interest. The increasing number of available canid reference genomes allows us to examine the impact the choice of reference genome makes with regard to reference genome quality and breed relatedness.

RESULTS

Here, we report two high quality de novo Basenji genome assemblies: a female, China (CanFam_Bas), and a male, Wags. We conduct pairwise comparisons and report structural variations between assembled genomes of three dog breeds: Basenji (CanFam_Bas), Boxer (CanFam3.1) and German Shepherd Dog (GSD) (CanFam_GSD). CanFam_Bas is superior to CanFam3.1 in terms of genome contiguity and comparable overall to the high quality CanFam_GSD assembly. By aligning short read data from 58 representative dog breeds to three reference genomes, we demonstrate how the choice of reference genome significantly impacts both read mapping and variant detection.

CONCLUSIONS

The growing number of high-quality canid reference genomes means the choice of reference genome is an increasingly critical decision in subsequent canid variant analyses. The basal position of the Basenji makes it suitable for variant analysis for targeted applications of specific dog breeds. However, we believe more comprehensive analyses across the entire family of canids is more suited to a pangenome approach. Collectively this work highlights the importance the choice of reference genome makes in all variation studies.

Detecting Causal Variants in Mendelian Disorders Using Whole-Genome Sequencing

Increasingly affordable sequencing technologies are revolutionizing the field of genomic medicine. It is now feasible to interrogate all major classes of variation in an individual across the entire genome for less than $1000 USD. While the generation of patient sequence information using these technologies has become routine, the analysis and interpretation of this data remains the greatest obstacle to widespread clinical implementation. This chapter summarizes the steps to identify, annotate, and prioritize variant information required for clinical report generation. We discuss methods to detect each variant class and describe strategies to increase the likelihood of detecting causal variant(s) in Mendelian disease. Lastly, we describe a sample workflow for synthesizing large amount of genetic information into concise clinical reports.

Detecting pathogenic variants in autoimmune diseases using high-throughput sequencing

Sequencing the first human genome in 2003 took 15 years and cost $2.7 billion. Advances in sequencing technologies have since decreased costs to the point where it is now feasible to resequence a whole human genome for $1000 in a single day. These advances have allowed the generation of huge volumes of high‐quality human sequence data used to construct increasingly large catalogs of both population‐level and disease‐causing variation. The existence of such databases, coupled with a high‐quality human reference genome, means we are able to interrogate and annotate all types of genetic variation and identify pathogenic variants for many diseases. Increasingly, sequencing‐based approaches are being used to elucidate the underlying genetic cause of autoimmune diseases, a group of roughly 80 polygenic diseases characterized by abnormal immune responses where healthy tissue is attacked. Although sequence data generation has become routine and affordable, significant challenges remain with no gold‐standard methodology to identify pathogenic variants currently available. This review examines the latest methodologies used to identify pathogenic variants in autoimmune diseases and considers available sequencing options and subsequent bioinformatic methodologies and strategies. The development of reliable and robust sequencing and analytic workflows to detect pathogenic variants is critical to realize the potential of precision medicine programs where patient variant information is used to inform clinical practice.

Absence of mucosal-associated invariant T cells in a person with a homozygous point mutation in MR1

The role unconventional T cells play in protective immunity in humans is unclear. Mucosal-associated invariant T (MAIT) cells are an unconventional T cell subset restricted to the antigen-presenting molecule MR1. Here, we report the discovery of a patient homozygous for a rare Arg31His (R9H in the mature protein) mutation in MR1 who has a history of difficult-to-treat viral and bacterial infections. MR1^R9H was unable to present the potent microbially derived MAIT cell stimulatory ligand. The MR1^R9H crystal structure revealed that the stimulatory ligand cannot bind due to the mutation lying within, and causing structural perturbation to, the ligand-binding domain of MR1. While MR1^R9H could bind and be up-regulated by a MAIT cell inhibitory ligand, the patient lacked circulating MAIT cells. This shows the importance of the stimulatory ligand for MAIT cell selection in humans. The patient had an expanded γδ T cell population, indicating a compensatory interplay between these unconventional T cell subsets.

Transcriptional profiling and immunophenotyping show sustained activation of blood monocytes in subpatent Plasmodium falciparum infection

OBJECTIVES

Malaria, caused by Plasmodium infection, remains a major global health problem. Monocytes are integral to the immune response, yet their transcriptional and functional responses in primary Plasmodium falciparum infection and in clinical malaria are poorly understood.

METHODS

The transcriptional and functional profiles of monocytes were examined in controlled human malaria infection with P. falciparum blood stages and in children and adults with acute malaria. Monocyte gene expression and functional phenotypes were examined by RNA sequencing and flow cytometry at peak infection and compared to pre-infection or at convalescence in acute malaria.

RESULTS

In subpatent primary infection, the monocyte transcriptional profile was dominated by an interferon (IFN) molecular signature. Pathways enriched included type I IFN signalling, innate immune response and cytokine-mediated signalling. Monocytes increased TNF and IL-12 production upon in vitro toll-like receptor stimulation and increased IL-10 production upon in vitro parasite restimulation. Longitudinal phenotypic analyses revealed sustained significant changes in the composition of monocytes following infection, with increased CD14⁺CD16^- and decreased CD14^-CD16⁺ subsets. In acute malaria, monocyte CD64/FcγRI expression was significantly increased in children and adults, while HLA-DR remained stable. Although children and adults showed a similar pattern of differentially expressed genes, the number and magnitude of gene expression change were greater in children.

CONCLUSIONS

Monocyte activation during subpatent malaria is driven by an IFN molecular signature with robust activation of genes enriched in pathogen detection, phagocytosis, antimicrobial activity and antigen presentation. The greater magnitude of transcriptional changes in children with acute malaria suggests monocyte phenotypes may change with age or exposure.

Human FCHO1 deficiency reveals role for clathrin-mediated endocytosis in development and function of T cells

Clathrin-mediated endocytosis (CME) is critical for internalisation of molecules across cell membranes. The FCH domain only 1 (FCHO1) protein is key molecule involved in the early stages of CME formation. The consequences of mutations in FCHO1 in humans were unknown. We identify ten unrelated patients with variable T and B cell lymphopenia, who are homozygous for six distinct mutations in FCHO1. We demonstrate that these mutations either lead to mislocalisation of the protein or prevent its interaction with binding partners. Live-cell imaging of cells expressing mutant variants of FCHO1 provide evidence of impaired formation of clathrin coated pits (CCP). Patient T cells are unresponsive to T cell receptor (TCR) triggering. Internalisation of the TCR receptor is severely perturbed in FCHO1-deficient Jurkat T cells but can be rescued by expression of wild-type FCHO1. Thus, we discovered a previously unrecognised critical role of FCHO1 and CME during T-cell development and function in humans.

Lymphoma Driver Mutations in the Pathogenic Evolution of an Iconic Human Autoantibody

Pathogenic autoantibodies arise in many autoimmune diseases, but it is not understood how the cells making them evade immune checkpoints. Here, single-cell multi-omics analysis demonstrates a shared mechanism with lymphoid malignancy in the formation of public rheumatoid factor autoantibodies responsible for mixed cryoglobulinemic vasculitis. By combining single-cell DNA and RNA sequencing with serum antibody peptide sequencing and antibody synthesis, rare circulating B lymphocytes making pathogenic autoantibodies were found to comprise clonal trees accumulating mutations. Lymphoma driver mutations in genes regulating B cell proliferation and V(D)J mutation (CARD11, TNFAIP3, CCND3, ID3, BTG2, and KLHL6) were present in rogue B cells producing the pathogenic autoantibody. Antibody V(D)J mutations conferred pathogenicity by causing the antigen-bound autoantibodies to undergo phase transition to insoluble aggregates at lower temperatures. These results reveal a pre-neoplastic stage in human lymphomagenesis and a cascade of somatic mutations leading to an iconic pathogenic autoantibody.

AMLVaran: a software approach to implement variant analysis of targeted NGS sequencing data in an oncological care setting

BACKGROUND

Next-Generation Sequencing (NGS) enables large-scale and cost-effective sequencing of genetic samples in order to detect genetic variants. After successful use in research-oriented projects, NGS is now entering clinical practice. Consequently, variant analysis is increasingly important to facilitate a better understanding of disease entities and prognoses. Furthermore, variant calling allows to adapt and optimize specific treatments of individual patients, and thus is an integral part of personalized medicine.However, the analysis of NGS data typically requires a number of complex bioinformatics processing steps. A flexible and reliable software that combines the variant analysis process with a simple, user-friendly interface is therefore highly desirable, but still lacking.

RESULTS

With AMLVaran (AML Variant Analyzer), we present a web-based software, that covers the complete variant analysis workflow of targeted NGS samples. The software provides a generic pipeline that allows free choice of variant calling tools and a flexible language (SSDL) for filtering variant lists. AMLVaran's interactive website presents comprehensive annotation data and includes curated information on relevant hotspot regions and driver mutations. A concise clinical report with rule-based diagnostic recommendations is generated.An AMLVaran configuration with eight variant calling tools and a complex scoring scheme, based on the somatic variant calling pipeline appreci8, was used to analyze three datasets from AML and MDS studies with 402 samples in total. Maximum sensitivity and positive predictive values were 1.0 and 0.96, respectively. The tool's usability was found to be satisfactory by medical professionals.

CONCLUSION

Coverage analysis, reproducible variant filtering and software usability are important for clinical assessment of variants. AMLVaran performs reliable NGS variant analyses and generates reports fulfilling the requirements of a clinical setting. Due to its generic design, the software can easily be adapted for use with different targeted panels for other tumor entities, or even for whole-exome data. AMLVaran has been deployed to a public web server and is distributed with Docker scripts for local use.

Recurrent miscalling of missense variation from short-read genome sequence data

BACKGROUND

Short-read resequencing of genomes produces abundant information of the genetic variation of individuals. Due to their numerous nature, these variants are rarely exhaustively validated. Furthermore, low levels of undetected variant miscalling will have a systematic and disproportionate impact on the interpretation of individual genome sequence information, especially should these also be carried through into in reference databases of genomic variation.

RESULTS

We find that sequence variation from short-read sequence data is subject to recurrent-yet-intermittent miscalling that occurs in a sequence intrinsic manner and is very sensitive to sequence read length. The miscalls arise from difficulties aligning short reads to redundant genomic regions, where the rate of sequencing error approaches the sequence diversity between redundant regions. We find the resultant miscalled variants to be sensitive to small sequence variations between genomes, and thereby are often intrinsic to an individual, pedigree, strain or human ethnic group. In human exome sequences, we identify 2-300 recurrent false positive variants per individual, almost all of which are present in public databases of human genomic variation. From the exomes of non-reference strains of inbred mice, we identify 3-5000 recurrent false positive variants per mouse - the number of which increasing with greater distance between an individual mouse strain and the reference C57BL6 mouse genome. We show that recurrently miscalled variants may be reproduced for a given genome from repeated simulation rounds of read resampling, realignment and recalling. As such, it is possible to identify more than two-thirds of false positive variation from only ten rounds of simulation.

CONCLUSION

Identification and removal of recurrent false positive variants from specific individual variant sets will improve overall data quality. Variant miscalls arising are highly sequence intrinsic and are often specific to an individual, pedigree or ethnicity. Further, read length is a strong determinant of whether given false variants will be called for any given genome - which has profound significance for cohort studies that pool datasets collected and sequenced at different points in time.

Comprehensive elaboration of database resources utilized in next-generation sequencing-based tumor somatic mutation detection

The rapid evolution of next-generation sequencing (NGS)-based tumor genomic profile detection and the emergence of molecularly targeted therapies have enabled precision oncology. In NGS-based analysis, various types of databases have been developed to perform different functions. However, many problems still exist when using these public databases. Therefore, it is important to better understand the characteristics and limitations of each database and have them complement each other to provide useful clinical evidence for NGS testing. In this review, we elaborate on the important role of databases and their concrete applications in NGS-based somatic mutation detection. We introduce the typically used databases for sequence alignment, variant filtration, and variant interpretation, and compare the differences between the databases with similar functions. Subsequently, we determine the limitations of each database and provide the corresponding solutions. Furthermore, we present an overview diagram to clearly illustrate the database used in the entire NGS-based somatic mutation detection pipeline.

Gene editing of the multi-copy H2A.B gene and its importance for fertility

BACKGROUND

Altering the biochemical makeup of chromatin by the incorporation of histone variants during development represents a key mechanism in regulating gene expression. The histone variant H2A.B, H2A.B.3 in mice, appeared late in evolution and is most highly expressed in the testis. In the mouse, it is encoded by three different genes. H2A.B expression is spatially and temporally regulated during spermatogenesis being most highly expressed in the haploid round spermatid stage. Active genes gain H2A.B where it directly interacts with polymerase II and RNA processing factors within splicing speckles. However, the importance of H2A.B for gene expression and fertility are unknown.

RESULTS

Here, we report the first mouse knockout of this histone variant and its effects on fertility, nuclear organization, and gene expression. In view of the controversy related to the generation of off-target mutations by gene editing approaches, we test the specificity of TALENs by disrupting the H2A.B multi-copy gene family using only one pair of TALENs. We show that TALENs do display a high level of specificity since no off-target mutations are detected by bioinformatics analyses of exome sequences obtained from three consecutive generations of knockout mice and by Sanger DNA sequencing. Male H2A.B.3 knockout mice are subfertile and display an increase in the proportion of abnormal sperm and clogged seminiferous tubules. Significantly, a loss of proper RNA Pol II targeting to distinct transcription-splicing territories and changes to pre-mRNA splicing are observed.

CONCLUSION

We have produced the first H2A.B knockout mouse using the TALEN approach.

Gain-of-function IKBKB mutation causes human combined immune deficiency

Genetic mutations account for many devastating early onset immune deficiencies. In contrast, less severe and later onset immune diseases, including in patients with no prior family history, remain poorly understood. Whole exome sequencing in two cohorts of such patients identified a novel heterozygous de novo IKBKB missense mutation (c.607G>A) in two separate kindreds in whom probands presented with immune dysregulation, combined T and B cell deficiency, inflammation, and epithelial defects. IKBKB encodes IKK2, which activates NF-κB signaling. IKK2^V203I results in enhanced NF-κB signaling, as well as T and B cell functional defects. IKK2^V203 is a highly conserved residue, and to prove causation, we generated an accurate mouse model by introducing the precise orthologous codon change in Ikbkb using CRISPR/Cas9. Mice and humans carrying this missense mutation exhibit remarkably similar cellular and biochemical phenotypes. Accurate mouse models engineered by CRISPR/Cas9 can help characterize novel syndromes arising from de novo germline mutations and yield insight into pathogenesis.

A homozygous KAT2B variant modulates the clinical phenotype of ADD3 deficiency in humans and flies

Recent evidence suggests that the presence of more than one pathogenic mutation in a single patient is more common than previously anticipated. One of the challenges hereby is to dissect the contribution of each gene mutation, for which animal models such as Drosophila can provide a valuable aid. Here, we identified three families with mutations in ADD3, encoding for adducin-γ, with intellectual disability, microcephaly, cataracts and skeletal defects. In one of the families with additional cardiomyopathy and steroid-resistant nephrotic syndrome (SRNS), we found a homozygous variant in KAT2B, encoding the lysine acetyltransferase 2B, with impact on KAT2B protein levels in patient fibroblasts, suggesting that this second mutation might contribute to the increased disease spectrum. In order to define the contribution of ADD3 and KAT2B mutations for the patient phenotype, we performed functional experiments in the Drosophila model. We found that both mutations were unable to fully rescue the viability of the respective null mutants of the Drosophila homologs, hts and Gcn5, suggesting that they are indeed pathogenic in flies. While the KAT2B/Gcn5 mutation additionally showed a significantly reduced ability to rescue morphological and functional defects of cardiomyocytes and nephrocytes (podocyte-like cells), this was not the case for the ADD3 mutant rescue. Yet, the simultaneous knockdown of KAT2B and ADD3 synergistically impaired kidney and heart function in flies as well as the adhesion and migration capacity of cultured human podocytes, indicating that mutations in both genes may be required for the full clinical manifestation. Altogether, our studies describe the expansion of the phenotypic spectrum in ADD3 deficiency associated with a homozygous likely pathogenic KAT2B variant and thereby identify KAT2B as a susceptibility gene for kidney and heart disease in ADD3-associated disorders.

Impact of the G84E variant on HOXB13 gene and protein expression in formalin-fixed, paraffin-embedded prostate tumours

The HOXB13 G84E variant is associated with risk of prostate cancer (PCa), however the role this variant plays in PCa development is unknown. This study examined 751 cases, 450 relatives and 355 controls to determine the contribution of this variant to PCa risk in Tasmania and investigated HOXB13 gene and protein expression in tumours from nine G84E heterozygote variant and 13 wild-type carriers. Quantitative PCR and immunohistochemistry showed that HOXB13 gene and protein expression did not differ between tumour samples from variant and wild-type carriers. Allele-specific transcription revealed that two of seven G84E carriers transcribed both the variant and wild-type allele, while five carriers transcribed the wild-type allele. Methylation of surrounding CpG sites was lower in the variant compared to the wild-type allele, however overall methylation across the region was very low. Notably, tumour characteristics were less aggressive in the two variant carriers that transcribed the variant allele compared to the five that did not. This study has shown that HOXB13 expression does not differ between tumour tissue of G84E variant carriers and non-carriers. Intriguingly, the G84E variant allele was rarely transcribed in carriers, suggesting that HOXB13 expression may be driven by the wild-type allele in the majority of carriers.

Plasmacytoid dendritic cells appear inactive during sub-microscopic Plasmodium falciparum blood-stage infection, yet retain their ability to respond to TLR stimulation

Plasmacytoid dendritic cells (pDC) are activators of innate and adaptive immune responses that express HLA-DR, toll-like receptor (TLR) 7, TLR9 and produce type I interferons. The role of human pDC in malaria remains poorly characterised. pDC activation and cytokine production were assessed in 59 malaria-naive volunteers during experimental infection with 150 or 1,800 P. falciparum-parasitized red blood cells. Using RNA sequencing, longitudinal changes in pDC gene expression were examined in five adults before and at peak-infection. pDC responsiveness to TLR7 and TLR9 stimulation was assessed in-vitro. Circulating pDC remained transcriptionally stable with gene expression altered for 8 genes (FDR < 0.07). There was no upregulation of co-stimulatory molecules CD86, CD80, CD40, and reduced surface expression of HLA-DR and CD123 (IL-3R-α). pDC loss from the circulation was associated with active caspase-3, suggesting pDC apoptosis during primary infection. pDC remained responsive to TLR stimulation, producing IFN-α and upregulating HLA-DR, CD86, CD123 at peak-infection. In clinical malaria, pDC retained HLA-DR but reduced CD123 expression compared to convalescence. These data demonstrate pDC retain function during a first blood-stage P. falciparum exposure despite sub-microscopic parasitaemia downregulating HLA-DR. The lack of evident pDC activation in both early infection and malaria suggests little response of circulating pDC to infection.

Rapid Identification of Pathogenic Variants in Two Cases of Charcot-Marie-Tooth Disease by Gene-Panel Sequencing

Charcot-Marie-Tooth disease (CMT) is a common inherited peripheral neuropathy affecting up to 1 in 1214 of the general population with more than 60 nuclear genes implicated in its pathogenesis. Traditional molecular diagnostic pathways based on relative prevalence and clinical phenotyping are limited by long turnaround time, population-specific prevalence of causative variants and inability to assess multiple co-existing variants. In this study, a CMT gene panel comprising 27 genes was used to uncover the pathogenic mutations in two index patients. The first patient is a 15-year-old boy, born of consanguineous parents, who has had frequent trips and falls since infancy, and was later found to have inverted champagne bottle appearance of bilateral legs and foot drop. His elder sister is similarly affected. The second patient is a 37-year-old woman referred for pre-pregnancy genetic diagnosis. During early adulthood, she developed progressive lower limb weakness, difficulties in tip-toe walking and thinning of calf muscles. Both patients are clinically compatible with CMT, have undergone multiple genetic testings and have not previously received a definitive genetic diagnosis. Patients 1 and 2 were found to have pathogenic homozygous HSPB1:NM_001540:c.250G>A (p.G84R) variant and heterozygous GDAP1:NM_018972:c.358C>T (p.R120W) variant, respectively. Advantages and limitations of the current approach are discussed.

Resequencing array for gene variant detection in malignant hyperthermia and butyrylcholinestherase deficiency

Malignant hyperthermia (MH) and butyrylcholinestherase (BCHE) deficiency are two relevant pharmacogenetic disorders in anesthetic practice linked with sequence variants, the former in the RyR1 and CACNA1S genes, the latter in the BCHE gene. Genotyping for known pathogenic variants in these genes is useful to help identify susceptible individuals, and others may exist but remain unknown, because full-length sequence of these genes is, in general, not investigated. To facilitate this task, we developed a resequencing DNA array, the perioperative patient safety (POPS) array, to be able to screen the entire coding sequences of the RyR1, CACNA1S and BCHE genes. MH-susceptible individuals (n = 121) identified with the in vitro contracture test, the standard diagnostic tool for MH susceptibility, were genotyped with the arrays. Compared with capillary sequencing, call rates with the arrays could achieve 100% at maximal sensitivity, although to reduce false positive rates, sensitivity was adjusted to 0.85, 0.87 and 0.66 for RyR1, CACNA1S and BCHE respectively, with overall base call specificity exceeding 99%. Detection of 29 predetermined RyR1 variants in 44 individuals was successful in 97% of the cases, among them all 16 variants of established diagnostic value. In a trial application of the arrays, 21 MH-susceptible subjects with no known RyR1 or CACNA1S variants were screened, resulting in the discovery of new variants, all confirmed by capillary sequencing. In conclusion, arrays offer an efficient high-throughput alternative for diagnostic genotyping of candidate genes affecting MH susceptibility, BCHE deficiency and other neuromuscular disorders, simultaneously enabling a comprehensive search for rare variants in these genes.

Mutations in sphingosine-1-phosphate lyase cause nephrosis with ichthyosis and adrenal insufficiency

Steroid-resistant nephrotic syndrome (SRNS) causes 15% of chronic kidney disease cases. A mutation in 1 of over 40 monogenic genes can be detected in approximately 30% of individuals with SRNS whose symptoms manifest before 25 years of age. However, in many patients, the genetic etiology remains unknown. Here, we have performed whole exome sequencing to identify recessive causes of SRNS. In 7 families with SRNS and facultative ichthyosis, adrenal insufficiency, immunodeficiency, and neurological defects, we identified 9 different recessive mutations in SGPL1, which encodes sphingosine-1-phosphate (S1P) lyase. All mutations resulted in reduced or absent SGPL1 protein and/or enzyme activity. Overexpression of cDNA representing SGPL1 mutations resulted in subcellular mislocalization of SGPL1. Furthermore, expression of WT human SGPL1 rescued growth of SGPL1-deficient dpl1Δ yeast strains, whereas expression of disease-associated variants did not. Immunofluorescence revealed SGPL1 expression in mouse podocytes and mesangial cells. Knockdown of Sgpl1 in rat mesangial cells inhibited cell migration, which was partially rescued by VPC23109, an S1P receptor antagonist. In Drosophila, Sply mutants, which lack SGPL1, displayed a phenotype reminiscent of nephrotic syndrome in nephrocytes. WT Sply, but not the disease-associated variants, rescued this phenotype. Together, these results indicate that SGPL1 mutations cause a syndromic form of SRNS.

Suitability of Different Mapping Algorithms for Genome-Wide Polymorphism Scans with Pool-Seq Data

The cost-effectiveness of sequencing pools of individuals (Pool-Seq) provides the basis for the popularity and widespread use of this method for many research questions, ranging from unraveling the genetic basis of complex traits, to the clonal evolution of cancer cells. Because the accuracy of Pool-Seq could be affected by many potential sources of error, several studies have determined, for example, the influence of sequencing technology, the library preparation protocol, and mapping parameters. Nevertheless, the impact of the mapping tools has not yet been evaluated. Using simulated and real Pool-Seq data, we demonstrate a substantial impact of the mapping tools, leading to characteristic false positives in genome-wide scans. The problem of false positives was particularly pronounced when data with different read lengths and insert sizes were compared. Out of 14 evaluated algorithms novoalign, bwa mem and clc4 are most suitable for mapping Pool-Seq data. Nevertheless, no single algorithm is sufficient for avoiding all false positives. We show that the intersection of the results of two mapping algorithms provides a simple, yet effective, strategy to eliminate false positives. We propose that the implementation of a consistent Pool-Seq bioinformatics pipeline, building on the recommendations of this study, can substantially increase the reliability of Pool-Seq results, in particular when libraries generated with different protocols are being compared.

DeepSNVMiner: a sequence analysis tool to detect emergent, rare mutations in subsets of cell populations

Background. Massively parallel sequencing technology is being used to sequence highly diverse populations of DNA such as that derived from heterogeneous cell mixtures containing both wild-type and disease-related states. At the core of such molecule tagging techniques is the tagging and identification of sequence reads derived from individual input DNA molecules, which must be first computationally disambiguated to generate read groups sharing common sequence tags, with each read group representing a single input DNA molecule. This disambiguation typically generates huge numbers of reads groups, each of which requires additional variant detection analysis steps to be run specific to each read group, thus representing a significant computational challenge. While sequencing technologies for producing these data are approaching maturity, the lack of available computational tools for analysing such heterogeneous sequence data represents an obstacle to the widespread adoption of this technology. Results. Using synthetic data we successfully detect unique variants at dilution levels of 1 in a 1,000,000 molecules, and find DeeepSNVMiner obtains significantly lower false positive and false negative rates compared to popular variant callers GATK, SAMTools, FreeBayes and LoFreq, particularly as the variant concentration levels decrease. In a dilution series with genomic DNA from two cells lines, we find DeepSNVMiner identifies a known somatic variant when present at concentrations of only 1 in 1,000 molecules in the input material, the lowest concentration amongst all variant callers tested. Conclusions. Here we present DeepSNVMiner; a tool to disambiguate tagged sequence groups and robustly identify sequence variants specific to subsets of starting DNA molecules that may indicate the presence of a disease. DeepSNVMiner is an automated workflow of custom sequence analysis utilities and open source tools able to differentiate somatic DNA variants from artefactual sequence variants that likely arose during DNA amplification. The workflow remains flexible such that it may be customised to variants of the data production protocol used, and supports reproducible analysis through detailed logging and reporting of results. DeepSNVMiner is available for academic non-commercial research purposes at https://github.com/mattmattmattmatt/DeepSNVMiner.

Clinical correlation of calpain-1 and glypican-3 expression with gallbladder carcinoma

Gallbladder carcinoma (GBC) possesses a poor prognosis, which is primarily attributed to the lack of early and timely surgical intervention. Calpain-1 and glypican-3 have been implicated in the progression of various types of cancer. The present study aimed to detect the expression of calpain-1 and glypican-3 in GBC, and analyzed whether the expression levels of these proteins correlated with any clinicopathological variables. A total of 100 patients with GBC and 30 patients with cholecystitis who accepted surgical treatment were enrolled in the present study. Pathological and clinical data were obtained from all patients. The expression of calpain-1 and glypican-3 was detected in paraffin-embedded tissues by immunohistochemistry. Calpain-1 expression was manually assessed with an immunohistochemical H-score with a slight modification. Glypican-3 expression was assessed as negative and positive. The correlations between protein expression and clinicopathological characteristics, and the associations between the proteins were analyzed. All patients exhibited positive expression of calpain-1. Notably, the high expression rate of calpain-1 was significantly increased in patients with GBC, compared with patients with cholecystitis (32.0 vs. 6.7%; χ²=7.668; P=0.006), suggesting that calpain-1 expression may be associated with progression from cholecystitis to GBC. In addition, the positive rate of glypican-3 expression was 53.0% in patients with GBC and 63.3% in patients with cholecystitis, with no significant difference (χ²=0.997; P=0.318). Furthermore, the expression of calpain-1 and glypican-3 had no significant correlation with gender, age, degree of tumor differentiation and tumor-node-metastasis classification in patients with GBC. Notably, the expression of calpain-1 and glypican-3 displayed a significant positive correlation in patients with GBC (r=0.517; P<0.01), but a significantly negative correlation (r=-0.856; P<0.01) in patients with cholecystitis. In conclusion, calpain-1 expression may be associated with progression from cholecystitis to GBC. Combined detection of calpain-1 and glypican-3 may be beneficial for prognosis assessment of GBC.