Computational identification and analysis of immune-associated nucleotide gene family in Arabidopsis thaliana

GTP-binding proteins represent a ubiquitous regulatory mechanism in controlling growth and development in eukaryotes under normal and stress conditions. The IAN/GIMAP proteins belong to a novel family of functionally uncharacterized GTP-binding proteins expressed in both plant and vertebrate cells during anti-pathogenic responses. To gain novel insights into their roles in plants, we did genome-wide analysis of the IAN/GIMAP gene family. We identified 13 Arabidopsis IAN/GIMAP genes, which share similar gene structures and mostly reside in a tandem cluster on chromosomes. Sequence comparison reveals that these genes encode 26–52 kDa proteins with one GTP-binding domain and a conserved box unique to the family. Phylogenetic analysis suggests that the IAN/GIMAP genes of angiosperms and vertebrates may have evolved by independent gene duplication events. GENEVESTIGATOR sources were mined for comprehensive and comparative Arabidopsis IAN/GIMAP gene family expression analysis. These data reveal that IAN/GIMAPs exhibit diverse expression patterns during development and in response to external stimuli, indicating that these paralogous genes are likely involved in complex biological processes in Arabidopsis. Our present findings provide a basis for elucidating the novel GTPase family protein-mediated regulatory mechanisms in the future.     

The Evolutionary History of MAPL (Mitochondria-Associated Protein Ligase) and Other Eukaryotic BAM/GIDE Domain Proteins

MAPL (mitochondria-associated protein ligase, also called MULAN/GIDE/MUL1) is a multifunctional mitochondrial outer membrane protein found in human cells that contains a unique BAM (beside a membrane) domain and a C-terminal RING-finger domain. MAPL has been implicated in several processes that occur in animal cells such as NF-kB activation, innate immunity and antiviral signaling, suppression of PINK1/parkin defects, mitophagy in skeletal muscle, and caspase-dependent apoptosis. Previous studies demonstrated that the BAM domain is present in diverse organisms in which most of these processes do not occur, including plants, archaea, and bacteria. Thus the conserved function of MAPL and its BAM domain remains an open question. In order to gain insight into its conserved function, we investigated the evolutionary origins of MAPL by searching for homologues in predicted proteomes of diverse eukaryotes. We show that MAPL proteins with a conserved BAM-RING architecture are present in most animals, protists closely related to animals, a single species of fungus, and several multicellular plants and related green algae. Phylogenetic analysis demonstrated that eukaryotic MAPL proteins originate from a common ancestor and not from independent horizontal gene transfers from bacteria. We also determined that two independent duplications of MAPL occurred, one at the base of multicellular plants and another at the base of vertebrates. Although no other eukaryote genome examined contained a verifiable MAPL orthologue, BAM domain-containing proteins were identified in the protists Bigelowiella natans and Ectocarpus siliculosis. Phylogenetic analyses demonstrated that these proteins are more closely related to prokaryotic BAM proteins and therefore likely arose from independent horizontal gene transfers from bacteria. We conclude that MAPL proteins with BAM-RING architectures have been present in the holozoan and viridiplantae lineages since their very beginnings. Our work paves the way for future studies into MAPL function in alternative model organisms like Capsaspora owczarzaki and Chlamydomonas reinhardtii that will help to answer the question of MAPL’s ancestral function in ways that cannot be answered by studying animal cells alone.     

Molecular characterization of 60 isolated wheat MYB genes and analysis of their expression during abiotic stress

The proteins of the MYB superfamily play central roles in developmental processes and defence responses in plants. Sixty unique wheat MYB genes that contain full-length cDNA sequences were isolated. These 60 genes were grouped into three categories, namely one R1R2R3-MYB, 22 R2R3-MYBs, and 37 MYB-related members. The sequence composition of the R2 and R3 repeats was conserved among the 22 wheat R2R3-MYB proteins. Phylogenetic comparison of the members of this superfamily among wheat, rice, and Arabidopsis revealed that the putative functions of some wheat MYB proteins were clustered into the Arabidopsis functional clades. Tissue-specific expression profiles showed that most of the wheat MYB genes were expressed in all of the tissues examined, suggesting that wheat MYB genes take part in multiple cellular processes. The expression analysis during abiotic stress identified a group of MYB genes that respond to one or more stress treatments. The overexpression of a salt-inducible gene, TaMYB32, enhanced the tolerance to salt stress in transgenic Arabidopsis. This study is the first comprehensive study of the MYB gene family in Triticeae.     

