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.

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.

Orientation Regulation of Class-switch Recombination in Human B Cells.

We developed a linear amplification-mediated high-throughput genome-wide translocation sequencing method to profile Ig class-switch recombination (CSR) in human B cells in an unbiased and quantitative manner. This enables us to characterize CSR junctions resulting from either deletional recombination or inversion for each Ig class/subclass. Our data showed that more than 90% of CSR junctions detected in peripheral blood in healthy control subjects were due to deletional recombination. We further identified two major CSR junction signatures/patterns in human B cells. Signature 1 consists of recombination junctions resulting from both IgG and IgA switching, with a dominance of Sµ-Sγ junctions (72%) and deletional recombination (87%). Signature 2 is contributed mainly by Sµ-Sα junctions (96%), and these junctions were almost all due to deletional recombination (99%) and were characterized by longer microhomologies. CSR junctions identified in healthy individuals can be assigned to both signatures but with a dominance of signature 1, whereas almost all CSR junctions found in patients with defects in DNA-PKcs or Artemis, two classical nonhomologous end joining (c-NHEJ) factors, align with signature 2. Thus, signature 1 may represent c-NHEJ activity during CSR, whereas signature 2 is associated with microhomology-mediated alternative end joining in the absence of the studied c-NHEJ factors. Our findings suggest that in human B cells, the efficiency of the c-NHEJ machinery and the features of switch regions are crucial for the regulation of CSR orientation. Finally, our high-throughput method can also be applied to study the mechanism of rare types of recombination, such as switching to IgD and locus suicide switching.

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.

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.

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.

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.

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.

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.

A model for postzygotic mosaicisms quantifies the allele fraction drift, mutation rate, and contribution to de novo mutations.

The allele fraction (AF) distribution, occurrence rate, and evolutionary contribution of postzygotic single-nucleotide mosaicisms (pSNMs) remain largely unknown. In this study, we developed a mathematical model to describe the accumulation and AF drift of pSNMs during the development of multicellular organisms. By applying the model, we quantitatively analyzed two large-scale data sets of pSNMs identified from human genomes. We found that the postzygotic mutation rate per cell division during early embryogenesis, especially during the first cell division, was higher than the average mutation rate in either male or female gametes. We estimated that the stochastic cell death rate per cell cleavage during human embryogenesis was ∼5%, and parental pSNMs occurring during the first three cell divisions contributed to ∼10% of the de novo mutations observed in children. We further demonstrated that the genomic profiles of pSNMs could be used to measure the divergence distance between tissues. Our results highlight the importance of pSNMs in estimating recurrence risk and clarified the quantitative relationship between postzygotic and de novo mutations.

Distinctive types of postzygotic single-nucleotide mosaicisms in healthy individuals revealed by genome-wide profiling of multiple organs.

Postzygotic single-nucleotide mosaicisms (pSNMs) have been extensively studied in tumors and are known to play critical roles in tumorigenesis. However, the patterns and origin of pSNMs in normal organs of healthy humans remain largely unknown. Using whole-genome sequencing and ultra-deep amplicon re-sequencing, we identified and validated 164 pSNMs from 27 postmortem organ samples obtained from five healthy donors. The mutant allele fractions ranged from 1.0% to 29.7%. Inter- and intra-organ comparison revealed two distinctive types of pSNMs, with about half originating during early embryogenesis (embryonic pSNMs) and the remaining more likely to result from clonal expansion events that had occurred more recently (clonal expansion pSNMs). Compared to clonal expansion pSNMs, embryonic pSNMs had higher proportion of C>T mutations with elevated mutation rate at CpG sites. We observed differences in replication timing between these two types of pSNMs, with embryonic and clonal expansion pSNMs enriched in early- and late-replicating regions, respectively. An increased number of embryonic pSNMs were located in open chromatin states and topologically associating domains that transcribed embryonically. Our findings provide new insights into the origin and spatial distribution of postzygotic mosaicism during normal human development.

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.

MosaicHunter: accurate detection of postzygotic single-nucleotide mosaicism through next-generation sequencing of unpaired, trio, and paired samples.

