IQCJ-SCHIP1, a novel fusion transcript encoding a calmodulin-binding IQ motif protein

Fast and flexible profiling of chromatin accessibility and total RNA expression in single nuclei using Microwell-seq3

Single cell chromatin accessibility profiling and transcriptome sequencing are the most widely used technologies for single-cell genomics. Here, we present Microwell-seq3, a high-throughput and facile platform for high-sensitivity single-nucleus chromatin accessibility or full-length transcriptome profiling. The method combines a preindexing strategy and a penetrable chip-in-a-tube for single nucleus loading and DNA amplification and therefore does not require specialized equipment. We used Microwell-seq3 to profile chromatin accessibility in more than 200,000 single nuclei and the full-length transcriptome in ~50,000 nuclei from multiple adult mouse tissues. Compared with the existing polyadenylated transcript capture methods, integrative analysis of cell type-specific regulatory elements and total RNA expression uncovered comprehensive cell type heterogeneity in the brain. Gene regulatory networks based on chromatin accessibility profiling provided an improved cell type communication model. Finally, we demonstrated that Microwell-seq3 can identify malignant cells and their specific regulons in spontaneous lung tumors of aged mice. We envision a broad application of Microwell-seq3 in many areas of research.

Sex-specific genetic architecture of late-life memory performance

BACKGROUND

Women demonstrate a memory advantage when cognitively healthy yet lose this advantage to men in Alzheimer's disease. However, the genetic underpinnings of this sex difference in memory performance remain unclear.

METHODS

We conducted the largest sex-aware genetic study on late-life memory to date (Nmales  = 11,942; Nfemales  = 15,641). Leveraging harmonized memory composite scores from four cohorts of cognitive aging and AD, we performed sex-stratified and sex-interaction genome-wide association studies in 24,216 non-Hispanic White and 3367 non-Hispanic Black participants.

RESULTS

We identified three sex-specific loci (rs67099044-CBLN2, rs719070-SCHIP1/IQCJ-SCHIP), including an X-chromosome locus (rs5935633-EGL6/TCEANC/OFD1), that associated with memory. Additionally, we identified heparan sulfate signaling as a sex-specific pathway and found sex-specific genetic correlations between memory and cardiovascular, immune, and education traits.

DISCUSSION

This study showed memory is highly and comparably heritable across sexes, as well as highlighted novel sex-specific genes, pathways, and genetic correlations that related to late-life memory.

HIGHLIGHTS

Demonstrated the heritable component of late-life memory is similar across sexes. Identified two genetic loci with a sex-interaction with baseline memory. Identified an X-chromosome locus associated with memory decline in females. Highlighted sex-specific candidate genes and pathways associated with memory. Revealed sex-specific shared genetic architecture between memory and complex traits.

Genomics and transcriptomics landscapes associated to changes in insulin sensitivity in response to endurance exercise training

Despite good adherence to supervised endurance exercise training (EET), some individuals experience no or little improvement in peripheral insulin sensitivity. The genetic and molecular mechanisms underlying this phenomenon are currently not understood. By investigating genome-wide variants associated with baseline and exercise-induced changes (∆) in insulin sensitivity index (S_i) in healthy volunteers, we have identified novel candidate genes whose mouse knockouts phenotypes were consistent with a causative effect on S_i. An integrative analysis of functional genomic and transcriptomic profiles suggests genetic variants have an aggregate effect on baseline S_i and ∆S_i, focused around cholinergic signalling, including downstream calcium and chemokine signalling. The identification of calcium regulated MEF2A transcription factor as the most statistically significant candidate driving the transcriptional signature associated to ∆S_i further strengthens the relevance of calcium signalling in EET mediated S_i response.

The plasma peptides of breast versus ovarian cancer

Background

