全基因组复制在动植物中普遍存在, 被认为是促进物种进化的重要动力之一。作为蕨类植物的单种科物种, 翼盖蕨(Didymochlaena trancatula)是真水龙骨类I的基部类群, 在蕨类中具有独特的演化地位。本研究基于高通量测序, 通过同义替换率(Ks)分析、相对定年分析揭示翼盖蕨的全基因组复制发生情况。Ks分析表明, 翼盖蕨至少经历了两次全基因组复制事件, 其中一次发生于59-62 million years ago (Mya), 另一次发生于90-94 Mya, 这两次全基因组复制事件分别和白垩纪第三纪的Cretaceous-Tertiary (C-T)大灭绝事件以及翼盖蕨的物种分化时间相吻合。进一步对两次全基因组复制保留的基因进行功能注释和富集分析, 结果显示与转录及代谢调控相关的基因优势被保留。翼盖蕨的全基因组复制事件可能促进了该物种的分化及其对极端环境的适应性。

The ongoing climate change is characterized by increased temperatures and altered precipitation patterns. In addition, there has been an increase in both the frequency and intensity of extreme climatic events such as drought. Episodes of drought induce a series of interconnected effects, all of which have the potential to alter the carbon balance of forest ecosystems profoundly at different scales of plant organization and ecosystem functioning. During recent years, considerable progress has been made in the understanding of how aboveground parts of trees respond to drought and how these responses affect carbon assimilation. In contrast, processes of belowground parts are relatively underrepresented in research on climate change. In this review, we describe current knowledge about responses of tree roots to drought. Tree roots are capable of responding to drought through a variety of strategies that enable them to avoid and tolerate stress. Responses include root biomass adjustments, anatomical alterations, and physiological acclimations. The molecular mechanisms underlying these responses are characterized to some extent, and involve stress signaling and the induction of numerous genes, leading to the activation of tolerance pathways. In addition, mycorrhizas seem to play important protective roles. The current knowledge compiled in this review supports the view that tree roots are well equipped to withstand drought situations and maintain morphological and physiological functions as long as possible. Further, the reviewed literature demonstrates the important role of tree roots in the functioning of forest ecosystems and highlights the need for more research in this emerging field.

The use of some multiple-sequence alignments in phylogenetic analysis, particularly those that are not very well conserved, requires the elimination of poorly aligned positions and divergent regions, since they may not be homologous or may have been saturated by multiple substitutions. A computerized method that eliminates such positions and at the same time tries to minimize the loss of informative sites is presented here. The method is based on the selection of blocks of positions that fulfill a simple set of requirements with respect to the number of contiguous conserved positions, lack of gaps, and high conservation of flanking positions, making the final alignment more suitable for phylogenetic analysis. To illustrate the efficiency of this method, alignments of 10 mitochondrial proteins from several completely sequenced mitochondrial genomes belonging to diverse eukaryotes were used as examples. The percentages of removed positions were higher in the most divergent alignments. After removing divergent segments, the amino acid composition of the different sequences was more uniform, and pairwise distances became much smaller. Phylogenetic trees show that topologies can be different after removing conserved blocks, particularly when there are several poorly resolved nodes. Strong support was found for the grouping of animals and fungi but not for the position of more basal eukaryotes. The use of a computerized method such as the one presented here reduces to a certain extent the necessity of manually editing multiple alignments, makes the automation of phylogenetic analysis of large data sets feasible, and facilitates the reproduction of the final alignment by other researchers.

Polyploids - organisms that have multiple sets of chromosomes - are common in certain plant and animal taxa, and can be surprisingly stable. The evidence that has emerged from genome analyses also indicates that many other eukaryotic genomes have a polyploid ancestry, suggesting that both humans and most other eukaryotes have either benefited from or endured polyploidy. Studies of polyploids soon after their formation have revealed genetic and epigenetic interactions between redundant genes. These interactions can be related to the phenotypes and evolutionary fates of polyploids. Here, I consider the advantages and challenges of polyploidy, and its evolutionary potential.

Nuclear DNA contents (4C) were estimated by Feulgen microdensitometry in 27 species of slipper orchids. These data and recent information concerning the molecular systematics of Cypripedioideae allow an interesting re-evaluation of karyotype and genome size variation among slipper orchids in a phylogenetic context. DNA amounts differed 5.7-fold, from 24.4 pg in Phragmipedium longifolium to 138.1 pg in Paphiopedilum wardii. The most derived clades of the conduplicate-leaved slipper orchids have undergone a radical process of genome fragmentation that is most parsimoniously explained by Robertsonian changes involving centric fission. This process seems to have occurred independently of genome size variation. However, it may reflect environmental or selective pressures favoring higher numbers of linkage groups in the karyotype.

