Large Set Of Genes Research Articles

Evolutionary Genomics of a Subdivided Species.

The ways in which genetic variation is distributed within and among populations is a key determinant of the evolutionary features of a species. However, most comprehensive studies of these features have been restricted to studies of subdivision in settings known to have been driven by local adaptation, leaving our understanding of the natural dispersion of allelic variation less than ideal. Here, we present a geographic population-genomic analysis of 10 populations of the freshwater microcrustacean Daphnia pulex, an emerging model system in evolutionary genomics. These populations exhibit a pattern of moderate isolation-by-distance, with an average migration rate of 0.6 individuals per generation, and average effective population sizes of ∼650,000 individuals. Most populations contain numerous private alleles, and genomic scans highlight the presence of islands of excessively high population subdivision for more common alleles. A large fraction of such islands of population divergence likely reflect historical neutral changes, including rare stochastic migration and hybridization events. The data do point to local adaptive divergence, although the precise nature of the relevant variation is diffuse and cannot be associated with particular loci, despite the very large sample sizes involved in this study. In contrast, an analysis of between-species divergence highlights positive selection operating on a large set of genes with functions nearly nonoverlapping with those involved in local adaptation, in particular ribosome structure, mitochondrial bioenergetics, light reception and response, detoxification, and gene regulation. These results set the stage for using D. pulex as a model for understanding the relationship between molecular and cellular evolution in the context of natural environments.

A transcriptional regulator that boosts grain yields and shortens the growth duration of rice.

Complex biological processes such as plant growth and development are often under the control of transcription factors that regulate the expression of large sets of genes and activate subordinate transcription factors in a cascade-like fashion. Here, by screening candidate photosynthesis-related transcription factors in rice, we identified a DREB (Dehydration Responsive Element Binding) family member, OsDREB1C, in which expression is induced by both light and low nitrogen status. We show that OsDREB1C drives functionally diverse transcriptional programs determining photosynthetic capacity, nitrogen utilization, and flowering time. Field trials with OsDREB1C-overexpressing rice revealed yield increases of 41.3 to 68.3% and, in addition, shortened growth duration, improved nitrogen use efficiency, and promoted efficient resource allocation, thus providing a strategy toward achieving much-needed increases in agricultural productivity.

Spruce giga-genomes: structurally similar yet distinctive with differentially expanding gene families and rapidly evolving genes.

Spruces (Picea spp.) are coniferous trees widespread in boreal and mountainous forests of the northern hemisphere, with large economic significance and enormous contributions to global carbon sequestration. Spruces harbor very large genomes with high repetitiveness, hampering their comparative analysis. Here, we present and compare the genomes of four different North American spruces: the genome assemblies for Engelmann spruce (Picea engelmannii) and Sitka spruce (Picea sitchensis) together with improved and more contiguous genome assemblies for white spruce (Picea glauca) and for a naturally occurring introgress of these three species known as interior spruce (P. engelmannii × glauca × sitchensis). The genomes were structurally similar, and a large part of scaffolds could be anchored to a genetic map. The composition of the interior spruce genome indicated asymmetric contributions from the three ancestral genomes. Phylogenetic analysis of the nuclear and organelle genomes revealed a topology indicative of ancient reticulation. Different patterns of expansion of gene families among genomes were observed and related with presumed diversifying ecological adaptations. We identified rapidly evolving genes that harbored high rates of non-synonymous polymorphisms relative to synonymous ones, indicative of positive selection and its hitchhiking effects. These gene sets were mostly distinct between the genomes of ecologically contrasted species, and signatures of convergent balancing selection were detected. Stress and stimulus response was identified as the most frequent function assigned to expanding gene families and rapidly evolving genes. These two aspects of genomic evolution were complementary in their contribution to divergent evolution of presumed adaptive nature. These more contiguous spruce giga-genome sequences should strengthen our understanding of conifer genome structure and evolution, as their comparison offers clues into the genetic basis of adaptation and ecology of conifers at the genomic level. They will also provide tools to better monitor natural genetic diversity and improve the management of conifer forests. The genomes of four closely related North American spruces indicate that their high similarity at the morphological level is paralleled by the high conservation of their physical genome structure. Yet, the evidence of divergent evolution is apparent in their rapidly evolving genomes, supported by differential expansion of key gene families and large sets of genes under positive selection, largely in relation to stimulus and environmental stress response.

