Ecological fitness, genomic islands and bacterial pathogenicity. A Darwinian view of the evolution of microbes

The compositions of bacterial genomes can be changed rapidly and dramatically through a variety of processes including horizontal gene transfer. This form of change is key to bacterial evolution, as it leads to ‘evolution in quantum leaps’. Horizontal gene transfer entails the incorporation of genetic elements transferred from another organism—perhaps in an earlier generation—directly into the genome, where they form ‘genomic islands’, i.e. blocks of DNA with signatures of mobile genetic elements. Genomic islands whose functions increase bacterial fitness, either directly or indirectly, have most likely been positively selected and can be termed ‘fitness islands’. Fitness islands can be divided into several subtypes: ‘ecological islands’ in environmental bacteria and ‘saprophytic islands’, ‘symbiosis islands’ or ‘pathogenicity islands’ (PAIs) in microorganisms that interact with living hosts. Here we discuss ways in which PAIs contribute to the pathogenic potency of bacteria, and the idea that genetic entities similar to genomic islands may also be present in the genomes of eukaryotes.     

HPV-CCDC106 integration alters local chromosome architecture and hijacks an enhancer by three-dimensional genome structure remodeling in cervical cancer

Integration of human papillomavirus (HPV) DNA into the human genome is a reputed key driver of cervical cancer. However, the effects of HPV integration on chromatin structural organization and gene expression are largely unknown. We studied a cohort of 61 samples and identified an integration hot spot in the CCDC106 gene on chromosome 19. We then selected fresh cancer tissue that contained the unique integration loci at CCDC106 with no HPV episomal DNA and performed whole-genome, RNA, chromatin immunoprecipitation and high-throughput chromosome conformation capture (Hi-C) sequencing to identify the mechanisms of HPV integration in cervical carcinogenesis. Molecular analyses indicated that chromosome 19 exhibited significant genomic variation and differential expression densities, with correlation found between three-dimensional (3D) structural change and gene expression. Importantly, HPV integration divided one topologically associated domain (TAD) into two smaller TADs and hijacked an enhancer from PEG3 to CCDC106, with a decrease in PEG3 expression and an increase in CCDC106 expression. This expression dysregulation was further confirmed using 10 samples from our cohort, which exhibited the same HPV-CCDC106 integration. In summary, we found that HPV-CCDC106 integration altered local chromosome architecture and hijacked an enhancer via 3D genome structure remodeling. Thus, this study provides insight into the 3D structural mechanism underlying HPV integration in cervical carcinogenesis.     

Distinct Gene Number-Genome Size Relationships for Eukaryotes and Non-Eukaryotes: Gene Content Estimation for Dinoflagellate Genomes

