Targeted sequencing for high-resolution evolutionary analyses following genome duplication in salmonid fish : Proof of concept for key components of the insulin-like growth factor axis

Acknowledgements This study was funded by a Natural Environment Research Council grant (NERC, project code: NBAF704). FML is funded by a NERC Doctoral Training Grant (Project Reference: NE/L50175X/1). RLS was an undergraduate student at the University of Aberdeen and benefitted from financial support from the School of Biological Sciences. DJM is indebted to Dr. Steven Weiss (University of Graz, Austria), Dr. Takashi Yada (National Research Institute of Fisheries Science, Japan), Dr. Robert Devlin (Fisheries and Oceans Canada, Canada), Prof. Samuel Martin (University of Aberdeen, UK), Mr. Neil Lincoln (Environment Agency, UK) and Prof. Colin Adams/Mr. Stuart Wilson (University of Glasgow, UK) for providing salmonid material or assisting with its sampling. We are grateful to staff at the Centre for Genomics Research (University of Liverpool, UK) (i.e. NERC Biomolecular Analysis Facility – Liverpool; NBAF-Liverpool) for performing sequence capture/Illumina sequencing and providing us with details on associated methods that were incorporated into the manuscript. Finally, we are grateful to the organizers of the Society of Experimental Biology Satellite meeting 'Genome-powered perspectives in integrative physiology and evolutionary biology' (held in Prague, July 2015) for inviting us to contribute to this special edition of Marine Genomics and hosting a really stimulating meeting.

Combined analysis of variation in core, accessory and regulatory genome regions provides a super-resolution view into the evolution of bacterial populations

The use of whole-genome phylogenetic analysis has revolutionized our understanding of the evolution and spread of many important bacterial pathogens due to the high resolution view it provides. However, the majority of such analyses do not consider the potential role of accessory genes when inferring evolutionary trajectories. Moreover, the recently discovered importance of the switching of gene regulatory elements suggests that an exhaustive analysis, combining information from core and accessory genes with regulatory elements could provide unparalleled detail of the evolution of a bacterial population. Here we demonstrate this principle by applying it to a worldwide multi-host sample of the important pathogenic E. coli lineage ST131. Our approach reveals the existence of multiple circulating subtypes of the major drug–resistant clade of ST131 and provides the first ever population level evidence of core genome substitutions in gene regulatory regions associated with the acquisition and maintenance of different accessory genome elements.

Genome-scale analysis of the non-cultivable Treponema pallidum reveals extensive within-patient genetic variation

Insights into the genomic adaptive traits of Treponema pallidum, the causative bacterium of syphilis, have long been hampered due to the absence of in vitro culture models and the constraints associated with its propagation in rabbits. Here, we have bypassed the culture bottleneck by means of a targeted strategy never applied to uncultivable bacterial human pathogens to directly capture whole-genome T. pallidum data in the context of human infection. This strategy has unveiled a scenario of discreet T. pallidum interstrain single-nucleotide-polymorphism-based microevolution, contrasting with a rampant within-patient genetic heterogeneity mainly targeting multiple phase-variable loci and a major antigen-coding gene (tprK). TprK demonstrated remarkable variability and redundancy, intra- and interpatient, suggesting ongoing parallel adaptive diversification during human infection. Some bacterial functions (for example, flagella- and chemotaxis-associated) were systematically targeted by both inter- and intrastrain single nucleotide polymorphisms, as well as by ongoing within-patient phase variation events. Finally, patient-derived genomes possess mutations targeting a penicillin-binding protein coding gene (mrcA) that had never been reported, unveiling it as a candidate target to investigate the impact on the susceptibility to penicillin. Our findings decode the major genetic mechanisms by which T. pallidum promotes immune evasion and survival, and demonstrate the exceptional power of characterizing evolving pathogen subpopulations during human infection.

Evidence that human chlamydia pneumoniae was zoonotically acquired

Zoonotic infections are a growing threat to global health. Chlamydia pneumoniae is a major human pathogen that is widespread in human populations, causing acute respiratory disease, and has been associated with chronic disease. C. pneumoniae was first identified solely in human populations; however, its host range now includes other mammals, marsupials, amphibians, and reptiles. Australian koalas (Phascolarctos cinereus) are widely infected with two species of Chlamydia, C. pecorum and C. pneumoniae. Transmission of C. pneumoniae between animals and humans has not been reported; however, two other chlamydial species, C. psittaci and C. abortus, are known zoonotic pathogens. We have sequenced the 1,241,024-bp chromosome and a 7.5-kb cryptic chlamydial plasmid of the koala strain of C. pneumoniae (LPCoLN) using the whole-genome shotgun method. Comparative genomic analysis, including pseudogene and single-nucleotide polymorphism (SNP) distribution, and phylogenetic analysis of conserved genes and SNPs against the human isolates of C. pneumoniae show that the LPCoLN isolate is basal to human isolates. Thus, we propose based on compelling genomic and phylogenetic evidence that humans were originally infected zoonotically by an animal isolate(s) of C. pneumoniae which adapted to humans primarily through the processes of gene decay and plasmid loss, to the point where the animal reservoir is no longer required for transmission.

