Cancer prognosis prediction using somatic point mutation and copy number variation data: a comparison of gene-level and pathway-based models

BACKGROUND

Genomic profiling of solid human tumors by projects such as The Cancer Genome Atlas (TCGA) has provided important information regarding the somatic alterations that drive cancer progression and patient survival. Although researchers have successfully leveraged TCGA data to build prognostic models, most efforts have focused on specific cancer types and a targeted set of gene-level predictors. Less is known about the prognostic ability of pathway-level variables in a pan-cancer setting. To address these limitations, we systematically evaluated and compared the prognostic ability of somatic point mutation (SPM) and copy number variation (CNV) data, gene-level and pathway-level models for a diverse set of TCGA cancer types and predictive modeling approaches.

RESULTS

We evaluated gene-level and pathway-level penalized Cox proportional hazards models using SPM and CNV data for 29 different TCGA cohorts. We measured predictive accuracy as the concordance index for predicting survival outcomes. Our comprehensive analysis suggests that the use of pathway-level predictors did not offer superior predictive power relative to gene-level models for all cancer types but had the advantages of robustness and parsimony. We identified a set of cohorts for which somatic alterations could not predict prognosis, and a unique cohort LGG, for which SPM data was more predictive than CNV data and the predictive accuracy is good for all model types. We found that the pathway-level predictors provide superior interpretative value and that there is often a serious collinearity issue for the gene-level models while pathway-level models avoided this issue.

CONCLUSION

Our comprehensive analysis suggests that when using somatic alterations data for cancer prognosis prediction, pathway-level models are more interpretable, stable and parsimonious compared to gene-level models. Pathway-level models also avoid the issue of collinearity, which can be serious for gene-level somatic alterations. The prognostic power of somatic alterations is highly variable across different cancer types and we have identified a set of cohorts for which somatic alterations could not predict prognosis. In general, CNV data predicts prognosis better than SPM data with the exception of the LGG cohort.

Seq2seq Fingerprint with Byte-Pair Encoding for Predicting Changes in Protein Stability upon Single Point Mutation

The engineering of stable proteins is crucial for various industrial purposes. Several machine learning methods have been developed to predict changes in the stability of proteins corresponding to single point mutations. To improve the prediction accuracy, we propose a new unsupervised descriptor for protein sequences, which is based on a sequence-to-sequence (seq2seq) neural network model combined with a sequence-compression method called byte-pair encoding (BPE). Our results demonstrate that BPE can encode a protein sequence into a sequence of shorter length, thereby enabling efficient training of the seq2seq model. Furthermore, we implement a basic predictor using the proposed descriptor, and our experimental results demonstrate that the predictor achieves state-of-the-art accuracy in tests for proteins that are not included in the training data.

Whole genome sequence-based haplotypes reveal a single origin of the 1393 bp HBB deletion

BACKGROUND

Mutations of HBB give rise to two prevalent haemoglobin disorders-sickle cell disease (SCD) and β-thalassaemia. While SCD is caused by a single base substitution, nearly 300 mutations that downregulate expression of HBB have been described. The vast majority of β-thalassaemia alleles are point mutations or small insertion/deletions within the HBB gene; deletions causing β-thalassaemia are very rare. We have identified three individuals with haemoglobin Sβ0-thalassaemia in which the β0-thalassaemia mutation is caused by a large deletion.

OBJECTIVE

To use whole genome sequence data to determine whether these deletions arose from a single origin.

METHODS

We used two approaches to confirm unrelatedness: pairwise comparison of SNPs and identity by descent analysis. Eagle, V.2.4, was used to generate phased haplotypes for the 683 individuals. The Neighbor-Net method implemented in SplitsTree V.4.13.1 was used to construct the network of haplotypes.

RESULTS

All three deletions involved 1393 bp, encompassing the β-promoter, exons 1 and 2, and part of intron 2, with identical breakpoints. The cases were confirmed to be unrelated. Haplotypes based on 29 SNPs in the HBB cluster showed that the three individuals harboured different βS haplotypes. In contrast, the haplotype harbouring the 1393 bp deletion was the same in all three individuals.

CONCLUSION

