Exome Sequencing and Optical Genome Mapping in Molecularly Unsolved Cases of Duchenne Muscular Dystrophy: Identification of a Causative X-Chromosomal Inversion Disrupting the DMD Gene.

Duchenne muscular dystrophy (DMD) is a severe progressive muscle disease that mainly affects boys due to X-linked recessive inheritance. In most affected individuals, MLPA or sequencing-based techniques detect deletions, duplications, or point mutations in the dystrophin-encoding DMD gene. However, in a small subset of patients clinically diagnosed with DMD, the molecular cause is not identified with these routine methods. Evaluation of the 60 DMD patients in our center revealed three cases without a known genetic cause. DNA samples of these patients were analyzed using whole-exome sequencing (WES) and, if unconclusive, optical genome mapping (OGM). WES led to a diagnosis in two cases: one patient was found to carry a splice mutation in the DMD gene that had not been identified during previous Sanger sequencing. In the second patient, we detected two variants in the fukutin gene (FKTN) that were presumed to be disease-causing. In the third patient, WES was unremarkable, but OGM identified an inversion disrupting the DMD gene (~1.28 Mb) that was subsequently confirmed with long-read sequencing. These results highlight the importance of reanalyzing unsolved cases using WES and demonstrate that OGM is a useful method for identifying large structural variants in cases with unremarkable exome sequencing.

Novel Filamin C Myofibrillar Myopathy Variants Cause Different Pathomechanisms and Alterations in Protein Quality Systems.

Myofibrillar myopathies (MFM) are a group of chronic muscle diseases pathophysiologically characterized by accumulation of protein aggregates and structural failure of muscle fibers. A subtype of MFM is caused by heterozygous mutations in the filamin C (FLNC) gene, exhibiting progressive muscle weakness, muscle structural alterations and intracellular protein accumulations. Here, we characterize in depth the pathogenicity of two novel truncating FLNc variants (p.Q1662X and p.Y2704X) and assess their distinct effect on FLNc stability and distribution as well as their impact on protein quality system (PQS) pathways. Both variants cause a slowly progressive myopathy with disease onset in adulthood, chronic myopathic alterations in muscle biopsy including the presence of intracellular protein aggregates. Our analyses revealed that p.Q1662X results in FLNc haploinsufficiency and p.Y2704X in a dominant-negative FLNc accumulation. Moreover, both protein-truncating variants cause different PQS alterations: p.Q1662X leads to an increase in expression of several genes involved in the ubiquitin-proteasome system (UPS) and the chaperone-assisted selective autophagy (CASA) system, whereas p.Y2704X results in increased abundance of proteins involved in UPS activation and autophagic buildup. We conclude that truncating FLNC variants might have different pathogenetic consequences and impair PQS function by diverse mechanisms and to varying extents. Further studies on a larger number of patients are necessary to confirm our observations.

Integrative functional genomics decodes herpes simplex virus 1.

The predicted 80 open reading frames (ORFs) of herpes simplex virus 1 (HSV-1) have been intensively studied for decades. Here, we unravel the complete viral transcriptome and translatome during lytic infection with base-pair resolution by computational integration of multi-omics data. We identify a total of 201 transcripts and 284 ORFs including all known and 46 novel large ORFs. This includes a so far unknown ORF in the locus deleted in the FDA-approved oncolytic virus Imlygic. Multiple transcript isoforms expressed from individual gene loci explain translation of the vast majority of ORFs as well as N-terminal extensions (NTEs) and truncations. We show that NTEs with non-canonical start codons govern the subcellular protein localization and packaging of key viral regulators and structural proteins. We extend the current nomenclature to include all viral gene products and provide a genome browser that visualizes all the obtained data from whole genome to single-nucleotide resolution.

A validation study revealed differences in design and performance of MEDLINE search filters for qualitative research.

OBJECTIVES: Several search filters exist to identify qualitative research, but so far none of them has been validated with an independent set of relevant references irrespective of a medical topic. The objective of this study was to provide a comparative overview of validation results for various MEDLINE search filters. STUDY DESIGN AND SETTING: Search filters were tested for plausibility. A relative recall approach was used to generate a gold standard based on an overview of systematic reviews of qualitative studies. For each review, the included qualitative studies were collected and checked for MEDLINE-indexing. The body of indexed articles yielded the gold standard. Validation tests were conducted to determine whether the references of the gold standard could be identified with the respective search filters. RESULTS: Thirteen search filters were validated in MEDLINE. One search filter by Wong et al. (2004) was found to be the most sensitive (93.63%). While medical subject heading "qualitative research" achieved the best precision (2.15%), sensitivity was the lowest (22.56%). University of Texas provided the best balanced search filter with a sensitivity of 81.96% and a precision of 0.80%. CONCLUSION: Search filters to identify qualitative research in MEDLINE differ greatly in design and performance. The selection of the appropriate search filter depends on project-specific demands and resources.

