Complex chromosomal rearrangements (CCRs) are rearrangements involving more than two chromosomes or more than two breakpoints. Whole genome sequencing (WGS) allows for outstanding high resolution characterization on the nucleotide level in unique sequences of such rearrangements, but problems remain for mapping breakpoints in repetitive regions of the genome, which are known to be prone to rearrangements. Hence, multiple complementary WGS experiments are sometimes needed to solve the structures of CCRs. We have studied three individuals with CCRs: Case 1 and Case 2 presented with de novo karyotypically balanced, complex interchromosomal rearrangements (46,XX,t(2;8;15)(q35;q24.1;q22) and 46,XY,t(1;10;5)(q32;p12;q31)), and Case 3 presented with a de novo, extremely complex intrachromosomal rearrangement on chromosome 1. Molecular cytogenetic investigation revealed cryptic deletions in the breakpoints of chromosome 2 and 8 in Case 1, and on chromosome 10 in Case 2, explaining their clinical symptoms. In Case 3, 26 breakpoints were identified using WGS, disrupting five known disease genes. All rearrangements were subsequently analyzed using optical maps, linked-read WGS, and short-read WGS. In conclusion, we present a case series of three unique de novo CCRs where we by combining the results from the different technologies fully solved the structure of each rearrangement. The power in combining short-read WGS with long-molecule sequencing or optical mapping in these unique de novo CCRs in a clinical setting is demonstrated.

The current human reference sequence (GRCh38) is a foundation for large-scale sequencing projects. However, recent studies have suggested that GRCh38 may be incomplete and give a suboptimal representation of specific population groups. Here, we performed a de novo assembly of two Swedish genomes that revealed over 10 Mb of sequences absent from the human GRCh38 reference in each individual. Around 6 Mb of these novel sequences (NS) are shared with a Chinese personal genome. The NS are highly repetitive, have an elevated GC-content, and are primarily located in centromeric or telomeric regions. Up to 1 Mb of NS can be assigned to chromosome Y, and large segments are also missing from GRCh38 at chromosomes 14, 17, and 21. Inclusion of NS into the GRCh38 reference radically improves the alignment and variant calling from short-read whole-genome sequencing data at several genomic loci. A re-analysis of a Swedish population-scale sequencing project yields > 75,000 putative novel single nucleotide variants (SNVs) and removes > 10,000 false positive SNV calls per individual, some of which are located in protein coding regions. Our results highlight that the GRCh38 reference is not yet complete and demonstrate that personal genome assemblies from local populations can improve the analysis of short-read whole-genome sequencing data.

Clustered copy number variants (CNVs) as detected by chromosomal microarray analysis (CMA) are often reported as germline chromothripsis. However, such cases might need further investigations by massive parallel whole genome sequencing (WGS) in order to accurately define the underlying complex rearrangement, predict the occurrence mechanisms and identify additional complexities. Here, we utilized WGS to delineate the rearrangement structure of 21 clustered CNV carriers first investigated by CMA and identified a total of 83 breakpoint junctions (BPJs). The rearrangements were further sub-classified depending on the patterns observed: I) Cases with only deletions (n = 8) often had additional structural rearrangements, such as insertions and inversions typical to chromothripsis; II) cases with only duplications (n = 7) or III) combinations of deletions and duplications (n = 6) demonstrated mostly interspersed duplications and BPJs enriched with microhomology. In two cases the rearrangement mutational signatures indicated both a breakage-fusion-bridge cycle process and haltered formation of a ring chromosome. Finally, we observed two cases with Alu- and LINE-mediated rearrangements as well as two unrelated individuals with seemingly identical clustered CNVs on 2p25.3, possibly a rare European founder rearrangement. In conclusion, through detailed characterization of the derivative chromosomes we show that multiple mechanisms are likely involved in the formation of clustered CNVs and add further evidence for chromoanagenesis mechanisms in both "simple" and highly complex chromosomal rearrangements. Finally, WGS characterization adds positional information, important for a correct clinical interpretation and deciphering mechanisms involved in the formation of these rearrangements.