We have implemented a high-performance computing (HPC) version of ProtTest that can be executed in parallel in multicore desktops and clusters. This version, called ProtTest 3, includes new features and extended capabilities.ProtTest 3 source code and binaries are freely available under GNU license for download from http://darwin.uvigo.es/software/prottest3, linked to a Mercurial repository at Bitbucket (https://bitbucket.org/).dposada@uvigo.esSupplementary data are available at Bioinformatics online.

Detection of protein families in large databases is one of the principal research objectives in structural and functional genomics. Protein family classification can significantly contribute to the delineation of functional diversity of homologous proteins, the prediction of function based on domain architecture or the presence of sequence motifs as well as comparative genomics, providing valuable evolutionary insights. We present a novel approach called TRIBE-MCL for rapid and accurate clustering of protein sequences into families. The method relies on the Markov cluster (MCL) algorithm for the assignment of proteins into families based on precomputed sequence similarity information. This novel approach does not suffer from the problems that normally hinder other protein sequence clustering algorithms, such as the presence of multi-domain proteins, promiscuous domains and fragmented proteins. The method has been rigorously tested and validated on a number of very large databases, including SwissProt, InterPro, SCOP and the draft human genome. Our results indicate that the method is ideally suited to the rapid and accurate detection of protein families on a large scale. The method has been used to detect and categorise protein families within the draft human genome and the resulting families have been used to annotate a large proportion of human proteins.

Paphiopedilum is an important genus of the orchid family Orchidaceae and has high horticultural value. The wild populations are under threat of extinction because of overcollection and habitat destruction. Mature seeds of most Paphiopedilum species are difficult to germinate, which severely restricts their germplasm conservation and commercial production. The factors inhibiting germination are largely unknown.In this study, large amounts of non-methylated lignin accumulated during seed maturation of Paphiopedilum armeniacum (P. armeniacum), which negatively correlates with the germination rate. The transcriptome profiles of P. armeniacum seed at different development stages were compared to explore the molecular clues for non-methylated lignin synthesis. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that a large number of genes associated with phenylpropanoid biosynthesis and phenylalanine metabolism during seed maturation were differentially expressed. Several key genes in the lignin biosynthetic pathway displayed different expression patterns during the lignification process. PAL, 4CL, HCT, and CSE upregulation was associated with C and H lignin accumulation. The expression of CCoAOMT, F5H, and COMT were maintained at a low level or down-regulated to inhibit the conversion to the typical G and S lignin. Quantitative real-time RT-PCR analysis confirmed the altered expression levels of these genes in seeds and vegetative tissues.This work demonstrated the plasticity of natural lignin polymer assembly in seed and provided a better understanding of the molecular mechanism of seed-specific lignification process.

A codon-based model for the evolution of protein-coding DNA sequences is presented for use in phylogenetic estimation. A Markov process is used to describe substitutions between codons. Transition/transversion rate bias and codon usage bias are allowed in the model, and selective restraints at the protein level are accommodated using physicochemical distances between the amino acids coded for by the codons. Analyses of two data sets suggest that the new codon-based model can provide a better fit to data than can nucleotide-based models and can produce more reliable estimates of certain biologically important measures such as the transition/transversion rate ratio and the synonymous/nonsynonymous substitution rate ratio.

