Single cell atlas of Xenoturbella bocki highlights limited cell-type complexity

Phylogenetic analyses over the last two decades have united a few small, and previously orphan clades, the nematodermatids, acoels and xenoturbelids, into the phylum Xenacoelomorpha. Some phylogenetic analyses support a sister relationship between Xenacoelomorpha and Ambulacraria (Xenambulacraria), while others suggest that Xenacoelomorpha may be sister to the rest of the Bilateria (Nephrozoa). An understanding of the cell type complements of Xenacoelomorphs is essential to assessing these alternatives as well as to our broader understanding of bilaterian cell type evolution. Employing whole organism single-cell RNA-seq in the marine xenacoelomorph worm Xenoturbella bocki, we show that Xenambulacrarian nerve nets share regulatory features and a peptidergic identity with those found in cnidarians and protostomes and more broadly share muscle and gland cell similarities with other metazoans. Taken together, these data are consistent with broad homologies of animal gland, muscle, and neurons as well as more specific affinities between Xenoturbella and acoel gut and epidermal tissues, consistent with the monophyly of Xenacoelomorpha.

Temporospatial hierarchy and allele-specific expression of zygotic genome activation revealed by distant interspecific urochordate hybrids

Zygotic genome activation (ZGA) is a universal process in early embryogenesis of metazoan, when the quiescent zygotic nucleus initiates global transcription. However, the mechanisms related to massive genome activation and allele-specific expression (ASE) remain not well understood. Here, we develop hybrids from two deeply diverged (120 Mya) ascidian species to symmetrically document the dynamics of ZGA. We identify two coordinated ZGA waves represent early developmental and housekeeping gene reactivation, respectively. Single-cell RNA sequencing reveals that the major expression wave exhibits spatial heterogeneity and significantly correlates with cell fate. Moreover, allele-specific expression occurs in a species- rather than parent-related manner, demonstrating the divergence of cis-regulatory elements between the two species. These findings provide insights into ZGA in chordates.

Key homeobox transcription factors regulate the development of the firefly's adult light organ and bioluminescence

Adult fireflies exhibit unique flashing courtship signals, emitted by specialized light organs, which develop mostly independently from larval light organs during the pupal stage. The mechanisms of adult light organ development have not been thoroughly studied until now. Here we show that key homeobox transcription factors AlABD-B and AlUNC-4 regulate the development of adult light organs and bioluminescence in the firefly Aquatica leii. Interference with the expression of AlAbd-B and AlUnc-4 genes results in undeveloped or non-luminescent adult light organs. AlABD-B regulates AlUnc-4, and they interact with each other. AlABD-B and AlUNC-4 activate the expression of the luciferase gene AlLuc1 and some peroxins. Four peroxins are involved in the import of AlLUC1 into peroxisomes. Our study provides key insights into the development of adult light organs and flash signal control in fireflies.

Beyond the model organism-Mechanistic analysis of protein post-translational modifications is ready for all challengers

A historic challenge for shotgun proteomics has been the requirement for high quality, simple, and nonredundant curated protein sequences in small .fasta text files. Due to the intrinsic informatic challenges and time required to assemble these files, proteomics has struggled to expand beyond the confines of a few model organisms. When considering post-translational modifications that may or may not be present on a specific peptide sequence, these factors inevitably compound. A study on how mangos continue to ripen on the shelf may not be the first thing you'd think of as proof of a scientific discipline shedding historic limitations. However, Bautiste-Valle et al., may be just that. These authors present a quantitative comparison of both peptide and glycopeptide alterations through the complexity of the fruit ripening process and in this we see the present state of a field that no longer needs to wait on genomics to obtain deep mechanistic insights.

Engineering the Soil Bacterium Pseudomonas synxantha 2-79 into a Ratiometric Bioreporter for Phosphorus Limitation

Microbial bioreporters hold promise for addressing challenges in medical and environmental applications. However, the difficulty in ensuring their stable persistence and function within the target environment remains a challenge. One strategy is to integrate information about the host strain and target environment into the design-build-test cycle of the bioreporter itself. Here, we present a case study for such an environmentally motivated design process by engineering the wheat commensal bacterium Pseudomonas synxantha 2-79 into a ratiometric bioreporter for phosphorus limitation. Comparative analysis showed that an exogenous P-responsive promoter outperformed its native counterparts. This reporter can selectively sense and report phosphorus limitation at plant-relevant concentrations of 25-100 μM without cross-activation from carbon or nitrogen limitation or high cell densities. Its performance is robust over a field-relevant pH range (5.8-8), and it responds only to inorganic phosphorus, even in the presence of common soil organic P. Finally, we used fluorescein-calibrated flow cytometry to assess whether the reporter's performance in shaken liquid culture predicts its performance in soil, finding that although the reporter is still functional at the bulk level, its variability in performance increases when grown in a soil slurry as compared to planktonic culture, with a fraction of the population not expressing the reporter proteins. Together, our environmentally aware design process provides an example of how laboratory bioengineering efforts can generate microbes with a greater promise to function reliably in their applied contexts.

Estimating admixture pedigrees of recent hybrids without a contiguous reference genome

The genome of recently admixed individuals or hybrids has characteristic genetic patterns that can be used to learn about their recent admixture history. One of these are patterns of interancestry heterozygosity, which can be inferred from SNP data from either called genotypes or genotype likelihoods, without the need for information on genomic location. This makes them applicable to a wide range of data that are often used in evolutionary and conservation genomic studies, such as low-depth sequencing mapped to scaffolds and reduced representation sequencing. Here we implement maximum likelihood estimation of interancestry heterozygosity patterns using two complementary models. We furthermore develop apoh (Admixture Pedigrees of Hybrids), a software that uses estimates of paired ancestry proportions to detect recently admixed individuals or hybrids, and to suggest possible admixture pedigrees. It furthermore calculates several hybrid indices that make it easier to identify and rank possible admixture pedigrees that could give rise to the estimated patterns. We implemented apoh both as a command line tool and as a Graphical User Interface that allows the user to automatically and interactively explore, rank and visualize compatible recent admixture pedigrees, and calculate the different summary indices. We validate the performance of the method using admixed family trios from the 1000 Genomes Project. In addition, we show its applicability on identifying recent hybrids from RAD-seq data of Grant's gazelle (Nanger granti and Nanger petersii) and whole genome low-depth data of waterbuck (Kobus ellipsiprymnus) which shows complex admixture of up to four populations.

The evolution of antimicrobial peptides in Chiroptera

High viral tolerance coupled with an extraordinary regulation of the immune response makes bats a great model to study host-pathogen evolution. Although many immune-related gene gains and losses have been previously reported in bats, important gene families such as antimicrobial peptides (AMPs) remain understudied. We built an exhaustive bioinformatic pipeline targeting the major gene families of defensins and cathelicidins to explore AMP diversity and analyze their evolution and distribution across six bat families. A combination of manual and automated procedures identified 29 AMP families across queried species, with α-, β-defensins, and cathelicidins representing around 10% of AMP diversity. Gene duplications were inferred in both α-defensins, which were absent in five species, and three β-defensin gene subfamilies, but cathelicidins did not show significant shifts in gene family size and were absent in Anoura caudifer and the pteropodids. Based on lineage-specific gains and losses, we propose diet and diet-related microbiome evolution may determine the evolution of α- and β-defensins gene families and subfamilies. These results highlight the importance of building species-specific libraries for genome annotation in non-model organisms and shed light on possible drivers responsible for the rapid evolution of AMPs. By focusing on these understudied defenses, we provide a robust framework for explaining bat responses to pathogens.

