Evolutionary conservation,diversity and specificity of LTR‐retrotransposons in flowering plants: insights from genome‐wide analysis and multi‐specific comparison

The availability of complete or nearly complete genome sequences from several plant species permits detailed discovery and cross‐species comparison of transposable elements (TEs) at the whole genome level. We initially investigated 510 long terminal repeat‐retrotransposon (LTR‐RT) families comprising 32 370 elements in soybean (Glycine max (L.) Merr.). Approximately 87% of these elements were located in recombination‐suppressed pericentromeric regions, where the ratio (1.26) of solo LTRs to intact elements (S/I) is significantly lower than that of chromosome arms (1.62). Further analysis revealed a significant positive correlation between S/I and LTR sizes, indicating that larger LTRs facilitate solo LTR formation. Phylogenetic analysis revealed seven Copia and five Gypsy evolutionary lineages that were present before the divergence of eudicot and monocot species, but the scales and timeframes within which they proliferated vary dramatically across families, lineages and species, and notably, a Copia lineage has been lost in soybean. Analysis of the physical association of LTR‐RTs with centromere satellite repeats identified two putative centromere retrotransposon (CR) families of soybean, which were grouped into the CR (e.g. CRR and CRM) lineage found in grasses, indicating that the ‘functional specification’ of CR pre‐dates the bifurcation of eudicots and monocots. However, a number of families of the CR lineage are not concentrated in centromeres, suggesting that their CR roles may now be defunct. Our data also suggest that the envelope‐like genes in the putative Copia retrovirus‐like family are probably derived from the Gypsy retrovirus‐like lineage, and thus we propose the hypothesis of a single ancient origin of envelope‐like genes in flowering plants.     

Evolutionary and ecological genomics of non‐model plants

Abstract Dissecting evolutionary dynamics of ecologically important traits is a long‐term challenge for biologists. Attempts to understand natural variation and molecular mechanisms have motivated a move from laboratory model systems to non‐model systems in diverse natural environments. Next generation sequencing methods, along with an expansion of genomic resources and tools, have fostered new links between diverse disciplines, including molecular biology, evolution, ecology, and genomics. Great progress has been made in a few non‐model wild plants, such as Arabidopsis relatives, monkey flowers, and wild sunflowers. Until recently, the lack of comprehensive genomic information has limited evolutionary and ecological studies to larger QTL (quantitative trait locus) regions rather than single gene resolution, and has hindered recognition of general patterns of natural variation and local adaptation. Further efforts in accumulating genomic data and developing bioinformatic and biostatistical tools are now poised to move this field forward. Integrative national and international collaborations and research communities are needed to facilitate development in the field of evolutionary and ecological genomics.     

稻属植物的基因组进化

水稻所在的稻属(Oryza)共有24个左右的物种。由于野生稻含有大量的优良农艺性状基因, 在水稻遗传学研究中日益受到重视。随着国际稻属基因组计划的开展, 越来越多的稻属基因组序列被测定, 稻属成为进行比较、功能和进化基因组学研究的模式系统。近期开展的一系列研究对稻属不同基因组区段以及全基因组序列的比较分析, 揭示了稻属在基因组大小、基因移动、多倍体进化、常染色质到异染色质的转化以及着丝粒区域的进化等方面的分子机制。转座子的活性以及转座子因非均等重组或非法重组而造成的删除, 对稻属基因组的扩增和收缩具有重要作用。DNA双链断裂修复介导的基因移动, 特别是非同源末端连接, 是稻属基因组非共线性基因形成的主要来源。稻属基因组从常染色质到异染色质的转换过程, 伴随着转座子的大量扩增、基因片段的区段性和串联重复以及从基因组其他位置不断捕获异染色质基因。对稻属不同物种间基因拷贝数、特异基因和重要农艺性状基因的进化等研究, 可揭示稻属不同物种间表型和适应性差异的分子基础, 将加速水稻的育种和改良。     

Spatio‐temporal patterns of genome evolution in allotetraploid species of the genus Oryza

Despite knowledge that polyploidy is widespread and a major evolutionary force in flowering plant diversification, detailed comparative molecular studies on polyploidy have been confined to only a few species and families. The genus Oryza is composed of 23 species that are classified into ten distinct ‘genome types’ (six diploid and four polyploid), and is emerging as a powerful new model system to study polyploidy. Here we report the identification, sequence and comprehensive comparative annotation of eight homoeologous genomes from a single orthologous region (Adh1–Adh2) from four allopolyploid species representing each of the known Oryza genome types (BC, CD, HJ and KL). Detailed comparative phylogenomic analyses of these regions within and across species and ploidy levels provided several insights into the spatio‐temporal dynamics of genome organization and evolution of this region in ‘natural’ polyploids of Oryza. The major findings of this study are that: (i) homoeologous genomic regions within the same nucleus experience both independent and parallel evolution, (ii) differential lineage‐specific selection pressures do not occur between polyploids and their diploid progenitors, (iii) there have been no dramatic structural changes relative to the diploid ancestors, (iv) a variation in the molecular evolutionary rate exists between the two genomes in the BC complex species even though the BC and CD polyploid species appear to have arisen <2 million years ago, and (v) there are no clear distinctions in the patterns of genome evolution in the diploid versus polyploid species.     

