Search

The rapid emergence of SARS-CoV-2 variants raised public health questions concerning the capability of diagnostic tests to detect new strains, the efficacy of vaccines, and how to map the geographical distribution of variants to understand transmission patterns and loads on healthcare resources. Next-generation sequencing (NGS) is the primary method for detecting and tracing new variants, but it is expensive, and it can take weeks before sequence data are available in public repositories. This article describes a customizable reverse transcription PCR (RT-PCR)-based genotyping approach which is significantly less expensive, accelerates reporting, and can be implemented in any lab that performs RT-PCR. Specific single-nucleotide polymorphisms (SNPs) and indels were identified which had high positive-percent agreement (PPA) and negative-percent agreement (NPA) compared to NGS for the major genotypes that circulated through September 11, 2021. Using a 48-marker panel, testing on 1,031 retrospective SARS-CoV-2 positive samples yielded a PPA and NPA ranging from 96.3 to 100% and 99.2 to 100%, respectively, for the top 10 most prevalent World Health Organization (WHO) lineages during that time. The effect of reducing the quantity of panel markers was explored, and a 16-marker panel was determined to be nearly as effective as the 48-marker panel at lineage assignment. Responding to the emergence of Omicron, a genotyping panel was developed which distinguishes Delta and Omicron using four highly specific SNPs. The results demonstrate the utility of the condensed panel to rapidly track the growing prevalence of Omicron across the US in December 2021 and January 2022.

The goal of this study is to generate and characterize a knock-in model of Pompe disease (PD) – a rare, progressive, fatal disorder primarily affecting the cardiac and musculoskeletal systems. While a murine model of PD exists, it bears a Cre/loxP induced exonic insertion of a neomycin cassette and does not completely recapitulate severe human PD - displaying nonfatal hypertrophic cardiomyopathy only late in its natural history. We therefore designed a CRISPR-Cas9 knock-in system targeting the Gaa gene to introduce the known pathogenic CRIM negative Gaa mutation c.1826insA (p.Y609*). Following optimization of our knock-in strategy in cultured murine myoblasts, we successfully generated a Gaac1826insA mouse model using a dual sgRNA with ssODN donor template approach. Whole genome sequencing and analysis of the Gaac1826insA murine model establishes that our system is highly specific for the Gaac1826 target locus and does not induce any off-target mutations or genomic rearrangements. Next, we examined GAA mRNA transcript, protein expression and enzymatic activity levels in our PD knock-in mice. Gaac1826insA mice display significantly reduced levels of GAA expression and enzymatic activity relative to wild-type mice. We performed echocardiography on Gaac1826insA mice to assess cardiac structure and function. Gaac1826insA mice exhibit early-onset, progressive cardiac hypertrophy as measured by significant increases in left ventricular wall thickness and mass index by 3 months of age. We also conducted functional tests – grip strength, inverted screen, gait analysis – on Gaac1826insA mice every 3 months to assess overall motor performance. Gaac1826insA mice display impaired motor strength and coordination relative to wild-type mice. Altogether, our results demonstrate that the Gaac1826insA murine model recapitulates human infantile-onset Pompe disease and is better suited for evaluation of therapeutic strategies such as genome correction.

The rapid emergence of new SARS-CoV-2 variants raises a number of public health questions including the capability of diagnostic tests to detect new strains, the efficacy of vaccines, and how to map the geographical distribution of variants to better understand patterns of transmission and possible load on healthcare resources. Next-Generation Sequencing (NGS) is the primary method for detecting and tracing the emergence of new variants, but it is expensive, and it can take weeks before sequence data is available in public repositories. Here, we describe a Polymerase Chain Reaction (PCR)-based genotyping approach that is significantly less expensive, accelerates reporting on SARS-CoV-2 variants, and can be implemented in any testing lab performing PCR.Specific Single Nucleotide Polymorphisms (SNPs) and indels are identified that have high positive percent agreement (PPA) and negative percent agreement (NPA) compared to NGS for the major genotypes that circulated in 2021. Using a 48-marker panel, testing on 1,128 retrospective samples yielded a PPA and NPA in the 96.3 to 100% and 99.2 to 100% range, respectively, for the top 10 most prevalent lineages. The effect on PPA and NPA of reducing the number of panel markers was also explored.In addition, with the emergence of Omicron, we also developed an Omicron genotyping panel that distinguishes the Delta and Omicron variants using four (4) highly specific SNPs. Data from testing demonstrates the capability to use the panel to rapidly track the growing prevalence of the Omicron variant in the United States in December 2021.

