Chromosome-level organization of the regulatory genome in the Drosophila nervous system.

Previous studies have identified topologically associating domains (TADs) as basic units of genome organization. We present evidence of a previously unreported level of genome folding, where distant TAD pairs, megabases apart, interact to form meta-domains. Within meta-domains, gene promoters and structural intergenic elements present in distant TADs are specifically paired. The associated genes encode neuronal determinants, including those engaged in axonal guidance and adhesion. These long-range associations occur in a large fraction of neurons but support transcription in only a subset of neurons. Meta-domains are formed by diverse transcription factors that are able to pair over long and flexible distances. We present evidence that two such factors, GAF and CTCF, play direct roles in this process. The relative simplicity of higher-order meta-domain interactions in Drosophila, compared with those previously described in mammals, allowed the demonstration that genomes can fold into highly specialized cell-type-specific scaffolds that enable megabase-scale regulatory associations.

Mechanopathology of biofilm-like Mycobacterium tuberculosis cords.

Mycobacterium tuberculosis (Mtb) cultured axenically without detergent forms biofilm-like cords, a clinical identifier of virulence. In lung-on-chip (LoC) and mouse models, cords in alveolar cells contribute to suppression of innate immune signaling via nuclear compression. Thereafter, extracellular cords cause contact-dependent phagocyte death but grow intercellularly between epithelial cells. The absence of these mechanopathological mechanisms explains the greater proportion of alveolar lesions with increased immune infiltration and dissemination defects in cording-deficient Mtb infections. Compression of Mtb lipid monolayers induces a phase transition that enables mechanical energy storage. Agent-based simulations demonstrate that the increased energy storage capacity is sufficient for the formation of cords that maintain structural integrity despite mechanical perturbation. Bacteria in cords remain translationally active despite antibiotic exposure and regrow rapidly upon cessation of treatment. This study provides a conceptual framework for the biophysics and function in tuberculosis infection and therapy of cord architectures independent of mechanisms ascribed to single bacteria.

DNA modifications impact natural transformation of Acinetobacter baumannii.

Acinetobacter baumannii is a dangerous nosocomial pathogen, especially due to its ability to rapidly acquire new genetic traits, including antibiotic resistance genes (ARG). In A. baumannii, natural competence for transformation, one of the primary modes of horizontal gene transfer (HGT), is thought to contribute to ARG acquisition and has therefore been intensively studied. However, knowledge regarding the potential role of epigenetic DNA modification(s) on this process remains lacking. Here, we demonstrate that the methylome pattern of diverse A. baumannii strains differs substantially and that these epigenetic marks influence the fate of transforming DNA. Specifically, we describe a methylome-dependent phenomenon that impacts intra- and inter-species DNA exchange by the competent A. baumannii strain A118. We go on to identify and characterize an A118-specific restriction-modification (RM) system that impairs transformation when the incoming DNA lacks a specific methylation signature. Collectively, our work contributes towards a more holistic understanding of HGT in this organism and may also aid future endeavors towards tackling the spread of novel ARGs. In particular, our results suggest that DNA exchanges between bacteria that share similar epigenomes are favored and could therefore guide future research into identifying the reservoir(s) of dangerous genetic traits for this multi-drug resistant pathogen.

tRNA-derived fragments in T lymphocyte-beta cell crosstalk and in type 1 diabetes pathogenesis in NOD mice.

