Structurally divergent and recurrently mutated regions of primate genomes.

We sequenced and assembled using multiple long-read sequencing technologies the genomes of chimpanzee, bonobo, gorilla, orangutan, gibbon, macaque, owl monkey, and marmoset. We identified 1,338,997 lineage-specific fixed structural variants (SVs) disrupting 1,561 protein-coding genes and 136,932 regulatory elements, including the most complete set of human-specific fixed differences. We estimate that 819.47 Mbp or ∼27% of the genome has been affected by SVs across primate evolution. We identify 1,607 structurally divergent regions wherein recurrent structural variation contributes to creating SV hotspots where genes are recurrently lost (e.g., CARD, C4, and OLAH gene families) and additional lineage-specific genes are generated (e.g., CKAP2, VPS36, ACBD7, and NEK5 paralogs), becoming targets of rapid chromosomal diversification and positive selection (e.g., RGPD gene family). High-fidelity long-read sequencing has made these dynamic regions of the genome accessible for sequence-level analyses within and between primate species.

Recurrent inversion polymorphisms in humans associate with genetic instability and genomic disorders.

Unlike copy number variants (CNVs), inversions remain an underexplored genetic variation class. By integrating multiple genomic technologies, we discover 729 inversions in 41 human genomes. Approximately 85% of inversions <2 kbp form by twin-priming during L1 retrotransposition; 80% of the larger inversions are balanced and affect twice as many nucleotides as CNVs. Balanced inversions show an excess of common variants, and 72% are flanked by segmental duplications (SDs) or retrotransposons. Since flanking repeats promote non-allelic homologous recombination, we developed complementary approaches to identify recurrent inversion formation. We describe 40 recurrent inversions encompassing 0.6% of the genome, showing inversion rates up to 2.7 × 10-4 per locus per generation. Recurrent inversions exhibit a sex-chromosomal bias and co-localize with genomic disorder critical regions. We propose that inversion recurrence results in an elevated number of heterozygous carriers and structural SD diversity, which increases mutability in the population and predisposes specific haplotypes to disease-causing CNVs.

Differential encoding in prefrontal cortex projection neuron classes across cognitive tasks.

Single-cell transcriptomics has been widely applied to classify neurons in the mammalian brain, while systems neuroscience has historically analyzed the encoding properties of cortical neurons without considering cell types. Here we examine how specific transcriptomic types of mouse prefrontal cortex (PFC) projection neurons relate to axonal projections and encoding properties across multiple cognitive tasks. We found that most types projected to multiple targets, and most targets received projections from multiple types, except PFC→PAG (periaqueductal gray). By comparing Ca2+ activity of the molecularly homogeneous PFC→PAG type against two heterogeneous classes in several two-alternative choice tasks in freely moving mice, we found that all task-related signals assayed were qualitatively present in all examined classes. However, PAG-projecting neurons most potently encoded choice in cued tasks, whereas contralateral PFC-projecting neurons most potently encoded reward context in an uncued task. Thus, task signals are organized redundantly, but with clear quantitative biases across cells of specific molecular-anatomical characteristics.

TBL1XR1 Mutations Drive Extranodal Lymphoma by Inducing a Pro-tumorigenic Memory Fate.

The most aggressive B cell lymphomas frequently manifest extranodal distribution and carry somatic mutations in the poorly characterized gene TBL1XR1. Here, we show that TBL1XR1 mutations skew the humoral immune response toward generating abnormal immature memory B cells (MB), while impairing plasma cell differentiation. At the molecular level, TBL1XR1 mutants co-opt SMRT/HDAC3 repressor complexes toward binding the MB cell transcription factor (TF) BACH2 at the expense of the germinal center (GC) TF BCL6, leading to pre-memory transcriptional reprogramming and cell-fate bias. Upon antigen recall, TBL1XR1 mutant MB cells fail to differentiate into plasma cells and instead preferentially reenter new GC reactions, providing evidence for a cyclic reentry lymphomagenesis mechanism. Ultimately, TBL1XR1 alterations lead to a striking extranodal immunoblastic lymphoma phenotype that mimics the human disease. Both human and murine lymphomas feature expanded MB-like cell populations, consistent with a MB-cell origin and delineating an unforeseen pathway for malignant transformation of the immune system.

