Search

Wolbachia is one of the most widespread intracellular bacteria on earth, estimated to infect between 40 and 66% of arthropod species in most ecosystems that have been surveyed. Their significance rests not only in their vast distribution, but also in their ability to modify the reproductive biology of their hosts, which can ultimately affect genetic diversity and speciation of infected populations. Wolbachia has yet to be formally identified in the fauna of New Zealand which has high levels of endemic biodiversity and this represents a gap in our understanding of the global biology of Wolbachia. Using High Throughput Sequencing (HTS) of host DNA in conjunction with traditional molecular techniques we identified six endemic Orthoptera species that were positive for Wolbachia infection. In addition, short-sequence amplification with Wolbachia specific primers applied to New Zealand and introduced invertebrates detected a further 153 individuals positive for Wolbachia. From these short-range DNA amplification products sequence data was obtained for the ftsZ gene region from 86 individuals representing 10 host species. Phylogenetic analysis using the sequences obtained in this study reveals that there are two distinct Wolbachia bacteria lineages in New Zealand hosts belonging to recognised Wolbachia supergroups (A and B). These represent the first described instances of Wolbachia in the New Zealand native fauna, including detection in putative parasitoids of infected Orthoptera suggesting a possible transmission path. Our detection of Wolbachia infections of New Zealand species provides the opportunity to study local transmission of Wolbachia and explore their role in the evolution of New Zealand invertebrates.

Phytoplasmas are cell wall-less plant pathogenic bacteria responsible for major crop losses throughout the world. In grapevine they cause grapevine yellows, a detrimental disease associated with a variety of symptoms. The high economic impact of this disease has sparked considerable interest among researchers to understand molecular mechanisms related to pathogenesis. Increasing evidence exist that a class of small non-coding endogenous RNAs, known as microRNAs (miRNAs), play an important role in post-transcriptional gene regulation during plant development and responses to biotic and abiotic stresses. Thus, we aimed to dissect complex high-throughput small RNA sequencing data for the genome-wide identification of known and novel differentially expressed miRNAs, using read libraries constructed from healthy and phytoplasma-infected Chardonnay leaf material. Furthermore, we utilised computational resources to predict putative miRNA targets to explore the involvement of possible pathogen response pathways. We identified multiple known miRNA sequence variants (isomiRs), likely generated through post-transcriptional modifications. Sequences of 13 known, canonical miRNAs were shown to be differentially expressed. A total of 175 novel miRNA precursor sequences, each derived from a unique genomic location, were predicted, of which 23 were differentially expressed. A homology search revealed that some of these novel miRNAs shared high sequence similarity with conserved miRNAs from other plant species, as well as known grapevine miRNAs. The relative expression of randomly selected known and novel miRNAs was determined with real-time RT-qPCR analysis, thereby validating the trend of expression seen in the normalised small RNA sequencing read count data. Among the putative miRNA targets, we identified genes involved in plant morphology, hormone signalling, nutrient homeostasis, as well as plant stress. Our results may assist in understanding the role that miRNA pathways play during plant pathogenesis, and may be crucial in understanding disease symptom development in aster yellows phytoplasma-infected grapevines.

Hypotheses of hybrid origin are common. Here we use next generation sequencing to test a hybrid hypothesis for a non-model insect with a large genome. We compared a putative hybrid triploid stick insect species (Acanthoxyla geisovii) with its putative paternal diploid taxon (Clitarchus hookeri), a relationship that provides clear predictions for the relative genetic diversity within each genome. The parental taxon is expected to have comparatively low allelic diversity that is nested within the diversity of the hybrid daughter genome. The scale of genome sequencing required was conveniently achieved by extracting mRNA and sequencing cDNA to examine expressed allelic diversity. This allowed us to test hybrid-progenitor relationships among non-model organisms with large genomes and different ploidy levels. Examination of thousands of independent loci avoids potential problems produced by the silencing of parts of one or other of the parental genomes, a phenomenon sometimes associated with the process of stabilisation of a hybrid genome. Transcript assembles were assessed for evidence of paralogs and/or alternative splice variants before proceeding. Comparison of transcript assemblies was not an appropriate measure of genetic variability, but by mapping reads back to clusters derived from each species we determined levels of allelic diversity. We found greater cDNA sequence diversity among alleles in the putative hybrid species (Acanthoxyla geisovii) than the non-hybrid. The allelic diversity within the putative paternal species (Clitachus hookeri) nested within the hybrid-daughter genome, supports the current view of a hybrid-progenitor relationship for these stick insect species. Next generation sequencing technology provides opportunities for testing evolutionary hypotheses with non-model organisms, including, as here, genomes that are large due to polyploidy.

