DHA Accumulation in the Polar Lipids of the Euryhaline Copepod Pseudodiaptomus inopinus and Its Transfer to Red Sea Bream Pagrus major Larvae

The euryhaline copepod Pseudodiaptomus inopinus play important roles in coastal waters as vectors of docosahexanoic (DHA) and eicosapentaenoic (EPA) acids for larval fish. While DHA and EPA in polar lipids (PLs) are more effective for fish larval development than non-polar lipid forms (NLs), there is little knowledge how much these lipids are accumulated in copepods from microalgae and are effective for early development of fish larvae. We report PLs fatty acid profiles of P. inopinus fed DHA-poor microalgae and evaluate its significance as a food source for larvae development of Pagrus major, compared with DHA-enriched rotifers. Copepods and rotifers were fed a mixed diet of three algal species (Phaeodactylum tricornutum, Tisochrysis lutea, and Pavlova lutheri), in addition of DHA-supplemented Super Fresh Chlorella (SFC) for rotifers. Compared with SFC, the algal mixture had higher EPA but lower DHA. Copepods had higher DHA and EPA in total lipids than rotifers fed each diet. Copepod PLs were specifically enriched with DHA and their contents were higher than both rotifers. On the other hand, PLs EPA contents were comparable between preys, indicating that copepods selectively fortified the PLs. Fish culture experiment showed that larvae fed copepods had higher growth than those fed SFC-enriched rotifers. Principal component analysis for each organism fatty acid composition emphasized trophic modification of DHA by copepods toward larval fish. This study highlighted that P. inopinus contribute to enhanced growth of coastal larval fish by efficiently transferring DHA via copepod fatty acid metabolism.

INTRODUCTION the importance of prey quality (i.e., nutrition) to explain fish stocks. Marine finfish lack the ability to synthesize n-3 highly unsaturated fatty acids (n-3 HUFAs), such as docosahexaenoic acid (DHA) and eicosapentaenoic acid (EPA), and therefore must supplement these essential fatty acids from diets for their growth and survival (Tocher, 2003). Indeed, larval fish reared under sufficient amounts of diets further enhance their growth, survival, and development when the diets are rich in n-3 HUFA (especially DHA; Watanabe et al., 1989;Nćss et al., 1995;Shields et al., 1999;Evjemo et al., 2003;Karlsen et al., 2015). Larvae inhabiting in the coastal areas, also exhibit higher growth rate when the ambient zooplankton contain higher n-3 HUFAs (Paulsen et al., 2014). These facts suggest the necessity to investigate how fatty acids are transferred toward larval fish for further understanding of fish recruitment.
Coastal and estuarine waters are known as a nursery area for various larval species and support their abundance (Borges et al., 2007). The euryhaline calanoid Pseudodiaptomus inopinus is a very common copepod species in such ecosystems around East Asia, Russia, and the United States (Sakaguchi and Ueda, 2011). The genus Pseudodiaptomus has a tolerance to a broad range of salinity (Beyrend-Dur et al., 2011), and high turbidity (Hwang et al., 2010), so that it is often the most abundant prey for larval fish. This copepod is reported to have high DHA and/or EPA in total lipids (Toledo et al., 1999;Rayner et al., 2015). Thus, P. inopinus with rich fatty acids can play important roles for larval development in coastal ecosystems.
On the other hand, in coastal and estuary waters, the dominant primary producers of n-3 HUFA are often diatoms (Kleppel, 1993;Trigueros and Orive, 2001;Ansotegui et al., 2003). They lack DHA despite having specifically high EPA (Kates and Volcani, 1966;Volkman et al., 1989). Indeed, DHA was relatively low in spring blooms where diatoms are dominated (Sargent et al., 1985). On the diatom-copepod-larval fish food chain, little is known how P. inopinus contributes to early development of larval fish by transferring DHA. To deal with this issue, fatty acid profiles, particularly polar lipids (PLs) of copepods, require further investigation (Olsen et al., 2014). For larval fish, dietary n-3 HUFA distributed to the PLs can be more efficiently absorbed than those distributed to the non-polar lipids (NLs; Salhi et al., 1999), resulting in enhanced larval growth and survival (Gisbert et al., 2005;Wold et al., 2009). However, the distribution of DHA and EPA in P. inopinus PLs is unknown.
Our objective is to evaluate the fatty acid profiles of PLs in P. inopinus cultivated using DHA-poor microalgae and their significance in early development of fish larvae. Rather than copepods, successful larviculture is well established using artificially DHA-enriched Brachionus rotifers as the diet (Koven et al., 1990;Castell et al., 2003;Thépot et al., 2016). In this study, we used such rotifers fed on a commercial DHA enrichment diet as a standard for evaluating copepod fatty acids. We also cultured P. major larvae with copepods, and compared the performances with those cultured using enriched rotifers. P. major is a widely distributed fish in coastal waters of Japan (Shoji et al., 2007), and is of commercial importance to Japan, Korea, and China (Hossain et al., 2016). We then discuss fatty acid trophic transfer in coastal ecosystems, and the significance of copepods on early development of coastal fish larvae.

