Search

Coffea arabica, an allotetraploid hybrid of Coffea eugenioidesand Coffea canephora, is the source of approximately 60% of coffee products worldwide, and its cultivated accessions have undergone several population bottlenecks. We present chromosome-level assemblies of a di-haploid C. arabicaaccession and modern representatives of its diploid progenitors, C. eugenioidesand C. canephora. The three species exhibit largely conserved genome structures between diploid parents and descendant subgenomes, with no obvious global subgenome dominance. We find evidence for a founding polyploidy event 350,000–610,000 years ago, followed by several pre-domestication bottlenecks, resulting in narrow genetic variation. A split between wild accessions and cultivar progenitors occurred ~30.5 thousand years ago, followed by a period of migration between the two populations. Analysis of modern varieties, including lines historically introgressed with C. canephora, highlights their breeding histories and loci that may contribute to pathogen resistance, laying the groundwork for future genomics-based breeding of C. arabica.

We present chromosome-level genome assemblies from representative species of three independently evolved seagrass lineages: Posidonia oceanica, Cymodocea nodosa, Thalassia testudinumand Zostera marina. We also include a draft genome of Potamogeton acutifolius, belonging to a freshwater sister lineage to Zosteraceae. All seagrass species share an ancient whole-genome triplication, while additional whole-genome duplications were uncovered for C. nodosa, Z. marinaand P. acutifolius. Comparative analysis of selected gene families suggests that the transition from submerged-freshwater to submerged-marine environments mainly involved fine-tuning of multiple processes (such as osmoregulation, salinity, light capture, carbon acquisition and temperature) that all had to happen in parallel, probably explaining why adaptation to a marine lifestyle has been exceedingly rare. Major gene losses related to stomata, volatiles, defence and lignification are probably a consequence of the return to the sea rather than the cause of it. These new genomes will accelerate functional studies and solutions, as continuing losses of the ‘savannahs of the sea’ are of major concern in times of climate change and loss of biodiversity.

ABSTRACTDiscus (Symphysodon spp.) is a popular freshwater ornamental fish with high market value due to its extraordinary skin pigmentation. The molecular mechanisms that govern skin colouration have yet to be identified. Six skin samples from wild discus were collected in the D1 (dorsal intense black bar stripe), D2 (ventral lighter bar stripes and orange carotenoid spots) and D3 (posterior without stripe nor spot) regions of wild discus for transcriptomic profiling and analysis. Differential analysis was conducted in 3 pair-wise comparisons, D1vsD3, D2vsD3 and D2vsD1. The results demonstrated that the melanin synthesis-associated genes (ERBB4, ADCY9, WNT2, 14-3-3proteins and TSPANgenes) were upregulated in D1 and D2 regions, whereas carotenoids metabolism-associated genes (PLIN2, AHCY, RETSAT, STARD10, TTC39Band RDH1genes) were uniquely expressed in D2 region. Amino acid metabolism-related pathways were enriched in both the D1 and D2 regions. MAPK pathway enrichment was detected merely in the D1 region, while various transporter activities were enriched in the D2 region. Further investigations speculated the key regulator roles of GNAQ, ADCY9, 14-3-3and TSPAN genes in black vertical bar formation. This is the first fundamental report on discus skin pigment genes that could be extended into more detailed studies in the future.

Currents are unique drivers of oceanic phylogeography and thus determine the distribution of marine coastal species, along with past glaciations and sea-level changes. Here we reconstruct the worldwide colonization history of eelgrass (Zostera marinaL.), the most widely distributed marine flowering plant or seagrass from its origin in the Northwest Pacific, based on nuclear and chloroplast genomes. We identified two divergent Pacific clades with evidence for admixture along the East Pacific coast. Two west-to-east (trans-Pacific) colonization events support the key role of the North Pacific Current. Time-calibrated nuclear and chloroplast phylogenies yielded concordant estimates of the arrival of Z. marinain the Atlantic through the Canadian Arctic, suggesting that eelgrass-based ecosystems, hotspots of biodiversity and carbon sequestration, have only been present there for ~243 ky (thousand years). Mediterranean populations were founded ~44 kya, while extant distributions along western and eastern Atlantic shores were founded at the end of the Last Glacial Maximum (~19 kya), with at least one major refuge being the North Carolina region. The recent colonization and five- to sevenfold lower genomic diversity of the Atlantic compared to the Pacific populations raises concern and opportunity about how Atlantic eelgrass might respond to rapidly warming coastal oceans.

Most fresh bananas belong to the Cavendish and Gros Michel subgroups. Here, we report chromosome-scale genome assemblies of Cavendish (1.48 Gb) and Gros Michel (1.33 Gb), defining three subgenomes, Ban, Dh and Ze, with Musa acuminatassp. banksii, malaccensisand zebrinaas their major ancestral contributors, respectively. The insertion of repeat sequences in the Fusarium oxysporumf. sp. cubense(Foc) tropical race 4 RGA2(resistance gene analog 2) promoter was identified in most diploid and triploid bananas. We found that the receptor-like protein (RLP) locus, including Focrace 1-resistant genes, is absent in the Gros Michel Ze subgenome. We identified two NAP (NAC-like, activated by apetala3/pistillata) transcription factor homologs specifically and highly expressed in fruit that directly bind to the promoters of many fruit ripening genes and may be key regulators of fruit ripening. Our genome data should facilitate the breeding and super-domestication of bananas.

