Complex balanced intrachromosomal rearrangement involving PITX2 identified as a cause of Axenfeld-Rieger Syndrome.

Axenfeld-Rieger Syndrome (ARS) type 1 is a rare autosomal dominant condition characterized by anterior chamber anomalies, umbilical defects, dental hypoplasia, and craniofacial anomalies, with Meckel's diverticulum in some individuals. Here, we describe a clinically ascertained female of childbearing age with ARS for whom clinical targeted sequencing and deletion/duplication analysis followed by clinical exome and genome sequencing resulted in no pathogenic variants or variants of unknown significance in PITX2 or FOXC1. Advanced bioinformatic analysis of the genome data identified a complex, balanced rearrangement disrupting PITX2. This case is the first reported intrachromosomal rearrangement leading to ARS, illustrating that for patients with compelling clinical phenotypes but negative genomic testing, additional bioinformatic analysis are essential to identify subtle genomic abnormalities in target genes.

Characterization of unusual iAMP21 B-lymphoblastic leukemia (iAMP21-ALL) from the Mayo Clinic and Children's Oncology Group.

Acute lymphoblastic leukemia (B-ALL) with intrachromosomal amplification of chromosome 21 (iAMP21-ALL) represents a recurrent high-risk cytogenetic abnormality and accurate identification is critical for appropriate clinical management. Identification of iAMP21-ALL has historically relied on fluorescence in situ hybridization (FISH) using a RUNX1 probe. Current classification requires ≥ five copies of RUNX1 per cell and ≥ three additional copies of RUNX1 on a single abnormal iAMP21-chromosome. We sought to evaluate the performance of the RUNX1 probe in the identification of iAMP21-ALL. This study was a retrospective evaluation of iAMP21-ALL in the Mayo Clinic and Children's Oncology Group cohorts. Of 207 cases of iAMP21-ALL, 188 (91%) were classified as "typical" iAMP21-ALL, while 19 (9%) cases were classified as "unusual" iAMP21-ALL. The "unusual" iAMP21 cases did not meet the current definition of iAMP21 by FISH but were confirmed to have iAMP21 by chromosomal microarray. Half of the "unusual" iAMP21-ALL cases had less than five RUNX1 signals, while the remainder had ≥ five RUNX1 signals with some located apart from the abnormal iAMP21-chromosome. Nine percent of iAMP21-ALL cases fail to meet the FISH definition of iAMP21-ALL demonstrating that laboratories are at risk of misidentification of iAMP21-ALL when relying only on the RUNX1 FISH probe. Incorporation of chromosomal microarray testing circumvents these risks.

Identification of EWSR1 rearrangements in patients with immature hematopoietic neoplasms: A case series and review of literature.

Rearrangement of the EWSR1 gene (22q12.2) is a well-recognized genetic lesion in bone and soft tissue tumors. However, few reports have suggested that EWSR1 rearrangements may also occur in the setting of hematopoietic tumors. We herein describe two cases of immature hematopoietic neoplasms presenting with EWSR1 rearrangements. The first occurred in a 41-year-old female diagnosed with mixed-phenotype acute leukemia, B/T/myeloid, in which conventional chromosome analysis revealed a t(2;22)(q35;q12). Further analysis with whole genome sequencing revealed that this rearrangement led to an EWSR1::FEV gene fusion. The second case was identified in an 18-year-old male with a high-grade B-cell lineage malignant neoplasm with immature features in which conventional chromosome analysis revealed a t(17;22)(q25;q12). Mate-pair sequencing, a next generation sequencing-based assay, was performed and revealed three in-frame chimeric gene fusions involving the EWSR1, TEF and STRADA gene regions. This report further expands the repertoire of hematopoietic neoplasms with EWSR1 fusions and partner genes involved in these rearrangements.

Increased complexity of t(11;14) rearrangements in plasma cell neoplasms compared with mantle cell lymphoma.

