Vertical distribution of megafauna on the Bering Sea slope based on ROV survey

Video surveys were carried out during the 75th cruise of the RV Akademik M.A. Lavrentyev (June 2016) along the northern slope of the Volcanologists Massif, in the south-western Bering Sea. The seafloor was explored using the ROV Comanche 18. Seven dives were performed in the depth range from 4,278 m to 349 m. Overall, about 180 species of megafauna were recognised. Fifteen types of megafauna communities corresponding to certain depth ranges were distinguished based on the most abundant taxa. Dominance changed with depth in the following order: the holothurian Kolga kamchatica at the maximum depth (4,277–4,278 m); the holothurian Scotoplanes kurilensis at 3,610–2,790 m; the ophiuroid Ophiura bathybia at 3,030–2,910 m; benthic shrimps of the family Crangonidae at 2,910–2,290 m; the holothurian Paelopatides solea at 2,650–2,290 m; benthic jellyfish from the family Rhopalonematidae at 2,470–2,130 m; the enteropneust Torquaratoridae at 2,290–1,830 m; the holothurian Synallactes chuni and the ophiuroid of the genera Ophiura and Ophiocantha at 1,830–1,750 m. At depths 1,750–720 m most of the megafauna was associated with live or dead colonies of the sponge Farrea spp. Depths 720–390 m were dominated by the coral Heteropolypus ritteri and/or Corallimorphus pilatus. At 390–350 m depth, the shallowest depth range, the dominant taxon was the zoantharian Epizoanthus sp. Soft sediment megafauna communities dominated by torquaratorid enteropneusts to our knowledge have not been observed before in the deep-sea, the same as communities with a dominance of benthopelagic rhopalonematid jellyfish. The depths of the largest community changes, or the largest turnover of dominant species, were revealed at ∼2,790 m between the bathyal and abyssal zones and ∼1,750 m and ∼720 m within the bathyal zone.