SUMMARY Nuclear factor κB (NF-κB) RelA is the potent transcriptional activator of inflammatory response genes. We stringently defined a list of direct RelA target genes by integrating physical (chromatin immunoprecipitation sequencing [ChIP-seq]) and functional (RNA sequencing [RNA-seq] in knockouts) datasets. We then dissected each gene’s regulatory strategy by testing RelA variants in a primary-cell genetic-complementation assay. All endogenous target genes require RelA to make DNA-base-specific contacts, and none are activatable by the DNA binding domain alone. However, endogenous target genes differ widely in how they employ the two transactivation domains. Through model-aided analysis of the dynamic time-course data, we reveal the gene-specific synergy and redundancy of TA1 and TA2. Given that post-translational modifications control TA1 activity and intrinsic affinity for coactivators determines TA2 activity, the differential TA logics suggests context-dependent versus context-independent control of endogenous RelA-target genes. Although some inflammatory initiators appear to require co-stimulatory TA1 activation, inflammatory resolvers are a part of the NF-κB RelA core response., Graphical Abstract, In Brief Ngo et al. developed a genetic complementation system for NF-κB RelA that reveals that NF-κB target-gene selection requires high-affinity RelA binding and transcriptional activation domains for gene induction. The synergistic and redundant functions of two transactivation domains define pro-inflammatory and inflammation-response genes.

Targeted deletion of BAFF causes severe deficiency of splenic B cells. BAFF-R is commonly thought to signal to nuclear factor κ-light-chain-enhancer of activated B cells (NF-κB)-inducing kinase dependent noncanonical NF-κB RelB. However, RelB-deficient mice have normal B-cell numbers. Recent studies showed that BAFF also signals to the canonical NF-κB pathway, and we found that both RelB and cRel are persistently activated, suggesting BAFF signaling coordinates both pathways to ensure robust B-cell development. Indeed, we report now that combined loss of these 2 NF-κB family members leads to impaired BAFF-mediated survival and development in vitro. Although single deletion of RelB and cRel was dispensable for normal B-cell development, double knockout mice displayed an early B-cell developmental blockade and decreased mature B cells. Despite disorganized splenic architecture in Relb(-/-)cRel(-/-) mice, generation of mixed-mouse chimeras established the developmental phenotype to be B-cell intrinsic. Together, our results indicate that BAFF signals coordinate both RelB and cRel activities to ensure survival during peripheral B-cell maturation.

Summary The regenerative capacity of the injured CNS in adult mammals is severely limited, yet axons in the peripheral nervous system (PNS) regrow, albeit to a limited extent, after injury. We reasoned that coordinate regulation of gene expression in injured neurons involving multiple pathways was central to PNS regenerative capacity. To provide a framework for revealing pathways involved in PNS axon regrowth after injury, we applied a comprehensive systems biology approach, starting with gene expression profiling of dorsal root ganglia (DRGs) combined with multi-level bioinformatic analyses and experimental validation of network predictions. We used this rubric to identify a drug that accelerates DRG neurite outgrowth in vitro and optic nerve outgrowth in vivo by inducing elements of the identified network. The work provides a functional genomics foundation for understanding neural repair and proof of the power of such approaches in tackling complex problems in nervous system biology.

The spliceosome catalyzes the removal of introns from pre-messenger RNA (mRNA) and subsequent pairing of exons with remarkable fidelity. Some exons are known to be skipped or included in the mature mRNA in a cell type- or context-dependent manner (cassette exons), thereby contributing to the diversification of the human proteome. Interestingly, splicing is initiated (and sometimes completed) co-transcriptionally. Here, we develop a kinetic mathematical modeling framework to investigate alternative co-transcriptional splicing (CTS) and, specifically, the control of cassette exons’ inclusion. We show that when splicing is co-transcriptional, default splice patterns of exon inclusion are more likely than when splicing is post-transcriptional, and that certain exons are more likely to be regulatable (i.e. cassette exons) than others, based on the exon–intron structure context. For such regulatable exons, transcriptional elongation rates may affect splicing outcomes. Within the CTS paradigm, we examine previously described hypotheses of co-operativity between splice sites of short introns (i.e. ‘intron definition’) or across short exons (i.e. ‘exon definition’), and find that models encoding these faithfully recapitulate observations in the fly and human genomes, respectively.

