Rare Risk Variants Identification by Identity-by-Descent Mapping and Whole-Exome Sequencing Implicates Neuronal Development Pathways in Schizophrenia and Bipolar Disorder

Schizophrenia (SCZ) and bipolar disorder (BPD) are highly heritable disorders with an estimated co-heritability of 68%. Hundreds of common alleles have been implicated, but recently a role for rare, high-penetrant variants has been also suggested in both disorders. This study investigated a familial cohort of SCZ and BPD patients from a closed population sample, where the high recurrence of the disorders and the homogenous genetic background indicate a possible enrichment in rare risk alleles. A total of 230 subjects (161 cases, 22 unaffected relatives, and 47 controls) were genetically investigated through an innovative strategy that integrates identity-by-descent (IBD) mapping and whole-exome sequencing (WES). IBD analysis allowed to track high-risk haplotypes (IBD_risk) shared exclusively by multiple patients from different families and possibly carrying the most penetrant alleles. A total of 444 non-synonymous sequence variants, of which 137 disruptive, were identified in IBD_risk haplotypes by WES. Interestingly, gene sets previously implicated in SCZ (i.e., post-synaptic density (PSD) proteins, voltage-gated calcium channels (VGCCs), and fragile X mental retardation protein (FMRP) targets) were found significantly enriched in genes carrying IBD_risk variants. Further, IBD_risk variants were preferentially affecting genes involved in the extracellular matrix (ECM) biology and axon guidance processes which appeared to be functionally connected in the pathway-derived meta-network analysis. Results thus confirm rare risk variants as key factors in SCZ and BPD pathogenesis and highlight a role for the development of neuronal connectivity in the etiology of both disorders.     

Whole Exome Sequencing in Patients with the Cuticular Drusen Subtype of Age-Related Macular Degeneration

Age-related macular degeneration (AMD) is the leading cause of irreversible blindness in elderly people worldwide. Cuticular drusen (CD) is a clinical subtype of AMD, which typically displays an earlier age at onset, and has a strong genetic component. Genetic studies support a role for rare sequence variants in CD susceptibility, and rare sequence variants in the CFH gene have been identified in 8.8% of CD cases. To further explore the role of rare variants in CD, we performed whole exome sequencing (WES) in 14 affected members of six families and 12 sporadic cases with CD. We detected rare sequence variants in CFH and FBLN5, which previously were shown to harbor rare variants in patients with CD. In addition, we detected heterozygous rare sequence variants in several genes encoding components of the extracellular matrix (ECM), including FBLN1, FBLN3/EFEMP1, FBLN5, FBLN6/HMCN1, FBN2, and COL15A1. Two rare pathogenic variants were identified in the COL15A1 gene: one in a sporadic case and another was found to segregate in a family with six affected individuals with CD. In addition, two rare pathogenic variants were identified in the FGL1 gene in three unrelated CD cases. These findings suggest that alterations in the ECM and in the coagulation pathway may play a role in the pathogenesis of CD. The identified candidate genes require further analyses in larger cohorts to confirm their role in the CD subtype of AMD. No evidence was found of rare sequence variants in a single gene that segregate with CD in the six families, suggesting that the disease is genetically heterogeneous.     

Network- and attribute-based classifiers can prioritize genes and pathways for autism spectrum disorders and intellectual disability

Autism spectrum disorders (ASD) are a group of related neurodevelopmental disorders with significant combined prevalence (～1%) and high heritability. Dozens of individually rare genes and loci associated with high-risk for ASD have been identified, which overlap extensively with genes for intellectual disability (ID). However, studies indicate that there may be hundreds of genes that remain to be identified. The advent of inexpensive massively parallel nucleotide sequencing can reveal the genetic underpinnings of heritable complex diseases, including ASD and ID. However, whole exome sequencing (WES) and whole genome sequencing (WGS) provides an embarrassment of riches, where many candidate variants emerge. It has been argued that genetic variation for ASD and ID will cluster in genes involved in distinct pathways and protein complexes. For this reason, computational methods that prioritize candidate genes based on additional functional information such as protein-protein interactions or association with specific canonical or empirical pathways, or other attributes, can be useful. In this study we applied several supervised learning approaches to prioritize ASD or ID disease gene candidates based on curated lists of known ASD and ID disease genes. We implemented two network-based classifiers and one attribute-based classifier to show that we can rank and classify known, and predict new, genes for these neurodevelopmental disorders. We also show that ID and ASD share common pathways that perturb an overlapping synaptic regulatory subnetwork. We also show that features relating to neuronal phenotypes in mouse knockouts can help in classifying neurodevelopmental genes. Our methods can be applied broadly to other diseases helping in prioritizing newly identified genetic variation that emerge from disease gene discovery based on WES and WGS.     