AIMS/HYPOTHESIS: tRNAs play a central role in protein synthesis. Besides this canonical function, they were recently found to generate non-coding RNA fragments (tRFs) regulating different cellular activities. The aim of this study was to assess the involvement of tRFs in the crosstalk between immune cells and beta cells and to investigate their contribution to the development of type 1 diabetes. METHODS: Global profiling of the tRFs present in pancreatic islets of 4- and 8-week-old NOD mice and in extracellular vesicles released by activated CD4+ T lymphocytes was performed by small RNA-seq. Changes in the level of specific fragments were confirmed by quantitative PCR. The transfer of tRFs from immune cells to beta cells occurring during insulitis was assessed using an RNA-tagging approach. The functional role of tRFs increasing in beta cells during the initial phases of type 1 diabetes was determined by overexpressing them in dissociated islet cells and by determining the impact on gene expression and beta cell apoptosis. RESULTS: We found that the tRF pool was altered in the islets of NOD mice during the initial phases of type 1 diabetes. Part of these changes were triggered by prolonged exposure of beta cells to proinflammatory cytokines (IL-1ß, TNF-α and IFN-γ) while others resulted from the delivery of tRFs produced by CD4+ T lymphocytes infiltrating the islets. Indeed, we identified several tRFs that were enriched in extracellular vesicles from CD4+/CD25- T cells and were transferred to beta cells upon adoptive transfer of these immune cells in NOD.SCID mice. The tRFs delivered to beta cells during the autoimmune reaction triggered gene expression changes that affected the immune regulatory capacity of insulin-secreting cells and rendered the cells more prone to apoptosis. CONCLUSIONS/INTERPRETATION: Our data point to tRFs as novel players in the crosstalk between the immune system and insulin-secreting cells and suggest a potential involvement of this novel class of non-coding RNAs in type 1 diabetes pathogenesis. DATA AVAILABILITY: Sequences are available from the Gene Expression Omnibus (GEO) with accession numbers GSE242568 and GSE256343.

Variants in the degron of AFF3 are associated with intellectual disability, mesomelic dysplasia, horseshoe kidney, and epileptic encephalopathy.

The ALF transcription factor paralogs, AFF1, AFF2, AFF3, and AFF4, are components of the transcriptional super elongation complex that regulates expression of genes involved in neurogenesis and development. We describe an autosomal dominant disorder associated with de novo missense variants in the degron of AFF3, a nine amino acid sequence important for its binding to ubiquitin ligase, or with de novo deletions of this region. The sixteen affected individuals we identified, along with two previously reported individuals, present with a recognizable pattern of anomalies, which we named KINSSHIP syndrome (KI for horseshoe kidney, NS for Nievergelt/Savarirayan type of mesomelic dysplasia, S for seizures, H for hypertrichosis, I for intellectual disability, and P for pulmonary involvement), partially overlapping the AFF4-associated CHOPS syndrome. Whereas homozygous Aff3 knockout mice display skeletal anomalies, kidney defects, brain malformations, and neurological anomalies, knockin animals modeling one of the microdeletions and the most common of the missense variants identified in affected individuals presented with lower mesomelic limb deformities like KINSSHIP-affected individuals and early lethality, respectively. Overexpression of AFF3 in zebrafish resulted in body axis anomalies, providing some support for the pathological effect of increased amount of AFF3. The only partial phenotypic overlap of AFF3- and AFF4-associated syndromes and the previously published transcriptome analyses of ALF transcription factors suggest that these factors are not redundant and each contributes uniquely to proper development.

Mutation-specific pathophysiological mechanisms define different neurodevelopmental disorders associated with SATB1 dysfunction.

Whereas large-scale statistical analyses can robustly identify disease-gene relationships, they do not accurately capture genotype-phenotype correlations or disease mechanisms. We use multiple lines of independent evidence to show that different variant types in a single gene, SATB1, cause clinically overlapping but distinct neurodevelopmental disorders. Clinical evaluation of 42 individuals carrying SATB1 variants identified overt genotype-phenotype relationships, associated with different pathophysiological mechanisms, established by functional assays. Missense variants in the CUT1 and CUT2 DNA-binding domains result in stronger chromatin binding, increased transcriptional repression, and a severe phenotype. In contrast, variants predicted to result in haploinsufficiency are associated with a milder clinical presentation. A similarly mild phenotype is observed for individuals with premature protein truncating variants that escape nonsense-mediated decay, which are transcriptionally active but mislocalized in the cell. Our results suggest that in-depth mutation-specific genotype-phenotype studies are essential to capture full disease complexity and to explain phenotypic variability.

SPICA: Swiss portal for immune cell analysis.