We suggest that all the reported cases of the 1393 bp HBB deletion, including the three cases here, are likely to be of the same ancestral origin.

Analyses of displacements resulting from a point mutation in proteins

The effects of a single residue substitution on the protein backbone are frequently quite small and there are many other potential sources of structural variation for protein. We present here a methodology considering different sources of distortions in order to isolate the very effect of the mutation. To validate our methodology, we consider a well-studied family with many single mutants: the human lysozyme. Most of the perturbations are expected to be at the very localisation of the mutation, but in many cases the effects are propagated at long range. We show that the distances between the mutated residue and the 5% most disturbed residues exponentially decreases. One third of the affected residues are in direct contact with the mutated position; the remaining two thirds are potential allosteric effects. We confirm the reliability of the residues identified as significantly perturbed by comparing our results to experimental studies. We confirm with the present method all the previously identified perturbations. This study shows that mutations have long-range impact on protein backbone that can be detected, although the displacement of the affected atoms is small.

A point mutation decouples the lipid transfer activities of microsomal triglyceride transfer protein

Apolipoprotein B-containing lipoproteins (B-lps) are essential for the transport of hydrophobic dietary and endogenous lipids through the circulation in vertebrates. Zebrafish embryos produce large numbers of B-lps in the yolk syncytial layer (YSL) to move lipids from yolk to growing tissues. Disruptions in B-lp production perturb yolk morphology, readily allowing for visual identification of mutants with altered B-lp metabolism. Here we report the discovery of a missense mutation in microsomal triglyceride transfer protein (Mtp), a protein that is essential for B-lp production. This mutation of a conserved glycine residue to valine (zebrafish G863V, human G865V) reduces B-lp production and results in yolk opacity due to aberrant accumulation of cytoplasmic lipid droplets in the YSL. However, this phenotype is milder than that of the previously reported L475P stalactite (stl) mutation. MTP transfers lipids, including triglycerides and phospholipids, to apolipoprotein B in the ER for B-lp assembly. In vitro lipid transfer assays reveal that while both MTP mutations eliminate triglyceride transfer activity, the G863V mutant protein unexpectedly retains ~80% of phospholipid transfer activity. This residual phospholipid transfer activity of the G863V mttp mutant protein is sufficient to support the secretion of small B-lps, which prevents intestinal fat malabsorption and growth defects observed in the mttpstl/stl mutant zebrafish. Modeling based on the recent crystal structure of the heterodimeric human MTP complex suggests the G865V mutation may block triglyceride entry into the lipid-binding cavity. Together, these data argue that selective inhibition of MTP triglyceride transfer activity may be a feasible therapeutic approach to treat dyslipidemia and provide structural insight for drug design. These data also highlight the power of yolk transport studies to identify proteins critical for B-lp biology.

Structural effects of the highly protective V127 polymorphism on human prion protein

Prion diseases, a group of incurable, lethal neurodegenerative disorders of mammals including humans, are caused by prions, assemblies of misfolded host prion protein (PrP). A single point mutation (G127V) in human PrP prevents prion disease, however the structural basis for its protective effect remains unknown. Here we show that the mutation alters and constrains the PrP backbone conformation preceding the PrP β-sheet, stabilising PrP dimer interactions by increasing intermolecular hydrogen bonding. It also markedly changes the solution dynamics of the β2-α2 loop, a region of PrP structure implicated in prion transmission and cross-species susceptibility. Both of these structural changes may affect access to protein conformers susceptible to prion formation and explain its profound effect on prion disease.

Structural Insights into the SPRED1-Neurofibromin-KRAS Complex and Disruption of SPRED1-Neurofibromin Interaction by Oncogenic EGFR

