Resolving haplotype variation and complex genetic architecture in the human immunoglobulin kappa chain locus in individuals of diverse ancestry.

Immunoglobulins (IGs), critical components of the human immune system, are composed of heavy and light protein chains encoded at three genomic loci. The IG Kappa (IGK) chain locus consists of two large, inverted segmental duplications. The complexity of the IG loci has hindered use of standard high-throughput methods for characterizing genetic variation within these regions. To overcome these limitations, we use long-read sequencing to create haplotype-resolved IGK assemblies in an ancestrally diverse cohort (n = 36), representing the first comprehensive description of IGK haplotype variation. We identify extensive locus polymorphism, including novel single nucleotide variants (SNVs) and novel structural variants harboring functional IGKV genes. Among 47 functional IGKV genes, we identify 145 alleles, 67 of which were not previously curated. We report inter-population differences in allele frequencies for 10 IGKV genes, including alleles unique to specific populations within this dataset. We identify haplotypes carrying signatures of gene conversion that associate with SNV enrichment in the IGK distal region, and a haplotype with an inversion spanning the proximal and distal regions. These data provide a critical resource of curated genomic reference information from diverse ancestries, laying a foundation for advancing our understanding of population-level genetic variation in the IGK locus.

Vaccination induces broadly neutralizing antibody precursors to HIV gp41.

A key barrier to the development of vaccines that induce broadly neutralizing antibodies (bnAbs) against human immunodeficiency virus (HIV) and other viruses of high antigenic diversity is the design of priming immunogens that induce rare bnAb-precursor B cells. The high neutralization breadth of the HIV bnAb 10E8 makes elicitation of 10E8-class bnAbs desirable; however, the recessed epitope within gp41 makes envelope trimers poor priming immunogens and requires that 10E8-class bnAbs possess a long heavy chain complementarity determining region 3 (HCDR3) with a specific binding motif. We developed germline-targeting epitope scaffolds with affinity for 10E8-class precursors and engineered nanoparticles for multivalent display. Scaffolds exhibited epitope structural mimicry and bound bnAb-precursor human naive B cells in ex vivo screens, protein nanoparticles induced bnAb-precursor responses in stringent mouse models and rhesus macaques, and mRNA-encoded nanoparticles triggered similar responses in mice. Thus, germline-targeting epitope scaffold nanoparticles can elicit rare bnAb-precursor B cells with predefined binding specificities and HCDR3 features.

Vaccine priming of rare HIV broadly neutralizing antibody precursors in nonhuman primates.

Germline-targeting immunogens hold promise for initiating the induction of broadly neutralizing antibodies (bnAbs) to HIV and other pathogens. However, antibody-antigen recognition is typically dominated by heavy chain complementarity determining region 3 (HCDR3) interactions, and vaccine priming of HCDR3-dominant bnAbs by germline-targeting immunogens has not been demonstrated in humans or outbred animals. In this work, immunization with N332-GT5, an HIV envelope trimer designed to target precursors of the HCDR3-dominant bnAb BG18, primed bnAb-precursor B cells in eight of eight rhesus macaques to substantial frequencies and with diverse lineages in germinal center and memory B cells. We confirmed bnAb-mimicking, HCDR3-dominant, trimer-binding interactions with cryo-electron microscopy. Our results demonstrate proof of principle for HCDR3-dominant bnAb-precursor priming in outbred animals and suggest that N332-GT5 holds promise for the induction of similar responses in humans.

Multiomics Analysis of PCB126's Effect on a Mouse Chronic-Binge Alcohol Feeding Model.

