Saccharina genomes provide novel insight into kelp biology

An improved approach for reconstructing consensus repeats from short sequence reads

BACKGROUND

Repeat elements are important components of most eukaryotic genomes. Most existing tools for repeat analysis rely either on high quality reference genomes or existing repeat libraries. Thus, it is still challenging to do repeat analysis for species with highly repetitive or complex genomes which often do not have good reference genomes or annotated repeat libraries. Recently we developed a computational method called REPdenovo that constructs consensus repeat sequences directly from short sequence reads, which outperforms an existing tool called RepARK. One major issue with REPdenovo is that it doesn't perform well for repeats with relatively high divergence rates or low copy numbers. In this paper, we present an improved approach for constructing consensus repeats directly from short reads. Comparing with the original REPdenovo, the improved approach uses more repeat-related k-mers and improves repeat assembly quality using a consensus-based k-mer processing method.

RESULTS

We compare the performance of the new method with REPdenovo and RepARK on Human, Arabidopsis thaliana and Drosophila melanogaster short sequencing data. And the new method fully constructs more repeats in Repbase than the original REPdenovo and RepARK, especially for repeats of higher divergence rates and lower copy number. We also apply our new method on Hummingbird data which doesn't have a known repeat library, and it constructs many repeat elements that can be validated using PacBio long reads.

CONCLUSION

We propose an improved method for reconstructing repeat elements directly from short sequence reads. The results show that our new method can assemble more complete repeats than REPdenovo (and also RepARK). Our new approach has been implemented as part of the REPdenovo software package, which is available for download at https://github.com/Reedwarbler/REPdenovo .

Phaeoviral Infections Are Present in Macrocystis, Ecklonia and Undaria (Laminariales) and Are Influenced by Wave Exposure in Ectocarpales

Two sister orders of the brown macroalgae (class Phaeophyceae), the morphologically complex Laminariales (commonly referred to as kelp) and the morphologically simple Ectocarpales are natural hosts for the dsDNA phaeoviruses (family Phycodnaviridae) that persist as proviruses in the genomes of their hosts. We have previously shown that the major capsid protein (MCP) and DNA polymerase concatenated gene phylogeny splits phaeoviruses into two subgroups, A and B (both infecting Ectocarpales), while MCP-based phylogeny suggests that the kelp phaeoviruses form a distinct third subgroup C. Here we used MCP to better understand the host range of phaeoviruses by screening a further 96 and 909 samples representing 11 and 3 species of kelp and Ectocarpales, respectively. Sporophyte kelp samples were collected from their various natural coastal habitats spanning five continents: Africa, Asia, Australia, Europe, and South America. Our phylogenetic analyses showed that while most of the kelp phaeoviruses, including one from Macrocystispyrifera, belonged to the previously designated subgroup C, new lineages of Phaeovirus in 3 kelp species, Ecklonia maxima, Ecklonia radiata, Undaria pinnatifida, grouped instead with subgroup A. In addition, we observed a prevalence of 26% and 63% in kelp and Ectocarpales, respectively. Although not common, multiple phaeoviral infections per individual were observed, with the Ectocarpales having both intra- and inter-subgroup phaeoviral infections. Only intra-subgroup phaeoviral infections were observed in kelp. Furthermore, prevalence of phaeoviral infections within the Ectocarpales is also linked to their exposure to waves. We conclude that phaeoviral infection is a widely occurring phenomenon in both lineages, and that phaeoviruses have diversified with their hosts at least since the divergence of the Laminariales and Ectocarpales.

Comparative analysis of plant MKK gene family reveals novel expansion mechanism of the members and sheds new light on functional conservation

Background

Mitogen-activated protein kinase (MAPK) cascades play critical functions in almost every aspect of plant growth and development, which regulates many physiological and biochemical processes. As a middle nodal point of the MAPK cascades, although evolutionary analysis of MKK from individual plant families had some reports, their evolutionary history in entire plants is still not clear.

Results

To better understand the evolution and function of plant MKKs, we performed systematical molecular evolutionary analysis of the MAPKK gene family and also surveyed their gene organizations, sequence features and expression patterns in different subfamilies. Phylogenetic analysis showed that plant MAPKK fall into five different groups (Group A–E). Majority orthology groups seemed to be a single or low-copy genes in all plant species analyzed in Group B, C and D, whereas group A MKKs undergo several duplication events, generating multiple gene copies. Further analysis showed that these duplication events were on account of whole genome duplications (WGDs) in plants and the duplicate genes maybe have undergone functional divergence. We also found that group E MKKs had mutation with one change of serine or theronine might lead to inactivity originated through the ancient tandem duplicates in monocots. Moreover, we also identified MKK3 integrated NTF2 domain that might have gradually lost the cytoplasmic-nuclear trafficking activity, which suggests that they may involve with the gene function more and more sophistication in the evolutionary process. Moreover, expression analyses indicated that plant MKK genes play probable roles in UV-B signaling.

Conclusion

In general, ancient gene and genome duplications are significantly conducive to the expansion of the plant MKK gene family. Our study reveals two distinct evolutionary patterns for plant MKK proteins and sheds new light on the functional evolution of this gene family.

Electronic supplementary material

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

Genomic analyses of unique carbohydrate and phytohormone metabolism in the macroalga Gracilariopsis lemaneiformis (Rhodophyta)

BACKGROUND

Red algae are economically valuable for food and in industry. However, their genomic information is limited, and the genomic data of only a few species of red algae have been sequenced and deposited recently. In this study, we annotated a draft genome of the macroalga Gracilariopsis lemaneiformis (Gracilariales, Rhodophyta).

RESULTS