P-436 EndoTime: Non-categorical timing estimates for luteal endometrium

Abstract Study question Can expression measurements of a small panel of genes be used to develop a continuous, non-categorical model for the improvement of endometrial biopsy timing accuracy? Summary answer Measuring expression levels of six genes (IL2RB, IGFBP1, CXCL14, DPP4, GPX3, and SLC15A2) is sufficient to obtain and assess substantially more accurate timing estimates. What is known already Commercially available endometrial timing approaches based on gene expression require much larger gene sets and use a categorical approach that classifies samples as pre-receptive, receptive, or post-receptive. Study design, size, duration Gene expression was measured by RT-qPCR in different sample sets, comprising a total of 664 endometrial biopsies obtained 4 to 12 days after a self-reported positive home ovulation test. A further 36 endometrial samples were profiled by RT-qPCR as well as RNA-sequencing. Participants/materials, setting, methods A computational procedure, named ‘EndoTime’, was established that models the temporal profile of each gene and estimates the timing of each sample. Iterating these steps, temporal profiles are gradually refined as sample timings are being updated, and confidence in timing estimates is increased. After convergence, the method reports updated timing estimates for each sample while preserving the overall distribution of time points. Main results and the role of chance The Wilcoxon rank-sum test was used to confirm that ordering samples by EndoTime estimates yields sharper temporal expression profiles for held-out genes (not used when determining sample timings) than ordering the same expression values by patient-reported times (GPX3 : p < 0.005; CXCL14 : p < 2.7e-6; DPP4 : p < 3.7e-13). Pearson correlation between EndoTime estimates for the same sample set but based on RT-qPCR or RNA-sequencing data showed high degree of congruency between the two (p = 8.6e-10, R-squared = 0.687). Estimated timings did not differ significantly between control subjects and patients with recurrent pregnancy loss or recurrent implantation failure (p > 0.05). Limitations, reasons for caution Timing estimates are informed by glandular gene expression and will only represent the temporal state of other endometrial cell types if in synchrony with the epithelium. Additionally, methods that estimate the day of ovulation are still required as these data are essential inputs in our method. Wider implications of the findings Our method is the first to assay the temporal state of endometrial samples on a continuous domain, enabling accurate temporal profiling of any gene in luteal phase biopsies for a wide range of research applications and, potentially, clinical use. It is freely available, open-source software including supporting data sets. Trial registration number Not applicable

Sequence and structural conservation reveal fingerprint residues in TRP channels.

Transient receptor potential (TRP) proteins are a large family of cation-selective channels, surpassed in variety only by voltage-gated potassium channels. Detailed molecular mechanisms governing how membrane voltage, ligand binding, or temperature can induce conformational changes promoting the open state in TRP channels are still a matter of debate. Aiming to unveil distinctive structural features common to the transmembrane domains within the TRP family, we performed phylogenetic reconstruction, sequence statistics, and structural analysis over a large set of TRP channel genes. Here, we report an exceptionally conserved set of residues. This fingerprint is composed of twelve residues localized at equivalent three-dimensional positions in TRP channels from the different subtypes. Moreover, these amino acids are arranged in three groups, connected by a set of aromatics located at the core of the transmembrane structure. We hypothesize that differences in the connectivity between these different groups of residues harbor the apparent differences in coupling strategies used by TRP subgroups.

Ethylene Promotes Expression of the Appressorium- and Pathogenicity-Related Genes via GPCR- and MAPK-Dependent Manners in Colletotrichum gloeosporioides.

