Y chromosome sequence and epigenomic reconstruction across human populations.

Recent advances in long-read sequencing technologies have allowed the generation and curation of more complete genome assemblies, enabling the analysis of traditionally neglected chromosomes, such as the human Y chromosome (chrY). Native DNA was sequenced on a MinION Oxford Nanopore Technologies sequencing device to generate genome assemblies for seven major chrY human haplogroups. We analyzed and compared the chrY enrichment of sequencing data obtained using two different selective sequencing approaches: adaptive sampling and flow cytometry chromosome sorting. We show that adaptive sampling can produce data to create assemblies comparable to chromosome sorting while being a less expensive and time-consuming technique. We also assessed haplogroup-specific structural variants, which would be otherwise difficult to study using short-read sequencing data only. Finally, we took advantage of this technology to detect and profile epigenetic modifications among the considered haplogroups. Altogether, we provide a framework to study complex genomic regions with a simple, fast, and affordable methodology that could be applied to larger population genomics datasets.

Genomic Adaptations to an Endoparasitic Lifestyle in the Morphologically Atypical Crustacean Sacculina carcini (Cirripedia: Rhizocephala).

The endoparasitic crustacean Sacculina carcini (Cirripedia: Rhizocephala) has a much simpler morphology than conventional filter-feeding barnacles, reflecting its parasitic lifestyle. To investigate the molecular basis of its refined developmental program, we produced a draft genome sequence for comparison with the genomes of nonparasitic barnacles and characterized the transcriptomes of internal and external tissues. The comparison of clusters of orthologous genes revealed the depletion of multiple gene families but also several unanticipated expansions compared to non-parasitic crustaceans. Transcriptomic analyses comparing interna and externa tissues revealed an unexpected variation of gene expression between rootlets sampled around host midgut and thoracic ganglia. Genes associated with lipid uptake were strongly expressed by the internal tissues. We identified candidate genes probably involved in host manipulation (suppression of ecdysis and gonad development) including those encoding crustacean neurohormones and the juvenile hormone binding protein. The evolution of Rhizocephala therefore appears to have involved a rapid turnover of genes (losses and expansions) as well as the fine tuning of gene expression.

In vitro and in vivo evidence that the switch from calcineurin to mTOR inhibitors may be a strategy for immunosuppression in Epstein-Barr virus-associated post-transplant lymphoproliferative disorder.

Post-transplant lymphoproliferative disorder is a life-threatening complication of immunosuppression following transplantation mediated by failure of T cells to control Epstein-Barr virus (EBV)-infected and transformed B cells. Typically, a modification or reduction of immunosuppression is recommended, but insufficiently defined thus far. In order to help delineate this, we characterized EBV-antigen-specific T cells and lymphoblastoid cell lines from healthy donors and in patients with a kidney transplant in the absence or presence of the standard immunosuppressants tacrolimus, cyclosporin A, prednisolone, rapamycin, and mycophenolic acid. Phenotypes of lymphoblastoid cell-lines and T cells, T cell-receptor-repertoire diversity, and T-cell reactivity upon co-culture with autologous lymphoblastoid cell lines were analyzed. Rapamycin and mycophenolic acid inhibited lymphoblastoid cell-line proliferation. T cells treated with prednisolone and rapamycin showed nearly normal cytokine production. Proliferation and the viability of T cells were decreased by mycophenolic acid, while tacrolimus and cyclosporin A were strong suppressors of T-cell function including their killing activity. Overall, our study provides a basis for the clinical decision for the modification and reduction of immunosuppression and adds information to the complex balance of maintaining anti-viral immunity while preventing acute rejection. Thus, an immunosuppressive regime based on mTOR inhibition and reduced or withdrawn calcineurin inhibitors could be a promising strategy for patients with increased risk of or manifested EBV-associated post-transplant lymphoproliferative disorder.

Alternative splicing of BUD13 determines the severity of a developmental disorder with lipodystrophy and progeroid features.