Bat teeth illuminate the diversification of mammalian tooth classes

Tooth classes are an innovation that has contributed to the evolutionary success of mammals. However, our understanding of the mechanisms by which tooth classes diversified remain limited. We use the evolutionary radiation of noctilionoid bats to show how the tooth developmental program evolved during the adaptation to new diet types. Combining morphological, developmental and mathematical modeling approaches, we demonstrate that tooth classes develop through independent developmental cascades that deviate from classical models. We show that the diversification of tooth number and size is driven by jaw growth rate modulation, explaining the rapid gain/loss of teeth in this clade. Finally, we mathematically model the successive appearance of tooth buds, supporting the hypothesis that growth acts as a key driver of the evolution of tooth number and size. Our work reveal how growth, by tinkering with reaction/diffusion processes, drives the diversification of tooth classes and other repeated structure during adaptive radiations.

Comprehensive RNA-Seq Analysis Pipeline for Non-Model Organisms and Its Application in Schmidtea mediterranea

RNA sequencing (RNA-seq) is a high-throughput technology that provides in-depth information on transcriptome. The advancement and dropping costs of RNA sequencing, accompanied by more available reference genomes for different species, make transcriptome analysis in non-model organisms possible. Current obstacles in analyzing RNA-seq data include a lack of functional annotation, which may complicate the process of linking genes to corresponding functions. Here, we provide a one-stop RNA-seq analysis pipeline, PipeOne-NM, for transcriptome functional annotation, non-coding RNA identification, and transcripts alternative splicing analysis of non-model organisms, intended for use with Illumina platform-based RNA-seq data. We performed PipeOne-NM on 237 Schmidtea mediterranea RNA-seq runs and assembled a transcriptome with 84,827 sequences from 49,320 genes, identifying 64,582 mRNA from 35,485 genes, 20,217 lncRNA from 17,084 genes, and 3481 circRNAs from 1103 genes. In addition, we performed a co-expression analysis of lncRNA and mRNA and identified that 1319 lncRNA co-express with at least one mRNA. Further analysis of samples from S. mediterranea sexual and asexual strains revealed the role of sexual reproduction in gene expression profiles. Samples from different parts of asexual S. mediterranea revealed that differential expression profiles of different body parts correlated with the function of conduction of nerve impulses. In conclusion, PipeOne-NM has the potential to provide comprehensive transcriptome information for non-model organisms on a single platform.

Epigenetics across the evolutionary tree: New paradigms from non-model animals

All animals have evolved solutions to manage their genomes, enabling the efficient organization of meters of DNA strands in the nucleus and allowing for nuanced regulation of gene expression while keeping transposable elements suppressed. Epigenetic modifications are central to accomplishing all these. Recent advances in sequencing technologies and the development of techniques that profile epigenetic marks and chromatin accessibility using reagents that can be used in any species has catapulted epigenomic studies in diverse animal species, shedding light on the multitude of epigenomic mechanisms utilized across the evolutionary tree. Now, comparative epigenomics is a rapidly growing field that is uncovering mechanistic aspects of epigenetic modifications and chromatin organization in non-model invertebrates, ranging from octopus to sponges. This review puts recent discoveries in the epigenetics of non-model invertebrates in historical context, and describes new insight into the patterning and functions of DNA methylation and other highly conserved epigenetic modifications.

gga-miRNOME, a microRNA-sequencing dataset from chick embryonic tissues

MicroRNAs (miRNAs) are small non-coding RNA molecules, with sizes ranging from 18 to 25 nucleotides, which are key players in gene expression regulation. These molecules play an important role in fine-tuning early vertebrate embryo development. However, there are scarce publicly available miRNA datasets from non-mammal embryos, such as the chicken (Gallus gallus), which is a classical model system to study vertebrate embryogenesis. Here, we performed microRNA-sequencing to characterize the early stages of trunk and limb development in the chick embryo. For this, we profiled three chick embryonic tissues, namely, Undetermined Presomitic Mesoderm (PSM_U), Determined Presomitic Mesoderm (PSM_D) and Forelimb Distal Cyclic Domain (DCD). We identified 926 known miRNAs, and 1,141 novel candidate miRNAs, which nearly duplicates the number of Gallus gallus entries in the miRBase database. These data will greatly benefit the avian research community, particularly by highlighting new miRNAs potentially involved in the regulation of early vertebrate embryo development, that can be prioritized for further experimental testing.

Diversifying Isoprenoid Platforms via Atypical Carbon Substrates and Non-model Microorganisms

Isoprenoid compounds are biologically ubiquitous, and their characteristic modularity has afforded products ranging from pharmaceuticals to biofuels. Isoprenoid production has been largely successful in Escherichia coli and Saccharomyces cerevisiae with metabolic engineering of the mevalonate (MVA) and methylerythritol phosphate (MEP) pathways coupled with the expression of heterologous terpene synthases. Yet conventional microbial chassis pose several major obstacles to successful commercialization including the affordability of sugar substrates at scale, precursor flux limitations, and intermediate feedback-inhibition. Now, recent studies have challenged typical isoprenoid paradigms by expanding the boundaries of terpene biosynthesis and using non-model organisms including those capable of metabolizing atypical C1 substrates. Conversely, investigations of non-model organisms have historically informed optimization in conventional microbes by tuning heterologous gene expression. Here, we review advances in isoprenoid biosynthesis with specific focus on the synergy between model and non-model organisms that may elevate the commercial viability of isoprenoid platforms by addressing the dichotomy between high titer production and inexpensive substrates.

Studying the genetic basis of masting

The mechanisms underlying mast seeding have traditionally been studied by collecting long-term observational data on seed crops and correlating seedfall with environmental variables. Significant progress in ecological genomics will improve our understanding of the evolution of masting by clarifying the genetic basis of masting traits and the role of natural selection in shaping those traits. Here, we summarize three important aspects in studying the evolution of masting at the genetic level: which traits govern masting, whether those traits are genetically regulated, and which taxa show wide variation in these traits. We then introduce recent studies on the molecular mechanisms of masting. Those studies measure seasonal changes in gene expression in natural conditions to quantify how multiple environmental factors combine to regulate floral initiation, which in many masting plant species is the single largest contributor to among-year variation in seed crops. We show that Fagaceae offers exceptional opportunities for evolutionary investigations because of its diversity at both the phenotypic and genetic levels and existing documented genome sequences. This article is part of the theme issue 'The ecology and evolution of synchronized seed production in plants'.

Characterization of germline development and identification of genes associated with germline specification in pineapple

