Genetic Adaptation and Neandertal Admixture Shaped the Immune System of Human Populations.

Humans differ in the outcome that follows exposure to life-threatening pathogens, yet the extent of population differences in immune responses and their genetic and evolutionary determinants remain undefined. Here, we characterized, using RNA sequencing, the transcriptional response of primary monocytes from Africans and Europeans to bacterial and viral stimuli-ligands activating Toll-like receptor pathways (TLR1/2, TLR4, and TLR7/8) and influenza virus-and mapped expression quantitative trait loci (eQTLs). We identify numerous cis-eQTLs that contribute to the marked differences in immune responses detected within and between populations and a strong trans-eQTL hotspot at TLR1 that decreases expression of pro-inflammatory genes in Europeans only. We find that immune-responsive regulatory variants are enriched in population-specific signals of natural selection and show that admixture with Neandertals introduced regulatory variants into European genomes, affecting preferentially responses to viral challenges. Together, our study uncovers evolutionarily important determinants of differences in host immune responsiveness between human populations.

The complete sequence and comparative analysis of ape sex chromosomes.

Apes possess two sex chromosomes-the male-specific Y chromosome and the X chromosome, which is present in both males and females. The Y chromosome is crucial for male reproduction, with deletions being linked to infertility1. The X chromosome is vital for reproduction and cognition2. Variation in mating patterns and brain function among apes suggests corresponding differences in their sex chromosomes. However, owing to their repetitive nature and incomplete reference assemblies, ape sex chromosomes have been challenging to study. Here, using the methodology developed for the telomere-to-telomere (T2T) human genome, we produced gapless assemblies of the X and Y chromosomes for five great apes (bonobo (Pan paniscus), chimpanzee (Pan troglodytes), western lowland gorilla (Gorilla gorilla gorilla), Bornean orangutan (Pongo pygmaeus) and Sumatran orangutan (Pongo abelii)) and a lesser ape (the siamang gibbon (Symphalangus syndactylus)), and untangled the intricacies of their evolution. Compared with the X chromosomes, the ape Y chromosomes vary greatly in size and have low alignability and high levels of structural rearrangements-owing to the accumulation of lineage-specific ampliconic regions, palindromes, transposable elements and satellites. Many Y chromosome genes expand in multi-copy families and some evolve under purifying selection. Thus, the Y chromosome exhibits dynamic evolution, whereas the X chromosome is more stable. Mapping short-read sequencing data to these assemblies revealed diversity and selection patterns on sex chromosomes of more than 100 individual great apes. These reference assemblies are expected to inform human evolution and conservation genetics of non-human apes, all of which are endangered species.

Single-Cell RNA-Seq Mapping of Human Thymopoiesis Reveals Lineage Specification Trajectories and a Commitment Spectrum in T Cell Development.

The challenges in recapitulating in vivo human T cell development in laboratory models have posed a barrier to understanding human thymopoiesis. Here, we used single-cell RNA sequencing (sRNA-seq) to interrogate the rare CD34+ progenitor and the more differentiated CD34- fractions in the human postnatal thymus. CD34+ thymic progenitors were comprised of a spectrum of specification and commitment states characterized by multilineage priming followed by gradual T cell commitment. The earliest progenitors in the differentiation trajectory were CD7- and expressed a stem-cell-like transcriptional profile, but had also initiated T cell priming. Clustering analysis identified a CD34+ subpopulation primed for the plasmacytoid dendritic lineage, suggesting an intrathymic dendritic specification pathway. CD2 expression defined T cell commitment stages where loss of B cell potential preceded that of myeloid potential. These datasets delineate gene expression profiles spanning key differentiation events in human thymopoiesis and provide a resource for the further study of human T cell development.

CRISPRi screen of long non-coding RNAs identifies LINC03045 regulating glioblastoma invasion.

