Intron Size Research Articles

Higher frequency of intron loss from the promoter proximally paused genes of Drosophila melanogaster

Although intron losses have been widely reported, it is not clear whether they are neutral and therefore random or driven by positive selection. Intron transcription and splicing are time-consuming and can delay the expression of its host gene. For genes that must be activated quickly to respond to physiological or stress signals, intron delay may be deleterious. Promoter proximally paused (PPP) genes are a group of rapidly expressed genes. To respond quickly to activation signals, they generally initiate transcription competently but stall after synthesizing a short RNA. In this study, performed in Drosophila melanogaster, the PPP genes were found to have a significantly higher rate of intron loss than control genes. However, further analysis did not find more significant shrinkage of intron size in PPP genes. Referring to previous studies on the rates of transcription and splicing and to the time saved by deletion of the introns from mouse gene Hes7, it is here suggested that transcription delay is comparable to splicing delay only when the intron is 28.5 kb or larger, which is greater in size than 95% of vertebrate introns, 99.5% of Drosophila introns, and all the annotated introns of Saccharomyces cerevisiae and Arabidopsis thaliana. Delays in intron splicing are probably a selective force, promoting intron loss from quickly expressed genes. In other genes, it may have been an exaptation during the emergency of developmental clocks.

Genome-Wide Analyses of a Plant-Specific LIM-Domain Gene Family Implicate Its Evolutionary Role in Plant Diversification

The Arabidopsis DA1 genes appear to have multiple functions in regulating organ size and abiotic stress response, but the biological roles of its closely related genes remain unknown. Evolutionary analyses might provide some clues to aid in an understanding of their functional diversification. In this work, we characterized the molecular evolution and expressional diversification of DA1-like genes. Surveying 354 sequenced genomes revealed 142 DA1-like genes only in plants, indicating plant-specificity of these genes. The DA1-like protein modular structure was composed of two UIMs (ubiquitin interaction motifs), one LIM-domain (from lin-11, isl-1, and mec-3), and a conserved C-terminal, and was distinguishable from the already defined three groups of LIM-domain proteins. We further found that the DA1-like genes diverged into Classes I and II at the ancestor of seed plants and acquired 13 clade-specific residues during their evolutionary history. Moreover, diverse intron size evolution was noted following the transition from size-expandable introns to minimal ones, accompanying the emergence and diversification of angiosperms. Functional diversification as it relates to gene expression was further investigated in soybean. Glycine max DA1 genes showed diverse tissues expression patterns during development and had substantially varied abiotic stress response expression. Thus, variations in the coding regions, intron size, and gene expression contributed to the functional diversification of this gene family. Our data suggest that the evolution of the DA1-like genes facilitated the development of diverse molecular and functional diversification patterns to accompany the successful radiation of plants into diverse environments during evolution.

Intron evolution in a phylogenetic perspective: Divergent trends in the two copies of the duplicated def gene in Impatiens L. (Balsaminaceae)

AbstractThe history of MADS box genes is well‐known in angiosperms. While duplication events and gene losses occur frequently, gene structure and intron positions are very conserved. We investigated all six introns in a duplicated MADS box gene (deficiens, def) in selected Impatiens taxa, thereby assessing intron features. For the first time, our study provides a comparison of molecular changes in all introns of a gene from a phylogenetic perspective. Interestingly, a uniform pattern of molecular evolution in the introns of each copy was not observed, but intron length increases, decreases, and size retention can be found in each copy. A tendency to accumulate long autapomorphic indels is also present, thus, a longer intron length does not reflect a higher number of parsimony‐informative characters. Substitution rates vary between introns of each gene copy. While four of the six introns of def1 exhibit a change in their substitution rate, five of the six def2 introns maintain their rates throughout the genus albeit at different levels. In MADS box genes several regulatory sequences are found residing in introns. Thus, presence of putative regulatory motifs was investigated. Most of them are not conserved in position and usually present in only one of the gene copies. In addition, the potential for phylogenetic reconstruction of introns in both def copies is shortly discussed.

