Genome-wide expression analysis of plant cell cycle modulated genes

Genome-wide expression analysis is rapidly becoming an essential tool for identifying and analysing genes involved in, or controlling, various biological processes ranging from development to responses to environmental cues. The control of cell division involves the temporal expression of different sets of genes, allowing the dividing cell to progress through the different phases of the cell cycle. A landmark study using DNA microarrays to follow the patterns of gene expression in synchronously dividing yeast cells has allowed the identification of several hundreds of genes that are involved in the cell cycle. Although DNA microarrays provide a convenient tool for genome-wide expression analysis, their use is limited to organisms for which the complete genome sequence or a large cDNA collection is available. For other organisms, including most plant species, DNA fragment analysis based methods, such as cDNA-AFLP, provide a more appropriate tool for genome-wide expression analysis. Furthermore, cDNA-AFLP exhibits properties that complement DNA microarrays and, hence, constitutes a useful tool for gene discovery.     

Genome-wide analysis of immunophilin FKBP genes and expression patterns in Zea mays

The receptors for the immunosuppression drugs FK506 and rapamycin are called FKBPs (FK506-binding proteins). FKBPs comprise a large family; they are found in many species, including bacteria, fungi, animals, and plants. As a class of peptidyl-prolyl cis-trans isomerase enzymes, the FKBP genes have been the focus of recent studies on plant stress tolerance and immunology. We identified and analyzed gene families encoding these proteins in maize using computational and molecular biology approaches. Thirty genes were found to encode putative FKBPs according to their FK506-binding domain. The FKBP genes can be classified into single domain and multiple domain members based on the number of the domains. By analysis of the physical locations, the 30 FKBP genes were found to be widely distributed on 10 chromosomes. After analysis of the FKBP phylogenetic tree in the maize genome, we found that the 30 genes revealed two major clades. Gene duplication played a major role in the evolution of FKBP genes, which suggests that the FKBP genes in maize have a pattern significantly different from that of these genes in rice. Based on semi-quantitative RT-PCR, we found that the 30 FKBPs were expressed differently in various tissues in maize, which suggests that FKBP genes play different roles in each tissue. Several FKBPs were expressed at higher levels in roots, indicating that these genes in maize may have similar or overlapping functions.     

Genome-wide aberrant DNA methylation of microRNA host genes in hepatocellular carcinoma

Mature microRNAs (miRNAs) are a class of small non-coding RNAs involved in posttranslational gene silencing. Previous studies found that downregulation of miRNAs is a common feature observed in solid tumors, including hepatocellular carcinoma (HCC). We employed a genome-wide approach to test the hypothesis that DNA methylation alterations in miRNA host genes may cause deregulated miRNA expression in HCC. We analyzed tumor and adjacent non-tumor tissues from 62 Taiwanese HCC cases using Infinium HumanMethylation27 DNA Analysis BeadChips that include 254 CpG sites covering 110 miRNAs from 64 host genes. Expression levels of three identified miRNAs (miR-10a, miR-10b and miR-196b) were measured in a subset of 37 HCC tumor and non-tumor tissues. After Bonferroni adjustment, a total of 54 CpG sites from 27 host genes significantly differed in DNA methylation levels between tumor and adjacent non-tumor tissues with 53 sites significantly hypermethylated in tumor tissues. Among the 54 significant CpG sites, 15 sites had more than 2-fold tumor/non-tumor changes, 17 sites had differences > 10%, and 10 sites had both features [including 8 significantly hypermethylated CpG sites in the host genes of miR-10a, miR-10b and miR-196b (HOXB4, HOXD4 and HOXA9, respectively)]. Significant downregulation of miR-10a was observed in tumor compared with non-tumor tissues (0.50 vs. 1.73, p = 0.031). The concordance for HOXB4 methylation alteration and dysregulation of miR-10a was 73.5%. No significant change was observed for miR-10b expression. Unexpectedly, miR-196b was significantly upregulated in tumor compared with non-tumor tissues (p = 0.0001). These data suggest that aberrant DNA methylation may lead to dysregulation of miR-10a in HCC tumor tissues.     

Genome-wide aberrant DNA methylation of microRNA host genes in hepatocellular carcinoma

