The simple fool's guide to population genomics via RNA ‐Seq: an introduction to high‐throughput sequencing data analysis

Magnetic orientation in juvenile Atlantic herring (Clupea harengus) could involve cryptochrome 4 as a potential magnetoreceptor

The Earth's magnetic field can provide reliable directional information, allowing migrating animals to orient themselves using a magnetic compass or estimate their position relative to a target using map-based orientation. Here we show for the first time that young, inexperienced herring (Clupea harengus, Ch) have a magnetic compass when they migrate hundreds of kilometres to their feeding grounds. In birds, such as the European robin (Erithacus rubecula), radical pair-based magnetoreception involving cryptochrome 4 (ErCRY4) was demonstrated; the molecular basis of magnetoreception in fish is still elusive. We show that cry4 expression in the eye of herring is upregulated during the migratory season, but not before, indicating a possible use for migration. The amino acid structure of herring ChCRY4 shows four tryptophans and a flavin adenine dinucleotide-binding site, a prerequisite for a magnetic receptor. Using homology modelling, we successfully reconstructed ChCRY4 of herring, DrCRY4 of zebrafish (Danio rerio) and StCRY4 of brown trout (Salmo trutta) and showed that ChCRY4, DrCRY4 and ErCRY4a, but not StCRY4, exhibit very comparable dynamic behaviour. The electron transfer could take place in ChCRY4 in a similar way to ErCRY4a. The combined behavioural, transcriptomic and simulation experiments provide evidence that CRY4 could act as a magnetoreceptor in Atlantic herring.

Complex dynamics of coral gene expression responses to low pH across species

Coral capacity to tolerate low pH affects coral community composition and, ultimately, reef ecosystem function. Low pH submarine discharges ('Ojo'; Yucatán, México) represent a natural laboratory to study plasticity and acclimatization to low pH in relation to ocean acidification. A previous >2-year coral transplant experiment to ambient and low pH common garden sites revealed differential survivorship across species and sites, providing a framework to compare mechanistic responses to differential pH exposures. Here, we examined gene expression responses of transplants of three species of reef-building corals (Porites astreoides, Porites porites and Siderastrea siderea) and their algal endosymbiont communities (Symbiodiniaceae) originating from low pH (Ojo) and ambient pH native origins (Lagoon or Reef). Transplant pH environment had the greatest effect on gene expression of Porites astreoides hosts and symbionts and P. porites hosts. Host P. astreoides Ojo natives transplanted to ambient pH showed a similar gene expression profile to Lagoon natives remaining in ambient pH, providing evidence of plasticity in response to ambient pH conditions. Although origin had a larger effect on host S. siderea gene expression due to differences in symbiont genera within Reef and Lagoon/Ojo natives, subtle effects of low pH on all origins demonstrated acclimatization potential. All corals responded to low pH by differentially expressing genes related to pH regulation, ion transport, calcification, cell adhesion and stress/immune response. This study demonstrates that the magnitude of coral gene expression responses to pH varies considerably among populations, species and holobionts, which could differentially affect acclimatization to and impacts of ocean acidification.

Diversity and distribution of benthic dinoflagellates in Tonga include the potentially harmful genera Gambierdiscus and Fukuyoa

Benthic dinoflagellates that can cause illness, such as ciguatera poisoning (CP), are prevalent around the Pacific but are poorly described in many locations. This study represents the first ecological assessment of benthic harmful algae species in the Kingdom of Tonga, a country where CP occurs regularly. Surveys were conducted in June 2016 in the Tongatapu island group, and in June 2017 across three island groups: Ha'apai, Vava'u, and Tongatapu. Shallow subtidal coastal habitats were investigated by measuring water quality parameters and conducting quadrat surveys. Microalgae samples were collected using either macrophyte collection or the artificial substrate method. Benthic dinoflagellates (Gambierdiscus and/or Fukuyoa, Ostreopsis, and Prorocentrum) were counted using light microscopy, followed by molecular analyses (real-time PCR in 2016 and high throughput sequencing (metabarcoding) in 2017) to identify Gambierdiscus and Fukuyoa to species level. Six species were detected from the Tongatapu island group in 2016 (G. australes, G. carpenteri, G. honu, G. pacificus, F. paulensis, and F. ruetzleri) using real-time PCR. Using the metabarcoding approach in 2017, a total of eight species (G. australes, G. carpenteri, G. honu, G. pacificus, G. cheloniae, G. lewisii, G. polynesiensis, and F. yasumotoi) were detected. Species were detected in mixed assemblages of up to six species, with G. pacificus and G. carpenteri being the most frequently observed. Ha'apai had the highest diversity with eight species detected, which identifies this area as a Gambierdiscus diversity 'hotspot'. Vava'u and Tongatapu had three and six species found respectively. Gambierdiscus polynesiensis, a described ciguatoxin producer and proposed causative agent of CP was found only in Ha'apai and Vava'u in 2017, but not in Tongatapu in either year. Ostreopsis spp. and Prorocentrum spp. were also frequently observed, with Prorocentrum most abundant at the majority of sites. In 2016, the highest number of Gambierdiscus and/or Fukuyoa cells were observed on seagrass (Halodule uninervis) from Sopu, Tongatapu. In 2017, the highest numbers of Gambierdiscus and/or Fukuyoa from artificial substrate samples were recorded in the Halimeda dominant habitat at Neiafu Tahi, Vava'u, a low energy site. This raised the question of the effect of wave motion or currents on abundance measurements from artificial substrates. Differences in detection were noticed between macrophytes and artificial substrates, with higher numbers of species found on artificial substrates. This study provides a baseline of benthic dinoflagellate distributions and diversity for Tonga that may be used for future studies and the development of monitoring programmes.

Simultaneous warming and acidification limit population fitness and reveal phenotype costs for a marine copepod

Phenotypic plasticity and evolutionary adaptation allow populations to cope with global change, but limits and costs to adaptation under multiple stressors are insufficiently understood. We reared a foundational copepod species, Acartia hudsonica, under ambient (AM), ocean warming (OW), ocean acidification (OA), and combined ocean warming and acidification (OWA) conditions for 11 generations (approx. 1 year) and measured population fitness (net reproductive rate) derived from six life-history traits (egg production, hatching success, survival, development time, body size and sex ratio). Copepods under OW and OWA exhibited an initial approximately 40% fitness decline relative to AM, but fully recovered within four generations, consistent with an adaptive response and demonstrating synergy between stressors. At generation 11, however, fitness was approximately 24% lower for OWA compared with the AM lineage, consistent with the cost of producing OWA-adapted phenotypes. Fitness of the OWA lineage was not affected by reversal to AM or low food environments, indicating sustained phenotypic plasticity. These results mimic those of a congener, Acartia tonsa, while additionally suggesting that synergistic effects of simultaneous stressors exert costs that limit fitness recovery but can sustain plasticity. Thus, even when closely related species experience similar stressors, species-specific costs shape their unique adaptive responses.

Differential gene expression analysis in the scallop Argopecten purpuratus exposed to altered pH and temperature conditions in an upwelling-influenced farming area

Increased carbon dioxide in the atmosphere and its absorption across the ocean surface will alter natural variations in pH and temperature levels, occurring in coastal upwelling ecosystems. The scallop Argopecten purpuratus, one of the most economically important species farmed in northern Chile, has been shown to be vulnerable to these environmental drivers. However, the regulatory responses at the gene-level of scallops to these climate stressors remain almost unknown. Consequently, we used an orthogonal experimental design and RNAseq approach to analyze the acute effects of variability in pH and temperature on gene expression in the muscle tissue of A. purpuratus. In respect to control conditions (pH ~ 8.0/ 14 °C), the influence of low pH (~ 7.7) and temperature (14 °C) induced the activation of several genes associated with apoptotic signaling pathways and protein localization to plasma membrane. Elevated temperature (18 °C) and pH (~8.0) conditions increased the expression of transcripts associated with the activation of muscle contraction, regulation, and sarcomere organization effects on muscle tissue. In scallops exposed to low pH and elevated temperature, the genes expressed were differentially associated with the oxidation-reduction process, signal translation, and positive regulation of GTPase activity. These results indicated that the differentially expressed genes under the experimental conditions tested are mainly related to the mitigation of cellular damage and homeostasis control. Our results add knowledge about the function of the adductor muscle in response to stressors in scallops. Furthermore, these results could help in the identification of molecular biomarkers of stress necessary to be integrated into the aquaculture programs for the mitigation of climate change.

Gene expression profile during seed development of Bixa orellana accessions varying in bixin pigment

