The Recurrent Atypical e8a2 BCR::ABL1 Transcript with Insertion of an Inverted 55 Base Pair ABL1 Intron 1b Sequence: A Detailed Molecular Analysis

Atypical BCR::ABL1 transcripts are found in approximately 2% of cases of chronic myeloid leukemia. It is important to detect them since affected patients also benefit from tyrosine kinase inhibitor therapy. In the rare e8a2 atypical BCR::ABL1 transcript, two out-of-frame exons are fused, thus, interposed nucleotides are usually found at the fusion site to restore the reading frame. In approximately half of previously reported e8a2 BCR::ABL1 cases, an inserted 55 bp sequence homologous to an inverted sequence from ABL1 intron 1b was detected. The generation of this recurrent transcript variant is not obvious. This work describes the molecular analysis of such an e8a2 BCR::ABL1 translocation from a CML patient. The genomic chromosomal breakpoint is identified, and the formation of this transcript is theoretically explained. The clinical course of the patient is reported, and recommendations are provided for the molecular analysis of future e8a2 BCR::ABL1 cases.

Lineage-Specific Variation in IR Boundary Shift Events, Inversions, and Substitution Rates among Caprifoliaceae s.l. (Dipsacales) Plastomes

Caprifoliaceae s.l. plastid genomes (plastomes) show that one inversion and two inverted repeat boundary shifts occurred in the common ancestor of this family, after which the plastomes are generally conserved. This study reports plastome sequences of five additional species, Fedia cornucopiae, Valeriana fauriei, and Valerianella locusta from the subfamily Valerianoideae, as well as Dipsacus japonicus and Scabiosa comosa from the subfamily Dipsacoideae. Combined with the published plastomes, these plastomes provide new insights into the structural evolution of plastomes within the family. Moreover, the three plastomes from the subfamily Valerianoideae exhibited accelerated nucleotide substitution rates, particularly at synonymous sites, across the family. The patterns of accD sequence divergence in the family are dynamic with structural changes, including interruption of the conserved domain and increases in nonsynonymous substitution rates. In particular, the Valeriana accD gene harbors a large insertion of amino acid repeat (AAR) motifs, and intraspecific polymorphism with a variable number of AARs in the Valeriana accD gene was detected. We found a correlation between intron losses and increased ratios of nonsynonymous to synonymous substitution rates in the clpP gene with intensified positive selection. In addition, two Dipsacoideae plastomes revealed the loss of the plastid-encoded rps15, and a potential functional gene transfer to the nucleus was confirmed.

Detection of a mosaic CDKL5 deletion and inversion by optical genome mapping ends an exhaustive diagnostic odyssey

BACKGROUND

Currently available structural variant (SV) detection methods do not span the complete spectrum of disease-causing SVs. Optical genome mapping (OGM), an emerging technology with the potential to resolve diagnostic dilemmas, was performed to investigate clinically-relevant SVs in a 4-year-old male with an epileptic encephalopathy of undiagnosed molecular origin.

METHODS

OGM was utilized to image long, megabase-size DNA molecules, fluorescently labeled at specific sequence motifs throughout the genome with high sensitivity for detection of SVs greater than 500 bp in size. OGM results were confirmed in a CLIA-certified laboratory via mate-pair sequencing.

RESULTS

OGM identified a mosaic, de novo 90 kb deletion and inversion on the X chromosome disrupting the CDKL5 gene. Detection of the mosaic deletion, which had been previously undetected by chromosomal microarray, an infantile epilepsy panel including exon-level microarray for CDKL5, exome sequencing as well as genome sequencing, resulted in a diagnosis of X-linked dominant early infantile epileptic encephalopathy-2.

CONCLUSION

OGM affords an effective technology for the detection of SVs, especially those that are mosaic, since these remain difficult to detect with current NGS technologies and with conventional chromosomal microarrays. Further research in undiagnosed populations with OGM is warranted.

Comparative genome analysis revealed gene inversions, boundary expansions and contractions, and gene loss in the Stemona sessilifolia (Miq.) Miq. chloroplast genome

Stemona sessilifolia (Miq.) Miq., commonly known as Baibu, is one of the most popular herbal medicines in Asia. In the Chinese Pharmacopoeia, Baibu has multiple authentic sources and there are many similar herbs sold as Baibu in herbal medicine markets. The existence of counterfeits of Baibu brings challenges to its identification. To assist in its accurate identification, we sequenced and analyzed the complete chloroplast genome of S. sessilifolia using next-generation sequencing technology. The genome was found to be 154,037 bp in length, possessing a typical quadripartite structure consisting of a pair of inverted repeats (IRs: 27,090 bp) separated by a large single copy (LSC: 81,949 bp) and a small single copy (SSC: 17,908 bp). A total of 112 unique genes were identified, including 80 protein-coding, 28 transfer RNA and four ribosomal RNA genes. In addition, 45 tandem, 27 forward, 23 palindromic and 104 simple sequence repeats were detected in the genome by repeated analysis. Compared with its counterfeits (Asparagus officinalis and Carludovica palmata) we found that IR expansion and SSC contraction events of S. sessilifolia resulted in two copies of the rpl22 gene in the IR regions and a partial duplication of the ndhF gene in the SSC region. An approximately 3-kb-long inversion was also identified in the LSC region, leading to the petA and cemA genes being presented in the complementary strand of the chloroplast DNA molecule. Comparative analysis revealed some highly variable regions, including trnF-GAA_ndhJ, atpB_rbcL, rps15_ycf1, trnG-UCC_trnR-UCU, ndhF_rpl32, accD_psaI, rps2_rpoC2, trnS-GCU_trnG-UCC, trnT-UGU_trnL-UAA and rps16_trnQ-UUG. Finally, gene loss events were investigated in the context of phylogenetic relationships. In summary, the complete plastome of S. sessilifolia will provide valuable information for the distinction between Baibu and its counterfeits and assist in elucidating the evolution of S. sessilifolia.