INTRODUCTION
Depth plays a crucial role in the distribution of species and communities in the World Ocean. It was shown that species generally occupy discrete depth ranges and progressively replace each other from the continental shelf to the abyssal depths (Carney, Haedrich & Rowe, 1983). Distinct bathymetric faunal zones have been identified for many marine taxa (Mironov, 1986;Billett, 1991;Cartes & Sarda, 1993) and for marine fauna in general (Rowe & Menzies, 1969;Menzies, George & Rowe, 1973;Haedrich, Rowe & Polloni, 1975). Very few studies have examined the depth-related distribution of deep-sea species within their depth ranges to understand not only critical depths of faunal change but to identify more subtle changes that are related to the abundance of species that can provide a better overall picture of faunal zonation. For example , Howell, Billett & Tyler (2002) showed that asteroid species have very narrow centres of distribution in which they are abundant, despite much wider total adult depth ranges. Depth-related faunistic boundaries occur at many locations worldwide indicating that important controlling factors are present at these depths and that they may occur globally. Causes of the vertical zonation are related to the biology and physiology of the marine organisms and controlled by change with depth of some pivotal environmental parameters such as pressure, light, temperature and food availability (Menzies, George & Rowe, 1973;Rex, 1976;Carney, 2005;Rex & Etter, 2010). Generally, the causes of change in species composition with depth are complex, and several factors might interact, leading to the certain patterns observed in the each specific area (Soltwedel et al., 2009).
Benthic megafauna may be defined as marine animals exceeding 0.5-1 cm in size and recognizable on seafloor images (Grassle et al., 1975;Rex, 1981). Information on abundance and distribution of the megafauna is important for deep-sea benthic ecology. Megafaunal organisms play an important role in benthic ecosystems through active recycling of sediment organic matter, bioturbation and food web linkages.
The number of megafaunal studies has increased with the growing use of deep-sea video and imagery for ecological purposes.
Investigation of the Bering Sea began in the 18th century by G.W. Steller during the V. Bering expedition (1,737-1,742). Various species from the Bering Sea were studied in the late 19th and early 20th century. The first quantitative results of the benthos distribution were obtained in 1,932-1,933 by K. Derugin (Ushakov & Kusakin, 1978). Further considerable contributions were made by expeditions of the RV Vityaz in 1,950-1,952 to the western part of the Bering Sea and by the expeditions of a number of research vessels in 1,958-1,960 to the eastern part of the Bering Sea (Ushakov, 1952;Vinogradova, 1954;Neyman, 1960;Neiman, 1963;Zenkevitch, 1963;Filatova & Barsanova, 1964;Kuznetzov, 1964;Semenov, 1964;Vinogradov & Neyman, 1965). Qualitative and quantitative studies of the distribution of the benthos and trophic structure were described mainly on the shelf and upper slope (Rowland, 1973;Stoker, 1973;Alton, 1974;Haflinger, 1978;Stoker, 1978;Jewett & Feder, 1981) while surveys of the Western (Komandorsky) and Central (Aleutian) Basins were extremely limited. Impact of climate change on the Bering Sea ecosystems is being actively studied; however, these investigations are restricted only to some economically important species and do not take into account benthic communities as a whole (Hamazaki et al., 2005;Grebmeier et al., 2006;Stevenson & Lauth, 2012). Available data on the deep basins of the Bering Sea were obtained using trawls and grabs. Image and video surveys of the deep-sea megafauna communities were not conducted before, with the exception of dives in the submersible Mir in 1990 to the Piip's Volcano. Based on these dives, vertical zonation of benthic fauna on the upper part of the northern and southern slope of the Volcanologists Massif was described (Galkin & Sagalevich, 2012).
It was mentioned about 30 species of megafauna in the depth range from 450 to 900 m and the vertical zonation in fauna distribution was revealed. Dominant taxa changed with depth with hexactinellid sponges in the lower part replaced by corals and unidentified small cnidarians in the upper part.
The Bering Sea is separated from the North Pacific Ocean by the Aleutian-Komandor Island systems. The channels between the islands often have a great depth, sometimes exceeding 4,000 m, permitting almost unrestricted exchange between the abyssal zone of the Bering Sea and the Pacific Ocean (Stoker, 1978). As a result, most of the deep-sea species of the Bering Sea are in common with the Pacific Ocean (Filatova & Barsanova, 1964). By contrast, the exchange with the Chukchi Sea and the Arctic Ocean is limited to the Bering Strait, less than 50 m deep. The biomass of benthos in the Komandorsky Basin is very high-up to 32 g/m 2 , average 14 g/m 2 (Belyev, 1960). It is explained by favourable nutritional conditions due to proximity to the slope, abundant flora and fauna of the sublittoral zone and currents carrying organic-rich water from bays and mixing with the deep layers. For comparison, in the Central Basin biomass is lower due to isolation from the shelf (average 6 g/m 2 ) (Filatova & Barsanova, 1964).
At least seven physical factors may influence the qualitative and quantitative distribution of the Bering Sea benthic fauna. These factors include sediment particle size, bottom temperature, salinity, depth, sedimentation rates, circulation intensity, and suspended particle content in the near-bottom water layer (Kuznetzov, 1964;Stoker, 1978). Usually species dominating faunal assemblages belong to several trophic groups (Neyman, 1960;Filatova & Barsanova, 1964;Stoker, 1973); however, one certain trophic type can prevail. The distribution of different trophic types depends on topography, distribution of organic carbon in the sediment, mobility of bottom waters and distribution of erosion and sedimented zones at the seafloor (Neiman, 1963;Rex, 1981;Gage & Tyler, 1991;Sokolova, 2000). For example, sessile suspension feeders are dominated in the Bering Sea at areas with high speed of currents, at areas with high concentrations of suspension particles in near bottom layers of water and on hard substrate (Kuznetzov, 1964). Deposit-feeders are dominated at areas with soft sediment enriched with organic carbon (Gage & Tyler, 1991).
The present study is the first detailed investigation on the megafauna communities of the slope of the Volcanologists Massif and adjacent area of the Komandorsky Basin based on video survey. Our main aim is to clarify patterns of the vertical distribution of megafauna communities on the slope of the Volcanologists Massif based on their composition and structure.