Single-cell transcriptomics allows the study of immune cell heterogeneity at an unprecedented level of resolution. The Swiss portal for immune cell analysis (SPICA) is a web resource dedicated to the exploration and analysis of single-cell RNA-seq data of immune cells. In contrast to other single-cell databases, SPICA hosts curated, cell type-specific reference atlases that describe immune cell states at high resolution, and published single-cell datasets analysed in the context of these atlases. Additionally, users can privately analyse their own data in the context of existing atlases and contribute to the SPICA database. SPICA is available at https://spica.unil.ch.

Variants in USP48 encoding ubiquitin hydrolase are associated with autosomal dominant non-syndromic hereditary hearing loss.

Non-Syndromic Hereditary Hearing Loss (NSHHL) is a genetically heterogeneous sensory disorder with about 120 genes already associated. Through exome sequencing (ES) and data aggregation, we identified a family with six affected individuals and one unrelated NSHHL patient with predicted-to-be deleterious missense variants in USP48. We also uncovered an eighth patient presenting unilateral cochlear nerve aplasia and a de novo splice variant in the same gene. USP48 encodes a ubiquitin carboxyl-terminal hydrolase under evolutionary constraint. Pathogenicity of the variants is supported by in vitro assays that showed that the mutated proteins are unable to hydrolyze tetra-ubiquitin. Correspondingly, three-dimensional representation of the protein containing the familial missense variant is situated in a loop that might influence the binding to ubiquitin. Consistent with a contribution of USP48 to auditory function, immunohistology showed that the encoded protein is expressed in the developing human inner ear, specifically in the spiral ganglion neurons, outer sulcus, interdental cells of the spiral limbus, stria vascularis, Reissner's membrane and in the transient Kolliker's organ that is essential for auditory development. Engineered zebrafish knocked-down for usp48, the USP48 ortholog, presented with a delayed development of primary motor neurons, less developed statoacoustic neurons innervating the ears, decreased swimming velocity and circling swimming behavior indicative of vestibular dysfunction and hearing impairment. Corroboratingly, acoustic startle response assays revealed a significant decrease of auditory response of zebrafish lacking usp48 at 600 and 800 Hz wavelengths. In conclusion, we describe a novel autosomal dominant NSHHL gene through a multipronged approach combining ES, animal modeling, immunohistology and molecular assays.

Biallelic variants in NSUN6 cause an autosomal recessive neurodevelopmental disorder.

PURPOSE: 5-methylcytosine RNA modifications are driven by NSUN methyltransferases. Although variants in NSUN2 and NSUN3 were associated with neurodevelopmental diseases, the physiological role of NSUN6 modifications on transfer RNAs and messenger RNAs remained elusive. METHODS: We combined exome sequencing of consanguineous families with functional characterization to identify a new neurodevelopmental disorder gene. RESULTS: We identified 3 unrelated consanguineous families with deleterious homozygous variants in NSUN6. Two of these variants are predicted to be loss-of-function. One maps to the first exon and is predicted to lead to the absence of NSUN6 via nonsense-mediated decay, whereas we showed that the other maps to the last exon and encodes a protein that does not fold correctly. Likewise, we demonstrated that the missense variant identified in the third family has lost its enzymatic activity and is unable to bind the methyl donor S-adenosyl-L-methionine. The affected individuals present with developmental delay, intellectual disability, motor delay, and behavioral anomalies. Homozygous ablation of the NSUN6 ortholog in Drosophila led to locomotion and learning impairment. CONCLUSION: Our data provide evidence that biallelic pathogenic variants in NSUN6 cause one form of autosomal recessive intellectual disability, establishing another link between RNA modification and cognition.

Soil protist function varies with elevation in the Swiss Alps.

Protists are abundant and play key trophic functions in soil. Documenting how their trophic contributions vary across large environmental gradients is essential to understand and predict how biogeochemical cycles will be impacted by global changes. Here, using amplicon sequencing of environmental DNA in open habitat soil from 161 locations spanning 2600 m of elevation in the Swiss Alps (from 400 to 3000 m), we found that, over the whole study area, soils are dominated by consumers, followed by parasites and phototrophs. In contrast, the proportion of these groups in local communities shows large variations in relation to elevation. While there is, on average, three times more consumers than parasites at low elevation (400-1000 m), this ratio increases to 12 at high elevation (2000-3000 m). This suggests that the decrease in protist host biomass and diversity toward mountains tops impact protist functional composition. Furthermore, the taxonomic composition of protists that infect animals was related to elevation while that of protists that infect plants or of protist consumers was related to soil pH. This study provides a first step to document and understand how soil protist functions vary along the elevational gradient.

