FUNC: a package for detecting significant associations between gene sets and ontological annotations

Nkx2.3 transcription factor is a key regulator of mucous cell identity in salivary glands

Saliva is vital to oral health, fulfilling multiple functions in the oral cavity. Three pairs of major salivary glands and hundreds of minor salivary glands contribute to saliva production. The secretory acinar cells within these glands include two distinct populations. Serous acinar cells secrete a watery saliva containing enzymes, while mucous acinar cells secrete a more viscous fluid containing highly glycosylated mucins. Despite their shared developmental origins, the parotid gland (PG) is comprised of only serous acinar cells, while the sublingual gland (SLG) contains predominantly mucous acinar cells. The instructive signals that govern the identity of serous versus mucous acinar cell phenotypes are not yet known. The homeobox transcription factor Nkx2.3 is uniquely expressed in the SLG. Disruption of the Nkx2.3 gene was reported to delay the maturation of SLG mucous acinar cells. To examine whether Nkx2.3 plays a role in directing the mucous cell phenotype, we analyzed SLG from Nkx2.3-/- mice using RNAseq, immunostaining and proteomic analysis of saliva. Our results indicate that Nkx2.3, most likely in concert with other transcription factors uniquely expressed in the SLG, is a key regulator of the molecular program that specifies the identity of mucous acinar cells.

Range of chromatin accessibility configurations are permissive of GABAergic fate acquisition in developing mouse brain

In recent single-cell -omics studies, both the differential activity of transcription factors regulating cell fate determination and differential genome activation have been tested for utility as descriptors of cell types. Naturally, genome accessibility and gene expression are interlinked. To understand the variability in genomic feature activation in the GABAergic neurons of different spatial origins, we have mapped accessible chromatin regions and mRNA expression in single cells derived from the developing mouse central nervous system (CNS). We first defined a reference set of open chromatin regions for scATAC-seq read quantitation across samples, allowing comparison of chromatin accessibility between brain regions and cell types directly. Second, we integrated the scATAC-seq and scRNA-seq data to form a unified resource of transcriptome and chromatin accessibility landscape for the cell types in di- and telencephalon, midbrain and anterior hindbrain of E14.5 mouse embryo. Importantly, we implemented resolution optimization at the clustering, and automatized the cell typing step. We show high level of concordance between the cell clustering based on the chromatin accessibility and the transcriptome in analyzed neuronal lineages, indicating that both genome and transcriptome features can be used for cell type definition. Hierarchical clustering by the similarity in accessible chromatin reveals that the genomic feature activation correlates with neurotransmitter phenotype, selector gene expression, cell differentiation stage and neuromere origins.

Lipid nanoparticle delivery limits antisense oligonucleotide activity and cellular distribution in the brain after intracerebroventricular injection

Antisense oligonucleotide (ASO) therapeutics are being investigated for a broad range of neurological diseases. While ASOs have been effective in the clinic, improving productive ASO internalization into target cells remains a key area of focus in the field. Here, we investigated how the delivery of ASO-loaded lipid nanoparticles (LNPs) affects ASO activity, subcellular trafficking, and distribution in the brain. We show that ASO-LNPs increase ASO activity up to 100-fold in cultured primary brain cells as compared to non-encapsulated ASO. However, in contrast to the widespread ASO uptake and activity observed following free ASO delivery in vivo, LNP-delivered ASOs did not downregulate mRNA levels throughout the brain after intracerebroventricular injection. This lack of activity was likely due to ASO accumulation in cells lining the ventricles and blood vessels. Furthermore, we reveal a formulation-dependent activation of the immune system post dosing, suggesting that LNP encapsulation cannot mask cellular ASO backbone-mediated toxicities. Together, these data provide insights into how LNP encapsulation affects ASO distribution as well as activity in the brain, and a foundation that enables future optimization of brain-targeting ASO-LNPs.

Novel insights into plant defensin ingestion induced metabolic responses in the polyphagous insect pest Helicoverpa armigera

Lepidopteran insect pest Helicoverpa armigera is one of the most destructive pests of crop plants and several biotechnological approaches are being developed for its control. Plant defensins are small cationic and cysteine-rich peptides that play a role in plant defense. Ingestion of a defensin from Capsicum annuum (CanDef-20) induced a dose-dependent reduction in larval and pupal mass, delayed metamorphosis and also severely reduced fecundity and fertility in H. armigera. To understand the molecular mechanisms of CanDef-20 ingestion-mediated antibiosis in H. armigera larvae, a comparative transcriptomics analysis was carried out. Predominant downregulation of GOs represents serine-type endopeptidases, structural constituents of ribosomes and integral membrane components and differential upregulation of ATP binding, nucleus and translation, while up-regulation of nucleic acid binding represented by transposable elements, were detected. Different isoforms of lipase, serine endopeptidase, glutathione S-transferase, cadherin, alkaline phosphatase and aminopeptidases were found to be upregulated as a compensatory response to CanDef-20 ingestion. In vitro enzyme assays and qPCR analysis of some representative genes associated with vital cellular processes like metamorphosis, food digestion and gut membrane indicated adaptive differential regulations in CanDef-20 fed H. armigera larvae. We conclude that CanDef-20 ingestion affects insect metabolism in a number of ways through its interaction with cell membrane, enzymes, cytoplasmic proteins and triggering transposon mobilization which are linked to growth retardation and adaptive strategies in H. armigera.

Mutation-class dependent signatures outweigh disease-associated processes in cystic fibrosis cells

BACKGROUND

The phenotypic heterogeneity observed in Cystic Fibrosis (CF) patients suggests the involvement of other genes, besides CFTR. Here, we combined transcriptome and proteome analysis to understand the global gene expression patterns associated with five prototypical CFTR mutations.

RESULTS

Evaluation of differentially expressed genes and proteins unveiled common and mutation-specific changes revealing functional signatures that are much more associated with the specific molecular defects associated with each mutation than to the CFTR loss-of-function phenotype. The combination of both datasets revealed that mutation-specific detected translated-transcripts (Dtt) have a high level of consistency.

CONCLUSIONS

This is the first combined transcriptomic and proteomic study focusing on prototypical CFTR mutations. Analysis of Dtt provides novel insight into the pathophysiology of CF, and the mechanisms through which each mutation class causes disease and will likely contribute to the identification of new therapeutic targets and/or biomarkers for CF.

Klarigi: Characteristic explanations for semantic biomedical data

Annotation of biomedical entities with ontology classes provides for formal semantic analysis and mobilisation of background knowledge in determining their relationships. To date, enrichment analysis has been routinely employed to identify classes that are over-represented in annotations across sets of groups, such as biosample gene expression profiles or patient phenotypes, and is useful for a range of tasks including differential diagnosis and causative variant prioritisation. These approaches, however, usually consider only univariate relationships, make limited use of the semantic features of ontologies, and provide limited information and evaluation of the explanatory power of both singular and grouped candidate classes. Moreover, they are not designed to solve the problem of deriving cohesive, characteristic, and discriminatory sets of classes for entity groups. We have developed a new tool, called Klarigi, which introduces multiple scoring heuristics for identification of classes that are both compositional and discriminatory for groups of entities annotated with ontology classes. The tool includes a novel algorithm for derivation of multivariable semantic explanations for entity groups, makes use of semantic inference through live use of an ontology reasoner, and includes a classification method for identifying the discriminatory power of candidate sets, in addition to significance testing apposite to traditional enrichment approaches. We describe the design and implementation of Klarigi, including its scoring and explanation determination methods, and evaluate its use in application to two test cases with clinical significance, comparing and contrasting methods and results with literature-based and enrichment analysis methods. We demonstrate that Klarigi produces characteristic and discriminatory explanations for groups of biomedical entities in two settings. We also show that these explanations recapitulate and extend the knowledge held in existing biomedical databases and literature for several diseases. We conclude that Klarigi provides a distinct and valuable perspective on biomedical datasets when compared with traditional enrichment methods, and therefore constitutes a new method by which biomedical datasets can be explored, contributing to improved insight into semantic data.

Integrated Omic Approaches Reveal Molecular Mechanisms of Tolerance during Soybean and Meloidogyne incognita Interactions

