Habitat suitability analysis of eelgrass Zostera marina L.in the subtidal zone of Xiaoheishan Island

We present a GIS-based habitat suitability index(HSI) model to identify suitable areas for Zostera marina L. restoration in the subtidal zone of Xiaoheishan Island. The controlling factors in the model,in order of importance,are Secchi depth,sediment composition,water temperature,salinity,current velocity,water depth and nutrient quality. Specific factor piecewise functions have been used to transform parameter values into normalized quality indexes. The weight of each factor was defined using expert knowledge and the analytic hierarchy process(AHP) method. All of the data thus obtained were interpolated using the inverse distance weighted(IDW) interpolation method to create maps for the entire region. In this study,the analysis of habitat suitability in the subtidal zone of Xiaoheishan Island was conducted for four seasons. According to the GIS-based HSI model,the optimal habitat of Z ostera marina L. appears in spring,although habitat remains suitable all year round. On the whole,the optimum site for eelgrass restoration is located in the eastern region,followed by the western and southern regions. We believe that the GIS-based HSI model could be a promising tool to select sites for Z ostera marina L. restoration and could also be applicable in other types of habitat evaluation.     

Identifying suitable habitats of three ungulates in Arjinshan National Nature Reserve,China

Arjinshan National Nature Reserve (ANNR) is one of 3 refuges of the endemic ungulates Tibetan wild ass (Equus kiang), Chiru (Pantholops hodgsonii) and wild yak (Bos mutus) that are endangered by natural and anthropogenic factors in the Qinghai-Tibetan Plateau (QTP). Identifying habitat suitability is critical for species protection management. We used a GIS-based niche model to delimit and classify habitat suitability using an integrated assessment system, which included 9 biotic and abiotic factors. We divided the research area into 4 habitat types for these ungulates: (1) high suitability habitat; (2) moderate suitability habitat; (3) low suitability habitat; and (4) unsuitable area. Results suggested that chiru have the most areas of high and moderate suitability habitats while Tibetan wild ass had the largest areas of low suitability habitat and unsuitable area. Wild yak had the largest area of moderate and low suitability habitat, but high suitability habitats for wild yak were smaller than those of other 2 ungulates. There was overlap of high-quality habitat for the three kinds of ungulates in the vicinity of Kardun inspection station, which could be regarded as the core area for the coexistence and conservation of these endangered ungulate populations.     

Seasonal variation in species composition and abundance of demersal fish and invertebrates in a Seagrass Natural Reserve on the eastern coast of the Shandong Peninsula,China

Seagrass habitats are structurally complex ecosystems, which support high productivity and biodiversity. In temperate systems the density of seagrass may change seasonally, and this may influence the associated fish and invertebrate community. Little is known about the role of seagrass beds as possible nursery areas for fish and invertebrates in China. To study the functioning of a seagrass habitat in northern China, demersal fish and invertebrates were collected monthly using traps, from February 2009 to January 2010. The density, leaf length and biomass of the dominant seagrass Zostera marina and water temperature were also measured. The study was conducted in a Seagrass Natural Reserve(SNR) on the eastern coast of the Shandong Peninsula, China. A total of 22 fish species and five invertebrate species were recorded over the year. The dominant fish species were Synechogobius ommaturus, Sebastes schlegelii, Pholis fangi, Pagrus major and Hexagrammos otakii and these species accounted for 87% of the total number of fish. The dominant invertebrate species were Charybdis japonica and Octopus variabilis and these accounted for 98% of the total abundance of invertebrates. There was high temporal variation in species composition and abundance. The peak number of fish species occurred in August–October 2009, while the number of individual fish and biomass was highest during November 2009. Invertebrate numbers and biomass was highest in March, April, July and September 2009. Temporal changes in species abundance of fishes and invertebrates corresponded with changes in the shoot density and leaf length of the seagrass, Zostera marina.     

Predicting the impact of climate change on the distribution of two threatened Himalayan medicinal plants of Liliaceae in Nepal

