From tides to nucleotides: Genomic signatures of adaptation to environmental heterogeneity in barnacles

The northern acorn barnacle (Semibalanus balanoides) is a robust system to study the genetic basis of adaptations to highly heterogeneous environments. Adult barnacles may be exposed to highly dissimilar levels of thermal stress depending on where they settle in the intertidal (i.e., closer to the upper or lower tidal boundary). For instance, barnacles near the upper tidal limit experience episodic summer temperatures above recorded heat coma levels. This differential stress at the microhabitat level is also dependent on the aspect of sun exposure. In the present study, we used pool-seq approaches to conduct a genome wide screen for loci responding to intertidal zonation across the North Atlantic basin (Maine, Rhode Island, and Norway). Our analysis discovered 382 genomic regions containing SNPs which are consistently zonated (i.e., SNPs whose frequencies vary depending on their position in the rocky intertidal) across all surveyed habitats. Notably, most zonated SNPs are young and private to the North Atlantic. These regions show high levels of genetic differentiation across ecologically extreme microhabitats concomitant with elevated levels of genetic variation and Tajima's D, suggesting the action of non-neutral processes. Overall, these findings support the hypothesis that spatially heterogeneous selection is a general and repeatable feature for this species, and that natural selection can maintain functional genetic variation in heterogeneous environments.publishedVersio

SARS-CoV-2 specific immune responses in overweight and obese COVID-19 patients

Obesity is a known risk factor for severe respiratory tract infections. In this prospective study, we assessed the impact of being obese or overweight on longitudinal SARS-CoV-2 humoral and cellular responses up to 18 months after infection. 274 patients provided blood samples at regular time intervals up to 18 months including obese (BMI ≥30, n=32), overweight (BMI 25-29.9, n=103) and normal body weight (BMI 18.5-24.9, n=134) SARS-CoV-2 patients. We determined SARS-CoV-2 spike-specific IgG, IgA, IgM levels by ELISA and neutralising antibody titres by neutralisation assay. RBD- and spike-specific memory B cells were investigated by ELISpot, spike- and non-spike-specific IFN-γ, IL-2 and IFN-γ/IL-2 secreting T cells by FluoroSpot and T cell receptor (TCR) sequencing was performed. Higher BMI correlated with increased COVID-19 severity. Humoral and cellular responses were stronger in overweight and obese patients than normal weight patients and associated with higher spike-specific IgG binding titres relative to neutralising antibody titres. Linear regression models demonstrated that BMI, age and COVID-19 severity correlated independently with higher SARS-CoV-2 immune responses. We found an increased proportion of unique SARS-CoV-2 specific T cell clonotypes after infection in overweight and obese patients. COVID-19 vaccination boosted humoral and cellular responses irrespective of BMI, although stronger immune boosting was observed in normal weight patients. Overall, our results highlight more severe disease and an over-reactivity of the immune system in overweight and obese patients after SARS-CoV-2 infection, underscoring the importance of recognizing overweight/obese individuals as a risk group for prioritisation for COVID-19 vaccination

evolutionPoster.pdf

Hybrid genome assembly of the northern acorn barnacle, S. balanoides

Identifying structural variants using linked-read sequencing data.

MOTIVATION:Structural variation, including large deletions, duplications, inversions, translocations and other rearrangements, is common in human and cancer genomes. A number of methods have been developed to identify structural variants from Illumina short-read sequencing data. However, reliable identification of structural variants remains challenging because many variants have breakpoints in repetitive regions of the genome and thus are difficult to identify with short reads. The recently developed linked-read sequencing technology from 10X Genomics combines a novel barcoding strategy with Illumina sequencing. This technology labels all reads that originate from a small number (∼5 to 10) DNA molecules ∼50 Kbp in length with the same molecular barcode. These barcoded reads contain long-range sequence information that is advantageous for identification of structural variants. RESULTS:We present Novel Adjacency Identification with Barcoded Reads (NAIBR), an algorithm to identify structural variants in linked-read sequencing data. NAIBR predicts novel adjacencies in an individual genome resulting from structural variants using a probabilistic model that combines multiple signals in barcoded reads. We show that NAIBR outperforms several existing methods for structural variant identification-including two recent methods that also analyze linked-reads-on simulated sequencing data and 10X whole-genome sequencing data from the NA12878 human genome and the HCC1954 breast cancer cell line. Several of the novel somatic structural variants identified in HCC1954 overlap known cancer genes. AVAILABILITY AND IMPLEMENTATION:Software is available at compbio.cs.brown.edu/software. SUPPLEMENTARY INFORMATION:Supplementary data are available at Bioinformatics online

