Search

Over the last six decades, populations of the bumblebees Bombus sylvicola and Bombus balteatus in Colorado have experienced decreases in tongue length, a trait important for plant-pollinator mutualisms. It has been hypothesized that this observation reflects selection resulting from shifts in floral composition under climate change. Here we used morphometrics and population genomics to determine whether morphological change is ongoing, investigate the genetic basis of morphological variation, and analyse population structure in these populations. We generated a genome assembly of B. balteatus. We then analysed whole-genome sequencing data and morphometric measurements of 580 samples of both species from seven high-altitude localities. Out of 281 samples originally identified as B. sylvicola, 67 formed a separate genetic cluster comprising a newly-discovered cryptic species ("incognitus"). However, an absence of genetic structure within species suggests that gene flow is common between mountains. We found a significant decrease in tongue length between bees collected between 2012-2014 and in 2017, indicating that morphological shifts are ongoing. We did not discover any genetic associations with tongue length, but a SNP related to production of a proteolytic digestive enzyme was implicated in body size variation. We identified evidence of covariance between kinship and both tongue length and body size, which is suggestive of a genetic component of these traits, although it is possible that shared environmental effects between colonies are responsible. Our results provide evidence for ongoing modification of a morphological trait important for pollination and indicate that this trait probably has a complex genetic and environmental basis.

Over the last six decades, populations of the bumblebees Bombus sylvicola and Bombus balteatus in Colorado have experienced decreases in tongue length, a trait important for plant-pollinator mutualisms. It has been hypothesized that this observation reflects selection resulting from shifts in floral composition under climate change. Here we used morphometrics and population genomics to determine whether morphological change is ongoing, investigate the genetic basis of morphological variation, and analyze population structure in these populations. We generated a genome assembly of B. balteatus. We then analyzed whole-genome sequencing data and morphometric measurements of 580 samples of both species from seven high-altitude localities. Out of 281 samples originally identified as B. sylvicola, 67 formed a separate genetic cluster comprising a newly-discovered cryptic species (“incognitus”). However, an absence of genetic structure within species suggests that gene flow is common between mountains. We found a significant decrease in tongue length between bees collected between 2012-2014 and in 2017, indicating that morphological shifts are ongoing. We did not discover any genetic associations with tongue length, but a SNP related to production of a proteolytic digestive enzyme was implicated in body size variation. We identified evidence of covariance between kinship and both tongue length and body size, which is suggestive of a genetic component of these traits, although it is possible that shared environmental effects between colonies are responsible. Our results provide evidence for ongoing modification of a morphological trait important for pollination and indicate that this trait likely has a complex genetic and environmental basis.

Background The vaginal microbiota has been reported to be associated with HPV infection and cervical cancer. This study was performed to compare the vaginal microbiota at two timepoints in women performing self-sampling and had a persistent or transient HPV16 infection. The women were tested for 12 high-risk HPV (hrHPV) types but only women with single type (HPV16) were included to reduce confounding variables. Methods In total 96 women were included in this study. Of these, 26 were single positive for HPV16 in the baseline test and HPV negative in the follow-up test and 38 were single positive for HPV16 in both tests and diagnosed with CIN2+ in histology. In addition, 32 women that were negative for all 12 HPV tested were included. The samples of vaginal fluid were analyzed with the Ion 16S™ Metagenomics Kit and Ion 16S™ metagenomics module within the Ion Reporter™ software. Results K-means clustering resulted in two Lactobacillus-dominated groups, one with Lactobacillus sp. and the other specifically with Lactobacillus iners. The two remaining clusters were dominated by a mixed non-Lactobacillus microbiota. HPV negative women had lower prevalence (28%) of the non-Lactobacill dominant cluster in the baseline test, as compared to women with HPV16 infection (42%) (p value = 0.0173). Transition between clusters were more frequent in women with persistent HPV16 infection (34%) as compared in women who cleared the HPV16 infection (19%) (p value = 0.036). Conclusions The vaginal microbiota showed a higher rate of transitioning between bacterial profiles in women with persistent HPV16 infection as compared to women with transient infection. This indicate an instability in the microenvironment in women with persistent HPV infection and development of CIN2+.