Cultivation of Rotifers
The L-type rotifer Brachionus plicatilis species complex (Obama strain) was used as control diet for larval fish. As an experimental control, artificially enriched rotifers were prepared by feeding Chlorella vulgaris containing DHA (Super Fresh Chlorella V-12, SFC; Chlorella industry, Tokyo, Japan). SFC is produced by heterotrophic cultivation of C. vulgaris with fish oil DHA as a carbon source (Hayashi et al., 2001). Following the procedure of rotifer enrichment described by Waqalevu et al. (2018), SFC was fed to rotifers at 46 × 10 3 cells per rotifer per a day. Seawater for the culture medium was filtered through 100 µm cartridge filters (Micro-Cilia; Roki Techno Co. Ltd., Japan); water temperature was maintained at 25 • C and salinity at 17 . Taurine, an essential amino acid for red sea bream larvae (Chen et al., 2004), was also enriched in rotifers using Aquaplus ET (Marubeni Nisshin Feed Co, Ltd., Japan) following manufacturer instructions (60 mg L −1 ). The enriched rotifers were prepared daily during larval rearing experiment (for 20 days) and harvested to directly supply to larvae. In addition, for fatty acid analysis, the rotifers were directly harvested from three randomly-selected cultures, washed with distilled water, and then freeze-dried and preserved at −80 • C until measurement. Starvation periods to make rotifer digestive tracts empty were not set before harvesting them.
To compare the fatty acid composition of rotifers and copepods fed the same diet, the rotifers were also fed same algal mixture as the copepods. Rotifers were fed the mixture daily, cultivated at 20 • C water temperature and salinity 17 . When rotifer density reached 800-1,000 individuals ml −1 they were harvested and preserve for later fatty acid analysis (these rotifers were not used for fish larval diet). Such rotifer samples were prepared in triplicate as SFC-enriched ones.

Cultivation of Microalgae
Copepods were fed on the diatom Phaeodactylum tricornutum, with the haptophytes Pavlova lutheri and Tisochrysis lutea (originated from Isochrysis affinis galbana Tahiti strain; Bendif et al., 2013) as a mixture. We supplemented the diatom with the two haptophytes for mass production of copepod, since monoalgal diatom is well reported to suppress copepod egg productivity and survival (Støttrup and Jensen, 1990;Schipp et al., 1999;Shields et al., 2005;Ohs et al., 2010). The strains of P. tricornutum and P. lutheri were supplied by National Research Institute of Aquaculture, Fisheries Research Agency (Mie, Japan); and T. lutea strain was from the Graduate School of Fisheries and Environmental Science, Nagasaki University, Japan. Each microalga was incubated in 5 L polycarbonate bottles containing a seawater culture medium at 20 • C with continuous 0.5 L min −1 aeration and irradiance of 120 µmol photons m −2 s −1 . Seawater for the microalgae culture medium was filtered through a GF/F filter (Whatman), diluted to salinity 17 , autoclaved, and enriched with KW21 medium (Daiichi Seimo Frontiers in Marine Science | www.frontiersin.org Co. Ltd., Japan). Fresh cultures for each species were prepared weekly and maintained for about 1 week. Cells of the three algal species were counted every day using a Thoma hematocytometer (0.1 mm depth; Sunlead Glass Co. Ltd., Tokyo, Japan). Cells were harvested for a diet and fatty acid analysis once they reached the end of logarithmic growth phase or the beginning of stationary phase. Microalgae were directly used as a diet without temporary storage. Three samples of each algal species for fatty acid analysis were twice washed with 0.5 M ammonium formate, and then freeze-dried and preserved at −80 • C.

