Mobile Element Insertions Research Articles

The Simons Genome Diversity Project: A Global Analysis of Mobile Element Diversity.

Ongoing retrotransposition of Alu, LINE-1, and SINE–VNTR–Alu elements generates diversity and variation among human populations. Previous analyses investigating the population genetics of mobile element insertions (MEIs) have been limited by population ascertainment bias or by relatively small numbers of populations and low sequencing coverage. Here, we use 296 individuals representing 142 global populations from the Simons Genome Diversity Project (SGDP) to discover and characterize MEI diversity from deeply sequenced whole-genome data. We report 5,742 MEIs not originally reported by the 1000 Genomes Project and show that high sampling diversity leads to a 4- to 7-fold increase in MEI discovery rates over the original 1000 Genomes Project data. As a result of negative selection, nonreference polymorphic MEIs are underrepresented within genes, and MEIs within genes are often found in the transcriptional orientation opposite that of the gene. Globally, 80% of Alu subfamilies predate the expansion of modern humans from Africa. Polymorphic MEIs show heterozygosity gradients that decrease from Africa to Eurasia to the Americas, and the number of MEIs found uniquely in a single individual are also distributed in this general pattern. The maximum fraction of MEI diversity partitioned among the seven major SGDP population groups (FST) is 7.4%, similar to, but slightly lower than, previous estimates and likely attributable to the diverse sampling strategy of the SGDP. Finally, we utilize these MEIs to extrapolate the primary Native American shared ancestry component to back to Asia and provide new evidence from genome-wide identical-by-descent genetic markers that add additional support for a southeastern Siberian origin for most Native Americans.

Mobile element insertion detection in 89,874 clinical exomes

PurposeExome sequencing (ES) is increasingly used for the diagnosis of rare genetic disease. However, some pathogenic sequence variants within the exome go undetected due to the technical difficulty of identifying them. Mobile element insertions (MEIs) are a known cause of genetic disease in humans but have been historically difficult to detect via ES and similar targeted sequencing methods. MethodsWe developed and applied a novel MEI detection method prospectively to samples received for clinical ES beginning in November 2017. Positive MEI findings were confirmed by an orthogonal method and reported back to the ordering provider. In this study, we examined 89,874 samples from 38,871 cases. ResultsDiagnostic MEIs were present in 0.03% (95% binomial test confidence interval: 0.02–0.06%) of all cases and account for 0.15% (95% binomial test confidence interval: 0.08–0.25%) of cases with a molecular diagnosis. One diagnostic MEI was a novel founder event. Most patients with pathogenic MEIs had prior genetic testing, three of whom had previous negative DNA sequencing analysis of the diagnostic gene. ConclusionMEI detection from ES is a valuable diagnostic tool, reveals molecular findings that may be undetected by other sequencing assays, and increases diagnostic yield by 0.15%.

Integrated Mobile Element Scanning (ME-Scan) method for identifying multiple types of polymorphic mobile element insertions

BackgroundMobile elements are ubiquitous components of mammalian genomes and constitute more than half of the human genome. Polymorphic mobile element insertions (pMEIs) are a major source of human genomic variation and are gaining research interest because of their involvement in gene expression regulation, genome integrity, and disease.ResultsBuilding on our previous Mobile Element Scanning (ME-Scan) protocols, we developed an integrated ME-Scan protocol to identify three major active families of human mobile elements, AluYb, L1HS, and SVA. This approach selectively amplifies insertion sites of currently active retrotransposons for Illumina sequencing. By pooling the libraries together, we can identify pMEIs from all three mobile element families in one sequencing run. To demonstrate the utility of the new ME-Scan protocol, we sequenced 12 human parent-offspring trios. Our results showed high sensitivity (> 90%) and accuracy (> 95%) of the protocol for identifying pMEIs in the human genome. In addition, we also tested the feasibility of identifying somatic insertions using the protocol.ConclusionsThe integrated ME-Scan protocol is a cost-effective way to identify novel pMEIs in the human genome. In addition, by developing the protocol to detect three mobile element families, we demonstrate the flexibility of the ME-Scan protocol. We present instructions for the library design, a sequencing protocol, and a computational pipeline for downstream analyses as a complete framework that will allow researchers to easily adapt the ME-Scan protocol to their own projects in other genomes.

