Sanitary quality of edible bivalve mollusks in Southeastern Brazil using an UV based depuration system

The increase in seafood production, especially in mariculture worldwide, has brought out the need of continued monitoring of shellfish production areas in order to ensure safety to human consumption. The purpose of this research was to evaluate pathogenic protozoa, viruses and bacteria contamination in oysters before and after UV depuration procedure, in brackish waters at all stages of cultivation and treatment steps and to enumerate microbiological indicators of fecal contamination from production site up to depuration site in an oyster cooperative located at the Southeastern estuarine area of Brazil. Oysters and brackish water were collected monthly from September 2009 to November 2010. Four sampling sites were selected for enteropathogens analysis: site 1- oyster growth, site 2- catchment water (before UV depuration procedure), site 3 - filtration stage of water treatment (only for protozoa analysis) and site 4- oyster's depuration tank. Three microbiological indicators ! were examined at sites 1, 2 and 4. The following pathogenic microorganisms were searched: Giardia cysts, Cryptosporidium oocysts, Human Adenovirus (HAdV), Hepatitis A virus (HAV), Human Norovirus (HnoV) (genogroups I and II), JC strain Polyomavirus (JCPyV) and Salmonella sp. Analysis consisted of molecular detection (qPCR) for viruses (oysters and water samples); immunomagnetic separation followed by direct immunofluorescence assay for Cryptosporidium oocysts and Giardia cysts and also molecular detection (PCR) for the latter (oysters and water samples); commercial kit (Reveal-Neogee (R)) for Salmonella analysis (oysters). Giardia was the most prevalent pathogen in all sites where it was detected: 36.3%, 18.1%, 36.3% and 27.2% of water from sites 1, 2, 3 and 4 respectively; 36.3% of oysters from site 1 and 54.5% of depurated oysters were harboring Giardia cysts. The huge majority of contaminated samples were classified as Giardia duodenalis. HAdv was detected in water and o! ysters from growth site and HnoV GI in two batches of oysters ! (site 1) in huge concentrations (2.11 x 10(13), 3.10 x 10(12) gc/g). In depuration tank site, Salmonella sp., HAV (4.84 x 10(3)) and HnoV GII (7.97 x 10(14)) were detected once in different batches of oysters. Cryptosporidium spp. oocysts were present in 9.0% of water samples from site four. These results reflect the contamination of oysters even when UV depuration procedures are employed in this shellfish treatment plant. Moreover, the molecular comprehension of the sources of contamination is necessary to develop an efficient management strategy allied to shellfish treatment improvement to prevent foodborne illnesses. (C) 2011 Elsevier Ltd. All rights reserved.

A molecular phylogeny of bivalve mollusks: ancient radiations and divergences as revealed by mitochondrial genes

The main scope of my PhD is the reconstruction of the large-scale bivalve phylogeny on the basis of four mitochondrial genes, with samples taken from all major groups of the class. To my knowledge, it is the first attempt of such a breadth in Bivalvia. I decided to focus on both ribosomal and protein coding DNA sequences (two ribosomal encoding genes -12s and 16s -, and two protein coding ones - cytochrome c oxidase I and cytochrome b), since either bibliography and my preliminary results confirmed the importance of combined gene signals in improving evolutionary pathways of the group. Moreover, I wanted to propose a methodological pipeline that proved to be useful to obtain robust results in bivalves phylogeny. Actually, best-performing taxon sampling and alignment strategies were tested, and several data partitioning and molecular evolution models were analyzed, thus demonstrating the importance of molding and implementing non-trivial evolutionary models. In the line of a more rigorous approach to data analysis, I also proposed a new method to assess taxon sampling, by developing Clarke and Warwick statistics: taxon sampling is a major concern in phylogenetic studies, and incomplete, biased, or improper taxon assemblies can lead to misleading results in reconstructing evolutionary trees. Theoretical methods are already available to optimize taxon choice in phylogenetic analyses, but most involve some knowledge about genetic relationships of the group of interest, or even a well-established phylogeny itself; these data are not always available in general phylogenetic applications. The method I proposed measures the "phylogenetic representativeness" of a given sample or set of samples and it is based entirely on the pre-existing available taxonomy of the ingroup, which is commonly known to investigators. Moreover, it also accounts for instability and discordance in taxonomies. A Python-based script suite, called PhyRe, has been developed to implement all analyses.

Evaluation of Microbial Contamination in Bivalve Mollusks: Epidemiology and Diagnosis