Nucleoid Size Scaling and Intracellular Organization of Translation across Bacteria.

The scaling of organelles with cell size is thought to be exclusive to eukaryotes. Here, we demonstrate that similar scaling relationships hold for the bacterial nucleoid. Despite the absence of a nuclear membrane, nucleoid size strongly correlates with cell size, independent of changes in DNA amount and across various nutrient conditions. This correlation is observed in diverse bacteria, revealing a near-constant ratio between nucleoid and cell size for a given species. As in eukaryotes, the nucleocytoplasmic ratio in bacteria varies greatly among species. This spectrum of nucleocytoplasmic ratios is independent of genome size, and instead it appears linked to the average population cell size. Bacteria with different nucleocytoplasmic ratios have a cytoplasm with different biophysical properties, impacting ribosome mobility and localization. Together, our findings identify new organizational principles and biophysical features of bacterial cells, implicating the nucleocytoplasmic ratio and cell size as determinants of the intracellular organization of translation.

An Integrative Model of Cellular States, Plasticity, and Genetics for Glioblastoma.

Diverse genetic, epigenetic, and developmental programs drive glioblastoma, an incurable and poorly understood tumor, but their precise characterization remains challenging. Here, we use an integrative approach spanning single-cell RNA-sequencing of 28 tumors, bulk genetic and expression analysis of 401 specimens from the The Cancer Genome Atlas (TCGA), functional approaches, and single-cell lineage tracing to derive a unified model of cellular states and genetic diversity in glioblastoma. We find that malignant cells in glioblastoma exist in four main cellular states that recapitulate distinct neural cell types, are influenced by the tumor microenvironment, and exhibit plasticity. The relative frequency of cells in each state varies between glioblastoma samples and is influenced by copy number amplifications of the CDK4, EGFR, and PDGFRA loci and by mutations in the NF1 locus, which each favor a defined state. Our work provides a blueprint for glioblastoma, integrating the malignant cell programs, their plasticity, and their modulation by genetic drivers.

Dissecting the Causal Mechanism of X-Linked Dystonia-Parkinsonism by Integrating Genome and Transcriptome Assembly.

X-linked Dystonia-Parkinsonism (XDP) is a Mendelian neurodegenerative disease that is endemic to the Philippines and is associated with a founder haplotype. We integrated multiple genome and transcriptome assembly technologies to narrow the causal mutation to the TAF1 locus, which included a SINE-VNTR-Alu (SVA) retrotransposition into intron 32 of the gene. Transcriptome analyses identified decreased expression of the canonical cTAF1 transcript among XDP probands, and de novo assembly across multiple pluripotent stem-cell-derived neuronal lineages discovered aberrant TAF1 transcription that involved alternative splicing and intron retention (IR) in proximity to the SVA that was anti-correlated with overall TAF1 expression. CRISPR/Cas9 excision of the SVA rescued this XDP-specific transcriptional signature and normalized TAF1 expression in probands. These data suggest an SVA-mediated aberrant transcriptional mechanism associated with XDP and may provide a roadmap for layered technologies and integrated assembly-based analyses for other unsolved Mendelian disorders.

Extrafollicular B cell responses correlate with neutralizing antibodies and morbidity in COVID-19.

A wide spectrum of clinical manifestations has become a hallmark of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) COVID-19 pandemic, although the immunological underpinnings of diverse disease outcomes remain to be defined. We performed detailed characterization of B cell responses through high-dimensional flow cytometry to reveal substantial heterogeneity in both effector and immature populations. More notably, critically ill patients displayed hallmarks of extrafollicular B cell activation and shared B cell repertoire features previously described in autoimmune settings. Extrafollicular activation correlated strongly with large antibody-secreting cell expansion and early production of high concentrations of SARS-CoV-2-specific neutralizing antibodies. Yet, these patients had severe disease with elevated inflammatory biomarkers, multiorgan failure and death. Overall, these findings strongly suggest a pathogenic role for immune activation in subsets of patients with COVID-19. Our study provides further evidence that targeted immunomodulatory therapy may be beneficial in specific patient subpopulations and can be informed by careful immune profiling.

Immune Monitoring of Trans-endothelial Transport by Kidney-Resident Macrophages.

