Molecular basis for differential Igk versus Igh V(D)J joining mechanisms.

In developing B cells, V(D)J recombination assembles exons encoding IgH and Igκ variable regions from hundreds of gene segments clustered across Igh and Igk loci. V, D and J gene segments are flanked by conserved recombination signal sequences (RSSs) that target RAG endonuclease1. RAG orchestrates Igh V(D)J recombination upon capturing a JH-RSS within the JH-RSS-based recombination centre1-3 (RC). JH-RSS orientation programmes RAG to scan upstream D- and VH-containing chromatin that is presented in a linear manner by cohesin-mediated loop extrusion4-7. During Igh scanning, RAG robustly utilizes only D-RSSs or VH-RSSs in convergent (deletional) orientation with JH-RSSs4-7. However, for Vκ-to-Jκ joining, RAG utilizes Vκ-RSSs from deletional- and inversional-oriented clusters8, inconsistent with linear scanning2. Here we characterize the Vκ-to-Jκ joining mechanism. Igk undergoes robust primary and secondary rearrangements9,10, which confounds scanning assays. We therefore engineered cells to undergo only primary Vκ-to-Jκ rearrangements and found that RAG scanning from the primary Jκ-RC terminates just 8 kb upstream within the CTCF-site-based Sis element11. Whereas Sis and the Jκ-RC barely interacted with the Vκ locus, the CTCF-site-based Cer element12 4 kb upstream of Sis interacted with various loop extrusion impediments across the locus. Similar to VH locus inversion7, DJH inversion abrogated VH-to-DJH joining; yet Vκ locus or Jκ inversion allowed robust Vκ-to-Jκ joining. Together, these experiments implicated loop extrusion in bringing Vκ segments near Cer for short-range diffusion-mediated capture by RC-based RAG. To identify key mechanistic elements for diffusional V(D)J recombination in Igk versus Igh, we assayed Vκ-to-JH and D-to-Jκ rearrangements in hybrid Igh-Igk loci generated by targeted chromosomal translocations, and pinpointed remarkably strong Vκ and Jκ RSSs. Indeed, RSS replacements in hybrid or normal Igk and Igh loci confirmed the ability of Igk-RSSs to promote robust diffusional joining compared with Igh-RSSs. We propose that Igk evolved strong RSSs to mediate diffusional Vκ-to-Jκ joining, whereas Igh evolved weaker RSSs requisite for modulating VH joining by RAG-scanning impediments.

TrieDedup: a fast trie-based deduplication algorithm to handle ambiguous bases in high-throughput sequencing.

BACKGROUND: High-throughput sequencing is a powerful tool that is extensively applied in biological studies. However, sequencers may produce low-quality bases, leading to ambiguous bases, 'N's. PCR duplicates introduced in library preparation are conventionally removed in genomics studies, and several deduplication tools have been developed for this purpose. Two identical reads may appear different due to ambiguous bases and the existing tools cannot address 'N's correctly or efficiently. RESULTS: Here we proposed and implemented TrieDedup, which uses the trie (prefix tree) data structure to compare and store sequences. TrieDedup can handle ambiguous base 'N's, and efficiently deduplicate at the level of raw sequences. We also reduced its memory usage by approximately 20% by implementing restrictedDict in Python. We benchmarked the performance of the algorithm and showed that TrieDedup can deduplicate reads up to 270-fold faster than pairwise comparison at a cost of 32-fold higher memory usage. CONCLUSIONS: The TrieDedup algorithm may facilitate PCR deduplication, barcode or UMI assignment, and repertoire diversity analysis of large-scale high-throughput sequencing datasets with its ultra-fast algorithm that can account for ambiguous bases due to sequencing errors.

The lipid globotriaosylceramide promotes germinal center B cell responses and antiviral immunity.

Influenza viruses escape immunity owing to rapid antigenic evolution, which requires vaccination strategies that allow for broadly protective antibody responses. We found that the lipid globotriaosylceramide (Gb3) expressed on germinal center (GC) B cells is essential for the production of high-affinity antibodies. Mechanistically, Gb3 bound and disengaged CD19 from its chaperone CD81, permitting CD19 to translocate to the B cell receptor complex to trigger signaling. Moreover, Gb3 regulated major histocompatibility complex class II expression to increase diversity of T follicular helper and GC B cells reactive with subdominant epitopes. In influenza infection, elevating Gb3, either endogenously or exogenously, promoted broadly reactive antibody responses and cross-protection. These data demonstrate that Gb3 determines the affinity and breadth of B cell immunity and has potential as a vaccine adjuvant.

