New electrocardiographic criteria for discriminating between Brugada types 2 and 3 patterns and incomplete right bundle branch block.

OBJECTIVES: The aim of this study was to evaluate new electrocardiographic (ECG) criteria for discriminating between incomplete right bundle branch block (RBBB) and the Brugada types 2 and 3 ECG patterns. BACKGROUND: Brugada syndrome can manifest as either type 2 or type 3 pattern. The latter should be distinguished from incomplete RBBB, present in 3% of the population. METHODS: Thirty-eight patients with either type 2 or type 3 Brugada pattern that were referred for an antiarrhythmic drug challenge (AAD) were included. Before AAD, 2 angles were measured from ECG leads V(1) and/or V(2) showing incomplete RBBB: 1) α, the angle between a vertical line and the downslope of the r'-wave, and 2) β, the angle between the upslope of the S-wave and the downslope of the r'-wave. Baseline angle values, alone or combined with QRS duration, were compared between patients with negative and positive results on AAD. Receiver-operating characteristic curves were constructed to identify optimal discriminative cutoff values. RESULTS: The mean β angle was significantly smaller in the 14 patients with negative results on AAD compared to the 24 patients with positive results on AAD (36 ± 20° vs. 62 ± 20°, p < 0.01). Its optimal cutoff value was 58°, which yielded a positive predictive value of 73% and a negative predictive value of 87% for conversion to type 1 pattern on AAD; α was slightly less sensitive and specific compared with β. When the angles were combined with QRS duration, it tended to improve discrimination. CONCLUSIONS: In patients with suspected Brugada syndrome, simple ECG criteria can enable discrimination between incomplete RBBB and types 2 and 3 Brugada patterns.

Block of the superior cervical ganglion, description of a novel ultrasound-guided technique in human cadavers.

OBJECTIVE.: Injection of opioids to the superior cervical ganglion (SCG) has been reported to provide pain relief in patients suffering from different kinds of neuropathic facial pain conditions, such as trigeminal neuralgia, postherpetic neuralgia, and atypical facial pain. The classic approach to the SCG is a transoral technique using a so-called "stopper" to prevent accidental carotid artery puncture. The main disadvantage of this technique is that the needle tip is positioned distant from the actual target, possibly impeding successful block of the SCG. A further limitation is that injection of local anesthetics due to potential carotid artery puncture is contraindicated. We hypothesized that the SCG can be identified and blocked using ultrasound imaging, potentially increasing precision of this technique. INTERVENTIONS.: In this pilot study, 20 US-guided simulated blocks of the SCG were performed in 10 human cadavers in order to determine the accuracy of this novel block technique. After injection of 0.1 mL of dye, the cadavers were dissected to evaluate the needle position and coloring of the SCG. RESULTS.: Nineteen of the 20 needle tips were located in or next to the SCG. This corresponded to a simulated block success rate of 95% (95% confidence interval 85-100%). In 17 cases, the SCG was completely colored, and in two cases, the caudal half of the SCG was colored with dye. CONCLUSIONS.: The anatomical dissections confirmed that our ultrasound-guided approach to the SCG is accurate. Ultrasound could become an attractive alternative to the "blind" transoral technique of SCG blocks.

Exome sequencing identifies DYNC2H1 mutations as a common cause of asphyxiating thoracic dystrophy (Jeune syndrome) without major polydactyly, renal or retinal involvement.