The root-knot nematode (RKN), Meloidogyne incognita, is a devastating soybean pathogen worldwide. The use of resistant cultivars is the most effective method to prevent economic losses caused by RKNs. To elucidate the mechanisms involved in resistance to RKN, we determined the proteome and transcriptome profiles from roots of susceptible (BRS133) and highly tolerant (PI 595099) Glycine max genotypes 4, 12, and 30 days after RKN infestation. After in silico analysis, we described major defense molecules and mechanisms considered constitutive responses to nematode infestation, such as mTOR, PI3K-Akt, relaxin, and thermogenesis. The integrated data allowed us to identify protein families and metabolic pathways exclusively regulated in tolerant soybean genotypes. Among them, we highlighted the phenylpropanoid pathway as an early, robust, and systemic defense process capable of controlling M. incognita reproduction. Associated with this metabolic pathway, 29 differentially expressed genes encoding 11 different enzymes were identified, mainly from the flavonoid and derivative pathways. Based on differential expression in transcriptomic and proteomic data, as well as in the expression profile by RT-qPCR, and previous studies, we selected and overexpressed the GmPR10 gene in transgenic tobacco to assess its protective effect against M. incognita. Transgenic plants of the T₂ generation showed up to 58% reduction in the M. incognita reproduction factor. Finally, data suggest that GmPR10 overexpression can be effective against the plant parasitic nematode M. incognita, but its mechanism of action remains unclear. These findings will help develop new engineered soybean genotypes with higher performance in response to RKN infections.

Multi-omic brain and behavioral correlates of cell-free fetal DNA methylation in macaque maternal obesity models

Maternal obesity during pregnancy is associated with neurodevelopmental disorder (NDD) risk. We utilized integrative multi-omics to examine maternal obesity effects on offspring neurodevelopment in rhesus macaques by comparison to lean controls and two interventions. Differentially methylated regions (DMRs) from longitudinal maternal blood-derived cell-free fetal DNA (cffDNA) significantly overlapped with DMRs from infant brain. The DMRs were enriched for neurodevelopmental functions, methylation-sensitive developmental transcription factor motifs, and human NDD DMRs identified from brain and placenta. Brain and cffDNA methylation levels from a large region overlapping mir-663 correlated with maternal obesity, metabolic and immune markers, and infant behavior. A DUX4 hippocampal co-methylation network correlated with maternal obesity, infant behavior, infant hippocampal lipidomic and metabolomic profiles, and maternal blood measurements of DUX4 cffDNA methylation, cytokines, and metabolites. We conclude that in this model, maternal obesity was associated with changes in the infant brain and behavior, and these differences were detectable in pregnancy through integrative analyses of cffDNA methylation with immune and metabolic factors.

Decoding Functional High-Density Lipoprotein Particle Surfaceome Interactions

High-density lipoprotein (HDL) is a mixture of complex particles mediating reverse cholesterol transport (RCT) and several cytoprotective activities. Despite its relevance for human health, many aspects of HDL-mediated lipid trafficking and cellular signaling remain elusive at the molecular level. During HDL’s journey throughout the body, its functions are mediated through interactions with cell surface receptors on different cell types. To characterize and better understand the functional interplay between HDL particles and tissue, we analyzed the surfaceome-residing receptor neighborhoods with which HDL potentially interacts. We applied a combination of chemoproteomic technologies including automated cell surface capturing (auto-CSC) and HATRIC-based ligand–receptor capturing (HATRIC-LRC) on four different cellular model systems mimicking tissues relevant for RCT. The surfaceome analysis of EA.hy926, HEPG2, foam cells, and human aortic endothelial cells (HAECs) revealed the main currently known HDL receptor scavenger receptor B1 (SCRB1), as well as 155 shared cell surface receptors representing potential HDL interaction candidates. Since vascular endothelial growth factor A (VEGF-A) was recently found as a regulatory factor of transendothelial transport of HDL, we next analyzed the VEGF-modulated surfaceome of HAEC using the auto-CSC technology. VEGF-A treatment led to the remodeling of the surfaceome of HAEC cells, including the previously reported higher surfaceome abundance of SCRB1. In total, 165 additional receptors were found on HAEC upon VEGF-A treatment representing SCRB1 co-regulated receptors potentially involved in HDL function. Using the HATRIC-LRC technology on human endothelial cells, we specifically aimed for the identification of other bona fide (co-)receptors of HDL beyond SCRB1. HATRIC-LRC enabled, next to SCRB1, the identification of the receptor tyrosine-protein kinase Mer (MERTK). Through RNA interference, we revealed its contribution to endothelial HDL binding and uptake. Furthermore, subsequent proximity ligation assays (PLAs) demonstrated the spatial vicinity of MERTK and SCRB1 on the endothelial cell surface. The data shown provide direct evidence for a complex and dynamic HDL receptome and that receptor nanoscale organization may influence binding and uptake of HDL.

Human follicular mites: Ectoparasites becoming symbionts

Most humans carry mites in the hair follicles of their skin for their entire lives. Follicular mites are the only metazoans tha continuously live on humans. We propose that Demodex folliculorum (Acari) represents a transitional stage from a host-injuring obligate parasite to an obligate symbiont. Here, we describe the profound impact of this transition on the genome and physiology of the mite. Genome sequencing revealed that the permanent host association of D. folliculorum led to an extensive genome reduction through relaxed selection and genetic drift, resulting in the smallest number of protein-coding genes yet identified among panarthropods. Confocal microscopy revealed that this gene loss coincided with an extreme reduction in the number of cells. Single uninucleate muscle cells are sufficient to operate each of the three segments that form each walking leg. While it has been assumed that the reduction of the cell number in parasites starts early in development, we identified a greater total number of cells in the last developmental stage (nymph) than in the terminal adult stage, suggesting that reduction starts at the adult or ultimate stage of development. This is the first evolutionary step in an arthropod species adopting a reductive, parasitic or endosymbiotic lifestyle. Somatic nuclei show underreplication at the diploid stage. Novel eye structures or photoreceptors as well as a unique human host melatonin-guided day/night rhythm are proposed for the first time. The loss of DNA repair genes coupled with extreme endogamy might have set this mite species on an evolutionary dead-end trajectory.

A Data-Driven Medical Decision Framework for Associating Adverse Drug Events with Drug-Drug Interaction Mechanisms

Adverse drug events (ADEs) occur when multiple drugs interact within an individual, thus causing effects that were not initially predicted. Such toxic interactions lead to morbidity and mortality. Contemporary research surrounding ADEs has tended to focus on the detection of potential ADEs without great concern for elucidating the associations of drug-drug interaction (DDI) mechanisms that can predict potential adverse drug reactions (ADRs). Such associations are of great practical importance for everyday pharmacovigilance efforts. This study presents a data-driven framework for conducting knowledge-driven data analysis that combines a semantic inference system and enrichment analysis in order to identify potential ADE mechanisms. The framework was used to rank mechanisms according to their relevance for DDIs and also to categorize ADEs based on the number of DDI mechanism associations identified through enrichment analysis. Its validity is demonstrated through using both commercial and publicly available DDI resources. The results of this study solidly prove the framework's effectiveness and highlight potential for future research by way of incorporating additional and broader data to deepen and expand its capabilities.

Placenta and fetal brain share a neurodevelopmental disorder DNA methylation profile in a mouse model of prenatal PCB exposure

Polychlorinated biphenyls (PCBs) are developmental neurotoxicants implicated as environmental risk factors for neurodevelopmental disorders (NDDs). Here, we report the effects of prenatal exposure to a human-relevant mixture of PCBs on the DNA methylation profiles of mouse placenta and fetal brain. Thousands of differentially methylated regions (DMRs) distinguish placenta and fetal brain from PCB-exposed mice from sex-matched vehicle controls. In both placenta and fetal brain, PCB-associated DMRs are enriched for functions related to neurodevelopment and cellular signaling and enriched within regions of bivalent chromatin. The placenta and brain PCB DMRs overlap significantly and map to a shared subset of genes enriched for Wnt signaling, Slit/Robo signaling, and genes differentially expressed in NDD models. The consensus PCB DMRs also significantly overlap with DMRs from human NDD brain and placenta. These results demonstrate that PCB-exposed placenta contains a subset of DMRs that overlap fetal brain DMRs relevant to an NDD.

Genomic basis for skin phenotype and cold adaptation in the extinct Steller's sea cow

