Short-read Whole-genome Sequencing Research Articles

Linked-Read Whole Genome Sequencing Solves a Double DMD Gene Rearrangement.

Next generation sequencing (NGS) has changed our approach to diagnosis of genetic disorders. Nowadays, the most comprehensive application of NGS is whole genome sequencing (WGS) that is able to detect virtually all DNA variations. However, even after accurate WGS, many genetic conditions remain unsolved. This may be due to the current NGS protocols, based on DNA fragmentation and short reads. To overcome these limitations, we applied a linked-read sequencing technology that combines single-molecule barcoding with short-read WGS. We were able to assemble haplotypes and distinguish between alleles along the genome. As an exemplary case, we studied the case of a female carrier of X-linked muscular dystrophy with an unsolved genetic status. A deletion of exons 16–29 in DMD gene was responsible for the disease in her family, but she showed a normal dosage of these exons by Multiplex Ligation-dependent Probe Amplification (MLPA) and array CGH. This situation is usually considered compatible with a “non-carrier” status. Unexpectedly, the girl also showed an increased dosage of flanking exons 1–15 and 30–34. Using linked-read WGS, we were able to distinguish between the two X chromosomes. In the first allele, we found the 16–29 deletion, while the second allele showed a 1–34 duplication: in both cases, linked-read WGS correctly mapped the borders at single-nucleotide resolution. This duplication in trans apparently restored the normal dosage of exons 16–29 seen by quantitative assays. This had a dramatic impact in genetic counselling, by converting a non-carrier into a double carrier status prediction. We conclude that linked-read WGS should be considered as a valuable option to improve our understanding of unsolved genetic conditions.

Vertical and horizontal dissemination of an IncC plasmid harbouring rmtB 16S rRNA methylase gene, conferring resistance to plazomicin, among invasive ST258 and ST16 KPC-producing Klebsiella pneumoniae

ObjectivesCarbapenem resistance in Klebsiella pneumoniae is a major clinical challenge. Aminoglycosides remain an important asset in the current therapeutic arsenal to treat these infections. We examined aminoglycoside resistance phenotypes and genomics in a collection of 100 invasive KPC-producing K. pneumoniae isolates sequentially collected in a Brazilian tertiary hospital between 2014 and 2016. MethodsAminoglycoside susceptibility testing was performed. We used a combined long-read (MinION) and short-read (Illumina) whole-genome sequencing strategy to provide a genomic picture of aminoglycoside resistance genes, with particular emphasis on 16S rRNA methyltransferases and related plasmids. Results68% of the strains were resistant to gentamicin and 42% to amikacin, with 35% resistant to both of these commonly used aminoglycosides. We identified the 16S rRNA methyltransferase gene rmtB in 30% of these isolates: 97% (29/30) belonged to sequence type 258 (ST258) and a single isolate to the emergent ST16 clone. In ST258 and ST16 the rmtB gene was located on large IncC plasmids of 177 kb and 174 kb, respectively, highly similar to a plasmid previously identified in Proteus mirabilis in the same hospital. Moreover, 99% of the isolates remained susceptible to the veterinary-approved drug apramycin, currently under clinical development for human medicine. ConclusionSuch findings in geographically and temporally related isolates suggest a combination of vertical clonal spread as well as horizontal interspecies and intraspecies plasmid transfer. This broad rmtB dissemination in an endemic setting for KPC-producing clones is worrisome since it provides resistance to most clinically available aminoglycosides, including the novel aminoglycoside-modifying enzyme-resistant plazomicin.

Targeted Long-Read Sequencing Identifies a Retrotransposon Insertion as a Cause of Altered GNAS Exon A/B Methylation in a Family With Autosomal Dominant Pseudohypoparathyroidism Type 1b (PHP1B).