Genomic analyses provide insights into the polyploidization-driven herbicide adaptation in Leptochloa weeds

Polyploidy confers a selective advantage under stress conditions; however, whether polyploidization mediates enhanced herbicide adaptation remains largely unknown. Tetraploid Leptochloa chinensis is a notorious weed in the rice ecosystem, causing severe yield loss in rice. In China, L. chinensis has only one sister species, the diploid L. panicea, whose damage is rarely reported. To gain insights into the effects of polyploidization on herbicide adaptation, we first assembled a high-quality genome of L. panicea and identified genome structure variations with L. chinensis. Moreover, we identified herbicide-resistance genes specifically expanded in L. chinensis, which may confer a greater herbicide adaptability in L. chinensis. Analysis of gene retention and loss showed that five herbicide target-site genes and several herbicide nontarget-site resistance gene families were retained during polyploidization. Notably, we identified three pairs of polyploidization-retained genes including LcABCC8, LcCYP76C1 and LcCYP76C4 that may enhance herbicide resistance. More importantly, we found that both copies of LcCYP76C4 were under herbicide selection during the spread of L. chinensis in China. Furthermore, we identified another gene potentially involved in herbicide resistance, LcCYP709B2, which is also retained during polyploidization and under selection. This study provides insights into the genomic basis of the enhanced herbicide adaptability of Leptochloa weeds during polyploidization and provides guidance for the precise and efficient control of polyploidy weeds.     

Dynamic virulence‐related regions of the plant pathogenic fungus Verticillium dahliae display enhanced sequence conservation

Plant pathogens continuously evolve to evade host immune responses. During host colonization, many fungal pathogens secrete effectors to perturb such responses, but these in turn may become recognized by host immune receptors. To facilitate the evolution of effector repertoires, such as the elimination of recognized effectors, effector genes often reside in genomic regions that display increased plasticity, a phenomenon that is captured in the two‐speed genome hypothesis. The genome of the vascular wilt fungus Verticillium dahliae displays regions with extensive presence/absence polymorphisms, so‐called lineage‐specific regions, that are enriched in in planta‐induced putative effector genes. As expected, comparative genomics reveals differential degrees of sequence divergence between lineage‐specific regions and the core genome. Unanticipated, lineage‐specific regions display markedly higher sequence conservation in coding as well as noncoding regions than the core genome. We provide evidence that disqualifies horizontal transfer to explain the observed sequence conservation and conclude that sequence divergence occurs at a slower pace in lineage‐specific regions of the V. dahliae genome. We hypothesize that differences in chromatin organisation may explain lower nucleotide substitution rates in the plastic, lineage‐specific regions of V. dahliae.     

Identifying co‐opted transposable elements using comparative epigenomics

The human genome gives rise to different epigenomic landscapes that define each cell type and can be deregulated in disease. Recent efforts by ENCODE, the NIH Roadmap and the International Human Epigenome Consortium (IHEC) have made significant advances towards assembling reference epigenomic maps of various tissues. Notably, these projects have found that approximately 80% of human DNA was biochemically active in at least one epigenomic assay while only approximately 10% of the sequence displayed signs of purifying selection. Given that transposable elements (TEs) make up at least 50% of the human genome and can be actively transcribed or act as regulatory elements either for their own purposes or be co‐opted for the benefit of their host; we are interested in exploring their overall contribution to the “functional” genome. Traditional methods used to identify functional DNA have relied on comparative genomics, conservation analysis and low throughput validation assays. To discover co‐opted TEs, and distinguish them from noisy genomic elements, we argue that comparative epigenomic methods will also be important.     

Population genomic analyses from low‐coverage RAD‐Seq data: a case study on the non‐model cucurbit bottle gourd