TypeTE: a tool to genotype mobile element insertions from whole genome resequencing data.

Alu retrotransposons account for more than 10% of the human genome, and insertions of these elements create structural variants segregating in human populations. Such polymorphic Alus are powerful markers to understand population structure, and they represent variants that can greatly impact genome function, including gene expression. Accurate genotyping of Alus and other mobile elements has been challenging. Indeed, we found that Alu genotypes previously called for the 1000 Genomes Project are sometimes erroneous, which poses significant problems for phasing these insertions with other variants that comprise the haplotype. To ameliorate this issue, we introduce a new pipeline – TypeTE – which genotypes Alu insertions from whole-genome sequencing data. Starting from a list of polymorphic Alus, TypeTE identifies the hallmarks (poly-A tail and target site duplication) and orientation of Alu insertions using local re-assembly to reconstruct presence and absence alleles. Genotype likelihoods are then computed after re-mapping sequencing reads to the reconstructed alleles. Using a high-quality set of PCR-based genotyping of >200 loci, we show that TypeTE improves genotype accuracy from 83% to 92% in the 1000 Genomes dataset. TypeTE can be readily adapted to other retrotransposon families and brings a valuable toolbox addition for population genomics.

Whole genome sequencing reveals the distribution of resistance and virulence genes of pathogenic Escherichia coli CCHTP from giant panda

Pathogenic Escherichia coli (E. coli) is the most common pathogen causing urinary tract infection in animals. We investigated the antibiotic resistance and virulence genes of pathogenic E. coli CCHTP derived from urine with occult blood of the giant panda by whole genome sequencing. The flanking sequencing of resistance and virulence genes in genomic islands were also analyzed. Our results demonstrate that E. coli CCHTP contains different families of antibiotic resistance genes, most of which are efflux pump related genes, including multiple drug resistance efflux pump genes mdfA, emrE, and mdtN. A total of 166 virulence factors and 563 virulence genes were identified, and the most virulence factors and related genes are involved in host cell attachment and invasion processes. Furthermore, sequence analysis of 19 genomic islands revealed that antibiotic and virulence genes are associated with mobile genetic elements (transposon and insertion sequence) in GIs011 and GIs017. These structures can mediate horizontal transfer of antibiotic and virulence genes. Our work described the distribution of antibiotic resistance genes and virulence genes in E. coli CCHTP, which may provide an important guidance for treatment and rational drug use of E. coli CCHTP infection in the giant panda.

Pedigree-based estimation of human mobile element retrotransposition rates.

Germline mutation rates in humans have been estimated for a variety of mutation types, including single-nucleotide and large structural variants. Here, we directly measure the germline retrotransposition rate for the three active retrotransposon elements: L1, Alu, and SVA. We used three tools for calling mobile element insertions (MEIs) (MELT, RUFUS, and TranSurVeyor) on blood-derived whole-genome sequence (WGS) data from 599 CEPH individuals, comprising 33 three-generation pedigrees. We identified 26 de novo MEIs in 437 births. The retrotransposition rate estimates for Alu elements, one in 40 births, is roughly half the rate estimated using phylogenetic analyses, a difference in magnitude similar to that observed for single-nucleotide variants. The L1 retrotransposition rate is one in 63 births and is within range of previous estimates (1:20–1:200 births). The SVA retrotransposition rate, one in 63 births, is much higher than the previous estimate of one in 900 births. Our large, three-generation pedigrees allowed us to assess parent-of-origin effects and the timing of insertion events in either gametogenesis or early embryonic development. We find a statistically significant paternal bias in Alu retrotransposition. Our study represents the first in-depth analysis of the rate and dynamics of human retrotransposition from WGS data in three-generation human pedigrees.