Pseudohypoparathyroidism type Ib (PHP1B) is characterized predominantly by resistance to parathyroid hormone (PTH) leading to hypocalcemia and hyperphosphatemia. These laboratory abnormalities are caused by maternal loss-of-methylation (LOM) at GNAS exon A/B, which reduces in cis expression of the stimulatory G protein α-subunit (Gsα). Paternal Gsα expression in proximal renal tubules is silenced through unknown mechanisms, hence LOM at exon A/B reduces further Gsα protein in this kidney portion, leading to PTH resistance. In a previously reported PHP1B family, affected members showed variable LOM at exon A/B, yet no genetic defect was found by whole-genome sequencing despite linkage to GNAS. Using targeted long-read sequencing (T-LRS), we discovered an approximately 2800-bp maternally inherited retrotransposon insertion nearly 1200 bp downstream of exon XL not found in public databases or in 13,675 DNA samples analyzed by short-read whole-genome sequencing. T-LRS data furthermore confirmed normal methylation at exons XL, AS, and NESP and showed that LOM comprising exon A/B is broader than previously thought. The retrotransposon most likely causes the observed epigenetic defect by impairing function of a maternally derived NESP transcript, consistent with findings in mice lacking full-length NESP mRNA and in PHP1B patients with deletion of exon NESP and adjacent intronic sequences. In addition to demonstrating that T-LRS is an effective strategy for identifying a small disease-causing variant that abolishes or severely reduces exon A/B methylation, our data demonstrate that this sequencing technology has major advantages for simultaneously identifying structural defects and altered methylation. © 2022 American Society for Bone and Mineral Research (ASBMR).

Extreme enrichment of VNTR-associated polymorphicity in human subtelomeres: genes with most VNTRs are predominantly expressed in the brain

The human genome harbors numerous structural variants (SVs) which, due to their repetitive nature, are currently underexplored in short-read whole-genome sequencing approaches. Using single-molecule, real-time (SMRT) long-read sequencing technology in combination with FALCON-Unzip, we generated a de novo assembly of the diploid genome of a 115-year-old Dutch cognitively healthy woman. We combined this assembly with two previously published haploid assemblies (CHM1 and CHM13) and the GRCh38 reference genome to create a compendium of SVs that occur across five independent human haplotypes using the graph-based multi-genome aligner REVEAL. Across these five haplotypes, we detected 31,680 euchromatic SVs (>50 bp). Of these, ~62% were comprised of repetitive sequences with ‘variable number tandem repeats’ (VNTRs), ~10% were mobile elements (Alu, L1, and SVA), while the remaining variants were inversions and indels. We observed that VNTRs with GC-content >60% and repeat patterns longer than 15 bp were 21-fold enriched in the subtelomeric regions (within 5 Mb of the ends of chromosome arms). VNTR lengths can expand to exceed a critical length which is associated with impaired gene transcription. The genes that contained most VNTRs, of which PTPRN2 and DLGAP2 are the most prominent examples, were found to be predominantly expressed in the brain and associated with a wide variety of neurological disorders. Repeat-induced variation represents a sizeable fraction of the genetic variation in human genomes and should be included in investigations of genetic factors associated with phenotypic traits, specifically those associated with neurological disorders. We make available the long and short-read sequence data of the supercentenarian genome, and a compendium of SVs as identified across 5 human haplotypes.

Cytogenetically visible inversions are formed by multiple molecular mechanisms

Cytogenetically detected inversions are generally assumed to be copy number and phenotypically neutral events. While nonallelic homologous recombination is thought to play a major role, recent data suggest the involvement of other molecular mechanisms in inversion formation. Using a combination of short‐read whole‐genome sequencing (WGS), 10X Genomics Chromium WGS, droplet digital polymerase chain reaction and array comparative genomic hybridization we investigated the genomic structure of 18 large unique cytogenetically detected chromosomal inversions and achieved nucleotide resolution of at least one chromosomal inversion junction for 13/18 (72%). Surprisingly, we observed that seemingly copy number neutral inversions can be accompanied by a copy‐number gain of up to 350 kb and local genomic complexities (3/18, 17%). In the resolved inversions, the mutational signatures are consistent with nonhomologous end‐joining (8/13, 62%) or microhomology‐mediated break‐induced replication (5/13, 38%). Our study indicates that short‐read 30x coverage WGS can detect a substantial fraction of chromosomal inversions. Moreover, replication‐based mechanisms are responsible for approximately 38% of those events leading to a significant proportion of inversions that are actually accompanied by additional copy‐number variation potentially contributing to the overall phenotypic presentation of those patients.