Restriction site‐associated DNA sequencing (RAD‐Seq), a next‐generation sequencing‐based genome ‘complexity reduction’ protocol, has been useful in population genomics in species with a reference genome. However, the application of this protocol to natural populations of genomically underinvestigated species, particularly under low‐to‐medium sequencing depth, has not been well justified. In this study, a Bayesian method was developed for calling genotypes from an F₂ population of bottle gourd [Lagenaria siceraria (Mol.) Standl.] to construct a high‐density genetic map. Low‐depth genome shotgun sequencing allowed the assembly of scaffolds/contigs comprising approximately 50% of the estimated genome, of which 922 were anchored for identifying syntenic regions between species. RAD‐Seq genotyping of a natural population comprising 80 accessions identified 3226 single nuclear polymorphisms (SNPs), based on which two sub‐gene pools were suggested for association with fruit shape. The two sub‐gene pools were moderately differentiated, as reflected by the Hudson's F_ST value of 0.14, and they represent regions on LG7 with strikingly elevated F_ST values. Seven‐fold reduction in heterozygosity and two times increase in LD (r²) were observed in the same region for the round‐fruited sub‐gene pool. Outlier test suggested the locus LX3405 on LG7 to be a candidate site under selection. Comparative genomic analysis revealed that the cucumber genome region syntenic to the high F_ST island on LG7 harbors an ortholog of the tomato fruit shape gene OVATE. Our results point to a bright future of applying RAD‐Seq to population genomic studies for non‐model species even under low‐to‐medium sequencing efforts. The genomic resources provide valuable information for cucurbit genome research.     

Kmasker plants – a tool for assessing complex sequence space in plant species

Many plant genomes display high levels of repetitive sequences. The assembly of these complex genomes using short high‐throughput sequence reads is still a challenging task. Underestimation or disregard of repeat complexity in these datasets can easily misguide downstream analysis. Detection of repetitive regions by k‐mer counting methods has proved to be reliable. Easy‐to‐use applications utilizing k‐mer counting are in high demand, especially in the domain of plants. We present Kmasker plants, a tool that uses k‐mer count information as an assistant throughout the analytical workflow of genome data that is provided as a command‐line and web‐based solution. Beside its core competence to screen and mask repetitive sequences, we have integrated features that enable comparative studies between different cultivars or closely related species and methods that estimate target specificity of guide RNAs for application of site‐directed mutagenesis using Cas9 endonuclease. In addition, we have set up a web service for Kmasker plants that maintains pre‐computed indices for 10 of the economically most important cultivated plants. Source code for Kmasker plants has been made publically available at https://github.com/tschmutzer/kmasker . The web service is accessible at https://kmasker.ipk-gatersleben.de .     

Evolutionary dynamics of duplicated genes in plants

Gene duplication, arising from region-specific duplication or genome-wide polyploidization, is a prominent feature in plant genome evolution. Understanding the mechanisms generating duplicate gene copies and the subsequent dynamics among gene duplicates is vital because these investigations shed light on regional and genome-wide aspects of evolutionary forces shaping intra- and interspecific genome contents, evolutionary relationships, and interactions. This review discusses recent gene duplication analyses in plants, focusing on the molecular and evolutionary dynamics occurring at three different timescales following duplication: (1). initial establishment and persistence of cytotypes, (2). interactions among duplicate gene copies, and (3). longer term differentiation between duplicated genes. These relative time points are presented in terms of their potential adaptive significance and impact on plant evolutionary genomics research.     

A cycad's non‐saturating response to carbon dioxide enrichment indicates Cenozoic carbon limitation in pre‐historic plants

Cycads were a dominant plant functional type during the Mesozoic Era when atmospheric carbon dioxide [CO₂] greatly exceeded current conditions. Cycads, now rare and endangered, are slow‐growing perennial gymnosperms that develop carbon‐rich structural biomass, such as sclerophyllous leaves, dense stems and massive reproductive cones. Is cycad carbon partitioning to specific organs a constraint of their high [CO₂] evolutionary history (CO₂ legacy hypothesis, CLH)? To explore changes in cycad growth, carbon partitioning and assimilation responses that could be expected during the CO₂ depletion of the Cenozoic Era, individuals of the cycad species Encephalartos villosus plants were grown at four CO₂ levels: 400, 550, 750 and 1000 μmol mol^?1. The CLH predicts that cycad biomass and growth rates would increase in elevated [CO₂] due to increased net assimilation rates, and that carbon‐dense structures would provide sufficient carbohydrate sinks to prevent photosynthetic down‐regulation even under super‐ambient [CO₂] of 1000 μmol mol^?1. Both hypotheses were confirmed, though the latter less strongly. Plant relative growth rates increased 23% and biomass accumulation increased 65% in 1000 μmol mol^?1relative to 400 μmol mol^?1 treatment groups. Mean net assimilation rates increased 130% at 1000 μmol mol^?1 relative to 400 μmol mol^?1 CO₂, though there was some down‐regulation of maximum rate of carboxylation (Vc_max). Assimilation rates, relative growth rates, biomass and mean leaf sugar content were linearly related to [CO₂] over the entire experimental range. Photosynthesis appears to be regulated by stomata at low CO₂ levels and by non‐stomatal (i.e. biochemical limitations) at greater concentrations. In general, our results suggest that growth and physiological performance of cycads have been severely compromised by declining [CO₂] during the Cenozoic Era, possibly contributing to the current rare and endangered status of this functional type.     

