Transcriptome analysis of different developmental stages of amphioxus reveals dynamic changes of distinct classes of genes during development

A transcriptomic examination of encased rotifer embryos reveals the developmental trajectory leading to long-term dormancy; are they "animal seeds"?

BACKGROUND

Organisms from many distinct evolutionary lineages acquired the capacity to enter a dormant state in response to environmental conditions incompatible with maintaining normal life activities. Most studied organisms exhibit seasonal or annual episodes of dormancy, but numerous less studied organisms enter long-term dormancy, lasting decades or even centuries. Intriguingly, many planktonic animals produce encased embryos known as resting eggs or cysts that, like plant seeds, may remain dormant for decades. Herein, we studied a rotifer Brachionus plicatilis as a model planktonic species that forms encased dormant embryos via sexual reproduction and non-dormant embryos via asexual reproduction and raised the following questions: Which genes are expressed at which time points during embryogenesis? How do temporal transcript abundance profiles differ between the two types of embryos? When does the cell cycle arrest? How do dormant embryos manage energy?

RESULTS

As the molecular developmental kinetics of encased embryos remain unknown, we employed single embryo RNA sequencing (CEL-seq) of samples collected during dormant and non-dormant embryogenesis. We identified comprehensive and temporal transcript abundance patterns of genes and their associated enriched functional pathways. Striking differences were uncovered between dormant and non-dormant embryos. In early development, the cell cycle-associated pathways were enriched in both embryo types but terminated with fewer nuclei in dormant embryos. As development progressed, the gene transcript abundance profiles became increasingly divergent between dormant and non-dormant embryos. Organogenesis was suspended in dormant embryos, concomitant with low transcript abundance of homeobox genes, and was replaced with an ATP-poor preparatory phase characterized by very high transcript abundance of genes encoding for hallmark dormancy proteins (e.g., LEA proteins, sHSP, and anti-ROS proteins, also found in plant seeds) and proteins involved in dormancy exit. Surprisingly, this period appeared analogous to the late maturation phase of plant seeds.

CONCLUSIONS

The study highlights novel divergent temporal transcript abundance patterns between dormant and non-dormant embryos. Remarkably, several convergent functional solutions appear during the development of resting eggs and plant seeds, suggesting a similar preparatory phase for long-term dormancy. This study accentuated the broad novel molecular features of long-term dormancy in encased animal embryos that behave like "animal seeds".

Three amphioxus reference genomes reveal gene and chromosome evolution of chordates

The slow-evolving invertebrate amphioxus has an irreplaceable role in advancing our understanding of the vertebrate origin and innovations. Here we resolve the nearly complete chromosomal genomes of three amphioxus species, one of which best recapitulates the 17 chordate ancestor linkage groups. We reconstruct the fusions, retention, or rearrangements between descendants of whole-genome duplications, which gave rise to the extant microchromosomes likely existed in the vertebrate ancestor. Similar to vertebrates, the amphioxus genome gradually establishes its three-dimensional chromatin architecture at the onset of zygotic activation and forms two topologically associated domains at the Hox gene cluster. We find that all three amphioxus species have ZW sex chromosomes with little sequence differentiation, and their putative sex-determining regions are nonhomologous to each other. Our results illuminate the unappreciated interspecific diversity and developmental dynamics of amphioxus genomes and provide high-quality references for understanding the mechanisms of chordate functional genome evolution.

The boom and bust of the aphid's essential amino acid metabolism across nymphal development