An integrated personal and population-based Egyptian genome reference

A small number of de novo assembled human genomes have been reported to date, and few have been complemented with population-based genetic variation, which is particularly important for North Africa, a region underrepresented in current genome-wide references. Here, we combine long- and short-read whole-genome sequencing data with recent assembly approaches into a de novo assembly of an Egyptian genome. The assembly demonstrates well-balanced quality metrics and is complemented with variant phasing via linked reads into haploblocks, which we associate with gene expression changes in blood. To construct an Egyptian genome reference, we identify genome-wide genetic variation within a cohort of 110 Egyptian individuals. We show that differences in allele frequencies and linkage disequilibrium between Egyptians and Europeans may compromise the transferability of European ancestry-based genetic disease risk and polygenic scores, substantiating the need for multi-ethnic genome references. Thus, the Egyptian genome reference will be a valuable resource for precision medicine.

IgCaller for reconstructing immunoglobulin gene rearrangements and oncogenic translocations from whole-genome sequencing in lymphoid neoplasms

Immunoglobulin (Ig) gene rearrangements and oncogenic translocations are routinely assessed during the characterization of B cell neoplasms and stratification of patients with distinct clinical and biological features, with the assessment done using Sanger sequencing, targeted next-generation sequencing, or fluorescence in situ hybridization (FISH). Currently, a complete Ig characterization cannot be extracted from whole-genome sequencing (WGS) data due to the inherent complexity of the Ig loci. Here, we introduce IgCaller, an algorithm designed to fully characterize Ig gene rearrangements and oncogenic translocations from short-read WGS data. Using a cohort of 404 patients comprising different subtypes of B cell neoplasms, we demonstrate that IgCaller identifies both heavy and light chain rearrangements to provide additional information on their functionality, somatic mutational status, class switch recombination, and oncogenic Ig translocations. Our data thus support IgCaller to be a reliable alternative to Sanger sequencing and FISH for studying the genetic properties of the Ig loci.

Gplas: a comprehensive tool for plasmid analysis using short-read graphs.

SummaryPlasmids can horizontally transmit genetic traits, enabling rapid bacterial adaptation to new environments and hosts. Short-read whole-genome sequencing data are often applied to large-scale bacterial comparative genomics projects but the reconstruction of plasmids from these data is facing severe limitations, such as the inability to distinguish plasmids from each other in a bacterial genome. We developed gplas, a new approach to reliably separate plasmid contigs into discrete components using sequence composition, coverage, assembly graph information and network partitioning based on a pruned network of plasmid unitigs. Gplas facilitates the analysis of large numbers of bacterial isolates and allows a detailed analysis of plasmid epidemiology based solely on short-read sequence data.Availability and implementationGplas is written in R, Bash and uses a Snakemake pipeline as a workflow management system. Gplas is available under the GNU General Public License v3.0 at https://gitlab.com/sirarredondo/gplas.git.Supplementary information Supplementary data are available at Bioinformatics online.

Combining short and long read sequencing to characterize antimicrobial resistance genes on plasmids applied to an unauthorized genetically modified Bacillus