BACKGROUND: Environmental pollutants, including polychlorinated biphenyls (PCBs) have been implicated in the pathogenesis of liver disease. Our group recently demonstrated that PCB126 promoted steatosis, hepatomegaly, and modulated intermediary metabolism in a rodent model of alcohol-associated liver disease (ALD). OBJECTIVE: To better understand how PCB126 promoted ALD in our previous model, the current study adopts multiple omics approaches to elucidate potential mechanistic hypotheses. METHODS: Briefly, male C57BL/6J mice were exposed to 0.2mg/kg polychlorinated biphenyl (PCB) 126 or corn oil vehicle prior to ethanol (EtOH) or control diet feeding in the chronic-binge alcohol feeding model. Liver tissues were collected and prepared for mRNA sequencing, phosphoproteomics, and inductively coupled plasma mass spectrometry for metals quantification. RESULTS: Principal component analysis showed that PCB126 uniquely modified the transcriptome in EtOH-fed mice. EtOH feeding alone resulted in >4,000 differentially expressed genes (DEGs), and PCB126 exposure resulted in more DEGs in the EtOH-fed group (907 DEGs) in comparison with the pair-fed group (503 DEGs). Top 20 significant gene ontology (GO) biological processes included "peptidyl tyrosine modifications," whereas top 25 significantly decreasing GO molecular functions included "metal/ion/zinc binding." Quantitative, label-free phosphoproteomics and western blot analysis revealed no major significant PCB126 effects on total phosphorylated tyrosine residues in EtOH-fed mice. Quantified hepatic essential metal levels were primarily significantly lower in EtOH-fed mice. PCB126-exposed mice had significantly lower magnesium, cobalt, and zinc levels in EtOH-fed mice. DISCUSSION: Previous work has demonstrated that PCB126 is a modifying factor in metabolic dysfunction-associated steatotic liver disease (MASLD), and our current work suggests that pollutants also modify ALD. PCB126 may, in part, be contributing to the malnutrition aspect of ALD, where metal deficiency is known to contribute and worsen prognosis. https://doi.org/10.1289/EHP14132.

Obesogenic polystyrene microplastic exposures disrupt the gut-liver-adipose axis.

Microplastics (MP) derived from the weathering of polymers, or synthesized in this size range, have become widespread environmental contaminants and have found their way into water supplies and the food chain. Despite this awareness, little is known about the health consequences of MP ingestion. We have previously shown that the consumption of polystyrene (PS) beads was associated with intestinal dysbiosis and diabetes and obesity in mice. To further evaluate the systemic metabolic effects of PS on the gut-liver-adipose tissue axis, we supplied C57BL/6J mice with normal water or that containing 2 sizes of PS beads (0.5 and 5 µm) at a concentration of 1 µg/ml. After 13 weeks, we evaluated indices of metabolism and liver function. As observed previously, mice drinking the PS-containing water had a potentiated weight gain and adipose expansion. Here we found that this was associated with an increased abundance of adipose F4/80+ macrophages. These exposures did not cause nonalcoholic fatty liver disease but were associated with decreased liver:body weight ratios and an enrichment in hepatic farnesoid X receptor and liver X receptor signaling. PS also increased hepatic cholesterol and altered both hepatic and cecal bile acids. Mice consuming PS beads and treated with the berry anthocyanin, delphinidin, demonstrated an attenuated weight gain compared with those mice receiving a control intervention and also exhibited a downregulation of cyclic adenosine monophosphate (cAMP) and peroxisome proliferator-activated receptor (PPAR) signaling pathways. This study highlights the obesogenic role of PS in perturbing the gut-liver-adipose axis and altering nuclear receptor signaling and intermediary metabolism. Dietary interventions may limit the adverse metabolic effects of PS consumption.

Nuclear RNA catabolism controls endogenous retroviruses, gene expression asymmetry, and dedifferentiation.

Endogenous retroviruses (ERVs) are remnants of ancient parasitic infections and comprise sizable portions of most genomes. Although epigenetic mechanisms silence most ERVs by generating a repressive environment that prevents their expression (heterochromatin), little is known about mechanisms silencing ERVs residing in open regions of the genome (euchromatin). This is particularly important during embryonic development, where induction and repression of distinct classes of ERVs occur in short temporal windows. Here, we demonstrate that transcription-associated RNA degradation by the nuclear RNA exosome and Integrator is a regulatory mechanism that controls the productive transcription of most genes and many ERVs involved in preimplantation development. Disrupting nuclear RNA catabolism promotes dedifferentiation to a totipotent-like state characterized by defects in RNAPII elongation and decreased expression of long genes (gene-length asymmetry). Our results indicate that RNA catabolism is a core regulatory module of gene networks that safeguards RNAPII activity, ERV expression, cell identity, and developmental potency.

Isoform-resolved transcriptome of the human preimplantation embryo.

Human preimplantation development involves extensive remodeling of RNA expression and splicing. However, its transcriptome has been compiled using short-read sequencing data, which fails to capture most full-length mRNAs. Here, we generate an isoform-resolved transcriptome of early human development by performing long- and short-read RNA sequencing on 73 embryos spanning the zygote to blastocyst stages. We identify 110,212 unannotated isoforms transcribed from known genes, including highly conserved protein-coding loci and key developmental regulators. We further identify 17,964 isoforms from 5,239 unannotated genes, which are largely non-coding, primate-specific, and highly associated with transposable elements. These isoforms are widely supported by the integration of published multi-omics datasets, including single-cell 8CLC and blastoid studies. Alternative splicing and gene co-expression network analyses further reveal that embryonic genome activation is associated with splicing disruption and transient upregulation of gene modules. Together, these findings show that the human embryo transcriptome is far more complex than currently known, and will act as a valuable resource to empower future studies exploring development.