Machine-learning predicts genomic determinants of meiosis-driven structural variation in a eukaryotic pathogen

Species harbor extensive structural variation underpinning recent adaptive evolution. However, the causality between genomic features and the induction of new rearrangements is poorly established. Here, we analyze a global set of telomere-to-telomere genome assemblies of a fungal pathogen of wheat to establish a nucleotide-level map of structural variation. We show that the recent emergence of pesticide resistance has been disproportionally driven by rearrangements. We use machine learning to train a model on structural variation events based on 30 chromosomal sequence features. We show that base composition and gene density are the major determinants of structural variation. Retrotransposons explain most inversion, indel and duplication events. We apply our model to Arabidopsis thaliana and show that our approach extends to more complex genomes. Finally, we analyze complete genomes of haploid offspring in a four-generation pedigree. Meiotic crossover locations are enriched for new rearrangements consistent with crossovers being mutational hotspots. The model trained on species-wide structural variation accurately predicts the position of >74% of newly generated variants along the pedigree. The predictive power highlights causality between specific sequence features and the induction of chromosomal rearrangements. Our work demonstrates that training sequence-derived models can accurately identify regions of intrinsic DNA instability in eukaryotic genomes.

Haplotype-resolved diverse human genomes and integrated analysis of structural variation

Long-read and strand-specific sequencing technologies together facilitate the de novo assembly of high-quality haplotype-resolved human genomes without parent-child trio data. We present 64 assembled haplotypes from 32 diverse human genomes. These highly contiguous haplotype assemblies (average minimum contig length needed to cover 50% of the genome: 26 million base pairs) integrate all forms of genetic variation, even across complex loci. We identified 107,590 structural variants (SVs), of which 68% were not discovered with short-read sequencing, and 278 SV hotspots (spanning megabases of gene-rich sequence). We characterized 130 of the most active mobile element source elements and found that 63% of all SVs arise through homology-mediated mechanisms. This resource enables reliable graph-based genotyping from short reads of up to 50,340 SVs, resulting in the identification of 1526 expression quantitative trait loci as well as SV candidates for adaptive selection within the human population.

Functional differences in the hypothalamic-pituitary-gonadal axis are associated with alternative reproductive tactics based on an inversion polymorphism

The evolution of social behavior depends on genetic changes, yet, how genomic variation manifests itself in behavioral diversity is still largely unresolved. Chromosomal inversions can play a pivotal role in producing distinct behavioral phenotypes, in particular, when inversion genes are functionally associated with hormone synthesis and signaling. Male ruffs exhibit alternative reproductive tactics (ARTs) with an autosomal inversion determining two alternative morphs with clear behavioral and hormonal differences from the ancestral morph. We investigated hormonal and transcriptomic differences in the pituitary and gonads. Using a GnRH challenge, we found that the ability to synthesize testosterone in inversion carriers is severely constrained, whereas the synthesis of androstenedione, a testosterone precursor, is not. Inversion morphs were able to produce a transient increase in androstenedione following the GnRH injection, supporting the view that pituitary sensitivity to GnRH is comparable to that of the ancestral morph. We then performed gene expression analyses in a second set of untreated birds and found no evidence of alterations to pituitary sensitivity, gonadotropin production or gonad sensitivity to luteinizing hormone or follicle-stimulating hormone across morphs. Inversion morphs also showed reduced progesterone receptor expression in the pituitary. Strikingly, in the gonads, inversion morphs over-expressed STAR, a gene that is located outside of the inversion and responsible for providing the cholesterol substrate required for the synthesis of sex hormones. In conclusion, our results suggest that the gonads determine morph-specific differences in hormonal regulation.

CRISPR-Cas9-mediated 75.5-Mb inversion in maize

CRISPR-Cas is a powerful double-strand-break technology with wide-ranging applications from gene discovery to commercial product development. Thus far, this tool has been almost exclusively used for gene knockouts and deletions, with a few examples of gene edits and targeted gene insertions. Here, we demonstrate the application of CRISPR-Cas9 technology to mediate targeted 75.5-Mb pericentric inversion in chromosome 2 in one of the elite maize inbred lines from Corteva Agriscience. This inversion unlocks a large chromosomal region containing substantial genetic variance for recombination, thus providing opportunities for the development of new maize varieties with improved phenotypes.

The barley pan-genome reveals the hidden legacy of mutation breeding

Genetic diversity is key to crop improvement. Owing to pervasive genomic structural variation, a single reference genome assembly cannot capture the full complement of sequence diversity of a crop species (known as the 'pan-genome'1). Multiple high-quality sequence assemblies are an indispensable component of a pan-genome infrastructure. Barley (Hordeum vulgare L.) is an important cereal crop with a long history of cultivation that is adapted to a wide range of agro-climatic conditions2. Here we report the construction of chromosome-scale sequence assemblies for the genotypes of 20 varieties of barley-comprising landraces, cultivars and a wild barley-that were selected as representatives of global barley diversity. We catalogued genomic presence/absence variants and explored the use of structural variants for quantitative genetic analysis through whole-genome shotgun sequencing of 300 gene bank accessions. We discovered abundant large inversion polymorphisms and analysed in detail two inversions that are frequently found in current elite barley germplasm; one is probably the product of mutation breeding and the other is tightly linked to a locus that is involved in the expansion of geographical range. This first-generation barley pan-genome makes previously hidden genetic variation accessible to genetic studies and breeding.