Whole Exome Sequencing of Distant Relatives in Multiplex Families Implicates Rare Variants in Candidate Genes for Oral Clefts

A dozen genes/regions have been confirmed as genetic risk factors for oral clefts in human association and linkage studies, and animal models argue even more genes may be involved. Genomic sequencing studies should identify specific causal variants and may reveal additional genes as influencing risk to oral clefts, which have a complex and heterogeneous etiology. We conducted a whole exome sequencing (WES) study to search for potentially causal variants using affected relatives drawn from multiplex cleft families. Two or three affected second, third, and higher degree relatives from 55 multiplex families were sequenced. We examined rare single nucleotide variants (SNVs) shared by affected relatives in 348 recognized candidate genes. Exact probabilities that affected relatives would share these rare variants were calculated, given pedigree structures, and corrected for the number of variants tested. Five novel and potentially damaging SNVs shared by affected distant relatives were found and confirmed by Sanger sequencing. One damaging SNV in CDH1, shared by three affected second cousins from a single family, attained statistical significance (P = 0.02 after correcting for multiple tests). Family-based designs such as the one used in this WES study offer important advantages for identifying genes likely to be causing complex and heterogeneous disorders.     

Multi‐genic pattern found in rare type of hypopituitarism: a whole‐exome sequencing study of Han Chinese with pituitary stalk interruption syndrome

Pituitary stalk interruption syndrome (PSIS) is a rare type of hypopituitarism manifesting various degrees of pituitary hormone deficiency. Although mutations have been identified in some familial cases, the underpinning mechanisms of sporadic patients with PSIS who are in a vast majority remain elusive, necessitating a comprehensive study using systemic approaches. We postulate that other genetic mechanisms may be responsible for the sporadic PSIS. To test this hypothesis, we conducted a study in 24 patients with PSIS of Han Chinese with no family history using whole‐exome sequencing (WES) and bioinformatic analysis. We identified a group of heterozygous mutations in 92% (22 of 24) of the patients, and these genes are mostly associated with Notch, Shh, Wnt signalling pathways. Importantly, 83% (20 of 24) of the patients had more than one mutation in those pathways suggesting synergy of compound mutations underpin the pathogenesis of sporadic PSIS.     

Lysoplex: An efficient toolkit to detect DNA sequence variations in the autophagy-lysosomal pathway

The autophagy-lysosomal pathway (ALP) regulates cell homeostasis and plays a crucial role in human diseases, such as lysosomal storage disorders (LSDs) and common neurodegenerative diseases. Therefore, the identification of DNA sequence variations in genes involved in this pathway and their association with human diseases would have a significant impact on health. To this aim, we developed Lysoplex, a targeted next-generation sequencing (NGS) approach, which allowed us to obtain a uniform and accurate coding sequence coverage of a comprehensive set of 891 genes involved in lysosomal, endocytic, and autophagic pathways. Lysoplex was successfully validated on 14 different types of LSDs and then used to analyze 48 mutation-unknown patients with a clinical phenotype of neuronal ceroid lipofuscinosis (NCL), a genetically heterogeneous subtype of LSD. Lysoplex allowed us to identify pathogenic mutations in 67% of patients, most of whom had been unsuccessfully analyzed by several sequencing approaches. In addition, in 3 patients, we found potential disease-causing variants in novel NCL candidate genes. We then compared the variant detection power of Lysoplex with data derived from public whole exome sequencing (WES) efforts. On average, a 50% higher number of validated amino acid changes and truncating variations per gene were identified. Overall, we identified 61 truncating sequence variations and 488 missense variations with a high probability to cause loss of function in a total of 316 genes. Interestingly, some loss-of-function variations of genes involved in the ALP pathway were found in homozygosity in the normal population, suggesting that their role is not essential. Thus, Lysoplex provided a comprehensive catalog of sequence variants in ALP genes and allows the assessment of their relevance in cell biology as well as their contribution to human disease.     

Family‐based genome scan for age at onset of late‐onset Alzheimer's disease in whole exome sequencing data