The lipid Gb3 promotes germinal center B cell responses and anti-viral immunity.

Influenza viruses escape immunity due to rapid antigenic evolution, which requires vaccination strategies that allow for broadly protective antibody responses. Here, we demonstrate that the lipid globotriaosylceramide (Gb3) expressed on germinal center (GC) B cells is essential for the production of high-affinity antibodies. Mechanistically, Gb3 binds and disengages CD19 from its chaperone CD81 for subsequent translocation to the B cell receptor (BCR) complex to trigger signaling. Abundance of Gb3 amplifies the PI3-kinase/Akt/Foxo1 pathway to drive affinity maturation. Moreover, this lipid regulates MHC-II expression to increase diversity of T follicular helper (Tfh) and GC B cells reactive with subdominant epitopes. In influenza infection, Gb3 promotes broadly reactive antibody responses and cross-protection. Thus, we show that Gb3 determines affinity as well as breadth in B cell immunity and propose this lipid as novel vaccine adjuvant against viral infection. One Sentence Summary: Gb3 abundance on GC B cells selects antibodies with high affinity and broad epitope reactivities, which are cross-protective against heterologous influenza infection.

Contribution of the IGCR1 regulatory element and the 3'Igh CTCF-binding elements to regulation of Igh V(D)J recombination.

Immunoglobulin heavy chain variable region exons are assembled in progenitor-B cells, from VH, D, and JH gene segments located in separate clusters across the Igh locus. RAG endonuclease initiates V(D)J recombination from a JH-based recombination center (RC). Cohesin-mediated extrusion of upstream chromatin past RC-bound RAG presents Ds for joining to JHs to form a DJH-RC. Igh has a provocative number and organization of CTCF-binding elements (CBEs) that can impede loop extrusion. Thus, Igh has two divergently oriented CBEs (CBE1 and CBE2) in the IGCR1 element between the VH and D/JH domains, over 100 CBEs across the VH domain convergent to CBE1, and 10 clustered 3'Igh-CBEs convergent to CBE2 and VH CBEs. IGCR1 CBEs segregate D/JH and VH domains by impeding loop extrusion-mediated RAG-scanning. Downregulation of WAPL, a cohesin unloader, in progenitor-B cells neutralizes CBEs, allowing DJH-RC-bound RAG to scan the VH domain and perform VH-to-DJH rearrangements. To elucidate potential roles of IGCR1-based CBEs and 3'Igh-CBEs in regulating RAG-scanning and elucidate the mechanism of the ordered transition from D-to-JH to VH-to-DJH recombination, we tested effects of inverting and/or deleting IGCR1 or 3'Igh-CBEs in mice and/or progenitor-B cell lines. These studies revealed that normal IGCR1 CBE orientation augments RAG-scanning impediment activity and suggest that 3'Igh-CBEs reinforce ability of the RC to function as a dynamic loop extrusion impediment to promote optimal RAG scanning activity. Finally, our findings indicate that ordered V(D)J recombination can be explained by a gradual WAPL downregulation mechanism in progenitor-B cells as opposed to a strict developmental switch.

Humanized V(D)J-rearranging and TdT-expressing mouse vaccine models with physiological HIV-1 broadly neutralizing antibody precursors.