Predicting the potential distribution of medicinal plants in response to climate change is essential for their conservation and management. Contributing to the management program, this study aimed to predict the distribution of two threatened medicinal plants, Fritillaria cirrhosa and Lilium nepalense. The location of focal species gathered from herbarium specimen housed in different herbaria and online databases were geo-referenced and checked for spatial autocorrelation. The predictive environmental variables were selected, and MaxEnt software was used to model the current and future distributions of focal species. Four Representative Concentration Pathway (RCP) trajectories of the BCC-CSM1.1 model were used as the future (2050) projection layer. The MaxEnt modelling delineated the potential distribution of F. cirrhosa and L. nepalense. The current suitability is projected towards Central and Eastern Hilly/Mountainous regions. Both species gain maximum suitability in RCP 4.5 which decline towards other trajectories for L. nepalense. Overall, both the focal species shift towards the north-west, losing their potential habitat in hilly and lower mountainous regions by 2050 across all trajectories. Our results highlight the impact of future climate change on two threatened and valuable species. The results can be further useful to initiate farming of these medicinally and economically important species based on climatically suitable zone and for designing a germplasm conservation strategy.     

Seasonal variations in the phytoplankton community and the relationship between environmental factors of the sea around Xiaoheishan Island in China

Seasonal variations in the phytoplankton community and the relationship between environmental factors of the sea area around Xiaoheishan Island are investigated in the present study. Xiaoheishan Island is located at 37°58′14″N and 120°38′46″E in Shandong Province, China. A total of 65 species of phytoplankton belonging to three phyla and 27 genera were identified, with Bacillariophyta having the largest number of species. The annual average chlorophyll a concentration for this area was 3.11 μg/L, and there occurs a Skeletonema costatum bloom in winter. The Shannon-Weaver indexes(log_2) of the phytoplankton from all stations were higher than 1, and the Pielou indexes were all higher than 0.3. The results of the canonical correspondence analysis(CCA) indicated that water temperature, PO_4~(3ˉ) and Cu were the environmental factors that had the greatest influence on the distribution of the phytoplankton community throughout the entire year. Although the concentration of heavy metal is well up to the state standards of the first grade of China(GB 3097-1997), these metals still have an impact on the phytoplankton community from this area.     

Impact of <Emphasis Type="Italic">Vibrio parahaemolyticus</Emphasis> and white spot syndrome virus (WSSV) co-infection on survival of penaeid shrimp <Emphasis Type="Italic">Litopenaeus vannamei</Emphasis>

White spot syndrome virus (WSSV) is an important viral pathogen that infects farmed penaeid shrimp, and the threat of Vibrio parahaemolyticus infection to shrimp farming has become increasingly severe. Viral and bacterial cross or superimposed infections may induce higher shrimp mortality. We used a feeding method to infect Litopenaeus vannamei with WSSV and then injected a low dose of V. parahaemolyticus (WSSV+Vp), or we first infected L. vannamei with a low-dose injection of V. parahaemolyticus and then fed the shrimp WSSV to achieve viral infection (Vp+WSSV). The eff ect of V. parahaemolyticus and WSSV co-infection on survival of L. vannamei was evaluated by comparing cumulative mortality rates between experimental and control groups. We also spread L. vannamei hemolymph on thiosulfate citrate bile salt sucrose agar plates to determine the number of Vibrio, and the WSSV copy number in L. vannamei gills was determined using an absolute quantitative polymerase chain reaction (PCR) method. LvMyD88 and Lvakt gene expression levels were detected in gills of L. vannamei by real-time PCR to determine the cause of the diff erent mortality rates. Our results show that (1) the cumulative mortality rate of L. vannamei in the WSSV+Vp group reached 100% on day 10 after WSSV infection, whereas the cumulative mortality rate of L. vannamei in the Vp+WSSV group and the WSSV-alone control group approached 100% on days 11 and 13 of infection; (2) the number of Vibrio in the L. vannamei group infected with V. parahaemolyticus alone declined gradually, whereas the other groups showed significant increases in the numbers of Vibrio (P<0.05); (3) the WSSV copy numbers in the gills of the WSSV+Vp, Vp+WSSV, and the WSSV-alone groups increased from 10⁵ to 10⁷ /mg tissue 72, 96, and 144 h after infection, respectively. These results suggest that V. parahaemolyticus infection accelerated proliferation of WSSV in L. vannamei and vice versa. The combined accelerated proliferation of both V. parahaemolyticus and WSSV led to massive death of L. vannamei.     

Mapping soil organic carbon stocks of northeastern China using expert knowledge and GIS-based methods

