Vertical distribution of bacterial and archaeal communities along discrete layers of a deep-sea cold sediment sample at the East Pacific Rise (similar to 13 degrees N)

The community structure and vertical distribution of prokaryotes in a deep-sea (ca. 3,191 m) cold sediment sample (ca. 43 cm long) collected at the East Pacific Rise (EPR) similar to 13 degrees N were studied with 16SrDNA-based molecular analyses. Total community DNA was extracted from each of four discrete layers EPRDS-1, -2, -3 and -4 (from top to bottom) and 16S rDNA were amplified by PCR. Cluster analysis of DGGE profiles revealed that the bacterial communities shifted sharply between EPRDS-1 and EPRDS-2 in similarity coefficient at merely 49%. Twenty-three sequences retrieved from DGGE bands fell into 11 groups based on BLAST and bootstrap analysis. The dominant groups in the bacterial communities were Chloroflexi, Gamma proteobacteria, Actinobacterium and unidentified bacteria, with their corresponding percentages varying along discrete layers. Pairwise Fst (F-statistics) values between the archaeal clone libraries indicated that the archaeal communities changed distinctly between EPRDS-2 and EPRDS-3. Sequences from the archaeal libraries were divided to eight groups. Crenarchaea Marine Group I (MGI) was prevalent in EPRDS-1 at 83%, while Uncultured Crenarchaea group II B (UCII B) abounded in EPRDS-4 at 61%. Our results revealed that the vertically stratified distribution of prokaryotic communities might be in response to the geochemical settings and suggested that the sampling area was influenced by hydrothermalism. The copresence of members related to hydrothermalism and cold deep-sea environments in the microbial community indicated that the area might be a transitional region from hydrothermal vents to cold deep-sea sediments.

Molecular cloning and expression of a novel Kazal-type serine proteinase inhibitor gene from Zhikong scallop Chlamys farreri, and the inhibitory activity of its recombinant domain

Serine proteinase inhibitors (SPIs) play important roles in host physiological and immunological processes in all multicellular organisms. A novel Kazal-type SPI gene was cloned from the Zhikong scallop Chlamys farreri (designated as CfKZSPI) by expressed sequence tag (EST) and rapid amplification of cDNA ends (RACE) approaches. The full-length cDNA of CfKZSPI was of 1788 nucleotides with a canonical polyadenylation signal sequence AATAAA and a polyA tail, and an open reading frame (ORF) encoding a polypeptide of 509 amino acids with a putative signal peptide of 22 amino acids. The deduced amino acid sequence of CfKZSPI contained 12 tandem Kazal domains with high similarity to other Kazal-type SPIs. The temporal expression of CfKZSPI in hemocytes after Vibrio anguillorum challenge was recorded by quantitative real-time RT-PCR. The relative mRNA expression level of CfKZSPI was up-regulated and reached 43.6-fold at 3 h post-challenge. After a decrease at 6 h, the expression Level increased again and reached 207.8-fold at 12 h post-challenge. The 12th Kazal domain of CfKZSPI was recombined into pET-32a(+) and expressed in Escherichia coli Rosetta-gami (DE3) to investigate its inhibitory activity. The purified recombinant protein (rCf KZSPI-1 2) showed significant inhibitory activity against trypsin but no activity against thrombin. When the molar ratio of inhibitor to trypsin reached 1:1, almost 90% of the enzyme activity could be inhibited, which suggested that one molecule of rCfKZSPI-12 was able to inhibit one molecule of trypsin. Kinetics analysis with Dixon plot showed that the inhibition constant (K-i) of rCfKZSPI-12 to trypsin was 173 nmol L-1. These results indicated that CfKZSPI was a novel Kazal-type SPI with significant inhibitory activity against trypsin, and was suspected to be involved in scallop immune response. (c) 2008 Elsevier Ltd. All rights reserved.

Lineage-specific domain fusion in the evolution of purine nucleotide cyclases in cyanobacteria

