Identifying Gene Interaction Enrichment for Gene Expression Data

Gene set analysis allows the inclusion of knowledge from established gene sets, such as gene pathways, and potentially improves the power of detecting differentially expressed genes. However, conventional methods of gene set analysis focus on gene marginal effects in a gene set, and ignore gene interactions which may contribute to complex human diseases. In this study, we propose a method of gene interaction enrichment analysis, which incorporates knowledge of predefined gene sets (e.g. gene pathways) to identify enriched gene interaction effects on a phenotype of interest. In our proposed method, we also discuss the reduction of irrelevant genes and the extraction of a core set of gene interactions for an identified gene set, which contribute to the statistical variation of a phenotype of interest. The utility of our method is demonstrated through analyses on two publicly available microarray datasets. The results show that our method can identify gene sets that show strong gene interaction enrichments. The enriched gene interactions identified by our method may provide clues to new gene regulation mechanisms related to the studied phenotypes. In summary, our method offers a powerful tool for researchers to exhaustively examine the large numbers of gene interactions associated with complex human diseases, and can be a useful complement to classical gene set analyses which only considers single genes in a gene set.     

The small MbtH-like protein encoded by an internal gene of the balhimycin biosynthetic gene cluster is not required for glycopeptide production

The balhimycin biosynthetic gene cluster of the glycopeptide producer Amycolatopsis balhimycina includes a gene (orf1) with unknown function. orf1 shows high similarity to the mbtH gene from Mycobacterium tuberculosis. In almost all nonribosomal peptide synthetase (NRPS) biosynthetic gene clusters, we could identify a small mbtH-like gene whose function in peptide biosynthesis is not known. The mbtH-like gene is always colocalized with the NRPS genes; however, it does not have a specific position in the gene cluster. In all glycopeptide biosynthetic gene clusters the orf1-like gene is always located downstream of the gene encoding the last module of the NRPS. We inactivated the orf1 gene in A. balhimycina by generating a deletion mutant. The balhimycin production is not affected in the orf1-deletion mutant and is indistinguishable from that of the wild type. For the first time, we show that the inactivation of an mbtH-like gene does not impair the biosynthesis of a nonribosomal peptide.     

From gene trees to species trees II: species tree inference by minimizing deep coalescence events

When gene copies are sampled from various species, the resulting gene tree might disagree with the containing species tree. The primary causes of gene tree and species tree discord include incomplete lineage sorting, horizontal gene transfer, and gene duplication and loss. Each of these events yields a different parsimony criterion for inferring the (containing) species tree from gene trees. With incomplete lineage sorting, species tree inference is to find the tree minimizing extra gene lineages that had to coexist along species lineages; with gene duplication, it becomes to find the tree minimizing gene duplications and/or losses. In this paper, we present the following results: 1) The deep coalescence cost is equal to the number of gene losses minus two times the gene duplication cost in the reconciliation of a uniquely leaf labeled gene tree and a species tree. The deep coalescence cost can be computed in linear time for any arbitrary gene tree and species tree. 2) The deep coalescence cost is always not less than the gene duplication cost in the reconciliation of an arbitrary gene tree and a species tree. 3) Species tree inference by minimizing deep coalescence events is NP-hard.     

人α-心钠素全基因的化学合成

用固相亚磷酰胺法(solid-phase phosphoramidite method)合成了人α-心钠素(简称α-hANP)的全基因。合成的α-hANP基因全长为96bp,包括编码α-hANP的结构基因、基因5＇端的谷氨酸(作为表达产物用内肽酶Glu-C专一裂解的位点)密码子GAA、基因3＇端的终止信号TGA及基因两端的接头部分。基因分7个寡聚核苷酸片段在DNA合成仪上合成,然后一次性酶促连接成为完整的α—hANP双链基因。化学合成的基因克隆到M13载体上,经分子杂交鉴定筛选出的α-hANP基因克隆株,用双脱氧链终止法进行序列分析,证明核苷酸顺序正确。     

The problem of the gene