The entire 88.98 Mb genome of Gp. lemaneiformis 981 was generated from 13,825 scaffolds (≥500 bp) with an N50 length of 30,590 bp, accounting for approximately 91% of this algal genome. A total of 38.73 Mb of scaffold sequences were repetitive, and 9281 protein-coding genes were predicted. A phylogenomic analysis of 20 genomes revealed the relationship among the Chromalveolata, Rhodophyta, Chlorophyta and higher plants. Homology analysis indicated phylogenetic proximity between Gp. lemaneiformis and Chondrus crispus. The number of enzymes related to the metabolism of carbohydrates, including agar, glycoside hydrolases, glycosyltransferases, was abundant. In addition, signaling pathways associated with phytohormones such as auxin, salicylic acid and jasmonates are reported for the first time for this alga.

CONCLUSION

We sequenced and analyzed a draft genome of the red alga Gp. lemaneiformis, and revealed its carbohydrate metabolism and phytohormone signaling characteristics. This work will be helpful in research on the functional and comparative genomics of the order Gracilariales and will enrich the genomic information on marine algae.

Functional genomics analysis reveals the biosynthesis pathways of important cellular components (alginate and fucoidan) of Saccharina

Although alginate and fucoidan are unique cellular components and have important biological significance in brown algae, and many possible involved genes are present in brown algal genomes, their functions and regulatory mechanisms have not been fully revealed. Both polysaccharides may play important roles in the evolution of multicellular brown algae, but specific and in-depth studies are still limited. In this study, a functional genomics analysis of alginate and fucoidan biosynthesis routes was conducted in Saccharina, and the key events in these pathways in brown algae were identified. First, genes from different sources, including eukaryotic hosts via endosymbiotic gene transfer and bacteria via horizontal gene transfer, were combined to build a complete pathway framework. Then, a critical event occurred to drive these pathways to have real function: one of the mannose-6-phosphate isomerase homologs that arose by gene duplication subsequently adopted the function of the mannose-1-phosphate guanylyltransferase (MGP) gene, which was absent in algal genomes. Further, downstream pathway genes proceeded with gene expansions and complex transcriptional mechanisms, which may be conducive to the synthesis of alginate and fucoidan with diverse structures and contents depending on the developmental stage, tissue structure, and environmental conditions. This study revealed the alginate and fucoidan synthesis pathways and all included genes from separate phylogenetic sources in brown algae. Enzyme assays confirmed the function of key genes and led to the determination of a substitute for the missing MPG. All gene families had constitutively expressed member(s) to maintain the basic synthesis; and the gene function differentiation, enzyme characterization and gene expression regulation differences separated brown algae from other algae lineages and were considered to be the major driving forces for sophisticated system evolution of brown algae.

Genome-wide searches and molecular analyses highlight the unique evolutionary path of flavone synthase I (FNSI) in Apiaceae

Flavone synthase is a key enzyme for flavone biosynthesis and is encoded by two gene families: flavone synthase I (FNSI) and flavone synthase II (FNSII). FNSII is widely distributed in plants, while FNSI has been reported in rice (Oryza sativa) and seven species of Apiaceae. FNSI has likely evolved from the duplication of flavanone 3β-hydroxylase (F3H). In this study, we used multiple bioinformatics tools to identify putative FNSI and F3H genes from 42 publicly available genome and transcriptome datasets. Results showed that rice FNSI does not share a common ancestral sequence with other known FNSI genes and that FNSI is absent from species outside of Apiaceae. Positive selection site identification analysis revealed that four sites within the FNSI tree branches of Apiaceae evolved under significant positive selection. The putative F3H genes identified in this study provide a valuable resource for further function analysis of flavone synthase.

Evolution and Expansion of the Prokaryote-Like Lipoxygenase Family in the Brown Alga Saccharina japonica

Lipoxygenase (LOX) plays important roles in fatty acid oxidation and lipid mediator biosynthesis. In this study, we give first insights into brown algal LOX evolution. Whole genome searches revealed four, three, and eleven LOXs in Ectocarpus siliculosus, Cladosiphon okamuranus, and Saccharina japonica, respectively. In phylogenetic analyses, LOXs from brown algae form a robust clade with those from prokaryotes, suggesting an ancestral origin and slow evolution. Brown algal LOXs were divided into two clades, C1 and C2 in a phylogenetic tree. Compared to the two species of Ectocarpales, LOX gene expansion occurred in the kelp S. japonica through tandem duplication and segmental duplication. Selection pressure analysis showed that LOX genes in brown algae have undergone strong purifying selection, while the selective constraint in the C2 clade was more relaxed than that in the C1 clade. Furthermore, within each clade, LOXs of S. japonica evolved under more relaxed selection constraints than E. siliculosus and C. okamuranus. Structural modeling showed that unlike LOXs of plants and animals, which contain a β barrel in the N-terminal part of the protein, LOXs in brown algae fold into a single domain. Analysis of previously published transcriptomic data showed that LOXs in E. siliculosus are responsive to hyposaline, hypersaline, oxidative, and copper stresses. Moreover, clear divergence of expression patterns was observed among different life stages, as well as between duplicate gene pairs. In E. siliculosus, all four LOXs are male-biased in immature gametophytes, and mature gametophytes showed significantly higher LOX mRNA levels than immature gametophytes and sporophytes. In S. japonica, however, our RNA-Seq data showed that most LOXs are highly expressed in sporophytes. Even the most recently duplicated gene pairs showed divergent expression patterns, suggesting that functional divergence has likely occurred since LOX genes duplicated, which potentially contributes to the production of various oxylipins in brown algae.

Detection of bacterial contaminants and hybrid sequences in the genome of the kelp Saccharina japonica using Taxoblast

Modern genome sequencing strategies are highly sensitive to contamination making the detection of foreign DNA sequences an important part of analysis pipelines. Here we use Taxoblast, a simple pipeline with a graphical user interface, for the post-assembly detection of contaminating sequences in the published genome of the kelp Saccharina japonica. Analyses were based on multiple blastn searches with short sequence fragments. They revealed a number of probable bacterial contaminations as well as hybrid scaffolds that contain both bacterial and algal sequences. This or similar types of analysis, in combination with manual curation, may thus constitute a useful complement to standard bioinformatics analyses prior to submission of genomic data to public repositories. Our analysis pipeline is open-source and freely available at http://sdittami.altervista.org/taxoblast and via SourceForge (https://sourceforge.net/projects/taxoblast).