Cyclic nucleotides (both cAMP and cGMP) play extremely important roles in cyanobacteria, such as regulating heterocyst formation, respiration, or gliding. Catalyzing the formation of cAMP and cGMP from ATP and GTP is a group of functionally important enzymes named adenylate cyclases and guanylate cyclases, respectively. To understand their evolutionary patterns, in this study, we presented a systematic analysis of all the cyclases in cyanobacterial genomes. We found that different cyanobacteria had various numbers of cyclases in view of their remarkable diversities in genome size and physiology. Most of these cyclases exhibited distinct domain architectures, which implies the versatile functions of cyanobacterial cyclases. Mapping the whole set of cyclase domain architectures from diverse prokaryotic organisms to their phylogenetic tree and detailed phylogenetic analysis of cyclase catalytic domains revealed that lineage-specific domain recruitment appeared to be the most prevailing pattern contributing to the great variability of cyanobacterial cyclase domain architectures. However, other scenarios, such as gene duplication, also occurred during the evolution of cyanobacterial cyclases. Sequence divergence seemed to contribute to the origin of putative guanylate cyclases which were found only in cyanobacteria. In conclusion, the comprehensive survey of cyclases in cyanobacteria provides novel insight into their potential evolutionary mechanisms and further functional implications.

A novel C1q-domain-containing protein from Zhikong scallop Chlamys farreri with lipopolysaccharide binding activity

The C1q-domain-containing (C1qDC) proteins are a family of proteins characterized by a globular C1q (gC1q) domain in their C-terminus. They are involved in various processes of vertebrates and supposed to be an important pattern recognition receptor in innate immunity of invertebrates. In this study, a novel member of C1q-domain-containing protein family was identified from Zhikong scallop Chlamys farreri (designated as CfC1qDC) by expressed sequence tag (EST) and rapid amplification of cDNA ends (RACE) approaches. The full-length cDNA of CfC1qDC was of 777 bp, consisting of a T-terminal untranslated region (UTR) of 62 bp and a 3' UTR of 178 bp with a polyadenylation signal sequence AATAAA and a poly (A) tail. The CfC1qDC cDNA encoded a polypeptide of 178 amino acids, including a signal peptide and a C1q-domain of 158 amino acids with the theoretical isoelectric point of 5.19 and the predicted molecular weight of 17.2 kDa. The C1q-domain in CfC1qDC exhibited homology with those in sialic acid binding lectin from mollusks and C1qDC proteins from higher vertebrates. The typical 10 beta-strand jelly-roll folding topology structure of C1q-domain and the residues essential for effective packing of the hydrophobic core were well conserved in CfC1qDC. By fluorescent quantitative real-time PCR, mRNA transcripts of CfC1qDC were mainly detected in kidney, mantle, adductor muscle and gill, and also marginally detectable in hemocytes. In the bacterial challenge experiment, after the scallops were challenged by Listonella anguillarum, there was a significant up-regulation in the relative expression level of CfC1qDC and at 6 h post-injection, the mRNA expression reached the maximum level and was 4.55-fold higher than that of control scallops. Similarly, the expression of CfC1qDC mRNA in mixed primary cultures of hemocytes stimulated by lipopolysaccharides (LPS) was up-regulated and reached the maximum level at 6 h post-stimulation, and then dropped back to the original level gradually. In order to investigate its function, the cDNA fragment encoding the mature peptide of CfC1qDC was recombined and expressed in Escherichia coli BL21 (DE3). The recombinant CfC1qDC protein displayed a significantly strong activity to bind LIDS from E. coli, although no obvious antibacterial or agglutinating activity toward Gram-negative bacteria E. coli JM109, L. anguillarum and Gram-positive bacteria Micrococcus luteus was observed. These results suggested that CfC1qDC was absolutely a novel member of the C1qDC protein family and was involved in the recognition of invading microorganisms probably as a pattern recognition molecule in mollusk. (c) 2008 Elsevier Ltd. All rights reserved.

A novel serine protease with clip domain from scallop Chlamys farreri

