全基因组复制事件的绝对定年揭示莲座蕨属植物的迟滞演化

全基因组复制在维管植物的物种形成过程中普遍存在, 被认为是物种适应极端环境的重要机制之一。确定全基因组复制事件的发生时间对理解生物的适应性演化具有重要意义。然而, 在维管植物, 特别是蕨类植物中, 全基因组复制事件的发生时间及其演化意义仍知之甚少。本研究以蕨类植物重要基部类群——福建莲座蕨(Angiopteris fokiensis)为例, 基于不同采样点(广东、广西、上海)的3个转录组学数据, 利用同义替换率(Ks)和绝对定年的方法分析全基因组复制事件的发生时间和物种单位时间内的分子演化速率, 并对事件发生后保留下的基因进行基因功能注释和富集分析。结果表明, 福建莲座蕨在159‒165 Mya发生了一次全基因组复制事件, 该复制事件优先保留的基因主要与营养代谢、信号传导、适应调节和组织结构生长相关。另外, 福建莲座蕨的分子演化速率为1.66 × 10^‒9 (同义替换/位点/年), 是除裸子植物外, 陆生植物中已知演化速率最缓慢的类群。综合以上研究结果, 我们推测福建莲座蕨全基因组复制的发生可能与裸子植物繁盛、核心被子植物集中兴起或托阿尔阶灭绝事件有关。而复制后显著保留基因可能促进了莲座蕨属(Angiopteris)植物的遗传和形态创新, 从而帮助其快速适应环境的剧烈变化。进一步对该类群植物演化速率缓慢的原因进行讨论, 推测莲座蕨属缓慢的演化速率可能与其本身世代周期长、基因组较大及其生长环境稳定有关。本研究通过分析福建莲座蕨的全基因组复制历史和复制基因的保留模式, 推测全基因组复制事件对促进演化速率较慢的植物适应极端环境变化具有重要意义, 可为理解其他陆生植物的适应性演化提供更多启发。     

基于转录组数据揭示4种兜兰的全基因组复制历史

多倍化或全基因组复制(WGD)是物种多样性发生的重要驱动力。目前, 在蕨类、菊科以及豆科等类群丰富的植物中已多次报道全基因组复制事件, 而兰科(Orchidaceae)全基因组复制事件报道极少, 与其丰富的物种多样性存在矛盾, 推测与前期样本量小但类群跨度大的研究策略有关。选取染色体数目变异丰富且多样性较高的兜兰属(Paphiopedilum)为兰科植物代表类群, 基于共享数据库中4种兜兰的转录组数据, 采用同义替换率(K_s)、系统发生基因组学以及相对定年的方法分析兜兰属植物是否发生过全基因组复制事件。结果表明, 在4种兜兰中均检测到3次全基因组复制事件, 分别发生在110.17-119.77 Mya (WGD1)、60.95-74.19 Mya (WGD2)和38.19-45.85 Mya (WGD3)。其中, WGD3为新检测到的全基因组复制事件, 推测其发生在杓兰亚科(Cypripedioideae)与姐妹类群分化后, 兜兰属与姐妹类群分化之前。此外, 3次全基因组复制事件发生后优先保留的基因拷贝在功能上多与当时的环境胁迫响应相关, 推测全基因组复制提高了兜兰属植物祖先对当时极端环境变化的适应性。     

葡萄、桃和可可基因组的进化分析

本研究以葡萄、桃和可可为研究对象,基于比较基因组学,利用基因同源共线性方法对基因组内的结构和基因组间同源信息进行比对分析,确定了物种基因组内和基因组间的同源片段。通过统计3个物种基因组间的同源共线基因的保留情况发现,葡萄基因组的保留情况最好,桃次之(为73.4%),可可最差(为68.9%),其丢失均可能是由于双子叶植物共有的三倍化导致基因组稳定性遭到破坏。另外,共线基因间的同义核苷酸替换率的频数分布证实,葡萄、桃和可可仅经历过一次古老的全基因组三倍化,并未经历最近的全基因组加倍,且可可基因组进化最快,葡萄基因组进化最保守;3个物种的分歧时间分别为:葡萄(~110 Mya)、可可(~90 Mya)、桃(~80 Mya)。本研究将为3个物种及双子叶植物基因组的结构、功能和进化等研究提供重要的理论依据。     

水稻和其他禾本科植物基因组多倍体起源的证据