Genes with a large intronic burden show greater evolutionary conservation on the protein level

BackgroundThe existence of introns in eukaryotic genes is believed to provide an evolutionary advantage by increasing protein diversity through exon shuffling and alternative splicing. However, this eukaryotic feature is associated with the necessity of exclusion of intronic sequences, which requires considerable energy expenditure and can lead to splicing errors. The relationship between intronic burden and evolution is poorly understood. The goal of this study was to analyze the relationship between the intronic burden and the level of evolutionary conservation of the gene.ResultsWe found a positive correlation between the level of evolutionary conservation of a gene and its intronic burden. The level of evolutionary conservation was estimated using the conservation index (CI). The CI value was determined on the basis of the most distant ortholog of the human protein sequence and ranged from 0 (the gene was unique to the human genome) to 9 (an ortholog of the human gene was detected in plants). In multivariable model, both the number of introns and total intron size remained significant predictors of CI. We also found that the number of alternative splice variants was positively correlated with CI.The expression level of a gene was negatively correlated with the number of introns and total size of intronic region. Genes with a greater intronic burden had lower density of missense and nonsense mutations in the coding regions of the gene, which suggests that they are under a stronger pressure from purifying selection.ConclusionsWe identified a positive association between intronic burden and CI. One of the possible explanations of this is the idea of a cost-benefits balance. Evolutionarily conserved (functionally important) genes can “afford” the negative consequences of maintaining multiple introns because these consequences are outweighed by the benefit of maintaining the gene. Evolutionarily conserved and functionally important genes may use introns to create novel splice variants to tune the gene function to developmental stage and tissue type.

PCR survey of 50 introns in animals: Cross-amplification of homologous EPIC loci in eight non-bilaterian, protostome and deuterostome phyla

Exon Primed Intron Crossing (EPIC) markers provide molecular tools that are susceptible to be variable within species while remaining amplifiable by PCR using potentially universal primers. In this study we tested the possibility of obtaining PCR products from 50 EPIC markers on 23 species belonging to seven different phyla (Porifera, Cnidaria, Arthropoda, Nematoda, Mollusca, Annelida, Echinodermata) using 70 new primer pairs. A previous study had identified and tested those loci in a dozen species, including another phylum, Urochordata (Chenuil et al., 2010). Results were contrasted among species. The best results were achieved with the oyster (Mollusca) where 28 loci provided amplicons susceptible to contain an intron according to their size. This was however not the case with the other mollusk Crepidula fornicata, which seems to have undergone a reduction in intron number or intron size. In the Porifera, 13 loci appeared susceptible to contain an intron, a surprisingly high number for this phylum considering its phylogenetic distance with genomic data used to design the primers. For two cnidarian species, numerous loci (24) were obtained. Ecdysozoan phyla (arthropods and nematodes) proved less successful than others as expected considering reports of their rapid rate of genome evolution and the worst results were obtained for several arthropods. Some general patterns among phyla arose, and we discuss how the results of this EPIC survey may give new insights into genome evolution of the study species. This work confirms that this set of EPIC loci provides an easy-to-use toolbox to identify genetic markers potentially useful for population genetics, phylogeography or phylogenetic studies for a large panel of metazoan species. We then argue that obtaining diploid sequence genotypes for these loci became simple and affordable owing to Next-Generation Sequencing development. Species surveyed in this study belong to several genera (Acanthaster, Alvinocaris, Aplysina, Aurelia, Crepidula, Eunicella, Hediste, Hemimysis, Litoditis, Lophelia, Mesopodopsis, Mya, Ophiocten, Ophioderma, Ostrea, Pelagia, Platynereis, Rhizostoma, Rimicaris), two of them, belonging to the family Vesicomydae and Eunicidae, could not be determined at the genus level.

Transcription factor binding to Caenorhabditis elegans first introns reveals lack of redundancy with gene promoters