PURPOSE: In this study we aimed to identify the molecular genetic cause of a progressive multisystem disease with prominent lipodystrophy. METHODS: In total, 5 affected individuals were investigated using exome sequencing. Dermal fibroblasts were characterized using RNA sequencing, proteomics, immunoblotting, immunostaining, and electron microscopy. Subcellular localization and rescue studies were performed. RESULTS: We identified a lipodystrophy phenotype with a typical facial appearance, corneal clouding, achalasia, progressive hearing loss, and variable severity. Although 3 individuals showed stunted growth, intellectual disability, and died within the first decade of life (A1, A2, and A3), 2 are adults with normal intellectual development (A4 and A5). All individuals harbored an identical homozygous nonsense variant affecting the retention and splicing complex component BUD13. The nucleotide substitution caused alternative splicing of BUD13 leading to a stable truncated protein whose expression positively correlated with disease expression and life expectancy. In dermal fibroblasts, we found elevated intron retention, a global reduction of spliceosomal proteins, and nuclei with multiple invaginations, which were more pronounced in A1, A2, and A3. Overexpression of both BUD13 isoforms normalized the nuclear morphology. CONCLUSION: Our results define a hitherto unknown syndrome and show that the alternative splice product converts a loss-of-function into a hypomorphic allele, thereby probably determining the severity of the disease and the survival of affected individuals.

Citizen-science reveals changes in the oral microbiome in Spain through age and lifestyle factors.

The relevance of the human oral microbiome to our understanding of human health has grown in recent years as microbiome studies continue to develop. Given the links of the oral cavity with the digestive, respiratory and circulatory systems, the composition of the oral microbiome is relevant beyond just oral health, impacting systemic processes across the body. However, we still have a very limited understanding about intrinsic and extrinsic factors that shape the composition of the healthy oral microbiome. Here, we followed a citizen-science approach to assess the relative impact on the oral microbiome of selected biological, social, and lifestyle factors in 1648 Spanish individuals. We found that the oral microbiome changes across age, with middle ages showing a more homogeneous composition, and older ages showing more diverse microbiomes with increased representation of typically low abundance taxa. By measuring differences within and between groups of individuals sharing a given parameter, we were able to assess the relative impact of different factors in driving specific microbial compositions. Chronic health disorders present in the analyzed population were the most impactful factors, followed by smoking and the presence of yeasts in the oral cavity. Finally, we corroborate findings in the literature that relatives tend to have more similar oral microbiomes, and show for the first time a similar effect for classmates. Multiple intrinsic and extrinsic factors jointly shape the oral microbiome. Comparative analysis of metabarcoding data from a large sample set allows us to disentangle the individual effects.

The Positive Rhinovirus/Enterovirus Detection and SARS-CoV-2 Persistence beyond the Acute Infection Phase: An Intra-Household Surveillance Study.

We aimed to assess the duration of nasopharyngeal severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA persistence in adults self-confined at home after acute infection; and to identify the associations of SARS-CoV-2 persistence with respiratory virus co-detection and infection transmission. A cross-sectional intra-household study was conducted in metropolitan Barcelona (Spain) during the time period of April to June 2020. Every adult who was the first family member reported as SARS-CoV-2-positive by reverse transcription polymerase chain reaction (RT-PCR) as well as their household child contacts had nasopharyngeal swabs tested by a targeted SARS-CoV-2 RT-PCR and a multiplex viral respiratory panel after a 15 day minimum time lag. Four-hundred and four households (404 adults and 708 children) were enrolled. SARS-CoV-2 RNA was detected in 137 (33.9%) adults and 84 (11.9%) children. Rhinovirus/Enterovirus (RV/EV) was commonly found (83.3%) in co-infection with SARS-CoV-2 in adults. The mean duration of SARS-CoV-2 RNA presence in adults' nasopharynx was 52 days (range 26-83 days). The persistence of SARS-CoV-2 was significantly associated with RV/EV co-infection (adjusted odds ratio (aOR) 9.31; 95% CI 2.57-33.80) and SARS-CoV-2 detection in child contacts (aOR 2.08; 95% CI 1.24-3.51). Prolonged nasopharyngeal SARS-CoV-2 RNA persistence beyond the acute infection phase was frequent in adults quarantined at home during the first epidemic wave; which was associated with RV/EV co-infection and could enhance intra-household infection transmission.

Citizen-science based study of the oral microbiome in Cystic fibrosis and matched controls reveals major differences in diversity and abundance of bacterial and fungal species.

