Genome-wide validation of Magnaporthe grisea gene structures based on transcription evidence

Accurate cDNA data is useful to validate gene structures in a genome. We sequenced 35 189 expressed sequence tags (ESTs) obtained from the highly destructive rice blast fungus, Magnaporthe grisea. Our custom-made computational programs mapped these ESTs on the M. grisea genome sequence, and reconstructed gene structures as well as protein-coding regions. As a result, we predicted 4480 protein-coding sequences, which were more accurate than ab initio predictions. Moreover, cross-species comparisons suggested that our predicted proteins were nearly complete. The cDNA clones obtained in this study will be important for further experimental studies. Our genome annotation is available at http://www.mg.dna.affrc.go.jp/.     

Complete sequence determination combined with analysis of transposition/site-specific recombination events to explain genetic organization of IncP-7 TOL plasmid pWW53 and related mobile genetic elements

Recent studies have indicated that the evolutionarily common catabolic gene clusters are loaded on structurally diverse toluene-catabolic (TOL) plasmids and their residing transposons. To elucidate the mechanisms supporting the diversification of catabolic plasmids and transposons, we determined here the complete 107,929 bp sequence of pWW53, a TOL plasmid from Pseudomonas putida MT53. pWW53 was found to belong to the IncP-7 incompatibility group that play important roles in the catabolism of several xenobiotics. pWW53 carried two distinct transposase-resolvase gene clusters (tnpAR modules), five short terminal inverted repeats (IRs), and three site-specific resolution (res) sites that are all typical of class II transposons. This organization of pWW53 suggested the four possible transposable regions, Tn4657 to Tn4660. The largest 86 kb region (Tn4657) spanned the three other regions, and Tn4657 and Tn4660 (62 kb) covered all of the 36 xyl genes for toluene catabolism. Our subsequent transposition experiments clarified that the three transposons, Tn4657 to Tn4659, indeed exhibit their transposability, and that pWW53 also generated another 37 kb toluene-catabolic transposon, Tn4656, which carried the two separated and inversely oriented segments of pWW53: the tnpRA-IR module of Tn4658 and a part of xyl gene clusters on Tn4657. The Tn4658 transposase was able to mediate the transposition of Tn4658, Tn4657, and Tn4656, while the Tn4659 transposase catalyzed only the transposition of Tn4659. Tn4656 was formed by the Tn4658 resolvase-mediated site-specific inversion between the two inversely oriented res sites on pWW53. These findings and comparison with other catabolic plasmids clearly indicate multiple copies of transposition-related genes and sites on one plasmid and their recombination activities contribute greatly to the diversification of plasmid structures as well as wide dissemination of the evolutionary common gene clusters in various plasmids.     

Comparative genomic analysis of dinucleotide repeats in Tritryps

The protozoans Trypanosoma cruzi, Trypanosoma brucei and Leishmania major (Tritryps), are evolutionarily ancient eukaryotes which cause worldwide human parasitosis. They present unique biological features. Indeed, canonical DNA/RNA cis-acting elements remain mostly elusive. Repetitive sequences, originally considered as selfish DNA, have been lately recognized as potentially important functional sequence elements in cell biology. In particular, the dinucleotide patterns have been related to genome compartmentalization, gene evolution and gene expression regulation. Thus, we perform a comparative analysis of the occurrence, length and location of dinucleotide repeats (DRs) in the Tritryp genomes and their putative associations with known biological processes. We observe that most types of DRs are more abundant than would be expected by chance. Complementary DRs usually display asymmetrical strand distribution, favoring TT and GT repeats in the coding strands. In addition, we find that GT repeats are among the longest DRs in the three genomes. We also show that specific DRs are non-uniformly distributed along the polycistronic unit, decreasing toward its boundaries. Distinctive non-uniform density patterns were also found in the intergenic regions, with predominance at the vicinity of the ORFs. These findings further support that DRs may control genome structure and gene expression.     

Identification and characterization of yak (Bos grunniens) b-Boule gene and its alternative splice variants