Plasma cell neoplasms (PCN) and mantle cell lymphoma (MCL) can both harbor t(11;14)(q13;q32) (CCND1/IGH), usually resulting in cyclin D1 overexpression. In some cases, particularly at low levels of disease, it can be morphologically challenging to distinguish between these entities in the bone marrow (BM) since PCN with t(11;14) are often CD20-positive with lymphoplasmacytic cytology, while MCL can rarely have plasmacytic differentiation. We compared the difference in CCND1/IGH by fluorescence in situ hybridization (FISH) in PCN and MCL to evaluate for possible differentiating characteristics. We identified 326 cases of MCL with t(11;14) and 279 cases of PCN with t(11;14) from either formalin-fixed, paraffin-embedded tissue or fresh BM specimens. The "typical," balanced CCND1/IGH FISH signal pattern was defined as three total CCND1 signals, three total IGH signals, and two total fusion signals. Any deviation from the "typical" pattern was defined as an "atypical" pattern, which was further stratified into "gain of fusion" vs "complex" patterns. There was a significantly higher proportion of cases that showed an atypical FISH pattern in PCN compared with MCL (53% vs 27%, P < .0001). There was also a significantly higher proportion of cases that showed a complex FISH pattern in PCN compared with MCL (47% vs 17%, P < .0001). We confirmed these findings using mate-pair sequencing of 25 PCN and MCL samples. PCN more often have a complex CCND1/IGH FISH pattern compared with MCL, suggesting possible differences in the genomic mechanisms underlying these rearrangements in plasma cells compared with B cells.

Identification of a novel KMT2A/GIMAP8 gene fusion in a pediatric patient with acute undifferentiated leukemia.

Acute undifferentiated leukemia (AUL) is a very rare hematologic neoplasm that expresses no markers specific for either myeloid or lymphoid lineages. While commonly observed in several acute leukemias, KMT2A rearrangements in AUL have been rarely reported in the literature. We report the third case to our knowledge of AUL harboring a KMT2A rearrangement. Furthermore, the KMT2A/GIMAP8 gene fusion identified in this case represents a novel KMT2A rearrangement.

Theragnostic chromosomal rearrangements in treatment-naive pancreatic ductal adenocarcinomas obtained via endoscopic ultrasound.

A crucial mutational mechanism in malignancy is structural variation, in which chromosomal rearrangements alter gene functions that drive cancer progression. Herein, the presence and pattern of structural variations were investigated in twelve prospectively acquired treatment-naïve pancreatic cancers specimens obtained via endoscopic ultrasound (EUS). In many patients, this diagnostic biopsy procedure and specimen is the only opportunity to identify somatic clinically relevant actionable alterations that may impact their care and outcome. Specialized mate pair sequencing (MPseq) provided genome-wide structural variance analysis (SVA) with a view to identifying prognostic markers and possible therapeutic targets. MPseq was successfully performed on all specimens, identifying highly rearranged genomes with complete SVA on all specimens with > 20% tumour content. SVA identified chimeric fusion proteins and potentially immunogenic readthrough transcripts, change of function truncations, gains and losses of key genes linked to tumour progression. Complex localized rearrangements, termed chromoanagenesis, with broad pattern heterogeneity were observed in 10 (83%) specimens, impacting multiple genes with diverse cellular functions that could influence theragnostic evaluation and responsiveness to immunotherapy regimens. This study indicates that genome-wide MPseq can be successfully performed on very limited clinically EUS obtained specimens for chromosomal rearrangement detection and potential theragnostic targets.

Clinical utility of next generation sequencing to detect IGH/IL3 rearrangements [t(5;14)(q31.1;q32.1)] in B-lymphoblastic leukemia/lymphoma.

The t(5;14)(q31.1;q32.1) associated with B-lymphoblastic leukemia/lymphoma (B-ALL/LBL) is a rare, recurrent genetic abnormality recognized as a distinct entity by the 2017 World Health Organization (WHO) classification. In these cases, the IGH enhancer region (14q32.1) is juxtaposed to the vicinity of the IL3 gene (5q31.1), resulting in increased production of interleukin-3 (IL3) and subsequently a characteristic reactive eosinophilia. B-ALL with t(5;14)(q31.1;q32.1) may have a low lymphoblast count that can complicate detection of t(5;14)(q31.1;q32.1) by conventional chromosome studies. We have identified four patients with IGH/IL3 rearrangements despite normal conventional chromosome studies in each case [one patient had a non-clonal t(5;14)(q31;q32) finding]. Fluorescence in situ hybridization utilizing a laboratory-developed IGH break-apart probe set identified IGH rearrangements in three of four cases, and a next generation sequencing (NGS) based assay, mate-pair sequencing (MPseq), was required to characterize the IGH/IL3 rearrangements in each case. Three patients demonstrated a balanced t(5;14)(q31.1;q32.1) while one patient had a cryptic insertion of the IL3 gene into the IGH region. These results demonstrate that NGS-based assays, such as MPseq, confer an advantage in the detection of IGH/IL3 rearrangements that are otherwise challenging to characterize by traditional cytogenetic methodologies.