Diverse morphological, cellular and physiological changes occur during seed maturation in Bixa orellana when the seed tissues form specialized cell glands that produce reddish latex with high bixin amounts. Transcriptomic profiling during seed development in three B. orellana accessions (P12, N4 and N5) with contrasting morphologic characteristics showed enrichment in pathways of triterpenes, sesquiterpenes, and cuticular wax biosynthesis. WGCNA allows groups of all identified genes in six modules the module turquoise, the largest and highly correlated with the bixin content. The high number of genes in this module suggests a diversification of regulatory mechanisms for bixin accumulation with the genes belonging to isoprene, triterpenes and carotene pathways, being more highly correlated with the bixin content. Analysis of key genes of the mevalonate (MVA) and the 2C-methyl-D-erythritol-4-phosphate (MEP) pathways revealed specific activities of orthologs of BoHMGR, BoFFP, BoDXS, and BoHDR. This suggests that isoprenoid production is necessary for compounds included in the reddish latex of developing seeds. The carotenoid-related genes BoPSY2, BoPDS1 and BoZDS displayed a high correlation with bixin production, consistent with the requirement for carotene precursors for apocarotenoid biosynthesis. The BoCCD gene member (BoCCD4-4) and some BoALDH (ALDH2B7.2 and ALDH3I1) and BoMET (BoSABATH1 and BoSABATH8) gene members were highly correlated to bixin in the final seed development stage. This suggested a contributing role for several genes in apocarotenoid production. The results revealed high genetic complexity in the biosynthesis of reddish latex and bixin in specialized seed cell glands in different accessions of B. orellana suggesting gene expression coordination between both metabolite biosynthesis processes.

Geographical distribution, molecular and toxin diversity of the dinoflagellate species Gambierdiscus honu in the Pacific region

An increase in cases of ciguatera poisoning (CP) and expansion of the causative species in the South Pacific region highlight the need for baseline data on toxic microalgal species to help identify new areas of risk and manage known hot spots. Gambierdiscus honu is a toxin producing and potential CP causing dinoflagellate species, first described in 2017. Currently no high-resolution geographical distribution, intraspecific genetic variation or toxin production diversity data is available for G. honu. This research aimed to further characterize G. honu by investigating its distribution using species-specific real-time polymerase chain reaction assays at 25 sites in an area spanning ∼8000 km of the Coral Sea/Pacific Ocean, and assessing intraspecific genetic variation, toxicity and toxin production of isolated strains. Assessment of genetic variation of the partial rRNA operon of isolates demonstrated no significant intraspecific population structure, in addition to a lack of adherence to isolation by distance (IBD) model of evolution. The detected distribution of G. honu in the Pacific region was within the expected tropical to temperate latitudinal ranges of 10° to -30° and extended from Australia to French Polynesia. In the lipophilic fractions, the neuroblastoma cell-based assay (CBA-N2a) showed no ciguatoxin (CTX)-like activity for nine of the 10 isolates, and an atypical pattern for CAWD233 isolate which showed cytotoxic activity in OV- and OV+ conditions. In the same way, liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis confirmed no Pacific-CTXs (CTX-3B, CTX-3C, CTX-4A, CTX-4B) were produced by the ten strains. The CBA-N2a assessment of the hydrophilic fractions showed moderate to high cytotoxicity in both OV- and OV+ condition for all the strains showing a cytotoxic profile similar to that of gambierone. Indeed, this study is the first to show the cytotoxic activity of gambierone on mouse neuroblastoma cells while no cytotoxicity was observed when 44-MG was analysed at the same concentrations using the CBA-N2a. Analysis of the hydrophilic via LC-MS/MS confirmed production of gambierone in all isolates, ranging from 2.1 to 38.1 pg/cell, with 44-methylgambierone (44-MG) also produced by eight of the isolates, ranging from 0.3 to 42.9 pg/cell. No maitotoxin-1 was detected in any of the isolates. Classification of the G. honu strains according to the quantities of gambierone produced aligned with the classification of their cytotoxicity using the CBA-N2a. Finally, no maitotoxin-1 (MTX) was detected in any of the isolates. This study shows G. honu is widely distributed within the Pacific region with no significant intraspecific population structure present. This aligns with the view of microalgal populations as global metapopulations, however more in-depth assessment with other genetic markers could detect further structure. Toxicity diversity across 10 isolates assessed did not display any geographical patterns.

Chemical and molecular tools to probe biological sex differences at multiple length scales

Biological sex differences are observed at multiple different length scales and across organ systems. Gaps in knowledge remain regarding our understanding of how molecular, cellular, and environmental factors contribute to physiological sex differences. Here, we provide our perspective on how chemical and molecular tools can be leveraged to explore sex differences in biology at the molecular, intracellular, extracellular, tissue, and organ length scales. We provide examples where chemical and molecular tools were used to explore sex differences in the cardiovascular, nervous, immune, and reproductive systems. We also provide a future outlook where chemical and molecular tools can be applied to continue investigating sex differences in biology, with the ultimate goal of addressing inequities in biomedical research and approaches to clinical treatments.

In planta transcriptomics reveals conflicts between pattern-triggered immunity and the AlgU sigma factor regulon

In previous work, we determined the transcriptomic impacts of flg22 pre-induced Pattern Triggered Immunity (PTI) in Arabidopsis thaliana on the pathogen Pseudomonas syringae pv. tomato DC3000 (Pto). During PTI exposure we observed expression patterns in Pto reminiscent of those previously observed in a Pto algU mutant. AlgU is a conserved extracytoplasmic function sigma factor which has been observed to regulate over 950 genes in Pto in growth media. We sought to identify the AlgU regulon when the bacteria are inside the plant host and which PTI-regulated genes overlapped with AlgU-regulated genes. In this study, we analyzed transcriptomic data from RNA-sequencing to identify the AlgU regulon (while in the host) and its relationship with PTI. Our results showed that the upregulation of 224 genes while inside the plant host require AlgU, while another 154 genes are downregulated dependent on AlgU in Arabidopsis during early infection. Both stress response and virulence-associated genes were upregulated in a manner dependent on AlgU, while the flagellar motility genes are downregulated in a manner dependent on AlgU. Under the pre-induced PTI condition, more than half of these AlgU-regulated genes have lost induction/suppression in contrast to mock treated plants, and almost all function groups regulated by AlgU were affected by PTI.

Population transcriptomics reveals the effect of gene flow on the evolution of range limits

One of the most important questions in evolutionary biology is how the spatial distribution of species is limited. Asymmetric gene flow from core populations is suggested to increase the number of poorly adapted immigrants in the populations at the range edge. Genetic load due to migration, i.e., migration load, should prevent adaptation to the local habitat, leading to decreases in distribution range via local extinction or the limiting range expansion. However, few experimental studies have examined the effects of immigration on fitness and natural selection within recipient populations. To investigate the influence of migration load on the evolution of distribution range, we performed field and laboratory observations as well as population transcriptomics for the common river snail, Semisulcospira reiniana. This species meets the conditions that migration from source populations can prevent local adaptation in a sink population because they inhabit the broader range of environments, including middle/upper reaches of a river and estuaries within a single river and they may be more vulnerable to being swept away by water currents due to lowered spontaneous (upward) locomotion activity. We found that river steepness was related to the lower distribution limit of S. reiniana, with a narrower distribution range in the steeper river. Population transcriptomic analysis showed that gene flow was heavily asymmetric from the upstream populations to downstream ones in the steep river, suggesting a greater migration load in the steep river. The number of genes putatively involved in adaptation to the local habitat was lower in the steep river than in the gentle river. Gene expression profiles suggested that individuals achieve better local adaptation in the gentle river. Laboratory experiments suggested that evolutionary differences in salinity tolerance among local populations were only found in the gentle river. Our results consistent with the hypothesis that migration load owing to asymmetric gene flow disturbs local adaptation and restricts the distribution range of river snails.

Protocol for using NoBadWordsCombiner to merge and minimize "bad words" from BLAST hits against multiple eukaryotic gene annotation databases

Annotating protein-coding genes can be challenging, especially when searching for the best hits against multiple functional databases. This is partly because of "bad words" appearing as top hits, such as hypothetical or uncharacterized proteins. To help alleviate some of these issues, we designed a bioinformatics tool called NoBadWordsCombiner, which efficiently merges the hits from various databases, strengthening gene definitions by minimizing functional descriptions containing "bad words." Unlike other available tools, NoBadWordsCombiner is user friendly, but it does require users to have some general bioinformatics skills, including a basic understanding of the BLAST package and dash shell in Linux/Unix environments. For complete details on the use and execution of this protocol, please refer to Zhang et al. (2021a).