Single-cell analysis of structural variations and complex rearrangements with tri-channel processing

Structural variation (SV), involving deletions, duplications, inversions and translocations of DNA segments, is a major source of genetic variability in somatic cells and can dysregulate cancer-related pathways. However, discovering somatic SVs in single cells has been challenging, with copy-number-neutral and complex variants typically escaping detection. Here we describe single-cell tri-channel processing (scTRIP), a computational framework that integrates read depth, template strand and haplotype phase to comprehensively discover SVs in individual cells. We surveyed SV landscapes of 565 single cells, including transformed epithelial cells and patient-derived leukemic samples, to discover abundant SV classes, including inversions, translocations and complex DNA rearrangements. Analysis of the leukemic samples revealed four times more somatic SVs than cytogenetic karyotyping, submicroscopic copy-number alterations, oncogenic copy-neutral rearrangements and a subclonal chromothripsis event. Advancing current methods, single-cell tri-channel processing can directly measure SV mutational processes in individual cells, such as breakage-fusion-bridge cycles, facilitating studies of clonal evolution, genetic mosaicism and SV formation mechanisms, which could improve disease classification for precision medicine.

Construction and comparison of three reference-quality genome assemblies for soybean

We report reference-quality genome assemblies and annotations for two accessions of soybean (Glycine max) and for one accession of Glycine soja, the closest wild relative of G. max. The G. max assemblies provided are for widely used US cultivars: the northern line Williams 82 (Wm82) and the southern line Lee. The Wm82 assembly improves the prior published assembly, and the Lee and G. soja assemblies are new for these accessions. Comparisons among the three accessions show generally high structural conservation, but nucleotide difference of 1.7 single-nucleotide polymorphisms (snps) per kb between Wm82 and Lee, and 4.7 snps per kb between these lines and G. soja. snp distributions and comparisons with genotypes of the Lee and Wm82 parents highlight patterns of introgression and haplotype structure. Comparisons against the US germplasm collection show placement of the sequenced accessions relative to global soybean diversity. Analysis of a pan-gene collection shows generally high conservation, with variation occurring primarily in genomically clustered gene families. We found approximately 40-42 inversions per chromosome between either Lee or Wm82v4 and G. soja, and approximately 32 inversions per chromosome between Wm82 and Lee. We also investigated five domestication loci. For each locus, we found two different alleles with functional differences between G. soja and the two domesticated accessions. The genome assemblies for multiple cultivated accessions and for the closest wild ancestor of soybean provides a valuable set of resources for identifying causal variants that underlie traits for the domestication and improvement of soybean, serving as a basis for future research and crop improvement efforts for this important crop species.

CRISPR/Cas9-Mediated Knockout and In Situ Inversion of the ORF57 Gene from All Copies of the Kaposi's Sarcoma-Associated Herpesvirus Genome in BCBL-1 Cells

Kaposi's sarcoma-associated herpesvirus (KSHV)-transformed primary effusion lymphoma cell lines contain ∼70 to 150 copies of episomal KSHV genomes per cell and have been widely used for studying the mechanisms of KSHV latency and lytic reactivation. Here, we report the first complete knockout (KO) of viral ORF57 gene from all ∼100 copies of KSHV genome per cell in BCBL-1 cells. This was achieved by a modified CRISPR/Cas9 technology to simultaneously express two guide RNAs (gRNAs) and Cas9 from a single expression vector in transfected cells in combination with multiple rounds of cell selection and single-cell cloning. CRISPR/Cas9-mediated genome engineering induces the targeted gene deletion and inversion in situ We found the inverted ORF57 gene in the targeted site in the KSHV genome in one of two characterized single cell clones. Knockout of ORF57 from the KSHV genome led to viral genome instability, thereby reducing viral genome copies and expression of viral lytic genes in BCBL-1-derived single-cell clones. The modified CRISPR/Cas9 technology was very efficient in knocking out the ORF57 gene in iSLK/Bac16 and HEK293/Bac36 cells, where each cell contains only a few copies of the KSHV genome. The ORF57 KO genome was stable in iSLK/Bac16 cells, and, upon lytic induction, was partially rescued by ectopic ORF57 to express viral lytic gene ORF59 and produce infectious virions. Together, the technology developed in this study has paved the way to express two separate gRNAs and the Cas9 enzyme simultaneously in the same cell and could be efficiently applied to any genetic alterations from various genomes, including those in extreme high copy numbers.IMPORTANCE This study provides the first evidence that CRISPR/Cas9 technology can be applied to knock out the ORF57 gene from all ∼100 copies of the KSHV genome in primary effusion lymphoma (PEL) cells by coexpressing two guide RNAs (gRNAs) and Cas9 from a single expression vector in combination with single-cell cloning. The gene knockout efficiency in this system was evaluated rapidly using a direct cell PCR screening. The current CRISPR/Cas9 technology also mediated ORF57 inversion in situ in the targeted site of the KSHV genome. The successful rescue of viral lytic gene expression and infectious virion production from the ORF57 knockout (KO) genome further reiterates the essential role of ORF57 in KSHV infection and multiplication. This modified technology should be useful for knocking out any viral genes from a genome to dissect functions of individual viral genes in the context of the virus genome and to understand their contributions to viral genetics and the virus life cycle.

scoreInvHap: Inversion genotyping for genome-wide association studies

