Mobile Element Insertions Research Articles

Insertion sequence elements and unique symmetrical genomic regions mediate chromosomal inversions in Streptococcus pyogenes.

Chromosomal inversions are a phenomenon in many bacterial species, often across the axis of replication. Inversions have been shown to alter gene expression, changing persistence of colonisation and infection following environmental stresses. In Streptococcus pyogenes, inversions have been reported. However, frequency and molecular markers of inversions have not been systematically examined. Here, 249 complete S.pyogenes genomes were analysed using a pangenomic core gene synteny framework to identify sequences associated with inversions. 47% of genomes (118/249) contained at least one inversion, from 23 unique inversion locations. Chromosomal locations enabling inversions were usually associated with mobile elements (insertion sequences n=9 and prophages n=7). Two insertion sequences, IS1548 and IS1239, accounted for>80% of insertion sequences and were the only insertion sequences associated with inversions. The most observed inversion location (n=104 genomes, 88% of genomes with an inversion) occurs between two conserved regions encoding rRNAs, tRNAsand sigma factor genes. The regions are symmetrically placed around the origin of replication forming a unique chromosomal structure inS. pyogenes, relative to other streptococci. Cataloging of thechromosomal location and frequency of inversions can direct dissection of phenotypic changes following chromosomal inversions. The framework used here can be transferred to other bacterial species to characterise chromosomal inversions.

Complex genetic variation in nearly complete human genomes.

Diverse sets of complete human genomes are required to construct a pangenome reference and to understand the extent of complex structural variation. Here, we sequence 65 diverse human genomes and build 130 haplotype-resolved assemblies (130 Mbp median continuity), closing 92% of all previous assembly gaps1,2 and reaching telomere-to-telomere (T2T) status for 39% of the chromosomes. We highlight complete sequence continuity of complex loci, including the major histocompatibility complex (MHC), SMN1/SMN2, NBPF8, and AMY1/AMY2, and fully resolve 1,852 complex structural variants (SVs). In addition, we completely assemble and validate 1,246 human centromeres. We find up to 30-fold variation in α-satellite high-order repeat (HOR) array length and characterize the pattern of mobile element insertions into α-satellite HOR arrays. While most centromeres predict a single site of kinetochore attachment, epigenetic analysis suggests the presence of two hypomethylated regions for 7% of centromeres. Combining our data with the draft pangenome reference1 significantly enhances genotyping accuracy from short-read data, enabling whole-genome inference3 to a median quality value (QV) of 45. Using this approach, 26,115 SVs per sample are detected, substantially increasing the number of SVs now amenable to downstream disease association studies.

Long-read genome sequencing and variant reanalysis increase diagnostic yield in neurodevelopmental disorders.

Variant detection from long-read genome sequencing (lrGS) has proven to be more accurate and comprehensive than variant detection from short-read genome sequencing (srGS). However, the rate at which lrGS can increase molecular diagnostic yield for rare disease is not yet precisely characterized. We performed lrGS using Pacific Biosciences HiFi technology on 96 short-read-negative probands with rare diseases that were suspected to be genetic. We generated hg38-aligned variants and de novo phased genome assemblies, and subsequently annotated, filtered, and curated variants using clinical standards. New disease-relevant or potentially relevant genetic findings were identified in 16/96 (16.7%) probands, nine of which (8/96, ~9.4%) harbored pathogenic or likely pathogenic variants. Nine probands (~9.4%) had variants that were accurately called in both srGS and lrGS and represent changes to clinical interpretation, mostly from recently published gene-disease associations. Seven cases included variants that were only correctly interpreted in lrGS, including copy-number variants, an inversion, a mobile element insertion, two low-complexity repeat expansions, and a 1 bp deletion. While evidence for each of these variants is, in retrospect, visible in srGS, they were either not called within srGS data, were represented by calls with incorrect sizes or structures, or failed quality-control and filtration. Thus, while reanalysis of older srGS data clearly increases diagnostic yield, we find that lrGS allows for substantial additional yield (7/96, 7.3%) beyond srGS. We anticipate that as lrGS analysis improves, and as lrGS datasets grow allowing for better variant frequency annotation, the additional lrGS-only rare disease yield will grow over time.