Unraveling the Distribution and Evolution of miR156 targeted SPLs in Plants by Phylogenetic Analysis

Squamosa promoter binding protein like genes (SPLs) are critical during plant development and mostly regulated by miR156. However, little is known about phylogenetic distribution and evolutionary patterns of miR156 targeted SPLs. In this study, 183 SPLs from nine genome sequenced species representing algae, bryophytes, lycophyte, monocots, and eudicots were computationally analyzed. Our results showed that miR156 responsive elements (MREs) on SPLs were present in land plants but absent from unicellular green algae. Phylogenetic analysis revealed that miR156 targeted SPLs only distributed in group II not group I of land plants, suggesting they originated from a common ancestor. In addition, group II were further divided into seven subgroups (IIa IIg) and miR156 targeted SPLs distributed in some specific members of SPLs from six subgroups except subgroup IId. Such distribution pattern was well elucidated by gene structure evolution of miR156 targeted SPLs based on the correlation of phylogenetic classification and gene structure. They could suffer from the exon loss events combined with MREs loss during evolution. Moreover, gene duplication contributed to the abundance of miR156 targeted SPLs, which had significantly increased after angiosperms and lower plants split. With Arabidopsis as the model species, we found segmental and tandem gene duplications predominated during miR156 targeted SPLs expansion. Taken together, these results provide better insights in understanding the function diversity and evolution of miR156 targeted SPLs in plants.     

miR156靶基因SPL在植物中的系统分布和进化模式分析

Squamosa promoterbinding proteinlike genes （SPLs）在植物发育过程中具有重要作用。很多SPLs被miR156调节,然而,对于它们在植物中的系统分布和进化模式还知之甚少。本文对9个测序物种（藻类,苔藓,石松,单子叶和双子叶植物）的183个SPLs进行了生物信息学分析。结果表明miR156应答元件（MREs）仅在陆生植物SPLs中发现,藻类中不存在。系统进化分析显示陆生植物SPLs分为两大分支：group I和group II。 MiR156靶基因仅分布于group II,表明它们有着共同的祖先。Group II进一步分为7个亚支（IIaIIg）,miR156靶基因分布在除IId外的其余6个亚支的特定SPLs。系统分类与基因结构的相关性反映了SPL靶基因结构上的变化。在进化过程中,它们可能发生外显子的丢失且伴随MRE的丢失。另外,基因重复对SPL靶基因的丰度变化影响很大,尤其是被子植物与低等植物分歧后它们数量明显增加。以拟南芥为模式植物分析发现串联重复和片段重复是SPL靶基因扩张的主要机制。     

Genomic organization and evolutionary conservation of plant D-type cyclins

Plants contain more genes encoding core cell cycle regulators than other organisms but it is unclear whether these represent distinct functions. D-type cyclins (CYCD) play key roles in the G1-to-S-phase transition, and Arabidopsis (Arabidopsis thaliana) contains 10 CYCD genes in seven defined subgroups, six of which are conserved in rice (Oryza sativa). Here, we identify 22 CYCD genes in the poplar (Populus trichocarpa) genome and confirm that these six CYCD subgroups are conserved across higher plants, suggesting subgroup-specific functions. Different subgroups show gene number increases, with CYCD3 having three members in Arabidopsis, six in poplar, and a single representative in rice. All three species contain a single CYCD7 gene. Despite low overall sequence homology, we find remarkable conservation of intron/exon boundaries, because in most CYCD genes of plants and mammals, the first exon ends in the conserved cyclin signature. Only CYCD3 genes contain the complete cyclin box in a single exon, and this structure is conserved across angiosperms, again suggesting an early origin for the subgroup. The single CYCD gene of moss has a gene structure closely related to those of higher plants, sharing an identical exon/intron structure with several higher plant subgroups. However, green algae have CYCD genes structurally unrelated to higher plants. Conservation is also observed in the location of potential cyclin-dependent kinase phosphorylation sites within CYCD proteins. Subgroup structure is supported by conserved regulatory elements, particularly in the eudicot species, including conserved E2F regulatory sites within CYCD3 promoters. Global expression correlation analysis further supports distinct expression patterns for CYCD subgroups.     

Beyond PKA: Evolutionary and structural insights that define a docking and dimerization domain superfamily

