A single-nucleus and spatial transcriptomic atlas of the COVID-19 liver reveals topological, functional, and regenerative organ disruption in patients.

The molecular underpinnings of organ dysfunction in acute COVID-19 and its potential long-term sequelae are under intense investigation. To shed light on these in the context of liver function, we performed single-nucleus RNA-seq and spatial transcriptomic profiling of livers from 17 COVID-19 decedents. We identified hepatocytes positive for SARS-CoV-2 RNA with an expression phenotype resembling infected lung epithelial cells. Integrated analysis and comparisons with healthy controls revealed extensive changes in the cellular composition and expression states in COVID-19 liver, reflecting hepatocellular injury, ductular reaction, pathologic vascular expansion, and fibrogenesis. We also observed Kupffer cell proliferation and erythrocyte progenitors for the first time in a human liver single-cell atlas, resembling similar responses in liver injury in mice and in sepsis, respectively. Despite the absence of a clinical acute liver injury phenotype, endothelial cell composition was dramatically impacted in COVID-19, concomitantly with extensive alterations and profibrogenic activation of reactive cholangiocytes and mesenchymal cells. Our atlas provides novel insights into liver physiology and pathology in COVID-19 and forms a foundational resource for its investigation and understanding.

SARS-CoV-2 infection produces chronic pulmonary epithelial and immune cell dysfunction with fibrosis in mice.

A subset of individuals who recover from coronavirus disease 2019 (COVID-19) develop post-acute sequelae of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) (PASC), but the mechanistic basis of PASC-associated lung abnormalities suffers from a lack of longitudinal tissue samples. The mouse-adapted SARS-CoV-2 strain MA10 produces an acute respiratory distress syndrome in mice similar to humans. To investigate PASC pathogenesis, studies of MA10-infected mice were extended from acute to clinical recovery phases. At 15 to 120 days after virus clearance, pulmonary histologic findings included subpleural lesions composed of collagen, proliferative fibroblasts, and chronic inflammation, including tertiary lymphoid structures. Longitudinal spatial transcriptional profiling identified global reparative and fibrotic pathways dysregulated in diseased regions, similar to human COVID-19. Populations of alveolar intermediate cells, coupled with focal up-regulation of profibrotic markers, were identified in persistently diseased regions. Early intervention with antiviral EIDD-2801 reduced chronic disease, and early antifibrotic agent (nintedanib) intervention modified early disease severity. This murine model provides opportunities to identify pathways associated with persistent SARS-CoV-2 pulmonary disease and test countermeasures to ameliorate PASC.

System-wide transcriptome damage and tissue identity loss in COVID-19 patients.

The molecular mechanisms underlying the clinical manifestations of coronavirus disease 2019 (COVID-19), and what distinguishes them from common seasonal influenza virus and other lung injury states such as acute respiratory distress syndrome, remain poorly understood. To address these challenges, we combine transcriptional profiling of 646 clinical nasopharyngeal swabs and 39 patient autopsy tissues to define body-wide transcriptome changes in response to COVID-19. We then match these data with spatial protein and expression profiling across 357 tissue sections from 16 representative patient lung samples and identify tissue-compartment-specific damage wrought by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, evident as a function of varying viral loads during the clinical course of infection and tissue-type-specific expression states. Overall, our findings reveal a systemic disruption of canonical cellular and transcriptional pathways across all tissues, which can inform subsequent studies to combat the mortality of COVID-19 and to better understand the molecular dynamics of lethal SARS-CoV-2 and other respiratory infections.

PSMA-targeting TGFß-insensitive armored CAR T cells in metastatic castration-resistant prostate cancer: a phase 1 trial.

Chimeric antigen receptor (CAR) T cells have demonstrated promising efficacy, particularly in hematologic malignancies. One challenge regarding CAR T cells in solid tumors is the immunosuppressive tumor microenvironment (TME), characterized by high levels of multiple inhibitory factors, including transforming growth factor (TGF)-ß. We report results from an in-human phase 1 trial of castration-resistant, prostate cancer-directed CAR T cells armored with a dominant-negative TGF-ß receptor (NCT03089203). Primary endpoints were safety and feasibility, while secondary objectives included assessment of CAR T cell distribution, bioactivity and disease response. All prespecified endpoints were met. Eighteen patients enrolled, and 13 subjects received therapy across four dose levels. Five of the 13 patients developed grade ≥2 cytokine release syndrome (CRS), including one patient who experienced a marked clonal CAR T cell expansion, >98% reduction in prostate-specific antigen (PSA) and death following grade 4 CRS with concurrent sepsis. Acute increases in inflammatory cytokines correlated with manageable high-grade CRS events. Three additional patients achieved a PSA reduction of ≥30%, with CAR T cell failure accompanied by upregulation of multiple TME-localized inhibitory molecules following adoptive cell transfer. CAR T cell kinetics revealed expansion in blood and tumor trafficking. Thus, clinical application of TGF-ß-resistant CAR T cells is feasible and generally safe. Future studies should use superior multipronged approaches against the TME to improve outcomes.

