Rare variant association on unrelated individuals in case-control studies using aggregation tests: existing methods and current limitations.

Over the past years, progress made in next-generation sequencing technologies and bioinformatics have sparked a surge in association studies. Especially, genome-wide association studies (GWASs) have demonstrated their effectiveness in identifying disease associations with common genetic variants. Yet, rare variants can contribute to additional disease risk or trait heterogeneity. Because GWASs are underpowered for detecting association with such variants, numerous statistical methods have been recently proposed. Aggregation tests collapse multiple rare variants within a genetic region (e.g. gene, gene set, genomic loci) to test for association. An increasing number of studies using such methods successfully identified trait-associated rare variants and led to a better understanding of the underlying disease mechanism. In this review, we compare existing aggregation tests, their statistical features and scope of application, splitting them into the five classical classes: burden, adaptive burden, variance-component, omnibus and other. Finally, we describe some limitations of current aggregation tests, highlighting potential direction for further investigations.

Epilepsy with faint capillary malformation or reticulated telangiectasia associated with mosaic AKT3 pathogenic variants.

Capillary malformations (CMs) are the most common type of vascular anomalies, affecting around 0.3% of newborns. They are usually caused by somatic pathogenic variants in GNAQ or GNA11. PIK3CA and PIK3R1, part of the phosphoinositide 3-kinase-protein kinase B-mammalian target of rapamycin pathway, are mutated in fainter CMs such as diffuse CM with overgrowth and megalencephaly CM. In this study, we present two young patients with a CM-like phenotype associated with cerebral anomalies and severe epilepsy. Pathogenic variants in PIK3CA and PIK3R1, as well as GNAQ and GNA11, were absent in affected cutaneous tissue biopsies. Instead, we identified two somatic pathogenic variants in the AKT3 gene. Subsequent analysis of the DNA obtained from surgically resected brain tissue of one of the two patients confirmed the presence of the AKT3 variant. Focal cortical dysplasia was also detected in this patient. Genetic analysis thus facilitated workup to reach a precise diagnosis for these patients, associating the vascular anomaly with the neurological symptoms. This study underscores the importance of searching for additional signs and symptoms to guide the diagnostic workup, especially in cases with atypical vascular malformations. In addition, it strongly emphasizes the significance of genotype-phenotype correlation studies in guiding clinicians' informed decision-making regarding patient care.

Excalibur: A new ensemble method based on an optimal combination of aggregation tests for rare-variant association testing for sequencing data.

The development of high-throughput next-generation sequencing technologies and large-scale genetic association studies produced numerous advances in the biostatistics field. Various aggregation tests, i.e. statistical methods that analyze associations of a trait with multiple markers within a genomic region, have produced a variety of novel discoveries. Notwithstanding their usefulness, there is no single test that fits all needs, each suffering from specific drawbacks. Selecting the right aggregation test, while considering an unknown underlying genetic model of the disease, remains an important challenge. Here we propose a new ensemble method, called Excalibur, based on an optimal combination of 36 aggregation tests created after an in-depth study of the limitations of each test and their impact on the quality of result. Our findings demonstrate the ability of our method to control type I error and illustrate that it offers the best average power across all scenarios. The proposed method allows for novel advances in Whole Exome/Genome sequencing association studies, able to handle a wide range of association models, providing researchers with an optimal aggregation analysis for the genetic regions of interest.

Ureteropelvic junction obstruction with primary lymphoedema associated with CELSR1 variants.

BACKGROUND: Primary lymphoedema (PL) is a chronic, debilitating disease caused by developmental and functional defects of the lymphatic system. It is marked by an accumulation of interstitial fluid, fat and tissue fibrosis. There is no cure. More than 50 genes and genetic loci have been linked to PL. We sought to study systematically cell polarity signalling protein Cadherin Epidermal Growth Factor Laminin G Seven-pass G-type Receptor 1 (CELSR1) variants linked to PL. METHODS: We investigated 742 index patients from our PL cohort using exome sequencing. RESULTS: We identified nine variants predicted to cause CELSR1 loss of function. Four of them were tested for nonsense-mediated mRNA decay, but none was observed. Most of the truncated CELSR1 proteins would lack the transmembrane domain, if produced. The affected individuals had puberty/late-onset PL on lower extremities. The variants had a statistically significant difference in penetrance between female patients (87%) and male patients (20%). Eight variant carriers had a kidney anomaly, mostly in the form of ureteropelvic junction obstruction, which has not been associated with CELSR1 before. CELSR1 is located in the 22q13.3 deletion locus of the Phelan-McDermid syndrome. As variable renal defects are often seen in patients with the Phelan-McDermid syndrome, CELSR1 may be the long-sought gene for the renal defects. CONCLUSION: PL associated with a renal anomaly suggests a CELSR1-related cause.