Genetic variation in the immunoglobulin heavy chain locus shapes the human antibody repertoire.

Variation in the antibody response has been linked to differential outcomes in disease, and suboptimal vaccine and therapeutic responsiveness, the determinants of which have not been fully elucidated. Countering models that presume antibodies are generated largely by stochastic processes, we demonstrate that polymorphisms within the immunoglobulin heavy chain locus (IGH) impact the naive and antigen-experienced antibody repertoire, indicating that genetics predisposes individuals to mount qualitatively and quantitatively different antibody responses. We pair recently developed long-read genomic sequencing methods with antibody repertoire profiling to comprehensively resolve IGH genetic variation, including novel structural variants, single nucleotide variants, and genes and alleles. We show that IGH germline variants determine the presence and frequency of antibody genes in the expressed repertoire, including those enriched in functional elements linked to V(D)J recombination, and overlapping disease-associated variants. These results illuminate the power of leveraging IGH genetics to better understand the regulation, function, and dynamics of the antibody response in disease.

HIV-Positive Liver Transplant Does not Alter the Latent Viral Reservoir in Recipients With Antiretroviral Therapy-Suppressed HIV.

The latent viral reservoir (LVR) remains a major barrier to HIV-1 curative strategies. It is unknown whether receiving a liver transplant from a donor with HIV might lead to an increase in the LVR because the liver is a large lymphoid organ. We found no differences in intact provirus, defective provirus, or the ratio of intact to defective provirus between recipients with ART-suppressed HIV who received a liver from a donor with (n = 19) or without HIV (n = 10). All measures remained stable from baseline by 1 year posttransplant. These data demonstrate that the LVR is stable after liver transplantation in people with HIV. Clinical Trials Registration. NCT02602262 and NCT03734393.

Persistent SARS-CoV-2 Infection in a Multiple Sclerosis Patient on Ocrelizumab: A Case Report.

A 44-year-old female patient with multiple sclerosis (MS) treated with ocrelizumab was hospitalized with SARS-CoV-2 pneumonia three times over the course of five months, eventually expiring. Viral sequencing of samples from her first and last admissions suggests a single persistent SARS-CoV-2 infection. We hypothesize that her immunocompromised state, due to MS treatment with an immunosuppressive monoclonal antibody, prevented her from achieving viral clearance.

FLAIRR-Seq: A Method for Single-Molecule Resolution of Near Full-Length Antibody H Chain Repertoires.

Current Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) using short-read sequencing strategies resolve expressed Ab transcripts with limited resolution of the C region. In this article, we present the near-full-length AIRR-seq (FLAIRR-seq) method that uses targeted amplification by 5' RACE, combined with single-molecule, real-time sequencing to generate highly accurate (99.99%) human Ab H chain transcripts. FLAIRR-seq was benchmarked by comparing H chain V (IGHV), D (IGHD), and J (IGHJ) gene usage, complementarity-determining region 3 length, and somatic hypermutation to matched datasets generated with standard 5' RACE AIRR-seq using short-read sequencing and full-length isoform sequencing. Together, these data demonstrate robust FLAIRR-seq performance using RNA samples derived from PBMCs, purified B cells, and whole blood, which recapitulated results generated by commonly used methods, while additionally resolving H chain gene features not documented in IMGT at the time of submission. FLAIRR-seq data provide, for the first time, to our knowledge, simultaneous single-molecule characterization of IGHV, IGHD, IGHJ, and IGHC region genes and alleles, allele-resolved subisotype definition, and high-resolution identification of class switch recombination within a clonal lineage. In conjunction with genomic sequencing and genotyping of IGHC genes, FLAIRR-seq of the IgM and IgG repertoires from 10 individuals resulted in the identification of 32 unique IGHC alleles, 28 (87%) of which were previously uncharacterized. Together, these data demonstrate the capabilities of FLAIRR-seq to characterize IGHV, IGHD, IGHJ, and IGHC gene diversity for the most comprehensive view of bulk-expressed Ab repertoires to date.

Characterization of the immunoglobulin lambda chain locus from diverse populations reveals extensive genetic variation.