Pervasive, Genome-Wide Transcription in the Organelle Genomes of Diverse Plastid-Bearing Protists

Organelle genomes are among the most sequenced kinds of chromosome. This is largely because they are small and widely used in molecular studies, but also because next-generation sequencing technologies made sequencing easier, faster, and cheaper. However, studies of organelle RNA have not kept pace with those of DNA, despite huge amounts of freely available eukaryotic RNA-sequencing (RNA-seq) data. Little is known about organelle transcription in nonmodel species, and most of the available eukaryotic RNA-seq data have not been mined for organelle transcripts. Here, we use publicly available RNA-seq experiments to investigate organelle transcription in 30 diverse plastid-bearing protists with varying organelle genomic architectures. Mapping RNA-seq data to organelle genomes revealed pervasive, genome-wide transcription, regardless of the taxonomic grouping, gene organization, or noncoding content. For every species analyzed, transcripts covered ≥85% of the mitochondrial and/or plastid genomes (all of which were ≤105 kb), indicating that most of the organelle DNA—coding and noncoding—is transcriptionally active. These results follow earlier studies of model species showing that organellar transcription is coupled and ubiquitous across the genome, requiring significant downstream processing of polycistronic transcripts. Our findings suggest that noncoding organelle DNA can be transcriptionally active, raising questions about the underlying function of these transcripts and underscoring the utility of publicly available RNA-seq data for recovering complete genome sequences. If pervasive transcription is also found in bigger organelle genomes (>105 kb) and across a broader range of eukaryotes, this could indicate that noncoding organelle RNAs are regulating fundamental processes within eukaryotic cells.

A review about brown algal cell walls and fucose-containing sulfated polysaccharides: Cell wall context, biomedical properties and key research challenges

Studies on brown algal cell walls have entered a new phase with the concomitant discovery of novel polysaccharides present in cell walls and the establishment of a comprehensive generic model for cell wall architecture. Brown algal cell walls are composites of structurally complex polysaccharides. In this review we discuss the most recent progress in the structural composition of brown algal cell walls, emphasizing the significance of extraction and screening techniques, and the biological activities of the corresponding polysaccharides, with a specific focus on the fucose-containing sulfated polysaccharides. They include valuable marine molecules that exert a broad range of pharmacological properties such as antioxidant and anti-inflammatory activities, functions in the regulation of immune responses and of haemostasis, anti-infectious and anticancer actions. We identify the key remaining challenges in this research field.

Acquisition and metabolism of carbon in the Ochrophyta other than diatoms

The acquisition and assimilation of inorganic C have been investigated in several of the 15 clades of the Ochrophyta other than diatoms, with biochemical, physiological and genomic data indicating significant mechanistic variation. Form ID Rubiscos in the Ochrophyta are characterized by a broad range of kinetics values. In spite of relatively high K0.5CO2 and low CO2 : O2 selectivity, diffusive entry of CO2 occurs in the Chrysophyceae and Synurophyceae. Eustigmatophyceae and Phaeophyceae, on the contrary, have CO2 concentrating mechanisms, usually involving the direct or indirect use of [Formula: see text] This variability is possibly due to the ecological contexts of the organism. In brown algae, C fixation generally takes place through a classical C3 metabolism, but there are some hints of the occurrence of C4 metabolism and low amplitude CAM in a few members of the Fucales. Genomic data show the presence of a number of potential C4 and CAM genes in Ochrophyta other than diatoms, but the other core functions of many of these genes give a very limited diagnostic value to their presence and are insufficient to conclude that C4 photosynthesis is present in these algae.This article is part of the themed issue 'The peculiar carbon metabolism in diatoms'.

The Algal Revolution

Algae are (mostly) photosynthetic eukaryotes that occupy multiple branches of the tree of life, and are vital for planet function and health. In this review, we highlight a transformative period in studies of the evolution and functioning of this extraordinary group of organisms and their potential for novel applications, wrought by high-throughput 'omic' and reverse genetic methods. We cover the origin and diversification of algal groups, explore advances in understanding the link between phenotype and genotype, consider algal sex determination, and review progress in understanding the roots of algal multicellularity. Experimental evolution studies to determine how algae evolve in changing environments are highlighted, as is their potential as production platforms for compounds of commercial interest, such as biofuel precursors, nutraceuticals, or therapeutics.

Finding a partner in the ocean: molecular and evolutionary bases of the response to sexual cues in a planktonic diatom

Microalgae play a major role as primary producers in aquatic ecosystems. Cell signalling regulates their interactions with the environment and other organisms, yet this process in phytoplankton is poorly defined. Using the marine planktonic diatom Pseudo-nitzschia multistriata, we investigated the cell response to cues released during sexual reproduction, an event that demands strong regulatory mechanisms and impacts on population dynamics. We sequenced the genome of P. multistriata and performed phylogenomic and transcriptomic analyses, which allowed the definition of gene gains and losses, horizontal gene transfers, conservation and evolutionary rate of sex-related genes. We also identified a small number of conserved noncoding elements. Sexual reproduction impacted on cell cycle progression and induced an asymmetric response of the opposite mating types. G protein-coupled receptors and cyclic guanosine monophosphate (cGMP) are implicated in the response to sexual cues, which overall entails a modulation of cell cycle, meiosis-related and nutrient transporter genes, suggesting a fine control of nutrient uptake even under nutrient-replete conditions. The controllable life cycle and the genome sequence of P. multistriata allow the reconstruction of changes occurring in diatoms in a key phase of their life cycle, providing hints on the evolution and putative function of their genes and empowering studies on sexual reproduction.

Multiple gene movements into and out of haploid sex chromosomes

BACKGROUND