RASA1 mosaic mutations in patients with capillary malformation-arteriovenous malformation.

BACKGROUND: Capillary malformation-arteriovenous malformation is an autosomal dominant disorder, characterised by capillary malformations and increased risk of fast-flow vascular malformations, caused by loss-of-function mutations in the RASA1 or EPHB4 genes. Around 25% of the patients do not seem to carry a germline mutation in either one of these two genes. Even if other genes could be involved, some individuals may have mutations in the known genes that escaped detection by less sensitive techniques. We tested the hypothesis that mosaic mutations could explain some of previously negative cases. METHODS: DNA was extracted from peripheral blood lymphocytes, saliva or vascular malformation tissues from four patients. RASA1 and EPHB4 coding regions and exon/intron boundaries were analysed by targeted custom gene panel sequencing. A second panel and/or Sanger sequencing were used to confirm the identified mutations. RESULTS: Four distinct mosaic RASA1 mutations, with an allele frequency ranging from 3% to 25%, were identified in four index patients with classical capillary malformation-arteriovenous malformation phenotype. Three mutations were known, one was novel. In one patient, a somatic second hit was also identified. One index case had three affected children, illustrating that the mosaicism was also present in the germline. CONCLUSION: This study shows that RASA1 mosaic mutations can cause capillary malformation-arteriovenous malformation. Thus, highly sensitive sequencing techniques should be considered as diagnostic tools, especially for patients with no family history. Even low-level mosaicism can cause the classical phenotype and increased risk for offspring. In addition, our study further supports the second-hit pathophysiological mechanism to explain the multifocality of vascular lesions in this disorder.

Tumor sequencing is useful to refine the analysis of germline variants in unexplained high-risk breast cancer families.

BACKGROUND: Multigene panels are routinely used to assess for predisposing germline mutations in families at high breast cancer risk. The number of variants of unknown significance thereby identified increases with the number of sequenced genes. We aimed to determine whether tumor sequencing can help refine the analysis of germline variants based on second somatic genetic events in the same gene. METHODS: Whole-exome sequencing (WES) was performed on whole blood DNA from 70 unrelated breast cancer patients referred for genetic testing and without a BRCA1, BRCA2, TP53, or CHEK2 mutation. Rare variants were retained in a list of 735 genes. WES was performed on matched tumor DNA to identify somatic second hits (copy number alterations (CNAs) or mutations) in the same genes. Distinct methods (among which immunohistochemistry, mutational signatures, homologous recombination deficiency, and tumor mutation burden analyses) were used to further study the role of the variants in tumor development, as appropriate. RESULTS: Sixty-eight patients (97%) carried at least one germline variant (4.7 ± 2.0 variants per patient). Of the 329 variants, 55 (17%) presented a second hit in paired tumor tissue. Of these, 53 were CNAs, resulting in tumor enrichment (28 variants) or depletion (25 variants) of the germline variant. Eleven patients received variant disclosure, with clinical measures for five of them. Seven variants in breast cancer-predisposing genes were considered not implicated in oncogenesis. One patient presented significant tumor enrichment of a germline variant in the oncogene ERBB2, in vitro expression of which caused downstream signaling pathway activation. CONCLUSION: Tumor sequencing is a powerful approach to refine variant interpretation in cancer-predisposing genes in high-risk breast cancer patients. In this series, the strategy provided clinically relevant information for 11 out of 70 patients (16%), adapted to the considered gene and the familial clinical phenotype.

SLC13A3 variants cause acute reversible leukoencephalopathy and α-ketoglutarate accumulation.

OBJECTIVE: SLC13A3 encodes the plasma membrane Na+ /dicarboxylate cotransporter 3, which imports inside the cell 4 to 6 carbon dicarboxylates as well as N-acetylaspartate (NAA). SLC13A3 is mainly expressed in kidney, in astrocytes, and in the choroid plexus. We describe two unrelated patients presenting with acute, reversible (and recurrent in one) neurological deterioration during a febrile illness. Both patients exhibited a reversible leukoencephalopathy and a urinary excretion of α-ketoglutarate (αKG) that was markedly increased and persisted over time. In one patient, increased concentrations of cerebrospinal fluid NAA and dicarboxylates (including αKG) were observed. Extensive workup was unsuccessful, and a genetic cause was suspected. METHODS: Whole exome sequencing (WES) was performed. Our teams were connected through GeneMatcher. RESULTS: WES analysis revealed variants in SLC13A3. A homozygous missense mutation (p.Ala254Asp) was found in the first patient. The second patient was heterozygous for another missense mutation (p.Gly548Ser) and an intronic mutation affecting splicing as demonstrated by reverse transcriptase polymerase chain reaction performed in muscle tissue (c.1016 + 3A > G). Mutations and segregation were confirmed by Sanger sequencing. Functional studies performed on HEK293T cells transiently transfected with wild-type and mutant SLC13A3 indicated that the missense mutations caused a marked reduction in the capacity to transport αKG, succinate, and NAA. INTERPRETATION: SLC13A3 deficiency causes acute and reversible leukoencephalopathy with marked accumulation of αKG. Urine organic acids (especially αKG and NAA) and SLC13A3 mutations should be screened in patients presenting with unexplained reversible leukoencephalopathy, for which SLC13A3 deficiency is a novel differential diagnosis. ANN NEUROL 2019;85:385-395.