Harbinger transposon insertion in ethylene signaling gene leads to emergence of new sexual forms in cucurbits

In flowering plants, the predominant sexual morph is hermaphroditism, and the emergence of unisexuality is poorly understood. Using Cucumis melo (melon) as a model system, we explore the mechanisms driving sexual forms. We identify a spontaneous mutant exhibiting a transition from bisexual to unisexual male flower, and identify the causal mutation as a Harbinger transposon impairing the expression of Ethylene Insensitive 2 (CmEIN2) gene. Genetics and transcriptomic analysis reveal a dual role of CmEIN2 in both sex determination and fruit shape formation. Upon expression of CmACS11, EIN2 is recruited to repress the expression of the carpel inhibitor, CmWIP1. Subsequently, EIN2 is recruited to mediate stamina inhibition. Following the sex determination phase, EIN2 promotes fruit shape elongation. Genome-wide analysis reveals that Harbinger transposon mobilization is triggered by environmental cues, and integrates preferentially in active chromatin, particularly within promoter regions. Characterization of a large collection of melon germplasm points to active transpositions in the wild, compared to cultivated accessions. Our study underscores the association between chromatin dynamics and the temporal aspects of mobile genetic element insertions, providing valuable insights into plant adaptation and crop genome evolution.

Analyzing imbalance of short tandem repeats for pancreatic cancer detection.

e16359 Background: Short tandem repeats (STRs), alternatively known as microsatellites, encompass repeat units of 1 to 6 base pairs (bp) in the genome and function as regulators of gene expression, impacting plasticity and complex traits. STRs are prone to hypermutation with alterations correlated to pathogenicity in multiple Mendelian and acquired disorders. STR variant classes include insertions, mobile element insertions (MEIs), deletions, and multi-allelic copy number variants (mCNVs). These can lead to repeat expansions, allelic imbalance (AI), and microsatellite instability (MSI) observed in diseases such as Huntington's disease, hereditary ataxia, and multiple cancers. In pancreatic cancer, AI and associated loss-of-heterozygosity (LOH) at tumor suppressor loci flanked by hypervariable STR regions correspond to a pathogenic phenotype with a risk of aggressive disease. Consequently, it is critical to develop accurate and robust methods to precisely identify and genotype pancreatic samples with varying specimen quality in a clinical setting to aid in diagnosis. Methods: Capillary electrophoresis (CE) and Sanger sequencing have been used to analyze STRs. The availability of newer tools such as 2nd generation short-read (NGS) and 3rd generation long-read (TGS) sequencers could allow for high throughput scaling, however, informatics challenges associated with analyzing repetitive regions remain. Here we developed amplification and enrichment methods for targeted analysis for detection by CE, NGS and TGS. Results: Amplification and enrichment strategies for detecting 17 short tandem repeat regions that play a role in pancreatic carcinogenesis were developed. Methods were optimized for detecting STR alterations in cell-free nucleic acids isolated from endoscope-guided pancreatic fine needle aspirate biopsies. Furthermore, informatics tools were developed to analyze short-read and long-read sequencing data to accurately analyze repeat regions. Conclusions: Here we compared 3 platforms (CE, NGS, and TGS) and found that CE remains the gold standard for analyzing AI/LOH with data comparable to long-read TGS, while short-read NGS is not always amenable to accurate analysis dependent on the region of interest. Additionally, we found that artifacts of amplification persist, irrespective of the 3 platforms analyzed.

Disentangling oncogenic amplicons in esophageal adenocarcinoma

