Genome-wide dissection of the chalcone synthase gene family in <Emphasis Type="Italic">Oryza sativa</Emphasis>

Enzymes of the chalcone synthase (CHS) family catalyze the generation of multiple secondary metabolites in fungi, plants, and bacteria. These metabolites have played key roles in antimicrobial activity, UV protection, flower pigmentation, and pollen fertility during the evolutionary process of land plants. We performed a genome-wide investigation about CHS genes in rice (Oryza sativa). The phylogenetic relationships, gene structures, chromosomal locations, and functional predictions of the family members were examined. Twenty-seven CHS family genes (OsCHS01–27) were identified in the rice genome and were found to cluster into six classes according to their phylogenetic relationships. The 27 OsCHS genes were unevenly distributed on six chromosomes, and 17 genes were found in the genome duplication zones with two segmental duplication and five tandem duplication events that may have played key roles in the expansion of the rice CHS gene family. In addition, the OsCHS genes exhibited diverse expression patterns under salicylic acid treatment. Our results revealed that the OsCHS genes exhibit both diversity and conservation in many aspects, which will contribute to further studies of the function of the rice CHS gene family and provide a reference for investigating this family in other plants.     

Genome-wide identification and resistance expression analysis of the <Emphasis Type="Italic">NBS</Emphasis> gene family in <Emphasis Type="Italic">Triticum urartu</Emphasis>

As the largest class of resistant genes, the nucleotide binding site (NBS) has been studied extensively at a genome-wide level in rice, sorghum, maize, barley and hexaploid wheat. However, no such comprehensive analysis has been conducted of the NBS gene family in Triticum urartu, the donor of the A genome to the common wheat. Using a bioinformatics method, 463 NBS genes were isolated from the whole genome of T. urartu, of which 461 had location information. The expansion pattern and evolution of the 461 NBS candidate proteins were analyzed, and 118 of them were duplicated. By calculating the lengths of the copies, it was inferred that the NBS resistance gene family of T. urartu has experienced at least two duplication events. Expression analysis based on RNA-seq data found that 6 genes were differentially expressed among Tu38, Tu138 and Tu158 in response to Blumeria graminis f.sp.tritici (Bgt). Following Bgt infection, the expression levels of these genes were up-regulated. These results provide critical references for further identification and analysis of NBS family genes with important functions.     

Ancestral and more recently acquired syntenic relationships of MADS-box genes uncovered by the <Emphasis Type="Italic">Physcomitrella patens</Emphasis> pseudochromosomal genome assembly

Key message

The Physcomitrella pseudochromosomal genome assembly revealed previously invisible synteny enabling realisation of the full potential of shared synteny as a tool for probing evolution of this plant’s MADS-box gene family.

Abstract

Assembly of the sequenced genome of Physcomitrella patens into 27 mega-scaffolds (pseudochromosomes) has confirmed the major predictions of our earlier model of expansion of the MADS-box gene family in the Physcomitrella lineage. Additionally, microsynteny has been conserved in the immediate vicinity of some recent duplicates of MADS-box genes. However, comparison of non-syntenic MIKC MADS-box genes and neighbouring genes indicates that chromosomal rearrangements and/or sequence degeneration have destroyed shared synteny over longer distances (macrosynteny) around MADS-box genes despite subsets comprising two or three MIKC genes having remained syntenic. In contrast, half of the type I MADS-box genes have been transposed creating new syntenic relations with MIKC genes. This implies that conservation of ancient ancestral synteny of MIKC genes and of more recently acquired synteny of type I and MIKC genes may be selectively advantageous. Our revised model predicts the birth rate of MIKC genes in Physcomitrella is higher than that of type I genes. However, this difference is attributable to an early tandem duplication and an early segmental duplication of MIKC genes prior to the two polyploidisations that account for most of the expansion of the MADS-box gene family in Physcomitrella. Furthermore, this early segmental duplication spawned two chromosomal lineages: one with a MIKC ^C gene, belonging to the PPM2 clade, in close proximity to one or a pair of MIKC* genes and another with a MIKC ^C gene, belonging to the PpMADS-S clade, characterised by greater separation from syntenic MIKC* genes. Our model has evolutionary implications for the Physcomitrella karyotype.

     