Phylogeographic reconstruction of a bacterial species with high levels of lateral gene transfer

Background Phylogeographic reconstruction of some bacterial populations is hindered by low diversity coupled with high levels of lateral gene transfer. A comparison of recombination levels and diversity at seven housekeeping genes for eleven bacterial species, most of which are commonly cited as having high levels of lateral gene transfer shows that the relative contributions of homologous recombination versus mutation for Burkholderia pseudomallei is over two times higher than for Streptococcus pneumoniae and is thus the highest value yet reported in bacteria. Despite the potential for homologous recombination to increase diversity, B. pseudomallei exhibits a relative lack of diversity at these loci. In these situations, whole genome genotyping of orthologous shared single nucleotide polymorphism loci, discovered using next generation sequencing technologies, can provide very large data sets capable of estimating core phylogenetic relationships. We compared and searched 43 whole genome sequences of B. pseudomallei and its closest relatives for single nucleotide polymorphisms in orthologous shared regions to use in phylogenetic reconstruction. Results Bayesian phylogenetic analyses of >14,000 single nucleotide polymorphisms yielded completely resolved trees for these 43 strains with high levels of statistical support. These results enable a better understanding of a separate analysis of population differentiation among >1,700 B. pseudomallei isolates as defined by sequence data from seven housekeeping genes. We analyzed this larger data set for population structure and allele sharing that can be attributed to lateral gene transfer. Our results suggest that despite an almost panmictic population, we can detect two distinct populations of B. pseudomallei that conform to biogeographic patterns found in many plant and animal species. That is, separation along Wallace's Line, a biogeographic boundary between Southeast Asia and Australia. Conclusion We describe an Australian origin for B. pseudomallei, characterized by a single introduction event into Southeast Asia during a recent glacial period, and variable levels of lateral gene transfer within populations. These patterns provide insights into mechanisms of genetic diversification in B. pseudomallei and its closest relatives, and provide a framework for integrating the traditionally separate fields of population genetics and phylogenetics for other bacterial species with high levels of lateral gene transfer.

Comparison of koala LPCoLN and human strains of Chlamydia pneumoniae highlights extended genetic diversity in the species

Background Chlamydia pneumoniae is a widespread pathogen causing upper and lower respiratory tract infections in addition to a range of other diseases in humans and animals. Previous whole genome analyses have focused on four essentially clonal (> 99% identity) C. pneumoniae human genomes (AR39, CWL029, J138 and TW183), providing relatively little insight into strain diversity and evolution of this species. Results We performed individual gene-by-gene comparisons of the recently sequenced C. pneumoniae koala genome and four C. pneumoniae human genomes to identify species-specific genes, and more importantly, to gain an insight into the genetic diversity and evolution of the species. We selected genes dispersed throughout the chromosome, representing genes that were specific to C. pneumoniae, genes with a demonstrated role in chlamydial biology and/or pathogenicity (n = 49), genes encoding nucleotide salvage or amino acid biosynthesis proteins (n = 6), and extrachromosomal elements (9 plasmid and 2 bacteriophage genes). Conclusions We have identified strain-specific differences and targets for detection of C. pneumoniae isolates from both human and animal origin. Such characterisation is necessary for an improved understanding of disease transmission and intervention.

Modulation of the Chlamydia trachomatis In vitro transcriptome response by the sex hormones estradiol and progesterone

BACKGROUND:Chlamydia trachomatis is a major cause of sexually transmitted disease in humans. Previous studies in both humans and animal models of chlamydial genital tract infection have suggested that the hormonal status of the genital tract epithelium at the time of exposure can influence the outcome of the chlamydial infection. We performed a whole genome transcriptional profiling study of C. trachomatis infection in ECC-1 cells under progesterone or estradiol treatment.RESULTS:Both hormone treatments caused a significant shift in the sub-set of genes expressed (25% of the transcriptome altered by more than 2-fold). Overall, estradiol treatment resulted in the down-regulation of 151 genes, including those associated with lipid and nucleotide metabolism. Of particular interest was the up-regulation in estradiol-supplemented cultures of six genes (omcB, trpB, cydA, cydB, pyk and yggV), which suggest a stress response similar to that reported previously in other models of chlamydial persistence. We also observed morphological changes consistent with a persistence response. By comparison, progesterone supplementation resulted in a general up-regulation of an energy utilising response.CONCLUSION:Our data shows for the first time, that the treatment of chlamydial host cells with key reproductive hormones such as progesterone and estradiol, results in significantly altered chlamydial gene expression profiles. It is likely that these chlamydial expression patterns are survival responses, evolved by the pathogen to enable it to overcome the host's innate immune response. The induction of chlamydial persistence is probably a key component of this survival response.