The serine proteases with clip domain are involved in various innate immune functions in invertebrate such as antimicrobial activity, cell adhesion, pattern recognition and regulation of the prophenoloxidase system. A serine protease with clip-domain cDNA (Cf SP) was obtained by Expressed sequence taggings (ESTs) method and rapid amplification of cDNA ends (RACE). The Cf SP full-length cDNA was of 1,152 bp, including a 5'-terminal untranslated region (UTR) of 63 bp, a 3'-terminal UTR of 81 bp with a canonical polyadenylation signal sequence AATAAA and a poly(A) tail, and an open reading frame of 1,008 bp encoding a polypeptide of 336 amino acids with a putative signal peptide of 19 amino acids. The deduced amino acid sequence of Cf SP contained an amino-terminal clip domain with three disulfide bonds formed six conserved Cys residues, a carboxyl-terminal trypsin-like domain with the conserved His-Asp-Ser catalytic triad, and a low complexity linker sequence. The Cf SP was strongly expressed in hemocytes and the mRNA expression of Cf SP was up-regulated and increased 3.2-fold and 2.6-fold at 16 h after injection of Vibrio anguillarum and Micrococcus luteus. The results suggested that Cf SP gene might be involved in immune response of Gram-negative and Gram-positive microbial infection in scallop.

Molecular cloning, characterization and expression of a serine protease with clip-domain homologue from scallop Chlamys farreri

Serine proteases play critical roles in a variety of invertebrate immune defense responses, including hemolymph coagulation, antimicrobial peptide synthesis, and melanization. The first mollusk serine protease with clip-domain (designated CFSP1) cDNA was obtained from the scallop Chlamys farreri challenged with Vibrio anguillarum by randomly sequencing a whole tissue cDNA library and rapid amplification of cDNA ends (RACE). The full-length cDNA of the C. farreri serine protease was 1211 bp, consisting of a 5-terminal untranslated region (UTR) of 72 bp, a 3'-terminal UTR of 77 bp with a canonical polyadenylation signal sequence AATAAA and a poly (A) tail, and an open reading frame of 1062 bp. The CFSP1 cDNA encoded a polypeptide of 354 amino acids with a putative signal peptide of 19 amino acids and a mature protein of 335 amino acids. The deduced amino acid sequence of CFSP1 contained an amino-terminal clip domain, a low complexity region, and a carboxyl-terminal serine protease domain. CFSP1 mRNA was mainly expressed constitutively in the hemocytes and was up-regulated and increased 2.9- and 1.9-fold at 16 h after injury and injection of bacteria. (c) 2006 Elsevier Ltd. All rights reserved.

Dissection and localization of the immunostimulating domain of Edwardsiella tarda FliC

Bacterial flagellin is known to induce potent immune response in vertebrate systems via the toll-like receptor (TLR) 5. As a result, flagellin has been studied extensively as a vaccine adjuvant. In a previous study, we examined the vaccine and adjuvant potentials of the flagellin (FliC) of the fish pathogen Edwardsiella tarda. We found that E. tarda FliC induced low protective immunity by itself but could function as a molecular adjuvant and potentiate the specific immune response induced by the E. tarda antigen Eta6. Since FliC is a large protein and organized into distinct structural domains, we wondered whether the immunostimulating effect observed with the full-length protein could be localized to a certain region. To investigate this question, we in the present study dissected the FliC protein into several segments according to its structural features: (i) N163, which consists of the conserved N-terminal 163 residues of FliC; (ii) M160, which consists of the variable middle 160 residues; (iii) C94, which consists of the conserved C-terminal 94 residues; (iv) NC257, which is an artificial fusion of N163 and C94. To examine the adjuvanticity of the FliC fragments, DNA vaccine plasmids expressing FliC fragments in fusion with Eta6 were constructed and used to immunize Japanese flounder. The results showed that N163 produced the best adjuvant effect, which, in respect to improvement in the relative percent survival of the vaccinated fish, was comparable to that of the full-length FliC. None of the other FliC fragments exhibited apparent immunopotentiating effect. Further analysis showed that N163 enhanced the production of serum specific antibodies and, like full-length FliC, significantly upregulated the expression of the genes that are possibly involved in innate and adaptive immunity. These results indicate that N163 is the immunodominant region of FliC and suggest that E. tarda FliC may induce immune responses in Japanese flounder via mechanisms alternative to that involving TLR5. (C) 2010 Elsevier Ltd. All rights reserved.

