Comprehensive Evolutionary and Expression Analysis of FCS-Like Zinc finger Gene Family Yields Insights into Their Origin,Expansion and Divergence

Plant evolution is characterized by frequent genome duplication events. Expansion of habitat resulted in the origin of many novel genes and genome duplication events which in turn resulted in the expansion of many regulatory gene families. The plant-specific FCS-Like Zinc finger (FLZ) gene family is characterized by the presence of a FCS-Like Zinc finger (FLZ) domain which mediates the protein-protein interaction. In this study, we identified that the expansion of FLZ gene family size in different species is correlated with ancestral and lineage-specific whole genome duplication events. The subsequent gene loss found to have a greater role in determining the size of this gene family in many species. However, genomic block duplications played the significant role in the expansion of FLZ gene family in some species. Comparison of Arabidopsis thaliana and Oryza sativa FLZ gene family revealed monocot and dicot specific evolutionary trends. The FLZ genes were found to be under high purifying selection. The spatiotemporal expression analyses of Arabidopsis thaliana FLZ gene family revealed that majority of the members are highly expressed in reproductive organs. FLZ genes were also found to be highly expressed during vegetative-to-reproductive phase transition which is correlated with the proposed role of this gene family in sugar signaling. The comparison of sequence, structural and expression features of duplicated genes identified lineage-specific redundancy and divergence. This extensive evolutionary analysis and expression analysis of Arabidopsis thaliana FLZ genes will pave the way for further functional analysis of FLZ genes.     

Genome-wide analysis and identification of stress-responsive genes of the CCCH zinc finger family in Solanum lycopersicum

Zinc finger genes comprise a large and diverse gene family. Based on their individual finger structures and spacing, zinc finger proteins are further divided into different families according to their specific molecular functions. Genes in the CCCH family encode zinc finger proteins containing a motif with three cysteines and one histidine. They play important roles in plant growth and development, and in response to biotic and abiotic stresses. However, the limited analysis of the genome sequence has meant that there is no detailed information concerning the CCCH zinc finger family in tomato (Solanum lycopersicum). Here, we identified 80 CCCH zinc finger protein genes in the tomato genome. A complete overview of this gene family in tomato was presented, including the chromosome locations, gene duplications, phylogeny, gene structures and protein motifs. Promoter sequences and expression profiles of putative stress-responsive members were also investigated. These results revealed that, with the exception of four genes, the 80 CCCH genes are distributed over all 12 chromosomes with different densities, and include six segmental duplication events. The CCCH family in tomato could be divided into 12 groups based on their different CCCH motifs and into eight subfamilies by phylogenetic analysis. Analysis showed that almost all CCCH genes contain putative stress-responsive cis-elements in their promoter regions. Nine CCCH genes chosen for further quantitative real-time PCR analysis showed differential expression patterns in three representative tomato tissues. In addition, their expression levels indicated that these genes are mostly involved in the response to mannitol, heat, salicylic acid, ethylene or methyl jasmonate treatments. To the best of our knowledge, this is the first report of a genome-wide analysis of the tomato CCCH zinc finger family. Our data provided valuable information on tomato CCCH proteins and form a foundation for future studies of these proteins, especially for those members that may play important roles in stress responses.     

新基因水稻 OsLSD1 的克隆及拟南芥和水稻类 LSD1 基因家族的生物信息学分析

类 LSD1 (LSD1-like) 基因家族是一类特殊的 C2C2 型锌指蛋白基因,编码植物特有的转录因子 . 目前已经研究的 2 个成员拟南芥 LSD1 (lesions stimulating disease resistance 1) 和 LOL1 (LSD-One-Like 1) 基因均参与植物细胞程序化死亡 (programmed cell death, PCD) 的调控 . 从水稻 cDNA 文库中克隆到 1 个类 LSD1 基因,命名为 OsLSD1. 该基因长 988 bp ,包含一个 432 bp 的开放阅读框,推导的氨基酸序列 (143 个氨基酸 ) 含有 3 个内部保守的锌指结构域 . DNA 印迹结果表明 OsLSD1 基因在水稻基因组中为单拷贝,且在根、茎和叶中表达 . 借助于生物信息学分析技术,从拟南芥和水稻数据库中各识别出 5 个和 7 个 ( 包括 OsLSD1) 类 LSD1 基因 . 分析了这些类 LSD1 基因的结构,蛋白质结构域组成 . 系统进化分析表明,无论基于编码区的核苷酸或氨基酸序列都可以将这些类 LSD1 基因分为 2 类 . 虽然不存在拟南芥或水稻特有的类 LSD1 蛋白,但有些结构域是水稻所特有的,也有些基因是来源于复制事件 .     