The ability to predict gene content is highly desirable for characterization of not-yet sequenced genomes like those of dinoflagellates. Using data from completely sequenced and annotated genomes from phylogenetically diverse lineages, we investigated the relationship between gene content and genome size using regression analyses. Distinct relationships between log₁₀-transformed protein-coding gene number (Y′) versus log₁₀-transformed genome size (X′, genome size in kbp) were found for eukaryotes and non-eukaryotes. Eukaryotes best fit a logarithmic model, Y′ = ln(-46.200+22.678X′, whereas non-eukaryotes a linear model, Y′ = 0.045+0.977X′, both with high significance (p<0.001, R²>0.91). Total gene number shows similar trends in both groups to their respective protein coding regressions. The distinct correlations reflect lower and decreasing gene-coding percentages as genome size increases in eukaryotes (82%–1%) compared to higher and relatively stable percentages in prokaryotes and viruses (97%–47%). The eukaryotic regression models project that the smallest dinoflagellate genome (3×10⁶ kbp) contains 38,188 protein-coding (40,086 total) genes and the largest (245×10⁶ kbp) 87,688 protein-coding (92,013 total) genes, corresponding to 1.8% and 0.05% gene-coding percentages. These estimates do not likely represent extraordinarily high functional diversity of the encoded proteome but rather highly redundant genomes as evidenced by high gene copy numbers documented for various dinoflagellate species.     

Local hopping mobile DNA implicated in pseudogene formation and reductive evolution in an obligate cyanobacteria-plant symbiosis

Background

Insertion sequences (ISs) are approximately 1 kbp long “jumping” genes found in prokaryotes. ISs encode the protein Transposase, which facilitates the excision and reinsertion of ISs in genomes, making these sequences a type of class I (“cut-and-paste”) Mobile Genetic Elements. ISs are proposed to be involved in the reductive evolution of symbiotic prokaryotes. Our previous sequencing of the genome of the cyanobacterium ‘Nostoc azollae’ 0708, living in a tight perpetual symbiotic association with a plant (the water fern Azolla), revealed the presence of an eroding genome, with a high number of insertion sequences (ISs) together with an unprecedented large proportion of pseudogenes. To investigate the role of ISs in the reductive evolution of ‘Nostoc azollae’ 0708, and potentially in the formation of pseudogenes, a bioinformatic investigation of the IS identities and positions in 47 cyanobacterial genomes was conducted. To widen the scope, the IS contents were analysed qualitatively and quantitatively in 20 other genomes representing both free-living and symbiotic bacteria.

Results

Insertion Sequences were not randomly distributed in the bacterial genomes and were found to transpose short distances from their original location (“local hopping”) and pseudogenes were enriched in the vicinity of IS elements. In general, symbiotic organisms showed higher densities of IS elements and pseudogenes than non-symbiotic bacteria. A total of 1108 distinct repeated sequences over 500 bp were identified in the 67 genomes investigated. In the genome of ‘Nostoc azollae’ 0708, IS elements were apparent at 970 locations (14.3%), with 428 being full-length. Morphologically complex cyanobacteria with large genomes showed higher frequencies of IS elements, irrespective of life style.

Conclusions

The apparent co-location of IS elements and pseudogenes found in prokaryotic genomes implies earlier IS transpositions into genes. As transpositions tend to be local rather than genome wide this likely explains the proximity between IS elements and pseudogenes. These findings suggest that ISs facilitate the reductive evolution in for instance in the symbiotic cyanobacterium ‘Nostoc azollae’ 0708 and in other obligate prokaryotic symbionts.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1386-7) contains supplementary material, which is available to authorized users.     

Human papillomavirus type 16 (HPV-16) genomes integrated in head and neck cancers and in HPV-16-immortalized human keratinocyte clones express chimeric virus-cell mRNAs similar to those found in cervical cancers

Many human papillomavirus (HPV)-positive high-grade lesions and cancers of the uterine cervix harbor integrated HPV genomes expressing the E6 and E7 oncogenes from chimeric virus-cell mRNAs, but less is known about HPV integration in head and neck cancer (HNC). Here we compared viral DNA status and E6-E7 mRNA sequences in HPV-16-positive HNC tumors to those in independent human keratinocyte cell clones derived from primary tonsillar or foreskin epithelia immortalized with HPV-16 genomes. Three of nine HNC tumors and epithelial clones containing unintegrated HPV-16 genomes expressed mRNAs spliced from HPV-16 SD880 to SA3358 and terminating at the viral early gene p(A) signal. In contrast, most integrated HPV genomes in six HNCs and a set of 31 keratinocyte clones expressed HPV-16 major early promoter (MEP)-initiated mRNAs spliced from viral SD880 directly to diverse cellular sequences, with a minority spliced to SA3358 followed by a cellular DNA junction. Sequence analysis of chimeric virus-cell mRNAs from HNC tumors and keratinocyte clones identified viral integration sites in a variety of chromosomes, with some located in or near growth control genes, including the c-myc protooncogene and the gene encoding FAP-1 phosphatase. Taken together, these findings support the hypothesis that HPV integration in cancers is a stochastic process resulting in clonal selection of aggressively expanding cells with altered gene expression of integrated HPV genomes and potential perturbations of cellular genes at or near viral integration sites. Furthermore, our results demonstrate that this selection also takes place and can be studied in primary human keratinocytes in culture.     