The emergence and mass utilization of high-throughput (HT) technologies, including sequencing technologies (genomics) and mass spectrometry (proteomics, metabolomics, lipids), has allowed geneticists, biologists, and biostatisticians to bridge the gap between genotype and phenotype on a massive scale. These new technologies have brought rapid advances in our understanding of cell biology, evolutionary history, microbial environments, and are increasingly providing new insights and applications towards clinical care and personalized medicine. Areas covered: The very success of this industry also translates into daunting big data challenges for researchers and institutions that extend beyond the traditional academic focus of algorithms and tools. The main obstacles revolve around analysis provenance, data management of massive datasets, ease of use of software, interpretability and reproducibility of results. Expert commentary: The authors review the challenges associated with implementing bioinformatics best practices in a large-scale setting, and highlight the opportunity for establishing bioinformatics pipelines that incorporate data tracking and auditing, enabling greater consistency and reproducibility for basic research, translational or clinical settings.

Genetics studies have revealed a genetic contribution to susceptibility for common or sporadic forms of neurodegenerative disease such as Alzheimer disease (AD), frontotemporal dementia (FTD), and progressive supra-nuclear palsy (PSP, a syndrome characterized by oculomotor and gait abnormalities). In AD, early genetic mapping approaches have identified rare variants in genes such as APP, PSEN1, and PSEN2that cause familial, early-onset forms.1 APOE was also pinpointed as a late-onset AD susceptibility gene.2 Genome-wide association studies3–5 (GWAS) targeted toward common variants in primarily European populations have identified many variants associated with AD, most clearly near APOE but also consistently near ABCA7, BIN1, CLU, CR1, PICALM, SORL1, and other genes. Next-generation sequencing approaches have also found rare variants with strong effect in the MAPT and TREM2 genes.6,7 In FTD, the most frequently observed mutations in familial cases occur in C9ORF72, GRN, MAPT, TARDBP, and other genes.8 In sporadic cases, a haplotype variant on the long arm of chromosome 17 has been repeatedly associated with PSP.9–11 In addition, GWAS have been performed for sporadic cases of FTD, identifying associated single-nucleotide polymorphisms (SNPs) near TMEM106B12 and BTNL2/HLA-DRA/HLA-DRB5 and RAB38/CTSC,13 as well as for PSP, identifying associated SNPs near MAPT, EIF2AK3, STX6, and MOBP.11 Despite progress in understanding the genetics of neurodegenerative diseases, known genetic risk factors cannot explain a large portion of the heritability of these diseases. For example, in AD, all common variants (including known and unknown risk variants) have been predicted to account for less than 25% of disease variance,14 and known high-penetrance rare variants account for few cases, collectively totaling only a fraction of the estimated 58% to 79% heritability of AD.15 Some of this missing heritability may be due to a blind spot in conventional genetic studies to date. A moderately rare variant with moderate effect size would be too uncommon to be tagged by a standard genotyping array and have too small of an effect to be detected by linkage or genome sequencing in practical sample sizes. The exome array bridges this gap by genotyping at low cost more than 200 000 coding variants identified through sequencing studies (Figure 1). This approach has been applied to phenotypes such as insulin homeostasis,16 bronchopulmonary dysplasia,17 and heart disease.18,19 For AD, Chung et al20 recently reported an exome array study in Korean participants that found an association with APOE, APOC1, and TOMM40 variants (near the APOE locus) but did not identify novel genetic variants. Herein, we report findings from the application of the exome array to the multiancestral Genetic Investigation in Frontotemporal Dementia and Alzheimer’s Disease (GIFT) Study cohort to determine the contribution of low-frequency coding variants to susceptibility to sporadic AD, PSP, and FTD. Figure 1 Comparison of the Exome Array and Related Genotyping and Sequencing Technologies

