Large-scale Genomic Rearrangements Research Articles

The chromosome-scale reference genome for the pinfish (Lagodon rhomboides) provides insights into their evolutionary and demographic history.

The pinfish (Lagodon rhomboides) is an ecologically, economically, and culturally relevant member of family Sparidae, playing crucial roles in marine food webs of the western Atlantic Ocean and Gulf of Mexico. Despite their high abundance and ecological importance, there is a scarcity of genomic resources for this species. We assembled and annotated a chromosome-scale genome for the pinfish, resulting in a highly contiguous 785 Mb assembly of 24 scaffolded chromosomes. The high-quality assembly contains 98.9% complete BUSCOs and shows strong synteny to other chromosome-scale genomes of fish in the family Sparidae, with a limited number of large-scale genomic rearrangements. Leveraging this new genomic resource, we found evidence of significant expansions of dietary gene families over the evolutionary history of the pinfish, which may be associated with an ontogenetic shift from carnivory to herbivory seen in this species. Estimates of historical patterns of population demography using this new reference genome identified several periods of population growth and contraction which were associated with ancient climatic shifts and sea level changes. This genome serves as a valuable reference for future studies of population genomics and differentiation and provides a much-needed genomic resource for this western Atlantic sparid.

Experimental evolution of extremotolerant and extremophilic fungi under osmotic stress.

Experimental evolution was carried out to investigate the adaptive responses of extremotolerant fungi to a stressful environment. For 12 cultivation cycles, the halotolerant black yeasts Aureobasidium pullulans and Aureobasidium subglaciale were grown at high NaCl or glycerol concentrations, and the halophilic basidiomycete Wallemia ichthyophaga was grown close to its lower NaCl growth limit. All evolved Aureobasidium spp. accelerated their growth at low water activity. Whole genomes of the evolved strains were sequenced. No aneuploidies were detected in any of the genomes, contrary to previous studies on experimental evolution at high salinity with other species. However, several hundred single-nucleotide polymorphisms were identified compared with the genomes of the progenitor strains. Two functional groups of genes were overrepresented among the genes presumably affected by single-nucleotide polymorphisms: voltage-gated potassium channels in A. pullulans at high NaCl concentration, and hydrophobins in W. ichthyophaga at low NaCl concentration. Both groups of genes were previously associated with adaptation to high salinity. Finally, most evolved Aureobasidium spp. strains were found to have increased intracellular and decreased extracellular glycerol concentrations at high salinity, suggesting that the strains have optimised their management of glycerol, their most important compatible solute. Experimental evolution therefore not only confirmed the role of potassium transport, glycerol management, and cell wall in survival at low water activity, but also demonstrated that fungi from extreme environments can further improve their growth rates under constant extreme conditions in a relatively short time and without large scale genomic rearrangements.

The number and pattern of viral genomic reassortments are not necessarily identifiable from segment trees.