INTRODUCTION: Glioblastoma (GBM) invasion studies have focused on coding genes, while few studies evaluate long non-coding RNAs (lncRNAs), transcripts without protein-coding potential, for role in GBM invasion. We leveraged CRISPR-interference (CRISPRi) to evaluate invasive function of GBM-associated lncRNAs in an unbiased functional screen, characterizing and exploring the mechanism of identified candidates. METHODS: We implemented a CRISPRi lncRNA loss-of-function screen evaluating association of lncRNA knockdown (KD) with invasion capacity in Matrigel. Top screen candidates were validated using CRISPRi and oligonucleotide(ASO)-mediated knockdown in three tumor lines. Clinical relevance of candidates was assessed via The Cancer Genome Atlas(TCGA) and Genotype-Tissue Expression(GTEx) survival analysis. Mediators of lncRNA effect were identified via differential expression analysis following lncRNA KD and assessed for tumor invasion using knockdown and rescue experiments. RESULTS: Forty-eight lncRNAs were significantly associated with 33-83% decrease in invasion (p<0.01) upon knockdown. The top candidate, LINC03045, identified from effect size and p-value, demonstrated 82.7% decrease in tumor cell invasion upon knockdown, while LINC03045 expression was significantly associated with patient survival and tumor grade(p<0.0001). RNAseq analysis of LINC03045 knockdown revealed that WASF3, previously implicated in tumor invasion studies, was highly correlated with lncRNA expression, while WASF3 KD was associated with significant decrease in invasion. Finally, WASF3 overexpression demonstrated rescue of invasive function lost with LINC03045 KD. CONCLUSION: CRISPRi screening identified LINC03045, a previously unannotated lncRNA, as critical to GBM invasion. Gene expression is significantly associated with tumor grade and survival. RNA-seq and mechanistic studies suggest that this novel lncRNA may regulate invasion via WASF3.

Negative feedback control of neuronal activity by microglia.

Microglia, the brain's resident macrophages, help to regulate brain function by removing dying neurons, pruning non-functional synapses, and producing ligands that support neuronal survival1. Here we show that microglia are also critical modulators of neuronal activity and associated behavioural responses in mice. Microglia respond to neuronal activation by suppressing neuronal activity, and ablation of microglia amplifies and synchronizes the activity of neurons, leading to seizures. Suppression of neuronal activation by microglia occurs in a highly region-specific fashion and depends on the ability of microglia to sense and catabolize extracellular ATP, which is released upon neuronal activation by neurons and astrocytes. ATP triggers the recruitment of microglial protrusions and is converted by the microglial ATP/ADP hydrolysing ectoenzyme CD39 into AMP; AMP is then converted into adenosine by CD73, which is expressed on microglia as well as other brain cells. Microglial sensing of ATP, the ensuing microglia-dependent production of adenosine, and the adenosine-mediated suppression of neuronal responses via the adenosine receptor A1R are essential for the regulation of neuronal activity and animal behaviour. Our findings suggest that this microglia-driven negative feedback mechanism operates similarly to inhibitory neurons and is essential for protecting the brain from excessive activation in health and disease.

Histone serotonylation is a permissive modification that enhances TFIID binding to H3K4me3.

Chemical modifications of histones can mediate diverse DNA-templated processes, including gene transcription1-3. Here we provide evidence for a class of histone post-translational modification, serotonylation of glutamine, which occurs at position 5 (Q5ser) on histone H3 in organisms that produce serotonin (also known as 5-hydroxytryptamine (5-HT)). We demonstrate that tissue transglutaminase 2 can serotonylate histone H3 tri-methylated lysine 4 (H3K4me3)-marked nucleosomes, resulting in the presence of combinatorial H3K4me3Q5ser in vivo. H3K4me3Q5ser displays a ubiquitous pattern of tissue expression in mammals, with enrichment observed in brain and gut, two organ systems responsible for the bulk of 5-HT production. Genome-wide analyses of human serotonergic neurons, developing mouse brain and cultured serotonergic cells indicate that H3K4me3Q5ser nucleosomes are enriched in euchromatin, are sensitive to cellular differentiation and correlate with permissive gene expression, phenomena that are linked to the potentiation of TFIID4-6 interactions with H3K4me3. Cells that ectopically express a H3 mutant that cannot be serotonylated display significantly altered expression of H3K4me3Q5ser-target loci, which leads to deficits in differentiation. Taken together, these data identify a direct role for 5-HT, independent from its contributions to neurotransmission and cellular signalling, in the mediation of permissive gene expression.

Human PD-1 agonist treatment alleviates neutrophilic asthma by reprogramming T cells.

