Effects of petrogenic pollutants on North Atlantic and Arctic Calanus copepods: From molecular mechanisms to population impacts.

Oil and gas industries in the Northern Atlantic Ocean have gradually moved closer to the Arctic areas, a process expected to be further facilitated by sea ice withdrawal caused by global warming. Copepods of the genus Calanus hold a key position in these cold-water food webs, providing an important energetic link between primary production and higher trophic levels. Due to their ecological importance, there is a concern about how accidental oil spills and produced water discharges may impact cold-water copepods. In this review, we summarize the current knowledge of the toxicity of petroleum on North Atlantic and Arctic Calanus copepods. We also review how recent development of high-quality transcriptomes from RNA-sequencing of copepods have identified genes regulating key biological processes, like molting, diapause and reproduction in Calanus copepods, to suggest linkages between exposure, molecular mechanisms and effects on higher levels of biological organization. We found that the available ecotoxicity threshold data for these copepods provide valuable information about their sensitivity to acute petrogenic exposures; however, there is still insufficient knowledge regarding underlying mechanisms of toxicity and the potential for long-term implications of relevance for copepod ecology and phenology. Copepod transcriptomics has expanded our understanding of how key biological processes are regulated in cold-water copepods. These advances can improve our understanding of how pollutants affect biological processes, and thus provide the basis for new knowledge frameworks spanning the effect continuum from molecular initiating events to adverse effects of regulatory relevance. Such efforts, guided by concepts such as adverse outcome pathways (AOPs), enable standardized and transparent characterization and evaluation of knowledge and identifies research gaps and priorities. This review suggests enhancing mechanistic understanding of exposure-effect relationships to better understand and link biomarker responses to adverse effects to improve risk assessments assessing ecological effects of pollutant mixtures, like crude oil, in Arctic areas.

De novo transcriptomes of six calanoid copepods (Crustacea): a resource for the discovery of novel genes.

This study presents eight new high-quality de novo transcriptomes from six co-occurring species of calanoid copepods, the first published for Neocalanus plumchrus, N. cristatus, Eucalanus bungii and Metridia pacifica and additional ones for N. flemingeri and Calanus marshallae. They are ecologically-important members of sub-arctic North Pacific marine zooplankton communities. 'Omics data for this diverse and numerous taxonomic group are sparse and difficult to obtain. Total RNA from single individuals was used to construct gene libraries that were sequenced on an Illumina Next-Seq platform. Quality filtered reads were assembled with Trinity software and validated using multiple criteria. The study's primary purpose is to provide a resource for gene expression studies. The integrated database can be used for quantitative inter- and intra-species comparisons of gene expression patterns across biological processes. An example of an additional use is provided for discovering novel and evolutionarily-significant proteins within the Calanoida. A workflow was designed to find and characterize unannotated transcripts with homologies across de novo assemblies that have also been shown to be eco-responsive.

Evidence for virus-mediated oncogenesis in bladder cancers arising in solid organ transplant recipients.

A small percentage of bladder cancers in the general population have been found to harbor DNA viruses. In contrast, up to 25% of tumors of solid organ transplant recipients, who are at an increased risk of developing bladder cancer and have an overall poorer outcomes, harbor BK polyomavirus (BKPyV). To better understand the biology of the tumors and the mechanisms of carcinogenesis from potential oncoviruses, we performed whole genome and transcriptome sequencing on bladder cancer specimens from 43 transplant patients. Nearly half of the tumors from this patient population contained viral sequences. The most common were from BKPyV (N=9, 21%), JC polyomavirus (N=7, 16%), carcinogenic human papillomaviruses (N=3, 7%), and torque teno viruses (N=5, 12%). Immunohistochemistry revealed variable Large T antigen expression in BKPyV-positive tumors ranging from 100% positive staining of tumor tissue to less than 1%. In most cases of BKPyV-positive tumors, the viral genome appeared to be clonally integrated into the host chromosome consistent with microhomology-mediated end joining and coincided with focal amplifications of the tumor genome similar to other virus-mediated cancers. Significant changes in host gene expression consistent with the functions of BKPyV Large T antigen were also observed in these tumors. Lastly, we identified four mutation signatures in our cases, with those attributable to APOBEC3 and SBS5 being the most abundant. Mutation signatures associated with an antiviral drug, ganciclovir, and aristolochic acid, a nephrotoxic compound found in some herbal medicines, were also observed. The results suggest multiple pathways to carcinogenesis in solid organ transplant recipients with a large fraction being virus-associated.