A model of persistent post SARS-CoV-2 induced lung disease for target identification and testing of therapeutic strategies.

COVID-19 survivors develop post-acute sequelae of SARS-CoV-2 (PASC), but the mechanistic basis of PASC-associated lung abnormalities suffers from a lack of longitudinal samples. Mouse-adapted SARS-CoV-2 MA10 produces an acute respiratory distress syndrome (ARDS) in mice similar to humans. To investigate PASC pathogenesis, studies of MA10-infected mice were extended from acute disease through clinical recovery. At 15-120 days post-virus clearance, histologic evaluation identified subpleural lesions containing collagen, proliferative fibroblasts, and chronic inflammation with tertiary lymphoid structures. Longitudinal spatial transcriptional profiling identified global reparative and fibrotic pathways dysregulated in diseased regions, similar to human COVID-19. Populations of alveolar intermediate cells, coupled with focal upregulation of pro-fibrotic markers, were identified in persistently diseased regions. Early intervention with antiviral EIDD-2801 reduced chronic disease, and early anti-fibrotic agent (nintedanib) intervention modified early disease severity. This murine model provides opportunities to identify pathways associated with persistent SARS-CoV-2 pulmonary disease and test countermeasures to ameliorate PASC.

The tumor immune microenvironment of primary and metastatic HER2- positive breast cancers utilizing gene expression and spatial proteomic profiling.

BACKGROUND: The characterization of the immune component of the tumor microenvironment (TME) of human epidermal growth factor receptor 2 positive (HER2+) breast cancer has been limited. Molecular and spatial characterization of HER2+ TME of primary, recurrent, and metastatic breast tumors has the potential to identify immune mediated mechanisms and biomarker targets that could be used to guide selection of therapies. METHODS: We examined 15 specimens from eight patients with HER2+ breast cancer: 10 primary breast tumors (PBT), two soft tissue, one lung, and two brain metastases (BM). Using molecular profiling by bulk gene expression TME signatures, including the Tumor Inflammation Signature (TIS) and PAM50 subtyping, as well as spatial characterization of immune hot, warm, and cold regions in the stroma and tumor epithelium using 64 protein targets on the GeoMx Digital Spatial Profiler. RESULTS: PBT had higher infiltration of immune cells relative to metastatic sites and higher protein and gene expression of immune activation markers when compared to metastatic sites. TIS scores were lower in metastases, particularly in BM. BM also had less immune infiltration overall, but in the stromal compartment with the highest density of immune infiltration had similar levels of T cells that were less activated than PBT stromal regions suggesting immune exclusion in the tumor epithelium. CONCLUSIONS: Our findings show stromal and tumor localized immune cells in the TME are more active in primary versus metastatic disease. This suggests patients with early HER2+ breast cancer could have more benefit from immune-targeting therapies than patients with advanced disease.

Visualizing in deceased COVID-19 patients how SARS-CoV-2 attacks the respiratory and olfactory mucosae but spares the olfactory bulb.

Anosmia, the loss of smell, is a common and often the sole symptom of COVID-19. The onset of the sequence of pathobiological events leading to olfactory dysfunction remains obscure. Here, we have developed a postmortem bedside surgical procedure to harvest endoscopically samples of respiratory and olfactory mucosae and whole olfactory bulbs. Our cohort of 85 cases included COVID-19 patients who died a few days after infection with SARS-CoV-2, enabling us to catch the virus while it was still replicating. We found that sustentacular cells are the major target cell type in the olfactory mucosa. We failed to find evidence for infection of olfactory sensory neurons, and the parenchyma of the olfactory bulb is spared as well. Thus, SARS-CoV-2 does not appear to be a neurotropic virus. We postulate that transient insufficient support from sustentacular cells triggers transient olfactory dysfunction in COVID-19. Olfactory sensory neurons would become affected without getting infected.