During the early 20th century the diverse practices of genetics were unified by the concept of the gene. This classical gene was simultaneously a unit of structure, function, mutation, and recombination. Starting in the 1940s, however, the classical gene began to fragment. Today when we speak of a gene for some malady, a regulatory gene, a structural gene, or a gene frequency, it is entirely possible that we are deploying different gene concepts even though we are using the same term. The problem of the gene addresses the fragmentation of the classical gene concept by asking to what extent a comprehensive and unifying gene concept is possible or desirable. Fully comprehensive gene concepts seem untenable today, but, within different disciplinary domains, unifying, but non-comprehensive, gene concepts can be epistemically worthwhile. The problem of the gene persists, however, not because of its epistemic value, but because of its political value. Using both the arguments for newly proposed gene concepts and the historical dispute over the classical gene, I argue that the desirability of gene concepts rests in part on the political ramifications of their deployment and contestation.     

利用套叠PCR技术进行基因突变和拼接

利用套叠PCR技术（又称重叠区扩增基因拼接法）对hGM-CSF基因内第28位氨基酸处的糖基化位点进行突变和进行人促性腺激素基因,腺苷酸激酶短肽与胰岛素样生长因子－基因三者之间的拼接,结果表明采用该技术能在体外实行有效的基因重组和定点突变,其成功率为100％,这一技术不需要内切酶消化和连接酶处理,技术操作员简单易行,在基因拼接,基因内部突变方面具有良好的应用价值。     

Reconstructing the evolution of genes along the species tree

A model and algorithm are proposed to infer the evolution of a gene family described by the corresponding gene tree, with respect to the species evolution described by the corresponding species tree. The model describes the evolution using the new concept of a nested tree. The algorithm performance is illustrated by the example of several orthologous protein groups. The considered evolutionary events are speciation, gene duplication and loss, and horizontal gene transfer retaining the original gene copy. The transfer event with the loss of the original gene copy is considered as a combination of gene transfer and loss. The model maps each evolutionary event onto the species phylogeny.     

A fusion gene in man: DNA sequence analysis of the abnormal globin gene of Hemoglobin Miyada

An abnormal globin gene from a patient heterozygous for Hemoglobin Miyada was cloned and sequenced. The results indicated that the 5′ flanking region and the 5′ side of the gene were identical to those of a β-globin gene and that the 3′ side was identical to that of a γ-globin gene. The part of the gene identical to a β-globin gene shifted to the part identical to the δ-globin gene somewhere in a homologous sequence region between the third nucleotide of the 17th codon and the second nucleotide of the 22nd codon of these two genes. Thus, results of analysis of the nucleotide sequence support the idea that the abnormal globin gene of Hemoglobin Miyada was generated as a fusion gene by unequal crossing over between a β- and a δ-globin gene.     

Molecular cloning, sequence analysis, and expression of the yeast alcohol acetyltransferase gene.

The ATF1 gene, which encodes alcohol acetyltransferase (AATase), was cloned from Saccharomyces cerevisiae and brewery lager yeast (Saccharomyces uvarum). The nucleotide sequence of the ATF1 gene isolated from S. cerevisiae was determined. The structural gene consists of a 1,575-bp open reading frame that encodes 525 amino acids with a calculated molecular weight of 61,059. Although the yeast AATase is considered a membrane-bound enzyme, the results of a hydrophobicity analysis suggested that this gene product does not have a membrane-spanning region that is significantly hydrophobic. A Southern analysis of the yeast genomes in which the ATF1 gene was used as a probe revealed that S. cerevisiae has one ATF1 gene, while brewery lager yeast has one ATF1 gene and another, homologous gene (Lg-ATF1). Transformants carrying multiple copies of the ATF1 gene or the Lg-ATF1 gene exhibited high AATase activity in static cultures and produced greater concentrations of acetate esters than the control.     

The lemA gene required for pathogenicity of Pseudomonas syringae pv. syringae on bean is a member of a family of two-component regulators.

The lemA gene of the plant pathogen Pseudomonas syringae pv. syringae is required for disease lesion formation on bean plants. Cosmid clones that complemented a lemA mutant in trans were isolated previously. The lemA gene was localized by subcloning and transposon mutagenesis. The lemA region and flanking DNA were sequenced, and an open reading frame of 2.7 kb was identified. The nucleotide and predicted amino acid sequences of the lemA gene showed sequence similarity to a family of prokaryotic two-component regulatory proteins. Unlike most of the previously described two-component systems, the lemA gene product contained homology to both components in one protein. Mutations introduced upstream and downstream of the lemA gene failed to locate a gene for a second protein component but identified the putative cysM gene of P. syringae pv. syringae. The cysM gene was located upstream of the lemA gene and was divergently transcribed. The lemA gene product was expressed at low levels in P. syringae pv. syringae and appeared to be positively auto-regulated.     