There is a need to demonstrate a proof of principle that proteomics has the capacity to analyze plasma from breast cancer versus other diseases and controls in a multisite clinical trial design. The peptides or proteins that show a high observation frequency, and/or precursor intensity, specific to breast cancer plasma might be discovered by comparison to other diseases and matched controls. The endogenous tryptic peptides of breast cancer plasma were compared to ovarian cancer, female normal, sepsis, heart attack, Alzheimer's and multiple sclerosis along with the institution-matched normal and control samples collected directly onto ice.

Methods

Endogenous tryptic peptides were extracted from individual breast cancer and control EDTA plasma samples in a step gradient of acetonitrile, and collected over preparative C18 for LC-ESI-MS/MS with a set of LTQ XL linear quadrupole ion traps working together in parallel to randomly and independently sample clinical populations. The MS/MS spectra were fit to fully tryptic peptides or phosphopeptides within proteins using the X!TANDEM algorithm. The protein observation frequency was counted using the SEQUEST algorithm after selecting the single best charge state and peptide sequence for each MS/MS spectra. The observation frequency was subsequently tested by Chi Square analysis. The log₁₀ precursor intensity was compared by ANOVA in the R statistical system.

Results

Peptides and/or phosphopeptides of common plasma proteins such as APOE, C4A, C4B, C3, APOA1, APOC2, APOC4, ITIH3 and ITIH4 showed increased observation frequency and/or precursor intensity in breast cancer. Many cellular proteins also showed large changes in frequency by Chi Square (χ² > 100, p < 0.0001) in the breast cancer samples such as CPEB1, LTBP4, HIF-1A, IGHE, RAB44, NEFM, C19orf82, SLC35B1, 1D12A, C8orf34, HIF1A, OCLN, EYA1, HLA-DRB1, LARS, PTPDC1, WWC1, ZNF562, PTMA, MGAT1, NDUFA1, NOGOC, OR1E1, OR1E2, CFI, HSA12, GCSH, ELTD1, TBX15, NR2C2, FLJ00045, PDLIM1, GALNT9, ASH2L, PPFIBP1, LRRC4B, SLCO3A1, BHMT2, CS, FAM188B2, LGALS7, SAT2, SFRS8, SLC22A12, WNT9B, SLC2A4, ZNF101, WT1, CCDC47, ERLIN1, SPFH1, EID2, THOC1, DDX47, MREG, PTPRE, EMILIN1, DKFZp779G1236 and MAP3K8 among others. The protein gene symbols with large Chi Square values were significantly enriched in proteins that showed a complex set of previously established functional and structural relationships by STRING analysis. An increase in mean precursor intensity of peptides was observed for QSER1 as well as SLC35B1, IQCJ-SCHIP1, MREG, BHMT2, LGALS7, THOC1, ANXA4, DHDDS, SAT2, PTMA and FYCO1 among others. In contrast, the QSER1 peptide QPKVKAEPPPK was apparently specific to ovarian cancer.

Conclusion

There was striking agreement between the breast cancer plasma peptides and proteins discovered by LC-ESI-MS/MS with previous biomarkers from tumors, cells lines or body fluids by genetic or biochemical methods. The results indicate that variation in plasma peptides from breast cancer versus ovarian cancer may be directly discovered by LC-ESI-MS/MS that will be a powerful tool for clinical research. It may be possible to use a battery of sensitive and robust linear quadrupole ion traps for random and independent sampling of plasma from a multisite clinical trial.

Fine mapping of genome-wide association study signals to identify genetic markers of the plasma triglyceride response to an omega-3 fatty acid supplementation

Background

Using a genome-wide association study (GWAS) approach, our group previously computed a genetic risk score (GRS) from single nucleotide polymorphisms (SNPs) of 10 loci that affect the plasma triglyceride (TG) response to an omega-3 (n-3) fatty acid (FA) supplementation.

Objectives

The objective was to compute a novel and more refined GRS using fine mapping to include a large number of genetic variants.

Methods