Ethylene (ET) represents a signal that can be sensed by plant pathogenic fungi to accelerate their spore germination and subsequent infection. However, the molecular mechanisms of responses to ET in fungi remain largely unclear. In this study, Colletotrichum gloeosporioides was investigated via transcriptomic analysis to reveal the genes that account for the ET-regulated fungal development and virulence. The results showed that ET promoted genes encoding for fungal melanin biosynthesis enzymes, extracellular hydrolases, and appressorium-associated structure proteins at 4 h after treatment. When the germination lasted until 24 h, ET induced multiple appressoria from every single spore, but downregulated most of the genes. Loss of selected ET responsive genes encoding for scytalone dehydratase (CgSCD1) and cerato-platanin virulence protein (CgCP1) were unable to alter ET sensitivity of C. gloeosporioides in vitro but attenuated the influence of ET on pathogenicity. Knockout of the G-protein-coupled receptors CgGPCR3-1/2 and the MAPK signaling pathway components CgMK1 and CgSte11 resulted in reduced ET sensitivity. Taken together, this study in C. gloeosporioides reports that ET can cause transcription changes in a large set of genes, which are mainly responsible for appressorium development and virulence expression, and these processes are dependent on the GPCR and MAPK pathways.

Cytokine and Chemokine Receptor Patterns of Human Malignant Melanoma Cell Lines.

Cytokine and chemokine receptors can promote tumor progression, invasion, and metastasis development by inducing different intracellular signaling pathways. The aim of this study was to determine the cytokine and chemokine receptor gene expression patterns in human melanoma cell lines. We found a large set of cytokine and chemokine receptor genes that were significantly differentially expressed between melanoma cell lines that originated from different subtypes of primary melanomas as well as cell lines that originated from melanoma metastases. The relative expressions of two receptor genes (CCR2 and TNFRSF11B) were positively correlated with the invasive potential of the cell lines, whereas a negative correlation was observed for the TNFRSF14 gene expression. We also found a small set of receptor genes that exhibited a significantly decreased expression in association with a BRAFV600E mutation. Based on our results, we assume that the analyzed cytokine and chemokine receptor collection may provide potential to distinguish the different subtypes of melanomas, helping us to understand the biological behavior of BRAFV600E-mutated melanoma cells.

Drought attenuates plant defence against bacterial pathogens by suppressing the expression of CBP60g/SARD1 during combined stress.

In nature, plants are frequently exposed to drought and bacterial pathogens simultaneously. However, information on how the drought and defence pathways interact and orchestrate global transcriptional regulation is limited. Here, we show that moderate drought stress enhances the susceptibility of Arabidopsis thaliana to Pseudomonas syringae pv. tomato DC3000. Using transcriptome meta-analysis, we found that drought and bacterial stress antagonistically modulate a large set of genes predominantly involved in salicylic acid (SA) and abscisic acid (ABA) signalling networks. We identified that the levels of SA and ABA are dynamically regulated during the course of stress. Importantly, under combined stress, drought through the ABA pathway downregulates the induction of Calmodulin-binding Protein 60 g (CBP60g) and Systemic Acquired Resistance Deficient 1 (SARD1), two transcription factors crucial for SA production upon bacterial infection. We also identified an important role of NPR1-LIKE PROTEIN 3 and 4 (NPR3/4) transcriptional repressors in the drought-mediated negative regulation of CBP60g/SARD1 expression. Using a genetic approach, we show that CBP60g/SARD1 expression is the key determinant of plant defence against bacterial pathogens under combined stress. Thus, these transcription factors act as critical nodes for the crosstalk between drought and bacterial stress signalling under combined stress in plants.

Genetic diversity of meat quality related genes in Argentinean pigs

Genetic influence on pork quality exists between breeds and within a breed. The variation is caused by a large set of genes, and pork quality traits have a multifactorial background. Research into the genetics of meat quality found causative mutations associated with marked effects on pig meat value. This study aimed to investigate the segregation of meat quality-related SNPs and compare their diversity and genetics in commercial and Creole pigs from different farms in the North-West of Argentina. A screen for SNPs in RYR1, PRKAG3, CAST, and SOX6 candidate genes and the differentiation of their genotypes by PCR–RFLP was conducted. All genes were characterized by a high level of polymorphism and heterozygosity, and populations showed no differences in the genetic structure for the analyzed SNPs. These results highlighted the role of pig genotypes as a source of basic variability potentially affecting processed meat products and fresh meat.