To improve our understanding of the origin and evolution of mycoheterotrophic plants, we here present the chromosome-scale genome assemblies of two sibling orchid species: partially mycoheterotrophic Platanthera zijinensisand holomycoheterotrophic Platanthera guangdongensis. Comparative analysis shows that mycoheterotrophy is associated with increased substitution rates and gene loss, and the deletion of most photoreceptor genes and auxin transporter genes might be linked to the unique phenotypes of fully mycoheterotrophic orchids. Conversely, trehalase genes that catalyse the conversion of trehalose into glucose have expanded in most sequenced orchids, in line with the fact that the germination of orchid non-endosperm seeds needs carbohydrates from fungi during the protocorm stage. We further show that the mature plant of P. guangdongensis, different from photosynthetic orchids, keeps expressing trehalase genes to hijack trehalose from fungi. Therefore, we propose that mycoheterotrophy in mature orchids is a continuation of the protocorm stage by sustaining the expression of trehalase genes. Our results shed light on the molecular mechanism underlying initial, partial and full mycoheterotrophy.

Cycads represent one of the most ancient lineages of living seed plants. Identifying genomic features uniquely shared by cycads and other extant seed plants, but not non-seed-producing plants, may shed light on the origin of key innovations, as well as the early diversification of seed plants. Here, we report the 10.5-Gb reference genome of Cycas panzhihuaensis, complemented by the transcriptomes of 339 cycad species. Nuclear and plastid phylogenomic analyses strongly suggest that cycads and Ginkgoform a clade sister to all other living gymnosperms, in contrast to mitochondrial data, which place cycads alone in this position. We found evidence for an ancient whole-genome duplication in the common ancestor of extant gymnosperms. The Cycasgenome contains four homologues of the fitDgene family that were likely acquired via horizontal gene transfer from fungi, and these genes confer herbivore resistance in cycads. The male-specific region of the Y chromosome of C. panzhihuaensiscontains a MADS-box transcription factor expressed exclusively in male cones that is similar to a system reported in Ginkgo, suggesting that a sex determination mechanism controlled by MADS-box genes may have originated in the common ancestor of cycads and Ginkgo. The C. panzhihuaensisgenome provides an important new resource of broad utility for biologists.

Euphyllophytes encompass almost all extant plants, including two sister clades, ferns and seed plants. Decoding genomes of ferns is the key to deep insight into the origin of euphyllophytes and the evolution of seed plants. Here we report a chromosome-level genome assembly of Adiantum capillus-venerisL., a model homosporous fern. This fern genome comprises 30 pseudochromosomes with a size of 4.8-gigabase and a contig N50 length of 16.22 Mb. Gene co-expression network analysis uncovered that homospore development in ferns has relatively high genetic similarities with that of the pollen in seed plants. Analysing fern defence response expands understanding of evolution and diversity in endogenous bioactive jasmonates in plants. Moreover, comparing fern genomes with those of other land plants reveals changes in gene families important for the evolutionary novelties within the euphyllophyte clade. These results lay a foundation for studies on fern genome evolution and function, as well as the origin and evolution of euphyllophytes.

Understanding how organisms adapt to extreme living conditions is central to evolutionary biology. Dark septate endophytes (DSEs) constitute an important component of the root mycobiome and they are often able to alleviate host abiotic stresses. Here, we investigated the molecular mechanisms underlying the beneficial association between the DSE Laburnicola rhizohalophilaand its host, the native halophyte Suaeda salsa, using population genomics. Based on genome-wide Fst (pairwise fixation index) and Vst analyses, which compared the variance in allele frequencies of single-nucleotide polymorphisms (SNPs) and copy number variants (CNVs), respectively, we found a high level of genetic differentiation between two populations. CNV patterns revealed population-specific expansions and contractions. Interestingly, we identified a ~20?kbp genomic island of high divergence with a strong sign of positive selection. This region contains a melanin-biosynthetic polyketide synthase gene cluster linked to six additional genes likely involved in biosynthesis, membrane trafficking, regulation, and localization of melanin. Differences in growth yield and melanin biosynthesis between the two populations grown under 2% NaCl stress suggested that this genomic island contributes to the observed differences in melanin accumulation. Our findings provide a better understanding of the genetic and evolutionary mechanisms underlying the adaptation to saline conditions of the L. rhizohalophila–S. salsasymbiosis.

Understanding how organisms adapt to extreme living conditions is central to evolutionary biology. Dark septate endophytes (DSEs) constitute an important component of the root mycobiome and they are often able to alleviate host abiotic stresses. Here, we investigated the molecular mechanisms underlying the beneficial association between the DSE Laburnicola rhizohalophilaand its host, the native halophyte Suaeda salsa, using population genomics. Based on genome-wide Fst (pairwise fixation index) and Vst analyses, which compared the variance in allele frequencies of single-nucleotide polymorphisms (SNPs) and copy number variants (CNVs), respectively, we found a high level of genetic differentiation between two populations. CNV patterns revealed population-specific expansions and contractions. Interestingly, we identified a ~20 kbp genomic island of high divergence with a strong sign of positive selection. This region contains a melanin-biosynthetic polyketide synthase gene cluster linked to six additional genes likely involved in biosynthesis, membrane trafficking, regulation, and localization of melanin. Differences in growth yield and melanin biosynthesis between the two populations grown under 2% NaCl stress suggested that this genomic island contributes to the observed differences in melanin accumulation. Our findings provide a better understanding of the genetic and evolutionary mechanisms underlying the adaptation to saline conditions of the L. rhizohalophila–S. salsasymbiosis.