Alzheimer's disease (AD) is a common and complex neurodegenerative disease. Age at onset (AAO) of AD is an important component phenotype with a genetic basis, and identification of genes in which variation affects AAO would contribute to identification of factors that affect timing of onset. Increase in AAO through prevention or therapeutic measures would have enormous benefits by delaying AD and its associated morbidities. In this paper, we performed a family‐based genome‐wide association study for AAO of late‐onset AD in whole exome sequence data generated in multigenerational families with multiple AD cases. We conducted single marker and gene‐based burden tests for common and rare variants, respectively. We combined association analyses with variance component linkage analysis, and with reference to prior studies, in order to enhance evidence of the identified genes. For variants and genes implicated by the association study, we performed a gene‐set enrichment analysis to identify potential novel pathways associated with AAO of AD. We found statistically significant association with AAO for three genes (WRN, NTN4 and LAMC3) with common associated variants, and for four genes (SLC8A3, SLC19A3, MADD and LRRK2) with multiple rare‐associated variants that have a plausible biological function related to AD. The genes we have identified are in pathways that are strong candidates for involvement in the development of AD pathology and may lead to a better understanding of AD pathogenesis.     

Differential Distribution of Elastic Tissue in Human Natural Skin and Tissue-Engineered Skin

Tissue-engineered skins (TES), manufactured by epidermal and dermal equivalents, are now being used in biological, pharmacotoxicological and clinical applications. It is thus interesting to know to what extent artificial organs are similar to natural counterparts. Elastic fibres are important constituents of the extracellular matrix of natural skin (NS). The aim of our study was to investigate the possible occurrence and distribution of elastic tissue in a model of human TES using different histochemical techniques, including classical Orcein and Fuchsin-Resorcin methods and immunohistochemistry, at both light and electron microscopical levels. Immunoperoxidase and high resolution immunogold methods were used. In NS, classical staining techniques and elastin-immunohistochemistry revealed a well-organized network of elastic fibres. High resolution immunocytochemistry revealed an intense labelling in the amorphous component of elastic fibres. Fibres of different diameters were immunostained. In TES, no stained elastic fibres were observed using classical staining techniques, and the interpretation of immunoperoxidase observations was not clear-cut. In contrast, immunogold staining at the electron microscopical level provided specific labelling of elastin-like immunoreactive material in the dermal equivalent. However, ultrastructural immunocytochemistry revealed that elastic tissue organization in TES was poor compared to that in NS. This study demonstrates that elastic fibres are a component of the extracellular matrix in this model of TES and suggests that fibroblasts of the dermal equivalent are engaged in matrix secretion. Nevertheless, the level of extracellular matrix organization in TES is low compared to NS. Moreover, this study also suggests that different models of bilayered TES may differ with respect to extracellular matrix organization. These aspects should be considered when TES is used in biological and pharmacotoxicological studies. A better understanding of the factors influencing extracellular matrix formation in TES is necessary to achieve further development of skin generation in vitro.     

Large-Scaled Metabolic Profiling of Human Dermal Fibroblasts Derived from Pseudoxanthoma Elasticum Patients and Healthy Controls

Mutations in the ABC transporter ABCC6 were recently identified as cause of Pseudoxanthoma elasticum (PXE), a rare genetic disorder characterized by progressive mineralization of elastic fibers. We used an untargeted metabolic approach to identify biochemical differences between human dermal fibroblasts from healthy controls and PXE patients in an attempt to find a link between ABCC6 deficiency, cellular metabolic alterations and disease pathogenesis. 358 compounds were identified by mass spectrometry covering lipids, amino acids, peptides, carbohydrates, nucleotides, vitamins and cofactors, xenobiotics and energy metabolites. We found substantial differences in glycerophospholipid composition, leucine dipeptides, and polypeptides as well as alterations in pantothenate and guanine metabolism to be significantly associated with PXE pathogenesis. These findings can be linked to extracellular matrix remodeling and increased oxidative stress, which reflect characteristic hallmarks of PXE. Our study could facilitate a better understanding of biochemical pathways involved in soft tissue mineralization.     

Pathways-Driven Sparse Regression Identifies Pathways and Genes Associated with High-Density Lipoprotein Cholesterol in Two Asian Cohorts