The roles of segmental and tandem gene duplication in the evolution of large gene families in <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>

Background

Most genes in Arabidopsis thaliana are members of gene families. How do the members of gene families arise, and how are gene family copy numbers maintained? Some gene families may evolve primarily through tandem duplication and high rates of birth and death in clusters, and others through infrequent polyploidy or large-scale segmental duplications and subsequent losses.

Results

Our approach to understanding the mechanisms of gene family evolution was to construct phylogenies for 50 large gene families in Arabidopsis thaliana, identify large internal segmental duplications in Arabidopsis, map gene duplications onto the segmental duplications, and use this information to identify which nodes in each phylogeny arose due to segmental or tandem duplication. Examples of six gene families exemplifying characteristic modes are described. Distributions of gene family sizes and patterns of duplication by genomic distance are also described in order to characterize patterns of local duplication and copy number for large gene families. Both gene family size and duplication by distance closely follow power-law distributions.

Conclusions

Combining information about genomic segmental duplications, gene family phylogenies, and gene positions provides a method to evaluate contributions of tandem duplication and segmental genome duplication in the generation and maintenance of gene families. These differences appear to correspond meaningfully to differences in functional roles of the members of the gene families.

     

Genome-wide analysis and characterization of molecular evolution of the <Emphasis Type="Italic">HCT</Emphasis> gene family in pear (<Emphasis Type="Italic">Pyrus bretschneideri</Emphasis>)

Lignin is a major component of stone cells in pear fruit, which significantly affects fruit quality. Hydroxycinnamoyl CoA: shikimate hydroxycinnamoyl transferase (HCT), a recently discovered enzyme in plants, is an important gene that participates in the formation of lignin. Although HCT gene cloning and expression patterns have been studied in several species, including pear, there is still no extensive genome-wide bioinformatics analysis on the whole gene family, and the evolutionary history of HCT gene family is still unknown. A total of 82 HCT genes were identified in pear, most of which have one or two exons, and all with the conserved HXXXD motif and transferase domains. Based on the structural characteristics and phylogenetic analysis of these sequences, the HCT gene family genes could be classified into four main groups. Structural analysis also revealed that 25 % of HCT genes share a MYB binding site. Expansion of the HCT gene family mostly occurred before the divergence between Arabidopsis and Rosaceae, with whole-genome duplication or segmental duplication events playing the most important role in the expansion of the HCT gene family in pear. At the same time, purifying selection also played a critical role in the evolution of HCT genes. Five of the 82 HCT genes were verified by qRT-PCR to correspond to the pattern of stone cell formation during pear fruit development. The genome-wide identification, chromosome localization, gene structures, synteny, and expression analyses of pear HCT genes provide an overall insight into HCT gene family and their potential involvement in growth and development of stone cells.     

Genome-wide bioinformatics analysis of Dof transcription factor gene family of chickpea and its comparative phylogenetic assessment with Arabidopsis and rice

The genome mining of chickpea (Cicer arietinum L.) revealed a total of 37 putative Dof genes using NCBI BLAST search against the genome with a highly conserved Dof domain. The translated Dof proteins possessed 150–493 amino acid residues with molecular weight ranging from 16.9 to 54.4 kD and pI varied from 4.98 to 9.64 as revealed by ExPASy server ProtParam. The exon–intron organization showed predominance of intronless Dof genes in chickpea. The predicted Dof genes were distributed among the eight chromosomes with a maximum of 9 Dof genes present on chromosome 7 and a single Dof gene was found on chromosome 8.The predominance of segmental gene duplication as compared to tandem duplication was observed which might be the prime cause of Dof gene family expansion in chickpea. The cis-regulatory element analysis revealed the presence of light-responsive, hormone-responsive, endosperm-specific, meristem-specific and stress-responsive elements. Comprehensive phylogenetic analyses of Dof genes of chickpea with Arabidopsis, rice, soybean and pigeonpea revealed several orthologs and paralogs assisting in understanding the putative functions of CaDof genes. The functional divergence and site-specific selective pressures of chickpea Dof genes have been investigated. The bioinformatics-based genome-wide assessment of Dof gene family of chickpea attempted in the present study could be a significant step for deciphering novel Dof genes based on genome-wide expression profiling.     