Background: The temporal origin and diversification of orchids (family Orchidaceae) has been subject to intense debate in the last decade. The description of the first reliable fossil in 2007 enabled a direct calibration of the orchid phylogeny, but little attention has been paid to the potential influence of dating methodology in obtaining reliable age estimates. Moreover, two new orchid fossils described in 2009 have not yet been incorporated in a molecular dating analysis. Here we compare the ages of major orchid clades estimated under two widely used methods, a Bayesian relaxed clock implemented in BEAST and Penalized Likelihood implemented in r8s. We then perform a new family-level analysis by integrating all 3 available fossils and using BEAST. To evaluate how the newly estimated ages may influence the evolutionary interpretation of a species-level phylogeny, we assess divergence times for the South American genus Hoffmannseggella (subfam. Epidendroideae), for which we present an almost complete phylogeny (40 out of 41 species sampled). Results: Our results provide additional support that all extant orchids shared a most recent common ancestor in the Late Cretaceous (similar to 77 million years ago, Ma). However, we estimate the crown age of the five orchid subfamilies to be generally (similar to 1-8 Ma) younger than previously calculated under the Penalized Likelihood algorithm and using a single internal fossil calibration. The crown age of Hoffmannseggella is estimated here at similar to 11 Ma, some 3 Ma more recently than estimated under Penalized Likelihood. Conclusions: Contrary to recent suggestions that orchid diversification began in a period of global warming, our results place the onset of diversification of the largest orchid subfamilies (Orchidoideae and Epidendroideae) in a period of global cooling subsequent to the Early Eocene Climatic Optimum. The diversification of Hoffmannseggella appears even more correlated to late Tertiary climatic fluctuations than previously suggested. With the incorporation of new fossils in the orchid phylogeny and the use of a method that is arguably more adequate given the present data, our results represent the most up-to-date estimate of divergence times in orchids.

Orthology assignment is ideally suited for functional inference. However, because predicting orthology is computationally intensive at large scale, and most pipelines are relatively inaccessible (e.g., new assignments only available through database updates), less precise homology-based functional transfer is still the default for (meta-)genome annotation. We, therefore, developed eggNOG-mapper, a tool for functional annotation of large sets of sequences based on fast orthology assignments using precomputed clusters and phylogenies from the eggNOG database. To validate our method, we benchmarked Gene Ontology (GO) predictions against two widely used homology-based approaches: BLAST and InterProScan. Orthology filters applied to BLAST results reduced the rate of false positive assignments by 11%, and increased the ratio of experimentally validated terms recovered over all terms assigned per protein by 15%. Compared with InterProScan, eggNOG-mapper achieved similar proteome coverage and precision while predicting, on average, 41 more terms per protein and increasing the rate of experimentally validated terms recovered over total term assignments per protein by 35%. EggNOG-mapper predictions scored within the top-5 methods in the three GO categories using the CAFA2 NK-partial benchmark. Finally, we evaluated eggNOG-mapper for functional annotation of metagenomics data, yielding better performance than interProScan. eggNOG-mapper runs ∼15× faster than BLAST and at least 2.5× faster than InterProScan. The tool is available standalone and as an online service at http://eggnog-mapper.embl.de.© The Author 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.

Whole-genome duplications are widespread across land plant phylogenies and particularly frequent in ferns and angiosperms. Genome duplications spurred the evolution of key innovations associated with diversification in many angiosperm clades and lineages. Such diversifications are not initiated by genome doubling per se. Rather, differentiation of the primary polyploid populations through a range of processes results in post-polyploid genome diploidization. Structural diploidization gradually reverts the polyploid genome to one functionally diploid-like through chromosomal rearrangements which frequently result in dysploid changes. Dysploidies may lead to reproductive isolation among post-polyploid offspring and significantly contribute to speciation and cladogenetic events.Copyright © 2018 Elsevier Ltd. All rights reserved.

Pineapple (Ananas comosus (L.) Merr.) is the most economically valuable crop possessing crassulacean acid metabolism (CAM), a photosynthetic carbon assimilation pathway with high water-use efficiency, and the second most important tropical fruit. We sequenced the genomes of pineapple varieties F153 and MD2 and a wild pineapple relative, Ananas bracteatus accession CB5. The pineapple genome has one fewer ancient whole-genome duplication event than sequenced grass genomes and a conserved karyotype with seven chromosomes from before the ρ duplication event. The pineapple lineage has transitioned from C3 photosynthesis to CAM, with CAM-related genes exhibiting a diel expression pattern in photosynthetic tissues. CAM pathway genes were enriched with cis-regulatory elements associated with the regulation of circadian clock genes, providing the first cis-regulatory link between CAM and circadian clock regulation. Pineapple CAM photosynthesis evolved by the reconfiguration of pathways in C3 plants, through the regulatory neofunctionalization of preexisting genes and not through the acquisition of neofunctionalized genes via whole-genome or tandem gene duplication.

