茉莉酸信号受体COI蛋白家族的分子进化与基因表达分析

茉莉酸(Jasmonic acid,JA)存在于所有高等植物中,是植物对病原微生物和虫害防御反应的关键激素。在茉莉酸信号转导中,COI1(COR-insensitive 1)作为茉莉酸信号受体蛋白在其中发挥关键作用。本研究采用生物信息学方法,从藻类、苔藓类、蕨类、裸子及单、双子叶植物多谱系对COI蛋白家族进行比较基因组学研究,并取得以下结果:(1)同源基因鉴定结果发现,在所选的7种陆生植物中一共鉴定了55个COIs同源基因,然而,在低等的水生植物包括绿藻类(Chlorophytes)、红藻类(Rhodophytes)、硅藻类(Bacillariophytes)、灰胞藻类(Glaucophytes)及褐藻类(Phaeophytes)等基因组中均未发现其同源基因;(2)系统进化树分析表明,植物COI蛋白家族可以分为4个保守的亚家族,且在陆生植物扩增的同时可能已发生功能分化;(3)基因结构分析显示,植物COI家族基因结构表现多样性,主要体现在内含子的数目和长度上;(4)基因表达数据提示,COI基因家族成员参与植物生长发育的多个时期,且在不同组织器官以及不同的胁迫应答反应中发挥不同的作用。以上结果将为植物COI基因家族的深入研究提供参考。     

Pairs of Mutually Compensatory Frameshifting Mutations Contribute to Protein Evolution

Insertions and deletions of lengths not divisible by 3 in protein-coding sequences cause frameshifts that usually induce premature stop codons and may carry a high fitness cost. However, this cost can be partially offset by a second compensatory indel restoring the reading frame. The role of such pairs of compensatory frameshifting mutations (pCFMs) in evolution has not been studied systematically. Here, we use whole-genome alignments of protein-coding genes of 100 vertebrate species, and of 122 insect species, studying the prevalence of pCFMs in their divergence. We detect a total of 624 candidate pCFM genes; six of them pass stringent quality filtering, including three human genes: RAB36, ARHGAP6, and NCR3LG1. In some instances, amino acid substitutions closely predating or following pCFMs restored the biochemical similarity of the frameshifted segment to the ancestral amino acid sequence, possibly reducing or negating the fitness cost of the pCFM. Typically, however, the biochemical similarity of the frameshifted sequence to the ancestral one was not higher than the similarity of a random sequence of a protein-coding gene to its frameshifted version, indicating that pCFMs can uncover radically novel regions of protein space. In total, pCFMs represent an appreciable and previously overlooked source of novel variation in amino acid sequences.     

Rethinking Opsins

Opsins, the protein moieties of animal visual photo-pigments, have emerged as moonlighting proteins with diverse, light-dependent and -independent physiological functions. This raises the need to revise some basic assumptions concerning opsin expression, structure, classification, and evolution.     

Remarkably Low KIR and HLA Diversity in Amerindians Reveals Signatures of Strong Purifying Selection Shaping the Centromeric KIR Region

The killer-cell immunoglobulin-like receptors (KIR) recognize human leukocyte antigen (HLA) molecules to regulate the cytotoxic and inflammatory responses of natural killer cells. KIR genes are encoded by a rapidly evolving gene family on chromosome 19 and present an unusual variation of presence and absence of genes and high allelic diversity. Although many studies have associated KIR polymorphism with susceptibility to several diseases over the last decades, the high-resolution allele-level haplotypes have only recently started to be described in populations. Here, we use a highly innovative custom next-generation sequencing method that provides a state-of-art characterization of KIR and HLA diversity in 706 individuals from eight unique South American populations: five Amerindian populations from Brazil (three Guarani and two Kaingang); one Amerindian population from Paraguay (Aché); and two urban populations from Southern Brazil (European and Japanese descendants from Curitiba). For the first time, we describe complete high-resolution KIR haplotypes in South American populations, exploring copy number, linkage disequilibrium, and KIR–HLA interactions. We show that all Amerindians analyzed to date exhibit the lowest numbers of KIR–HLA interactions among all described worldwide populations, and that 83–97% of their KIR–HLA interactions rely on a few HLA-C molecules. Using multiple approaches, we found signatures of strong purifying selection on the KIR centromeric region, which codes for the strongest NK cell educator receptors, possibly driven by the limited HLA diversity in these populations. Our study expands the current knowledge of KIR genetic diversity in populations to understand KIR–HLA coevolution and its impact on human health and survival.     

Lower Statistical Support with Larger Data Sets: Insights from the Ochrophyta Radiation