Boule is responsible for meiotic arrest of sperms and male sterility during mammalian spermatogenesis. In the present study, we first identified yak b-Boule gene and its two alternative splice variants. The full length coding region of yak b-Boule is 888 bp and encodes a 295-amino acid protein with a typical RNA-recognition motif (RRM) and a Deleted in Azoospermia (DAZ) repetitive sequence motif. Two alternative splice variants of yak b-Boule were generated following the consensus “GT-AG” rule and named b-Boule1 (36 bp deletion in exon 3) and b-Boule2 (deletion of integral exon 7), respectively. In male yak, b-Boule, b-Boule1 and b-Boule2 were found to be exclusively expressed in the testes at a ratio of 81:0.1:1. Intriguingly, the mRNA expression levels of b-Boule and b-Boule1 in yak testis were significantly higher than those in cattle–yak, although no significant difference was observed for b-Boule2 expression between the yak and cattle–yak. These results suggest that b-Boule gene, which is partially regulated by alternative splicing, may be involved in the process of yak spermatogenesis.     

Identification and characterization of a differentially expressed protein (CAPZB) in skeletal muscle between Meishan and Large White pigs

Actin capping protein beta (CAPZB) protein was identified with considerable differences in the longissimus dorsi muscle between Large White and Meishan pigs using proteomics approach. However, in pigs, the information on CAPZB is very limited. In this study, we cloned and characterized the porcine actin capping protein beta (CAPZB) gene. In addition, we present two novel porcine CAPZB splice variants CAPZB1 and CAPZB2. CAPZB1 was expressed in all twenty tissues. However, CAPZB2 was predominantly expressed in the skeletal muscle and heart. In addition, the two isoforms had different expression profiles during the skeletal muscle development and between breeds. Moreover, the SNP T394G was identified in the coding region of the CAPZB gene, which was significantly associated with the carcass traits including the LFW, CFW, SFT and LEA. Data presented in our study suggests that the CAPZB gene may be a candidate gene of meat production trait and provides useful information for further studies on its roles in porcine skeletal muscle.     

Molecular characterization and expression patterns of the actinin-associated LIM protein (ALP) subfamily genes in porcine skeletal muscle

The actinin-associated LIM protein (ALP) subfamily has important functions in cell signal transduction, cell proliferation, and integration of cytoskeletal architecture. To detect their functions in pig skeletal muscle, we cloned and characterized the pig ALP subfamily genes, drew their genomic structure maps, and detected their tissue expression patterns. We identified a new spliced variant of PDLIM3 in pig skeletal muscle and named it as PDLIM3-4, which was only expressed in the heart and skeletal muscle. Our results showed that PDLIM3-4 was expressed in adult pig skeletal muscle with the highest expression level, and both PDLIM3-4 isoform and PDLIM4 had different expression profiles during the prenatal and postnatal stages of skeletal muscle development among the three pig breeds. These studies provide useful information for further research on the functions of pig ALP subfamily genes in skeletal muscle development.     

Inactivation of the MSLNL gene encoding mesothelin-like protein during African great ape evolution

Loss of gene function is implicated in the emergence of novel phenotypes during organism evolution. Here, we report the inactivation of the MSLNL gene encoding mesothelin-like protein in African great ape evolution. Human MSLNL has a nonsense mutation in exon 10 and two polymorphic mutations: a frameshift in exon 3 and a nonsense codon in exon 8. The gorilla gene also shows multiple deleterious mutations, including a premature stop codon, a deletion, and a splice site mutation. Molecular evolutionary analysis indicated relaxed selection pressure on MSLNL in African great ape lineages, which suggested that MSLNL might have become inactivated before the divergence of human, chimpanzee and gorilla. The mouse Mslnl gene is highly expressed in olfactory epithelium and moderately expressed in several other tissues. We propose that the loss of MSLNL may be associated with the evolution of the olfactory system in African great apes including human.     

Polymorphisms in the promoter region of ESR2 gene and susceptibility to ovarian cancer

Susceptibility to ovarian cancer might be affected by genetic variations in genes involved in estrogen biosynthesis, metabolism or signal transduction. In this study we tested the hypothesis that single nucleotide polymorphisms (SNPs) in the promoter of human ESR2 gene, coding for estrogen receptor β, may be associated with increased risk for ovarian cancer. Three SNPs in the promoter region of human ESR2 gene were genotyped by means of allele-specific tetra-primer PCR. A total of 184 ovarian cancer cases and the same numbers of controls were included in the study. With regard to homozygous analysis, the AA genotype of SNP rs3020449 was found to be significantly more frequent in ovarian cancer cases staged as FIGO III + IV than in cases staged as I + II (OR 2.717, p = 0.027). With regard to allele frequency, the G allele of this SNP was less frequent in FIGO I + II cases than in cases with higher FIGO stages (OR 1.739, p = 0.018). With regard to genotype frequency, allele frequency, allele positivity or haplotype frequency of SNPs rs2987983, rs3020449 and rs3020450 we did not observe a significant difference between the cancer and the control group. Our data suggest that SNPs in the promoter region of ESR2 gene do not affect susceptibility to ovarian cancer, but SNP rs3020449 might affect progression of this disease.     