Spatially Resolved Immunometabolism to Understand Infectious Disease Progression

Infectious diseases, including those of viral, bacterial, fungal, and parasitic origin are often characterized by focal inflammation occurring in one or more distinct tissues. Tissue-specific outcomes of infection are also evident in many infectious diseases, suggesting that the local microenvironment may instruct complex and diverse innate and adaptive cellular responses resulting in locally distinct molecular signatures. In turn, these molecular signatures may both drive and be responsive to local metabolic changes in immune as well as non-immune cells, ultimately shaping the outcome of infection. Given the spatial complexity of immune and inflammatory responses during infection, it is evident that understanding the spatial organization of transcripts, proteins, lipids, and metabolites is pivotal to delineating the underlying regulation of local immunity. Molecular imaging techniques like mass spectrometry imaging and spatially resolved, highly multiplexed immunohistochemistry and transcriptomics can define detailed metabolic signatures at the microenvironmental level. Moreover, a successful complementation of these two imaging techniques would allow multi-omics analyses of inflammatory microenvironments to facilitate understanding of disease pathogenesis and identify novel targets for therapeutic intervention. Here, we describe strategies for downstream data analysis of spatially resolved multi-omics data and, using leishmaniasis as an exemplar, describe how such analysis can be applied in a disease-specific context.

Elucidating the biochemical basis of trans-16:1 fatty acid change in leaves during cold acclimation in wheat

In plant cells, phosphatidylglycerol (PG) in the chloroplast has a characteristic trans-∆3-hexadecenoic acid (t16:1) at the sn-2 position. The t16:1 content in wheat leaf tissues decreases during cold treatment, but the significance of this fatty acid compositional change and the underlying biochemical mechanism remains poorly understood. Using a large collection of wheat cultivars displaying a varying capacity of freezing tolerance, we show for the first time under field conditions that this low temperature induced t16:1 change is associated with winter hardiness. To explore the metabolic mechanism responsible for the reduction of t16:1, we performed detailed lipid analysis and comparative transcriptome study with four selected wheat lines under cold acclimation. Our results show that wheat leaf tissues experience a gradual decrease in chloroplast lipid pathway activity during cold acclimation and that the decline in chloroplast lipid synthesis manifests itself in the decrease of t16:1 in PG. Comparative RNA-seq analyses with leaf tissues further reveal concerted transcriptome shifts indicating a rebalancing of chloroplast and cytosolic lipid synthesis during cold acclimation. Our study, thus, provides mechanistic understanding on chloroplast lipid adjustments as a "molecular ideotype" and the t16:1 change as a specific metabolite marker for screening freezing tolerance in wheat.

Finding Nemo's clock reveals switch from nocturnal to diurnal activity

Timing mechanisms play a key role in the biology of coral reef fish. Typically, fish larvae leave their reef after hatching, stay for a period in the open ocean before returning to the reef for settlement. During this dispersal, larvae use a time-compensated sun compass for orientation. However, the timing of settlement and how coral reef fish keep track of time via endogenous timing mechanisms is poorly understood. Here, we have studied the behavioural and genetic basis of diel rhythms in the clown anemonefish Amphiprion ocellaris. We document a behavioural shift from nocturnal larvae to diurnal adults, while juveniles show an intermediate pattern of activity which potentially indicates flexibility in the timing of settlement on a host anemone. qRTPCR analysis of six core circadian clock genes (bmal1, clocka, cry1b, per1b, per2, per3) reveals rhythmic gene expression patterns that are comparable in larvae and juveniles, and so do not reflect the corresponding activity changes. By establishing an embryonic cell line, we demonstrate that clown anemonefish possess an endogenous clock with similar properties to that of the zebrafish circadian clock. Furthermore, our study provides a first basis to study the multi-layered interaction of clocks from fish, anemones and their zooxanthellae endosymbionts.

Systems Wide Analysis of CCM Signaling Complex Alterations in CCM-Deficient Models Using Omics Approaches

Omics research has garnered popularity recently to integrate in-depth analysis of alterations at the molecular level to elucidate observable phenotypes resulting from knockdown/knockout models. Genomics, performed through RNA-seq, allows the user to evaluate alterations at the transcription level, oftentimes more sensitive than other types of analysis, especially when attempting to understand lack of observation of an expected phenotype. Proteomics facilitates an understanding of mechanisms being altered at the translational level allowing for an understanding of multiple layers of regulation occurring, elucidating discrepancies between what is seen at the RNA level compared to what is translated to a functional protein. Here we describe the methods currently being used to evaluate CCM-deficient strains in human brain microvascular endothelial cells (HBMVEC), zebrafish embryos as well as in vivo mouse model to evaluate impacts on various signaling cascades resulting from deficiencies in KRIT1 (CCM1), MGC4607 (CCM2), and PDCD10 (CCM3). The integration of data from genomics and proteomics analysis allows for the composition of interactomes, elucidating systems wide impacts resulting from disruption of the CCM signaling complex (CSC).

Protein Complexes Form a Basis for Complex Hybrid Incompatibility

Proteins are the workhorses of the cell and execute many of their functions by interacting with other proteins forming protein complexes. Multi-protein complexes are an admixture of subunits, change their interaction partners, and modulate their functions and cellular physiology in response to environmental changes. When two species mate, the hybrid offspring are usually inviable or sterile because of large-scale differences in the genetic makeup between the two parents causing incompatible genetic interactions. Such reciprocal-sign epistasis between inter-specific alleles is not limited to incompatible interactions between just one gene pair; and, usually involves multiple genes. Many of these multi-locus incompatibilities show visible defects, only in the presence of all the interactions, making it hard to characterize. Understanding the dynamics of protein-protein interactions (PPIs) leading to multi-protein complexes is better suited to characterize multi-locus incompatibilities, compared to studying them with traditional approaches of genetics and molecular biology. The advances in omics technologies, which includes genomics, transcriptomics, and proteomics can help achieve this end. This is especially relevant when studying non-model organisms. Here, we discuss the recent progress in the understanding of hybrid genetic incompatibility; omics technologies, and how together they have helped in characterizing protein complexes and in turn multi-locus incompatibilities. We also review advances in bioinformatic techniques suitable for this purpose and propose directions for leveraging the knowledge gained from model-organisms to identify genetic incompatibilities in non-model organisms.

Signals of selection beyond bottlenecks between exotic populations of the bull-headed dung beetle, Onthophagus taurus

Colonization of new environments can lead to population bottlenecks and rapid phenotypic evolution that could be due to neutral and selective processes. Exotic populations of the bull-headed dung beetle (Onthophagus taurus) have differentiated in opposite directions from native beetles in male horn-to-body size allometry and female fecundity. Here we test for genetic and transcriptional differences among two exotic and one native O. taurus populations after three generations in common garden conditions. We sequenced RNA from 24 individuals for each of the three populations including both sexes, and spanning four developmental stages for the two exotic, differentiated populations. Identifying 270,400 high-quality single nucleotide polymorphisms, we revealed a strong signal of genetic differentiation between the three populations, and evidence of recent bottlenecks within and an excess of outlier loci between exotic populations. Differences in gene expression between populations were greatest in prepupae and early adult life stages, stages during which differences in male horn development and female fecundity manifest. Finally, genes differentially expressed between exotic populations also had greater genetic differentiation and performed functions related to chitin biosynthesis and nutrient sensing, possibly underlying allometry and fecundity trait divergences. Our results suggest that beyond bottlenecks, recent introductions have led to genetic and transcriptional differences in genes correlated with observed phenotypic differences.

Sequence analysis of the Petunia inflata S-locus region containing 17 S-Locus F-Box genes and the S-RNase gene involved in self-incompatibility

Self-incompatibility in Petunia is controlled by the polymorphic S-locus, which contains S-RNase encoding the pistil determinant and 16-20 S-locus F-box (SLF) genes collectively encoding the pollen determinant. Here we sequenced and assembled approximately 3.1 Mb of the S₂ -haplotype of the S-locus in Petunia inflata using bacterial artificial chromosome clones collectively containing all 17 SLF genes, SLFLike1, and S-RNase. Two SLF pseudogenes and 28 potential protein-coding genes were identified, 20 of which were also found at the S-loci of both the S_6a -haplotype of P. inflata and the S_N -haplotype of self-compatible Petunia axillaris, but not in the S-locus remnants of self-compatible potato (Solanum tuberosum) and tomato (Solanum lycopersicum). Comparative analyses of S-locus sequences of these three S-haplotypes revealed potential genetic exchange in the flanking regions of SLF genes, resulting in highly similar flanking regions between different types of SLF and between alleles of the same type of SLF of different S-haplotypes. The high degree of sequence similarity in the flanking regions could often be explained by the presence of similar long terminal repeat retroelements, which were enriched at the S-loci of all three S-haplotypes and in the flanking regions of all S-locus genes examined. We also found evidence of the association of transposable elements with SLF pseudogenes. Based on the hypothesis that SLF genes were derived by retrotransposition, we identified 10 F-box genes as putative SLF parent genes. Our results shed light on the importance of non-coding sequences in the evolution of the S-locus, and on possible evolutionary mechanisms of generation, proliferation, and deletion of SLF genes.