Standard approaches to data analysis in genome-wide association studies (GWAS) ignore any potential functional relationships between gene variants. In contrast gene pathways analysis uses prior information on functional structure within the genome to identify pathways associated with a trait of interest. In a second step, important single nucleotide polymorphisms (SNPs) or genes may be identified within associated pathways. The pathways approach is motivated by the fact that genes do not act alone, but instead have effects that are likely to be mediated through their interaction in gene pathways. Where this is the case, pathways approaches may reveal aspects of a trait''s genetic architecture that would otherwise be missed when considering SNPs in isolation. Most pathways methods begin by testing SNPs one at a time, and so fail to capitalise on the potential advantages inherent in a multi-SNP, joint modelling approach. Here, we describe a dual-level, sparse regression model for the simultaneous identification of pathways and genes associated with a quantitative trait. Our method takes account of various factors specific to the joint modelling of pathways with genome-wide data, including widespread correlation between genetic predictors, and the fact that variants may overlap multiple pathways. We use a resampling strategy that exploits finite sample variability to provide robust rankings for pathways and genes. We test our method through simulation, and use it to perform pathways-driven gene selection in a search for pathways and genes associated with variation in serum high-density lipoprotein cholesterol levels in two separate GWAS cohorts of Asian adults. By comparing results from both cohorts we identify a number of candidate pathways including those associated with cardiomyopathy, and T cell receptor and PPAR signalling. Highlighted genes include those associated with the L-type calcium channel, adenylate cyclase, integrin, laminin, MAPK signalling and immune function.     

Integrated Model of De Novo and Inherited Genetic Variants Yields Greater Power to Identify Risk Genes

De novo mutations affect risk for many diseases and disorders, especially those with early-onset. An example is autism spectrum disorders (ASD). Four recent whole-exome sequencing (WES) studies of ASD families revealed a handful of novel risk genes, based on independent de novo loss-of-function (LoF) mutations falling in the same gene, and found that de novo LoF mutations occurred at a twofold higher rate than expected by chance. However successful these studies were, they used only a small fraction of the data, excluding other types of de novo mutations and inherited rare variants. Moreover, such analyses cannot readily incorporate data from case-control studies. An important research challenge in gene discovery, therefore, is to develop statistical methods that accommodate a broader class of rare variation. We develop methods that can incorporate WES data regarding de novo mutations, inherited variants present, and variants identified within cases and controls. TADA, for Transmission And De novo Association, integrates these data by a gene-based likelihood model involving parameters for allele frequencies and gene-specific penetrances. Inference is based on a Hierarchical Bayes strategy that borrows information across all genes to infer parameters that would be difficult to estimate for individual genes. In addition to theoretical development we validated TADA using realistic simulations mimicking rare, large-effect mutations affecting risk for ASD and show it has dramatically better power than other common methods of analysis. Thus TADA''s integration of various kinds of WES data can be a highly effective means of identifying novel risk genes. Indeed, application of TADA to WES data from subjects with ASD and their families, as well as from a study of ASD subjects and controls, revealed several novel and promising ASD candidate genes with strong statistical support.     

A heart for fibrillin: spatial arrangement in adult wild-type murine myocardial tissue

Fibrillins are major constituents of microfibrils, which are essential components of the extracellular matrix of connective tissues where they contribute to the tissue homeostasis. Although it is known that microfibrils are abundantly expressed in the left ventricle of the heart, limited data are available about the presence of microfibrils in the other parts of the myocardial tissue and whether there are age or sex-related differences in the spatial arrangement of the microfibrils. This basic knowledge is essential to better understand the impact of fibrillin-1 pathogenic variants on the myocardial tissue as seen in Marfan related cardiomyopathy. We performed histological analyses on wild-type male and female murine myocardial tissue collected at different time-points (1, 3 and 6 months). Fibrillin-1 and -2 immunofluorescence stainings were performed on cross-sections at the level of the apex, the mid-ventricles and the atria. In addition, other myocardial matrix components such as collagen and elastin were also investigated. Fibrillin-1 presented as long fibres in the apex, mid-ventricles and atria. The spatial arrangement differed between the investigated regions, but not between age groups or sexes. Collagen had a similar broad spatial arrangement to that of fibrillin-1, whereas elastic fibres were primarily present in the atria and the vessels. In contrast to fibrillin-1, limited amounts of fibrillin-2 were observed. Fibrillin-rich fibres contribute to the architecture of the myocardial tissue in a region-dependent manner in wild-type murine hearts. This knowledge is helpful for future experimental set-ups of studies evaluating the impact of fibrillin-1 pathogenic variants on the myocardial tissue.     

Matrix Gla protein accumulates at the border of regions of calcification and normal tissue in the media of the arterial vessel wall.