Homozygosity for Mobile Element Insertions Associated with WBSCR17 Could Predict Success in Assistance Dog Training Programs.

Assistance dog training programs can see as many as 60% of their trainees dismissed. Many training programs utilize behavioral assays prior to admittance to identify likely successful candidates, yet such assays can be insconsistent. Recently, four canine retrotransposon mobile element insertions (MEIs) in or near genes WBSCR17 (Cfa6.6 and Cfa6.7), GTF2I (Cfa6.66) and POM121 (Cfa6.83) were identified in domestic dogs and gray wolves. Variations in these MEIs were significantly associated with a heightened propensity to initiate prolonged social contact or hypersociability. Using our dataset of 837 dogs, 228 of which had paired survey-based behavioral data, we discovered that one of the insertions in WBSCR17 is the most important predictor of dog sociable behaviors related to human proximity, measured by the Canine Behavioral Assessment Research Questionnaire (C-BARQ©). We found a positive correlation between insertions at Cfa6.6 and dog separation distress in the form of restlessness when about to be left alone by the owner. Lastly, assistance dogs showed significant heterozygosity deficiency at locus Cfa6.6 and higher frequency of insertions at Cfa6.6 and Cfa6.7. We suggest that training programs could utilize this genetic survey to screen for MEIs at WBSCR17 to identify dogs with sociable traits compatible with successful assistance dog performance.

Mobile genetic elements explain size variation in the mitochondrial genomes of four closely-related Armillaria species

BackgroundSpecies in the genus Armillaria (fungi, basidiomycota) are well-known as saprophytes and pathogens on plants. Many of them cause white-rot root disease in diverse woody plants worldwide. Mitochondrial genomes (mitogenomes) are widely used in evolutionary and population studies, but despite the importance and wide distribution of Armillaria, the complete mitogenomes have not previously been reported for this genus. Meanwhile, the well-supported phylogeny of Armillaria species provides an excellent framework in which to study variation in mitogenomes and how they have evolved over time.ResultsHere we completely sequenced, assembled, and annotated the circular mitogenomes of four species: A. borealis, A. gallica, A. sinapina, and A. solidipes (116,443, 98,896, 103,563, and 122,167 bp, respectively). The variation in mitogenome size can be explained by variable numbers of mobile genetic elements, introns, and plasmid-related sequences. Most Armillaria introns contained open reading frames (ORFs) that are related to homing endonucleases of the LAGLIDADG and GIY-YIG families. Insertions of mobile elements were also evident as fragments of plasmid-related sequences in Armillaria mitogenomes. We also found several truncated gene duplications in all four mitogenomes.ConclusionsOur study showed that fungal mitogenomes have a high degree of variation in size, gene content, and genomic organization even among closely related species of Armillara. We suggest that mobile genetic elements invading introns and intergenic sequences in the Armillaria mitogenomes have played a significant role in shaping their genome structure. The mitogenome changes we describe here are consistent with widely accepted phylogenetic relationships among the four species.

RMETL: sensitive mobile element insertion detection with long read realignment.

Mobile element insertion (MEI) is a major category of structure variations (SVs). The rapid development of long read sequencing technologies provides the opportunity to detect MEIs sensitively. However, the signals of MEI implied by noisy long reads are highly complex due to the repetitiveness of mobile elements as well as the high sequencing error rates. Herein, we propose the Realignment-based Mobile Element insertion detection Tool for Long read (rMETL). Benchmarking results of simulated and real datasets demonstrate that rMETL enables to handle the complex signals to discover MEIs sensitively. It is suited to produce high-quality MEI callsets in many genomics studies. rMETL is available from https://github.com/hitbc/rMETL. Supplementary data are available at Bioinformatics online.