Cultivation of Copepods
Copepods (P. inopinus) were collected using a plankton net (mesh 63 µm) from the surface of Yakugachi River (Kagoshima, Japan), in waters of 19.7 • C and salinity 11.5 . Seawater for the copepod culture medium was pumped from Kagoshima Bay (Kagoshima, Japan), filtered through a graded series of 100, 25, and 10 µm pore sized cartridge filters (Micro-Cilia; Roki Techno Co. Ltd., Japan), before being subjected to ultraviolet irradiance (UVF-1000; Iwaki Co. Ltd., Japan). Copepods were cultivated at 20 • C, salinity 17 , and photoperiod of 12L:12D, initially using 5-L plastic beakers or 50-L polycarbonate tanks as a container. As the copepod population grew, the cultures were transferred to larger polycarbonate tanks (100 and 200 L). Aeration was gently performed via a glass tube at 1 ml min −1 for <50 L of medium, 3 ml min −1 for <100 L of medium, and 6 ml min −1 for <200 L of medium. Microalgae were supplied to copepods as a mixture every 2 days to achieve the following dry mass-based concentrations (Phaeodactylum, 31 µg ml −1 , Pavlova, 20 µg ml −1 , and Tisochrysis 12 µg ml −1 , at the ratio of 5:3:2), where cell concentrations were 400 × 10 3 , 200 × 10 3 , and 400 × 10 3 cells ml −1 , respectively. The medium for copepod culture was fully exchanged every 2-3 days. Generations of the copepods were repeatably incubated for a year before the present experiment, and the species was confirmed to be P. inopinus from the subsamples collected monthly. Except for the 200-L culture of copepods used as a diet for larval fish, three 50-L culture tanks were prepared for copepod fatty acid analysis. Protocols for copepod harvest, wash, and preservation were the same as those for rotifers.

Larval Rearing
Naturally fertilized and spawned red sea bream eggs were obtained from Ogata Suisan Inc. (Amakusa, Japan). The eggs were placed into black 100-L polyethylene tanks at a stocking density of 2,000 eggs per tank. The number of hatched larvae was estimated using a volumetric method described by Kotani et al. (2017), to be 1,900 individuals in each tank. Rearing tanks (100 L) were maintained under 12L:12D light conditions, with water quality monitored twice daily; pH and dissolved oxygen (DO) were measured by portable pH (D-51; Horiba Co. Ltd., Kyoto, Japan) and DO (FDO925; WTW, Bavaria, Germany) meters. Temperature was maintained at 21.8 ± 1.5 • C via a 1 kW titanium heater (Nittokizai Co, Ltd., Saitama, Japan); salinity was maintained at 33 ± 1 . Rearing tank aeration was initially supplied at 15 ml min −1 , and this gradually increased with larval growth. The first water exchange occurred 5 days post hatching (dph), after which water was exchanged at 3.5 ml s −1 via a siphoning system.
Larvae were reared until 20 dph following two prey conditions: SFC-enriched rotifers (control), and algal mixture-fed copepods (Figure 1). Both preys were sufficiently fed to fish larvae based on the dry mass, twice a day (at 08:00 and 16:00) from 3 dph until 19 dph. On critical days 3 and 4 post-hatching, to achieve a successful first feeding under copepod supply, the number of nauplii fed to fish larvae was increased to that of rotifers (ca. 5 inds ml −1 ). Larval rearing was terminated at 20 dph morning before feeding. Twenty fish larvae in each treatment were sampled at 0, 5, 10, 15, and 20 dph. All fish larvae surviving in control and copepod treatments at 20 dph morning were sampled, washed with distilled water, and separately pooled as 1 sample to detect their fatty acids. Larvae for fatty acid analysis were freeze-dried and preserved at −80 • C until lipid measurement.