The main aim of this paper was to calculate soil organic carbon stock (SOCS) with consideration of the pedogenetic horizons using expert knowledge and GIS-based methods in northeastern China. A novel prediction process was presented and was referred to as model-then-calculate with respect to the variable thicknesses of soil horizons (MCV). The model-then-calculate with fixed-thickness (MCF), soil profile statistics (SPS), pedological professional knowledge-based (PKB) and vegetation type-based (Veg) methods were carried out for comparison. With respect to the similar pedological information, nine common layers from topsoil to bedrock were grouped in the MCV. Validation results suggested that the MCV method generated better performance than the other methods considered. For the comparison of polygon based approaches, the Veg method generated better accuracy than both SPS and PKB, as limited soil data were incorporated. Additional prediction of the pedogenetic horizons within MCV benefitted the regional SOCS estimation and provided information for future soil classification and understanding of soil functions. The intermediate product, that is, horizon thickness maps were fluctuant enough and reflected many details in space. The linear mixed model indicated that mean annual air temperature (MAAT) was the most important predictor for the SOCS simulation. The minimal residual of the linear mixed models was achieved in the vegetation type-based model, whereas the maximal residual was fitted in the soil type-based model. About 95% of SOCS could be found in Argosols, Cambosols and Isohumosols. The largest SOCS was found in the croplands with vegetation of Triticum aestivum L., Sorghum bicolor (L.) Moench, Glycine max (L.) Merr., Zea mays L. and Setaria italica (L.) P. Beauv.     

Effects of different phosphorus concentrations and N/P ratios on the growth and photosynthetic characteristics of <Emphasis Type="Italic">Skeletonema costatum</Emphasis> and <Emphasis Type="Italic">Prorocentrum donghaiense</Emphasis>

The effects of different phosphorus (P) concentrations (0.36, 3.6, and 36 μmol/L corresponding to low-, middle-, and high-P concentration groups, respectively) and nitrogen (N)/P ratios on the growth and photosynthetic characteristics of Skeletonema costatum and Prorocentrum donghaiense were studied. For both species, the high-P (HP) concentration group showed the greatest algal density and highest specific growth rate. Changes in the maximum efficiency of photosystem II (F _v /F _m ) were monitored under the various P and N/P conditions. The largest decrease in F _v /F _m was in the low-P (LP) group in S. costatum and in the HP group in P. donghaiense. There were high rapid light curves and photochemical quantum yields (Φ _PSII) for S. costatum in the HP group, while the actual photosynthetic capacity was higher in P. donghaiense than in S. costatum in the MP group. Under eutrophic but relatively P-restricted conditions, P. donghaiense had higher photosynthetic activity and potential, which could cause this dinoflagellate to increasingly dominate the phytoplankton community in these conditions. Under the same P concentration and N/P ratio, P. donghaiense had a larger relative maximum rate of electron transport and higher Φ _{PS II} values than those of S. costatum. These differences between P. donghaiense and S. costatum may explain the interaction and succession patterns of these two species in the Changjiang (Yangtze) River estuary from a photosynthesis perspective.     

Identification of three common Loliginidae squid species in the South China Sea by analyzing hard tissues with geometric outline method

The hard tissues of squid can provide important information for species identification. In this study, we used statolith and beak to identify three squid species including Uroteuthis duvaucelii, Loliolus beka, and U. edulis in the South China Sea. Because of the highly overlapping habitat and similar body morphology of the three squid species, we explored four different ways to identify them, by using statolith, upper beak, lower beak and a combination of statolith and beak. An outline geometric morphometric method and stepwise discriminant analysis were used to evaluate the most suitable method for the identification. We found that the combination of statolith and beak had the highest cross validation rate that was 75.0%, 87.5% and 88.7% for U. duvaucelii, L. beka and U. edulis, respectively. Using two beaks had similar results and the lowest cross validation rate was 60.0%, 50.0%, and 73.7% for the upper beak, 46.9%, 58.5% and 75.3% for the lower beak of U. duvaucelii, L. beka and U. edulis, respectively. Analyzing with the statolith had moderate cross validation which was 72.2%, 80.0%, and 87.7% for U. duvaucelii, L. beka and U. edulis, respectively. From the results it is suggested when the entire body of a squid is available, a combination of statolith and beak should be used for the identification. When only one hard tissue is available, species identification can be subjected to large errors.     

Biodeposition,respiration, and excretion rates of an introduced clam <Emphasis Type="Italic">Mercenaria mercenaria</Emphasis> in ponds with implications for potential competition with the native clam <Emphasis Type="Italic">Meretrix meretrix</Emphasis> in Shuangtaizi estuary,China