We describe here the reconstruction of the genome of the most recent common ancestor (MRCA) of modern monocots and eudicots, accounting for 95% of extant angiosperms, with its potential repertoire of 22,899 ancestral genes conserved in present-day crops. The MRCA provides a starting point for deciphering the reticulated evolutionary plasticity between species (rapidly versus slowly evolving lineages), subgenomes (pre-versus post-duplication blocks), genomic compartments (stable versus labile loci), genes (ancestral versus species-specific genes) and functions (gained versus lost ontologies), the key mutational forces driving the success of polyploidy in crops. The estimation of the timing of angiosperm evolution, based on MRCA genes, suggested that this group emerged 214 million years ago during the late Triassic era, before the oldest recorded fossil. Finally, the MRCA constitutes a unique resource for scientists to dissect major agronomic traits in translational genomics studies extending from model species to crops.

Whole-genome duplication (WGD) in angiosperms has been hypothesized to be advantageous in unstable environments and/or to increase diversification rates, leading to radiations. Under the first hypothesis, floras in stable environments are predicted to have lower proportions of polyploids than highly, recently disturbed floras, whereas species-rich floras would be expected to have higher than expected proportions of polyploids under the second. The South African Cape flora is used to discriminate between these two hypotheses because it features a hyperdiverse flora predominantly generated by a limited number of radiations (Cape clades), against a backdrop of climatic and geological stability.We compiled all known chromosome counts for species in 21 clades present in the Cape (1653 species, including 24 Cape clades), inferred ploidy levels for these species by inspection or derived from the primary literature, and compared Cape to non-Cape ploidy levels in these clades (17,520 species) using G tests.The Cape flora has anomalously low proportions of polyploids compared with global levels. This pattern is consistently observed across nearly half the clades and across global latitudinal gradients, although individual lineages seem to be following different paths to low levels of WGD and to differing degrees.This pattern shows that the diversity of the Cape flora is the outcome of primarily diploid radiations and supports the hypothesis that WGD may be rare in stable environments.© 2016 Botanical Society of America.

Gene duplications provide evolutionary potentials for generating novel functions, while polyploidization or whole genome duplication (WGD) doubles the chromosomes initially and results in hundreds to thousands of retained duplicates. WGDs are strongly supported by evidence commonly found in many species-rich lineages of eukaryotes, and thus are considered as a major driving force in species diversification. We performed comparative genomic and phylogenomic analyses of 59 public genomes/transcriptomes and 46 newly sequenced transcriptomes covering major lineages of angiosperms to detect large-scale gene duplication events by surveying tens of thousands of gene family trees. These analyses confirmed most of the previously reported WGDs and provided strong evidence for novel ones in many lineages. The detected WGDs supported a model of exponential gene loss during evolution with an estimated half-life of approximately 21.6 million years, and were correlated with both the emergence of lineages with high degrees of diversification and periods of global climate changes. The new datasets and analyses detected many novel WGDs widely spread during angiosperm evolution, uncovered preferential retention of gene functions in essential cellular metabolisms, and provided clues for the roles of WGD in promoting angiosperm radiation and enhancing their adaptation to environmental changes.Copyright © 2018 The Author. Published by Elsevier Inc. All rights reserved.

Using plant EST collections, we obtained 1392 potential gene duplicates across 8 plant species: Zea mays, Oryza sativa, Sorghum bicolor, Hordeum vulgare, Solanum tuberosum, Lycopersicon esculentum, Medicago truncatula, and Glycine max. We estimated the synonymous and nonsynonymous distances between each gene pair and identified two to three mixtures of normal distributions corresponding to one to three rounds of genome duplication in each species. Within the Poaceae, we found a conserved duplication event among all four species that occurred approximately 50-60 million years ago (Mya); an event that probably occurred before the major radiation of the grasses. In the Solanaceae, we found evidence for a conserved duplication event approximately 50-52 Mya. A duplication in soybean occurred approximately 44 Mya and a duplication in Medicago about 58 Mya. Comparing synonymous and nonsynonymous distances allowed us to determine that most duplicate gene pairs are under purifying, negative selection. We calculated Pearson's correlation coefficients to provide us with a measure of how gene expression patterns have changed between duplicate pairs, and compared this across evolutionary distances. This analysis showed that some duplicates seemed to retain expression patterns between pairs, whereas others showed uncorrelated expression.

Finite mixture models are being used increasingly to model a wide variety of random phenomena for clustering, classification and density estimation. is a powerful and popular package which allows modelling of data as a Gaussian finite mixture with different covariance structures and different numbers of mixture components, for a variety of purposes of analysis. Recently, version 5 of the package has been made available on CRAN. This updated version adds new covariance structures, dimension reduction capabilities for visualisation, model selection criteria, initialisation strategies for the EM algorithm, and bootstrap-based inference, making it a full-featured R package for data analysis via finite mixture modelling.