Within long-term symbioses, animals integrate their physiology and development with their symbiont. In a model nutritional mutualism, aphids harbor the endosymbiont, Buchnera, within specialized bacteriocyte cells. Buchnera synthesizes essential amino acids (EAAs) and vitamins for their host, which are lacking from the aphid's plant sap diet. It is unclear if the aphid host differentially expresses aphid EAA metabolism pathways and genes that collaborate with Buchnera for the production of EAA and vitamins throughout nymphal development when feeding on plants. It is also unclear if aphid bacteriocytes are differentially methylated throughout aphid development as DNA methylation may play a role in gene regulation. By analyzing aphid gene expression, we determined that the bacteriocyte is metabolically more active in metabolizing Buchnera's EAAs and vitamins early in nymphal development compared to intermediate or later immature and adult lifestages. The largest changes in aphid bacteriocyte gene expression, especially for aphid genes that collaborate with Buchnera, occurred during the 3rd to 4th instar transition. During this transition, there is a huge shift in the bacteriocyte from a high energy "nutrient-consuming state" to a "recovery and growth state" where patterning and signaling genes and pathways are upregulated and differentially methylated, and de novo methylation is reduced as evidenced by homogenous DNA methylation profiles after the 2nd instar. Moreover, bacteriocyte number increased and Buchnera's titer decreased throughout aphid nymphal development. These data suggest in combination that bacteriocytes of older nymphal and adult lifestages depend less on the nutritional symbiosis compared to early nymphal lifestages.

An Updated Staging System for Cephalochordate Development: One Table Suits Them All

Chordates are divided into three subphyla: Vertebrata, Tunicata, and Cephalochordata. Phylogenetically, the Cephalochordata, more commonly known as lancelets or amphioxus, constitute the sister group of Vertebrata and Tunicata. Lancelets are small, benthic, marine filter feeders, and their roughly three dozen described species are divided into three genera: Branchiostoma, Epigonichthys, and Asymmetron. Due to their phylogenetic position and their stereotypical chordate morphology and genome architecture, lancelets are key models for understanding the evolutionary history of chordates. Lancelets have thus been studied by generations of scientists, with the first descriptions of adult anatomy and developmental morphology dating back to the 19th century. Today, several different lancelet species are used as laboratory models, predominantly for developmental, molecular and genomic studies. Surprisingly, however, a universal staging system and an unambiguous nomenclature for developing lancelets have not yet been adopted by the scientific community. In this work, we characterized the development of the European lancelet (Branchiostoma lanceolatum) using confocal microscopy and compiled a streamlined developmental staging system, from fertilization through larval life, including an unambiguous stage nomenclature. By tracing growth curves of the European lancelet reared at different temperatures, we were able to show that our staging system permitted an easy conversion of any developmental time into a specific stage name. Furthermore, comparisons of embryos and larvae from the European lancelet (B. lanceolatum), the Florida lancelet (Branchiostoma floridae), two Asian lancelets (Branchiostoma belcheri and Branchiostoma japonicum), and the Bahamas lancelet (Asymmetron lucayanum) demonstrated that our staging system could readily be applied to other lancelet species. Although the detailed staging description was carried out on developing B. lanceolatum, the comparisons with other lancelet species thus strongly suggested that both staging and nomenclature are applicable to all extant lancelets. We conclude that this description of embryonic and larval development will be of great use for the scientific community and that it should be adopted as the new standard for defining and naming developing lancelets. More generally, we anticipate that this work will facilitate future studies comparing representatives from different chordate lineages.

Endodermal Maternal Transcription Factors Establish Super-Enhancers during Zygotic Genome Activation

Elucidation of the sequence of events underlying the dynamic interaction between transcription factors and chromatin states is essential. Maternal transcription factors function at the top of the regulatory hierarchy to specify the primary germ layers at the onset of zygotic genome activation (ZGA). We focus on the formation of endoderm progenitor cells and examine the interactions between maternal transcription factors and chromatin state changes underlying the cell specification process. Endoderm-specific factors Otx1 and Vegt together with Foxh1 orchestrate endoderm formation by coordinated binding to select regulatory regions. These interactions occur before the deposition of enhancer histone marks around the regulatory regions, and these TFs recruit RNA polymerase II, regulate enhancer activity, and establish super-enhancers associated with important endodermal genes. Therefore, maternal transcription factors Otx1, Vegt, and Foxh1 combinatorially regulate the activity of super-enhancers, which in turn activate key lineage-specifying genes during ZGA.