Optimizing clinical cytology touch preparations for next generation sequencing.

Intraoperative diagnosis is routinely performed on cytology touch preparations (TPs) from core needle biopsies (CNBs). Current interest promotes their utility as an important source of patient tissue for clinical genomic testing. Herein we present whole genome structural variant analysis (SVA) from mate-pair sequencing (MPseq) and whole exome sequencing (WES) mutation calling in DNA directly whole genome amplified (WGA) from TPs. Chromosomal copy changes and somatic DNA junction detection from MPseq of TPs were highly consistent with associated CNBs and bulk resected tissues in all cases. While increased frequency coverage noise from limitations of amplification of limited sample input was significant, this was effectively compensated by natural tumor enrichment during the TP process, which also enhanced variant detection and loss of heterozygosity evaluations from WES. This novel TP methodology enables expanded utility of frequently limited CNB for both clinical and research genomic testing.

Cryptic and atypical KMT2A-USP2 and KMT2A-USP8 rearrangements identified by mate pair sequencing in infant and childhood leukemia.

Infant leukemias are a rare group of neoplasms that are clinically and biologically distinct from their pediatric and adult counterparts. Unlike leukemia in older children where survival rates are generally favorable, infants with leukemia have a 5-year event-free survival rate of <50%. The majority of infant leukemias are characterized by KMT2A (MLL) rearrangements (~70 to 80% in acute lymphoblastic leukemia), which appear to be drivers of early leukemogenesis. In this report, we describe three cases: a 9-month-old female infant with B-acute lymphoblastic leukemia (B-ALL), an 8-month-old female presenting with B/myeloid mixed phenotype acute leukemia (MPAL), and a 16-month-old male with B-ALL. The first case had a normal karyotype and B-ALL FISH results consistent with an atypical KMT2A rearrangement. The second case had trisomy 10 as the sole chromosomal abnormality and a normal KMT2A FISH result. Case 3 had trisomy 8 and a t(11;15)(q23;q21), an atypical KMT2A rearrangement by FISH studies, and a focal deletion of 15q with a breakpoint within the USP8 gene by chromosomal microarray. Mate pair sequencing was performed on all three cases and identified a KMT2A-USP2 rearrangement (cases 1 and 2) or a KMT2A-USP8 rearrangement (case 3). These recently characterized KMT2A fusions have been described exclusively in infant and pediatric leukemia cases where the incidence varies vary according to leukemia subtype, are considered high-risk, with a high incidence of central nervous system involvement, poor response to initial prednisone treatment, and poor event free survival. Additionally, approximately half of cases are unable to be resolved using standard cytogenetic approaches and are likely under recognized. Therefore, targeted molecular approaches are suggested in genetically unresolved infant leukemia cases to characterize these prognostically relevant clones.

Detection of cryptic CCND1 rearrangements in mantle cell lymphoma by next generation sequencing.

The accurate detection of recurrent genetic abnormalities for most hematologic neoplasms is critical for diagnosis, prognosis and/or treatment. Rearrangements involving CCND1 are observed in a subset of mature B-cell neoplasms and can be reliably detected by fluorescence in situ hybridization (FISH) in most cases. However, cryptic and complex chromosomal rearrangements may pose a technical challenge for accurate diagnosis. Herein, we describe two patients with suspected mantle cell lymphoma that lacked obvious CCND1 rearrangements by FISH studies. A next generation sequencing (NGS) based assay, mate-pair sequencing (MPseq), was utilized in each case to investigate potential cryptic CCND1 rearrangements and revealed cryptic insertional events resulting in CCND1/IGH and CCND1/IGK rearrangements. These cases demonstrate that NGS-based assays, including MPseq, are a powerful approach to identify cryptic rearrangements of clinical importance that are not detected by current clinical genomics evaluation.