A total of 208 participants of the Fatty Acid Sensor (FAS) Study received 5 g fish oil/d, containing 1.9-2.2 g eicosapentaenoic acid and 1.1 g docosahexanoic acid, for 6 wk. Plasma TG concentrations were measured before and after supplementation. Dense genotyping and genotype imputation were used to refine mapping around GWAS hits. A GRS was computed by summing the number of at-risk alleles of tagging SNPs. Analyses were replicated in samples of the FINGEN study.

Results

A total of 31 tagging SNPs associated with the TG response were used for GRS calculation in the FAS study. In a general linear model adjusted for age, sex, and body mass index, the GRS explained 49.73% of TG response variance (P < 0.0001). Nonresponders to the n-3 FA supplementation had a higher GRS than did responders. In the FINGEN replication study, the GRS explained 3.67% of TG response variance (P = 0.0006).

Conclusions

Fine mapping proved to be effective to refine the previous GRS. Carrying increasing numbers of at-risk alleles of 31 SNPs confers a higher risk of being nonresponsive to n-3 FAs. The genetic profile therefore appears to be an important determinant of the plasma TG response to an n-3 FA supplementation and could be used to target those most likely to gain clinical benefit. This trial was registered at http://www.clinicaltrials.gov as NCT01343342.

Homozygous nonsense mutation in SCHIP1/IQCJ-SCHIP1 causes a neurodevelopmental brain malformation syndrome

We report a consanguineous Arab family with 3 affected siblings who display a disorder of global developmental delay, learning difficulties, facial dysmorphism, hearing impairments, and cataract. The clinical phenotype was associated with characteristic brain magnetic resonance imaging (MRI) features of axonal guidance defects involving anterior commissure agenesis as well as scattered areas of polymicrogyria-cobblestone complex. Whole genome sequencing revealed a novel nonsense mutation (159609921C>T) that segregated in the family consistent in an autosomal recessive pattern. This mutation located in the C-terminal region shared by the Schwanomin-Interacting Protein1 (SCHIP1) isoforms including the IQCJ-SCHIP1. The in vitro expression of SCHIP1 and IQCJ-SCHIP1 truncated mutant isoforms (NM_001197109.1; p.R209* and NM_001197114.1; p.R501*, respectively) were markedly reduced as compared to their full-length versions suggesting protein stability/folding impairment. The pathogenic nature of this mutation is supported by a previously reported mouse knockout of Schip1 isoforms, which phenocopied the human axon guidance abnormality. This is the first report of a SCHIP1/IQCJ-SCHIP1 point mutation in humans associated with a neurological-developmental phenotype.

Plasma Triglyceride Levels May Be Modulated by Gene Expression of IQCJ, NXPH1, PHF17 and MYB in Humans

A genome-wide association study (GWAS) by our group identified loci associated with the plasma triglyceride (TG) response to ω-3 fatty acid (FA) supplementation in IQCJ, NXPH1, PHF17 and MYB. Our aim is to investigate potential mechanisms underlying the associations between single nucleotide polymorphisms (SNPs) in the four genes and TG levels following ω-3 FA supplementation. 208 subjects received 3 g/day of ω-3 FA (1.9–2.2 g of EPA and 1.1 g of docosahexaenoic acid (DHA)) for six weeks. Plasma TG were measured before and after the intervention. 67 SNPs were selected to increase the density of markers near GWAS hits. Genome-wide expression and methylation analyses were conducted on respectively 30 and 35 participants’ blood sample together with in silico analyses. Two SNPs of IQCJ showed different affinities to splice sites depending on alleles. Expression levels were influenced by genotype for one SNP in NXPH1 and one in MYB. Associations between 12 tagged SNPs of IQCJ, 26 of NXPH1, seven of PHF17 and four of MYB and gene-specific CpG site methylation levels were found. The response of plasma TG to ω-3 FA supplementation may be modulated by the effect of DNA methylation on expression levels of genes revealed by GWAS.