Sprouty-related, EVH1 domain-containing (SPRED) proteins negatively regulate RAS/mitogen-activated protein kinase (MAPK) signaling following growth factor stimulation. This inhibition of RAS is thought to occur primarily through SPRED1 binding and recruitment of neurofibromin, a RasGAP, to the plasma membrane. Here, we report the structure of neurofibromin (GTPase-activating protein [GAP]-related domain) complexed with SPRED1 (EVH1 domain) and KRAS. The structure provides insight into how the membrane targeting of neurofibromin by SPRED1 allows simultaneous interaction with activated KRAS. SPRED1 and NF1 loss-of-function mutations occur across multiple cancer types and developmental diseases. Analysis of the neurofibromin-SPRED1 interface provides a rationale for mutations observed in Legius syndrome and suggests why SPRED1 can bind to neurofibromin but no other RasGAPs. We show that oncogenic EGFR(L858R) signaling leads to the phosphorylation of SPRED1 on serine 105, disrupting the SPRED1-neurofibromin complex. The structural, biochemical, and biological results provide new mechanistic insights about how SPRED1 interacts with neurofibromin and regulates active KRAS levels in normal and pathologic conditions.

A tailored Cln3Q352X mouse model for testing therapeutic interventions in CLN3 Batten disease

CLN3 Batten disease (CLN3 disease) is a pediatric lysosomal storage disorder that presents with progressive blindness, motor and cognitive decline, seizures, and premature death. CLN3 disease results from mutations in CLN3 with the most prevalent mutation, a 966 bp deletion spanning exons 7-8, affecting ~ 75% of patients. Mouse models with complete Cln3 deletion or Cln3Δex7/8 mutation have been invaluable for learning about both the basic biology of CLN3 and the underlying pathological changes associated with CLN3 disease. These models, however, vary in their disease presentation and are limited in their utility for studying the role of nonsense mediated decay, and as a consequence, in testing nonsense suppression therapies and read-through compounds. In order to develop a model containing a disease-causing nonsense point mutation, here we describe a first-of-its-kind Cln3Q352X mouse model containing a c.1054C > T (p.Gln352Ter) point mutation. Similar to previously characterized Cln3 mutant mouse lines, this novel model shows pathological deficits throughout the CNS including accumulation of lysosomal storage material and glial activation, and has limited perturbation in behavioral measures. Thus, at the molecular and cellular level, this mouse line provides a valuable tool for testing nonsense suppression therapies or read through compounds in CLN3 disease in the future.

Functional Features of TREHALOSE-6-PHOSPHATE SYNTHASE1, an Essential Enzyme in Arabidopsis

In Arabidopsis (Arabidopsis thaliana), TREHALOSE-6-PHOSPHATE SYNTHASE1 (TPS1) catalyzes the synthesis of the sucrose-signaling metabolite trehalose 6-phosphate (Tre6P) and is essential for embryogenesis and normal postembryonic growth and development. To understand its molecular functions, we transformed the embryo-lethal tps1-1 null mutant with various forms of TPS1 and with a heterologous TPS (OtsA) from Escherichia coli, under the control of the TPS1 promoter, and tested for complementation. TPS1 protein localized predominantly in the phloem-loading zone and guard cells in leaves, root vasculature, and shoot apical meristem, implicating it in both local and systemic signaling of Suc status. The protein is targeted mainly to the nucleus. Restoring Tre6P synthesis was both necessary and sufficient to rescue the tps1-1 mutant through embryogenesis. However, postembryonic growth and the sucrose-Tre6P relationship were disrupted in some complementation lines. A point mutation (A119W) in the catalytic domain or truncating the C-terminal domain of TPS1 severely compromised growth. Despite having high Tre6P levels, these plants never flowered, possibly because Tre6P signaling was disrupted by two unidentified disaccharide-monophosphates that appeared in these plants. The noncatalytic domains of TPS1 ensure its targeting to the correct subcellular compartment and its catalytic fidelity and are required for appropriate signaling of Suc status by Tre6P.

Temperature-sensitive spinal muscular atrophy-causing point mutations lead to SMN instability, locomotor defects and premature lethality in Drosophila