Acute leukemias harboring KMT2A/MLLT10 fusion: a 10-year experience from a single genomics laboratory.

The MLLT10 (formerly AF10) gene is the fourth most common KMT2A fusion partner across all acute leukemias and requires at least 3 breaks to form an in-frame KMT2A/MLLT10 fusion due to the opposite orientation of each gene. A 10-year retrospective review was performed to identify individuals from all age groups that harbor KMT2A/MLLT10 fusion obtained by our KMT2A/MLLT10 dual-color dual-fusion fluorescence in situ hybridization (D-FISH) assay. Of the 60 unique individuals identified, 31 were male and 29 were female (M:F ratio, 1.1:1) with ages ranging from 3 days to 86 years (mean 21.5 years, median 5.5 years). The diagnoses included acute myeloid leukemia (AML) (49 patients, 82%), B- or T-lymphoblastic leukemia/lymphoma (7 patients, 12%), myeloid sarcoma (3 patients, 5%), and a single case (2%) of undifferentiated leukemia. Twenty-seven of 49 patients (55%) with AML were in the infant or pediatric age group. Fifty-three of 60 patients (88%) had KMT2A/MLLT10 D-FISH signal patterns mostly consisting of single fusions. In addition, 10 (26%) of 38 patients with conventional chromosome studies had "normal" (5 patients) or abnormal (5 patients) chromosome studies that lacked structural or numeric abnormalities involving chromosomes 10 or 11, implying cryptic cytogenetic mechanisms for KMT2A/MLLT10 fusion. Lastly, mate-pair sequencing was performed on 4 AML cases, 2 of which had "normal" chromosome studies and cryptic KMT2A/MLLT10 fusion as detected by KMT2A/MLLT10 D-FISH studies, and verified the multiple breaks required to generate KMT2A/MLLT10 fusion.

Genome U-Plot: a whole genome visualization.

Motivation: The ability to produce and analyze whole genome sequencing (WGS) data from samples with structural variations (SV) generated the need to visualize such abnormalities in simplified plots. Conventional two-dimensional representations of WGS data frequently use either circular or linear layouts. There are several diverse advantages regarding both these representations, but their major disadvantage is that they do not use the two-dimensional space very efficiently. We propose a layout, termed the Genome U-Plot, which spreads the chromosomes on a two-dimensional surface and essentially quadruples the spatial resolution. We present the Genome U-Plot for producing clear and intuitive graphs that allows researchers to generate novel insights and hypotheses by visualizing SVs such as deletions, amplifications, and chromoanagenesis events. The main features of the Genome U-Plot are its layered layout, its high spatial resolution and its improved aesthetic qualities. We compare conventional visualization schemas with the Genome U-Plot using visualization metrics such as number of line crossings and crossing angle resolution measures. Based on our metrics, we improve the readability of the resulting graph by at least 2-fold, making apparent important features and making it easy to identify important genomic changes. Results: A whole genome visualization tool with high spatial resolution and improved aesthetic qualities. Availability and implementation: An implementation and documentation of the Genome U-Plot is publicly available at https://github.com/gaitat/GenomeUPlot. Contact: vasmatzis.george@mayo.edu. Supplementary information: Supplementary data are available at Bioinformatics online.

Mate pair sequencing improves detection of genomic abnormalities in acute myeloid leukemia.