Alignment quality may have as much impact on phylogenetic reconstruction as the phylogenetic methods used. Not only the alignment algorithm, but also the method used to deal with the most problematic alignment regions, may have a critical effect on the final tree. Although some authors remove such problematic regions, either manually or using automatic methods, in order to improve phylogenetic performance, others prefer to keep such regions to avoid losing any information. Our aim in the present work was to examine whether phylogenetic reconstruction improves after alignment cleaning or not. Using simulated protein alignments with gaps, we tested the relative performance in diverse phylogenetic analyses of the whole alignments versus the alignments with problematic regions removed with our previously developed Gblocks program. We also tested the performance of more or less stringent conditions in the selection of blocks. Alignments constructed with different alignment methods (ClustalW, Mafft, and Probcons) were used to estimate phylogenetic trees by maximum likelihood, neighbor joining, and parsimony. We show that, in most alignment conditions, and for alignments that are not too short, removal of blocks leads to better trees. That is, despite losing some information, there is an increase in the actual phylogenetic signal. Overall, the best trees are obtained by maximum-likelihood reconstruction of alignments cleaned by Gblocks. In general, a relaxed selection of blocks is better for short alignment, whereas a stringent selection is more adequate for longer ones. Finally, we show that cleaned alignments produce better topologies although, paradoxically, with lower bootstrap. This indicates that divergent and problematic alignment regions may lead, when present, to apparently better supported although, in fact, more biased topologies.

The phylogeny and biogeography of the genus were evaluated by using phylogenetic trees derived from analysis of nuclear ribosomal internal transcribed spacer (ITS) sequences, the plastid L intron, the L-F spacer, and the B-L spacer. This genus was divided into three subgenera:,, and. Each of them is monophyletic with high bootstrap supports according to the highly resolved phylogenetic tree reconstructed by combined sequences. There are five sections within the subgenus, including,,,, and. The subgenus is phylogenetic basal, which suggesting that is comprising more ancestral characters than other subgenera. The evolutionary trend of genus was deduced based on the maximum likelihood (ML) tree and Bayesian Evolutionary Analysis Sampling Trees (BEAST). Reconstruct Ancestral State in Phylogenies (RASP) analyses based on the combined sequence data. The biogeographic analysis indicates that species were firstly derived in Southern China and Southeast Asia, subsequently dispersed into the Southeast Asian archipelagoes. The subgenera was likely derived after these historical dispersals and vicariance events. Our research reveals the relevance of the differentiation of in Southeast Asia and geological history. Moreover, the biogeographic analysis explains that the significant evolutionary hotspots of these orchids in the Sundaland and Wallacea might be attributed to repeated migration and isolation events between the south-eastern Asia mainland and the Sunda Super Islands.Copyright © 2020 Tsai, Liao, Ko, Chen and Chiang.

Stress acclimating plants respond to abiotic and biotic stress by remodeling membrane fluidity and by releasing alpha-linolenic (18:3) from membrane lipids. The modification of membrane fluidity is mediated by changes in unsaturated fatty acid levels, a function provided in part by the regulated activity of fatty acid desaturases. Adjustment of membrane fluidity maintains an environment suitable for the function of critical integral proteins during stress. alpha-Linolenic acid, released from membrane lipid by regulated lipase activity, is the precursor molecule for phyto-oxylipin biosynthesis. The modulation of chloroplast oleic acid (18:1) levels is central to the normal expression of defense responses to pathogens in Arabidopsis. Oleic (18:1) and linolenic (18:2) acid levels, in part, regulate development, seed colonization, and mycotoxin production by Aspergillus spp.

Abiotic stress is one of the severe stresses of environment that lowers the growth and yield of any crop even on irrigated land throughout the world. A major phytohormone abscisic acid (ABA) plays an essential part in acting toward varied range of stresses like heavy metal stress, drought, thermal or heat stress, high level of salinity, low temperature, and radiation stress. Its role is also elaborated in various developmental processes including seed germination, seed dormancy, and closure of stomata. ABA acts by modifying the expression level of gene and subsequent analysis of cis- and trans-acting regulatory elements of responsive promoters. It also interacts with the signaling molecules of processes involved in stress response and development of seeds. On the whole, the stress to a plant can be susceptible or tolerant by taking into account the coordinated activities of various stress-responsive genes. Numbers of transcription factor are involved in regulating the expression of ABA responsive genes by acting together with their respective cis- acting elements. Hence, for improvement in stress-tolerance capacity of plants, it is necessary to understand the mechanism behind it. On this ground, this article enlightens the importance and role of ABA signaling with regard to various stresses as well as regulation of ABA biosynthetic pathway along with the transcription factors for stress tolerance.