It is commonly assumed that increasing the number of characters has the potential to resolve evolutionary radiations. Here, we studied photosynthetic stramenopiles (Ochrophyta) using alignments of heterogeneous origin mitochondrion, plastid, and nucleus. Surprisingly while statistical support for the relationships between the six major Ochrophyta lineages increases when comparing the mitochondrion (6,762 sites) and plastid (21,692 sites) trees, it decreases in the nuclear (209,105 sites) tree. Statistical support is not simply related to the data set size but also to the quantity of phylogenetic signal available at each position and our ability to extract it. Here, we show that this ability for current phylogenetic methods is limited, because conflicting results were obtained when varying taxon sampling. Even though the use of a better fitting model improved signal extraction and reduced the observed conflicts, the plastid data set provided higher statistical support for the ochrophyte radiation than the larger nucleus data set. We propose that the higher support observed in the plastid tree is due to an acceleration of the evolutionary rate in one short deep internal branch, implying that more phylogenetic signal per position is available to resolve the Ochrophyta radiation in the plastid than in the nuclear data set. Our work therefore suggests that, in order to resolve radiations, beyond the obvious use of data sets with more positions, we need to continue developing models of sequence evolution that better extract the phylogenetic signal and design methods to search for genes/characters that contain more signal specifically for short internal branches.     

Variation in Release Factor Abundance Is Not Needed to Explain Trends in Bacterial Stop Codon Usage

In bacteria stop codons are recognized by one of two class I release factors (RF1) recognizing TAG, RF2 recognizing TGA, and TAA being recognized by both. Variation across bacteria in the relative abundance of RF1 and RF2 is thus hypothesized to select for different TGA/TAG usage. This has been supported by correlations between TAG:TGA ratios and RF1:RF2 ratios across multiple bacterial species, potentially also explaining why TAG usage is approximately constant despite extensive variation in GC content. It is, however, possible that stop codon trends are determined by other forces and that RF ratios adapt to stop codon usage, rather than vice versa. Here, we determine which direction of the causal arrow is the more parsimonious. Our results support the notion that RF1/RF2 ratios become adapted to stop codon usage as the same trends, notably the anomalous TAG behavior, are seen in contexts where RF1:RF2 ratios cannot be, or are unlikely to be, causative, that is, at 3′untranslated sites never used for translation termination, in intragenomic analyses, and across archaeal species (that possess only one RF1). We conclude that specifics of RF biology are unlikely to fully explain TGA/TAG relative usage. We discuss why the causal relationships for the evolution of synonymous stop codon usage might be different from those affecting synonymous sense codon usage, noting that transitions between TGA and TAG require two-point mutations one of which is likely to be deleterious.     

Current biogeographical roles of the Kunlun Mountains

Large‐scale patterns of biodiversity and formation have garnered increasing attention in biogeography and macroecology. The Qinghai‐Tibet Plateau (QTP) is an ideal area for exploring these issues. However, the QTP consists of multiple geographic subunits, which are understudied. The Kunlun Mountains is a geographical subunit situated in the northern edge of the QTP, in northwest China. The diversity pattern, community phylogenetic structures, and biogeographical roles of the current flora of the Kunlun Mountains were analyzed by collecting and integrating plant distribution, regional geological evolution, and phylogeography. A total of 1911 species, 397 genera, and 75 families present on the Kunlun Mountains, of which 29.8% of the seed plants were endemic to China. The mean divergence time (MDT) of the Kunlun Mountains flora was in the early Miocene (19.40 Ma). Analysis of plant diversity and MDT indicated that the eastern regions of the Kunlun Mountains were the center of species richness, endemic taxa, and ancient taxa. Geographical origins analysis showed that the Kunlun Mountains flora was diverse and that numerous clades were from East Asia and Tethyan. Analysis of geographical origins and geological history together highlighted that the extant biodiversity on the Kunlun Mountains appeared through species recolonization after climatic fluctuations and glaciations during the Quaternary. The nearest taxon index speculated that habitat filtering was the most important driving force for biodiversity patterns. These results suggest that the biogeographical roles of the Kunlun Mountains are corridor and sink, and the corresponding key processes are species extinction and immigration. The Kunlun Mountains also form a barrier, representing a boundary among multiple floras, and convert the Qinghai‐Tibet Plateau into a relatively closed geographical unit.     

Hundreds of Out-of-Frame Remodeled Gene Families in the Escherichia coli Pangenome

All genomes include gene families with very limited taxonomic distributions that potentially represent new genes and innovations in protein-coding sequence, raising questions on the origins of such genes. Some of these genes are hypothesized to have formed de novo, from noncoding sequences, and recent work has begun to elucidate the processes by which de novo gene formation can occur. A special case of de novo gene formation, overprinting, describes the origin of new genes from noncoding alternative reading frames of existing open reading frames (ORFs). We argue that additionally, out-of-frame gene fission/fusion events of alternative reading frames of ORFs and out-of-frame lateral gene transfers could contribute to the origin of new gene families. To demonstrate this, we developed an original pattern-search in sequence similarity networks, enhancing the use of these graphs, commonly used to detect in-frame remodeled genes. We applied this approach to gene families in 524 complete genomes of Escherichia coli. We identified 767 gene families whose evolutionary history likely included at least one out-of-frame remodeling event. These genes with out-of-frame components represent ∼2.5% of all genes in the E. coli pangenome, suggesting that alternative reading frames of existing ORFs can contribute to a significant proportion of de novo genes in bacteria.