Rapid Adaptation to Temperature via a Potential Genomic Island of Divergence in the Invasive Green Crab, Carcinus maenas

Widespread species often adapt easily to novel conditions – both those found in new habitats and those generated by climate change. However, rapid adaptation may be hindered in the marine realm, where long-distance dispersal and consequently high gene flow are predicted to limit potential for local adaptation. Here, we use a highly dispersive invasive marine crab to test the nature and speed of adaptation to temperature in the sea. Using single nucleotide polymorphisms (SNPs) generated from cardiac transcriptome sequencing, we characterized six populations of the European green crab (Carcinus maenas) located across parallel thermal gradients in their native and invasive ranges. We compared SNP frequencies with local temperatures and previously generated data on cardiac heat and cold tolerance to identify candidate markers associated with population-level differences in thermal physiology. Of 10,790 SNPs, 104 were identified as frequency outliers, a signal that was strongly driven by association with temperature and/or cold tolerance. Seventy-two of these outlier markers, representing 28 different genes, were in a cluster of SNPs identified as a potential inversion polymorphism using linkage disequilibrium network analysis. This SNP cluster was unique in the data set, which was otherwise characterized by low levels of linkage disequilibrium, and markers in this cluster showed a significant enrichment of coding substitutions relative to the full SNP set. These 72 outlier SNPs appear to be transmitted as a unit, and represent a putative genomic island of divergence which varied in frequency with organismal cold tolerance. This relationship was strikingly similar across both native and invasive populations, all of which showed a very strong correlation with cold tolerance (R2 = 0.96 over all six populations). Notably, three of these populations have diverged recently (<100 years) and show little to no neutral divergence, suggesting that this genomic region may be responding to temperature on a relatively short time scale. This relationship indicates adaptation to temperature based on the action of a putative genomic island of divergence, perhaps partially explaining the extraordinary invasive ability of this species.

Rapid evolution of genome-wide gene expression and plasticity during saline to freshwater invasions by the copepod Eurytemora affinis species complex

Saline migrants into freshwater habitats constitute among the most destructive invaders in aquatic ecosystems throughout the globe. However, the evolutionary and physiological mechanisms underlying such habitat transitions remain poorly understood. To explore the mechanisms of freshwater adaptation and distinguish between adaptive (evolutionary) and acclimatory (plastic) responses to salinity change, we examined genome-wide patterns of gene expression between ancestral saline and derived freshwater populations of the Eurytemora affinis species complex, reared under two different common-garden conditions (0 versus 15 PSU). We found that evolutionary shifts in gene expression (between saline and freshwater inbred lines) showed far greater changes and were more widespread than acclimatory responses to salinity (0 versus 15 PSU). Most notably, 30-40 genes showing evolutionary shifts in gene expression across the salinity boundary were associated with ion transport function, with inorganic cation transmembrane transport forming the largest Gene Ontology category. Of particular interest was the sodium transporter, the Na⁺ /H⁺ antiporter (NHA) gene family, which was discovered in animals relatively recently. Thirty key ion regulatory genes, such as NHA paralogue #7, demonstrated concordant evolutionary and plastic shifts in gene expression, suggesting the evolution of ion transporter function and plasticity during rapid invasions into novel salinities. Moreover, freshwater invasions were associated with the evolution of reduced plasticity in the freshwater population, again for the same key ion transporters, consistent with the predicted evolution of canalization following adaptation to stressful conditions. Our results have important implications for understanding evolutionary and physiological mechanisms of range expansions by some of the most widespread invaders in aquatic habitats.

Adaptive divergence, neutral panmixia, and algal symbiont population structure in the temperate coral Astrangia poculata along the Mid-Atlantic United States

Astrangia poculata is a temperate scleractinian coral that exists in facultative symbiosis with the dinoflagellate alga Breviolum psygmophilum across a range spanning the Gulf of Mexico to Cape Cod, Massachusetts. Our previous work on metabolic thermal performance of Virginia (VA) and Rhode Island (RI) populations of A. poculata revealed physiological signatures of cold (RI) and warm (VA) adaptation of these populations to their respective local thermal environments. Here, we used whole-transcriptome sequencing (mRNA-Seq) to evaluate genetic differences and identify potential loci involved in the adaptive signature of VA and RI populations. Sequencing data from 40 A. poculata individuals, including 10 colonies from each population and symbiotic state (VA-white, VA-brown, RI-white, and RI-brown), yielded a total of 1,808 host-associated and 59 algal symbiont-associated single nucleotide polymorphisms (SNPs) post filtration. Fst outlier analysis identified 66 putative high outlier SNPs in the coral host and 4 in the algal symbiont. Differentiation of VA and RI populations in the coral host was driven by putatively adaptive loci, not neutral divergence (Fst = 0.16, p = 0.001 and Fst = 0.002, p = 0.269 for outlier and neutral SNPs respectively). In contrast, we found evidence of neutral population differentiation in B. psygmophilum (Fst = 0.093, p = 0.001). Several putatively adaptive host loci occur on genes previously associated with the coral stress response. In the symbiont, three of four putatively adaptive loci are associated with photosystem proteins. The opposing pattern of neutral differentiation in B. psygmophilum, but not the A. poculata host, reflects the contrasting dynamics of coral host and algal symbiont population connectivity, dispersal, and gene by environment interactions.

An 85K SNP Array Uncovers Inbreeding and Cryptic Relatedness in an Antarctic Fur Seal Breeding Colony

High density single nucleotide polymorphism (SNP) arrays allow large numbers of individuals to be rapidly and cost-effectively genotyped at large numbers of genetic markers. However, despite being widely used in studies of humans and domesticated plants and animals, SNP arrays are lacking for most wild organisms. We developed a custom 85K Affymetrix Axiom array for an intensively studied pinniped, the Antarctic fur seal (Arctocephalus gazella). SNPs were discovered from a combination of genomic and transcriptomic resources and filtered according to strict criteria. Out of a total of 85,359 SNPs tiled on the array, 75,601 (88.6%) successfully converted and were polymorphic in 270 animals from a breeding colony at Bird Island in South Georgia. Evidence was found for inbreeding, with three genomic inbreeding coefficients being strongly intercorrelated and the proportion of the genome in runs of homozygosity being non-zero in all individuals. Furthermore, analysis of genomic relatedness coefficients identified previously unknown first-degree relatives and multiple second-degree relatives among a sample of ostensibly unrelated individuals. Such “cryptic relatedness” within fur seal breeding colonies may increase the likelihood of consanguineous matings and could therefore have implications for understanding fitness variation and mate choice. Finally, we demonstrate the cross-amplification potential of the array in three related pinniped species. Overall, our SNP array will facilitate future studies of Antarctic fur seals and has the potential to serve as a more general resource for the wider pinniped research community.

Genomic, Transcriptomic and Epigenomic Tools to Study the Domestication of Plants and Animals: A Field Guide for Beginners

In the last decade, genomics and the related fields of transcriptomics and epigenomics have revolutionized the study of the domestication process in plants and animals, leading to new discoveries and new unresolved questions. Given that some domesticated taxa have been more studied than others, the extent of genomic data can range from vast to nonexistent, depending on the domesticated taxon of interest. This review is meant as a rough guide for students and academics that want to start a domestication research project using modern genomic tools, as well as for researchers already conducting domestication studies that are interested in following a genomic approach and looking for alternate strategies (cheaper or more efficient) and future directions. We summarize the theoretical and technical background needed to carry out domestication genomics, starting from the acquisition of a reference genome and genome assembly, to the sampling design for population genomics, paleogenomics, transcriptomics, epigenomics and experimental validation of domestication-related genes. We also describe some examples of the aforementioned approaches and the relevant discoveries they made to understand the domestication of the studied taxa.

Molecular and biochemical responses of vitellogenin in the mussel Mytilus galloprovincialis exposed to the glyphosate-based herbicide Roundup® Power 2.0