Polymorphic inversions contribute to adaptation and phenotypic variation. However, large multi-centric association studies of inversions remain challenging. We present scoreInvHap, a method to genotype inversions from SNP data for genome-wide association studies (GWASs), overcoming important limitations of current methods and outperforming them in accuracy and applicability. scoreInvHap calls individual inversion-genotypes from a similarity score to the SNPs of experimentally validated references. It can be used on different sources of SNP data, including those with low SNP coverage such as exome sequencing, and is easily adaptable to genotype new inversions, either in humans or in other species. We present 20 human inversions that can be reliably and easily genotyped with scoreInvHap to discover their role in complex human traits, and illustrate a first genome-wide association study of experimentally-validated human inversions. scoreInvHap is implemented in R and it is freely available from Bioconductor.

Chromosomal inversions are structural variants consisting on an orientation change of a chromosome segment. Inversions have been linked to some phenotypic differences between individuals and to genetic divergence. However, their overall contribution to complex diseases is largely underdetermined as there are no high-throughput methods to call inversion-genotypes in large cohort studies. Here, we propose a new method, scoreInvHap, to call individual inversion genotypes from their haplotype similarity. We show that scoreInvHap has a high performance when analyzing heterogeneous sources of SNP data. Our current implementation contains 20 human inversions that can be readily genotyped in existing GWAS datasets. We exemplify the utility of scoreInvHap by running the first-genome wide association of experimentally validated inversions and a multi-centric inversion association study. All in all, scoreInvHap can substantially contribute to increase our knowledge of the role of chromosomal inversions in complex diseases by re-analyzing data from existing genetic association studies.

Mobile Elements Shape Plastome Evolution in Ferns

Plastid genomes display remarkable organizational stability over evolutionary time. From green algae to angiosperms, most plastid genomes are largely collinear, with only a few cases of inversion, gene loss, or, in extremely rare cases, gene addition. These plastome insertions are mostly clade-specific and are typically of nuclear or mitochondrial origin. Here, we expand on these findings and present the first family-level survey of plastome evolution in ferns, revealing a novel suite of dynamic mobile elements. Comparative plastome analyses of the Pteridaceae expose several mobile open reading frames that vary in sequence length, insertion site, and configuration among sampled taxa. Even between close relatives, the presence and location of these elements is widely variable when viewed in a phylogenetic context. We characterize these elements and refer to them collectively as Mobile Open Reading Frames in Fern Organelles (MORFFO). We further note that the presence of MORFFO is not restricted to Pteridaceae, but is found across ferns and other plant clades. MORFFO elements are regularly associated with inversions, intergenic expansions, and changes to the inverted repeats. They likewise appear to be present in mitochondrial and nuclear genomes of ferns, indicating that they can move between genomic compartments with relative ease. The origins and functions of these mobile elements are unknown, but MORFFO appears to be a major driver of structural genome evolution in the plastomes of ferns, and possibly other groups of plants.

Precise and Predictable CRISPR Chromosomal Rearrangements Reveal Principles of Cas9-Mediated Nucleotide Insertion

Chromosomal rearrangements including large DNA-fragment inversions, deletions, and duplications by Cas9 with paired sgRNAs are important to investigate genome structural variations and developmental gene regulation, but little is known about the underlying mechanisms. Here, we report that disrupting CtIP or FANCD2, which have roles in alternative non-homologous end joining, enhances precise DNA-fragment deletion. By analyzing the inserted nucleotides at the junctions of DNA-fragment editing of deletions, inversions, and duplications and characterizing the cleaved products, we find that Cas9 endonucleolytically cleaves the noncomplementary strand with a flexible scissile profile upstream of the -3 position of the PAM site in vivo and in vitro, generating double-strand break ends with 5' overhangs of 1-3 nucleotides. Moreover, we find that engineered Cas9 nucleases have distinct cleavage profiles. Finally, Cas9-mediated nucleotide insertions are nonrandom and are equal to the combined sequences upstream of both PAM sites with predicted frequencies. Thus, precise and predictable DNA-fragment editing could be achieved by perturbing DNA repair genes and using appropriate PAM configurations.

Multiple large inversions and breakpoint rewiring of gene expression in the evolution of the fire ant social supergene

Supergenes consist of co-adapted loci that segregate together and are associated with adaptive traits. In the fire ant Solenopsis invicta, two 'social' supergene variants regulate differences in colony queen number and other traits. Suppressed recombination in this system is maintained, in part, by a greater than 9 Mb inversion, but the supergene is larger. Has the supergene in S. invicta undergone multiple large inversions? The initial gene content of the inverted allele of a supergene would be the same as that of the wild-type allele. So, how did the inversion increase in frequency? To address these questions, we cloned one extreme breakpoint in the fire ant supergene. In doing so, we found a second large (greater than 800 Kb) rearrangement. Furthermore, we determined the temporal order of the two big inversions based on the translocation pattern of a third small fragment. Because the S. invicta supergene lacks evolutionary strata, our finding of multiple inversions may support an introgression model of the supergene. Finally, we showed that one of the inversions swapped the promoter of a breakpoint-adjacent gene, which might have conferred a selective advantage relative to the non-inverted allele. Our findings provide a rare example of gene alterations arising directly from an inversion event.

Detecting Small Inversions Using SRinversion

Rapid development of next generation sequencing (NGS) technology has substantially improved our ability to detect genomic variations. However, unlike other variations, such as point mutations, insertions, and deletions, which can be identified in high sensitivities and specificities based on NGS reads, most of inversions, especially those shorter than 1 kb, remain difficult to detect. Here we introduce a new framework, SRinversion, which was developed specifically for detection of inversions shorter than 1 kb by splitting and realigning poorly mapped or unmapped reads of the NGS data.

