Next-generation sequencing: A frameshift in skeletal dysplasia gene discovery

In the last decade, huge breakthroughs in genetics - driven by new technology and different statistical approaches - have resulted in a plethora of new disease genes identified for both common and rare diseases. Massive parallel sequencing, commonly known as next-generation sequencing, is the latest advance in genetics, and has already facilitated the discovery of the molecular cause of many monogenic disorders. This article describes this new technology and reviews how this approach has been used successfully in patients with skeletal dysplasias. Moreover, this article illustrates how the study of rare diseases can inform understanding and therapeutic developments for common diseases such as osteoporosis. © International Osteoporosis Foundation and National Osteoporosis Foundation 2013.

Whole-exome re-sequencing in a family quartet identifies POP1 mutations as the cause of a novel skeletal dysplasia

Recent advances in DNA sequencing have enabled mapping of genes for monogenic traits in families with small pedigrees and even in unrelated cases. We report the identification of disease-causing mutations in a rare, severe, skeletal dysplasia, studying a family of two healthy unrelated parents and two affected children using whole-exome sequencing. The two affected daughters have clinical and radiographic features suggestive of anauxetic dysplasia (OMIM 607095), a rare form of dwarfism caused by mutations of RMRP. However, mutations of RMRP were excluded in this family by direct sequencing. Our studies identified two novel compound heterozygous loss-of-function mutations in POP1, which encodes a core component of the RNase mitochondrial RNA processing (RNase MRP) complex that directly interacts with the RMRP RNA domains that are affected in anauxetic dysplasia. We demonstrate that these mutations impair the integrity and activity of this complex and that they impair cell proliferation, providing likely molecular and cellular mechanisms by which POP1 mutations cause this severe skeletal dysplasia. © 2011 Glazov et al.

Whole-genome sequencing identifies EN1 as a determinant of bone density and fracture