Immunoglobulins (IGs), crucial components of the adaptive immune system, are encoded by three genomic loci. However, the complexity of the IG loci severely limits the effective use of short read sequencing, limiting our knowledge of population diversity in these loci. We leveraged existing long read whole-genome sequencing (WGS) data, fosmid technology, and IG targeted single-molecule, real-time (SMRT) long-read sequencing (IG-Cap) to create haplotype-resolved assemblies of the IG Lambda (IGL) locus from 6 ethnically diverse individuals. In addition, we generated 10 diploid assemblies of IGL from a diverse cohort of individuals utilizing IG-Cap. From these 16 individuals, we identified significant allelic diversity, including 36 novel IGLV alleles. In addition, we observed highly elevated single nucleotide variation (SNV) in IGLV genes relative to IGL intergenic and genomic background SNV density. By comparing SNV calls between our high quality assemblies and existing short read datasets from the same individuals, we show a high propensity for false-positives in the short read datasets. Finally, for the first time, we nucleotide-resolved common 5-10 Kb duplications in the IGLC region that contain functional IGLJ and IGLC genes. Together these data represent a significant advancement in our understanding of genetic variation and population diversity in the IGL locus.

The Assembled Genome of the Stroke-Prone Spontaneously Hypertensive Rat.

BACKGROUND: We report the creation and evaluation of a de novo assembly of the genome of the spontaneously hypertensive rat, the most widely used model of human cardiovascular disease. METHODS: The genome is assembled from long read sequencing (PacBio HiFi and continuous long read data [CLR]) and scaffolded with long-range structural information obtained from Bionano optical maps and proximity ligation sequencing proximity analysis of the genome. The genome assembly was polished with Illumina short reads. Completeness of the assembly was investigated using Benchmarking Universal Single Copy Orthologs analysis. The genome assembly was also evaluated with the rat reference gene set, using NCBI automated protocols. We also generated orthogonal single molecule transcript sequence reads (Iso-Seq) from 8 tissues and used them to validate the coding assembly, to annotate the assembly with RNA transcripts representing unique full length transcript isoforms for each gene and to determine whether divergences between RefSeq sequences and the assembly were attributable to assembly errors or polymorphisms. RESULTS: The assembly analysis indicates that this assembly is comparable in contiguity and completeness to the current rat reference assembly, while the use of HiFi sequencing yields an assembly that is more correct at the single base level. Synteny analysis was performed to uncover the extent of synteny and the presence and distribution of chromosomal rearrangements between the reference and this assembly. CONCLUSION: The resulting genome assembly is reference quality and captures significant structural variation.

The effect of induction immunosuppression for kidney transplant on the latent HIV reservoir.

The HIV latent viral reservoir (LVR) remains a major challenge in the effort to find a cure for HIV. There is interest in lymphocyte-depleting agents, used in solid organ and bone marrow transplantation to reduce the LVR. This study evaluated the LVR and T cell receptor repertoire in HIV-infected kidney transplant recipients using intact proviral DNA assay and T cell receptor sequencing in patients receiving lymphocyte-depleting or lymphocyte-nondepleting immunosuppression induction therapy. CD4+ T cells and intact and defective provirus frequencies decreased following lymphocyte-depleting induction therapy but rebounded to near baseline levels within 1 year after induction. In contrast, these biomarkers were relatively stable over time in the lymphocyte-nondepleting group. The lymphocyte-depleting group had early TCRß repertoire turnover and newly detected and expanded clones compared with the lymphocyte-nondepleting group. No differences were observed in TCRß clonality and repertoire richness between groups. These findings suggest that, even with significant decreases in the overall size of the circulating LVR, the reservoir can be reconstituted in a relatively short period of time. These results, while from a relatively unique population, suggest that curative strategies aimed at depleting the HIV LVR will need to achieve specific and durable levels of HIV-infected T cell depletion.

Reemergence of pathogenic, autoantibody-producing B cell clones in myasthenia gravis following B cell depletion therapy.