The CD155/TIGIT axis promotes and maintains immune evasion in neoantigen-expressing pancreatic cancer.

The CD155/TIGIT axis can be co-opted during immune evasion in chronic viral infections and cancer. Pancreatic adenocarcinoma (PDAC) is a highly lethal malignancy, and immune-based strategies to combat this disease have been largely unsuccessful to date. We corroborate prior reports that a substantial portion of PDAC harbors predicted high-affinity MHC class I-restricted neoepitopes and extend these findings to advanced/metastatic disease. Using multiple preclinical models of neoantigen-expressing PDAC, we demonstrate that intratumoral neoantigen-specific CD8+ T cells adopt multiple states of dysfunction, resembling those in tumor-infiltrating lymphocytes of PDAC patients. Mechanistically, genetic and/or pharmacologic modulation of the CD155/TIGIT axis was sufficient to promote immune evasion in autochthonous neoantigen-expressing PDAC. Finally, we demonstrate that the CD155/TIGIT axis is critical in maintaining immune evasion in PDAC and uncover a combination immunotherapy (TIGIT/PD-1 co-blockade plus CD40 agonism) that elicits profound anti-tumor responses in preclinical models, now poised for clinical evaluation.

Systemic Tissue and Cellular Disruption from SARS-CoV-2 Infection revealed in COVID-19 Autopsies and Spatial Omics Tissue Maps.

The Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) virus has infected over 115 million people and caused over 2.5 million deaths worldwide. Yet, the molecular mechanisms underlying the clinical manifestations of COVID-19, as well as what distinguishes them from common seasonal influenza virus and other lung injury states such as Acute Respiratory Distress Syndrome (ARDS), remains poorly understood. To address these challenges, we combined transcriptional profiling of 646 clinical nasopharyngeal swabs and 39 patient autopsy tissues, matched with spatial protein and expression profiling (GeoMx) across 357 tissue sections. These results define both body-wide and tissue-specific (heart, liver, lung, kidney, and lymph nodes) damage wrought by the SARS-CoV-2 infection, evident as a function of varying viral load (high vs. low) during the course of infection and specific, transcriptional dysregulation in splicing isoforms, T cell receptor expression, and cellular expression states. In particular, cardiac and lung tissues revealed the largest degree of splicing isoform switching and cell expression state loss. Overall, these findings reveal a systemic disruption of cellular and transcriptional pathways from COVID-19 across all tissues, which can inform subsequent studies to combat the mortality of COVID-19, as well to better understand the molecular dynamics of lethal SARS-CoV-2 infection and other viruses.

The spatial landscape of lung pathology during COVID-19 progression.

Recent studies have provided insights into the pathology of and immune response to COVID-191-8. However, a thorough investigation of the interplay between infected cells and the immune system at sites of infection has been lacking. Here we use high-parameter imaging mass cytometry9 that targets the expression of 36 proteins to investigate the cellular composition and spatial architecture of acute lung injury in humans (including injuries derived from SARS-CoV-2 infection) at single-cell resolution. These spatially resolved single-cell data unravel the disordered structure of the infected and injured lung, alongside the distribution of extensive immune infiltration. Neutrophil and macrophage infiltration are hallmarks of bacterial pneumonia and COVID-19, respectively. We provide evidence that SARS-CoV-2 infects predominantly alveolar epithelial cells and induces a localized hyperinflammatory cell state that is associated with lung damage. We leverage the temporal range of fatal outcomes of COVID-19 in relation to the onset of symptoms, which reveals increased macrophage extravasation and increased numbers of mesenchymal cells and fibroblasts concomitant with increased proximity between these cell types as the disease progresses-possibly as a result of attempts to repair the damaged lung tissue. Our data enable us to develop a biologically interpretable landscape of lung pathology from a structural, immunological and clinical standpoint. We use this landscape to characterize the pathophysiology of the human lung from its macroscopic presentation to the single-cell level, which provides an important basis for understanding COVID-19 and lung pathology in general.

Shotgun transcriptome, spatial omics, and isothermal profiling of SARS-CoV-2 infection reveals unique host responses, viral diversification, and drug interactions.