Protein-interaction domains can create unique macromolecular complexes that drive evolutionary innovation. By combining bioinformatic and phylogenetic analyses with structural approaches, we have discovered that the docking and dimerization (D/D) domain of the PKA regulatory subunit is an ancient and conserved protein fold. An archetypal function of this module is to interact with A-kinase-anchoring proteins (AKAPs) that facilitate compartmentalization of this key cell-signaling enzyme. Homology searching reveals that D/D domain proteins comprise a superfamily with 18 members that function in a variety of molecular and cellular contexts. Further in silico analyses indicate that D/D domains segregate into subgroups on the basis of their similarity to type I or type II PKA regulatory subunits. The sperm autoantigenic protein 17 (SPA17) is a prototype of the type II or R2D2 subgroup that is conserved across metazoan phyla. We determined the crystal structure of an extended D/D domain from SPA17 (amino acids 1–75) at 1.72 Å resolution. This revealed a four-helix bundle-like configuration featuring terminal β-strands that can mediate higher order oligomerization. In solution, SPA17 forms both homodimers and tetramers and displays a weak affinity for AKAP18. Quantitative approaches reveal that AKAP18 binding occurs at nanomolar affinity when SPA17 heterodimerizes with the ropporin-1-like D/D protein. These findings expand the role of the D/D fold as a versatile protein-interaction element that maintains the integrity of macromolecular architectures within organelles such as motile cilia.     

miR156靶基因SPL在植物中的系统分布和进化模式分析

Squamosa promoter binding protein like genes (SPLs)在植物发育过程中具有重要作用。很多SPLs被miR156调节,然而,对于它们在植物中的系统分布和进化模式还知之甚少。本文对9个测序物种（藻类,苔藓,石松,单子叶和双子叶植物）的183个SPLs进行了生物信息学分析。结果表明miR156应答元件（MREs）仅在陆生植物SPLs中发现,藻类中不存在。系统进化分析显示陆生植物SPLs分为两大分支：group I和group II。 MiR156靶基因仅分布于group II,表明它们有着共同的祖先。Group II进一步分为7个亚支(IIa IIg),miR156靶基因分布在除IId外的其余6个亚支的特定SPLs。系统分类与基因结构的相关性反映了SPL靶基因结构上的变化。在进化过程中,它们可能发生外显子的丢失且伴随MRE的丢失。另外,基因重复对SPL靶基因的丰度变化影响很大,尤其是被子植物与低等植物分歧后它们数量明显增加。以拟南芥为模式植物分析发现串联重复和片段重复是SPL靶基因扩张的主要机制。     

Genome-wide identification and analysis of membrane-bound O-acyltransferase (MBOAT) gene family in plants

Membrane bound O-acyl transferase (MBOAT) family is composed of gene members encoding a variety of acyltransferase enzymes, which play important roles in plant acyl lipid metabolism. Here, we present the first genome-enabled identification and analysis of MBOAT gene models in plants. In total, we identified 136 plant MBOAT sequences from 14 plant species with complete genomes. Phylogenetic relationship analyses suggested the plant MBOAT gene models fell into four major groups, two of which likely encode enzymes of diacylglycerol acyltransferase 1 (DGAT1) and lysophospholipid acyltransferase (LPLAT), respectively, with one–three copies of paralogs present in each of the most plant species. A group of gene sequences, which are homologous to Saccharomyces cerevisiae glycerol uptake proteins (GUP), was identified in plants; copy numbers were conserved, with only one copy represented in each of the most plant species; analyses showed that residues essential for acyltransferases were more prone to be conserved than vertebrate orthologs. Among four groups, one was inferred to emerge in land plants and experience a rapid expansion in genomes of angiosperms, which suggested their important roles in adaptation of plants in lands. Sequence and phylogeny analyses indicated that genes in all four groups encode enzymes with acyltransferases. Comprehensive sequence identification of MBOAT family members and investigation into classification provide a complete picture of the MBOAT gene family in plants, and could shed light into enzymatic functions of different MBOAT genes in plants.     

Two new subfamilies of DNA mismatch repair proteins (MutS) specifically abundant in the marine environment