Spinal muscular atrophy (SMA) is the leading genetic cause of death in young children, arising from homozygous deletion or mutation of the survival motor neuron 1 (SMN1) gene. SMN protein expressed from a paralogous gene, SMN2, is the primary genetic modifier of SMA; small changes in overall SMN levels cause dramatic changes in disease severity. Thus, deeper insight into mechanisms that regulate SMN protein stability should lead to better therapeutic outcomes. Here, we show that SMA patient-derived missense mutations in the Drosophila SMN Tudor domain exhibit a pronounced temperature sensitivity that affects organismal viability, larval locomotor function and adult longevity. These disease-related phenotypes are domain specific and result from decreased SMN stability at elevated temperature. This system was utilized to manipulate SMN levels during various stages of Drosophila development. Owing to a large maternal contribution of mRNA and protein, Smn is not expressed zygotically during embryogenesis. Interestingly, we find that only baseline levels of SMN are required during larval stages, whereas high levels of the protein are required during pupation. This previously uncharacterized period of elevated SMN expression, during which the majority of adult tissues are formed and differentiated, could be an important and translationally relevant developmental stage in which to study SMN function. Taken together, these findings illustrate a novel in vivo role for the SMN Tudor domain in maintaining SMN homeostasis and highlight the necessity for high SMN levels at crucial developmental time points that are conserved from Drosophila to humans.

Mutational signatures are jointly shaped by DNA damage and repair

Cells possess an armamentarium of DNA repair pathways to counter DNA damage and prevent mutation. Here we use C. elegans whole genome sequencing to systematically quantify the contributions of these factors to mutational signatures. We analyse 2,717 genomes from wild-type and 53 DNA repair defective backgrounds, exposed to 11 genotoxins, including UV-B and ionizing radiation, alkylating compounds, aristolochic acid, aflatoxin B1, and cisplatin. Combined genotoxic exposure and DNA repair deficiency alters mutation rates or signatures in 41% of experiments, revealing how different DNA alterations induced by the same genotoxin are mended by separate repair pathways. Error-prone translesion synthesis causes the majority of genotoxin-induced base substitutions, but averts larger deletions. Nucleotide excision repair prevents up to 99% of point mutations, almost uniformly across the mutation spectrum. Our data show that mutational signatures are joint products of DNA damage and repair and suggest that multiple factors underlie signatures observed in cancer genomes.

Antagonism of PP2A is an independent and conserved function of HIV-1 Vif and causes cell cycle arrest

The seminal description of the cellular restriction factor APOBEC3G and its antagonism by HIV-1 Vif has underpinned two decades of research on the host-virus interaction. We recently reported that HIV-1 Vif is also able to degrade the PPP2R5 family of regulatory subunits of key cellular phosphatase PP2A (PPP2R5A-E; Greenwood et al., 2016; Naamati et al., 2019). We now identify amino acid polymorphisms at positions 31 and 128 of HIV-1 Vif which selectively regulate the degradation of PPP2R5 family proteins. These residues covary across HIV-1 viruses in vivo, favouring depletion of PPP2R5A-E. Through analysis of point mutants and naturally occurring Vif variants, we further show that degradation of PPP2R5 family subunits is both necessary and sufficient for Vif-dependent G2/M cell cycle arrest. Antagonism of PP2A by HIV-1 Vif is therefore independent of APOBEC3 family proteins, and regulates cell cycle progression in HIV-infected cells.

A Single Amino Acid Substitution in STKc_GSK3 Kinase Conferring Semispherical Grains and Its Implications for the Origin of Triticum sphaerococcum

Six subspecies of hexaploid wheat (Triticum aestivum) have been identified, but the origin of Indian dwarf wheat (Triticum sphaerococcum), the only subspecies with round grains, is currently unknown. Here, we isolated the grain-shape gene Tasg-D1 in T sphaerococcum via positional cloning. Tasg-D1 encodes a Ser/Thr protein kinase glycogen synthase kinase3 (STKc_GSK3) that negatively regulates brassinosteroid signaling. Expression of TaSG-D1 and the mutant form Tasg-D1 in Arabidopsis (Arabidopsis thaliana) suggested that a single amino acid substitution in the Thr-283-Arg-284-Glu-285-Glu-286 domain of TaSG-D1 enhances protein stability in response to brassinosteroids, likely leading to formation of round grains in wheat. This gain-of-function mutation has pleiotropic effects on plant architecture and exhibits incomplete dominance. Haplotype analysis of 898 wheat accessions indicated that the origin of T sphaerococcum in ancient India involved at least two independent mutations of TaSG-D1 Our results demonstrate that modest genetic changes in a single gene can induce dramatic phenotypic changes.

A point mutation in LTT1 enhances cold tolerance at the booting stage in rice