How the avian sex chromosomes first evolved from autosomes remains elusive as 100 million years (Myr) of divergence and degeneration obscure their evolutionary history. Sylvioidea songbirds is an emerging model for understanding avian sex chromosome evolution because a unique chromosome fusion event ∼24 Myr ago has formed enlarged “neo-sex chromosomes” consisting of an added (new) and an ancestral (old) part. Here, we report the female genome (ZW) of one Sylvioidea species, the great reed warbler (Acrocephalus arundinaceus). We confirm that the added region has been translocated to both Z and W, and show that the added-W part has been heavily reorganised within itself and with the ancestral-W. Next, we show that recombination between Z and W continued after the fusion event, and that recombination suppression took ∼10 Myr to be completed and arose locally and non-linearly along the sex chromosomes. This pattern is inconsistent with that of large inversions and instead suggests gradual and mosaic recombination suppression. We find that the added-W mirrors the ancestral-W in terms of repeat accumulation, loss of genetic variation, and gene degeneration. Lastly, we show that genes being maintained on W are slowly evolving and dosage sensitive, and that highly conserved genes across broad taxonomic groups, regardless of sex-linkage, evolve slower on both Z and W. This study reveals complex expansion of recombination suppression along avian sex chromosomes, and that the evolutionary trajectory of sex-linked genes is highly predictable and governed partly by sex-linkage per se, partly by their functional properties.

Structural variation (SV) accounts for a substantial part of genetic mutations segregating across eukaryotic genomes with important medical and evolutionary implications. Here, we characterized SV across evolutionary time scales in the songbird genus Corvus using de novo assembly and read mapping approaches. Combining information from short-read (N = 127) and long-read re-sequencing data (N = 31) as well as from optical maps (N = 16) revealed a total of 201,738 insertions, deletions and inversions. Population genetic analysis of SV in the Eurasian crow speciation model revealed an evolutionary young (~530,000 years) cis-acting 2.25-kb retrotransposon insertion reducing expression of the NDP gene with consequences for premating isolation. Our results attest to the wealth of SV segregating in natural populations and demonstrate its evolutionary significance.

To further characterize the function of the Borrelia burgdorferi C-terminal protease CtpA, we used site-directed mutagenesis to alter the putative CtpA cleavage site of one of its known substrates, the outer membrane (OM) porin P13. These mutations resulted in only partial blockage of P13 processing. Ectopic expression of a C-terminally truncated P13 in B. burgdorferi indicated that the C-terminal peptide functions as a safeguard against misfolding or mislocalization prior to its proteolytic removal by CtpA. In a parallel study of Borrelia burgdorferi lipoprotein sorting mechanisms, we observed a lower-molecular-weight variant of surface lipoprotein OspC that was particularly prominent with OspC mutants that mislocalized to the periplasm or contained C-terminal epitope tags. Further investigation revealed that the variant resulted from C-terminal proteolysis by CtpA. Together, these findings indicate that CtpA rather promiscuously targets polypeptides that lack structurally constrained C termini, as proteolysis appears to occur independently of a specific peptide recognition sequence. Low-level processing of surface lipoproteins such as OspC suggests the presence of a CtpA-dependent quality control mechanism that may sense proper translocation of integral outer membrane proteins and surface lipoproteins by detecting the release of C-terminal peptides.

Borrelia burgdorferi has an unusual genomic structure containing 21 plasmids. These plasmids carry genes that are essential for infectivity and survival of the spirochetes in vivo. Several plasmids are lost during cultivation in vitro, which might lead to a heterogeneous population after multiple passages and loss of infectivity in laboratory animals. Herein, we present a simple and inexpensive multiplex PCR method that detects the complete plasmid profile of B. burgdorferi B31 in just two PCR tubes.

Relapsing fever is a worldwide, endemic disease caused by several spirochetal species belonging to the genus Borrelia . During the recurring fever peaks, borreliae proliferate remarkably quickly compared to the slow dissemination of Lyme disease Borrelia and therefore require efficient nutrient uptake from the blood of their hosts. This study describes the identification and characterization of the first relapsing fever porin, which is present in the outer membranes of B. duttonii , B. hermsii , B. recurrentis , and B. turicatae . The pore-forming protein was purified by hydroxyapatite chromatography and designated Oms38, for outer membrane-spanning protein of 38 kDa. Biophysical characterization of Oms38 was done by using the black lipid bilayer method, demonstrating that Oms38 forms small, water-filled channels of 80 pS in 1 M KCl that did not exhibit voltage-dependent closure. The Oms38 channel is slightly selective for anions and shows a ratio of permeability for cations over anions of 0.41 in KCl. Analysis of the deduced amino acid sequences demonstrated that Oms38 contains an N-terminal signal sequence which is processed under in vivo conditions. Oms38 is highly conserved within the four studied relapsing fever species, sharing an overall amino acid identity of 58% and with a strong indication for the presence of amphipathic β-sheets.