Unexpected post-glacial colonisation route explains the white colour of barn owls (Tyto alba) from the British Isles.

The climate fluctuations of the Quaternary shaped the movement of species in and out of glacial refugia. In Europe, the majority of species followed one of the described traditional postglacial recolonization routes from the southern peninsulas towards the north. Like most organisms, barn owls are assumed to have colonized the British Isles by crossing over Doggerland, a land bridge that connected Britain to northern Europe. However, while they are dark rufous in northern Europe, barn owls in the British Isles are conspicuously white, a contrast that could suggest selective forces are at play on the islands. Yet, our analysis of known candidate genes involved in coloration found no signature of selection. Instead, using whole genome sequences and species distribution modelling, we found that owls colonised the British Isles soon after the last glaciation, directly from a white coloured refugium in the Iberian Peninsula, before colonising northern Europe. They would have followed a hitherto unknown post-glacial colonization route to the Isles over a westwards path of suitable habitat in now submerged land in the Bay of Biscay, thus not crossing Doggerland. As such, they inherited the white colour of their Iberian founders and maintained it through low gene flow with the mainland that prevents the import of rufous alleles. Thus, we contend that neutral processes probably explain this contrasting white colour compared to continental owls. With the barn owl being a top predator, we expect future research will show this unanticipated route was used by other species from its paleo community.

Genomic consequences of colonisation, migration and genetic drift in barn owl insular populations of the eastern Mediterranean.

The study of insular populations was key in the development of evolutionary theory. The successful colonisation of an island depends on the geographic context, and specific characteristics of the organism and the island, but also on stochastic processes. As a result, apparently identical islands may harbour populations with contrasting histories. Here, we use whole genome sequences of 65 barn owls to investigate the patterns of inbreeding and genetic diversity of insular populations in the eastern Mediterranean Sea. We focus on Crete and Cyprus, islands with similar size, climate and distance to mainland, that provide natural replicates for a comparative analysis of the impacts of microevolutionary processes on isolated populations. We show that barn owl populations from each island have a separate origin, Crete being genetically more similar to other Greek islands and mainland Greece, and Cyprus more similar to the Levant. Further, our data show that their respective demographic histories following colonisation were also distinct. On the one hand, Crete harbours a small population and maintains very low levels of gene flow with neighbouring populations. This has resulted in low genetic diversity, strong genetic drift, increased relatedness in the population and remote inbreeding. Cyprus, on the other hand, appears to maintain enough gene flow with the mainland to avoid such an outcome. Our study provides a comparative population genomic analysis of the effects of neutral processes on a classical island-mainland model system. It provides empirical evidence for the role of stochastic processes in determining the fate of diverging isolated populations.

Anti-adipogenic signals at the onset of obesity-related inflammation in white adipose tissue.

Chronic inflammation that affects primarily metabolic organs, such as white adipose tissue (WAT), is considered as a major cause of human obesity-associated co-morbidities. However, the molecular mechanisms initiating this inflammation in WAT are poorly understood. By combining transcriptomics, ChIP-seq and modeling approaches, we studied the global early and late responses to a high-fat diet (HFD) in visceral (vWAT) and subcutaneous (scWAT) AT, the first being more prone to obesity-induced inflammation. HFD rapidly triggers proliferation of adipocyte precursors within vWAT. However, concomitant antiadipogenic signals limit vWAT hyperplastic expansion by interfering with the differentiation of proliferating adipocyte precursors. Conversely, in scWAT, residing beige adipocytes lose their oxidizing properties and allow storage of excessive fatty acids. This phase is followed by tissue hyperplastic growth and increased angiogenic signals, which further enable scWAT expansion without generating inflammation. Our data indicate that scWAT and vWAT differential ability to modulate adipocyte number and differentiation in response to obesogenic stimuli has a crucial impact on the different susceptibility to obesity-related inflammation of these adipose tissue depots.