Understanding germline specification in plants could be advantageous for agricultural applications. In recent decades, substantial efforts have been made to understand germline specification in several plant species, including Arabidopsis, rice, and maize. However, our knowledge of germline specification in many agronomically important plant species remains obscure. Here, we characterized the female germline specification and subsequent female gametophyte development in pineapple using callose staining, cytological, and whole-mount immunolocalization analyses. We also determined the male germline specification and gametophyte developmental timeline and observed male meiotic behavior using chromosome spreading assays. Furthermore, we identified 229 genes that are preferentially expressed at the megaspore mother cell (MMC) stage during ovule development and 478 genes that are preferentially expressed at the pollen mother cell (PMC) stage of anther development using comparative transcriptomic analysis. The biological functions, associated regulatory pathways and expression patterns of these genes were also analyzed. Our study provides a convenient cytological reference for exploring pineapple germline development and a molecular basis for the future functional analysis of germline specification in related plant species.

Integration of transcriptomic and proteomic approaches for snake venom profiling

INTRODUCTION

Snake venoms contain many protein and peptide isoforms with high levels of sequence variation, even within a single species.

AREAS COVERED

In this review, we highlight several examples, from both published and unpublished work in our lab, demonstrating how a combined venom gland transcriptome and proteome methodology allows for comprehensive characterization of venoms, including those from understudied rear-fanged snake species, and we provide recommendations for using these approaches.

EXPERT OPINION

When characterizing venoms, peptide mass fingerprinting using databases built predominately from protein sequences originating from model organisms can be disadvantageous, especially when the intention is to document protein diversity. Therefore, the use of species-specific venom gland transcriptomes corrects for the absence of these unique peptide sequences in databases. The integration of transcriptomics and proteomics improves the accuracy of either approach alone for venom profiling.

Morphological studies of rose prickles provide new insights

Prickles are common structures in plants that play a key role in defense against herbivores. In the Rosa genus, prickles are widely present with great diversity in terms of form and density. For cut rose production, prickles represent an important issue, as they can damage the flower and injure workers. Our objectives were to precisely describe the types of prickles that exist in roses, their tissues of origin and their development. We performed a detailed histological analysis of prickle initiation and development in a rose F1 population. Based on the prickle investigation of 110 roses, we proposed the first categorization of prickles in the Rosa genus. They are mainly divided into two categories, nonglandular prickles (NGPs) and glandular prickles (GPs), and subcategories were defined based on the presence/absence of hairs and branches. We demonstrated that NGPs and GPs both originate from multiple cells of the ground meristem beneath the protoderm. For GPs, the gland cells originate from the protoderm of the GP at the early developmental stage. Our findings clearly demonstrate that prickles are not modified trichomes (which originate from the protoderm). These conclusions are different from the current mainstream hypothesis. These results provide a foundation for further studies on prickle initiation and development in plants.

Development of the first axillary in vitro shoot multiplication protocol for coconut palms

The coconut palm or "tree of life" is one of nature's most useful plants and the demand for its fruit is increasing. However, coconut production is threatened by ageing plantations, pests and diseases. Currently, the palm is exclusively propagated via seeds, limiting the amount of planting material. A novel micropropagation method is presented, based on axillary shoot formation. Apical meristems of in vitro coconut seedlings are cultured onto Y3 medium containing 1 µM TDZ. This induces the apical meristem to proliferate through axillary shoots in ~ 27% of the initiated explants. These axillary shoots are seen as white clumps of proliferating tissue and can be multiplied at a large scale or regenerated into rooted in vitro plantlets. This innovative micropropagation method will enable the production of disease-free, high quality in vitro plantlets, which will solve the worldwide scarcity of coconut planting material.

Morphological, phenological, and transcriptional analyses provide insight into the diverse flowering traits of a mutant of the relic woody plant Liriodendron chinense

Flowering is crucial to plant reproduction and controlled by multiple factors. However, the mechanisms underlying the regulation of flowering in perennial plants are still largely unknown. Here, we first report a super long blooming 1 (slb1) mutant of the relict tree Liriodendron chinense possessing a prolonged blooming period of more than 5 months, in contrast to the 1 month blooming period in the wild type (WT). Phenotypic characterization showed that earlier maturation of lateral shoots was caused by accelerated axillary bud fate, leading to the phenotype of continuous flowering in slb1 mutants. The transcriptional activity of genes related to hormone signaling (auxin, cytokinin, and strigolactone), nutrient availability, and oxidative stress relief further indicated active outgrowth of lateral buds in slb1 mutants. Interestingly, we discovered a unique FT splicing variant with intron retention specific to slb1 mutants, representing a potential causal mutation in the slb1 mutants. Surprisingly, most slb1 inbred offspring flowered precociously with shorter juvenility (~4 months) than that (usually 8-10 years) required in WT plants, indicating heritable variation underlying continuous flowering in slb1 mutants. This study reports an example of a perennial tree mutant that flowers continuously, providing a rare resource for both breeding and genetic research.

Sugars promote graft union development in the heterograft of cucumber onto pumpkin

The use of heterografts is widely applied for the production of several important commercial crops, but the molecular mechanism of graft union formation remains poorly understood. Here, cucumber grafted onto pumpkin was used to study graft union development, and genome-wide tempo-spatial gene expression at the graft interface was comprehensively investigated. Histological analysis suggested that resumption of the rootstock growth occurred after both phloem and xylem reconnection, and the scion showed evident callus production compared with the rootstock 3 days after grafting. Consistently, transcriptome data revealed specific responses between the scion and rootstock in the expression of genes related to cambium development, the cell cycle, and sugar metabolism during both vascular reconnection and healing, indicating distinct mechanisms. Additionally, lower levels of sugars and significantly changed sugar enzyme activities at the graft junction were observed during vascular reconnection. Next, we found that the healing process of grafted etiolated seedlings was significantly delayed, and graft success, xylem reconnection, and the growth of grafted plants were enhanced by exogenous glucose. This demonstrates that graft union formation requires the correct sugar content. Furthermore, we also found that graft union formation was delayed with a lower energy charge by the target of rapamycin (TOR) inhibitor AZD-8055, and xylem reconnection and the growth of grafted plants were enhanced under AZD-8055 with exogenous glucose treatment. Taken together, our results reveal that sugars play a positive role in graft union formation by promoting the growth of cucumber/pumpkin and provide useful information for understanding graft union healing and the application of heterografting in the future.

Single-cell RNA-seq analysis reveals penaeid shrimp hemocyte subpopulations and cell differentiation process

Crustacean aquaculture is expected to be a major source of fishery commodities in the near future. Hemocytes are key players of the immune system in shrimps; however, their classification, maturation, and differentiation are still under debate. To date, only discrete and inconsistent information on the classification of shrimp hemocytes has been reported, showing that the morphological characteristics are not sufficient to resolve their actual roles. Our present study using single-cell RNA sequencing revealed six types of hemocytes of Marsupenaeus japonicus based on their transcriptional profiles. We identified markers of each subpopulation and predicted the differentiation pathways involved in their maturation. We also predicted cell growth factors that might play crucial roles in hemocyte differentiation. Different immune roles among these subpopulations were suggested from the analysis of differentially expressed immune-related genes. These results provide a unified classification of shrimp hemocytes, which improves the understanding of its immune system.

Comparative analysis of phenotypic plasticity sheds light on the evolution and molecular underpinnings of locust phase polyphenism

