Characterization and mapping of expressed sequence tags isolated from a subtracted cDNA library of Litopenaeus vannamei injected with white spot syndrome virus

A large number of polymorphic simple sequence repeats (SSRs) or microsatellites are needed to develop a genetic map for shrimp. However, developing an SSR map is very time-consuming, expensive, and most SSRs are not specifically linked to gene loci of immediate interest. We report here on our strategy to develop polymorphic markers using expressed sequence tags (ESTs) by designing primers flanking single or multiple SSRs with three or more repeats. A subtracted cDNA library was prepared using RNA from specific pathogen-free (SPF) Litopenaeus vannamei juveniles (similar to 1 g) collected before (0) and after (48 h) inoculation with the China isolate of white spot syndrome virus (WSSV). A total of 224 clones were sequenced, 194 of which were useful for homology comparisons against annotated genes in NCBI nonredundant (nr) and protein databases, providing 179 sequences encoded by nuclear DNA, 4 mitochondrial DNA, and 11 were similar to portions of WSSV genome. The nuclear sequences clustered in 43 groups, 11 of which were homologous to various ESTs of unknown function, 4 had no homology to any sequence, and 28 showed similarities to known genes of invertebrates and vertebrates, representatives of cellular metabolic processes such as calcium ion balance, cytoskeleton mRNAs, and protein synthesis. A few sequences were homologous to immune system-related (allergens) genes and two were similar to motifs of the sex-lethal gene of Drosophila. A large number of EST sequences were similar to domains of the EF-hand superfamily (Ca2+ binding motif and FRQ protein domain of myosin light chains). Single or multiple SSRs with three or more repeats were found in approximately 61 % of the 179 nuclear sequences. Primer sets were designed from 28 sequences representing 19 known or putative genes and tested for polymorphism (EST-SSR marker) in a small test panel containing 16 individuals. Ten (53%) of the 19 putative or unknown function genes were polymorphic, 4 monomorphic, and 3 either failed to satisfactorily amplify genomic DNA or the allele amplification conditions need to be further optimized. Five polymorphic ESTs were genotyped with the entire reference mapping family, two of them (actin, accession #CX535973 and shrimp allergen arginine kinase, accession #CX535999) did not amplify with all offspring of the IRMF panel suggesting presence of null alleles, and three of them amplified in most of the IRM F offspring and were used for linkage analysis. EF-hand motif of myosin light chain (accession #CX535935) was placed in ShrimpMap's linkage group 7, whereas ribosomal protein S5 (accession #CX535957) and troponin I (accession #CX535976) remained unassigned. Results indicate that (a) a large number of ESTs isolated from this cDNA library are similar to cytoskeleton mRNAs and may reflect a normal pathway of the cellular response after im infection with WSSV, and (b) primers flanking single or multiple SSRs with three or more repeats from shrimp ESTs could be an efficient approach to develop polymorphic markers useful for linkage mapping. Work is underway to map additional SSR-containing ESTs from this and other cDNA libraries as a plausible strategy to increase marker density in ShrimpMap.

Some palaemonid shrimps (Crustacea : Decapoda) from northern South China Sea, with descripitions of two new species

Thirty-four palaemonid species from the northern area of South China Sea are reported in the present paper, including two new species, Paranchistus lini, new species, and Periclimenes chacei, new species, and eight other species are reported for the first time from South China Sea: Leander tenuicornis (Say, 1818), Leandrites stenopus Holthuis, 1950, Palaemon debilis Dana, 1852, Urocaridella antonbrurnii (Bruce, 1967), Kemponia platycheles (Holthuis, 1952), Paraclimenes franklini (Bruce, 1990b), Paranchistus spondylis Suzuki, 1971, Periclimenes hirsutus Bruce, 1971.

Two trichodinid ectoparasites from marine molluscs in the Yellow Sea, off China, with the description of Trichodina caecellae n. sp (Protozoa : Ciliophora : Peritrichia)

The morphology and infraciliature of two ectoparasitic ciliates, Trichodina caecellae n. sp. and T. ruditapicis Xu, Song & Warren, 2000, parasitising the gills of marine molluscs from the Shandong coast of the Yellow Sea, China, were investigated following wet silver nitrate and protargol impregnation. T. caecellae was found on the small marine sand clam Caecella chinensis Deshayes and is distinguished mainly by the acute triangle-like blade, the very delicate central part and the needle-shaped ray. T. ruditapicis was studied based on four populations from three clams: two populations from Ruditapes philippinarum (Adams) and one each from Saxidomus purpuratus (Sowerby) and Solen grandis Dunker. All four populations fell within the range of morphometry and agreed closely in the overall appearance of the adhesive disc. However, variability was found in the denticle structure, especially in populations from different host clams. Our observations suggest that denticle morphology may be more or less variable between and within populations, and that such minor differences should not be overestimated. It should be emphasised that, except for the denticle morphology, the bright granules or circles in the centre of the adhesive disc represent another important feature facilitating the identification of this trichodinid species.

东喜马拉雅构造结新生代碰撞变形与岩石圈演化

The Eastern Himalayan Syntaxis (EHS) is one of the strongest deformation area along the Himalayan belt resulted from the collision between Indian plate and the Eurasian Plate since the 50～60Ma, and has sensitivity tracked and preserved the whole collisional processes. It should depend on the detail geological investigations to establish the deformational accommodate mode, and the uplift history, to elucidate the deep structure and the crust-mantle interaction of the Tibet Plateau of the EHS. The deep-seated (Main Mantle Thrusts) structures were exhumed in the EHS. The MMT juxtapose the Gangdese metamorphic basement and some relic of Gangdese mantle on the high Himalayan crystalline series. The Namjagbawa group which is 1200～1500Ma dated by U/Pb age of zircon and the Namla group which is 550Ma dated by U/Pb age of zircon is belong to High Himalayan crystalline series and Gangdese basement respectively. There is some ophiolitic relic along the MMT, such as metamorphic ocean mantle peridotite and metamorphic tholeiite of the upper part of ocean-crust. The metamorphic ocean mantle peridotites (spinel-orthopyroxene peridotite) show U type REE patterns. The ~(87)Sr/~(86)Sr ratios were, 0.709314～0.720788, and the ~(143)Nd/~(144)Nd ratios were 0.512073～0.512395, plotting in the forth quadrant on the ~(87)Sr/~(86)Sr-~(143)Nd/~(144)Nd isotope diagram. Some metamorphic basalt (garnet amphibolite) enclosures have been found in the HP garnet-kynite granulite. The garnet amphibolites can be divided two groups, the first group is deplete of LREE, and the second group is flat or rich LREE, and their ~(87)Sr/~(86)Sr, ~(143)Nd/~(144)Nd ratios were 0.70563～0.705381 and 0.512468～0.51263 respectively. Trace element and isotopic characteristics of the garnet amphibolites display that they formed in the E-MORB environment. Some phlogolite amphibole harzburgites, which exhibit extensive replacement by Phl, Amp, Tc and Dol etc, were exhumed along the MMT. The Phl-Amp harzburgites are rich in LREE and LILE, such as Rb, K etc, and depletes Eu (Eu~* = 0.36 ～ 0.68) and HFSE, such as Nb, Ta, Zr, Hf, P, Ti etc. The trace element indicate that the Phl-Amp harzburgites have island arc signature. Their ~(87)Sr/~(86)Sr are varied from 0.708912 to 0.879839, ~(143)Nd/~(144)Nd from 0.511993 to 0.512164, ε Nd from- 9.2 to - 12.6. Rb/Sr isochrone age of the phlogolite amphibole harzburgite shows the metasomatism took place at 41Ma, and the Amp ~(40)Ar/~(39)Ar cooling age indcate the Phl-Amp harzburgite raising at 16Ma. There is an intense crust shortening resulted from the thrust faults and folds in the Cayu block which is shortened more 120km than that of the Lasha block in 35～90Ma. With the NE corner of the India plate squash into the Gangdese arc, the sinistral Pai shear fault and the dextral Aniqiao shear fault on the both sides of the Great bent of Yalun Zangbu river come into active in 21～26Ma. On the other hand, the right-lateral Gongrigabu strike-slip faults come into activity at the same period, a lower age bound for the Gongrigabu strike-slip fault is estimated to be 23～24Ma from zircon of ion-probe U/Pb thermochronology. The Gongrigabu strike-slip faults connect with the Lhari strike-slip fault in the northwestern direction and with the Saganing strike-slip at the southeastern direction. Another important structure in the EHS is the Gangdese detachment fault system (GDS) which occurs between the sedimental cover and the metamorphic basement. The lower age of the GDS is to be 16Ma from the preliminary 40Ar/39Ar thermochronology of white mica. The GDS is thought to be related to the reverse of the subducted Indian crust and the fast uplift of the EHS. Structural and thermochronology investigation of the EHS suggest that the eastern Tibet and the western Yunnan rotated clockwise around the EHS in the period of 35～60Ma. Later, the large-scale strike-slip faults (RRD, Gaoligong and Saganing fault) prolongate into the EHS, and connect with the Guyu fault and Gongrigabu fault, which suggest that the Indianchia block escape along these faults. Two kind of magmatic rocks in the EHS have been investigated, one is the mantle-derived amphibole gabbro, dioposide diorite and amphibole diorite, another is crust origin biotit-garnet adamellite, biotit-garnet granodiorite and garnet-amphibole-biotite granite. The amphibole gabbro dioposite diorite and amphibole diorite are rich in LREE, and LILE, such as Ba, Rb, Th, K, Sr etc, depleted in HFSE, such as Nb, Ta, Zr, Hf, Ti etc. The ratio of ~(87)Sr/~(86)Sr are from 0.7044 to 0.7048, ~(143)Nd/~(144)Nd are from 0.5126 to 0.5127. The age of the mantle origin magamatic rocks, which result from the partial melt of the raising and decompression anthenosphere, is 8Ma by ~(40)Ar/~(39)Ar dating of amphibole from the diorite. The later crust origin biotite-garnet adamellite, biotite-garnet granodiorite and garnet-amphibole-biotite granite are characterized by aboudance in LREE, and strong depletion of Eu. The ratios of ~(87)Sr-~(86)Sr are from 0.795035 to 0.812028, ~(143)Nd/~(144)Nd from 0.51187 to 0.511901. The ~(40)Ar/~(39)Ar plateau age of the amphibole from the garnet-amphibole-biotite granite is 17.5±0.3Ma, and the isochrone age is 16.8±0.6Ma. Their geochemical characteristics show that the crust-derived magmatic rocks formed from partial melting of the lower curst in the post-collisional environment. A group of high-pressure kaynite-garnet granulites and enclave of high-pressure garnet-clinopyroxene grnulites and calc-silicate grnulites are outcroped along the MMT. The peak metamorphic condition of the high-pressure granulites yields T=800～960 ℃, P=1.4～1.8Gpa, corresponding the condition of 60km depth. The retrograde assemblages of the high-pressure grnulites occur at the condition of T=772.3～803.3 ℃, P=0.63～0.64Gpa. The age of the peak metamorphic assemblages are 45 ～ 69Ma indicated by the zircon U/Pb ion-plobe thermochronology, and the retrograde assemblage ages are 13～26Ma by U/Pb, ~(40)Ar/~(39)Ar thermochronology. The ITD paths of the high-pressure granulites show that they were generated during the tectonic thickening and more rapid tectonic exhumation caused by the subducting of the Indian plate and subsequent break-off of the subducted slab. A great deal of apatite, zircon and sphene fission-track ages, isotopic thermochronology of the rocks in the EHS show that its rapid raising processes of the EHS can be divided into three main periods. There are 35～60Ma, 13～25Ma, 0～3Ma. 3Ma is a turn in the course of raising in the EHS which is characterized by abruptly acceleration of uplifting. The uplift ratios are lower than 1mm .a~(-1) before 3Ma, and higher than 1mm .a~(-1) with a maximum ratio of 30mm .a~(-1) since 3Ma. The bottom (knick point) of the partial anneal belt is 3.8km above sea level in the EHS, and correspond to age of 3Ma determined by fission-track age of apatite. The average uplift ratio is about 1.4 mm .a~(-1) below the knick point. The EHS has raised 4.3km from the surface of 2.36km above sea level since 3Ma estimated by the fossil partial anneal belt of the EHS. We propose a two-stage subduction model (B+A model) basing on Structural, thermochronological, magmatical, metamorphic and geophysical investigations of the EHS. The first stage is the subduction of the Indian continental margin following after the subduction of the Tethys Ocean crust and subsequent collision with the Gangdese arc, and the second stage is the Indian crust injecting into the lower crust and upper mantle of the Tibet plateau. Slab break-off seems to be occurred between these two stages.

小波域地震资料的处理方法（WSSP）研究

The modeling formula based on seismic wavelet can well simulate zero - phase wavelet and hybrid-phase wavelet, and approximate maximal - phase and minimal - phase wavelet in a certain sense. The modeling wavelet can be used as wavelet function after suitable modification item added to meet some conditions. On the basis of the modified Morlet wavelet, the derivative wavelet function has been derived. As a basic wavelet, it can be sued for high resolution frequency - division processing and instantaneous feature extraction, in acoordance with the signal expanding characters in time and scale domains by each wavelet structured. Finally, an application example proves the effectiveness and reasonability of the method. Based on the analysis of SVD (Singular Value Decomposition) filter, by taking wavelet as basic wavelet and combining SVD filter and wavelet transform, a new de - noising method, which is Based on multi - dimension and multi-space de - noising method, is proposed. The implementation of this method is discussed the detail. Theoretical analysis and modeling show that the method has strong capacity of de - noising and keeping attributes of effective wave. It is a good tool for de - noising when the S/N ratio is poor. To give prominence to high frequency information of reflection event of important layer and to take account of other frequency information under processing seismic data, it is difficult for deconvolution filter to realize this goal. A filter from Fourier Transform has some problems for realizing the goal. In this paper, a new method is put forward, that is a method of processing seismic data in frequency division from wavelet transform and reconstruction. In ordinary seismic processing methods for resolution improvement, deconvolution operator has poor part characteristics, thus influencing the operator frequency. In wavelet transform, wavelet function has very good part characteristics. Frequency - division data processing in wavelet transform also brings quite good high resolution data, but it needs more time than deconvolution method does. On the basis of frequency - division processing method in wavelet domain, a new technique is put forward, which involves 1) designing filter operators equivalent to deconvolution operator in time and frequency domains in wavelet transform, 2) obtaining derivative wavelet function that is suitable to high - resolution seismic data processing, and 3) processing high resolution seismic data by deconvolution method in time domain. In the method of producing some instantaneous characteristic signals by using Hilbert transform, Hilbert transform is very sensitive to high - frequency random noise. As a result, even though there exist weak high - frequency noises in seismic signals, the obtained instantaneous characteristics of seismic signals may be still submerged by the noises. One method for having instantaneous characteristics of seismic signals in wavelet domain is put forward, which obtains directly the instantaneous characteristics of seismic signals by taking the characteristics of both the real part (real signals, namely seismic signals) and the imaginary part (the Hilbert transfom of real signals) of wavelet transform. The method has the functions of frequency division and noise removal. What is more, the weak wave whose frequency is lower than that of high - frequency random noise is retained in the obtained instantaneous characteristics of seismic signals, and the weak wave may be seen in instantaneous characteristic sections (such as instantaneous frequency, instantaneous phase and instantaneous amplitude). Impedance inversion is one of tools in the description of oil reservoir. one of methods in impedance inversion is Generalized Linear Inversion. This method has higher precision of inversion. But, this method is sensitive to noise of seismic data, so that error results are got. The description of oil reservoir in researching important geological layer, in order to give prominence to geological characteristics of the important layer, not only high frequency impedance to research thin sand layer, but other frequency impedance are needed. It is difficult for some impedance inversion method to realize the goal. Wavelet transform is very good in denoising and processing in frequency division. Therefore, in the paper, a method of impedance inversion is put forward based on wavelet transform, that is impedance inversion in frequency division from wavelet transform and reconstruction. in this paper, based on wavelet transform, methods of time - frequency analysis is given. Fanally, methods above are in application on real oil field - Sansan oil field.