Detecting correlation between sequence and expression divergences in a comparative analysis of human serpin genes

Physiological functions and characteristic structures of the serpin gene superfamily have been studied extensively, yet the evolution of the serpin genes remains unclear. Gene duplication in this superfamily may shed light on this issue. Two models are used to predict the preservation of duplicated genes: the classical model and the duplication–degeneration–complementation (DDC) model. In this study, we analyzed the phylogenetic relationships of 33 human serpin genes and the expression data of some members of the serpin superfamily from a DNA microarray of human leukemia U937 cells with stably inducible expression of the leukemia-related AML1-ETO gene. We then determined the utility of the DDC model by mapping serpin superfamily expression data to the phylogenetic tree. The correlation between sequence and expression divergences as measured by the Pearson correlation coefficient indicated that human serpin genes evolved under the DDC model. Our study provides a new strategy for comparative analysis of gene sequences and microarray data.     

Detecting correlation between characters in a comparative analysis with uncertain phylogeny

Abstract.— The importance of accommodating the phylogenetic history of a group when performing a comparative analysis is now widely recognized. The typical approaches either assume the tree is known without error, or they base inferences on a collection of well-supported trees or on a collection of trees generated under a stochastic model of cladogenesis. However, these approaches do not adequately account for the uncertainty of phylogenetic trees in a comparative analysis, especially when data relevant to the phylogeny of a group are available. Here, we develop a method for performing comparative analyses that is based on an extension of Felsenstein's independent contrasts method. Uncertainties in the phylogeny, branch lengths, and other parameters are accommodated by averaging over all possible trees, weighting each by the probability that the tree is correct. We do this in a Bayesian framework and use Markov chain Monte Carlo to perform the high-dimensional summations and integrations required by the analysis. We illustrate the method using comparative characters sampled from Anolis lizards.     

A comparative analysis of serpin genes in the silkworm genome

Serine protease inhibitors (serpins) are a superfamily of proteins, most of which control protease-mediated processes by inhibiting their cognate enzymes. Sequencing of the silkworm genome provides an opportunity to investigate serpin structure, function, and evolution at the genome level. There are thirty-four serpin genes in Bombyx mori. Six are highly similar to their Manduca sexta orthologs that regulate innate immunity. Three alternative exons in serpin1 gene and four in serpin28 encode a variable region including the reactive site loop. Splicing of serpin2 pre-mRNA yields variations in serpin2A, 2A′ and 2B. Sequence similarity and intron positions reveal the evolutionary pathway of seven serpin genes in group C. RT-PCR indicates an increase in the mRNA levels of serpin1, 3, 5, 6, 9, 12, 13, 25, 27, 32 and 34 in fat body and hemocytes of larvae injected with bacteria. These results suggest that the silkworm serpins play regulatory roles in defense responses.     

Variation in synonymous substitution rates among mammalian genes and the correlation between synonymous and nonsynonymous divergences

Using mammalian gene sequences, the variances in the numbers of synonymous and nonsynonymous substitutions among genes were estimated together with the correlation coefficient between the two. The expected correlation coefficient can be obtained under the neutral theory using these estimated values of the variances. The expected coefficient is found to often be one-half to two-thirds of the observed value. Possible causes for the disagreement were discussed, such as correlated selective constraints on the two types of substitutions and excess doublet mutations. The variance of mutation rate and that of selective constraint were also estimated. The results show that the coefficient of variation of the former is 0.2–0.3, whereas that of the latter is 0.7–0.9. Correspondence to: T. Ohta     

Conservation of a pair of serpin 2 genes and their expression in Amphiesmenoptera

Silk secreted by the larvae of Hydropsyche angustipennis (Trichoptera) contains serpins HaSerp2A and HaSerp2B that are homologous to serpin 2 known from several lepidopterans and some other insects. The gene HaSerp2A is 2684 bp downstream from the HaSerp2B gene. The genes possess identical exon/intron segmentation (9 exons) and their sequences are nearly identical: only 8 out of 1203 nt differ in the coding region, 4 out of 567 nt in the introns and 2 out of 52 nt in 3' UTR. Both genes are highly expressed in the silk glands whereas expression in larval carcass devoid of the silk glands is hard to detect. Translation products of the genes consist of 401 amino acids, are 98.8% identical, and are secreted as 45 kDa proteins into silk. Homologous genes in similar tandem arrangement occur on chromosome 15 of Bombyx mori (Lepidoptera). The upstream gene BmSerp2B is modified in several exons and does not seem to produce functional mRNA. The gene BmSerp2A contains two copies of exon 9, of which only the second one is used. One kind of mRNA does and the other does not include exon 1, which encodes a signal peptide. The mRNA yielding secreted BmSerp2A is expressed in the posterior, and that encoding the cytoplasmic BmSerp2A in the middle silk gland region; both kinds are strongly expressed in the anterior region. The data indicate that (1) A duplication of serpin 2 gene occurred either before Trichoptera and Lepidoptera diverged as separate orders or independently in early phylogeny of either order; (2) In the caddisfly H. angustipennis, both genes are expressed specifically in the silk glands and generate proteins deposited in the silk; (3) Only one gene seems to be functional in B. mori and is expressed in a cytoplasmic and in a secreted forms in diverse organs, including the silk glands.     