Mounting evidence suggests that terrestrialization of plants started in streptophyte green algae, favoured by their dual existence in freshwater and subaerial/terrestrial environments. Here, we present the genomes of Mesostigma virideand Chlorokybus atmophyticus, two sister taxa in the earliest-diverging clade of streptophyte algae dwelling in freshwater and subaerial/terrestrial environments, respectively. We provide evidence that the common ancestor of M. virideand C. atmophyticus(and thus of streptophytes) had already developed traits associated with a subaerial/terrestrial environment, such as embryophyte-type photorespiration, canonical plant phytochrome, several phytohormones and transcription factors involved in responses to environmental stresses, and evolution of cellulose synthase and cellulose synthase-like genes characteristic of embryophytes. Both genomes differed markedly in genome size and structure, and in gene family composition, revealing their dynamic nature, presumably in response to adaptations to their contrasting environments. The ancestor of M. viridepossibly lost several genomic traits associated with a subaerial/terrestrial environment following transition to a freshwater habitat.

Hornworts, liverworts and mosses are three early diverging clades of land plants, and together comprise the bryophytes. Here, we report the draft genome sequence of the hornwort Anthoceros angustus. Phylogenomic inferences confirm the monophyly of bryophytes, with hornworts sister to liverworts and mosses. The simple morphology of hornworts correlates with low genetic redundancy in plant body plan, while the basic transcriptional regulation toolkit for plant development has already been established in this early land plant lineage. Although the Anthocerosgenome is small and characterized by minimal redundancy, expansions are observed in gene families related to RNA editing, UV protection and desiccation tolerance. The genome of A. angustusbears the signatures of horizontally transferred genes from bacteria and fungi, in particular of genes operating in stress-response and metabolic pathways. Our study provides insight into the unique features of hornworts and their molecular adaptations to live on land.

Water lilies belong to the angiosperm order Nymphaeales. Amborellales, Nymphaeales and Austrobaileyales together form the so-called ANA-grade of angiosperms, which are extant representatives of lineages that diverged the earliest from the lineage leading to the extant mesangiosperms1–3. Here we report the 409-megabase genome sequence of the blue-petal water lily (Nymphaea colorata). Our phylogenomic analyses support Amborellales and Nymphaeales as successive sister lineages to all other extant angiosperms. The N. coloratagenome and 19 other water lily transcriptomes reveal a Nymphaealean whole-genome duplication event, which is shared by Nymphaeaceae and possibly Cabombaceae. Among the genes retained from this whole-genome duplication are homologues of genes that regulate flowering transition and flower development. The broad expression of homologues of floral ABCE genes in N. coloratamight support a similarly broadly active ancestral ABCE model of floral organ determination in early angiosperms. Water lilies have evolved attractive floral scents and colours, which are features shared with mesangiosperms, and we identified their putative biosynthetic genes in N. colorata. The chemical compounds and biosynthetic genes behind floral scents suggest that they have evolved in parallel to those in mesangiosperms. Because of its unique phylogenetic position, the N. coloratagenome sheds light on the early evolution of angiosperms.

Starmerella (Candida) bombicolais the biosurfactant-producing species that caught the greatest deal of attention in the academic and industrial world due to its ability of producing large amounts of sophorolipids. Despite its high economic potential, the biochemistry behind the sophorolipid biosynthesis is still poorly understood. Here we present the first proteomic characterization of S. bombicolafor which we created a lys1Δ mutant to allow the use of SILAC for quantitative analysis. To characterize the processes behind the production of these biosurfactants, we compared the proteome of sophorolipid producing (early stationary phase) and nonproducing cells (exponential phase). We report the simultaneous production of all known enzymes involved in sophorolipid biosynthesis including a predicted sophorolipid transporter. In addition, we identified the heme binding protein Dap1 as a possible regulator for Cyp52M1. Our results further indicate that ammonium and phosphate limitation are not the sole limiting factors inducing sophorolipid biosynthesis.