Introduction: Cystic fibrosis (CF) is an autosomal genetic disease, associated with the production of excessively thick mucosa and with life-threatening chronic lung infections. The microbiota of the oral cavity can act as a reservoir or as a barrier for infectious microorganisms that can colonize the lungs. However, the specific composition of the oral microbiome in CF is poorly understood.Methods: In collaboration with CF associations in Spain, we collected oral rinse samples from 31 CF persons (age range 7-47) and matched controls, and then performed 16S rRNA metabarcoding and high-throughput sequencing, combined with culture and proteomics-based identification of fungi to survey the bacterial and fungal oral microbiome.Results: We found that CF is associated with less diverse oral microbiomes, which were characterized by higher prevalence of Candida albicans and differential abundances of a number of bacterial taxa that have implications in both the connection to lung infections in CF, as well as potential oral health concerns, particularly periodontitis and dental caries.Conclusion: Overall, our study provides a first global snapshot of the oral microbiome in CF. Future studies are required to establish the relationships between the composition of the oral and lung microbiomes in CF.

E2F6 initiates stable epigenetic silencing of germline genes during embryonic development.

In mouse development, long-term silencing by CpG island DNA methylation is specifically targeted to germline genes; however, the molecular mechanisms of this specificity remain unclear. Here, we demonstrate that the transcription factor E2F6, a member of the polycomb repressive complex 1.6 (PRC1.6), is critical to target and initiate epigenetic silencing at germline genes in early embryogenesis. Genome-wide, E2F6 binds preferentially to CpG islands in embryonic cells. E2F6 cooperates with MGA to silence a subgroup of germline genes in mouse embryonic stem cells and in embryos, a function that critically depends on the E2F6 marked box domain. Inactivation of E2f6 leads to a failure to deposit CpG island DNA methylation at these genes during implantation. Furthermore, E2F6 is required to initiate epigenetic silencing in early embryonic cells but becomes dispensable for the maintenance in differentiated cells. Our findings elucidate the mechanisms of epigenetic targeting of germline genes and provide a paradigm for how transient repression signals by DNA-binding factors in early embryonic cells are translated into long-term epigenetic silencing during mouse development.

A global metagenomic map of urban microbiomes and antimicrobial resistance.

We present a global atlas of 4,728 metagenomic samples from mass-transit systems in 60 cities over 3 years, representing the first systematic, worldwide catalog of the urban microbial ecosystem. This atlas provides an annotated, geospatial profile of microbial strains, functional characteristics, antimicrobial resistance (AMR) markers, and genetic elements, including 10,928 viruses, 1,302 bacteria, 2 archaea, and 838,532 CRISPR arrays not found in reference databases. We identified 4,246 known species of urban microorganisms and a consistent set of 31 species found in 97% of samples that were distinct from human commensal organisms. Profiles of AMR genes varied widely in type and density across cities. Cities showed distinct microbial taxonomic signatures that were driven by climate and geographic differences. These results constitute a high-resolution global metagenomic atlas that enables discovery of organisms and genes, highlights potential public health and forensic applications, and provides a culture-independent view of AMR burden in cities.

The intratumoral CXCR3 chemokine system is predictive of chemotherapy response in human bladder cancer.

Chemotherapy has direct toxic effects on cancer cells; however, long-term cancer control and complete remission are likely to involve CD8+ T cell immune responses. To study the role of CD8+ T cell infiltration in the success of chemotherapy, we examined patients with muscle invasive bladder cancer (MIBC) who were categorized on the basis of the response to neoadjuvant chemotherapy (NAC). We identified the intratumoral CXCR3 chemokine system (ligands and receptor splice variants) as a critical component for tumor eradication upon NAC in MIBC. Through characterization of CD8+ T cells, we found that stem-like T cell subpopulations with abundant CXCR3alt, a variant form of the CXCL11 receptor, responded to CXCL11 in culture as demonstrated by migration and enhanced effector function. In tumor biopsies of patients with MIBC accessed before treatment, CXCL11 abundance correlated with high numbers of tumor-infiltrating T cells and response to NAC. The presence of CXCR3alt and CXCL11 was associated with improved overall survival in MIBC. Evaluation of both CXCR3alt and CXCL11 enabled discrimination between responder and nonresponder patients with MIBC before treatment. We validated the prognostic role of the CXCR3-CXCL11 chemokine system in an independent cohort of chemotherapy-treated and chemotherapy-naïve patients with MIBC from data in TCGA. In summary, our data revealed stimulatory activity of the CXCR3alt-CXCL11 chemokine system on CD8+ T cells that is predictive of chemotherapy responsiveness in MIBC. This may offer immunotherapeutic options for targeted activation of intratumoral stem-like T cells in solid tumors.