netNMF-sc: leveraging gene-gene interactions for imputation and dimensionality reduction in single-cell expression analysis.

Single-cell RNA-sequencing (scRNA-seq) enables high-throughput measurement of RNA expression in single cells. However, because of technical limitations, scRNA-seq data often contain zero counts for many transcripts in individual cells. These zero counts, or dropout events, complicate the analysis of scRNA-seq data using standard methods developed for bulk RNA-seq data. Current scRNA-seq analysis methods typically overcome dropout by combining information across cells in a lower-dimensional space, leveraging the observation that cells generally occupy a small number of RNA expression states. We introduce netNMF-sc, an algorithm for scRNA-seq analysis that leverages information across both cells and genes. netNMF-sc learns a low-dimensional representation of scRNA-seq transcript counts using network-regularized non-negative matrix factorization. The network regularization takes advantage of prior knowledge of gene-gene interactions, encouraging pairs of genes with known interactions to be nearby each other in the low-dimensional representation. The resulting matrix factorization imputes gene abundance for both zero and nonzero counts and can be used to cluster cells into meaningful subpopulations. We show that netNMF-sc outperforms existing methods at clustering cells and estimating gene-gene covariance using both simulated and real scRNA-seq data, with increasing advantages at higher dropout rates (e.g., >60%). We also show that the results from netNMF-sc are robust to variation in the input network, with more representative networks leading to greater performance gains

Symptom burden and immune dynamics 6 to 18 months following mild SARS-CoV-2 infection -a case-control study

Background The burden and duration of persistent symptoms after nonsevere coronavirus disease 2019 (COVID-19) remains uncertain. This study aimed to assess postinfection symptom trajectories in home-isolated COVID-19 cases compared with age- and time- matched seronegative controls, and investigate immunological correlates of long COVID. Methods A prospective case-control study included home-isolated COVID-19 cases between February 28 and April 4, 2020, and followed for 12 (n = 233) to 18 (n = 149) months, and 189 age-matched severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-naive controls. We collected clinical data at baseline, 6, 12, and 18 months postinfection, and blood samples at 2, 4, 6, and 12 months for analysis of SARS-CoV-2-specific humoral and cellular responses. Results Overall, 46% (108/233) had persisting symptoms 12 months after COVID-19. Compared with controls, adult cases had a high risk of fatigue (27% excess risk, sex, and comorbidity adjusted odds ratio [aOR] 5.86; 95% confidence interval [CI], 3.27–10.5), memory problems (21% excess risk; aOR 7.42; CI, 3.51–15.67), concentration problems (20% excess risk; aOR 8.88; 95% CI, 3.88–20.35), and dyspnea (10% excess risk; aOR 2.66; 95% CI, 1.22–5.79). The prevalence of memory problems increased overall from 6 to 18 months (excess risk 11.5%; 95% CI, 1.5–21.5; P = .024) and among women (excess risk 18.7%; 95% CI, 4.4–32.9; P = .010). Longitudinal spike immunoglobulin G was significantly associated with dyspnea at 12 months. The spike-specific clonal CD4+ T-cell receptor β depth was significantly associated with both dyspnea and number of symptoms at 12 months. Conclusions This study documents a high burden of persisting symptoms after mild COVID-19 and suggests that infection induced SARS-CoV-2-specific immune responses may influence long-term symptoms.publishedVersio

From tides to nucleotides: Genomic signatures of adaptation to environmental heterogeneity in barnacles