Биоинформатический алгоритм обработки данных высокопроизводительного секвенирования для прецизионной идентификации соматических инсерций ретроэлементов

Retroelements are considered as one of the important sources of genomic variability in modern humans. It is known that transposition activity of retroelements in germline cells generates new insertions in various genomic loci and sometimes results in genetic diseases. Retroelements activity in somatic cells is restricted by different cellular mechanisms; however, there is an evidence for it in some tissue types. Somatic insertions can trigger tumorigenesis or participate in normal functioning such as generation of neurons' plasticity. In spite of the rapid development of high-throughput sequencing methods a confident detection of somatic insertions is still quite a challenging task. That, in part, is due to the absence of adequate bioinformatic tools for the analysis of sequencing data. Here, we propose an advanced computational pipeline for the identification of somatic insertions in datasets generated by selective amplification and high-throughput sequencing of genomic regions flanking insertions of AluYa5. Particular attention is paid for the identification of various artifacts arising in course of library preparation and the parameters for their filtration. Pipeline sensitivity is confirmed by in silico experiments with artificial datasets. Using the proposed pipeline we remove at least 80% of artifacts and preserve 75% of potentially somatic insertions. The approaches used in this work can be applied for the study of other mobile elements insertion variability.

Antibiotics, Resistome and Resistance Mechanisms: A Bacterial Perspective.

History of mankind is regarded as struggle against infectious diseases. Rather than observing the withering away of bacterial diseases, antibiotic resistance has emerged as a serious global health concern. Medium of antibiotic resistance in bacteria varies greatly and comprises of target protection, target substitution, antibiotic detoxification and block of intracellular antibiotic accumulation. Further aggravation to prevailing situation arose on observing bacteria gradually becoming resistant to different classes of antibiotics through acquisition of resistance genes from same and different genera of bacteria. Attributing bacteria with feature of better adaptability, dispersal of antibiotic resistance genes to minimize effects of antibiotics by various means including horizontal gene transfer (conjugation, transformation, and transduction), Mobile genetic elements (plasmids, transposons, insertion sequences, integrons, and integrative-conjugative elements) and bacterial toxin-antitoxin system led to speedy bloom of antibiotic resistance amongst bacteria. Proficiency of bacteria to obtain resistance genes generated an unpleasant situation; a grave, but a lot unacknowledged, feature of resistance gene transfer.

Activity of Genes with Functions in Human Williams-Beuren Syndrome Is Impacted by Mobile Element Insertions in the Gray Wolf Genome.

In canines, transposon dynamics have been associated with a hyper-social behavioral syndrome, although the functional mechanism has yet to be described. We investigate the epigenetic and transcriptional consequences of these behavior-associated mobile element insertions (MEIs) in dogs and Yellowstone gray wolves. We posit that the transposons themselves may not be the causative feature; rather, their transcriptional regulation may exert the functional impact. We survey four outlier transposons associated with hyper-sociability, with the expectation that they are targeted for epigenetic silencing. We predict hyper-methylation of MEIs, suggestive that the epigenetic silencing of and not the MEIs themselves may be driving dysregulation of nearby genes. We found that transposon-derived sequences are significantly hyper-methylated, regardless of their copy number or species. Further, we have assessed transcriptome sequence data and found evidence that MEIs impact the expression levels of six genes (WBSCR17, LIMK1, GTF2I, WBSCR27, BAZ1B, and BCL7B), all of which have known roles in human Williams–Beuren syndrome due to changes in copy number, typically hemizygosity. Although further evidence is needed, our results suggest that a few insertions alter local expression at multiple genes, likely through a cis-regulatory mechanism that excludes proximal methylation.

Co-evolution with recombination affects the stability of mobile genetic element insertions within gene families of Salmonella.