Understanding the functions of multi-cellular organs in terms of the molecular networks within each cell is an important step in the quest to predict phenotype from genotype. B-lymphocyte population dynamics, which are predictive of immune response and vaccine effectiveness, are determined by individual cells undergoing division or death seemingly stochastically. Based on tracking single-cell time-lapse trajectories of hundreds of B cells, single-cell transcriptome, and immunofluorescence analyses, we constructed an agent-based multi-modular computational model to simulate lymphocyte population dynamics in terms of the molecular networks that control NF-κB signaling, the cell cycle, and apoptosis. Combining modeling and experimentation, we found that NF-κBcRel enforces the execution of a cellular decision between mutually exclusive fates by promoting survival in growing cells. But as cRel deficiency causes growing B cells to die at similar rates to non-growing cells, our analysis reveals that the phenomenological decision model of wild-type cells is rooted in a biased race of cell fates. We show that a multi-scale modeling approach allows for the prediction of dynamic organ-level physiology in terms of intra-cellular molecular networks.

When messenger RNA splicing occurs co-transcriptionally, the potential for kinetic control based on transcription dynamics is widely recognized. Indeed, perturbation studies have reported that when transcription kinetics are perturbed genetically or pharmacologically splice patterns may change. However, whether kinetic control is contributing to the control of splicing within the normal range of physiological conditions remains unknown. We examined if the kinetic determinants for co-transcriptional splicing (CTS) might be reflected in the structure and expression patterns of the genome and epigenome. To identify and then quantitatively relate multiple, simultaneous CTS determinants, we constructed a scalable mathematical model of the kinetic interplay of RNA synthesis and CTS and parameterized it with diverse next generation sequencing (NGS) data. We thus found a variety of CTS determinants encoded in vertebrate genomes and epigenomes, and that these combine variously for different groups of genes such as housekeeping versus regulated genes. Together, our findings indicate that the kinetic basis of splicing is functionally and physiologically relevant, and may meaningfully inform the analysis of genomic and epigenomic data to provide insights that are missed when relying on statistical approaches alone.

SUMMARY BAFF, an activator of the noncanonical NFκB pathway, provides critical survival signals during B cell maturation and contributes to B cell proliferation. We found that the NFκB family member RelB is required ex vivo for B cell maturation, but cRel is required for proliferation. Combined molecular network modeling and experimentation revealed Nfkb2 p100 as a pathway switch; at moderate p100 synthesis rates in maturing B cells, BAFF fully utilizes p100 to generate the RelB:p52 dimer, whereas at high synthesis rates, p100 assembles into multimeric IκBsome complexes, which BAFF neutralizes in order to potentiate cRel activity and B cell expansion. Indeed, moderation of p100 expression or disruption of IκBsome assembly circumvented the BAFF requirement for full B cell expansion. Our studies emphasize the importance of p100 in determining distinct NFκB network states during B cell biology, which causes BAFF to have context-dependent functional consequences., Graphical abstract

The transcription factor NF-κB is a regulator of inflammatory and adaptive immune responses, yet only IκBα was shown to limit NF-κB activation and inflammatory responses. We investigated another negative feedback regulator, IκBε, in the regulation of B cell proliferation and survival. Loss of IκBε resulted in increased B cell proliferation and survival in response to both antigenic and innate stimulation. NF-κB activity was elevated during late-phase activation, but the dimer composition was stimulus specific. In response to IgM, cRel dimers were elevated in IκBε-deficient cells, yet in response to LPS, RelA dimers also were elevated. The corresponding dimer-specific sequences were found in the promoters of hyperactivated genes. Using a mathematical model of the NF-κB–signaling system in B cells, we demonstrated that kinetic considerations of IκB kinase–signaling input and IκBε’s interactions with RelA- and cRel-specific dimers could account for this stimulus specificity. cRel is known to be the key regulator of B cell expansion. We found that the RelA-specific phenotype in LPS-stimulated cells was physiologically relevant: unbiased transcriptome profiling revealed that the inflammatory cytokine IL-6 was hyperactivated in IκBε−/− B cells. When IL-6R was blocked, LPS-responsive IκBε−/− B cell proliferation was reduced to near wild-type levels. Our results provide novel evidence for a critical role for immune-response functions of IκBε in B cells; it regulates proliferative capacity via at least two mechanisms involving cRel- and RelA-containing NF-κB dimers. This study illustrates the importance of kinetic considerations in understanding the functional specificity of negative-feedback regulators.