SARS-CoV-2 Seroprevalence and Antibody Kinetics Among Health Care Workers in a Spanish Hospital After 3 Months of Follow-up.

BACKGROUND: At the COVID-19 spring 2020 pandemic peak in Spain, prevalence of SARS-CoV-2 infection in a cohort of 578 randomly selected health care workers (HCWs) from Hospital Clínic de Barcelona was 11.2%. METHODS: A follow-up survey 1 month later (April-May 2020) measured infection by rRT-PCR and IgM, IgA, and IgG to the receptor-binding domain of the spike protein by Luminex. Antibody kinetics, including IgG subclasses, was assessed until month 3. RESULTS: At month 1, the prevalence of infection measured by rRT-PCR and serology was 14.9% (84/565) and seroprevalence 14.5% (82/565). We found 25 (5%) new infections in 501 participants without previous evidence of infection. IgM, IgG, and IgA levels declined in 3 months (antibody decay rates 0.15 [95% CI, .11-.19], 0.66 [95% CI, .54-.82], and 0.12 [95% CI, .09-.16], respectively), and 68.33% of HCWs had seroreverted for IgM, 3.08% for IgG, and 24.29% for IgA. The most frequent subclass responses were IgG1 (highest levels) and IgG2, followed by IgG3, and only IgA1 but no IgA2 was detected. CONCLUSIONS: Continuous and improved surveillance of SARS-CoV-2 infections in HCWs remains critical, particularly in high-risk groups. The observed fast decay of IgA and IgM levels has implications for seroprevalence studies using these isotypes.

Seroprevalence of antibodies against SARS-CoV-2 among health care workers in a large Spanish reference hospital.

Health care workers (HCW) are a high-risk population to acquire SARS-CoV-2 infection from patients or other fellow HCW. This study aims at estimating the seroprevalence against SARS-CoV-2 in a random sample of HCW from a large hospital in Spain. Of the 578 participants recruited from 28 March to 9 April 2020, 54 (9.3%, 95% CI: 7.1-12.0) were seropositive for IgM and/or IgG and/or IgA against SARS-CoV-2. The cumulative prevalence of SARS-CoV-2 infection (presence of antibodies or past or current positive rRT-PCR) was 11.2% (65/578, 95% CI: 8.8-14.1). Among those with evidence of past or current infection, 40.0% (26/65) had not been previously diagnosed with COVID-19. Here we report a relatively low seroprevalence of antibodies among HCW at the peak of the COVID-19 epidemic in Spain. A large proportion of HCW with past or present infection had not been previously diagnosed with COVID-19, which calls for active periodic rRT-PCR testing in hospital settings.

Oral microbiome in down syndrome and its implications on oral health.

Introduction: The oral cavity harbors an abundant and diverse microbial community (i.e. the microbiome), whose composition and roles in health and disease have been the focus of intense research. Down syndrome (DS) is associated with particular characteristics in the oral cavity, and with a lower incidence of caries and higher incidence of periodontitis and gingivitis compared to control populations. However, the overall composition of the oral microbiome in DS and how it varies with diverse factors like host age or the pH within the mouth are still poorly understood. Methods: Using a Citizen-Science approach in collaboration with DS associations in Spain, we performed 16S rRNA metabarcoding and high-throughput sequencing, combined with culture and proteomics-based identification of fungi to survey the bacterial and fungal oral microbiome in 27 DS persons (age range 7-55) and control samples matched by geographical distribution, age range, and gender. Results: We found that DS is associated with low salivary pH and less diverse oral microbiomes, which were characterized by lower levels of Alloprevotella, Atopobium, Candidatus Saccharimonas, and higher amounts of Kingella, Staphylococcus, Gemella, Cardiobacterium, Rothia, Actinobacillus, and greater prevalence of Candida. Conclusion: Altogether, our study provides a first global snapshot of the oral microbiome in DS. Future studies are required to establish whether the observed differences are related to differential pathology in the oral cavity in DS.