Highly penetrant cancer syndromes account for less than 5% of all cases with familial colorectal cancer (CRC), and other genetic contribution explains the majority of the genetic contribution to CRC. A CRC susceptibility locus on chromosome 9q has been suggested. In this study, families where risk of CRC was linked to the region, were used to search for predisposing mutations in all genes in the region. No disease-causing mutation was found. Next, haplotype association studies were performed in the region, comparing Swedish CRC cases (2664) and controls (4782). Two overlapping haplotypes were suggested. One 10-SNP haplotype was indicated in familial CRC (OR 1.4, p = 0.00005) and one 25-SNP haplotype was indicated in sporadic CRC (OR 2.2, p = 0.0000012). The allele frequencies of the 10-SNP and the 25-SNP haplotypes were 13.7% and 2.5% respectively and both included one RNA, RP11-332M4.1 and RP11-l80l4.2, in the non-overlapping regions. The sporadic 25-SNP haplotype could not be studied further, but the familial 10-SNP haplotype was analyzed in 61 additional CRC families, and 6 of them were informative for all markers and had the risk haplotype. Targeted sequencing of the 10-SNP region in the linked families identified one variant in RP11-332M4.1, suggestive to confer the increased CRC risk on this haplotype. Our results support the presence of two loci at 9q22.32, each with one RNA as the putative cause of increased CRC risk. These RNAs could exert their effect through the same, or different, genes/pathways, possibly through the regulation of neighboring genes, such as PTCH1, FANCC, DKFZP434H0512, ERCC6L2 or the processed transcript LINC00046.

In Kinetoplastids, protein-coding genes are transcribed polycistronically by RNA polymerase II. Individual mature mRNAs are generated from polycistronic precursors by 5' trans splicing of a 39-nt capped leader RNA and 3' polyadenylation. It was previously known that trans splicing generally occurs at an AG dinucleotide downstream of a polypyrimidine tract, and that polyadenylation is coupled to downstream trans splicing. The few polyadenylation sites that had been examined were 100-400 nt upstream of the polypyrimidine tract which marked the adjacent trans splice site. We wished to define the sequence requirements for trypanosome mRNA processing more tightly and to generate a predictive algorithm. By scanning all available Trypanosoma brucei cDNAs for splicing and polyadenylation sites, we found that trans splicing generally occurs at the first AG following a polypyrimidine tract of 8-25 nt, giving rise to 5'-UTRs of a median length of 68 nt. We also found that in general, polyadenylation occurs at a position with one or more A residues located between 80 and 140 nt from the downstream polypyrimidine tract. These data were used to calibrate free parameters in a grammar model with distance constraints, enabling prediction of polyadenylation and trans splice sites for most protein-coding genes in the trypanosome genome. The data from the genome analysis and the program are available from: .

Temporal lobe resection for epilepsy involves a risk of damaging the anterior part of the optic radiation, Meyer's loop, causing a contralateral upper quadrant visual field defect. This study aims to assess the intersubject variability in the course of Meyer's loop in vivo by diffusion tensor imaging and tractography.

Trypanosomatids are widespread parasites that cause three major tropical diseases. In trypanosomatids, as in most other organisms, acetylation is a common protein modification that is important in multiple, diverse processes. This paper describes a new member of the Trypanosoma cruzi acetyltransferase family. The gene is single copy and orthologs are also present in the other two sequenced trypanosomatids, Trypanosoma brucei and Leishmania major. This protein (TcAT-1) has the essential motifs present in members of the GCN5-related acetyltransferase (GNAT) family, as well as an additional motif also found in some enzymes from plant and animal species. The protein is evolutionarily more closely related to this group of enzymes than to histone acetyltransferases. The native protein has a cytosolic cellular location and is present in all three life-cycle stages of the parasite. The recombinant protein was shown to have autoacetylation enzymatic activity.