AiC1qDC-1, a novel gC1q-domain-containing protein from bay scallop Argopecten irradians with fungi agglutinating activity

The globular C1q-domain-containing (C1qDC) proteins are a family of versatile pattern recognition receptors via their globular C1q (gC1q) domain to bind various ligands including several PAMPs on pathogens. In this study, a new gC1q-domain-containing protein (AiC1qDC-1) gene was cloned from Argopecten irradians by rapid amplification of cDNA ends (RACE) approaches and expressed sequence tag (EST) analysis. The full-length cDNA of AiC1qDC-1 was composed of 733 bp, encoding a signal peptide of 19 residues and a typical gC1q domain of 137 residues containing all eight invariant amino acids in human C1qDC proteins and seven aromatic residues essential for effective packing of the hydrophobic core of AiC1qDC-1. The gC1q domain of AiC1qDC-1, which possessed the typical 10-stranded beta-sandwich fold with a jelly-roll topology common to all C1q family members, showed high homology not only to those of Cl qDC proteins in mollusk but also to those of C1qDC proteins in human. The AiC1qDC-1 transcripts were mainly detected in the tissue of hepatopancreas and also marginally detectable in adductor, heart, mantle, gill and hemocytes by fluorescent quantitative real-time PCR. In the microbial challenge experiment, there was a significant up-regulation in the relative expression level of AiC1qDC-1 in hepatopancreas and hemocytes of the scallops challenged by fungi Pichia pastoris GS115, Gram-positive bacteria Micrococcus luteus and Gram-negative bacteria Listonella anguillarum. The recombinant AiC1qDC-1 (rAiC1qDC-1) protein displayed no obvious agglutination against M. luteus and L. anguillarum, but it aggregated P. pastoris remarkably. This agglutination could be inhibited by D-mannose and PGN but not by LPS, glucan or D-galactose. These results indicated that AiC1qDC-1 functioned as a pattern recognition receptor in the immune defense of scallops against pathogens and provided clues for illuminating the evolution of the complement classical pathway. (C) 2010 Elsevier Ltd. All rights reserved.

SmC2P1, a C2 domain protein from Scophthalmus maximus that binds bacterial pathogen-infected lymphocytes and reduces bacterial survival

C2 domains are protein structural modules found in many eukaryotic proteins involved in signal transduction, membrane trafficking, and immune defense. Most of the studied C2 domain-containing proteins are multi-domained in structure, in which the C2 domain is an independently folded motif and plays an essential role in calcium-dependent membrane-targeting. Although C2 domains isolated from intact proteins have been studied for biological functions, no study on natural proteins containing C2 domain only has been documented. In this study, we identified a Scophthalmus maximus protein SmC2P1 that is comprised of a single C2 domain and lacks any other apparent domain structures. The deduced amino acid sequence of SmC2P1 contains 129 residues and shares 36-38% identities with the C2 domains of the perforins of several fish species. Like typical C2 domains, SmC2P1 is predicted to organize into eight beta-strands with a Ca2+-binding site located in inter-strand loops. SmC2P1 expression was detected, in deceasing order, in liver, spleen, blood, brain, muscle, kidney, gill, and heart. Experimental challenge of turbot with a bacterial pathogen significantly upregulated SmC2P1 expression in kidney in a time-dependent manner. Recombinant SmC2P1 purified from yeast exhibits no hemolytic activity but binds to pathogen-infected kidney lymphocytes in the presence of calcium. Furthermore, interaction of recombinant SmC2P1 with bacterium-infected lymphocytes reduced bacterial survival. These results indicate that SmC2P1 is a functional protein that is involved in host immune defense against bacterial infection. (C) 2010 Elsevier Ltd. All rights reserved.