Shellfish are filter-feeding organisms that can accumulate many bacteria and viruses. Considering that depuration procedures are not effective in removal of certain microorganisms, shellfish-borne diseases are frequent in many parts of the world, and their control must rely primarily on investigation of prevalence of human pathogens in shellfish and water environment. However, the diffusion of enteric viruses and Vibrio bacteria is not known in many geographical areas, for example in Sardinia, Italy. A survey aimed at investigating the prevalence of Norovirus (NoV), hepatitis A virus (HAV), V. parahaemolyticus, V. cholerae and V. vulnificus was carried out, analyzing both local and imported purified, non-purified and retail shellfish from North Italy and Sardinia. Shellfish from both areas were found contaminated by NoVs, HAV and Vibrio, including retail and purified animals. Molecular analysis evidenced different NoV genogroups and genotypes, including bovine NoVs, as well as pathogenic Vibrio strains, underlining the risk for shellfish consumers. However, also other approaches are needed to control the diffusion of shellfish-borne diseases. It was originally thought that enteric viruses are passively accumulated by shellfish. Recently, it was proven that NoVs bind to specific carbohydrate ligands in oysters, and various NoV strains are characterized by a different bioaccumulation pattern. To deepen the knowledge on this argument, a study was carried out, analyzing bioaccumulation of up to 8 different NoV strains in four different species of shellfish. Different bioaccumulation patterns were observed for each shellfish species and NoV strain used, potentially important in setting up effective shellfish purification protocols. Finally, a novel study of evaluation of viral contamination in shellfish from the French Atlantic coast was carried out following the passage of Xynthia tempest over Western Europe which caused massive destruction. Different enteric viruses were found over a one month period, evidencing the potential of these events of contaminating shellfish.

Bivalve mollusks in the sediments of the Anapa Bay Bar (Black Sea) in October 2010

Sandy beaches of the Anapa Bay Bar are a unique natural resource, but they are gradually being degrade under both natural and anthropogenic factors. Emissions of sand and shelly ground from the adjacent sea bottom partly compensate for this process. Concentration of carbonates may reach up to 50% in beach sands, and most of these carbonates are of mollusk origin. The major deposit formation role belongs to the key bivalve species: Chamelea gallina (Linnaeus, 1758). Average biomass of this mollusk species reaches up to 450 g/m**2 at depths 5-10 m. The other two subdominating mollusk species, bivalve Donax trunculus (Linnaeus, 1758) and gastropod Rapana venosa (Valenciennes, 1846), may impact as 16 g/m**2 and 6 g/m**2, respectively. Annually, 350 kg of shelly ground per running meter are newly deposited on the Anapa beach.

The influence of temperature and seawater carbonate saturation state on 13C-18O bond ordering in bivalve mollusks

The shells of marine mollusks are widely used archives of past climate and ocean chemistry. Whilst the measurement of mollusk delta 18O to develop records of past climate change is a commonly used approach, it has proven challenging to develop reliable independent paleothermometers that can be used to deconvolve the contributions of temperature and fluid composition on molluscan oxygen isotope compositions. Here we investigate the temperature dependence of 13C-18O bond abundance, denoted by the measured parameter Delta 47, in shell carbonates of bivalve mollusks and assess its potential to be a useful paleothermometer. We report measurements on cultured specimens spanning a range in water temperatures of 5 to 25 °C, and field collected specimens spanning a range of -1 to 29 °C. In addition we investigate the potential influence of carbonate saturation state on bivalve stable isotope compositions by making measurements on both calcitic and aragonitic specimens that have been cultured in seawater that is either supersaturated or undersaturated with respect to aragonite. We find a robust relationship between Delta 47 and growth temperature. We also find that the slope of a linear regression through all the Delta 47 data for bivalves plotted against seawater temperature is significantly shallower than previously published inorganic and biogenic carbonate calibration studies produced in our laboratory and go on to discuss the possible sources of this difference. We find that changing seawater saturation state does not have significant effect on the Delta 47 of bivalve shell carbonate in two taxa that we examined, and we do not observe significant differences between Delta 47-temperature relationships between calcitic and aragonitic taxa.

Evaluation of tropical water sources and mollusks in southern Brazil using microbiological, biochemical, and chemical parameters