How do maternal transcription factors interact with chromatin regions to coordinate the endodermal gene regulatory program? Paraiso et al. demonstrate that combinatorial binding of maternal Otx1, Vegt, and Foxh1 to select enhancers and super-enhancers in the genome controls endodermal cell fate specification during zygotic gene activation.

Transcriptome sequencing and characterization of Astragalus membranaceus var. mongholicus root reveals key genes involved in flavonoids biosynthesis

BACKGROUND

Astragalus membranaceus (Fisch.) Bge. var. mongholicus (Bge.) Hsiao is a traditional medicinal herb of Leguminosae since it contains bioactive compounds such as flavonoids, which have significant pharmacological effects on immunity and antioxidant. However, the scanty genomic and transcriptome resources of Astragalus membranaceus have hindered further exploration of its biosynthesis and accumulation mechanism.

OBJECTIVE

This project aim to further improve our understanding of the relationship between transcriptional behavior and flavonoids content of A. mongholicus.

METHODS

The accumulation of flavonoids and related gene expression in five different developmental stages (A: vegetative, B: florescence, C: fruiting, D: fruit ripening and E: defoliating stages) of A. mongholicus root were studied by combining UV spectrophotometry and transcriptomic techniques. The de novo assembly, annotation and functional evaluation of the contigs were performed with bioinformatics tools.

RESULTS

After screening and assembling the raw data, there were a total of 158,123 unigenes with an average length of 644.89 bp were finally obtained, which has 8362 unigenes could be jointly annotated by NR, SwissProt, eggNOG, GO, KEGG and Pfam databases. KEGG enrichment analysis was performed on differentially expressed genes(DEGs)in the four groups (A vs. B, B vs. C, C vs. D, D vs. E). The results showed that many DEGs in each group were significantly enriched to flavonoids biosynthesis related pathways. Among them, a number of 86 were involved in the biosynthesis of isoflavonoid (12), flavonoid (5) and phenylpropanoid (69). Further analysis of these DEGs revealed that the expression levels of key genes such as PAL, 4CL, CCR, COMT, DFR, etc. were all down-regulated at the fruiting stage, and then raised at the fruit ripening stage. This expression pattern was similar to the accumulation trend of total flavonoids content.

CONCLUSIONS

In summary, this comprehensive transcriptome dataset allowed the identification of genes associated with flavonoids metabolic pathways. The results laid a foundation for the biosynthesis and regulation of flavonoids. It also provided a scientific basis for the most suitable harvest time and resource utilization of A. mongholicus.

Evolutionary History of GLIS Genes Illuminates Their Roles in Cell Reprograming and Ciliogenesis

The GLIS family transcription factors, GLIS1 and GLIS3, potentiate generation of induced pluripotent stem cells (iPSCs). In contrast, another GLIS family member, GLIS2, suppresses cell reprograming. To understand how these disparate roles arose, we examined evolutionary origins and genomic organization of GLIS genes. Comprehensive phylogenetic analysis shows that GLIS1 and GLIS3 originated during vertebrate whole genome duplication, whereas GLIS2 is a sister group to the GLIS1/3 and GLI families. This result is consistent with their opposing functions in cell reprograming. Glis1 evolved faster than Glis3, losing many protein-interacting motifs. This suggests that Glis1 acquired new functions under weakened evolutionary constraints. In fact, GLIS1 induces induced pluripotent stem cells more strongly. Transcriptomic data from various animal embryos demonstrate that glis1 is maternally expressed in some tetrapods, whereas vertebrate glis3 and invertebrate glis1/3 genes are rarely expressed in oocytes, suggesting that vertebrate (or tetrapod) Glis1 acquired a new expression domain and function as a maternal factor. Furthermore, comparative genomic analysis reveals that glis1/3 is part of a bilaterian-specific gene cluster, together with rfx3, ndc1, hspb11, and lrrc42. Because known functions of these genes are related to cilia formation and function, the last common ancestor of bilaterians may have acquired this cluster by shuffling gene order to establish more sophisticated epithelial tissues involving cilia. This evolutionary study highlights the significance of GLIS1/3 for cell reprograming, development, and diseases in ciliated organs such as lung, kidney, and pancreas.