Locusts exhibit one of nature's most spectacular examples of complex phenotypic plasticity, in which changes in density cause solitary and cryptic individuals to transform into gregarious and conspicuous locusts forming large migrating swarms. We investigated how these coordinated alternative phenotypes might have evolved by studying the Central American locust and three closely related non-swarming grasshoppers in a comparative framework. By experimentally isolating and crowding during nymphal development, we induced density-dependent phenotypic plasticity and quantified the resulting behavioural, morphological, and molecular reaction norms. All four species exhibited clear plasticity, but the individual reaction norms varied among species and showed different magnitudes. Transcriptomic responses were species-specific, but density-responsive genes were functionally similar across species. There were modules of co-expressed genes that were highly correlated with plastic reaction norms, revealing a potential molecular basis of density-dependent phenotypic plasticity. These findings collectively highlight the importance of studying multiple reaction norms from a comparative perspective.

Mutation of SlARC6 leads to tissue-specific defects in chloroplast development in tomato

The proliferation and development of chloroplasts are important for maintaining the normal chloroplast population in plant tissues. Most studies have focused on chloroplast maintenance in leaves. In this study, we identified a spontaneous mutation in a tomato mutant named suffulta (su), in which the stems appeared albinic while the leaves remained normal. Map-based cloning showed that Su encodes a DnaJ heat shock protein that is a homolog of the Arabidopsis gene AtARC6, which is involved in chloroplast division. Knockdown and knockout of SlARC6 in wild-type tomato inhibit chloroplast division, indicating the conserved function of SlARC6. In su mutants, most mesophyll cells contain only one or two giant chloroplasts, while no chloroplasts are visible in 60% of stem cells, resulting in the albinic phenotype. Compared with mature tissues, the meristem of su mutants suggested that chloroplasts could partially divide in meristematic cells, suggesting the existence of an alternative mechanism in those dividing cells. Interestingly, the adaxial petiole cells of su mutants contain more chloroplasts than the abaxial cells. In addition, prolonged lighting can partially rescue the albinic phenotypes in su mutants, implying that light may promote SlACR6-independent chloroplast development. Our results verify the role of SlACR6 in chloroplast division in tomato and uncover the tissue-specific regulation of chloroplast development.

CsTCPs regulate shoot tip development and catechin biosynthesis in tea plant (Camellia sinensis)

The growth of leaves and biosynthesis of characteristic secondary metabolites are critically important for tea production and quality control. However, little is known about the coordinated regulation of leaf development and catechin biosynthesis in tea plants. Here, we reported that TCP TFs are involved in both catechin biosynthesis and leaf development. An integrated analysis of catechin profiling and CsTCP expression in different tissues of plants under various environmental conditions at different developmental stages indicated significant correlations between the transcript levels of CIN-type TCPs and catechin production. CIN-type CsTCP3 and CsTCP4 and PCF-type CsTCP14 interacted with the MYB-bHLH-WD40 repeat (MBW) complex by forming a CsTCP3-CsTT8 heterodimer and modulating the transactivation activity of the promoters of anthocyanin synthase (CsANS1) and anthocyanidin reductase (CsANR1). Four types of microRNA/target modules, miR319b/CsTCP3-4, miR164b/CsCUC, miR396/CsGRF-GIF, and miR165b/HD-ZIPIII ones, were also identified and characterized for their functions in the regulation of the development of tea plant shoot tips and leaf shape. The results of these modules were reflected by their different expression patterns in developing buds and leaves that had distinctly different morphologies in three different tea plant varieties. Their roles in the regulation of catechin biosynthesis were also further verified by manipulation of microRNA319b (miR319b), which targets the transcripts of CsTCP3 and CsTCP4. Thus, CsTCPs represent at least one of these important groups of TFs that can integrate tea plant leaf development together with secondary metabolite biosynthesis. Our study provides new insight into shoot tip development and catechin production in tea plants and lays a foundation for further mechanistic understanding of the regulation of tea plant leaf development and secondary metabolism.

Sensing the world and its dangers: An evolutionary perspective in neuroimmunology

Detecting danger is key to the survival and success of all species. Animal nervous and immune systems cooperate to optimize danger detection. Preceding studies have highlighted the benefits of bringing neurons into the defense game, including regulation of immune responses, wound healing, pathogen control, and survival. Here, we summarize the body of knowledge in neuroimmune communication and assert that neuronal participation in the immune response is deeply beneficial in each step of combating infection, from inception to resolution. Despite the documented tight association between the immune and nervous systems in mammals or invertebrate model organisms, interdependence of these two systems is largely unexplored across metazoans. This review brings a phylogenetic perspective of the nervous and immune systems in the context of danger detection and advocates for the use of non-model organisms to diversify the field of neuroimmunology. We identify key taxa that are ripe for investigation due to the emergence of key evolutionary innovations in their immune and nervous systems. This novel perspective will help define the primordial principles that govern neuroimmune communication across taxa.

Identification and Validation of Reference Genes in Clostridium beijerinckii NRRL B-598 for RT-qPCR Using RNA-Seq Data

Gene expression analysis through reverse transcription-quantitative real-time polymerase chain reaction (RT-qPCR) depends on correct data normalization by reference genes with stable expression. Although Clostridium beijerinckii NRRL B-598 is a promising Gram-positive bacterium for the industrial production of biobutanol, validated reference genes have not yet been reported. In this study, we selected 160 genes with stable expression based on an RNA sequencing (RNA-Seq) data analysis, and among them, seven genes (zmp, rpoB1, rsmB, greA, rpoB2, topB2, and rimO) were selected for experimental validation by RT-qPCR and gene ontology (GO) enrichment analysis. According to statistical analyses, zmp and greA were the most stable and suitable reference genes for RT-qPCR normalization. Furthermore, our methodology can be useful for selection of the reference genes in other strains of C. beijerinckii and it also suggests that the RNA-Seq data can be used for the initial selection of novel reference genes, however, their validation is required.

PbEIL1 acts upstream of PbCysp1 to regulate ovule senescence in seedless pear

Numerous environmental and endogenous signals control the highly orchestrated and intricate process of plant senescence. Ethylene, a well-known inducer of senescence, has long been considered a key endogenous regulator of leaf and flower senescence, but the molecular mechanism of ethylene-induced ovule senescence has not yet been elucidated. In this study, we found that blockage of fertilization caused ovule abortion in the pear cultivar '1913'. According to transcriptome and phytohormone content data, ethylene biosynthesis was activated by pollination. At the same time, ethylene overaccumulated in ovules, where cells were sensitive to ethylene signals in the absence of fertilization. We identified a transcription factor in the ethylene signal response, ethylene-insensitive 3-like (EIL1), as a likely participant in ovule senescence. Overexpression of PbEIL1 in tomato caused precocious onset of ovule senescence. We further found that EIL1 could directly bind to the promoter of the SENESCENCE-ASSOCIATED CYSTEINE PROTEINASE 1 (PbCysp1) gene and act upstream of senescence. Yeast one-hybrid and dual-luciferase assays revealed the interaction of the transcription factor and the promoter DNA sequence and demonstrated that PbEIL1 enhanced the action of PbCysp1. Collectively, our results provide new insights into how ethylene promotes the progression of unfertilized ovule senescence.

Biosynthetic labeling with 3-O-propargylcaffeyl alcohol reveals in vivo cell-specific patterned lignification in loquat fruits during development and postharvest storage