OBJECTIVE: Acute myeloid leukemia (AML) can be subtyped based on recurrent cytogenetic and molecular genetic abnormalities with diagnostic and prognostic significance. Although cytogenetic characterization classically involves conventional chromosome and/or fluorescence in situ hybridization (FISH) assays, limitations of these techniques include poor resolution and the inability to precisely identify breakpoints. METHOD: We evaluated whether an NGS-based methodology that detects structural abnormalities and copy number changes using mate pair sequencing (MPseq) can enhance the diagnostic yield for patients with AML. RESULTS: Using 68 known abnormal and 20 karyotypically normal AML samples, each recurrent primary AML-specific abnormality previously identified in the abnormal samples was confirmed using MPseq. Importantly, in eight cases with abnormalities that could not be resolved by conventional cytogenetic studies, MPseq was utilized to molecularly define eight recurrent AML-fusion events. In addition, MPseq uncovered two cryptic abnormalities that were missed by conventional cytogenetic studies. Thus, MPseq improved the diagnostic yield in the detection of AML-specific structural rearrangements in 10/88 (11%) of cases analyzed. CONCLUSION: Utilization of MPseq represents a precise, molecular-based technique that can be used as an alternative to conventional cytogenetic studies for newly diagnosed AML patients with the potential to revolutionize the diagnosis of hematologic malignancies.

Copy number variant analysis using genome-wide mate-pair sequencing.

Copy number variation (CNV) is a common form of structural variation detected in human genomes, occurring as both constitutional and somatic events. Cytogenetic techniques like chromosomal microarray (CMA) are widely used in analyzing CNVs. However, CMA techniques cannot resolve the full nature of these structural variations (i.e. the orientation and location of associated breakpoint junctions) and must be combined with other cytogenetic techniques, such as karyotyping or FISH, to do so. This makes the development of a next-generation sequencing (NGS) approach capable of resolving both CNVs and breakpoint junctions desirable. Mate-pair sequencing (MPseq) is a NGS technology designed to find large structural rearrangements across the entire genome. Here we present an algorithm capable of performing copy number analysis from mate-pair sequencing data. The algorithm uses a step-wise procedure involving normalization, segmentation, and classification of the sequencing data. The segmentation technique combines both read depth and discordant mate-pair reads to increase the sensitivity and resolution of CNV calls. The method is particularly suited to MPseq, which is designed to detect breakpoint junctions at high resolution. This allows for the classification step to accurately calculate copy number levels at the relatively low read depth of MPseq. Here we compare results for a series of hematological cancer samples that were tested with CMA and MPseq. We demonstrate comparable sensitivity to the state-of-the-art CMA technology, with the benefit of improved breakpoint resolution. The algorithm provides a powerful analytical tool for the analysis of MPseq results in cancer.

Integrated mate-pair and RNA sequencing identifies novel, targetable gene fusions in peripheral T-cell lymphoma.

Peripheral T-cell lymphomas (PTCLs) represent a heterogeneous group of T-cell malignancies that generally demonstrate aggressive clinical behavior, often are refractory to standard therapy, and remain significantly understudied. The most common World Health Organization subtype is PTCL, not otherwise specified (NOS), essentially a "wastebasket" category because of inadequate understanding to assign cases to a more specific diagnostic entity. Identification of novel fusion genes has contributed significantly to improving the classification, biologic understanding, and therapeutic targeting of PTCLs. Here, we integrated mate-pair DNA and RNA next-generation sequencing to identify chromosomal rearrangements encoding expressed fusion transcripts in PTCL, NOS. Two of 11 cases had novel fusions involving VAV1, encoding a truncated form of the VAV1 guanine nucleotide exchange factor important in T-cell receptor signaling. Fluorescence in situ hybridization studies identified VAV1 rearrangements in 10 of 148 PTCLs (7%). These were observed exclusively in PTCL, NOS (11%) and anaplastic large cell lymphoma (11%). In vitro, ectopic expression of a VAV1 fusion promoted cell growth and migration in a RAC1-dependent manner. This growth was inhibited by azathioprine, a clinically available RAC1 inhibitor. We also identified novel kinase gene fusions, ITK-FER and IKZF2-ERBB4, as candidate therapeutic targets that show similarities to known recurrent oncogenic ITK-SYK fusions and ERBB4 transcript variants in PTCLs, respectively. Additional novel and potentially clinically relevant fusions also were discovered. Together, these findings identify VAV1 fusions as recurrent and targetable events in PTCLs and highlight the potential for clinical sequencing to guide individualized therapy approaches for this group of aggressive malignancies.