The cold tolerance of rice at the booting stage is a main factor determining sustainability and regional adaptability. However, relatively few cold tolerance genes have been identified that can be effectively used in breeding programmes. Here, we show that a point mutation in the low-temperature tolerance 1 (LTT1) gene improves cold tolerance by maintaining tapetum degradation and pollen development, by activation of systems that metabolize reactive oxygen species (ROS). Cold-induced ROS accumulation is therefore prevented in the anthers of the ltt1 mutants allowing correct development. In contrast, exposure to cold stress dramatically increases ROS accumulation in the wild type anthers, together with the expression of genes encoding proteins associated with programmed cell death and with the accelerated degradation of the tapetum that ultimately leads to pollen abortion. These results demonstrate that appropriate ROS management is critical for the cold tolerance of rice at the booting stage. Hence, the ltt1 mutation can significantly improve the seed setting ability of cold-sensitive rice varieties under low-temperature stress conditions, with little yield penalty under optimal temperature conditions. This study highlights the importance of a valuable genetic resource that may be applied in rice breeding programmes to enhance cold tolerance.

Clinical significance of EGFR mutation types in lung adenocarcinoma: A multi-centre Korean study

Adenocarcinoma is the most common type of non-small cell lung cancer. Some causative genomic alterations in epidermal growth factor receptor (EGFR), including deletions in exon 19 (E19 dels) and a point mutation in E21, are known to have favourable prognoses due to sensitivity to tyrosine kinase inhibitors; however, the prognoses of other uncommon mutations are unclear. This study analysed the clinical significance of EGFR mutation types in lung adenocarcinoma. We retrospectively reviewed 1,020 subjects (mean age: 66.8 years, female: 41.7%) who were diagnosed with advanced lung adenocarcinoma, had EGFR mutation data, and did not undergo surgery from five medical institutes between 2010 and 2016. Subjects were classified according to EGFR mutation status, particularly for exon-specific mutations. EGFR positivity was defined as the presence of mutation and EGFR negativity was defined as wild-type EGFR. EGFR positivity was 38.0%, with the incidence of mutations in E18, E19, E20, and E21 was 3.6%, 51.0%, 3.4%, and 42.0%, respectively. The EGFR positive group survived significantly longer than the negative group (p<0.001), and there was a significant difference in survival among the four EGFR mutation sites (p = 0.003); E19 dels were the only significant factor that lowered mortality (HR: 0.678, p = 0.002), while an E21 mutation was the prognostic factor associated with the most increased mortality (HR: 1.365, p = 0.015). Amongst EGFR positive subjects, the proportion of E19 dels in TKI-responders was significantly higher and that of E21 mutations significantly lower, compared with non-responders. In TKI treatment, mutations in E18 and E20 were not worse factors than the E21 L858R mutation. In conclusion, the presence of EGFR mutations in advanced lung adenocarcinoma can predict a good prognosis; E19 dels prospect to have a better prognosis than other mutations, while an E21 mutation is expected to increase mortality.

Deep2Full: Evaluating strategies for selecting the minimal mutational experiments for optimal computational predictions of deep mutational scan outcomes

Performing a complete deep mutational scan with all single point mutations may not be practical, and may not even be required, especially if predictive computational models can be developed. Computational models are however naive to cellular response in the myriads of assay-conditions. In a realistic paradigm of assay context-aware predictive hybrid models that combine minimal experimental data from deep mutational scans with structure, sequence information and computational models, we define and evaluate different strategies for choosing this minimal set. We evaluated the trivial strategy of a systematic reduction in the number of mutational studies from 85% to 15%, along with several others about the choice of the types of mutations such as random versus site-directed with the same 15% data completeness. Interestingly, the predictive capabilities by training on a random set of mutations and using a systematic substitution of all amino acids to alanine, asparagine and histidine (ANH) were comparable. Another strategy we explored, augmenting the training data with measurements of the same mutants at multiple assay conditions, did not improve the prediction quality. For the six proteins we analyzed, the bin-wise error in prediction is optimal when 50-100 mutations per bin are used in training the computational model, suggesting that good prediction quality may be achieved with a library of 500-1000 mutations.