Reassortment is an evolutionary process common in viruses with segmented genomes. These viruses can swap whole genomic segments during cellular co-infection, giving rise to novel progeny formed from the mixture of parental segments. Because large-scale genome rearrangements have the potential to generate new phenotypes, reassortment is important to both evolutionary biology and public health research. However, statistical inference of the pattern of reassortment events from phylogenetic data is exceptionally difficult, potentially involving inference of general graphs in which individual segment trees are embedded. In this paper, we argue that, in general, the number and pattern of reassortment events are not identifiable from segment trees alone, even with theoretically ideal data. We call this fact the fundamental problem of reassortment, which we illustrate using the concept of the `first-infection tree', a typically but not always counterfactual genealogy that would have been observed in the segment trees had no reassortment occurred. Further, we illustrate four additional problems that can arise logically in the inference of reassortment events and show, using simulated data, that these problems are not rare and can potentially distort our perception of reassortment even in small data sets. Finally, we discuss how existing methods can be augmented or adapted to account for not only the fundamental problem of reassortment but also the four additional situations that can complicate the inference of reassortment.

Abstract 405: Translocation detection in cancer using low-pass pore-c sequencing

Abstract Complex chromosomal rearrangements, including translocations, play a critical role in oncogenesis and are often identified as recurrent genetic aberrations in hematologic malignancies and solid tumors[1]. Translocation detection by karyotyping of an individual’s metaphase chromosomes remains challenging for detection of balanced rearrangements, which have no gain or loss of genetic material. Rearrangement detection using whole-genome sequencing is a viable alternative, but it relies on breakpoint-spanning reads and requires high sequencing depth (30-60x) with long reads to achieve high sensitivity, especially in repetitive regions of the genome[2]. Here, we describe a workflow for translocation detection at low sequencing depth with Oxford Nanopore’s long read chromatin conformation capture technique, Pore-C[3]. Importantly, for translocation calling, Pore-C does not rely on breakpoint-spanning reads but rather the high intrachromosomal interaction frequency of genomic regions around the breakpoint, thus largely lowering sequencing depth requirements and reducing mapping issues in low complexity regions. We prepared, barcoded, and sequenced Pore-C libraries of 3 cancer cell lines including lung cancer (A549), Acute Monocytic Leukemia (THP-1), melanoma (COLO 829) and its matching normal (COLO 829 BL) on a single Q20+ MinION flow cell (FLO-MIN114). Genome-wide contact maps of low-pass data (<1.5 Gbps per sample) were compared against higher depth (>30 Gbps per sample) Hi-C maps [4], revealing that low-pass Pore-C successfully captured large-scale genomic rearrangement. For example, chromosomal pairs chr8-chr11 and chr15-chr19 in A549, chr9-chr11 and chr1-chr20 in THP1, and chr7-chr15 in Colo829 were detected at <0.5X depth of sequencing coverage. These translocation events were further validated by breakpoint analysis using adaptive sampling, where targeted regions were previously confirmed by independent studies[4],[5]. Low-pass Pore-C detects translocations in an unbiased manner and does not require prior knowledge of the translocation structure. Combined with its simple sample-preparation workflow and the capability to provide genome-wide copy-number information in a single experiment, it can serve as a cost-effective and comprehensive tool for cancer genomic studies. [1] Chromosomal translocations in human cancer, Nature, 372, 143 (1994) [2] Hi-C as a tool for precise detection and characterisation of chromosomal rearrangements and copy number variation in human tumors, Genome Biology, 18, 125 (2017) [3] Identifying synergistic high-order 3D chromatin conformations from genome-scale nanopore concatemer sequencing, Nature Biotechnology, 40, 1488(2022) [4] Chromosomal translocations detection in cancer cells using chromosomal conformation capture data, Genes (Basel), 13, 7 (2022) [5] A multi-platform reference for somatic structural variation detection, Cell Genomics, 2, 6 (2022) Citation Format: Scott Hickey, Xiaoguang Dai, Sergey Aganezov, John Beaulaurier, Eoghan Harrington, Sissel Juul. Translocation detection in cancer using low-pass pore-c sequencing [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 405.

Genomic consequences associated with Agrobacterium-mediated transformation of plants.

Attenuated strains of the naturally occurring plant pathogen Agrobacterium tumefaciens can transfer virtually any DNA sequence of interest to model plants and crops. This has made Agrobacterium-mediated transformation (AMT) one of the most commonly used tools in agricultural biotechnology. Understanding AMT, and its functional consequences, is of fundamental importance given that it sits at the intersection of many fundamental fields of study, including plant-microbe interactions, DNA repair/genome stability, and epigenetic regulation of gene expression. Despite extensive research and use of AMT over the last 40 years, the extent of genomic disruption associated with integrating exogenous DNA into plant genomes using this method remains underappreciated. However, new technologies like long-read sequencing make this disruption more apparent, complementing previous findings from multiple research groups that have tackled this question in the past. In this review, we cover progress on the molecular mechanisms involved in Agrobacterium-mediated DNA integration into plant genomes. We also discuss localized mutations at the site of insertion and describe the structure of these DNA insertions, which can range from single copy insertions to large concatemers, consisting of complex DNA originating from different sources. Finally, we discuss the prevalence of large-scale genomic rearrangements associated with the integration of DNA during AMT with examples. Understanding the intended and unintended effects of AMT on genome stability is critical to all plant researchers who use this methodology to generate new genetic variants.

Whole-genome and genome-wide association studies improve key agricultural traits of safflower for industrial and medicinal use.

Safflower (Carthamus tinctorius) is widely cultivated around the world for its seeds and flowers. The presence of linoleic acid (LA) in its seeds and hydroxysafflor yellow A (HSYA) in its flowers are the crucial traits that enable safflower to be used for industrial and medicinal purposes. Understanding the genetic control of these traits is essential for optimizing the quality of safflower and its breeding. To further this research, we present a chromosome-scale assembly of the genome of the safflower variety 'Chuanhonghua 1', which was achieved using an integrated strategy combining Illumina, Oxford Nanopore, and Hi-C sequencing. We obtained a 1.17-Gb assembly with a contig N50 of 1.08Mb, and all assembled sequences were assigned to 12 pseudochromosomes. Safflower's evolution involved the core eudicot γ-triplication event and a whole-genome duplication event, which led to large-scale genomic rearrangements. Extensive genomic shuffling has occurred since the divergence of the ancestor of dicotyledons. We conducted metabolite and transcriptome profiles with time- and part-dependent changes and screened candidate genes that significantly contribute to seed lipid biosynthesis. We also analyzed key gene families that participate in LA and HSYA biosynthesis. Additionally, we re-sequenced 220 safflower lines and carried out a genome-wide association study using high-quality SNP data for eight agronomic traits. We identified SNPs related to important traits in safflower. Besides, the candidate gene HH_034464 (CtCGT1) was shown to be involved in the biosynthesis of HSYA. Overall, we provide a high-quality reference genome and elucidate the genetic basis of LA and HSYA biosynthesis in safflower. This vast amount of data will benefit further research for functional gene mining and breeding in safflower.

Bedaquiline and clofazimine resistance in Mycobacterium tuberculosis: an in-vitro and in-silico data analysis

Bedaquiline is a core drug for the treatment of multidrug-resistant tuberculosis; however, the understanding of resistance mechanisms is poor, which is hampering rapid molecular diagnostics. Some bedaquiline-resistant mutants are also cross-resistant to clofazimine. To decipher bedaquiline and clofazimine resistance determinants, we combined experimental evolution, protein modelling, genome sequencing, and phenotypic data. For this in-vitro and in-silico data analysis, we used a novel in-vitro evolutionary model using subinhibitory drug concentrations to select bedaquiline-resistant and clofazimine-resistant mutants. We determined bedaquiline and clofazimine minimum inhibitory concentrations and did Illumina and PacBio sequencing to characterise selected mutants and establish a mutation catalogue. This catalogue also includes phenotypic and genotypic data of a global collection of more than 14 000 clinical Mycobacterium tuberculosis complex isolates, and publicly available data. We investigated variants implicated in bedaquiline resistance by protein modelling and dynamic simulations. We discerned 265 genomic variants implicated in bedaquiline resistance, with 250 (94%) variants affecting the transcriptional repressor (Rv0678) of the MmpS5-MmpL5 efflux system. We identified 40 new variants in vitro, and a new bedaquiline resistance mechanism caused by a large-scale genomic rearrangement. Additionally, we identified in vitro 15 (7%) of 208 mutations found in clinical bedaquiline-resistant isolates. From our in-vitro work, we detected 14 (16%) of 88 mutations so far identified as being associated with clofazimine resistance and also seen in clinically resistant strains, and catalogued 35 new mutations. Structural modelling of Rv0678 showed four major mechanisms of bedaquiline resistance: impaired DNA binding, reduction in protein stability, disruption of protein dimerisation, and alteration in affinity for its fatty acid ligand. Our findings advance the understanding of drug resistance mechanisms in M tuberculosis complex strains. We have established an extended mutation catalogue, comprising variants implicated in resistance and susceptibility to bedaquiline and clofazimine. Our data emphasise that genotypic testing can delineate clinical isolates with borderline phenotypes, which is essential for the design of effective treatments. Leibniz ScienceCampus Evolutionary Medicine of the Lung, Deutsche Forschungsgemeinschaft, Research Training Group 2501 TransEvo, Rhodes Trust, Stanford University Medical Scientist Training Program, National Institute for Health and Care Research Oxford Biomedical Research Centre, Oxford University Hospitals NHS Foundation Trust, Bill & Melinda Gates Foundation, Wellcome Trust, and Marie Skłodowska-Curie Actions.

Genomic analysis of Leishmania turanica strains from different regions of Central Asia.

The evolution in Leishmania is governed by the opposite forces of clonality and sexual reproduction, with vicariance being an important factor. As such, Leishmania spp. populations may be monospecific or mixed. Leishmania turanica in Central Asia is a good model to compare these two types. In most areas, populations of L. turanica are mixed with L. gerbilli and L. major. Notably, co-infection with L. turanica in great gerbils helps L. major to withstand a break in the transmission cycle. Conversely, the populations of L. turanica in Mongolia are monospecific and geographically isolated. In this work, we compare genomes of several well-characterized strains of L. turanica originated from monospecific and mixed populations in Central Asia in order to shed light on genetic factors, which may drive evolution of these parasites in different settings. Our results illustrate that evolutionary differences between mixed and monospecific populations of L. turanica are not dramatic. On the level of large-scale genomic rearrangements, we confirmed that different genomic loci and different types of rearrangements may differentiate strains originated from mixed and monospecific populations, with genome translocations being the most prominent example. Our data suggests that L. turanica has a significantly higher level of chromosomal copy number variation between the strains compared to its sister species L. major with only one supernumerary chromosome. This suggests that L. turanica (in contrast to L. major) is in the active phase of evolutionary adaptation.

Integration of Hi-C with short and long-read genome sequencing reveals the structure of germline rearranged genomes

Structural variants are a common cause of disease and contribute to a large extent to inter-individual variability, but their detection and interpretation remain a challenge. Here, we investigate 11 individuals with complex genomic rearrangements including germline chromothripsis by combining short- and long-read genome sequencing (GS) with Hi-C. Large-scale genomic rearrangements are identified in Hi-C interaction maps, allowing for an independent assessment of breakpoint calls derived from the GS methods, resulting in >300 genomic junctions. Based on a comprehensive breakpoint detection and Hi-C, we achieve a reconstruction of whole rearranged chromosomes. Integrating information on the three-dimensional organization of chromatin, we observe that breakpoints occur more frequently than expected in lamina-associated domains (LADs) and that a majority reshuffle topologically associating domains (TADs). By applying phased RNA-seq, we observe an enrichment of genes showing allelic imbalanced expression (AIG) within 100 kb around the breakpoints. Interestingly, the AIGs hit by a breakpoint (19/22) display both up- and downregulation, thereby suggesting different mechanisms at play, such as gene disruption and rearrangements of regulatory information. However, the majority of interpretable genes located 200 kb around a breakpoint do not show significant expression changes. Thus, there is an overall robustness in the genome towards large-scale chromosome rearrangements.

APOBEC mutagenesis, kataegis, chromothripsis in EGFR-mutant osimertinib-resistant lung adenocarcinomas

Studies of targeted therapy resistance in lung cancer have primarily focused on single-gene alterations. Based on prior work implicating apolipoprotein b mRNA-editing enzyme, catalytic polypeptide-like (APOBEC) mutagenesis in histological transformation of epidermal growth factor receptor (EGFR)-mutant lung cancers, we hypothesized that mutational signature analysis may help elucidate acquired resistance to targeted therapies. APOBEC mutational signatures derived from an Food and Drug Administration-cleared multigene panel [Memorial Sloan Kettering Cancer Center Integrated Mutation Profiling of Actionable Cancer Targets (MSK-IMPACT)] using the Signature Multivariate Analysis (SigMA) algorithm were validated against the gold standard of mutational signatures derived from whole-exome sequencing. Mutational signatures were decomposed in 3276 unique lung adenocarcinomas (LUADs), including 93 paired osimertinib-naïve and -resistant EGFR-mutant tumors. Associations between APOBEC and mechanisms of resistance to osimertinib were investigated. Whole-genome sequencing was carried out on available EGFR-mutant lung cancer samples (10 paired, 17 unpaired) to investigate large-scale genomic alterations potentially contributing to osimertinib resistance. APOBEC mutational signatures were more frequent in receptor tyrosine kinase (RTK)-driven lung cancers (EGFR, ALK, RET, and ROS1; 25%) compared to LUADs at large (20%, P < 0.001); across all subtypes, APOBEC mutational signatures were enriched in subclonal mutations (P < 0.001). In EGFR-mutant lung cancers, osimertinib-resistant samples more frequently displayed an APOBEC-dominant mutational signature compared to osimertinib-naïve samples (28% versus 14%, P= 0.03). Specifically, mutations detected in osimertinib-resistant tumors but not in pre-treatment samples significantly more frequently displayed an APOBEC-dominant mutational signature (44% versus 23%, P < 0.001). EGFR-mutant samples with APOBEC-dominant signatures had enrichment of large-scale genomic rearrangements (P= 0.01) and kataegis (P= 0.03) in areas of APOBEC mutagenesis. APOBEC mutational signatures are frequent in RTK-driven LUADs and increase under the selective pressure of osimertinib in EGFR-mutant lung cancer. APOBEC mutational signature enrichment in subclonal mutations, private mutations acquired after osimertinib treatment, and areas of large-scale genomic rearrangements highlights a potentially fundamental role for APOBEC mutagenesis in the development of resistance to targeted therapies, which may be potentially exploited to overcome such resistance.

Phylogenetic analysis based on single-copy orthologous proteins in highly variable chloroplast genomes of Corydalis

Corydalis is one of the few lineages that have been reported to have extensive large-scale chloroplast genome (cp-genome) rearrangements. In this study, novel cp-genome rearrangements of Corydalis pinnata, C. mucronate, and C. sheareri are described. C. pinnata is a narrow endemic species only distributed at Qingcheng Mountain in southwest China. Two independent relocations of the same four genes (trnM-CAU-rbcL) were found relocated from the typically posterior part of the large single-copy region to the front of it. A uniform inversion of an 11–14-kb segment (ndhB-trnR-ACG) was found in the inverted repeat region; and extensive losses of accD, clpP, and trnV-UAC genes were detected in all cp-genomes of all three species of Corydalis. In addition, a phylogenetic tree was reconstructed based on 31 single-copy orthologous proteins in 27 cp-genomes. This study provides insights into the evolution of cp-genomes throughout the genus Corydalis and also provides a reference for further studies on the taxonomy, identification, phylogeny, and genetic transformation of other lineages with extensive rearrangements in cp-genomes.

Multiparametric High-Content Cell Painting Identifies Copper Ionophores as Selective Modulators of Esophageal Cancer Phenotypes.

Esophageal adenocarcinoma is of increasing global concern due to increasing incidence, a lack of effective treatments, and poor prognosis. Therapeutic target discovery and clinical trials have been hindered by the heterogeneity of the disease, the lack of “druggable” driver mutations, and the dominance of large-scale genomic rearrangements. We have previously undertaken a comprehensive small-molecule phenotypic screen using the high-content Cell Painting assay to quantify the morphological response to a total of 19,555 small molecules across a panel of genetically distinct human esophageal cell lines to identify new therapeutic targets and small molecules for the treatment of esophageal adenocarcinoma. In this current study, we report for the first time the dose–response validation studies for the 72 screening hits from the target-annotated LOPAC and Prestwick FDA-approved compound libraries and the full list of 51 validated esophageal adenocarcinoma-selective small molecules (71% validation rate). We then focus on the most potent and selective hit molecules, elesclomol, disulfiram, and ammonium pyrrolidinedithiocarbamate. Using a multipronged, multitechnology approach, we uncover a unified mechanism of action and a vulnerability in esophageal adenocarcinoma toward copper-dependent cell death that could be targeted in the future.

Identification of novel, functional, long noncoding RNAs involved in programmed, large-scale genome rearrangements.

Noncoding RNAs (ncRNAs) make up to ∼98% percent of the transcriptome of a given organism. In recent years, one relatively new class of ncRNAs, long noncoding RNAs (lncRNAs), were shown to be more than mere by-products of gene expression and regulation. The unicellular eukaryote Paramecium tetraurelia is a member of the ciliate phylum, an extremely heterogeneous group of organisms found in most bodies of water across the globe. A hallmark of ciliate genetics is nuclear dimorphism and programmed elimination of transposons and transposon-derived DNA elements, the latter of which is essential for the maintenance of the somatic genome. Paramecium and ciliates in general harbor a plethora of different ncRNA species, some of which drive the process of large-scale genome rearrangements, including DNA elimination, during sexual development. Here, we identify and validate the first known functional lncRNAs in ciliates to date. Using deep-sequencing and subsequent bioinformatic processing and experimental validation, we show that Paramecium expresses at least 15 lncRNAs. These candidates were predicted by a highly conservative pipeline, and informatic analyses hint at differential expression during development. Depletion of two lncRNAs, lnc1 and lnc15, resulted in clear phenotypes, decreased survival, morphological impairment, and a global effect on DNA elimination.

P-OGC32 Copper-dependent cell death is a cancer specific vulnerability in oesophageal adenocarcinoma and provides the opportunity for future biomarker-stratified clinical trials

Abstract Background Oesophageal adenocarcinoma (OAC) is of increasing global concern due to increasing incidence, lack of effective treatments, and poor prognosis. Therapeutic target discovery and clinical trials have been hindered by the heterogeneity of the disease, dominance of large-scale genomic rearrangements, and lack of driver mutations. We have profiled small-molecule compounds using an innovative high content imaging assay in a panel of transformed and non-transformed oesophageal cell lines to identify OAC-specific cytotoxic compounds, new therapeutic targets, potential drug repurposing opportunities, and chemical starting points for the treatment of OAC. Methods We have comprehensively profiled 19,555 small-molecule compounds using an innovative high content assay to quantify 1000’s of subcellular imaging features to capture important phenotypes missed by standard approaches. Prioritised molecules then underwent functional, transcriptomic, and metabolomic characterisation across panels of oesophageal cell lines and patient-derived organoids for the identification of OAC-specific drug mechanisms. Results We identified 72 lead compounds as exhibiting OAC-specific cytotoxicity and characterised three of the most potent and selective compounds in depth, each of different proposed classes and chemical structures. Using several orthogonal methods we uncovered a unified mechanism of action and a targetable vulnerability in OAC involving copper-dependent cell death. Strikingly no phenotypic effects or changes in gene-expression were observed following treatment with these compounds in non-transformed oesophageal cell lines or normal gastric organoids providing support for this mechanism as a cancer-specific phenomenon. Conclusions We have applied high content imaging, transcriptomic and metabolomic analyses to reveal a unique vulnerability in OAC. We have defined a unified mechanism of OAC-specific copper-dependent cell death for the three highly potent compounds. Finally, through the integration of transcriptomic and metabolomics analyses we gained insight into drug sensitivity and provide the basis for a future biomarker-stratified clinical trial of these drugs in OAC.

The potential of family-free rearrangements towards gene orthology inference.

Recently, we proposed an efficient ILP formulation [Rubert DP, Martinez FV, Braga MDV, Natural family-free genomic distance, Algorithms Mol Biol16:4, 2021] for exactly computing the rearrangement distance of two genomes in a family-free setting. In such a setting, neither prior classification of genes into families, nor further restrictions on the genomes are imposed. Given two genomes, the mentioned ILP computes an optimal matching of the genes taking into account simultaneously local mutations, given by gene similarities, and large-scale genome rearrangements. Here, we explore the potential of using this ILP for inferring groups of orthologs across several species. More precisely, given a set of genomes, our method first computes all pairwise optimal gene matchings, which are then integrated into gene families in the second step. Our approach is implemented into a pipeline incorporating the pre-computation of gene similarities. It can be downloaded from gitlab.ub.uni-bielefeld.de/gi/FFGC. We obtained promising results with experiments on both simulated and real data.

Large scale genomic rearrangements in selected Arabidopsis thaliana T-DNA lines are caused by T-DNA insertion mutagenesis

BackgroundExperimental proof of gene function assignments in plants is based on mutant analyses. T-DNA insertion lines provided an invaluable resource of mutants and enabled systematic reverse genetics-based investigation of the functions of Arabidopsis thaliana genes during the last decades.ResultsWe sequenced the genomes of 14 A. thaliana GABI-Kat T-DNA insertion lines, which eluded flanking sequence tag-based attempts to characterize their insertion loci, with Oxford Nanopore Technologies (ONT) long reads. Complex T-DNA insertions were resolved and 11 previously unknown T-DNA loci identified, resulting in about 2 T-DNA insertions per line and suggesting that this number was previously underestimated. T-DNA mutagenesis caused fusions of chromosomes along with compensating translocations to keep the gene set complete throughout meiosis. Also, an inverted duplication of 800 kbp was detected. About 10 % of GABI-Kat lines might be affected by chromosomal rearrangements, some of which do not involve T-DNA. Local assembly of selected reads was shown to be a computationally effective method to resolve the structure of T-DNA insertion loci. We developed an automated workflow to support investigation of long read data from T-DNA insertion lines. All steps from DNA extraction to assembly of T-DNA loci can be completed within days.ConclusionsLong read sequencing was demonstrated to be an effective way to resolve complex T-DNA insertions and chromosome fusions. Many T-DNA insertions comprise not just a single T-DNA, but complex arrays of multiple T-DNAs. It is becoming obvious that T-DNA insertion alleles must be characterized by exact identification of both T-DNA::genome junctions to generate clear genotype-to-phenotype relations.

Plasmodium vinckei genomes provide insights into the pan-genome and evolution of rodent malaria parasites

BackgroundRodent malaria parasites (RMPs) serve as tractable tools to study malaria parasite biology and host-parasite-vector interactions. Among the four RMPs originally collected from wild thicket rats in sub-Saharan Central Africa and adapted to laboratory mice, Plasmodium vinckei is the most geographically widespread with isolates collected from five separate locations. However, there is a lack of extensive phenotype and genotype data associated with this species, thus hindering its use in experimental studies.ResultsWe have generated a comprehensive genetic resource for P. vinckei comprising of five reference-quality genomes, one for each of its subspecies, blood-stage RNA sequencing data for five P. vinckei isolates, and genotypes and growth phenotypes for ten isolates. Additionally, we sequenced seven isolates of the RMP species Plasmodium chabaudi and Plasmodium yoelii, thus extending genotypic information for four additional subspecies enabling a re-evaluation of the genotypic diversity and evolutionary history of RMPs.The five subspecies of P. vinckei have diverged widely from their common ancestor and have undergone large-scale genome rearrangements. Comparing P. vinckei genotypes reveals region-specific selection pressures particularly on genes involved in mosquito transmission. Using phylogenetic analyses, we show that RMP multigene families have evolved differently across the vinckei and berghei groups of RMPs and that family-specific expansions in P. chabaudi and P. vinckei occurred in the common vinckei group ancestor prior to speciation. The erythrocyte membrane antigen 1 and fam-c families in particular show considerable expansions among the lowland forest-dwelling P. vinckei parasites. The subspecies from the highland forests of Katanga, P. v. vinckei, has a uniquely smaller genome, a reduced multigene family repertoire and is also amenable to transfection making it an ideal parasite for reverse genetics. We also show that P. vinckei parasites are amenable to genetic crosses.ConclusionsPlasmodium vinckei isolates display a large degree of phenotypic and genotypic diversity and could serve as a resource to study parasite virulence and immunogenicity. Inclusion of P. vinckei genomes provide new insights into the evolution of RMPs and their multigene families. Amenability to genetic crossing and transfection make them also suitable for classical and functional genetics to study Plasmodium biology.

Comparative genome analysis of multidrug-resistant and susceptible Glaesserella parasuis strains isolated from diseased pigs in China

Glaesserella parasuis (G. parasuis) is a respiratory pathogen of swine and the etiological agent of Glässer’s disease. Although the emergence of multidrug-resistant (MDR) G. parasuis is a critical problem in the swine industry, there are few publications on the genetic basis of antimicrobial resistance of G. parasuis. In this study, comparative genome analyses were used to identify genomic differences between two phenotypically distinct isolates, an MDR isolate (HPS-1) and a susceptible isolate (HPS-2), from diseased swines in China. These isolates were both serovar 4, which is predominant in cases of Glässer’s disease and is the most prevalent serovar in China. Based on clusters of orthologous group (COG) annotations, genes assigned to the extracellular structure category were only detected in HPS-1 and genes related to cell motility were more abundant in HPS-1 than in HPS-2. A comparative genomic analysis showed that these two isolates are closely related, although there was a large-scale genomic rearrangement. Eighteen percent of the genome consisted of strain-specific accessory genes of HPS-1. Notably, only the two genes aac(6′)-Ie-aph(2′')-Ia and blaROB-1 on a plasmid were specific to HPS-1. We also detected 30,599 single nucleotide polymorphisms (SNPs), including nonsynonymous SNPs in the aminoglycoside resistance gene aph(3′')-Ib, the fusidic acid resistance gene fusA, and the two rifampicin resistance genes rpoC and rpoB in HPS-1. These findings improve our understanding of the differences between MDR and susceptible isolates and will aid the development of treatment strategies to decrease the prevalence and disease burden caused by G. parasuis.

A mapping platform for mitotic crossover by single-cell multi-omics.

Mitotic crossovers have the potential to cause large-scale genome rearrangements. Here, we describe high-throughput, single-cell, whole-genome sequencing methods for mapping crossovers genome-wide at scale. The methods are generalizable to various eukaryotes and to other end points requiring high-throughput, high-coverage single cell sequencing.

Unraveling Genome Evolution Throughout Visual Analysis: The XCout Portal.

Due to major breakthroughs in sequencing technologies throughout the last decades, the time and cost per sequencing experiment have reduced drastically, overcoming the data generation barrier during the early genomic era. Such a shift has encouraged the scientific community to develop new computational methods that are able to compare large genomic sequences, thus enabling large-scale studies of genome evolution. The field of comparative genomics has proven itself invaluable for studying the evolutionary mechanisms and the forces driving genome evolution. In this line, a full genome comparison study between 2 species requires a quadratic number of comparisons in terms of the number of sequences (around 400 chromosome comparisons in the case of mammalian genomes); however, when studying conserved syntenies or evolutionary rearrangements, many sequence comparisons can be skipped for not all will contain significant signals. Subsequently, the scientific community has developed fast heuristics to perform multiple pairwise comparisons between large sequences to determine whether significant sets of conserved similarities exist. The data generation problem is no longer an issue, yet the limitations have shifted toward the analysis of such massive data. Therefore, we present XCout, a Web-based visual analytics application for multiple genome comparisons designed to improve the analysis of large-scale evolutionary studies using novel techniques in Web visualization. XCout enables to work on hundreds of comparisons at once, thus reducing the time of the analysis by identifying significant signals between chromosomes across multiple species. Among others, XCout introduces several techniques to aid in the analysis of large-scale genome rearrangements, particularly (1) an interactive heatmap interface to display comparisons using automatic color scales based on similarity thresholds to ease detection at first sight, (2) an overlay system to detect individual signal contributions between chromosomes, (3) a tracking tool to trace conserved blocks across different species to perform evolutionary studies, and (4) a search engine to search annotations throughout different species.