Lignification is a major cell wall modification that often results in the formation of sophisticated subcellular patterns during plant development or in response to environmental stresses. Precise localization of the spatiotemporal deposition of lignin is of great importance for revealing the lignification regulatory mechanism of individual cells. In loquat fruits, lignification typically increases the flesh lignin content and firmness, reducing their edibility and processing quality. However, the precise localization of the spatiotemporal active zones of lignification inside loquat fruit flesh remains poorly understood, and little is known about the contribution of patterned lignification to cell wall structure dynamics and the subsequent fruit-quality deterioration. Here, we performed an emerging bioorthogonal chemistry imaging technique to trace the in vivo patterned lignification dynamics in cells of loquat fruit flesh during development and storage. In developing fruits, lignified cells (LCs) and vascular bundles (VBs) were the zones of active lignification, and ring-like LCs deposited lignin at both the inner wall layer and the outer periphery sides. The domino effect of the generation of LCs was preliminarily visualized. In mature fruits, the newly formed lignin in the flesh of fruits during storage was specifically deposited in the corners and middle lamellae of parenchyma cells surrounding the VBs, resulting in the development of a reticular structure. Based on the findings, distinct spatiotemporal patterned lignification modes for different flesh cells in loquat fruits were proposed. These findings provide loquat lignification dynamics together with spatiotemporal data that can improve our understanding of the lignification process in planta.

Reannotation of the cultivated strawberry genome and establishment of a strawberry genome database

Cultivated strawberry (Fragaria × ananassa) is an important fruit crop species whose fruits are enjoyed by many worldwide. An octoploid of hybrid origin, the complex genome of this species was recently sequenced, serving as a key reference genome for cultivated strawberry and related species of the Rosaceae family. The current annotation of the F. ananassa genome mainly relies on ab initio predictions and, to a lesser extent, transcriptome data. Here, we present the structure and functional reannotation of the F. ananassa genome based on one PacBio full-length RNA library and ninety-two Illumina RNA-Seq libraries. This improved annotation of the F. ananassa genome, v1.0.a2, comprises a total of 108,447 gene models, with 97.85% complete BUSCOs. The models of 19,174 genes were modified, 360 new genes were identified, and 11,044 genes were found to have alternatively spliced isoforms. Additionally, we constructed a strawberry genome database (SGD) for strawberry gene homolog searching and annotation downloading. Finally, the transcriptome of the receptacles and achenes of F. ananassa at four developmental stages were reanalyzed and qualified, and the expression profiles of all the genes in this annotation are also provided. Together, this study provides an updated annotation of the F. ananassa genome, which will facilitate genomic analyses across the Rosaceae family and gene functional studies in cultivated strawberry.

Hydrogen sulfide promotes flowering in heading Chinese cabbage by S-sulfhydration of BraFLCs

Heading Chinese cabbage (Brassica rapa L. syn. B. campestris L. ssp. chinensis Makino var. pekinensis (Rupr.) J. Cao et Sh. Cao) is a cruciferous Brassica vegetable that has a triplicate genome, owing to an ancient genome duplication event. It is unclear whether the duplicated homologs have conserved or diversified functions. Hydrogen sulfide (H₂S) is a plant gasotransmitter that plays important physiological roles in growth, development, and responses to environmental stresses. The modification of cysteines through S-sulfhydration is an important mechanism of H₂S, which regulates protein functions. H₂S promotes flowering in Arabidopsis and heading Chinese cabbage. Here we investigated the molecular mechanisms of H₂S used to promote flowering in the latter. Four, five, and four BraFLC, BraSOC I, and BraFT homologs were identified in heading Chinese cabbage. Different BraFLC proteins were bound to different CArG boxes in the promoter regions of the BraSOC I and BraFT homologs, producing different binding patterns. Thus, there may be functionally diverse BraFLC homologs in heading Chinese cabbage. Exogenous H₂S at 100 μmol L^-1 significantly promoted flowering by compensating for insufficient vernalization. BraFLC 1 and BraFLC 3 underwent S-sulfhydration by H₂S, after which their abilities to bind most BraSOC I or BraFT promoter probes weakened or even disappeared. These changes in binding ability were consistent with the expression pattern of the BraFT and BraSOC I homologs in seedlings treated with H₂S. These results indicated that H₂S signaling regulates flowering time. In summary, H₂S signaling promoted plant flowering by weakening or eliminating the binding abilities of BraFLCs to downstream promoters through S-sulfhydration.

Melon short internode (CmSi) encodes an ERECTA-like receptor kinase regulating stem elongation through auxin signaling

Plant height is one of the most important agronomic traits that directly determines plant architecture, and compact or dwarf plants can allow for increased planting density and land utilization as well as increased lodging resistance and economic yield. At least four dwarf/semidwarf genes have been identified in different melon varieties, but none of them have been cloned, and little is known about the molecular mechanisms underlying internode elongation in melon. Here, we report map-based cloning and functional characterization of the first semidwarf gene short internode (Cmsi) in melon, which encodes an ERECTA-like receptor kinase regulating internode elongation. Spatial-temporal expression analyses revealed that CmSI exhibited high expression in the vascular bundle of the main stem during internode elongation. The expression level of CmSI was positively correlated with stem length in the different melon varieties examined. Ectopic expression of CmSI in Arabidopsis and cucumber suggested CmSI as a positive regulator of internode elongation in both species. Phytohormone quantitation and transcriptome analysis showed that the auxin content and the expression levels of a number of genes involved in the auxin signaling pathway were altered in the semidwarf mutant, including several well-known auxin transporters, such as members of the ABCB family and PIN-FORMED genes. A melon polar auxin transport protein CmPIN2 was identified by protein-protein interaction assay as physically interacting with CmSI to modulate auxin signaling. Thus, CmSI functions in an auxin-dependent regulatory pathway to control internode elongation in melon. Our findings revealed that the ERECTA family gene CmSI regulates stem elongation in melon through auxin signaling, which can directly affect polar auxin transport.

Genomic and transcriptomic variation defines the chromosome-scale assembly of Haemonchus contortus, a model gastrointestinal worm

Haemonchus contortus is a globally distributed and economically important gastrointestinal pathogen of small ruminants and has become a key nematode model for studying anthelmintic resistance and other parasite-specific traits among a wider group of parasites including major human pathogens. Here, we report using PacBio long-read and OpGen and 10X Genomics long-molecule methods to generate a highly contiguous 283.4 Mbp chromosome-scale genome assembly including a resolved sex chromosome for the MHco3(ISE).N1 isolate. We show a remarkable pattern of conservation of chromosome content with Caenorhabditis elegans, but almost no conservation of gene order. Short and long-read transcriptome sequencing allowed us to define coordinated transcriptional regulation throughout the parasite's life cycle and refine our understanding of cis- and trans-splicing. Finally, we provide a comprehensive picture of chromosome-wide genetic diversity both within a single isolate and globally. These data provide a high-quality comparison for understanding the evolution and genomics of Caenorhabditis and other nematodes and extend the experimental tractability of this model parasitic nematode in understanding helminth biology, drug discovery and vaccine development, as well as important adaptive traits such as drug resistance.