EndoTime: non-categorical timing estimates for luteal endometrium.

STUDY QUESTIONCan the accuracy of timing of luteal phase endometrial biopsies based on urinary ovulation testing be improved by measuring the expression of a small number of genes and a continuous, non-categorical modelling approach?SUMMARY ANSWERMeasuring the expression levels of six genes (IL2RB, IGFBP1, CXCL14, DPP4, GPX3 and SLC15A2) is sufficient to obtain substantially more accurate timing estimates and to assess the reliability of timing estimates for each sample.WHAT IS KNOWN ALREADYCommercially available endometrial timing approaches based on gene expression require large gene sets and use a categorical approach that classifies samples as pre-receptive, receptive or post-receptive.STUDY DESIGN, SIZE, DURATIONGene expression was measured by RTq-PCR in different sample sets, comprising a total of 664 endometrial biopsies obtained 4–12 days after a self-reported positive home ovulation test. A further 36 endometrial samples were profiled by RTq-PCR as well as RNA-sequencing.PARTICIPANTS/MATERIALS, SETTING, METHODSA computational procedure, named ‘EndoTime’, was established that models the temporal profile of each gene and estimates the timing of each sample. Iterating these steps, temporal profiles are gradually refined as sample timings are being updated, and confidence in timing estimates is increased. After convergence, the method reports updated timing estimates for each sample while preserving the overall distribution of time points.MAIN RESULTS AND THE ROLE OF CHANCEThe Wilcoxon rank-sum test was used to confirm that ordering samples by EndoTime estimates yields sharper temporal expression profiles for held-out genes (not used when determining sample timings) than ordering the same expression values by patient-reported times (GPX3: P < 0.005; CXCL14: P < 2.7e−6; DPP4: P < 3.7e−13). Pearson correlation between EndoTime estimates for the same sample set but based on RTq-PCR or RNA-sequencing data showed a high degree of congruency between the two (P = 8.6e−10, R2 = 0.687). Estimated timings did not differ significantly between control subjects and patients with recurrent pregnancy loss or recurrent implantation failure (P > 0.05).LARGE SCALE DATAThe RTq-PCR data files are available via the GitHub repository for the EndoTime software at https://github.com/AE-Mitchell/EndoTime, as is the code used for pre-processing of RTq-PCR data. The RNA-sequencing data are available on GEO (accession GSE180485).LIMITATIONS, REASONS FOR CAUTIONTiming estimates are informed by glandular gene expression and will only represent the temporal state of other endometrial cell types if in synchrony with the epithelium. Methods that estimate the day of ovulation are still required as these data are essential inputs in our method. Our approach, in its current iteration, performs batch correction such that larger sample batches impart greater accuracy to timing estimations. In theory, our method requires endometrial samples obtained at different days in the luteal phase. In practice, however, this is not a concern as timings based on urinary ovulation testing are associated with a sufficient level of noise to ensure that a variety of time points will be sampled.WIDER IMPLICATIONS OF THE FINDINGSOur method is the first to assay the temporal state of luteal-phase endometrial samples on a continuous domain. It is freely available with fully shared data and open-source software. EndoTime enables accurate temporal profiling of any gene in luteal endometrial samples for a wide range of research applications and, potentially, clinical use.STUDY FUNDING/COMPETING INTEREST(S)This study was supported by a Wellcome Trust Investigator Award (Grant/Award Number: 212233/Z/18/Z) and the Tommy's National Miscarriage Research Centre. None of the authors have any competing interests. J.L. was funded by the Biotechnology and Biological Sciences Research Council (UK) through the Midlands Integrative Biology Training Partnership (MIBTP, BB/M01116X/1).