BACKGROUND: Neutrophilic asthma is associated with disease severity and corticosteroid insensitivity. Novel therapies are required to manage this life-threatening asthma phenotype. Programmed cell death protein-1 (PD-1) is a key homeostatic modulator of the immune response for T-cell effector functions. OBJECTIVE: We sought to investigate the role of PD-1 in the regulation of acute neutrophilic inflammation in a murine model of airway hyperreactivity (AHR). METHODS: House dust mite was used to induce and compare neutrophilic AHR in wild-type and PD-1 knockout mice. Then, the therapeutic potential of a human PD-1 agonist was tested in a humanized mouse model in which the PD-1 extracellular domain is entirely humanized. Single-cell RNA sequencing and flow cytometry were mainly used to investigate molecular and cellular mechanisms. RESULTS: PD-1 was highly induced on pulmonary T cells in our inflammatory model. PD-1 deficiency was associated with an increased neutrophilic AHR and high recruitment of inflammatory cells to the lungs. Consistently, PD-1 agonist treatment dampened AHR, decreased neutrophil recruitment, and modulated cytokine production in a humanized PD-1 mouse model. Mechanistically, we demonstrated at the transcriptional and protein levels that the inhibitory effect of PD-1 agonist is associated with the reprogramming of pulmonary effector T cells that showed decreased number and activation. CONCLUSIONS: PD-1 agonist treatment is efficient in dampening neutrophilic AHR and lung inflammation in a preclinical humanized mouse model.

ENL links histone acetylation to oncogenic gene expression in acute myeloid leukaemia.

Cancer cells are characterized by aberrant epigenetic landscapes and often exploit chromatin machinery to activate oncogenic gene expression programs. Recognition of modified histones by 'reader' proteins constitutes a key mechanism underlying these processes; therefore, targeting such pathways holds clinical promise, as exemplified by the development of bromodomain and extra-terminal (BET) inhibitors. We recently identified the YEATS domain as an acetyl-lysine-binding module, but its functional importance in human cancer remains unknown. Here we show that the YEATS domain-containing protein ENL, but not its paralogue AF9, is required for disease maintenance in acute myeloid leukaemia. CRISPR-Cas9-mediated depletion of ENL led to anti-leukaemic effects, including increased terminal myeloid differentiation and suppression of leukaemia growth in vitro and in vivo. Biochemical and crystal structural studies and chromatin-immunoprecipitation followed by sequencing analyses revealed that ENL binds to acetylated histone H3, and co-localizes with H3K27ac and H3K9ac on the promoters of actively transcribed genes that are essential for leukaemia. Disrupting the interaction between the YEATS domain and histone acetylation via structure-based mutagenesis reduced the recruitment of RNA polymerase II to ENL-target genes, leading to the suppression of oncogenic gene expression programs. Notably, disrupting the functionality of ENL further sensitized leukaemia cells to BET inhibitors. Together, our data identify ENL as a histone acetylation reader that regulates oncogenic transcriptional programs in acute myeloid leukaemia, and suggest that displacement of ENL from chromatin may be a promising epigenetic therapy, alone or in combination with BET inhibitors, for aggressive leukaemia.

Dynamic changes in cis-regulatory occupancy by Six1 and its cooperative interactions with distinct cofactors drive lineage-specific gene expression programs during progressive differentiation of the auditory sensory epithelium.

The transcription factor Six1 is essential for induction of sensory cell fate and formation of auditory sensory epithelium, but how it activates gene expression programs to generate distinct cell-types remains unknown. Here, we perform genome-wide characterization of Six1 binding at different stages of auditory sensory epithelium development and find that Six1-binding to cis-regulatory elements changes dramatically at cell-state transitions. Intriguingly, Six1 pre-occupies enhancers of cell-type-specific regulators and effectors before their expression. We demonstrate in-vivo cell-type-specific activity of Six1-bound novel enhancers of Pbx1, Fgf8, Dusp6, Vangl2, the hair-cell master regulator Atoh1 and a cascade of Atoh1's downstream factors, including Pou4f3 and Gfi1. A subset of Six1-bound sites carry consensus-sequences for its downstream factors, including Atoh1, Gfi1, Pou4f3, Gata3 and Pbx1, all of which physically interact with Six1. Motif analysis identifies RFX/X-box as one of the most significantly enriched motifs in Six1-bound sites, and we demonstrate that Six1-RFX proteins cooperatively regulate gene expression through binding to SIX:RFX-motifs. Six1 targets a wide range of hair-bundle regulators and late Six1 deletion disrupts hair-bundle polarity. This study provides a mechanistic understanding of how Six1 cooperates with distinct cofactors in feedforward loops to control lineage-specific gene expression programs during progressive differentiation of the auditory sensory epithelium.