Antibody heavy chain (HC) and light chain (LC) variable region exons are assembled by V(D)J recombination. V(D)J junctional regions encode complementarity-determining-region 3 (CDR3), an antigen-contact region immensely diversified through nontemplated nucleotide additions ("N-regions") by terminal deoxynucleotidyl transferase (TdT). HIV-1 vaccine strategies seek to elicit human HIV-1 broadly neutralizing antibodies (bnAbs), such as the potent CD4-binding site VRC01-class bnAbs. Mice with primary B cells that express receptors (BCRs) representing bnAb precursors are used as vaccination models. VRC01-class bnAbs uniformly use human HC VH1-2 and commonly use human LCs Vκ3-20 or Vκ1-33 associated with an exceptionally short 5-amino-acid (5-aa) CDR3. Prior VRC01-class models had nonphysiological precursor levels and/or limited precursor diversity. Here, we describe VRC01-class rearranging mice that generate more physiological primary VRC01-class BCR repertoires via rearrangement of VH1-2, as well as Vκ1-33 and/or Vκ3-20 in association with diverse CDR3s. Human-like TdT expression in mouse precursor B cells increased LC CDR3 length and diversity and also promoted the generation of shorter LC CDR3s via N-region suppression of dominant microhomology-mediated Vκ-to-Jκ joins. Priming immunization with eOD-GT8 60mer, which strongly engages VRC01 precursors, induced robust VRC01-class germinal center B cell responses. Vκ3-20-based responses were enhanced by N-region addition, which generates Vκ3-20-to-Jκ junctional sequence combinations that encode VRC01-class 5-aa CDR3s with a critical E residue. VRC01-class-rearranging models should facilitate further evaluation of VRC01-class prime and boost immunogens. These new VRC01-class mouse models establish a prototype for the generation of vaccine-testing mouse models for other HIV-1 bnAb lineages that employ different HC or LC Vs.

RJAfinder: An automated tool for quantification of responding to joint attention behaviors in autism spectrum disorder using eye tracking data.

Deficits in responding to joint attention (RJA) are early symptoms of autism spectrum disorder (ASD). Currently, no automated tools exist for identifying and quantifying RJA behaviors. A few eye tracking studies have investigated RJA in ASD children but have produced conflicting results. In addition, little is known about the trajectory of RJA development through developmental age. Here, a new video was designed including 12 clips of an actor pointing to or looking at an object. Eye tracking technology was used to monitor RJA in three groups: 143 ASD children assessed with the Autism Diagnostic Interview-Revised (ADI-R) and the Autism Diagnostic Observation Schedule (ADOS) (4-7 years old), 113 age- and gender-matched typically developing children (TDC), and 43 typically developing adults (TDA) (19-32 years old). RJAfinder was developed in R and MATLAB to quantify RJA events from the eye tracking data. RJA events were compared among the three groups. Spearman correlation coefficients between total number of RJA events in ASD and the Social Responsiveness Scale (SRS) scores were calculated. A logistic regression model was built using the average valid sampling rate and the total number of RJA events as two predictive variables to classify ASD and TDC groups. ASD children displayed statistically significantly less RJA events than the TDC and TDA groups with medium-to-large-sized effects. ASD and TDC children both displayed more RJA events in response to pointing stimuli than to looking stimuli. Our logistic regression model predicted ASD tendency with 0.76 accuracy in the testing set. RJA ability improved more slowly between the ages of 4-7 years old in the ASD group than in the TDC group. In ASD children, RJA ability showed negative correlation with SRS total T-score as well as the scores of five subdomains. Our study provides an automated tool for quantifying RJA and insights for the study of RJA in ASD children, which may help improve ASD screening, subtyping, and behavior interventions.

Scratch-AID, a deep learning-based system for automatic detection of mouse scratching behavior with high accuracy.

Mice are the most commonly used model animals for itch research and for development of anti-itch drugs. Most laboratories manually quantify mouse scratching behavior to assess itch intensity. This process is labor-intensive and limits large-scale genetic or drug screenings. In this study, we developed a new system, Scratch-AID (Automatic Itch Detection), which could automatically identify and quantify mouse scratching behavior with high accuracy. Our system included a custom-designed videotaping box to ensure high-quality and replicable mouse behavior recording and a convolutional recurrent neural network trained with frame-labeled mouse scratching behavior videos, induced by nape injection of chloroquine. The best trained network achieved 97.6% recall and 96.9% precision on previously unseen test videos. Remarkably, Scratch-AID could reliably identify scratching behavior in other major mouse itch models, including the acute cheek model, the histaminergic model, and a chronic itch model. Moreover, our system detected significant differences in scratching behavior between control and mice treated with an anti-itch drug. Taken together, we have established a novel deep learning-based system that could replace manual quantification for mouse scratching behavior in different itch models and for drug screening.