Noninvasive detection of F8 int22h-related inversions and sequence variants in maternal plasma of hemophilia carriers

Direct detection of F8 and F9 sequence variants in maternal plasma of hemophilia carriers has been demonstrated by microfluidics digital PCR. Noninvasive prenatal assessment of the most clinically relevant group of sequence variants among patients with hemophilia, namely, those involving int22h-related inversions disrupting the F8 gene, poses additional challenges because of its molecular complexity. We investigated the use of droplet digital PCR (ddPCR) and targeted massively parallel sequencing (MPS) for maternal plasma DNA analysis to noninvasively determine fetal mutational status in pregnancies at risk for hemophilia. We designed family-specific ddPCR assays to detect causative sequence variants scattered across the F8 and F9 genes. A haplotype-based approach coupled with targeted MPS was applied to deduce fetal genotype by capturing a 7.6-Mb region spanning the F8 gene in carriers with int22h-related inversions. The ddPCR analysis correctly determined fetal hemophilia status in 15 at-risk pregnancies in samples obtained from 8 to 42 weeks of gestation. There were 3 unclassified samples, but no misclassification. Detailed fetal haplotype maps of the F8 gene region involving int22h-related inversions obtained through targeted MPS enabled correct diagnoses of fetal mutational status in 3 hemophilia families. Our data suggest it is feasible to apply targeted MPS to interrogate maternally inherited F8 int22h-related inversions, whereas ddPCR represents an affordable approach for the identification of F8 and F9 sequence variants in maternal plasma. These advancements may bring benefits for the pregnancy management for carriers of hemophilia sequence variants; in particular, the common F8 int22h-related inversions, associated with the most severe clinical phenotype.

Expression, purification and crystallization of a protein resulting from the inversion of the amino-acid sequence of a helical bundle

Earlier studies have found that the occurrence of inverse sequence identity in proteins is not indicative of three-dimensional similarity, but rather leads to different folds or unfolded proteins. Short helices, however, frequently keep their conformations when their sequences are inverted. To explore the impact of sequence inversion on long helices, revRM6, with the inverse amino-acid sequence relative to RM6, a highly stable variant of the ColE1 Rop protein, was engineered. RM6 is a highly regular four-α-helical bundle that serves as a model system for protein-folding studies. Here, the crystallization and preliminary crystallographic characterization of revRM6 are reported. The protein was overexpressed in Escherichia coli, purified to homogeneity and crystallized. The crystals belonged to space group P41212, with unit-cell parameters a = b = 44.98, c = 159.74 Å, and diffracted to a resolution of 3.45 Å.

The Structural Features of Thousands of T-DNA Insertion Sites Are Consistent with a Double-Strand Break Repair-Based Insertion Mechanism

Transformation by Agrobacterium tumefaciens, an important tool in modern plant research, involves the integration of T-DNA initially present on a plasmid in agrobacteria into the genome of plant cells. The process of attachment of the agrobacteria to plant cells and the transport of T-DNA into the cell and further to the nucleus has been well described. However, the exact mechanism of integration into the host's DNA is still unclear, although several models have been proposed. During confirmation of T-DNA insertion alleles from the GABI-Kat collection of Arabidopsis thaliana mutants, we have generated about 34,000 sequences from the junctions between inserted T-DNA and adjacent genome regions. Here, we describe the evaluation of this dataset with regard to existing models for T-DNA integration. The results suggest that integration into the plant genome is mainly mediated by the endogenous plant DNA repair machinery. The observed integration events showed characteristics highly similar to those of repair sites of double-strand breaks with respect to microhomology and deletion sizes. In addition, we describe unexpected integration events, such as large deletions and inversions at the integration site that are relevant for correct interpretation of results from T-DNA insertion mutants in reverse genetics experiments.

Formation of linear amplicons with inverted duplications in Leishmania requires the MRE11 nuclease

Extrachromosomal DNA amplification is frequent in the protozoan parasite Leishmania selected for drug resistance. The extrachromosomal amplified DNA is either circular or linear, and is formed at the level of direct or inverted homologous repeated sequences that abound in the Leishmania genome. The RAD51 recombinase plays an important role in circular amplicons formation, but the mechanism by which linear amplicons are formed is unknown. We hypothesized that the Leishmania infantum DNA repair protein MRE11 is required for linear amplicons following rearrangements at the level of inverted repeats. The purified LiMRE11 protein showed both DNA binding and exonuclease activities. Inactivation of the LiMRE11 gene led to parasites with enhanced sensitivity to DNA damaging agents. The MRE11^−/− parasites had a reduced capacity to form linear amplicons after drug selection, and the reintroduction of an MRE11 allele led to parasites regaining their capacity to generate linear amplicons, but only when MRE11 had an active nuclease activity. These results highlight a novel MRE11-dependent pathway used by Leishmania to amplify portions of its genome to respond to a changing environment.

Extrachromosomal DNA amplification is frequent in the human protozoan parasite Leishmania when challenged with drug or other stressful conditions. DNA amplicons, either circular or linear, are formed by recombination between direct or inverted repeats spread throughout the genome of the parasite. The recombinase RAD51 is involved in the formation of circular amplicons, but the mechanism by which linear amplicons are formed is still unknown in this parasite. Studies in other organisms have provided some evidence that a DNA break is required for linear amplifications, and that the DNA repair protein MRE11 can be involved in this process. In this work, we present our biochemical, cellular and molecular characterization of the Leishmania infantum MRE11 orthologue and provide evidence that this nuclease is involved in the formation of linear amplicons in Leishmania. Our results highlight a novel MRE11-dependent pathway used by Leishmania to amplify portions of its genome to respond to a changing environment.