Long-term evolution of sex chromosomes is a dynamic process shaped by gene gain and gene loss. Sex chromosome gene traffic has been studied in XY and ZW systems but no detailed analyses have been carried out for haploid phase UV sex chromosomes. Here, we explore sex-specific sequences of seven brown algal species to understand the dynamics of the sex-determining region (SDR) gene content across 100 million years of evolution.

RESULTS

A core set of sex-linked genes is conserved across all the species investigated, but we also identify modifications of both the U and the V SDRs that occurred in a lineage-specific fashion. These modifications involve gene loss, gene gain and relocation of genes from the SDR to autosomes. Evolutionary analyses suggest that the SDR genes are evolving rapidly and that this is due to relaxed purifying selection. Expression analysis indicates that genes that were acquired from the autosomes have been retained in the SDR because they confer a sex-specific role in reproduction. By examining retroposed genes in Saccharina japonica, we demonstrate that UV sex chromosomes have generated a disproportionate number of functional orphan retrogenes compared with autosomes. Movement of genes out of the UV sex chromosome could be a means to compensate for gene loss from the non-recombining region, as has been suggested for Y-derived retrogenes in XY sexual systems.

CONCLUSION

This study provides the first analysis of gene traffic in a haploid UV system and identifies several features of general relevance to the evolution of sex chromosomes.

Sequencing of complete mitochondrial genome of brown algal Saccharina sp. ye-C

The complete sequence (37,609 bp) of the mitochondrial DNA (mtDNA) of the Saccharina sp. ye-C was determined using Illumina sequencing data. The genome contains 38 protein-coding genes (PCG), 3 ribosomal RNA (rRNA), 25 transfer RNA (tRNA) genes that are typical of Saccharina mtDNA. The phylogenetic analysis based on the mitochondrial genomes of brown algae showed that Saccharina sp. ye-C formed a robust clade with Saccharina coriacea, which strongly supports their close phylogenetic relationship.

Genome-wide computational analysis of the secretome of brown algae (Phaeophyceae)

Brown algae have evolved complex multicellularity in the heterokont lineage. They are phylogenetically distant to land plants, fungi and animals. Especially, the members of Laminariales (so-called kelps) have developed highly differentiated tissues. Extracellular matrix (ECM) plays pivotal roles in a number of essential processes in multicellular organisms, such as cell adhesion, cell and tissue differentiations, cell-to-cell communication, and responses to environmental stimuli. In these processes, a set of extracellular secreted proteins called the secretome operates remodeling of the physicochemical nature of ECM and signal transduction by interacting with cell surface proteins and signaling molecules. Characterization of the secretome is a critical step to clarify the contributions of ECM to the multicellularity of brown algae. However, the identity of the brown algal secretome has been poorly understood. In order to reveal the repertory of the brown algal secretome and its involvement in the evolution of Laminariales, we conducted a genome-wide analysis of the brown algal secretome utilizing the published complete genome data of Ectocarpus siliculosus and Saccharina japonica as well as newly obtained RNA-seq data of seven laminarialean species (Agarum clathratum, Alaria crassifolia, Aureophycus aleuticus, Costaria costata, Pseudochorda nagaii, Saccharina angustata and Undaria pinnatifida) largely covering the laminarialean families. We established the in silico pipeline to systematically and accurately detect the secretome by combining multiple prediction algorithms for the N-terminal signal peptide and transmembrane domain within the protein sequence. From 16,189 proteins of E. siliculosus and 18,733 proteins of S. japonica, 552 and 964 proteins respectively were predicted to be classified as the secretome. Conserved domain analysis showed that the domain repertory were very similar to each other, and that of the brown algal secretome was partially common with that of the secretome of other multicellular organisms (land plants, fungi and animals). In the laminarialean species, it was estimated that the gene abundance and the domain architecture of putative ECM remodeling-related proteins were altered compared with those of E. siliculosus, and that the alteration started from the basal group of Laminariales. These results suggested that brown algae have developed their own secretome, and its functions became more elaborated in the more derived members in Laminariales.

Re-annotation, improved large-scale assembly and establishment of a catalogue of noncoding loci for the genome of the model brown alga Ectocarpus

The genome of the filamentous brown alga Ectocarpus was the first to be completely sequenced from within the brown algal group and has served as a key reference genome both for this lineage and for the stramenopiles. We present a complete structural and functional reannotation of the Ectocarpus genome. The large-scale assembly of the Ectocarpus genome was significantly improved and genome-wide gene re-annotation using extensive RNA-seq data improved the structure of 11 108 existing protein-coding genes and added 2030 new loci. A genome-wide analysis of splicing isoforms identified an average of 1.6 transcripts per locus. A large number of previously undescribed noncoding genes were identified and annotated, including 717 loci that produce long noncoding RNAs. Conservation of lncRNAs between Ectocarpus and another brown alga, the kelp Saccharina japonica, suggests that at least a proportion of these loci serve a function. Finally, a large collection of single nucleotide polymorphism-based markers was developed for genetic analyses. These resources are available through an updated and improved genome database. This study significantly improves the utility of the Ectocarpus genome as a high-quality reference for the study of many important aspects of brown algal biology and as a reference for genomic analyses across the stramenopiles.

High-density genetic map and identification of QTLs for responses to temperature and salinity stresses in the model brown alga Ectocarpus

Deciphering the genetic architecture of adaptation of brown algae to environmental stresses such as temperature and salinity is of evolutionary as well as of practical interest. The filamentous brown alga Ectocarpus sp. is a model for the brown algae and its genome has been sequenced. As sessile organisms, brown algae need to be capable of resisting the various abiotic stressors that act in the intertidal zone (e.g. osmotic pressure, temperature, salinity, UV radiation) and previous studies have shown that an important proportion of the expressed genes is regulated in response to hyposaline, hypersaline or oxidative stress conditions. Using the double digest RAD sequencing method, we constructed a dense genetic map with 3,588 SNP markers and identified 39 QTLs for growth-related traits and their plasticity under different temperature and salinity conditions (tolerance to high temperature and low salinity). GO enrichment tests within QTL intervals highlighted membrane transport processes such as ion transporters. Our study represents a significant step towards deciphering the genetic basis of adaptation of Ectocarpus sp. to stress conditions and provides a substantial resource to the increasing list of tools generated for the species.

