Frequent genes in rare diseases: panel‐based next generation sequencing to disclose causal mutations in hereditary neuropathies

Hereditary neuropathies comprise a wide variety of chronic diseases associated to more than 80 genes identified to date. We herein examined 612 index patients with either a Charcot‐Marie‐Tooth phenotype, hereditary sensory neuropathy, familial amyloid neuropathy, or small fiber neuropathy using a customized multigene panel based on the next generation sequencing technique. In 121 cases (19.8%), we identified at least one putative pathogenic mutation. Of these, 54.4% showed an autosomal dominant, 33.9% an autosomal recessive, and 11.6% an X‐linked inheritance. The most frequently affected genes were PMP22 (16.4%), GJB1 (10.7%), MPZ, and SH3TC2 (both 9.9%), and MFN2 (8.3%). We further detected likely or known pathogenic variants in HINT1, HSPB1, NEFL, PRX, IGHMBP2, NDRG1, TTR, EGR2, FIG4, GDAP1, LMNA, LRSAM1, POLG, TRPV4, AARS, BIC2, DHTKD1, FGD4, HK1, INF2, KIF5A, PDK3, REEP1, SBF1, SBF2, SCN9A, and SPTLC2 with a declining frequency. Thirty‐four novel variants were considered likely pathogenic not having previously been described in association with any disorder in the literature. In one patient, two homozygous mutations in HK1 were detected in the multigene panel, but not by whole exome sequencing. A novel missense mutation in KIF5A was considered pathogenic because of the highly compatible phenotype. In one patient, the plasma sphingolipid profile could functionally prove the pathogenicity of a mutation in SPTLC2. One pathogenic mutation in MPZ was identified after being previously missed by Sanger sequencing. We conclude that panel based next generation sequencing is a useful, time‐ and cost‐effective approach to assist clinicians in identifying the correct diagnosis and enable causative treatment considerations.

     

Anistropically varying conductivity in irreversible electroporation simulations

Background

One recent area of cancer research is irreversible electroporation (IRE). Irreversible electroporation is a minimally invasive procedure where needle electrodes are inserted into the body to ablate tumor cells with electricity. The aim of this paper is to propose a mathematical model that incorporates a tissue’s conductivity increasing more in the direction of the electrical field as this has been shown to occur in experiments.

Method

It was necessary to mathematically derive a valid form of the conductivity tensor such that it is dependent on the electrical field direction and can be easily implemented into numerical software. The derivation of a conductivity tensor that can take arbitrary functions for the conductivity in the directions tangent and normal to the electrical field is the main contribution of this paper. Numerical simulations were performed for isotropic-varying and anisotropic-varying conductivities to evaluate the importance of including the electrical field’s direction in the formulation for conductivity.

Results

By starting from previously published experimental results, this paper derived a general formulation for an anistropic-varying tensor for implementation into irreversible electroporation modeling software. The anistropic-varying tensor formulation allows the conductivity to take into consideration both electrical field direction and magnitude, as opposed to previous published works that only took into account electrical field magnitude.The anisotropic formulation predicts roughly a five percent decrease in ablation size for the monopolar simulation and approximately a ten percent decrease in ablation size for the bipolar simulations. This is a positive result as previously reported results found the isotropic formulation to overpredict ablation size for both monopolar and bipolar simulations. Furthermore, it was also reported that the isotropic formulation overpredicts the ablation size more for the bipolar case than the monopolar case. Thus, our results are following the experimental trend by having a larger percentage change in volume for the bipolar case than the monopolar case.

Conclusions

The predicted volume of ablated cells decreased, and could be a possible explanation for the slight over-prediction seen by isotropic-varying formulations.

     

Comprehensive proteomic analysis of Penicillium verrucosum

Mass spectrometric identification of proteins in species lacking validated sequence information is a major problem in veterinary science. In the present study, we used ochratoxin A producing Penicillium verrucosum to identify and quantitatively analyze proteins of an organism with yet no protein information available. The work presented here aimed to provide a comprehensive protein identification of P. verrucosum using shotgun proteomics. We were able to identify 3631 proteins in an “ab initio” translated database from DNA sequences of P. verrucosum. Additionally, a sequential window acquisition of all theoretical fragment‐ion spectra analysis was done to find differentially regulated proteins at two different time points of the growth curve. We compared the proteins at the beginning (day 3) and at the end of the log phase (day 12).     

Systemic acquired resistance: the elusive signal(s)

Systemic acquired resistance (SAR) is a form of inducible resistance that is triggered in systemic healthy tissues of locally infected plants. The nature of the mobile signal that travels through the phloem from the site of infection to establish systemic immunity has been sought after for decades. Several candidate signaling molecules have emerged in the past two years, including the methylated derivative of a well-known defense hormone (methyl salicylate), the defense hormone jasmonic acid, a yet undefined glycerolipid-derived factor, and a group of peptides that is involved in cell-to-cell basal defense signaling. Systemic SAR signal amplification increasingly appears to parallel salicylic acid-dependent defense responses, and is concomitantly fine-tuned by auxin.     

Biochemical and mutagenic analysis of I-CreII reveals distinct but important roles for both the H-N-H and GIY-YIG motifs

