Case-only exome sequencing and complex disease susceptibility gene discovery: study design considerations

Case report: Genetic analysis of a novel intronic inversion variant in the SPTB gene associated with hereditary spherocytosis

Background: Hereditary spherocytosis (HS) is a congenital haemolytic anaemia attributed to dysregulation or abnormal quantities of erythrocyte membrane proteins. Currently, the most common erythrocytic gene, spectrin β (SPTB), variants are located in exons and give rise to mRNA defects. However, the genetic characteristics and pathogenic mechanisms of SPTB intronic variants are not completely understood. This study aimed to analyse a rare intronic inversion variant in the SPTB gene associated with HS, and explore the impact of the variant on SPTB mRNA splicing. Method: The clinical manifestations of the patient were summarised and analysed for spherocytosis phenotype diagnosis. The pathogenic variant was identified in the proband using targeted next-generation and Sanger sequencing. RNA sequencing was performed to analyse whether SPTB gene splicing and expression were affected. Results: Targeted next-generation sequencing identified a novel disease-associated intronic inversion variant of the SPTB gene in the proband. The inversion variant was located between intron 19 and 20, and contained the entire exon 20 and partial sequences of adjacent introns. Sanger sequencing confirmed that the intronic inversion variant only appeared in the genome of the proband, not in his parents. RNA sequencing revealed that the variant could result in the skipping of exon 20 and reduced expression of SPTB mRNA. Conclusion: This study identifies a rare intronic inversion variant in the SPTB gene associated with hereditary spherocytosis. The pathogenic variant can lead to exon 20 skipping and decreased SPTB gene expression. This finding has not been previously reported in any literature. This study can expand the intronic variant spectrum of the SPTB gene, deepen our understanding of HS pathogenesis, and contribute to the genetic diagnosis and clinical management of patients.

Citrullinemia type II accompanied by mental derangement combined with multidrug resistance 3 decrease, case report

Background

Here we report a rare case of citrullinemia type II (CTLN2) accompanied by mental derangement with a deficiency of multidrug resistance 3 (MDR3) in the liver.

Case presentation

The clinical data of a 17-year-old girl were collected. Liver puncture was performed, and hepatic expression of MDR3 was determined by immunohistochemistry. Serum amino acids of the patient and her parents wwere determined by a chemical isotope labeling liquid chromatography-mass spectrometry (CIL LC-MS). Genetic mutations of ABCB4 and SLC25A13 were screened by whole-exome sequencing. Immunohistochemical analysis showed a remarkably lower expression of MDR3. Mutation in ABCB4 gene was not found and whole-exome sequencing revealed the SLC25A13 mutation 852-855 del. Elevated serum levels of citrulline, homocitrulline, and homoarginine in the patient and her mother were found.

Conclusions

We reported a rare case of CTLN2 combined with MDR3 deficiency, without mutation of ABCB4. The link between MDR3 down-expression and CTLN2 warrants further investigation. Meanwhile, clinicians need to further rule out the possibility of CTLN2 if MDR3 decreases in adolescent patients with mental disorders and abnormal liver function.

Uncovering the Contribution of Moderate-Penetrance Susceptibility Genes to Breast Cancer by Whole-Exome Sequencing and Targeted Enrichment Sequencing of Candidate Genes in Women of European Ancestry

Rare variants in at least 10 genes, including BRCA1, BRCA2, PALB2, ATM, and CHEK2, are associated with increased risk of breast cancer; however, these variants, in combination with common variants identified through genome-wide association studies, explain only a fraction of the familial aggregation of the disease. To identify further susceptibility genes, we performed a two-stage whole-exome sequencing study. In the discovery stage, samples from 1528 breast cancer cases enriched for breast cancer susceptibility and 3733 geographically matched unaffected controls were sequenced. Using five different filtering and gene prioritization strategies, 198 genes were selected for further validation. These genes, and a panel of 32 known or suspected breast cancer susceptibility genes, were assessed in a validation set of 6211 cases and 6019 controls for their association with risk of breast cancer overall, and by estrogen receptor (ER) disease subtypes, using gene burden tests applied to loss-of-function and rare missense variants. Twenty genes showed nominal evidence of association (p-value < 0.05) with either overall or subtype-specific breast cancer. Our study had the statistical power to detect susceptibility genes with effect sizes similar to ATM, CHEK2, and PALB2, however, it was underpowered to identify genes in which susceptibility variants are rarer or confer smaller effect sizes. Larger sample sizes would be required in order to identify such genes.