Small immune complexes cause type III hypersensitivity reactions that frequently result in tissue injury. The responsible mechanisms, however, remain unclear and differ depending on target organs. Here, we identify a kidney-specific anatomical and functional unit, formed by resident macrophages and peritubular capillary endothelial cells, which monitors the transport of proteins and particles ranging from 20 to 700 kDa or 10 to 200 nm into the kidney interstitium. Kidney-resident macrophages detect and scavenge circulating immune complexes "pumped" into the interstitium via trans-endothelial transport and trigger a FcγRIV-dependent inflammatory response and the recruitment of monocytes and neutrophils. In addition, FcγRIV and TLR pathways synergistically "super-activate" kidney macrophages when immune complexes contain a nucleic acid. These data identify a physiological function of tissue-resident kidney macrophages and a basic mechanism by which they initiate the inflammatory response to small immune complexes in the kidney.

Cellulose degradation by polysaccharide monooxygenases.

Polysaccharide monooxygenases (PMOs), also known as lytic PMOs (LPMOs), enhance the depolymerization of recalcitrant polysaccharides by hydrolytic enzymes and are found in the majority of cellulolytic fungi and actinomycete bacteria. For more than a decade, PMOs were incorrectly annotated as family 61 glycoside hydrolases (GH61s) or family 33 carbohydrate-binding modules (CBM33s). PMOs have an unusual surface-exposed active site with a tightly bound Cu(II) ion that catalyzes the regioselective hydroxylation of crystalline cellulose, leading to glycosidic bond cleavage. The genomes of some cellulolytic fungi contain more than 20 genes encoding cellulose-active PMOs, suggesting a diversity of biological activities. PMOs show great promise in reducing the cost of conversion of lignocellulosic biomass to fermentable sugars; however, many questions remain about their reaction mechanism and biological function. This review addresses, in depth, the structural and mechanistic aspects of oxidative depolymerization of cellulose by PMOs and considers their biological function and phylogenetic diversity.

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.

CHD4 is recruited by GATA4 and NKX2-5 to repress noncardiac gene programs in the developing heart.

The nucleosome remodeling and deacetylase (NuRD) complex is one of the central chromatin remodeling complexes that mediates gene repression. NuRD is essential for numerous developmental events, including heart development. Clinical and genetic studies have provided direct evidence for the role of chromodomain helicase DNA-binding protein 4 (CHD4), the catalytic component of NuRD, in congenital heart disease (CHD), including atrial and ventricular septal defects. Furthermore, it has been demonstrated that CHD4 is essential for mammalian cardiomyocyte formation and function. A key unresolved question is how CHD4/NuRD is localized to specific cardiac target genes, as neither CHD4 nor NuRD can directly bind DNA. Here, we coupled a bioinformatics-based approach with mass spectrometry analyses to demonstrate that CHD4 interacts with the core cardiac transcription factors GATA4, NKX2-5, and TBX5 during embryonic heart development. Using transcriptomics and genome-wide occupancy data, we characterized the genomic landscape of GATA4, NKX2-5, and TBX5 repression and defined the direct cardiac gene targets of the GATA4-CHD4, NKX2-5-CHD4, and TBX5-CHD4 complexes. These data were used to identify putative cis-regulatory elements controlled by these complexes. We genetically interrogated two of these silencers in vivo: Acta1 and Myh11 We show that deletion of these silencers leads to inappropriate skeletal and smooth muscle gene misexpression, respectively, in the embryonic heart. These results delineate how CHD4/NuRD is localized to specific cardiac loci and explicates how mutations in the broadly expressed CHD4 protein lead to cardiac-specific disease states.

A Hyper-IgM Syndrome Mutation in Activation-Induced Cytidine Deaminase Disrupts G-Quadruplex Binding and Genome-wide Chromatin Localization.

Precise targeting of activation-induced cytidine deaminase (AID) to immunoglobulin (Ig) loci promotes antibody class switch recombination (CSR) and somatic hypermutation (SHM), whereas AID targeting of non-Ig loci can generate oncogenic DNA lesions. Here, we examined the contribution of G-quadruplex (G4) nucleic acid structures to AID targeting in vivo. Mice bearing a mutation in Aicda (AIDG133V) that disrupts AID-G4 binding modeled the pathology of hyper-IgM syndrome patients with an orthologous mutation, lacked CSR and SHM, and had broad defects in genome-wide AIDG133V chromatin localization. Genome-wide analyses also revealed that wild-type AID localized to MHCII genes, and AID expression correlated with decreased MHCII expression in germinal center B cells and diffuse large B cell lymphoma. Our findings indicate a crucial role for G4 binding in AID targeting and suggest that AID activity may extend beyond Ig loci to regulate the expression of genes relevant to the physiology and pathology of activated B cells.