Genomic organization and comparative sequence analysis of the mouse and human FRS2, FRS3 genes

The signaling adapter proteins FRS2 and FRS3 are implicated in the transmission of extracellular signals from nerve growth factor (NGF) or fibroblast growth factor (FGF) receptors to the Ras/mitogen-activated protein kinase signaling cascade. This study presents the genomic sequence and exon-intron organization of the mouse FRS2 and FRS3 loci as well as their evolutionary conservation with their human counterparts. Both FRS2 and FRS3 contain 5 coding exons spanning over 7 kb of genomic sequence with similar exon sizes and organization. Comparative genomic sequence analyses show a highly conserved genomic organization between mouse and human in both FRS2 and FRS3 genes. Non-coding sequences, highly conserved between mouse and human, were identified in the FRS3 introns that may potentially function as regulatory elements. To assay potential differences in their patterns of expression, RT-PCR analysis was used to assay FRS2 and FRS3 expression in the developing embryo and neural tube (NT) during the time of neurogenesis.     

The porcine tumor necrosis factor-encoding genes: sequence and comparative analysis

We have cloned and sequenced a 10.22-kb fragment of the genomic locus of the porcine tumor necrosis factor-encoding genes, TNF-alpha and TNF-beta. A liver genomic DNA library, partially digested with Sau3AI, was cloned into the phage lambda EMBL4 and screened with a porcine TNF-alpha cDNA probe. Analysis showed that both the TNF-alpha and TNF-beta genes were present on the cloned fragment. In addition, the cloned fragment contained about 2 kb of repetitive sequences 5' to the TNF-beta gene. The TNF genes are arranged in a tandem repeat, as is the case for the human, mouse and rabbit TNF genes. The comparison of both genes with their human homologues displayed a considerable degree of conservation (80%), suggesting an equal evolution rate.     

Comparative sequence analysis of imprinted genes between human and mouse to reveal imprinting signatures

We performed a comparative genomic sequence analysis between human and mouse for 24 imprinted genes on human chromosomes 1, 6, 7, 11, 13, 14, 15, 18, 19, and 20. The MEME program was used to search for motifs within conserved sequences among the imprinted genes and we then used the MAST program to analyze for the presence or absence of motifs in the imprinted genes and 128 nonimprinted genes. Our analysis identified 15 motifs that were significantly enriched in the imprinted genes. We generated a logistic regression model by combining multiple motifs as input variables and the 24 imprinted genes and the 128 nonimprinted genes as a training set. The accuracy, sensitivity, and specificity of our model were 98, 92, and 99%, respectively. The model was further validated by an open test on 12 additional imprinted genes. The motifs identified in this study are novel imprinting signatures, which should improve our understanding of genomic imprinting and the role of genomic imprinting in human diseases.     

The correlation between GC3 and hydropathy in human genes

A positive correlation holds between the GC level of third codon positions of human genes (GC(3)) and hydropathy of the encoded proteins. This correlation may appear counterintuitive, since it links a physical property of proteins to the base composition of 'synonymous' sites. We here establish the nontriviality of the correlation, which has recently been contested. In particular, the correlation cannot simply be a consequence of an analogous correlation for first and second codon positions, since no such correlation exists. More generally, for any explanation via two chained correlations, the intermediate property would need to be strongly correlated with hydrophobicity and/or GC(3).     

Isolation and comparative expression analysis of six MBD genes in wheat

The 5-methylcytosines (m5C) play critical roles in epigenetic control, often being recognized by proteins containing an MBD. In this study, we isolated six wheat cDNAs with open reading frame encoding putative methyl-binding domain proteins, which were designated as TaMBD1-TaMBD6, respectively. BLASTX searches and phylogenetic analysis suggested that the six TaMBD genes belonged to four (I, II, III and VIII) of the eight subclasses of MBD family. Genomic analysis showed that a 1386 bp intron was included in TaMBD1 and a 12-bp intron was found in TaMBD4. The expression profiles of the six TaMBDs were studied via Q-RT-PCR and the results indicated that the TaMBDs were differentially expressed in detected wheat tissues. It was interesting to note that 3 TaMBDs were highly expressed in dry seeds and endosperms. Moreover, the differential expression patterns of TaMBDs were observed in leaves and roots under water-stress. We concluded that multiple wheat MBD genes were present and they might play important roles in wheat growth and development, as well as in the water-stress response.     

Detecting complexes from edge-weighted PPI networks via genes expression analysis

Background

Identifying complexes from PPI networks has become a key problem to elucidate protein functions and identify signal and biological processes in a cell. Proteins binding as complexes are important roles of life activity. Accurate determination of complexes in PPI networks is crucial for understanding principles of cellular organization.