Vitamin K-dependent matrix Gla protein (MGP) has been suggested to play a role in the inhibition of soft-tissue calcification. Here we report the expression of recombinant prokaryotic MGP as part of a fusion protein and the preparation of two antibodies that specifically recognize MGP. Monoclonal antibodies were raised against synthetic peptides homologous to the sequences 3-15 and 63-75 of human MGP. Both antibodies recognize recombinant and synthetic human MGP. Immunohistochemical analysis showed that MGP was associated with the extracellular matrix of noncalcified bone and with chondrocytes in cartilage. In the healthy human arterial vessel wall, MGP antigen was demonstrated in association with smooth muscle cells and elastic laminae of the tunica media and with the extracellular matrix of the adventitia. Colocalization with the elastic laminae was lost at sites of medial calcification; in both human and rat arteries, high amounts of MGP were found in the extracellular matrix at borders of intimal and medial calcification. Our data demonstrate the close association between MGP and calcification. It is suggested that undercarboxylated MGP is biologically inactive and that poor vascular vitamin K status may form a risk factor for vascular calcification.     

Wnt and RUNX2 mediate cartilage breakdown by osteoarthritis synovial fibroblast‐derived ADAMTS‐7 and ‐12

Failure of therapeutic approaches for the treatment of osteoarthritis (OA) based on the inhibition of metalloproteinases, might be because of their constitutive expression in homeostasis, together with their network complexity. The knowledge of this network would contribute to selective target pathological conditions. In this sense, blockade of mediators produced by neighbouring joint cells, such as synovial fibroblasts (SF), would prevent cartilage damage. Thus, we studied the contribution of ADAMTS‐7 and ‐12 from SF to cartilage oligomeric matrix protein (COMP) degradation, and the signalling pathways involved in their expression. We report for the first time in SF, the involvement of ERK‐Runx2 axis and Wnt/β‐catenin signalling in ADAMTS‐12 and ADAMTS‐7 expressions, respectively, with the subsequent consequences in COMP degradation from cartilage extracellular matrix. After stimulation with IL‐1β or fibronectin fragments, we showed that ERK inhibition decreased Runx2 activation and ADAMTS‐12 expression in OA‐SF, also reducing Fn‐fs‐induced COMP degradation. Blockage of Wnt signalling by DKK1 reduced ADAMTS‐7 and COMP degradation in OA‐SF as well. In addition, Wnt7B expression was induced by IL‐1β and by itself, also increasing ADAMTS‐7. Our results could contribute to the development of disease‐modifying OA drugs targeting ADAMTS‐7 and ‐12 for the prevention of extracellular matrix components degradation like COMP.     

Coordinated activation of metabolic pathways for antioxidants and defence compounds by jasmonates and their roles in stress tolerance in Arabidopsis

Jasmonic acid (JA) and methyl jasmonate (MeJA), collectively termed jasmonates, are ubiquitous plant signalling compounds. Several types of stress conditions, such as wounding and pathogen infection, cause endogenous JA accumulation and the expression of jasmonate-responsive genes. Although jasmonates are important signalling components for the stress response in plants, the mechanism by which jasmonate signalling contributes to stress tolerance has not been clearly defined. A comprehensive analysis of jasmonate-regulated metabolic pathways in Arabidopsis was performed using cDNA macroarrays containing 13516 expressed sequence tags (ESTs) covering 8384 loci. The results showed that jasmonates activate the coordinated gene expression of factors involved in nine metabolic pathways belonging to two functionally related groups: (i) ascorbate and glutathione metabolic pathways, which are important in defence responses to oxidative stress, and (ii) biosynthesis of indole glucosinolate, which is a defence compound occurring in the Brassicaceae family. We confirmed that JA induces the accumulation of ascorbate, glutathione and cysteine and increases the activity of dehydroascorbate reductase, an enzyme in the ascorbate recycling pathway. These antioxidant metabolic pathways are known to be activated under oxidative stress conditions. Ozone (O3) exposure, a representative oxidative stress, is known to cause activation of antioxidant metabolism. We showed that O3 exposure caused the induction of several genes involved in antioxidant metabolism in the wild type. However, in jasmonate-deficient Arabidopsis 12-oxophytodienoate reductase 3 (opr3) mutants, the induction of antioxidant genes was abolished. Compared with the wild type, opr3 mutants were more sensitive to O3 exposure. These results suggest that the coordinated activation of the metabolic pathways mediated by jasmonates provides resistance to environmental stresses.