Increased presence and differential molecular imprinting of transit amplifying cells in psoriasis.

Psoriasis is a very common chronic inflammatory skin disease characterized by epidermal thickening and scaling resulting from keratinocyte hyperproliferation and impaired differentiation. Pathomechanistic studies in psoriasis are often limited by using whole skin tissue biopsies, neglecting their stratification and cellular diversity. This study aimed at characterizing epidermal alterations in psoriasis at the level of keratinocyte populations. Epidermal cell populations were purified from skin biopsies of psoriasis patients and healthy donors using a novel cell type-specific approach. Molecular characterization of the transit-amplifying cells (TAC), the key players of epidermal renewal, was performed using immunocytofluorescence-technique and integrated multiscale-omics analyses. Already TAC from non-lesional psoriatic skin showed altered methylation and differential expression in 1.7% and 1.0% of all protein-coding genes, respectively. In psoriatic lesions, TAC were strongly expanded showing further increased differentially methylated (10-fold) and expressed (22-fold) genes numbers. Importantly, 17.2% of differentially expressed genes were associated with respective gene methylations. Compared with non-lesional TAC, pathway analyses revealed metabolic alterations as one feature predominantly changed in TAC derived from active psoriatic lesions. Overall, our study showed stage-specific molecular alterations, allows new insights into the pathogenesis, and implies the involvement of epigenetic mechanisms in lesion development in psoriasis. KEY MESSAGES: Transit amplifying cell (TAC) numbers are highly increased in psoriatic lesions Psoriatic TAC show profound molecular alterations & stage-specific identity TAC from unaffected areas already show first signs of molecular alterations Lesional TAC show a preference in metabolic-related alterations.

Computational Processing and Quality Control of Hi-C, Capture Hi-C and Capture-C Data.

Hi-C, capture Hi-C (CHC) and Capture-C have contributed greatly to our present understanding of the three-dimensional organization of genomes in the context of transcriptional regulation by characterizing the roles of topological associated domains, enhancer promoter loops and other three-dimensional genomic interactions. The analysis is based on counts of chimeric read pairs that map to interacting regions of the genome. However, the processing and quality control presents a number of unique challenges. We review here the experimental and computational foundations and explain how the characteristics of restriction digests, sonication fragments and read pairs can be exploited to distinguish technical artefacts from valid read pairs originating from true chromatin interactions.

The Role of Pre-existing Cross-Reactive Central Memory CD4 T-Cells in Vaccination With Previously Unseen Influenza Strains.

Influenza vaccination is a common approach to prevent seasonal and pandemic influenza. Pre-existing antibodies against close viral strains might impair antibody formation against previously unseen strains-a process called original antigenic sin. The role of this pre-existing cellular immunity in this process is, despite some hints from animal models, not clear. Here, we analyzed cellular and humoral immunity in healthy individuals before and after vaccination with seasonal influenza vaccine. Based on influenza-specific hemagglutination inhibiting (HI) titers, vaccinees were grouped into HI-negative and -positive cohorts followed by in-depth cytometric and TCR repertoire analysis. Both serological groups revealed cross-reactive T-cell memory to the vaccine strains at baseline that gave rise to the majority of vaccine-specific T-cells post vaccination. On the contrary, very limited number of vaccine-specific T-cell clones was recruited from the naive pool. Furthermore, baseline quantity of vaccine-specific central memory helper T-cells and clonotype richness of this population directly correlated with the vaccination efficacy. Our findings suggest that the deliberate recruitment of pre-existing cross-reactive cellular memory might help to improve vaccination outcome.

BKV Clearance Time Correlates With Exhaustion State and T-Cell Receptor Repertoire Shape of BKV-Specific T-Cells in Renal Transplant Patients.