Well-ordered HIV-1 envelope glycoprotein (Env) trimers are prioritized for clinical evaluation, and there is a need for an improved understanding about how elicited B cell responses evolve following immunization. To accomplish this, we prime-boosted rhesus macaques with clade C NFL trimers and identified 180 unique Ab lineages from ∼1,000 single-sorted Env-specific memory B cells. We traced all lineages in high-throughput heavy chain (HC) repertoire (Rep-seq) data generated from multiple immune compartments and time points and expressed several as monoclonal Abs (mAbs). Our results revealed broad dissemination and high levels of somatic hypermutation (SHM) of most lineages, including tier 2 virus neutralizing lineages, following boosting. SHM was highest in the Ab complementarity determining regions (CDRs) but also surprisingly high in the framework regions (FRs), especially FR3. Our results demonstrate the capacity of the immune system to affinity-mature large numbers of Env-specific B cell lineages simultaneously, supporting the use of regimens consisting of repeated boosts to improve each Ab, even those belonging to less expanded lineages.

Exome and genome sequencing is becoming the method of choice for rare disease diagnostics. One of the key challenges remaining is distinguishing the disease causing variants from the benign background variation. After analysis and annotation of the sequencing data there are typically thousands of candidate variants requiring further investigation. One of the most effective and least biased ways to reduce this number is to assess the rarity of a variant in any population. Currently, there are a number of reliable sources of information for major population frequencies when considering single nucleotide variants (SNVs) and small insertion and deletions (INDELs), with gnomAD as the most prominent public resource available. However, local variation or frequencies in sub-populations may be underrepresented in these public resources. In contrast, for structural variation (SV), the background frequency in the general population is more or less unknown mostly due to challenges in calling SVs in a consistent way. Keeping track of local variation is one way to overcome these problems and significantly reduce the number of potential disease causing variants retained for manual inspection, both for SNVs and SVs.

Analysis of cell-free tumour DNA, a liquid biopsy, is a promising biomarker for cancer. We have performed a proof-of principle study to test the applicability in the clinical setting, analysing copy number alterations (CNAs) in plasma and tumour tissue from 44 patients with gastro-oesophageal cancer.

Ataxia channelopathies share common features such as slow motor progression and variable degrees of cognitive dysfunction. Mutations in potassium voltage-gated channel subfamily D member 3 (KCND3), encoding the K+ channel, Kv4.3, are associated with spinocerebellar ataxia (SCA) 19, allelic with SCA22. Mutations in potassium voltage-gated channel subfamily C member 3 (KCNC3), encoding another K+ channel, Kv3.3, cause SCA13. First, a comprehensive phenotype assessment was carried out in a family with autosomal dominant ataxia harboring 2 genetic variants in KCNC3 and KCND3. To evaluate the physiological impact of these variants on channel currents, in vitro studies were performed.