Differential Responses of CO2 Assimilation,Carbohydrate Allocation and Gene Expression to NaCl Stress in Perennial Ryegrass with Different Salt Tolerance

Little is known about the effects of NaCl stress on perennial ryegrass (Lolium perenne L.) photosynthesis and carbohydrate flux. The objective of this study was to understand the carbohydrate metabolism and identify the gene expression affected by salinity stress. Seventy-four days old seedlings of two perennial ryegrass accessions (salt-sensitive ‘PI 538976’ and salt-tolerant ‘Overdrive’) were subjected to three levels of salinity stress for 5 days. Turf quality in all tissues (leaves, stems and roots) of both grass accessions negatively and significantly correlated with GFS (Glu+Fru+Suc) content, except for ‘Overdrive’ stems. Relative growth rate (RGR) in leaves negatively and significantly correlated with GFS content in ‘Overdrive’ (P<0.01) and ‘PI 538976’ (P<0.05) under salt stress. ‘Overdrive’ had higher CO₂ assimilation and F_v/F_m than ‘PI 538976’. Intercellular CO₂ concentration, however, was higher in ‘PI 538976’ treated with 400 mM NaCl relative to that with 200 mM NaCl. GFS content negatively and significantly correlated with RGR in ‘Overdrive’ and ‘PI 538976’ leaves and in ‘PI 538976’ stems and roots under salt stress. In leaves, carbohydrate allocation negatively and significantly correlated with RGR (r² = 0.83, P<0.01) and turf quality (r² = 0.88, P<0.01) in salt-tolerant ‘Overdrive’, however, the opposite trend for salt-sensitive ‘PI 538976’ (r² = 0.71, P<0.05 for RGR; r² = 0.62, P>0.05 for turf quality). A greater up-regulation in the expression of SPS, SS, SI, 6-SFT gene was observed in ‘Overdrive’ than ‘PI 538976’. A higher level of SPS and SS expression in leaves was found in ‘PI 538976’ relative to ‘Overdrive’. Accumulation of hexoses in roots, stems and leaves can induce a feedback repression to photosynthesis in salt-stressed perennial ryegrass and the salt tolerance may be changed with the carbohydrate allocation in leaves and stems.     

In silico prediction of HIV-1-host molecular interactions and their directionality

Human immunodeficiency virus type 1 (HIV-1) continues to be a major cause of disease and premature death. As with all viruses, HIV-1 exploits a host cell to replicate. Improving our understanding of the molecular interactions between virus and human host proteins is crucial for a mechanistic understanding of virus biology, infection and host antiviral activities. This knowledge will potentially permit the identification of host molecules for targeting by drugs with antiviral properties. Here, we propose a data-driven approach for the analysis and prediction of the HIV-1 interacting proteins (VIPs) with a focus on the directionality of the interaction: host-dependency versus antiviral factors. Using support vector machine learning models and features encompassing genetic, proteomic and network properties, our results reveal some significant differences between the VIPs and non-HIV-1 interacting human proteins (non-VIPs). As assessed by comparison with the HIV-1 infection pathway data in the Reactome database (sensitivity > 90%, threshold = 0.5), we demonstrate these models have good generalization properties. We find that the ‘direction’ of the HIV-1-host molecular interactions is also predictable due to different characteristics of ‘forward’/pro-viral versus ‘backward’/pro-host proteins. Additionally, we infer the previously unknown direction of the interactions between HIV-1 and 1351 human host proteins. A web server for performing predictions is available at http://hivpre.cvr.gla.ac.uk/.     

sPARTA: a parallelized pipeline for integrated analysis of plant miRNA and cleaved mRNA data sets,including new miRNA target-identification software