This study aimed to evaluate the potential impacts of an introduced clam Mercenaria mercenaria on estuarine ecosystem, and implications for the niche competition with a native clam Meretrix meretrix. The biodeposition, respiration, and excretion rates of M. mercenaria were determined seasonally using a sediment trap and a closed respirator in field. The biodeposition rates of M. mercenaria were 0.06–0.37 g/ (ind.·d), and the respiration rates were 0.31–14.66 mg/(ind.·d). The ammonia and phosphate excretion rates were 0.18–36.70 and 1.44–14.87 μg/(ind.·d), respectively. The hard clam M. mercenaria may discharge dry deposits up to 2.1×10⁵ t, contribute 18.3 t ammonia and 9.0 t phosphate to culture ponds, and consume 7.9×10³ t O₂ from ponds annually. It suggested that the hard clam M. mercenaria might play an important role in pelagic-benthic coupling in pond ecosystem through biodeposition and excretion. A comparison of the key physiological parameters of the introduced clam M. mercenaria and the native clam Meretrix meretrix suggested that M. mercenaria had a niche similar to that of Meretrix meretrix in Shuangtaizi estuary and might have a potential competition with Meretrix meretrix for habitat and food if M. mercenaria species escaped from the culture pond or artificially released in estuarine ecosystem.     

Loss of preferred habitat and pollen limitation threatens reproduction in a rare mountain <Emphasis Type="Italic">Paeonia delavayi</Emphasis>

Paeonia delavayi is a wild tree peony species endemic to high-altitude regions in southwestern China. Recent agricultural land expansion, however, not only causes its dwindling population size, but also poses a severe threat to its long-term persistence. Since our knowledge of the reproductive biology of P. delavayi is very limited, and some management misconceptions have further exacerbated the already worrisome situation, the aim of the present study is to provide some scientific evidence regarding habitat preference and breeding system of P. delavayi, and to help correct some management misconceptions. Specifically, morphological traits of P. delavayi populations growing in two different habitats (i.e., the open area versus forest understory), including basal stem diameter, plant height and flower numbers per plant, were measured and compared with each other. A range of floral manipulation experiments was further conducted to study the extent of self-compatibility and pollen limitation. The results showed that P. delavayi preferred the open habitat over forest understory as it generally grew bigger with thicker stems, and produced a higher number of flowers in the former habitat. Therefore, the open habitat should receive immediate priority for conservation action instead of being converted to agricultural land; also, P. delavayi is pollen limited, as evidenced by the fact that experimental flowers receiving supplemental non-self pollen generally displayed elevated seed production than flowers subjected to other treatments (e.g., non-supplemental self pollen, non-supplemental non-self pollen and supplemental self pollen), which is consistent with reported observations that alpine plants are more likely to be pollen-limited than lowland plants. We suggest that human intervention might be necessary to guarantee the long-term persistence of P. delavayi as harsh alpine environment, intense competition for pollinators and different anthropogenic perturbations co-limit its reproductive success.     

<Emphasis Type="Italic">GST</Emphasis> (<Emphasis Type="Italic">phi</Emphasis>) gene from Macrophyte <Emphasis Type="Italic">Lemna minor</Emphasis> is involved in cadmium exposure responses

Reactive oxygen species (ROS) scavengers, including ascorbate peroxidase, superoxide dismutase, catalase and peroxidase, are the most commonly used biomarkers in assessing an organisms’ response to many biotic and abiotic stresses. In this study, we cloned an 866 bp GST (phi) gene in Lemna minor and investigated its characteristics, expression and enzymatic activities under 75 μmol/L cadmium concentrations in comparison with other ROS scavengers. GST (phi) gene expression patterns were similar to those of other scavengers of ROS. This suggests that GST (phi) might be involved in responding to heavy metal (cadmium) stress and that its expression level could be used as a bio-indicator in monitoring cadmium pollution.     

Effects of nitrogen and phosphorus on the growth of <Emphasis Type="Italic">Levanderina fissa</Emphasis>: How it blooms in Pearl River Estuary

Effects of nitrogen (N) and phosphorus (P) from different sources and at different concentrations on the growth of Levanderina fissa (= Gyrodinium instriatum) were studied in laboratory conditions. The findings might explain the recurrent blooms of this species in Pearl River Estuary, China. Results showed that nutrient limitation significantly inhibited the growth of L. fissa. The values of specific growth rate (μ _max) and half-saturation nutrient concentration (K _S) were 0.37 divisions/d and 8.49 μmol L^?1 for N, and 0.39 divisions/d and 1.99 μmol L^?1 for P, respectively. Based on K _S values, dissolved inorganic N level in PRE was sufficient to support the high proliferation of L. fissa, while dissolved inorganic P concentration was far lower than the minimum requirement for its effective growth. L. fissa was not able to utilize dissolved organic N (DON) compounds such as urea, amino acids, and uric acid. However, it grew well by using a wide variety of dissolved organic P (DOP) sources like nucleotides, glycerophosphate, and 4-nitrophenylphosphate. The results from this study suggested that the ability in DOP utilization of L. fissa offers this species a competitive advantage in phytoplankton communities. The high level and continuous supply of DIN, enrichment of DOP, together with warm climate and low salinity in the Pearl River Estuary provided a suitable nutrient niche for the growth of L. fissa, and resulted in the recurrent blooms in the estuary.     