Larval fish counteract ram and suction to capture evasive prey.

A simple hydrodynamic model of predator-prey interactions between larval clownfish and copepod prey is used to elucidate how larval fish capture highly evasive copepods. Fish larvae are considered to be suction feeders; however, video observations revealed that successful captures by clownfish larvae were preceded by rapidly accelerating lunges (ram), while the role of suction to draw prey into the fish's mouth was less clear. Simulations were made of the fish's strike, varying strengths of ram and suction to characterize optimal strategies for copepod capture given known evasive capabilities. Our results suggest that, contrary to expectations, suction feeding is dominant only in older larvae, whereas ram feeding is the dominant mode for early larvae. Despite the relatively weak suction produced by smaller larvae, it still plays a crucial role in prey capture through hydrodynamic stealth. Escape-triggering water deformations from the strike can be cancelled through controlled suction. Experimental data obtained from larval clownfish agree with model results, suggesting that the primary role of suction in early larvae is providing hydrodynamic stealth rather than capture.

Physiological acclimatization in high-latitude zooplankton.

How individual organisms adapt to nonoptimal conditions through physiological acclimatization is central to predicting the consequences of unusual abiotic and biotic conditions such as those produced by marine heat waves. The Northeast Pacific, including the Gulf of Alaska, experienced an extreme warming event (2014-2016, "The Blob") that affected all trophic levels and led to large-scale changes in the community. The marine copepod Neocalanus flemingeri is a key member of the subarctic Pacific pelagic ecosystem. During the spring phytoplankton bloom this copepod builds substantial lipid stores as it prepares for its nonfeeding adult phase. A 3-year comparison of gene expression profiles of copepods collected in Prince William Sound in the Gulf of Alaska between 2015 and 2017 included two high-temperature years (2015 and 2016) and one year with very low phytoplankton abundances (2016). The largest differences in gene expression were between high and low chlorophyll years, and not between warm and cool years. The observed gene expression patterns were indicative of physiological acclimatization. The predominant signal in 2016 was the down-regulation of genes involved in glycolysis and its incoming pathways, consistent with the modulation of metabolic rates in response to prolonged low food conditions. Despite the down-regulation of genes involved in metabolism, there was no evidence of suppression of protein synthesis based on gene expression or behavioural activity. Genes involved in muscle function were up-regulated, and the copepods were actively swimming and responsive to stimuli at collection. However, genes involved in fatty acid metabolism were down-regulated in 2016, suggesting reduced lipid accumulation.

Characterizing the tumor microenvironment in rare renal cancer histological types.

The tumor microenvironment (TME), including immune cells, cancer-associated fibroblasts, endothelial cells, adjacent normal cells, and others, plays a crucial role in influencing tumor behavior and progression. Here, we characterized the TME in 83 primary renal tumors and matched metastatic or recurrence tissue samples (n = 15) from papillary renal cell carcinoma (pRCC) types 1 (n = 20) and 2 (n = 49), collecting duct carcinomas (CDC; n = 14), and high-grade urothelial carcinomas (HGUC; n = 5). We investigated 10 different markers of immune infiltration, vasculature, cell proliferation, and epithelial-to-mesenchymal transition by using machine learning image analysis in conjunction with immunohistochemistry. Marker expression was compared by Mann-Whitney and Kruskal-Wallis tests and correlations across markers using Spearman's rank correlation coefficient. Multivariable Poisson regression analysis was used to compare marker expression between histological types, while accounting for variation in tissue size. Several immune markers showed different rates of expression across histological types of renal carcinoma. Using pRCC1 as reference, the incidence rate ratio (IRR) of CD3+ T cells (IRR [95% confidence interval, CI] = 2.48 [1.53-4.01]) and CD20+ B cells (IRR [95% CI] = 4.38 [1.22-5.58]) was statistically significantly higher in CDC. In contrast, CD68+ macrophages predominated in pRCC1 (IRR [95% CI] = 2.35 [1.42-3.9]). Spatial analysis revealed CD3+ T-cell and CD20+ B-cell expressions in CDC to be higher at the proximal (p < 0.0001) and distal (p < 0.0001) tumor periphery than within the central tumor core. In contrast, expression of CD68+ macrophages in pRCC2 was higher in the tumor center compared to the proximal (p = 0.0451) tumor periphery and pRCC1 showed a distance-dependent reduction, from the central tumor, in CD68+ macrophages with the lowest expression of CD68 marker at the distal tumor periphery (p = 0.004). This study provides novel insights into the TME of rare kidney cancer types, which are often understudied. Our findings of differences in marker expression and localization by histological subtype could have implications for tumor progression and response to immunotherapies or other targeted therapies.