The Complete Genome of Chelonus insularis Reveals Dynamic Arrangement of Genome Components in Parasitoid Wasps That Produce Bracoviruses.

Bracoviruses (BVs) are endogenized nudiviruses in parasitoid wasps of the microgastroid complex (family Braconidae). Microgastroid wasps have coopted nudivirus genes to produce replication-defective virions that females use to transfer virulence genes to parasitized hosts. The microgastroid complex further consists of six subfamilies and ∼50,000 species but current understanding of BV gene inventories and organization primarily derives from analysis of two wasp species in the subfamily Microgastrinae (Microplitis demolitor and Cotesia congregata) that produce M. demolitor BV (MdBV) and C. congregata BV (CcBV). Notably, several genomic features of MdBV and CcBV remain conserved since divergence of M. demolitor and C. congregata ∼53 million years ago (MYA). However, it is unknown whether these conserved traits more broadly reflect BV evolution, because no complete genomes exist for any microgastroid wasps outside the Microgastrinae. In this regard, the subfamily Cheloninae is of greatest interest because it diverged earliest from the Microgastrinae (∼85 MYA) after endogenization of the nudivirus ancestor. Here, we present the complete genome of Chelonus insularis, which is an egg-larval parasitoid in the Cheloninae that produces C. insularis BV (CinsBV). We report that the inventory of nudivirus genes in C. insularis is conserved but are dissimilarly organized compared to M. demolitor and C. congregata. Reciprocally, CinsBV proviral segments share organizational features with MdBV and CcBV but virulence gene inventories exhibit almost no overlap. Altogether, our results point to the functional importance of a conserved inventory of nudivirus genes and a dynamic set of virulence genes for the successful parasitism of hosts. Our results also suggest organizational features previously identified in MdBV and CcBV are likely not essential for BV virion formation. IMPORTANCE Bracoviruses are a remarkable example of virus endogenization, because large sets of genes from a nudivirus ancestor continue to produce virions that thousands of wasp species rely upon to parasitize hosts. Understanding how these genes interact and have been coopted by wasps for novel functions is of broad interest in the study of virus evolution. This work characterizes bracovirus genome components in the parasitoid wasp Chelonus insularis, which together with existing wasp genomes captures a large portion of the diversity among wasp species that produce bracoviruses. Results provide new information about how bracovirus genome components are organized in different wasps while also providing additional insights on key features required for function.

Uncontrolled CD21low age-associated and B1 B cell accumulation caused by failure of an EGR2/3 tolerance checkpoint.

CD21low age-associated or atypical memory B cells are autoantibody enriched and poised for plasma cell differentiation. These cells overaccumulate in chronic infections, autoimmune disease, and immunodeficiency, posing the question of what checkpoints normally oppose their accumulation. Here, we reveal a critical role for paralogous calcium-NFAT-regulated transcription factors EGR2 and EGR3 that are induced in self-reactive B cells. CD21low and B1 B cells lacking EGR2 and EGR3 accumulate and circulate in young mice in numbers 10- to 20-fold greater than normal and overexpress a large set of EGR2 ChIP-seq target genes, including known drivers of plasma cell differentiation. Most follicular B cells constitutively express Egr2 proportionally to surface IgM downregulation by self-antigens, and EGR2/3 deficiency abolishes this cardinal feature of B cell anergy. These results explain the cardinal features of B cell anergy, define a key transcriptional checkpoint repressing CD21low B cell formation, and inform how NFATC1 or EGR2 mutations promote B1 cell-derived chronic lymphocytic leukemias.

Genomic and Metabolomic Analyses of the Marine Fungus Emericellopsis cladophorae: Insights into Saltwater Adaptability Mechanisms and Its Biosynthetic Potential.