BACKGROUND: Jeune asphyxiating thoracic dystrophy (JATD) is a rare, often lethal, recessively inherited chondrodysplasia characterised by shortened ribs and long bones, sometimes accompanied by polydactyly, and renal, liver and retinal disease. Mutations in intraflagellar transport (IFT) genes cause JATD, including the IFT dynein-2 motor subunit gene DYNC2H1. Genetic heterogeneity and the large DYNC2H1 gene size have hindered JATD genetic diagnosis. AIMS AND METHODS: To determine the contribution to JATD we screened DYNC2H1 in 71 JATD patients JATD patients combining SNP mapping, Sanger sequencing and exome sequencing. RESULTS AND CONCLUSIONS: We detected 34 DYNC2H1 mutations in 29/71 (41%) patients from 19/57 families (33%), showing it as a major cause of JATD especially in Northern European patients. This included 13 early protein termination mutations (nonsense/frameshift, deletion, splice site) but no patients carried these in combination, suggesting the human phenotype is at least partly hypomorphic. In addition, 21 missense mutations were distributed across DYNC2H1 and these showed some clustering to functional domains, especially the ATP motor domain. DYNC2H1 patients largely lacked significant extra-skeletal involvement, demonstrating an important genotype-phenotype correlation in JATD. Significant variability exists in the course and severity of the thoracic phenotype, both between affected siblings with identical DYNC2H1 alleles and among individuals with different alleles, which suggests the DYNC2H1 phenotype might be subject to modifier alleles, non-genetic or epigenetic factors. Assessment of fibroblasts from patients showed accumulation of anterograde IFT proteins in the ciliary tips, confirming defects similar to patients with other retrograde IFT machinery mutations, which may be of undervalued potential for diagnostic purposes.

Probabilistic base calling of Solexa sequencing data.

BACKGROUND: Solexa/Illumina short-read ultra-high throughput DNA sequencing technology produces millions of short tags (up to 36 bases) by parallel sequencing-by-synthesis of DNA colonies. The processing and statistical analysis of such high-throughput data poses new challenges; currently a fair proportion of the tags are routinely discarded due to an inability to match them to a reference sequence, thereby reducing the effective throughput of the technology. RESULTS: We propose a novel base calling algorithm using model-based clustering and probability theory to identify ambiguous bases and code them with IUPAC symbols. We also select optimal sub-tags using a score based on information content to remove uncertain bases towards the ends of the reads. CONCLUSION: We show that the method improves genome coverage and number of usable tags as compared with Solexa's data processing pipeline by an average of 15%. An R package is provided which allows fast and accurate base calling of Solexa's fluorescence intensity files and the production of informative diagnostic plots.

Exome sequencing of index patients with retinal dystrophies as a tool for molecular diagnosis.

BACKGROUND: Retinal dystrophies (RD) are a group of hereditary diseases that lead to debilitating visual impairment and are usually transmitted as a Mendelian trait. Pathogenic mutations can occur in any of the 100 or more disease genes identified so far, making molecular diagnosis a rather laborious process. In this work we explored the use of whole exome sequencing (WES) as a tool for identification of RD mutations, with the aim of assessing its applicability in a diagnostic context. METHODOLOGY/PRINCIPAL FINDINGS: We ascertained 12 Spanish families with seemingly recessive RD. All of the index patients underwent mutational pre-screening by chip-based sequence hybridization and resulted to be negative for known RD mutations. With the exception of one pedigree, to simulate a standard diagnostic scenario we processed by WES only the DNA from the index patient of each family, followed by in silico data analysis. We successfully identified causative mutations in patients from 10 different families, which were later verified by Sanger sequencing and co-segregation analyses. Specifically, we detected pathogenic DNA variants (∼50% novel mutations) in the genes RP1, USH2A, CNGB3, NMNAT1, CHM, and ABCA4, responsible for retinitis pigmentosa, Usher syndrome, achromatopsia, Leber congenital amaurosis, choroideremia, or recessive Stargardt/cone-rod dystrophy cases. CONCLUSIONS/SIGNIFICANCE: Despite the absence of genetic information from other family members that could help excluding nonpathogenic DNA variants, we could detect causative mutations in a variety of genes known to represent a wide spectrum of clinical phenotypes in 83% of the patients analyzed. Considering the constant drop in costs for human exome sequencing and the relative simplicity of the analyses made, this technique could represent a valuable tool for molecular diagnostics or genetic research, even in cases for which no genotypes from family members are available.

Sequencing and characterizing the genome of Estrella lausannensis as an undergraduate project: training students and biological insights.