Reactivation of the BK polyomavirus is known to lead to severe complications in kidney transplant patients. The current treatment strategy relies on decreasing the immunosuppression to allow the immune system to clear the virus. Recently, we demonstrated a clear association between the resolution of BKV reactivation and reconstitution of BKV-specific CD4+ T-cells. However, which factors determine the duration of viral infection clearance remains so far unclear. Here we apply a combination of in-depth multi-parametric flow cytometry and NGS-based CDR3 beta chain receptor repertoire analysis of BKV-specific T-cells to a cohort of 7 kidney transplant patients during the clinical course of BKV reactivation. This way we followed TCR repertoires at single clone levels and functional activity of BKV-specific T-cells during the resolution of BKV infection. The duration of BKV clearance did not depend on the number of peripheral blood BKV-specific T-cells nor on a few immunodominant BKV-specific T-cell clones. Rather, the T-cell receptor repertoire diversity and exhaustion status of BKV-specific T-cells affected the duration of viral clearance: high clonotype diversity and lack of PD1 and TIM3 exhaustion markers on BKV-specific T-cells was associated with short clearance time. Our data thus demonstrate how the diversity and the exhaustion state of the T-cells can determine the clinical course of BKV infection.

PAX5 biallelic genomic alterations define a novel subgroup of B-cell precursor acute lymphoblastic leukemia.

Chromosomal rearrangements and specific aneuploidy patterns are initiating events and define subgroups in B-cell precursor acute lymphoblastic leukemia (BCP-ALL). Here we analyzed 250 BCP-ALL cases and identified a novel subgroup ('PAX5-plus', n = 19) by distinct DNA methylation and gene expression profiles. All patients in this subgroup harbored mutations in the B-lineage transcription factor PAX5, with p.P80R as hotspot. Mutations either affected two independent codons, consistent with compound heterozygosity, or suffered LOH predominantly through chromosome 9p aberrations. These biallelic events resulted in disruption of PAX5 transcriptional programs regulating B-cell differentiation and tumor suppressor functions. Homozygous CDKN2A/B deletions and RAS-activating hotspot mutations were highly enriched as cooperating events in the genomic profile of PAX5-plus ALL. Together, this defined a specific pattern of triple alterations, exclusive to the novel subgroup. PAX5-plus ALL was observed in pediatric and adult patients. Although restricted by the limited sample size, a tendency for more favorable clinical outcome was observed, with 10 of 12 adult PAX5-plus patients achieving long-term survival. PAX5-plus represents the first BCP-ALL subgroup defined by sequence alterations in contrast to gross chromosomal events and exemplifies how deregulated differentiation (PAX5), impaired cell cycle control (CDKN2A/B) and sustained proliferative signaling (RAS) cooperatively drive leukemogenesis.

GOPHER: Generator Of Probes for capture Hi-C Experiments at high Resolution.

BACKGROUND: Target enrichment combined with chromosome conformation capturing methodologies such as capture Hi-C (CHC) can be used to investigate spatial layouts of genomic regions with high resolution and at scalable costs. A common application of CHC is the investigation of regulatory elements that are in contact with promoters, but CHC can be used for a range of other applications. Therefore, probe design for CHC needs to be adapted to experimental needs, but no flexible tool is currently available for this purpose. RESULTS: We present a Java desktop application called GOPHER (Generator Of Probes for capture Hi-C Experiments at high Resolution) that implements three strategies for CHC probe design. GOPHER's simple approach is similar to the probe design of previous approaches that employ CHC to investigate all promoters, with one probe being placed at each margin of a single digest that overlaps the transcription start site (TSS) of each promoter. GOPHER's simple-patched approach extends this methodology with a heuristic that improves coverage of viewpoints in which the TSS is located near to one of the boundaries of the digest. GOPHER's extended approach is intended mainly for focused investigations of smaller gene sets. GOPHER can also be used to design probes for regions other than TSS such as GWAS hits or large blocks of genomic sequence. GOPHER additionally provides a number of features that allow users to visualize and edit viewpoints, and outputs a range of files useful for documentation, ordering probes, and downstream analysis. CONCLUSION: GOPHER is an easy-to-use and robust desktop application for CHC probe design. Source code and a precompiled executable can be downloaded from the GOPHER GitHub page at https://github.com/TheJacksonLaboratory/Gopher .