Rapid Evolution of microRNA Loci in the Brown Algae

Stringent searches for microRNAs (miRNAs) have so far only identified these molecules in animals, land plants, chlorophyte green algae, slime molds and brown algae. The identification of miRNAs in brown algae was based on the analysis of a single species, the filamentous brown alga Ectocarpus sp. Here, we have used deep sequencing of small RNAs and a recently published genome sequence to identify miRNAs in a second brown alga, the kelp Saccharina japonica. S. japonica possesses a large number of miRNAs (117) and these miRNAs are highly diverse, falling into 98 different families. Surprisingly, none of the S. japonica miRNAs share significant sequence similarity with the Ectocarpus sp. miRNAs. However, the miRNA repertoires of the two species share a number of structural and genomic features indicating that they were generated by similar evolutionary processes and therefore probably evolved within the context of a common, ancestral miRNA system. This lack of sequence similarity suggests that miRNAs evolve rapidly in the brown algae (the two species are separated by ∼95 Myr of evolution). The sets of predicted targets of miRNAs in the two species were also very different suggesting that the divergence of the miRNAs may have had significant consequences for miRNA function.

Effect of domestication on the genetic diversity and structure of Saccharina japonica populations in China

Saccharina japonica is a commercially and ecologically important seaweed and is an excellent system for understanding the effects of domestication on marine crops. In this study, we used 19 selected simple sequence repeat (SSR) markers to investigate the influence of domestication on the genetic diversity and structure of S. japonica populations. Wild kelp populations exhibited higher genetic diversity than cultivated populations based on total NA, HE, HO, NP and AR. Discriminant analysis of principal components (DAPC), a neighbour-joining (NJ) tree and STRUCTURE analyses indicated that S. japonica populations could be divided into two groups (a cultivated/introduced group and a wild indigenous group) with significant genetic differentiation (P < 0.0001). Divergent selection, continuous inbreeding and inter-specific hybridization have caused the divergence of these two genetically separate gene pools. The significant genetic differentiation between northern and southern cultivated populations appears to be due to inter-specific hybridization and wild germplasm introduction during the domestication process. In addition, the cultivation of S. japonica has not resulted in any serious genetic disturbance of wild introduced S. japonica populations. An understanding of the genetic diversity and genetic structure of domesticated S. japonica will be necessary for further genetic improvement and effective use of germplasm.

The Ectocarpus IMMEDIATE UPRIGHT gene encodes a member of a novel family of cysteine-rich proteins with an unusual distribution across the eukaryotes

The sporophyte generation of the brown alga Ectocarpus sp. exhibits an unusual pattern of development compared with the majority of brown algae. The first cell division is symmetrical and the apical-basal axis is established late in development. In the immediate upright (imm) mutant, the initial cell undergoes an asymmetric division to immediately establish the apical-basal axis. We provide evidence which suggests that this phenotype corresponds to the ancestral state of the sporophyte. The IMM gene encodes a protein of unknown function that contains a repeated motif also found in the EsV-1-7 gene of the Ectocarpus virus EsV-1. Brown algae possess large families of EsV-1-7 domain genes but these genes are rare in other stramenopiles, suggesting that the expansion of this family might have been linked with the emergence of multicellular complexity. EsV-1-7 domain genes have a patchy distribution across eukaryotic supergroups and occur in several viral genomes, suggesting possible horizontal transfer during eukaryote evolution.

Meneco, a Topology-Based Gap-Filling Tool Applicable to Degraded Genome-Wide Metabolic Networks

Increasing amounts of sequence data are becoming available for a wide range of non-model organisms. Investigating and modelling the metabolic behaviour of those organisms is highly relevant to understand their biology and ecology. As sequences are often incomplete and poorly annotated, draft networks of their metabolism largely suffer from incompleteness. Appropriate gap-filling methods to identify and add missing reactions are therefore required to address this issue. However, current tools rely on phenotypic or taxonomic information, or are very sensitive to the stoichiometric balance of metabolic reactions, especially concerning the co-factors. This type of information is often not available or at least prone to errors for newly-explored organisms. Here we introduce Meneco, a tool dedicated to the topological gap-filling of genome-scale draft metabolic networks. Meneco reformulates gap-filling as a qualitative combinatorial optimization problem, omitting constraints raised by the stoichiometry of a metabolic network considered in other methods, and solves this problem using Answer Set Programming. Run on several artificial test sets gathering 10,800 degraded Escherichia coli networks Meneco was able to efficiently identify essential reactions missing in networks at high degradation rates, outperforming the stoichiometry-based tools in scalability. To demonstrate the utility of Meneco we applied it to two case studies. Its application to recent metabolic networks reconstructed for the brown algal model Ectocarpus siliculosus and an associated bacterium Candidatus Phaeomarinobacter ectocarpi revealed several candidate metabolic pathways for algal-bacterial interactions. Then Meneco was used to reconstruct, from transcriptomic and metabolomic data, the first metabolic network for the microalga Euglena mutabilis. These two case studies show that Meneco is a versatile tool to complete draft genome-scale metabolic networks produced from heterogeneous data, and to suggest relevant reactions that explain the metabolic capacity of a biological system.

MetaSort untangles metagenome assembly by reducing microbial community complexity

Most current approaches to analyse metagenomic data rely on reference genomes. Novel microbial communities extend far beyond the coverage of reference databases and de novo metagenome assembly from complex microbial communities remains a great challenge. Here we present a novel experimental and bioinformatic framework, metaSort, for effective construction of bacterial genomes from metagenomic samples. MetaSort provides a sorted mini-metagenome approach based on flow cytometry and single-cell sequencing methodologies, and employs new computational algorithms to efficiently recover high-quality genomes from the sorted mini-metagenome by the complementary of the original metagenome. Through extensive evaluations, we demonstrated that metaSort has an excellent and unbiased performance on genome recovery and assembly. Furthermore, we applied metaSort to an unexplored microflora colonized on the surface of marine kelp and successfully recovered 75 high-quality genomes at one time. This approach will greatly improve access to microbial genomes from complex or novel communities.