Development of vascular disease models to explore disease causation and pathomechanisms of rare vascular diseases

As the field of medicine is striving forward heralded by a new era of next-generation sequencing (NGS) and integrated technologies such as bioprinting and biological material development, the utility of rare monogenetic vascular disease modeling in this landscape is starting to emerge. With their genetic simplicity and broader applicability, these patient-specific models are at the forefront of modern personalized medicine. As a collective, rare diseases are a significant burden on global healthcare systems, and rare vascular diseases make up a significant proportion of this. High costs are due to a lengthy diagnostic process, affecting all ages from infants to adults, as well as the severity and chronic nature of the disease. Their complex nature requires sophisticated disease models and integrated approaches involving multidisciplinary teams. Here, we review these emerging vascular disease models, how they contribute to our understanding of the pathomechanisms in rare vascular diseases and provide useful platforms for therapeutic discovery.

Variant Enrichment Analysis to Explore Pathways Functionality in Complex Autoinflammatory Skin Disorders through Whole Exome Sequencing Analysis

The challenge of unravelling the molecular basis of multifactorial disorders nowadays cannot rely just on association studies searching for potential causative variants shared by groups of patients and not present in healthy individuals; indeed, association studies have as a main limitation the lack of information on the interactions between the disease-causing variants. Thus, new genomic analysis tools focusing on disrupted pathways rather than associated gene variants are required to better understand the complexity of a disease. Therefore, we developed the Variant Enrichment Analysis (VEA) workflow, a tool applicable for whole exome sequencing data, able to find differences between the numbers of genetic variants in a given pathway in comparison with a reference dataset. In this study, we applied VEA to discover novel pathways altered in patients with complex autoinflammatory skin disorders, namely PASH (n = 9), 3 of whom are overlapping with SAPHO) and PAPASH (n = 3). With this approach we have been able to identify pathways related to neutrophil and endothelial cells homeostasis/activations, as disrupted in our patients. We hypothesized that unregulated neutrophil transendothelial migration could elicit increased neutrophil infiltration and tissue damage. Based on our findings, VEA, in our experimental dataset, allowed us to predict novel pathways impaired in subjects with autoinflammatory skin disorders.

Coding Variants Coupled With Rapid Modeling in Zebrafish Implicate Dynein Genes, dnaaf1 and zmynd10, as Adolescent Idiopathic Scoliosis Candidate Genes

Adolescent idiopathic scoliosis (AIS) is the most common pediatric spine disorder affecting ∼3% of children worldwide. Human genetic studies suggest a complex polygenic disease model for AIS with large genetic and phenotypic heterogeneity. However, the overall genetic etiology of AIS remains poorly understood. To identify additional AIS susceptibility loci, we performed whole-exome sequencing (WES) on a cohort of 195 Southern Chinese AIS patients. Bioinformatics analysis identified 237 novel rare variants associated with AIS, located in 232 new susceptibility loci. Enrichment analysis of these variants revealed 10 gene families associated with our AIS cohort. We screened these gene families by comparing our candidate gene list with IS candidate genes in the Human Phenotype Ontology (HPO) database and previous reported studies. Two candidate gene families, axonemal dynein and axonemal dynein assembly factors, were retained for their associations with ciliary architecture and function. The damaging effects of candidate variants in dynein genes dnali1, dnah1, dnaaf, and zmynd10, as well as in one fibrillin-related gene tns1, were functionally analyzed in zebrafish using targeted CRISPR/Cas9 screening. Knockout of two candidate genes, dnaaf1 or zmynd10, recapitulated scoliosis in viable adult zebrafish. Altogether, our results suggest that the disruption of one or more dynein-associated factors may correlate with AIS susceptibility in the Southern Chinese population.