An antibody from single human VH-rearranging mouse neutralizes all SARS-CoV-2 variants through BA.5 by inhibiting membrane fusion.

SARS-CoV-2 Omicron subvariants have generated a worldwide health crisis due to resistance to most approved SARS-CoV-2 neutralizing antibodies and evasion of vaccination-induced antibodies. To manage Omicron subvariants and prepare for new ones, additional means of isolating broad and potent humanized SARS-CoV-2 neutralizing antibodies are desirable. Here, we describe a mouse model in which the primary B cell receptor (BCR) repertoire is generated solely through V(D)J recombination of a human VH1-2 heavy chain (HC) and, substantially, a human Vκ1-33 light chain (LC). Thus, primary humanized BCR repertoire diversity in these mice derives from immensely diverse HC and LC antigen-contact CDR3 sequences generated by nontemplated junctional modifications during V(D)J recombination. Immunizing this mouse model with SARS-CoV-2 (Wuhan-Hu-1) spike protein immunogens elicited several VH1-2/Vκ1-33-based neutralizing antibodies that bound RBD in a different mode from each other and from those of many prior patient-derived VH1-2-based neutralizing antibodies. Of these, SP1-77 potently and broadly neutralized all SARS-CoV-2 variants through BA.5. Cryo-EM studies revealed that SP1-77 bound RBD away from the receptor-binding motif via a CDR3-dominated recognition mode. Lattice light-sheet microscopy-based studies showed that SP1-77 did not block ACE2-mediated viral attachment or endocytosis but rather blocked viral-host membrane fusion. The broad and potent SP1-77 neutralization activity and nontraditional mechanism of action suggest that it might have therapeutic potential. Likewise, the SP1-77 binding epitope may inform vaccine strategies. Last, the type of humanized mouse models that we have described may contribute to identifying therapeutic antibodies against future SARS-CoV-2 variants and other pathogens.

Generation of a live attenuated influenza A vaccine by proteolysis targeting.

The usefulness of live attenuated virus vaccines has been limited by suboptimal immunogenicity, safety concerns or cumbersome manufacturing processes and techniques. Here we describe the generation of a live attenuated influenza A virus vaccine using proteolysis-targeting chimeric (PROTAC) technology to degrade viral proteins via the endogenous ubiquitin-proteasome system of host cells. We engineered the genome of influenza A viruses in stable cell lines engineered for virus production to introduce a conditionally removable proteasome-targeting domain, generating fully infective PROTAC viruses that were live attenuated by the host protein degradation machinery upon infection. In mouse and ferret models, PROTAC viruses were highly attenuated and able to elicit robust and broad humoral, mucosal and cellular immunity against homologous and heterologous virus challenges. PROTAC-mediated attenuation of viruses may be broadly applicable for generating live attenuated vaccines.

Mycobacterial and Plasmodium ovale-associated destruction of the jaw bones.

OBJECTIVES: The project aims were to identify infectious mechanisms responsible for an extreme form of mandibular osteonecrosis and osteomyelitis in West African populations and test the hypothesis that Mycobacterium tuberculosis plays a pivotal role. MATERIALS AND METHODS: DNA was extracted from mandibular fragments of 9 of 19 patients previously included in a prospective study leading to the mycobacterial hypothesis. Amplified DNAs were used for preparing libraries suitable for next-generation sequencing. For comparison of the whole-genome sequencing data of the 9 patients with DNAs of both microbiota and human tissues, DIAMOND v0.9.26 was used to align sequencing reads to NCBI-nr database and MEGAN 6 for taxonomy binning and identification of Mycobacterium tuberculosis strains. RESULTS: The data show that mandibular bone fragments of all 9 patients not only contain Homo sapiens and Mycobacterium tuberculosis DNAs; they also contain DNAs of Plasmodium ovale wallikeri, Staphylococcus aureus, Staphylococcus hominis, and Prevotella P3-120/intermedia; as well as large numbers of DNAs from other infectious components. CONCLUSIONS: The data obtained provide direct evidence to support the conclusion that combinations of Mycobacterium tuberculosis, Plasmodium ovale wallikeri, and other oral bacteria are involved in this particular type of mandibular destruction in West African individuals of many ages.