Mammary collagen architecture and its association with mammographic density and lesion severity among women undergoing image-guided breast biopsy.

BACKGROUND: Elevated mammographic breast density is a strong breast cancer risk factor with poorly understood etiology. Increased deposition of collagen, one of the main fibrous proteins present in breast stroma, has been associated with increased mammographic density. Collagen fiber architecture has been linked to poor outcomes in breast cancer. However, relationships of quantitative collagen fiber features assessed in diagnostic biopsies with mammographic density and lesion severity are not well-established. METHODS: Clinically indicated breast biopsies from 65 in situ or invasive breast cancer cases and 73 frequency matched-controls with a benign biopsy result were used to measure collagen fiber features (length, straightness, width, alignment, orientation and density (fibers/µm2)) using second harmonic generation microscopy in up to three regions of interest (ROIs) per biopsy: normal, benign breast disease, and cancer. Local and global mammographic density volumes were quantified in the ipsilateral breast in pre-biopsy full-field digital mammograms. Associations of fibrillar collagen features with mammographic density and severity of biopsy diagnosis were evaluated using generalized estimating equation models with an independent correlation structure to account for multiple ROIs within each biopsy section. RESULTS: Collagen fiber density was positively associated with the proportion of stroma on the biopsy slide (p < 0.001) and with local percent mammographic density volume at both the biopsy target (p = 0.035) and within a 2 mm perilesional ring (p = 0.02), but not with global mammographic density measures. As severity of the breast biopsy diagnosis increased at the ROI level, collagen fibers tended to be less dense, shorter, straighter, thinner, and more aligned with one another (p < 0.05). CONCLUSIONS: Collagen fiber density was positively associated with local, but not global, mammographic density, suggesting that collagen microarchitecture may not translate into macroscopic mammographic features. However, collagen fiber features may be markers of cancer risk and/or progression among women referred for biopsy based on abnormal breast imaging.

Genomic and evolutionary classification of lung cancer in never smokers.

Lung cancer in never smokers (LCINS) is a common cause of cancer mortality but its genomic landscape is poorly characterized. Here high-coverage whole-genome sequencing of 232 LCINS showed 3 subtypes defined by copy number aberrations. The dominant subtype (piano), which is rare in lung cancer in smokers, features somatic UBA1 mutations, germline AR variants and stem cell-like properties, including low mutational burden, high intratumor heterogeneity, long telomeres, frequent KRAS mutations and slow growth, as suggested by the occurrence of cancer drivers' progenitor cells many years before tumor diagnosis. The other subtypes are characterized by specific amplifications and EGFR mutations (mezzo-forte) and whole-genome doubling (forte). No strong tobacco smoking signatures were detected, even in cases with exposure to secondhand tobacco smoke. Genes within the receptor tyrosine kinase-Ras pathway had distinct impacts on survival; five genomic alterations independently doubled mortality. These findings create avenues for personalized treatment in LCINS.

Post-diapause transcriptomic restarts: insight from a high-latitude copepod.