Evolutionary and ecological genomics of non-model plants

Dissecting evolutionary dynamics of ecologically important traits is a long-term challenge for biologists.Attempts to understand natural variation and molecular mechanisms have motivated a move from laboratory model systems to non-model systems in diverse natural environments.Next generation sequencing methods,along with an expansion of genomic resources and tools,have fostered new links between diverse disciplines,including molecular biology,evolution,ecology,and genomics.Great progress has been made in a few non-model wild plants,such as Arabidopsis relatives,monkey flowers,and wild sunflowers.Until recently,the lack of comprehensive genomic information has limited evolutionary and ecological studies to larger QTL (quantitative trait locus) regions rather than single gene resolution,and has hindered recognition of general patterns of natural variation and local adaptation.Further efforts in accumulating genomic data and developing bioinformatic and biostatistical tools are now poised to move this field forward.Integrative national and international collaborations and research communities are needed to facilitate development in the field of evolutionary and ecological genomics.     

Unexpected side effects in biocrust after treating non‐native plants using carbon addition

Carbon addition has been proposed as an alternative to herbicide and manual removal methods to treat non‐native plants and reduce non‐target effects of treatments (e.g. impacts on native plants; surface disturbance). On Mojave Desert pavement and biocrust substrates after experimental soil disturbance and carbon addition (1,263 g C/m² as sucrose), we observed declines in lichens and moss cover in sucrose‐treated plots. To further explore this unforeseen potential side effect of using carbon addition as a non‐native plant treatment, we conducted biocrust surveys 5 and 7 years after treatments, sampled surface soils to observe if treatments additionally affected soil filamentous cyanobacteria, and conducted laboratory trials testing the effects of different levels of sucrose on cyanobacteria and desert mosses. Sucrose addition to biocrust plots reduced lichen and moss cover by 33–78% and species richness by 40–80%. Sucrose reduced biocrust cover in biocrust plots to levels similarly detected in pavement plots (<1%). While cyanobacteria in the field did not appear to be affected by sucrose, laboratory tests showed negative effects of sucrose on both cyanobacteria and mosses. Cyanobacteria declined by 41% 1 month after exposure to 5.4 g C/m² equivalent solutions. We detected injury to photosynthesis in mosses after 96 hour exposure to 79–316 g C/m² equivalent solutions. Caution is warranted when using carbon addition, at least in the form and concentration of sucrose, as a treatment for reducing non‐native plants on sites where conserving biocrust is a goal.     

Evolution of trace amine associated receptor (TAAR) gene family in vertebrates: lineage-specific expansions and degradations of a second class of vertebrate chemosensory receptors expressed in the olfactory epithelium

The trace amine-associated receptors (TAARs) form a specific family of G protein-coupled receptors in vertebrates. TAARs were initially considered neurotransmitter receptors, but recent study showed that mouse TAARs function as chemosensory receptors in the olfactory epithelium. To clarify the evolutionary dynamics of the TAAR gene family in vertebrates, near-complete repertoires of TAAR genes and pseudogenes were identified from the genomic assemblies of 4 teleost fishes (zebrafish, fugu, stickleback, and medaka), western clawed frogs, chickens, 3 mammals (humans, mice, and opossum), and sea lampreys. Database searches revealed that fishes had many putatively functional TAAR genes (13-109 genes), whereas relatively small numbers of TAAR genes (3-22 genes) were identified in tetrapods. Phylogenetic analysis of these genes indicated that the TAAR gene family was subdivided into 5 subfamilies that diverged before the divergence of ray-finned fishes and tetrapods. In tetrapods, virtually all TAAR genes were located in 1 specific region of their genomes as a gene cluster; however, in fishes, TAAR genes were scattered throughout more than 2 genomic locations. This possibly reflects a whole-genome duplication that occurred in the common ancestor of ray-finned fishes. Expression analysis of zebrafish and stickleback TAAR genes revealed that many TAARs in these fishes were expressed in the olfactory organ, suggesting the relatively high importance of TAARs as chemosensory receptors in fishes. A possible evolutionary history of the vertebrate TAAR gene family was inferred from the phylogenetic and comparative genomic analyses.