Monozygotic Twins Concordant for Common Variable Immunodeficiency: Strikingly Similar Clinical and Immune Profile Associated With a Polygenic Burden

Monozygotic twins provide a unique opportunity to better understand complex genetic diseases and the relative contribution of heritable factors in shaping the immune system throughout life. Common Variable Immunodeficiency Disorders (CVID) are primary antibody defects displaying wide phenotypic and genetic heterogeneity, with monogenic transmission accounting for only a minority of the cases. Here, we report a pair of monozygotic twins concordant for CVID without a family history of primary immunodeficiency. They featured a remarkably similar profile of clinical manifestations and immunological alterations at diagnosis (established at age 37) and along the subsequent 15 years of follow-up. Interestingly, whole-exome sequencing failed to identify a monogenic cause for CVID, but unraveled a combination of heterozygous variants, with a predicted deleterious impact. These variants were found in genes involved in relevant immunological pathways, such as JUN, PTPRC, TLR1, ICAM1, and JAK3. The potential for combinatorial effects translating into the observed disease phenotype is inferred from their roles in immune pathways, namely in T and B cell activation. The combination of these genetic variants is also likely to impose a significant constraint on environmental influences, resulting in a similar immunological phenotype in both twins, despite exposure to different living conditions. Overall, these cases stress the importance of integrating NGS data with clinical and immunological phenotypes at the single-cell level, as provided by multi-dimensional flow-cytometry, in order to understand the complex genetic landscape underlying the vast majority of patients with CVID, as well as those with other immunodeficiencies.

Association score testing for rare variants and binary traits in family data with shared controls

Genome-wide association studies have been an important approach used to localize trait loci, with primary focus on common variants. The multiple rare variant-common disease hypothesis may explain the missing heritability remaining after accounting for identified common variants. Advances of sequencing technologies with their decreasing costs, coupled with methodological advances in the context of association studies in large samples, now make the study of rare variants at a genome-wide scale feasible. The resurgence of family-based association designs because of their advantage in studying rare variants has also stimulated more methods development, mainly based on linear mixed models (LMMs). Other tests such as score tests can have advantages over the LMMs, but to date have mainly been proposed for single-marker association tests. In this article, we extend several score tests (χcorrected2, WQLS, and SKAT) to the multiple variant association framework. We evaluate and compare their statistical performances relative with the LMM. Moreover, we show that three tests can be cast as the difference between marker allele frequencies (AFs) estimated in each of the group of affected and unaffected subjects. We show that these tests are flexible, as they can be based on related, unrelated or both related and unrelated subjects. They also make feasible an increasingly common design that only sequences a subset of affected subjects (related or unrelated) and uses for comparison publicly available AFs estimated in a group of healthy subjects. Finally, we show the great impact of linkage disequilibrium on the performance of all these tests.

Genetic and experimental evidence for the involvement of the CD6 lymphocyte receptor in psoriasis

Psoriasis is a chronic inflammatory skin disease with a strong genetic background and is triggered by environmental factors. Available evidence supports CD6, a lymphocyte surface receptor mostly expressed by T cells, as a putative target in autoimmunity. Accordingly, a humanized anti-CD6 antibody has been assayed for the treatment of certain autoimmune disorders, including psoriasis. Here, we present novel evidence in mice and humans for a direct involvement of CD6 in psoriasis pathophysiology. First, an attenuated form of imiquimod-induced psoriasis-like skin inflammation was demonstrated in CD6-deficient mice, as deduced from lower epidermal thickness and local reduced production of pro-inflammatory cytokines, namely, interleukin-17A. Thus, isolated CD4⁺CD62L⁺ T cells from CD6-deficient mice displayed decreased in vitro T-helper type 17 polarization. Second, a statistically significant association between CD6 single-nucleotide polymorphisms (rs17824933, rs11230563 and rs12360861) and more severe forms of psoriasis was demonstrated in a cohort of 304 patients at three public hospitals from the metropolitan area of Barcelona. Taken together, these results provide new supportive evidence of the contribution of the CD6 lymphocyte receptor in psoriasis at both experimental and clinical levels.