Gene expression is controlled through the binding of transcription factors (TFs) to regulatory genomic regions. First introns are longer than other introns in multiple eukaryotic species and are under selective constraint. Here we explore the importance of first introns in TF binding in the nematode Caenorhabditis elegans by combining computational predictions and experimentally derived TF–DNA interaction data. We found that first introns of C. elegans genes, particularly those for families enriched in long first introns, are more conserved in length, have more conserved predicted TF interactions and are bound by more TFs than other introns. We detected a significant positive correlation between first intron size and the number of TF interactions obtained from chromatin immunoprecipitation assays or determined by yeast one-hybrid assays. TFs that bind first introns are largely different from those binding promoters, suggesting that the different interactions are complementary rather than redundant. By combining first intron and promoter interactions, we found that genes that share a large fraction of TF interactions are more likely to be co-expressed than when only TF interactions with promoters are considered. Altogether, our data suggest that C. elegans gene regulation may be additive through the combined effects of multiple regulatory regions.

Systematic analysis of intron size and abundance parameters in diverse lineages

All eukaryotic genomes have genes with introns in variable sizes. As far as spliceosomal introns are concerned, there are at least three basic parameters to stratify introns across diverse eukaryotic taxa: size, number, and sequence context. The number parameter is highly variable in lower eukaryotes, especially among protozoan and fungal species, which ranges from less than 4% to 78% of the genes. Over greater evolutionary time scales, the number parameter undoubtedly increases as observed in higher plants and higher vertebrates, reaching greater than 12.5 exons per gene in average among mammalian genomes. The size parameter is more complex, where multiple modes appear at work. Aside from intronless genes, there are three other types of intron-containing genes: half-sized, minimal, and size-expandable introns. The half-sized introns have only been found in a limited number of genomes among protozoan and fungal lineages and the other two types are prevalent in all animal and plant genomes. Among the size-expandable introns, the sizes of plant introns are expansion-limited in that the large introns exceeding 1000 bp are fewer in numbers and transposon-free as compared to the large introns among animals, where the larger introns are filled with transposable elements and appear expansion-flexible, reaching several kilobasepairs (kbp) and even thousands of kbp in size. Most of the intron parameters can be studied as signatures of the specific splicing machineries of different eukaryotic lineages and are highly relevant to the regulation of gene expression and functionality. In particular, the transcription-splicing-export coupling of eukaryotic intron dispensing leads to a working hypothesis that all intron parameters are evolved to be efficient and function-related in processing and routing the spliced transcripts.

Phylogenetics and Gene Structure Dynamics of Polygalacturonase Genes in Aspergillus and Neurospora crassa.

Polygalacturonase (PG) gene is a typical gene family present in eukaryotes. Forty-nine PGs were mined from the genomes of Neurospora crassa and five Aspergillus species. The PGs were classified into 3 clades such as clade 1 for rhamno-PGs, clade 2 for exo-PGs and clade 3 for exo- and endo-PGs, which were further grouped into 13 sub-clades based on the polypeptide sequence similarity. In gene structure analysis, a total of 124 introns were present in 44 genes and five genes lacked introns to give an average of 2.5 introns per gene. Intron phase distribution was 64.5% for phase 0, 21.8% for phase 1, and 13.7% for phase 2, respectively. The introns varied in their sequences and their lengths ranged from 20 bp to 424 bp with an average of 65.9 bp, which is approximately half the size of introns in other fungal genes. There were 29 homologous intron blocks and 26 of those were sub-clade specific. Intron losses were counted in 18 introns in which no obvious phase preference for intron loss was observed. Eighteen introns were placed at novel positions, which is considerably higher than those of plant PGs. In an evolutionary sense both intron loss and gain must have taken place for shaping the current PGs in these fungi. Together with the small intron size, low conservation of homologous intron blocks and higher number of novel introns, PGs of fungal species seem to have recently undergone highly dynamic evolution.