Inhibition of G-protein signalling in cardiac dysfunction of intellectual developmental disorder with cardiac arrhythmia (IDDCA) syndrome.

BACKGROUND: Pathogenic variants of GNB5 encoding the ß5 subunit of the guanine nucleotide-binding protein cause IDDCA syndrome, an autosomal recessive neurodevelopmental disorder associated with cognitive disability and cardiac arrhythmia, particularly severe bradycardia. METHODS: We used echocardiography and telemetric ECG recordings to investigate consequences of Gnb5 loss in mouse. RESULTS: We delineated a key role of Gnb5 in heart sinus conduction and showed that Gnb5-inhibitory signalling is essential for parasympathetic control of heart rate (HR) and maintenance of the sympathovagal balance. Gnb5-/- mice were smaller and had a smaller heart than Gnb5+/+ and Gnb5+/- , but exhibited better cardiac function. Lower autonomic nervous system modulation through diminished parasympathetic control and greater sympathetic regulation resulted in a higher baseline HR in Gnb5-/- mice. In contrast, Gnb5-/- mice exhibited profound bradycardia on treatment with carbachol, while sympathetic modulation of the cardiac stimulation was not altered. Concordantly, transcriptome study pinpointed altered expression of genes involved in cardiac muscle contractility in atria and ventricles of knocked-out mice. Homozygous Gnb5 loss resulted in significantly higher frequencies of sinus arrhythmias. Moreover, we described 13 affected individuals, increasing the IDDCA cohort to 44 patients. CONCLUSIONS: Our data demonstrate that loss of negative regulation of the inhibitory G-protein signalling causes HR perturbations in Gnb5-/- mice, an effect mainly driven by impaired parasympathetic activity. We anticipate that unravelling the mechanism of Gnb5 signalling in the autonomic control of the heart will pave the way for future drug screening.

KIAA1109 Variants Are Associated with a Severe Disorder of Brain Development and Arthrogryposis.

Whole-exome and targeted sequencing of 13 individuals from 10 unrelated families with overlapping clinical manifestations identified loss-of-function and missense variants in KIAA1109 allowing delineation of an autosomal-recessive multi-system syndrome, which we suggest to name Alkuraya-Kucinskas syndrome (MIM 617822). Shared phenotypic features representing the cardinal characteristics of this syndrome combine brain atrophy with clubfoot and arthrogryposis. Affected individuals present with cerebral parenchymal underdevelopment, ranging from major cerebral parenchymal thinning with lissencephalic aspect to moderate parenchymal rarefaction, severe to mild ventriculomegaly, cerebellar hypoplasia with brainstem dysgenesis, and cardiac and ophthalmologic anomalies, such as microphthalmia and cataract. Severe loss-of-function cases were incompatible with life, whereas those individuals with milder missense variants presented with severe global developmental delay, syndactyly of 2nd and 3rd toes, and severe muscle hypotonia resulting in incapacity to stand without support. Consistent with a causative role for KIAA1109 loss-of-function/hypomorphic variants in this syndrome, knockdowns of the zebrafish orthologous gene resulted in embryos with hydrocephaly and abnormally curved notochords and overall body shape, whereas published knockouts of the fruit fly and mouse orthologous genes resulted in lethality or severe neurological defects reminiscent of the probands' features.

PamgeneAnalyzeR: open and reproducible pipeline for kinase profiling.

Protein phosphorylation--catalyzed by protein kinases-is the most common post-translational modification. It increases the functional diversity of the proteome and influences various aspects of normal physiology and can be altered in disease states. High throughput profiling of kinases is becoming an essential experimental approach to investigate their activity and this can be achieved using technologies such as PamChip® arrays provided by PamGene for kinase activity measurement. Here, we present 'pamgeneAnalyzeR', an R package developed as an alternative to the manual steps necessary to extract the data from PamChip® peptide microarrays images in a reproducible and robust manner. The extracted data can be directly used for downstream analysis. AVAILABILITY AND IMPLEMENTATION: PamgeneAnalyzeR is implemented in R and can be obtained from https://github.com/amelbek/pamgeneAnalyzeR.