Schwannomin-interacting Protein 1 Isoform IQCJ-SCHIP1 Is a Multipartner Ankyrin- and Spectrin-binding Protein Involved in the Organization of Nodes of Ranvier

The nodes of Ranvier are essential regions for action potential conduction in myelinated fibers. They are enriched in multimolecular complexes composed of voltage-gated Nav and Kv7 channels associated with cell adhesion molecules. Cytoskeletal proteins ankyrin-G (AnkG) and βIV-spectrin control the organization of these complexes and provide mechanical support to the plasma membrane. IQCJ-SCHIP1 is a cytoplasmic protein present in axon initial segments and nodes of Ranvier. It interacts with AnkG and is absent from nodes and axon initial segments of βIV-spectrin and AnkG mutant mice. Here, we show that IQCJ-SCHIP1 also interacts with βIV-spectrin and Kv7.2/3 channels and self-associates, suggesting a scaffolding role in organizing nodal proteins. IQCJ-SCHIP1 binding requires a βIV-spectrin-specific domain and Kv7 channel 1-5-10 calmodulin-binding motifs. We then investigate the role of IQCJ-SCHIP1 in vivo by studying peripheral myelinated fibers in Schip1 knock-out mutant mice. The major nodal proteins are normally enriched at nodes in these mice, indicating that IQCJ-SCHIP1 is not required for their nodal accumulation. However, morphometric and ultrastructural analyses show an altered shape of nodes similar to that observed in βIV-spectrin mutant mice, revealing that IQCJ-SCHIP1 contributes to nodal membrane-associated cytoskeleton organization, likely through its interactions with the AnkG/βIV-spectrin network. Our work reveals that IQCJ-SCHIP1 interacts with several major nodal proteins, and we suggest that it contributes to a higher organizational level of the AnkG/βIV-spectrin network critical for node integrity.

Novel Genetic Loci Associated with the Plasma Triglyceride Response to an Omega-3 Fatty Acid Supplementation

BACKGROUND

A recent genome-wide association study (GWAS) by our group identified 13 loci associated with the plasma triglyceride (TG) response to omega-3 (n-3) fatty acid (FA) supplementation. This study aimed to test whether single-nucleotide polymorphisms (SNPs) within the IQCJ, NXPH1, PHF17 and MYB genes are associated with the plasma TG response to an n-3 FA supplementation.

METHODS

A total of 208 subjects followed a 6-week n-3 FA supplementation of 5 g/day of fish oil (1.9-2.2 g of eicosapentaenoic acid and 1.1 g of docosahexaenoic acid). Measurements of plasma lipids were made before and after the supplementation. Sixty-seven tagged SNPs were selected to increase the density of markers near GWAS hits.

RESULTS

In a repeated model, independent effects of the genotype and the gene-supplementation interaction were associated with plasma TG. Genotype effects were observed with two SNPs of NXPH1, and gene-diet interactions were observed with ten SNPs of IQCJ, four SNPs of NXPH1 and three SNPs of MYB. Positive and negative responders showed different genotype frequencies with nine SNPs of IQCJ, two SNPs of NXPH1 and two SNPs of MYB.

CONCLUSION

Fine mapping in GWAS-associated loci allowed the identification of SNPs partly explaining the large interindividual variability observed in plasma TG levels in response to an n-3 FA supplementation.

CK2-regulated schwannomin-interacting protein IQCJ-SCHIP-1 association with AnkG contributes to the maintenance of the axon initial segment