Length and Dry Mass of Rotifers and Copepods
Lengths of rotifers and copepods were measured using a microruler under microscopy (SMZ800; Nikon Co, Ltd., Tokyo, Japan) after fixation in acidic Lugol's solution. The length of rotifer lorica (from crown spines to the lower extremity; Fu et al., 1991) was measured. Four categories of copepods were measured: nauplius (stages I-VI), copepodite (stages I-V), and adult male and female (copepodite stage VI). Adult males were distinguished from other stages by their asymmetrically developed first antennae; females were deemed adult when they possessed egg sacs. Body length was measured for nauplii (n = 40), and prosome length for copepodites (n = 30) and adults (n = 20).
Rotifer dry mass was measured by rinsing samples from three cultures with distilled water, freeze-dried, and weighed using a semi-micro balance (readability of 0.01 mg; GH-202; A and D, Japan). Dry mass of individual copepod was estimated from body length (µm) and prosome length (µm) using the equation adapted by Uye et al. (1983), with a factor (0.45) to convert carbon content to dry mass (Omori and Ikeda, 1984):

Larval Growth and Gut Contents
Fish standard length (SL) was measured after fixation in 10% formalin-seawater solution. Fifteen larvae sampled 2 h after feeding (5, 10, and 15 dph) were used for gut content analysis. In gut contents, rotifers and copepods at each developmental stage (nauplius, copepodite, and adult) were identified by microscopy (SMZ800; Nikon Co, Ltd., Tokyo, Japan) and counted. Total dry mass (µg larva −1 ) and DHA mass (ng larva −1 ) of ingested preys were calculated from the number of prey items in larval guts FIGURE 1 | Temporal change in prey items of larval red sea bream Pagrus major in experimental treatments (rotifers and copepods). Values in parentheses present densities based on individual number (ind. mL −1 ) and dry mass (µg mL −1 ). Rotifers were prepared by feeding docosahexaenoic acid-supplemented Super Fresh Chlorella (SFC); and the copepods fed algal mixture were used. Rotifer: Brachionus plicatilis sp. complex. Copepod: Pseudodiaptomus inopinus. multiplied by their dry mass (µg ind −1 ) and DHA content (mg g −1 ).

Fatty Acid Analysis
After homogenization of dried samples, total lipids were extracted following the method of Folch et al. (1957). Extracted lipids were separated into NLs and PLs fractions by column chromatography using a Sep-Pak silica cartridge (Waters, SA, United States; Juaneda and Rocquelin, 1985), with a chloroformmethanol mixture (98:2, v/v) for NLs and methanol for PLs. Methyl-esterification of fatty acids for total lipids, NLs, and PLs, was conducted as described by Matsui et al. (2019). Gas chromatography (GC-2010 Plus; Shimadzu, Kyoto, Japan) equipped with a hydrogen flame ionizing type detector (260 • C) and Omegawax capillary column (30 m length, 0.32 mm internal diameter, 0.25 µm film; Supelco, PA, United States) was used for detection of each fatty acid. Each fatty acid species was determined according to the equivalent chain length (ECL), and its proportion was calculated from the peak area. Quantification of fatty acid content was carried out based on the content of nonadecanoic acid (C19:0) analyzed together as an internal standard. Except for microalgae, rotifers, and copepods (n = 3), replicates for lipid analysis are not made for larval fish, due to difficulties in mass production of copepods used for later larviculture. Since fatty acid compositions for Phaeodactylum, Tisochrysis, and Pavlva microalgae have already been well reported (Volkman et al., 1989;Reitan et al., 1994a), each species was analyzed once (n = 1) and added to this study only as supplementary data for algal mixture.