Florianopolis, a city located in the Santa Catarina State in southern Brazil, is the national leading producer of bivalve mollusks. The quality of bivalve mollusks is closely related to the sanitary conditions of surrounding waters where they are cultivated. Presently, cultivation areas receive large amounts of effluents derived mainly from treated and non-treated domestic, rural, and urban sewage. This contributes to the contamination of mollusks with trace metals, pesticides, other organic compounds, and human pathogens such as viruses, bacteria, and protozoan. The aim of this study was to perform a thorough diagnosis of the shellfish growing areas in Florianopolis, on the coast of Santa Catarina. The contamination levels of seawater, sediments, and oysters were evaluated for their microbiological, biochemical, and chemical parameters at five sea sites in Florianopolis, namely three regular oyster cultivation areas (Sites 1, 2, and oyster supplier), a polluted site (Site 3), and a heavily polluted site (Site 4). Samples were evaluated at day zero and after 14 days. Seawater and sediment samples were collected just once, at the end of the experiment. Antioxidant defenses, which may occur in contaminated environments in response to the increased production of reactive oxygen species (ROS) by organisms, were analyzed in oysters, as well as organic compounds (in oysters and sediment samples) and microbiological contamination (in oysters and seawater samples). The results showed the presence of the following contaminants: fecal coliforms in seawater samples (four sites), human adenovirus (all sites), human noroviruses GI and GII (two sites), Hepatitis A viruses (one site), JC Polyomavirus in an oyster sample from the oyster supplier, Giardia duodenalis cysts, and Cryptosporidium sp oocysts (one site). Among organochlorine pesticides, only DDT (dichlorodiphenyltrichloroethane) and HCH (hexachlorocyclohexane) were detected in some sediment and oysters samples in very low levels; site 4 had the highest concentrations of total aliphatic hydrocarbons. PAHs, and linear alkylbenzenes (LABs) found either in oysters or in sediment samples. The major concentration of fecal sterol coprostanol was found at site 4, followed by site 3. After 14 days of allocation in the four selected sites, there was a significant difference in the enzymes analyzed at the monitored spots. The detection of different contaminants in oysters, seawater, and sediment samples in the present study shows the impact untreated or inadequately treated effluents have on coastal areas. These results highlight the need for public investment in adequate wastewater treatment and adequate treatment of oysters, ensuring safe areas for shellfish production as well as healthier bivalve mollusks for consumption.

Calibration study of bivalve shells : implications for paleoenvironmental reconstruction

The accretionary shells of bivalve mollusks can provide environmental information, such as water temperature, precipitation, freshwater fluxes, primary productivity and anthropogenic activities in the form of variable growth rates and variable geochemical properties, such as stable oxygen and carbon isotopes. However, paleoenvironmental reconstructions are constrained by uncertainties about isotopic equilibrium fractionation during shell formation, which is generally acknowledged as a reasonable assumption for bivalves, but it has been disputed in several species. Furthermore, the variation in shell growth rates is accepted to rely on multiple environmental variables, such as temperature, food availability and salinity, but can differ from species to species. Therefore, it is necessary to perform species-specific calibration studies for both isotope proxies and shell growth rates before they can be used with confidence for environmental interpretations of the past. Accordingly, the principal objective of this Ph.D research is to examine the reliability of selected bivalve species, the long-lived Eurhomalea exalbida (Dillwyn), the short-lived and fast growing species Paphia undulata (Born 1778), and the freshwater mussel Margaritifera falcata (Gould 1850), as paleoenvironmental proxy archives.rnThe first part is focused on δ18Oshell and shell growth history of live-collected E. exalbida from the Falkland Islands. The most remarkable finding, however, is that E. exalbida formed its shell with an offset of -0.48‰ to -1.91‰ from the expected oxygen isotopic equilibrium with the ambient water. If this remained unnoticed, paleotemperature estimates would overestimate actual water temperatures by 2.1-8.3°C. With increasing ontogenetic age, the discrepancy between measured and reconstructed temperatures increased exponentially, irrespective of the seasonally varying shell growth rates. This study clearly demonstrates that, when the disequilibrium fractionation effect is taken into account, E. exalbida can serve as a high-resolution paleoclimate archive for the southern South America. The species therefore provides quantifiable temperature estimates, which yields new insights into long-term paleoclimate dynamics for mid to high latitudes on the southern hemisphere.rnThe stable carbon isotope of biogenic carbonates is generally considered to be useful for reconstruction of seawater dissolved inorganic carbon. The δ13Cshell composition of E. exalbida was therefore, investigated in the second part of this study. This chapter focuses on inter-annual and intra-annual variations in δ13Cshell. Environmental records in δ13Cshell are found to be strongly obscured by changes in shell growth rates, even if removing the ontogenetic decreasing trend. This suggests that δ13Cshell in E. exalbida may not be useful as an environmental proxy, but a potential tool for ecological investigations. rnIn addition to long-lived bivalve species, short-lived species that secrete their shells extremely fast, can also be useful for environmental reconstructions, especially as a high-resolution recorder. Therefore, P. undulata from Daya Bay, South China Sea was utilized in Chapter 4 to evaluate and establish a potential proxy archive for past variations of the East Asian monsoon on shorter time-scales. The δ18Oshell can provide qualitative estimates of the amount of monsoonal rain and terrestrial runoff and the δ13Cshell likely reflect the relative amount of isotopically light terrestrial carbon that reaches the ocean during the summer monsoon season. Therefore, shells of P. undulata can provide serviceable proxy archives to reconstruct the frequency of exceptional summer monsoons in the past. The relative strength of monsoon-related precipitation and associated changes in ocean salinity and the δ13C ratios of the dissolved inorganic carbon signature (δ13CDIC) can be estimated from the δ18Oshell and δ13Cshell values as well as shell growth patterns. rnIn the final part, the freshwater pearl shell M. falcata from four rivers in British Columbia, Canada was preliminarily studied concerning the lifespans and the shell growth rates. Two groups separated by the Georgia Strait can be clearly distinguished. Specimens from the western group exhibit a shorter lifespan, while the eastern group live longer. Moreover, the average lifespan seems to decrease from south to north. The computed growth equations from the eastern and western groups differ as well. The western group exhibits a lower growth rate, while bivalves from the eastern group grow faster. The land use history seems to be responsible for the differences in lifespans of the specimens from the two groups. Differences in growth rate may be induced by differences in water temperature or nutrient input also related to the land use activities.