Genome-wide characterization of the pectate lyase-like (PLL) genes in Brassica rapa

Pectate lyases (PL) depolymerize demethylated pectin (pectate, EC 4.2.2.2) by catalyzing the eliminative cleavage of α-1,4-glycosidic linked galacturonan. Pectate Lyase-like (PLL) genes are one of the largest and most complex families in plants. However, studies on the phylogeny, gene structure, and expression of PLL genes are limited. To understand the potential functions of PLL genes in plants, we characterized their intron–exon structure, phylogenetic relationships, and protein structures, and measured their expression patterns in various tissues, specifically the reproductive tissues in Brassica rapa. Sequence alignments revealed two characteristic motifs in PLL genes. The chromosome location analysis indicated that 18 of the 46 PLL genes were located in the least fractionated sub-genome (LF) of B. rapa, while 16 were located in the medium fractionated sub-genome (MF1) and 12 in the more fractionated sub-genome (MF2). Quantitative RT-PCR analysis showed that BrPLL genes were expressed in various tissues, with most of them being expressed in flowers. Detailed qRT-PCR analysis identified 11 pollen specific PLL genes and several other genes with unique spatial expression patterns. In addition, some duplicated genes showed similar expression patterns. The phylogenetic analysis identified three PLL gene subfamilies in plants, among which subfamily II might have evolved from gene neofunctionalization or subfunctionalization. Therefore, this study opens the possibility for exploring the roles of PLL genes during plant development.     

Identification of the mulberry genes involved in ethylene biosynthesis and signaling pathways and the expression of MaERF-B2-1 and MaERF-B2-2 in the response to flooding stress

The phytohormone ethylene is essential to plant growth and development. It plays crucial roles in responses to biotic and abiotic stress. The mulberry tree is an important crop plant in countries in which people rear silkworms for silk production. The availability of the mulberry genome has made it possible to identify mulberry genes involved in ethylene biosynthesis and signal pathways. A total of 145 mulberry genes were identified by both homology-based and hidden Markov model (HMM) search, including 29 genes associated with ethylene biosynthesis and 116 genes in the AP2/ERF family. Studies on gene structure have provided a genetic basis for understanding the functions of these genes. The differences in gene expression were also observed in different tissues. The expression of two mulberry genes in the AP2/ERF family, MaERF-B2-1 and MaERF-B2-2, was found to be associated with the response to flooding stress.     

Fine mapping and epistatic interactions of the vernalization gene VRN-D4 in hexaploid wheat

Wheat vernalization requirement is mainly controlled by the VRN1, VRN2, VRN3, and VRN4 genes. The first three have been cloned and have homoeologs in all three genomes. VRN4 has been found only in the D genome (VRN-D4) and has not been cloned. We constructed a high-density genetic map of the VRN-D4 region and mapped VRN-D4 within a 0.09 cM interval in the centromeric region of chromosome 5D. Using telocentric 5D chromosomes generated from the VRN-D4 donor Triple Dirk F, we determined that VRN-D4 is located on the short arm. The VRN-D4 candidate region is colinear with a 2.24 Mb region on Brachypodium distachyon chromosome 4, which includes 127 predicted genes. Ten of these genes have predicted roles in development but we detected no functional polymorphisms associated to VRN-D4. Two recombination events separated VRN-D4 from TaVIL-D1, the wheat homolog of Arabidopsis vernalization gene VIL1, confirming that this gene is not a candidate for VRN-D4. We detected significant interactions between VRN-D4 and other four genes controlling vernalization requirement (Vrn-A1, Vrn-B1, Vrn-D1, and Vrn-B3), which confirmed that VRN-D4 is part of the vernalization pathway and that it is either upstream or is part of the regulatory feedback loop involving VRN1, VRN2 and VRN3 genes. The precise mapping of VRN-D4 and the characterization of its interactions with other vernalization genes provide valuable information for the utilization of VRN-D4 in wheat improvement and for our current efforts to clone this vernalization gene.     