Statistical Analysis
Statistical analysis was undertaken using Sigma-Plot ver 11.0 software (Systat Software Inc., CA, United States) and JMP version 11.0 software (SAS Institute Inc., NC, United States), p = 0.05. The length and dry mass of copepods at different stages were non-parametrically analyzed by Kruskal-Wallis test, followed by Steel-Dwass test, since Shapiro-Wilk's normality test and Bartlett's test both failed. A Student's t-test was conducted to compare total lipid fatty acid species of SFC and algal mixtures. Between rotifer fed on SFC and algal mixture, and copepods, the fatty acid proportions and contents for total lipids, NLs, and PLs were compared by one-way analysis of variance (ANOVA), then subjected to a Tukey pairwise multiple comparison post hoc test if significant differences were observed. An arcsine transformation was performed for proportions of fatty acids in each lipid before being subjected to a one-way ANOVA or Student's t-test. Larval growth rate was compared between feeding treatments (control and copepod) and between rearing times using twoway ANOVA. Total lipid fatty acid proportions of microalgae, rotifers, copepods, and larvae were analyzed using a principal component analysis (PCA) to characterize fatty acid profiles of each trophic level organism.

Length and Dry Mass of Copepod Stages
The length of copepods increased ontogenetically (Figure 2A). Nauplii showed similar dry mass but larger body length than those for rotifers (p < 0.05; Figure 2B). Length and dry mass were significantly higher in copepodites and adults than nauplii (p < 0.05). The length and dry mass were larger for adult females than adult males (p < 0.05).
those of SFC (p < 0.05). EPA was higher in the algal mixture (p < 0.05) because of the high proportion of fatty acids in P. lutheri (24.44%) and P. tricornutum (23.74%). While the algal mixture also contained DHA, mainly due to the fatty acids supplied from T. lutea (11.85%) and P. lutheri (9.62%), the DHA was lower than SFC (p < 0.05). The sum of the EPA and DHA (n-3 HUFAs) were occupied at approximately 20% of total fatty acids in both SFC and algal mixture. However, there was remarkable difference in DHA/EPA ratio between SFC (2.21) and algal mixture (0.35; p < 0.05). Arachidonic acid (ARA) was scarce (<1%) in all microalgae.

Fatty Acid Profiles of Rotifers and Copepods
The total lipid fatty acid composition of rotifers fed an algal mixture was compared with that of control diet (rotifers fed SFC; Table 2). The fatty acids C16:0, C18:1 (n-9 + n-7), and C18:2n-6 were high in total lipids of rotifers fed SFC. Rotifers fed SFC also exhibited EPA and DHA at 4.66 and 6.40%, respectively, resulting in DHA/EPA ratio of 1.38. When rotifers were fed an algal mixture, the C14:0, C16:1n-7, and C18:4n-3 fatty acid proportions increased significantly, and tended to have more EPA compared with those fed SFC (although the difference was not significant). On the contrary, the 18:2n-6 fatty acid was drastically decreased in rotifers fed an algal mixture (p < 0.05). The DHA was almost half as much as that of rotifers fed SFC. The fatty acid compositions of the two rotifer groups (those fed SFC, and those fed an algal mixture) corresponded with fatty acid peaks in SFC and the algal mixture ( Table 1).