Results

We propose a novel method to identify complexes on PPI networks, based on different co-expression information. First, we use Markov Cluster Algorithm with an edge-weighting scheme to calculate complexes on PPI networks. Then, we propose some significant features, such as graph information and gene expression analysis, to filter and modify complexes predicted by Markov Cluster Algorithm. To evaluate our method, we test on two experimental yeast PPI networks.

Conclusions

On DIP network, our method has Precision and F-Measure values of 0.6004 and 0.5528. On MIPS network, our method has F-Measure and S _n values of 0.3774 and 0.3453. Comparing to existing methods, our method improves Precision value by at least 0.1752, F-Measure value by at least 0.0448, S _n value by at least 0.0771. Experiments show that our method achieves better results than some state-of-the-art methods for identifying complexes on PPI networks, with the prediction quality improved in terms of evaluation criteria.

     

Genomic sequence around butterfly wing development genes: annotation and comparative analysis

Background

Analysis of genomic sequence allows characterization of genome content and organization, and access beyond gene-coding regions for identification of functional elements. BAC libraries, where relatively large genomic regions are made readily available, are especially useful for species without a fully sequenced genome and can increase genomic coverage of phylogenetic and biological diversity. For example, no butterfly genome is yet available despite the unique genetic and biological properties of this group, such as diversified wing color patterns. The evolution and development of these patterns is being studied in a few target species, including Bicyclus anynana, where a whole-genome BAC library allows targeted access to large genomic regions.

Methodology/Principal Findings

We characterize ∼1.3 Mb of genomic sequence around 11 selected genes expressed in B. anynana developing wings. Extensive manual curation of in silico predictions, also making use of a large dataset of expressed genes for this species, identified repetitive elements and protein coding sequence, and highlighted an expansion of Alcohol dehydrogenase genes. Comparative analysis with orthologous regions of the lepidopteran reference genome allowed assessment of conservation of fine-scale synteny (with detection of new inversions and translocations) and of DNA sequence (with detection of high levels of conservation of non-coding regions around some, but not all, developmental genes).

Conclusions

The general properties and organization of the available B. anynana genomic sequence are similar to the lepidopteran reference, despite the more than 140 MY divergence. Our results lay the groundwork for further studies of new interesting findings in relation to both coding and non-coding sequence: 1) the Alcohol dehydrogenase expansion with higher similarity between the five tandemly-repeated B. anynana paralogs than with the corresponding B. mori orthologs, and 2) the high conservation of non-coding sequence around the genes wingless and Ecdysone receptor, both involved in multiple developmental processes including wing pattern formation.     

Radiation hybrid mapping and comparative sequence analysis of bovine RIG-I and MAVS genes.

Retinoic acid inducible gene I (RIG-I) and mitochondrial antiviral signaling (MAVS) proteins have recently been found to operate in a pathway for the detection and subsequent elimination of replicating viral genomes. Because of this innate immunity role, RIG-I and MAVS are candidates for studies of disease resistance. The objectives of this work were to (1) radiation hybrid (RH) map bovine RIG-I and MAVS and (2) perform comparative sequence analysis of partial genomic sequence from each gene. Using a bovine 5000(rad) RH panel, RIG-I was localized to BTA08 (LOD > 12) and MAVS was localized to BTA13 (LOD > 12). RIG-I exon 14 and partial MAVS exon five were sequenced in nine breeds and compared with available sequence from the Bovine Genome Project. RIG-I exon 14 and partial MAYS exon five were conserved in all samples examined. One T-A transversion SNP was found in intronic sequence downstream of RIG-I exon 14.     

Entire sequence of a mouse chromosomal segment containing the gene Rhced and a comparative analysis of the homologous human sequence

The mouse genomic sequence of the region containing the gene Rhced, the orthologue to the human gene RH30, was determined to elucidate the structure of Rhced and its flanking regions and to compare these with the corresponding human genomic region. Two genes, Smp1 and AK003528 (an orthologue of FLJ10747), flank Rhced. Neither sequences homologous to the characteristic nucleotide elements flanking the RHD gene in humans (rhesus boxes) nor an additional Rh gene were found within the mouse region sequenced. This result and that of a previous report demonstrate that this chromosomal region of the mouse comprises five genes (FLJ10747-RHCE-SMP1-NPD014-P29) that exhibit syntenic homology with the corresponding human region, which suggests that the RHD gene and rhesus boxes were inserted later. Evaluations of tissue distribution and subcellular localization of these genes indicate that the SMP1 orthologue has a ubiquitous tissue distribution and cytoplasmic localization, whereas AK003528 is expressed slightly higher in testis with a strong subcellular localization in the nucleus. Despite the steady improvements in the draft sequence of the human genome, this study demonstrates the continuing benefits of comparative genetic analyses in increasing our understanding of human genomic structure.