Massive habitat-specific genomic response in D. melanogaster populations during experimental evolution in hot and cold environments

Experimental evolution in combination with whole-genome sequencing (evolve and resequence [E&R]) is a promising approach to define the genotype-phenotype map and to understand adaptation in evolving populations. Many previous studies have identified a large number of putative selected sites (i.e., candidate loci), but it remains unclear to what extent these loci are genuine targets of selection or experimental noise. To address this question, we exposed the same founder population to two different selection regimes-a hot environment and a cold environment-and quantified the genomic response in each. We detected large numbers of putative selected loci in both environments, albeit with little overlap between the two sets of candidates, indicating that most resulted from habitat-specific selection. By quantifying changes across multiple independent biological replicates, we demonstrate that most of the candidate SNPs were false positives that were linked to selected sites over distances much larger than the typical linkage disequilibrium range of Drosophila melanogaster. We show that many of these mid- to long-range associations were attributable to large segregating inversions and confirm by computer simulations that such patterns could be readily replicated when strong selection acts on rare haplotypes. In light of our findings, we outline recommendations to improve the performance of future Drosophila E&R studies which include using species with negligible inversion loads, such as D. mauritiana and D. simulans, instead of D. melanogaster.

InvFEST, a database integrating information of polymorphic inversions in the human genome