A systems genetics resource and analysis of sleep regulation in the mouse.

Sleep is essential for optimal brain functioning and health, but the biological substrates through which sleep delivers these beneficial effects remain largely unknown. We used a systems genetics approach in the BXD genetic reference population (GRP) of mice and assembled a comprehensive experimental knowledge base comprising a deep "sleep-wake" phenome, central and peripheral transcriptomes, and plasma metabolome data, collected under undisturbed baseline conditions and after sleep deprivation (SD). We present analytical tools to interactively interrogate the database, visualize the molecular networks altered by sleep loss, and prioritize candidate genes. We found that a one-time, short disruption of sleep already extensively reshaped the systems genetics landscape by altering 60%-78% of the transcriptomes and the metabolome, with numerous genetic loci affecting the magnitude and direction of change. Systems genetics integrative analyses drawing on all levels of organization imply α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptor trafficking and fatty acid turnover as substrates of the negative effects of insufficient sleep. Our analyses demonstrate that genetic heterogeneity and the effects of insufficient sleep itself on the transcriptome and metabolome are far more widespread than previously reported.

Differential regulation of RNA polymerase III genes during liver regeneration.

Mouse liver regeneration after partial hepatectomy involves cells in the remaining tissue synchronously entering the cell division cycle. We have used this system and H3K4me3, Pol II and Pol III profiling to characterize adaptations in Pol III transcription. Our results broadly define a class of genes close to H3K4me3 and Pol II peaks, whose Pol III occupancy is high and stable, and another class, distant from Pol II peaks, whose Pol III occupancy strongly increases after partial hepatectomy. Pol III regulation in the liver thus entails both highly expressed housekeeping genes and genes whose expression can adapt to increased demand.

Greater topoclimatic control of above- versus below-ground communities.

Assessing the degree to which climate explains the spatial distributions of different taxonomic and functional groups is essential for anticipating the effects of climate change on ecosystems. Most effort so far has focused on above-ground organisms, which offer only a partial view on the response of biodiversity to environmental gradients. Here including both above- and below-ground organisms, we quantified the degree of topoclimatic control on the occurrence patterns of >1,500 taxa and phylotypes along a c. 3,000 m elevation gradient, by fitting species distribution models. Higher model performances for animals and plants than for soil microbes (fungi, bacteria and protists) suggest that the direct influence of topoclimate is stronger on above-ground species than on below-ground microorganisms. Accordingly, direct climate change effects are predicted to be stronger for above-ground than for below-ground taxa, whereas factors expressing local soil microclimate and geochemistry are likely more important to explain and forecast the occurrence patterns of soil microbiota. Detailed mapping and future scenarios of soil microclimate and microhabitats, together with comparative studies of interacting and ecologically dependent above- and below-ground biota, are thus needed to understand and realistically forecast the future distribution of ecosystems.

GNB5 Mutations Cause an Autosomal-Recessive Multisystem Syndrome with Sinus Bradycardia and Cognitive Disability.

GNB5 encodes the G protein ß subunit 5 and is involved in inhibitory G protein signaling. Here, we report mutations in GNB5 that are associated with heart-rate disturbance, eye disease, intellectual disability, gastric problems, hypotonia, and seizures in nine individuals from six families. We observed an association between the nature of the variants and clinical severity; individuals with loss-of-function alleles had more severe symptoms, including substantial developmental delay, speech defects, severe hypotonia, pathological gastro-esophageal reflux, retinal disease, and sinus-node dysfunction, whereas related heterozygotes harboring missense variants presented with a clinically milder phenotype. Zebrafish gnb5 knockouts recapitulated the phenotypic spectrum of affected individuals, including cardiac, neurological, and ophthalmological abnormalities, supporting a direct role of GNB5 in the control of heart rate, hypotonia, and vision.