Study area and video surveys
Video surveys were carried out during the 75th cruise of the RV Akademik M. A. Lavrentyev (June 2016) in the northern slope of the Volcanologists Massif, the south-west Bering Sea (Fig. 1). The seafloor was explored using the ROV Comanche 18 (Galkin & Ivin, 2019). Seven dives were performed in the depth range from 4,278 m to 349 m. For technical reasons, dives tracks were directed upwards along the slope. No specific permissions were  Table 1.
Full-size DOI: 10.7717/peerj.8628/ fig-1 required for the investigated location because it is the territorial waters of the Russian Federation and it is not the protected area.
Visual observations were conducted in a continuous mode, total duration is 40 hrs 44 min. Observations were accompanied by video recording and sampling of individual organisms by ROV manipulator. The field studies did not involve endangered or protected species. The ROV was equipped with the Kongsberg Underwater HDTV colour camera OE14-502 and with a ten-centimetre laser scale. A sonardyne USBL Fusion underwater navigation system was used. In total 22 hrs 29 min of video footage were analysed. Dive details are given in Table 1.
The present investigation was focused on soft sediment epibenthic megafauna. To estimate the relative abundance of the soft sediment megafauna, a series of video transects were performed at depths 4,278-1,370 m. Transects for quantitative analysis were selected . Sediment cover on these transects exceeded 70% of the total seafloor square. The number of transects selected for the analysis from each dive depended on the quality of video records. The quality in turn was subject to conditions of dive (speed and altitude) influenced by the seafloor topography, currents and technical aspects. Transects were selected if they complied with the following characteristics: camera was looking almost vertically down, average speed was about 0.5 knots, altitude was approximately 1.5 m, there were no sediment clouds and illumination was satisfactory (Table 1). Altogether 33 quantitative transects were analysed with total time of 2 hrs 25 min. The time of individual transects varied from 1.5 min to 8 min (about 5 min in most cases). The fauna of the hard substrate on all dives was described only qualitatively due to difficulties in adequate identification to a low taxonomic level and counting of animals on hard substrate with complex topology. On two dives (1 and 9) at depths 1,370-350 m, quantitative transects could not be selected owing to the predominance of hard substrata and extremely rough topography. However, records from these transects were used for qualitative analysis of benthic communities.

Video analysis
Video records were analysed on land using the QuickTime software. Visible megafauna was registered and taxonomic experts provided identification of the observed organisms to the lowest possible taxonomic level. For identification of some species, specimens were collected using the ROV manipulator. Complete lists of taxa identified on video and in collected material are given in Table S1. Images of identified taxa on video are given in Data S1. Association with hard or soft substrate was noted for all identified species. Abundance of all identified species on hard and soft substrate was evaluated visually using four-point scale based on video (Table S1).
All specimens of each visible soft sediment megafauna species on quantitative transects were counted (Table S2). The following taxa/organisms were excluded from the statistical analysis because they could not be adequately identified and/or counted on video: highly mobile taxa (isopods, polychaetes Macellicephalinae, amphipods, shrimps, fishes), infauna represented only by Lebensspuren, gelatinous zooplankton, small-size organisms (<one cm) and organisms that could not be identified at least to the phylum level. For shrimps and jellyfishes, exceptions were made for near bottom species numerous in some areas.
Because the densities of some species on soft sediment were relatively high it was most optimal to extract individual images (print screen) from the quantitative transects with an interval of 30 s to determine densities. In total 365 digital images were examined using the Photoshop software. A ten-centimetre laser scale was used for calculation of the seafloor surface area on images.