Bacteria can have multiple copies of a gene at separate locations on the same chromosome. Some of these gene families, including tuf (translation elongation factor EF-Tu) and rrl (ribosomal RNA), encode functions critically important for bacterial fitness. Genes within these families are known to evolve in concert using homologous recombination to transfer genetic information from one gene to another. This mechanism can counteract the detrimental effects of nucleotide sequence divergence over time. Whether such mechanisms can also protect against the potentially lethal effects of mobile genetic element insertion is not well understood. To address this we constructed two different length insertion cassettes to mimic mobile genetic elements and inserted these into various positions of the tuf and rrl genes. Wemeasured rates of recombinational repair that removed the inserted cassette and studied the underlying mechanism. Our results indicate that homologous recombination can protect the tuf and rrl genes from inactivation by mobile genetic elements, but forinsertions within shorter gene sequences the efficiency of repair is very low. Intriguingly, we found that physical distance separating genes on the chromosome directly affects the rate of recombinational repair suggesting that relative location will influence the ability of homologous recombination to maintain homogeneity.

Comparative sequence alignment reveals River Buffalo genomic structural differences compared with cattle

This study sought to characterize differences in gene content, regulation and structure between taurine cattle and river buffalo (one subspecies of domestic water buffalo) using the extensively annotated UMD3.1 cattle reference genome as a basis for comparisons. We identified 127 deletion CNV regions in river buffalo representing 5 annotated cattle genes. We also characterized 583 merged mobile element insertion (MEI) events within the upstream regions of annotated cattle genes. Transcriptome analysis in various tissue types on river buffalo confirmed the absence of four cattle genes. Four genes which may be related to phenotypic differences in meat quality and color, had upstream MEI predictions and were found to have significantly elevated expression in river buffalo compared with cattle. Our comparative alignment approach and gene expression analyses suggested a functional role for many genomic structural variations, which may contribute to the unique phenotypes of river buffalo.

Molecular characteristic of stable and unstable white gene alleles in highly mutable lines from natural Drosophila melanogaster populations

Mutations in the white locus emerged in highly mutable isofemale Drosophila melanogaster lines from the populations of Novosibirsk 2013 (NS3 line), Nalchik 2014 (N119 line), and Sakhalin Island 2014 (S46 line). A single white-eyed male found in the NS3 line was sterile. Phenotypically mutant derivatives (white gene alleles) differing in eye color (pure white, different shades of yellow (honey), orange (apricot), cherry, and red (wild type)) emerged during the N119 and S46 line breeding in the laboratory. Molecular genetic study of the structure of wild type white locus in initial lines and white-mutant derivatives de novo emerging from them, as well as other white lines from the fund of the Laboratory of Population Genetics of the Institute of Cytology and Genetics (Siberian Branch, Russian Academy of Sciences), was conducted. The pairs of primers flanking different white gene regions were selected. Six such pairs overlapped the coding part of the gene. Molecular genetic analysis demonstrated that most DNA defects were limited to the region which includes the first exon (34 lines). Among them, four mutant events were accompanied by an insertion of DNA fragments of approximately 800 bp; one mutation event was accompanied by a deletion of approximately 200 bp; in 29 cases, no PCR product was obtained (this can indicate that as a minimum one of the primer binding sites is damaged). The inserted DNA fragments have no homology with known D. melanogaster sequences presented in the NCBI database. The complete white gene deletion with the manifestation of mutant “white eyes” phenotype was registered in four cases (and only in the N119 line derivatives). Normal PCR product was obtained in 22 cases for all six DNA fragments. Among them, there are both alleles phenotypically mutant by the eye color (white, cherry, or orange) and revertants to the wild type (red). The abundance of defects in the beginning of the gene can indicate a multiplicity of mobile genetic element insertion sites in this part of the white gene in D. melanogaster.

The Mobile Element Locator Tool (MELT): population-scale mobile element discovery and biology.

