角蟾科（Megophryidae; Anura）蝌蚪的形态多样性和生态适应

角蟾科（Megophryidae）是以角蟾属（Megophrys Kuhl and Van Hasselt, 1822）为模式属而建立的，隶于无尾目（Anura），变凹型亚目（Anomocoela）。角蟾科包括2 亚科11 属142 种，分布于东洋界，从巴基斯坦、中国西部向东直到菲律宾和苏达群岛；中国有9 属75 种分布于华中和华南地区。角蟾科被认为是原始的两栖动物之一，其分类学、系统学、生态学、动物地理学的研究均深受中外科学家的瞩目。近年来，通过形态学、古生物学、细胞学、生态学、支序系统学的研究，角蟾科的分类与系统学研究取得了较大进展。与成体形态和分子系统学研究结果相比较，蝌蚪的研究存在更多的问题和挑战，尚需深入研究：（1）角蟾科蝌蚪的形态多样性分析；（2）角蟾科的系统发育关系与蝌蚪的演化，以及口漏斗的起源；（3）角蟾科蝌蚪表型分化与栖息环境和觅食行为的适应演化。针对上述问题，本文对角蟾科9 属30 种蝌蚪的形态特征，包括外部宏观形态和口器外部结构特征、口器内部显微结构、唇齿和角质颌的亚显微结构作了深入细致、多层次的比较研究；通过12s rRNA 和cytochrome b 基因构建最大简约树，采用贝叶斯系统发育进行分析，蝌蚪型的演化采用祖先性状的重建方法分析；得到如下结论：1）初步将角蟾科蝌蚪分为4 种类型；并且建立了2 种新的角蟾科蝌蚪类型。A 型：拟髭蟾型蝌蚪，该型蝌蚪包括拟髭蟾属、髭蟾属、齿蟾属和齿突蟾属的物种；B 型：新类型，掌突蟾型蝌蚪，该型蝌蚪在本文中包括掌突蟾属、小臂蟾属的物种；C 型：新类型，短腿蟾型蝌蚪，一种特化类型，该型蝌蚪在本文中仅包括短腿蟾属的物种；D 型：角蟾型蝌蚪，该型蝌蚪在本文中包括无耳蟾属、小口拟角蟾属和异角蟾属的物种。2）对角蟾科的分类进行了修订：（1）支持角蟾科两个亚科的分类系统；（2）角蟾亚科包括拟角蟾属、异角蟾属、无耳蟾属和短腿蟾属；该亚科形态差异小，系统学关系比较复杂，暂不作族级分类的再划分；（3）拟髭蟾亚科分为2 个族：拟髭蟾族，该族物种具有类型A 的蝌蚪，包括4 个属：拟髭蟾属、髭蟾属、齿蟾属、齿突蟾属；掌突蟾族，该族物种具有类型B 的蝌蚪，包括2 个属：掌突蟾属和小臂蟾属。3）结合分子系统进化关系探讨了4 种蝌蚪类型的演化。（1）角蟾科蝌蚪的最近共同祖先来自于一类具有拟髭蟾型蝌蚪性状的蝌蚪；（2）掌突蟾型蝌蚪和角蟾亚科的蝌蚪是由具有拟髭蟾型蝌蚪性状的祖先蝌蚪分别演化而来；（3）短腿蟾型蝌蚪是角蟾型蝌蚪的一种特化类型；（4）外群蝌蚪具有与拟髭蟾型蝌蚪相似的性状，进一步印证了类拟髭蟾型蝌蚪是角蟾科蝌蚪的最近共同祖先的假说；（5）具有口漏斗的蝌蚪类型是由不具口漏斗的蝌蚪类型演化而来，在角蟾科中口漏斗是一种衍生性状。4）分析了角蟾科四种蝌蚪类型与栖息环境的适应演化。（1）角蟾科蝌蚪的口部和体形的变化反映了该科蝌蚪由缓流向类似静水生境的回水凼的渐变式适应，角蟾科蝌蚪的形态显示了多方面的适应变化；（2）随着蝌蚪类型由A 向D的演化，当水速较大时，拟髭蟾型的蝌蚪营流水攀吸型生活方式；当水速递减时，掌突蟾型蝌蚪营流水附着型生活方式；当水速进一步递减时，具有较小口漏斗的短腿蟾型蝌蚪和具有大漏斗的角蟾型蝌蚪营流水浮泳型生活。角蟾科蝌蚪对于水流递减的适应演化说明蝌蚪的生态学适应是具有进化意义的；（3）蝌蚪口器内部结构的分化揭示了蝌蚪和食性的适应关系，蝌蚪以口部的唇齿与角质颌刮取或吞吸水中的物质，然后，通过口乳突有选择地过滤进入口腔中食物。拟髭蟾亚科蝌蚪的唇齿多而窄，唇齿间距宽，颌鞘粗而稀，反映了其植食性为主的特点；它们的舌前乳突一般为指状，在口腔入口处所占面积小，其机械过滤的作用很多被唇齿和角质颌分担了；而角蟾亚科的蝌蚪，其角质颌弱，其舌前乳突一般为匙状，几乎填满了口腔入口处，因此舌前乳突起了主要的机械过滤作用。The family Megophryidae is the largest and most diverse families inArchaeobatrachia, and most of its species occur in India, Pakistan, and eastward intoChina, Southeast Asia, Borneo and the Philippines to the Sunda Islands. Currently thefamily includes 142 species have been grouped into two subfamilies, Megophryinaeand Leptobrachiinae. The mountains of central and southern China are rich in speciesof Megophryidae, 75 species belong to 9 genera and two subfamilies.The family was supposed to be ideal materials of studies in many fields of biology,such as taxonomy, evolution, systematics, ecology, and biogeography. Recently, therehave a great development in taxonomy and systematics of megophryids throughstudied by morphology, paleontology, cytology, ecology, and cladistics. However,larvae of megophryids were generally unknown, although the tadpoles might be veryimportant for above studies.In this paper, we examined the evolutionary scenario of the tadpoles’ morphologyin the context of a phylogenetic framework. Our objectives are (1) to evaluate thedivergence of larval body shape and oral discs in the family Megophryidae, (2) toexplore the evolutionary trends of the larvae in megophryidae, and test if thefunnel-shaped oral disc is apomorphic, and (3) to explore the relationship of the larvalstructure, diet and microhabitat.We examined larval morphology of 30 megophryid species, the larval body shape,oral discs, the buccopharyngeal cavity, and jaw sheaths and denticles of the Chinesemegophryid frogs were re-examined. We constructed a phylogeny of the species on thebasis of published mitochondrial cytochrome b and 16S rRNA gene segments usingpartitioned Bayesian analyses. Furthermore, hypothetical changes of larval morphologywere inferred using parsimony principle on the phylogeny. The results showed that:1) Four tadpole types in Megophryidae. The larval morphological charactersseries in Chinese megophryids fall into four general categories according to the bodyshape and oral discs: (A) Leptobrachiini type, species from genera Leptobrachium,Oreolalax, Scutiger and, Vibrissaphora share this type of tadpoles. (B) Leptolalax type,species of genus Leptolalax have this type of tadpoles. (C) Brachytarsophrys type,species of the genus Brachytarsophrys have this type of tadpoles. (D) Megophryinitype, species of the genera Atympanophrys, Ophryophryne, and Xenophrys share this type of tadpoles. Of which B and C are two novel types.2）Taxonomic implications. The present study leads us to reconsider the generalclassification of tribes attributed to members of Megophryidae. More specifically,concerning the phylogenetic relationships and the two novel tadpole types describedherein, we propose a provisional taxonomy for the family but suggest that further taxasampling of other megophryids be performed to confirm this taxonomic change. TheMegophryidae is composed of two subfamilies (Leptobrachiinae and Megophryinae).The Leptobrachiinae was recogonized the two tribes: (1) tribe Leptobrachiini sensuDubois, corresponding to the tadpole of type A, including four genera, i.e.,Leptobrachium, Oreolalax, Scutiger and, Vibrissaphora; (2) tribe Leptolalaxini,corresponding to the tadpole of novel type B, including two genera, i.e., Leptolalaxand Leptobrachella. However, the relationships among the genera of Megophryinaewere largely unresolved, they recognized no monophyletic groups above the generalevel. A more thorough sampling will likely foster a better taxonomic solution.3) The larval evolutionary scenario in Megophryidae.Type A is characteristicof normal-mouthed with multiple tooth rows, representing the tadpole type of theMRCA of Chinese megophryids. Type B is characteristic of normal-mouthed withreduced tooth rows, prolonging labium, and integumetary glands. Type C ischaracteristic of no labial teeth and smaller umbeliform oral disc. Type D ischaracteristic of no labial teeth, enlarged umbeliform oral disc, representing the tadpoleof the MRCA of subfamily Megophryinae. A previous hypothesis, referring tofunnel-shaped oral discs as an apomorphy, is supported.4) The larval adaptation to habitats in Megophryidae. Tadpoles generallyadhere to substrates using their mouths, and the microhabitat that the tadpoles occupyreflects the degree of adhesion and oral complexity. The morphological changes inmegophryid tadpoles virtually allow a progressive adaptation to a changing habitatfrom faster water to slower water. Within the tadpoles of Type A to type D, the TOTbecomes smaller and smaller, and the oral disc orientates from anteroventral toumbelliform upturned, and eye position orientates from dorsal to lateral, and the trunkis more and more depressed and tail becomes relatively longer and slender. Within therunning water, the normal-mouthed with multiple tooth rows of Leptobrachiini tadpoles are correlated with lotic-suctorial, benthic feeders with anteroventral oraldisc and the largest body. With the water’s velocity decreasing, the lotic-adherentfeeders of Leptolalax tadpoles have tube-shaped labium with reduced tooth rows andintegumetary glands. And then, the smaller umbeliform in Brachytarsophrys tadpolesand the enlarged umbeliform oral disc in the Megophryini tadpoles are inhabitmicrohabitats of non-flowing backwaters of rivers, indicative of adaptive traits oflotic-neustonic surface feeders. The scheme of megophryid tadpoles andmicrohabitats provided the first clear evidence which congruent with the hypothesis ofAltig and Johnston (1989). The ecological divergence plays a general role in thedivergence and evolution of megophrid larvae. There is a definite correlation amongthe buccopharyngeal cavity, diet and feeding mechanisms, the tadpole graze orswallow the food particles, then through papillae which like a sieve and sort out foodparticles to the oesophagus. The tadpole of Leptobrachiinae possess multiple toothrows, wide intertooth distance as well as thick and sparse jaw sheath, these tadpolesinhabit bottom of the streams and graze on epiphyton or major detritus of organicmatter on the substrates, their prelingual papillae like single finger, the mechanicalpurpose of papillae served share in by tooth and jaw. The tadpoles of Megophryinaeoccur near the water surface of small streams and are the filter feeder, their dietincludes plankton and organic debris floating on the water surface, those tadpolepossess weak jaw, their prelingual papillae like spoon, the mechanical purpose ofpapillae served mostly for sieve.