Efficient depletion of ribosomal RNA for RNA sequencing in planarians.

BACKGROUND: The astounding regenerative abilities of planarian flatworms prompt steadily growing interest in examining their molecular foundation. Planarian regeneration was found to require hundreds of genes and is hence a complex process. Thus, RNA interference followed by transcriptome-wide gene expression analysis by RNA-seq is a popular technique to study the impact of any particular planarian gene on regeneration. Typically, the removal of ribosomal RNA (rRNA) is the first step of all RNA-seq library preparation protocols. To date, rRNA removal in planarians was primarily achieved by the enrichment of polyadenylated (poly(A)) transcripts. However, to better reflect transcriptome dynamics and to cover also non-poly(A) transcripts, a procedure for the targeted removal of rRNA in planarians is needed. RESULTS: In this study, we describe a workflow for the efficient depletion of rRNA in the planarian model species S. mediterranea. Our protocol is based on subtractive hybridization using organism-specific probes. Importantly, the designed probes also deplete rRNA of other freshwater triclad families, a fact that considerably broadens the applicability of our protocol. We tested our approach on total RNA isolated from stem cells (termed neoblasts) of S. mediterranea and compared ribodepleted libraries with publicly available poly(A)-enriched ones. Overall, mRNA levels after ribodepletion were consistent with poly(A) libraries. However, ribodepleted libraries revealed higher transcript levels for transposable elements and histone mRNAs that remained underrepresented in poly(A) libraries. As neoblasts experience high transposon activity this suggests that ribodepleted libraries better reflect the transcriptional dynamics of planarian stem cells. Furthermore, the presented ribodepletion procedure was successfully expanded to the removal of ribosomal RNA from the gram-negative bacterium Salmonella typhimurium. CONCLUSIONS: The ribodepletion protocol presented here ensures the efficient rRNA removal from low input total planarian RNA, which can be further processed for RNA-seq applications. Resulting libraries contain less than 2% rRNA. Moreover, for a cost-effective and efficient removal of rRNA prior to sequencing applications our procedure might be adapted to any prokaryotic or eukaryotic species of choice.

GPD1 Specifically Marks Dormant Glioma Stem Cells with a Distinct Metabolic Profile.

Brain tumor stem cells (BTSCs) are a chemoresistant population that can drive tumor growth and relapse, but the lack of BTSC-specific markers prevents selective targeting that spares resident stem cells. Through a ribosome-profiling analysis of mouse neural stem cells (NSCs) and BTSCs, we find glycerol-3-phosphate dehydrogenase 1 (GPD1) expression specifically in BTSCs and not in NSCs. GPD1 expression is present in the dormant BTSC population, which is enriched at tumor borders and drives tumor relapse after chemotherapy. GPD1 inhibition prolongs survival in mouse models of glioblastoma in part through altering cellular metabolism and protein translation, compromising BTSC maintenance. Metabolomic and lipidomic analyses confirm that GPD1+ BTSCs have a profile distinct from that of NSCs, which is dependent on GPD1 expression. Similar GPD1 expression patterns and prognostic associations are observed in human gliomas. This study provides an attractive therapeutic target for treating brain tumors and new insights into mechanisms regulating BTSC dormancy.

Selective ribosome profiling to study interactions of translating ribosomes in yeast.

A number of enzymes, targeting factors and chaperones engage ribosomes to support fundamental steps of nascent protein maturation, including enzymatic processing, membrane targeting and co-translational folding. The selective ribosome profiling (SeRP) method is a new tool for studying the co-translational activity of maturation factors that provides proteome-wide information on a factor's nascent interactome, the onset and duration of binding and the mechanisms controlling factor engagement. SeRP is based on the combination of two ribosome-profiling (RP) experiments, sequencing the ribosome-protected mRNA fragments from all ribosomes (total translatome) and the ribosome subpopulation engaged by the factor of interest (factor-bound translatome). We provide a detailed SeRP protocol, exemplified for the yeast Hsp70 chaperone Ssb (stress 70 B), for studying factor interactions with nascent proteins that is readily adaptable to identifying nascent interactomes of other co-translationally acting eukaryotic factors. The protocol provides general guidance for experimental design and optimization, as well as detailed instructions for cell growth and harvest, the isolation of (factor-engaged) monosomes, the generation of a cDNA library and data analysis. Experience in biochemistry and RNA handling, as well as basic programing knowledge, is necessary to perform SeRP. Execution of a SeRP experiment takes 8-10 working days, and initial data analysis can be completed within 1-2 d. This protocol is an extension of the originally developed protocol describing SeRP in bacteria.