With the widespread availability of high-throughput sequencing technologies, sequencing projects have become pervasive in the molecular life sciences. The huge bulk of data generated daily must be analyzed further by biologists with skills in bioinformatics and by "embedded bioinformaticians," i.e., bioinformaticians integrated in wet lab research groups. Thus, students interested in molecular life sciences must be trained in the main steps of genomics: sequencing, assembly, annotation and analysis. To reach that goal, a practical course has been set up for master students at the University of Lausanne: the "Sequence a genome" class. At the beginning of the academic year, a few bacterial species whose genome is unknown are provided to the students, who sequence and assemble the genome(s) and perform manual annotation. Here, we report the progress of the first class from September 2010 to June 2011 and the results obtained by seven master students who specifically assembled and annotated the genome of Estrella lausannensis, an obligate intracellular bacterium related to Chlamydia. The draft genome of Estrella is composed of 29 scaffolds encompassing 2,819,825 bp that encode for 2233 putative proteins. Estrella also possesses a 9136 bp plasmid that encodes for 14 genes, among which we found an integrase and a toxin/antitoxin module. Like all other members of the Chlamydiales order, Estrella possesses a highly conserved type III secretion system, considered as a key virulence factor. The annotation of the Estrella genome also allowed the characterization of the metabolic abilities of this strictly intracellular bacterium. Altogether, the students provided the scientific community with the Estrella genome sequence and a preliminary understanding of the biology of this recently-discovered bacterial genus, while learning to use cutting-edge technologies for sequencing and to perform bioinformatics analyses.

The early Miocene rotation of the Corso-Sardinian block. New paleomagnetic constraints for the end of the motion

The paleomagnetic investigations carried out in the 70's on Oligo-Miocene volcanics of Sardinia have demonstrated that the island was turned by 35-30 degrees clockwise from 33 Ma up to 3-1-20.5 Ma and rotated counterclockwise in a few million years [De Jong et al., 1969, 1973; Bobier et Coulon, 1970; Coulon et al., 1974; Manzoni, 1974, 1975; Bellon rr nl.. 1977: Edel et Lortscher, 1977; Edel, 1979, 1980]. Since then, the end of the rotation fixed at 19 Ma by Montigny er al. [1981] was the subject of discussions and several studies associating paleomagnetism and radiometric dating were undertaken [Assorgia er al., 1994: Vigliotti et Langenheim, 1995: Deino et al., 1997; Gattacceca rt Deino, 1999]. This is a contribution to this debate that is hampered by thr important secular variation recorded in the volcanics. The only way to get our of this problem is to sample series of successive flows as completely as possible, and to reduce the effect of secular variation by the calculation of means. Sampling was performed north of Bonorva in 5 pyroclastic flows that belong to the upper ignimbritic series SI2 according to Coulon rr nl. [1974] or LBLS according to Assorgia et al. [1997] (fig. I). Ar-40/Ar-39 dating of biotites from the debris flow (MDF) has yielded an age or 18.35 +/- 0.03 Ma [Dubois, 2000]. Five of the investigated sites are located beneath the debris flow ITV, TVB, TVD, SPM85, SPM86), one site was cured in the matrix of the debris flow (MDF) and one in 4 metric blocks included in the flow (DFC). Another site was sampled in the upper ash flow (PDM) that marks the end of the pyroclastic activity, just before the marine transgression. According to micropaleontological and radiometric dating this transgression has occurred between 18.35 and 17.6 Ma [Dubois, 2000]. After removal of a soft viscous component, the thermal demagnetization generally shows a univectorial behaviour of the remanent magnetization (fig. 2a). The maximum unblocking temperatures of 580-620 degrees (tab. I) and a rapid saturation below 100 mT (fig. 3) indicate that the carrier of the characteristic magnetization is magnetite. The exception comes: from the upper site PDM in which were found two characteristic components, one with a normal polarity and low unblocking temperatures up to 350 degreesC and one with a reversed polarity and maximum unblocking temperatures at 580-600 degreesC of magnetite. After calculation of a mean direction for each flow, the mean << Al >> direction 4 degrees /57 degrees (alpha (95) = 13 degrees) computed with the mean directions for the 5 flows may be considered as weakly affected by secular variation. But the results require a more careful examination. The declinations are N to NNW beneath the debris flow. NNW in the debris flow. and NNE (or SSW) above the debris flow, The elongated distribution of the directions obtained at sites TVB and TVD. scattered from the mean direction of TV to the mean direction of MDF is interpreted as due to partial overprinting during the debris How volcanic episode, The low temperature component PDMa is likely related to the alteration seen on thin sections and is also viewed as an overprint. As NNE/SSW directions occur as well below (mean direction << B >> : 5 degrees /58 degrees) as above the debris flow (PDMb : 200 degrees/-58 degrees). the NNW directions (<< C >> : 337 degrees /64 degrees) associated with the debris flow volcanism may be interpreted as resulting from a magnetic field excursion. According to the polarity scale of Cande and Kent [1992, 1995] and the radiometric age of MDF, the directions with normal polarity (TV, TVB, TVD, SPM85. SPM86a. MDF. DFC) may represent the period 5En. while the directions with reversed polarity PDMb and SPM86b were likely acquired during the period 5Dr. Using the mean << Al >> direction, the mean << B >>, or the PDM direction (tab. I). the deviation in declination with the direction of stable Europe 6.4 degrees /58.7 degrees (alpha (95) = 8 degrees) for a selection of 4 middle Tertiary poles by Besse et Courtillot [1991] or 7 degrees /56 degrees (alpha (95) = 3 degrees) for 19 poles listed by Edel [1980] can be considered as negligible. Using the results from the uppermost ignimbritic layer of Anglona also emplaced around 18.3 Ma [Odin rt al.. 1994]. the mean direction << E >> (3 degrees /51.5 degrees) leads to the same conclusion. On the contrary, when taking into account all dated results available for the period 5En (mean direction << D >> 353 degrees /56 degrees for 45 sites) (tab. II). the deviation 13 degrees is much more significant. As the rotation of Sardinia started around 21-20.5 Ma. the assumption of a constant velocity of rotation and the deviations of the Sardinia directions with respect to the stable Europe direction locate the end of the motion between 18.3 and 17.2 or 16.7 Ma (fig. 4). During the interval 18.35-17.5 Ma, the marine transgression took place. At the same period a NE-SW shortening interpreted as resulting from the collision of Sardinia with Apulia affected different parts of the island [Letouzey et al., 1982]. Consequently, the new paleomagnetic results and the tectono-sedimentary evolution are in favour of an end of the rotation at 17.5-18 Ma.