Myasthenia gravis (MG) is an autoantibody-mediated autoimmune disorder of the neuromuscular junction. A small subset of patients (<10%) with MG, have autoantibodies targeting muscle-specific tyrosine kinase (MuSK). MuSK MG patients respond well to CD20-mediated B cell depletion therapy (BCDT); most achieve complete stable remission. However, relapse often occurs. To further understand the immunomechanisms underlying relapse, we studied autoantibody-producing B cells over the course of BCDT. We developed a fluorescently labeled antigen to enrich for MuSK-specific B cells, which was validated with a novel Nalm6 cell line engineered to express a human MuSK-specific B cell receptor. B cells (≅ 2.6 million) from 12 different samples collected from nine MuSK MG patients were screened for MuSK specificity. We successfully isolated two MuSK-specific IgG4 subclass-expressing plasmablasts from two of these patients, who were experiencing a relapse after a BCDT-induced remission. Human recombinant MuSK mAbs were then generated to validate binding specificity and characterize their molecular properties. Both mAbs were strong MuSK binders, they recognized the Ig1-like domain of MuSK, and showed pathogenic capacity when tested in an acetylcholine receptor (AChR) clustering assay. The presence of persistent clonal relatives of these MuSK-specific B cell clones was investigated through B cell receptor repertoire tracing of 63,977 unique clones derived from longitudinal samples collected from these two patients. Clonal variants were detected at multiple timepoints spanning more than five years and reemerged after BCDT-mediated remission, predating disease relapse by several months. These findings demonstrate that a reservoir of rare pathogenic MuSK autoantibody-expressing B cell clones survive BCDT and reemerge into circulation prior to manifestation of clinical relapse. Overall, this study provides both a mechanistic understanding of MuSK MG relapse and a valuable candidate biomarker for relapse prediction.

Reference long-read isoform-aware transcriptomes of 4 human peripheral blood lymphocyte subsets.

Long-read sequencing technologies such as isoform sequencing can generate highly accurate sequences of full-length mRNA transcript isoforms. Such long-read transcriptomics may be especially useful in investigations of lymphocyte functional plasticity as it relates to human health and disease. However, no long-read isoform-aware reference transcriptomes of human circulating lymphocytes are readily available despite being valuable as benchmarks in a variety of transcriptomic studies. To begin to fill this gap, we purified 4 lymphocyte populations (CD4+ T, CD8+ T, NK, and Pan B cells) from the peripheral blood of a healthy male donor and obtained high-quality RNA (RIN > 8) for isoform sequencing and parallel RNA-Seq analyses. Many novel polyadenylated transcript isoforms, supported by both isoform sequencing and RNA-Seq data, were identified within each sample. The datasets met several metrics of high quality and have been deposited to the Gene Expression Omnibus database (GSE202327, GSE202328, GSE202329) as both raw and processed files to serve as long-read reference transcriptomes for future studies of human circulating lymphocytes.

Deep oncopanel sequencing reveals within block position-dependent quality degradation in FFPE processed samples.

BACKGROUND: Clinical laboratories routinely use formalin-fixed paraffin-embedded (FFPE) tissue or cell block cytology samples in oncology panel sequencing to identify mutations that can predict patient response to targeted therapy. To understand the technical error due to FFPE processing, a robustly characterized diploid cell line was used to create FFPE samples with four different pre-tissue processing formalin fixation times. A total of 96 FFPE sections were then distributed to different laboratories for targeted sequencing analysis by four oncopanels, and variants resulting from technical error were identified. RESULTS: Tissue sections that fail more frequently show low cellularity, lower than recommended library preparation DNA input, or target sequencing depth. Importantly, sections from block surfaces are more likely to show FFPE-specific errors, akin to "edge effects" seen in histology, while the inner samples display no quality degradation related to fixation time. CONCLUSIONS: To assure reliable results, we recommend avoiding the block surface portion and restricting mutation detection to genomic regions of high confidence.

mRatBN7.2: familiar and unfamiliar features of a new rat genome reference assembly.

Rat genomic tools have been slower to emerge than for those of humans and mice and have remained less thorough and comprehensive. The arrival of a new and improved rat reference genome, mRatBN7.2, in late 2020 is a welcome event. This assembly, like predecessor rat reference assemblies, is derived from an inbred Brown Norway rat. In this "user" survey we hope to provide other users of this assembly some insight into its characteristics and some assessment of its improvements as well as a few caveats that arise from the unique aspects of this assembly. mRatBN7.2 was generated by the Wellcome Sanger Institute as part of the large Vertebrate Genomes Project. This rat assembly has now joined human, mouse, chicken, and zebrafish in the National Center for Biotechnology Information (NCBI)'s Genome Reference Consortium, which provides ongoing curation of the assembly. Here we examine the technical procedures by which the assembly was created and assess how this assembly constitutes an improvement over its predecessor. We also indicate the technical limitations affecting the assembly, providing illustrations of how these limitations arise and the impact that results for this reference assembly.