Steller's sea cow, an extinct sirenian and one of the largest Quaternary mammals, was described by Georg Steller in 1741 and eradicated by humans within 27 years. Here, we complement Steller's descriptions with paleogenomic data from 12 individuals. We identified convergent evolution between Steller's sea cow and cetaceans but not extant sirenians, suggesting a role of several genes in adaptation to cold aquatic (or marine) environments. Among these are inactivations of lipoxygenase genes, which in humans and mouse models cause ichthyosis, a skin disease characterized by a thick, hyperkeratotic epidermis that recapitulates Steller's sea cows' reportedly bark-like skin. We also found that Steller's sea cows' abundance was continuously declining for tens of thousands of years before their description, implying that environmental changes also contributed to their extinction.

Transcriptome Responses of Wild Arachis to UV-C Exposure Reveal Genes Involved in General Plant Defense and Priming

Stress priming is an important strategy for enhancing plant defense capacity to deal with environmental challenges and involves reprogrammed transcriptional responses. Although ultraviolet (UV) light exposure is a widely adopted approach to elicit stress memory and tolerance in plants, the molecular mechanisms underlying UV-mediated plant priming tolerance are not fully understood. Here, we investigated the changes in the global transcriptome profile of wild Arachis stenosperma leaves in response to UV-C exposure. A total of 5751 differentially expressed genes (DEGs) were identified, with the majority associated with cell signaling, protein dynamics, hormonal and transcriptional regulation, and secondary metabolic pathways. The expression profiles of DEGs known as indicators of priming state, such as transcription factors, transcriptional regulators and protein kinases, were further characterized. A meta-analysis, followed by qRT-PCR validation, identified 18 metaDEGs as being commonly regulated in response to UV and other primary stresses. These genes are involved in secondary metabolism, basal immunity, cell wall structure and integrity, and may constitute important players in the general defense processes and establishment of a priming state in A. stenosperma. Our findings contribute to a better understanding of transcriptional dynamics involved in wild Arachis adaptation to stressful conditions of their natural habitats.

RNA sequencing of whole blood reveals early alterations in immune cells and gene expression in Parkinson's disease

Changes in the blood-based RNA transcriptome have the potential to inform biomarkers of Parkinson's disease (PD) progression. Here we sequenced a discovery set of whole-blood RNA species in 4,871 longitudinally collected samples from 1,570 clinically phenotyped individuals from the Parkinson's Progression Marker Initiative (PPMI) cohort. Samples were sequenced to an average of 100 million read pairs to create a high-quality transcriptome. Participants with PD in the PPMI had significantly altered RNA expression (>2,000 differentially expressed genes), including an early and persistent increase in neutrophil gene expression, with a concomitant decrease in lymphocyte cell counts. This was validated in a cohort from the Parkinson's Disease Biomarkers Program (PDBP) consisting of 1,599 participants and by alterations in immune cell subtypes. This publicly available transcriptomic dataset, coupled with available detailed clinical data, provides new insights into PD biological processes impacting whole blood and new paths for developing diagnostic and prognostic PD biomarkers.

An ancestral recombination graph of human, Neanderthal, and Denisovan genomes

Many humans carry genes from Neanderthals, a legacy of past admixture. Existing methods detect this archaic hominin ancestry within human genomes using patterns of linkage disequilibrium or direct comparison to Neanderthal genomes. Each of these methods is limited in sensitivity and scalability. We describe a new ancestral recombination graph inference algorithm that scales to large genome-wide datasets and demonstrate its accuracy on real and simulated data. We then generate a genome-wide ancestral recombination graph including human and archaic hominin genomes. From this, we generate a map within human genomes of archaic ancestry and of genomic regions not shared with archaic hominins either by admixture or incomplete lineage sorting. We find that only 1.5 to 7% of the modern human genome is uniquely human. We also find evidence of multiple bursts of adaptive changes specific to modern humans within the past 600,000 years involving genes related to brain development and function.

Genome Expression Dynamics Reveal the Parasitism Regulatory Landscape of the Root-Knot Nematode Meloidogyne incognita and a Promoter Motif Associated with Effector Genes

Root-knot nematodes (genus Meloidogyne) are the major contributor to crop losses caused by nematodes. These nematodes secrete effector proteins into the plant, derived from two sets of pharyngeal gland cells, to manipulate host physiology and immunity. Successful completion of the life cycle, involving successive molts from egg to adult, covers morphologically and functionally distinct stages and will require precise control of gene expression, including effector genes. The details of how root-knot nematodes regulate transcription remain sparse. Here, we report a life stage-specific transcriptome of Meloidogyne incognita. Combined with an available annotated genome, we explore the spatio-temporal regulation of gene expression. We reveal gene expression clusters and predicted functions that accompany the major developmental transitions. Focusing on effectors, we identify a putative cis-regulatory motif associated with expression in the dorsal glands, providing an insight into effector regulation. We combine the presence of this motif with several other criteria to predict a novel set of putative dorsal gland effectors. Finally, we show this motif, and thereby its utility, is broadly conserved across the Meloidogyne genus, and we name it Mel-DOG. Taken together, we provide the first genome-wide analysis of spatio-temporal gene expression in a root-knot nematode and identify a new set of candidate effector genes that will guide future functional analyses.

Reintroduction of the archaic variant of NOVA1 in cortical organoids alters neurodevelopment

The evolutionarily conserved splicing regulator neuro-oncological ventral antigen 1 (NOVA1) plays a key role in neural development and function. NOVA1 also includes a protein-coding difference between the modern human genome and Neanderthal and Denisovan genomes. To investigate the functional importance of an amino acid change in humans, we reintroduced the archaic allele into human induced pluripotent cells using genome editing and then followed their neural development through cortical organoids. This modification promoted slower development and higher surface complexity in cortical organoids with the archaic version of NOVA1 Moreover, levels of synaptic markers and synaptic protein coassociations correlated with altered electrophysiological properties in organoids expressing the archaic variant. Our results suggest that the human-specific substitution in NOVA1, which is exclusive to modern humans since divergence from Neanderthals, may have had functional consequences for our species' evolution.

Transcriptome Analysis and Knockdown of the Juvenile Hormone Esterase Gene Reveal Abnormal Feeding Behavior in the Sugarcane Giant Borer

The sugarcane giant borer (SGB), Telchin licus licus, is a pest that has strong economic relevance for sugarcane producers. Due to the endophytic behavior of the larva, current methods of management are inefficient. A promising biotechnological management option has been proposed based on RNA interference (RNAi), a process that uses molecules of double-stranded RNA (dsRNA) to specifically knock down essential genes and reduce insect survival. The selection of suitable target genes is often supported by omic sciences. Studies have shown that genes related to feeding adaptation processes are good candidates to be targeted by RNAi for pest management. Among those genes, esterases are highlighted because of their impact on insect development. In this study, the objective was to evaluate the transcriptome responses of the SGB’s gut in order to provide curated data of genes that could be used for pest management by RNAi in future studies. Further, we validated the function of an esterase-coding gene and its potential as a target for RNAi-based control. We sequenced the gut transcriptome of SGB larvae by Illumina HiSeq and evaluated its gene expression profiles in response to different diets (sugarcane stalk and artificial diet). We obtained differentially expressed genes (DEGs) involved in detoxification, digestion, and transport, which suggest a generalist mechanism of adaptation in SGB larvae. Among the DEGs, was identified and characterized a candidate juvenile hormone esterase gene (Tljhe). We knocked down the Tljhe gene by oral delivery of dsRNA molecules and evaluated gene expression in the gut. The survival and nutritional parameters of the larvae were measured along the developmental cycle of treated insects. We found that the gene Tljhe acts as a regulator of feeding behavior. The knockdown of Tljhe triggered a forced starvation state in late larval instars that significantly reduced the fitness of the larvae. However, the mechanism of action of this gene remains unclear, and the correlation between the expression of Tljhe and the levels of juvenile hormone (JH) metabolites in the hemolymph of the SGB must be assessed in future research.

Comparative Genomics Reveals Novel Target Genes towards Specific Control of Plant-Parasitic Nematodes