Evidence from a familial case suggests maternal inheritance of primary biliary cholangitis

Primary biliary cholangitis (PBC) is an idiopathic autoimmune liver disease characterized by chronic cholestasis and destruction of the intrahepatic bile ducts. Similar to other autoimmune diseases, the pathogenesis of PBC is considered to be a complex etiologic phenomenon involving the interaction of genetic and environmental factors. Although a number of common variants associated with PBC have been reported from genome-wide association studies, a precise genetic mechanism underlying PBC has yet to be identified. Here, we describe a family with four sisters who were diagnosed with PBC. After the diagnosis of the index patient who was in an advanced stage of PBC, one sister presented with acute hepatitis, and two sisters were subsequently diagnosed with PBC. Notably, one half-sister with a different mother exhibited no evidence of PBC following clinical investigation. Our report suggests the possibility of a maternal inheritance of PBC susceptibility. Moreover, judging from the high-penetrance of the disease observed in this family, we inferred that a pathogenic genetic variant might be the cause of PBC development. We describe a family that exhibited diverse clinical presentations of PBC that included asymptomatic stages with mildly increased liver enzyme levels and symptomatic stages with acute hepatitis or advanced liver fibrosis. Additional studies are needed to investigate the role of genetic factors in the pathogenesis of this rare autoimmune disease.

Germline CDKN2A/P16INK4A mutations contribute to genetic determinism of sarcoma

BACKGROUND

Sarcomas are rare mesenchymal malignancies whose pathogenesis is poorly understood; both environmental and genetic risk factors could contribute to their aetiology.

METHODS AND RESULTS

We performed whole-exome sequencing (WES) in a familial aggregation of three individuals affected with soft-tissue sarcoma (STS) without TP53 mutation (Li-Fraumeni-like, LFL) and found a shared pathogenic mutation in CDKN2A tumour suppressor gene. We searched for individuals with sarcoma among 474 melanoma-prone families with a CDKN2A-/+ genotype and for CDKN2A mutations in 190 TP53-negative LFL families where the index case was a sarcoma. Including the initial family, eight independent sarcoma cases carried a germline mutation in the CDKN2A/p16INK4A gene. In five out of seven formalin-fixed paraffin-embedded sarcomas, heterozygosity was lost at germline CDKN2A mutations sites demonstrating complete loss of function. As sarcomas are rare in CDKN2A/p16INK4A carriers, we searched in constitutional WES of nine carriers for potential modifying rare variants and identified three in platelet-derived growth factor receptor (PDGFRA) gene. Molecular modelling showed that two never-described variants could impact the PDGFRA extracellular domain structure.

CONCLUSION

Germline mutations in CDKN2A/P16INK4A, a gene known to predispose to hereditary melanoma, pancreatic cancer and tobacco-related cancers, account also for a subset of hereditary sarcoma. In addition, we identified PDGFRA as a candidate modifier gene.

Discovery of Variants Underlying Host Susceptibility to Virus Infection Using Whole-Exome Sequencing

The clinical course of any viral infection greatly differs in individuals. This variation results from various viral, host, and environmental factors. The identification of host genetic factors influencing inter-individual variation in susceptibility to several pathogenic viruses has tremendously increased our understanding of the mechanisms and pathways required for immunity. Next-generation sequencing of whole exomes represents a powerful tool in biomedical research. In this chapter, we briefly introduce whole-exome sequencing in the context of genetic approaches to identify host susceptibility genes to viral infections. We then describe general aspects of the workflow for whole-exome sequence analysis together with the tools and online resources that can be used to identify and annotate variant calls, and then prioritize them for their potential association to phenotypes of interest.

Systematic evaluation of underlying defects in DNA repair as an approach to case-only assessment of familial prostate cancer