Esophageal adenocarcinoma is a prominent example of cancer characterized by frequent amplifications in oncogenes. However, the mechanisms leading to amplicons that involve breakage-fusion-bridge cycles and extrachromosomal DNA are poorly understood. Here, we use 710 esophageal adenocarcinoma cases with matched samples and patient-derived organoids to disentangle complex amplicons and their associated mechanisms. Short-read sequencing identifies ERBB2, MYC, MDM2, and HMGA2 as the most frequent oncogenes amplified in extrachromosomal DNAs. We resolve complex extrachromosomal DNA and breakage-fusion-bridge cycles amplicons by integrating of de-novo assemblies and DNA methylation in nine long-read sequenced cases. Complex amplicons shared between precancerous biopsy and late-stage tumor, an enrichment of putative enhancer elements and mobile element insertions are potential drivers of complex amplicons’ origin. We find that patient-derived organoids recapitulate extrachromosomal DNA observed in the primary tumors and single-cell DNA sequencing capture extrachromosomal DNA-driven clonal dynamics across passages. Prospectively, long-read and single-cell DNA sequencing technologies can lead to better prediction of clonal evolution in esophageal adenocarcinoma.

Mapping and functional characterization of structural variation in 1060 pig genomes

BackgroundStructural variations (SVs) have significant impacts on complex phenotypes by rearranging large amounts of DNA sequence.ResultsWe present a comprehensive SV catalog based on the whole-genome sequence of 1060 pigs (Sus scrofa) representing 101 breeds, covering 9.6% of the pig genome. This catalog includes 42,487 deletions, 37,913 mobile element insertions, 3308 duplications, 1664 inversions, and 45,184 break ends. Estimates of breed ancestry and hybridization using genotyped SVs align well with those from single nucleotide polymorphisms. Geographically stratified deletions are observed, along with known duplications of the KIT gene, responsible for white coat color in European pigs. Additionally, we identify a recent SINE element insertion in MYO5A transcripts of European pigs, potentially influencing alternative splicing patterns and coat color alterations. Furthermore, a Yorkshire-specific copy number gain within ABCG2 is found, impacting chromatin interactions and gene expression across multiple tissues over a stretch of genomic region of ~200 kb. Preliminary investigations into SV’s impact on gene expression and traits using the Pig Genotype-Tissue Expression (PigGTEx) data reveal SV associations with regulatory variants and gene-trait pairs. For instance, a 51-bp deletion is linked to the lead eQTL of the lipid metabolism regulating gene FADS3, whose expression in embryo may affect loin muscle area, as revealed by our transcriptome-wide association studies.ConclusionsThis SV catalog serves as a valuable resource for studying diversity, evolutionary history, and functional shaping of the pig genome by processes like domestication, trait-based breeding, and adaptive evolution.

Long-read sequencing improves the genetic diagnosis of retinitis pigmentosa by identifying an Alu retrotransposon insertion in the EYS gene

BackgroundBiallelic variants in EYS are the major cause of autosomal recessive retinitis pigmentosa (arRP) in certain populations, a clinically and genetically heterogeneous disease that may lead to legal blindness. EYS is one of the largest genes (~ 2 Mb) expressed in the retina, in which structural variants (SVs) represent a common cause of disease. However, their identification using short-read sequencing (SRS) is not always feasible. Here, we conducted targeted long-read sequencing (T-LRS) using adaptive sampling of EYS on the MinION sequencing platform (Oxford Nanopore Technologies) to definitively diagnose an arRP family, whose affected individuals (n = 3) carried the heterozygous pathogenic deletion of exons 32–33 in the EYS gene. As this was a recurrent variant identified in three additional families in our cohort, we also aimed to characterize the known deletion at the nucleotide level to assess a possible founder effect.ResultsT-LRS in family A unveiled a heterozygous AluYa5 insertion in the coding exon 43 of EYS (chr6(GRCh37):g.64430524_64430525ins352), which segregated with the disease in compound heterozygosity with the previously identified deletion. Visual inspection of previous SRS alignments using IGV revealed several reads containing soft-clipped bases, accompanied by a slight drop in coverage at the Alu insertion site. This prompted us to develop a simplified program using grep command to investigate the recurrence of this variant in our cohort from SRS data. Moreover, LRS also allowed the characterization of the CNV as a ~ 56.4kb deletion spanning exons 32–33 of EYS (chr6(GRCh37):g.64764235_64820592del). The results of further characterization by Sanger sequencing and linkage analysis in the four families were consistent with a founder variant.ConclusionsTo our knowledge, this is the first report of a mobile element insertion into the coding sequence of EYS, as a likely cause of arRP in a family. Our study highlights the value of LRS technology in characterizing and identifying hidden pathogenic SVs, such as retrotransposon insertions, whose contribution to the etiopathogenesis of rare diseases may be underestimated.