ColIb plasmid genes that inhibit the replication of T5 and T7 bacteriophage

The colicin Ib (ColIb) plasmid genes that inhibit the replication of the T5-like and T7 bacteriophage have been cloned on an approximately 7200-bp ClaI fragment and their sites relative to each other and to the colicin immunity (imm) gene have been mapped. The inhibition of wild-type T7 by the clone is shown to be caused by the same gene or genes (pic) that cause the inhibition of T7 kinase-negative mutants and is a different gene than the one that causes inhibition of T5 (ibf or abi). The pic gene does not hybridize to the pif genes of the F plasmid that also cause the replication of T7 to be inhibited. The abi gene and the pic gene map very closely together but are under the control of different promoters. The abi gene has a maximum size of 900 bp and lies approximately 3000 bp away from the immunity gene, distal to the colicin gene. A site which maps in or near the gene binds very tightly to Escherichia coli RNA polymerase. The pic gene or genes lie between the abi gene and the imm gene and are contiguous with abi. Promoters for pic have been mapped and hypotheses to explain the inhibition of T7 by a cloned gene but not the whole ColIb plasmid are presented.     

Organ-specific and dosage-dependent expression of a leaf/stem specific gene from potato after tagging and transfer into potato and tobacco plants.

ST-LS1, a single copy gene from potato displaying a leaf/stem specific gene expression, was tagged by an exon modification and introduced into both potato and tobacco cells using Agrobacterium vectors. After regeneration of whole plants, the expression of the tagged gene was analyzed with respect to its organ specificity and compared to the expression of the corresponding resident gene. The expression of the transferred gene in transgenic plants closely followed the expression of the resident gene. No marked influence of the plant species serving as host was observed. The level of expression of the introduced gene varied by a factor of at least 100 in independent transformants when normalized to the expression of the resident gene. Southern analysis performed on the transformed plants indicated a correlation between copy number of the introduced gene and its expression level. The activity of the tagged gene as well as of the resident gene was significantly inhibited by treatment of the transgenic plants with the herbicide norfluorazon, indicating that this gene activity is dependent on the presence of functional chloroplasts in the leaves.     

Molecular cloning of the Candida maltosa ADE1 gene.

The structural gene (ADE1) encoding phosphoribosyl-aminoimidazole-succinocarboxamide synthetase (SAICAR synthetase; EC 6.3.2.6) in Candida maltosa has been isolated by functional complementation of an ade1 strain of Saccharomyces cerevisiae. The gene was localized on a 2.5-kb BamHI DNA fragment. Nucleotide sequence analysis of the cloned gene has revealed an open reading frame encoding a protein (SAICAR synthetase) with an Mr of 32,751. The codon bias index, 0.68, indicates that the ADE1 gene is a moderately highly expressed gene. The cloned gene shows 63.5% nt identity and 65.2% deduced amino acid identity with the S. cerevisiae ADE1 gene which encodes the same enzymatic activity. The gene may be used as a convenient genetic marker for construction of a new host-vector system for C. maltosa.     

The evolution of the novel Sdic gene cluster in Drosophila melanogaster

The origin of new genes and of new functions for existing genes are fundamental processes in molecular evolution. Sdic is a newly evolved gene that arose recently in the D. melanogaster lineage. The gene encodes a novel sperm motility protein. It is a chimeric gene formed by duplication of two other genes followed by multiple deletions and other sequence rearrangements. The Sdic gene exists in several copies in the X chromosome, and is presumed to have undergone several duplications to form a tandemly arrayed gene cluster. Given the very recent origin of the gene and the gene cluster, the analysis of the composition of this gene cluster represents an excellent opportunity to study the origin and evolution of new gene functions and the fate of gene duplications. We have analyzed the nucleotide sequence of this region and reconstructed the evolutionary history of this gene cluster. We found that the cluster is composed by four tandem copies of Sdic; these duplicates are very similar but can be distinguished by the unique pattern of insertions, deletions, and point mutations in each copy. The oldest gene copy in the array has a 3' exon that has undergone accelerated diversification, and also shows divergent regulatory sequences. Moreover, there is evidence that this might be the only gene copy in the tandem array that is transcribed at a significant level, expressing a novel sperm-specific protein. There is also a retrotransposon located at the 3' end of each Sdic gene copy. We argue that this gene cluster was formed in the last two million years by at least three tandem duplications and one retrotransposition event.