应用基因芯片分析甘蓝型油菜柱头特异表达基因

以甘蓝型油菜(Brassica napus)野生型(宁油10号)及其柱头授粉功能缺失突变体FS-M1为材料,使用油菜基因表达谱芯片筛选甘蓝型油菜柱头特异表达基因。在含有16 540个基因的油菜基因表达谱芯片中(43 803探针),获得了4 410条差异表达探针,选择部分差异表达基因进行实时定量PCR,所得结果与芯片检测结果相吻合。其中,野生型较FS-M1显著上调且获得209个功能注释的探针,对应198个基因,这些特异表达的基因主要富集在水解酶、转移酶、氧化还原酶和转录因子中;涉及较大的基因家族包括:细胞色素P450基因、GDSL脂肪酶/水解酶基因、ABC转运蛋白基因、myb转录因子基因、bHLH转录因子基因、过氧化物酶家族和受体激酶基因等。推测这些基因与甘蓝型油菜柱头发育及授粉功能有关。     

NADP-Dependent enzymes. II: Evolution of the mono- and dinucleotide binding domains

Nicotinamide adenine dinucleotides [NAD and NADP with both referred to as NAD(P)] are among the more diffuse redox cofactors. Despite their stereochemical similarity where the only difference is a phosphomonoester on the ribose near the adenine of NADP, they show different biochemical reactivities with NAD behaving as an oxidant and NADP as a reductant. NAD(P)-dependent enzymes generally share a common open α/β fold with few exceptions only recently structurally characterized. This study of the molecular evolution of the NAD(P) binding domains, possible given the large number of known molecular structures, addresses two main questions: 1) can a common fold exist in different biological systems (divergent evolution) and 2) does a relationship exist among similar biological systems that display different folds (convergent evolution)? Both the structures of mono- and dinucleotide binding domains have been classified by cluster analysis based on the similarity evaluated by their main chain Cα superposition. Moreover, the cofactor conformations and the stereochemical characteristics of their pockets have also been classified by analogous methods on the basis of the published tertiary structures. Two primary results appear: 1) the classification of the mononucleotide binding domains is different from that of the dinucleotide binding folds and 2) both divergent and convergent evolutionary pathways can be hypothesized, the latter less frequently observed and less pronounced but nevertheless evident. The generally accepted hypothesis that dinucleotide binding domains have evolved by gene duplication of primordial genes coding for the smaller mononucleotide binding domains is acceptable but the two halves of the resulting dinucleotide binding domains are evolutionarly uncorrelated. The NH₂-terminal mononucleotide binding domain is less variable than the COOH-terminal half, probably because it involves the binding of the ADP moiety of NAD(P) invariant in all examined systems. There is evidence to postulate that evolutionary pathways for NAD(P)-dependent enzymes are both divergent and convergent. In fact, nearly all combinations of similarity/dissimilarity in overall fold, cofactor conformation, and cofactor binding pocket structural characteristics for each enzyme pair examined are possible. The NAD(P)-dependent enzymes apparently provide a canonical example of an evolutionary principle that “anything goes.” © 1997 Wiley-Liss Inc.     

First trimester trophoblasts forming endothelial-like tubes in vitro emulate a ‘blood vessel development’ gene expression profile

Extravillous cytotrophoblasts isolated from first trimester placenta, and immortalised cell lines derived from them, have the intrinsic ability to form endothelial-like tubes when cultured on Matrigel™ extracellular matrix. This in vitro tube formation may model placental angiogenesis and/or endovascular differentiation by trophoblasts. To interpret the relevance of this phenomenon to placental development, we used a gene expression microarray approach to identify which genes and pathways are associated with the tube-forming phenotype of HTR8/SVneo first trimester trophoblasts (HTR8-M), compared with HTR8/SVneo not forming tubes on plastic culture surface (HTR8-P). Furthermore, we used weighted gene co-expression network analysis (WGCNA) of microarray data to identify modules of co-expressed genes underlying the biological processes. There were 481 genes differentially expressed between HTR8-M and HTR8-P and these were significantly enriched for blood vessel development and related gene ontologies. WGCNA clustered the genes into 9 co-expression modules. One module was significantly associated with HTR8-M (p = 1.15E-05) and contained genes involved in actin cytoskeleton organization, cell migration and blood vessel development, consistent with tube formation on Matrigel. Another module was significantly associated with HTR8-P (p = 1.94E-05) and was enriched for genes involved in mitosis, consistent with proliferation by cells on plastic which do not differentiate. Up-regulation of angiogenesis and vascular development pathways in endovascular trophoblasts in vivo could underpin spiral artery remodelling processes, which are defective in preeclamptic pregnancies.     

Arabidopsis thaliana nicotinate/nicotinamide mononucleotide adenyltransferase (AtNMNAT) is required for pollen tube growth

While mammals and fungi possess nicotinate/nicotinamide mononucleotide adenyltransferase (NMNAT) isoforms, Arabidopsis thaliana only contains a single NMNAT gene, AtNMNAT (At5g55810). We analyzed the enzymatic activity of the AtNMNAT-encoded protein to determine the role of AtNMNAT in plant development. AtNMNAT catalyzed the synthesis of nicotinate adenine dinucleotide (NaAD) from nicotinate mononucleotide (NaMN) in the Preiss-Handler-dependent pathway, and of nicotinamide adenine dinucleotide (NAD) from nicotiamide mononucleotide (NMN) in the Preiss-Handler-independent pathway. Prominent AtNMNAT expression was detected in the male gametophyte. Moreover, AtNMNAT expression was spatio-temporally regulated during microspore development and pollen tube growth. Disruption of the AtNMNAT gene (atnmnat mutant) was characterized by a decrease in NAD content in pollen. Cytological examinations revealed that the atnmnat mutant was gametophytically impaired in in vivo and in vitro pollen tube growth. Our results suggest that metabolic fulfillment via the NAD pathway is indispensable for normal pollen growth and subsequent normal seed production.     