MutS proteins are ubiquitous in cellular organisms and have important roles in DNA mismatch repair or recombination. In the virus world, the amoeba-infecting Mimivirus, as well as the recently sequenced Cafeteria roenbergensis virus are known to encode a MutS related to the homologs found in octocorals and ɛ-proteobacteria. To explore the presence of MutS proteins in other viral genomes, we performed a genomic survey of four giant viruses (‘giruses'') (Pyramimonas orientalis virus (PoV), Phaeocystis pouchetii virus (PpV), Chrysochromulina ericina virus (CeV) and Heterocapsa circularisquama DNA virus (HcDNAV)) that infect unicellular marine algae. Our analysis revealed the presence of a close homolog of Mimivirus MutS in all the analyzed giruses. These viral homologs possess a specific domain structure, including a C-terminal HNH-endonuclease domain, defining the new MutS7 subfamily. We confirmed the presence of conserved mismatch recognition residues in all members of the MutS7 subfamily, suggesting their role in DNA mismatch repair rather than DNA recombination. PoV and PpV were found to contain an additional type of MutS, which we propose to call MutS8. The MutS8 proteins in PoV and PpV were found to be closely related to homologs from ‘Candidatus Amoebophilus asiaticus'', an obligate intracellular amoeba-symbiont belonging to the Bacteroidetes. Furthermore, our analysis revealed that MutS7 and MutS8 are abundant in marine microbial metagenomes and that a vast majority of these environmental sequences are likely of girus origin. Giruses thus seem to represent a major source of the underexplored diversity of the MutS family in the microbial world.     

The complete chloroplast genomes of three Hamamelidaceae species: Comparative and phylogenetic analyses

Hamamelidaceae is an important group that represents the origin and early evolution of angiosperms. Its plants have many uses, such as timber, medical, spice, and ornamental uses. In this study, the complete chloroplast genomes of Loropetalum chinense (R. Br.) Oliver, Corylopsis glandulifera Hemsl., and Corylopsis velutina Hand.‐Mazz. were sequenced using the Illumina NovaSeq 6000 platform. The sizes of the three chloroplast genomes were 159,402 bp (C. glandulifera), 159,414 bp (C. velutina), and 159,444 bp (L. chinense), respectively. These chloroplast genomes contained typical quadripartite structures with a pair of inverted repeat (IR) regions (26,283, 26,283, and 26,257 bp), a large single‐copy (LSC) region (88,134, 88,146, and 88,160 bp), and a small single‐copy (SSC) region (18,702, 18,702, and 18,770 bp). The chloroplast genomes encoded 132–133 genes, including 85–87 protein‐coding genes, 37–38 tRNA genes, and 8 rRNA genes. The coding regions were composed of 26,797, 26,574, and 26,415 codons, respectively, most of which ended in A/U. A total of 37–43 long repeats and 175–178 simple sequence repeats (SSRs) were identified, and the SSRs contained a higher number of A + T than G + C bases. The genome comparison showed that the IR regions were more conserved than the LSC or SSC regions, while the noncoding regions contained higher variability than the gene coding regions. Phylogenetic analyses revealed that species in the same genus tended to cluster together. Chunia Hung T. Chang, Mytilaria Lecomte, and Disanthus Maxim. may have diverged early and Corylopsis Siebold & Zucc. was closely related to Loropetalum R. Br. This study provides valuable information for further species identification, evolution, and phylogenetic studies of Hamamelidaceae plants.     

Chromosome‐level genome assembly for the horned‐gall aphid provides insights into interactions between gall‐making insect and its host plant

The aphid Schlechtendalia chinensis is an economically important insect that can induce horned galls, which are valuable for the medicinal and chemical industries. Up to now, more than twenty aphid genomes have been reported. Most of the sequenced genomes are derived from free‐living aphids. Here, we generated a high‐quality genome assembly from a galling aphid. The final genome assembly is 271.52 Mb, representing one of the smallest sequenced genomes of aphids. The genome assembly is based on contig and scaffold N50 values of the genome sequence are 3.77 Mb and 20.41 Mb, respectively. Nine‐seven percent of the assembled sequences was anchored onto 13 chromosomes. Based on BUSCO analysis, the assembly involved 96.9% of conserved arthropod and 98.5% of the conserved Hemiptera single‐copy orthologous genes. A total of 14,089 protein‐coding genes were predicted. Phylogenetic analysis revealed that S. chinensis diverged from the common ancestor of Eriosoma lanigerum approximately 57 million years ago (MYA). In addition, 35 genes encoding salivary gland proteins showed differentially when S. chinensis forms a gall, suggesting they have potential roles in gall formation and plant defense suppression. Taken together, this high‐quality S. chinensis genome assembly and annotation provide a solid genetic foundation for future research to reveal the mechanism of gall formation and to explore the interaction between aphids and their host plants.