基因加倍(Gene duplication)被认为是进化的加速器。古老的基因组加倍事件已经在多个物种中被确定,包括酵母、脊椎动物以及拟南芥等。本研究发现水稻基因组同样存在全基因组加倍事件,大概发生在禾谷类作物分化之前,距今约7000万年。在水稻基因组中,共找到117个加倍区段(Duplicated block),分布在水稻的全部12条染色体,覆盖约60％的水稻基因组。在加倍区段,大约有20％的基因保留了加倍后的姊妹基因对(Duplicated pairs)。与此形成鲜明对照的是加倍区段的转录因子保留了60％的姊妹基因。禾本科植物全基因组加倍事件的确定对研究禾本科植物基因组的进化具有重要影响,暗示了多倍体化及随后的基因丢失、染色体重排等在禾谷类物种分化中扮演了重要角色。     

基因重复研究进展

基因复制是基因通过不等交换,反转录转座或由全基因组复制等途径产生一个与原基因相似的基因或碱基序列,它与生物体基因组大小的进化、新基因的起源、物种的分化以及基因抗突变的能力大小等都密切相关。本文综述了复制基因的产生和保留机制、选择作用、分化的途径以及复制基因进化速率等方面的相关研究,揭示了基因复制对于生物进化的重要性,以引起大家对该领域的了解与关注。关键词：基因复制;复制基因;不等交换;反转录转座;全基因组复制     

全基因组与串联复制后白菜基因的保留

全基因组复制与串联复制是两种重要的基因扩增途径,在生物进化过程中普遍存在.这两种复制方式相互关系的研究在拟南芥中已经取得很多成果.白菜(Brassica rapa)属于十字花科(Brassicaceae)芸薹属(Brassca),是一类重要的经济作物,也是研究基因组多倍化和形态演化的模式植物.白菜基因组的测序与组装工作已经取得了重大成就,运用比较基因组学的方法,通过比较白菜与模式植物拟南芥,可以清晰鉴定白菜基因组经历的全基因组三倍化事件.同时,白菜与拟南芥同属于十字花科,有较近的起源关系和良好的基因组共线性关系.因此,拟南芥可以作为外群研究白菜全基因组三倍化以及串联重复之后基因的偏向性保留.结果发现,在白菜中存在物种特有的偏向性保留基因,即与环境刺激相关的基因和与激素相关的基因.     

假基因的组成、分布及其分子进化

假基因(pseudogene)是指基因组中与正常基因序列相似, 但是缺乏功能的DNA 序列。通过序列同源性搜索, 可以收集基因组中假基因的群体特性、染色体分布和同源家族等特性。假基因很好地保留了数百万年前基因组中祖先基因的分子记录, 被视为“基因化石”, 因此假基因在进化和比较基因组学中是重要的资源。应用假基因和基因比较体系, 可以探究生物基因的进化史和基因组稳定性。如: 用Ka/Ks比值确定假基因的自然选择压、物种亲缘关系和进化距离, 分析假基因自身的进化趋势, 探讨DNA 突变的成因等。     

槲蕨属叶绿体基因组密码子偏好性分析

叶绿体基因组密码子偏好性使用模式往往影响基因表达效率,为促进药用植物叶绿体基因工程的发展,提高槲蕨药用品质,该研究以川滇槲蕨、栎叶槲蕨和槲蕨三个近缘药用植物为材料,使用CodonW、CUSP和SPSS等软件分析其叶绿体基因组编码基因密码子使用偏好性,筛选出三个物种的最优密码子。结果表明:川滇槲蕨、栎叶槲蕨和槲蕨叶绿体基因组的有效密码子数(ENC)范围分别为40.10~61、40.33~61和40.15~61,其密码子偏好性较弱;ndhE、rpl22、rpl14、rpl20、ccsA、rps4和rpl16编码基因的ENC值差异较大,表明近缘物种中,部分基因的密码子偏好性存在一定差异;三个物种编码基因的ENC频数集中于-0.1~0.1之间,说明槲蕨属基因密码子偏好性主要受到突变的影响。川滇槲蕨12个最优密码子有6个和栎叶槲蕨相同,分别是UCU、ACU、GCU、CAA、AAA和GAU;栎叶槲蕨10个最优密码子有2个和槲蕨相同,分别是UUA和AUU,而川滇槲蕨与槲蕨无相同的最优密码子。该研究结果可为槲蕨属药用植物基因工程中外源基因的改良及其表达奠定基础。     