Multiplying messages LRRK beneath Parkinson disease.

Parkinson Disease (PD) is a progressive neurodegenerative disorder with limited therapeutic options. In this issue of Cell, Martin et al. link PD protein leucine-rich repeat kinase 2 (LRRK2) to abnormalities of translational control, a pathogenic mechanism implicated in an increasing number of CNS neurodegenerative diseases, as well as in normal aging.

Reconstructing and reprogramming the tumor-propagating potential of glioblastoma stem-like cells.

Developmental fate decisions are dictated by master transcription factors (TFs) that interact with cis-regulatory elements to direct transcriptional programs. Certain malignant tumors may also depend on cellular hierarchies reminiscent of normal development but superimposed on underlying genetic aberrations. In glioblastoma (GBM), a subset of stem-like tumor-propagating cells (TPCs) appears to drive tumor progression and underlie therapeutic resistance yet remain poorly understood. Here, we identify a core set of neurodevelopmental TFs (POU3F2, SOX2, SALL2, and OLIG2) essential for GBM propagation. These TFs coordinately bind and activate TPC-specific regulatory elements and are sufficient to fully reprogram differentiated GBM cells to "induced" TPCs, recapitulating the epigenetic landscape and phenotype of native TPCs. We reconstruct a network model that highlights critical interactions and identifies candidate therapeutic targets for eliminating TPCs. Our study establishes the epigenetic basis of a developmental hierarchy in GBM, provides detailed insight into underlying gene regulatory programs, and suggests attendant therapeutic strategies. PAPERCLIP:

Increased mutation and gene conversion within human segmental duplications.

Single-nucleotide variants (SNVs) in segmental duplications (SDs) have not been systematically assessed because of the limitations of mapping short-read sequencing data1,2. Here we constructed 1:1 unambiguous alignments spanning high-identity SDs across 102 human haplotypes and compared the pattern of SNVs between unique and duplicated regions3,4. We find that human SNVs are elevated 60% in SDs compared to unique regions and estimate that at least 23% of this increase is due to interlocus gene conversion (IGC) with up to 4.3 megabase pairs of SD sequence converted on average per human haplotype. We develop a genome-wide map of IGC donors and acceptors, including 498 acceptor and 454 donor hotspots affecting the exons of about 800 protein-coding genes. These include 171 genes that have 'relocated' on average 1.61 megabase pairs in a subset of human haplotypes. Using a coalescent framework, we show that SD regions are slightly evolutionarily older when compared to unique sequences, probably owing to IGC. SNVs in SDs, however, show a distinct mutational spectrum: a 27.1% increase in transversions that convert cytosine to guanine or the reverse across all triplet contexts and a 7.6% reduction in the frequency of CpG-associated mutations when compared to unique DNA. We reason that these distinct mutational properties help to maintain an overall higher GC content of SD DNA compared to that of unique DNA, probably driven by GC-biased conversion between paralogous sequences5,6.

A median fin derived from the lateral plate mesoderm and the origin of paired fins.

The development of paired appendages was a key innovation during evolution and facilitated the aquatic to terrestrial transition of vertebrates. Largely derived from the lateral plate mesoderm (LPM), one hypothesis for the evolution of paired fins invokes derivation from unpaired median fins via a pair of lateral fin folds located between pectoral and pelvic fin territories1. Whilst unpaired and paired fins exhibit similar structural and molecular characteristics, no definitive evidence exists for paired lateral fin folds in larvae or adults of any extant or extinct species. As unpaired fin core components are regarded as exclusively derived from paraxial mesoderm, any transition presumes both co-option of a fin developmental programme to the LPM and bilateral duplication2. Here, we identify that the larval zebrafish unpaired pre-anal fin fold (PAFF) is derived from the LPM and thus may represent a developmental intermediate between median and paired fins. We trace the contribution of LPM to the PAFF in both cyclostomes and gnathostomes, supporting the notion that this is an ancient trait of vertebrates. Finally, we observe that the PAFF can be bifurcated by increasing bone morphogenetic protein signalling, generating LPM-derived paired fin folds. Our work provides evidence that lateral fin folds may have existed as embryonic anlage for elaboration to paired fins.