Data analysis
Statistical analysis was performed in Primer V6 (Clarke & Warwick, 2001).
In order to discern whether depth was driving changes in community composition a cluster analysis was performed based on quantitative transects on the soft substrate. Each quantitative transect was treated as a separate sample. The relative contribution (%) of each species to total abundance on each transect was calculated. Relative abundance data were square-root transformed to reduce the dominance of the most abundant species. The similarity between transects was estimated using the quantitative Bray-Curtis index. Similarity matrixes was used for cluster analyses. Clusters were generated using group averaged linkage. The resulting divisions were tested using SIMPROF permutation tests, for looking statistically significant evidence of genuine clusters in samples which are a priori unstructured. Revealed upper and lower boundaries of depth ranges based on the analysis of quantitative transects were adjusted based on the visual observations of dominance of the habitat-forming taxa. It was done because there were significant sections of the video that didn't fall on the quantitative transects. Average density per depth range (±standard deviation) of the most abundant soft sediment taxa was calculated in order to quantitatively describe differences between detected zones.
Boundaries of the depth ranges within 1,370-350 m with >90% hard substrate were estimated visually based on video observations of the dominant taxa. At the depths 1,750-1,370 m with 70% hard substrate we separately analysed areas of the soft substrate that were used in the cluster analysis and areas with a hard substrate at which dominant taxa was evaluated only visually based on video observations.
The bathymetric distribution of dominant species was then plotted to visually determine the approximate depths where the largest turnover of dominant species occurred. These are the depths where species dominant over a wide depth range are replaced by other dominants. This information was then used to assess the number of identified species, aggregated into several major taxonomic groups, that occur between determined depths.

RESULTS
Overall about 170 species of megafauna were recognised on images. 30% of these species were collected using ROV and identified in the laboratory by taxonomic experts (Table S1). Two species of Hexactinellida (genus Farrea), three species of Actiniaria (genus Sicyonis), one species of Corallimorpharia (Corallimorphidae gen. sp.) and one species of Holothuroidea (Zygothuria sp.1) collected using ROV appeared to be new to science and they will be described separately. In addition, three species are assumed to also be new to science based on images and/or video records (one species each of ascidian, antipatharian and benthic siphonophorae), but these cannot be verified without physical samples.

The vertical distribution of communities
The non-metric cluster analysis revealed ten distinct groups of soft-sediment megabenthic communities distinguished by depth of occurrence (Fig. 2) at the similarity level of 56% (SIMPROF: R = 7,438, p = 0.001). Transects that cannot be significantly differentiated based on SIMPROF tests were connected on Fig. 2 by red dotted lines. In two cases transects differentiated based on SIMPROF tests were located in a very close depth range, which may indicate possible other factors affecting their community structure.
We distinguished separately the depth zone 1,750-1,370 m with 70% hard substrate where sections of the soft and hard substrate were analysed independently. This zone was Characteristics of communities in each depth zone (dominant species on the soft substrate and common species on hard and soft substrate) are given in Table 2 (upper and lower boundaries of depth ranges are given with adjusted based on the visual observations of dominance of the habitat-forming taxa in the whole video footage).
The dominant taxa changed with depth in the following order: the holothurian Kolga kamchatica Rogacheva, 2012 at the maximum depth (4,277-4,278 m); the holothurian Scotoplanes kurilensis Gebruk, 1983 was dominant or abundant at depths 3,610-2,790 m; the ophiuroid Ophiura bathybia Clark, 1911 sharply dominated at depths 3,030-2,910 m; benthic crangonid shrimps dominated at 2,910-2,290 m with a peak at 2,790-2,470 m; the holothurian Paelopatides solea Baranova, 1955 at 2,650-2,290 m; benthic jellyfish at 2,470-2,130 m; enteropneusts at 2,290-1,830 m; the holothurian Synallactes chuni and ophiuroid of the genera Ophiura and Ophiocantha at 1,830-1,370 m; the glass sponge Farrea spp. at 1,750-1,370 m. At the depths 1,370-720 m almost all benthic epifauna was associated with live or dead colonies of Farrea spp., except of areas of the soft substrate (10% of the total area) dominated by the ophiuroids Ophiura and Ophiocantha. Depths 720-390 m were dominated by the coral Heteropolypus ritteri Nutting, 1909. The corallimorpharian   At the minimum in our study depth 390-350 m, the zoantharian Epizoanthus sp. was the dominant taxon. Contribution to the abundance (in %) per depth zone of the most abundant taxa, based on analysis of quantitative transects, is presented in Fig. 3 and Table 3. Densities of the most abundant taxa are given in Table 3. Characteristic images of all distinguished depth zones with dominant species are presented in Figs. 4 and 5.