Larviculture Performance of Red Sea Bream
Fish larvae fed rotifers and copepods exhibited exponential growth curves (p < 0.05; Figure 4A). SL at 20 dph was greater in fish larvae fed copepods (p < 0.05; Figure 4B). The growth coefficient of larvae fed copepods was higher than that for rotifer feeding treatment (p < 0.05; Table 5).
Based on gut contents at 5 dph (Figures 5A,B), fish larvae fed rotifers ingested an average of 2.0 ± 1.3 rotifers larva −1 , while those fed copepods ingested an average of 1.0 ± 1.2 nauplii larva −1 . Larval ingestion rates increased ontogenetically in both rotifer-and copepod-feeding treatments. Copepodites occurred in the guts of fish larvae from 10 dph. The total dry mass of ingested prey (µg larva −1 ) was estimated from the number of preys in the gut by multiplying each dry mass ( Figure 5C). Total dry mass ingested did not differ between fish larvae fed rotifers or copepods at 5 dph, although the difference in total dry matter was greater in fish larvae fed copepods at both 10 dph (1.7 ± 1.8 µg larva −1 ) and 15 dph (5.1 ± 3.6 µg larva −1 ; p < 0.05). Estimates of the total DHA mass ingested by fish larvae (ng larva −1 ) based on their dry mass-based contents in preys ( Figure 5D) showed that copepod-fed fish larvae ingested more PLs-DHA at 5 dph (0.77 ± 0.89 ng larva −1 ), 10 dph (5.85 ± 6.20 ng larva −1 ), and 15 dph (17.87 ± 12.47 ng larva −1 ) than fish larvae fed rotifers. The NLs-DHA mass ingested was comparable between fish larvae fed copepods and rotifers irrespective of their ages.
on PC1 were found for SFC and the scores of rotifers fed SFC and larvae fed those rotifers were close to the SFC scores. On the contrary, algal mixture had the most positive score on PC1. The scores of rotifers fed algal mixture were close to zero, and copepods fed algal mixture had positive scores between algal mixture and the rotifers. Larvae fed those copepods had relatively higher score than larvae fed SFC-enriched rotifers, although their scores were commonly negative. The PC2 mainly separated larvae from the lower trophic level organisms. Irrespective of the diets, larvae were commonly positioned at the most negative scores. On the contrary, SFC and algal mixture had the most positive scores. The scores of rotifers fed each diet little varied within the FIGURE 3 | Dry mass-based contents of eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) for non-polar lipids (NLs) and polar lipids (PLs) in enriched preys. SFC-Rotifer represents rotifers fed docosahexaenoic acid-supplemented Super Freshwater Chlorella (SFC). AM-Rotifer and AM-Copepod represent rotifers and copepods fed algal mixture, respectively. Values represent means ± standard deviation (n = 3). Letters above bars indicate significant differences (Tukey test, p < 0.05, a > b > c). positive side, unlike those for PC1. Compared with rotifers, copepods modified the PC2 scores of microalgae toward those of larvae.

Accumulation of DHA in PLs of Copepods
The result of PCA for total lipid fatty acid proportions of each trophic level revealed the followings; (1) fatty acid compositions of microalgae (algal mixture and SFC) were commonly the most different from those for larval fish, and (2) compared with microalgae, P. inopinus and rotifers had closer characteristics of fatty acids toward larvae based on DHA and the others (except for C16:0), respectively, (Figure 6). In particular, P. inopinus had PLs specifically enriched with DHA despite its diet containing low DHA relative to EPA (Tables 1, 4). Such selection of fatty acids to accumulate in the PLs can support the function of P. inopinus as a vector of DHA to larval fish. Fraser et al. (1989) reported wild-caught Calanus finmarchicus, Pseudocalanus sp., and Temora longicornis copepods to have abundant DHA in PLs (30.9, 31.8, and 31.9%, FIGURE 4 | Growth curve for Pagrus major larvae reared under different feeding regimes (A) and a histogram of larval standard length at the end of rearing periods (20 days post hatching) (B). The curve for larvae fed rotifers (SFC-Rotifer; solid line) and those fed copepods (broken line) is fitted based on data from 20 samples obtained at 0, 5, 10, 15, and 20 dph. respectively), which is consistent with our copepod species, P. inopinus (Table 4). Each copepod contained little DHA in triacylglycerol and wax esters, the main components of NLs (Fraser et al., 1989). It stresses the importance of DHA as a preferential component in the PLs of copepods. PLs, including phospholipids, mainly form lipid bilayer of plasma and organelle membranes in animals (Spector and Yorek, 1985). An abundance of DHA in C. finmarchicus eggs and nauplii was reported by Sargent and Falk-Petersen (1988). Egg production and hatching success for T. longicornis are enhanced with increased DHA concentration in the diet (Evjemo et al., 2008;Jónasdóttir et al., 2009). Considering specific abundance of DHA in PLs, the DHA in copepods probably functions as a key membrane component in the gonad and embryo during development.
On the other hand, proportions of DHA in PLs of artificially enriched rotifers were lower than in C. finmarchicus (Bell et al., 2003;Olsen et al., 2014), which is consistent with our comparison ( Table 4). Rotifers fed SFC had DHA enriched in NLs but not in PLs (Figure 3). There was no significant difference in PLs-DHA of rotifers fed SFC and the algal mixture irrespective of differences in their microalgae DHA (Table 1). Thus, rotifers possibly have PLs with relatively limited capacities for DHA (Li et al., 2015). Rotifers are reported to use DHA mainly as an energy source (Rainuzzo et al., 1994), which is a different DHA metabolic strategy to copepods.
In this study, it was unclear how P. inopinus could retain high concentrations of DHA in their PLs. In mammals and fish, phospholipase A2 and phospholipid acyltransferase function to modify fatty acid composition belonging to PLs (Sun et al., 1979;Balsinde, 2002;Tocher et al., 2008). Such enzymatic activities can also be investigated in copepods to reveal how DHA is preferentially accumulated in their PLs. In addition, copepods can biosynthesize DHA via fatty acid desaturation and elongation even when they cannot ingest diets including DHA (Monroig et al., 2013;Nielsen et al., 2019;Lee et al., 2020). In fact, DHA is detected in Pseudodiaptomus annandalei adults and nauplii fed non-DHA possessing microalga Tetraselmis chuii (Rayner et al., 2017). It indicates that both activities (i.e., accumulation of DHA in the PLs, and fatty acid biosynthesis) can support copepods to maintain DHA in their PLs.