Occurrence and distribution of Vibrio parahaemolyticus in retail oysters in Sao Paulo State, Brazil

Vibrio parahaemolyticus is a potentially pathogenic bacterium that occurs naturally in estuarine environments worldwide, and is often associated with gastroenteritis in humans following consumption of raw bivalve mollusks, especially raw oysters. The occurrence of total and pathogenic V. parahaemolyticus in 74 samples of raw oysters collected in restaurants, supermarkets, groceries and beach huts in Sao Paulo State, was monitored between February 2006 and January 2007. Enumeration of V. parahaemolyticus was performed according to the most probable number (MPN) procedure. Five to ten typical colonies were selected from thiosulfate-citrate-bile salts-sucrose (TCBS) agar plates for confirmation by the presence of the species-specific gene tlh and the virulence genes tdh and trh by multiplex PCR. V. parahaemolyticus was detected in 100% of samples. The densities of total V. parahaemolyticus varied from 1.78 to 6.04 logio (MPN/g), with higher densities being detected in fall and summer, and lower densities in winter (P < 0.05). There was no statistical difference among densities of V parahaemolyticus regarding the site of collection. None of the 1943 V parahaemolyticus isolates contained tdh and/or trh. These data provide information for the assessment of exposure to V. parahaemolyticus in oysters consumed in Sao Paulo, State, Brazil. (C) 2010 Elsevier Ltd. All rights reserved.

A trophic analysis of target species of macrobenthos in a subtropical coast community: a taxa relationship essay

Studies on the feeding habits of aquatic organisms are a requirement for the management and sustainable use of marine ecosystems. The aim of the present research was to analyze the habits and trophic similarities of decapods, starfish and fish in order to propose trophic relationships between taxa, using Hennigian methods of phylogenetic systematics. This new grouping hypothesis, based on shared and exclusive food items and food types, corresponds to the broad taxonomic groups used in the analysis. Our results indicate that algae, Mollusca, Polychaeta, Crustacea, Echinodermata and Actinopterygii are the most exploited common resources among the species studied. Starfish were differentiated from other organisms for being stenophagic, and were grouped for feeding on bivalve mollusks. A larger group of fish and crustaceans shares algae and mainly crustaceans as food items. A third group united all eight species of Actinopterygii. This largest subgroup of fish is typically carnivorous, feeding on Anthozoa and a great quantity of Crustacea. Synodus foetens has a special position among fishes, due to its unique feeding on nematodes. A Euclidean distance dendrogram obtained in a previous publication grouped S. foetens with starfish. That result was based on a few non-exclusive shared similarities in feeding modes, as well as on shared absences of items, which are not an adequate grouping factor. Starfish are stenophagic, eating bivalves almost exclusively. Synodus foetens and Isopisthus parvipinnis have restricted food items, and are thus intermediary in relation to starfish, decapods, and other fish, which are euryphagous. The trophic cladogram displays details of food items, whether or not shared by all species. The resulting trophic analysis is consistent with known historical relationships.