Glyphosate-based herbicides (GBHs) occur in aquatic ecosystems at concentrations of hundreds of micrograms per liter. As formulation adjuvants are suspected to be endocrine-disrupting chemicals, we assessed the effects of the recent GBH formulation Roundup® Power 2.0 on vitellogenin (VTG) in Mytilus galloprovincialis. Mussels were exposed for 7, 14, and 21 days to two concentrations of the commercial formulation, corresponding to 100 and 1000 μg/L of glyphosate. The expression of the vtg gene in gonads of females and males, as well as the levels of alkali labile phosphates (ALP) in gonads and non-gonadal tissues from the two sexes were measured. No significant alterations were observed in vtg expression values during the exposure. Conversely, a significant reduction in gonadal ALP levels was observed in females exposed for 21 days and in males exposed for 7 days. In addition, ALP levels increased significantly in gonads from males exposed for 21 days to the two concentrations of Roundup®. As for non-gonadal tissues, ALP levels did not change significantly in females, whereas ALP levels decreased significantly in non-gonadal tissues from males exposed for 21 days to the lowest concentration tested. An overall statistically significant difference in ALP levels was found between females and males. Although preliminary, our study suggests that GBH can affect reproduction-related parameters in mussels.

Biogeochemical and physical controls on methane fluxes from two ferruginous meromictic lakes

Meromictic lakes with anoxic bottom waters often have active methane cycles whereby methane is generally produced biogenically under anoxic conditions and oxidized in oxic surface waters prior to reaching the atmosphere. Lakes that contain dissolved ferrous iron in their deep waters (i.e., ferruginous) are rare, but valuable, as geochemical analogues of the conditions that dominated the Earth's oceans during the Precambrian when interactions between the iron and methane cycles could have shaped the greenhouse regulation of the planet's climate. Here, we explored controls on the methane fluxes from Brownie Lake and Canyon Lake, two ferruginous meromictic lakes that contain similar concentrations (max. >1 mM) of dissolved methane in their bottom waters. The order Methanobacteriales was the dominant methanogen detected in both lakes. At Brownie Lake, methanogen abundance, an increase in methane concentration with respect to depths closer to the sediment, and isotopic data suggest methanogenesis is an active process in the anoxic water column. At Canyon Lake, methanogenesis occurred primarily in the sediment. The most abundant aerobic methane-oxidizing bacteria present in both water columns were associated with the Gammaproteobacteria, with little evidence of anaerobic methane oxidizing organisms being present or active. Direct measurements at the surface revealed a methane flux from Brownie Lake that was two orders of magnitude greater than the flux from Canyon Lake. Comparison of measured versus calculated turbulent diffusive fluxes indicates that most of the methane flux at Brownie Lake was non-diffusive. Although the turbulent diffusive methane flux at Canyon Lake was attenuated by methane oxidizing bacteria, dissolved methane was detected in the epilimnion, suggestive of lateral transport of methane from littoral sediments. These results highlight the importance of direct measurements in estimating the total methane flux from water columns, and that non-diffusive transport of methane may be important to consider from other ferruginous systems.

Transcriptome Characterization of Reverse Development in Turritopsis dohrnii (Hydrozoa, Cnidaria)

Medusae of Turritopsis dohrnii undergo reverse development in response to physical damage, adverse environmental conditions, or aging. Senescent, weakened or damaged medusae transform into a cluster of poorly differentiated cells (known as the cyst stage), which metamorphose back into a preceding life cycle stage, the polyp. During the metamorphosis, cell transdifferentiation occurs. The cyst represents the intermediate stage between a reverting medusa and a healthy polyp, during which cell transdifferentiation and tissue reorganization take place. Here we characterize and compare the transcriptomes of the polyp and newborn medusa stages of T. dohrnii with that of the cyst, to identify biological networks potentially involved in the reverse development and transdifferentiation processes. The polyp, medusa and cyst of T. dohrnii were sequenced through Illumina RNA-sequencing and assembled using a de novo approach, resulting in 92,569, 74,639 and 86,373 contigs, respectively. The transcriptomes were annotated and comparative analyses among the stages identified biological networks that were significantly over-and under-expressed in the cyst as compared to the polyp and medusa stages. Biological processes that occur at the cyst stage such as telomerase activity, regulation of transposable elements and DNA repair systems, and suppression of cell signaling pathways, mitotic cell division and cellular differentiation and development may be involved in T. dohrnii's reverse development and transdifferentiation. Our results are the first attempt to understand T. dohrnii's life-cycle reversal at the genetic level, and indicate possible avenues of future research on developmental strategies, cell transdifferentiation, and aging using T. dohrnii as a non-traditional in vivo system.

Transcriptome analyses provide the first insight into the molecular basis of cold tolerance in Larimichthys polyactis

Larimichthys polyactis is one of the most economically important marine fish species that have become newly cultured in China in recent years. The gene expression changes that L. polyactis experiences in cold toleranceis still unknown, limiting the expansion of its cultivation, fast growth, and high yield. To investigate the molecular mechanism behind L. polyactis's cold tolerance and to provide a resource for conducting genetic research on L. polyactis, transcriptome sequencing using RNA-seq was performed on individuals that survived cold stress at 4 °C (cold tolerant, CT), and individuals that barely survived 4 °C (cold sensitive, CS), which was considered as the control. A number of 387,607,550 clean reads were obtained from the transcriptomes, and comparative transcriptomic analysis identified 141 differently expressed genes (DEGs), of which 67 were up-regulated and 74 were down-regulated in CT compared to CS under cold stress. Furthermore, ten differently expressed genes were selected from the RNA-Seq analysis to be further validated by real-time PCR. Functional network analysis indicated that L. polyactis adapted to cold stress by employing a series of mechanisms to minimize damages caused by exposure to cold temperatures. The molecular mechanisms identified through RNA-Seq included Extracellular matrix (ECM) receptor interaction, glycerolipid metabolism, regulation of autophagy and focal adhesion pathway as playing vital roles in cold tolerance in L. polyactis. This study may help elucidate how L. polyactis tolerates cold, which is of value for breeding cold-tolerant L. polyactis stocks for cultivation.

Effects of a male meiotic driver on male and female transcriptomes in the house mouse

Not all genetic loci follow Mendel's rules, and the evolutionary consequences of this are not yet fully known. Genomic conflict involving multiple loci is a likely outcome, as restoration of Mendelian inheritance patterns will be selected for, and sexual conflict may also arise when sexes are differentially affected. Here, we investigate effects of the t haplotype, an autosomal male meiotic driver in house mice, on genome-wide gene expression patterns in males and females. We analysed gonads, liver and brain in adult same-sex sibling pairs differing in genotype, allowing us to identify t-associated differences in gene regulation. In testes, only 40% of differentially expressed genes mapped to the approximately 708 annotated genes comprising the t haplotype. Thus, much of the activity of the t haplotype occurs in trans, and as upregulation. Sperm maturation functions were enriched among both cis and trans acting t haplotype genes. Within the t haplotype, we observed more downregulation and differential exon usage. In ovaries, liver and brain, the majority of expression differences mapped to the t haplotype, and were largely independent of the differences seen in the testis. Overall, we found widespread transcriptional effects of this male meiotic driver in the house mouse genome.

Sex identification from distinctive gene expression patterns in Antarctic krill (Euphausia superba)

Antarctic krill (Euphausia superba) is a highly abundant keystone species of the Southern Ocean ecosystem, directly connecting primary producers to high-trophic level predators. Sex ratios of krill vary remarkably between swarms and this phenomenon is poorly understood, as identification of krill sex relies on external morphological differences that appear late during development. Sex determination mechanisms in krill are unknown, but could include genetic, environmental or parasitic mechanisms. Similarly, virtually nothing is known about molecular sex differentiation. The krill genome has to date not been sequenced, and due to its enormous size and large amount of repetitive elements, it is currently not feasible to develop sex-specific DNA markers. To produce a reliable molecular marker for sex in krill and to investigate molecular sex differentiation we therefore focused on identifying sex-specific transcriptomic differences. Through transcriptomic analysis, we found large gene expression differences between testes and ovaries and identified three genes exclusively expressed in female whole krill from early juvenile stages onwards. The sex-specific expression of these three genes persisted through sexual regression, although our regressed samples originated from a krill aquarium and may differ from wild-regressed krill. Two slightly male-biased genes did not display sufficient expression differences to clearly differentiate sexes. Based on the expression of the three female-specific genes we developed a molecular test that for the first time allows the unambiguous sex determination of krill samples lacking external sex-specific features from juvenile stages onwards, including the sexually regressed krill we examined.