LanceletDB: an integrated genome database for lancelet, comparing domain types and combination in orthologues among lancelet and other species

Lancelet (amphioxus) represents the most basally divergent extant chordate (cephalochordates) that diverged from the other two chordate lineages (urochordates and vertebrates) more than half a billion years ago. As it occupies a key position in evolution, it is considered as one of the best proxies for understanding the chordate ancestral state. Thus, the construction of a database with multiple lancelet genomes and gene annotation data, including protein domains, is urgently needed to investigate the loss and gain of domains in orthologues among species, especially ancient domain types (non-vertebrate-specific domains) and novel domain combination, which is helpful for providing new insight into the chordate ancestral state and vertebrate evolution. Here, we present an integrated genome database for lancelet, LanceletDB, which provides reference haploid genome sequence and annotation data for lancelet (Branchiostoma belcheri), including gene models and annotation, protein domain types, gene expression pattern in embryogenesis, different expression sequence tag sets and alternative polyadenylation (APA) sites profiled by the sequencing APA sites method. Especially, LanceletDB allows comparison of domain types and combination in orthologues among type species so as to decode the ancient domain types and novel domain combination during evolution. We also integrated the released diploid lancelet genome annotation data (Branchiostoma floridae) to expand LanceletDB and extend its usefulness. These data are available through the search and analysis page, basic local alignment search tool page and genome browser to provide an integrated display.

Nature-inspired engineering of an F-type lectin for increased binding strength

Individual lectin-carbohydrate interactions are usually of low affinity. However, high avidity is frequently attained by the multivalent presentation of glycans on biological surfaces coupled with the occurrence of high order lectin oligomers or tandem repeats of lectin domains in the polypeptide. F-type lectins are l-fucose binding lectins with a typical sequence motif, HX(26)RXDX(4)R/K, whose residues participate in l-fucose binding. We previously reported the presence of a few eukaryotic F-type lectin domains with partial sequence duplication that results in the presence of two l-fucose-binding sequence motifs. We hypothesized that such partial sequence duplication would result in greater avidity of lectin-ligand interactions. Inspired by this example from Nature, we attempted to engineer a bacterial F-type lectin domain from Streptosporangium roseum to attain avid binding by mimicking partial duplication. The engineered lectin demonstrated 12-fold greater binding strength than the wild-type lectin to multivalent fucosylated glycoconjugates. However, the affinity to the monosaccharide l-fucose in solution was similar and partial sequence duplication did not result in an additional functional l-fucose binding site. We also cloned, expressed and purified a Branchiostoma floridae F-type lectin domain with naturally occurring partial sequence duplication and confirmed that the duplicated region with the F-type lectin sequence motif did not participate in l-fucose binding. We found that the greater binding strength of the engineered lectin from S. roseum was instead due to increased oligomerization. We believe that this Nature-inspired strategy might be useful for engineering lectins to improve binding strength in various applications.

Comparative transcriptome analysis provides global insight into gene expression differences between two orchid cultivars