Analysing ambiguities in trypanosomatids taxonomy by barcoding.

BACKGROUND: Biodiversity screens and phylogenetic studies are dependent on reliable DNA sequences in public databases. Biological collections possess vouchered specimens with a traceable history. Therefore, DNA sequencing of samples available at institutional collections can greatly contribute to taxonomy, and studies on evolution and biodiversity. METHODS: We sequenced part of the glycosomal glyceraldehyde phosphate dehydrogenase (gGAPDH) and the SSU rRNA (V7/V8) genes from 102 trypanosomatid cultures, which are available on request at www.colprot.fiocruz.br. OBJECTIVE: The main objective of this work was to use phylogenetic inferences, using the obtained DNA sequences and those from representatives of all Trypanosomatidae genera, to generate phylogenetic trees that can simplify new isolates screenings. FINDINGS: A DNA sequence is provided for the first time for several isolates, the phylogenetic analysis allowed the classification or reclassification of several specimens, identification of candidates for new genera and species, as well as the taxonomic validation of several deposits. MAIN CONCLUSIONS: This survey aimed at presenting a list of validated species and their associated DNA sequences combined with a short historical overview of each isolate, which can support taxonomic and biodiversity research and promote culture collections.

Loss of ADAMTS3 activity causes Hennekam lymphangiectasia-lymphedema syndrome 3.

Primary lymphedema is due to developmental and/or functional defects in the lymphatic system. It may affect any part of the body, with predominance for the lower extremities. Twenty-seven genes have already been linked to primary lymphedema, either isolated, or as part of a syndrome. The proteins that they encode are involved in VEGFR3 receptor signaling. They account for about one third of all primary lymphedema cases, underscoring the existence of additional genetic factors. We used whole-exome sequencing to investigate the underlying cause in a non-consanguineous family with two children affected by lymphedema, lymphangiectasia and distinct facial features. We discovered bi-allelic missense mutations in ADAMTS3. Both were predicted to be highly damaging. These amino acid substitutions affect well-conserved residues in the prodomain and in the peptidase domain of ADAMTS3. In vitro, the mutant proteins were abnormally processed and sequestered within cells, which abolished proteolytic activation of pro-VEGFC. VEGFC processing is also affected by CCBE1 mutations that cause the Hennekam lymphangiectasia-lymphedema syndrome syndrome type1. Our data identifies ADAMTS3 as a novel gene that can be mutated in individuals affected by the Hennekam syndrome. These patients have distinctive facial features similar to those with mutations in CCBE1. Our results corroborate the recent in vitro and murine data that suggest a close functional interaction between ADAMTS3 and CCBE1 in triggering VEGFR3 signaling, a cornerstone for the differentiation and function of lymphatic endothelial cells.

PDGFRB gain-of-function mutations in sporadic infantile myofibromatosis.

Infantile myofibromatosis is one of the most prevalent soft tissue tumors of infancy and childhood. Multifocal nodules with visceral lesions are associated with a poor prognosis. A few familial cases have been linked to mutations in various genes including PDGFRB. In this study, we sequenced PDGFRB, which encodes a receptor tyrosine kinase, in 16 cases of myofibromatosis or solitary myofibroma. Mutations in the coding sequence of PDGFRB were identified in 6 out of 8 patients with the sporadic multicentric form of the disease and in 1 out of 8 patients with isolated myofibroma. Two patients had the same mutation in multiple separated lesions. By contrast, a third patient had three different PDGFRB mutations in the three nodules analyzed. Mutations were located in the transmembrane, juxtamembrane and kinase domains of the receptor. We showed that these mutations activated receptor signaling in the absence of ligand and transformed fibroblasts. In one case, a weakly-activating germline variant was associated with a stronger somatic mutation, suggesting a two-hit model for familial myofibromatosis. Furthermore, the mutant receptors were sensitive to the tyrosine kinase inhibitor imatinib, except D850V, which was inhibited by dasatinib and ponatinib, suggesting a targeted therapy for severe myofibromatosis. In conclusion, we identified gain-of-function PDGFRB mutations in the majority of multifocal infantile myofibromatosis cases, shedding light on the mechanism of disease development, which is reminiscent of multifocal venous malformations induced by TIE2 mutations. Our results provide a genetic test to facilitate diagnosis, and preclinical data for development of molecular therapies.