Protein Thermal Stability Engineering Using HoTMuSiC

The rational design of enzymes is a challenging research field, which plays an important role in the optimization of a wide series of biotechnological processes. Computational approaches allow screening all possible amino acid substitutions in a target protein and to identify a subset likely to have the desired properties. They can thus be used to guide and restrict the huge, time-consuming search in sequence space to reach protein optimality. Here we present HoTMuSiC, a tool that predicts the impact of point mutations on the protein melting temperature, which uses the experimental or modeled protein structure as sole input and is available at the dezyme.com website. Its main advantages include accuracy and speed, which makes it a perfect instrument for thermal stability engineering projects aiming at designing new proteins that feature increased heat resistance or remain active and stable in nonphysiological conditions. We set up a HoTMuSiC-based pipeline, which uses additional information to avoid mutations of functionally important residues, identified as being too well conserved among homologous proteins or too close to annotated functional sites. The efficiency of this pipeline is successfully demonstrated on Rhizomucor miehei lipase.

OsCIPK7 point-mutation leads to conformation and kinase-activity change for sensing cold response

Calcineurin B-like interacting protein kinases (CIPKs) play important roles via environmental stress. However, less is known how to sense the stress in molecular structure conformation level. Here, an OsCIPK7 mutant via TILLING procedure with a point mutation in the kinase domain showed increased chilling tolerance, which could be potentially used in the molecular breeding. We found that this point mutation of OsCIPK7 led to a conformational change in the activation loop of the kinase domain, subsequently with an increase of protein kinase activity, thus conferred an increased tolerance to chilling stress.

A single point mutation in the C-terminal extension of wheat Rubisco activase dramatically reduces ADP inhibition via enhanced ATP binding affinity

Ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) activase (Rca) is a AAA+ enzyme that uses ATP to remove inhibitors from the active site of Rubisco, the central carboxylation enzyme of photosynthesis. Rca α and β isoforms exist in most higher plant species, with the α isoform being identical to the β form but having an additional 25-45 amino acids at the Rca C terminus, known as the C-terminal extension (CTE). Rca is inhibited by ADP, and the extent of ADP sensitivity of the Rca complex can be modulated by the CTE of the α isoform, particularly in relation to a disulfide bond structure that is specifically reduced by the redox-regulatory enzyme thioredoxin-f. Here, we introduced single point mutations of Lys-428 in the CTE of Rca-α from wheat (Triticum aestivum) (TaRca2-α). Substitution of Lys-428 with Arg dramatically altered ADP inhibition, independently of thioredoxin-f regulation. We determined that the reduction in ADP inhibition in the K428R variant is not due to a change in ADP affinity, as the apparent constant for ADP binding was not altered by the K428R substitution. Rather, we observed that the K428R substitution strongly increased ATP substrate affinity and ATP-dependent catalytic velocity. These results suggest that the Lys-428 residue is involved in interacting with the γ-phosphate of ATP. Considering that nucleotide-dependent Rca activity regulates Rubisco and thus photosynthesis during fluctuating irradiance, the K428R substitution could potentially provide a mechanism for boosting the performance of wheat grown in the dynamic light environments of the field.

Predictive Value of Serum Antibodies and Point Mutations of AQP4, AQP1 and MOG in A Cohort of Spanish Patients with Neuromyelitis Optica Spectrum Disorders

The detection of IgG aquaporin-4 antibodies in the serum of patients with Neuromyelitis optica (NMO) has dramatically improved the diagnosis of this disease and its distinction from multiple sclerosis. Recently, a group of patients have been described who have an NMO spectrum disorder (NMOsd) and who are seronegative for AQP4 antibodies but positive for IgG aquaporin-1 (AQP1) or myelin oligodendrocyte glycoprotein (MOG) antibodies. The purpose of this study was to determine whether AQP1 and MOG could be considered new biomarkers of this disease; and if point mutations in the gDNA of AQP4, AQP1 and MOG genes could be associated with the etiology of NMOsd. We evaluated the diagnostic capability of ELISA and cell-based assays (CBA), and analyzed their reliability, specificity, and sensitivity in detecting antibodies against these three proteins. The results showed that both assays can recognize these antigen proteins under appropriate conditions, but only anti-AQP4 antibodies, and not AQP1 or MOG, appears to be a clear biomarker for NMOsd. CBA is the best method for detecting these antibodies; and serum levels of AQP4 antibodies do not correlate with the progression of this disease. So far, the sequencing analysis has not revealed a genetic basis for the etiology of NMOsd, but a more extensive analysis is required before definitive conclusions can be drawn.