Large Scale Full-Length cDNA Sequencing Reveals a Unique Genomic Landscape in a Lepidopteran Model Insect, Bombyx mori

The establishment of a complete genomic sequence of silkworm, the model species of Lepidoptera, laid a foundation for its functional genomics. A more complete annotation of the genome will benefit functional and comparative studies and accelerate extensive industrial applications for this insect. To realize these goals, we embarked upon a large-scale full-length cDNA collection from 21 full-length cDNA libraries derived from 14 tissues of the domesticated silkworm and performed full sequencing by primer walking for 11,104 full-length cDNAs. The large average intron size was 1904 bp, resulting from a high accumulation of transposons. Using gene models predicted by GLEAN and published mRNAs, we identified 16,823 gene loci on the silkworm genome assembly. Orthology analysis of 153 species, including 11 insects, revealed that among three Lepidoptera including Monarch and Heliconius butterflies, the 403 largest silkworm-specific genes were composed mainly of protective immunity, hormone-related, and characteristic structural proteins. Analysis of testis-/ovary-specific genes revealed distinctive features of sexual dimorphism, including depletion of ovary-specific genes on the Z chromosome in contrast to an enrichment of testis-specific genes. More than 40% of genes expressed in specific tissues mapped in tissue-specific chromosomal clusters. The newly obtained FL-cDNA sequences enabled us to annotate the genome of this lepidopteran model insect more accurately, enhancing genomic and functional studies of Lepidoptera and comparative analyses with other insect orders, and yielding new insights into the evolution and organization of lepidopteran-specific genes.

First Extensive Examination of Genome Size in Phylum Brachiopoda (Lamp Shells) Collected from Japan

Genome size (C value) is a fundamental characteristic of every species and is very important for the progress of cytogenetic, genomic, and phylogenic studies. However, information on the C value of phylum Brachiopoda is scarce. In this study, we collected 8 brachiopod species from Japan, including 5 from class Articulata, and determined their C values by flow cytometry. The mean C values for these 8 species-namely, Lingula anatina, Lingula reevii, Discradisca stella, Terebratulina crossei, Laqueus rubellus, Laqueus blanfordi, Terebratalia coreanica, and Platidia japonica—were 0.41 pg, 0.41 pg, 0.46 pg, 0.33 pg, 0.44 pg, 0.42 pg, 0.37 pg, and 0.31 pg, respectively. Although the C values were examined across the various taxa within Brachiopoda, we detected very little variation (approximately 1.5-fold, from 0.31 pg in P. japonica to 0.46 pg in D. stella). This low C value variation suggests that the occurrences of evolutionary events such as changing the number of transposable elements and intron size, polyploidy, and DNA loss—all of which could cause a change in genome size—are fewer than in other taxa of the animal kingdom.

The Conserved Intronic Cleavage and Polyadenylation Site of CstF-77 Gene Imparts Control of 3′ End Processing Activity through Feedback Autoregulation and by U1 snRNP

The human gene encoding the cleavage/polyadenylation (C/P) factor CstF-77 contains 21 exons. However, intron 3 (In3) accounts for nearly half of the gene region, and contains a C/P site (pA) with medium strength, leading to short mRNA isoforms with no apparent protein products. This intron contains a weak 5′ splice site (5′SS), opposite to the general trend for large introns in the human genome. Importantly, the intron size and strengths of 5′SS and pA are all highly conserved across vertebrates, and perturbation of these parameters drastically alters intronic C/P. We found that the usage of In3 pA is responsive to the expression level of CstF-77 as well as several other C/P factors, indicating it attenuates the expression of CstF-77 via a negative feedback mechanism. Significantly, intronic C/P of CstF-77 pre-mRNA correlates with global 3′UTR length across cells and tissues. In addition, inhibition of U1 snRNP also leads to regulation of the usage of In3 pA, suggesting that the C/P activity in the cell can be cross-regulated by splicing, leading to coordination between these two processes. Importantly, perturbation of CstF-77 expression leads to widespread alternative cleavage and polyadenylation (APA) and disturbance of cell proliferation and differentiation. Thus, the conserved intronic pA of the CstF-77 gene may function as a sensor for cellular C/P and splicing activities, controlling the homeostasis of CstF-77 and C/P activity and impacting cell proliferation and differentiation.