Antimicrobial resistance (AMR) is a major public health threat. Plasmids are able to transfer AMR genes among bacterial isolates. Whole genome sequencing (WGS) is a powerful tool to monitor AMR determinants. However, plasmids are difficult to reconstruct from WGS data. This study aimed to improve the characterization, including the localization of AMR genes using short and long read WGS strategies. We used a genetically modified (GM) Bacillus subtilis isolated as unexpected contamination in a feed additive, and therefore considered unauthorized (RASFF 2014.1249), as a case study. In GM organisms, AMR genes are used as selection markers. Because of the concern of spread of these AMR genes when present on mobile genetic elements, it is crucial to characterize their location. Our approach resulted in an assembly of one chromosome and one plasmid, each with several AMR determinants of which five are against critically important antibiotics. Interestingly, we found several plasmids, containing AMR genes, integrated in the chromosome in a repetitive region of at least 53 kb. Our findings would have been impossible using short reads only. We illustrated the added value of long read sequencing in addressing the challenges of plasmid reconstruction within the context of evaluating the risk of AMR spread.

New genomic features of the polled intersex syndrome variant in goats unraveled by long-read whole-genome sequencing.

In domestic goats, the polled intersex syndrome (PIS) refers to XX female-to-male sex reversal associated with the absence of horn growth (polled). The causal variant was previously reported as a 11.7kb deletion at approximately 129Mb on chromosome 1 that affects the transcription of both FOXL2 and several long non-coding RNAs. In the meantime the presence of different versions of the PIS deletion was postulated and trials to establish genetic testing with the existing molecular genetic information failed. Therefore, we revisited this variant by long-read whole-genome sequencing of two genetically female (XX) goats, a PIS-affected and a horned control. This revealed the presence of a more complex structural variant consisting of a deletion with a total length of 10159bp and an inversely inserted approximately 480kb-sized duplicated segment of a region located approximately 21Mb further downstream on chromosome 1 containing two genes, KCNJ15 and ERG. Publicly available short-read whole-genome sequencing data, Sanger sequencing of the breakpoints and FISH using BAC clones corresponding to both involved genome regions confirmed this structural variant. A diagnostic PCR was developed for simultaneous genotyping of carriers for this variant and determination of their genetic sex. We showed that the variant allele was present in all 334 genotyped polled goats of diverse breeds and that all analyzed 15 PIS-affected XX goats were homozygous. Our findings enable for the first time a precise genetic diagnosis for polledness and PIS in goats and add a further genomic feature to the complexity of the PIS phenomenon.

Discordant bioinformatic predictions of antimicrobial resistance from whole-genome sequencing data of bacterial isolates: an inter-laboratory study.

Antimicrobial resistance (AMR) poses a threat to public health. Clinical microbiology laboratories typically rely on culturing bacteria for antimicrobial-susceptibility testing (AST). As the implementation costs and technical barriers fall, whole-genome sequencing (WGS) has emerged as a ‘one-stop’ test for epidemiological and predictive AST results. Few published comparisons exist for the myriad analytical pipelines used for predicting AMR. To address this, we performed an inter-laboratory study providing sets of participating researchers with identical short-read WGS data from clinical isolates, allowing us to assess the reproducibility of the bioinformatic prediction of AMR between participants, and identify problem cases and factors that lead to discordant results. We produced ten WGS datasets of varying quality from cultured carbapenem-resistant organisms obtained from clinical samples sequenced on either an Illumina NextSeq or HiSeq instrument. Nine participating teams (‘participants’) were provided these sequence data without any other contextual information. Each participant used their choice of pipeline to determine the species, the presence of resistance-associated genes, and to predict susceptibility or resistance to amikacin, gentamicin, ciprofloxacin and cefotaxime. We found participants predicted different numbers of AMR-associated genes and different gene variants from the same clinical samples. The quality of the sequence data, choice of bioinformatic pipeline and interpretation of the results all contributed to discordance between participants. Although much of the inaccurate gene variant annotation did not affect genotypic resistance predictions, we observed low specificity when compared to phenotypic AST results, but this improved in samples with higher read depths. Had the results been used to predict AST and guide treatment, a different antibiotic would have been recommended for each isolate by at least one participant. These challenges, at the final analytical stage of using WGS to predict AMR, suggest the need for refinements when using this technology in clinical settings. Comprehensive public resistance sequence databases, full recommendations on sequence data quality and standardization in the comparisons between genotype and resistance phenotypes will all play a fundamental role in the successful implementation of AST prediction using WGS in clinical microbiology laboratories.