Generation of large-scale vortex dislocations in a three-dimensional wake-type flow

Numerical study of three-dimensional evolution of wake-type flow and vortex dislocations is performed by using a compact finite diffenence-Fourier spectral method to solve 3-D incompressible Navier-Stokes equations. A local spanwise nonuniformity in momentum defect is imposed on the incoming wake-type flow. The present numerical results have shown that the flow instability leads to three-dimensional vortex streets, whose frequency, phase as well as the strength vary with the span caused by the local nonuniformity. The vortex dislocations are generated in the nonuniform region and the large-scale chain-like vortex linkage structures in the dislocations are shown. The generation and the characteristics of the vortex dislocations are described in detail.

Structure of Vortex Dislocations in a Wake-type Flow

Structure and dynamical processes of vortex dislocations in a kind of wake-type flow are described clearly by vortex lines, which are directly constructed from data of three-dimensional direct numerical simulations of the flow evolution.

Investigation of Mn-Implanted n-Type Ge

Mn+ irons were implanted to n-type Ge(1 1 1) single crystal at room temperature with an energy of 100 keV and a dose of 3 x 10(16) cm(-2). Subsequently annealing was performed at 400degreesC for 1 h under flowing nitrogen gas. X-ray diffraction measurements show that as-implanted sample is amorphous and the structure of crystal is restored after annealing. Polycrystalline germanium is formed in annealed sample. There are no new phases found except germanium. The samples surface morphologies indicate that annealed sample has island-like feature while there is no such kind of characteristic in as-implanted sample. The elemental composition of annealed sample was analyzed by Auger electron spectroscopy. It shows that manganese ions are deeply implanted into germanium substrate and the highest manganese atomic concentration is 8% at the depth of 120 nm. The magnetic properties of samples were investigated by an alternating gradient magnetometer. The annealed sample shows ferromagnetic behavior at room temperature.