Spearman Correlation Identifies Statistically Significant Gene Expression Clusters in Spinal Cord Development and Injury

An important problem in the analysis of large-scale gene expression data is the validation of gene expression clusters. By examining the temporal expression patterns of 74 genes expressed in rat spinal cord under three different experimental conditions, we have found evidence that some genes cluster together under multiple conditions. Using RT-PCR data from spinal cord development and two sets of microarray data from spinal injury, we applied Spearman correlation to identify clusters and to assign P values to pairs of genes with highly similar temporal expression patterns. We found that 15% of genes occurred in statistically significant pairs in all three experimental conditions, providing both statistical and experimental support for the idea that genes that cluster together are co-regulated. In addition, we demonstrated that DNA microarray and RT-PCR data are comparable, and can be combined to confirm gene expression relationships.     

Molecular dissection of the response of a rice leucine-rich repeat receptor-like kinase (LRR-RLK) gene to abiotic stresses

Leucine-rich repeat (LRR) receptor-like kinase (RLK) proteins play key roles in a variety of biological pathways. In a previous study, we analyzed the members of the rice LRR-RLK gene family using in silico analysis. A total of 23 LRR-RLK genes were selected based on the expression patterns of a genome-wide dataset of microarrays. The Oryza sativa gamma-ray induced LRR-RLK1 (OsGIRL1) gene was highly induced by gamma irradiation. Therefore, we studied its expression pattern in response to various different abiotic and phytohormone treatments. OsGIRL1 was induced on exposure to abiotic stresses such as salt, osmotic, and heat, salicylic acid (SA), and abscisic acid (ABA), but exhibited downregulation in response to jasmonic acid (JA) treatment. The OsGIRL1 protein was clearly localized at the plasma membrane. The truncated proteins harboring juxtamembrane and kinase domains (or only harboring a kinase domain) exhibited strong autophosphorylation. The biological function of OsGIRL1 was investigated via heterologous overexpression of this gene in Arabidopsis plants subjected to gamma-ray irradiation, salt stress, osmotic stress, and heat stress. A hypersensitive response was observed in response to salt stress and heat stress, whereas a hyposensitive response was observed in response to gamma-ray treatment and osmotic stress. These results provide critical insights into the molecular functions of the rice LRR-RLK genes as receptors of external signals.     

EDISA: extracting biclusters from multiple time-series of gene expression profiles

Background  

Cells dynamically adapt their gene expression patterns in response to various stimuli. This response is orchestrated into a number of gene expression modules consisting of co-regulated genes. A growing pool of publicly available microarray datasets allows the identification of modules by monitoring expression changes over time. These time-series datasets can be searched for gene expression modules by one of the many clustering methods published to date. For an integrative analysis, several time-series datasets can be joined into a three-dimensional gene-condition-time dataset, to which standard clustering or biclustering methods are, however, not applicable. We thus devise a probabilistic clustering algorithm for gene-condition-time datasets.     

Highly interconnected genes in disease-specific networks are enriched for disease-associated polymorphisms

Background

Complex diseases are associated with altered interactions between thousands of genes. We developed a novel method to identify and prioritize disease genes, which was generally applicable to complex diseases.

Results

We identified modules of highly interconnected genes in disease-specific networks derived from integrating gene-expression and protein interaction data. We examined if those modules were enriched for disease-associated SNPs, and could be used to find novel genes for functional studies. First, we analyzed publicly available gene expression microarray and genome-wide association study (GWAS) data from 13, highly diverse, complex diseases. In each disease, highly interconnected genes formed modules, which were significantly enriched for genes harboring disease-associated SNPs. To test if such modules could be used to find novel genes for functional studies, we repeated the analyses using our own gene expression microarray and GWAS data from seasonal allergic rhinitis. We identified a novel gene, FGF2, whose relevance was supported by functional studies using combined small interfering RNA-mediated knock-down and gene expression microarrays. The modules in the 13 complex diseases analyzed here tended to overlap and were enriched for pathways related to oncological, metabolic and inflammatory diseases. This suggested that this union of the modules would be associated with a general increase in susceptibility for complex diseases. Indeed, we found that this union was enriched with GWAS genes for 145 other complex diseases.

Conclusions

Modules of highly interconnected complex disease genes were enriched for disease-associated SNPs, and could be used to find novel genes for functional studies.     

LinkNMF: Identification of histone modification modules in the human genome using nonnegative matrix factorization

Histone modifications are ubiquitous processes involved in various cellular mechanisms. Systemic analysis of multiple chromatin modifications has been used to characterize various chromatin states associated with functional DNA elements, gene expression, and specific biological functions. However, identification of modular modification patterns is still required to understand the functional associations between histone modification patterns and specific chromatin/DNA binding factors. To recognize modular modification patterns, we developed a novel algorithm that combines nonnegative matrix factorization (NMF) and a clique-detection algorithm. We applied it, called LinkNMF, to generate a comprehensive modification map in human CD4 + T cell promoter regions. Initially, we identified 11 modules not recognized by conventional approaches. The modules were grouped into two major classes: gene activation and repression. We found that genes targeted by each module were enriched with distinguishable biological functions, suggesting that each modular pattern plays a unique functional role. To explain the formation of modular patterns, we investigated the module-specific binding patterns of chromatin regulators. Application of LinkNMF to histone modification maps of diverse cells and developmental stages will be helpful for understanding how histone modifications regulate gene expression. The algorithm is available on our website at biodb.kaist.ac.kr/LinkNMF.