BACKGROUND: Diapause is a seasonal dormancy that allows organisms to survive unfavorable conditions and optimizes the timing of reproduction and growth. Emergence from diapause reverses the state of arrested development and metabolic suppression returning the organism to an active state. The physiological mechanisms that regulate the transition from diapause to post-diapause are still unknown. In this study, this transition has been characterized for the sub-arctic calanoid copepod Neocalanus flemingeri, a key crustacean zooplankter that supports the highly productive North Pacific fisheries. Transcriptional profiling of females, determined over a two-week time series starting with diapausing females collected from > 400 m depth, characterized the molecular mechanisms that regulate the post-diapause trajectory. RESULTS: A complex set of transitions in relative gene expression defined the transcriptomic changes from diapause to post-diapause. Despite low temperatures (5-6 °C), the switch from a "diapause" to a "post-diapause" transcriptional profile occurred within 12 h of the termination stimulus. Transcriptional changes signaling the end of diapause were activated within one-hour post collection and included the up-regulation of genes involved in the 20E cascade pathway, the TCA cycle and RNA metabolism in combination with the down-regulation of genes associated with chromatin silencing. By 12 h, females exhibited a post-diapause phenotype characterized by the up-regulation of genes involved in cell division, cell differentiation and multiple developmental processes. By seven days post collection, the reproductive program was fully activated as indicated by up-regulation of genes involved in oogenesis and energy metabolism, processes that were enriched among the differentially expressed genes. CONCLUSIONS: The analysis revealed a finely structured, precisely orchestrated sequence of transcriptional changes that led to rapid changes in the activation of biological processes paving the way to the successful completion of the reproductive program. Our findings lead to new hypotheses related to potentially universal mechanisms that terminate diapause before an organism can resume its developmental program.

Diapause vs. reproductive programs: transcriptional phenotypes in a keystone copepod.

Many arthropods undergo a seasonal dormancy termed "diapause" to optimize timing of reproduction in highly seasonal environments. In the North Atlantic, the copepod Calanus finmarchicus completes one to three generations annually with some individuals maturing into adults, while others interrupt their development to enter diapause. It is unknown which, why and when individuals enter the diapause program. Transcriptomic data from copepods on known programs were analyzed using dimensionality reduction of gene expression and functional analyses to identify program-specific genes and biological processes. These analyses elucidated physiological differences and established protocols that distinguish between programs. Differences in gene expression were associated with maturation of individuals on the reproductive program, while those on the diapause program showed little change over time. Only two of six filters effectively separated copepods by developmental program. The first one included all genes annotated to RNA metabolism and this was confirmed using differential gene expression analysis. The second filter identified 54 differentially expressed genes that were consistently up-regulated in individuals on the diapause program in comparison with those on the reproductive program. Annotated to oogenesis, RNA metabolism and fatty acid biosynthesis, these genes are both indicators for diapause preparation and good candidates for functional studies.

Diesel exhaust and bladder cancer risk by pathologic stage and grade subtypes.

BACKGROUND: The International Agency for Research on Cancer (IARC) classifies diesel engine exhaust as carcinogenic to humans based on sufficient evidence for lung cancer. IARC noted, however, an increased risk of bladder cancer (based on limited evidence). OBJECTIVE: To evaluate the association between quantitative, lifetime occupational diesel exhaust exposure and risk of urothelial cell carcinoma of the bladder (UBC) overall and according to pathological subtypes. METHODS: Data from personal interviews with 1944 UBC cases, as well as formalin-fixed paraffin-embedded tumor tissue blocks, and 2135 controls were pooled from two case-control studies conducted in the U.S. and Spain. Lifetime occupational histories combined with exposure-oriented questions were used to estimate cumulative exposure to respirable elemental carbon (REC), a primary surrogate for diesel exhaust. Unconditional logistic regression and two-stage polytomous logistic regression were used to calculate odds ratios (ORs) and 95% confidence intervals (CIs), adjusting for smoking and other risk factors. RESULTS: Exposure to cumulative REC was associated with an increased risk of UBC; workers with cumulative REC >396 µg/m3-years had an OR of 1.61 (95% CI, 1.08-2.40). At this level of cumulative exposure, similar results were observed in the U.S. and Spain, OR = 1.75 (95% CI, 0.97-3.15) and OR = 1.54 (95% CI, 0.89-2.68), respectively. In lagged analysis, we also observed a consistent increased risk among workers with cumulative REC >396 µg/m3-years (range of ORs = 1.52-1.93) for all lag intervals evaluated (5-40 years). When we accounted for tumor subtypes defined by stage and grade, a significant association between diesel exhaust exposure and UBC was apparent (global test for association p = 0.0019). CONCLUSIONS: Combining data from two large epidemiologic studies, our results provide further evidence that diesel exhaust exposure increases the risk of UBC.