The genus Emericellopsis is found in terrestrial, but mainly in marine, environments with a worldwide distribution. Although Emericellopsis has been recognized as an important source of bioactive compounds, the range of metabolites expressed by the species of this genus, as well as the genes involved in their production are still poorly known. Untargeted metabolomics, using UPLC- QToF–MS/MS, and genome sequencing (Illumina HiSeq) was performed to unlock E. cladophorae MUM 19.33 chemical diversity. The genome of E. cladophorae is 26.9 Mb and encodes 8572 genes. A large set of genes encoding carbohydrate-active enzymes (CAZymes), secreted proteins, transporters, and secondary metabolite biosynthetic gene clusters were identified. Our analysis also revealed genomic signatures that may reflect a certain fungal adaptability to the marine environment, such as genes encoding for (1) the high-osmolarity glycerol pathway; (2) osmolytes’ biosynthetic processes; (3) ion transport systems, and (4) CAZymes classes allowing the utilization of marine polysaccharides. The fungal crude extract library constructed revealed a promising source of antifungal (e.g., 9,12,13-Trihydroxyoctadec-10-enoic acid, hymeglusin), antibacterial (e.g., NovobiocinA), anticancer (e.g., daunomycinone, isoreserpin, flavopiridol), and anti-inflammatory (e.g., 2’-O-Galloylhyperin) metabolites. We also detected unknown compounds with no structural match in the databases used. The metabolites’ profiles of E. cladophorae MUM 19.33 fermentations were salt dependent. The results of this study contribute to unravel aspects of the biology and ecology of this marine fungus. The genome and metabolome data are relevant for future biotechnological exploitation of the species.

Exercise Preconditioning Blunts Early Atrogenes Expression and Atrophy in Gastrocnemius Muscle of Hindlimb Unloaded Mice.

A large set of FoxOs-dependent genes play a primary role in controlling muscle mass during hindlimb unloading. Mitochondrial dysfunction can modulate such a process. We hypothesized that endurance exercise before disuse can protect against disuse-induced muscle atrophy by enhancing peroxisome proliferator-activated receptor-γ coactivator-1α (PGC1α) expression and preventing mitochondrial dysfunction and energy-sensing AMP-activated protein kinase (AMPK) activation. We studied cross sectional area (CSA) of muscle fibers of gastrocnemius muscle by histochemistry following 1, 3, 7, and 14 days of hindlimb unloading (HU). We used Western blotting and qRT-PCR to study mitochondrial dynamics and FoxOs-dependent atrogenes’ expression at 1 and 3 days after HU. Preconditioned animals were submitted to moderate treadmill exercise for 7 days before disuse. Exercise preconditioning protected the gastrocnemius from disuse atrophy until 7 days of HU. It blunted alterations in mitochondrial dynamics up to 3 days after HU and the expression of most atrogenes at 1 day after disuse. In preconditioned mice, the activation of atrogenes resumed 3 days after HU when mitochondrial dynamics, assessed by profusion and pro-fission markers (mitofusin 1, MFN1, mitofusin 2, MFN2, optic atrophy 1, OPA1, dynamin related protein 1, DRP1 and fission 1, FIS1), PGC1α levels, and AMPK activation were at a basal level. Therefore, the normalization of mitochondrial dynamics and function was not sufficient to prevent atrogenes activation just a few days after HU. The time course of sirtuin 1 (SIRT1) expression and content paralleled the time course of atrogenes’ expression. In conclusion, seven days of endurance exercise counteracted alterations of mitochondrial dynamics and the activation of atrogenes early into disuse. Despite the normalization of mitochondrial dynamics, the effect on atrogenes’ suppression died away within 3 days of HU. Interestingly, muscle protection lasted until 7 days of HU. A longer or more intense exercise preconditioning may prolong atrogenes suppression and muscle protection.

Changes in antennal gene expression underlying sensory system maturation in Rhodnius prolixus