The axon initial segment (AIS) plays a central role in electrogenesis and in the maintenance of neuronal polarity. Its molecular organization is dependent on the scaffolding protein ankyrin (Ank) G and is regulated by kinases. For example, the phosphorylation of voltage-gated sodium channels by the protein kinase CK2 regulates their interaction with AnkG and, consequently, their accumulation at the AIS. We previously showed that IQ motif containing J-Schwannomin-Interacting Protein 1 (IQCJ-SCHIP-1), an isoform of the SCHIP-1, accumulated at the AIS in vivo. Here, we analyzed the molecular mechanisms involved in IQCJ-SCHIP-1-specific axonal location. We showed that IQCJ-SCHIP-1 accumulation in the AIS of cultured hippocampal neurons depended on AnkG expression. Pull-down assays and surface plasmon resonance analysis demonstrated that AnkG binds to CK2-phosphorylated IQCJ-SCHIP-1 but not to the non-phosphorylated protein. Surface plasmon resonance approaches using IQCJ-SCHIP-1, SCHIP-1a, another SCHIP-1 isoform, and their C-terminus tail mutants revealed that a segment including multiple CK2-phosphorylatable sites was directly involved in the interaction with AnkG. Pharmacological inhibition of CK2 diminished both IQCJ-SCHIP-1 and AnkG accumulation in the AIS. Silencing SCHIP-1 expression reduced AnkG cluster at the AIS. Finally, over-expression of IQCJ-SCHIP-1 decreased AnkG concentration at the AIS, whereas a mutant deleted of the CK2-regulated AnkG interaction site did not. Our study reveals that CK2-regulated IQJC-SCHIP-1 association with AnkG contributes to AIS maintenance. The axon initial segment (AIS) organization depends on ankyrin (Ank) G and kinases. Here we showed that AnkG binds to CK2-phosphorylated IQCJ-SCHIP-1, in a segment including 12 CK2-phosphorylatable sites. In cultured neurons, either pharmacological inhibition of CK2 or IQCJ-SCHIP-1 silencing reduced AnkG clustering. Overexpressed IQCJ-SCHIP-1 decreased AnkG concentration at the AIS whereas a mutant deleted of the CK2-regulated AnkG interaction site did not. Thus, CK2-regulated IQJC-SCHIP-1 association with AnkG contributes to AIS maintenance.

Deep RNA sequencing analysis of readthrough gene fusions in human prostate adenocarcinoma and reference samples

BACKGROUND

Readthrough fusions across adjacent genes in the genome, or transcription-induced chimeras (TICs), have been estimated using expressed sequence tag (EST) libraries to involve 4-6% of all genes. Deep transcriptional sequencing (RNA-Seq) now makes it possible to study the occurrence and expression levels of TICs in individual samples across the genome.

METHODS

We performed single-end RNA-Seq on three human prostate adenocarcinoma samples and their corresponding normal tissues, as well as brain and universal reference samples. We developed two bioinformatics methods to specifically identify TIC events: a targeted alignment method using artificial exon-exon junctions within 200,000 bp from adjacent genes, and genomic alignment allowing splicing within individual reads. We performed further experimental verification and characterization of selected TIC and fusion events using quantitative RT-PCR and comparative genomic hybridization microarrays.

RESULTS

Targeted alignment against artificial exon-exon junctions yielded 339 distinct TIC events, including 32 gene pairs with multiple isoforms. The false discovery rate was estimated to be 1.5%. Spliced alignment to the genome was less sensitive, finding only 18% of those found by targeted alignment in 33-nt reads and 59% of those in 50-nt reads. However, spliced alignment revealed 30 cases of TICs with intervening exons, in addition to distant inversions, scrambled genes, and translocations. Our findings increase the catalog of observed TIC gene pairs by 66%.We verified 6 of 6 predicted TICs in all prostate samples, and 2 of 5 predicted novel distant gene fusions, both private events among 54 prostate tumor samples tested. Expression of TICs correlates with that of the upstream gene, which can explain the prostate-specific pattern of some TIC events and the restriction of the SLC45A3-ELK4 e4-e2 TIC to ERG-negative prostate samples, as confirmed in 20 matched prostate tumor and normal samples and 9 lung cancer cell lines.

CONCLUSIONS