Significance of Copepod Fatty Acid Profiles to Fish Larvae
Larvae fed rotifers and copepods had extremely high DHA relative to EPA (at the ratio of 4.50 and 3.10; Table 6), indicating the preference of DHA for larvae. In particular, PLs in larval fish fed both copepods and rotifers were commonly enriched with DHA regardless of dietary differences in total lipid DHA ( Table 2). DHA is highly concentrated in the larval fish brain; it is a key membrane component, and functions in development of the nervous system (Furuita et al., 1998;Masuda et al., 1999). In this case the DHA is often paired with C16:0 (Farkas et al., 2000), explaining why C16:0 was also presentative fatty acid for larvae (Figure 6). Because of the limited capacity of n-3 HUFA biosynthesis in marine fish (Tocher, 2003), red sea bream might preferentially retain DHA obtained from its diet in its PLs. In this study, although there are no replicates in the fatty acid analysis for larvae, NLs-DHA in larvae fed copepods is quite high and therefore non-negligible ( Table 6). Dietary DHA is efficiently incorporated into marine finfish larvae when the FIGURE 5 | Gut contents of larval Pagrus major based on the number of preys (A,B), total dry mass (C), and docosahexaenoic acid (DHA) mass (D). Ingested mass of DHA is presented for non-polar (NLs-DHA) and polar (PLs-DHA) lipid forms. Values represent means ± standard deviation (n = 15). Asterisks depict significant differences in dry and DHA masses between rotifers fed docosahexaenoic acid-supplemented Super Fresh Chlorella treatment (SFC-Rotifer) and copepod treatments (Student's t-test, p < 0.05).
fatty acid is formed in the PLs (Salhi et al., 1999;Gisbert et al., 2005;Wold et al., 2009). In fact, P. inopinus had high PLs-DHA, while rotifers fed SFC had high DHA in NLs but not in PLs (Figure 3). Larvae also ingested considerably more DHA mass from copepod PLs (Figure 5D). Our results suggest that accumulation of DHA in copepod PLs functions to efficiently transfer DHA to fish larvae.
rotifers and copepods. Larval Atlantic cod Gadus morhua growth was positively correlated with proportions of DHA (Cutts et al., 2006). Thus, higher DHA in P. inopinusfed larvae (especially in those NLs) can be linked to the enhanced growth.
As the other n-3 HUFA, EPA is essential fatty acid to further enhanced growth and survival of larvae (Ganuza et al., 2008;Eryalçın et al., 2013). In fact, an increase in concentrations of dietary n-3 HUFA resulted in higher larval growth as long as the DHA/EPA ratio is appropriate (Rodríguez et al., 1998). On the contrary, excess EPA over DHA (i.e., low DHA/EPA ratio) negatively affects larval growth and development (Reitan et al., 1994b;Rodríguez et al., 1997). P. inopinus, despite extremely low DHA/EPA in their diet (Table 1), increased that ratio of total lipids until larvae predated the preys ( Table 2). In particular, P. inopinus accumulated excess EPA in the NLs (Figure 4), which might be consumed as an energy source like NLs DHA for rotifers. They indicate that P. inopinus contribute to larval growth by transferring DHA and/or excluding EPA from microalgae.

Roles of P. inopinus in Coastal Regions in Transferring n-3 HUFA to Larvae
In the debate on fish recruitment, contribution of microalgae fatty acids to early development of larval fish has been discussed (Bell and Sargent, 1996). In marine ecosystems, specific fatty acids in microalgae are sufficiently transferred via mainly copepods to wild fish larvae so that they are used as biomarker for microalgae (Klungsøyr et al., 1989;Rossi et al., 2006). Especially in coastal and estuarine regions like P. inopinus habitats, microalgae diatoms are often dominated (Ubertini et al., 2012) and such blooms had higher proportions of EPA than DHA (Sargent et al., 1985). Sediments, which can be fed by the genus Pseudodiaptomus copepods (Jacoby and Greenwood, 1988), also include diatoms as a main FIGURE 6 | Principal component analysis (PCA) of total lipid fatty acid compositions of microalgae (circle), zooplankton (rotifers, diamond; copepods, cross), and larval fish (square). Green: docosahexaenoic acid-supplemented Super Fresh Chlorella (SFC). Brown: algal mixture. The PCA was calculated using proportions of C14:0, C16:0, C16:1, C16:2, C16:3, C18:2, C18:3, C18:4, EPA, C22:5, and DHA, which were significant differences of one or more percent between SFC and algal mixture. The PC1 scores were correlated with C14:0 (r = 0.913), C16:1 (r = 0.973), C16:2 (r = 0.708), C16:3 (r = 0.944), C18:2 (r = 0.860), C18:3 (r = 0.733), C18:4 (r = 0.909), EPA (r = 0.892), C22:5 (r = 0.892); and the PC2 ones were correlated with C16:0 (r = 0.783) and DHA (r = 0.742). component (mainly benthic microalgae; Grippo et al., 2010). However, larvae experienced on diatom-fed copepods showed lower growth than those on copepods fed dinoflagellate with higher DHA and lower EPA (St. John et al., 2001). They can hypothesize that fatty acid distribution in coastal microalgae is not suited for larval growth. According to the result of PCA on each fatty acid species for three trophic levels, however, two kinds of microalgae with distinct fatty acid compositions did not regulate larval DHA (Figure 6). Compared with rotifers, P. inopinus shifted algal fatty acid compositions much more toward DHA peak where larvae were positioned, and furthermore, enhanced larval growth (Figure 4). Our results suggest that what feeds on microalgae is more important factors than microalgae to determine larval growth. On the classical food chain in coastal waters, P. inopinus can function as a buffer of microalgae fatty acids for enhanced larval recruitment.

CONCLUSION
We artificially recreated a series of trophic levels in the coastal ecosystem using the cultured copepod P. inopinus, its prey (DHApoor algae), and predatory fish larvae. Fatty acid profiles of P. inopinus and fish larvae fed these copepods were compared with those fed DHA-enriched rotifers. Copepods upgraded DHA levels compared with the dietary fatty acids and selectively fortify their PLs. Concentrations of PLs-DHA were considerably higher in copepods compared with rotifers. Red sea bream larvae fed copepods had higher trends in NLs-DHA and enhanced the growth performance. We suggest that P. inopinus play an important ecological role as a carrier of the essential fatty acid DHA in PLs to coastal finfish larvae.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.

ETHICS STATEMENT
The animal study was reviewed and approved by Animal Experimental Committee, Kagoshima University.

AUTHOR CONTRIBUTIONS
HM, TS, and TKot designed the study and experiments. HM and TS conducted the experiments and collected the data. TKob identified copepod species. HM prepared a draft of the article and revised it with TKob and VW. MI provided the methodologies for GC analysis. HM performed data analyses. TKot supervised the studies. All the authors read and approved the final manuscript.

FUNDING
This work was supported by the Grant-in-Aid for Scientific Research (C) from Japan Society for the Promotion of Science (#24580275). Also, we received Functional Enhancement Expenses in Kagoshima University Operating Cost Grant "Advanced research infrastructure of biology and physiology of eel species inhabiting in Southern Kyushu, Japan, Taiwan, and South East Asia and Creation of Industry-academia-government network to induce the social implementation of those research result."