Differential expression of chemokine and matrix re-modelling genes is associated with contrasting schistosome-induced hepatopathology in murine models

The pathological outcomes of schistosomiasis are largely dependent on the molecular and cellular mechanisms of the host immune response. In this study, we investigated the contribution of variations in host gene expression to the contrasting hepatic pathology observed between two inbred mouse strains following Schistosoma japonicum infection. Whole genome microarray analysis was employed in conjunction with histological and immunohistochemical analysis to define and compare the hepatic gene expression profiles and cellular composition associated with the hepatopathology observed in S. japonicum-infected BALB/c and CBA mice. We show that the transcriptional profiles differ significantly between the two mouse strains with high statistical confidence. We identified specific genes correlating with the more severe pathology associated with CBA mice, as well as genes which may confer the milder degree of pathology associated with BALB/c mice. In BALB/c mice, neutrophil genes exhibited striking increases in expression, which coincided with the significantly greater accumulation of neutrophils at granulomatous regions seen in histological sections of hepatic tissue. In contrast, up-regulated expression of the eosinophil chemokine CCL24 in CBA mice paralleled the cellular influx of eosinophils to the hepatic granulomas. Additionally, there was greater down-regulation of genes involved in metabolic processes in CBA mice, reflecting the more pronounced hepatic damage in these mice. Profibrotic genes showed similar levels of expression in both mouse strains, as did genes associated with Th1 and Th2 responses. However, imbalances in expression of matrix metalloproteinases (e.g. MMP12, MMP13) and tissue inhibitors of metalloproteinases (TIMP1) may contribute to the contrasting pathology observed in the two strains. Overall, these results provide a more complete picture of the molecular and cellular mechanisms which govern the pathological outcome of hepatic schistosomiasis. This improved understanding of the immunopathogenesis in the murine model schistosomiasis provides the basis for a better appreciation of the complexities associated with chronic human schistosomiasis.

Gene expression profile of primary prostate epithelial and stromal cells in response to sulforaphane or iberin exposure

BACKGROUND: Broccoli consumption has been associated with a reduced risk of prostate cancer. Isothiocyanates (ITCs) derived from glucosinolates that accumulate in broccoli are dietary compounds that may mediate these health effects. Sulforaphane (SF, 4-methylsulphinylbutyl ITC) derives from heading broccoli (calabrese) and iberin (IB, 3-methylsulphinypropyl ITC) from sprouting broccoli. While there are many studies regarding the biological activity of SF, mainly undertaken with cancerous cells, there are few studies associated with IB. METHODS: Primary epithelial and stromal cells were derived from benign prostatic hyperplasia tissue. Affymetrix U133 Plus 2.0 whole genome arrays were used to compare global gene expression between these cells, and to quantify changes in gene expression following exposure to physiologically appropriate concentrations of SF and IB. Ontology and pathway analyses were used to interpret results. Changes in expression of a subset of genes were confirmed by real-time RT-PCR. RESULTS: Global gene expression profiling identified epithelial and stromal-specific gene expression profiles. SF induced more changes in epithelial cells, whereas IB was more effective in stromal cells. Although IB and SF induced different changes in gene expression in both epithelial and stromal cells, these were associated with similar pathways, such as cell cycle and detoxification. Both ITCs increased expression of PLAGL1, a tumor suppressor gene, in stromal cells and suppressed expression of the putative tumor promoting genes IFITM1, CSPG2, and VIM in epithelial cells. CONCLUSION: These data suggest that IB and SF both alter genes associated with cancer prevention, and IB should be investigated further as a potential chemopreventative agent.

Lineage-specific biomarkers predict response to FGFR inhibition

In this issue of Cancer Discovery, Guagnano and colleagues use a large and diverse annotated collection of cancer cell lines, the Cancer Cell Line Encyclopedia, to correlate whole-genome expression and genomic alteration datasets with cell line sensitivity data to the novel pan-fibroblast growth factor receptor (FGFR) inhibitor NVP-BGJ398. Their findings underscore not only the preclinical use of such cell line panels in identifying predictive biomarkers, but also the emergence of the FGFRs as valid therapeutic targets, across an increasingly broad range of malignancies.

Chlamydia pneumoniae : modern insights into an ancient pathogen