Deep transcriptional sequencing and analysis with targeted and spliced alignment methods can effectively identify TIC events across the genome in individual tissues. Prostate and reference samples exhibit a wide range of TIC events, involving more genes than estimated previously using ESTs. Tissue specificity of TIC events is correlated with expression patterns of the upstream gene. Some TIC events, such as MSMB-NCOA4, may play functional roles in cancer.

[Axon and Schwann cells... so far away, so close]

Myelination was a major step in the evolution of the nervous system. Appearing first in jaw fish, myelination allows the fast and secure propagation of action potentials at a low energetic cost, and without exaggerated increase in axonal diameter. In the peripheral nervous system of mammals, myelination results from the tight interactions between Schwann cells and axons, leading to the formation of highly differentiated domains along the axon. The molecular determinants of these interactions are starting to be well identified. Their understanding provides a precise framework to interpret the defects, which occur in pathological circumstances. This review summarizes the present state of knowledge concerning axoglial interactions in peripheral nerves.

Schwannomin-interacting protein-1 isoform IQCJ-SCHIP-1 is a late component of nodes of Ranvier and axon initial segments

Axon initial segments (AISs) and nodes of Ranvier (NRs) are essential regions for saltatory conduction of the action potential along the axon. These two domains are enriched in similar multimolecular complexes, which include voltage-gated sodium channels (Na(v)), NF186 (neurofascin 186), NrCAM (neuron glia-related cell adhesion molecule), and cytoskeleton linkers ankyrin G (AnkG) and betaIV-spectrin. Identification of novel members of these complexes is critical to better understand their formation, function, and maintenance. Here we report that IQCJ-SCHIP-1, a recently identified isoform of schwannomin-interacting protein-1 (SCHIP-1), is a novel component of both AISs and NRs in the central and peripheral nervous systems. We show that IQCJ-SCHIP-1 binds calmodulin in the absence of Ca(2+) and is highly enriched at AISs and NRs. IQCJ-SCHIP-1 accumulation at AISs and NRs is a late event, suggesting that IQCJ-SCHIP-1 is likely to play a role in mature AISs and NRs rather than during their formation. IQCJ-SCHIP-1 was not detected at AISs in the absence of AnkG and interacted in vitro with this protein. IQCJ-SCHIP-1 was also absent from central NRs and AISs of quivering mice, which have a mutation of betaIV-spectrin. We suggest that IQCJ-SCHIP-1 might participate, along with AnkG and betaIV-spectrin, in the stabilization or function of the multimolecular complexes of AISs and NRs, possibly by participating in Ca(2+)-mediated responses.

The merlin interacting proteins reveal multiple targets for NF2 therapy

The neurofibromatosis 2 (NF2) tumor suppressor protein merlin is commonly mutated in human benign brain tumors. The gene altered in NF2 was located on human chromosome 22q12 in 1993 and the encoded protein named merlin and schwannomin. Merlin has homology to ERM family proteins, ezrin, radixin, and moesin, within the protein 4.1 superfamily. In efforts to determine merlin function several groups have discovered 34 merlin interacting proteins, including ezrin, radixin, moesin, CD44, layilin, paxillin, actin, N-WASP, betaII-spectrin, microtubules, TRBP, eIF3c, PIKE, NHERF, MAP, RalGDS, RhoGDI, EG1/magicin, HEI10, HRS, syntenin, caspr/paranodin, DCC, NGB, CRM1/exportin, SCHIP1, MYPT-1-PP1delta, RIbeta, PKA, PAK (three types), calpain and Drosophila expanded. Many of the proteins that interact with the merlin N-terminal domain also bind ezrin, while other merlin interacting proteins do not bind other members of the ERM family. Merlin also interacts with itself. This review describes these proteins, their possible roles in NF2, and the resultant hypothesized merlin functions. Review of all of the merlin interacting proteins and functional consequences of losses of these interactions reveals multiple merlin actions in PI3-kinase, MAP kinase and small GTPase signaling pathways that might be targeted to inhibit the proliferation of NF2 tumors.