Identification, expression, function and localization of a DUF985 domain-containing hypothetical gene from amphioxus Branchiostoma belcheri

The progress in genome sequencing has led to an increasing submission of uncharacterized hypothetical genes with the domain of unknown function, DUF985, in GenBank, and none of these genes is related to a known protein. We therefore underwent an experimental study to identify the function of a DUF985 domain-containing hypothetical gene BbDUF985 (GenBank Accession No. AY273818) isolated from amphioxus Branchiostoma belcheri (B. belcheri). BbDUF985 was successfully expressed in both prokaryotic and eukaryotic systems, and its recombinant proteins expressed in both systems definitely exhibited an activity of phosphoglucose isomerase (PGI). Both tissue-section in situ hybridization and immunohistochemistry demonstrated that BbDUF985 was expressed in a tissue-specific manner, with most abundant levels in the hepatic caecum and ovary. In CHO cells transfected with the expression plasmid pEGFP-N1/BbDUF985, the fusion protein was targeted in the cytoplasm of CHO cells, suggesting that BbDUF985 is a cytosolic protein. In contrast, Western blotting indicated that BbDUF985 was also present in amphioxus humoral fluids, suggesting that it exists as a secreted protein as well. Our study provided a framework for further understanding the biochemical properties and physiological function of DUF985-containing hypothetical proteins in other species. (c) 2008 Elsevier Inc. All rights reserved.

jλ变换与离散振荡指标Md在海洋机器人控制系统动态综合中的应用

本文给出由连续系统的S域,到离散领域?的转换的通式。并给出了离散传递函数W(?)的基本属性、离散振荡指标Md的基本概念、?变换的稳定判据、二阶无静差数字随动系统的动态综合、Md综合在海洋机器人航向控制中的应用。

Fourier有限差分法深度偏移与正演模拟

With the development of seismic exploration, the target becomes more and more complex, which leads to a higher demand for the accuracy and efficiency in 3D exploration. Fourier finite-difference (FFD) method is one of the most valuable methods in complex structure exploration, which keeps the ability of finite-differenc method in dealing with laterally varing media and inherits the predominance of the phase-screen method in stablility and efficiency. In this thesis, the accuracy of the FFD operator is highly improved by using simulated annealing algorithm. This method takes the extrapolation step and band width into account, which is more suitable to various band width and discrete scale than the commonely-used optimized method based on velocity contrast alone. In this thesis, the FFD method is extended to viscoacoustic modeling. Based on one-way wave equation, the presented method is implemented in frequency domain; thus, it is more efficient than two-way methods, and is more convenient than time domain methods in handling attenuation and dispersion effects. The proposed method can handle large velocity contrast and has a high efficiency, which is helpful to further research on earth absorption and seismic resolution. Starting from the frequency dispersion of the acoustic VTI wave equation, this thesis extends the FFD migration method to the acoustic VTI media. Compared with the convetional FFD method, the presented method has a similar computational efficiency, and keeps the abilities of dealing with large velocity contrasts and steep dips. The numerical experiments based on the SEG salt model show that the presented method is a practical migration method for complex acoustical VTI media, because it can handle both large velocity contrasts and large anisotropy variations, and its accuracy is relatively high even in strong anisotropic media. In 3D case, the two-way splitting technique of FFD operator causes artificial azimuthal anisotropy. These artifacts become apparent with increasing dip angles and velocity contrasts, which prevent the application of the FFD method in 3D complex media. The current methods proposed to reduce the azimuthal anisotropy significantly increase the computational cost. In this thesis, the alternating-direction-implicit plus interpolation scheme is incorporated into the 3D FFD method to reduce the azimuthal anisotropy. By subtly utilizing the Fourier based scheme of the FFD method, the improved fast algorithm takes approximately no extra computation time. The resulting operator keeps both the accuracy and the efficiency of the FFD method, which is helpful to the inhancements of both the accuracy and the efficiency for prestack depth migration. The general comparison is presented between the FFD operator and the generalized-screen operator, which is valuable to choose the suitable method in practice. The percentage relative error curves and migration impulse responses show that the generalized-screen operator is much sensiutive to the velocity contrasts than the FFD operator. The FFD operator can handle various velocity contrasts, while the generalized-screen operator can only handle some range of the velocity contrasts. Both in large and weak velocity contrasts, the higher order term of the generalized-screen operator has little effect on improving accuracy. The FFD operator is more suitable to large velocity contrasts, while the generalized-screen operator is more suitable to middle velocity contrasts. Both the one-way implicit finite-difference migration and the two-way explicit finite-differenc modeling have been implemented, and then they are compared with the corresponding FFD methods respectively. This work gives a reference to the choosen of proper method. The FFD migration is illustrated to be more attractive in accuracy, efficiency and frequency dispertion than the widely-used implicit finite-difference migration. The FFD modeling can handle relatively coarse grids than the commonly-used explicit finite-differenc modeling, thus it is much faster in 3D modeling, especially for large-scale complex media.