Risk assessment for prostate cancer is challenging due to its genetic heterogeneity. In this study, our goal was to develop an operational framework to select and evaluate gene variants that may contribute to familial prostate cancer risk. Drawing on orthogonal sources, we developed a candidate list of genes relevant to prostate cancer, then analyzed germline exomes from 12 case-only prostate cancer patients from high-risk families to identify patterns of protein-damaging gene variants. We described an average of 5 potentially disruptive variants in each individual and annotated them in the context of public databases representing human variation. Novel damaging variants were found in several genes of relevance to prostate cancer. Almost all patients had variants associated with defects in DNA damage response. Many also had variants linked to androgen signaling. Treatment of primary T-lymphocytes from these prostate cancer patients versus controls with DNA damaging agents showed elevated levels of the DNA double strand break (DSB) marker γH2AX (p < 0.05), supporting the idea of an underlying defect in DNA repair. This work suggests the value of focusing on underlying defects in DNA damage in familial prostate cancer risk assessment and demonstrates an operational framework for exome sequencing in case-only prostate cancer genetic evaluation.

Identification of novel hereditary cancer genes by whole exome sequencing

Whole exome sequencing (WES) provides a powerful tool for medical genetic research. Several dozens of WES studies involving patients with hereditary cancer syndromes have already been reported. WES led to breakthrough in understanding of the genetic basis of some exceptionally rare syndromes; for example, identification of germ-line SMARCA4 mutations in patients with ovarian hypercalcemic small cell carcinomas indeed explains a noticeable share of familial aggregation of this disease. However, studies on common cancer types turned out to be more difficult. In particular, there is almost a dozen of reports describing WES analysis of breast cancer patients, but none of them yet succeeded to reveal a gene responsible for the significant share of missing heritability. Virtually all components of WES studies require substantial improvement, e.g. technical performance of WES, interpretation of WES results, mode of patient selection, etc. Most of contemporary investigations focus on genes with autosomal dominant mechanism of inheritance; however, recessive and oligogenic models of transmission of cancer susceptibility also need to be considered. It is expected that the list of medically relevant tumor-predisposing genes will be rapidly expanding in the next few years.

Pharmacologic Therapy of Diabetes and Overall Cancer Risk and Mortality: A Meta-Analysis of 265 Studies

Different anti-diabetic medications (ADMs) may modify cancer risk and mortality in patients with diabetes. We conducted a systematic review and meta-analysis to estimate the magnitude of association and quality of supporting evidence for each ADM. A total of 265 studies (44 cohort studies, 39 case-control studies, and 182 randomized controlled trials (RCT)) were identified, involving approximately 7.6 million and 137,540 patients with diabetes for observational studies and RCTs, respectively. The risk of bias overall was moderate. Meta-analysis demonstrated that the use of metformin or thiazolidinediones was associated with a lower risk of cancer incidence (RR = 0.86, 95% CI 0.83-0.90, I² = 88.61%; RR = 0.93, 95% CI 0.91-0.96, I² = 0.00% respectively). On the other hand, insulin, sulfonylureas and alpha glucosidase inhibitor use was associated with an increased risk of cancer incidence (RR = 1.21, 95% CI 1.08-1.36, I² = 96.31%; RR = 1.20, 95% CI 1.13-1.27, I² = 95.02%; RR = 1.10, 95% CI 1.05-1.15, I² = 0.00% respectively). Use of other types of ADMs was not significantly associated with cancer risk. This study indicates that some ADMs may modify the risk of cancer in individuals with diabetes. Knowledge of this risk may affect the choice of ADM in individuals concerned about cancer or at increased risk for cancer.

Exome-Wide Association Study Identifies New Low-Frequency and Rare UGT1A1 Coding Variants and UGT1A6 Coding Variants Influencing Serum Bilirubin in Elderly Subjects: A Strobe Compliant Article