In less than nine months, the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) killed over a million people, including >25,000 in New York City (NYC) alone. The COVID-19 pandemic caused by SARS-CoV-2 highlights clinical needs to detect infection, track strain evolution, and identify biomarkers of disease course. To address these challenges, we designed a fast (30-minute) colorimetric test (LAMP) for SARS-CoV-2 infection from naso/oropharyngeal swabs and a large-scale shotgun metatranscriptomics platform (total-RNA-seq) for host, viral, and microbial profiling. We applied these methods to clinical specimens gathered from 669 patients in New York City during the first two months of the outbreak, yielding a broad molecular portrait of the emerging COVID-19 disease. We find significant enrichment of a NYC-distinctive clade of the virus (20C), as well as host responses in interferon, ACE, hematological, and olfaction pathways. In addition, we use 50,821 patient records to find that renin-angiotensin-aldosterone system inhibitors have a protective effect for severe COVID-19 outcomes, unlike similar drugs. Finally, spatial transcriptomic data from COVID-19 patient autopsy tissues reveal distinct ACE2 expression loci, with macrophage and neutrophil infiltration in the lungs. These findings can inform public health and may help develop and drive SARS-CoV-2 diagnostic, prevention, and treatment strategies.

A single-cell and spatial atlas of autopsy tissues reveals pathology and cellular targets of SARS-CoV-2.

The SARS-CoV-2 pandemic has caused over 1 million deaths globally, mostly due to acute lung injury and acute respiratory distress syndrome, or direct complications resulting in multiple-organ failures. Little is known about the host tissue immune and cellular responses associated with COVID-19 infection, symptoms, and lethality. To address this, we collected tissues from 11 organs during the clinical autopsy of 17 individuals who succumbed to COVID-19, resulting in a tissue bank of approximately 420 specimens. We generated comprehensive cellular maps capturing COVID-19 biology related to patients' demise through single-cell and single-nucleus RNA-Seq of lung, kidney, liver and heart tissues, and further contextualized our findings through spatial RNA profiling of distinct lung regions. We developed a computational framework that incorporates removal of ambient RNA and automated cell type annotation to facilitate comparison with other healthy and diseased tissue atlases. In the lung, we uncovered significantly altered transcriptional programs within the epithelial, immune, and stromal compartments and cell intrinsic changes in multiple cell types relative to lung tissue from healthy controls. We observed evidence of: alveolar type 2 (AT2) differentiation replacing depleted alveolar type 1 (AT1) lung epithelial cells, as previously seen in fibrosis; a concomitant increase in myofibroblasts reflective of defective tissue repair; and, putative TP63+ intrapulmonary basal-like progenitor (IPBLP) cells, similar to cells identified in H1N1 influenza, that may serve as an emergency cellular reserve for severely damaged alveoli. Together, these findings suggest the activation and failure of multiple avenues for regeneration of the epithelium in these terminal lungs. SARS-CoV-2 RNA reads were enriched in lung mononuclear phagocytic cells and endothelial cells, and these cells expressed distinct host response transcriptional programs. We corroborated the compositional and transcriptional changes in lung tissue through spatial analysis of RNA profiles in situ and distinguished unique tissue host responses between regions with and without viral RNA, and in COVID-19 donor tissues relative to healthy lung. Finally, we analyzed genetic regions implicated in COVID-19 GWAS with transcriptomic data to implicate specific cell types and genes associated with disease severity. Overall, our COVID-19 cell atlas is a foundational dataset to better understand the biological impact of SARS-CoV-2 infection across the human body and empowers the identification of new therapeutic interventions and prevention strategies.

Minimal residual disease undetectable by next-generation sequencing predicts improved outcome in CLL after chemoimmunotherapy.