Background The late-gestation fetal sheep responds to hypoxia with physiological, neuroendocrine, and cellular responses that aid in fetal survival. The response of the fetus to hypoxia represents a coordinated effort to maximize oxygen transfer from the mother and minimize wasteful oxygen consumption by the fetus. While there have been many studies aimed at investigating the coordinated physiological and endocrine responses to hypoxia, and while immunohistochemical or in situ hybridization studies have revealed pathways supporting the endocrine function of the pituitary, there is little known about the coordinated cellular response of the pituitary to the hypoxia. Results Thirty min hypoxia (from 17.0±1.7 to 8.0±0.8 mm Hg, followed by 30 min normoxia) upregulated 595 and downregulated 790 genes in fetal pituitary (123-132 days' gestation; term = 147 days). Network inference of up- and down- regulated genes revealed a high degree of functional relatedness amongst the gene sets. Gene ontology analysis revealed upregulation of cellular metabolic processes (e.g., RNA synthesis, response to estrogens) and downregulation of protein phosphorylation, protein metabolism, and mitosis. Genes found to be at the center of the network of upregulated genes included genes important for purine binding and signaling. At the center of the downregulated network were genes involved in mRNA processing, DNA repair, sumoylation, and vesicular trafficking. Transcription factor analysis revealed that both up- and down-regulated gene sets are enriched for control by several transcription factors (e.g., SP1, MAZ, LEF1, NRF1, ELK1, NFAT, E12, PAX4) but not for HIF-1, which is known to be an important controller of genomic responses to hypoxia. Conclusions The multiple analytical approaches used in this study suggests that the acute response to 30 min of transient hypoxia in the late-gestation fetus results in reduced cellular metabolism and a pattern of gene expression that is consistent with cellular oxygen and ATP starvation. In this early time point, we see a vigorous gene response. But, like the hypothalamus, the transcriptomic response is not consistent with mediation by HIF-1. If HIF-1 is a significant controller of gene expression in the fetal pituitary after hypoxia, it must be at a later time.

Rickettsiae are associated with a diverse range of invertebrate hosts. Of these, mosquitoes could emerge as one of the most important vectors because of their ability to transmit significant numbers of pathogens and parasites throughout the world. Recent studies have implicated Anopheles gambiae as a potential vector of Rickettsia felis. Herein we report that a metagenome sequencing study identified rickettsial sequence reads in culicine mosquitoes from the Republic of Korea. The detected rickettsiae were characterized by a genus-specific quantitative real-time PCR assay and sequencing of rrs, gltA, 17kDa, ompB, and sca4 genes. Three novel rickettsial genotypes were detected (Rickettsia sp. A12.2646, Rickettsia sp. A12.2638 and Rickettsia sp. A12.3271), from Mansonia uniformis, Culex pipiens, and Aedes esoensis, respectively. The results underscore the need to determine the Rickettsia species diversity associated with mosquitoes, their evolution, distribution and pathogenic potential.

BACKGROUND AND AIMS: The methylene-tetrahydrofolate dehydrogenase (NADP+ dependent) 1-like (MTHFD1L) gene is involved in mitochondrial tetrahydrofolate metabolism. Polymorphisms in MTHFD1L, including rs6922269, have been implicated in risk for coronary artery disease (CAD). We investigated the association between rs6922269 and known metabolic risk factors and survival in two independent cohorts of coronary heart disease patients. METHODS AND RESULTS: DNA and plasma from 1940 patients with acute coronary syndromes were collected a median of 32 days after index hospital admission (Coronary Disease Cohort Study, CDCS). Samples from a validation cohort of 842 patients post-myocardial infarction (PMI) were taken 24-96 hours after hospitalization. DNA samples were genotyped for rs6922269, using a TaqMan assay. Homocysteine and active vitamin B12 were measured by immunoassay in baseline CDCS plasma samples, but not PMI plasma. All cause mortality was documented over follow-up of 4.1 (CDCS) and 8.8 (PMI) years, respectively. rs6922269 genotype frequencies were AA n = 135, 7.0%; GA n = 785, 40.5% and GG n = 1020, 52.5% in the CDCS and similar in the PMI cohort. CDCS patients with AA genotype for rs6922269 had lower levels of co-variate adjusted baseline plasma active vitamin B12 (p = 0.017) and poorer survival than patients with GG or GA genotype (mortality: AA 19.6%, GA 12.0%, GG 11.6%; p = 0.007). In multivariate analysis, rs6922269 genotype predicted survival, independent of established covariate predictors (p = 0.03). However the association between genotype and survival was not validated in the PMI cohort. CONCLUSION: MTHFD1L rs6922269 genotype is associated with active vitamin B12 levels at baseline and may be a marker of prognostic risk in patients with established coronary heart disease.

Genome-wide association (GWA) studies have been limited by the reliance on common variants present on microarrays or imputable from the HapMap Project data. More recently, the completion of the 1000 Genomes Project has provided variant and haplotype information for several million variants derived from sequencing over 1,000 individuals. To help understand the extent to which more variants (including low frequency (1% ≤ MAF

Soybean cyst nematode (SCN; Heterodera glycines Ichinohe) reproduces on the roots of common bean (Phaseolus vulgaris L.) and can cause reductions in plant growth and seed yield. The molecular changes in common bean roots caused by SCN infection are unknown. Identification of genetic factors associated with SCN resistance could help in development of improved bean varieties with high SCN resistance. Gene expression profiling was conducted on common bean roots infected by SCN HG type 0 using next generation RNA sequencing technology. Two pinto bean genotypes, PI533561 and GTS-900, resistant and susceptible to SCN infection, respectively, were used as RNA sources eight days post inoculation. Total reads generated ranged between ~ 3.2 and 5.7 million per library and were mapped to the common bean reference genome. Approximately 70-90% of filtered RNA-seq reads uniquely mapped to the reference genome. In the inoculated roots of resistant genotype PI533561, a total of 353 genes were differentially expressed with 154 up-regulated genes and 199 down-regulated genes when compared to the transcriptome of non- inoculated roots. On the other hand, 990 genes were differentially expressed in SCN-inoculated roots of susceptible genotype GTS-900 with 406 up-regulated and 584 down-regulated genes when compared to non-inoculated roots. Genes encoding nucleotide-binding site leucine-rich repeat resistance (NLR) proteins, WRKY transcription factors, pathogenesis-related (PR) proteins and heat shock proteins involved in diverse biological processes were differentially expressed in both resistant and susceptible genotypes. Overall, suppression of the photosystem was observed in both the responses. Furthermore, RNA-seq results were validated through quantitative real time PCR. This is the first report describing genes/transcripts involved in SCN-common bean interaction and the results will have important implications for further characterization of SCN resistance genes in common bean.

Many genomes have been sequenced to high-quality draft status using Sanger capillary electrophoresis and/or newer short-read sequence data and whole genome assembly techniques. However, even the best draft genomes contain gaps and other imperfections due to limitations in the input data and the techniques used to build draft assemblies. Sequencing biases, repetitive genomic features, genomic polymorphism, and other complicating factors all come together to make some regions difficult or impossible to assemble. Traditionally, draft genomes were upgraded to “phase 3 finished” status using time-consuming and expensive Sanger-based manual finishing processes. For more facile assembly and automated finishing of draft genomes, we present here an automated approach to finishing using long-reads from the Pacific Biosciences RS (PacBio) platform. Our algorithm and associated software tool, PBJelly, (publicly available at https://sourceforge.net/projects/pb-jelly/) automates the finishing process using long sequence reads in a reference-guided assembly process. PBJelly also provides “lift-over” co-ordinate tables to easily port existing annotations to the upgraded assembly. Using PBJelly and long PacBio reads, we upgraded the draft genome sequences of a simulated Drosophila melanogaster, the version 2 draft Drosophila pseudoobscura, an assembly of the Assemblathon 2.0 budgerigar dataset, and a preliminary assembly of the Sooty mangabey. With 24× mapped coverage of PacBio long-reads, we addressed 99% of gaps and were able to close 69% and improve 12% of all gaps in D. pseudoobscura. With 4× mapped coverage of PacBio long-reads we saw reads address 63% of gaps in our budgerigar assembly, of which 32% were closed and 63% improved. With 6.8× mapped coverage of mangabey PacBio long-reads we addressed 97% of gaps and closed 66% of addressed gaps and improved 19%. The accuracy of gap closure was validated by comparison to Sanger sequencing on gaps from the original D. pseudoobscura draft assembly and shown to be dependent on initial reference quality.

The Hedgehog (Hh) signaling pathway directs a multitude of cellular responses during embryogenesis and adult tissue homeostasis. Stimulation of the pathway results in activation of Hh target genes by the transcription factor Ci/Gli, which binds to specific motifs in genomic enhancers. In Drosophila, only a few enhancers (patched, decapentaplegic, wingless, stripe, knot, hairy, orthodenticle) have been shown by in vivo functional assays to depend on direct Ci/Gli regulation. All but one (orthodenticle) contain more than one Ci/Gli site, prompting us to directly test whether homotypic clustering of Ci/Gli binding sites is sufficient to define a Hh-regulated enhancer. We therefore developed a computational algorithm to identify Ci/Gli clusters that are enriched over random expectation, within a given region of the genome. Candidate genomic regions containing Ci/Gli clusters were functionally tested in chicken neural tube electroporation assays and in transgenic flies. Of the 22 Ci/Gli clusters tested, seven novel enhancers (and the previously known patched enhancer) were identified as Hh-responsive and Ci/Gli-dependent in one or both of these assays, including: Cuticular protein 100A (Cpr100A); invected (inv), which encodes an engrailed-related transcription factor expressed at the anterior/posterior wing disc boundary; roadkill (rdx), the fly homolog of vertebrate Spop; the segment polarity gene gooseberry (gsb); and two previously untested regions of the Hh receptor-encoding patched (ptc) gene. We conclude that homotypic Ci/Gli clustering is not sufficient information to ensure Hh-responsiveness; however, it can provide a clue for enhancer recognition within putative Hedgehog target gene loci.

Malaria remains a significant global health burden. The development of an effective malaria vaccine remains as a major challenge with the potential to significantly reduce morbidity and mortality. While Plasmodium spp. have been shown to contain a large number of intrinsically disordered proteins (IDPs) or disordered protein regions, the relationship of protein structure to subcellular localisation and adaptive immune responses remains unclear. In this study, we employed several computational prediction algorithms to identify IDPs at the proteome level of six Plasmodium spp. and to investigate the potential impact of protein disorder on adaptive immunity against P. falciparum parasites. IDPs were shown to be particularly enriched within nuclear proteins, apical proteins, exported proteins and proteins localised to the parasitophorous vacuole. Furthermore, several leading vaccine candidates, and proteins with known roles in host-cell invasion, have extensive regions of disorder. Presentation of peptides by MHC molecules plays an important role in adaptive immune responses, and we show that IDP regions are predicted to contain relatively few MHC class I and II binding peptides owing to inherent differences in amino acid composition compared to structured domains. In contrast, linear B-cell epitopes were predicted to be enriched in IDPs. Tandem repeat regions and non-synonymous single nucleotide polymorphisms were found to be strongly associated with regions of disorder. In summary, immune responses against IDPs appear to have characteristics distinct from those against structured protein domains, with increased antibody recognition of linear epitopes but some constraints for MHC presentation and issues of polymorphisms. These findings have major implications for vaccine design, and understanding immunity to malaria.

FasL expression on human immune cells and cancer cells plays important roles in immune homeostasis and in cancer development. Our previous study suggests that polymorphisms in the FasL promoter can significantly affect the gene expression in human cells. In addition to the functional FasL SNP -844C>T (rs763110), three other SNPs (SNP -756A>G or rs2021837, SNP -478A>T or rs41309790, and SNP -205 C>G or rs74124371) exist in the proximal FasL promoter. In the current study, we established three major FasL hyplotypes in humans. Interestingly, a transcription motif search revealed that the FasL promoter possessed two consensus T-cell factor (TCF/LEF1) binding elements (TBEs), which is either polymorphic (SNP -205C>G) or close to the functional SNP -844C>T. Subsequently, we demonstrate that both FasL TBEs formed complexes with the TCF-4 and β-catenin transcription factors in vitro and in vivo. Co-transfection of LEF-1 and β-catenin transcription factors significantly increased FasL promoter activities, suggesting that FasL is a target gene of the β-catenin/T-cell factor pathway. More importantly, we found that the rare allele (-205G) of the polymorphic FasL TBE (SNP -205C>G) failed to bind the TCF-4 transcription factor and that SNP -205 C>G significantly affected the promoter activity. Furthermore, promoter reporter assays revealed that FasL SNP haplotypes influenced promoter activities in human colon cancer cells and in human T cells. Finally, β-catenin knockdown significantly decreased the FasL expression in human SW480 colon cancer cells. Collectively, our data suggest that β-catenin may be involved in FasL gene regulation and that FasL expression is influenced by FasL SNP haplotypes, which may have significant implications in immune response and tumorigenesis.

Glucocorticoids (GCs) are key mediators of stress response and are widely used as pharmacological agents to treat immune diseases, such as asthma and inflammatory bowel disease, and certain types of cancer. GCs act mainly by activating the GC receptor (GR), which interacts with other transcription factors to regulate gene expression. Here, we combined different functional genomics approaches to gain molecular insights into the mechanisms of action of GC. By profiling the transcriptional response to GC over time in 4 Yoruba (YRI) and 4 Tuscans (TSI) lymphoblastoid cell lines (LCLs), we suggest that the transcriptional response to GC is variable not only in time, but also in direction (positive or negative) depending on the presence of specific interacting transcription factors. Accordingly, when we performed ChIP-seq for GR and NF-κB in two YRI LCLs treated with GC or with vehicle control, we observed that features of GR binding sites differ for up- and down-regulated genes. Finally, we show that eQTLs that affect expression patterns only in the presence of GC are 1.9-fold more likely to occur in GR binding sites, compared to eQTLs that affect expression only in its absence. Our results indicate that genetic variation at GR and interacting transcription factors binding sites influences variability in gene expression, and attest to the power of combining different functional genomic approaches.

The investigation of associations between rare genetic variants and diseases or phenotypes has two goals. Firstly, the identification of which genes or genomic regions are associated, and secondly, discrimination of associated variants from background noise within each region. Over the last few years, many new methods have been developed which associate genomic regions with phenotypes. However, classical methods for high-dimensional data have received little attention. Here we investigate whether several classical statistical methods for high-dimensional data: ridge regression (RR), principal components regression (PCR), partial least squares regression (PLS), a sparse version of PLS (SPLS), and the LASSO are able to detect associations with rare genetic variants. These approaches have been extensively used in statistics to identify the true associations in data sets containing many predictor variables. Using genetic variants identified in three genes that were Sanger sequenced in 1998 individuals, we simulated continuous phenotypes under several different models, and we show that these feature selection and feature extraction methods can substantially outperform several popular methods for rare variant analysis. Furthermore, these approaches can identify which variants are contributing most to the model fit, and therefore both goals of rare variant analysis can be achieved simultaneously with the use of regression regularization methods. These methods are briefly illustrated with an analysis of adiponectin levels and variants in the ADIPOQ gene.

International audience; Streptococcus pyogenes, is an important human pathogen classified within the pyogenic group of streptococci, exclusively adapted to the human host. Our goal was to employ a comparative evolutionary approach to better understand the genomic events concomitant with S. pyogenes human adaptation. As part of ascertaining these events, we sequenced the genome of one of the potential sister species, the agricultural pathogen S. canis, and combined it in a comparative genomics reconciliation analysis with two other closely related species, Streptococcus dysgalactiae and Streptococcus equi, to determine the genes that were gained and lost during S. pyogenes evolution. Genome wide phylogenetic analyses involving 15 Streptococcus species provided convincing support for a clade of S. equi, S. pyogenes, S. dysgalactiae, and S. canis and suggested that the most likely S. pyogenes sister species was S. dysgalactiae. The reconciliation analysis identified 113 genes that were gained on the lineage leading to S. pyogenes. Almost half (46%) of these gained genes were phage associated and 14 showed significant matches to experimentally verified bacteria virulence factors. Subsequent to the origin of S. pyogenes, over half of the phage associated genes were involved in 90 different LGT events, mostly involving different strains of S. pyogenes, but with a high proportion involving the horse specific pathogen S. equi subsp. equi, with the directionality almost exclusively (86%) in the S. pyogenes to S. equi direction. Streptococcus agalactiae appears to have played an important role in the evolution of S. pyogenes with a high proportion of LGTs originating from this species. Overall the analysis suggests that S. pyogenes adaptation to the human host was achieved in part by (i) the integration of new virulence factors (e.g. speB, and the sal locus) and (ii) the construction of new regulation networks (e.g. rgg, and to some extent speB).