Homing endonucleases typically contain one of four conserved catalytic motifs, and other elements that confer tight DNA binding. I-CreII, which catalyzes homing of the Cr.psbA4 intron, is unusual in containing two potential catalytic motifs, H-N-H and GIY-YIG. Previously, we showed that cleavage by I-CreII leaves ends (2-nt 3′ overhangs) that are characteristic of GIY-YIG endonucleases, yet it has a relaxed metal requirement like H-N-H enzymes. Here we show that I-CreII can bind DNA without an added metal ion, and that it binds as a monomer, akin to GIY-YIG enzymes. Moreover, cleavage of supercoiled DNA, and estimates of strand-specific cleavage rates, suggest that I-CreII uses a sequential cleavage mechanism. Alanine substitution of a number of residues in the GIY-YIG motif, however, did not block cleavage activity, although DNA binding was substantially reduced in several variants. Substitution of conserved histidines in the H-N-H motif resulted in variants that did not promote DNA cleavage, but retained high-affinity DNA binding—thus identifying it as the catalytic motif. Unlike the non-specific H-N-H colicins, however; substitution of the conserved asparagine substantially reduced DNA binding (though not the ability to promote cleavage). These results indicate that, in I-CreII, two catalytic motifs have evolved to play important roles in specific DNA binding. The data also indicate that only the H-N-H motif has retained catalytic ability.     

The effect of thermal treatment on antibacterial properties of nanostructured TiO<Subscript>2</Subscript>(N) films illuminated with visible light

This work focuses on the photocatalytic performances and antibacterial activity of nitrogen doped TiO₂ nanosystems with three and five layers obtained by a sol-gel route, followed by thermal treatment in oxygen or ammonia atmosphere at temperatures between 400 and 1000°C. Subsequently, the antibacterial activity of the obtained nanosystems on the Escherichia coli cells are determined and discussed. The obtained results show a significant dependence of the functional performances on the system’s composition. In particular, the antimicrobial activity of nitrogen-doped TiO₂ films is correlated with the temperature of thermal treatment and illumination time with visible artificial light.     

Mutation-selection equilibrium in games with multiple strategies

In evolutionary games the fitness of individuals is not constant but depends on the relative abundance of the various strategies in the population. Here we study general games among n strategies in populations of large but finite size. We explore stochastic evolutionary dynamics under weak selection, but for any mutation rate. We analyze the frequency dependent Moran process in well-mixed populations, but almost identical results are found for the Wright-Fisher and Pairwise Comparison processes. Surprisingly simple conditions specify whether a strategy is more abundant on average than 1/n, or than another strategy, in the mutation-selection equilibrium. We find one condition that holds for low mutation rate and another condition that holds for high mutation rate. A linear combination of these two conditions holds for any mutation rate. Our results allow a complete characterization of n×n games in the limit of weak selection.     

Early life relict feature in peptide mass distribution

The molecular mass of a biomolecule is characterized by the monoisitopic mass M_mono and the average isotopic mass M_av. We found that tryptic peptide masses mapped on a plane made by two parameters derived from M_mono and M_av form a peculiar feature in the form of a ‘band gap’ stretching across the whole ‘peptide galaxy’, with a narrow line in the centre. The purpose of this study was to investigate possible reasons for the emergence of such a feature, provided it is not a random occurrence. The a priori probability of such a feature to emerge by chance was found to be less than 1:100. Peptides contributing to the central line have elemental compositions following the rules S = 0; Z = C - (N + H)/2 = 0, which nine out of 20 amino acid residues satisfy. The relative abundances of amino acids in the peptides contributing to the central line correlate with the consensus order of emergence of these amino acids, with ancient amino acids being overrepresented in on-Line peptides. Since linear correlation between M_av and M_mono reduces the complexity of polypeptide molecules, and the turnover rate of less complex molecules should be faster in non-equilibrium abiotic synthesis, we hypothesize that the line could be a signature of abiotic production of primordial biopolymers. The linear dependence between the average isotopic masses and monoisotopic masses may have influenced the selection of amino acid residues for terrestrial life.     

Identification of key amino acid residues that determine the ability of high risk HPV16-E7 to dysregulate major histocompatibility complex class I expression

High risk human Papillomavirus (HPV) types are the major causative agents of cervical cancer. Reduced expression of major histocompatibility complex class I (MHC I) on HPV-infected cells might be responsible for insufficient T cell response and contribute to HPV-associated malignancy. The viral gene product required for subversion of MHC I synthesis is the E7 oncoprotein. Although it has been suggested that high and low risk HPVs diverge in their ability to dysregulate MHC I expression, it is not known what sequence determinants of HPV-E7 are responsible for this important functional difference. To investigate this, we analyzed the capability to affect MHC I of a set of chimeric E7 variants containing sequence elements from either high risk HPV16 or low risk HPV11. HPV16-E7, but not HPV11-E7, causes significant diminution of mRNA synthesis and surface presentation of MHC I, which depend on histone deacetylase activity. Our experiments demonstrate that the C-terminal region within the zinc finger domain of HPV-E7 is responsible for the contrasting effects of HPV11- and HPV16-E7 on MHC I. By using different loss- and gain-of-function mutants of HPV11- and HPV16-E7, we identify for the first time a residue variation at position 88 that is highly critical for HPV16-E7-mediated suppression of MHC I. Furthermore, our studies suggest that residues at position 78, 80, and 88 build a minimal functional unit within HPV16-E7 required for binding and histone deacetylase recruitment to the MHC I promoter. Taken together, our data provide new insights into how high risk HPV16-E7 dysregulates MHC I for immune evasion.