基于孢子形态和分子证据探讨鳞盖蕨属(碗蕨科)系统分类

孢粉学是解决植物分类中疑难类群物种微形态分化的重要方法, 随着分子系统学的发展, 结合这两门学科的优势可以更加有效地解决疑难类群的分类学问题。鳞盖蕨属(Microlepia)是一个分类困难的疑难类群, 采用孢粉学与分子系统学一一对应的方法, 以及居群取样方式, 选取280份样本, 联合4个叶绿体片段(rbcL、trnL-F、psbA-trnH和rps4), 采用最大似然法和贝叶斯法构建该属的系统发生关系, 在此基础上对凭证标本中100份材料的孢子进行观察和分析。综合分子系统学和孢粉学的研究结果, 得出结论: (1) 在形态学研究中广泛被接受的15个物种得到了单系支持, 并厘清了分类困难的复合群; (2) 发现边缘鳞盖蕨(M. marginata)可能存在隐性种; (3) 建议恢复过去归并处理为异名的瑶山鳞盖蕨(M. yaoshanica)、罗浮鳞盖蕨(M. lofoushanensis)、四川鳞盖蕨(M. szechuanica)以及滇西鳞盖蕨(M. subspeluncae); (4) 提出鳞盖蕨属可能存在杂交现象; (5) 提出鳞盖蕨属完整的属下分类建议。     

甘蓝型油菜MADS-box基因家族的鉴定与系统进化分析

MADS-box基因家族参与调控开花时间、花器官分化、根系生长、分生组织分化、子房和配子发育、果实膨大及衰老等植物生长发育的重要过程。基于甘蓝型油菜(Brassica napus)基因组测序数据,利用生物信息学方法对甘蓝型油菜MADS-box基因家族进行鉴定和注释及基因结构与系统进化分析。结果显示,在甘蓝型油菜中鉴定出307个MADS-box基因家族成员,根据进化关系可将其分为两大类型,I型(M-type)包含α、β、γ三个亚家族,II型(MIKC-type)包括MIKCC和MIKC*两个亚家族,MIKCC可进一步分为13个小类;甘蓝型油菜A基因组染色体上分布的MADS-box基因多于C基因组。在基因结构上,MIKC-type亚家族基因序列普遍比M-type长且含有较多的外显子;M-type亚家族蛋白序列中的motif数量为2–5个,MIKC-type亚家族蛋白序列中平均含有7个motif。拟南芥(Arabidopsis thaliana)与甘蓝型油菜MADS-box基因共线性分析结果显示,全基因组复制事件对MADS-box基因家族尤其是MIKC亚家族的扩张起重要作用;MIKC亚家族基因在进化过程中受到的选择压力约为M-type的2倍,这表明MIKC-type亚家族在进化过程中被选择性保留。     

Recurrent genome duplication events likely contributed to both the ancient and recent rise of ferns

Ferns, the second largest group of vascular plants, originated ～400 mil ion years ago(Mya). They became dominant in the ancient Earth landscape before the angiosperms and are stil important in current ecosystems.Many ferns have exceptional y high chromosome numbers,possibly resulting from whole-genome duplications(WGDs).However, WGDs have not been investigated molecularly across fern diversity. Here we detected and dated fern WGDs using a phylogenomic approach and by calculating synonymous substitution rates(Ks). We also investigated a possible correlation between proposed WGDs and shifts in species diversification rates. We identified 19 WGDs: three ancient events along the fern phylogenetic backbone that are shared by 66%–97% of extant ferns, with additional lineage-specific WGDs for eight orders, providing strongevidence for recurring genome duplications across fern evolutionary history. We also observed similar Ks peak values for more than half of these WGDs, with multiple WGDs occurring close to the Cretaceous(～145–66 Mya). Despite the repeated WGD events, the biodiversity of ferns declined during the Cretaceous, implying that other factors probably contributed to the floristic turnover from ferns to angiosperms. This study provides molecular evidence for recurring WGDs in ferns and offers important clues to the genomic evolutionary history of ferns.     

Evolution of the Relaxin/Insulin-like Gene Family in Placental Mammals: Implications for Its Early Evolution