Diel variation in metabolism and ammonia excretion of <Emphasis Type="Italic">Marphysa sanguinea</Emphasis> (Polychaeta: Eunicidae)

Polychaetes provide an excellent food resource for fish and represent the dominant zoobenthos in marine ecosystems. Diel variation in the rates of metabolism and ammonia-N excretion of Marphysa sanguinea were studied. The worms were grouped according to their wet body weight into small (S; 1.24±0.06 g), medium (M; 4.00±0.30 g), and large (L; 8.54±1.08 g) categories. Their weight-specific metabolic rates, based on aerobic respiration (R), were measured at 16°C (±0.2°C) and classed as either routine (R _R) or standard (R _S) rates. Both respiration types decreased with increasing body weight. Respiration was described by R = a W ^b, where b was -0.400 9 and -0.532 0 for R _R and R _S, respectively. Diurnal changes in R _S for each group was relatively flat, with a slightly increasing trend with time, but was relatively stable as a whole. R _R of the diurnal variation of worms was higher than R S, but both had similar overall trends. The peak values of specific dynamic action (SDA) (R _SDA) in the S, M, and L groups were 2.704, 1.149, and 0.682 mg/(g?h), respectively. The durations of SDA were 13, 6, and 6 h, respectively and the energy expenditures of SDA were 377.98, 117.34, and 74.94 J/g, respectively. These data indicate that the metabolic rates were higher in smaller individuals, which is advantageous for their rapid growth.     

Effects of dietary genistein on GH/IGF-I axis of Nile tilapia <Emphasis Type="Italic">Oreochromis niloticus</Emphasis>

There is considerable concern that isoflavones, such as genistein in fish feed composed of soybean protein, aff ects somatic growth in fish. Our previous works demonstrated that 30 and 300 μg/g dietary genistein had no significant eff ect on growth performance in Nile tilapia (Oreochromis niloticus), but the higher level of genistein (3 000 μg/g) significantly depressed growth. This study was conducted to further examine the eff ects of dietary genistein on the endocrine disruption on growth hormone/insulin-like growth factor-I (GH/IGF-I) axis in Nile tilapia (O. niloticus). Juvenile fish were fed by hand twice daily to satiation with one of four isonitrogenous and isoenergetic diets, each containing either 0, 30, 300 or 3 000 μg/g genistein. Following an 8-week feeding period, plasma GH and IGF-I levels were investigated by radioimmunoassay and gene expression levels of gh, ghrelin, gnrhs, ghr, npy, npyrs, pacap, ghrs, igf-I, igf-Ir, and igfbp3 were examined by real-time PCR. The results show that no significant change in plasma GH and IGF-I levels in fish fed with diets containing 30 μg/g and 300 μg/g genistein. mRNA expression of genes along the GH/IGF-I axis remained unaff ected, except for igf-Ir, which was stimulated by the 300 μg/g genistein diet. While in fish fed the 3 000 μg/g genistein diet, the plasma GH and IGF-I levels decreased, and mRNA expression of gh, ghr2, npyr1, igf-I, and igf-Ir were also significantly depressed. In contrast, npy and igfbp3 mRNA expression were enhanced. This study provides convincing evidence for growth impediment by genistein by disturbing the GH/IGF-I axis in Nile tilapia O. niloticus.     

Immunogenicity and protective role of antigenic regions from five outer membrane proteins of <Emphasis Type="Italic">Flavobacterium columnare</Emphasis> in grass carp <Emphasis Type="Italic">Ctenopharyngodon idella</Emphasis>