WUSCHEL-related homeobox1 (WOX1) regulates vein patterning and leaf size in Cucumis sativus

In plants, WUSCHEL-related homeobox1 (WOX1) homologs promote lamina mediolateral outgrowth. However, the downstream components linking WOX1 and lamina development remain unclear. In this study, we revealed the roles of WOX1 in palmate leaf expansion in cucumber (Cucumis sativus). A cucumber mango fruit (mf) mutant, resulting from truncation of a WOX1-type protein (CsWOX1), displayed abnormal lamina growth and defects in the development of secondary and smaller veins. CsWOX1 was expressed in the middle mesophyll and leaf margins and rescued defects of the Arabidopsis wox1 prs double mutant. Transcriptomic analysis revealed that genes involved in auxin polar transport and auxin response were highly associated with leaf development. Analysis of the cucumber mf rl (round leaf) double mutant revealed that CsWOX1 functioned in vein development via PINOID (CsPID1)-controlled auxin transport. Overexpression of CsWOX1 in cucumber (CsWOX1-OE) affected vein patterning and produced 'butterfly-shaped' leaves. CsWOX1 physically interacted with CsTCP4a, which may account for the abnormal lamina development in the mf mutant line and the smaller leaves in the CsWOX1-OE plants. Our findings demonstrated that CsWOX1 regulates cucumber leaf vein development by modulating auxin polar transport; moreover, CsWOX1 regulates leaf size by controlling CIN-TCP genes.

Using automated reasoning to explore the metabolism of unconventional organisms: a first step to explore host-microbial interactions

Systems modelled in the context of molecular and cellular biology are difficult to represent with a single calibrated numerical model. Flux optimisation hypotheses have shown tremendous promise to accurately predict bacterial metabolism but they require a precise understanding of metabolic reactions occurring in the considered species. Unfortunately, this information may not be available for more complex organisms or non-cultured microorganisms such as those evidenced in microbiomes with metagenomic techniques. In both cases, flux optimisation techniques may not be applicable to elucidate systems functioning. In this context, we describe how automatic reasoning allows relevant features of an unconventional biological system to be identified despite a lack of data. A particular focus is put on the use of Answer Set Programming, a logic programming paradigm with combinatorial optimisation functionalities. We describe its usage to over-approximate metabolic responses of biological systems and solve gap-filling problems. In this review, we compare steady-states and Boolean abstractions of metabolic models and illustrate their complementarity via applications to the metabolic analysis of macro-algae. Ongoing applications of this formalism explore the emerging field of systems ecology, notably elucidating interactions between a consortium of microbes and a host organism. As the first step in this field, we will illustrate how the reduction in microbiotas according to expected metabolic phenotypes can be addressed with gap-filling problems.

A Cost Reduced Variant of Epi-Genotyping by Sequencing for Studying DNA Methylation in Non-model Organisms

Reference-free reduced representation bisulfite sequencing uses enzymatic digestion for reducing genome complexity and allows detection of markers to study DNA methylation of a high number of individuals in natural populations of non-model organisms. Current methods like epiGBS enquire the use of a higher number of methylated DNA oligos with a significant cost (especially for small labs and first pilot studies). In this paper, we present a modification of this epiGBS protocol that requires the use of only one hemimethylated P2 (common) adapter, which is combined with unmethylated barcoded adapters. The unmethylated cytosines of one chain of the barcoded adapter are replaced by methylated cytosines using nick translation with methylated cytosines in dNTP solution. The basic version of our technique uses only one restriction enzyme, and as a result, genomic fragments are integrated into two orientations with respect to the adapter sequences. Comparing the sequences of two chain orientations makes it possible to reconstruct the original sequence before bisulfite treatment with the help of standard software and newly developed software written in C and described here. We provide a proof of concept via data obtained from almond (Prunus dulcis). Example data and a detailed description of the complete software pipeline starting from the raw reads up until the final differentially methylated cytosines are given in Supplementary Material making this technique accessible to non-expert computer users. The adapter design showed in this paper should allow the use of a two restriction enzyme approach with minor changes in software parameters.

Proteogenomics-Guided Evaluation of RNA-Seq Assembly and Protein Database Construction for Emergent Model Organisms

Proteogenomics is gaining momentum as, today, genomics, transcriptomics, and proteomics can be readily performed on any new species. This approach allows key alterations to molecular pathways to be identified when comparing conditions. For animals and plants, RNA-seq-informed proteomics is the most popular means of interpreting tandem mass spectrometry spectra acquired for species for which the genome has not yet been sequenced. It relies on high-performance de novo RNA-seq assembly and optimized translation strategies. Here, several pre-treatments for Illumina RNA-seq reads before assembly are explored to translate the resulting contigs into useful polypeptide sequences. Experimental transcriptomics and proteomics datasets acquired for individual Gammarus fossarum freshwater crustaceans are used, the most relevant procedure is defined by the ratio of MS/MS spectra assigned to peptide sequences. Removing reads with a mean quality score of less than 17-which represents a single probable nucleotide error on 150-bp reads-prior to assembly, increases the proteomics outcome. The best translation using Transdecoder is achieved with a minimal open reading frame length of 50 amino acids and systematic selection of ORFs longer than 900 nucleotides. Using these parameters, transcriptome assembly and translation informed by proteomics pave the way to further improvements in proteogenomics.

Near-future ocean warming and acidification alter foraging behaviour, locomotion, and metabolic rate in a keystone marine mollusc

Environmentally-induced changes in fitness are mediated by direct effects on physiology and behaviour, which are tightly linked. We investigated how predicted ocean warming (OW) and acidification (OA) affect key ecological behaviours (locomotion speed and foraging success) and metabolic rate of a keystone marine mollusc, the sea hare Stylocheilus striatus, a specialist grazer of the toxic cyanobacterium Lyngbya majuscula. We acclimated sea hares to OW and/or OA across three developmental stages (metamorphic, juvenile, and adult) or as adults only, and compare these to sea hares maintained under current-day conditions. Generally, locomotion speed and time to locate food were reduced ~1.5- to 2-fold when the stressors (OW or OA) were experienced in isolation, but reduced ~3-fold when combined. Decision-making was also severely altered, with correct foraging choice nearly 40% lower under combined stressors. Metabolic rate appeared to acclimate to the stressors in isolation, but was significantly elevated under combined stressors. Overall, sea hares that developed under OW and/or OA exhibited a less severe impact, indicating beneficial phenotypic plasticity. Reduced foraging success coupled with increased metabolic demands may impact fitness in this species and highlight potentially large ecological consequences under unabated OW and OA, namely in regulating toxic cyanobacteria blooms on coral reefs.

Too Many False Targets for MicroRNAs: Challenges and Pitfalls in Prediction of miRNA Targets and Their Gene Ontology in Model and Non-model Organisms

Short ("seed") or extended base pairing between microRNAs (miRNAs) and their target RNAs enables post-transcriptional silencing in many organisms. These interactions allow the computational prediction of potential targets. In model organisms, predicted targets are frequently validated experimentally; hence meaningful miRNA-regulated processes are reported. However, in non-models, these reports mostly rely on computational prediction alone. Many times, further bioinformatic analyses such as Gene Ontology (GO) enrichment are based on these in silico projections. Here such approaches are reviewed, their caveats are highlighted and the ease of picking false targets from predicted lists is demonstrated. Discoveries that shed new light on how miRNAs evolved to regulate targets in various phyletic groups are discussed, in addition to the pitfalls of target identification in non-model organisms. The goal is to prevent the misuse of bioinformatic tools, as they cannot bypass the biological understanding of miRNA-target regulation.