基于匹配追踪方法的地震信号时频分析研究

Seismic signal is a typical non-stationary signal, whose frequency is continuously changing with time and is determined by the bandwidth of seismic source and the absorption characteristic of the media underground. The most interesting target of seismic signal’s processing and explaining is to know about the local frequency’s abrupt changing with the time, since this kind of abrupt changing is indicating the changing of the physical attributes of the media underground. As to the seismic signal’s instantaneous attributes taken from time-frequency domain, the key target is to search a effective, non-negative and fast algorithm time-frequency distribution, and transform the seismic signal into this time-frequency domain to get its instantaneous power spectrum density, and then use the process of weighted adding and average etc. to get the instantaneous attributes of seismic signal. Time-frequency analysis as a powerful tool to deal with time variant non-stationary signal is becoming a hot researching spot of modern signal processing, and also is an important method to make seismic signal’s attributes analysis. This kind of method provides joint distribution message about time domain and frequency domain, and it clearly plots the correlation of signal’s frequency changing with the time. The spectrum decomposition technique makes seismic signal’s resolving rate reach its theoretical level, and by the method of all frequency scanning and imaging the three dimensional seismic data in frequency domain, it improves and promotes the resolving abilities of seismic signal vs. geological abnormal objects. Matching pursuits method is an important way to realize signal’s self-adaptive decomposition. Its main thought is that any signal can be expressed by a series of time-frequency atoms’ linear composition. By decomposition the signal within an over completed library, the time-frequency atoms which stand for the signal itself are selected neatly and self-adaptively according to the signal’s characteristics. This method has excellent sparse decomposition characteristics, and is widely used in signal de-noising, signal coding and pattern recognizing processing and is also adaptive to seismic signal’s decomposition and attributes analysis. This paper takes matching pursuits method as the key research object. As introducing the principle and implementation techniques of matching pursuits method systematically, it researches deeply the pivotal problems of atom type’s selection, the atom dictionary’s discrete, and the most matching atom’s searching algorithm, and at the same time, applying this matching pursuits method into seismic signal’s processing by picking-up correlative instantaneous messages from time-frequency analysis and spectrum decomposition to the seismic signal. Based on the research of the theory and its correlative model examination of the adaptively signal decomposition with matching pursuit method, this paper proposes a fast optimal matching time-frequency atom’s searching algorithm aimed at seismic signal’s decomposition by frequency-dominated pursuit method and this makes the MP method pertinence to seismic signal’s processing. Upon the research of optimal Gabor atom’s fast searching and matching algorithm, this paper proposes global optimal searching method using Simulated Annealing Algorithm, Genetic Algorithm and composed Simulated Annealing and Genetic Algorithm, so as to provide another way to implement fast matching pursuit method. At the same time, aimed at the characteristics of seismic signal, this paper proposes a fast matching atom’s searching algorithm by means of designating the max energy points of complex seismic signal, searching for the most optimal atom in the neighbor area of these points according to its instantaneous frequency and instantaneous phase, and this promotes the calculating efficiency of seismic signal’s matching pursuit algorithm. According to these methods proposed above, this paper implements them by programmed calculation, compares them with some open algorithm and proves this paper’s conclusions. It also testifies the active results of various methods by the processing of actual signals. The problems need to be solved further and the aftertime researching targets are as follows: continuously seeking for more efficient fast matching pursuit algorithm and expanding its application range, and also study the actual usage of matching pursuit method.