The depths with the largest community changes
Based on the vertical distribution of dominant megafaunal taxa, the depths of the largest community changes, or the largest turnover of dominant species, were revealed (Fig. 6). The first pronounced change occurred at ∼2,790 m depth and was identified by the replacement of the holothurian Scotoplanes kurilensis by other species of holothurians, enteropneusts, benthic jellyfish and shrimps. The second pronounced change occurred at ∼1,750 m: here the glass sponge Farrea spp. took the role of the landscape determining species. The third change at ∼720 m was associated with the replacement of glass sponges by soft corals. We also suggest that a community change may have occurred at depths between 3,650-4,200 m (but a lack of observations between 4,277-3,650 m limited more precise zoning). This change was indicated by the replacement of the holothurian Kolga kamchatica by the holothurian Scotoplanes kurilensis. The number of identified species, aggregated into several major taxonomic groups, that occur between depths with the largest community changes is presented in Fig. 7. The highest number of species was found at depths 3,650-2,790 m, the lowest-at 4,277-4,278 m. The taxon with the highest species richness at the maximum depth 4,277-4,278 m was Holothuroidea (eight species). The relatively high number of holothurian species was also at 3,650-2,790 m (nine species). The Order Pennatulacea was recorded only at 3,650-2,790 m (two species). The highest diversity of actiniarians was at depths 2,790-1,750 Table 3 The mean densities (ind m −2 ) ± standard deviations and the contribution to total abundance (in %) of the most abundant taxa in different depth zones. The dominant taxa are shaded grey.  m (14 species). Alcyonacean corals showed the highest species richness at 1,750-720 m (not less than 9 species). The highest species richness of hexactinellid sponges was at 1,750-720 m (not less than ten species).

H e t e r o p o l y p u s r i t t e r i C o r a l l i m o r p h u s p i l a t u s E p i z o a n t h u s s p .
Depth, m dominant species common species rare species

DISCUSSION
In the present study, the vertical distribution of benthic megafaunal communities was described quantitatively for the first time on the slope of the Volcanologists Massif, in the south-west Bering Sea and adjacent Komandorsky Basin in the depth range 4,278-349 m. We were able to identify approximately 170 species, based on the photo and video surveys together with samples taken by ROV. Fifteen depth zones were revealed. Seven species collected using ROV appeared new to science and will be described separately. In addition, three species are assumed to be new to science based on images and/or video records, but further research and sampling using ROV in this area are required for taxonomic clarification.