Characterization and identification of the integrin family in silkworm, Bombyx mori

As an important economic insect, Bombyx mori is also a useful model organism for lepidopteran insect. Integrins are evolutionarily conserved from sponges to humans, and play vital roles in many physiological and pathological processes. To explore their diverse functions of integrins in insect, eleven integrins including six α and five β subunits were cloned and characterized from silkworm. Our results showed that integrins from silkworm own more family members compared to other invertebrates. Among those α subunits, integrins α1, α2, and the other four subunits belong to PS1, PS2, and PS3 groups, respectively. The β subunits mainly gather in the insect βν group except the β1 subunit which belongs to the insect β group. Expression profiles demonstrated that the integrins exhibited distinct patterns, but were mainly expressed in hemocytes. α1 and β2 subunits are the predominant ones either in the embryogenesis or larva stages. Interestingly, integrins were significantly up-regulated after stimulated by 20-hydroxyecdysone (20-E) in vivo. These results indicate that integrins perform diverse functions in hemocytes of silkworm. Overall, our results provide a new insight into the functional and evolutionary features of integrins.     

Molecular cloning and characterization of the β-catenin gene from fine-wool sheep

β-Catenin is an evolutionarily conserved molecule that functions as a crucial effector in both cell-to-cell adhesion and Wnt signaling. To gain a better understanding of its role in the development of hair follicles, we cloned the cDNA sequence of the β-catenin gene from the skin of Aohan fine-wool sheep and performed a variety of bioinformatics analyses. We obtained the full-length sequence, which was 4573-bp long and contained a 2346-bp open reading frame encoding a protein of 781 amino acids. The protein had a predicted molecular weight of 85.4 kDa and a theoretical isoelectric point of 5.57. Domain architecture analysis of the β-catenin protein revealed an armadillo repeat region, which is a common feature of β-catenin in other species. The ovine β-catenin gene shares 97.91%, 94.25%, 94.59%, 83.89%, and 89.39% sequence identity with its homologs in Bos taurus, Homo sapiens, Sus scrofa, Gallus gallus, and Mus musculus, respectively, while the amino acid sequence is more than 99% identical with each of these species. The expression of β-catenin mRNA was detected in the heart, liver, spleen, lung, kidney, skin, muscle, and adipose tissue. Expression levels were maximal in the lung and minimal in the muscle, and the difference in expression in these tissues was significant (P < 0.01). Western blot analysis revealed the presence of the β-catenin protein in all tissues examined; expression was lowest in the skin and adipose tissues.     

A shotgun approach to discovering and reconstructing consensus retrotransposons ex novo from dense contigs of short sequences derived from Genbank Genome Survey Sequence database records

Retrotransposons constitute the majority of pseudogenic protein coding regions of most eukaryotic genomes. Most genomes carry tens to thousands of retrotransposon copies derived from dozens of distinct families, but most if not all of these copies are non-functional and contain disabling mutations, including large numbers of indels. Until recently, most regions rich in these elements were virtually ignored in all but the most complete genome sequencing projects, and the full extent of their impact on the structure and function of the genomes of higher eukaryotes was under-appreciated. Even when new retrotransposons are encountered and annotated by automated gene finding programs and similarity searches, coding regions are treated as exons and invariably and not surprisingly mistranslated because of numerous frameshift mutations and large indels. Very few functional retrotransposons contain introns, as in silico annotations imply. While many repetitive DNA consensus sequences have been assembled from collections of largely full-length copies using full-length templates, we have shown that repetitive DNA consensus sequence contigs representing long, moderately high copy-number elements can also be generated ex novo in the absence of templates from very short overlapping sequences. We have devised an in silico strategy to recover and reconstruct consensus sequences of elements up to 20,000 bp by building dense contigs of hundreds of overlapping 400 to 900-bp records found in the Genbank Genome Survey Sequence database. The results are hypothetical ancestral sequences that encode elements that appear to be fully functional with intact open reading frames and other conserved features.     