Cotranslational assembly of protein complexes in eukaryotes revealed by ribosome profiling.

The folding of newly synthesized proteins to the native state is a major challenge within the crowded cellular environment, as non-productive interactions can lead to misfolding, aggregation and degradation1. Cells cope with this challenge by coupling synthesis with polypeptide folding and by using molecular chaperones to safeguard folding cotranslationally2. However, although most of the cellular proteome forms oligomeric assemblies3, little is known about the final step of folding: the assembly of polypeptides into complexes. In prokaryotes, a proof-of-concept study showed that the assembly of heterodimeric luciferase is an organized cotranslational process that is facilitated by spatially confined translation of the subunits encoded on a polycistronic mRNA4. In eukaryotes, however, fundamental differences-such as the rarity of polycistronic mRNAs and different chaperone constellations-raise the question of whether assembly is also coordinated with translation. Here we provide a systematic and mechanistic analysis of the assembly of protein complexes in eukaryotes using ribosome profiling. We determined the in vivo interactions of the nascent subunits from twelve hetero-oligomeric protein complexes of Saccharomyces cerevisiae at near-residue resolution. We find nine complexes assemble cotranslationally; the three complexes that do not show cotranslational interactions are regulated by dedicated assembly chaperones5-7. Cotranslational assembly often occurs uni-directionally, with one fully synthesized subunit engaging its nascent partner subunit, thereby counteracting its propensity for aggregation. The onset of cotranslational subunit association coincides directly with the full exposure of the nascent interaction domain at the ribosomal tunnel exit. The action of the ribosome-associated Hsp70 chaperone Ssb8 is coordinated with assembly. Ssb transiently engages partially synthesized interaction domains and then dissociates before the onset of partner subunit association, presumably to prevent premature assembly interactions. Our study shows that cotranslational subunit association is a prevalent mechanism for the assembly of hetero-oligomers in yeast and indicates that translation, folding and the assembly of protein complexes are integrated processes in eukaryotes.

Role for ribosome-associated complex and stress-seventy subfamily B (RAC-Ssb) in integral membrane protein translation.

Targeting of most integral membrane proteins to the endoplasmic reticulum is controlled by the signal recognition particle, which recognizes a hydrophobic signal sequence near the protein N terminus. Proper folding of these proteins is monitored by the unfolded protein response and involves protein degradation pathways to ensure quality control. Here, we identify a new pathway for quality control of major facilitator superfamily transporters that occurs before the first transmembrane helix, the signal sequence recognized by the signal recognition particle, is made by the ribosome. Increased rates of translation elongation of the N-terminal sequence of these integral membrane proteins can divert the nascent protein chains to the ribosome-associated complex and stress-seventy subfamily B chaperones. We also show that quality control of integral membrane proteins by ribosome-associated complex-stress-seventy subfamily B couples translation rate to the unfolded protein response, which has implications for understanding mechanisms underlying human disease and protein production in biotechnology.

Profiling Ssb-Nascent Chain Interactions Reveals Principles of Hsp70-Assisted Folding.

The yeast Hsp70 chaperone Ssb interacts with ribosomes and nascent polypeptides to assist protein folding. To reveal its working principle, we determined the nascent chain-binding pattern of Ssb at near-residue resolution by in vivo selective ribosome profiling. Ssb associates broadly with cytosolic, nuclear, and hitherto unknown substrate classes of mitochondrial and endoplasmic reticulum (ER) nascent proteins, supporting its general chaperone function. Ssb engages most substrates by multiple binding-release cycles to a degenerate sequence enriched in positively charged and aromatic amino acids. Timely association with this motif upon emergence at the ribosomal tunnel exit requires ribosome-associated complex (RAC) but not nascent polypeptide-associated complex (NAC). Ribosome footprint densities along orfs reveal faster translation at times of Ssb binding, mainly imposed by biases in mRNA secondary structure, codon usage, and Ssb action. Ssb thus employs substrate-tailored dynamic nascent chain associations to coordinate co-translational protein folding, facilitate accelerated translation, and support membrane targeting of organellar proteins.

SecA Cotranslationally Interacts with Nascent Substrate Proteins In Vivo.