Triatomine bugs are the blood feeding insect vectors transmitting Chagas disease to humans, a neglected tropical disease that affects over 8 million people, mainly in Latin America. The behavioral responses to host cues and bug signals in Rhodnius prolixus are state dependent, i.e., they vary as a function of post-ecdysis age. At the molecular level, these changes in behavior are probably due to a modulation of peripheral and central processes. In the present study, we report a significant modulation of the expression of a large set of sensory-related genes. Results were generated by means of antennal transcriptomes of 5th instar larvae along the first week (days 0, 2, 4, 6 and 8) after ecdysis sequenced using the Illumina HiSeq platform. Significant age-induced changes in transcript abundance were established for more than 6120 genes (54,7% of 11,186 genes expressed) in the antenna of R. prolixus. This was especially true between the first two days after ecdysis when more than 2500 genes had their expression significantly altered. In contrast, expression profiles were almost identical between day 6 and 8, with only a few genes showing significant modulation of their expression. A total of 86 sensory receptors, odorant carriers and odorant degrading enzymes were significantly modulated across age points and clustered into three distinct expression profiles. The set of sensory genes whose expression increased with age (profile 3) may include candidates underlying the increased responsiveness to host cues shown by R. prolixus during the first days after molting. For the first time, we describe the maturation process undergone at the molecular level by the peripheral sensory system of a hemimetabolous insect.

StartLink and StartLink+: Prediction of Gene Starts in Prokaryotic Genomes.

State-of-the-art algorithms of ab initio gene prediction for prokaryotic genomes were shown to be sufficiently accurate. A pair of algorithms would agree on predictions of gene 3′ends. Nonetheless, predictions of gene starts would not match for 15–25% of genes in a genome. This discrepancy is a serious issue that is difficult to be resolved due to the absence of sufficiently large sets of genes with experimentally verified starts. We have introduced StartLink that infers gene starts from conservation patterns revealed by multiple alignments of homologous nucleotide sequences. We also have introduced StartLink+ combining both ab initio and alignment-based methods. The ability of StartLink to predict the start of a given gene is restricted by the availability of homologs in a database. We observed that StartLink made predictions for 85% of genes per genome on average. The StartLink+ accuracy was shown to be 98–99% on the sets of genes with experimentally verified starts. In comparison with database annotations, we observed that the annotated gene starts deviated from the StartLink+ predictions for ∼5% of genes in AT-rich genomes and for 10–15% of genes in GC-rich genomes on average. The use of StartLink+ has a potential to significantly improve gene start annotation in genomic databases.

Transcriptomic Response to Perkinsus marinus in Two Crassostrea Oysters Reveals Evolutionary Dynamics of Host-Parasite Interactions.

Infectious disease outbreaks are causing widespread declines of marine invertebrates including corals, sea stars, shrimps, and molluscs. Dermo is a lethal infectious disease of the eastern oyster Crassostrea virginica caused by the protist Perkinsus marinus. The Pacific oyster Crassostrea gigas is resistant to Dermo due to differences in the host-parasite interaction that is not well understood. We compared transcriptomic responses to P. marinus challenge in the two oysters at early and late infection stages. Dynamic and orchestrated regulation of large sets of innate immune response genes were observed in both species with remarkably similar patterns for most orthologs, although responses in C. virginica were stronger, suggesting strong or over-reacting immune response could be a cause of host mortality. Between the two species, several key immune response gene families differed in their expansion, sequence variation and/or transcriptional response to P. marinus, reflecting evolutionary divergence in host-parasite interaction. Of note, significant upregulation of inhibitors of apoptosis (IAPs) was observed in resistant C. gigas but not in susceptible C. virginica, suggesting upregulation of IAPs is an active defense mechanism, not a passive response orchestrated by P. marinus. Compared with C. gigas, C. virginica exhibited greater expansion of toll-like receptors (TLRs) and positive selection in P. marinus responsive TLRs. The C1q domain containing proteins (C1qDCs) with the galactose-binding lectin domain that is involved in P. marinus recognition, were only present and significantly upregulated in C. virginica. These results point to previously undescribed differences in host defense genes between the two oyster species that may account for the difference in susceptibility, providing an expanded portrait of the evolutionary dynamics of host-parasite interaction in lophotrochozoans that lack adaptive immunity. Our findings suggest that C. virginica and P. marinus have a history of coevolution and the recent outbreaks may be due to increased virulence of the parasite.