Mature microRNAs (miRNAs) are a class of small non-coding RNAs involved in posttranslational gene silencing. Previous studies found that downregulation of miRNAs is a common feature observed in solid tumors, including hepatocellular carcinoma (HCC). We employed a genome-wide approach to test the hypothesis that DNA methylation alterations in miRNA host genes may cause deregulated miRNA expression in HCC. We analyzed tumor and adjacent non-tumor tissues from 62 Taiwanese HCC cases using Infinium HumanMethylation27 DNA Analysis BeadChips that include 254 CpG sites covering 110 miRNAs from 64 host genes. Expression levels of three identified miRNAs (miR-10a, miR-10b and miR-196b) were measured in a subset of 37 HCC tumor and non-tumor tissues. After Bonferroni adjustment, a total of 54 CpG sites from 27 host genes significantly differed in DNA methylation levels between tumor and adjacent non-tumor tissues with 53 sites significantly hypermethylated in tumor tissues. Among the 54 significant CpG sites, 15 sites had more than 2-fold tumor/non-tumor changes, 17 sites had differences > 10%, and 10 sites had both features [including 8 significantly hypermethylated CpG sites in the host genes of miR-10a, miR-10b and miR-196b (HOXB4, HOXD4 and HOXA9, respectively)]. Significant downregulation of miR-10a was observed in tumor compared with non-tumor tissues (0.50 vs. 1.73, p = 0.031). The concordance for HOXB4 methylation alteration and dysregulation of miR-10a was 73.5%. No significant change was observed for miR-10b expression. Unexpectedly, miR-196b was significantly upregulated in tumor compared with non-tumor tissues (p = 0.0001). These data suggest that aberrant DNA methylation may lead to dysregulation of miR-10a in HCC tumor tissues.     

Comparative genomic analysis reveals evolutionary characteristics and patterns of microRNA clusters in vertebrates

MicroRNAs (miRNAs) are a class of small non-coding RNAs that can play important regulatory roles in many important biological processes. Although clustering patterns of miRNA clusters have been uncovered in animals, the origin and evolution of miRNA clusters in vertebrates are still poorly understood. Here, we performed comparative genomic analyses to construct 51 sets of orthologous miRNA clusters (SOMCs) across seven test vertebrate species, a collection of miRNA clusters from two or more species that are likely to have evolved from a common ancestral miRNA cluster, and used these to systematically examine the evolutionary characteristics and patterns of miRNA clusters in vertebrates. We found that miRNA clusters are continuously generated, and most of them tend to be conserved and maintained in vertebrate genomes, although some adaptive gains and losses of miRNA cluster have occurred during evolution. Furthermore, miRNA clusters appeared relatively early in the evolutionary history might suffer from more complicated adaptive gain-and-loss than those young miRNA clusters. Detailed analysis showed that genomic duplication events of ancestral miRNAs or miRNA clusters are likely to be major driving force and apparently contribute to origin and evolution of miRNA clusters. Comparison of conserved with lineage-specific miRNA clusters revealed that the contribution of duplication events for the formation of miRNA cluster appears to be more important for conserved miRNA clusters than lineage-specific. Our study provides novel insights for further exploring the origins and evolution of miRNA clusters in vertebrates at a genome scale.     

Genome-wide analysis of DNA methylation patterns

Cytosine methylation is the most common covalent modification of DNA in eukaryotes. DNA methylation has an important role in many aspects of biology, including development and disease. Methylation can be detected using bisulfite conversion, methylation-sensitive restriction enzymes, methyl-binding proteins and anti-methylcytosine antibodies. Combining these techniques with DNA microarrays and high-throughput sequencing has made the mapping of DNA methylation feasible on a genome-wide scale. Here we discuss recent developments and future directions for identifying and mapping methylation, in an effort to help colleagues to identify the approaches that best serve their research interests.     

Genome-wide analysis of plant glutaredoxin systems