Assembly of 43 human Y chromosomes reveals extensive complexity and variation.

The prevalence of highly repetitive sequences within the human Y chromosome has prevented its complete assembly to date1 and led to its systematic omission from genomic analyses. Here we present de novo assemblies of 43 Y chromosomes spanning 182,900 years of human evolution and report considerable diversity in size and structure. Half of the male-specific euchromatic region is subject to large inversions with a greater than twofold higher recurrence rate compared with all other chromosomes2. Ampliconic sequences associated with these inversions show differing mutation rates that are sequence context dependent, and some ampliconic genes exhibit evidence for concerted evolution with the acquisition and purging of lineage-specific pseudogenes. The largest heterochromatic region in the human genome, Yq12, is composed of alternating repeat arrays that show extensive variation in the number, size and distribution, but retain a 1:1 copy-number ratio. Finally, our data suggest that the boundary between the recombining pseudoautosomal region 1 and the non-recombining portions of the X and Y chromosomes lies 500 kb away from the currently established1 boundary. The availability of fully sequence-resolved Y chromosomes from multiple individuals provides a unique opportunity for identifying new associations of traits with specific Y-chromosomal variants and garnering insights into the evolution and function of complex regions of the human genome.

FXR inhibition may protect from SARS-CoV-2 infection by reducing ACE2.

Preventing SARS-CoV-2 infection by modulating viral host receptors, such as angiotensin-converting enzyme 2 (ACE2)1, could represent a new chemoprophylactic approach for COVID-19 that complements vaccination2,3. However, the mechanisms that control the expression of ACE2 remain unclear. Here we show that the farnesoid X receptor (FXR) is a direct regulator of ACE2 transcription in several tissues affected by COVID-19, including the gastrointestinal and respiratory systems. We then use the over-the-counter compound z-guggulsterone and the off-patent drug ursodeoxycholic acid (UDCA) to reduce FXR signalling and downregulate ACE2 in human lung, cholangiocyte and intestinal organoids and in the corresponding tissues in mice and hamsters. We show that the UDCA-mediated downregulation of ACE2 reduces susceptibility to SARS-CoV-2 infection in vitro, in vivo and in human lungs and livers perfused ex situ. Furthermore, we reveal that UDCA reduces the expression of ACE2 in the nasal epithelium in humans. Finally, we identify a correlation between UDCA treatment and positive clinical outcomes after SARS-CoV-2 infection using retrospective registry data, and confirm these findings in an independent validation cohort of recipients of liver transplants. In conclusion, we show that FXR has a role in controlling ACE2 expression and provide evidence that modulation of this pathway could be beneficial for reducing SARS-CoV-2 infection, paving the way for future clinical trials.

Design and testing of a humanized porcine donor for xenotransplantation.

Recent human decedent model studies1,2 and compassionate xenograft use3 have explored the promise of porcine organs for human transplantation. To proceed to human studies, a clinically ready porcine donor must be engineered and its xenograft successfully tested in nonhuman primates. Here we describe the design, creation and long-term life-supporting function of kidney grafts from a genetically engineered porcine donor transplanted into a cynomolgus monkey model. The porcine donor was engineered to carry 69 genomic edits, eliminating glycan antigens, overexpressing human transgenes and inactivating porcine endogenous retroviruses. In vitro functional analyses showed that the edited kidney endothelial cells modulated inflammation to an extent that was indistinguishable from that of human endothelial cells, suggesting that these edited cells acquired a high level of human immune compatibility. When transplanted into cynomolgus monkeys, the kidneys with three glycan antigen knockouts alone experienced poor graft survival, whereas those with glycan antigen knockouts and human transgene expression demonstrated significantly longer survival time, suggesting the benefit of human transgene expression in vivo. These results show that preclinical studies of renal xenotransplantation could be successfully conducted in nonhuman primates and bring us closer to clinical trials of genetically engineered porcine renal grafts.