A gene-edited mouse model of limb-girdle muscular dystrophy 2C for testing exon skipping

Limb-girdle muscular dystrophy type 2C is caused by autosomal recessive mutations in the γ-sarcoglycan (SGCG) gene. The most common SGCG mutation is a single nucleotide deletion from a stretch of five thymine residues in SGCG exon 6 (521ΔT). This founder mutation disrupts the transcript reading frame, abolishing protein expression. An antisense oligonucleotide exon-skipping method to reframe the human 521ΔT transcript requires skipping four exons to generate a functional, internally truncated protein. In vivo evaluation of this multi-exon skipping, antisense-mediated therapy requires a genetically appropriate mouse model. The human and mouse γ-sarcoglycan genes are highly homologous in sequence and gene structure, including the exon 6 region harboring the founder mutation. Herein, we describe a new mouse model of this form of limb-girdle muscular dystrophy generated using CRISPR/Cas9-mediated gene editing to introduce a single thymine deletion in murine exon 6, recreating the 521ΔT point mutation in Sgcg These mice express the 521ΔT transcript, lack γ-sarcoglycan protein and exhibit a severe dystrophic phenotype. Phenotypic characterization demonstrated reduced muscle mass, increased sarcolemmal leak and fragility, and decreased muscle function, consistent with the human pathological findings. Furthermore, we showed that intramuscular administration of a murine-specific multiple exon-directed antisense oligonucleotide cocktail effectively corrected the 521ΔT reading frame. These data demonstrate a molecularly and pathologically suitable model for in vivo testing of a multi-exon skipping strategy to advance preclinical development of this genetic correction approach.

Multi-individual microsatellite identification: A multiple genome approach to microsatellite design (MiMi)