Soft sediment communities
The greatest depths of the Komandorsky Basin (4,278 m) were dominated by holothurians (eight species) with the prevailing of the elpidiid Kolga kamchatica. Aggregations of small elpidiid holothurians Kolga in the abyssal zone is a phenomenon observed in several areas of the ocean. For example, Kolga nana Théel, 1879 aggregations were found in the Charlie-Gibbs Fracture Zone (Gebruk & Krylova, 2013) with densities of up to 76 ind/m 2 and in the Porcupine Seabight in the North Atlantic (Billett & Hansen, 1982) with a mean density of 34 ind/m 2 at depths at around 4,000 m. Specimens of K. nana were observed in the Charlie-Gibbs Fracture Zone feeding on the soft sediment mainly in flat and gentle slope areas (Rogacheva, 2012). Another species of this genus, K. hyalina Danielssen & Koren, 1879, was abundant in the Central Arctic Basins with densities significantly dependent on seasonal ice-algae supply to the seafloor (Rybakova et al., 2019). The density of K. kamchatica in our study (17 ind/m 2 ) was lower than in the mentioned reports, although it was still relatively high. Another elpidiid holothurian, Scotoplanes kurilensis, with lower mean density 1-4 ind/m 2 , was abundant at depths 3,610-2,790 m. In general, elpidiid holothurians can reach high abundances, especially in seafloor depressions in the abyssal zone. Rowe (1971) regarded aggregations of elpidiid holothurians as a typical feature of canyons which trap organic matter that attracts deposit-feeders. However, aggregations of elpidiid holothurians are not exclusive to seafloor depressions: small elpidiids commonly aggregate in areas of accumulation of organic matter on the seafloor (Gebruk, 1990;Billett, 1991;Roberts et al., 2000;Gebruk et al., 2013). Amaro et al. (2010) showed on the example of infaunal molpadiid holothurian Molpadia musculus Risso, 1826 from the Nazare Canyon (NE Atlantic) the key functional role of holothurians in the degradation of organic matter from the sediments and its redistribution and availability for other fauna. Generally, depths greater 3,500 m in our study were characterized by a high diversity of holothurians (11 species), with one species usually dominant at any given depth. Depositfeeding holothurians are known to have different feeding strategies and activity (Roberts et al., 2000). Feeding-guilds and niche separation optimize the use of food resources. Studies in various regions showed that holothurians, especially elpidiids, are able to react rapidly to the changes in food supply (Billett et al., 2001;Billet et al., 2010;Gutt et al., 2011;Kuhnz et al., 2014;Huffard et al., 2016).
At 3,610-3,450 m depth, a high diversity of suspension feeders associated with soft substrate was notable, including the hexactinellid sponge Holascus sp., some unidentified species of ceriantharians and two species of Pennatulacea. It can be suggested that this is related to peculiarities of near-bottom hydrodynamics in this area. The growth of the relative importance of sessile suspension feeders (sponges and ceriantharians) deeper than 3,000 m was mentioned by Kuznetzov (1964) based on limited data from the Central and Komandorsky Basins of the Bering Sea: these species developed specific adaptations to the soft-clay habitats, such as raising their bodies above the sediment and the ability to utilize food from the near bottom water layers.
The ophiuroid Ophiura bathybia formed aggregations at 3,030-2,910 m depth with a very high mean density (74 ind/m 2 ). Many ophiuroid taxa are predators (Warner, 1982). Perhaps local sedimentation or hydrological conditions may affect food concentration for them in this area.
At 2,790-2,650 m depth the species diversity and abundance of megafauna were low. Unidentified shrimps of the family Crangonidae were the most notable and abundant softsediment megafauna. The hard substrate together with the rare synallactid holothurians and corals at this depth were covered with a considerable sediment layer. We hypothesize that this was an effect of active re-deposition or rapid sedimentation causing decrease in abundance and diversity of megafauna.
On the lower slope, at depths from 3,000 m to 2,150 m, the synallactid holothurian, Paelopatides solea, played a significant role. Numerically they were not dominant, however owing to big size, their contribution to the biomass must have been significant. It means that they may play an important role in ecosystem functioning at that depth range utilizing organic matter from the sediment.
Communities with numerical dominance of benthopelagic rhopalonematid jellyfish at depths >1,500 m (2,470-2,130 m in the present study) to our knowledge have not been reported before. Traditional methods of benthos investigations (trawls and grabs) do not allow the detection of aggregations of gelatinous organisms at the seafloor, because such organisms are usually destroyed in catches. Visual seafloor inspection can resolve this problem. Different observations at depths of up to one and a half thousand meters suggest that deep-sea medusae can be abundant near the bottom (Larson et al., 1992;Matsumoto, Baxter & Chen, 1997;Grange et al., 2017). It can be suggested that aggregations of jellyfish at the seafloor are related to the local near-bottom current regime, apparently increasing food availability through resuspension of organic matter (Larson et al., 1992). Soft sediment megafauna communities dominated by torquaratorid enteropneusts (at the depths 2,130-1,830 m) have not been reported before in the deep-sea either. Deep-sea enteropneusts are also usually gelatinous and easily destroyed in trawl and grab catches. With increased use of ROVs, enteropneusts have been recorded in different deep ocean basins . Deep-sea enteropneusts are still poorly known, however their function in benthic ecosystems can be important, especially in organic matter processing and surficial bioturbation (Jones et al., 2013). Deep-water enteropneusts feed on the sediment surface, collecting detritus from the substrate, apparently with little or no selectivity as they crawl forward leaving behind a faecal trail of undigestible presumably mineralised material (Holland et al., 2009). Some species have demonstrated the ability to use near-bottom currents to move between feeding grounds in a controlled manner through changes in body posture Jones et al., 2013). Estimation of the area of their traces in the northern part of the Mid-Atlantic Ridge at around 2,500 m depth revealed that enteropneusts contribute significantly to the surficial deposit feeding (Jones et al., 2013). This was unexpected because local deposit feeding megafaunal density was dominated by holothurians, with only a relatively small proportion of enteropneusts (Jones et al., 2013). According to our observations, the mean density of enteropneusts was unusually high reaching 12 ind m −2 at depths 2,130-1,830 m . The highest density ever previously recorded (for a mixture of Australian torquaratorid species) was 29 worms per 1,000 m 2 (Anderson, Przeslawski & Tran, 2011). The highest density found during a 15-year study of Tergivelum baldwinae Holland et al., 2009 from the Pacific Ocean was two worms per 1,000 m 2 (Smith, Holland & Ruhl, 2005). In our study area, this taxon definitely must play a major role in the utilization of organic matter at the seafloor.
Small areas of soft sediment (10% of the total area) at depths 1,830-1,370 m were inhabited mainly by the holothurians Synallactes chuni and Pannychia aff. moseleyi Théel, 1882 and some species of ophiuroids. Notably, these species of holothurians and ophiuroids occurred also on hard substrate in that area which may indicate that they are ubiquitous, with no specific adaptations to soft sediment.