A Glimpse of Nucleo-Cytoplasmic Large DNA Virus Biodiversity through the Eukaryotic Genomics Window

The nucleocytoplasmic large DNA viruses (NCLDV) are a group of extremely complex double-stranded DNA viruses, which are major parasites of a variety of eukaryotes. Recent studies showed that certain eukaryotes contain fragments of NCLDV DNA integrated in their genome, when surprisingly many of these organisms were not previously shown to be infected by NCLDVs. We performed an update survey of NCLDV genes hidden in eukaryotic sequences to measure the incidence of this phenomenon in common public sequence databases. A total of 66 eukaryotic genomic or transcriptomic datasets-many of which are from algae and aquatic protists-contained at least one of the five most consistently conserved NCLDV core genes. Phylogenetic study of the eukaryotic NCLDV-like sequences identified putative new members of already recognized viral families, as well as members of as yet unknown viral clades. Genomic evidence suggested that most of these sequences resulted from viral DNA integrations rather than contaminating viruses. Furthermore, the nature of the inserted viral genes helped predicting original functional capacities of the donor viruses. These insights confirm that genomic insertions of NCLDV DNA are common in eukaryotes and can be exploited to delineate the contours of NCLDV biodiversity.

Identification of Genes under Positive Selection Reveals Differences in Evolutionary Adaptation between Brown-Algal Species

Brown algae are an important taxonomic group in coastal ecosystems. The model brown algal species Ectocarpus siliculosus and Saccharina japonica are closely related lineages. Despite their close phylogenetic relationship, they vary greatly in morphology and physiology. To obtain further insights into the evolutionary forces driving divergence in brown algae, we analyzed 3,909 orthologs from both species to identify Genes Under Positive Selection (GUPS). About 12% of the orthologs in each species were considered to be under positive selection. Many GUPS are involved in membrane transport, regulation of homeostasis, and sexual reproduction in the small sporophyte of E. siliculosus, which is known to have a complex life cycle and to occupy a wide range of habitats. Genes involved in photosynthesis and cell division dominated the group of GUPS in the large kelp of S. japonica, which might explain why this alga has evolved the ability to grow very rapidly and to form some of the largest sporophytes. A significant number of molecular chaperones (e.g., heat-shock proteins) involved in stress responses were identified to be under positive selection in both species, potentially indicating their important roles for macroalgae to cope with the relatively variable environment of coastal ecosystems. Moreover, analysis of previously published microarray data of E. siliculosus showed that many GUPS in E. siliculosus were responsive to stress conditions, such as oxidative and hyposaline stress, whereas our RNA-seq data of S. japonica showed that GUPS in this species were most highly expressed in large sporophytes, which supports the suggestion that selection largely acts on different sets of genes in both marcoalgal species, potentially reflecting their adaptation to different ecological niches.

A draft genome of the brown alga, Cladosiphon okamuranus, S-strain: a platform for future studies of 'mozuku' biology

The brown alga, Cladosiphon okamuranus (Okinawa mozuku), is economically one of the most important edible seaweeds, and is cultivated for market primarily in Okinawa, Japan. C. okamuranus constitutes a significant source of fucoidan, which has various physiological and biological activities. To facilitate studies of seaweed biology, we decoded the draft genome of C. okamuranus S-strain. The genome size of C. okamuranus was estimated as ∼140 Mbp, smaller than genomes of two other brown algae, Ectocarpus siliculosus and Saccharina japonica Sequencing with ∼100× coverage yielded an assembly of 541 scaffolds with N50 = 416 kbp. Together with transcriptomic data, we estimated that the C. okamuranus genome contains 13,640 protein-coding genes, approximately 94% of which have been confirmed with corresponding mRNAs. Comparisons with the E. siliculosus genome identified a set of C. okamuranus genes that encode enzymes involved in biosynthetic pathways for sulfated fucans and alginate biosynthesis. In addition, we identified C. okamuranus genes for enzymes involved in phlorotannin biosynthesis. The present decoding of the Cladosiphon okamuranus genome provides a platform for future studies of mozuku biology.

A novel codon-based de Bruijn graph algorithm for gene construction from unassembled transcriptomes

Most gene prediction methods detect coding sequences from transcriptome assemblies in the absence of closely related reference genomes. Such methods are of limited application due to high transcript fragmentation and extensive assembly errors, which may lead to redundant or false coding sequence predictions. We present inGAP-CDG, which can construct full-length and non-redundant coding sequences from unassembled transcriptomes by using a codon-based de Bruijn graph to simplify the assembly process and a machine learning-based approach to filter false positives. Compared with other methods, inGAP-CDG exhibits a significant increase in predicted coding sequence length and robustness to sequencing errors and varied read length.

Dynamics of cell wall assembly during early embryogenesis in the brown alga Fucus

Zygotes from Fucus species have been used extensively to study cell polarization and rhizoid outgrowth, and in this model system cell wall deposition aligns with the establishment of polarity. Monoclonal antibodies are essential tools for the in situ analysis of cell wall glycans, and here we report the characteristics of six monoclonal antibodies to alginates (BAM6-BAM11). The use of these, in conjunction with monoclonal antibodies to brown algal sulfated fucans, has enabled the study of the developmental dynamics of the Fucus zygote cell walls. Young zygotes are spherical and all alginate epitopes are deposited uniformly following cellulose deposition. At germination, sulfated fucans are secreted in the growing rhizoid wall. The redistribution of cell wall epitopes was investigated during treatments that cause reorientation of the growth axis (change in light direction) or disrupt rhizoid development (arabinogalactan-protein-reactive Yariv reagent). Alginate modeling was drastically impaired in the latter, and both treatments cause a redistribution of highly sulfated fucan epitopes. The dynamics of cell wall glycans in this system have been visualized in situ for the first time, leading to an enhanced understanding of the early developmental mechanisms of Fucus species. These sets of monoclonal antibodies significantly extend the available molecular tools for brown algal cell wall studies.