Parallel analysis of RNA ends (PARE) is a technique utilizing high-throughput sequencing to profile uncapped, mRNA cleavage or decay products on a genome-wide basis. Tools currently available to validate miRNA targets using PARE data employ only annotated genes, whereas important targets may be found in unannotated genomic regions. To handle such cases and to scale to the growing availability of PARE data and genomes, we developed a new tool, ‘sPARTA’ (small RNA-PARE target analyzer) that utilizes a built-in, plant-focused target prediction module (aka ‘miRferno’). sPARTA not only exhibits an unprecedented gain in speed but also it shows greater predictive power by validating more targets, compared to a popular alternative. In addition, the novel ‘seed-free’ mode, optimized to find targets irrespective of complementarity in the seed-region, identifies novel intergenic targets. To fully capitalize on the novelty and strengths of sPARTA, we developed a web resource, ‘comPARE’, for plant miRNA target analysis; this facilitates the systematic identification and analysis of miRNA-target interactions across multiple species, integrated with visualization tools. This collation of high-throughput small RNA and PARE datasets from different genomes further facilitates re-evaluation of existing miRNA annotations, resulting in a ‘cleaner’ set of microRNAs.     

Conjugative plasmids: vessels of the communal gene pool

Comparative whole-genome analyses have demonstrated that horizontal gene transfer (HGT) provides a significant contribution to prokaryotic genome innovation. The evolution of specific prokaryotes is therefore tightly linked to the environment in which they live and the communal pool of genes available within that environment. Here we use the term supergenome to describe the set of all genes that a prokaryotic ‘individual’ can draw on within a particular environmental setting. Conjugative plasmids can be considered particularly successful entities within the communal pool, which have enabled HGT over large taxonomic distances. These plasmids are collections of discrete regions of genes that function as ‘backbone modules’ to undertake different aspects of overall plasmid maintenance and propagation. Conjugative plasmids often carry suites of ‘accessory elements’ that contribute adaptive traits to the hosts and, potentially, other resident prokaryotes within specific environmental niches. Insight into the evolution of plasmid modules therefore contributes to our knowledge of gene dissemination and evolution within prokaryotic communities. This communal pool provides the prokaryotes with an important mechanistic framework for obtaining adaptability and functional diversity that alleviates the need for large genomes of specialized ‘private genes’.     

Recombinant viruses as transformation vectors of marine macroalgae

The large dsDNA viruses that are known to infect eukaryotic algae show promise as genetic vectors for algal biotechnology. The large size (150–330 kbp) of these viral genomes may permit insertion of large sequences of foreign DNA. The viruses infecting filamentous marine brown algae appear to be integrated into the genomes of their hosts, and may provide integration mechanisms that can be used for directing insertion of foreign genes into algal chromosomes.     

High-throughput,cost-effective verification of structural DNA assembly

DNA ‘assembly’ from ‘building blocks’ remains a cornerstone in synthetic biology, whether it be for gene synthesis (∼1 kb), pathway engineering (∼10 kb) or synthetic genomes (>100 kb). Despite numerous advances in the techniques used for DNA assembly, verification of the assembly is still a necessity, which becomes cost-prohibitive and a logistical challenge with increasing scale. Here we describe for the first time a comprehensive, high-throughput solution for structural DNA assembly verification by restriction digest using exhaustive in silico enzyme screening, rolling circle amplification of plasmid DNA, capillary electrophoresis and automated digest pattern recognition. This low-cost and robust methodology has been successfully used to screen over 31 000 clones of DNA constructs at <$1 per sample.     

Recognizing the pseudogenes in bacterial genomes