Genome-Wide Identification,Classification, and Expression Analysis of 14-3-3 Gene Family in Populus

Background

In plants, 14-3-3 proteins are encoded by a large multigene family and are involved in signaling pathways to regulate plant development and protection from stress. Although twelve Populus 14-3-3s were identified based on the Populus trichocarpa genome V1.1 in a previous study, no systematic analysis including genome organization, gene structure, duplication relationship, evolutionary analysis and expression compendium has been conducted in Populus based on the latest P. trichocarpa genome V3.0.

Principal Findings

Here, a comprehensive analysis of Populus 14-3-3 family is presented. Two new 14-3-3 genes were identified based on the latest P. trichocarpa genome. In P. trichocarpa, fourteen 14-3-3 genes were grouped into ε and non-ε group. Exon-intron organizations of Populus 14-3-3s are highly conserved within the same group. Genomic organization analysis indicated that purifying selection plays a pivotal role in the retention and maintenance of Populus 14-3-3 family. Protein conformational analysis indicated that Populus 14-3-3 consists of a bundle of nine α-helices (α1-α9); the first four are essential for formation of the dimer, while α3, α5, α7, and α9 form a conserved peptide-binding groove. In addition, α1, α3, α5, α7, and α9 were evolving at a lower rate, while α2, α4, and α6 were evolving at a relatively faster rate. Microarray analyses showed that most Populus 14-3-3s are differentially expressed across tissues and upon exposure to various stresses.

Conclusions

The gene structures and their coding protein structures of Populus 14-3-3s are highly conserved among group members, suggesting that members of the same group might also have conserved functions. Microarray and qRT-PCR analyses showed that most Populus 14-3-3s were differentially expressed in various tissues and were induced by various stresses. Our investigation provided a better understanding of the complexity of the 14-3-3 gene family in poplars.     

Whole genome identification and analysis of FK506-binding protein family genes in grapevine (Vitis vinifera L.)

In plant and animal species FK506-binding protein (FKBP) family genes are important conserved genes and it is defined as the receptors of FK506 and rapamycin, where they work as PPIase and protein folding chaperones. FKBP have been isolated from Arabidopsis thaliana, Oryza sativa, and Zea mays. In grape, twenty-three genes containing the FK506-binding domain (FKBP_C) were first time identified by HMMER and blast research, they were classified into three groups and 17 out of the 23 genes were located on 11 chromosomes (Chr1, 3, 5, 7, 8, 14, 15, 16, 17, 18, and 19). The predicted gene expression pattern and semi-quantitative RT-PCR results revealed that five VvFKBPs were expressed in all tissues, while seven VvFKBPs were expressed only in some of the tissues, and the remaining VvFKBPs were not expressed in leaf, stem, inflorescences, flowers, and a mixture of fruit tissues (small, medium and big-sized fruits). Most of the VvFKBPs in grapevine ‘Summer Black’ were similar to those predicted one in ‘Pinot Noir’ except for VvFKBP16-4 and VvFKBPa. VvFKBP12, FaFKBP12 and PpFKBP12 were cloned from ‘Summer Black’, ‘Sweet Charlie’ and ‘Xiahui 6’. Protein structure analysis confirmed that homologous genes have some differences during the process of protein structure construction. In this study, we characterized and verified 23 FKBP family genes in grapevine (Vitis vinifera L.) as well as their sub-cellular and chromosome location. The successful cloning of CDS regions and protein structural analysis of VvFKBP12, FaFKBP12, and PpFKBP12 can provide useful information for further study.     

Genomics and expression analysis of DHHC-cysteine-rich domain S-acyl transferase protein family in apple