The orchids GL and YL are two cultivars of Cymbidium longibracteatum. YL displays an obviously yellowing rhizome and yellow leaves, while GL ('Longchangsu') shows dark green leaves and greenish rhizome. But the molecular mechanism for the differences between the two cultivars is poorly understood. In the present study, we showed that the structure of chloroplasts was significantly damaged in YL. Biochemical analysis uncovered the contents of chlorophyll a, chlorophyll b, total chlorophyll and carotenoid were notably decreased in YL. Using RNA-Seq technology, more than 38 million clean reads were generated in each pool, and 116,422 unigenes were assembled de novo. 6,660 unigenes with differential expression patterns (FDR≤0.01 and |log2 ratio|≥1) were totally identified between the two cultivars. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of differentially expressed unigenes (DEGs) suggested 33 KEGG pathways were notably enriched, including biological processes such as “phenylpropanoid biosynthesis”, “phagosome”, “starch and sucrose metabolism”, “drug metabolism—cytochrome P450”, “fatty acid elongation”, and “flavone and flavonol biosynthesis”. Further analysis revealed that chlorophyll degeneration related unigene (c48794_g1) and flavonoid biosynthesis related unigenes (c16388_g1, c48963_g1, c63571_g1, c4492_g1, c52282_g1, c78740_g1, c4645_g1) were up-regulated while carotenoid biosynthesis related unigene (c7212_g1) were down-regulated in YL. Additionally, six of NAC, R2R3-MYB, bHLH transcription factors (c42861_g1, c105949_g1, c61265_g1, c42659_g1, c82171_g1, c19158_g1) might be involved in regulation of pigment biosynthesis. The chlorophyll degeneration and the flavonoid biosynthesis related unigenes up-regulation together with the carotenoid biosynthesis related unigenes down-regulation may contribute to the yellowing phenotype of YL.

Hagfish and lamprey Hox genes reveal conservation of temporal colinearity in vertebrates

Hox genes exert fundamental roles for proper regional specification along the main rostro-caudal axis of animal embryos. They are generally expressed in restricted spatial domains according to their position in the cluster (spatial colinearity)-a feature that is conserved across bilaterians. In jawed vertebrates (gnathostomes), the position in the cluster also determines the onset of expression of Hox genes (a feature known as whole-cluster temporal colinearity (WTC)), while in invertebrates this phenomenon is displayed as a subcluster-level temporal colinearity. However, little is known about the expression profile of Hox genes in jawless vertebrates (cyclostomes); therefore, the evolutionary origin of WTC, as seen in gnathostomes, remains a mystery. Here, we show that Hox genes in cyclostomes are expressed according to WTC during development. We investigated the Hox repertoire and Hox gene expression profiles in three different species-a hagfish, a lamprey and a shark-encompassing the two major groups of vertebrates, and found that these are expressed following a whole-cluster, temporally staggered pattern, indicating that WTC has been conserved during the past 500 million years despite drastically different genome evolution and morphological outputs between jawless and jawed vertebrates.

Transcriptome analysis of different growth stages of Aspergillus oryzae reveals dynamic changes of distinct classes of genes during growth

Background

The gene expression profile and metabolic pathways of Aspergillus oryzae underlying the anatomical and morphological differentiation across different growth stages have not been fully characterized. The rapid development of next-generation sequencing technologies provides advanced knowledge of the genomic organization of A. oryzae.

Results

In this study, we characterized the growth and development of A. oryzae at different growth stages, including the adaptive phase, logarithmic phase, and stationary phase. Our results revealed that A. oryzae undergoes physiological and morphological differentiation across the different stages. RNA-seq was employed to analyze the three stages of A. oryzae, which generated more than 27 million high-quality reads per sample. The analysis of differential gene expression showed more genes expressed differentially upon transition from the adaptive phase to the logarithmic and stationary phases, while relatively steady trend was observed during the transition from the logarithmic phase to the stationary phase. GO classification of the differentially expressed genes among different growth stages revealed that most of these genes were enriched for single-organism process, metabolic process, and catalytic activity. These genes were then subjected to a clustering analysis. The results showed that the cluster with the majority of genes with increased expression upon transition from the adaptive phase to the logarithmic phase, and steady expression from the logarithmic phase to the stationary phase was mainly involved in the carbohydrate and amino acid metabolism.

Conclusion

Our results provide a foundation for identifying developmentally important genes and understanding the biological processes across various growth stages.

Electronic supplementary material

The online version of this article (10.1186/s12866-018-1158-z) contains supplementary material, which is available to authorized users.