地震波场模拟及叠前深度偏移的李群方法

The processes of seismic wave propagation in phase space and one way wave extrapolation in frequency-space domain, if without dissipation, are essentially transformation under the action of one parameter Lie groups. Consequently, the numerical calculation methods of the propagation ought to be Lie group transformation too, which is known as Lie group method. After a fruitful study on the fast methods in matrix inversion, some of the Lie group methods in seismic numerical modeling and depth migration are presented here. Firstly the Lie group description and method of seismic wave propagation in phase space is proposed, which is, in other words, symplectic group description and method for seismic wave propagation, since symplectic group is a Lie subgroup and symplectic method is a special Lie group method. Under the frame of Hamiltonian, the propagation of seismic wave is a symplectic group transformation with one parameter and consequently, the numerical calculation methods of the propagation ought to be symplectic method. After discrete the wave field in time and phase space, many explicit, implicit and leap-frog symplectic schemes are deduced for numerical modeling. Compared to symplectic schemes, Finite difference (FD) method is an approximate of symplectic method. Consequently, explicit, implicit and leap-frog symplectic schemes and FD method are applied in the same conditions to get a wave field in constant velocity model, a synthetic model and Marmousi model. The result illustrates the potential power of the symplectic methods. As an application, symplectic method is employed to give synthetic seismic record of Qinghai foothills model. Another application is the development of Ray+symplectic reverse-time migration method. To make a reasonable balance between the computational efficiency and accuracy, we combine the multi-valued wave field & Green function algorithm with symplectic reverse time migration and thus develop a new ray+wave equation prestack depth migration method. Marmousi model data and Qinghai foothills model data are processed here. The result shows that our method is a better alternative to ray migration for complex structure imaging. Similarly, the extrapolation of one way wave in frequency-space domain is a Lie group transformation with one parameter Z and consequently, the numerical calculation methods of the extrapolation ought to be Lie group methods. After discrete the wave field in depth and space, the Lie group transformation has the form of matrix exponential and each approximation of it gives a Lie group algorithm. Though Pade symmetrical series approximation of matrix exponential gives a extrapolation method which is traditionally regarded as implicit FD migration, it benefits the theoretic and applying study of seismic imaging for it represent the depth extrapolation and migration method in a entirely different way. While, the technique of coordinates of second kind for the approximation of the matrix exponential begins a new way to develop migration operator. The inversion of matrix plays a vital role in the numerical migration method given by Pade symmetrical series approximation. The matrix has a Toepelitz structure with a helical boundary condition and is easy to inverse with LU decomposition. A efficient LU decomposition method is spectral factorization. That is, after the minimum phase correlative function of each array of matrix had be given by a spectral factorization method, all of the functions are arranged in a position according to its former location to get a lower triangular matrix. The major merit of LU decomposition with spectral factorization (SF Decomposition) is its efficiency in dealing with a large number of matrixes. After the setup of a table of the spectral factorization results of each array of matrix, the SF decomposition can give the lower triangular matrix by reading the table. However, the relationship among arrays is ignored in this method, which brings errors in decomposition method. Especially for numerical calculation in complex model, the errors is fatal. Direct elimination method can give the exact LU decomposition But even it is simplified in our case, the large number of decomposition cost unendurable computer time. A hybrid method is proposed here, which combines spectral factorization with direct elimination. Its decomposition errors is 10 times little than that of spectral factorization, and its decomposition speed is quite faster than that of direct elimination, especially in dealing with a large number of matrix. With the hybrid method, the 3D implicit migration can be expected to apply on real seismic data. Finally, the impulse response of 3D implicit migration operator is presented.