Differential gene expression during substrate probing in larvae of the Caribbean coral Porites astreoides

The transition from larva to adult is a critical step in the life history strategy of most marine animals. However, the genetic basis of this life history change remains poorly understood in many taxa, including most coral species. Recent evidence suggests that coral planula larvae undergo significant changes at the physiological and molecular levels throughout the development. To investigate this, we characterized differential gene expression (DGE) during the transition from planula to adult polyp in the abundant Caribbean reef-building coral Porites astreoides, that is from nonprobing to actively substrate-probing larva, a stage required for colony initiation. This period is crucial for the coral, because it demonstrates preparedness to locate appropriate substrata for settlement based on vital environmental cues. Through RNA-Seq, we identified 860 differentially expressed holobiont genes between probing and nonprobing larvae (p ≤ .01), the majority of which were upregulated in probing larvae. Surprisingly, differentially expressed genes of endosymbiotic dinoflagellate origin greatly outnumbered coral genes, compared with a nearly 1:1 ratio of coral-to-dinoflagellate gene representation in the holobiont transcriptome. This unanticipated result suggests that dinoflagellate endosymbionts may play a significant role in the transition from nonprobing to probing behaviour in dinoflagellate-rich larvae. Putative holobiont genes were largely involved in protein and nucleotide binding, metabolism and transport. Genes were also linked to environmental sensing and response and integral signalling pathways. Our results thus provide detailed insight into molecular changes prior to larval settlement and highlight the complex physiological and biochemical changes that occur in early transition stages from pelagic to benthic stages in corals, and perhaps more importantly, in their endosymbionts.

Landscape of ROD9 Island: Functional annotations and biological network of hypothetical proteins in Salmonella enterica

Salmonella, an Enterobacteria is a therapeutically important pathogen for the host. The advancement of genome sequencing of S. enterica serovar Enteritidis have identified a distinct ROD9 pathogenic island, imparting virulence. The occurrence of 17 ROD9 hypothetical proteins, necessitates subsequent bioinformatics approach for structural and functional aspects of protein-protein relations or networks in different pathogenic phenotypes express. A collective analysis using predictive bioinformatics tools that includes NCBI-BLASTp and BLAST2GO annotated the motif patterns and functional significance. The VFDB identified 10 virulence proteins at both genomic and metagenomic level. Phylogenetic analysis revealed a divergent and convergent relationship between 17 ROD9 and 41 SP-1 proteins. Here, combining a comprehensive approach from sequence based, motif recognitions, domain identification, virulence ability to structural modelling provides a precise function to ROD9 proteins biological network, for which no experimental information is available.

Integrated Genomic and Transcriptomic Analysis of the Peridinin Dinoflagellate Amphidinium carterae Plastid

The plastid genomes of peridinin-containing dinoflagellates are highly unusual, possessing very few genes, which are located on small chromosomal elements termed "minicircles". These minicircles may contain genes, or no recognisable coding information. Transcripts produced from minicircles may undergo unusual processing events, such as the addition of a 3' poly(U) tail. To date, little is known about the genetic or transcriptional diversity of non-coding sequences in peridinin dinoflagellate plastids. These sequences include empty minicircles, and regions of non-coding DNA in coding minicircles. Here, we present an integrated plastid genome and transcriptome for the model peridinin dinoflagellate Amphidinium carterae, identifying a previously undescribed minicircle. We also profile transcripts covering non-coding regions of the psbA and petB/atpA minicircles. We present evidence that antisense transcripts are produced within the A. carterae plastid, but show that these transcripts undergo different end cleavage events from sense transcripts, and do not receive 3' poly(U) tails. The difference in processing events between sense and antisense transcripts may enable the removal of non-coding transcripts from peridinin dinoflagellate plastid transcript pools.

Transcriptomic analysis of Pacific white shrimp (Litopenaeus vannamei, Boone 1931) in response to acute hepatopancreatic necrosis disease caused by Vibrio parahaemolyticus

Acute hepatopancreatic necrosis disease (AHPND), caused by marine bacteria Vibrio Parahaemolyticus, is a huge problem in shrimp farms. The V. parahaemolyticus infecting material is contained in a plasmid which encodes for the lethal toxins PirABVp, whose primary target tissue is the hepatopancreas, causing sloughing of epithelial cells, necrosis, and massive hemocyte infiltration. To get a better understanding of the hepatopancreas response during AHPND, juvenile shrimp Litopenaeus vannamei were infected by immersion with V. parahaemolyticus. We performed transcriptomic mRNA sequencing of infected shrimp hepatopancreas, at 24 hours post-infection, to identify novel differentially expressed genes a total of 174,098 transcripts were examined of which 915 transcripts were found differentially expressed after comparative transcriptomic analysis: 442 up-regulated and 473 down-regulated transcripts. Gene Ontology term enrichment analysis for up-regulated transcripts includes metabolic process, regulation of programmed cell death, carbohydrate metabolic process, and biological adhesion, whereas for down-regulated transcripts include, microtubule-based process, cell activation, and chitin metabolic process. The analysis of protein- protein network between up and down-regulated genes indicates that the first gene interactions are connected to oxidation-processes and sarcomere organization. Additionally, protein-protein networks analysis identified 20-top highly connected hub nodes. Based on their immunological or metabolic function, ten candidate transcripts were selected to measure their mRNA relative expression levels in AHPND infected shrimp hepatopancreas by RT-qPCR. Our results indicate a close connection between the immune and metabolism systems during AHPND infection. Our RNA-Seq and RT-qPCR data provide the possible immunological and physiological scenario as well as the molecular pathways that take place in the shrimp hepatopancreas in response to an infectious disease.

Small- and large-scale heterogeneity in genetic variation across the collard flycatcher genome: implications for estimating genetic diversity in nonmodel organisms

Population genetic studies in nonmodel organisms are often hampered by a lack of reference genomes that are essential for whole-genome resequencing. In the light of this, genotyping methods have been developed to effectively eliminate the need for a reference genome, such as genotyping by sequencing or restriction site-associated DNA sequencing (RAD-seq). However, what remains relatively poorly studied is how accurately these methods capture both average and variation in genetic diversity across an organism's genome. In this issue of Molecular Ecology Resources, Dutoit et al. (2016) use whole-genome resequencing data from the collard flycatcher to assess what factors drive heterogeneity in nucleotide diversity across the genome. Using these data, they then simulate how well different sequencing designs, including RAD sequencing, could capture most of the variation in genetic diversity. They conclude that for evolutionary and conservation-related studies focused on the estimating genomic diversity, researchers should emphasize the number of loci analysed over the number of individuals sequenced.

Identification and quantification of domestic feline cytochrome P450 transcriptome across multiple tissues

Understanding of cytochrome P450 (CYP) isoform distribution and function in the domestic feline is limited. Only a few studies have defined individual CYP isoforms across metabolically relevant tissues, hampering the ability to predict drug metabolism and potential drug-drug interactions. Using RNA sequencing (RNA-seq), transcriptomes from the 99 Lives Cat Genome Sequencing Initiative databank combined with experimentally acquired whole transcriptome sequencing of healthy, adult male (n = 2) and female (n = 2) domestic felines, expression of 42 CYP isoforms were identified in 20 different tissues. Thirty-seven of these isoforms had not been previously reported in cats. Depending on the tissue, three to twenty-nine CYP isoform transcripts were expressed. The feline genome annotations did not differentiate CYP2E1 and 2E2 genes, demonstrating poor annotation for this gene using the reference genome. As the majority of the sequences are based on automated pipelines, complete cDNA sequences for translation into CYP protein sequences could not be determined. This study is the first to identify and characterize 37 additional CYP isoforms in feline tissues, increasing the number of identified CYP from the previously reported seven isoforms to 42 across 20 tissues.

Nonrandom RNAseq gene expression associated with RNAlater and flash freezing storage methods

RNA sequencing is a popular next-generation sequencing technique for assaying genome-wide gene expression profiles. Nonetheless, it is susceptible to biases that are introduced by sample handling prior gene expression measurements. Two of the most common methods for preserving samples in both field-based and laboratory conditions are submersion in RNAlater and flash freezing in liquid nitrogen. Flash freezing in liquid nitrogen can be impractical, particularly for field collections. RNAlater is a solution for stabilizing tissue for longer-term storage as it rapidly permeates tissue to protect cellular RNA. In this study, we assessed genome-wide expression patterns in 30-day-old fry collected from the same brood at the same time point that were flash-frozen in liquid nitrogen and stored at -80°C or submerged and stored in RNAlater at room temperature, simulating conditions of fieldwork. We show that sample storage is a significant factor influencing observed differential gene expression. In particular, genes with elevated GC content exhibit higher observed expression levels in liquid nitrogen flash-freezing relative to RNAlater storage. Further, genes with higher expression in RNAlater relative to liquid nitrogen experience disproportionate enrichment for functional categories, many of which are involved in RNA processing. This suggests that RNAlater may elicit a physiological response that has the potential to bias biological interpretations of expression studies. The biases introduced to observed gene expression arising from mimicking many field-based studies are substantial and should not be ignored.