Mobile element insertions (MEIs) represent ∼25% of all structural variants in human genomes. Moreover, when they disrupt genes, MEIs can influence human traits and diseases. Therefore, MEIs should be fully discovered along with other forms of genetic variation in whole genome sequencing (WGS) projects involving population genetics, human diseases, and clinical genomics. Here, we describe the Mobile Element Locator Tool (MELT), which was developed as part of the 1000 Genomes Project to perform MEI discovery on a population scale. Using both Illumina WGS data and simulations, we demonstrate that MELT outperforms existing MEI discovery tools in terms of speed, scalability, specificity, and sensitivity, while also detecting a broader spectrum of MEI-associated features. Several run modes were developed to perform MEI discovery on local and cloud systems. In addition to using MELT to discover MEIs in modern humans as part of the 1000 Genomes Project, we also used it to discover MEIs in chimpanzees and ancient (Neanderthal and Denisovan) hominids. We detected diverse patterns of MEI stratification across these populations that likely were caused by (1) diverse rates of MEI production from source elements, (2) diverse patterns of MEI inheritance, and (3) the introgression of ancient MEIs into modern human genomes. Overall, our study provides the most comprehensive map of MEIs to date spanning chimpanzees, ancient hominids, and modern humans and reveals new aspects of MEI biology in these lineages. We also demonstrate that MELT is a robust platform for MEI discovery and analysis in a variety of experimental settings.

Evolution of ribosomal internal transcribed spacers in Deuterostomia

Sequences of ribosomal internal transcribed spacers (ITSs) are of great importance to molecular phylogenetics and DNA barcoding, but remain unstudied in some large taxa of Deuterostomia. We have analyzed complete ITS1 and ITS2 sequences in 62 species from 16 Deuterostomia classes, with ITS sequences in 24 species from 11 classes initially obtained using unannotated contigs and raw read sequences. A general tendency for both ITS length and GC-content increase from interior to superior Deuterostomia taxa, a uniform GC-content in both ITSs within the same species, thymine content decrease in sense DNA sequences of both ITSs are shown. A possible role of GC-based gene conversion in Deuterostomia ITS evolutionary changes is hypothesized. The first example of non-LTR retrotransposon insertion into ITS sequence in Deuterostomia is described in turtle Geochelone nigra. The roles of mobile genetic element insertions in the evolution of ITS sequences in some Sauropsida taxa are discussed as well.

Evolutionary loss of melanogenesis in the tunicate Molgula occulta

BackgroundAnalyzing close species with diverse developmental modes is instrumental for investigating the evolutionary significance of physiological, anatomical and behavioral features at a molecular level. Many examples of trait loss are known in metazoan populations living in dark environments. Tunicates are the closest living relatives of vertebrates and typically present a lifecycle with distinct motile larval and sessile adult stages. The nervous system of the motile larva contains melanized cells associated with geotactic and light-sensing organs. It has been suggested that these are homologous to vertebrate neural crest-derived melanocytes. Probably due to ecological adaptation to distinct habitats, several species of tunicates in the Molgulidae family have tailless (anural) larvae that fail to develop sensory organ-associated melanocytes. Here we studied the evolution of Tyrosinase family genes, indispensible for melanogenesis, in the anural, unpigmented Molgula occulta and in the tailed, pigmented Molgula oculata by using phylogenetic, developmental and molecular approaches.ResultsWe performed an evolutionary reconstruction of the tunicate Tyrosinase gene family: in particular, we found that M. oculata possesses genes predicted to encode one Tyrosinase (Tyr) and three Tyrosinase-related proteins (Tyrps) while M. occulta has only Tyr and Tyrp.a pseudogenes that are not likely to encode functional proteins. Analysis of Tyr sequences from various M. occulta individuals indicates that different alleles independently acquired frameshifting short indels and/or larger mobile genetic element insertions, resulting in pseudogenization of the Tyr locus. In M. oculata, Tyr is expressed in presumptive pigment cell precursors as in the model tunicate Ciona robusta. Furthermore, a M. oculata Tyr reporter gene construct was active in the pigment cell precursors of C. robusta embryos, hinting at conservation of the regulatory network underlying Tyr expression in tunicates. In contrast, we did not observe any expression of the Tyr pseudogene in M. occulta embryos. Similarly, M. occulta Tyr allele expression was not rescued in pigmented interspecific M. occulta × M. oculata hybrid embryos, suggesting deleterious mutations also to its cis-regulatory sequences. However, in situ hybridization for transcripts from the M. occulta Tyrp.a pseudogene revealed its expression in vestigial pigment cell precursors in this species.ConclusionsWe reveal a complex evolutionary history of the melanogenesis pathway in tunicates, characterized by distinct gene duplication and loss events. Our expression and molecular data support a tight correlation between pseudogenization of Tyrosinase family members and the absence of pigmentation in the immotile larvae of M. occulta. These results suggest that relaxation of purifying selection has resulted in the loss of sensory organ-associated melanocytes and core genes in the melanogenesis biosynthetic pathway in M. occulta.