Clinician-Driven Reanalysis of Exome Sequencing Data From Patients With Inherited Retinal Diseases

Despite advances in next-generation sequencing (NGS), a significant proportion of patients with inherited retinal disease (IRD) remain undiagnosed after initial genetic testing. Exome sequencing (ES) reanalysis in the clinical setting has been suggested as one method for improving diagnosis of IRD. To investigate the association of clinician-led reanalysis of ES data, which incorporates updated clinical information and comprehensive bioinformatic analysis, with the diagnostic yield in a cohort of patients with IRDs in Korea. This was a multicenter prospective cohort study involving 264 unrelated patients with IRDs, conducted in Korea between March 2018 and February 2020. Comprehensive ophthalmologic examinations and ES analyses were performed, and ES data were reanalyzed by an IRD specialist for single nucleotide variants, copy number variants, mobile element insertions, and mitochondrial variants. Data were analyzed from March to July 2023. Diagnostic rate of conventional bioinformatic analysis and clinician-driven ES reanalysis. A total of 264 participants (151 [57.2%] male; mean [SD] age at genetic testing, 33.6 [18.9] years) were enrolled, including 129 patients (48.9%) with retinitis pigmentosa and 26 patients (9.8%) with Stargardt disease or macular dystrophy. Initial bioinformatic analysis diagnosed 166 patients (62.9%). Clinician-driven reanalysis identified the molecular cause of diseases in an additional 22 patients, corresponding to an 8.3-percentage point increase in diagnostic rate. Key factors associated with new molecular diagnoses included clinical phenotype updates (4 patients) and detection of previously overlooked variation, such as structural variants (9 patients), mitochondrial variants (3 patients), filtered or not captured variants (4 patients), and noncanonical splicing variants (2 patients). Among the 22 patients, variants in 7 patients (31.8%) were observed in the initial analysis but not reported to patients, while those in the remaining 15 patients (68.2%) were newly detected by the ES reanalysis. In this cohort study, clinician-centered reanalysis of ES data was associated with improved molecular diagnostic yields in patients with IRD. This approach is important for uncovering missed genetic causes of retinal disease.

Molecular characterization of lugdunin inactivation mechanisms and their association with Staphylococcus lugdunensis genetic types

Background and purposeOur previous studies showed that lugdunin activities are associated with Staphylococcus lugdunensis genotypes, and most isolates do not exhibit lugdunin activity. As a continuation of our previous analysis, we focused on the reasons for defects in lugdunin production in S. lugdunensis clinical isolates. MethodsA comparative analysis of 36 S. lugdunensis whole genome sequencing data revealed three major mutation types, unknown deletion mechanism that caused most of lug operon genes lost, mobile genetic element (MGE) insertion, and nonsense mutations, which potentially damaged lugdunin production. A total of 152 S. lugdunensis clinical isolates belonging to lugdunin nonproducers were further examined for the above three mutation types. PCR products were sequenced to examine these variations. ResultsForty-six of the 152 isolates were CRISPR-Cas IIC isolates, including 26 ST27, 14 ST4, and 6 ST29 isolates; further investigation confirmed that all of their lug operons had lost almost all lug operon genes except lugM. An IS256 insertion in lugA was identified in 16 isolates, and most isolates (15 over 16) belonged to ST3. In addition, three nonsense mutations caused by single nucleotide substitutions (an adenine deletion in lugB at the 361th and 1219th nucleotides and an adenine deletion in lugC at the 1612nd nucleotide) that were frequently observed among 36 S. lugdunensis whole genome sequencing data were further observed in our clinical isolates. These three nonsense mutations were frequently found in most of CRISPR-Cas IIIA strains, especially in ST6 isolates. ConclusionOur findings suggest that the mechanisms affecting lugdunin production are associated with S. lugdunensis molecular types.

Genome sequencing detects a wide range of clinically relevant copy-number variants and other genomic alterations