Whole genome sequencing in patients with retinitis pigmentosa reveals pathogenic DNA structural changes and NEK2 as a new disease gene.

We performed whole genome sequencing in 16 unrelated patients with autosomal recessive retinitis pigmentosa (ARRP), a disease characterized by progressive retinal degeneration and caused by mutations in over 50 genes, in search of pathogenic DNA variants. Eight patients were from North America, whereas eight were Japanese, a population for which ARRP seems to have different genetic drivers. Using a specific workflow, we assessed both the coding and noncoding regions of the human genome, including the evaluation of highly polymorphic SNPs, structural and copy number variations, as well as 69 control genomes sequenced by the same procedures. We detected homozygous or compound heterozygous mutations in 7 genes associated with ARRP (USH2A, RDH12, CNGB1, EYS, PDE6B, DFNB31, and CERKL) in eight patients, three Japanese and five Americans. Fourteen of the 16 mutant alleles identified were previously unknown. Among these, there was a 2.3-kb deletion in USH2A and an inverted duplication of ∼446 kb in EYS, which would have likely escaped conventional screening techniques or exome sequencing. Moreover, in another Japanese patient, we identified a homozygous frameshift (p.L206fs), absent in more than 2,500 chromosomes from ethnically matched controls, in the ciliary gene NEK2, encoding a serine/threonine-protein kinase. Inactivation of this gene in zebrafish induced retinal photoreceptor defects that were rescued by human NEK2 mRNA. In addition to identifying a previously undescribed ARRP gene, our study highlights the importance of rare structural DNA variations in Mendelian diseases and advocates the need for screening approaches that transcend the analysis of the coding sequences of the human genome.

Association study of 182 candidate genes in anorexia nervosa.