Sequencing of complete mitochondrial genome of brown algal Saccharina sp. ye-C6

The complete mitochondrial sequence (37 657 bp) of the Saccharina sp. ye-C6 was obtained from the assembled and annotated genome data sequenced by Illumina sequencing technology (KT271760). The circular genome contains 38 protein-coding genes (PCG), three ribosomal RNA (rRNA), and 25 transfer RNA (tRNA) genes. Base composition of genome is A (28.41%), T (36.29%), C (14.72%), and G (20.58%). Through a phylogenetic analysis based on the mitochondrial genomes of brown algae, we found that Saccharina sp. ye-C6 and Saccharina japonica are the most closely related species and strongly supports their close phylogenetic affinity.

Sequencing of complete mitochondrial genome of brown algal Saccharina sp. ye-F

The complete sequence (37 657 bp) of the mitochondrial DNA (mtDNA) of the Saccharina sp. ye-F was determined using Illumina sequencing data (Illumina Inc., San Diego, CA). The genome contains 38 protein-coding genes (PCG), three ribosomal RNA (rRNA), and 25 transfer RNA (tRNA) genes that are typical of Saccharina mtDNA. A phylogenetic analysis based on the mitochondrial genomes of brown algae indicated that Saccharina sp. ye-F and Saccharina longissima, Saccharina japonica are the most closely related species, which strongly supports their close phylogenetic affinity.

Sequencing of complete mitochondrial genome of brown algal Saccharina sp. ye-W

The complete sequence (37 657 bp) of the mitochondrial DNA (mtDNA) of the Saccharina sp. ye-W was determined using Illumina sequencing data. The genome contains 38 protein-coding genes (PCG), three ribosomal RNA (rRNA), 25 transfer RNA (tRNA) genes that are typical of Saccharina mtDNA. Phylogenetic analysis based on the mitochondrial genomes of brown algae indicated that Saccharina sp. ye-W and Saccharina longissima, Saccharina japonica are the most closely related species, which strongly supports their close phylogenetic affinity.

Comparative characterization of two GDP-mannose dehydrogenase genes from Saccharina japonica (Laminariales, Phaeophyceae)

BACKGROUND

Saccharina japonica is an important commercial brown seaweed, its main product is alginate, which is used in food, textile and by the cosmetic and pharmaceutical industries. GDP-mannose dehydrogenase (GMD) is the key enzyme involved in the synthesis of alginate. However, little is known about GMD in S. japonica. Here we report comparative biochemical analysis of two GMD genes in S. japonica.

RESULTS

Two GMD genes from S. japonica (Sjgmd1, Sjgmd2) were cloned. The open reading frame lengths of Sjgmd1, Sjgmd2 are 963 bp and 948 bp, respectively. Alignment analysis showed that the two SjGMD sequences shared 79.38 % identity. Both proteins possess the GGxCLPKDV and GxGxVG sequence motifs characteristic of the short-chain dehydrogenase/reductase superfamily. The optimum temperatures for SjGMDs were 30 °C (SjGMD1) and 20 °C (SjGMD2), and the optimum pH values were 8.0 (SjGMD1) and 8.25 (SjGMD2). Kinetic analysis demonstrated the Km values for the substrate GDP-mannose were 289 μM (SjGMD1) and 177 μM (SjGMD2), and the Km values for the cofactor NAD(+) were 139 μM (SjGMD1) and 195 μM (SjGMD2). The metal iron Zn(2+) is a potent inhibitor of SjGMD1 and SjGMD2. Real-time PCR analysis showed that heat and desiccation treatments resulted in a significant increase in Sjgmd1 and Sjgmd2 transcript abundance, suggesting that the SjGMDs are directly involved in the acclimitisation of S. japonica to abiotic stresses.

CONCLUSION

Our work identified two novel genes encoding GMD in S. japonica, comparatively characterized their structural characteristics and enzyme kinetics, and revealed the function of GMD in the stress adaptability of S. japonica. The knowledge obtained here enriched our understanding of the alginate synthesis mechanism in S. japonica, and may promote further research on functional differences between GMD genes.

A General Strategy for Facile Synthesis and In Situ Screening of Self-Assembled Polymer-Peptide Nanomaterials

A universal strategy for efficient, mild, and purification-free synthesis and in situ screening of functional polymer-peptide nanomaterials is described. More than 1000 polymer-peptide conjugates (PPCs) with various chemical structures, compositions, and therapeutic efficacy are created. According to this strategy, the structure-function relationship of the PPCs is revealed, and the antitumor efficacies of the top performing PPCs are evaluated in vivo.

Improving seedless kelp (Saccharina japonica) during its domestication by hybridizing gametophytes and seedling-raising from sporophytes

Dongfang no.7 (Saccharina japonica) was bred and maintained by hybridizing gametophytes, self-crossing the best individuals, selecting the best self-crossing line and seedling-raising from yearly reconstructed sporophytes. It increased the air dry yield by 43.2% in average over 2 widely farmed controls. Dongfang no.7 was seedling-raised from bulked sporophytes reconstructed from its representative gametophyte clones. Such strategy ensured it against variety contamination due to possible cross fertilization and occasional mixing and inbred depletion due to self-crossing number-limited sporophytes year after year. It derived from an intraspecific hybrid through 4 rounds of self-crossing and selection and retained a certain degree of genetic heterozygosity, thus being immune to inbred depletion due to purification of unknown detrimental alleles. Most importantly, it can be farmed in currently available system as the seedlings for large scale culture can be raised from reconstructed Dongfang no.7 sporophytes. Breeding and maintaining Dongfang no.7 provided a model that other varieties of kelp (S. japonica) and brown algae may follow during their domestication.