PurposeCopy-number variants (CNVs) and other non–single nucleotide variant/indel variant types contribute an important proportion of diagnoses in individuals with suspected genetic disease. This study describes the range of such variants detected by genome sequencing (GS). MethodsFor a pediatric cohort of 1032 participants undergoing clinical GS, we characterize the CNVs and other non–single nucleotide variant/indel variant types that were reported, including aneuploidies, mobile element insertions, and uniparental disomies, and we describe the bioinformatic pipeline used to detect these variants. ResultsTogether, these genetic alterations accounted for 15.8% of reported variants. Notably, 67.9% of these were deletions, 32.9% of which overlapped a single gene, and many deletions were reported together with a second variant in the same gene in cases of recessive disease. A retrospective medical record review in a subset of this cohort revealed that up to 6 additional genetic tests were ordered in 68% (26/38) of cases, some of which failed to report the CNVs/rare variants reported on GS. ConclusionGS detected a broad range of reported variant types, including CNVs ranging in size from 1 Kb to 46 Mb.

Mobile element insertions in rare diseases: a comparative benchmark and reanalysis of 60,000 exome samples

Mobile element insertions (MEIs) are a known cause of genetic disease but have been underexplored due to technical limitations of genetic testing methods. Various bioinformatic tools have been developed to identify MEIs in Next Generation Sequencing data. However, most tools have been developed specifically for genome sequencing (GS) data rather than exome sequencing (ES) data, which remains more widely used for routine diagnostic testing. In this study, we benchmarked six MEI detection tools (ERVcaller, MELT, Mobster, SCRAMble, TEMP2 and xTea) on ES data and on GS data from publicly available genomic samples (HG002, NA12878). For all the tools we evaluated sensitivity and precision of different filtering strategies. Results show that there were substantial differences in tool performance between ES and GS data. MELT performed best with ES data and its combination with SCRAMble increased substantially the detection rate of MEIs. By applying both tools to 10,890 ES samples from Solve-RD and 52,624 samples from Radboudumc we were able to diagnose 10 patients who had remained undiagnosed by conventional ES analysis until now. Our study shows that MELT and SCRAMble can be used reliably to identify clinically relevant MEIs in ES data. This may lead to an additional diagnosis for 1 in 3000 to 4000 patients in routine clinical ES.

The landscape of human SVA retrotransposons.

SINE-VNTR-Alu (SVA) retrotransposons are evolutionarily young and still-active transposable elements (TEs) in the human genome. Several pathogenic SVA insertions have been identified that directly mutate host genes to cause neurodegenerative and other types of diseases. However, due to their sequence heterogeneity and complex structures as well as limitations in sequencing techniques and analysis, SVA insertions have been less well studied compared to other mobile element insertions. Here, we identified polymorphic SVA insertions from 3646 whole-genome sequencing (WGS) samples of>150 diverse populations and constructed a polymorphic SVA insertion reference catalog. Using 20 long-read samples, we also assembled reference and polymorphic SVA sequences and characterized the internal hexamer/variable-number-tandem-repeat (VNTR) expansions as well as differing SVA activity for SVA subfamilies and human populations. In addition, we developed a module to annotate both reference and polymorphic SVA copies. By characterizing the landscape of both reference and polymorphic SVA retrotransposons, our study enables more accurate genotyping of these elements and facilitate the discovery of pathogenic SVA insertions.

Emerging Opportunities to Study Mobile Element Insertions and Their Source Elements in an Expanding Universe of Sequenced Human Genomes.

Three mobile element classes, namely Alu, LINE-1 (L1), and SVA elements, remain actively mobile in human genomes and continue to produce new mobile element insertions (MEIs). Historically, MEIs have been discovered and studied using several methods, including: (1) Southern blots, (2) PCR (including PCR display), and (3) the detection of MEI copies from young subfamilies. We are now entering a new phase of MEI discovery where these methods are being replaced by whole genome sequencing and bioinformatics analysis to discover novel MEIs. We expect that the universe of sequenced human genomes will continue to expand rapidly over the next several years, both with short-read and long-read technologies. These resources will provide unprecedented opportunities to discover MEIs and study their impact on human traits and diseases. They also will allow the MEI community to discover and study the source elements that produce these new MEIs, which will facilitate our ability to study source element regulation in various tissue contexts and disease states. This, in turn, will allow us to better understand MEI mutagenesis in humans and the impact of this mutagenesis on human biology.