Tandem duplications of receptor-like kinase genes reshaped the homologous chromosome 1 regions in <Emphasis Type="Italic">Oryza sativa</Emphasis> ssp. <Emphasis Type="Italic">indica</Emphasis> and <Emphasis Type="Italic">japonica</Emphasis> and their possible ancestral genomes

Plant receptor-like kinase (Rlk) genes form a large family, each encoding a protein with a signal motif, a single transmembrane region, and a cytoplasmic kinase domain. Various gene duplications have contributed to the establishment and expansion of the family. Here, we characterized the formation and evolution of the Rlk gene family in cultivated rice and their possible progenitors. Using wheat Rlk gene sequences, we identified orthologs from the genomes of domesticated rice subspecies Oryza sativa ssp. japonica and ssp. indica and their putative progenitors O. glaberrima and O. rufipogon. The four chromosome 1 orthologous regions ranged from 103 to 281 kb comprising 181 syntenic blocks with 75 to 100% sequence identity. These regions contained 11–19 Triticum aestivum kinases (Taks) and 10–15 Lr10 receptor-like kinases (Lrks) organized in clusters and 3–12 transposable elements (TEs). Dot plot analyses showed that the 4 regions had 21–37 conserved catalytic domains, mainly in protein kinases (PKs) and tyrosine kinases (TyrKs) in coupling state. Over 50% of the sequences of glaberrima/rufipogon and japonica/indica pairs were colinear, while japonica/indica displayed a marked sequence expansion with duplicated genes and TEs. A total of 2312 single nucleotide polymorphisms (SNPs) and insertion-deletions (INDELs) were identified between japonica and indica. Duplication of the Rlk genes in O. glaberrima and O. rufipogon occurred after the grass species radiation and before the divergence of O. rufipogon from O. glaberrima; the orthologous Rlk genes from O. japonica and O. indica duplicated after O. sativa separated from O. rufipogon; paralogs, obtained through extensive duplication, happened after the separation of rice from maize. Tandem duplication was the major factor contributing to the gene copy number variation and genome size expansion.     

Low molecular weight glutenin subunit gene composition at <Emphasis Type="Italic">Glu</Emphasis>-<Emphasis Type="Italic">D3</Emphasis> loci of <Emphasis Type="Italic">Aegilops tauschii</Emphasis> and common wheat and a further view of wheat evolution

Key message

A comprehensive comparison of LMW-GS genes between Ae. tauschii and its progeny common wheat.

Abstract

Low molecular weight glutenin subunits (LMW-GSs) are determinant of wheat flour processing quality. However, the LMW-GS gene composition in Aegilops tauschii, the wheat D genome progenitor, has not been comprehensively elucidated and the impact of allohexaploidization on the Glu-D3 locus remains elusive. In this work, using the LMW-GS gene molecular marker system and the full-length gene-cloning method, LMW-GS genes at the Glu-D3 loci of 218 Ae. tauschii and 173 common wheat (Triticum aestivum L.) were characterized. Each Ae. tauschii contained 11 LMW-GS genes, and the whole collection was divided into 25 haplotypes (AeH01–AeH25). The Glu-D3 locus in common wheat lacked the LMW-GS genes D3-417, D3-507 and D3-552, but shared eight genes of identical open reading frame (ORF) sequences when compared to that of Ae. tauschii. Therefore, the allohexaploidization induces deletions, but exerts no influence on LMW-GS gene coding sequences at the Glu-D3 locus. 92.17% Ae. tauschii had 7-9 LMW-GSs, more than the six subunits in common wheat. The haplotypes AeH16, AeH20 and AeH23 of Ae. tauschii ssp. strangulate distributed in southeastern Caspian Iran were the main putative D genome donor of common wheat. These results facilitate the utilization of the Ae. tauschii glutenin gene resources and the understanding of wheat evolution.