The relaxin (RLN) and insulin-like (INSL) gene family is a group of genes involved in a variety of physiological roles that includes bone formation, testicular descent, trophoblast development, and cell differentiation. This family appears to have expanded in vertebrates relative to non-vertebrate chordates, but the relative contribution of whole genome duplications (WGDs) and tandem duplications to the observed diversity of genes is still an open question. Results from our comparative analyses favor a model of divergence post vertebrate WGDs in which a single-copy progenitor found in the last common ancestor of vertebrates experienced two rounds of WGDs before the functional differentiation that gave rise to the RLN and INSL genes. One of the resulting paralogs was subsequently lost, resulting in three proto-RLN/INSL genes on three separate chromosomes. Subsequent rounds of tandem gene duplication and divergence originated the set of paralogs found on a given cluster in extant vertebrates. Our study supports the hypothesis that differentiation of the RLN and INSL genes took place independently in each RLN/INSL cluster after the two WGDs during the evolutionary history of vertebrates. In addition, we show that INSL4 represents a relatively old gene that has been apparently lost independently in all Euarchontoglires other than apes and Old World monkeys, and that RLN2 derives from an ape-specific duplication.     

Tangled up in two: a burst of genome duplications at the end of the Cretaceous and the consequences for plant evolution

Genome sequencing has demonstrated that besides frequent small-scale duplications, large-scale duplication events such as whole genome duplications (WGDs) are found on many branches of the evolutionary tree of life. Especially in the plant lineage, there is evidence for recurrent WGDs, and the ancestor of all angiosperms was in fact most likely a polyploid species. The number of WGDs found in sequenced plant genomes allows us to investigate questions about the roles of WGDs that were hitherto impossible to address. An intriguing observation is that many plant WGDs seem associated with periods of increased environmental stress and/or fluctuations, a trend that is evident for both present-day polyploids and palaeopolyploids formed around the Cretaceous–Palaeogene (K–Pg) extinction at 66 Ma. Here, we revisit the WGDs in plants that mark the K–Pg boundary, and discuss some specific examples of biological innovations and/or diversifications that may be linked to these WGDs. We review evidence for the processes that could have contributed to increased polyploid establishment at the K–Pg boundary, and discuss the implications on subsequent plant evolution in the Cenozoic.     

Preferential retention of the slowly evolving gene in pairs of duplicates in angiosperm genomes

Gene duplication provides raw material for functional innovation, but gene duplicability varies considerably. Previous studies have found widespread asymmetrical sequence evolution between paralogs. However, it remains unknown whether the rate of evolution among paralogs affects their propensity of being retained after another round of whole-genome duplication (WGD). In this study, we investigated gene groups that have experienced two successive WGDs to determine which of two older duplicates with different evolutionary rates was more likely to retain both younger duplicates. To uncouple the measurement of evolutionary rates from any assignment of duplicate or singleton status, we measured the evolutionary rates of singleton genes in out-lineages but classified these singleton genes according to whether they are retained or not in a crown group of species. We found that genes that retained younger duplicates in the crown group of genomes were more constrained prior to the younger duplication event than those that failed to leave duplicates. In addition, we also found that the retained clades have more genes in out-lineages. Subsequent analyses showed that genes in the retained clades were expressed more broadly and highly than genes in the singleton clades. We concluded that the set of repeatedly retained genes after two WGDs is biased toward slowly evolving genes in angiosperms, suggesting that the potential of genes for both functional conservation and divergence likely affects their propensity of being retained after WGD in angiosperms.     

Cultivated hawthorn (Crataegus pinnatifida var. major) genome sheds light on the evolution of Maleae (apple tribe)

Cultivated hawthorn (Crataegus pinnatifida var. major) is an important medicinal and edible plant with a long history of use for health protection in China. Herein, we provide a de novo chromosome-level genome sequence of the hawthorn cultivar “Qiu Jinxing.” We assembled an 823.41 Mb genome encoding 40 571 genes and further anchored the 779.24 Mb sequence into 17 pseudo-chromosomes, which account for 94.64% of the assembled genome. Phylogenomic analyses revealed that cultivated hawthorn diverged from other species within the Maleae (apple tribe) at approximately 35.4 Mya. Notably, genes involved in the flavonoid and triterpenoid biosynthetic pathways have been significantly amplified in the hawthorn genome. In addition, our results indicated that the Maleae share a unique ancient tetraploidization event; however, no recent independent whole-genome duplication event was specifically detected in hawthorn. The amplification of non-specific long terminal repeat retrotransposons contributed the most to the expansion of the hawthorn genome. Furthermore, we identified two paleo-sub-genomes in extant species of Maleae and found that these two sub-genomes showed different rearrangement mechanisms. We also reconstructed the ancestral chromosomes of Rosaceae and discussed two possible paleo-polyploid origin patterns (autopolyploidization or allopolyploidization) of Maleae. Overall, our study provides an improved context for understanding the evolution of Maleae species, and this new high-quality reference genome provides a useful resource for the horticultural improvement of hawthorn.     