Relapsing fever caused by Borrelia crocidurae and B. duttonii in Togo may be misdiagnosed., Given the prevalence of relapsing fever (RF) in Senegal, this disease may cause illness and death in other areas of West Africa. We performed a cross-sectional, clinic-based study to investigate the presence of RF in Togo during 2002–2004. Blood samples from patients with fever were examined for RF spirochetes by microscopy, PCR, and DNA sequencing of amplicons and for antibodies to the glycerophosphodiester phosphodiesterase antigen. Although no spirochetes were seen in blood smears, ≈10% of the patients were positive by PCR and ≈13% were seropositive for spirochetes. DNA sequencing demonstrated that Borrelia crocidurae and B. duttonii were present. Most patients were treated for malaria whether or not plasmodia were observed. Thus, many RF patients originally had a misdiagnosis of malaria, which resulted in ineffective treatment. The inability of microscopic analysis to detect spirochetes compared with PCR demonstrates the need for tests with greater sensitivity.

The main Borrelia species causing Lyme borreliosis in Europe and Asia are Borrelia afzelii, B. garinii, B. burgdorferi and B. bavariensis. This is in contrast to the United States, where infections are exclusively caused by B. burgdorferi. Until to date the genome sequences of four B. afzelii strains, of which only two include the numerous plasmids, are available. In order to further assess the genetic diversity of B. afzelii, the most common species in Europe, responsible for the large variety of clinical manifestations of Lyme borreliosis, we have determined the full genome sequence of the B. afzelii strain K78, a clinical isolate from Austria. The K78 genome contains a linear chromosome (905,949 bp) and 13 plasmids (8 linear and 5 circular) together presenting 1,309 open reading frames of which 496 are located on plasmids. With the exception of lp28-8, all linear replicons in their full length including their telomeres have been sequenced. The comparison with the genomes of the four other B. afzelii strains, ACA-1, PKo, HLJ01 and Tom3107, as well as the one of B. burgdorferi strain B31, confirmed a high degree of conservation within the linear chromosome of B. afzelii, whereas plasmid encoded genes showed a much larger diversity. Since some plasmids present in B. burgdorferi are missing in the B. afzelii genomes, the corresponding virulence factors of B. burgdorferi are found in B. afzelii on other unrelated plasmids. In addition, we have identified a species specific region in the circular plasmid, cp26, which could be used for species determination. Different non-coding RNAs have been located on the B. afzelii K78 genome, which have not previously been annotated in any of the published Borrelia genomes.

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.

Infections by HIV increase the risk of acquiring secondary viral and bacterial infections and methods are needed to determine the spectrum of co-infections for proper treatment. We used rolling circle amplification (RCA) and Ion Proton sequencing to investigate the vaginal microbiome of 20 HIV positive women from South Africa. A total of 46 different human papillomavirus (HPV) types were found, many of which are not detected by existing genotyping assays. Moreover, the complete genomes of two novel HPV types were determined. Abundance of HPV infections was highly correlated with real-time PCR estimates, indicating that the RCA-Proton method can be used for quantification of individual pathogens. We also identified a large number of other viral, bacterial and parasitic co-infections and the spectrum of these co-infections varied widely between individuals. Our method provides rapid detection of a broad range of pathogens and the ability to reconstruct complete genomes of novel infectious agents.

Dekkera bruxellensis can outcompete Saccharomyces cerevisiae in environments with low sugar concentrations. It is usually regarded as a spoilage yeast but has lately been identified as an alternative ethanol production organism. In this study, global gene expression in the industrial isolate D. bruxellensis CBS 11270 under oxygen and glucose limitation was investigated by whole transcriptome sequencing using the AB SOLiD technology. Among other observations, we noted expression of respiratory complex I NADH-ubiquinone reductase although D. bruxellensis is a Crabtree positive yeast. The observed higher expression of NADH-generating enzymes compared to NAD(+)-generating enzymes might be the reason for the previously observed NADH imbalance and resulting Custer effect in D. bruxellensis. Low expression of genes involved in glycerol production is probably the molecular basis for high efficiency of D. bruxellensis metabolism under nutrient limitation. No D. bruxellensis homologs to the genes involved in the final reactions of glycerol biosynthesis were detected. A high number of expressed sugar transporter genes is consistent with the hypothesis that the competitiveness of D. bruxellensis is due to a higher affinity for the limiting substrate.