Hard substrate communities
The bottom of the Komandorsky Basin was 100% sediment covered, with rare occurrences of colonies of the hexactinellid sponge, Farrea sp. 1, with associated fauna. From the base of the slope up to 1,830 m depth, the area covered by sediment smoothly decreased to 50%, while the area of hard substrate increased. Accordingly, the diversity of species living on rocky outcrops increased, including sessile forms of corals, actiniarians, sponges, ascidians and others. Hard substrate provides more heterogeneous habitat as compared to soft sediment, additionally some organisms use sessile megafaunal species (e.g., sponges, corals) as substrates for attachment (e.g., ophiuroids living on gorgonian corals) (Levin et al., 2001;BuhlMortensen et al., 2010). Some corals and sponges modify habitats through their physical presence. They provide living space, alter physical conditions and affect biological interactions, considerably enhancing diversity (Felley, Vecchione & Wilson, 2008;Miller et al., 2012;Gebruk & Krylova, 2013).
At 1,830 to 1,370 m depth, the area of the hard substrate increased up to 90%. Rocky outcrops alternated with small areas of soft sediment. The deposit feeding holothurians, Synallactes chuni and Pannychia aff. moseleyi, and some species of ophiuroids occurred together with highly diverse sessile suspension feeding fauna on rocks covered by a thin layer of sediment. The ratio of sessile suspension feeders is known to increase on rocks owing to increased active water circulation (Kuznetzov, 1964;Witman & Dayton, 2001;Tyler, 2003;Miller & Etter, 2008;Miller & Etter, 2011). Starting from a depth of 1,750 m, the hard substrates were covered by live and dead colonies of the hexactinellid sponge, Farrea spp., with associated fauna. At 720 m this community was sharply replaced by an assemblage of cnidarians, including Heteropolypus ritteri, Corallimorphus pilatus, Epizoanthus sp. and various actiniarians. From a depth of about 1,750 m, the slope became steeper, possibly affecting water circulation and providing a habitat favourable for suspension feeders. The exact reasons of replacement of sponges by cnidarians remain unknown.
It is important to note that, based on the video or photo images, species could not be equally identified across all taxa, thus limiting the comparison of species diversity at different depths.