t-Distributed Stochastic Neighbor Embedding (t-SNE): A tool for eco-physiological transcriptomic analysis.

High-throughput RNA sequencing (RNA-Seq) has transformed the ecophysiological assessment of individual plankton species and communities. However, the technology generates complex data consisting of millions of short-read sequences that can be difficult to analyze and interpret. New bioinformatics workflows are needed to guide experimentation, environmental sampling, and to develop and test hypotheses. One complexity-reducing tool that has been used successfully in other fields is "t-distributed Stochastic Neighbor Embedding" (t-SNE). Its application to transcriptomic data from marine pelagic and benthic systems has yet to be explored. The present study demonstrates an application for evaluating RNA-Seq data using previously published, conventionally analyzed studies on the copepods Calanus finmarchicus and Neocalanus flemingeri. In one application, gene expression profiles were compared among different developmental stages. In another, they were compared among experimental conditions. In a third, they were compared among environmental samples from different locations. The profile categories identified by t-SNE were validated by reference to published results using differential gene expression and Gene Ontology (GO) analyses. The analyses demonstrate how individual samples can be evaluated for differences in global gene expression, as well as differences in expression related to specific biological processes, such as lipid metabolism and responses to stress. As RNA-Seq data from plankton species and communities become more common, t-SNE analysis should provide a powerful tool for determining trends and classifying samples into groups with similar transcriptional physiology, independent of collection site or time.

Diapause within the Context of Life-History Strategies in Calanid Copepods (Calanoida: Crustacea).

Post-embryonic diapause in copepods is an adaptation that allows species in the copepod family Calanidae to thrive in high-latitude environments by transforming a short spring phytoplankton bloom into large numbers of lipid-rich individuals capable of surviving a long period of starvation. The copepods, with their high-energy lipid reservoirs, are a critical food source for higher trophic levels, making the Calanidae a key component of high-latitude marine ecosystems. The physiological ecology of the developmental program remains poorly understood. However, new studies using high-throughput RNA sequencing approaches are giving detailed access to physiological status by generating gene expression profiles for both field-collected and laboratory-incubated individuals. These are beginning to characterize the diapause phenotype, elucidate the transcriptional and physiological progression through the diapause program, and illustrate the effects of organism-environment interactions. This paper reviews gene expression profiling studies on the life cycle and diapause program of Neocalanus flemingeri Miller (1988) that were conducted as part of a long-term observation program in the northern Gulf of Alaska. It summarizes recent findings and relates them to the ecology of this species and to that of other calanids.

Regional heterogeneity impacts gene expression in the subarctic zooplankter Neocalanus flemingeri in the northern Gulf of Alaska.

Marine pelagic species are being increasingly challenged by environmental change. Their ability to persist will depend on their capacity for physiological acclimatization. Little is known about limits of physiological plasticity in key species at the base of the food web. Here we investigate the capacity for acclimatization in the copepod Neocalanus flemingeri, which inhabits the Gulf of Alaska, a heterogeneous and highly seasonal environment. RNA-Seq analysis of field-collected pre-adults identified large regional differences in expression of genes involved in metabolic and developmental processes and response to stressors. We found that lipid synthesis genes were up-regulated in individuals from Prince William Sound and down-regulated in the Gulf of Alaska. Up-regulation of lipid catabolic genes in offshore individuals suggests they are experiencing nutritional deficits. The expression differences demonstrate physiological plasticity in response to a steep gradient in food availability. Our transcriptional analysis reveals mechanisms of acclimatization that likely contribute to the observed resilience of this population.

Going with the flow: hydrodynamic cues trigger directed escapes from a stalking predator.