Low-Dimensional Dynamical Systems for a Wake-Type Shear Flow with Vortex Dislocations

Low-dimensional systems are constructed to investigate dynamics of vortex dislocations in a wake-type shear flow. High-resolution direct numerical simulations are employed to obtain flow snapshots from which the most energetic modes are extracted using proper orthogonal decomposition (POD). The first 10 modes are classified into two groups. One represents the general characteristics of two-dimensional wake-type shear flow, and the other is related to the three-dimensional properties or non-uniform characteristics along the span. Vortex dislocations are generated by these two kinds of coherent structures. The results from the first 20 three-dimensional POD modes show that the low- dimensional systems have captured the basic properties of the wake-type shear flow with vortex dislocation, such as two incommensurable frequencies and their beat frequency.

Magnetic Properties of Mn-Implanted n-Type Ge

Mn+ ions were implanted into n-type Ge(111) single crystal at room temperature at an energy of 100 keV with a dose of 3 x 1016 cm-2. Subsequent annealing was performed on the samples at 400 °C and 600 °C in a flowing nitrogen atmosphere. The magnetic properties of the samples have been investigated by alternating gradient magnetometer at room temperature. The compositional properties of the annealed samples were studied by Auger electron spectroscopy and the structural properties were analyzed by X-ray diffraction measurements. Magnetization measurements reveal room-temperature ferromagnetism for the annealed samples. The magnetic analysis supported by compositional and structural properties indicates that forming the diluted magnetic semiconductor (DMS) MnxGe1-x after annealing may account for the ferromagnetic behavior in the annealed samples.