MicroTom has a short growth cycle and high transformation efficiency, and is a prospective model plant for studying organ development, metabolism, and plant–microbe interactions. Here, with a newly assembled reference genome for this tomato cultivar and abundant RNA-seq data derived from tissues of different organs/developmental stages/treatments, we constructed multiple gene co-expression networks, which will provide valuable clues for the identification of important genes involved in diverse regulatory pathways during plant growth, e.g. arbuscular mycorrhizal symbiosis and fruit development. Additionally, non-coding RNAs, including miRNAs, lncRNAs, and circRNAs were also identified, together with their potential targets. Interacting networks between different types of non-coding RNAs (miRNA-lncRNA), and non-coding RNAs and genes (miRNA-mRNA and lncRNA-mRNA) were constructed as well. Our results and data will provide valuable information for the study of organ differentiation and development of this important fruit. Lastly, we established a database (http://eplant.njau.edu.cn/microTomBase/) with genomic and transcriptomic data, as well as details of gene co-expression and interacting networks on MicroTom, and this database should be of great value to those who want to adopt MicroTom as a model plant for research.

Ferns are the closest sister group to all seed plants, yet little is known about their genomes other than that they are generally colossal. Here, we report on the genomes of Azolla filiculoidesand Salvinia cucullata(Salviniales) and present evidence for episodic whole-genome duplication in ferns—one at the base of ‘core leptosporangiates’ and one specific to Azolla. One fern-specific gene that we identified, recently shown to confer high insect resistance, seems to have been derived from bacteria through horizontal gene transfer. Azollacoexists in a unique symbiosis with N2-fixing cyanobacteria, and we demonstrate a clear pattern of cospeciation between the two partners. Furthermore, the Azollagenome lacks genes that are common to arbuscular mycorrhizal and root nodule symbioses, and we identify several putative transporter genes specific to Azolla–cyanobacterial symbiosis. These genomic resources will help in exploring the biotechnological potential of Azollaand address fundamental questions in the evolution of plant life.

Although classical evolutionary theory, i.e., population genetics and the Modern Synthesis, was already implicitly ‘gene-centred’, the organism was, in practice, still generally regarded as the individual unit of which a population is composed. The gene-centred approach to evolution only reached a logical conclusion with the advent of the gene-selectionist or gene’s eye view in the 1960s and 1970s. Whereas classical evolutionary theory can only work with (genotypically represented) fitness differences between individual organisms, gene-selectionism is capable of working with fitness differences among genes within the same organism and genome. Here, we explore the explanatory potential of ‘intra-organismic’ and ‘intra-genomic’ gene-selectionism, i.e., of a behavioural-ecological ‘gene’s eye view’ on genetic, genomic and organismal evolution. First, we give a general outline of the framework and how it complements the—to some extent—still ‘organism-centred’ approach of classical evolutionary theory. Secondly, we give a more in-depth assessment of its explanatory potential for biological evolution, i.e., for Darwin’s ‘common descent with modification’ or, more specifically, for ‘historical continuity or homology with modular evolutionary change’ as it has been studied by evolutionary developmental biology (evo-devo) during the last few decades. In contrast with classical evolutionary theory, evo-devo focuses on ‘within-organism’ developmental processes. Given the capacity of gene-selectionism to adopt an intra-organismal gene’s eye view, we outline the relevance of the latter model for evo-devo. Overall, we aim for the conceptual integration between the gene’s eye view on the one hand, and more organism-centred evolutionary models (both classical evolutionary theory and evo-devo) on the other.

Gnetophytes are an enigmatic gymnosperm lineage comprising three genera, Gnetum, Welwitschiaand Ephedra, which are morphologically distinct from all other seed plants. Their distinctiveness has triggered much debate as to their origin, evolution and phylogenetic placement among seed plants. To increase our understanding of the evolution of gnetophytes, and their relation to other seed plants, we report here a high-quality draft genome sequence for Gnetum montanum, the first for any gnetophyte. By using a novel genome assembly strategy to deal with high levels of heterozygosity, we assembled >4 Gb of sequence encoding 27,491 protein-coding genes. Comparative analysis of the G. montanumgenome with other gymnosperm genomes unveiled some remarkable and distinctive genomic features, such as a diverse assemblage of retrotransposons with evidence for elevated frequencies of elimination rather than accumulation, considerable differences in intron architecture, including both length distribution and proportions of (retro) transposon elements, and distinctive patterns of proliferation of functional protein domains. Furthermore, a few gene families showed Gnetum-specific copy number expansions (for example, cellulose synthase) or contractions (for example, Late Embryogenesis Abundant protein), which could be connected with Gnetum’s distinctive morphological innovations associated with their adaptation to warm, mesic environments. Overall, the G. montanumgenome enables a better resolution of ancestral genomic features within seed plants, and the identification of genomic characters that distinguish Gnetumfrom other gymnosperms. The evolution and phylogenetic placement of gnetophytes have remained elusive. The draft genome for a gnetophyte, Gnetum montanum, is now reported. Comparative analyses identify genomic features distinguishing gnetophytes from other gymnosperms.

The genomes of Cercospora zeina, Fusarium pininemorale, Hawksworthiomyces lignivorus, Huntiella decipiens, and Ophiostoma ips arepresented in this genome announcement. Three of these genomes are from plant pathogens and otherwise economically important fungal species. Fusarium pininemoraleand H. decipiensare not known to cause significant disease but are closely related to species of economic importance. The genome sizes range from 25.99 Mb in the case of O. ipsto 4.82 Mb for H. lignivorus.These genomes include the first reports of a genome from the genus Hawksworthiomyces.The availability of these genome data will allow the resolution of longstanding questions regarding the taxonomy of these species. In addition these genome sequences through comparative studies with closely related organisms will increase our understanding of how these species or close relatives cause disease.

WebComparing the whole genome sequence of Apostasia shenzhenicawith transcriptome and genome data from five orchid subfamilies permits the reconstruction of an ancestral gene toolkit, providing insight into orchid origins, evolution and diversification.

Polyploidy, or the duplication of entire genomes, has been observed in prokaryotic and eukaryotic organisms, and in somatic and germ cells. The consequences of polyploidization are complex and variable, and they differ greatly between systems (clonal or non-clonal) and species, but the process has often been considered to be an evolutionary 'dead end'. Here, we review the accumulating evidence that correlates polyploidization with environmental change or stress, and that has led to an increased recognition of its short-term adaptive potential. In addition, we discuss how, once polyploidy has been established, the unique retention profile of duplicated genes following whole-genome duplication might explain key longer-term evolutionary transitions and a general increase in biological complexity.

ABSTRACTMost of the proteins that are specifically turned over by selective autophagy are recognized by the presence of short Atg8 interacting motifs (AIMs) that facilitate their association with the autophagy apparatus. Such AIMs can be identified by bioinformatics methods based on their defined degenerate consensus F/W/Y-X-X-L/I/V sequences in which X represents any amino acid. Achieving reliability and/or fidelity of the prediction of such AIMs on a genome-wide scale represents a major challenge. Here, we present a bioinformatics approach, high fidelity AIM (hfAIM), which uses additional sequence requirements—the presence of acidic amino acids and the absence of positively charged amino acids in certain positions—to reliably identify AIMs in proteins. We demonstrate that the use of the hfAIM method allows for in silico high fidelity prediction of AIMs in AIM-containing proteins (ACPs) on a genome-wide scale in various organisms. Furthermore, by using hfAIM to identify putative AIMs in the Arabidopsisproteome, we illustrate a potential contribution of selective autophagy to various biological processes. More specifically, we identified 9 peroxisomal PEX proteins that contain hfAIM motifs, among which AtPEX1, AtPEX6 and AtPEX10 possess evolutionary-conserved AIMs. Bimolecular fluorescence complementation (BiFC) results verified that AtPEX6 and AtPEX10 indeed interact with Atg8 in planta. In addition, we show that mutations occurring within or nearby hfAIMs in PEX1, PEX6 and PEX10 caused defects in the growth and development of various organisms. Taken together, the above results suggest that the hfAIM tool can be used to effectively perform genome-wide in silico screens of proteins that are potentially regulated by selective autophagy. The hfAIM system is a web tool that can be accessed at link: http://bioinformatics.psb.ugent.be/hfAIM/.

Seagrasses colonized the sea on at least three independent occasions to form the basis of one of the most productive and widespread coastal ecosystems on the planet. Here we report the genome of Zostera marina (L.), the first, to our knowledge, marine angiosperm to be fully sequenced. This reveals unique insights into the genomic losses and gains involved in achieving the structural and physiological adaptations required for its marine lifestyle, arguably the most severe habitat shift ever accomplished by flowering plants. Key angiosperm innovations that were lost include the entire repertoire of stomatal genes, genes involved in the synthesis of terpenoids and ethylene signalling, and genes for ultraviolet protection and phytochromes for far-red sensing. Seagrasses have also regained functions enabling them to adjust to full salinity. Their cell walls contain all of the polysaccharides typical of land plants, but also contain polyanionic, low-methylated pectins and sulfated galactans, a feature shared with the cell walls of all macroalgae and that is important for ion homoeostasis, nutrient uptake and O2/CO2exchange through leaf epidermal cells. The Z. marina genome resource will markedly advance a wide range of functional ecological studies from adaptation of marine ecosystems under climate warming, to unravelling the mechanisms of osmoregulation under high salinities that may further inform our understanding of the evolution of salt tolerance in crop plants.

The genome of a magnoliid, Aristolochia fimbriata, shows no evidence of whole-genome duplication, a feature uniquely shared with Amborella trichopoda, the sister species to all other angiosperms. This ancestral structure of the Aristolochiagenome offers great opportunities for comparative genomics.

Hibiscus hamabois a semi-mangrove species with strong tolerance to salt and waterlogging stress. However, the molecular basis and mechanisms that underlie this strong adaptability to harsh environments remain poorly understood. Here, we assembled a high-quality, chromosome-level genome of this semi-mangrove plant and analyzed its transcriptome under different stress treatments to reveal regulatory responses and mechanisms. Our analyses suggested that H. hamabohas undergone two recent successive polyploidy events, a whole-genome duplication followed by a whole-genome triplication, resulting in an unusually large gene number (107 309 genes). Comparison of the H. hamabogenome with that of its close relative Hibiscus cannabinus, which has not experienced a recent WGT, indicated that genes associated with high stress resistance have been preferentially preserved in the H. hamabogenome, suggesting an underlying association between polyploidy and stronger stress resistance. Transcriptomic data indicated that genes in the roots and leaves responded differently to stress. In roots, genes that regulate ion channels involved in biosynthetic and metabolic processes responded quickly to adjust the ion concentration and provide metabolic products to protect root cells, whereas no such rapid response was observed from genes in leaves. Using co-expression networks, potential stress resistance genes were identified for use in future functional investigations. The genome sequence, along with several transcriptome datasets, provide insights into genome evolution and the mechanism of salt and waterlogging tolerance in H. hamabo, suggesting the importance of polyploidization for environmental adaptation.

Ginger (Zingiber officinale) is one of the most valued spice plants worldwide; it is prized for its culinary and folk medicinal applications and is therefore of high economic and cultural importance. Here, we present a haplotype-resolved, chromosome-scale assembly for diploid ginger anchored to 11 pseudochromosome pairs with a total length of 3.1 Gb. Remarkable structural variation was identified between haplotypes, and two inversions larger than 15 Mb on chromosome 4 may be associated with ginger infertility. We performed a comprehensive, spatiotemporal, genome-wide analysis of allelic expression patterns, revealing that most alleles are coordinately expressed. The alleles that exhibited the largest differences in expression showed closer proximity to transposable elements, greater coding sequence divergence, more relaxed selection pressure, and more transcription factor binding site differences. We also predicted the transcription factors potentially regulating 6-gingerol biosynthesis. Our allele-aware assembly provides a powerful platform for future functional genomics, molecular breeding, and genome editing in ginger.

Polyploidization plays a key role in plant evolution, but the forces driving the fate of homoeologs in polyploid genomes, i.e., paralogs resulting from a whole-genome duplication (WGD) event, remain to be elucidated. Here, we present a chromosome-scale genome assembly of tetraploid scarlet sage (Salvia splendens), one of the most diverse ornamental plants. We found evidence for three WGD events following an older WGD event shared by most eudicots (the γ event). A comprehensive, spatiotemporal, genome-wide analysis of homoeologs from the most recent WGD unveiled expression asymmetries, which could be associated with genomic rearrangements, transposable element proximity discrepancies, coding sequence variation, selection pressure, and transcription factor binding site differences. The observed differences between homoeologs may reflect the first step toward sub- and/or neofunctionalization. This assembly provides a powerful tool for understanding WGD and gene and genome evolution and is useful in developing functional genomics and genetic engineering strategies for scarlet sage and other Lamiaceae species.

The marvelously diverse Orchidaceae constitutes the largest family of angiosperms. The genus Cymbidiumin Orchidaceae is well known for its unique vegetation, floral morphology, and flower scent traits. Here, a chromosome-scale assembly of the genome of Cymbidium ensifolium(Jianlan) is presented. Comparative genomic analysis showed that C. ensifoliumhas experienced two whole-genome duplication (WGD) events, the most recent of which was shared by all orchids, while the older event was the τ event shared by most monocots. The results of MADS-box genes analysis provided support for establishing a unique gene model of orchid flower development regulation, and flower shape mutations in C. ensifoliumwere shown to be associated with the abnormal expression of MADS-box genes. The most abundant floral scent components identified included methyl jasmonate, acacia alcohol and linalool, and the genes involved in the floral scent component network of C. ensifoliumwere determined. Furthermore, the decreased expression of photosynthesis-antennae and photosynthesis metabolic pathway genes in leaves was shown to result in colorful striped leaves, while the increased expression of MADS-box genes in leaves led to perianth-like leaves. Our results provide fundamental insights into orchid evolution and diversification.

Inferring comprehensive regulatory networks from high-throughput data is one of the foremost challenges of modern computational biology. As high-throughput expression profiling experiments have gained common ground in many laboratories, different techniques have been proposed to infer transcriptional regulatory networks from them. Furthermore, with the advent of diverse types of high-throughput data, the research in network inference has received a new impulse. The use of diverse types of data, together with the increasing tendency of building the inference on biologically plausible simplifications, allows a more reliable and more complete description of networks. Here, we discuss how the research focus in the field of network inference is increasingly shifting from methods trying to reconstruct networks from a single data type towards integrative approaches dealing with several data sources simultaneously to infer regulatory modules.

Inferring comprehensive regulatory networks from high-throughput data is one of the foremost challenges of modern computational biology. As high-throughput expression profiling experiments have gained common ground in many laboratories, different techniques have been proposed to infer transcriptional regulatory networks from them. Furthermore, with the advent of diverse types of high-throughput data, the research in network inference has received a new impulse. The use of diverse types of data, together with the increasing tendency of building the inference on biologically plausible simplifications, allows a more reliable and more complete description of networks. Here, we discuss how the research focus in the field of network inference is increasingly shifting from methods trying to reconstruct networks from a single data type towards integrative approaches dealing with several data sources simultaneously to infer regulatory modules.

It is expected that one of the merits of comparative genomics lies in the transfer of structural and functional information from one genome to another. This is based on the observation that, although the number of chromosomal rearrangements that occur in genomes is extensive, different species still exhibit a certain degree of conservation regarding gene content and gene order. It is in this respect that we have developed a new software tool for the Automatic Detection of Homologous Regions (ADHoRe). ADHoRe was primarily developed to find large regions of microcolinearity, taking into account different types of microrearrangements such as tandem duplications, gene loss and translocations, and inversions. Such rearrangements often complicate the detection of colinearity, in particular when comparing more anciently diverged species. Application of ADHoRe to the complete genome of Arabidopsis and a large collection of concatenated rice BACs yields more than 20 regions showing statistically significant microcolinearity between both plant species. These regions comprise from 4 up to 11 conserved homologous gene pairs. We predict the number of homologous regions and the extent of microcolinearity to increase significantly once better annotations of the rice genome become available.

Motivation: The large amount of available annotated Arabidopsis thaliana sequences allows the induction of splice site prediction models with supervised learning algorithms (see Haussler (1998) for a review and references). These algorithms need information sources or features from which the models can be computed. For splice site prediction, the features we consider in this study are the presence or absence of certain nucleotides in close proximity to the splice site. Since it is not known how many and which nucleotides are relevant for splice site prediction, the set of features is chosen large enough such that the probability that all relevant information sources are in the set is very high. Using only those features that are relevant for constructing a splice site prediction system might improve the system and might also provide us with useful biological knowledge. Using fewer features will of course also improve the prediction speed of the system.Results: A wrapper-based feature subset selection algorithm using a support vector machine or a naive Bayes prediction method was evaluated against the traditional method for selecting features relevant for splice site prediction. Our results show that this wrapper approach selects features that improve the performance against the use of all features and against the use of the features selected by the traditional method.Availability: The data and additional interactive graphs on the selected feature subsets are available at http://www.psb.rug.ac.be/gpsContact: svgro@gengenp.rug.ac.beyvdp@gengenp.rug.ac.be

The large amount of available annotated Arabidopsis thaliana sequences allows the induction of splice site prediction models with supervised learning algorithms (see Haussler (1998) for a review and references). These algorithms need information sources or features from which the models can be computed. For splice site prediction, the features we consider in this study are the presence or absence of certain nucleotides in close proximity to the splice site. Since it is not known how many and which nucleotides are relevant for splice site prediction, the set of features is chosen large enough such that the probability that all relevant information sources are in the set is very high. Using only those features that are relevant for constructing a splice site prediction system might improve the system and might also provide us with useful biological knowledge. Using fewer features will of course also improve the prediction speed of the system.

There are approximately 25 000 species in the division Teleostei and most are believed to have arisen during a relatively short period of time ca.200 Myr ago. The discovery of ‘extra’ Hoxgene clusters in zebrafish (Danio rerio), medaka (Oryzias latipes), and pufferfish (Fugu rubripes), has led to the hypothesis that genome duplication provided the genetic raw material necessary for the teleost radiation. We identified 27 groups of orthologous genes which included one gene from man, mouse and chicken, one or two genes from tetraploid Xenopusand two genes from zebrafish. A genome duplication in the ancestor of teleost fishes is the most parsimonious explanation for the observations that for 15 of these genes, the two zebrafish orthologues are sister sequences in phylogenies that otherwise match the expected organismal tree, the zebrafish gene pairs appear to have been formed at approximately the same time, and are unlinked. Phylogenies of nine genes differ a little from the tree predicted by the fish–specific genome duplication hypothesis: one tree shows a sister sequence relationship for the zebrafish genes but differs slightly from the expected organismal tree and in eight trees, one zebrafish gene is the sister sequence to a clade which includes the second zebrafish gene and orthologues from Xenopus, chicken, mouse and man. For these nine gene trees, deviations from the predictions of the fish–specific genome duplication hypothesis are poorly supported. The two zebrafish orthologues for each of the three remaining genes are tightly linked and are, therefore, unlikely to have been formed during a genome duplication event. We estimated that the unlinked duplicated zebrafish genes are between 300 and 450 Myr. Thus, genome duplication could have provided the genetic raw material for teleost radiation. Alternatively, the loss of different duplicates in different populations (i.e. ‘divergent resolution’) may have promoted speciation in ancient teleost populations.

Chlamydomyxa labyrinthuloides is a heterokont alga known since the last century. It lives on Sphagnum and other water plants as aplanospores or plasmodia. We have investigated the taxonomic position of Chlamydomyxa labyrinthuloides by combining results from morphological studies, pigment analyses and a molecular phylogenetic analysis of the small subunit (SSU) rRNA gene. Chlamydomyxa labyrinthuloides shares morphological features with xanthophytes and chrysophytes, whereas pigment composition indicates a grouping with the phaeophytes, raphidophytes and chrysophytes. The sequence of the SSU rRNA gene and its phylogenetic reconstruction unambiguously demonstrate that Chlamydomyxa labyrinthuloides is related to the chrysophytes.

Evolutionary trees based on partial small ribosomal subunit RNA sequences of 22 metazoa species have been published [(1988) Science 239, 748-753]. In these trees, cnidarians (Radiata) seemed to have evolved independently from the Bilateria, which is in contradiction with the general evolutionary view. In order to further investigate this problem, the complete srRNA sequence of the sea anemone Ammonia sulcatawas determined and evolutionary trees were constructed using a matrix optimization method. In the tree thus obtained the sea anemone and Bilateria together form a monophyletic cluster, with the sea anemone forming the first line of descent of the metazoan group.

Evolutionary trees were constructed, by distance methods, from an alignment of 225 complete large subunit (LSU) rRNA sequences, representing Eucarya, Archaea, Bacteria, plastids, and mitochondria. A comparison was made with trees based on sets of small subunit (SSU) rRNA sequences. Trees constructed on the set of 172 species and organelles for which the sequences of both molecules are known had a very similar topology, at least with respect to the divergence order of large taxa such as the eukaryotic kingdoms and the bacterial divisions. However, since there are more than ten times as many SSU as LSU rRNA sequences, it is possible to select many SSU rRNA sequence sets of equivalent size but different species composition. The topologies of these trees showed considerable differences according to the particular species set selected.

Summary A tree was constructed from a structurally conserved area in an alignment of 83 small ribosomal subunit sequences of eukaryotic, archaebacterial, eubacterial, plastidial, and mitochondrial origin. The algorithm involved computation and optimization of a dissimilarity matrix. According to the tree, only plant mitochondria belong to the eubacterial primary kingdom, whereas animal, fungal, algal, and ciliate mitochondria branch off from an internal node situated between the three primary kingdoms. This result is at variance with a parsimony tree of similar size published by Cedergren et al. (J Mol Evol 28∶98–112, 1988), which postulates the mitochondria to be monophyletic and to belong to the eubacterial primary kingdom. The discrepancy does not follow from the use of conflicting sequence alignments, hence it must be due to the use of different treeing algorithms. We tested our algorithm on a set of sequences resulting from a simulated evolution and found it capable of faith-fully reconstructing a branching topology that involved very unequal evolutionary rates. The use of more limited or more extended areas of the complete sequence alignment, comprising only very conserved or also more variable portions of the small ribosomal subunit structure, does have some influence on the tree topology. In all cases, however, the nonplant mitochondria seem to branch off before the emergence of eubacteria, and the differences are limited to the branching pattern among different types of mitochondria.

We have begun to reconstruct the ancient history of the nematode phylum based on cytochrome c and globin amino acid sequences. The data suggest that the nematode ancestor diverged from a line leading to mammals about 1 billion years ago and that the most recent common ancestor of the extant species Caenorhabditis elegans, Trichostrongylus colubriformis, Nippostrongylus brasiliensis, Ascaris suum,and Pseudoterranova decipienslived about 550 MY ago. The rhabditids and strongylids emerged as one offshoot of this ancestor, the ascarids as another. Rhabditids and strongylids diverged some 400 MY ago, whereas the genera Trichostrongylusand Nippostrongylusdiverged slightly over 200 MY ago. A gene duplication event in the strongylid branch is predicted to have occurred around 250-335 MY ago. There are two globin genes in Nippostrongylus, expressed in anatomically distinct compartments (body and cuticle), and the single sequence from Trichostrongylusis most like the Nippostrongylusbody globin gene. A strikingly different duplication event occurred within the same period in the line leading to the extant ascarid genera, creating a single polypeptide containing two globin domains. The genera Ascarisand Pseudoterranovadiverged some 150-250 MY ago. Interestingly, the second globin repeat evolved at a faster rate in both species examined. This is possibly related to the acquisition of an unusual carboxyterminal extension, composed of alternating positively and negatively charged residues, that is necessary for the assembly of several monomers into the native polymeric molecules.

The complete nucleotide sequence of the small ribosomal subunit RNA of the gastropod, Limicolaria kambeul, was determined and used to infer a secondary structure model. In order to clarify the phylogenetic position of the Mollusca among the Metazoa, an evolutionary tree was constructed by neighbor-joining, starting from an alignment of small ribosomal subunit RNA sequences. The Mollusca appear to be a monophyletic group, related to Arthropoda and Chordata in an unresolved trichotomy.

Abstract.: In this study we constructed a bootstrapped distance tree of 500 small subunit ribosomal RNA sequences from organisms belonging to the so-called crown of eukaryote evolution. Taking into account the substitution rate of the individual nucleotides of the rRNA sequence alignment, our results suggest that (1) animals, true fungi, and choanoflagellates share a common origin: The branch joining these taxa is highly supported by bootstrap analysis (bootstrap support [BS] > 90%), (2) stramenopiles and alveolates are sister groups (BS = 75%), (3) within the alveolates, dinoflagellates and apicomplexans share a common ancestor BS > 95%), while in turn they both share a common origin with the ciliates (BS > 80%), and (4) within the stramenopiles, heterokont algae, hyphochytriomycetes, and oomycetes form a monophyletic grouping well supported by bootstrap analysis (BS > 85%), preceded by the well-supported successive divergence of labyrinthulomycetes and bicosoecids. On the other hand, many evolutionary relationships between crown taxa are still obscure on the basis of 18S rRNA. The branching order between the animal-fungal-choanoflagellates clade and the chlorobionts, the alveolates and stramenopiles, red algae, and several smaller groups of organisms remains largely unresolved. When among-site rate variation is not considered, the inferred tree topologies are inferior to those where the substitution rate spectrum for the 18S rRNA is taken into account. This is primarily indicated by the erroneous branching of fast-evolving sequences. Moreover, when different substitution rates among sites are not considered, the animals no longer appear as a monophyletic grouping in most distance trees.

Bacterial diversity in the water column of a freshwater lake in the Netherlands was investigated by analysis of 16S rRNA gene sequences recovered through PCR amplification from total community DNA. Among 23 unique cloned sequences, two appeared to belong to the recently described bacterial division Verrucomicrobiales. One of the two sequences was most similar to a group of environmental clones that form a distinct lineage within the division. The other sequence was divergent (less than 85% similarity) from all 16S rRNA gene sequences, both from cultivated species and from environmental clones, known in this division to date. Analysis by denaturing gradient gel electrophoresis (DGGE) and sequencing of DNA recovered through excision from the DGGE gel showed that the two sequence types were present in the lake throughout the year.