SecA is an essential component of the Sec machinery in bacteria, which is responsible for transporting proteins across the cytoplasmic membrane. Recent work from our laboratory indicates that SecA binds to ribosomes. Here, we used two different approaches to demonstrate that SecA also interacts with nascent polypeptides in vivo and that these polypeptides are Sec substrates. First, we photo-cross-linked SecA to ribosomes in vivo and identified mRNAs that copurify with SecA. Microarray analysis of the copurifying mRNAs indicated a strong enrichment for proteins containing Sec-targeting sequences. Second, we used a 2-dimensional (2-D) gel approach to analyze radioactively labeled nascent polypeptides that copurify with SecA, including maltose binding protein, a well-characterized SecA substrate. The interaction of SecA with nascent chains was not strongly affected in cells lacking SecB or trigger factor, both of which also interact with nascent Sec substrates. Indeed, the ability of SecB to interact with nascent chains was disrupted in strains in which the interaction between SecA and the ribosome was defective. Analysis of the interaction of SecA with purified ribosomes containing arrested nascent chains in vitro indicates that SecA can begin to interact with a variety of nascent chains when they reach a length of ∼110 amino acids, which is considerably shorter than the length required for interaction with SecB. Our results suggest that SecA cotranslationally recognizes nascent Sec substrates and that this recognition could be required for the efficient delivery of these proteins to the membrane-embedded Sec machinery. IMPORTANCE: SecA is an ATPase that provides the energy for the translocation of proteins across the cytoplasmic membrane by the Sec machinery in bacteria. The translocation of most of these proteins is uncoupled from protein synthesis and is frequently described as "posttranslational." Here, we show that SecA interacts with nascent Sec substrates. This interaction is not dependent on SecB or trigger factor, which also interact with nascent Sec substrates. Moreover, the interaction of SecB with nascent polypeptides is dependent on the interaction of SecA with the ribosome, suggesting that interaction of the nascent chain with SecA precedes interaction with SecB. Our results suggest that SecA could recognize substrate proteins cotranslationally in order to efficiently target them for uncoupled protein translocation.

Structural analysis of a signal peptide inside the ribosome tunnel by DNP MAS NMR.

Proteins are synthesized in cells by ribosomes and, in parallel, prepared for folding or targeting. While ribosomal protein synthesis is progressing, the nascent chain exposes amino-terminal signal sequences or transmembrane domains that mediate interactions with specific interaction partners, such as the signal recognition particle (SRP), the SecA-adenosine triphosphatase, or the trigger factor. These binding events can set the course for folding in the cytoplasm and translocation across or insertion into membranes. A distinction of the respective pathways depends largely on the hydrophobicity of the recognition sequence. Hydrophobic transmembrane domains stabilize SRP binding, whereas less hydrophobic signal sequences, typical for periplasmic and outer membrane proteins, stimulate SecA binding and disfavor SRP interactions. In this context, the formation of helical structures of signal peptides within the ribosome was considered to be an important factor. We applied dynamic nuclear polarization magic-angle spinning nuclear magnetic resonance to investigate the conformational states of the disulfide oxidoreductase A (DsbA) signal peptide stalled within the exit tunnel of the ribosome. Our results suggest that the nascent chain comprising the DsbA signal sequence adopts an extended structure in the ribosome with only minor populations of helical structure.

Complex folding and misfolding effects of deer-specific amino acid substitutions in the ß2-α2 loop of murine prion protein.

The ß2-α2 loop of PrP(C) is a key modulator of disease-associated prion protein misfolding. Amino acids that differentiate mouse (Ser169, Asn173) and deer (Asn169, Thr173) PrP(C) appear to confer dramatically different structural properties in this region and it has been suggested that amino acid sequences associated with structural rigidity of the loop also confer susceptibility to prion disease. Using mouse recombinant PrP, we show that mutating residue 173 from Asn to Thr alters protein stability and misfolding only subtly, whilst changing Ser to Asn at codon 169 causes instability in the protein, promotes oligomer formation and dramatically potentiates fibril formation. The doubly mutated protein exhibits more complex folding and misfolding behaviour than either single mutant, suggestive of differential effects of the ß2-α2 loop sequence on both protein stability and on specific misfolding pathways. Molecular dynamics simulation of protein structure suggests a key role for the solvent accessibility of Tyr168 in promoting molecular interactions that may lead to prion protein misfolding. Thus, we conclude that 'rigidity' in the ß2-α2 loop region of the normal conformer of PrP has less effect on misfolding than other sequence-related effects in this region.