We performed association studies with 5,151 SNPs that were judged as likely candidate genetic variations conferring susceptibility to anorexia nervosa (AN) based on location under reported linkage peaks, previous results in the literature (182 candidate genes), brain expression, biological plausibility, and estrogen responsivity. We employed a case-control design that tested each SNP individually as well as haplotypes derived from these SNPs in 1,085 case individuals with AN diagnoses and 677 control individuals. We also performed separate association analyses using three increasingly restrictive case definitions for AN: all individuals with any subtype of AN (All AN: n = 1,085); individuals with AN with no binge eating behavior (AN with No Binge Eating: n = 687); and individuals with the restricting subtype of AN (Restricting AN: n = 421). After accounting for multiple comparisons, there were no statistically significant associations for any individual SNP or haplotype block with any definition of illness. These results underscore the importance of large samples to yield appropriate power to detect genotypic differences in individuals with AN and also motivate complementary approaches involving Genome-Wide Association (GWA) studies, Copy Number Variation (CNV) analyses, sequencing-based rare variant discovery assays, and pathway-based analysis in order to make up for deficiencies in traditional candidate gene approaches to AN.

Probing the T-cell receptor repertoire with deep sequencing.

PURPOSE OF REVIEW: To review major findings on the T-cell receptor (TCR) repertoire diversity in response to several viral infections based on conventional methods of PCR, cloning and sequencing and to discuss their limitations in light of the recent methodological advances in deep sequencing.¦RECENT FINDINGS: Direct sequencing of TCR expressed by Ag-specific T cells isolated ex vivo has revealed that the TCR repertoire is not as restricted as previously estimated. Furthermore, analyses performed independently of the T-cell clonal hierarchy have brought to light an unexpected diversity. The choice of methods is critical to characterize the complexity of the repertoire. Recent advances in deep sequencing have uncovered the diversity of the TCR repertoire and shown that the size of the repertoire in naive and Ag-experienced memory T cells is three-fold to 15-fold larger than formerly estimated. Interestingly, the TCR complementary determining region 3 sequences are not randomly selected and a certain degree of shared TCR repertoire has been observed between different individuals.¦SUMMARY: Deep sequencing is a major methodological advance allowing more accurate molecular characterization of the TCR repertoire. In the near future, such technologies will further contribute to delineate the complexity of pathogen-specific T-cell response and help defining correlates of a protective immunity.

Exome sequencing identifies recurrent somatic MAP2K1 and MAP2K2 mutations in melanoma.

We performed exome sequencing to detect somatic mutations in protein-coding regions in seven melanoma cell lines and donor-matched germline cells. All melanoma samples had high numbers of somatic mutations, which showed the hallmark of UV-induced DNA repair. Such a hallmark was absent in tumor sample-specific mutations in two metastases derived from the same individual. Two melanomas with non-canonical BRAF mutations harbored gain-of-function MAP2K1 and MAP2K2 (MEK1 and MEK2, respectively) mutations, resulting in constitutive ERK phosphorylation and higher resistance to MEK inhibitors. Screening a larger cohort of individuals with melanoma revealed the presence of recurring somatic MAP2K1 and MAP2K2 mutations, which occurred at an overall frequency of 8%. Furthermore, missense and nonsense somatic mutations were frequently found in three candidate melanoma genes, FAT4, LRP1B and DSC1.

IROme, a New High-Throughput Molecular Tool for the Diagnosis of Inherited Retinal Dystrophies-A Price Comparison with Sanger Sequencing.

The molecular diagnosis of retinal dystrophies (RD) is difficult because of genetic and clinical heterogeneity. Previously, the molecular screening of genes was done one by one, sometimes in a scheme based on the frequency of sequence variants and the number of exons/length of the candidate genes. Payment for these procedures was complicated and the sequential billing of several genes created endless paperwork. We therefore evaluated the costs of generating and sequencing a hybridization-based DNA library enriched for the 64 most frequently mutated genes in RD, called IROme, and compared them to the costs of amplifying and sequencing these genes by the Sanger method. The production cost generated by the high-throughput (HT) sequencing of IROme was established at CHF 2,875.75 per case. Sanger sequencing of the same exons cost CHF 69,399.02. Turnaround time of the analysis was 3 days for IROme. For Sanger sequencing, it could only be estimated, as we never sequenced all 64 genes in one single patient. Sale cost for IROme calculated on the basis of the sale cost of one exon by Sanger sequencing is CHF 8,445.88, which corresponds to the sale price of 40 exons. In conclusion, IROme is cheaper and faster than Sanger sequencing and therefore represents a sound approach for the diagnosis of RD, both scientifically and economically. As a drop in the costs of HT sequencing is anticipated, target resequencing might become the new gold standard in the molecular diagnosis of RD.