Patients with chronic lymphocytic leukemia (CLL) who achieve blood or bone marrow (BM) undetectable minimal residual disease (U-MRD) status after first-line fludarabine, cyclophosphamide, and rituximab (FCR) have prolonged progression-free survival (PFS), when assessed by an assay with sensitivity 10-4 (MRD4). Despite reaching U-MRD4, many patients, especially those with unmutated IGHV, subsequently relapse, suggesting residual disease <10-4 threshold and the need for more sensitive MRD evaluation. MRD evaluation by next-generation sequencing (NGS) has a sensitivity of 10-6 (MRD6). To better assess the depth of remission following first-line FCR treatment, we used NGS (Adaptive Biotechnologies Corporation) to assess MRD in 62 patients, all of whom had BM U-MRD by multicolor flow cytometry (sensitivity 10-4) at end-of-FCR treatment. Samples from these patients included 57 BM samples, 29 peripheral blood mononuclear cell (PBMC) samples, and 32 plasma samples. Only 27.4% of the 62 patients had U-MRD by NGS. Rate of U-MRD by NGS was lowest in BM (25%), compared with PBMC (55%) or plasma (75%). No patient with U-MRD by NGS in BM or PBMC was MRD+ in plasma. Patients with mutated IGHV were more likely to have U-MRD by NGS at the end of treatment (EOT; 41% vs 13%, P = .02) than those with unmutated IGHV. Median follow-up was 81.6 months. Patients with U-MRD at EOT had superior PFS vs MRD+ patients, regardless of sample type assessed (BM, P = .02, median not reached [NR] vs 67 months; PBMC, P = .02, median NR vs 74 months). More sensitive MRD6 testing increases prognostic discrimination over MRD4 testing.

Identifying and Tracking Low-Frequency Virus-Specific TCR Clonotypes Using High-Throughput Sequencing.

Tracking antigen-specific T cell responses over time within individuals is difficult because of lack of knowledge of antigen-specific TCR sequences, limitations in sample size, and assay sensitivities. We hypothesized that analyses of high-throughput sequencing of TCR clonotypes could provide functional readouts of individuals' immunological histories. Using high-throughput TCR sequencing, we develop a database of TCRß sequences from large cohorts of mice before (naive) and after smallpox vaccination. We computationally identify 315 vaccine-associated TCR sequences (VATS) that are used to train a diagnostic classifier that distinguishes naive from vaccinated samples in mice up to 9 months post-vaccination with >99% accuracy. We determine that the VATS library contains virus-responsive TCRs by in vitro expansion assays and virus-specific tetramer sorting. These data outline a platform for advancing our capabilities to identify pathogen-specific TCR sequences, which can be used to identify and quantitate low-frequency pathogen-specific TCR sequences in circulation over time with exceptional sensitivity.

Correction: Cryptic chytridiomycosis linked to climate and genetic variation in amphibian populations of the southeastern United States.

[This corrects the article DOI: 10.1371/journal.pone.0175843.].

Cryptic chytridiomycosis linked to climate and genetic variation in amphibian populations of the southeastern United States.

North American amphibians have recently been impacted by two major emerging pathogens, the fungus Batrachochytrium dendrobatidis (Bd) and iridoviruses in the genus Ranavirus (Rv). Environmental factors and host genetics may play important roles in disease dynamics, but few studies incorporate both of these components into their analyses. Here, we investigated the role of environmental and genetic factors in driving Bd and Rv infection prevalence and severity in a biodiversity hot spot, the southeastern United States. We used quantitative PCR to characterize Bd and Rv dynamics in natural populations of three amphibian species: Notophthalmus perstriatus, Hyla squirella and Pseudacris ornata. We combined pathogen data, genetic diversity metrics generated from neutral markers, and environmental variables into general linear models to evaluate how these factors impact infectious disease dynamics. Occurrence, prevalence and intensity of Bd and Rv varied across species and populations, but only one species, Pseudacris ornata, harbored high Bd intensities in the majority of sampled populations. Genetic diversity and climate variables both predicted Bd prevalence, whereas climatic variables alone predicted infection intensity. We conclude that Bd is more abundant in the southeastern United States than previously thought and that genetic and environmental factors are both important for predicting amphibian pathogen dynamics. Incorporating both genetic and environmental information into conservation plans for amphibians is necessary for the development of more effective management strategies to mitigate the impact of emerging infectious diseases.

Genetic regulatory network motifs constrain adaptation through curvature in the landscape of mutational (co)variance.

Systems biology is accumulating a wealth of understanding about the structure of genetic regulatory networks, leading to a more complete picture of the complex genotype-phenotype relationship. However, models of multivariate phenotypic evolution based on quantitative genetics have largely not incorporated a network-based view of genetic variation. Here we model a set of two-node, two-phenotype genetic network motifs, covering a full range of regulatory interactions. We find that network interactions result in different patterns of mutational (co)variance at the phenotypic level (the M-matrix), not only across network motifs but also across phenotypic space within single motifs. This effect is due almost entirely to mutational input of additive genetic (co)variance. Variation in M has the effect of stretching and bending phenotypic space with respect to evolvability, analogous to the curvature of space-time under general relativity, and similar mathematical tools may apply in each case. We explored the consequences of curvature in mutational variation by simulating adaptation under divergent selection with gene flow. Both standing genetic variation (the G-matrix) and rate of adaptation are constrained by M, so that G and adaptive trajectories are curved across phenotypic space. Under weak selection the phenotypic mean at migration-selection balance also depends on M.