In vivo CRISPR screening for novel noncoding RNA functional targets in glioblastoma models.

CRISPR (clustered regularly interspaced short palindromic repeat)-based genetic screens offer unbiased and powerful tools for systematic and specific evaluation of phenotypes associated with specific target genes. CRISPR screens have been utilized heavily in vitro to identify functional coding and noncoding genes in a large number of cell types, including glioblastoma (GB), though no prior study has described the evaluation of CRISPR screening in GB in vivo. Here, we describe a protocol for targeting and transcriptionally repressing GB-specific long noncoding RNAs (lncRNAs) by CRISPR interference (CRISPRi) system in vivo, with tumor growth in the mouse cerebral cortex. Given the target-specific parameters of each individual screen, we list general steps involved in transducing guide RNA libraries into GB tumor lines, maintaining sufficient coverage, as well as cortically injecting and subsequently isolating transduced screen tumor cell populations for analysis. Finally, in order to demonstrate the use of this technique to discern an essential lncRNA, HOTAIR, from a nonessential lncRNA, we injected a 1:1 (HOTAIR:control nonessential lncRNA knockdown) mixture of fluorescently tagged U87 GB cells into the cortex of eight mice, evaluating selective depletion of HOTAIR-tagged cells at 2 weeks of growth. Fluorescently tagged populations were analyzed via flow cytometry for hiBFP (control knockdown) and green fluorescent protein (HOTAIR knockdown), revealing 17% (p = 0.007) decrease in fluorescence associated with HOTAIR knockdown relative to control. The described in vivo CRISPR screening methodology thus appears to be an effective option for identifying noncoding (and coding) genes affecting GB growth within the mouse cortex.

Reconstituting the genus Mycobacterium.

The definition of a genus has wide-ranging implications both in terms of binomial species names and also evolutionary relationships. In recent years, the definition of the genus Mycobacterium has been debated due to the proposed split of this genus into five new genera (Mycolicibacterium, Mycolicibacter, Mycolicibacillus, Mycobacteroides and an emended Mycobacterium). Since this group of species contains many important obligate and opportunistic pathogens, it is important that any renaming of species does not cause confusion in clinical treatment as outlined by the nomen periculosum rule (56a) of the Prokaryotic Code. In this study, we evaluated the proposed and original genus boundaries for the mycobacteria, to determine if the split into five genera was warranted. By combining multiple approaches for defining genus boundaries (16S rRNA gene similarity, amino acid identity index, average nucleotide identity, alignment fraction and percentage of conserved proteins) we show that the original genus Mycobacterium is strongly supported over the proposed five-way split. Thus, we propose that the original genus label be reapplied to all species within this group, with the proposed five genera potentially used as sub-genus complex names.

Aberrant H3.3 dynamics in NAc promote vulnerability to depressive-like behavior.

Human major depressive disorder (MDD), along with related mood disorders, is among the world's greatest public health concerns; however, its pathophysiology remains poorly understood. Persistent changes in gene expression are known to promote physiological aberrations implicated in MDD. More recently, histone mechanisms affecting cell type- and regional-specific chromatin structures have also been shown to contribute to transcriptional programs related to depressive behaviors, as well as responses to antidepressants. Although much emphasis has been placed in recent years on roles for histone posttranslational modifications and chromatin-remodeling events in the etiology of MDD, it has become increasingly clear that replication-independent histone variants (e.g., H3.3), which differ in primary amino acid sequence from their canonical counterparts, similarly play critical roles in the regulation of activity-dependent neuronal transcription, synaptic connectivity, and behavioral plasticity. Here, we demonstrate a role for increased H3.3 dynamics in the nucleus accumbens (NAc)-a key limbic brain reward region-in the regulation of aberrant social stress-mediated gene expression and the precipitation of depressive-like behaviors in mice. We find that molecular blockade of these dynamics promotes resilience to chronic social stress and results in a partial renormalization of stress-associated transcriptional patterns in the NAc. In sum, our findings establish H3.3 dynamics as a critical, and previously undocumented, regulator of mood and suggest that future therapies aimed at modulating striatal histone dynamics may potentiate beneficial behavioral adaptations to negative emotional stimuli.

