Search

SARS-CoV-2 infection causing the novel coronavirus disease 2019 (COVID–19) has been responsible for more than 2.8 million deaths and nearly 125 million infections worldwide as of March 2021. In March 2020, the World Health Organization determined that the COVID–19 outbreak is a global pandemic. The urgency and magnitude of this pandemic demanded immediate action and coordination between local, regional, national, and international actors. In that mission, researchers require access to high-quality biological materials and data from SARS-CoV-2 infected and uninfected patients, covering the spectrum of disease manifestations. The “Biobanque québécoise de la COVID-19” (BQC19) is a pan–provincial initiative undertaken in Québec, Canada to enable the collection, storage and sharing of samples and data related to the COVID-19 crisis. As a disease-oriented biobank based on high-quality biosamples and clinical data of hospitalized and non-hospitalized SARS-CoV-2 PCR positive and negative individuals. The BQC19 follows a legal and ethical management framework approved by local health authorities. The biosamples include plasma, serum, peripheral blood mononuclear cells and DNA and RNA isolated from whole blood. In addition to the clinical variables, BQC19 will provide in-depth analytical data derived from the biosamples including whole genome and transcriptome sequencing, proteome and metabolome analyses, multiplex measurements of key circulating markers as well as anti-SARS-CoV-2 antibody responses. BQC19 will provide the scientific and medical communities access to data and samples to better understand, manage and ultimately limit, the impact of COVID-19. In this paper we present BQC19, describe the process according to which it is governed and organized, and address opportunities for future research collaborations. BQC19 aims to be a part of a global communal effort addressing the challenges of COVID–19.

Hybridization capture with in-solution oligonucleotide probes has quickly become the preferred method for enriching specific DNA loci from degraded or ancient samples prior to high-throughput sequencing (HTS). Several companies synthesize sets of probes for in-solution hybridization capture, but these commercial reagents are usually expensive. Methods for economical in-house probe synthesis have been described, but they do not directly address one of the major advantages of commercially synthesised probes: that probe sequences matching many species can be synthesised in parallel and pooled. The ability to make "phylogenetically diverse" probes increases the cost-effectiveness of commercial probe sets, as they can be used across multiple projects (or for projects involving multiple species). However, it is labour-intensive to replicate this with in-house methods, as template molecules must first be generated for each species of interest. While it has been observed that probes can be used to enrich for phylogenetically distant targets, the ability of this effect to compensate for the lack of phylogenetically diverse probes in in-house synthesised probe sets has not been tested. In this study, we present a refined protocol for in-house RNA probe synthesis and evaluated the ability of probes generated using this method from a single species to successfully enrich for the target locus in phylogenetically distant species. We demonstrated that probes synthesized using long-range PCR products from a placental mammal mitochondrion (Bison spp.) could be used to enrich for mitochondrial DNA in birds and marsupials (but not plants). Importantly, our results were obtained for approximately a third of the cost of similar commercially available reagents.

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.

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 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.

Background There is little evidence on the accuracy of psychosis relapse prediction models. Our objective was to undertake a systematic review of relapse prediction models in psychosis. Method We conducted a literature search including studies that developed and/or validated psychosis relapse prediction models, with or without external model validation. Models had to target people with psychosis and predict relapse. The key databases searched were; Embase, Medline, Medline In-Process Citations & Daily Update, PsychINFO, BIOSIS Citation Index, CINAHL, and Science Citation Index, from inception to September 2016. Prediction modelling studies were assessed for risk of bias and applicability using the PROBAST tool. Results There were two eligible studies, which included 33,088 participants. One developed a model using prodromal symptoms and illness-related variables, which explained 14% of relapse variance but was at high risk of bias. The second developed a model using administrative data which was moderately discriminative (C = 0.631) and associated with relapse (OR 1.11 95% CI 1.10, 1.12) and achieved moderately discriminative capacity when validated (C = 0.630). The risk of bias was low. Conclusions Due to a lack of high quality evidence it is not possible to make any specific recommendations about the predictors that should be included in a prognostic model for relapse. For instance, it is unclear whether prodromal symptoms are useful for predicting relapse. The use of routine data to develop prediction models may be a more promising approach, although we could not empirically compare the two included studies.

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.

Campylobacter jejuni is a major cause of bacterial gastroenteritis worldwide, and the capsular polysaccharide (CPS) of this organism is required for persistence and disease. C. jejuni produces over 47 different capsular structures, including a unique O-methyl phosphoramidate (MeOPN) modification present on most C. jejuni isolates. Although the MeOPN structure is rare in nature it has structural similarity to some synthetic pesticides. In this study, we have demonstrated, by whole genome comparisons and high resolution magic angle spinning NMR, that MeOPN modifications are common to several Campylobacter species. Using MeOPN biosynthesis and transferase mutants generated in C. jejuni strain 81-176, we observed that loss of MeOPN from the cell surface correlated with increased invasion of Caco-2 epithelial cells and reduced resistance to killing by human serum. In C. jejuni, the observed serum mediated killing was determined to result primarily from activation of the classical complement pathway. The C. jejuni MeOPN transferase mutant showed similar levels of colonization relative to the wild-type in chickens, but showed a five-fold drop in colonization when co-infected with the wild-type in piglets. In Galleria mellonella waxmoth larvae, the MeOPN transferase mutant was able to kill the insects at wild-type levels. Furthermore, injection of the larvae with MeOPN-linked monosaccharides or CPS purified from the wild-type strain did not result in larval killing, indicating that MeOPN does not have inherent insecticidal activity.