Goldilocks Meets Santa Rosalia: An Ephemeral Speciation Model Explains Patterns of Diversification Across Time Scales.

Understanding the rate at which new species form is a key question in studying the evolution of life on earth. Here we review our current understanding of speciation rates, focusing on studies based on the fossil record, phylogenies, and mathematical models. We find that speciation rates estimated from these different studies can be dramatically different: some studies find that new species form quickly and often, while others find that new species form much less frequently. We suggest that instead of being contradictory, differences in speciation rates across different scales can be reconciled by a common model. Under the "ephemeral speciation model", speciation is very common and very rapid but the new species produced almost never persist. Evolutionary studies should therefore focus on not only the formation but also the persistence of new species.

Fat frogs, mobile genes: unexpected phylogeographic patterns for the ornate chorus frog (Pseudacris ornata).

The southeastern coastal plain of the United States is a region marked by extraordinary phylogeographic congruence that is frequently attributed to the changing sea levels that occurred during the glacial-interglacial cycles of the Pleistocene epoch. A phylogeographic break corresponding to the Apalachicola River has been suggested for many species studied to date that are endemic to this region. Here, we used this pattern of phylogeographic congruence to develop and test explicit hypotheses about the genetic structure in the ornate chorus frog (Pseudacris ornata). Using 1299 bp of mtDNA sequence and seven nuclear microsatellite markers in 13 natural populations of P. ornata, we found three clades corresponding to geographically distinct regions; one spans the Apalachicola River (Southern Clade), one encompasses Georgia and South Carolina (Central Clade) and a third comprises more northerly individuals (Northern Clade). However, it does not appear that typical phylogeographic barriers demarcate these clades. Instead, isolation by distance across the range of the entire species explained the pattern of genetic variation that we observed. We propose that P. ornata was historically widespread in the southeastern United States, and that a balance between genetic drift and migration was the root of the genetic divergence among populations. Additionally, we investigated fine-scale patterns of genetic structure and found the spatial scale at which there was significant genetic structure varied among the regions studied. Furthermore, we discuss our results in light of other phylogeographic studies of southeastern coastal plain organisms and in relation to amphibian conservation and management.

Permanent Genetic Resources added to Molecular Ecology Resources Database 1 August 2009-30 September 2009.

This article documents the addition of 238 microsatellite marker loci and 72 pairs of Single Nucleotide Polymorphism (SNP) sequencing primers to the Molecular Ecology Resources Database. Loci were developed for the following species: Adelges tsugae, Artemisia tridentata, Astroides calycularis, Azorella selago, Botryllus schlosseri, Botrylloides violaceus, Cardiocrinum cordatum var. glehnii, Campylopterus curvipennis, Colocasia esculenta, Cynomys ludovicianus, Cynomys leucurus, Cynomys gunnisoni, Epinephelus coioides, Eunicella singularis, Gammarus pulex, Homoeosoma nebulella, Hyla squirella, Lateolabrax japonicus, Mastomys erythroleucus, Pararge aegeria, Pardosa sierra, Phoenicopterus ruber ruber and Silene latifolia. These loci were cross-tested on the following species: Adelges abietis, Adelges cooleyi, Adelges piceae, Pineus pini, Pineus strobi, Tubastrea micrantha, three other Tubastrea species, Botrylloides fuscus, Botrylloides simodensis, Campylopterus hemileucurus, Campylopterus rufus, Campylopterus largipennis, Campylopterus villaviscensio, Phaethornis longuemareus, Florisuga mellivora, Lampornis amethystinus, Amazilia cyanocephala, Archilochus colubris, Epinephelus lanceolatus, Epinephelus fuscoguttatus, Symbiodinium temperate-A clade, Gammarus fossarum, Gammarus roeselii, Dikerogammarus villosus and Limnomysis benedeni. This article also documents the addition of 72 sequencing primer pairs and 52 allele specific primers for Neophocaena phocaenoides.