The newest genomic advances have uncovered an unprecedented degree of structural variation throughout genomes, with great amounts of data accumulating rapidly. Here we introduce InvFEST (http://invfestdb.uab.cat), a database combining multiple sources of information to generate a complete catalogue of non-redundant human polymorphic inversions. Due to the complexity of this type of changes and the underlying high false-positive discovery rate, it is necessary to integrate all the available data to get a reliable estimate of the real number of inversions. InvFEST automatically merges predictions into different inversions, refines the breakpoint locations, and finds associations with genes and segmental duplications. In addition, it includes data on experimental validation, population frequency, functional effects and evolutionary history. All this information is readily accessible through a complete and user-friendly web report for each inversion. In its current version, InvFEST combines information from 34 different studies and contains 1092 candidate inversions, which are categorized based on internal scores and manual curation. Therefore, InvFEST aims to represent the most reliable set of human inversions and become a central repository to share information, guide future studies and contribute to the analysis of the functional and evolutionary impact of inversions on the human genome.

The plastid genome of mycoheterotrophic monocot Petrosavia stellaris exhibits both gene losses and multiple rearrangements

Plastid genomes of nonphotosynthetic plants represent a perfect model for studying evolution under relaxed selection pressure. However, the information on their sequences is still limited. We sequenced and assembled plastid genome of Petrosavia stellaris, a rare mycoheterotrophic monocot plant. After orchids, Petrosavia represents only the second family of nonphotosynthetic monocots to have its plastid genome examined. Several unusual features were found: retention of the ATP synthase genes and rbcL gene; extensive gene order rearrangement despite a relative lack of repeat sequences; an unusually short inverted repeat region that excludes most of the rDNA operon; and a lack of evidence for accelerated sequence evolution. Plastome of photosynthetic relative of P. stellaris, Japonolirion osense, has standard gene order and does not have the predisposition to inversions. Thus, the rearrangements in the P. stellaris plastome are the most likely associated with transition to heterotrophic way of life.

Novel mutational mechanism in the thyroglobulin gene: imperfect DNA inversion as a cause for hereditary hypothyroidism

The objective of this study was to perform genetic analysis in three brothers of Turkish origin born from consanguineus parents and affected by congenital hypothyroidism, goiter and low levels of serum TG. The combination of sequencing of DNA, PCR mapping, quantitative real-time PCR, inverse-PCR (I-PCR), multiplex PCR and bioinformatics analysis were used in order to detect TG mutations. We demonstrated that the three affected siblings are homozygous for a DNA inversion of 16,962bp in the TG gene associated with two deleted regions at both sides of the inversion limits. The inversion region includes the first 9bp of exon 48, 1015bp of intron 47, 191bp of exon 47, 1523bp of intron 46, 135bp of exon 46 and the last 14,089bp of intron 45. The proximal deletion corresponds to 27bp of TG intron 45, while the distal deletion spans the last 230bp of TG exon 48 and the first 588bp of intergenic region downstream TG end. The parents were heterozygous carriers of the complex rearrangement. In conclusion, a novel large imperfect DNA inversion within the TG gene was identified by the strategy of I-PCR. This aberration was not detectable by normal sequencing of the exons and exon/intron boundaries. Remarkably, the finding represents the first description of a TG deficiency disease caused by a DNA inversion.

Generalized bacterial genome editing using mobile group II introns and Cre-lox

Efficient bacterial genetic engineering approaches with broad-host applicability are rare. We combine two systems, mobile group II introns ('targetrons') and Cre/lox, which function efficiently in many different organisms, into a versatile platform we call GETR (Genome Editing via Targetrons and Recombinases). The introns deliver lox sites to specific genomic loci, enabling genomic manipulations. Efficiency is enhanced by adding flexibility to the RNA hairpins formed by the lox sites. We use the system for insertions, deletions, inversions, and one-step cut-and-paste operations. We demonstrate insertion of a 12-kb polyketide synthase operon into the lacZ gene of Escherichia coli, multiple simultaneous and sequential deletions of up to 120 kb in E. coli and Staphylococcus aureus, inversions of up to 1.2 Mb in E. coli and Bacillus subtilis, and one-step cut-and-pastes for translocating 120 kb of genomic sequence to a site 1.5 Mb away. We also demonstrate the simultaneous delivery of lox sites into multiple loci in the Shewanella oneidensis genome. No selectable markers need to be placed in the genome, and the efficiency of Cre-mediated manipulations typically approaches 100%.

A model-based clustering method for genomic structural variant prediction and genotyping using paired-end sequencing data

Structural variation (SV) has been reported to be associated with numerous diseases such as cancer. With the advent of next generation sequencing (NGS) technologies, various types of SV can be potentially identified. We propose a model based clustering approach utilizing a set of features defined for each type of SV events. Our method, termed SVMiner, not only provides a probability score for each candidate, but also predicts the heterozygosity of genomic deletions. Extensive experiments on genome-wide deep sequencing data have demonstrated that SVMiner is robust against the variability of a single cluster feature, and it significantly outperforms several commonly used SV detection programs. SVMiner can be downloaded from http://cbc.case.edu/svminer/.

A possible mechanism for Inv22-related F8 large deletions in severe hemophilia A patients with high responding factor VIII inhibitors

BACKGROUND

Intron 22 inversion (Inv22) of the coagulation factor (F)VIII gene (F8) is a frequent cause of severe hemophilia A. In addition to Inv22, a variety of F8 mutations (1492 unique mutations) causing hemophilia A have been reported, of which 171 involve deletions of over 50 bp (HAMSTeRs database; http://hadb.org.uk/). However, only 10% of these large deletions have been fully characterized at the nucleotide level.

PATIENTS AND METHODS

We investigated gene abnormalities in three unrelated severe hemophilia A patients with high titer FVIII inhibitors. They had previously been shown to carry large deletions of the F8, but the precise gene abnormalities remain to be elucidated.

RESULTS

Inverse shifting-PCR (IS-PCR) Inv22 diagnostic tests revealed that these patients carried either type I or II Inv22. However, they showed a wild-type (WT) pattern in the IS-PCR Inv22 complementary tests. We further analyzed their X chromosomes to account for the puzzling results, and found that they had different centromeric breakpoints in the Inv22 X chromosomes, adjacent to the palindromic regions containing int22h-2 or -3, and their spacer region, respectively. The connections appeared to be shifted towards the telomere of the WT F8 Xq28, resulting in a new telomere with an additional intact int22h copy.

CONCLUSIONS

These gene rearrangements might result from double-strand breaks in the most distal regions of the long arms of the Inv22 X chromosomes, followed by DNA restorations using the WT F8 Xq28 by non-homologous end joining or break-induced replication; thus leading to large F8 deletions in severe hemophilia A patients.

Application of a novel strategy of engineering conditional alleles to a single exon gene, Sox2

BACKGROUND

The Conditional by Inversion (COIN) method for engineering conditional alleles relies on an invertible optimized gene trap-like element, the COIN module, for imparting conditionality. The COIN module contains an optimized 3' splice site-polyadenylation signal pair, but is inserted antisense to the target gene and therefore does not alter transcription, until it is inverted by Cre recombinase. In order to make COIN applicable to all protein-coding genes, the COIN module has been engineered within an artificial intron, enabling insertion into an exon.

METHODOLOGY/PRINCIPAL FINDINGS

Therefore, theoretically, the COIN method should be applicable to single exon genes, and to test this idea we engineered a COIN allele of Sox2. This single exon gene presents additional design challenges, in that its proximal promoter and coding region are entirely contained within a CpG island, and are also spanned by an overlapping transcript, Sox2Ot, which contains mmu-miR1897. Here, we show that despite disruption of the CpG island by the COIN module intron, the COIN allele of Sox2 (Sox2(COIN)) is phenotypically wild type, and also does not interfere with expression of Sox2Ot and miR1897. Furthermore, the inverted COIN allele of Sox2, Sox2(INV) is functionally null, as homozygotes recapitulate the phenotype of Sox2(βgeo/βgeo) mice, a well-characterized Sox2 null. Lastly, the benefit of the eGFP marker embedded in the COIN allele is demonstrated as it mirrors the expression pattern of Sox2.

CONCLUSIONS/SIGNIFICANCE

Our results demonstrate the applicability of the COIN technology as a method of choice for targeting single exon genes.

A single promoter inversion switches Photorhabdus between pathogenic and mutualistic states

Microbial populations stochastically generate variants with strikingly different properties, such as virulence or avirulence and antibiotic tolerance or sensitivity. Photorhabdus luminescens bacteria have a variable life history in which they alternate between pathogens to a wide variety of insects and mutualists to their specific host nematodes. Here, we show that the P. luminescens pathogenic variant (P form) switches to a smaller-cell variant (M form) to initiate mutualism in host nematode intestines. A stochastic promoter inversion causes the switch between the two distinct forms. M-form cells are much smaller (one-seventh the volume), slower growing, and less bioluminescent than P-form cells; they are also avirulent and produce fewer secondary metabolites. Observations of form switching by individual cells in nematodes revealed that the M form persisted in maternal nematode intestines, were the first cells to colonize infective juvenile (IJ) offspring, and then switched to P form in the IJ intestine, which armed these nematodes for the next cycle of insect infection.

A strong deletion bias in nonallelic gene conversion

Gene conversion is the unidirectional transfer of genetic information between orthologous (allelic) or paralogous (nonallelic) genomic segments. Though a number of studies have examined nucleotide replacements, little is known about length difference mutations produced by gene conversion. Here, we investigate insertions and deletions produced by nonallelic gene conversion in 338 Drosophila and 10,149 primate paralogs. Using a direct phylogenetic approach, we identify 179 insertions and 614 deletions in Drosophila paralogs, and 132 insertions and 455 deletions in primate paralogs. Thus, nonallelic gene conversion is strongly deletion-biased in both lineages, with almost 3.5 times as many conversion-induced deletions as insertions. In primates, the deletion bias is considerably stronger for long indels and, in both lineages, the per-site rate of gene conversion is orders of magnitudes higher than that of ordinary mutation. Due to this high rate, deletion-biased nonallelic gene conversion plays a key role in genome size evolution, leading to the cooperative shrinkage and eventual disappearance of selectively neutral paralogs.

Gene conversion is a process whereby a DNA sequence is copied from one segment of the genome (donor) to another (recipient), resulting in the replacement, insertion, or deletion of a DNA sequence in the recipient. This exchange is facilitated by the high sequence similarity of the two segments, which is due to their evolutionary relationship. Here, we study insertions and deletions produced by gene conversion between paralogs, segments related by DNA duplication events. By comparing paralog sequences in multiple species of fruit flies and primates, we find that deletions occur more than three times as frequently as insertions. We also discover that the rate of gene conversion between paralogs is quite high. The deletion bias and high rate of this process causes paralogs to shrink cooperatively and eventually be eliminated from the genome. Because of the abundance of paralogs in animal genomes, this phenomenon can lead to a significant reduction in genome size. Therefore, our finding enhances our understanding of the forces that lead to changes in genome size during evolution.

Eighteen years of molecular genotyping the hemophilia inversion hotspot: from southern blot to inverse shifting-PCR

The factor VIII gene (F8) intron 22 inversion (Inv22) is a paradigmatic duplicon-mediated rearrangement, found in about one half of patients with severe hemophilia A worldwide. The identification of this prevalent cause of hemophilia was delayed for nine years after the F8 characterization in 1984. The aim of this review is to present the wide diversity of practical approaches that have been developed for genotyping the Inv22 (and related int22h rearrangements) since discovery in 1993. The sequence- Southern blot, long distance-PCR and inverse shifting-PCR-for Inv22 genotyping is an interesting example of scientific ingenuity and evolution in order to resolve challenging molecular diagnostic problems.

Evolution of the rpoB-psbZ region in fern plastid genomes: notable structural rearrangements and highly variable intergenic spacers

BACKGROUND

The rpoB-psbZ (BZ) region of some fern plastid genomes (plastomes) has been noted to go through considerable genomic changes. Unraveling its evolutionary dynamics across all fern lineages will lead to clarify the fundamental process shaping fern plastome structure and organization.

RESULTS

A total of 24 fern BZ sequences were investigated with taxon sampling covering all the extant fern orders. We found that: (i) a tree fern Plagiogyria japonica contained a novel gene order that can be generated from either the ancestral Angiopteris type or the derived Adiantum type via a single inversion; (ii) the trnY-trnE intergenic spacer (IGS) of the filmy fern Vandenboschia radicans was expanded 3-fold due to the tandem 27-bp repeats which showed strong sequence similarity with the anticodon domain of trnY; (iii) the trnY-trnE IGSs of two horsetail ferns Equisetum ramosissimum and E. arvense underwent an unprecedented 5-kb long expansion, more than a quarter of which was consisted of a single type of direct repeats also relevant to the trnY anticodon domain; and (iv) ycf66 has independently lost at least four times in ferns.

CONCLUSIONS

Our results provided fresh insights into the evolutionary process of fern BZ regions. The intermediate BZ gene order was not detected, supporting that the Adiantum type was generated by two inversions occurring in pairs. The occurrence of Vandenboschia 27-bp repeats represents the first evidence of partial tRNA gene duplication in fern plastomes. Repeats potentially forming a stem-loop structure play major roles in the expansion of the trnY-trnE IGS.

Insertion sequence inversions mediated by ectopic recombination between terminal inverted repeats

Transposable elements are widely distributed and diverse in both eukaryotes and prokaryotes, as exemplified by DNA transposons. As a result, they represent a considerable source of genomic variation, for example through ectopic (i.e. non-allelic homologous) recombination events between transposable element copies, resulting in genomic rearrangements. Ectopic recombination may also take place between homologous sequences located within transposable element sequences. DNA transposons are typically bounded by terminal inverted repeats (TIRs). Ectopic recombination between TIRs is expected to result in DNA transposon inversions. However, such inversions have barely been documented. In this study, we report natural inversions of the most common prokaryotic DNA transposons: insertion sequences (IS). We identified natural TIR-TIR recombination-mediated inversions in 9% of IS insertion loci investigated in Wolbachia bacteria, which suggests that recombination between IS TIRs may be a quite common, albeit largely overlooked, source of genomic diversity in bacteria. We suggest that inversions may impede IS survival and proliferation in the host genome by altering transpositional activity. They may also alter genomic instability by modulating the outcome of ectopic recombination events between IS copies in various orientations. This study represents the first report of TIR-TIR recombination within bacterial IS elements and it thereby uncovers a novel mechanism of structural variation for this class of prokaryotic transposable elements.

[Mobile genetic elements in the epidemiology of bacterial resistance to antibiotics]

The study of the role of mobile elements and mobilization of resistance genes is crucial for understanding the epidemiology of antibiotic resistance. This review summarizes recent data on the insertion sequences, transposons, integrons and plasmids that are involved in the mobilization of bacterial antibiotic resistance genes.