Ku70 suppresses alternative end joining in G1-arrested progenitor B cells.

Classical nonhomologous end joining (C-NHEJ) repairs DNA double-strand breaks (DSBs) throughout interphase but predominates in G1 phase when homologous recombination is unavailable. Complexes containing the Ku70/80 ("Ku") and XRCC4/ligase IV (Lig4) core C-NHEJ factors are required, respectively, for sensing and joining DSBs. While XRCC4/Lig4 are absolutely required for joining RAG1/2 endonuclease ("RAG")-initiated DSBs during V(D)J recombination in G1-phase progenitor lymphocytes, cycling cells deficient for XRCC4/Lig4 also can join chromosomal DSBs by alternative end-joining (A-EJ) pathways. Restriction of V(D)J recombination by XRCC4/Lig4-mediated joining has been attributed to RAG shepherding V(D)J DSBs exclusively into the C-NHEJ pathway. Here, we report that A-EJ of DSB ends generated by RAG1/2, Cas9:gRNA, and Zinc finger endonucleases in Lig4-deficient G1-arrested progenitor B cell lines is suppressed by Ku. Thus, while diverse DSBs remain largely as free broken ends in Lig4-deficient G1-arrested progenitor B cells, deletion of Ku70 increases DSB rejoining and translocation levels to those observed in Ku70-deficient counterparts. Correspondingly, while RAG-initiated V(D)J DSB joining is abrogated in Lig4-deficient G1-arrested progenitor B cell lines, joining of RAG-generated DSBs in Ku70-deficient and Ku70/Lig4 double-deficient lines occurs through a translocation-like A-EJ mechanism. Thus, in G1-arrested, Lig4-deficient progenitor B cells are functionally end-joining suppressed due to Ku-dependent blockage of A-EJ, potentially in association with G1-phase down-regulation of Lig1. Finally, we suggest that differential impacts of Ku deficiency versus Lig4 deficiency on V(D)J recombination, neuronal apoptosis, and embryonic development results from Ku-mediated inhibition of A-EJ in the G1 cell cycle phase in Lig4-deficient developing lymphocyte and neuronal cells.

Loop extrusion mediates physiological Igh locus contraction for RAG scanning.

RAG endonuclease initiates Igh V(D)J recombination in progenitor B cells by binding a JH-recombination signal sequence (RSS) within a recombination centre (RC) and then linearly scanning upstream chromatin, presented by loop extrusion mediated by cohesin, for convergent D-RSSs1,2. The utilization of convergently oriented RSSs and cryptic RSSs is intrinsic to long-range RAG scanning3. Scanning of RAG from the DJH-RC-RSS to upstream convergent VH-RSSs is impeded by D-proximal CTCF-binding elements (CBEs)2-5. Primary progenitor B cells undergo a mechanistically undefined contraction of the VH locus that is proposed to provide distal VHs access to the DJH-RC6-9. Here we report that an inversion of the entire 2.4-Mb VH locus in mouse primary progenitor B cells abrogates rearrangement of both VH-RSSs and normally convergent cryptic RSSs, even though locus contraction still occurs. In addition, this inversion activated both the utilization of cryptic VH-RSSs that are normally in opposite orientation and RAG scanning beyond the VH locus through several convergent CBE domains to the telomere. Together, these findings imply that broad deregulation of CBE impediments in primary progenitor B cells promotes RAG scanning of the VH locus mediated by loop extrusion. We further found that the expression of wings apart-like protein homologue (WAPL)10, a cohesin-unloading factor, was low in primary progenitor B cells compared with v-Abl-transformed progenitor B cell lines that lacked contraction and RAG scanning of the VH locus. Correspondingly, depletion of WAPL in v-Abl-transformed lines activated both processes, further implicating loop extrusion in the locus contraction mechanism.

MosaicBase: A Knowledgebase of Postzygotic Mosaic Variants in Noncancer Disease-related and Healthy Human Individuals.