Plant-parasitic nematodes cause extensive annual yield losses to worldwide agricultural production. Most cultivated plants have no known resistance against nematodes and the few bearing a resistance gene can be overcome by certain species. Chemical methods that have been deployed to control nematodes have largely been banned from use due to their poor specificity and high toxicity. Hence, there is an urgent need for the development of cleaner and more specific control methods. Recent advances in nematode genomics, including in phytoparasitic species, provide an unprecedented opportunity to identify genes and functions specific to these pests. Using phylogenomics, we compared 61 nematode genomes, including 16 for plant-parasitic species and identified more than 24,000 protein families specific to these parasites. In the genome of Meloidogyne incognita, one of the most devastating plant parasites, we found ca. 10,000 proteins with orthologs restricted only to phytoparasitic species and no further homology in protein databases. Among these phytoparasite-specific proteins, ca. 1000 shared the same properties as known secreted effectors involved in essential parasitic functions. Of these, 68 were novel and showed strong expression during the endophytic phase of the nematode life cycle, based on both RNA-seq and RT-qPCR analyses. Besides effector candidates, transcription-related and neuro-perception functions were enriched in phytoparasite-specific proteins, revealing interesting targets for nematode control methods. This phylogenomics analysis constitutes a unique resource for the further understanding of the genetic basis of nematode adaptation to phytoparasitism and for the development of more efficient control methods.

Transcriptome Profiling-Based Analysis of Carbohydrate-Active Enzymes in Aspergillus terreus Involved in Plant Biomass Degradation

Given the global abundance of plant biomass residues, potential exists in biorefinery-based applications with lignocellulolytic fungi. Frequently isolated from agricultural cellulosic materials, Aspergillus terreus is a fungus efficient in secretion of commercial enzymes such as cellulases, xylanases and phytases. In the context of biomass saccharification, lignocellulolytic enzyme secretion was analyzed in a strain of A. terreus following liquid culture with sugarcane bagasse (SB) (1% w/v) and soybean hulls (SH) (1% w/v) as sole carbon source, in comparison to glucose (G) (1% w/v). Analysis of the fungal secretome revealed a maximum of 1.017 UI.mL^–1 xylanases after growth in minimal medium with SB, and 1.019 UI.mL^–1 after incubation with SH as carbon source. The fungal transcriptome was characterized on SB and SH, with gene expression examined in comparison to equivalent growth on G as carbon source. Over 8000 genes were identified, including numerous encoding enzymes and transcription factors involved in the degradation of the plant cell wall, with significant expression modulation according to carbon source. Eighty-nine carbohydrate-active enzyme (CAZyme)-encoding genes were identified following growth on SB, of which 77 were differentially expressed. These comprised 78% glycoside hydrolases, 8% carbohydrate esterases, 2.5% polysaccharide lyases, and 11.5% auxiliary activities. Analysis of the glycoside hydrolase family revealed significant up-regulation for genes encoding 25 different GH family proteins, with predominance for families GH3, 5, 7, 10, and 43. For SH, from a total of 91 CAZyme-encoding genes, 83 were also significantly up-regulated in comparison to G. These comprised 80% glycoside hydrolases, 7% carbohydrate esterases, 5% polysaccharide lyases, 7% auxiliary activities (AA), and 1% glycosyltransferases. Similarly, within the glycoside hydrolases, significant up-regulation was observed for genes encoding 26 different GH family proteins, with predominance again for families GH3, 5, 10, 31, and 43. A. terreus is a promising species for production of enzymes involved in the degradation of plant biomass. Given that this fungus is also able to produce thermophilic enzymes, this first global analysis of the transcriptome following cultivation on lignocellulosic carbon sources offers considerable potential for the application of candidate genes in biorefinery applications.

Low-pass whole genome bisulfite sequencing of neonatal dried blood spots identifies a role for RUNX1 in Down syndrome DNA methylation profiles

Neonatal dried blood spots (NDBS) are a widely banked sample source that enables retrospective investigation into early life molecular events. Here, we performed low-pass whole genome bisulfite sequencing (WGBS) of 86 NDBS DNA to examine early life Down syndrome (DS) DNA methylation profiles. DS represents an example of genetics shaping epigenetics, as multiple array-based studies have demonstrated that trisomy 21 is characterized by genome-wide alterations to DNA methylation. By assaying over 24 million CpG sites, thousands of genome-wide significant (q < 0.05) differentially methylated regions (DMRs) that distinguished DS from typical development and idiopathic developmental delay were identified. Machine learning feature selection refined these DMRs to 22 loci. The DS DMRs mapped to genes involved in neurodevelopment, metabolism, and transcriptional regulation. Based on comparisons with previous DS methylation studies and reference epigenomes, the hypermethylated DS DMRs were significantly (q < 0.05) enriched across tissues while the hypomethylated DS DMRs were significantly (q < 0.05) enriched for blood-specific chromatin states. A ~28 kb block of hypermethylation was observed on chromosome 21 in the RUNX1 locus, which encodes a hematopoietic transcription factor whose binding motif was the most significantly enriched (q < 0.05) overall and specifically within the hypomethylated DMRs. Finally, we also identified DMRs that distinguished DS NDBS based on the presence or absence of congenital heart disease (CHD). Together, these results not only demonstrate the utility of low-pass WGBS on NDBS samples for epigenome-wide association studies, but also provide new insights into the early life mechanisms of epigenomic dysregulation resulting from trisomy 21.

Comparative gut transcriptome analysis of Diatraea saccharalis in response to the dietary source

The sugarcane borer (Diatraea saccharalis, Fabricius, 1794) is a devastating pest that causes millions of dollars of losses each year to sugarcane producers by reducing sugar and ethanol yields. The control of this pest is difficult due to its endophytic behavior and rapid development. Pest management through biotechnological approaches has emerged in recent years as an alternative to currently applied methods. Genetic information about the target pests is often required to perform biotechnology-based management. The genomic and transcriptomic data for D. saccharalis are very limited. Herein, we report a tissue-specific transcriptome of D. saccharalis larvae and a differential expression analysis highlighting the physiological characteristics of this pest in response to two different diets: sugarcane and an artificial diet. Sequencing was performed on the Illumina HiSeq 2000 platform, and a de novo assembly was generated. A total of 27,626 protein-coding unigenes were identified, among which 1,934 sequences were differentially expressed between treatments. Processes such as defence, digestion, detoxification, signaling, and transport were highly represented among the differentially expressed genes (DEGs). Furthermore, seven aminopeptidase genes were identified as candidates to encode receptors of Cry proteins, which are toxins of Bacillus thuringiensis used to control lepidopteran pests. Since plant-insect interactions have produced a considerable number of adaptive responses in hosts and herbivorous insects, the success of phytophagous insects relies on their ability to overcome challenges such as the response to plant defences and the intake of nutrients. In this study, we identified metabolic pathways and specific genes involved in these processes. Thus, our data strongly contribute to the knowledge advancement of insect transcripts, which can be a source of target genes for pest management.

A high-coverage Neandertal genome from Chagyrskaya Cave

We sequenced the genome of a Neandertal from Chagyrskaya Cave in the Altai Mountains, Russia, to 27-fold genomic coverage. We show that this Neandertal was a female and that she was more related to Neandertals in western Eurasia [Prüfer et al., Science 358, 655-658 (2017); Hajdinjak et al., Nature 555, 652-656 (2018)] than to Neandertals who lived earlier in Denisova Cave [Prüfer et al., Nature 505, 43-49 (2014)], which is located about 100 km away. About 12.9% of the Chagyrskaya genome is spanned by homozygous regions that are between 2.5 and 10 centiMorgans (cM) long. This is consistent with the fact that Siberian Neandertals lived in relatively isolated populations of less than 60 individuals. In contrast, a Neandertal from Europe, a Denisovan from the Altai Mountains, and ancient modern humans seem to have lived in populations of larger sizes. The availability of three Neandertal genomes of high quality allows a view of genetic features that were unique to Neandertals and that are likely to have been at high frequency among them. We find that genes highly expressed in the striatum in the basal ganglia of the brain carry more amino-acid-changing substitutions than genes expressed elsewhere in the brain, suggesting that the striatum may have evolved unique functions in Neandertals.

Genetic Architecture of Gene Expression in European and African Americans: An eQTL Mapping Study in GENOA