Faraway, so close. The comparative method and the potential of non-model animals in mitochondrial research

Inference from model organisms has been the engine for many discoveries in life science, but indiscriminate generalization leads to oversimplifications and misconceptions. Model organisms and inductive reasoning are irreplaceable: there is no other way to tackle the complexity of living systems. At the same time, it is not advisable to infer general patterns from a restricted number of species, which are very far from being representative of the diversity of life. Not all models are equal. Some organisms are suitable to find similarities across species, other highly specialized organisms can be used to focus on differences. In this opinion piece, we discuss the dominance of the mechanistic/reductionist approach in life sciences and make a case for an enhanced application of the comparative approach to study processes in all their various forms across different organisms. We also enlist some rising animal models in mitochondrial research, to exemplify how non-model organisms can be chosen in a comparative framework. These taxa often do not possess implemented tools and dedicated methods/resources. However, because of specific features, they have the potential to address still unanswered biological questions. Finally, we discuss future perspectives and caveats of the comparative method in the age of 'big data'. This article is part of the theme issue 'Linking the mitochondrial genotype to phenotype: a complex endeavour'.

Impacts of the synthetic androgen Trenbolone on gonad differentiation and development - comparisons between three deeply diverged anuran families

Using a recently developed approach for testing endocrine disruptive chemicals (EDCs) in amphibians, comprising synchronized tadpole exposure plus genetic and histological sexing of metamorphs in a flow-through-system, we tested the effects of 17β-Trenbolone (Tb), a widely used growth promoter in cattle farming, in three deeply diverged anuran families: the amphibian model species Xenopus laevis (Pipidae) and the non-models Bufo(tes) viridis (Bufonidae) and Hyla arborea (Hylidae). Trenbolone was applied in three environmentally and/or physiologically relevant concentrations (0.027 µg/L (10-10 M), 0.27 µg/L (10-9 M), 2.7 µg/L (10-8 M)). In none of the species, Tb caused sex reversals or masculinization of gonads but had negative species-specific impacts on gonad morphology and differentiation after the completion of metamorphosis, independently of genetic sex. In H. arborea and B. viridis, mounting Tb-concentration correlated positively with anatomical abnormalities at 27 µg/L (10-9 M) and 2.7 µg/L (10-8 M), occurring in X. laevis only at the highest Tb concentration. Despite anatomical aberrations, histologically all gonadal tissues differentiated seemingly normally when examined at the histological level but at various rates. Tb-concentration caused various species-specific mortalities (low in Xenopus, uncertain in Bufo). Our data suggest that deep phylogenetic divergence modifies EDC-vulnerability, as previously demonstrated for Bisphenol A (BPA) and Ethinylestradiol (EE2).

Scanning RNA-Seq and RAD-Seq approach to develop SNP markers closely linked to MALE STERILITY 1 (MS1) in Cryptomeria japonica D. Don

Cryptomeria japonica is a major forestry tree species in Japan. Male sterility of the species is caused by a recessive gene, which shows dysfunction of pollen development and results in no dispersed pollen. Because the pollen of C. japonica induces pollinosis, breeding of pollen-free C. japonica is desired. In this study, single nucleotide polymorphism (SNP) markers located at 1.78 and 0.58 cM to a male sterility locus (MS1) were identified from an analysis of RNA-Seq and RAD-Seq, respectively. SNPs closely linked to MS1 were first scanned by a method similar to MutMap, where a type of index was calculated to measure the strength of the linkage between a marker sequence and MS1. Linkage analysis of selected SNP markers confirmed a higher efficiency of the current method to construct a partial map around MS1. Allele-specific PCR primer pair for the most closely linked SNP with MS1 was developed as a codominant marker, and visualization of the PCR products on an agarose gel enabled rapid screening of male sterile C. japonica. The allele-specific primers developed in this study would be useful for establishing the selection of male sterile C. japonica.

Tissue Mechanical Forces and Evolutionary Developmental Changes Act Through Space and Time to Shape Tooth Morphology and Function

Efforts from diverse disciplines, including evolutionary studies and biomechanical experiments, have yielded new insights into the genetic, signaling, and mechanical control of tooth formation and functions. Evidence from fossils and non-model organisms has revealed that a common set of genes underlie tooth-forming potential of epithelia, and changes in signaling environments subsequently result in specialized dentitions, maintenance of dental stem cells, and other phenotypic adaptations. In addition to chemical signaling, tissue forces generated through epithelial contraction, differential growth, and skeletal constraints act in parallel to shape the tooth throughout development. Here recent advances in understanding dental development from these studies are reviewed and important gaps that can be filled through continued application of evolutionary and biomechanical approaches are discussed.

Melanomacrophage Centers As a Histological Indicator of Immune Function in Fish and Other Poikilotherms

Melanomacrophage centers (MMCs) are aggregates of highly pigmented phagocytes found primarily in the head kidney and spleen, and occasionally the liver of many vertebrates. Preliminary histological analyses suggested that MMCs are structurally similar to the mammalian germinal center (GC), leading to the hypothesis that the MMC plays a role in the humoral adaptive immune response. For this reason, MMCs are frequently described in the literature as “primitive GCs” or the “evolutionary precursors” to the mammalian GC. However, we argue that this designation may be premature, having been pieced together from mainly descriptive studies in numerous distinct species. This review provides a comprehensive overview of the MMC literature, including a phylogenetic analysis of MMC distribution across vertebrate species. Here, we discuss the current understanding of the MMCs function in immunity and lingering questions. We suggest additional experiments needed to confirm that MMCs serve a GC-like role in fish immunity. Finally, we address the utility of the MMC as a broadly applicable histological indicator of the fish (as well as amphibian and reptilian) immune response in both laboratory and wild populations of both model and non-model vertebrates. We highlight the factors (sex, pollution exposure, stress, stocking density, etc.) that should be considered when using MMCs to study immunity in non-model vertebrates in wild populations.

Genomic Methods Take the Plunge: Recent Advances in High-Throughput Sequencing of Marine Mammals

The dramatic increase in the application of genomic techniques to non-model organisms (NMOs) over the past decade has yielded numerous valuable contributions to evolutionary biology and ecology, many of which would not have been possible with traditional genetic markers. We review this recent progression with a particular focus on genomic studies of marine mammals, a group of taxa that represent key macroevolutionary transitions from terrestrial to marine environments and for which available genomic resources have recently undergone notable rapid growth. Genomic studies of NMOs utilize an expanding range of approaches, including whole genome sequencing, restriction site-associated DNA sequencing, array-based sequencing of single nucleotide polymorphisms and target sequence probes (e.g., exomes), and transcriptome sequencing. These approaches generate different types and quantities of data, and many can be applied with limited or no prior genomic resources, thus overcoming one traditional limitation of research on NMOs. Within marine mammals, such studies have thus far yielded significant contributions to the fields of phylogenomics and comparative genomics, as well as enabled investigations of fitness, demography, and population structure. Here we review the primary options for generating genomic data, introduce several emerging techniques, and discuss the suitability of each approach for different applications in the study of NMOs.