Ocurrence of Vibrio spp., positive coagulase staphylococci and enteric bacteria in oysters (Crassostrea gigas) harvested in the south bay of Santa Catarina island, Brazil

The aim of this study was to assess the contamination of oysters (Crassostrea gigas), harvested in six different regions of the South Bay of Santa Catarina Island, with Coliforms at 45 ºC, Escherichia coli, Vibrio spp., positive coagulase staphylococci, and Salmonella sp. over a period of one year. One hundred eighty oyster samples were collected directly from their culture sites and analyzed. Each sample consisted of a pool of 12 oysters. All of the samples analyzed showed absence of Salmonella, 18 (10%) samples showed presence of Escherichia coli, 15 (8.3%) samples were positive for V. alginolyticus, and Vibriocholerae was detected in 4 samples (2.2%). The counts of positive-coagulase staphylococci varied from <10 to 1.9 x 102 CFU.g-1, whereas the counts of Coliforms at 45 ºC and E. coli ranged from <3 to 1.5 x 102 MPN.g-1 and <3 and 4.3 x 10 MPN.g-1, respectively. Counts of V. parahaemolyticus and V. vulnificus ranged between <3 and 7 MPN.g-1, for both microorganisms. This suggests the need for monitoring these Vibrios contamination in oysters. Based on the results of the microbiological assays, the samples analyzed showed acceptable bacteriological quality, i.e., they were within the parameters established by Brazilian Legislation.

Salinity influences glutathione S-transferase activity and lipid peroxidation responses in the Crassostrea gigas oyster exposed to diesel oil