Most existing expression quantitative trait locus (eQTL) mapping studies have been focused on individuals of European ancestry and are underrepresented in other populations including populations with African ancestry. Lack of large-scale well-powered eQTL mapping studies in populations with African ancestry can both impede the dissemination of eQTL mapping results that would otherwise benefit individuals with African ancestry and hinder the comparable analysis for understanding how gene regulation is shaped through evolution. We fill this critical knowledge gap by performing a large-scale in-depth eQTL mapping study on 1,032 African Americans (AA) and 801 European Americans (EA) in the GENOA cohort. We identified a total of 354,931 eSNPs in AA and 371,309 eSNPs in EA, with 112,316 eSNPs overlapped between the two. We found that eQTL harboring genes (eGenes) are enriched in metabolic pathways and tend to have higher SNP heritability compared to non-eGenes. We found that eGenes that are common in the two populations tend to be less conserved than eGenes that are unique to one population, which are less conserved than non-eGenes. Through conditional analysis, we found that eGenes in AA tend to harbor more independent eQTLs than eGenes in EA, suggesting potentially diverse genetic architecture underlying expression variation in the two populations. Finally, the large sample sizes in GENOA allow us to construct accurate expression prediction models in both AA and EA, facilitating powerful transcriptome-wide association studies. Overall, our results represent an important step toward revealing the genetic architecture underlying expression variation in African Americans.

Over 1000 genetic loci influencing blood pressure with multiple systems and tissues implicated

High blood pressure (BP) remains the major heritable and modifiable risk factor for cardiovascular disease. Persistent high BP, or hypertension, is a complex trait with both genetic and environmental interactions. Despite swift advances in genomics, translating new discoveries to further our understanding of the underlying molecular mechanisms remains a challenge. More than 500 loci implicated in the regulation of BP have been revealed by genome-wide association studies (GWAS) in 2018 alone, taking the total number of BP genetic loci to over 1000. Even with the large number of loci now associated to BP, the genetic variance explained by all loci together remains low (~5.7%). These genetic associations have elucidated mechanisms and pathways regulating BP, highlighting potential new therapeutic and drug repurposing targets. A large proportion of the BP loci were discovered and reported simultaneously by multiple research groups, creating a knowledge gap, where the reported loci to date have not been investigated in a harmonious way. Here, we review the BP-associated genetic variants reported across GWAS studies and investigate their potential impact on the biological systems using in silico enrichment analyses for pathways, tissues, gene ontology and genetic pleiotropy.

Transcriptome and gene expression analysis of three developmental stages of the coffee berry borer, Hypothenemus hampei

Coffee production is a global industry valued at approximately 173 billion US dollars. One of the main challenges facing coffee production is the management of the coffee berry borer (CBB), Hypothenemus hampei, which is considered the primary arthropod pest of coffee worldwide. Current control strategies are inefficient for CBB management. Although biotechnological alternatives, including RNA interference (RNAi), have been proposed in recent years to control insect pests, characterizing the genetics of the target pest is essential for the successful application of these emerging technologies. In this study, we employed RNA-seq to obtain the transcriptome of three developmental stages of the CBB (larva, female and male) to increase our understanding of the CBB life cycle in relation to molecular features. The CBB transcriptome was sequenced using Illumina Hiseq and assembled de novo. Differential gene expression analysis was performed across the developmental stages. The final assembly produced 29,434 unigenes, of which 4,664 transcripts were differentially expressed. Genes linked to crucial physiological functions, such as digestion and detoxification, were determined to be tightly regulated between the reproductive and nonreproductive stages of CBB. The data obtained in this study help to elucidate the critical roles that several genes play as regulatory elements in CBB development.

Genetic Variation in Pan Species Is Shaped by Demographic History and Harbors Lineage-Specific Functions

Chimpanzees (Pan troglodytes) and bonobos (Pan paniscus) are the closest living relatives of humans, but the two species show distinct behavioral and physiological differences, particularly regarding female reproduction. Despite their recent rapid decline, the demographic histories of the two species have been different during the past 1–2 Myr, likely having an impact on their genomic diversity. Here, we analyze the inferred functional consequences of genetic variation across 69 individuals, making use of the most complete data set of genomes in the Pan clade to date. We test to which extent the demographic history influences the efficacy of purifying selection in these species. We find that small historical effective population sizes (N_e) correlate not only with low levels of genetic diversity but also with a larger number of deleterious alleles in homozygosity and an increased proportion of deleterious changes at low frequencies. To investigate the putative genetic basis for phenotypic differences between chimpanzees and bonobos, we exploit the catalog of putatively deleterious protein-coding changes in each lineage. We show that bonobo-specific nonsynonymous changes are enriched in genes related to age at menarche in humans, suggesting that the prominent physiological differences in the female reproductive system between chimpanzees and bonobos might be explained, in part, by putatively adaptive changes on the bonobo lineage.

A catalog of single nucleotide changes distinguishing modern humans from archaic hominins

Throughout the past decade, studying ancient genomes has provided unique insights into human prehistory, and differences between modern humans and other branches like Neanderthals can enrich our understanding of the molecular basis of unique modern human traits. Modern human variation and the interactions between different hominin lineages are now well studied, making it reasonable to go beyond fixed genetic changes and explore changes that are observed at high frequency in present-day humans. Here, we identify 571 genes with non-synonymous changes at high frequency. We suggest that molecular mechanisms in cell division and networks affecting cellular features of neurons were prominently modified by these changes. Complex phenotypes in brain growth trajectory and cognitive traits are likely influenced by these networks and other non-coding changes presented here. We propose that at least some of these changes contributed to uniquely human traits, and should be prioritized for experimental validation.

Genetics of adaptation in modern chicken

We carried out whole genome resequencing of 127 chicken including red jungle fowl and multiple populations of commercial broilers and layers to perform a systematic screening of adaptive changes in modern chicken (Gallus gallus domesticus). We uncovered >21 million high quality SNPs of which 34% are newly detected variants. This panel comprises >115,000 predicted amino-acid altering substitutions as well as 1,100 SNPs predicted to be stop-gain or -loss, several of which reach high frequencies. Signatures of selection were investigated both through analyses of fixation and differentiation to reveal selective sweeps that may have had prominent roles during domestication and breed development. Contrasting wild and domestic chicken we confirmed selection at the BCO2 and TSHR loci and identified 34 putative sweeps co-localized with ALX1, KITLG, EPGR, IGF1, DLK1, JPT2, CRAMP1, and GLI3, among others. Analysis of enrichment between groups of wild vs. commercials and broilers vs. layers revealed a further panel of candidate genes including CORIN, SKIV2L2 implicated in pigmentation and LEPR, MEGF10 and SPEF2, suggestive of production-oriented selection. SNPs with marked allele frequency differences between wild and domestic chicken showed a highly significant deficiency in the proportion of amino-acid altering mutations (P<2.5×10-6). The results contribute to the understanding of major genetic changes that took place during the evolution of modern chickens and in poultry breeding.

Nail abnormalities identified in an ageing study of 30 inbred mouse strains

In a large-scale ageing study, 30 inbred mouse strains were systematically screened for histologic evidence of lesions in all organ systems. Ten strains were diagnosed with similar nail abnormalities. The highest frequency was noted in NON/ShiLtJ mice. Lesions identified fell into two main categories: acute to chronic penetration of the third phalangeal bone through the hyponychium with associated inflammation and bone remodelling or metaplasia of the nail matrix and nail bed associated with severe orthokeratotic hyperkeratosis replacing the nail plate. Penetration of the distal phalanx through the hyponychium appeared to be the initiating feature resulting in nail abnormalities. The accompanying acute to subacute inflammatory response was associated with osteolysis of the distal phalanx. Evaluation of young NON/ShiLtJ mice revealed that these lesions were not often found, or affected only one digit. The only other nail unit abnormality identified was sporadic subungual epidermoid inclusion cysts which closely resembled similar lesions in human patients. These abnormalities, being age-related developments, may have contributed to weight loss due to impacts upon feeding and should be a consideration for future research due to the potential to interact with other experimental factors in ageing studies using the affected strains of mice.

Sexual dimorphisms in the transcriptomes of murine salivary glands

Transcriptional profiling identified 933 sexually dimorphic genes out of the 14 371 protein‐coding genes expressed in the three major murine salivary glands: parotid, sublingual, and submandibular. Most (89%) sex‐specific genes were enriched in a single gland, while only 0.5% of the sexually dimorphic genes were enriched in all glands. The sublingual gland displayed a strong male sex bias (94% of sex‐enriched genes), while a sex preference was not obvious in the parotid or submandibular glands. A subset of transcription factor genes was correlated with the expression of gland‐specific, sex‐enriched genes. Higher expression of Cftr chloride and Scnn1 sodium channels in the male submandibular correlated with greater NaCl reabsorption. In conclusion, adult salivary glands display sex‐ and gland‐specific differences in gene expression that reflect their unique functional properties.

Quantitative evaluation of ontology design patterns for combining pathology and anatomy ontologies