Mosaic variants resulting from postzygotic mutations are prevalent in the human genome and play important roles in human diseases. However, except for cancer-related variants, there is no collection of postzygotic mosaic variants in noncancer disease-related and healthy individuals. Here, we present MosaicBase, a comprehensive database that includes 6698 mosaic variants related to 266 noncancer diseases and 27,991 mosaic variants identified in 422 healthy individuals. Genomic and phenotypic information of each variant was manually extracted and curated from 383 publications. MosaicBase supports the query of variants with Online Mendelian Inheritance in Man (OMIM) entries, genomic coordinates, gene symbols, or Entrez IDs. We also provide an integrated genome browser for users to easily access mosaic variants and their related annotations for any genomic region. By analyzing the variants collected in MosaicBase, we find that mosaic variants that directly contribute to disease phenotype show features distinct from those of variants in individuals with mild or no phenotypes, in terms of their genomic distribution, mutation signatures, and fraction of mutant cells. MosaicBase will not only assist clinicians in genetic counseling and diagnosis but also provide a useful resource to understand the genomic baseline of postzygotic mutations in the general human population. MosaicBase is publicly available at http://mosaicbase.com/ or http://49.4.21.8:8000.

CTCF orchestrates long-range cohesin-driven V(D)J recombinational scanning.

The RAG endonuclease initiates Igh locus V(D)J recombination in progenitor (pro)-B cells1. Upon binding a recombination centre-based JH, RAG scans upstream chromatin via loop extrusion, potentially mediated by cohesin, to locate Ds and assemble a DJH-based recombination centre2. CTCF looping factor-bound elements (CBEs) within IGCR1 upstream of Ds impede RAG scanning3-5; however, their inactivation allows scanning to proximal VHs, where additional CBEs activate rearrangement and impede scanning any further upstream5. Distal VH utilization is thought to involve diffusional access to the recombination centre following large-scale Igh locus contraction6-8. Here we test the potential of linear RAG scanning to mediate distal VH usage in G1-arrested v-Abl pro-B cell lines9, which undergo robust D-to-JH but little VH-to-DJH rearrangements, presumably owing to lack of locus contraction2,5. Through an auxin-inducible approach10, we degraded the cohesin component RAD2110-12 or CTCF12,13 in these G1-arrested lines. Degradation of RAD21 eliminated all V(D)J recombination and interactions associated with RAG scanning, except for reecombination centre-located DQ52-to-JH joining, in which synapsis occurs by diffusion2. Remarkably, while degradation of CTCF suppressed most CBE-based chromatin interactions, it promoted robust recombination centre interactions with, and robust VH-to-DJH joining of, distal VHs, with patterns similar to those of 'locus-contracted' primary pro-B cells. Thus, downmodulation of CTCF-bound scanning-impediment activity promotes cohesin-driven RAG scanning across the 2.7-Mb Igh locus.

BCR selection and affinity maturation in Peyer's patch germinal centres.

The antigen-binding variable regions of the B cell receptor (BCR) and of antibodies are encoded by exons that are assembled in developing B cells by V(D)J recombination1. The BCR repertoires of primary B cells are vast owing to mechanisms that create diversity at the junctions of V(D)J gene segments that contribute to complementarity-determining region 3 (CDR3), the region that binds antigen1. Primary B cells undergo antigen-driven BCR affinity maturation through somatic hypermutation and cellular selection in germinal centres (GCs)2,3. Although most GCs are transient3, those in intestinal Peyer's patches (PPs)-which depend on the gut microbiota-are chronic4, and little is known about their BCR repertoires or patterns of somatic hypermutation. Here, using a high-throughput assay that analyses both V(D)J segment usage and somatic hypermutation profiles, we elucidate physiological BCR repertoires in mouse PP GCs. PP GCs from different mice expand public BCR clonotypes (clonotypes that are shared between many mice) that often have canonical CDR3s in the immunoglobulin heavy chain that, owing to junctional biases during V(D)J recombination, appear much more frequently than predicted in naive B cell repertoires. Some public clonotypes are dependent on the gut microbiota and encode antibodies that are reactive to bacterial glycans, whereas others are independent of gut bacteria. Transfer of faeces from specific-pathogen-free mice to germ-free mice restored germ-dependent clonotypes, directly implicating BCR selection. We identified somatic hypermutations that were recurrently selected in such public clonotypes, indicating that affinity maturation occurs in mouse PP GCs under homeostatic conditions. Thus, persistent gut antigens select recurrent BCR clonotypes to seed chronic PP GC responses.