The northern acorn barnacle (Semibalanus balanoides) is a robust system to study the genetic basis of adaptations to highly heterogeneous environments. Adult barnacles may be exposed to highly dissimilar levels of thermal stress depending on where they settle in the intertidal (i.e., closer to the upper or lower tidal boundary). For instance, barnacles near the upper tidal limit experience episodic summer temperatures above recorded heat coma levels. This differential stress at the microhabitat level is also dependent on the aspect of sun exposure. In the present study, we used pool-seq approaches to conduct a genome wide screen for loci responding to intertidal zonation across the North Atlantic basin (Maine, Rhode Island, and Norway). Our analysis discovered 382 genomic regions containing SNPs which are consistently zonated (i.e., SNPs whose frequencies vary depending on their position in the rocky intertidal) across all surveyed habitats. Notably, most zonated SNPs are young and private to the North Atlantic. These regions show high levels of genetic differentiation across ecologically extreme microhabitats concomitant with elevated levels of genetic variation and Tajima's D, suggesting the action of non-neutral processes. Overall, these findings support the hypothesis that spatially heterogeneous selection is a general and repeatable feature for this species, and that natural selection can maintain functional genetic variation in heterogeneous environments

Ecological load and balancing selection in circumboreal barnacles

Acorn barnacle adults experience environmental heterogeneity at various spatial scales of their circumboreal habitat, raising the question of how adaptation to high environmental variability is maintained in the face of strong juvenile dispersal and mortality. Here, we show that 4% of genes in the barnacle genome experience balancing selection across the entire range of the species. Many of these genes harbor mutations maintained across 2 My of evolution between the Pacific and Atlantic oceans. These genes are involved in ion regulation, pain reception, and heat tolerance, functions which are essential in highly variable ecosystems. The data also reveal complex population structure within and between basins, driven by the trans-Arctic interchange and the last glaciation. Divergence between Atlantic and Pacific populations is high, foreshadowing the onset of allopatric speciation, and suggesting that balancing selection is strong enough to maintain functional variation for millions of years in the face of complex demography

Durable immune responses after BNT162b2 vaccination in home-dwelling old adults

Objectives: Elderly are an understudied, high-risk group vulnerable to severe COVID-19. We comprehensively analyzed the durability of humoral and cellular immune responses after BNT162b2 vaccination and SARS-CoV-2 infection in elderly and younger adults. Methods: Home-dwelling old (n = 100, median 86 years) and younger adults (n = 449, median 38 years) were vaccinated with two doses of BNT162b2 vaccine at 3-week intervals and followed for 9-months. Vaccine-induced responses were compared to home-isolated COVID-19 patients (n = 183, median 47 years). Our analysis included neutralizing antibodies, spike-specific IgG, memory B-cells, IFN-γ and IL-2 secreting T-cells and sequencing of the T-cell receptor (TCR) repertoire. Results: Spike-specific breadth and depth of the CD4+ and CD8+ TCR repertoires were significantly lower in the elderly after one and two vaccinations. Both vaccinations boosted IFN-γ and IL-2 secreting spike-specific T-cells responses, with 96 % of the elderly and 100 % of the younger adults responding after the second dose, although responses were not maintained at 9-months. In contrast, T-cell responses persisted up to 12-months in infected patients. Spike-specific memory B-cells were induced after the first dose in 87 % of the younger adults compared to 38 % of the elderly, which increased to 83 % after the second dose. Memory B-cells were maintained at 9-months post-vaccination in both vaccination groups. Neutralizing antibody titers were estimated to last for 1-year in younger adults but only 6-months in the older vaccinees. Interestingly, infected older patients (n = 15, median 75 years) had more durable neutralizing titers estimated to last 14-months, 8-months longer than the older vaccinees. Conclusions: Vaccine-induced spike-specific IgG and neutralizing antibodies were consistently lower in the older than younger vaccinees. Overall, our data provide valuable insights into the kinetics of the humoral and cellular immune response in the elderly after SARS-CoV-2 vaccination or infection, highlighting the need for two doses, which can guide future vaccine design.Clinical trials.gov; NCT04706390