Gene expression in the phenotypically plastic Arctic charr (Salvelinus alpinus): A focus on growth and ossification at early stages of development

Gene expression during development shapes the phenotypes of individuals. Although embryonic gene expression can have lasting effects on developmental trajectories, few studies consider the role of maternal effects, such as egg size, on gene expression. Using qPCR, we characterize relative expression of 14 growth and/or skeletal promoting genes across embryonic development in Arctic charr (Salvelinus alpinus). We test to what extent their relative expression is correlated with egg size and size at early life‐stages within the study population. We predict smaller individuals to have higher expression of growth and skeletal promoting genes, due to less maternal resources (i.e., yolk) and prioritization of energy toward ossification. We found expression levels to vary across developmental stages and only three genes (Mmp9, Star, and Sgk1) correlated with individual size at a given developmental stage. Contrary to our hypothesis, expression of Mmp9 and Star showed a non‐linear relationship with size (at post fertilization and hatching, respectively), whilst Sgk1 was higher in larger embryos at hatching. Interestingly, these genes are also associated with craniofacial divergence of Arctic charr morphs. Our results indicate that early life‐stage variation in gene expression, concomitant to maternal effects, can influence developmental plasticity and potentially the evolution of resource polymorphism in fishes.

Diet-dependent gene expression highlights the importance of Cytochrome P450 in detoxification of algal secondary metabolites in a marine isopod

Isopods of the genus Idotea have an unusual ability to feed on algae containing high amounts of chemical defense molecules, such as species of the genera Fucus and Ulva. In this study, we compared gene expression patterns of Idotea balthica individuals fed with Fucus vesiculosus to individuals fed with Ulva lactuca. We generated the first-ever transcriptome assembly for this species, and found 3,233 differentially expressed genes across feeding regimes. However, only a handful of biological functions were enriched with regard to differentially expressed genes, the most notable being "alkaloid metabolic process". Within this category, we found eight differentially expressed cytochrome P450 (CYP) unigenes, all of which had a higher expression in the U. lactuca diet treatment. A phylogenetic analysis showed that the differentially expressed CYP genes are closely related to a CYP gene described from the hepatopancreas of the spiny lobster Panulirus argus, and we hypothesize that these transcripts are involved in metabolite detoxification. This is a first step in the understanding of this algae-grazer interaction, and will form a basis for future work to characterize cytochrome P450 functioning in marine crustaceans.

Integrated transcriptome and hormone profiling highlight the role of multiple phytohormone pathways in wheat resistance against fusarium head blight

Fusarium head blight (FHB or scab) caused by Fusarium spp. is a destructive disease of wheat. Since the most effective sources of FHB resistance are typically associated with unfavorable agronomic traits, breeding commercial cultivars that combine desired agronomic traits and a high level of FHB resistance remains a considerable challenge. A better understanding of the molecular mechanisms governing FHB resistance will help to design more efficient and precise breeding strategies. Here, multiple molecular tools and assays were deployed to compare the resistant variety Sumai3 with three regionally adapted Canadian cultivars. Macroscopic and microscopic disease evaluation established the relative level of Type II FHB resistance of the four varieties and revealed that the F. graminearum infection process displayed substantial temporal differences among organs. The rachis was found to play a critical role in preventing F. graminearum spread within spikes. Large-scale, organ-specific RNA-seq at different times after F. graminearum infection demonstrated that diverse defense mechanisms were expressed faster and more intensely in the spikelet of resistant varieties. The roles of plant hormones during the interaction of wheat with F. graminearum was inferred based on the transcriptomic data obtained and the quantification of the major plant hormones. Salicylic acid and jasmonic acid were found to play predominantly positive roles in FHB resistance, whereas auxin and ABA were associated with susceptibility, and ethylene appeared to play a dual role during the interaction with F graminearum.

Microbial metabolism of methanol and methylamine in the Gulf of Mexico: insight into marine carbon and nitrogen cycling

One carbon (C1) metabolism plays an important role in marine carbon cycling but the dynamics and modes of C1 transformations are not fully understood. We made contemporaneous measurements of methylamine and methanol metabolism to elucidate the role of C1 compounds as sources of carbon, energy and nitrogen. Methanol and methylamine were predominantly used as an energy source in offshore waters (oxidation rate constant: k_methanol : 0.02-0.10 day^-1 ; k_methylamine : 0.01-0.18 day^-1 ), but were also important sources of biomass carbon in coastal waters (assimilation rate constant: k_methanol : 0.04-0.10 day^-1 ; k_methylamine : 0.01-0.05 day^-1 ). The relative extent of assimilation versus oxidation for these substrates correlated positively with chlorophyll, nutrients and heterotrophic bacterial production. Methanol oxidation and assimilation were stimulated significantly by nutrient addition. In contrast, methylamine metabolism was inhibited by ammonium or nitrate, suggesting that methylamine served as a nitrogen source. A preliminary metagenomic survey revealed a diverse population of putative C1-utilizing microorganisms. These results show that the remineralization of methylamine could provide both C and N sources for microbes. Both methanol and methylamine contribute to microbial energetic and carbon substrate demands with a distinctly different signature in nearshore versus offshore environments.

Transcriptome dynamics associated with resistance and susceptibility against fusarium head blight in four wheat genotypes

Background

Fusarium head blight (FHB) of wheat in North America is caused mostly by the fungal pathogen Fusarium graminearum (Fg). Upon exposure to Fg, wheat initiates a series of cellular responses involving massive transcriptional reprogramming. In this study, we analyzed transcriptomics data of four wheat genotypes (Nyubai, Wuhan 1, HC374, and Shaw), at 2 and 4 days post inoculation (dpi) with Fg, using RNA-seq technology.

Results

A total of 37,772 differentially expressed genes (DEGs) were identified, 28,961 from wheat and 8811 from the pathogen. The susceptible genotype Shaw exhibited the highest number of host and pathogen DEGs, including 2270 DEGs associating with FHB susceptibility. Protein serine/threonine kinases and LRR-RK were associated with susceptibility at 2 dpi, while several ethylene-responsive, WRKY, Myb, bZIP and NAC-domain containing transcription factors were associated with susceptibility at 4 dpi. In the three resistant genotypes, 220 DEGs were associated with resistance. Glutathione S-transferase (GST), membrane proteins and distinct LRR-RKs were associated with FHB resistance across the three genotypes. Genes with unique, high up-regulation by Fg in Wuhan 1 were mostly transiently expressed at 2 dpi, while many defense-associated genes were up-regulated at both 2 and 4 dpi in Nyubai; the majority of unique genes up-regulated in HC374 were detected at 4 dpi only. In the pathogen, most genes showed increased expression between 2 and 4 dpi in all genotypes, with stronger levels in the susceptible host; however two pectate lyases and a hydrolase were expressed higher at 2 dpi, and acetyltransferase activity was highly enriched at 4 dpi.

Conclusions

There was an early up-regulation of LRR-RKs, different between susceptible and resistant genotypes; subsequently, distinct sets of genes associated with defense response were up-regulated. Differences in expression profiles among the resistant genotypes indicate genotype-specific defense mechanisms. This study also shows a greater resemblance in transcriptomics of HC374 to Nyubai, consistent with their sharing of two FHB resistance QTLs on 3BS and 5AS, compared to Wuhan 1 which carries one QTL on 2DL in common with HC374.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5012-3) contains supplementary material, which is available to authorized users.

Transcriptome dynamics associated with resistance and susceptibility against fusarium head blight in four wheat genotypes

BACKGROUND

Fusarium head blight (FHB) of wheat in North America is caused mostly by the fungal pathogen Fusarium graminearum (Fg). Upon exposure to Fg, wheat initiates a series of cellular responses involving massive transcriptional reprogramming. In this study, we analyzed transcriptomics data of four wheat genotypes (Nyubai, Wuhan 1, HC374, and Shaw), at 2 and 4 days post inoculation (dpi) with Fg, using RNA-seq technology.

RESULTS