Data are increasingly annotated with multiple ontologies to capture rich information about the features of the subject under investigation. Analysis may be performed over each ontology separately, but recently there has been a move to combine multiple ontologies to provide more powerful analytical possibilities. However, it is often not clear how to combine ontologies or how to assess or evaluate the potential design patterns available. Here we use a large and well-characterized dataset of anatomic pathology descriptions from a major study of aging mice. We show how different design patterns based on the MPATH and MA ontologies provide orthogonal axes of analysis, and perform differently in over-representation and semantic similarity applications. We discuss how such a data-driven approach might be used generally to generate and evaluate ontology design patterns.

Spring Is Coming: Genetic Analyses of the Bud Break Date Locus Reveal Candidate Genes From the Cold Perception Pathway to Dormancy Release in Apple (Malus × domestica Borkh.)

Chilling requirement (CR) for bud dormancy completion determines the time of bud break in apple (Malus × domestica Borkh.). The molecular control of bud dormancy is highly heritable, suggesting a strong genetic control of the trait. An available Infinium II SNP platform for genotyping containing 8,788 single nucleotide polymorphic markers was employed, and linkage maps were constructed in a F₁ cross from the low CR M13/91 and the moderate CR cv. Fred Hough. These maps were used to identify quantitative trait loci (QTL) for bud break date as a trait related to dormancy release. A major QTL for bud break was detected at the beginning of linkage group 9 (LG9). This QTL remained stable during seven seasons in two different growing sites. To increase mapping efficiency in detecting contributing genes underlying this QTL, 182 additional SNP markers located at the locus for bud break were used. Combining linkage mapping and structural characterization of the region, the high proportion of the phenotypic variance in the trait explained by the QTL is related to the coincident positioning of Arabidopsis orthologs for ICE1, FLC, and PRE1 protein-coding genes. The proximity of these genes from the most explanatory markers of this QTL for bud break suggests potential genetic additive effects, reinforcing the hypothesis of inter-dependent mechanisms controlling dormancy induction and release in apple trees.

A genome-wide scan shows evidence for local adaptation in a widespread keystone Neotropical forest tree

The role of natural selection in shaping patterns of diversity is still poorly understood in the Neotropics. We carried out the first genome-wide population genomics study in a Neotropical tree, Handroanthus impetiginosus (Bignoniaceae), sampling 75,838 SNPs by sequence capture in 128 individuals across 13 populations. We found evidences for local adaptation using Bayesian correlations of allele frequency and environmental variables (32 loci in 27 genes) complemented by an analysis of selective sweeps and genetic hitchhiking events using SweepFinder2 (81 loci in 47 genes). Fifteen genes were identified by both approaches. By accounting for population genetic structure, we also found 14 loci with selection signal in a STRUCTURE-defined lineage comprising individuals from five populations, using Outflank. All approaches pinpointed highly diverse and structurally conserved genes affecting plant development and primary metabolic processes. Spatial interpolation forecasted differences in the expected allele frequencies at loci under selection over time, suggesting that H. impetiginosus may track its habitat during climate changes. However, local adaptation through natural selection may also take place, allowing species persistence due to niche evolution. A high genetic differentiation was seen among the H. impetiginosus populations, which, together with the limited power of the experiment, constrains the improved detection of other types of soft selective forces, such as background, balanced, and purifying selection. Small differences in allele frequency distribution among widespread populations and the low number of loci with detectable adaptive sweeps advocate for a polygenic model of adaptation involving a potentially large number of small genome-wide effects.

GOTrapper: a tool to navigate through branches of gene ontology hierarchy

Background

Gene Ontology (GO) is a useful resource of controlled vocabulary that provides information about annotated genes. Based on such resource, finding the biological function is useful for biologists to come up with different hypotheses and help further investigations of an experiment. The biological function for desired genes and gene associations is picked up from a randomly chosen list or through the analysis of differential gene expression. Many tools have been developed to utilize GO knowledge and cluster genes according to relevant biological functions. The retrieved GO terms include both specific and non-specific terms, which is not user-friendly in terms of data analysis. Thus one approach is still missing, which allows navigating through different levels of GO hierarchy manually.

Result

We developed a tool, GOTrapper, which allows moving up or down to the very bottom of the GO hierarchy. This is performed manually by the user, based on an assigned threshold. This tool grabs the shared terms by the desired set of input genes of Homo sapiens. Here, two inputs are possible. “Within” is to find associated terms within one gene list, and “Between” is to find associated terms between two lists. The tool also provides the option to return the terms with the pre-selected evidence codes.

Conclusion

GOTrapper is a user-friendly Java tool that helps the user move up and down the ontology tree, which leads to new hypotheses and devising new association of the input genes. It also allows returning terms of associated genes based on selected evidence codes. This tool can be accessed and is freely available at https://github.com/BioGeneTools/GOTrapper.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2581-8) contains supplementary material, which is available to authorized users.

Immune Gene Diversity in Archaic and Present-day Humans

Genome-wide analyses of two Neandertals and a Denisovan have shown that these archaic humans had lower genetic heterozygosity than present-day people. A similar reduction in genetic diversity of protein-coding genes (gene diversity) was found in exome sequences of three Neandertals. Reduced gene diversity, particularly in genes involved in immunity, may have important functional consequences. In fact, it has been suggested that reduced diversity in immune genes may have contributed to Neandertal extinction. We therefore explored gene diversity in different human groups, and at different time points on the Neandertal lineage, with a particular focus on the diversity of genes involved in innate immunity and genes of the Major Histocompatibility Complex (MHC).

We find that the two Neandertals and a Denisovan have similar gene diversity, all significantly lower than any present-day human. This is true across gene categories, with no gene set showing an excess decrease in diversity compared with the genome-wide average. Innate immune-related genes show a similar reduction in diversity to other genes, both in present-day and archaic humans. There is also no observable decrease in gene diversity over time in Neandertals, suggesting that there may have been no ongoing reduction in gene diversity in later Neandertals, although this needs confirmation with a larger sample size. In both archaic and present-day humans, genes with the highest levels of diversity are enriched for MHC-related functions. In fact, in archaic humans the MHC genes show evidence of having retained more diversity than genes involved only in the innate immune system.

Identification of copper (Cu) stress-responsive grapevine microRNAs and their target genes by high-throughput sequencing

MicroRNAs (miRNAs) are a class of single-stranded non-coding small RNAs (sRNAs) that are 20-24 nucleotides (nt) in length. Extensive studies have indicated that miRNAs play important roles in plant growth, development and stress responses. With more copper (Cu) and copper containing compounds used as bactericides and fungicides in plants, Cu stress has become one of the serious environmental problems that affect plant growth and development. In order to uncover the hidden response mechanisms of Cu stress, two small RNA libraries were constructed from Cu-treated and water-treated (Control) leaves of 'Summer Black' grapevine. Following high-throughput sequencing and filtering, a total of 158 known and 98 putative novel miRNAs were identified in the two libraries. Among these, 100 known and 47 novel miRNAs were identified as differentially expressed under Cu stress. Subsequently, the expression patterns of nine Cu-responsive miRNAs were validated by quantitative real-time PCR (qRT-PCR). There existed some consistency in expression levels of Cu-responsive miRNAs between high throughput sequencing and qRT-PCR assays. The targets prediction of miRNAs indicates that miRNA may regulate some transcription factors, including AP2, SBP, NAC, MYB and ARF during Cu stress. The target genes for two known and two novel miRNAs showed specific cleavage sites at the 10th and/or 11th nucleotide from the 5'-end of the miRNA corresponding to their miRNA complementary sequences. The findings will lay the foundation for exploring the role of the regulation of miRNAs in response to Cu stress and provide valuable gene information for breeding some Cu-tolerant grapevine cultivars.

Analysis of the Trichuris suis excretory/secretory proteins as a function of life cycle stage and their immunomodulatory properties