Glioma is the most frequent malignant brain tumor in adults. Platelet-derived growth factor (PDGF) signaling is commonly activated in glioma. We have used a retrovirus-driven PDGFB-induced murine glioma model that causes tumors that closely resemble human gliomas of various grades. Knowing that retroviruses have a capacity to induce insertional mutagenesis, we have employed whole genome sequencing to identify potential genes that, together with PDGFB, drive glioma development. Gliomas were induced by RCAS virus injection into the brains of mice expressing the RCAS retroviral receptor from specific promoters. Genomic DNA from tumor cell lines was probed for retroviral tags and sequenced to identify genomic targets of the retrovirus. A streamlined analysis pipeline was developed for retrovirus integration detection and mapping to the reference mouse genome. Integration sites were analyzed and a common integration site (CIS) label was assigned to a gene, given that it was either tagged by a retrovirus more than once within a discovery set or found within the Retroviral Tagged Cancer Gene Database (RTCGD). In a small discovery subset of 15 murine gliomas, we have identified 40 CIS, of which 37 were validated by Sanger sequencing. When compared with previously identified CIS in RTCGD, 5.5% of them were shared with our older screen, where we overexpressed PDGFB from another retrovirus in order to induce glioma. Less CIS genes were shared with other published tumor models induced by viruses driven by other cancer genes/viruses. The majority of genes identified in our screen were tagged twice. However, Nfic, Cuecd1, Thra, Foxj1 and Nrxn1 were tagged three times, Ppfibp1 and Rhbg four times, and Mir29a/29b-1 seven times. As compared to control tumor lines, two top candidate genes, Mir29a and Ppfibp1, demonstrated significantly increased expression in tumor lines in were they were respectively tagged. Mir29a is often found downregulated in human tumors including gliomas, still high levels of Mir29a are sometimes found in certain aggressive cancers and in metastases. Interestingly, we found that specific PDGFR inhibition negatively regulates Mir29a, indicating a possible role for PDGF signaling in Mir29a regulation. Ppfibp1 has not been extensively studied in cancer. However, Ppfibp1 seems to have a subgroup-specific expression in human glioblastoma, making it an interesting candidate for further analysis. Here we present a new screening method that can be employed to identify genes involved in PDGFB-driven gliomagenesis. So far, we have identified 37 candidate genes by whole genome sequencing. Two of the most frequently tagged candidates, Mir29a and Ppfibp1 were upregulated as a consequence of retroviral mutagenesis. Their precise role in driving glioma formation in collaboration with PDGF is currently explored. Citation Format: Matko Cancer, Holger Weishaupt, Gabriela Rosen, Ignas Bunikis, Yiwen Jiang, Smitha Sreedharan, Sara Bolin, Ulf Gyllensten, Oren J. Becher, Lene Uhrbom, Adam Ameur, Fredrik J. Swartling. A forward genetics screen of murine brain tumors identifies novel candidate genes involved in gliomagenesis. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 2688.

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., Author Summary Lyme disease is caused by infection with the spirochete Borrelia burgdorferi. These spirochetes cycle between Ixodes ticks and vertebrate reservoirs, mainly rodents, but also birds. Previous studies have revealed major differences in the B. burgdorferi cell envelope structure and membrane composition compared to those of other bacteria. Proteins embedded in the bacterial membranes fulfill a number of tasks that are crucial for bacterial cells, such as solute and protein transport, as well as signal transduction, and interaction with other cells. Microorganisms have evolved mechanisms to protect themselves against harmful substances and secrete these through efflux pumps. So far, little is known about mechanisms of drug efflux systems in borreliae. Herein we identified an outer membrane channel forming protein important for B. burgdorferi to cause infection in mice and that also is involved in antibiotic resistance. We believe that this work will be helpful to understand the mechanisms underlying borreliae infection biology as well as in developing new therapeutic agents aiming to block this multi drug efflux system.

The genus Borrelia is the cause of the two human diseases: Lyme disease (LD) and relapsing fever (RF). Both LD and RF Borrelia species are obligate parasites and are dependent on nutrients provided by their hosts. The first step of nutrient uptake across the outer membrane of these Gram-negative bacteria is accomplished by water-filled channels, so-called porins. The knowledge of the porin composition in the outer membranes of the different pathogenic Borrelia species is limited. Only one porin has been described in relapsing fever spirochetes to date, whereas four porins are known to be present in Lyme disease agents. From these, the Borrelia burgdorferi outer membrane channel P66 is known to act as an adhesin and was well studied as a porin. To investigate if P66 porins are expressed and similarly capable of pore formation in other Borrelia causing Lyme disease or relapsing fever three LD species (B. burgdorferi, B. afzelii, B. garinii) and three RF species (B. duttonii, B. recurrentis and B. hermsii) were investigated for outer membrane proteins homologous to P66. A search in current published RF genomes, comprising the ones of B. duttonii, B. recurrentis and B. hermsii, indicated that they all contained P66 homologues. The P66 homologues of the six Borrelia species were purified to homogeneity and their pore-forming abilities as well as the biophysical properties of the pores were analyzed using the black lipid bilayer assay.