A total of 37,772 differentially expressed genes (DEGs) were identified, 28,961 from wheat and 8811 from the pathogen. The susceptible genotype Shaw exhibited the highest number of host and pathogen DEGs, including 2270 DEGs associating with FHB susceptibility. Protein serine/threonine kinases and LRR-RK were associated with susceptibility at 2 dpi, while several ethylene-responsive, WRKY, Myb, bZIP and NAC-domain containing transcription factors were associated with susceptibility at 4 dpi. In the three resistant genotypes, 220 DEGs were associated with resistance. Glutathione S-transferase (GST), membrane proteins and distinct LRR-RKs were associated with FHB resistance across the three genotypes. Genes with unique, high up-regulation by Fg in Wuhan 1 were mostly transiently expressed at 2 dpi, while many defense-associated genes were up-regulated at both 2 and 4 dpi in Nyubai; the majority of unique genes up-regulated in HC374 were detected at 4 dpi only. In the pathogen, most genes showed increased expression between 2 and 4 dpi in all genotypes, with stronger levels in the susceptible host; however two pectate lyases and a hydrolase were expressed higher at 2 dpi, and acetyltransferase activity was highly enriched at 4 dpi.

CONCLUSIONS

There was an early up-regulation of LRR-RKs, different between susceptible and resistant genotypes; subsequently, distinct sets of genes associated with defense response were up-regulated. Differences in expression profiles among the resistant genotypes indicate genotype-specific defense mechanisms. This study also shows a greater resemblance in transcriptomics of HC374 to Nyubai, consistent with their sharing of two FHB resistance QTLs on 3BS and 5AS, compared to Wuhan 1 which carries one QTL on 2DL in common with HC374.

Comparative analysis of whole flower transcriptomes in the Zingiberales

The advancement of next generation sequencing technologies (NGS) has revolutionized our ability to generate large quantities of data at a genomic scale. Despite great challenges, these new sequencing technologies have empowered scientists to explore various relevant biological questions on non-model organisms, even in the absence of a complete sequenced reference genome. Here, we analyzed whole flower transcriptome libraries from exemplar species across the monocot order Zingiberales, using a comparative approach in order to gain insight into the evolution of the molecular mechanisms underlying flower development in the group. We identified 4,153 coding genes shared by all floral transcriptomes analyzed, and 1,748 genes that are only retrieved in the Zingiberales. We also identified 666 genes that are unique to the ginger lineage, and 2,001 that are only found in the banana group, while in the outgroup species Dichorisandra thyrsiflora J.C. Mikan (Commelinaceae) we retrieved 2,686 unique genes. It is possible that some of these genes underlie lineage-specific molecular mechanisms of floral diversification. We further discuss the nature of these lineage-specific datasets, emphasizing conserved and unique molecular processes with special emphasis in the Zingiberales. We also briefly discuss the strengths and shortcomings of de novo assembly for the study of developmental processes across divergent taxa from a particular order. Although this comparison is based exclusively on coding genes, with particular emphasis in transcription factors, we believe that the careful study of other regulatory mechanisms, such as non-coding RNAs, might reveal new levels of complexity, which were not explored in this work.

Genetic and genomic monitoring with minimally invasive sampling methods

The decreasing cost and increasing scope and power of emerging genomic technologies are reshaping the field of molecular ecology. However, many modern genomic approaches (e.g., RAD-seq) require large amounts of high-quality template DNA. This poses a problem for an active branch of conservation biology: genetic monitoring using minimally invasive sampling (MIS) methods. Without handling or even observing an animal, MIS methods (e.g., collection of hair, skin, faeces) can provide genetic information on individuals or populations. Such samples typically yield low-quality and/or quantities of DNA, restricting the type of molecular methods that can be used. Despite this limitation, genetic monitoring using MIS is an effective tool for estimating population demographic parameters and monitoring genetic diversity in natural populations. Genetic monitoring is likely to become more important in the future as many natural populations are undergoing anthropogenically driven declines, which are unlikely to abate without intensive adaptive management efforts that often include MIS approaches. Here, we profile the expanding suite of genomic methods and platforms compatible with producing genotypes from MIS, considering factors such as development costs and error rates. We evaluate how powerful new approaches will enhance our ability to investigate questions typically answered using genetic monitoring, such as estimating abundance, genetic structure and relatedness. As the field is in a period of unusually rapid transition, we also highlight the importance of legacy data sets and recommend how to address the challenges of moving between traditional and next-generation genetic monitoring platforms. Finally, we consider how genetic monitoring could move beyond genotypes in the future. For example, assessing microbiomes or epigenetic markers could provide a greater understanding of the relationship between individuals and their environment.

RNA sequencing, de novo assembly and differential analysis of the gill transcriptome of freshwater climbing perch Anabas testudineus after 6 days of seawater exposure

To obtain transcriptomic insights into branchial responses to salinity challenge in Anabas testudineus, this study employed RNA sequencing (RNA-Seq) to analyse the gill transcriptome of A. testudineus exposed to seawater (SW) for 6 days compared with the freshwater (FW) control group. A combined FW and SW gill transcriptome was de novo assembled from 169.9 million 101 bp paired-end reads. In silico validation employing 17 A. testudineus Sanger full-length coding sequences showed that 15/17 of them had greater than 80% of their sequences aligned to the de novo assembled contigs where 5/17 had their full-length (100%) aligned and 9/17 had greater than 90% of their sequences aligned. The combined FW and SW gill transcriptome was mapped to 13,780 unique human identifiers at E-value ≤1.0E-20 while 952 and 886 identifiers were determined as up and down-regulated by 1.5 fold, respectively, in the gills of A. testudineus in SW when compared with FW. These genes were found to be associated with at least 23 biological processes. A larger proportion of genes encoding enzymes and transporters associated with molecular transport, energy production, metabolisms were up-regulated, while a larger proportion of genes encoding transmembrane receptors, G-protein coupled receptors, kinases and transcription regulators associated with cell cycle, growth, development, signalling, morphology and gene expression were relatively lower in the gills of A. testudineus in SW when compared with FW. High correlation (R = 0.99) was observed between RNA-Seq data and real-time quantitative PCR validation for 13 selected genes. The transcriptomic sequence information will facilitate development of molecular resources and tools while the findings will provide insights for future studies into branchial iono-osmoregulation and related cellular processes in A. testudineus.

A Multi-Omics Analysis Pipeline for the Metabolic Pathway Reconstruction in the Orphan Species Quercus ilex

Holm oak (Quercus ilex) is the most important and representative species of the Mediterranean forest and of the Spanish agrosilvo-pastoral "dehesa" ecosystem. Despite its environmental and economic interest, Holm oak is an orphan species whose biology is very little known, especially at the molecular level. In order to increase the knowledge on the chemical composition and metabolism of this tree species, the employment of a holistic and multi-omics approach, in the Systems Biology direction would be necessary. However, for orphan and recalcitrant plant species, specific analytical and bioinformatics tools have to be developed in order to obtain adequate quality and data-density before to coping with the study of its biology. By using a plant sample consisting of a pool generated by mixing equal amounts of homogenized tissue from acorn embryo, leaves, and roots, protocols for transcriptome (NGS-Illumina), proteome (shotgun LC-MS/MS), and metabolome (GC-MS) studies have been optimized. These analyses resulted in the identification of around 62629 transcripts, 2380 protein species, and 62 metabolites. Data are compared with those reported for model plant species, whose genome has been sequenced and is well annotated, including Arabidopsis, japonica rice, poplar, and eucalyptus. RNA and protein sequencing favored each other, increasing the number and confidence of the proteins identified and correcting erroneous RNA sequences. The integration of the large amount of data reported using bioinformatics tools allows the Holm oak metabolic network to be partially reconstructed: from the 127 metabolic pathways reported in KEGG pathway database, 123 metabolic pathways can be visualized when using the described methodology. They included: carbohydrate and energy metabolism, amino acid metabolism, lipid metabolism, nucleotide metabolism, and biosynthesis of secondary metabolites. The TCA cycle was the pathway most represented with 5 out of 10 metabolites, 6 out of 8 protein enzymes, and 8 out of 8 enzyme transcripts. On the other hand, gaps, missed pathways, included metabolism of terpenoids and polyketides and lipid metabolism. The multi-omics resource generated in this work will set the basis for ongoing and future studies, bringing the Holm oak closer to model species, to obtain a better understanding of the molecular mechanisms underlying phenotypes of interest (productive, tolerant to environmental cues, nutraceutical value) and to select elite genotypes to be used in restoration and reforestation programs, especially in a future climate change scenario.