Isolating pure microbial cultures and cultivating them in the laboratory on defined media is used to more fully characterize the metabolism and physiology of organisms. However, identifying an appropriate growth medium for a novel isolate remains a challenging task. Even organisms with sequenced and annotated genomes can be difficult to grow, despite our ability to build genome-scale metabolic networks that connect genomic data with metabolic function. The scientific literature is scattered with information about defined growth media used successfully for cultivating a wide variety of organisms, but to date there exists no centralized repository to inform efforts to cultivate less characterized organisms by bridging the gap between genomic data and compound composition for growth media. Here we present MediaDB, a manually curated database of defined media that have been used for cultivating organisms with sequenced genomes, with an emphasis on organisms with metabolic network models. The database is accessible online, can be queried by keyword searches or downloaded in its entirety, and can generate exportable individual media formulation files. The data assembled in MediaDB facilitate comparative studies of organism growth media, serve as a starting point for formulating novel growth media, and contribute to formulating media for in silico investigation of metabolic networks. MediaDB is freely available for public use at https://mediadb.systemsbiology.net.

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

Background Crocodilians exhibit a spectrum of rostral shape from long snouted (longirostrine), through to short snouted (brevirostrine) morphologies. The proportional length of the mandibular symphysis correlates consistently with rostral shape, forming as much as 50% of the mandible’s length in longirostrine forms, but 10% in brevirostrine crocodilians. Here we analyse the structural consequences of an elongate mandibular symphysis in relation to feeding behaviours. Methods/Principal Findings Simple beam and high resolution Finite Element (FE) models of seven species of crocodile were analysed under loads simulating biting, shaking and twisting. Using beam theory, we statistically compared multiple hypotheses of which morphological variables should control the biomechanical response. Brevi- and mesorostrine morphologies were found to consistently outperform longirostrine types when subject to equivalent biting, shaking and twisting loads. The best predictors of performance for biting and twisting loads in FE models were overall length and symphyseal length respectively; for shaking loads symphyseal length and a multivariate measurement of shape (PC1– which is strongly but not exclusively correlated with symphyseal length) were equally good predictors. Linear measurements were better predictors than multivariate measurements of shape in biting and twisting loads. For both biting and shaking loads but not for twisting, simple beam models agree with best performance predictors in FE models. Conclusions/Significance Combining beam and FE modelling allows a priori hypotheses about the importance of morphological traits on biomechanics to be statistically tested. Short mandibular symphyses perform well under loads used for feeding upon large prey, but elongate symphyses incur high strains under equivalent loads, underlining the structural constraints to prey size in the longirostrine morphotype. The biomechanics of the crocodilian mandible are largely consistent with beam theory and can be predicted from simple morphological measurements, suggesting that crocodilians are a useful model for investigating the palaeobiomechanics of other aquatic tetrapods.

UNLABELLEDPyruvate dehydrogenase (PDH) complex (PDC) deficiency is an inborn error of pyruvate metabolism causing a variety of neurologic manifestations. Systematic analyses of development of affected brain structures and the cellular processes responsible for their impairment have not been performed due to the lack of an animal model for PDC deficiency.METHODSIn the present study we investigated a murine model of systemic PDC deficiency by interrupting the X-linked Pdha1 gene encoding the α subunit of PDH to study its role on brain development and behavioral studies.RESULTSMale embryos died prenatally but heterozygous females were born. PDC activity was reduced in the brain and other tissues in female progeny compared to age-matched control females. Immunohistochemical analysis of several brain regions showed that approximately 40% of cells were PDH(-). The oxidation of glucose to CO2 and incorporation of glucose-carbon into fatty acids were reduced in brain slices from 15 day-old PDC-deficient females. Histological analyses showed alterations in several structures in white and gray matters in 35 day-old PDC-deficient females. Reduction in total cell number and reduced dendritic arbors in Purkinje neurons were observed in PDC-deficient females. Furthermore, cell proliferation, migration and differentiation into neurons by newly generated cells were reduced in the affected females during pre- and postnatal periods. PDC-deficient mice had normal locomotor activity in a novel environment but displayed decreased startle responses to loud noises and there was evidence of abnormal pre-pulse inhibition of the startle reflex.CONCLUSIONSThe results show that a reduction in glucose metabolism resulting in deficit in energy production and fatty acid biosynthesis impairs cellular differentiation and brain development in PDC-deficient mice.

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.

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.

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.

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 Comparative Toxicogenomics Database (CTD; http://ctdbase.org/) is a public resource that curates interactions between environmental chemicals and gene products, and their relationships to diseases, as a means of understanding the effects of environmental chemicals on human health. CTD provides a triad of core information in the form of chemical-gene, chemical-disease, and gene-disease interactions that are manually curated from scientific articles. To increase the efficiency, productivity, and data coverage of manual curation, we have leveraged text mining to help rank and prioritize the triaged literature. Here, we describe our text-mining process that computes and assigns each article a document relevancy score (DRS), wherein a high DRS suggests that an article is more likely to be relevant for curation at CTD. We evaluated our process by first text mining a corpus of 14,904 articles triaged for seven heavy metals (cadmium, cobalt, copper, lead, manganese, mercury, and nickel). Based upon initial analysis, a representative subset corpus of 3,583 articles was then selected from the 14,094 articles and sent to five CTD biocurators for review. The resulting curation of these 3,583 articles was analyzed for a variety of parameters, including article relevancy, novel data content, interaction yield rate, mean average precision, and biological and toxicological interpretability. We show that for all measured parameters, the DRS is an effective indicator for scoring and improving the ranking of literature for the curation of chemical-gene-disease information at CTD. Here, we demonstrate how fully incorporating text mining-based DRS scoring into our curation pipeline enhances manual curation by prioritizing more relevant articles, thereby increasing data content, productivity, and efficiency.

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).