Associations of genome-wide structural variations with phenotypic differences in cross-bred Eurasian pigs

BackgroundDuring approximately 10,000 years of domestication and selection, a large number of structural variations (SVs) have emerged in the genome of pig breeds, profoundly influencing their phenotypes and the ability to adapt to the local environment. SVs (≥ 50 bp) are widely distributed in the genome, mainly in the form of insertion (INS), mobile element insertion (MEI), deletion (DEL), duplication (DUP), inversion (INV), and translocation (TRA). While studies have investigated the SVs in pig genomes, genome-wide association studies (GWAS)-based on SVs have been rarely conducted.ResultsHere, we obtained a high-quality SV map containing 123,151 SVs from 15 Large White and 15 Min pigs through integrating the power of several SV tools, with 53.95% of the SVs being reported for the first time. These high-quality SVs were used to recover the population genetic structure, confirming the accuracy of genotyping. Potential functional SV loci were then identified based on positional effects and breed stratification. Finally, GWAS were performed for 36 traits by genotyping the screened potential causal loci in the F2 population according to their corresponding genomic positions. We identified a large number of loci involved in 8 carcass traits and 6 skeletal traits on chromosome 7, with FKBP5 containing the most significant SV locus for almost all traits. In addition, we found several significant loci in intramuscular fat, abdominal circumference, heart weight, and liver weight, etc.ConclusionsWe constructed a high-quality SV map using high-coverage sequencing data and then analyzed them by performing GWAS for 25 carcass traits, 7 skeletal traits, and 4 meat quality traits to determine that SVs may affect body size between European and Chinese pig breeds.

Structural Variations Contribute to the Genetic Etiology of Autism Spectrum Disorder and Language Impairments.

Autism spectrum disorder (ASD) is a neurodevelopmental disorder characterized by restrictive interests and/or repetitive behaviors and deficits in social interaction and communication. ASD is a multifactorial disease with a complex polygenic genetic architecture. Its genetic contributing factors are not yet fully understood, especially large structural variations (SVs). In this study, we aimed to assess the contribution of SVs, including copy number variants (CNVs), insertions, deletions, duplications, and mobile element insertions, to ASD and related language impairments in the New Jersey Language and Autism Genetics Study (NJLAGS) cohort. Within the cohort, ~77% of the families contain SVs that followed expected segregation or de novo patterns and passed our filtering criteria. These SVs affected 344 brain-expressed genes and can potentially contribute to the genetic etiology of the disorders. Gene Ontology and protein-protein interaction network analysis suggested several clusters of genes in different functional categories, such as neuronal development and histone modification machinery. Genes and biological processes identified in this study contribute to the understanding of ASD and related neurodevelopment disorders.

Extensive Independent Amplification of Platy-1 Retroposons in Tamarins, Genus Saguinus.

Platy-1 retroposons are short interspersed elements (SINEs) unique to platyrrhine primates. Discovered in the common marmoset (Callithrix jacchus) genome in 2016, these 100 bp mobile element insertions (MEIs) appeared to be novel drivers of platyrrhine evolution, with over 2200 full-length members across 62 different subfamilies, and strong evidence of ongoing proliferation in C. jacchus. Subsequent characterization of Platy-1 elements in Aotus, Saimiri and Cebus genera, suggested that the widespread mobilization detected in marmoset (family Callithrichidae) was perhaps an anomaly. Two additional Callithrichidae genomes are now available, a scaffold level genome assembly for Saguinus imperator (tamarin; SagImp_v1) and a chromosome-level assembly for Saguinus midas (Midas tamarin; ASM2_v1). Here, we report that each tamarin genome contains over 11,000 full-length Platy-1 insertions, about 1150 are shared by both Saguinus tamarins, 7511 are unique to S. imperator, and another 8187 are unique to S. midas. Roughly 325 are shared among the three callithrichids. We identified six new Platy-1 subfamilies derived from Platy-1-8, with the youngest new subfamily, Platy-1-8c_Saguinus, being the primary source of the Saguinus amplification burst. This constitutes the largest expansion of Platy-1 MEIs reported to date and the most extensive independent SINE amplification between two closely related species.