Somatic LINE-1 retrotransposition in cortical neurons and non-brain tissues of Rett patients and healthy individuals.

Mounting evidence supports that LINE-1 (L1) retrotransposition can occur postzygotically in healthy and diseased human tissues, contributing to genomic mosaicism in the brain and other somatic tissues of an individual. However, the genomic distribution of somatic human-specific LINE-1 (L1Hs) insertions and their potential impact on carrier cells remain unclear. Here, using a PCR-based targeted bulk sequencing approach, we profiled 9,181 somatic insertions from 20 postmortem tissues from five Rett patients and their matched healthy controls. We identified and validated somatic L1Hs insertions in both cortical neurons and non-brain tissues. In Rett patients, somatic insertions were significantly depleted in exons-mainly contributed by long genes-than healthy controls, implying that cells carrying MECP2 mutations might be defenseless against a second exonic L1Hs insertion. We observed a significant increase of somatic L1Hs insertions in the brain compared with non-brain tissues from the same individual. Compared to germline insertions, somatic insertions were less sense-depleted to transcripts, indicating that they underwent weaker selective pressure on the orientation of insertion. Our observations demonstrate that somatic L1Hs insertions contribute to genomic diversity and MeCP2 dysfunction alters their genomic patterns in Rett patients.

GPA: A Microbial Genetic Polymorphisms Assignments Tool in Metagenomic Analysis by Bayesian Estimation.

Identifying antimicrobial resistant (AMR) bacteria in metagenomics samples is essential for public health and food safety. Next-generation sequencing (NGS) technology has provided a powerful tool in identifying the genetic variation and constructing the correlations between genotype and phenotype in humans and other species. However, for complex bacterial samples, there lacks a powerful bioinformatic tool to identify genetic polymorphisms or copy number variations (CNVs) for given genes. Here we provide a Bayesian framework for genotype estimation for mixtures of multiple bacteria, named as Genetic Polymorphisms Assignments (GPA). Simulation results showed that GPA has reduced the false discovery rate (FDR) and mean absolute error (MAE) in CNV and single nucleotide variant (SNV) identification. This framework was validated by whole-genome sequencing and Pool-seq data from Klebsiella pneumoniae with multiple bacteria mixture models, and showed the high accuracy in the allele fraction detections of CNVs and SNVs in AMR genes between two populations. The quantitative study on the changes of AMR genes fraction between two samples showed a good consistency with the AMR pattern observed in the individual strains. Also, the framework together with the genome annotation and population comparison tools has been integrated into an application, which could provide a complete solution for AMR gene identification and quantification in unculturable clinical samples. The GPA package is available at https://github.com/IID-DTH/GPA-package.

Targeted resequencing of 358 candidate genes for autism spectrum disorder in a Chinese cohort reveals diagnostic potential and genotype-phenotype correlations.

Autism spectrum disorder (ASD) is a childhood neuropsychiatric disorder with a complex genetic architecture. The diagnostic potential of a targeted panel of ASD genes has only been evaluated in small cohorts to date and is especially understudied in the Chinese population. Here, we designed a capture panel with 358 genes (111 syndromic and 247 nonsyndromic) for ASD and sequenced a Chinese cohort of 539 cases evaluated with the Autism Diagnostic Interview-Revised (ADI-R) and the Autism Diagnostic Observation Schedule (ADOS) as well as 512 controls. ASD cases were found to carry significantly more ultra-rare functional variants than controls. A subset of 78 syndromic and 54 nonsyndromic genes was the most significantly associated and should be given high priority in the future screening of ASD patients. Pathogenic and likely pathogenic variants were detected in 9.5% of cases. Variants in SHANK3 and SHANK2 were the most frequent, especially in females, and occurred in 1.2% of cases. Duplications of 15q11-13 were detected in 0.8% of cases. Variants in CNTNAP2 and MEF2C were correlated with epilepsy/tics in cases. Our findings reveal the diagnostic potential of ASD genetic panel testing and new insights regarding the variant spectrum. Genotype-phenotype correlations may facilitate the diagnosis and management of ASD.