Generalized patterns of the vertical zonation
Depths with the largest community changes identified in our study as depths of the largest turnover of dominant taxa, in general correspond to the patterns of vertical zonation recognized by other authors (Carney, Haedrich & Rowe, 1983;Mironov, 1986;Howell, Billett & Tyler, 2002). Major faunal changes occurred at the transition between bathyal and abyssal zones, at ∼2,790 m depth. At this depth, the abyssal community dominated by the deposit-feeding holothurian Scotoplanes was replaced by the bathyal/slope communities of different feeding types (with several dominant species characterized by different feeding types) that may reflect changes in the quality and quantity of food resources at the bottom. The boundary between the abyssal and bathyal zones at a depth of ∼3,000 m was commonly noticed in previous works, however some authors recognized this boundary closer to 2,000 m (Gage & Tyler, 1991) or 4,000 m (Thurman, 1985). In the scheme of the vertical zonation suggested by Belyaev et al. (1959), the bathyal is defined as a depth zone between 200 m and 3,000 m, with a zone between 2,500 m and 3,500 m considered as the transition between the bathyal and abyssal zones. The boundary between the abyssal and bathyal zones can vary around the world based among other factors on geomorphological features of the seabed (Menzies, George & Rowe, 1973). Also, the depth at which this boundary stands out significantly depends on the methodology used. Within the abyssal zone we also identified considerable community change between the Scotoplanes dominated and Kolga dominated communities at depths between 3,650 m and 4,277 m.
The zoning of the bathyal zone also differs in published literature. Many authors distinguish the lower, middle and upper bathyal/slope based on the distribution of species or communities, however the depth boundaries of these horizons differ (Howell, Billett & Tyler, 2002). In the present study we distinguished a variety of communities dominated by soft or hard substrate fauna and deposit-feeding or suspension-feeding species at depths from ∼2,790 m to ∼1,750 m which may correspond to a lower bathyal zone. The mid bathyal zone we recognized between ∼1,750 m and ∼720 m, dominated by suspensionfeeding hexactinellid sponges associated with the increased availability of hard substrate. The upper bathyal zone appeared to occur between ∼720 m and 350 m, dominated by suspension-feeding cnidarians on hard substrate. In distinguishing the sub-zones within the bathyal zone, we followed the approach of Belyaev et al. (1959).
The relatively high species richness in our study occurred between the depths of 3,650-720 compared to depth ranges of 4,277-4,278 m and 720-350 m. The peak in species richness occurs at different depths in different World Ocean regions (Gebruk, Budaeva & King, 2010). Our results are concordant with several other authors who observed a peak at around 1,500-3,000 m (Rex, 1981;Howell, Billett & Tyler, 2002;Rex & Etter, 2010).
Thus, the vertical distribution of megafauna communities on the northern slope of the Volcanologists Massif (in the south-western Bering Sea) based on the video survey was examined for the first time. We revealed a distinct pattern of vertical zonation. Fifteen types of megafaunal communities, corelated with particular depth ranges, were distinguished. The largest community changes occurred at depths of ∼2,790 m, ∼1,750 m and ∼720 m. These depths correspond to transitions from the abyssal to the bathyal zone, and to transitions between the lower, middle and upper bathyal zones, as shown in a range of publications by different authors. Our results contribute towards better understanding of the regional patterns of vertical faunal zonation in the World Ocean.