Precise characterization of somatic complex structural variations from tumor/control paired long-read sequencing data with nanomonsv.

We present our novel software, nanomonsv, for detecting somatic structural variations (SVs) using tumor and matched control long-read sequencing data with a single-base resolution. The current version of nanomonsv includes two detection modules, Canonical SV module, and Single breakend SV module. Using tumor/control paired long-read sequencing data from three cancer and their matched lymphoblastoid lines, we demonstrate that Canonical SV module can identify somatic SVs that can be captured by short-read technologies with higher precision and recall than existing methods. In addition, we have developed a workflow to classify mobile element insertions while elucidating their in-depth properties, such as 5' truncations, internal inversions, as well as source sites for 3' transductions. Furthermore, Single breakend SV module enables the detection of complex SVs that can only be identified by long-reads, such as SVs involving highly-repetitive centromeric sequences, and LINE1- and virus-mediated rearrangements. In summary, our approaches applied to cancer long-read sequencing data can reveal various features of somatic SVs and will lead to a better understanding of mutational processes and functional consequences of somatic SVs.

Analysis of an IncR plasmid carried by carbapenem-resistant Klebsiella pneumoniae: A survey of swine Klebsiella pneumoniae in Jilin Province

This study was conducted in Jilin Province to investigate the mechanism involved in the antibiotic resistance and pathogenicity of Klebsiella pneumoniae. Lung samples were collected from large-scale pig farms in Jilin Province. Antimicrobial susceptibility and mouse lethality assays were carried out. K. pneumoniae isolate JP20, with high virulence and antibiotic resistance, was chosen for whole-genome sequencing. The complete sequence of its genome was annotated, and the virulence and antibiotic resistance mechanism were analysed. A total of 32 K. pneumoniae strains were isolated and tested for antibiotic resistance and pathogenicity. Among them, the JP20 strain showed high levels of resistance to all tested antimicrobial agents and strong pathogenicity in mice (lethal dose of 1.35×1011 CFU/mL). Sequencing of the multidrug-resistant and highly virulent K. pneumoniae JP20 strain revealed that the antibiotic resistance genes were mainly carried by an IncR plasmid. We speculate that extended-spectrum β-lactamases and loss of outer membrane porin OmpK36 play an important role in carbapenem antibiotic resistance. This plasmid contains a mosaic structure consisting of a large number of mobile elements. Through genome-wide analysis, we found that an lncR plasmid carried by the JP20 strain may have evolved in pig farms, possibly leading to multidrug resistance in the JP20 strain. It is speculated that the antibiotic resistance of K. pneumoniae in pig farms is mainly mediated by mobile elements (insertion sequences, transposons, and plasmids). These data provide a basis for monitoring the antibiotic resistance of K. pneumoniae and lay a foundation for an improved understanding of the genomic characteristics and antibiotic resistance mechanism of K. pneumoniae.

Detection of trait-associated structural variations using short-read sequencing

Genomic structural variation (SV) affects genetic and phenotypic characteristics in diverse organisms, but the lack of reliable methods to detect SV has hindered genetic analysis. We developed a computational algorithm (MOPline) that includes missing call recovery combined with high-confidence SV call selection and genotyping using short-read whole-genome sequencing (WGS) data. Using 3,672 high-coverage WGS datasets, MOPline stably detected ∼16,000 SVs per individual, which is over ∼1.7-3.3-fold higher than previous large-scale projects while exhibiting a comparable level of statistical quality metrics. We imputed SVs from 181,622 Japanese individuals for 42 diseases and 60 quantitative traits. A genome-wide association study with the imputed SVs revealed 41 top-ranked or nearly top-ranked genome-wide significant SVs, including 8 exonic SVs with 5 novel associations and enriched mobile element insertions. This study demonstrates that short-read WGS data can be used to identify rare and common SVs associated with a variety of traits.