From cytogenetics to cytogenomics: whole-genome sequencing as a first-line test comprehensively captures the diverse spectrum of disease-causing genetic variation underlying intellectual disability

BackgroundSince different types of genetic variants, from single nucleotide variants (SNVs) to large chromosomal rearrangements, underlie intellectual disability, we evaluated the use of whole-genome sequencing (WGS) rather than chromosomal microarray analysis (CMA) as a first-line genetic diagnostic test.MethodsWe analyzed three cohorts with short-read WGS: (i) a retrospective cohort with validated copy number variants (CNVs) (cohort 1, n = 68), (ii) individuals referred for monogenic multi-gene panels (cohort 2, n = 156), and (iii) 100 prospective, consecutive cases referred to our center for CMA (cohort 3). Bioinformatic tools developed include FindSV, SVDB, Rhocall, Rhoviz, and vcf2cytosure.ResultsFirst, we validated our structural variant (SV)-calling pipeline on cohort 1, consisting of three trisomies and 79 deletions and duplications with a median size of 850 kb (min 500 bp, max 155 Mb). All variants were detected. Second, we utilized the same pipeline in cohort 2 and analyzed with monogenic WGS panels, increasing the diagnostic yield to 8%. Next, cohort 3 was analyzed by both CMA and WGS. The WGS data was processed for large (> 10 kb) SVs genome-wide and for exonic SVs and SNVs in a panel of 887 genes linked to intellectual disability as well as genes matched to patient-specific Human Phenotype Ontology (HPO) phenotypes. This yielded a total of 25 pathogenic variants (SNVs or SVs), of which 12 were detected by CMA as well. We also applied short tandem repeat (STR) expansion detection and discovered one pathologic expansion in ATXN7. Finally, a case of Prader-Willi syndrome with uniparental disomy (UPD) was validated in the WGS data.Important positional information was obtained in all cohorts. Remarkably, 7% of the analyzed cases harbored complex structural variants, as exemplified by a ring chromosome and two duplications found to be an insertional translocation and part of a cryptic unbalanced translocation, respectively.ConclusionThe overall diagnostic rate of 27% was more than doubled compared to clinical microarray (12%). Using WGS, we detected a wide range of SVs with high accuracy. Since the WGS data also allowed for analysis of SNVs, UPD, and STRs, it represents a powerful comprehensive genetic test in a clinical diagnostic laboratory setting.

Plasmid Reconstruction from Next-Gen Data: A Detailed Protocol for the Use of PLACNETw for the Reconstruction of Plasmids from WGS Datasets.

Mobile Genetic Elements (MGE) play essential roles in adaptive bacterial evolution, facilitating genetic exchange for extrachromosomal DNA, especially antibiotic resistance genes and virulence factors. For this reason, high-throughput next-generation sequencing of bacteria is of great relevance, especially for clinical pathogenic bacteria. Accurate identification of MGE from whole-genome sequencing (WGS) datasets is one of the major challenges, still hindered by methodological limitations and high sequencing costs.This chapter encompasses the protocol used for plasmid reconstruction by applying the PLACNETw methodology, from raw reads to assembled plasmids and chromosome. PLACNETw is a graphical user-friendly interface to visualize and reconstruct MGE from short-read WGS datasets. No bioinformatic background or sophisticated computational resources are required and high precision and sensitivity are achieved.

Draft Genome Sequences of Acinetobacter baumannii Isolates Recovered from Sewage Water from a Poultry Slaughterhouse in Germany