In the coevolution of predator and prey, different and less well-understood rules for threat assessment apply to freely suspended organisms than to substrate-dwelling ones. Particularly vulnerable are small prey carried with the bulk movement of a surrounding fluid and thus deprived of sensory information within the bow waves of approaching predators. Some planktonic prey have solved this apparent problem, however. We quantified cues generated by the slow approach of larval clownfish ( Amphiprion ocellaris) that triggered a calanoid copepod ( Bestiolina similis) to escape before the fish could strike. To estimate water deformation around the copepod immediately preceding its jump, we represented the body of the fish as a rigid sphere in a hydrodynamic model that we parametrized with measurements of fish size, approach speed and distance to the copepod. Copepods of various developmental stages (CII-CVI) were sensitive to the water flow caused by the live predator, at deformation rates as low as 0.04 s-1. This rate is far lower than that predicted from experiments that used artificial predator-mimics. Additionally, copepods localized the source, with 87% of escapes directed away (greater than or equal to 90°) from the predator. Thus, copepods' survival in life-threatening situations relied on their detection of small nonlinear signals within an environment of locally linear deformation.

Predatory posture and performance in a precocious larval fish targeting evasive copepods.

Predatory fishes avoid detection by prey through a stealthy approach, followed by a rapid and precise fast-start strike. Although many first-feeding fish larvae strike at non-evasive prey using an S-start, the clownfish Amphiprion ocellaris feeds on highly evasive calanoid copepods from a J-shaped position, beginning 1 day post-hatch (dph). We quantified this unique strike posture by observing successful predatory interactions between larval clownfish (1 to 14 dph) and three developmental stages of the calanoid copepod Bestiolina similis The J-shaped posture of clownfish became less tightly curled (more L-shaped) during larval development. Larvae were also less tightly curled when targeting adult copepods, which are more evasive than younger copepod stages. Strike performance measured as time to capture and as peak speed improved only slightly with larval age. Therefore, the J-posture may allow first-feeding larvae to minimize disturbance during their approach of sensitive prey, and may represent an alternative predatory strategy to the prototypical S-start.

Physiological characterization of the emergence from diapause: A transcriptomics approach.

Organisms inhabiting high-latitude environments have evolved adaptations, such as diapause to time reproduction and growth to optimize their survival. However, the physiological regulation of the timing of complex life histories is poorly understood, particularly for marine copepods, that diapause at depth. A member of the pelagic community of the sub-Arctic Pacific Ocean, Neocalanus flemingeri enters diapause in June. Egg production occurs in winter/spring. In order to characterize the transition from diapause to egg release, females were collected in late September from 400-700 m depth, incubated in the dark at 4-5 °C and sampled for RNASeq at weekly intervals. The diapause phenotype showed down-regulation of protein turnover and up-regulation of stress genes. Activation of the reproductive program was marked by the up-regulation of genes involved in germline development. Thereafter, progress through phases of oocyte development could be linked to changes in gene expression. At 5 weeks, females showed up-regulation of spermatogenesis, indicating that stored sperm had been in a quiescent stage and completed their maturation inside the female. Gene expression profiles provide a framework to stage field-collected females. The 7-week progression from diapause to late oogenesis suggests that females typically spawning in January initiated the reproductive program in November.

The effects of food stoichiometry and temperature on copepods are mediated by ontogeny.

Climate change is warming the oceans, increasing carbon dioxide partial pressure and reducing nutrient recycling from deep layers. This will affect carbon (C) and phosphorus (P) availability in the oceans, thus, altering the balance between the nutrient content of consumers and their food resource. The combined effects of food quality and temperature have been investigated for adult copepods; however, nauplii, the early developmental stages of copepods, often far outnumber adults, grow more rapidly and have a higher phosphorus body content and demand than later life stages. Consequently, ontogeny may affect how copepods respond to the combined stressors of increasing temperature and altered food stoichiometry. We conducted temperature-controlled experiments (24, 28 and 32 °C) where Parvocalanus crassirostris was fed either a P-replete or a P-limited phytoplankton food source. Reduced survival of nauplii and copepodites at the highest temperature was ameliorated when fed P-replete food. At higher temperatures, copepodite growth remained stable, but internal C:P stoichiometry diverged in the direction of phytoplankton C:P, suggesting that increased temperature affected copepodite stoichiometric homeostasis. In contrast, naupliar P content increased with temperature and naupliar growth was P limited, suggesting nauplii required additional phosphorus at higher temperatures. We conclude that resource stoichiometry plays a key role in how copepod survival and growth are impacted by temperature, and that ontogeny mediates these responses. Our results suggest that as the extent of warming oceans and phytoplankton nutrient limitation increase, copepod survival and the growth of early life stages may decline.