Genomic organization and expression analysis of a farnesyl diphosphate synthase gene (FPPS2) in apples (Malus domestica Borkh.)

A farnesyl diphosphate synthase gene (FPPS2), which contains 11 introns and 12 exons, was isolated from the apple cultivar “White Winter Pearmain”. When it was compared to our previously reported FPPS1, its each intron size was different, its each exon size was the same as that of FPPS1 gene, 30 nucleotide differences were found in its coding sequence. Based on these nucleotide differences, specific primers were designed to perform expression analysis; the results showed that it expressed in both fruit and leaf, its expression level was obviously lower than that of FPPS1 gene in fruit which was stored at 4°C for 5weeks. This is the first report concerning two FPPS genes and their expression comparison in apples.

A simple model to explain evolutionary trends of eukaryotic gene architecture and expression

Enormous phylogenetic variation exists in the number and sizes of introns in protein-coding genes. Although some consideration has been given to the underlying role of the population-genetic environment in defining such patterns, the influence of the intracellular environment remains virtually unexplored. Drawing from observations on interactions between co-transcriptional processes involved in splicing and mRNA 3'-end formation, a mechanistic model is proposed for splice-site recognition that challenges the commonly accepted intron- and exon-definition models. Under the suggested model, splicing factors that outcompete 3'-end processing factors for access to intronic binding sites concurrently favor the recruitment of 3'-end processing factors at the pre-mRNA tail. This hypothesis sheds new light on observations such as the intron-mediated enhancement of gene expression and the negative correlation between intron length and levels of gene expression.

OLego: fast and sensitive mapping of spliced mRNA-Seq reads using small seeds

A crucial step in analyzing mRNA-Seq data is to accurately and efficiently map hundreds of millions of reads to the reference genome and exon junctions. Here we present OLego, an algorithm specifically designed for de novo mapping of spliced mRNA-Seq reads. OLego adopts a multiple-seed-and-extend scheme, and does not rely on a separate external aligner. It achieves high sensitivity of junction detection by strategic searches with small seeds (∼14 nt for mammalian genomes). To improve accuracy and resolve ambiguous mapping at junctions, OLego uses a built-in statistical model to score exon junctions by splice-site strength and intron size. Burrows–Wheeler transform is used in multiple steps of the algorithm to efficiently map seeds, locate junctions and identify small exons. OLego is implemented in C++ with fully multithreaded execution, and allows fast processing of large-scale data. We systematically evaluated the performance of OLego in comparison with published tools using both simulated and real data. OLego demonstrated better sensitivity, higher or comparable accuracy and substantially improved speed. OLego also identified hundreds of novel micro-exons (<30 nt) in the mouse transcriptome, many of which are phylogenetically conserved and can be validated experimentally in vivo. OLego is freely available at http://zhanglab.c2b2.columbia.edu/index.php/OLego.

A DNA Microarray for the Detection of Point Mutations and Copy Number Variation Causing Familial Hypercholesterolemia in Europe

To facilitate genetic cascade screening for familial hypercholesterolemia (FH) in Europe, two versions (7 and 9) of a DNA microarray were developed to detect the most frequent point mutations in the low-density lipoprotein receptor (LDLR), apolipoprotein B (APOB), and proprotein convertase subtilisin/kexin 9 (PCSK9) genes. The design of these microarrays is based on LIPOchip, version 4, which detects 191 LDLR and APOB mutations identified in Spanish patients with FH. A major improvement of LIPOchip, versions 7 and 9, is the ability to detect copy number variation (deletions or duplications of entire exons) in LDLR, thus abolishing the need to perform multiplex ligase-dependent probe amplification in patients with FH. The aim of this study was to validate a tool capable of detecting point mutations and copy number variations simultaneously and to evaluate its use and the newly developed software for analysis in clinical practice by reanalysis of several patients with known mutations causing FH. With the help of these validations, several aspects were analyzed, improved, and implemented in a newer version, which was evaluated through an internal validation.