A temperature-sensitive mutation in the WD repeat-containing protein Smu1 is related to maintenance of chromosome integrity

Temperature-sensitive CHO-K1 mutant cell line tsTM18 exhibits chromosomal instability and cell cycle arrest at S and G2 phases with decreased DNA synthesis at the nonpermissive temperature, 39 degrees C. To identify the causative mutation, we fused tsTM18 cells with normal human cells to generate hybrids carrying fragments of human chromosomes. Analysis of chromosome content of temperature-resistant transformants and introduction of a bacterial artificial chromosome containing part of human chromosome 9 led to isolation of the human SMU1 gene. Comparison of sequences of the Smu1 gene from wild-type and mutant cells revealed that the mutant phenotype is caused by a G-to-A transition that yields a gly-to-arg substitution at position 489 in hamster Smu1. The substituted glycine is located in the WD-repeat domain of Smu1. Single-stranded DNA accumulated in the nuclei of mutant cells at 39 degrees C. Furthermore, cdc2 kinase was not activated during G2 phase, and there was no chromosome segregation due to incomplete assembly of the spindle during M phase. Thus, Smu1 appears to be involved directly or indirectly in DNA replication, activation of cdc2 kinase, spindle assembly, and maintenance of chromosome integrity, reflecting the important roles of Smu1 in cellular function.     

Common variant in FUT2 gene is associated with levels of vitamin B12 in Indian population

Vitamin B₁₂ is an essential micronutrient synthesized by microorganisms. Mammals including humans have evolved ways for transport and absorption of this vitamin. Deficiency of vitamin B₁₂ (either due to low intake or polymorphism in genes involved in absorption and intracellular transport of this vitamin) has been associated with various complex diseases. Genome-wide association studies have recently identified several common single nucleotide polymorphisms (SNPs) in fucosyl transferase 2 gene (FUT2) to be associated with levels of vitamin B₁₂—the strongest association was with a non-synonymous SNP rs602662 in this gene. In the present study, we attempted to replicate the association of this SNP (rs602662) in an Indian population since a significant proportion has been reported to have low levels of vitamin B₁₂ in this population. A total of 1146 individuals were genotyped for this SNP using a single base extension method and association with levels of vitamin B₁₂ was assessed in these individuals. Regression analysis was performed to analyze the association considering various confounding factors like for age, sex, diet, hypertension, diabetes mellitus and coronary artery disease status. We found that the SNP rs602662 was significantly associated with the levels of vitamin B₁₂ (p value < 0.0001). We also found that individuals adhering to a vegetarian diet with GG (homozygous major genotype) have significantly lower levels of vitamin B₁₂ in these individuals. Thus, our study reveals that vegetarian diet along with polymorphism in the FUT2 gene may contribute significantly to the high prevalence of vitamin B₁₂ deficiency in India.     

Molecular characterization of two y-type high molecular weight glutenin subunit alleles 1Ay12 and 1Ay8 from cultivated einkorn wheat (Triticum monococcum ssp. monococcum)

Two y-type high molecular weight glutenin subunits (HMW-GSs) 1Ay12^? and 1Ay8^? from the two accessions PI560720 and PI345186 of cultivated einkorn wheat (Triticum monococcum ssp. monococcum, AA, 2n = 2x = 14), were identified by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE). The mobility of 1Ay12^? and 1Ay8^? was similar to that of 1Dy12 and 1By8 from common wheat Chinese Spring, respectively. Their ORFs respectively consisted of 1812 bp and 1935 bp, encoding 602 and 643 amino acid residues with the four typical structural domains of HMW-GS including signal peptide, conserved N-, and C-terminal and central repetitive domains. Compared with the most similar active 1Ay alleles previous published, there were a total of 15 SNPs and 2 InDels in them. Their encoding functions were confirmed by successful heterogeneous expression. The two novel 1Ay alleles were named as 1Ay12^? and 1Ay8^? with the accession No. JQ318694 and JQ318695 in GenBank, respectively. The two alleles were classed into the two distinct groups, Phe-type and Cys-type, which might be relevant to the differentiation of Glu-A1-2 alleles. Of which, 1Ay8^? belonged to Cys-type group, and its protein possessed an additional conserved cysteine residue in central repetitive region besides the six common ones in N- and C-terminal regions of Phe-type group, and was the second longest in all the known active 1Ay alleles. These results suggested that the subunit 1Ay8^? of cultivated einkorn wheat accession PI345186 might have a potential ability to strengthen the gluten polymer interactions and be a valuable genetic resource for wheat quality improvement.     