Biochemical responses in bivalve mollusks are commonly employed in environmental studies as biomarkers of aquatic contamination. The present study evaluated the possible influence of salinity (35, 25,15 and 9 ppt) in the biomarker responses of Crassostrea gigas oysters exposed to diesel at different nominal concentrations (0.01, 0.1 and 1 mLL(-1)) using a semi-static exposure system. Salinity alone did not resulted in major changes in the gill`s catalase activity (CAT), glutathione S-transferase activity (GST) and lipid peroxidation levels (measured as malondialdehyde. MDA), but influenced diesel related responses. At 25 ppt salinity, but not at the other salinity levels, oysters exposed to diesel showed a strikingly positive concentration-dependent GST response. At 25 ppt and 1 mLL(-1) diesel, the GST activity in the gills remained elevated, even after one week of depuration in clean water. The increased MDA levels in the oysters exposed to diesel comparing to control groups at 9, 15 and 35 ppt salinities suggest the occurrence of lipid peroxidation in those salinities, but not at 25 ppt salinity. The MDA quickly returned to basal levels after 24 h of depuration. CAT activity was unaltered by the treatments employed. High toxicity for 1 mLL(-1) diesel was observed only at 35 ppt salinity, but not in the other salinities. Results from this study strongly suggest that salinity influences the diesel related biomarker responses and toxicity in C. gigas, and that some of those responses remain altered even after depuration. (C) 2011 Elsevier B.V. All rights reserved.

Shell beds as paleoecological puzzles: A case study from the Upper Permian of the Parana Basin, Brazil

Microstratigraphic, sedimentological, and taphonomic features of the Ferraz Shell Bed, from the Upper Permian (Kazanian-Tatarian?) Corumbatai Formation of Rio Claro Region (the Parana Basin, Brazil), indicate that the bed consists of four distinct microstratigraphic units. They include, from bottom to top, a lag concentration (Unit 1), a partly reworked storm deposit (Unit 2), a rapidly deposited sandstone unit with three thin horizons recording episodes of reworking (Unit 3), and a shell-rich horizon generated by reworking/winnowing that was subsequently buried by storm-induced obrution deposit (Unit 4). The bioclasts of the Ferraz Shell Bed represent exclusively bivalve mollusks. Pinzonella illusa and Terraia aequilateralis are the dominant species. Taphonomic analysis indicates that mollusks are heavily time-averaged (except for some parts of Unit 3). Moreover, different species are time-averaged to a different degree (disharmonious time-averaging). The units differ statistically from one another in their taxonomic and ecological composition, in their taphonomic pattern, and in the size-frequency distributions of the two most common species. Other Permian shell beds of the Parana Basin are similar to the Ferraz Shell Bed in their faunal composition (they typically contain similar sets of 5 to 10 bivalve species) and in their taphonomic, sedimentologic, and microstratigraphic characteristics. However, rare shell beds that include 2-3 species only and are dominated by articulated shells preserved in life position also occur. Diversity levels in the Permian benthic associations of the Parana Basin were very low, with the point diversity of 2-3 species and with the within-habitat and basin-wide (alpha and gamma) diversities of 10 species, at most. The Parana Basin benthic communities may have thus been analogous to low-diversity bivalve-dominated associations of the present-day Baltic Sea. The 'Ferraz-type' shell beds of the Parana Basin represent genetically complex and highly heterogeneous sources of paleontological data. They are cumulative records of spectra of benthic ecosystems time-averaged over long periods of time (10(2)-10(4) years judging from actualistic research). Detailed biostratinomic reconstructions of shell beds can not only offer useful insights into their depositional histories, but may also allow paleoecologists to optimize their sampling designs, and consequently, refine paleoecological and paleoenvironmental interpretations.

Salinity influences glutathione S-transferase activity and lipid peroxidation responses in the Crassostrea gigas oyster exposed to diesel oil

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Drill holes in shells of Permian benthic invertebrates

Newly discovered benthic fossils and specimens illustrated in the paleontological literature indicate that drilling predators (or parasites) were present in the Permian. New field data from southern Brazil document the first drill holes ever reported for Permian bivalve mollusks. In addition, a literature review revealed drill holes in shells of articulate brachiopods from Russia, Greece, and West Texas. Holes range in size from 0.1 to 5.8 mm and are typically round, cylindrical, singular penetrations perpendicular to the valve surface. Incomplete, healed, and multiple holes are absent. Drilling frequency, a proxy for predation intensity, is very low: less than 1 percent (this estimate may be seriously affected by taphonomic and monographic biases). Literature data suggest that frequency of drilled specimens varied significantly among higher brachiopod taxa. The geography and stratigraphy of drilled specimens indicate that drilling organisms were worldwide in their occurrence and continuously present in marine ecosystems throughout the Permian. This report is consistent with other recent studies indicating that although drillers were continuously present throughout the Phanerozoic, drilling intensity was lower in the Late Paleozoic and early Mesozoic.

Comparative analysis of drilling frequencies in recent brachiopod-mollusk associations from the Southern Brazilian Shelf

Over 14,000 specimens-5,204 brachiopods, 9,137 bivalves, and 178 gastropods-acquired from 30 collecting stations (0 to 45 m depth) in the Ubatuba and Picinguaba bays, southern Brazil, were compared for drilling frequencies. Beveled (countersunk) circular-to-subcircular borings (Oichnus-like drill holes) were found in diverse bivalves but also in the rhynchonelliform brachiopod Bouchardia rosea-a small, semi-infaunal to epifaunal, free-lying species that dominates the brachiopod fauna of the southern Brazilian shelf. Drill holes in bivalve mollusks and brachiopods are comparable in their morphology, average diameter, and diameter range, indicating attacks by a single type of drilling organism. Drill holes in brachiopods were rare (0.4%) and found only at five sampling sites. Drillings in bivalves were over 10 times as frequent as in brachiopods, but the average drilling frequency was still low (5.6%) compared to typical boring frequencies of Cenozoic mollusks. Some common bivalve species, however, were drilled at frequencies up to 50 times higher than those observed for shells of B. rosea from the same samples. Due to scarcity of drilled brachiopods, it is not possible to evaluate if the driller displayed a nonrandom (stereotyped) site, size, or valve preference. Drilled brachiopods may record (1) naticid or muricid predation, (2) predation by other drillers, (3) parasitic drillings, and (4) mistaken or opportunistic attacks. Low drilling frequency in brachiopods is consistent with recent reports on ancient and modern examples. The scarcity of drilling in brachiopods, coupled with much higher drilling frequencies observed in sympatric bivalves, suggests that drilling in brachiopods may have been due to facultative or erroneous attacks. The drilling frequencies observed here for the brachiopod-bivalve assemblages are remarkably similar to those reported for Permian brachiopod-bivalves associations. This report adds to the growing evidence for an intriguing macroecological stasis: multiple meta-analytical surveys of present-day and fossil rhynchonelliform brachiopods conducted in recent years also point to persistent scarcity and low intensity of biotic interactions between brachiopods and drilling organisms throughout their evolutionary history.