Gene make-up: rapid and massive intron gains after horizontal transfer of a bacterial α-amylase gene to Basidiomycetes

BackgroundIncreasing genome data show that introns, a hallmark of eukaryotes, already existed at a high density in the last common ancestor of extant eukaryotes. However, intron content is highly variable among species. The tempo of intron gains and losses has been irregular and several factors may explain why some genomes are intron-poor whereas other are intron-rich.ResultsWe studied the dynamics of intron gains and losses in an α-amylase gene, whose product breaks down starch and other polysaccharides. It was transferred from an Actinobacterium to an ancestor of Agaricomycotina. This gene underwent further duplications in several species. The results indicate a high rate of intron insertions soon after the gene settled in the fungal genome. A number of these oldest introns, regularly scattered along the gene, remained conserved. Subsequent gains and losses were lineage dependent, with a majority of losses. Moreover, a few species exhibited a high number of both specific intron gains and losses in recent periods. There was little sequence conservation around insertion sites, then probably little information for splicing, whereas splicing sites, inside introns, showed typical and conserved patterns. There was little variation of intron size.ConclusionsSince most Basidiomycetes have intron-rich genomes and this richness was ancestral in Fungi, long before the transfer event, we suggest that the new gene was shaped to comply with requirements of the splicing machinery, such as short exon and intron sizes, in order to be correctly processed.

Spliceosomal introns in the 5′ untranslated region of plant BTL RING-H2 ubiquitin ligases are evolutionary conserved and required for gene expression

BackgroundIntrons located close to the 5′ end of a gene or in the 5′ untranslated region often exert positive effects on gene expression. This effect, known as intron-mediated enhancement (IME), has been observed in diverse eukaryotic organisms, including plants. The sequences involved in IME seem to be spread across the intron and function in an additive manner. The IMEter algorithm was developed to predict plant introns that may enhance gene expression. We have identified several plant members of the BTL class of E3s, which may have orthologs across eukaryotes, that contain a 5′UTR intron. The RING finger E3 ligases are key enzymes of the ubiquitination system that mediate the transfer of ubiquitin to substrates.ResultsIn this study, we retrieved BTL sequences from several angiosperm species and found that 5′UTR introns showing a strong IMEter score were predicted, suggesting that they may be conserved by lineage. Promoter-GUS fusion lines were used to confirm the IME effect of these 5′UTR introns on gene expression. IMEter scores of BTLs were compared with the 5′UTR introns of two gene families MHX and polyubiquitin genes.ConclusionsAnalysis performed in two Arabidopsis BTL E3 ligases genes indicated that the 5′UTR introns were essential for gene expression in all the tissues tested. Comparison of the average 5′UTR intron size on three gene families in ten angiosperm species suggests that a prevalent size for a 5′UTR intron is in the range of 600 nucleotides, and that the overall IMEter score within a gene family is preserved across several angiosperms. Our results indicated that gene expression dependent on a 5′UTR intron is an efficient regulatory mechanism in BTL E3 ligases that has been preserved throughout plant evolution.

Clinical utility gene card for: Rothmund–Thomson syndrome