Acinetobacter baumannii is an important human pathogen usually associated with severe hospital-acquired infections. Here, we announce the draft genome sequences of two livestock-associated isolates recovered from sewage water from a poultry slaughterhouse in Germany. Short-read whole-genome sequencing was conducted to determine the genetic basis of their antimicrobial resistance phenotype.

Misassembly of long reads undermines de novo-assembled ethnicity-specific genomes: validation in a Chinese Han population.

An ethnicity is characterized by genomic fragments, single nucleotide polymorphisms (SNPs), and structural variations specific to it. However, the widely used 'standard human reference genome' GRCh37/38 is based on Caucasians. Therefore, de novo-assembled reference genomes for specific ethnicities would have advantages for genetics and precision medicine applications, especially with the long-read sequencing techniques that facilitate genome assembly. In this study, we assessed the de novo-assembled Chinese Han reference genome HX1 vis-à-vis the standard GRCh38 for improving the quality of assembly and for ethnicity-specific applications. Surprisingly, all genomic sequencing datasets mapped better to GRCh38 than to HX1, even for the datasets of the Chinese Han population. This gap was mainly due to the massive structural misassembly of the HX1 reference genome rather than the SNPs between the ethnicities, and this misassembly could not be corrected by short-read whole-genome sequencing (WGS). For example, HX1 and the other de novo-assembled personal genomes failed to assemble the mitochondrial genome as a contig. We mapped 97.1% of dbSNP, 98.8% of ClinVar, and 97.2% of COSMIC variants to HX1. HX1-absent, non-synonymous ClinVar SNPs were involved in 140 genes and many important functions in various diseases, most of which were due to the assembly failure of essential exons. In contrast, the HX1-specific regions were scantly expressible, as shown in the cell lines and clinical samples of Chinese patients. Our results demonstrated that the de novo-assembled individual genome such as HX1 did not have advantages against the standard GRCh38 genome due to insufficient assembly quality, and that it is, therefore, not recommended for common use.

ConFindr: rapid detection of intraspecies and cross-species contamination in bacterial whole-genome sequence data.

Whole-genome sequencing (WGS) of bacterial pathogens is currently widely used to support public-health investigations. The ability to assess WGS data quality is critical to underpin the reliability of downstream analyses. Sequence contamination is a quality issue that could potentially impact WGS-based findings; however, existing tools do not readily identify contamination from closely-related organisms. To address this gap, we have developed a computational pipeline, ConFindr, for detection of intraspecies contamination. ConFindr determines the presence of contaminating sequences based on the identification of multiple alleles of core, single-copy, ribosomal-protein genes in raw sequencing reads. The performance of this tool was assessed using simulated and lab-generated Illumina short-read WGS data with varying levels of contamination (0–20% of reads) and varying genetic distance between the designated target and contaminant strains. Intraspecies and cross-species contamination was reliably detected in datasets containing 5% or more reads from a second, unrelated strain. ConFindr detected intraspecies contamination with higher sensitivity than existing tools, while also being able to automatically detect cross-species contamination with similar sensitivity. The implementation of ConFindr in quality-control pipelines will help to improve the reliability of WGS databases as well as the accuracy of downstream analyses. ConFindr is written in Python, and is freely available under the MIT License at github.com/OLC-Bioinformatics/ConFindr.

Resolving the full spectrum of human genome variation using Linked-Reads.