Quantifying Cell Binding Kinetics Mediated By Surface-Bound Blood Type B Antigen To Immobilized Antibodies

Cell adhesion is crucial to many biological processes, such as inflammatory responses, tumor metastasis and thrombosis formation. Recently a commercial surface plasmon resonance (SPR)-based BIAcore biosensor has been extended to determine cell binding mediated by surface-bound biomolecular interactions. How such cell binding is quantitatively governed by kinetic rates and regulating factors, however, has been poorly understood. Here we developed a novel assay to determine the binding kinetics of surface-bound biomolecular interactions using a commercial BIAcore 3000 biosensor. Human red blood cells (RBCs) presenting blood group B antigen and CM5 chip bearing immobilized anti-B monoclonal antibody (mAb) were used to obtain the time courses of response unit, or sensorgrams, when flowing RBCs over the chip surface. A cellular kinetic model was proposed to correlate the sensorgrams with kinetic rates. Impacts of regulating factors, such as cell concentration, flow duration and rate, antibody-presenting level, as well as pH value and osmotic pressure of suspending medium were tested systematically, which imparted the confidence that the approach can be applied to kinetic measurements of cell adhesion mediated by surface-bound biomolecular interactions. These results provided a new insight into quantifying cell binding using a commercial SPR-based BIAcore biosensor.