1.5 Mutational spectrum Biallelic mutations in the RECQL4 gene are associated with Rothmund–Thomson (RTS) and two additional recessive disorders: RAPADILINO (RAdial hypoplasia, PAtellae hypoplasia and cleft or arched PAlate, DIarrhea and DIslocated joints, LIttle size and LImb malformation, slender NOse and NOrmal intelligence) and Baller– Gerold syndrome (BGS). More than 60 disease-causing mutations have been reported, of which at least 40 have been detected in RTS patients.1,2 The types of observed mutations are as follows, in order of decreasing prevalence: nonsense or frameshift mutations; splicing alterations, including substitutions at canonical splice junctions or at splice-site consensus sequences and subtle intronic deletions that reduce intron size below the threshold (o80 bp) required for correct splicing;3,4 and missense mutations. There are a few recurrent mutations, among which the most common, exon 9 c.1573delT (p.Cys525AlafsX33), has been detected in patients with all three RECQL4-associated diseases. This truncating mutation accounts for approximately one-third of RTS mutations and has only been found in compound heterozygous patients from multiple ethnic backgrounds. A deletion of the entire gene has never been identified.

Cloning and expression of two HSC70 genes in walking catfishClarias macrocephalus(Günther, 1864) challenged withAeromonas hydrophila

Two different cDNAs encoding heat shock cognate protein 70 genes were isolated from the liver of walking catfish. The complete cDNA sequences for walking catfish HSC70-1 and HSC70-2 were of identical length (2278 bp), with an open reading frame of 1950 bp and a predicted 649 amino acid protein. Genomic sequences of the walking catfish HSC70-1 and HSC70-2 genes were composed of eight exons and seven introns, but exon and intron sizes were different. Amino acid sequences of both proteins shared 94% similarity with 38 substitutions. The walking catfish HSC70-1 and HSC70-2 proteins shared 82–95% identity of amino acids with those of other teleosts, chicken and human. Phylogenetic analysis revealed that walking catfish HSC70-1 and HSC70-2 clustered with channel catfish HSC71. Under normal conditions, HSC70-1 and HSC70-2 transcripts were expressed at similar levels in liver, gills, brain and skeletal muscle of walking catfish. Bacterial infection by Aeromonas hydrophila differentially induced the expression of HSC70s in these tissues during 48 h. HSC70-1 was constitutively expressed at low levels in most tissues, whereas HSC70-2 was up-regulated at moderate and high levels in liver and skeletal muscle respectively. The significant increase in the expression of HSC70-2 in these tissues may relate to the role of HSC70-2 in the immune response of walking catfish.

Genomic organization, polymorphisms and molecular evolution of the goose-type lysozyme gene from Zhikong scallop Chlamys farreri

Lysozyme is a ubiquitous hydrolase that plays an important role in protecting the host against pathogenic infection. In the present study, the genomic DNA sequence of an invertebrate goose type (G type) lysozyme (designated CfLysG) was cloned from Zhikong scallop Chlamys farreri by genome walking technique. The full-length DNA of CfLysG gene was of 9455bp, and the fragment from the transcription site to the polyadenylation site was of 8217bp. A total of 104 SNPs and 29 ins–del polymorphisms were identified in the genomic sequence of CfLysG, and most of them were located in the promoter and intron regions, except for three SNPs located in the exon regions. Some putative transcription factor binding sites in the promoter region suggested the involvement of CfLysG in immune responses. There were six exons of 55, 60, 90, 113, 145 and 140bp interrupted by five relatively large introns in the genomic DNA sequence. CfLysG exhibited a unique exon–intron organization which was different from both its vertebrate and invertebrate homologues. Though some introns were lost in the urochordate homologues, four of the five introns in CfLysG DNA shared homologous positions with vertebrate G type lysozyme genes, which indicated the existence of these introns in the ancestry of G type lysozyme. With respect to the number and size of both exons and introns, the gene organization of CfLysG was more similar to that of vertebrate G type lysozyme, but its amino acid sequence shared higher similarity with that of other invertebrate G type lysozymes. In the phylogenic tree, G type lysozymes from mollusk were clustered together and formed a sister clade to the urochordate and vertebrate G type lysozymes. G type lysozyme was separated from C and I type lysozymes and closely matched to the phage-type lysozyme. The results suggested that G type lysozyme might be evolutionarily closest to the lysozyme ancestor, and they would be helpful to understand the evolution of lysozyme genes.