The extent to which low-frequency (minor allele frequency (MAF) between 1-5%) and rare (MAF sequencing (n = 2,882 from UK10K (ref. 10); a population-based genome sequencing consortium), whole-exome sequencing (n = 3,549), deep imputation of genotyped samples using a combined UK10K/1000 Genomes reference panel (n = 26,534), and de novo replication genotyping (n = 20,271). We identified a low-frequency non-coding variant near a novel locus, EN1, with an effect size fourfold larger than the mean of previously reported common variants for lumbar spine BMD (rs11692564(T), MAF = 1.6%, replication effect size = +0.20 s.d., Pmeta = 2 x 10(-14)), which was also associated with a decreased risk of fracture (odds ratio = 0.85; P = 2 x 10(-11); ncases = 98,742 and ncontrols = 409,511). Using an En1(cre/flox) mouse model, we observed that conditional loss of En1 results in low bone mass, probably as a consequence of high bone turnover. We also identified a novel low-frequency non-coding variant with large effects on BMD near WNT16 (rs148771817(T), MAF = 1.2%, replication effect size = +0.41 s.d., Pmeta = 1 x 10(-11)). In general, there was an excess of association signals arising from deleterious coding and conserved non-coding variants. These findings provide evidence that low-frequency non-coding variants have large effects on BMD and fracture, thereby providing rationale for whole-genome sequencing and improved imputation reference panels to study the genetic architecture of complex traits and disease in the general population.

Characterization of the restriction enzyme-like endonuclease encoded by the Entamoeba histolytica non-long terminal repeat retrotransposon EhLINE1

The genome of the human pathogen Entamoeba histolytica, a primitive protist, contains non-long terminal repeat retrotransposable elements called EhLINEs. These encode reverse transcriptase and endonuclease required for retrotransposition. The endonuclease shows sequence similarity with bacterial restriction endonucleases. Here we report the salient enzymatic features of one such endonuclease. The kinetics of an EhLINE1-encoded endonuclease catalyzed reaction, determined under steady-state and single-turnover conditions, revealed a significant burst phase followed by a slower steady-state phase, indicating that release of product could be the slower step in this reaction. For circular supercoiled DNA the K-m was 2.6 x 10-8 m and the k(cat) was 1.6 x 10-2 sec-1. For linear E. histolytica DNA substrate the K-m and k(cat) values were 1.3 x 10-8 m and 2.2 x 10-4 sec-1 respectively. Single-turnover reaction kinetics suggested a noncooperative mode of hydrolysis. The enzyme behaved as a monomer. While Mg2+ was required for activity, 60% activity was seen with Mn2+ and none with other divalent metal ions. Substitution of PDX12-14D (a metal-binding motif) with PAX(12-14)D caused local conformational change in the protein tertiary structure, which could contribute to reduced enzyme activity in the mutated protein. The protein underwent conformational change upon the addition of DNA, which is consistent with the known behavior of restriction endonucleases. The similarities with bacterial restriction endonucleases suggest that the EhLINE1-encoded endonuclease was possibly acquired from bacteria through horizontal gene transfer. The loss of strict sequence specificity for nicking may have been subsequently selected to facilitate spread of the retrotransposon to intergenic regions of the E. histolytica genome.

Next generation sequencing identifies novel CACNA1A gene mutations in Episodic Ataxia type 2

Episodic Ataxia type 2 (EA2) is a rare autosomal dominantly inherited neurological disorder characterized by recurrent disabling imbalance, vertigo and episodes of ataxia lasting minutes to hours. EA2 is caused most often by loss of function mutations of the calcium channel gene CACNA1A. In addition to EA2, mutations in CACNA1A are responsible for two other allelic disorders: familial hemiplegic migraine type1 (FHM1) and spinocerebellar ataxia type 6 (SCA6). Herein, we have utilised Next Generation Sequencing (NGS) to screen the coding sequence, exon-intron boundaries and UTRs of five genes where mutation is known to produce symptoms related to EA2, including CACNA1A. We performed this screening in a group of 31 unrelated patients with EA2 symptoms. Both novel and known mutations were detected through NGS technology, and confirmed through Sanger sequencing. Genetic testing showed in total 15 mutation bearing patients (48%), of which 9 were novel mutations (6 missense and 3 small frameshift deletion mutations) and six known mutations (4 missense and 2 nonsense).These results demonstrate the efficiency of our NGS-panel for detecting known and novel mutations for EA2 in the CACNA1A gene, also identifying a novel missense mutation in ATP1A2 which is not a normal target for EA2 screening.

Mycobacterium tuberculosis FtsZ requires at least one arginine residue at the C-terminal end for polymerization in vitro

We examined whether C-terminal residues of soluble recombinant FtsZ of Mycobacterium tuberculosis (MtFtsZ) have any role in MtFtsZ polymerization in vitro. MtFtsZ-delta C1, which lacks C-terminal extreme Arg residue (underlined in the C-terminal extreme stretch of 13 residues, DDDDVDVPPFMRR), but retaining the penultimate Arg residue (DDDDVDVPPFMR), polymerizes like full-length MtFtsZ in vitro. However, MtFtsZ-delta C2 that lacks both the Arg residues at the C-terminus (DDDDVDVPPFM), neither polymerizes at pH 6.5 nor forms even single- or double-stranded filaments at pH 7.7 in the presence of 10 mM CaCl2. Neither replacement of the penultimate Arg residue, in the C-terminal Arg deletion mutant DDDDVDVPPFMR, with Lys or His or Ala or Asp (DDDDVDVPPFMK/H/A/D) enabled polymerization. Although MtFtsZ-delta C2 showed secondary and tertiary structural changes, which might have affected polymerization, GTPase activity of MtFtsZ-delta C2 was comparable to that of MtFtsZ. These data suggest that MtFtsZ requires an Arg residue as the extreme C-terminal residue for polymerization in vitro. The polypeptide segment containing C-terminal 67 residues, whose coordinates were absent from MtFtsZ crystal structure, was modeled on tubulin and MtFtsZ dimers. Possibilities for the influence of the C-terminal Arg residues on the stability of the dimer and thereby on MtFtsZ polymerization have been discussed.

Characterization of single-stranded DNA-binding proteins from Mycobacteria. The carboxyl-terminal of domain of SSB is essential for stable association with its cognate RecA protein

Single-stranded DNA-binding proteins (SSB) play an important role in most aspects of DNA metabolism including DNA replication, repair, and recombination. We report here the identification and characterization of SSB proteins of Mycobacterium smegmatis and Mycobacterium tuberculosis. Sequence comparison of M. smegmatis SSB revealed that it is homologous to M. tuberculosis SSB, except for a small spacer connecting the larger amino-terminal domain with the extreme carboxyl-terminal tail. The purified SSB proteins of mycobacteria bound single-stranded DNA with high affinity, and the association and dissociation constants were similar to that of the prototype SSB. The proteolytic signatures of free and bound forms of SSB proteins disclosed that DNA binding was associated with structural changes at the carboxyl-terminal domain. Significantly, SSB proteins from mycobacteria displayed high affinity for cognate RecA, whereas Escherichia coli SSB did not under comparable experimental conditions. Accordingly, SSB and RecA were coimmunoprecipitated from cell lysates, further supporting an interaction between these proteins in vivo. The carboxyl-terminal domain of M. smegmatis SSB, which is not essential for interaction with ssDNA, is the site of binding of its cognate RecA. These studies provide the first evidence for stable association of eubacterial SSB proteins with their cognate RecA, suggesting that these two proteins might function together during DNA repair and/or recombination.

Ultra high throughput sequencing excludes MDH1 as candidate gene for RP28-linked retinitis pigmentosa

Purpose: Mutations in IDH3B, an enzyme participating in the Krebs cycle, have recently been found to cause autosomal recessive retinitis pigmentosa (arRP). The MDH1 gene maps within the RP28 arRP linkage interval and encodes cytoplasmic malate dehydrogenase, an enzyme functionally related to IDH3B. As a proof of concept for candidate gene screening to be routinely performed by ultra high throughput sequencing (UHTs), we analyzed MDH1 in a patient from each of the two families described so far to show linkage between arRP and RP28. Methods: With genomic long-range PCR, we amplified all introns and exons of the MDH1 gene (23.4 kb). PCR products were then sequenced by short-read UHTs with no further processing. Computer-based mapping of the reads and mutation detection were performed by three independent software packages. Results: Despite the intrinsic complexity of human genome sequences, reads were easily mapped and analyzed, and all algorithms used provided the same results. The two patients were homozygous for all DNA variants identified in the region, which confirms previous linkage and homozygosity mapping results, but had different haplotypes, indicating genetic or allelic heterogeneity. None of the DNA changes detected could be associated with the disease. Conclusions: The MDH1 gene is not the cause of RP28-linked arRP. Our experimental strategy shows that long-range genomic PCR followed by UHTs provides an excellent system to perform a thorough screening of candidate genes for hereditary retinal degeneration.

Peptide mimics for structural features in proteins. Crystal structures of three heptapeptide helices with a C-terminal 6-->1 hydrogen bond.

The crystal structure determination of three heptapeptides containing alpha-aminoisobutyryl (Aib) residues as a means of helix stabilization provides a high-resolution characterization of 6-->1 hydrogen-bonded conformations, reminiscent of helix-terminating structural features in proteins. The crystal parameters for the three peptides, Boc-Val-Aib-X-Aib-Ala-Aib-Y-OMe, where X and Y are Phe, Leu (I), Leu, Phe (II) and Leu, Leu (III) are: (I) space group P1, Z = 1, a = 9.903 A, b = 10.709 A, c = 11.969 A, alpha = 102.94 degrees, beta = 103.41 degrees, gamma = 92.72 degrees, R = 4.55%; (II) space group P21, Z = 2, a = 10.052 A, b = 17.653 A, c = 13.510 A, beta = 108.45 degrees, R = 4.49%; (III) space group P1, Z = 2 (two independent molecules IIIa and IIIb in the asymmetric unit), a = 10.833 A, b = 13.850 A, c = 16.928 A, alpha = 99.77 degrees, beta = 105.90 degrees, gamma = 90.64 degrees, R = 8.54%. In all cases the helices form 3(10)/alpha-helical (or 3(10)helical) structures, with helical columns formed by head-to-tail hydrogen bonding. The helices assemble in an all-parallel motif in crystals I and III and in an antiparallel motif in II. In the four crystallographically characterized molecules, I, II, IIIa and IIIb, Aib(6) adopts a left-handed helical (hL) conformation with positive phi, psi values, resulting in 6-->1 hydrogen-bond formation between Aib(2) CO and Leu(7)/Phe(7) NH groups. In addition a 4-->1 hydrogen bond is seen between Aib(3) CO and Aib(6) NH groups. This pattern of hydrogen bonding is often observed at the C-terminus of helices proteins, with the terminal pi-type turn being formed by four residues adopting the hRhRhRhL conformation.

Sequencing of T-superfamily conotoxins from Conus virgo: Pyroglutamic acid identification and disulfide arrangement by MALDI mass spectrometry

De novo mass spectrometric sequencing of two Conus peptides, Vi1359 and Vi1361, from the vermivorous cone snail Conus virgo, found off the southern Indian coast, is presented. The peptides, whose masses differ only by 2 Da, possess two disulfide bonds and an amidated C-terminus. Simple chemical modifications and enzymatic cleavage coupled with matrix assisted laser desorption ionization (MALDI) mass spectrometric analysis aided in establishing the sequences of Vi1359, ZCCITIPECCRI-NH2, and Vi1361, ZCCPTMPECCRI-NH2, Which differ only at residues 4 and 6 (Z = pyroglutamic acid). The presence of the pyroglutamyl residue at the N-terminus was unambiguously identified by chemical hydrolysis of the cyclic amide, followed by esterification. The presence of Ile residues in both the peptides was confirmed from high-energy collision induced dissociation (CID) studies, using the observation Of W-n- and d(n)-ions as a diagnostic. Differential cysteine labeling, in conjunction with MALDI-MS/MS, permitted establishment of disulfide connectivity in both peptides as Cys2-Cys9 and Cys3-Cys10. The cysteine pattern clearly reveals that the peptides belong to the class of T-superfamily conotoxins, in particular the T-1 superfamily.

Amino acid sequence of the winged bean (Psophocarpus tetragonolobus) basic lectin. Adenine binding and identification of the active-site tryptophan residue

The complete amino acid sequence of winged bean basic agglutinin (WBA I) was obtained by a combination of manual and gas-phase sequencing methods. Peptide fragments for sequence analyses were obtained by enzymatic cleavages using trypsin and Staphylococcus aureus V8 endoproteinase and by chemical cleavages using iodosobenzoic acid, hydroxylamine, and formic acid. COOH-terminal sequence analysis of WBA I and other peptides was performed using carboxypeptidase Y. The primary structure of WBA I was homologous to those of other legume lectins and more so to Erythrina corallodendron. Interestingly, the sequence shows remarkable identities in the regions involved in the association of the two monomers of E. corallodendron lectin. Other conserved regions are the double metal-binding site and residues contributing to the formation of the hydrophobic cavity and the carbohydrate-binding site. Chemical modification studies both in the presence and absence of N-acetylgalactosamine together with sequence analyses of tryptophan-containing tryptic peptides demonstrate that tryptophan 133 is involved in the binding of carbohydrate ligands by the lectin. The location of tryptophan 133 at the active center of WBA I for the first time subserves to explain a role for one of the most conserved residues in legume lectins.

Ultra high throughput sequencing excludes MDH1 as candidate gene for RP28-linked retinitis pigmentosa

Purpose: Mutations in IDH3B, an enzyme participating in the Krebs cycle, have recently been found to cause autosomal recessive retinitis pigmentosa (arRP). The MDH1 gene maps within the RP28 arRP linkage interval and encodes cytoplasmic malate dehydrogenase, an enzyme functionally related to IDH3B. As a proof of concept for candidate gene screening to be routinely performed by ultra high throughput sequencing (UHTs), we analyzed MDH1 in a patient from each of the two families described so far to show linkage between arRP and RP28. Methods: With genomic long-range PCR, we amplified all introns and exons of the MDH1 gene (23.4 kb). PCR products were then sequenced by short-read UHTs with no further processing. Computer-based mapping of the reads and mutation detection were performed by three independent software packages. Results: Despite the intrinsic complexity of human genome sequences, reads were easily mapped and analyzed, and all algorithms used provided the same results. The two patients were homozygous for all DNA variants identified in the region, which confirms previous linkage and homozygosity mapping results, but had different haplotypes, indicating genetic or allelic heterogeneity. None of the DNA changes detected could be associated with the disease.

The transient response of a two-terminal solion to a step-current input

The transition time associated with the time-variation of the voltage across a two-terminal diaphragm-less solion in response to a step-current stimulus has been studied experimentally. A theoretical analysis has also been made by solving the diffusion problem under the appropriate initial and boundary conditions. The behaviour of the theoretically predicted transition times is in agreement with the observed behaviour. The systems under study have been shown to be different from those used hitherto in thin-layer chronopotentiometry.

The transient response of a two-terminal solion to a stepcrrent input

The transition time associated with the time-variation of the voltage across a two-terminal diaphragm-less solion in response to a step-current stimulus has been studied experimentally. A theoretical analysis has also been made by solving the diffusion problem under the appropriate initial and boundary conditions. The behaviour of the theoretically predicted transition times is in agreement with the observed behaviour. The systems under study have been shown to be different from those used hitherto in thin-layer chronopotentiometry.