The NF-κB protein RelB controls dendritic cell (DC) maturation and may be targeted therapeutically to manipulate T cell responses in disease. Here we report that RelB promoted DC activation not as the expected RelB-p52 effector of the non-canonical NF-κB pathway, but as a RelB-p50 dimer regulated by canonical IκBs, IκBα and IκBε. IκB control of RelB minimized spontaneous maturation but enabled rapid pathogen-responsive maturation. Computational modeling of the NF-κB signaling module identified control points of this unexpected cell-type-specific regulation. Fibroblasts that were engineered accordingly showed DC-like RelB control. Canonical pathway control of RelB regulated pathogen-responsive gene expression programs. This work illustrates the potential utility of systems analyses in guiding the development of combination therapeutics for modulating DC-dependent T cell responses.

Understanding human-specific patterns of brain gene expression and regulation can provide key insights into human brain evolution and speciation. Here, we use next-generation sequencing, and Illumina and Affymetrix microarray platforms, to compare the transcriptome of human, chimpanzee, and macaque telencephalon. Our analysis reveals a predominance of genes differentially expressed within human frontal lobe and a striking increase in transcriptional complexity specific to the human lineage in the frontal lobe. In contrast, caudate nucleus gene expression is highly conserved. We also identify gene coexpression signatures related to either neuronal processes or neuropsychiatric diseases, including a human-specific module with CLOCK as its hub gene and another module enriched for neuronal morphological processes and genes coexpressed with FOXP2, a gene important for language evolution. These data demonstrate that transcriptional networks have undergone evolutionary remodeling even within a given brain region, providing a window through which to view the foundation of uniquely human cognitive capacities. © 2012 Elsevier Inc.

An enigmatic feature of age-related neurodegenerative diseases is that they seldom, if ever, are fully manifested in nonhuman species under natural conditions. The neurodegenerative tauopathies are typified by the intracellular aggregation of hyperphosphorylated microtubule-associated protein tau (MAPT) and the dysfunction and death of affected neurons. We document the first case of tauopathy with paired helical filaments in an aged chimpanzee (Pan troglodytes). Pathologic forms of tau in neuronal somata, neuropil threads, and plaque-like clusters of neurites were histologically identified throughout the neocortex and, to a lesser degree, in allocortical and subcortical structures. Ultrastructurally, the neurofibrillary tangles consisted of tau-immunoreactive paired helical filaments with a diameter and helical periodicity indistinguishable from those seen in Alzheimer's disease. A moderate degree of Abeta deposition was present in the cerebral vasculature and, less frequently, in senile plaques. Sequencing of the exons and flanking intronic regions in the genomic MAPT locus disclosed no mutations that are associated with the known human hereditary tauopathies, nor any polymorphisms of obvious functional significance. Although the lesion profile in this chimpanzee differed somewhat from that in Alzheimer's disease, the copresence of paired helical filaments and Abeta-amyloidosis indicates that the molecular mechanisms for the pathogenesis of the two canonical Alzheimer lesions--neurofibrillary tangles and senile plaques--are present in aged chimpanzees.

SummaryGenome-wide transcriptional profiling was used to characterize the molecular underpinnings of neocortical organization in rhesus macaque, including cortical areal specialization and laminar cell-type diversity. Microarray analysis of individual cortical layers across sensorimotor and association cortices identified robust and specific molecular signatures for individual cortical layers and areas, prominently involving genes associated with specialized neuronal function. Overall, transcriptome-based relationships were related to spatial proximity, being strongest between neighboring cortical areas and between proximal layers. Primary visual cortex (V1) displayed the most distinctive gene expression compared to other cortical regions in rhesus and human, both in the specialized layer 4 as well as other layers. Laminar patterns were more similar between macaque and human compared to mouse, as was the unique V1 profile that was not observed in mouse. These data provide a unique resource detailing neocortical transcription patterns in a nonhuman primate with great similarity in gene expression to human.

RNA-seq technologies are now replacing microarrays for profiling gene expression. Here we describe a robust RNA-seq strategy for multiplex analysis of RNA samples based on deep sequencing. First, an oligo-dT linked to an adaptor sequence is used to prime cDNA synthesis. Upon solid phase selection, second strand synthesis is initiated using a random primer linked to another adaptor sequence. Finally, the library is released from the beads and amplified using a bar-coded primer together with a common primer. This method, referred to as Multiplex Analysis of PolyA-linked Sequences (MAPS), preserves strand information, permits rapid identification of potentially new polyadenylation sites, and profiles gene expression in a highly cost effective manner. We have applied this technology to determine the transcriptome response to knockdown of the RNA binding protein TLS, and compared the result to current microarray technology, demonstrating the ability of MAPS to robustly detect regulated gene expression.