Polyploidy is a prominent process in plants and has been significant in the evolutionary history of vertebrates and other eukaryotes. In plants, interdisciplinary approaches combining phylogenetic and molecular genetic perspectives have enhanced our awareness of the myriad genetic interactions made possible by polyploidy. Here, processes and mechanisms of gene and genome evolution in polyploids are reviewed. Genes duplicated by polyploidy may retain their original or similar function, undergo diversification in protein function or regulation, or one copy may become silenced through mutational or epigenetic means. Duplicated genes also may interact through inter-locus recombination, gene conversion, or concerted evolution. Recent experiments have illuminated important processes in polyploids that operate above the organizational level of duplicated genes. These include inter-genomic chromosomal exchanges, saltational, non-Mendelian genomic evolution in nascent polyploids, inter-genomic invasion, and cytonuclear stabilization. Notwithstanding many recent insights, much remains to be learned about many aspects of polyploid evolution, including: the role of transposable elements in structural and regulatory gene evolution; processes and significance of epigenetic silencing; underlying controls of chromosome pairing; mechanisms and functional significance of rapid genome changes; cytonuclear accommodation; and coordination of regulatory factors contributed by two, sometimes divergent progenitor genomes. Continued application of molecular genetic approaches to questions of polyploid genome evolution holds promise for producing lasting insight into processes by which novel genotypes are generated and ultimately into how polyploidy facilitates evolution and adaptation.

We present the 1.06 Gb sequenced genome of Gastrodia elata, an obligate mycoheterotrophic plant, which contains 18,969 protein-coding genes. Many genes conserved in other plant species have been deleted from the G. elata genome, including most of those for photosynthesis. Additional evidence of the influence of genome plasticity in the adaptation of this mycoheterotrophic lifestyle is evident in the large number of gene families that are expanded in G. elata, including glycoside hydrolases and urease that likely facilitate the digestion of hyphae are expanded, as are genes associated with strigolactone signaling, and ATPases that may contribute to the atypical energy metabolism. We also find that the plastid genome of G. elata is markedly smaller than that of green plant species while its mitochondrial genome is one of the largest observed to date. Our report establishes a foundation for studying adaptation to a mycoheterotrophic lifestyle.

Since 65 million years ago (Ma), Earth's climate has undergone a significant and complex evolution, the finer details of which are now coming to light through investigations of deep-sea sediment cores. This evolution includes gradual trends of warming and cooling driven by tectonic processes on time scales of 10(5) to 10(7) years, rhythmic or periodic cycles driven by orbital processes with 10(4)- to 10(6)-year cyclicity, and rare rapid aberrant shifts and extreme climate transients with durations of 10(3) to 10(5) years. Here, recent progress in defining the evolution of global climate over the Cenozoic Era is reviewed. We focus primarily on the periodic and anomalous components of variability over the early portion of this era, as constrained by the latest generation of deep-sea isotope records. We also consider how this improved perspective has led to the recognition of previously unforeseen mechanisms for altering climate.

Ancient whole-genome duplications (WGDs) have been uncovered in almost all major lineages of life on Earth and the search for traces or remnants of such events has become standard practice in most genome analyses. This is especially true for plants, where ancient WGDs are abundant. Common approaches to find evidence for ancient WGDs include the construction of KS distributions and the analysis of intragenomic colinearity. Despite the increased interest in WGDs and the acknowledgment of their evolutionary importance, user-friendly and comprehensive tools for their analysis are lacking. Here, we present an easy to use command-line tool for KS distribution construction named wgd. The wgd suite provides commonly used KS and colinearity analysis workflows together with tools for modeling and visualization, rendering these analyses accessible to genomics researchers in a convenient manner.wgd is free and open source software implemented in Python and is available at https://github.com/arzwa/wgd.Supplementary data are available at Bioinformatics online.© The Author(s) 2018. Published by Oxford University Press.