Sequencing of complete mitochondrial genome of brown algal Saccharina sp. ye-C5

The complete genome (37 635 bp) mitochondrial DNA (mtDNA) of the Saccharina sp. ye-C5 was determined, which contains 38 protein-coding genes (PCG), three ribosomal RNA (rRNA) and 25 transfer RNA (tRNA) genes that are consistent in Saccharina genus. The phylogenetic tree that was established based on the mitochondrial genomes of brown algae, which indicated that Saccharina sp. ye-C5 and Saccharina longissima are the most closely related species.

Sequencing of complete mitochondrial genome of brown algal Saccharina sp. ye-C2

The complete sequence (37 657 bp) of the mitochondrial DNA (mtDNA) of the Saccharina sp. ye-C2 was determined. About 38 protein-coding genes (PCG), 3 ribosomal RNAs (rRNA) and 25 transfer RNA (tRNA) genes were annotated in the genome. The phylogenetic analysis strongly supports the close phylogenetic affinity of Saccharina sp. ye-C2 and Saccharina japonica based on the mitochondrial genomes of other brown algae.

Sequencing of complete mitochondrial genome of brown algal Saccharina sp. ye- C12

The complete sequence (37 654 bp) of the mitochondrial DNA (mtDNA) of the Saccharina sp. ye-C12 was determined using Illumina sequencing data. The genome contains 38 protein-coding genes (PCG), three ribosomal RNA (rRNA), 25 transfer RNA (tRNA) genes and the phylogenetic analysis based on the mitochondrial genomes of brown algae indicated that Saccharina sp. ye-C12 and Saccharina japonica are the most closely related species, which strongly supports their close phylogenetic affinity.

Construction of a high-density genetic map and mapping of a sex-linked locus for the brown alga Undaria pinnatifida (Phaeophyceae) based on large scale marker development by specific length amplified fragment (SLAF) sequencing

BACKGROUND

Undaria pinnatifida is an important economic brown alga in East Asian countries. However, its genetic and genomic information is very scarce, which hinders further research in this species. A high-density genetic map is a basic tool for fundamental and applied research such as discovery of functional genes and mapping of quantitative trait loci (QTL). In this study the recently developed specific length amplified fragment sequencing (SLAF-seq) technology was employed to construct a high-density genetic linkage map and locate a sex determining locus for U. pinnatifida.

RESULTS

A total of 28.06 Gb data including 140.31 M pair-end reads was obtained. After linkage analysis 4626 SLAF markers were mapped onto the genetic map. After adding the sex linked simple sequence repeat (SSR) marker [GenBank:AY738602.1], the final genetic map was 1816.28 cM long, consisting of 30 linkage groups with an average distance of 0.39 cM between adjacent markers. The length of LGs ranged from 20.12 to 106.95 cM. A major sex associated QTL was mapped to LG22 within a window starting at 29.01 cM and ending at 68.81 cM with a total of 68 SLAF markers. The SSR marker and five SLAF markers (Marker6556, 19020, 43089, 60771 and 26359) were identified as tightly sex-linked markers, as indicated by the absence of recombination between them and the sex phenotype. These markers were located at the position of 59.50 cM, which was supposed to be the sex determining region.

CONCLUSIONS

A high-density genetic linkage map was constructed using SLAF-seq technique and F1 gametophyte population for the first time in the economically important brown alga U. pinnatifida. For the first time, a major sex associated QTL suggesting a sex determining region was mapped to a single LG. This map will facilitate the further fundamental and applied research such as QTL mapping and map-based gene clone in U. pinnatifida and provide a reference for studies in other kelp species.

The Pseudoautosomal Regions of the U/V Sex Chromosomes of the Brown Alga Ectocarpus Exhibit Unusual Features

The recombining regions of sex chromosomes (pseudoautosomal regions, PARs) are predicted to exhibit unusual features due to their being genetically linked to the nonrecombining, sex-determining region. This phenomenon is expected to occur in both diploid (XY, ZW) and haploid (UV) sexual systems, with slightly different consequences for UV sexual systems because of the absence of masking during the haploid phase (when sex is expressed) and because there is no homozygous sex in these systems. Despite a considerable amount of theoretical work on PAR genetics and evolution, these genomic regions have remained poorly characterized empirically. We show here that although the PARs of the U/V sex chromosomes of the brown alga Ectocarpus recombine at a similar rate to autosomal regions of the genome, they exhibit many genomic features typical of nonrecombining regions. The PARs were enriched in clusters of genes that are preferentially, and often exclusively, expressed during the sporophyte generation of the life cycle, and many of these genes appear to have evolved since the Ectocarpales diverged from other brown algal lineages. A modeling-based approach was used to investigate possible evolutionary mechanisms underlying this enrichment in sporophyte-biased genes. Our results are consistent with the evolution of the PAR in haploid systems being influenced by differential selection pressures in males and females acting on alleles that are advantageous during the sporophyte generation of the life cycle.

Recent Major Advances of Biotechnology and Sustainable Aquaculture in China

Background:

Global aquaculture production has increased continuously over the last five decades, and particularly in China. Its aquaculture has become the fastest growing and most efficient agri-sector, with production accounting for more than 70% of the world’s aquaculture output. In the new century, with serious challenges regarding population, resources and the environment, China has been working to develop high-quality, effective, healthy, and sustainable blue agriculture through the application of modern biotechnology. Sound knowledge related to the biology and ecology of aquatic organisms has laid a solid foundation and provided the innovation and technology for rapid development of the aquaculture industry. Marine biotechnology, which is enabling solutions for ocean productivity and sustainability, has been promoted since the last decades of the 20th Century in China.

Objective:

In this article, priority areas of research, mainly genetic breeding, omics studies, novel production systems, biosecurity, bioprocesses and biorefinery, as well as the major progress of marine biotechnology R&D in China are reviewed.

Conclusion:

Current innovative achievements in China are not enough and the level and frequency of academic advancements must be improved. International cooperation and assistance remain crucial for the success of marine biotechnology.