Candidate gene discovery in autoimmunity by using extreme phenotypes, next generation sequencing and whole exome capture

Large-scale whole-exome sequencing analyses identified protein-coding variants associated with immune-mediated diseases in 350,770 adults

The genetic contribution of protein-coding variants to immune-mediated diseases (IMDs) remains underexplored. Through whole exome sequencing of 40 IMDs in 350,770 UK Biobank participants, we identified 162 unique genes in 35 IMDs, among which 124 were novel genes. Several genes, including FLG which is associated with atopic dermatitis and asthma, showed converging evidence from both rare and common variants. 91 genes exerted significant effects on longitudinal outcomes (interquartile range of Hazard Ratio: 1.12-5.89). Mendelian randomization identified five causal genes, of which four were approved drug targets (CDSN, DDR1, LTA, and IL18BP). Proteomic analysis indicated that mutations associated with specific IMDs might also affect protein expression in other IMDs. For example, DXO (celiac disease-related gene) and PSMB9 (alopecia areata-related gene) could modulate CDSN (autoimmune hypothyroidism-, psoriasis-, asthma-, and Graves' disease-related gene) expression. Identified genes predominantly impact immune and biochemical processes, and can be clustered into pathways of immune-related, urate metabolism, and antigen processing. Our findings identified protein-coding variants which are the key to IMDs pathogenesis and provided new insights into tailored innovative therapies.

Multilayer concept of autoimmune mechanisms and manifestations in inborn errors of immunity: Relevance for precision therapy

Autoimmunity in inborn errors of immunity (IEIs) has a multifactorial pathogenesis and develops subsequent to a genetic predisposition in conjunction with gene regulation, environmental modifiers, and infectious triggers. On the basis of incremental data availability owing to upfront application of omics technologies, a more granular and dynamic view of mechanisms and manifestations is warranted. Here, we present a comprehensive novel concept of autoimmunity in IEIs that considers multiple layers of interdependent elements and connects 101 causative genes or deletions according to the quality of the allelic variants with 47 molecular pathways and 22 immune effector mechanisms. Furthermore, we list 50 resulting manifestations together with the corresponding Human Phenotype Ontology terms and review the types and frequencies of the most relevant clinical presentations. When all of its elements are taken together, this concept (1) extends the historical anatomic view of central versus peripheral tolerance toward multiple interdependent mechanisms of immune tolerance, (2) delineates the mechanisms underlying the protean clinical manifestations, and thereby, (3) points toward the most suitable precision therapy for autoimmunity in IEIs. The multilayer concept of autoimmune mechanisms and manifestations in IEIs will facilitate research design and provide clinical guidance on the use of precision medicine irrespective of the data depth available in each health care scenario.

Multiple autoimmune syndrome: Clinical, immunological and genotypic characterization

INTRODUCTION

The existence of subphenotypes common to several autoimmune diseases (AIDs) suggests a shared physiopathology - autoimmune tautology. Multiple Autoimmune Syndrome (MAS) - the coexistence of three or more AIDs in one person-, best illustrates that polyautoimmunity is more than a coincidence.

OBJECTIVES

Characterize and compare the monoautoimmune and MAS patients. Understand if clustering of AIDs leads to differences in disease severity, autoantibodies expression or genetic polymorphisms that could be markers for polyautoimmunity.

METHODS

Currently adult patients were selected from unit cohort. MAS was assumed when ≥3 AIDs were present. 343 patients were included after exclusion criteria: having two AIDs or undetermined diagnosis. Clinical and immunological data were collected from medical files. HLA-DRB1 was genotyped by PCR-SSP methodology and PTPN22(rs2476601) polymorphisms by TaqMan Real Time PCR. Data were analysed using Chi-Square, Fisher's exact tests and logistic regression. Odds ratios (OR) and 95% confidence intervals were calculated.

RESULTS

In comparison with control population: ELEVATED FREQUENCIES: HLA-DRB1*03 in study cohort (OR=3.68,p<0.001) and in monoautoimmune SLE (OR=2.79,p<0.001) and SjS (OR=8.27,p<0.001); HLA-DRB1*15 in monoautoimmune SjS (OR=2.39,p = 0.011); HLA-DRB1*16 in MAS SLE (OR=2.67,p = 0.031); PTPN22_T in all groups except monoautoimmune SjS and triple positive systemic MAS.

DIMINISHED FREQUENCIES

HLA-DRB1*11 in study cohort (OR=0.57,p = 0.013), in MAS SLE (OR=0.39,p = 0.031) and monoautoimmune SjS (OR=0.10,p = 0.005); HLA-DRB1*13 in study cohort (OR=0.52,p = 0.001) and in monoautoimmune SLE (OR=0.53,p = 0.009) and SjS (OR=0.38,p = 0.031); HLA-DRB1*14 in study cohort (OR=0.32,p = 0.013) and monoautoimmune SLE (OR=0.21,p = 0.021); SLE group: HLA-DRB1*07 frequency was higher in monoautoimmune patients (OR=0.43,p = 0.023). MAS patients had significantly more NPSLE (OR=2.99,p<0.001), subacute cutaneous lesions (OR=2.30,p = 0.037), muscle&tendon (OR=2.00,p = 0.045), and haematological (OR=3.18,p = 0.006) involvement and Raynaud's (OR=2.94,p<0.001). SjS group: MAS patients had more frequently cryoglobulins (OR=2.96,p = 0.030), low complement (OR=2.43,p = 0.030) and Raynaud's (OR=4.38,p<0.001); monoautoimmune patients had more parotid enlargement (OR=0.12,p<0.001). APS group: MAS patients had more non-thrombotic manifestations (OR=4.69,p = 0.020) and Raynaud's (OR=9.12,p<0.001). Triple positive systemic MAS (SLE+SjS+APS) had more frequently severe kidney involvement (OR=11.67,p = 0.021) and CNS thrombosis (OR=4.44,p = 0.009). Anti-U1RNP increased frequency was transversally attributable to MAS.

CONCLUSIONS

The coexistence of AIDs contributes to a more severe disease course. We confirmed previously established genetic risk and protection factors and suggest a new protective one - HLA-DRB1*14. HLA-DRB1*07 and anti-U1RNP could be markers for mono and polyautoimmunity, respectively; HLA-DRB1*13 could be a predictor for vascular risk in patients with multiple AIDs. PTPN22(rs2476601) polymorphism could be associated with less severe disease.

Muscle and bone characteristics of a Chinese family with spinal muscular atrophy, lower extremity predominant 1 (SMALED1) caused by a novel missense DYNC1H1 mutation

BACKGROUND

Spinal muscular atrophy, lower extremity predominant (SMALED) is a type of non-5q spinal muscular atrophy characterised by weakness and atrophy of lower limb muscles without sensory abnormalities. SMALED1 can be caused by dynein cytoplasmic 1 heavy chain 1 (DYNC1H1) gene variants. However, the phenotype and genotype of SMALED1 may overlap with those of other neuromuscular diseases, making it difficult to diagnose clinically. Additionally, bone metabolism and bone mineral density (BMD) in patients with SMALED1 have never been reported.

METHODS

We investigated a Chinese family in which 5 individuals from 3 generations had lower limb muscle atrophy and foot deformities. Clinical manifestations and biochemical and radiographic indices were analysed, and mutational analysis was performed by whole-exome sequencing (WES) and Sanger sequencing.

RESULTS

A novel mutation in exon 4 of the DYNC1H1 gene (c.587T > C, p.Leu196Ser) was identified in the proband and his affected mother by WES. Sanger sequencing confirmed that the proband and 3 affected family members were carriers of this mutation. As leucine is a hydrophobic amino acid and serine is hydrophilic, the hydrophobic interaction resulting from mutation of amino acid residue 196 could influence the stability of the DYNC1H1 protein. Leg muscle magnetic resonance imaging of the proband revealed severe atrophy and fatty infiltration, and electromyographic recordings showed chronic neurogenic impairment of the lower extremities. Bone metabolism markers and BMD of the proband were all within normal ranges. None of the 4 patients had experienced fragility fractures.

CONCLUSION

This study identified a novel DYNC1H1 mutation and expands the spectrum of phenotypes and genotypes of DYNC1H1-related disorders. This is the first report of bone metabolism and BMD in patients with SMALED1.

Structural Protein Effects Underpinning Cognitive Developmental Delay of the PURA p.Phe233del Mutation Modelled by Artificial Intelligence and the Hybrid Quantum Mechanics–Molecular Mechanics Framework

A whole-exome capture and next-generation sequencing was applied to an 11 y/o patient with a clinical history of congenital hypotonia, generalized motor and cognitive neurodevelopmental delay, and severe cognitive deficit, and without any identifiable Syndromic pattern, and to her parents, we disclosed a de novo heterozygous pathogenic mutation, c.697_699del p.Phe233del (rs786204835)(ACMG classification PS2, PM1, PM2, PP5), harbored in the PURA gene (MIM*600473) (5q31.3), associated with Autosomal Dominant Mental Retardation 31 (MIM # 616158). We used the significant improvement in the accuracy of protein structure prediction recently implemented in AlphaFold that incorporates novel neural network architectures and training procedures based on the evolutionary, physical, and geometric constraints of protein structures. The wild-type (WT) sequence and the mutated sequence, missing the Phe233, were reconstructed. The predicted local Distance Difference Test (lDDT) for the PURAwt and the PURA–Phe233del showed that the occurrence of the Phe233del affects between 220–320 amino acids. The distortion in the PURA structural conformation in the ~5 Å surrounding area after the p.Phe233del produces a conspicuous disruption of the repeat III, where the DNA and RNA helix unwinding capability occurs. PURA Protein–DNA docking corroborated these results in an in silico analysis that showed a loss of the contact of the PURA–Phe233del III repeat domain model with the DNA. Together, (i) the energetic and stereochemical, (ii) the hydropathic indexes and polarity surfaces, and (iii) the hybrid Quantum Mechanics–Molecular Mechanics (QM–MM) analyses of the PURA molecular models demarcate, at the atomic resolution, the specific surrounding region affected by these mutations and pave the way for future cell-based functional analysis. To the best of our knowledge, this is the first report of a de novo mutation underpinning a PURA syndrome in a Latin American patient and highlights the importance of predicting the molecular effects in protein structure using artificial intelligence algorithms and molecular and atomic resolution stereochemical analyses.

Ataxia with Ocular Apraxia Type 1 (AOA1) (APTX, W279* Mutation): Neurological, Neuropsychological, and Molecular Outlining of a Heterogenous Phenotype in Four Colombian Siblings

Hereditary ataxias are a group of devastating neurological disorders that affect coordination of gait and are often associated with poor coordination of hands, speech, and eye movements. Ataxia with ocular apraxia type 1 (AOA1) (OMIM: 606,350.0006) is characterized by slowly progressive symptoms of childhood-onset and pathogenic mutations in APTX; the only known cause underpinning AOA1. APTX encodes the protein aprataxin, composed of three domains sharing homology with proteins involved in DNA damage, signaling, and repair. We present four siblings from an endogamic family in a rural, isolated town of Colombia with ataxia and ocular apraxia of childhood-onset and confirmed molecular diagnosis of AOA1, homozygous for the W279* p.Trp279Ter mutation. We predicted the mutated APTX with AlphaFold to demonstrate the effects of this stop-gain mutation that deletes three beta helices encoded by amino acid 270 to 339 rescinding the C2H2-type zinc fingers (Znf) (C2H2 Znf) DNA-binding, the DNA-repair domain, and the whole 3D structure of APTX. All siblings exhibited different ages of onset (4, 6, 8, and 11 years old) and heterogeneous patterns of dysarthria (ranging from absence to mild-moderate dysarthria). Neuropsychological evaluation showed no neurocognitive impairment in three siblings, but one sibling showed temporospatial disorientation, semantic and phonologic fluency impairment, episodic memory affection, constructional apraxia, moderate anomia, low executive function, and symptoms of depression. To our knowledge, this report represents the most extensive series of siblings affected with AOA1 in Latin America, and the genetic analysis completed adds important knowledge to outline this family's disease and general complex phenotype of hereditary ataxias.

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.

GMEPS: a fast and efficient likelihood approach for genome-wide mediation analysis under extreme phenotype sequencing

Due to many advantages such as higher statistical power of detecting the association of genetic variants in human disorders and cost saving, extreme phenotype sequencing (EPS) is a rapidly emerging study design in epidemiological and clinical studies investigating how genetic variations associate with complex phenotypes. However, the investigation of the mediation effect of genetic variants on phenotypes is strictly restrictive under the EPS design because existing methods cannot well accommodate the non-random extreme tails sampling process incurred by the EPS design. In this paper, we propose a likelihood approach for testing the mediation effect of genetic variants through continuous and binary mediators on a continuous phenotype under the EPS design (GMEPS). Besides implementing in EPS design, it can also be utilized as a general mediation analysis procedure. Extensive simulations and two real data applications of a genome-wide association study of benign ethnic neutropenia under EPS design and a candidate-gene study of neurocognitive performance in patients with sickle cell disease under random sampling design demonstrate the superiority of GMEPS under the EPS design over widely used mediation analysis procedures, while demonstrating compatible capabilities under the general random sampling framework.

Whole-Exome Sequencing in Patients Affected by Stevens-Johnson Syndrome and Toxic Epidermal Necrolysis Reveals New Variants Potentially Contributing to the Phenotype

Background

Adverse drug reactions (ADRs) are frequent occurring events that can essentially be defined as harmful or unpleasant symptoms secondary to the use of a medicinal product. ADRs involve a wide spectrum of clinical manifestations ranging from minor itching and rash to life-threatening reactions. Stevens-Johnson syndrome (SJS) and toxic epidermal necrolysis (TEN) are rare ADRs. SJS-TEN may be considered a polygenic pathology due to additive/epistatic effects caused by sequence variants in numerous genes. Next-generation sequencing (NGS) represents a potentially interesting exploration tool in such scenario as it facilitates the simultaneous analysis of large genomic regions and genes at affordable cost.

Methods

The present study has involved using whole-exome sequencing (WES) for the first time on SJS-TEN patients. It involved robust and innovative multistep bioinformatics analysis focusing on 313 candidate genes potentially participating in the disease’s aetiology, specific drugs’ metabolism and gene regulation.

Results

We identified combinations of frequently occurring and rare variants that may contribute to the disease’s pathogenesis. Depending on the specific drug being taken, different variants (and alleles) in NAT2, CYP2D8, CYP2B6, ABCC2, UGT2B7 and TCF3 were identified as coherent candidates representing potential future markers for SJS-TEN.

Conclusion

The present study proposed and has described (for the first time) a large-scale genomic analysis of patients affected by SJS-TEN. The genes and variants identified represent relevant candidates potentially participating in the disease’s pathogenesis. Corroborating that proposed by others, we found that complex combinations of frequently occurring and rare variants participating in particular drug metabolism molecular cascades could be associated with the phenotype. TCF3 TF may be considered a coherent candidate for SJS-TEN that should be analysed in new cohorts of patients having ADRs.

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.

Identifying new potential genetic biomarkers for HELLP syndrome using massive parallel sequencing

BACKGROUND

Preeclampsia (PE) is a frequently occurring multisystemic disease affecting ~5% of pregnancies. PE patients may develop HELLP syndrome (haemolysis, elevated liver enzymes, and low platelet), a mother and foetus life-threatening condition. Research into HELLP's genetic origin has been relatively unsuccessful, mainly because normal placental function and blood pressure regulation involve the fine-regulation of hundreds of genes.

OBJECTIVE

To identify new genes and mutations constituting potential biomarkers for HELLP syndrome.

STUDY DESIGN

The present case-control study involved whole-exome sequencing of 79 unrelated HELLP women. Candidate variants were screened in a control population constituted by 176 individuals. Stringent bioinformatics filters were used for selecting potentially etiological sequence variants in a subset of 487 genes. We used robust in silico mutation modelling for predicting the potential effect on protein structure.

RESULTS

We identified numerous sequence variants in genes related to angiogenesis/coagulation/blood pressure regulation, cell differentiation/communication/adhesion, cell cycle and transcriptional gene regulation, extracellular matrix biology, lipid metabolism and immunological response. Five sequence variants generated premature stop codons in genes playing an essential role in placental physiology (STOX1, PDGFD, IGF2, MMP1 and DNAH11). Six variants (ERAP1- p.Ile915Thr, ERAP2- p.Leu837Ser, COMT-p.His192Gln, CSAD-p.Pro418Ser, CDH1- p.Ala298Thr and CCR2-p.Met249Lys) led to destabilisation of protein structure as they had significant energy and residue interaction-related changes. We identified at least two mutations in 57% of patients, arguing in favour of a polygenic origin for the HELLP syndrome.

CONCLUSION

Our results provide novel evidence regarding PE/HELLP's genetic origin, leading to new biomarkers, having potential clinical usefulness, being proposed.

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.

Familial Alzheimer's Disease and Recessive Modifiers

Alzheimer’s disease (AD) is progressive brain disorder that affects ~ 50 million people worldwide and has no current effective treatment. AD age of onset (ADAOO) has shown to be critical for the identification of genes that modify the appearance of AD signs and symptoms in a specific population. We clinically characterized and whole-exome genotyped 71 individuals with AD from the Paisa genetic isolate, segregating the (PSEN1) E280A dominant fully penetrant mutation, and analyzed the potential recessive effects of ~ 50,000 common functional genomic variants to the ADAOO. Standard quality control and filtering procedures were applied, and recessive single- and multi-locus linear mixed-effects models were used. We identified genetic variants in the SLC9C1, CSN1S1, and LOXL4 acting recessively to delay ADAOO up to ~ 11, ~ 6, and ~ 9 years on average, respectively. In contrast, the CC recessive genotype in marker DHRS4L2-rs2273946 accelerates ADAOO by ~ 8 years. This study, reports new recessive variants modifying ADAOO in PSEN1 E280A mutation carriers. This set of genes are implicated in important biological processes and molecular functions commonly affected by genes associated with the etiology of AD such as APP, APOE, and CLU. Future functional studies using modern techniques such as induced pluripotent stem cells will allow a better understanding of the over expression and down regulation of these recessive modifier variants and hence the pathogenesis of AD. These results are important for prediction of AD and ultimately, substantial to develop new therapeutic strategies for individuals at risk or affected by AD.

New insights into the genetics and epigenetics of systemic sclerosis

Systemic sclerosis (SSc) is a severe autoimmune disease that is characterized by vascular abnormalities, immunological alterations and fibrosis of the skin and internal organs. The results of genetic studies in patients with SSc have revealed statistically significant genetic associations with disease manifestations and progression. Nevertheless, genetic susceptibility to SSc is moderate, and the functional consequences of genetic associations remain only partially characterized. A current hypothesis is that, in genetically susceptible individuals, epigenetic modifications constitute the driving force for disease initiation. As epigenetic alterations can occur years before fibrosis appears, these changes could represent a potential link between inflammation and tissue fibrosis. Epigenetics is a fast-growing discipline, and a considerable number of important epigenetic studies in SSc have been published in the past few years that span histone post-translational modifications, DNA methylation, microRNAs and long non-coding RNAs. This Review describes the latest insights into genetic and epigenetic contributions to the pathogenesis of SSc and aims to provide an improved understanding of the molecular pathways that link inflammation and fibrosis. This knowledge will be of paramount importance for the development of medicines that are effective in treating or even reversing tissue fibrosis.

Genetic Variation Underpinning ADHD Risk in a Caribbean Community

Attention Deficit Hyperactivity Disorder (ADHD) is a highly heritable and prevalent neurodevelopmental disorder that frequently persists into adulthood. Strong evidence from genetic studies indicates that single nucleotide polymorphisms (SNPs) harboured in the ADGRL3 (LPHN3), SNAP25, FGF1, DRD4, and SLC6A2 genes are associated with ADHD. We genotyped 26 SNPs harboured in genes previously reported to be associated with ADHD and evaluated their potential association in 386 individuals belonging to 113 nuclear families from a Caribbean community in Barranquilla, Colombia, using family-based association tests. SNPs rs362990-SNAP25 (T allele; p = 2.46 × 10⁻⁴), rs2282794-FGF1 (A allele; p = 1.33 × 10⁻²), rs2122642-ADGRL3 (C allele, p = 3.5 × 10⁻²), and ADGRL3 haplotype CCC (markers rs1565902-rs10001410-rs2122642, OR = 1.74, P_permuted = 0.021) were significantly associated with ADHD. Our results confirm the susceptibility to ADHD conferred by SNAP25, FGF1, and ADGRL3 variants in a community with a significant African American component, and provide evidence supporting the existence of specific patterns of genetic stratification underpinning the susceptibility to ADHD. Knowledge of population genetics is crucial to define risk and predict susceptibility to disease.

Variants at potential loci associated with Sjogren's syndrome in Koreans: A genetic association study

Sjogren's syndrome (SS), a chronic autoimmune disease, typically causes or involves inflammation in the salivary and lacrimal glands. Although recent genetic association studies have contributed to the discovery of SS susceptible genes, few studies have reported on the Korean population. Here, we did a genetic association study of SS in Korean patients using whole-exome sequencing data of 15 patients and 100 healthy controls. In addition to confirming previously described SS susceptibility loci MSH5 (p = 1.67 × 10-5) and RELN (p = 4.91 × 10-6), we also validated PRAMEF13 (p = 2.28 × 10-5), TARBP1 (p = 1.87 × 10-5), UGT2B28 (p = 1.33 × 10-5), TRBV5-6 (p = 2.27 × 10-5) and NAPB (p = 3.73 × 10-5) as novel susceptibility loci for SS. Furthermore, we identified UGT2B28, TARBP1 and PRAMEF13 as associated with human immune function. These findings may provide useful insight into to the pathways and pathogenesis contributing to SS susceptibility in the Korean population.

Whole-Transcriptome Sequencing: a Powerful Tool for Vascular Tissue Engineering and Endothelial Mechanobiology

Among applicable high-throughput techniques in cardiovascular biology, whole-transcriptome sequencing is of particular use. By utilizing RNA that is isolated from virtually all cells and tissues, the entire transcriptome can be evaluated. In comparison with other high-throughput approaches, RNA sequencing is characterized by a relatively low-cost and large data output, which permits a comprehensive analysis of spatiotemporal variation in the gene expression profile. Both shear stress and cyclic strain exert hemodynamic force upon the arterial endothelium and are considered to be crucial determinants of endothelial physiology. Laminar blood flow results in a high shear stress that promotes atheroresistant endothelial phenotype, while a turbulent, oscillatory flow yields a pathologically low shear stress that disturbs endothelial homeostasis, making respective arterial segments prone to atherosclerosis. Severe atherosclerosis significantly impairs blood supply to the organs and frequently requires bypass surgery or an arterial replacement surgery that requires tissue-engineered vascular grafts. To provide insight into patterns of gene expression in endothelial cells in native or bioartificial arteries under different biomechanical conditions, this article discusses applications of whole-transcriptome sequencing in endothelial mechanobiology and vascular tissue engineering.

Recent advances in understanding autoimmune thyroid disease: the tallest tree in the forest of polyautoimmunity

Autoimmune thyroid disease (AITD) is often observed together with other autoimmune diseases. The coexistence of two or more autoimmune diseases in the same patient is referred to as polyautoimmunity, and AITD is the autoimmune disease most frequently involved. The occurrence of polyautoimmunity has led to the hypothesis that the affected patients suffer from a generalized dysregulation of their immune system. The present review summarizes recent discoveries unravelling the immunological mechanisms involved in autoimmunity, ranging from natural autoimmunity to disease-specific autoimmunity. Furthermore, the clinical grounds for considering AITD in a setting of polyautoimmunity are explored. A better understanding of these may pave the way for designing new treatment modalities targeting the underlying immune dysregulation when AITD appears in the context of polyautoimmunity.

A complex intragenic rearrangement of ERCC8 in Chinese siblings with Cockayne syndrome

Cockayne syndrome is an autosomal recessive disorder principally characterized by postnatal growth failure and progressive neurological dysfunction, due primarily to mutations in ERCC6 and ERCC8. Here, we report our diagnostic experience for two patients in a Chinese family suspected on clinical grounds to have Cockayne syndrome. Using multiple molecular techniques, including whole exome sequencing, array comparative genomic hybridization and quantitative polymerase chain reaction, we identified compound heterozygosity for a maternal splicing variant (chr5:60195556, NM_000082:c.618-2A > G) and a paternal complex deletion/inversion/deletion rearrangement removing exon 4 of ERCC8, confirming the suspected pathogenesis in these two subjects. Microhomology (TAA and AGCT) at the breakpoints indicated that microhomology-mediated FoSTeS events were involved in this complex ERCC8 rearrangement. This diagnostic experience illustrates the value of high-throughput genomic technologies combined with detailed phenotypic assessment in clinical genetic diagnosis.

Mutations modifying sporadic Alzheimer's disease age of onset

The identification of mutations modifying the age of onset (AOO) in Alzheimer's disease (AD) is crucial for understanding the natural history of AD and, therefore, for early interventions. Patients with sporadic AD (sAD) from a genetic isolate in the extremes of the AOO distribution were whole-exome genotyped. Single- and multi-locus linear mixed-effects models were used to identify functional variants modifying AOO. A posteriori enrichment and bioinformatic analyses were applied to evaluate the non-random clustering of the associate variants to physiopathological pathways involved in AD. We identified more than 20 pathogenic, genome-wide statistically significant mutations of major modifier effect on the AOO. These variants are harbored in genes implicated in neuron apoptosis, neurogenesis, inflammatory processes linked to AD, oligodendrocyte differentiation, and memory processes. This set of new genes harboring these mutations could be of importance for prediction, follow-up and eventually as therapeutical targets of AD. © 2016 Wiley Periodicals, Inc.

APOE*E2 allele delays age of onset in PSEN1 E280A Alzheimer's disease

Alzheimer's disease (AD) age of onset (ADAOO) varies greatly between individuals, with unique causal mutations suggesting the role of modifying genetic and environmental interactions. We analyzed ~50 000 common and rare functional genomic variants from 71 individuals of the 'Paisa' pedigree, the world's largest pedigree segregating a severe form of early-onset AD, who were affected carriers of the fully penetrant E280A mutation in the presenilin-1 (PSEN1) gene. Affected carriers with ages at the extremes of the ADAOO distribution (30s-70s age range), and linear mixed-effects models were used to build single-locus regression models outlining the ADAOO. We identified the rs7412 (APOE*E2 allele) as a whole exome-wide ADAOO modifier that delays ADAOO by ~12 years (β=11.74, 95% confidence interval (CI): 8.07-15.41, P=6.31 × 10(-8), PFDR=2.48 × 10(-3)). Subsequently, to evaluate comprehensively the APOE (apolipoprotein E) haplotype variants (E1/E2/E3/E4), the markers rs7412 and rs429358 were genotyped in 93 AD affected carriers of the E280A mutation. We found that the APOE*E2 allele, and not APOE*E4, modifies ADAOO in carriers of the E280A mutation (β=8.24, 95% CI: 4.45-12.01, P=3.84 × 10(-5)). Exploratory linear mixed-effects multilocus analysis suggested that other functional variants harbored in genes involved in cell proliferation, protein degradation, apoptotic and immune dysregulation processes (i.e., GPR20, TRIM22, FCRL5, AOAH, PINLYP, IFI16, RC3H1 and DFNA5) might interact with the APOE*E2 allele. Interestingly, suggestive evidence as an ADAOO modifier was found for one of these variants (GPR20) in a set of patients with sporadic AD from the Paisa genetic isolate. This is the first study demonstrating that the APOE*E2 allele modifies the natural history of AD typified by the age of onset in E280A mutation carriers. To the best of our knowledge, this is the largest analyzed sample of patients with a unique mutation sharing uniform environment. Formal replication of our results in other populations and in other forms of AD will be crucial for prediction, follow-up and presumably developing new therapeutic strategies for patients either at risk or affected by AD.

Definition of mutations in polyautoimmunity

OBJECTIVES

Familial autoimmunity and polyautoimmunity represent extreme phenotypes ideal for identifying major genomic variants contributing to autoimmunity. Whole exome sequencing (WES) and linkage analysis are well suited for this purpose due to its strong resolution upon familial segregation patterns of functional protein coding and splice variants. The primary objective of this study was to identify potentially autoimmune causative variants using WES data from extreme pedigrees segregating polyautoimmunity phenotypes.

METHODS

DNA of 47 individuals across 10 extreme pedigrees, ascertained from probands affected with polyautoimmunity and familial autoimmunity, were selected for WES. Variant calls were obtained through Genome Analysis Toolkit. Filtration and prioritization framework to identify mutation(s) were applied, and later implemented for genetic linkage analysis. Sanger sequencing corroborated variants with significant linkage.

RESULTS

Novel and mostly rare variants harbored in SRA1, MLL4, ABCB8, DHX34 and PLAUR showed significant linkage (LOD scores are >3.0). The strongest signal was in SRA1, with a LOD score of 5.48. Network analyses indicated that SRA1, PLAUR and ABCB8 contribute to regulation of apoptotic processes.

CONCLUSIONS

Novel and rare variants in genetic linkage with polyautoimmunity were identified throughout WES. Genes harboring these variants might be major players of autoimmunity.

Identification of a de novo DYNC1H1 mutation via WES according to published guidelines

De novo mutations that contribute to rare Mendelian diseases, including neurological disorders, have been recently identified. Whole-exome sequencing (WES) has become a powerful tool for the identification of inherited and de novo mutations in Mendelian diseases. Two important guidelines were recently published regarding the investigation of causality of sequence variant in human disease and the interpretation of novel variants identified in human genome sequences. In this study, a family with supposed movement disorders was sequenced via WES (including the proband and her unaffected parents), and a standard investigation and interpretation of the identified variants was performed according to the published guidelines. We identified a novel de novo mutation (c.2327C > T, p.P776L) in DYNC1H1 gene and confirmed that it was the causal variant. The phenotype of the affected twins included delayed motor milestones, pes cavus, lower limb weakness and atrophy, and a waddling gait. Electromyographic (EMG) recordings revealed typical signs of chronic denervation. Our study demonstrates the power of WES to discover the de novo mutations associated with a neurological disease on the whole exome scale, and guidelines to conduct WES studies and interpret of identified variants are a preferable option for the exploration of the pathogenesis of rare neurological disorders.

Novel and rare functional genomic variants in multiple autoimmune syndrome and Sjögren's syndrome

Background

Multiple autoimmune syndrome (MAS), an extreme phenotype of autoimmune disorders, is a very well suited trait to tackle genomic variants of these conditions. Whole exome sequencing (WES) is a widely used strategy for detection of protein coding and splicing variants associated with inherited diseases.

Methods

The DNA of eight patients affected by MAS [all of whom presenting with Sjögren’s syndrome (SS)], four patients affected by SS alone and 38 unaffected individuals, were subject to WES. Filters to identify novel and rare functional (pathogenic–deleterious) homozygous and/or compound heterozygous variants in these patients and controls were applied. Bioinformatics tools such as the Human gene connectome as well as pathway and network analysis were applied to test overrepresentation of genes harbouring these variants in critical pathways and networks involved in autoimmunity.

Results

Eleven novel and rare functional variants were identified in cases but not in controls, harboured in: MACF1, KIAA0754, DUSP12, ICA1, CELA1, LRP1/STAT6, GRIN3B, ANKLE1, TMEM161A, and FKRP. These were subsequently subject to network analysis and their functional relatedness to genes already associated with autoimmunity was evaluated. Notably, the LRP1/STAT6 novel mutation was homozygous in one MAS affected patient and heterozygous in another. LRP1/STAT6 disclosed the strongest plausibility for autoimmunity. LRP1/STAT6 are involved in extracellular and intracellular anti-inflammatory pathways that play key roles in maintaining the homeostasis of the immune system. Further; networks, pathways, and interaction analyses showed that LRP1 is functionally related to the HLA-B and IL10 genes and it has a substantial impact within immunological pathways and/or reaction to bacterial and other foreign proteins (phagocytosis, regulation of phospholipase A2 activity, negative regulation of apoptosis and response to lipopolysaccharides). Further, ICA1 and STAT6 were also closely related to AIRE and IRF5, two very well known autoimmunity genes.

Conclusions

Novel and rare exonic mutations that may account for autoimmunity were identified. Among those, the LRP1/STAT6 novel mutation has the strongest case for being categorised as potentially causative of MAS given the presence of intriguing patterns of functional interaction with other major genes shaping autoimmunity.

Electronic supplementary material

The online version of this article (doi:10.1186/s12967-015-0525-x) contains supplementary material, which is available to authorized users.