bacterial infection drives the expression dynamics of microRNAs and their isomiRs.

The optimal coordination of the transcriptional response of host cells to infection is essential for establishing appropriate immunological outcomes. In this context, the role of microRNAs (miRNAs)--important epigenetic regulators of gene expression--in regulating mammalian immune systems is increasingly well recognised. However, the expression dynamics of miRNAs, and that of their isoforms, in response to infection remains largely unexplored. Here, we characterized the genome-wide miRNA transcriptional responses of human dendritic cells, over time, to various mycobacteria differing in their virulence as well as to other bacteria outside the genus Mycobacterium, using small RNA-sequencing. We detected the presence of a core temporal response to infection, shared across bacteria, comprising 49 miRNAs, highlighting a set of miRNAs that may play an essential role in the regulation of basic cellular responses to stress. Despite such broadly shared expression dynamics, we identified specific elements of variation in the miRNA response to infection across bacteria, including a virulence-dependent induction of the miR-132/212 family in response to mycobacterial infections. We also found that infection has a strong impact on both the relative abundance of the miRNA hairpin arms and the expression dynamics of miRNA isoforms. That we observed broadly consistent changes in relative arm expression and isomiR distribution across bacteria suggests that this additional, internal layer of variability in miRNA responses represents an additional source of subtle miRNA-mediated regulation upon infection. Collectively, this study increases our understanding of the dynamism and role of miRNAs in response to bacterial infection, revealing novel features of their internal variability and identifying candidate miRNAs that may contribute to differences in the pathogenicity of mycobacterial infections.

Origins of shared genetic variation in African cichlids.

Cichlid fishes have evolved tremendous morphological and behavioral diversity in the waters of East Africa. Within each of the Great Lakes Tanganyika, Malawi, and Victoria, the phenomena of hybridization and retention of ancestral polymorphism explain allele sharing across species. Here, we explore the sharing of single nucleotide polymorphisms (SNPs) between the major East African cichlid assemblages. A set of approximately 200 genic and nongenic SNPs was ascertained in five Lake Malawi species and genotyped in a diverse collection of ~160 species from across Africa. We observed segregating polymorphism outside of the Malawi lineage for more than 50% of these loci; this holds similarly for genic versus nongenic SNPs, as well as for SNPs at putative CpG versus non-CpG sites. Bayesian and principal component analyses of genetic structure in the data demonstrate that the Lake Malawi endemic flock is not monophyletic and that river species have likely contributed significantly to Malawi genomes. Coalescent simulations support the hypothesis that river cichlids have transported polymorphism between lake assemblages. We observed strong genetic differentiation between Malawi lineages for approximately 8% of loci, with contributions from both genic and nongenic SNPs. Notably, more than half of these outlier loci between Malawi groups are polymorphic outside of the lake. Cichlid fishes have evolved diversity in Lake Malawi as new mutations combined with standing genetic variation shared across East Africa.

Brain diversity evolves via differences in patterning.

Differences in brain region size among species are thought to arise late in development via adaptive control over neurogenesis, as cells of previously patterned compartments proliferate, die, and/or differentiate into neurons. Here we investigate comparative brain development in ecologically distinct cichlid fishes from Lake Malawi and demonstrate that brains vary among recently evolved lineages because of early patterning. Divergence among rock-dwellers and sand-dwellers in the relative size of the telencephalon versus the thalamus is correlated with gene expression variation in a regulatory circuit (composed of six3, fezf2, shh, irx1b, and wnt1) known from model organisms to specify anterior-posterior (AP) brain polarity and position the shh-positive signaling boundary zona limitans intrathalamica (ZLI) in the forebrain. To confirm that changes in this coexpression network are sufficient to produce the differences we observe, we manipulated WNT signaling in vivo by treating rock-dwelling cichlid embryos with temporally precise doses of LiCl. Chemically treated rock-dwellers develop gene expression patterns, ZLIs, and forebrains distinct from controls and untreated conspecifics, but strongly resembling those of sand-dwellers. Notably, endemic Malawi rock- and sand-dwelling lineages are alternately fixed for an SNP in irx1b, a mediator of WNT signaling required for proper thalamus and ZLI. Together, these natural experiments in neuroanatomy, development, and genomics suggest that evolutionary changes in AP patterning establish ecologically relevant differences in the elaboration of cichlid forebrain compartments. In general, variation in developmental patterning might lay the foundations on which neurogenesis erects diverse brain architectures.