Identification of downstream effectors of retinoic acid specifying the zebrafish pancreas by integrative genomics

Retinoic acid (RA) is a key signal for the specification of the pancreas. Still, the gene regulatory cascade triggered by RA in the endoderm remains poorly characterized. In this study, we investigated this regulatory network in zebrafish by combining RNA-seq, RAR ChIP-seq and ATAC-seq assays. By analysing the effect of RA and of the RA receptor (RAR) inverse-agonist BMS493 on the transcriptome and on the chromatin accessibility of endodermal cells, we identified a large set of genes and regulatory regions regulated by RA signalling. RAR ChIP-seq further defined the direct RAR target genes in zebrafish, including hox genes as well as several pancreatic regulators like mnx1, insm1b, hnf1ba and gata6. Comparison of zebrafish and murine RAR ChIP-seq data highlighted the conserved direct target genes and revealed that some RAR sites are under strong evolutionary constraints. Among them, a novel highly conserved RAR-induced enhancer was identified downstream of the HoxB locus and driving expression in the nervous system and in the gut in a RA-dependent manner. Finally, ATAC-seq data unveiled the role of the RAR-direct targets Hnf1ba and Gata6 in opening chromatin at many regulatory loci upon RA treatment.

High Genomic Identity between Clinical and Environmental Strains of Herbaspirillum frisingense Suggests Pre-Adaptation to Different Hosts and Intrinsic Resistance to Multiple Drugs.

The genus Herbaspirillum is widely studied for its ability to associate with grasses and to perform biological nitrogen fixation. However, the bacteria of the Herbaspirillum genus have frequently been isolated from clinical samples. Understanding the genomic characteristics that allow these bacteria to switch environments and become able to colonize human hosts is essential for monitoring emerging pathogens and predicting outbreaks. In this work, we describe the sequencing, assembly, and annotation of the genome of H. frisingense AU14559 isolated from the sputum of patients with cystic fibrosis, and its comparison with the genomes of the uropathogenic strain VT-16–41 and the environmental strains GSF30, BH-1, IAC152, and SG826. The genes responsible for biological nitrogen fixation were absent from all strains except for GSF30. On the other hand, genes encoding virulence and host interaction factors were mostly shared with environmental strains. We also identified a large set of intrinsic antibiotic resistance genes that were shared across all strains. Unlike other strains, in addition to unique genomic islands, AU14559 has a mutation that renders the biosynthesis of rhamnose and its incorporation into the exopolysaccharide unfeasible. These data suggest that H. frisingense has characteristics that provide it with the metabolic diversity needed to infect and colonize human hosts.

Transcriptome Analysis of Insulin Signaling-Associated Transcription Factors in C. elegans Reveal Their Genome-Wide Target Genes Specificity and Complexity.

Insulin/IGF-1-like signaling (IIS) plays a crucial, conserved role in development, growth, reproduction, stress tolerance, and longevity. In Caenorhabditis elegans, the enhanced longevity under reduced insulin signaling (rIIS) is primarily regulated by the transcription factors (TFs) DAF-16/FOXO, SKN-1/Nrf-1, and HSF1/HSF-1. The specific and coordinated regulation of gene expression by these TFs under rIIS has not been comprehensively elucidated. Here, using RNA-sequencing analysis, we report a systematic study of the complexity of TF-dependent target gene interactions during rIIS under analogous genetic and experimental conditions. We found that DAF-16 regulates only a fraction of the C. elegans transcriptome but controls a large set of genes under rIIS; SKN-1 and HSF-1 show the opposite trend. Both of the latter TFs function as activators and repressors to a similar extent, while DAF-16 is predominantly an activator. For expression of the genes commonly regulated by TFs under rIIS conditions, DAF-16 is the principal determining factor, dominating over the other two TFs, irrespective of whether they activate or repress these genes. The functional annotations and regulatory networks presented in this study provide novel insights into the complexity of the gene regulatory networks downstream of the IIS pathway that controls diverse phenotypes, including longevity.