S-acylation is one of a group of lipid modifications that occurs on eukaryotic proteins, mediated by DHHC-CRD-containing proteins, which plays an important role in regulating the membrane association, trafficking and function of target proteins. Although genome-wide identification of PAT family has been carried out in yeast, mice, humans and Arabidopsis, little is known about apple PAT genes. In this study, a total of 33 putative apple PAT proteins, containing DHHC-CRD by domain analysis, have been identified, and were classified into three groups according to the phylogenetic analysis of PAT proteins in apple and Arabidopsis. More complex TMDs in the most MdPATs revealed the PM location of the gene family. Gene structure, gene chromosomal location and paralogs analysis of MdPAT genes within the apple genome demonstrated that tandem and segmental duplications, as well as whole genome duplications, have likely contributed to the expansion and evolution of the PAT gene family in apple. According to the microarray and expressed sequence tag (ESTs) analysis, the different expression patterns indicate that they may play different roles during fruit development and rootstock-scion interactions process. Moreover, PATs were performed expression profile analyses in different tissues, indicating that the PATs are involved in various aspects of physiological and developmental processes of apple. To our knowledge, this is the first report of a genome-wide analysis of the apple PAT gene family, and this genomic analysis of apple DHHC-CRD PAT genes provides the first step towards a functional study of this gene family in apple.     

Genome-wide identification,characterization and expression analysis of the TLP gene family in melon (Cucumis melo L.)

Thaumatin-like proteins (TLPs), which belong to pathogenesis-related (PR) protein family 5 (PR5), are involved in plant host defense and various developmental processes. The functions of the TLP family have been extensively discussed in multiple organisms, whereas the detailed information of this family in melon has not been reported yet. In this study, we identified 28 TLP genes in the melon genome and a N-terminal signal peptide was found highly conserved within each member of this family. Phylogeny analysis indicated that TLPs from melon and other plant species were clustered into ten groups. Twelve segmental and seven tandem duplication gene pairs that underwent purifying selection were identified. TLP genes expressed differentially in different tissues/organs, and were significantly induced after Podosphaera xanthii infection. TLPs in breeding line MR-1 tend to express early after pathogen infection compared with cultivar Top Mark. Our study provides a comprehensive understanding of the melon TLP family and demonstrates their potential roles in disease resistance, therefore provides more reference for further research.     

Fine mapping of Co-x,an anthracnose resistance gene to a highly virulent strain of Colletotrichum lindemuthianum in common bean

Key message

The Co - x anthracnose R gene of common bean was fine-mapped into a 58 kb region at one end of chromosome 1, where no canonical NB-LRR-encoding genes are present in G19833 genome sequence.

Abstract

Anthracnose, caused by the phytopathogenic fungus Colletotrichum lindemuthianum, is one of the most damaging diseases of common bean, Phaseolus vulgaris. Various resistance (R) genes, named Co-, conferring race-specific resistance to different strains of C. lindemuthianum have been identified. The Andean cultivar JaloEEP558 was reported to carry Co-x on chromosome 1, conferring resistance to the highly virulent strain 100. To fine map Co-x, 181 recombinant inbred lines derived from the cross between JaloEEP558 and BAT93 were genotyped with polymerase chain reaction (PCR)-based markers developed using the genome sequence of the Andean genotype G19833. Analysis of RILs carrying key recombination events positioned Co-x at one end of chromosome 1 to a 58 kb region of the G19833 genome sequence. Annotation of this target region revealed eight genes: three phosphoinositide-specific phospholipases C (PI-PLC), one zinc finger protein and four kinases, suggesting that Co-x is not a classical nucleotide-binding leucine-rich encoding gene. In addition, we identified and characterized the seven members of common bean PI-PLC gene family distributed into two clusters located at the ends of chromosomes 1 and 8. Co-x is not a member of Co-1 allelic series since these two genes are separated by at least 190 kb. Comparative analysis between soybean and common bean revealed that the Co-x syntenic region, located at one end of Glycine max chromosome 18, carries Rhg1, a major QTL contributing to soybean cyst nematode resistance. The PCR-based markers generated in this study should be useful in marker-assisted selection for pyramiding Co-x with other R genes.