Chromoanasynthesis is a common mechanism that leads to ERBB2 amplifications in a cohort of early stage HER2+ breast cancer samples.

BACKGROUND: HER2 positive (HER2+) breast cancers involve chromosomal structural alterations that act as oncogenic driver events. METHODS: We interrogated the genomic structure of 18 clinically-defined HER2+ breast tumors through integrated analysis of whole genome and transcriptome sequencing, coupled with clinical information. RESULTS: ERBB2 overexpression in 15 of these tumors was associated with ERBB2 amplification due to chromoanasynthesis with six of them containing single events and the other nine exhibiting multiple events. Two of the more complex cases had adverse clinical outcomes. Chromosomes 8 was commonly involved in the same chromoanasynthesis with 17. In ten cases where chromosome 8 was involved we observed NRG1 fusions (two cases), NRG1 amplification (one case), FGFR1 amplification and ADAM32 or ADAM5 fusions. ERBB3 over-expression was associated with NRG1 fusions and EGFR and ERBB3 expressions were anti-correlated. Of the remaining three cases, one had a small duplication fully encompassing ERBB2 and was accompanied with a pathogenic mutation. CONCLUSION: Chromoanasynthesis involving chromosome 17 can lead to ERBB2 amplifications in HER2+ breast cancer. However, additional large genomic alterations contribute to a high level of genomic complexity, generating the hypothesis that worse outcome could be associated with multiple chromoanasynthetic events.

Princeton_TIGRESS 2.0: High refinement consistency and net gains through support vector machines and molecular dynamics in double-blind predictions during the CASP11 experiment.

Protein structure refinement is the challenging problem of operating on any protein structure prediction to improve its accuracy with respect to the native structure in a blind fashion. Although many approaches have been developed and tested during the last four CASP experiments, a majority of the methods continue to degrade models rather than improve them. Princeton_TIGRESS (Khoury et al., Proteins 2014;82:794-814) was developed previously and utilizes separate sampling and selection stages involving Monte Carlo and molecular dynamics simulations and classification using an SVM predictor. The initial implementation was shown to consistently refine protein structures 76% of the time in our own internal benchmarking on CASP 7-10 targets. In this work, we improved the sampling and selection stages and tested the method in blind predictions during CASP11. We added a decomposition of physics-based and hybrid energy functions, as well as a coordinate-free representation of the protein structure through distance-binning Cα-Cα distances to capture fine-grained movements. We performed parameter estimation to optimize the adjustable SVM parameters to maximize precision while balancing sensitivity and specificity across all cross-validated data sets, finding enrichment in our ability to select models from the populations of similar decoys generated for targets in CASPs 7-10. The MD stage was enhanced such that larger structures could be further refined. Among refinement methods that are currently implemented as web-servers, Princeton_TIGRESS 2.0 demonstrated the most consistent and most substantial net refinement in blind predictions during CASP11. The enhanced refinement protocol Princeton_TIGRESS 2.0 is freely available as a web server at http://atlas.engr.tamu.edu/refinement/. Proteins 2017; 85:1078-1098. © 2017 Wiley Periodicals, Inc.

Genomic rearrangements in sporadic lymphangioleiomyomatosis: an evolving genetic story.