Quantifying the Number of Independent Organelle DNA Insertions in Genome Evolution and Human Health.

Fragments of organelle genomes are often found as insertions in nuclear DNA. These fragments of mitochondrial DNA (numts) and plastid DNA (nupts) are ubiquitous components of eukaryotic genomes. They are, however, often edited out during the genome assembly process, leading to systematic underestimation of their frequency. Numts and nupts, once inserted, can become further fragmented through subsequent insertion of mobile elements or other recombinational events that disrupt the continuity of the inserted sequence relative to the genuine organelle DNA copy. Because numts and nupts are typically identified through sequence comparison tools such as BLAST, disruption of insertions into smaller fragments can lead to systematic overestimation of numt and nupt frequencies. Accurate identification of numts and nupts is important, however, both for better understanding of their role during evolution, and for monitoring their increasingly evident role in human disease. Human populations are polymorphic for 141 numt loci, five numts are causal to genetic disease, and cancer genomic studies are revealing an abundance of numts associated with tumor progression. Here, we report investigation of salient parameters involved in obtaining accurate estimates of numt and nupt numbers in genome sequence data. Numts and nupts from 44 sequenced eukaryotic genomes reveal lineage-specific differences in the number, relative age and frequency of insertional events as well as lineage-specific dynamics of their postinsertional fragmentation. Our findings outline the main technical parameters influencing accurate identification and frequency estimation of numts in genomic studies pertinent to both evolution and human health.

Characterization and pathogenicity of acute hepatopancreatic necrosis disease natural mutants, pirABvp (‐) V. parahaemolyticus, and pirABvp (+) V. campbellii strains

Two Vibrio parahaemolyticus virulence genes, pirAvp and pirBvp, are known to encode a binary photorhabdus insect-related (Pir) toxin that causes acute hepatopancreatic necrosis disease (AHPND) in shrimp. These genes are flanked with repeats of a mobile element (insertion sequence) in a large plasmid. This insertion sequence is closely (92%) related to the known insertion sequence ISVal1. The pirABvp genes and the flanking ISVal1 forms a 5535-bp composite transposon, designated Tn6264. There are pirABvp gene deletions in some strains of V. parahaemolyticus. During 2013–2016, we found 2 types of pirABvp deletion mutants from AHPND-affected farms. The type I mutants included 3 strains with deletions (4.4-kb or 6.0-kb) of entire pirABvp genes and the downstream ISVal1, and these mutants were named pirABvp(‐). The type II mutants included 3 strains with smaller deletions (1.5-kb or 1.6-kb) including a pirAvp gene and a partial pirBvp gene, and were named pirAvp(‐). In laboratory bioassays, these were not pathogenic to shrimp confirming that both pirAvp and pirBvp are required for AHPND pathogenicity. During 2016, we also isolated 4 V. campbellii strains carrying pirABvp genes from diseased shrimp. These V. campbellii stains were found, through laboratory bioassays and histological evaluation, to cause AHPND.