Flavobacterium columnare causes columnaris disease in freshwater fish. In the present study, the antigenic regions of five outer membrane proteins (OMPs), including zinc metalloprotease, prolyl oligopeptidase, thermolysin, collagenase and chondroitin AC lyase, were bioinformatically analyzed, fused together, and then expressed as a recombinant fusion protein in Escherichia coli. The expressed protein of 95.6 kDa, as estimated by 10% sodium dodecyl sulfate-polyacrylamide gel electrophoresis, was consistent with the molecular weight deduced from the amino acid sequence. The purified recombinant protein was used to vaccinate the grass carp, Ctenopharyngodon idella. Following vaccination of the fish their IgM antibody levels were examined, as was the expression of IgM, IgD and IgZ immunoglobulin genes and other genes such as MHC Iα and MHC IIβ, which are also involved in adaptive immunity. Interleukin genes (IL), including IL-1β, IL-8 and IL-10, and type I and type II interferon (IFN) genes were also examined. At 3 and 4 weeks post-vaccination (wpv), significant increases in IgM antibody levels were observed in the fish vaccinated with the recombinant fusion protein, and an increase in the expression levels of IgM, IgD and IgZ genes was also detected following the vaccinations, thus indicating that an adaptive immune response was induced by the vaccinations. Early increases in the expression levels of IL and IFN genes were also observed in the vaccinated fish. At four wpv, the fish were challenged with F. columnare, and the vaccinated fish showed a good level of protection against this pathogen, with 39% relative percent survival (RPS) compared with the control group. It can be concluded, therefore, that the five OMPs, in the form of a recombinant fusion protein vaccine, induced an immune response in fish and protection against F. columnare.     

A review of interaction between neon flying squid (<Emphasis Type="Italic">Ommastrephes bartramii</Emphasis>) and oceanographic variability in the North Pacific Ocean

The neon flying squid (Ommastrephes bartramii) is a short-lived opportunistic species widely distributed in subtropical and temperate waters in the North Pacific Ocean. The life cycle of O. bartramii from planktonic eggs to nektonic adults is closely linked to oceanographic conditions. The fluctuations in O. bartramii abundance and distribution tend to increase and widen continuously due to the heavy influences of ocean-climate events on various spatio-temporal scales. In this study, we reviewed the interaction between O. bartramii and oceanography variability in the North Pacific with respect to large-scale climatic-oceanic phenomena including El Niño, La Niña, Kuroshio, Oyashio and Pacific Decadal Oscillation (PDO), as well as regional environmental variables such as sea surface temperature (SST), sea surface height (SSH), sea surface salinity (SSS), chlorophyll-a (Chl-a) concentration, and plankton density. The population dynamics of O. bartramii is mediated mainly by meso- and large-scale climatic-oceanic events (e.g., Kuroshio and Oyashio Currents) rather than other local environmental conditions (e.g., SST and Chl-a concentration), because all of the oceanographic influences are imposed on the context of large-scale climate changes (e.g., PDO). An unstructured-grid finite-volume coastal ocean model coupled with an individual-based model is proposed to simulate relevant physical-biological oceanographic processes for identifying ocean-climate influence and predicting O. bartramii distribution and abundance in the North Pacific. Future research needs to be focused on improving the knowledge about early life history of O. bartramii and evaluating the relationship between marine physical environment and two separate passive drifting life stages of O. bartramii including free-floating eggs and planktonic paralarvae.     

Inbreeding depression on growth and survival of full-sib family of Manila clam (<Emphasis Type="Italic">Ruditapes philippinarum</Emphasis>)

In present study, the inbreeding depression (ID) of growth and survival of Manila clam (Ruditapes philippinarum) was investigated at larval and juvenile stages. Nine inbred families (A ₂, B ₂, C ₂, D ₂, E ₂, F ₂, G ₂, H ₂ and I ₂) were established by mating within nine full-sib families with expected inbreeding coefficient of 0.25. Inbred families showed significant differences in shell length and hatching rate of D-larvae (straight-hinged larvae). The larvae of the nine inbred families grew slower than those of control group (CG), and their ID value ranged from 0.81% ± 6.09% to 16.10% ± 1.49%. The ID value of larval survival rate varied between 27.47% ± 9.36% and 70.50% ± 13.66%. The ID was also detected for juvenile growth in A ₂, B ₂, C ₂, and D ₂, which ranged from 4.60 ± 2.21 to 17.71 ± 7.73. The A ₂ family maintained the highest juvenile survival rate, whereas the other inbred families exhibited ID values varying between 62.79% ± 4.54% and 96.14% ± 0.87%. The linear relationship of estimated ID between growth and survival was negatively correlated (R = ?0.434, P < 0.05). The results of this study suggested that the ID of growth was common at the larval stage but was less prevalent at juvenile stage. In contrast, the ID of survival increased from larval to juvenile stage. A better understanding of the effect of inbreeding may aid to selective breeding of Manila clam.