Recent advances in sequencing technology have created unprecedented opportunities for biological research. However, the increasing throughput of these technologies has created many challenges for data management and analysis. As the demand for sophisticated analyses increases, the development time of software and algorithms is outpacing the speed of traditional publication. As technologies continue to be developed, methods change rapidly, making publications less relevant for users. The SEQanswers wiki (SEQwiki) is a wiki database that is actively edited and updated by the members of the SEQanswers community (http://SEQanswers.com/). The wiki provides an extensive catalogue of tools, technologies and tutorials for high-throughput sequencing (HTS), including information about HTS service providers. It has been implemented in MediaWiki with the Semantic MediaWiki and Semantic Forms extensions to collect structured data, providing powerful navigation and reporting features. Within 2 years, the community has created pages for over 500 tools, with approximately 400 literature references and 600 web links. This collaborative effort has made SEQwiki the most comprehensive database of HTS tools anywhere on the web. The wiki includes task-focused mini-reviews of commonly used tools, and a growing collection of more than 100 HTS service providers. SEQwiki is available at: http://wiki.SEQanswers.com/.

The pathway for RNA interference is widespread in metazoans and participates in numerous cellular tasks, from gene silencing to chromatin remodeling and protection against retrotransposition. The unicellular eukaryote Trypanosoma cruzi is missing the canonical RNAi pathway and is unable to induce RNAi-related processes. To further understand alternative RNA pathways operating in this organism, we have performed deep sequencing and genome-wide analyses of a size-fractioned cDNA library (16-61 nt) from the epimastigote life stage. Deep sequencing generated 582,243 short sequences of which 91% could be aligned with the genome sequence. About 95-98% of the aligned data (depending on the haplotype) corresponded to small RNAs derived from tRNAs, rRNAs, snRNAs and snoRNAs. The largest class consisted of tRNA-derived small RNAs which primarily originated from the 3' end of tRNAs, followed by small RNAs derived from rRNA. The remaining sequences revealed the presence of 92 novel transcribed loci, of which 79 did not show homology to known RNA classes.

There is tremendous potential for genome sequencing to improve clinical diagnosis and care once it becomes routinely accessible, but this will require formalizing research methods into clinical best practices in the areas of sequence data generation, analysis, interpretation and reporting. The CLARITY Challenge was designed to spur convergence in methods for diagnosing genetic disease starting from clinical case history and genome sequencing data. DNA samples were obtained from three families with heritable genetic disorders and genomic sequence data were donated by sequencing platform vendors. The challenge was to analyze and interpret these data with the goals of identifying disease-causing variants and reporting the findings in a clinically useful format. Participating contestant groups were solicited broadly, and an independent panel of judges evaluated their performance.

Familial hemophagocytic lymphohistiocytosis (FHL) is an often-fatal hyperinflammatory disorder caused by autosomal recessive mutations in PRF1, UNC13D, STX11, and STXBP2. We identified a homozygous STX11 mutation, c.173T > C (p.L58P), in three patients presenting clinically with hemophagocytic lymphohistiocytosis from unrelated Pakistani families. The mutation yields an amino acid substitution in the N-terminal Habc domain of syntaxin-11 and resulted in defective natural killer cell degranulation. Notably, syntaxin-11 expression was decreased in patient cells. However, in an ectopic expression system, syntaxin-11 L58P was expressed at levels comparable to wild-type syntaxin-11, but did not bind Munc18-2. Moreover, another N-terminal syntaxin-11 mutant, R4A, also did not bind Munc18-2. Thus, we have identified a novel missense STX11 mutation causative of FHL type 4. The syntaxin-11 R4A and L58P mutations reveal that both the N-terminus and Habc domain of syntaxin-11 are required for binding to Munc18-2, implying similarity to the dynamic binary binding of neuronal syntaxin-1 to Munc18-1.

The recently launched veterinary anthelmintic drench for sheep (Novartis Animal Health Inc., Switzerland) containing the nematocide monepantel represents a new class of anthelmintics: the amino-acetonitrile derivatives (AADs), much needed in view of widespread resistance to the classical drugs. Recently, it was shown that the ACR-23 protein in Caenorhabditis elegans and a homologous protein, MPTL-1 in Haemonchus contortus, are potential targets for AAD action. Both proteins belong to the DEG-3 subfamily of acetylcholine receptors, which are thought to be nematode-specific, and different from those targeted by the imidazothiazoles (e.g. levamisole). Here we provide further evidence that Cel-ACR-23 and Hco-MPTL-1-like subunits are involved in the monepantel-sensitive phenotype. We performed comparative genomics of ligand-gated ion channel genes from several nematodes and subsequently assessed their sensitivity to anthelmintics. The nematode species in the Caenorhabditis genus, equipped with ACR-23/MPTL-1-like receptor subunits, are sensitive to monepantel (EC(50)<1.25 µM), whereas the related nematodes Pristionchus pacificus and Strongyloides ratti, which lack an ACR-23/MPTL-1 homolog, are insensitive (EC(50)>43 µM). Genome sequence information has long been used to identify putative targets for therapeutic intervention. We show how comparative genomics can be applied to predict drug sensitivity when molecular targets of a compound are known or suspected.

Peroxisome proliferator-activated receptors (PPARs) are key transcription factors that regulate adipose development and function, and the conversion of white into brown-like adipocytes. Here we investigated whether PPARα and PPARγ activation synergize to induce the browning of white fat.

The detection of recurrent somatic chromosomal rearrangements is standard of care for most leukemia types. Even though karyotype analysis-a low-resolution genome-wide chromosome analysis-is still the gold standard, it often needs to be complemented with other methods to increase resolution. To evaluate the feasibility and applicability of mate pair whole genome sequencing (MP-WGS) to detect structural chromosomal rearrangements in the diagnostic setting, we sequenced ten bone marrow samples from leukemia patients with recurrent rearrangements. Samples were selected based on cytogenetic and FISH results at leukemia diagnosis to include common rearrangements of prognostic relevance. Using MP-WGS and in-house bioinformatic analysis all sought rearrangements were successfully detected. In addition, unexpected complexity or additional, previously undetected rearrangements was unraveled in three samples. Finally, the MP-WGS analysis pinpointed the location of chromosome junctions at high resolution and we were able to identify the exact exons involved in the resulting fusion genes in all samples and the specific junction at the nucleotide level in half of the samples. The results show that our approach combines the screening character from karyotype analysis with the specificity and resolution of cytogenetic and molecular methods. As a result of the straightforward analysis and high-resolution detection of clinically relevant rearrangements, we conclude that MP-WGS is a feasible method for routine leukemia diagnostics of structural chromosomal rearrangements.

Duplicative horizontal gene transfer may bring two previously separated homologous genes together, which may raise questions about the interplay between the gene products. One such gene pair is the "native" PgiC1 and "foreign" PgiC2 in the perennial grass Festuca ovina. Both PgiC1 and PgiC2 encode cytosolic phosphoglucose isomerase, a dimeric enzyme whose proper binding is functionally essential. Here, we use biophysical simulations to explore the inter-monomer binding of the two homodimers and the heterodimer that can be produced by PgiC1 and PgiC2 in F. ovina. Using simulated native-state ensembles, we examine the structural properties and binding tightness of the dimers. In addition, we investigate their ability to withstand dissociation when pulled by a force. Our results suggest that the inter-monomer binding is tighter in the PgiC2 than the PgiC1 homodimer, which could explain the more frequent occurrence of the foreign PgiC2 homodimer in dry habitats. We further find that the PgiC1 and PgiC2 monomers are compatible with heterodimer formation; the computed binding tightness is comparable to that of the PgiC1 homodimer. Enhanced homodimer stability and capability of heterodimer formation with PgiC1 are properties of PgiC2 that may contribute to the retaining of the otherwise redundant PgiC2 gene.

Spastic ataxias are rare neurogenetic disorders involving spinocerebellar and pyramidal tracts. Many genes are involved. Among them, CAPN1, when mutated, is responsible for a complex inherited form of spastic paraplegia (SPG76). We report the largest published series of 21 novel patients with nine new CAPN1 disease-causing variants and their clinical characteristics from two European university hospitals (Paris and Stockholm). After a formal clinical examination, causative variants were identified by next-generation sequencing and confirmed by Sanger sequencing. CAPN1 variants are a rare cause (~ 1.4%) of young-adult-onset spastic ataxia; however, together with all published cases, they allowed us to better describe the clinical and genetic spectra of this form. Truncating variants are the most frequent, and missense variants lead to earlier age at onset in favor of an additional deleterious effect. Cerebellar ataxia with cerebellar atrophy, dysarthria and lower limb weakness are often associated with spasticity. We also suggest that cognitive impairment and depression should be assessed specifically in the follow-up of SPG76 cases.