IS elements in Aliivibrio salmonicida LFI1238: Occurrence,variability and impact on adaptability

Insertion sequence (IS) elements are short, self-replicating DNA sequences that are capable of efficiently spreading over the host genome. Possessing varied integration specificity IS elements are capable of the irreversible inactivation of genes, which diversifies the pool of intact genetic determinants in host populations. In the current study, we performed a complex analysis of IS elements (Vsa IS) in the previously sequenced genome of Aliivibrio salmonicida LFI1238 and proposed a model of the spread of the Vsa IS elements over the genome of this microorganism. Along with the prediction of the integration sites for Vsa IS elements, the current study provides an overview of the properties of A. salmonicida IS elements, as well as information regarding their occurrence in different bacterial classes. An analysis of individual alleles of the IS elements has allowed us to depict a history of the accumulation of mutations and to describe distinctive microevolution lines for actively transposing Vsa IS elements in the genome of A. salmonicida LFI1238. Our results demonstrate the high importance of the dead end microevolution of actively transposing Vsa IS elements for the inactivation of genes in A. salmonicida LFI1238.     

Ankyrin-G in skeletal muscle: tissue-specific alternative splicing contributes to the complexity of the sarcolemmal cytoskeleton

Ankyrins are versatile adaptor proteins that join the spectrin-based cytoskeleton to transmembrane proteins, and have roles in organizing the microstructure of cell membranes. Molecular diversity of ankyrins in mammals arises from extensive alternative splicing of the products of three genes. There has been no systematic analysis of the diversity of expression of ankyrins-G, the widely expressed Ank3 gene products, in a complex tissue. We previously described Ank(G107), the first muscle-specific ankyrin-G. Here, we combined cDNA and database analyses to gain novel insight into the ankyrins-G of skeletal muscle. We find: (i) that Ank3 is composed of at least 53 exons, many of which are subject to tissue-specific splicing; (ii) five novel full-length cDNAs encoding two canonical (Ank(G197), Ank(G217)) and three small isoforms (Ank(G109), Ank(G128), Ank(G130)) bring to six the number of ankyrins-G expressed in skeletal muscle; (iii) a 76-residue insert in the C-terminal domain is a 'signature' for muscle ankyrins; (iv) variably spliced sequences of 17/18 and 195 residues increase diversity in the C-terminal domains. Comparison of endogenous ankyrins-G with in vitro translated cDNAs revealed that small ankyrins account for the majority of the immunoreactivity for ankyrin-G in soleus muscle. The small ankyrins, when expressed in vivo in the rat muscle, are all targeted to sarcolemmal costameres. Our results demonstrate the tissue-dependent alternative splicing of Ank3 in skeletal muscle and point to novel functions of small ankyrins-G in organizing microdomains of the plasma membrane.     

Optimizing high-resolution melting analysis for the detection of mutations of GPR30/GPER-1 in breast cancer

G protein-coupled receptor 30/G protein estrogen receptor-1 (GPR30/GPER-1) is a novel membrane receptor for estrogen whose mRNA is expressed at high levels in estrogen-dependent cells such as breast cancer cell lines. However, mutations in GRP30 related to diseases remain unreported. To detect unknown mutations in the GPR30 open reading frame (ORF) quickly, the experimental conditions for high-resolution melting (HRM) analysis were examined for PCR primers, Taq polymerases, saturation DNA binding dyes, Mg(2+) concentration, and normalized temperatures. Nine known SNPs and 13 artificial point mutations within the GPR30 ORF, as well as single nucleotide variants in DNA extracted from subjects with breast cancers were tested under the optimal experimental conditions. The combination of Expand High Fidelity(PLUS) and SYTO9 in the presence of 2.0 mM MgCl(2) produced the best separation in melting curves of mutations in all regions of the GPR30 ORF. Under these experimental conditions, the mutations were clearly detected in both heterozygotes and homozygotes. HRM analysis of GPR30 using genomic DNA from subjects with breast cancers showed a novel single nucleotide variant, 111C>T in GPR30 and 4 known SNPs. The experimental conditions determined in this study for HRM analysis are useful for high throughput assays to detect unknown mutations within the GPR30 ORF.