The Complete Sequence and Comparative Analysis of Ape Sex Chromosomes.

Apes possess two sex chromosomes-the male-specific Y and the X shared by males and females. The Y chromosome is crucial for male reproduction, with deletions linked to infertility. The X chromosome carries genes vital for reproduction and cognition. Variation in mating patterns and brain function among great apes suggests corresponding differences in their sex chromosome structure and evolution. However, due to their highly repetitive nature and incomplete reference assemblies, ape sex chromosomes have been challenging to study. Here, using the state-of-the-art experimental and computational methods developed for the telomere-to-telomere (T2T) human genome, we produced gapless, complete assemblies of the X and Y chromosomes for five great apes (chimpanzee, bonobo, gorilla, Bornean and Sumatran orangutans) and a lesser ape, the siamang gibbon. These assemblies completely resolved ampliconic, palindromic, and satellite sequences, including the entire centromeres, allowing us to untangle the intricacies of ape sex chromosome evolution. We found that, compared to the X, ape Y chromosomes vary greatly in size and have low alignability and high levels of structural rearrangements. This divergence on the Y arises from the accumulation of lineage-specific ampliconic regions and palindromes (which are shared more broadly among species on the X) and from the abundance of transposable elements and satellites (which have a lower representation on the X). Our analysis of Y chromosome genes revealed lineage-specific expansions of multi-copy gene families and signatures of purifying selection. In summary, the Y exhibits dynamic evolution, while the X is more stable. Finally, mapping short-read sequencing data from >100 great ape individuals revealed the patterns of diversity and selection on their sex chromosomes, demonstrating the utility of these reference assemblies for studies of great ape evolution. These complete sex chromosome assemblies are expected to further inform conservation genetics of nonhuman apes, all of which are endangered species.

Decoding the IGF1 signaling gene regulatory network behind alveologenesis from a mouse model of bronchopulmonary dysplasia.

Lung development is precisely controlled by underlying gene regulatory networks (GRN). Disruption of genes in the network can interrupt normal development and cause diseases such as bronchopulmonary dysplasia (BPD) - a chronic lung disease in preterm infants with morbid and sometimes lethal consequences characterized by lung immaturity and reduced alveolarization. Here, we generated a transgenic mouse exhibiting a moderate severity BPD phenotype by blocking IGF1 signaling in secondary crest myofibroblasts (SCMF) at the onset of alveologenesis. Using approaches mirroring the construction of the model GRN in sea urchin's development, we constructed the IGF1 signaling network underlying alveologenesis using this mouse model that phenocopies BPD. The constructed GRN, consisting of 43 genes, provides a bird's eye view of how the genes downstream of IGF1 are regulatorily connected. The GRN also reveals a mechanistic interpretation of how the effects of IGF1 signaling are transduced within SCMF from its specification genes to its effector genes and then from SCMF to its neighboring alveolar epithelial cells with WNT5A and FGF10 signaling as the bridge. Consistently, blocking WNT5A signaling in mice phenocopies BPD as inferred by the network. A comparative study on human samples suggests that a GRN of similar components and wiring underlies human BPD. Our network view of alveologenesis is transforming our perspective to understand and treat BPD. This new perspective calls for the construction of the full signaling GRN underlying alveologenesis, upon which targeted therapies for this neonatal chronic lung disease can be viably developed.

Divergence in cis-regulatory sequences surrounding the opsin gene arrays of African cichlid fishes.