Parasitic worms have a remarkable ability to modulate host immune responses through several mechanisms including excreted/secreted proteins (ESP), yet the exact nature of these proteins and their targets often remains elusive. Here, we performed mass spectrometry analyses of ESP (TsESP) from larval and adult stages of the pig whipworm Trichuris suis (Ts) and identified ~350 proteins. Transcriptomic analyses revealed large subsets of differentially expressed genes in the various life cycle stages of the parasite. Exposure of bone marrow-derived macrophages and dendritic cells to TsESP markedly diminished secretion of the pro-inflammatory cytokines TNFα and IL-12p70. Conversely, TsESP exposure strongly induced release of the anti-inflammatory cytokine IL-10, and also induced high levels of nitric oxide (NO) and upregulated arginase activity in macrophages. Interestingly, TsESP failed to directly induce CD4+ CD25+ FoxP3+ regulatory T cells (Treg cells), while OVA-pulsed TsESP-treated dendritic cells suppressed antigen-specific OT-II CD4+ T cell proliferation. Fractionation of TsESP identified a subset of proteins that promoted anti-inflammatory functions, an activity that was recapitulated using recombinant T. suis triosephosphate isomerase (TPI) and nucleoside diphosphate kinase (NDK). Our study helps illuminate the intricate balance that is characteristic of parasite-host interactions at the immunological interface, and further establishes the principle that specific parasite-derived proteins can modulate immune cell functions.

Analysis of the Transcriptome in Aspergillus tamarii During Enzymatic Degradation of Sugarcane Bagasse

The production of bioethanol from non-food agricultural residues represents an alternative energy source to fossil fuels for incorporation into the world's economy. Within the context of bioconversion of plant biomass into renewable energy using improved enzymatic cocktails, Illumina RNA-seq transcriptome profiling was conducted on a strain of Aspergillus tamarii, efficient in biomass polysaccharide degradation, in order to identify genes encoding proteins involved in plant biomass saccharification. Enzyme production and gene expression was compared following growth in liquid and semi-solid culture with steam-exploded sugarcane bagasse (SB) (1% w/v) and glucose (1% w/v) employed as contrasting sole carbon sources. Enzyme production following growth in liquid minimum medium supplemented with SB resulted in 0.626 and 0.711 UI.mL^-1 xylanases after 24 and 48 h incubation, respectively. Transcriptome profiling revealed expression of over 7120 genes, with groups of genes modulated according to solid or semi-solid culture, as well as according to carbon source. Gene ontology analysis of genes expressed following SB hydrolysis revealed enrichment in xyloglucan metabolic process and xylan, pectin and glucan catabolic process, indicating up-regulation of genes involved in xylanase secretion. According to carbohydrate-active enzyme (CAZy) classification, 209 CAZyme-encoding genes were identified with significant differential expression on liquid or semi-solid SB, in comparison to equivalent growth on glucose as carbon source. Up-regulated CAZyme-encoding genes related to cellulases (CelA, CelB, CelC, CelD) and hemicellulases (XynG1, XynG2, XynF1, XylA, AxeA, arabinofuranosidase) showed up to a 10-fold log2FoldChange in expression levels. Five genes from the AA9 (GH61) family, related to lytic polysaccharide monooxygenase (LPMO), were also identified with significant expression up-regulation. The transcription factor gene XlnR, involved in induction of hemicellulases, showed up-regulation on liquid and semi-solid SB culture. Similarly, the gene ClrA, responsible for regulation of cellulases, showed increased expression on liquid SB culture. Over 150 potential transporter genes were also identified with increased expression on liquid and semi-solid SB culture. This first comprehensive analysis of the transcriptome of A. tamarii contributes to our understanding of genes and regulatory systems involved in cellulose and hemicellulose degradation in this fungus, offering potential for application in improved enzymatic cocktail development for plant biomass degradation in biorefinery applications.

Comparative root transcriptome of wild Arachis reveals NBS-LRR genes related to nematode resistance

BACKGROUND

The Root-Knot Nematode (RKN), Meloidogyne arenaria, significantly reduces peanut grain quality and yield worldwide. Whilst the cultivated species has low levels of resistance to RKN and other pests and diseases, peanut wild relatives (Arachis spp.) show rich genetic diversity and harbor high levels of resistance to many pathogens and environmental constraints. Comparative transcriptome analysis can be applied to identify candidate resistance genes.

RESULTS

Transcriptome analysis during the early stages of RKN infection of two peanut wild relatives, the highly RKN resistant Arachis stenosperma and the moderately susceptible A. duranensis, revealed genes related to plant immunity with contrasting expression profiles. These included genes involved in hormone signaling and secondary metabolites production and also members of the NBS-LRR class of plant disease resistance (R) genes. From 345 NBS-LRRs identified in A.duranensis reference genome, 52 were differentially expressed between inoculated and control samples, with the majority occurring in physical clusters unevenly distributed on eight chromosomes with preferential tandem duplication. The majority of these NBS-LRR genes showed contrasting expression behaviour between A. duranensis and A. stenosperma, particularly at 6 days after nematode inoculation, coinciding with the onset of the Hypersensitive Response in the resistant species. The physical clustering of some of these NBS-LRR genes correlated with their expression patterns in the contrasting genotypes. Four NBS-LRR genes exclusively expressed in A. stenosperma are located within clusters on chromosome Aradu. A09, which harbors a QTL for RKN resistance, suggesting a functional role for their physical arrangement and their potential involvement in this defense response.

CONCLUSION

The identification of functional novel R genes in wild Arachis species responsible for triggering effective defense cascades can contribute to the crop genetic improvement and enhance peanut resilience to RKN.

Species-Specific Changes in a Primate Transcription Factor Network Provide Insights into the Molecular Evolution of the Primate Prefrontal Cortex

The human prefrontal cortex (PFC) differs from that of other primates with respect to size, histology, and functional abilities. Here, we analyzed genome-wide expression data of humans, chimpanzees, and rhesus macaques to discover evolutionary changes in transcription factor (TF) networks that may underlie these phenotypic differences. We determined the co-expression networks of all TFs with species-specific expression including their potential target genes and interaction partners in the PFC of all three species. Integrating these networks allowed us inferring an ancestral network for all three species. This ancestral network as well as the networks for each species is enriched for genes involved in forebrain development, axonogenesis, and synaptic transmission. Our analysis allows us to directly compare the networks of each species to determine which links have been gained or lost during evolution. Interestingly, we detected that most links were gained on the human lineage, indicating increase TF cooperativity in humans. By comparing network changes between different tissues, we discovered that in brain tissues, but not in the other tissues, the human networks always had the highest connectivity. To pinpoint molecular changes underlying species-specific phenotypes, we analyzed the sub-networks of TFs derived only from genes with species-specific expression changes in the PFC. These sub-networks differed significantly in structure and function between the human and chimpanzee. For example, the human-specific sub-network is enriched for TFs implicated in cognitive disorders and for genes involved in synaptic plasticity and cognitive functions. Our results suggest evolutionary changes in TF networks that might have shaped morphological and functional differences between primate brains, in particular in the human PFC.

Early responses to dehydration in contrasting wild Arachis species

Wild peanut relatives (Arachis spp.) are genetically diverse and were selected throughout evolution to a range of environments constituting, therefore, an important source of allelic diversity for abiotic stress tolerance. In particular, A. duranensis and A. stenosperma, the parents of the reference Arachis A-genome genetic map, show contrasting transpiration behavior under limited water conditions. This study aimed to build a comprehensive gene expression profile of these two wild species under dehydration stress caused by the withdrawal of hydroponic nutrient solution. For this purpose, roots of both genotypes were collected at seven time-points during the early stages of dehydration and used to construct cDNA paired-end libraries. Physiological analyses indicated initial differences in gas exchange parameters between the drought-tolerant genotype of A. duranensis and the drought-sensitive genotype of A. stenosperma. High-quality Illumina reads were mapped against the A. duranensis reference genome and resulted in the identification of 1,235 and 799 Differentially Expressed Genes (DEGs) that responded to the stress treatment in roots of A. duranensis and A. stenosperma, respectively. Further analysis, including functional annotation and identification of biological pathways represented by these DEGs confirmed the distinct gene expression behavior of the two contrasting Arachis species genotypes under dehydration stress. Some species-exclusive and common DEGs were then selected for qRT-PCR analysis, which corroborated the in silico expression profiling. These included genes coding for regulators and effectors involved in drought tolerance responses, such as activation of osmosensing molecular cascades, control of hormone and osmolyte content, and protection of macromolecules. This dataset of transcripts induced during the dehydration process in two wild Arachis genotypes constitute new tools for the understanding of the distinct gene regulation processes in these closely related species but with contrasting drought responsiveness. In addition, our findings provide insights into the nature of drought tolerance in wild germoplasm, which might be explored as novel sources of diversity and useful wild alleles to develop climate-resilient crop varieties.

