70 research outputs found

    Combination of short-read, long-read, and optical mapping assemblies reveals large-scale tandem repeat arrays with population genetic implications

    Get PDF
    Accurate and contiguous genome assembly is key to a comprehensive understanding of the processes shaping genomic diversity and evolution. Yet, it is frequently constrained by constitutive heterochromatin, usually characterized by highly repetitive DNA. As a key feature of genome architecture associated with centromeric and subtelomeric regions, it locally influences meiotic recombination. In this study, we assess the impact of large tandem repeat arrays on the recombination rate landscape in an avian speciation model, the Eurasian crow. We assembled two high-quality genome references using single-molecule real-time sequencing (long-read assembly [LR]) and single-molecule optical maps (optical map assembly [OM]). A three-way comparison including the published short-read assembly (SR) constructed for the same individual allowed assessing assembly properties and pinpointing misassemblies. By combining information from all three assemblies, we characterized 36 previously unidentified large repetitive regions in the proximity of sequence assembly breakpoints, the majority of which contained complex arrays of a 14-kb satellite repeat or its 1.2-kb subunit. Using whole-genome population resequencing data, we estimated the population-scaled recombination rate (ρ) and found it to be significantly reduced in these regions. These findings are consistent with an effect of low recombination in regions adjacent to centromeric or subtelomeric heterochromatin and add to our understanding of the processes generating widespread heterogeneity in genetic diversity and differentiation along the genome. By combining three different technologies, our results highlight the importance of adding a layer of information on genome structure that is inaccessible to each approach independently

    Genomic analyses of the Linum distyly supergene reveal convergent evolution at the molecular level

    Get PDF
    Supergenes govern multi-trait-balanced polymorphisms in a wide range of systems; however, our understanding of their origins and evolution remains incomplete. The reciprocal placement of stigmas and anthers in pin and thrum floral morphs of distylous species constitutes an iconic example of a balanced polymorphism governed by a supergene, the distyly S-locus. Recent studies have shown that the Primula and Turnera distyly supergenes are both hemizygous in thrums, but it remains unknown whether hemizygosity is pervasive among distyly S-loci. As hemizygosity has major consequences for supergene evolution and loss, clarifying whether this genetic architecture is shared among distylous species is critical. Here, we have characterized the genetic architecture and evolution of the distyly supergene in Linum by generating a chromosome-level genome assembly of Linum tenue, followed by the identification of the S-locus using population genomic data. We show that hemizygosity and thrum-specific expression of S-linked genes, including a pistil-expressed candidate gene for style length, are major features of the Linum S-locus. Structural variation is likely instrumental for recombination suppression, and although the non-recombining dominant haplotype has accumulated transposable elements, S-linked genes are not under relaxed purifying selection. Our findings reveal remarkable convergence in the genetic architecture and evolution of independently derived distyly supergenes, provide a counterexample to classic inversion-based supergenes, and shed new light on the origin and maintenance of an iconic floral polymorphism.European Research Council (ERC) 757451Swedish Research Council 2019-04452, 2018-0597

    Identifying the causes and consequences of assembly gaps using a multiplatform genome assembly of a bird-of-paradise

    Get PDF
    Genome assemblies are currently being produced at an impressive rate by consortia and individual laboratories. The low costs and increasing efficiency of sequencing technologies now enable assembling genomes at unprecedented quality and contiguity. However, the difficulty in assembling repeat-rich and GC-rich regions (genomic “dark matter”) limits insights into the evolution of genome structure and regulatory networks. Here, we compare the efficiency of currently available sequencing technologies (short/linked/long reads and proximity ligation maps) and combinations thereof in assembling genomic dark matter. By adopting different de novo assembly strategies, we compare individual draft assemblies to a curated multiplatform reference assembly and identify the genomic features that cause gaps within each assembly. We show that a multiplatform assembly implementing long-read, linked-read and proximity sequencing technologies performs best at recovering transposable elements, multicopy MHC genes, GC-rich microchromosomes and the repeat-rich W chromosome. Telomere-to-telomere assemblies are not a reality yet for most organisms, but by leveraging technology choice it is now possible to minimize genome assembly gaps for downstream analysis. We provide a roadmap to tailor sequencing projects for optimized completeness of both the coding and noncoding parts of nonmodel genomes

    Genetic barriers to historical gene flow between cryptic species of alpine bumblebees revealed by comparative population genomics

    Get PDF
    Evidence is accumulating that gene flow commonly occurs between recently diverged species, despite the existence of barriers to gene flow in their genomes. However, we still know little about what regions of the genome become barriers to gene flow and how such barriers form. Here, we compare genetic differentiation across the genomes of bumblebee species living in sympatry and allopatry to reveal the potential impact of gene flow during species divergence and uncover genetic barrier loci. We first compared the genomes of the alpine bumblebee Bombus sylvicola and a previously unidentified sister species living in sympatry in the Rocky Mountains, revealing prominent islands of elevated genetic divergence in the genome that colocalize with centromeres and regions of low recombination. This same pattern is observed between the genomes of another pair of closely related species living in allopatry (B. bifarius and B. vancouverensis). Strikingly however, the genomic islands exhibit significantly elevated absolute divergence (dXY) in the sympatric, but not the allopatric, comparison indicating that they contain loci that have acted as barriers to historical gene flow in sympatry. Our results suggest that intrinsic barriers to gene flow between species may often accumulate in regions of low recombination and near centromeres through processes such as genetic hitchhiking, and that divergence in these regions is accentuated in the presence of gene flow

    An RND-Type Efflux System in Borrelia burgdorferi Is Involved in Virulence and Resistance to Antimicrobial Compounds

    Get PDF
    Borrelia burgdorferi is remarkable for its ability to thrive in widely different environments due to its ability to infect various organisms. In comparison to enteric Gram-negative bacteria, these spirochetes have only a few transmembrane proteins some of which are thought to play a role in solute and nutrient uptake and excretion of toxic substances. Here, we have identified an outer membrane protein, BesC, which is part of a putative export system comprising the components BesA, BesB and BesC. We show that BesC, a TolC homolog, forms channels in planar lipid bilayers and is involved in antibiotic resistance. A besC knockout was unable to establish infection in mice, signifying the importance of this outer membrane channel in the mammalian host. The biophysical properties of BesC could be explained by a model based on the channel-tunnel structure. We have also generated a structural model of the efflux apparatus showing the putative spatial orientation of BesC with respect to the AcrAB homologs BesAB. We believe that our findings will be helpful in unraveling the pathogenic mechanisms of borreliae as well as in developing novel therapeutic agents aiming to block the function of this secretion apparatus

    Borrelia channel-forming proteins : structure and function

    No full text
    Borrelia is a Gram-negative, corkscrew-shaped bacterium transmitted by infected ticks or lice. Borreliae are subdivided into pathogens of two diseases: Lyme disease, caused mainly by B. burgdorferi, B. afzelii and B. garinii; and relapsing fever caused primarily by B. duttonii, B. hermsii, B. recurrentis or B. crocidurae. Both diseases differ in their manifestations, duration times and dissemination patterns. Antibiotics are the major therapeutics, although unfortunately antibiotic treatment is not always beneficial. To date, drug resistance mechanisms in B. burgdorferi are unknown. Transporters of the resistance-nodulation-division (RND) family appear to be involved in drug resistance, especially in Gram-negative bacteria. They consist of three components: a cytoplasmic membrane export system, a membrane fusion protein (MFP), and an outer membrane factor (OMP). The major antibiotic efflux activity of this type in Escherichia coli is mediated by the tripartite multidrug resistance pump AcrAB-TolC. Based on the sequence homology we conclude that the besA (bb0140), besB (bb0141) and besC (bb0142) genes code for a similar efflux system in B. burgdorferi. We created a deletion mutant of besC. The minimal inhibitory concentration (MIC) values of B. burgdorferi carrying an inactive besC gene were 4- to 8-fold lower than in the wild type strain. Animal experiments showed that the besC mutant was unable to infect mice. Black lipid bilayer experiments were carried out to determine the biophysical properties of purified BesC. This study showed the importance of BesC protein for B. burgdorferi pathogenicity and resistance to antibiotics, although its importance in clinical isolates is not known. Due to its small genome, Borrelia is metabolically and biosynthetically deficient, thereby making it highly dependent on nutrients provided by their hosts. The uptake of nutrients by Borrelia is not yet completely understood. We describe the purification and characterization of a 36-kDa protein that functions as a putative dicarboxylate-specific porin in the outer membrane of Borrelia. The protein was designated as DipA, for dicarboxylate-specific porin A. DipA was biophysically characterized using the black lipid bilayer assay. The permeation of KCl through the channel could be partly blocked by titrating the DipA-mediated membrane conductance with increasing concentrations of different organic dicarboxylic anions. The obtained results imply that DipA does not form a general diffusion pore, but a porin with a binding site specific for dicarboxylates which play important key roles in the deficient metabolic and biosynthetic pathways of Borrelia species. The presence of porin P66 has been shown in both Lyme disease and relapsing fever spirochetes. In our study, purified P66 homologues from Lyme disease species B. burgdorferi, B. afzelii and B. garinii and relapsing fever species B. duttonii, B. recurrentis and B. hermsii were compared and their biophysical properties were further characterized in black lipid bilayer assay. Subsequently, the channel diameter of B. burgdorferi P66 was investigated in more detail. For this study, different nonelectrolytes with known hydrodynamic radii were used. This allowed us to determine the effective diameter of the P66 channel lumen. Furthermore, the blockage of the channel after addition of nonelectrolytes revealed seven subconducting states and indicated a heptameric structure of the P66 channel. These results may give more insight into the functional properties of this important porin

    The Etiological Agent of Lyme Disease, Borrelia burgdorferi, Appears To Contain Only a Few Small RNA Molecules

    No full text
    Small regulatory RNAs (sRNAs) have recently been shown to be the main controllers of several regulatory pathways. The function of sRNAs depends in many cases on the RNA-binding protein Hfq, especially for sRNAs with an antisense function. In this study, the genome of Borrelia burgdorferi was subjected to different searches for sRNAs, including direct homology and comparative genomics searches and ortholog- and annotation-based search strategies. Two new sRNAs were found, one of which showed complementarity to the rpoS region, which it possibly controls by an antisense mechanism. The role of the other sRNA is unknown, although observed complementarities against particular mRNA sequences suggest an antisense mechanism. We suggest that the low level of sRNAs observed in B. burgdorferi is at least partly due to the presumed lack of both functional Hfq protein and RNase E activity
    corecore