Large-scale population analyses coupled with advances in technology have demonstrated that the human genome is more diverse than originally thought. To date, this diversity has largely been uncovered using short-read whole-genome sequencing. However, these short-read approaches fail to give a complete picture of a genome. They struggle to identify structural events, cannot access repetitive regions, and fail to resolve the human genome into haplotypes. Here, we describe an approach that retains long range information while maintaining the advantages of short reads. Starting from ∼1 ng of high molecular weight DNA, we produce barcoded short-read libraries. Novel informatic approaches allow for the barcoded short reads to be associated with their original long molecules producing a novel data type known as “Linked-Reads”. This approach allows for simultaneous detection of small and large variants from a single library. In this manuscript, we show the advantages of Linked-Reads over standard short-read approaches for reference-based analysis. Linked-Reads allow mapping to 38 Mb of sequence not accessible to short reads, adding sequence in 423 difficult-to-sequence genes including disease-relevant genes STRC, SMN1, and SMN2. Both Linked-Read whole-genome and whole-exome sequencing identify complex structural variations, including balanced events and single exon deletions and duplications. Further, Linked-Reads extend the region of high-confidence calls by 68.9 Mb. The data presented here show that Linked-Reads provide a scalable approach for comprehensive genome analysis that is not possible using short reads alone.

Comprehensive structural variation genome map of individuals carrying complex chromosomal rearrangements.

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.

Complex structural variants in Mendelian disorders: identification and breakpoint resolution using short- and long-read genome sequencing

BackgroundStudies have shown that complex structural variants (cxSVs) contribute to human genomic variation and can cause Mendelian disease. We aimed to identify cxSVs relevant to Mendelian disease using short-read whole-genome sequencing (WGS), resolve the precise variant configuration and investigate possible mechanisms of cxSV formation.MethodsWe performed short-read WGS and analysis of breakpoint junctions to identify cxSVs in a cohort of 1324 undiagnosed rare disease patients. Long-read WGS and gene expression analysis were used to resolve one case.ResultsWe identified three pathogenic cxSVs: a de novo duplication-inversion-inversion-deletion affecting ARID1B, a de novo deletion-inversion-duplication affecting HNRNPU and a homozygous deletion-inversion-deletion affecting CEP78. Additionally, a de novo duplication-inversion-duplication overlapping CDKL5 was resolved by long-read WGS demonstrating the presence of both a disrupted and an intact copy of CDKL5 on the same allele, and gene expression analysis showed both parental alleles of CDKL5 were expressed. Breakpoint analysis in all the cxSVs revealed both microhomology and longer repetitive elements.ConclusionsOur results corroborate that cxSVs cause Mendelian disease, and we recommend their consideration during clinical investigations. We show that resolution of breakpoints can be critical to interpret pathogenicity and present evidence of replication-based mechanisms in cxSV formation.

TETyper: a bioinformatic pipeline for classifying variation and genetic contexts of transposable elements from short-read whole-genome sequencing data.

Much of the worldwide dissemination of antibiotic resistance has been driven by resistance gene associations with mobile genetic elements (MGEs), such as plasmids and transposons. Although increasing, our understanding of resistance spread remains relatively limited, as methods for tracking mobile resistance genes through multiple species, strains and plasmids are lacking. We have developed a bioinformatic pipeline for tracking variation within, and mobility of, specific transposable elements (TEs), such as transposons carrying antibiotic-resistance genes. TETyper takes short-read whole-genome sequencing data as input and identifies single-nucleotide mutations and deletions within the TE of interest, to enable tracking of specific sequence variants, as well as the surrounding genetic context(s), to enable identification of transposition events. A major advantage of TETyper over previous methods is that it does not require a genome reference. To investigate global dissemination of Klebsiella pneumoniae carbapenemase (KPC) and its associated transposon Tn4401, we applied TETyper to a collection of over 3000 publicly available Illumina datasets containing blaKPC. This revealed surprising diversity, with over 200 distinct flanking genetic contexts for Tn4401, indicating high levels of transposition. Integration of sample metadata revealed insights into associations between geographic locations, host species, Tn4401 sequence variants and flanking genetic contexts. To demonstrate the ability of TETyper to cope with high-copy-number TEs and to track specific short-term evolutionary changes, we also applied it to the insertion sequence IS26 within a defined K. pneumoniae outbreak. TETyper is implemented in python and is freely available at https://github.com/aesheppard/TETyper.