Screening synteny blocks in pairwise genome comparisons through integer programming

Background  

It is difficult to accurately interpret chromosomal correspondences such as true orthology and paralogy due to significant divergence of genomes from a common ancestor. Analyses are particularly problematic among lineages that have repeatedly experienced whole genome duplication (WGD) events. To compare multiple "subgenomes" derived from genome duplications, we need to relax the traditional requirements of "one-to-one" syntenic matchings of genomic regions in order to reflect "one-to-many" or more generally "many-to-many" matchings. However this relaxation may result in the identification of synteny blocks that are derived from ancient shared WGDs that are not of interest. For many downstream analyses, we need to eliminate weak, low scoring alignments from pairwise genome comparisons. Our goal is to objectively select subset of synteny blocks whose total scores are maximized while respecting the duplication history of the genomes in comparison. We call this "quota-based" screening of synteny blocks in order to appropriately fill a quota of syntenic relationships within one genome or between two genomes having WGD events.     

Identification of paralogous HERV-K LTRs on human chromosomes 3, 4, 7 and 11 in regions containing clusters of olfactory receptor genes

A locus harboring a human endogenous retroviral LTR (long terminal repeat) was mapped on the short arm of human chromosome 7 (7p22), and its evolutionary history was investigated. Sequences of two human genome fragments that were homologous to the LTR-flanking sequences were found in human genome databases: (1) an LTR-containing DNA fragment from region 3p13 of the human genome, which includes clusters of olfactory receptor genes and pseudogenes; and (2) a fragment of region 21q22.1 lacking LTR sequences. PCR analysis demonstrated that LTRs with highly homologous flanking sequences could be found in the genomes of human, chimp, gorilla, and orangutan, but were absent from the genomes of gibbon and New World monkeys. A PCR assay with a primer set corresponding to the sequence from human Chr 3 allowed us to detect LTR-containing paralogous sequences on human chromosomes 3, 4, 7, and 11. The divergence times for the LTR-flanking sequences on chromosomes 3 and 7, and the paralogous sequence on chromosome 21, were evaluated and used to reconstruct the order of duplication events and retroviral insertions. (1) An initial duplication event that occurred 14-17 Mya and before LTR insertion - produced two loci, one corresponding to that located on Chr 21, while the second was the ancestor of the loci on chromosomes 3 and 7. (2) Insertion of the LTR (most probably as a provirus) into this ancestral locus took place 13 Mya. (3) Duplication of the LTR-containing ancestral locus occurred 11 Mya, forming the paralogous modern loci on Chr 3 and 7.     

Mining EST databases to resolve evolutionary events in major crop species.

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.     

Screening of duplicated loci reveals hidden divergence patterns in a complex salmonid genome

A whole‐genome duplication (WGD) doubles the entire genomic content of a species and is thought to have catalysed adaptive radiation in some polyploid‐origin lineages. However, little is known about general consequences of a WGD because gene duplicates (i.e., paralogs) are commonly filtered in genomic studies; such filtering may remove substantial portions of the genome in data sets from polyploid‐origin species. We demonstrate a new method that enables genome‐wide scans for signatures of selection at both nonduplicated and duplicated loci by taking locus‐specific copy number into account. We apply this method to RAD sequence data from different ecotypes of a polyploid‐origin salmonid (Oncorhynchus nerka) and reveal signatures of divergent selection that would have been missed if duplicated loci were filtered. We also find conserved signatures of elevated divergence at pairs of homeologous chromosomes with residual tetrasomic inheritance, suggesting that joint evolution of some nondiverged gene duplicates may affect the adaptive potential of these genes. These findings illustrate that including duplicated loci in genomic analyses enables novel insights into the evolutionary consequences of WGDs and local segmental gene duplications.