Pseudogenes are now known to be a regular feature of bacterial genomes and are found in particularly high numbers within the genomes of recently emerged bacterial pathogens. As most pseudogenes are recognized by sequence alignments, we use newly available genomic sequences to identify the pseudogenes in 11 genomes from 4 bacterial genera, each of which contains at least 1 human pathogen. The numbers of pseudogenes range from 27 in Staphylococcus aureus MW2 to 337 in Yersinia pestis CO92 (e.g. 1–8% of the annotated genes in the genome). Most pseudogenes are formed by small frameshifting indels, but because stop codons are A + T-rich, the two low-G + C Gram-positive taxa (Streptococcus and Staphylococcus) have relatively high fractions of pseudogenes generated by nonsense mutations when compared with more G + C-rich genomes. Over half of the pseudogenes are produced from genes whose original functions were annotated as ‘hypothetical’ or ‘unknown’; however, several broadly distributed genes involved in nucleotide processing, repair or replication have become pseudogenes in one of the sequenced Vibrio vulnificus genomes. Although many of our comparisons involved closely related strains with broadly overlapping gene inventories, each genome contains a largely unique set of pseudogenes, suggesting that pseudogenes are formed and eliminated relatively rapidly from most bacterial genomes.     

Bacteriophage Mu integration in yeast and mammalian genomes

Genomic parasites have evolved distinctive lifestyles to optimize replication in the context of the genomes they inhabit. Here, we introduced new DNA into eukaryotic cells using bacteriophage Mu DNA transposition complexes, termed ‘transpososomes’. Following electroporation of transpososomes and selection for marker gene expression, efficient integration was verified in yeast, mouse and human genomes. Although Mu has evolved in prokaryotes, strong biases were seen in the target site distributions in eukaryotic genomes, and these biases differed between yeast and mammals. In Saccharomyces cerevisiae transposons accumulated outside of genes, consistent with selection against gene disruption. In mouse and human cells, transposons accumulated within genes, which previous work suggests is a favorable location for efficient expression of selectable markers. Naturally occurring transposons and viruses in yeast and mammals show related, but more extreme, targeting biases, suggesting that they are responding to the same pressures. These data help clarify the constraints exerted by genome structure on genomic parasites, and illustrate the wide utility of the Mu transpososome technology for gene transfer in eukaryotic cells.     

Integration sites of human papillomavirus 18 DNA sequences on HeLa cell chromosomes

Human papillomaviruses (HPV) 16 and 18 are closely linked with human genital cancer. In most cervical carcinomas, viral sequences are integrated into the host genome. HeLa, a cervical carcinoma cell line, has multiple copies of integrated HPV 18 DNA. In this study, in situ chromosome hybridization was used to assign the integration sites of HPV 18 DNA sequences on HeLa cell chromosomes. Four sites of hybridization were identified at 8q23----q24, 9q31----q34, p11----p13 on an abnormal chromosome 5, and q12----q13 on an abnormal 22. Three of these sites correspond with the locations of MYC, ABL, and SIS protooncogenes, and are at or in close proximity to fragile sites. The chromosomal localization of HPV 18 DNA may be useful in assessing the role of viral integration in the development of this malignancy.     

A new approach for the identification and cloning of genes: the pBACwich system using Cre/lox site-specific recombination

With current plant transformation methods (Agrobacterium, biolistics and protoplast fusion), insertion of DNA into the genome occurs randomly and in many instances at multiple sites. Associated position effects, copy number differences and multigene interactions can make gene expression experiments difficult to interpret and plant phenotypes less predictable. An alternative approach to random integration of large DNA fragments into plants is to utilize one of several site-specific recombination (SSR) systems, such as Cre/lox. Cre has been shown in numerous instances to mediate lox site-specific recombination in animal and plant cells. By incorporating the Cre/lox SSR system into a bacterial artificial chromosome (BAC) vector, a more precise evaluation of large DNA inserts for genetic complementation should be possible. Site-specific insertion of DNA into predefined sites in the genome may eliminate unwanted ‘position effects’ caused by the random integration of exogenously introduced DNA. In an effort to make the Cre/lox system an effective tool for site-directed integration of large DNAs, we constructed and tested a new vector potentially capable of integrating large DNA inserts into plant and fungal genomes. In this study, we present the construction of a new BAC vector, pBACwich, for the system and the use of this vector to demonstrate SSR of large DNA inserts (up to 230 kb) into plant and fungal genomes.