Lattice type of fracture model for concrete

Concrete is usually described as a three-phase material, where matrix, aggregate and interface zones are distinguished. The beam lattice model has been applied widely by many investigators to simulate fracture processes in concrete. Due to the extremely large computational effort, however, the beam lattice model faces practical difficulties. In our investigation, a new lattice called generalized beam (GB) lattice is developed to reduce computational effort. Numerical experiments conducted on a panel subjected to uniaxial tension show that the GB lattice model can reproduce the load-displacement curves and crack patterns in agreement to what are observed in tests. Moreover, the effects of the particle overlay on the fracture process are discussed in detail. (C) 2007 Elsevier Ltd. All rights reserved.

The mixed-type integral equations for transient elastodynamic plane crack problems

In the present paper, by use of the boundary integral equation method and the techniques of Green fundamental solution and singularity analysis, the dynamic infinite plane crack problem is investigated. For the first time, the problem is reduced to solving a system of mixed-typed integral equations in Laplace transform domain. The equations consist of ordinary boundary integral equations along the outer boundary and Cauchy singular integral equations along the crack line. The equations obtained are strictly proved to be equivalent with the dual integral equations obtained by Sih in the special case of dynamic Griffith crack problem. The mixed-type integral equations can be solved by combining the numerical method of singular integral equation with the ordinary boundary element method. Further use the numerical method for Laplace transform, several typical examples are calculated and their dynamic stress intensity factors are obtained. The results show that the method proposed is successful and can be used to solve more complicated problems.

A new type of boundary integral-equation for plane problems of elasticity including rotational forces

Based on the idea proposed by Hu [Scientia Sinica Series A XXX, 385-390 (1987)], a new type of boundary integral equation for plane problems of elasticity including rotational forces is derived and its boundary element formulation is presented. Numerical results for a rotating hollow disk are given to demonstrate the accuracy of the new type of boundary integral equation.

Vortex dislocation in wake-type flow caused by local spanwise nonuniformity

DNS of spatiotemporal evolution of a wake-type flow is performed. In the incoming flow, a local spanwise nonuniformity in momentum defect is initially imposed. Results show that the spanwise nonuniformity leads to a series of symmetric twist vortex dislocation in downstream of the flow. Vortex line variations and substantial transition of vorticity from spanwise to the streamwise and vertical directions clearly feature the generation of a vortex dislocation and the real vortex linking in the dislocation. Dynamical process and the mechanism responsible for the vortex dislocation are described.

A simple method to simulate shrinkage-induced cracking in cement-based composites by lattice-type modeling

A new numerical procedure is proposed to investigate cracking behaviors induced by mismatch between the matrix phase and aggregates due to matrix shrinkage in cement-based composites. This kind of failure processes is simplified in this investigation as a purely spontaneous mechanical problem, therefore, one main difficulty during simulating the phenomenon lies that no explicit external load serves as the drive to propel development of this physical process. As a result, it is different from classical mechanical problems and seems hard to be solved by using directly the classical finite element method (FEM), a typical kind of "load -> medium -> response" procedures. As a solution, the actual mismatch deformation field is decomposed into two virtual fields, both of which can be obtained by the classical FEM. Then the actual response is obtained by adding together the two virtual displacement fields based on the principle of superposition. Then, critical elements are detected successively by the event-by-event technique. The micro-structure of composites is implemented by employing the generalized beam (GB) lattice model. Numerical examples are given to show the effectiveness of the method, and detailed discussions are conducted on influences of material properties.