Integrated studies of organismal plasticity through physiological and transcriptomic approaches: examples from marine polar regions

Transcriptomic methods are now widely used in functional genomic research. The vast amount of information received from these studies comes along with the challenge of developing a precise picture of the functional consequences and the characteristic regulatory mechanisms. Here we assess recent studies in marine species and their adaptation to polar (and seasonal) cold and explore how they have been able to draw reliable conclusions from transcriptomic patterns on functional consequences in the organisms. Our analysis indicates that the interpretation of transcriptomic data suffers from insufficient understanding of the consequences for whole organism performance and fitness and comes with the risk of supporting only preliminary and superficial statements.We propose that the functional understanding of transcriptomic data may be improved by their tighter integration into overarching physiological concepts that support the more specific interpretation of the 'omics' data and, at the same time, can be developed further through embedding the transcriptomic phenomena observed. Such possibilities have not been fully exploited.In the context of thermal adaptation and limitation, we explore preliminary evidence that the concept of oxygen and capacity limited thermal tolerance (OCLTT) may provide sufficient complexity to guide the integration of such data and the development of associated functional hypotheses. At the same time, we identify a lack of methodological approaches linking genes and function to higher levels of integration, in terms of organism and ecosystem functioning, at temporal and geographical scales, to support more reliable conclusions and be predictive with respect to the effects of global changes.

Immunoglobulin detection in wild birds: effectiveness of three secondary anti-avian IgY antibodies in direct ELISAs in 41 avian species

Immunological reagents for wild, non-model species are limited or often non-existent for many species.In this study, we compare the reactivity of a new anti-passerine IgY secondary antibody with existing secondary antibodies developed for use with birds. Samples from 41 species from the following six avian orders were analysed: Anseriformes (1 family, 1 species), Columbiformes (1 family, 2 species), Galliformes (1 family, 1 species), Passeriformes (16 families, 34 species), Piciformes (1 family, 2 species) and Suliformes (1 family, 1 species). Direct ELISAs were performed to detect total IgY using goat anti-passerine IgY, goat anti-chicken IgY or goat anti-bird IgY secondary antibodies.The anti-passerine antibody exhibited significantly higher IgY reactivity compared to the anti-chicken and/or anti-bird antibodies in 80% of the passerine families tested. Birds in the order Piciformes (woodpeckers) and order Suliformes (cormorants) were poorly detected by all three secondary antibodies. A comparison of serum and plasma IgY levels was made within the same individuals for two passerine species (house finch and white-crowned sparrow), and serum exhibited significantly more IgY than the plasma for all three secondary antibodies. This result indicates that serum may be preferred to plasma when measuring total antibody levels in blood.This study indicates that the anti-passerine IgY secondary antibody can effectively be used in immunological assays to detect passerine IgY for species in most passerine families and is preferred over anti-chicken and anti-bird secondary antibodies for the majority of passerine species. This anti-passerine antibody will allow for more accurate detection and quantification of IgY in more wild bird species than was possible with previously available secondary antibodies.

Differential DNA Methylation Analysis without a Reference Genome

Genome-wide DNA methylation mapping uncovers epigenetic changes associated with animal development, environmental adaptation, and species evolution. To address the lack of high-throughput methods for DNA methylation analysis in non-model organisms, we developed an integrated approach for studying DNA methylation differences independent of a reference genome. Experimentally, our method relies on an optimized 96-well protocol for reduced representation bisulfite sequencing (RRBS), which we have validated in nine species (human, mouse, rat, cow, dog, chicken, carp, sea bass, and zebrafish). Bioinformatically, we developed the RefFreeDMA software to deduce ad hoc genomes directly from RRBS reads and to pinpoint differentially methylated regions between samples or groups of individuals (http://RefFreeDMA.computational-epigenetics.org). The identified regions are interpreted using motif enrichment analysis and/or cross-mapping to annotated genomes. We validated our method by reference-free analysis of cell-type-specific DNA methylation in the blood of human, cow, and carp. In summary, we present a cost-effective method for epigenome analysis in ecology and evolution, which enables epigenome-wide association studies in natural populations and species without a reference genome.

Characterization of MHC class II B polymorphism in multiple populations of wild gorillas using non-invasive samples and next-generation sequencing

Genes encoded by the major histocompatibility complex (MHC) are crucial for the recognition and presentation of antigens to the immune system. In contrast to their closest relatives, chimpanzees and humans, much less is known about variation in gorillas at these loci. This study explored the exon 2 variation of -DPB1, -DQB1, and -DRB genes in 46 gorillas from four populations while simultaneously evaluating the feasibility of using fecal samples for high-throughput MHC genotyping. By applying strict similarity- and frequency-based analysis, we found, despite our modest sample size, a total of 18 alleles that have not been described previously, thereby illustrating the potential for efficient and highly accurate MHC genotyping from non-invasive DNA samples. We emphasize the importance of controlling for multiple potential sources of error when applying this massively parallel short-read sequencing technology to PCR products generated from low concentration DNA extracts. We observed pronounced differences in MHC variation between species, subspecies and populations that are consistent with both the ancient and recent demographic histories experienced by gorillas.

Mechanisms of bacterial morphogenesis: evolutionary cell biology approaches provide new insights

How Darwin's "endless forms most beautiful" have evolved remains one of the most exciting questions in biology. The significant variety of bacterial shapes is most likely due to the specific advantages they confer with respect to the diverse environments they occupy. While our understanding of the mechanisms generating relatively simple shapes has improved tremendously in the last few years, the molecular mechanisms underlying the generation of complex shapes and the evolution of shape diversity are largely unknown. The emerging field of bacterial evolutionary cell biology provides a novel strategy to answer this question in a comparative phylogenetic framework. This relatively novel approach provides hypotheses and insights into cell biological mechanisms, such as morphogenesis, and their evolution that would have been difficult to obtain by studying only model organisms. We discuss the necessary steps, challenges, and impact of integrating "evolutionary thinking" into bacterial cell biology in the genomic era.

SNPMeta: SNP annotation and SNP metadata collection without a reference genome

The increase in availability of resequencing data is greatly accelerating SNP discovery and has facilitated the development of SNP genotyping assays. This, in turn, is increasing interest in annotation of individual SNPs. Currently, these data are only available through curation, or comparison to a reference genome. Many species lack a reference genome, but are still important genetic models or are significant species in agricultural production or natural ecosystems. For these species, it is possible to annotate SNPs through comparison with cDNA, or data from well-annotated genes in public repositories. We present SNPMeta, a tool which gathers information about SNPs by comparison with sequences present in GenBank databases. SNPMeta is able to annotate SNPs from contextual sequence in SNP assay designs, and SNPs discovered through genotyping by sequencing (GBS) approaches. However, SNPs discovered through GBS occur throughout the genome, rather than only in gene space, and therefore do not annotate at high rates. SNPMeta can therefore be used to annotate SNPs in nonmodel species or species that lack a reference genome. Annotations generated by SNPMeta are highly concordant with annotations that would be obtained from a reference genome.