Transcriptional profiling reveals gland-specific differential expression in the three major salivary glands of the adult mouse

RNA-Seq was used to better understand the molecular nature of the biological differences among the three major exocrine salivary glands in mammals. Transcriptional profiling found that the adult murine parotid, submandibular, and sublingual salivary glands express greater than 14,300 protein-coding genes, and nearly 2,000 of these genes were differentially expressed. Principle component analysis of the differentially expressed genes revealed three distinct clusters according to gland type. The three salivary gland transcriptomes were dominated by a relatively few number of highly expressed genes (6.3%) that accounted for more than 90% of transcriptional output. Of the 912 transcription factors expressed in the major salivary glands, greater than 90% of them were detected in all three glands, while expression for ~2% of them was enriched in an individual gland. Expression of these unique transcription factors correlated with sublingual and parotid specific subsets of both highly expressed and differentially expressed genes. Gene ontology analyses revealed that the highly expressed genes common to all glands were associated with global functions, while many of the genes expressed in a single gland play a major role in the function of that gland. In summary, transcriptional profiling of the three murine major salivary glands identified a limited number of highly expressed genes, differentially expressed genes, and unique transcription factors that represent the transcriptional signatures underlying gland-specific biological properties.

The Transcriptomes of Xiphinema index and Longidorus elongatus Suggest Independent Acquisition of Some Plant Parasitism Genes by Horizontal Gene Transfer in Early-Branching Nematodes

Nematodes have evolved the ability to parasitize plants on at least four independent occasions, with plant parasites present in Clades 1, 2, 10 and 12 of the phylum. In the case of Clades 10 and 12, horizontal gene transfer of plant cell wall degrading enzymes from bacteria and fungi has been implicated in the evolution of plant parasitism. We have used ribonucleic acid sequencing (RNAseq) to generate reference transcriptomes for two economically important nematode species, Xiphinema index and Longidorus elongatus, representative of two genera within the early-branching Clade 2 of the phylum Nematoda. We used a transcriptome-wide analysis to identify putative horizontal gene transfer events. This represents the first in-depth transcriptome analysis from any plant-parasitic nematode of this clade. For each species, we assembled ~30 million Illumina reads into a reference transcriptome. We identified 62 and 104 transcripts, from X. index and L. elongatus, respectively, that were putatively acquired via horizontal gene transfer. By cross-referencing horizontal gene transfer prediction with a phylum-wide analysis of Pfam domains, we identified Clade 2-specific events. Of these, a GH12 cellulase from X. index was analysed phylogenetically and biochemically, revealing a likely bacterial origin and canonical enzymatic function. Horizontal gene transfer was previously shown to be a phenomenon that has contributed to the evolution of plant parasitism among nematodes. Our findings underline the importance and the extensiveness of this phenomenon in the evolution of plant-parasitic life styles in this speciose and widespread animal phylum.

Detecting ancient positive selection in humans using extended lineage sorting

Natural selection that affected modern humans early in their evolution has likely shaped some of the traits that set present-day humans apart from their closest extinct and living relatives. The ability to detect ancient natural selection in the human genome could provide insights into the molecular basis for these human-specific traits. Here, we introduce a method for detecting ancient selective sweeps by scanning for extended genomic regions where our closest extinct relatives, Neandertals and Denisovans, fall outside of the present-day human variation. Regions that are unusually long indicate the presence of lineages that reached fixation in the human population faster than expected under neutral evolution. Using simulations, we show that the method is able to detect ancient events of positive selection and that it can differentiate those from background selection. Applying our method to the 1000 Genomes data set, we find evidence for ancient selective sweeps favoring regulatory changes and present a list of genomic regions that are predicted to underlie positively selected human specific traits.

Using the Neandertal genome to study the evolution of small insertions and deletions in modern humans

Background

Small insertions and deletions occur in humans at a lower rate compared to nucleotide changes, but evolve under more constraint than nucleotide changes. While the evolution of insertions and deletions have been investigated using ape outgroups, the now available genome of a Neandertal can shed light on the evolution of indels in more recent times.

Results

We used the Neandertal genome together with several primate outgroup genomes to differentiate between human insertion/deletion changes that likely occurred before the split from Neandertals and those that likely arose later. Changes that pre-date the split from Neandertals show a smaller proportion of deletions than those that occurred later. The presence of a Neandertal-shared allele in Europeans or Asians but the absence in Africans was used to detect putatively introgressed indels in Europeans and Asians. A larger proportion of these variants reside in intergenic regions compared to other modern human variants, and some variants are linked to SNPs that have been associated with traits in modern humans.

Conclusions

Our results are in agreement with earlier results that suggested that deletions evolve under more constraint than insertions. When considering Neandertal introgressed variants, we find some evidence that negative selection affected these variants more than other variants segregating in modern humans. Among introgressed variants we also identify indels that may influence the phenotype of their carriers. In particular an introgressed deletion associated with a decrease in the time to menarche may constitute an example of a former Neandertal-specific trait contributing to modern human phenotypic diversity.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1018-8) contains supplementary material, which is available to authorized users.

Gene expression analysis in Musa acuminata during compatible interactions with Meloidogyne incognita

BACKGROUND AND AIMS

Endoparasitic root-knot nematodes (RKNs) ( Meloidogyne spp.) cause considerable losses in banana ( Musa spp.), with Meloidogyne incognita a predominant species in Cavendish sub-group bananas. This study investigates the root transcriptome in Musa acuminata genotypes 4297-06 (AA) and Cavendish Grande Naine (CAV; AAA) during early compatible interactions with M. incognita .

METHODS

Roots were analysed by brightfield light microscopy over a 35 d period to examine nematode penetration and morphological cell transformation. RNA samples were extracted 3, 7 and 10 days after inoculation (DAI) with nematode J2 juveniles, and cDNA libraries were sequenced using lllumina HiSeq technology. Sequences were mapped to the M. acuminata ssp. malaccensis var. Pahang genome sequence, differentially expressed genes (DEGs) identified and transcript representation determined by gene set enrichment and pathway mapping.

KEY RESULTS

Microscopic analysis revealed a life cycle of M. incognita completing in 24 d in CAV and 27 d in 4279-06. Comparable numbers of DEGs were up- and downregulated in each genotype, with potential involvement of many in early host defence responses involving reactive oxygen species and jasmonate/ethylene signalling. DEGs revealed concomitant auxin metabolism and cell wall modification processes likely to be involved in giant cell formation. Notable transcripts related to host defence included those coding for leucine-rich repeat receptor-like serine/threonine-protein kinases, peroxidases, thaumatin-like pathogenesis-related proteins, and DREB, ERF, MYB, NAC and WRKY transcription factors. Transcripts related to giant cell development included indole acetic acid-amido synthetase GH3.8 genes, involved in auxin metabolism, as well as genes encoding expansins and hydrolases, involved in cell wall modification.

CONCLUSIONS

Expression analysis in M. acuminata during compatible interactions with RKNs provides insights into genes modulated during infection and giant cell formation. Increased understanding of both defence responses to limit parasitism during compatible interactions and effector-targeted host genes in this complex interaction will facilitate the development of genetic improvement measures for RKNs.

Improved genome assembly of American alligator genome reveals conserved architecture of estrogen signaling

The American alligator, Alligator mississippiensis, like all crocodilians, has temperature-dependent sex determination, in which the sex of an embryo is determined by the incubation temperature of the egg during a critical period of development. The lack of genetic differences between male and female alligators leaves open the question of how the genes responsible for sex determination and differentiation are regulated. Insight into this question comes from the fact that exposing an embryo incubated at male-producing temperature to estrogen causes it to develop ovaries. Because estrogen response elements are known to regulate genes over long distances, a contiguous genome assembly is crucial for predicting and understanding their impact. We present an improved assembly of the American alligator genome, scaffolded with in vitro proximity ligation (Chicago) data. We use this assembly to scaffold two other crocodilian genomes based on synteny. We perform RNA sequencing of tissues from American alligator embryos to find genes that are differentially expressed between embryos incubated at male- versus female-producing temperature. Finally, we use the improved contiguity of our assembly along with the current model of CTCF-mediated chromatin looping to predict regions of the genome likely to contain estrogen-responsive genes. We find that these regions are significantly enriched for genes with female-biased expression in developing gonads after the critical period during which sex is determined by incubation temperature. We thus conclude that estrogen signaling is a major driver of female-biased gene expression in the post-temperature sensitive period gonads.