Bespoke microsatellite marker panels are increasingly affordable and tractable to researchers and conservationists. The rate of microsatellite discovery is very high within a shotgun genomic data set, but extensive laboratory testing of markers is required for confirmation of amplification and polymorphism. By incorporating shotgun next-generation sequencing data sets from multiple individuals of the same species, we have developed a new method for the optimal design of microsatellite markers. This new tool allows us to increase the rate at which suitable candidate markers are selected by 58% in direct comparisons and facilitate an estimated 16% reduction in costs associated with producing a novel microsatellite panel. Our method enables the visualisation of each microsatellite locus in a multiple sequence alignment allowing several important quality checks to be made. Polymorphic loci can be identified and prioritised. Loci containing fragment-length-altering mutations in the flanking regions, which may invalidate assumptions regarding the model of evolution underlying variation at the microsatellite, can be avoided. Priming regions containing point mutations can be detected and avoided, helping to reduce sample-site-marker specificity arising from genetic isolation, and the likelihood of null alleles occurring. We demonstrate the utility of this new approach in two species: an echinoderm and a bird. Our method makes a valuable contribution towards minimising genotyping errors and reducing costs associated with developing a novel marker panel. The Python script to perform our method of multi-individual microsatellite identification (MiMi) is freely available from GitHub (https://github.com/graemefox/mimi).

Resurgent and Gating Pore Currents Induced by De Novo SCN2A Epilepsy Mutations

Over 150 mutations in the SCN2A gene, which encodes the neuronal Nav1.2 protein, have been implicated in human epilepsy cases. Of these, R1882Q and R853Q are two of the most commonly reported mutations. This study utilized voltage-clamp electrophysiology to characterize the biophysical effects of the R1882Q and R853Q mutations on the hNav1.2 channel, including their effects on resurgent current and gating pore current, which are not typically investigated in the study of Nav1.2 channel mutations. HEK cells transiently transfected with DNA encoding either wild-type (WT) or mutant hNav1.2 revealed that the R1882Q mutation induced a gain-of-function phenotype, including slowed fast inactivation, depolarization of the voltage dependence of inactivation, and increased persistent current. In this model system, the R853Q mutation primarily produced loss-of-function effects, including reduced transient current amplitude and density, hyperpolarization of the voltage dependence of inactivation, and decreased persistent current. The presence of a Navβ4 peptide (KKLITFILKKTREK-OH) in the pipette solution induced resurgent currents, which were increased by the R1882Q mutation and decreased by the R853Q mutation. Further study of the R853Q mutation in Xenopus oocytes indicated a reduced surface expression and revealed a robust gating pore current at negative membrane potentials, a function absent in the WT channel. This not only shows that different epileptogenic point mutations in hNav1.2 have distinct biophysical effects on the channel, but also illustrates that individual mutations can have complex consequences that are difficult to identify using conventional analyses. Distinct mutations may, therefore, require tailored pharmacotherapies in order to eliminate seizures.

A case-comparison study of pregnant women with mitochondrial disease - what to expect?

OBJECTIVE

Mitochondrial disease is a disorder of energy metabolism that affects 1 in 4300 adults in the UK. Pregnancy is associated with physiological demands that have implications for energy metabolism. We were interested to know how pregnancy was affected in women with mitochondrial disease, particularly those with the most common pathogenic mutation m.3243A>G.

DESIGN

Retrospective case-comparison study.

POPULATION/SETTING

Sixty-seven women with genetically confirmed mitochondrial disease from the UK Mitochondrial Diseases Cohort and 69 unaffected women participated.

METHODS

Participants answered questionnaires regarding each of their pregnancies. Patients were divided into two groups according to genetic mutation, with those harbouring m.3243A>G comprising a single group.

MAIN OUTCOME MEASURES

Pregnancy-related complications, mode of delivery, gestational age and birthweight of newborns.

RESULTS

Of 139 live births in the comparison group, 62 were in the m.3243A>G group and 87 were in the 'all other mutations' group. Pregnancies of women with the m.3243A>G mutation had significantly more gestational diabetes (odds ratio [OR] = 8.2, 95% CI 1.3-50.1), breathing difficulties (OR = 7.8, 95% CI 1.0-59.1) and hypertension (OR = 8.2, 95% CI 3.1-21.5) than the comparison group. Only half of the pregnancies in the m.3243A>G group had normal vaginal delivery, with emergency caesarean section accounting for 24.2% of deliveries. Babies were born significantly earlier to mothers harbouring m.3243A>G with 53.3% of them preterm (<37 weeks). These babies were also more likely to require resuscitation and admission.

CONCLUSION

Women who carried the m.3243A>G mutation appeared to be at higher risk of complications during pregnancies, caesarean section and preterm delivery than the unaffected women or those with other forms of mitochondrial disease.

TWEETABLE ABSTRACT

Pregnant women with mitochondrial disease - m.3243A>G mutation - are at greatly increased risk of complications and preterm delivery.

Conserved and novel roles of miR164-CUC2 regulatory module in specifying leaf and floral organ morphology in strawberry

MicroRNAs (miRNAs) are a kind of short noncoding RNA (20-24 nt), playing versatile roles in plant growth and development. Strawberry generates leaves and flowers with unique features. However, few miRNAs have been functionally characterised in strawberry, especially for their developmental regulation. Here, we identified one ethyl methanesulfonate (EMS) mutant, deeply serrated (des), in the woodland strawberry Fragaria vesca that has wrinkled leaves with deeper serrations, serrated petals and deformed carpels. The causative mutation occurs in the 19^th nucleotide of the FvemiR164a mature sequence. Overexpressing FveMIR164A rescued the phenotypes of des/fvemir164a except the petal serrations. Furthermore, we identified two allelic mutants of FveCUC2a, one target of FvemiR164a, which developed leaves with smooth margins and fused leaflets. Phenotypes of the double mutant fvemir164a fvecuc2a indicated that the two genes act linearly in leaf and carpel development, but synergistically in the development of other floral organs and inflorescence architecture. This work demonstrates the conserved and novel roles of the miR164-CUC2 module in leaf and flower development in different plant species, and reveals that the 19^th nucleotide of FvemiR164a is important for its processing.