Genome-wide association studies (GWASs) have identified loci contributing to total serum bilirubin level. However, no exome-wide approaches have been performed to address this question. Using exome-wide approach, we assessed the influence of protein-coding variants on unconjugated, conjugated, and total serum bilirubin levels in a well-characterized cohort of 773 ambulatory elderly subjects from Italy. Coding variants were replicated in 227 elderly subjects from the same area. We identified 4 missense rare (minor allele frequency, MAF < 0.5%) and low-frequency (MAF, 0.5%-5%) coding variants located in the first exon of the UGT1A1 gene, which encodes for the substrate-binding domain (rs4148323 [MAF = 0.06%; p.Gly71Arg], rs144398951 [MAF = 0.06%; p.Ile215Val], rs35003977 [MAF = 0.78%; p.Val225Gly], and rs57307513 [MAF = 0.06%; p.Ser250Pro]). These variants were in strong linkage disequilibrium with 3 intronic UGT1A1 variants (rs887829, rs4148325, rs6742078), which were significantly associated with total bilirubin level (P = 2.34 × 10(-34), P = 7.02 × 10(-34), and P = 8.27 × 10(-34)), as well as unconjugated, and conjugated bilirubin levels. We also identified UGT1A6 variants in association with total (rs6759892, p.Ser7Ala, P = 1.98 × 10(-26); rs2070959, p.Thr181Ala, P = 2.87 × 10(-27); and rs1105879, p.Arg184Ser, P = 3.27 × 10(-29)), unconjugated, and conjugated bilirubin levels. All UGT1A1 intronic variants (rs887829, rs6742078, and rs4148325) and UGT1A6 coding variants (rs6759892, rs2070959, and rs1105879) were significantly associated with gallstone-related cholecystectomy risk. The UGT1A6 variant rs2070959 (p.Thr181Ala) was associated with the highest risk of gallstone-related cholecystectomy (OR, 4.58; 95% CI, 1.58-13.28; P = 3.21 × 10(-3)). Using an exome-wide approach we identified coding variants on UGT1A1 and UGT1A6 genes in association with serum bilirubin level and hyperbilirubinemia risk in elderly subjects. UGT1A1 intronic single-nucleotide polymorphisms (SNPs) (rs6742078, rs887829, rs4148324) serve as proxy markers for the low-frequency and rare UGT1A1 variants, thereby providing mechanistic explanation to the relationship between UGT1A1 intronic SNPs and the UGT1A1 enzyme activity. UGT1A1 and UGT1A6 variants might be potentially associated with gallstone-related cholecystectomy risk.

Zinc transporter genes and urological cancers: integrated analysis suggests a role for ZIP11 in bladder cancer

Although zinc transporters were shown to play roles in the development of prostate, bladder, and renal cancer, no study has evaluated the genetic variants in zinc transporter genes with risk of urological cancers. A candidate gene association study using genome-wide association study (GWAS) datasets was conducted for variants in 24 zinc transporter genes. Genotypes were analyzed using logistic regression models adjusted for covariates. The function of identified variants was assessed by using the Encyclopedia of DNA Elements (ENCODE). We further evaluated tumors for somatic change of the implicated gene(s) and the associations between identified variants and patient survival from data in The Cancer Genome Atlas (TCGA). A ZIP11 variant, rs8081059, was significantly associated with increased risk of renal cell carcinoma (odds ratios (OR) = 1.28, 95 % confidence intervals (CI) (1.13-1.45), p = 0.049). No zinc transporter variants were associated with prostate cancer risk. Four variants within ZIP11 were significantly associated with bladder cancer risk: rs11871756 (OR = 1.43, 95 % CI (1.24-1.63), p = 0.0002), rs11077654 (OR = 0.76, 95 % CI (0.68-0.85), p = 0.001), rs9913017 (OR = 0.76, 95 % CI (0.68-0.85), p = 0.002), and rs4969054 (OR = 0.78, 95 % CI (0.69-0.88), p = 0.02); the three protective variants were co-located and highly correlated. These variants were located within predicted transcribed or enhancer regions. Among the 253 bladder cancer patients in TCGA, two had tumors that contained deleterious missense mutations in ZIP11. Moreover, rs11077654 was significantly associated with survival of bladder cancer patients (p = 0.046). In conclusion, zinc transporter gene, ZIP11, may play an important role in bladder cancer. Further studies of the gene are warranted.