BACKGROUND: Divergence within cis-regulatory sequences may contribute to the adaptive evolution of gene expression, but functional alleles in these regions are difficult to identify without abundant genomic resources. Among African cichlid fishes, the differential expression of seven opsin genes has produced adaptive differences in visual sensitivity. Quantitative genetic analysis suggests that cis-regulatory alleles near the SWS2-LWS opsins may contribute to this variation. Here, we sequence BACs containing the opsin genes of two cichlids, Oreochromis niloticus and Metriaclima zebra. We use phylogenetic footprinting and shadowing to examine divergence in conserved non-coding elements, promoter sequences, and 3'-UTRs surrounding each opsin in search of candidate cis-regulatory sequences that influence cichlid opsin expression. RESULTS: We identified 20 conserved non-coding elements surrounding the opsins of cichlids and other teleosts, including one known enhancer and a retinal microRNA. Most conserved elements contained computationally-predicted binding sites that correspond to transcription factors that function in vertebrate opsin expression; O. niloticus and M. zebra were significantly divergent in two of these. Similarly, we found a large number of relevant transcription factor binding sites within each opsin's proximal promoter, and identified five opsins that were considerably divergent in both expression and the number of transcription factor binding sites shared between O. niloticus and M. zebra. We also found several microRNA target sites within the 3'-UTR of each opsin, including two 3'-UTRs that differ significantly between O. niloticus and M. zebra. Finally, we examined interspecific divergence among 18 phenotypically diverse cichlids from Lake Malawi for one conserved non-coding element, two 3'-UTRs, and five opsin proximal promoters. We found that all regions were highly conserved with some evidence of CRX transcription factor binding site turnover. We also found three SNPs within two opsin promoters and one non-coding element that had weak association with cichlid opsin expression. CONCLUSIONS: This study is the first to systematically search the opsins of cichlids for putative cis-regulatory sequences. Although many putative regulatory regions are highly conserved across a large number of phenotypically diverse cichlids, we found at least nine divergent sequences that could contribute to opsin expression differences in cis and stand out as candidates for future functional analyses.

Mechanistic basis for chromosomal translocations at the E2A gene and its broader relevance to human B cell malignancies.

Analysis of translocation breakpoints in human B cell malignancies reveals that DNA double-strand breaks at oncogenes most frequently occur at CpG sites located within 20-600 bp fragile zones and depend on activation-induced deaminase (AID). AID requires single-stranded DNA (ssDNA) to act, but it has been unclear why or how this region transiently acquires a ssDNA state. Here, we demonstrate the ssDNA state in the 23 bp E2A fragile zone using several methods, including native bisulfite DNA structural analysis in live human pre-B cells. AID deamination within the E2A fragile zone does not require but is increased upon transcription. High C-string density, nascent RNA tails, and direct DNA sequence repeats prolong the ssDNA state of the E2A fragile zone and increase AID deamination at overlapping AID hotspots that contain the CpG sites at which breaks occur in patients. These features provide key insights into lymphoid fragile zones generally.

Creation of a Single Cell RNASeq Meta-Atlas to Define Human Liver Immune Homeostasis.

The liver is unique in both its ability to maintain immune homeostasis and in its potential for immune tolerance following solid organ transplantation. Single-cell RNA sequencing (scRNA seq) is a powerful approach to generate highly dimensional transcriptome data to understand cellular phenotypes. However, when scRNA data is produced by different groups, with different data models, different standards, and samples processed in different ways, it can be challenging to draw meaningful conclusions from the aggregated data. The goal of this study was to establish a method to combine 'human liver' scRNA seq datasets by 1) characterizing the heterogeneity between studies and 2) using the meta-atlas to define the dominant phenotypes across immune cell subpopulations in healthy human liver. Publicly available scRNA seq data generated from liver samples obtained from a combined total of 17 patients and ~32,000 cells were analyzed. Liver-specific immune cells (CD45+) were extracted from each dataset, and immune cell subpopulations (myeloid cells, NK and T cells, plasma cells, and B cells) were examined using dimensionality reduction (UMAP), differential gene expression, and ingenuity pathway analysis. All datasets co-clustered, but cell proportions differed between studies. Gene expression correlation demonstrated similarity across all studies, and canonical pathways that differed between datasets were related to cell stress and oxidative phosphorylation rather than immune-related function. Next, a meta-atlas was generated via data integration and compared against PBMC data to define gene signatures for each hepatic immune subpopulation. This analysis defined key features of hepatic immune homeostasis, with decreased expression across immunologic pathways and enhancement of pathways involved with cell death. This method for meta-analysis of scRNA seq data provides a novel approach to broadly define the features of human liver immune homeostasis. Specific pathways and cellular phenotypes described in this human liver immune meta-atlas provide a critical reference point for further study of immune mediated disease processes within the liver.