The recent release of the first tree genome (Populus trichocarpa) has allowed a comparison to be made of the multigenic glutaredoxin (Grx) and glutathione reductase (GR) families of this tree with those of other sequenced organisms and especially of the two other fully sequenced plant species, Arabidopsis thaliana and Oryza sativa. Grxs are small proteins involved in disulphide bridge or protein-glutathione adduct reduction, and they are maintained in a reduced form using glutathione and an NADPH-dependent GR. While the P. trichocarpa and O. sativa genomes are nearly five times larger than that of A. thaliana, they contain approximately 45 000 and 37 500 genes compared with the 25 500 genes of A. thaliana. On the one hand, the GR gene composition varies little between species and the gene structures are relatively conserved. On the other hand, the Grx gene family can be divided into three subgroups and the gene content is larger in P. trichocarpa (36 genes) compared with A. thaliana and O. sativa (31 and 27 genes, respectively). This could be partly explained by the occurrence of more duplication events, and this is especially true for one of the three identified Grx subgroups (subgroup III). The expression of most of these genes was confirmed by analysing expressed sequence tags present in various databases. In addition, the expression of Grx of subgroups I and II was examined by RT-PCR in various poplar organs. A complete classification based essentially on gene structure and sequence identity is proposed.     

Genome-wide analysis of NBS-LRR-encoding genes in Arabidopsis

The Arabidopsis genome contains approximately 200 genes that encode proteins with similarity to the nucleotide binding site and other domains characteristic of plant resistance proteins. Through a reiterative process of sequence analysis and reannotation, we identified 149 NBS-LRR-encoding genes in the Arabidopsis (ecotype Columbia) genomic sequence. Fifty-six of these genes were corrected from earlier annotations. At least 12 are predicted to be pseudogenes. As described previously, two distinct groups of sequences were identified: those that encoded an N-terminal domain with Toll/Interleukin-1 Receptor homology (TIR-NBS-LRR, or TNL), and those that encoded an N-terminal coiled-coil motif (CC-NBS-LRR, or CNL). The encoded proteins are distinct from the 58 predicted adapter proteins in the previously described TIR-X, TIR-NBS, and CC-NBS groups. Classification based on protein domains, intron positions, sequence conservation, and genome distribution defined four subgroups of CNL proteins, eight subgroups of TNL proteins, and a pair of divergent NL proteins that lack a defined N-terminal motif. CNL proteins generally were encoded in single exons, although two subclasses were identified that contained introns in unique positions. TNL proteins were encoded in modular exons, with conserved intron positions separating distinct protein domains. Conserved motifs were identified in the LRRs of both CNL and TNL proteins. In contrast to CNL proteins, TNL proteins contained large and variable C-terminal domains. The extant distribution and diversity of the NBS-LRR sequences has been generated by extensive duplication and ectopic rearrangements that involved segmental duplications as well as microscale events. The observed diversity of these NBS-LRR proteins indicates the variety of recognition molecules available in an individual genotype to detect diverse biotic challenges.     

人源microRNA前体的全基因组预测

microRNA（miRNA）是一类不编码蛋白的调控小分子RNA,在真核生物中发挥着广泛而重要的调控功能.由于miRNA的表达具有时空特异性,因而通过计算方法预测miRNA而后有针对性的实验验证是miRNA发现的一条重要途径.降低假阳性率是miRNA预测方法面临的重要挑战.本研究采用集成学习方法构建预测miRNA前体的分类器SVMbagging,对训练集、测试集和独立测试集的结果表明,本研究的方法性能稳健、假阳性率低,具有很好的泛化能力,尤其是当阈值取0.9时,特异性高达99.90%,敏感性在26%以上,适合于全基因组预测.采用SVMbagging在人全基因组中预测miRNA前体,当取阈值0.9时,得到14933个可能的miRNA前体.通过与高通量小RNA测序数据的比较,发现其中4481个miRNA前体具有完全匹配的小RNA序列,与理论估计的真阳性数值非常接近.最后,对32个可能的miRNA进行实验验证,确定其中2条为真实的miRNA.     

Genome-wide linkage analysis of Arabidopsis genes required for leaf development

In most crop species, primary productivity depends mainly on the leaf. However, the genes that contribute to the making of plant leaves remain largely unknown. With a view to identifying the genes involved in leaf development in Arabidopsis thaliana, we previously isolated EMS-induced mutants with abnormally shaped leaves and demonstrated that they fall into 94 complementation groups. We present here the map positions of 76 of these genes, which have been obtained using a high-throughput genetic mapping method, based on the simultaneous coamplification by PCR of 21 polymorphic microsatellites and the semiautomated fluorescent detection of the products. The map positions and F2 mapping populations obtained in this work will be instrumental in the positional cloning of these genes, which are essential for leaf development.