Sporadic lymphangioleiomyomatosis is a progressive pulmonary cystic disease resulting from the infiltration of smooth muscle-like lymphangioleiomyomatosis cells into the lung. The migratory/metastasizing properties of the lymphangioleiomyomatosis cell together with the presence of somatic mutations, primarily in the tuberous sclerosis complex gene (TSC2), lead many to consider this a low-grade malignancy. As malignant tumors characteristically accumulate somatic structural variations, which have not been well studied in sporadic lymphangioleiomyomatosis, we utilized mate pair sequencing to define structural variations within laser capture microdissected enriched lymphangioleiomyomatosis cell populations from five sporadic lymphangioleiomyomatosis patients. Lymphangioleiomyomatosis cells were confirmed in each tissue by hematoxylin eosin stain review and by HMB-45 immunohistochemistry in four cases. A mutation panel demonstrated characteristic TSC2 driver mutations in three cases. Genomic profiles demonstrated normal diploid coverage across all chromosomes, with no aneuploidy or detectable gains/losses of whole chromosomal arms typical of neoplastic diseases. However, somatic rearrangements and smaller deletions were validated in the two cases which lacked TSC2 driver mutations. Most significantly, one of these sporadic lymphangioleiomyomatosis cases contained two different size deletions encompassing the entire TSC1 locus. The detection of a homozygous deletion of TSC1 driving a predicted case of sporadic lymphangioleiomyomatosis, consistent with the common two-hit TSC2 mutation model, has never been reported for sporadic lymphangioleiomyomatosis. However, while no evidence of the hereditary tuberous sclerosis complex disease was reported for this patient, the potential for mosaicism and sub-clinical phenotype cannot be ruled out. Nevertheless, this study demonstrates that somatic structural rearrangements are present in lymphangioleiomyomatosis disease and provides a novel method of genomic characterization of sporadic lymphangioleiomyomatosis cells, aiding in defining cases with no detected mutations by conventional methodologies. These structural rearrangements could represent additional pathogenic mechanisms in sporadic lymphangioleiomyomatosis disease, potentially affecting response to therapy and adding to the complex genetic story of this rare disease.

conSSert: Consensus SVM Model for Accurate Prediction of Ordered Secondary Structure.

Accurate prediction of protein secondary structure remains a crucial step in most approaches to the protein-folding problem, yet the prediction of ordered secondary structure, specifically beta-strands, remains a challenge. We developed a consensus secondary structure prediction method, conSSert, which is based on support vector machines (SVM) and provides exceptional accuracy for the prediction of beta-strands with QE accuracy of over 0.82 and a Q2-EH of 0.86. conSSert uses as input probabilities for the three types of secondary structure (helix, strand, and coil) that are predicted by four top performing methods: PSSpred, PSIPRED, SPINE-X, and RAPTOR. conSSert was trained/tested using 4261 protein chains from PDBSelect25, and 8632 chains from PISCES. Further validation was performed using targets from CASP9, CASP10, and CASP11. Our data suggest that poor performance in strand prediction is likely a result of training bias and not solely due to the nonlocal nature of beta-sheet contacts. conSSert is freely available for noncommercial use as a webservice: http://ares.tamu.edu/conSSert/.

De novo design and experimental characterization of ultrashort self-associating peptides.

Self-association is a common phenomenon in biology and one that can have positive and negative impacts, from the construction of the architectural cytoskeleton of cells to the formation of fibrils in amyloid diseases. Understanding the nature and mechanisms of self-association is important for modulating these systems and in creating biologically-inspired materials. Here, we present a two-stage de novo peptide design framework that can generate novel self-associating peptide systems. The first stage uses a simulated multimeric template structure as input into the optimization-based Sequence Selection to generate low potential energy sequences. The second stage is a computational validation procedure that calculates Fold Specificity and/or Approximate Association Affinity (K*association) based on metrics that we have devised for multimeric systems. This framework was applied to the design of self-associating tripeptides using the known self-associating tripeptide, Ac-IVD, as a structural template. Six computationally predicted tripeptides (Ac-LVE, Ac-YYD, Ac-LLE, Ac-YLD, Ac-MYD, Ac-VIE) were chosen for experimental validation in order to illustrate the self-association outcomes predicted by the three metrics. Self-association and electron microscopy studies revealed that Ac-LLE formed bead-like microstructures, Ac-LVE and Ac-YYD formed fibrillar aggregates, Ac-VIE and Ac-MYD formed hydrogels, and Ac-YLD crystallized under ambient conditions. An X-ray crystallographic study was carried out on a single crystal of Ac-YLD, which revealed that each molecule adopts a ß-strand conformation that stack together to form parallel ß-sheets. As an additional validation of the approach, the hydrogel-forming sequences of Ac-MYD and Ac-VIE were shuffled. The shuffled sequences were computationally predicted to have lower K*association values and were experimentally verified to not form hydrogels. This illustrates the robustness of the framework in predicting self-associating tripeptides. We expect that this enhanced multimeric de novo peptide design framework will find future application in creating novel self-associating peptides based on unnatural amino acids, and inhibitor peptides of detrimental self-aggregating biological proteins.