Chlamydia pneumoniae is an enigmatic human and animal pathogen. Originally discovered in association with acute human respiratory disease, it is now associated with a remarkably wide range of chronic diseases as well as having a cosmopolitan distribution within the animal kingdom. Molecular typing studies suggest that animal strains are ancestral to human strains and that C. pneumoniae crossed from animals to humans as the result of at least one relatively recent zoonotic event. Whole genome analyses appear to support this concept – the human strains are highly conserved whereas the single animal strain that has been fully sequenced has a larger genome with several notable differences. When compared to the other, better known chlamydial species that is implicated in human infection, Chlamydia trachomatis, C. pneumoniae demonstrates pertinent differences in its cell biology, development, and genome structure. Here, we examine the characteristic facets of C. pneumoniae biology, offering insights into the diversity and evolution of this silent and ancient pathogen.

Mapping of multiple quantitative trait loci affecting bovine spongiform encephalopathy

A whole-genome scan was conducted to map quantitative trait loci (QTL) for BSE resistance or susceptibility. Cows from four half-sib families were included and 173 microsatellite markers were used to construct a 2835-cM (Kosambi) linkage map covering 29 autosomes and the pseudoautosomal region of the sex chromosome. Interval mapping by linear regression was applied and extended to a multiple-QTL analysis approach that used identified QTL on other chromosomes as cofactors to increase mapping power. In the multiple-QTL analysis, two genome-wide significant QTL (BTA17 and X/Y ps) and four genome-wide suggestive QTL (BTA1, 6, 13, and 19) were revealed. The QTL identified here using linkage analysis do not overlap with regions previously identified using TDT analysis. One factor that may explain the disparity between the results is that a more extensive data set was used in the present study. Furthermore, methodological differences between TDT and linkage analyses may affect the power of these approaches.

Elucidating the host–pathogen interaction between human colorectal cells and invading Enterovirus 71 using transcriptomics profiling

Enterovirus 71 (EV71) is one of the main etiological agents for Hand, Foot and Mouth Disease (HFMD) and has been shown to be associated with severe clinical manifestation. Currently, there is no antiviral therapeutic for the treatment of HFMD patients owing to a lack of understanding of EV71 pathogenesis. This study seeks to elucidate the transcriptomic changes that result from EV71 infection. Human whole genome microarray was employed to monitor changes in genomic profiles between infected and uninfected cells. The results reveal altered expression of human genes involved in critical pathways including the immune response and the stress response. Together, data from this study provide valuable insights into the host–pathogen interaction between human colorectal cells and EV71.

An ancient duplication of apple MYB transcription factors is responsible for novel red fruit-flesh phenotyp

Anthocyanin accumulation is coordinated in plants by a number of conserved transcription factors. In apple (Malus × domestica), an R2R3 MYB transcription factor has been shown to control fruit flesh and foliage anthocyanin pigmentation (MYB10) and fruit skin color (MYB1). However, the pattern of expression and allelic variation at these loci does not explain all anthocyanin-related apple phenotypes. One such example is an open-pollinated seedling of cv Sangrado that has green foliage and develops red flesh in the fruit cortex late in maturity. We used methods that combine plant breeding, molecular biology, and genomics to identify duplicated MYB transcription factors that could control this phenotype. We then demonstrated that the red-flesh cortex phenotype is associated with enhanced expression of MYB110a, a paralog of MYB10. Functional characterization of MYB110a showed that it was able to up-regulate anthocyanin biosynthesis in tobacco (Nicotiana tabacum). The chromosomal location of MYB110a is consistent with a whole-genome duplication event that occurred during the evolution of apple within the Maloideae family. Both MYB10 and MYB110a have conserved function in some cultivars, but they differ in their expression pattern and response to fruit maturity.

Effect of 5′UTR introns on gene expression in Arabidopsis thaliana

Background The majority of introns in gene transcripts are found within the coding sequences (CDSs). A small but significant fraction of introns are also found to reside within the untranslated regions (5′UTRs and 3′UTRs) of expressed sequences. Alignment of the whole genome and expressed sequence tags (ESTs) of the model plant Arabidopsis thaliana has identified introns residing in both coding and non-coding regions of the genome. Results A bioinformatic analysis revealed some interesting observations: (1) the density of introns in 5′UTRs is similar to that in CDSs but much higher than that in 3′UTRs; (2) the 5′UTR introns are preferentially located close to the initiating ATG codon; (3) introns in the 5′UTRs are, on average, longer than introns in the CDSs and 3′UTRs; and (4) 5′UTR introns have a different nucleotide composition to that of CDs and 3′UTR introns. Furthermore, we show that the 5′UTR intron of the A. thaliana EFIα-A3 gene affects the gene expression and the size of the 5′UTR intron influences the level of gene expression. Conclusion Introns within the 5′UTR show specific features that distinguish them from introns that reside within the coding sequence and the 3′UTR. In the EFIα-A3 gene, the presence of a long intron in the 5′UTR is sufficient to enhance gene expression in plants in a size dependent manner.