Genomic mosaicism arising from postzygotic mutations has long been associated with cancer and more recently with non-cancer diseases. It has also been detected in healthy individuals including healthy parents of children affected with genetic disorders, highlighting its critical role in the origin of genetic mutations. However, most existing software for the genome-wide identification of single-nucleotide mosaicisms (SNMs) requires a paired control tissue obtained from the same individual which is often unavailable for non-cancer individuals and sometimes missing in cancer studies. Here, we present MosaicHunter (http://mosaichunter.cbi.pku.edu.cn), a bioinformatics tool that can identify SNMs in whole-genome and whole-exome sequencing data of unpaired samples without matched controls using Bayesian genotypers. We evaluate the accuracy of MosaicHunter on both simulated and real data and demonstrate that it has improved performance compared with other somatic mutation callers. We further demonstrate that incorporating sequencing data of the parents can be an effective approach to significantly improve the accuracy of detecting SNMs in an individual when a matched control sample is unavailable. Finally, MosaicHunter also has a paired mode that can take advantage of matched control samples when available, making it a useful tool for detecting SNMs in both non-cancer and cancer studies.

Postzygotic single-nucleotide mosaicisms contribute to the etiology of autism spectrum disorder and autistic traits and the origin of mutations.

The roles and characteristics of postzygotic single-nucleotide mosaicisms (pSNMs) in autism spectrum disorders (ASDs) remain unclear. In this study of the whole exomes of 2,361 families in the Simons Simplex Collection, we identified 1,248 putative pSNMs in children and 285 de novo SNPs in children with detectable parental mosaicism. Ultra-deep amplicon resequencing suggested a validation rate of 51%. Analyses of validated pSNMs revealed that missense/loss-of-function (LoF) pSNMs with a high mutant allele fraction (MAF≥ 0.2) contributed to ASD diagnoses (P = 0.022, odds ratio [OR] = 5.25), whereas missense/LoF pSNMs with a low MAF (MAF<0.2) contributed to autistic traits in male non-ASD siblings (P = 0.033). LoF pSNMs in parents were less likely to be transmitted to offspring than neutral pSNMs (P = 0.037), and missense/LoF pSNMs in parents with a low MAF were transmitted more to probands than to siblings (P = 0.016, OR = 1.45). We estimated that pSNMs in probands or de novo mutations inherited from parental pSNMs increased the risk of ASD by approximately 6%. Adding pSNMs into the transmission and de novo association test model revealed 13 new ASD risk genes. These results expand the existing repertoire of genes involved in ASD and shed new light on the contribution of genomic mosaicisms to ASD diagnoses and autistic traits.

Amplicon Resequencing Identified Parental Mosaicism for Approximately 10% of "de novo" SCN1A Mutations in Children with Dravet Syndrome.

The majority of children with Dravet syndrome (DS) are caused by de novo SCN1A mutations. To investigate the origin of the mutations, we developed and applied a new method that combined deep amplicon resequencing with a Bayesian model to detect and quantify allelic fractions with improved sensitivity. Of 174 SCN1A mutations in DS probands which were considered "de novo" by Sanger sequencing, we identified 15 cases (8.6%) of parental mosaicism. We identified another five cases of parental mosaicism that were also detectable by Sanger sequencing. Fraction of mutant alleles in the 20 cases of parental mosaicism ranged from 1.1% to 32.6%. Thirteen (65% of 20) mutations originated paternally and seven (35% of 20) maternally. Twelve (60% of 20) mosaic parents did not have any epileptic symptoms. Their mutant allelic fractions were significantly lower than those in mosaic parents with epileptic symptoms (P = 0.016). We identified mosaicism with varied allelic fractions in blood, saliva, urine, hair follicle, oral epithelium, and semen, demonstrating that postzygotic mutations could affect multiple somatic cells as well as germ cells. Our results suggest that more sensitive tools for detecting low-level mosaicism in parents of families with seemingly "de novo" mutations will allow for better informed genetic counseling.

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.

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.