Cerebral organoids display dynamic clonal growth and tunable tissue replenishment.

During brain development, neural progenitors expand through symmetric divisions before giving rise to differentiating cell types via asymmetric divisions. Transition between those modes varies among individual neural stem cells, resulting in clones of different sizes. Imaging-based lineage tracing allows for lineage analysis at high cellular resolution but systematic approaches to analyse clonal behaviour of entire tissues are currently lacking. Here we implement whole-tissue lineage tracing by genomic DNA barcoding in 3D human cerebral organoids, to show that individual stem cell clones produce progeny on a vastly variable scale. By using stochastic modelling we find that variable lineage sizes arise because a subpopulation of lineages retains symmetrically dividing cells. We show that lineage sizes can adjust to tissue demands after growth perturbation via chemical ablation or genetic restriction of a subset of cells in chimeric organoids. Our data suggest that adaptive plasticity of stem cell populations ensures robustness of development in human brain organoids.

Neutral competition explains the clonal composition of neural organoids.

Neural organoids model the development of the human brain and are an indispensable tool for studying neurodevelopment. Whole-organoid lineage tracing has revealed the number of progenies arising from each initial stem cell to be highly diverse, with lineage sizes ranging from one to more than 20,000 cells. This high variability exceeds what can be explained by existing stochastic models of corticogenesis and indicates the existence of an additional source of stochasticity. To explain this variability, we introduce the SAN model which distinguishes Symmetrically diving, Asymmetrically dividing, and Non-proliferating cells. In the SAN model, the additional source of stochasticity is the survival time of a lineage's pool of symmetrically dividing cells. These survival times result from neutral competition within the sub-population of all symmetrically dividing cells. We demonstrate that our model explains the experimentally observed variability of lineage sizes and derive the quantitative relationship between survival time and lineage size. We also show that our model implies the existence of a regulatory mechanism which keeps the size of the symmetrically dividing cell population constant. Our results provide quantitative insight into the clonal composition of neural organoids and how it arises. This is relevant for many applications of neural organoids, and similar processes may occur in other developing tissues both in vitro and in vivo.

Learning From an Artificial Neural Network in Phylogenetics.

We show that an iterative ansatz of deep learning and human intelligence guided simplification may lead to surprisingly simple solutions for a difficult problem in phylogenetics. Distinguishing Farris and Felsenstein trees is a longstanding problem in phylogenetic tree reconstruction. The Artificial Neural Network F-zoneNN solves this problem for 4-taxon alignments evolved under the Jukes-Cantor model. It distinguishes between Farris and Felsenstein trees, but owing to its complexity, lacks transparency in its mechanism of discernment. Based on the simplification of F-zoneNN and alignment properties we constructed the function FarFelDiscerner. In contrast to F-zoneNN, FarFelDiscerner's decision process is understandable. Moreover, FarFelDiscerner is significantly simpler than F-zoneNN. Despite its simplicity this function infers the tree-type almost perfectly on noise-free data, and also performs well on simulated noisy alignments of finite length. We applied FarFelDiscerner to the historical Holometabola alignments where it places Strepsiptera with beetles, concordant with the current scientific view.

Critical Growth of Cerebral Tissue in Organoids: Theory and Experiments.

We develop a Fokker-Planck theory of tissue growth with three types of cells (symmetrically dividing, asymmetrically dividing, and nondividing) as main agents to study the growth dynamics of human cerebral organoids. Fitting the theory to lineage tracing data obtained in next generation sequencing experiments, we show that the growth of cerebral organoids is a critical process. We derive analytical expressions describing the time evolution of clonal lineage sizes and show how power-law distributions arise in the limit of long times due to the vanishing of a characteristic growth scale. We discuss that the independence of critical growth on initial conditions could be biologically advantageous.

ModelRevelator: Fast phylogenetic model estimation via deep learning.

Selecting the best model of sequence evolution for a multiple-sequence-alignment (MSA) constitutes the first step of phylogenetic tree reconstruction. Common approaches for inferring nucleotide models typically apply maximum likelihood (ML) methods, with discrimination between models determined by one of several information criteria. This requires tree reconstruction and optimisation which can be computationally expensive. We demonstrate that neural networks can be used to perform model selection, without the need to reconstruct trees, optimise parameters, or calculate likelihoods. We introduce ModelRevelator, a model selection tool underpinned by two deep neural networks. The first neural network, NNmodelfind, recommends one of six commonly used models of sequence evolution, ranging in complexity from Jukes and Cantor to General Time Reversible. The second, NNalphafind, recommends whether or not a Γ-distributed rate heterogeneous model should be incorporated, and if so, provides an estimate of the shape parameter, ɑ. Users can simply input an MSA into ModelRevelator, and swiftly receive output recommending the evolutionary model, inclusive of the presence or absence of rate heterogeneity, and an estimate of ɑ. We show that ModelRevelator performs comparably with likelihood-based methods and the recently published machine learning method ModelTeller over a wide range of parameter settings, with significant potential savings in computational effort. Further, we show that this performance is not restricted to the alignments on which the networks were trained, but is maintained even on unseen empirical data. We expect that ModelRevelator will provide a valuable alternative for phylogeneticists, especially where traditional methods of model selection are computationally prohibitive.

SVhound: detection of regions that harbor yet undetected structural variation.

BACKGROUND: Recent population studies are ever growing in number of samples to investigate the diversity of a population or species. These studies reveal new polymorphism that lead to important insights into the mechanisms of evolution, but are also important for the interpretation of these variations. Nevertheless, while the full catalog of variations across entire species remains unknown, we can predict which regions harbor additional not yet detected variations and investigate their properties, thereby enhancing the analysis for potentially missed variants. RESULTS: To achieve this we developed SVhound ( https://github.com/lfpaulin/SVhound ), which based on a population level SVs dataset can predict regions that harbor unseen SV alleles. We tested SVhound using subsets of the 1000 genomes project data and showed that its correlation (average correlation of 2800 tests r = 0.7136) is high to the full data set. Next, we utilized SVhound to investigate potentially missed or understudied regions across 1KGP and CCDG. Lastly we also apply SVhound on a small and novel SV call set for rhesus macaque (Macaca mulatta) and discuss the impact and choice of parameters for SVhound. CONCLUSIONS: SVhound is a unique method to identify potential regions that harbor hidden diversity in model and non model organisms and can also be potentially used to ensure high quality of SV call sets.

Molecular archaeology of human cognitive traits.

The brains and minds of our human ancestors remain inaccessible for experimental exploration. Therefore, we reconstructed human cognitive evolution by projecting nonsynonymous/synonymous rate ratios (ω values) in mammalian phylogeny onto the anatomically modern human (AMH) brain. This atlas retraces human neurogenetic selection and allows imputation of ancestral evolution in task-related functional networks (FNs). Adaptive evolution (high ω values) is associated with excitatory neurons and synaptic function. It shifted from FNs for motor control in anthropoid ancestry (60-41 mya) to attention in ancient hominoids (26-19 mya) and hominids (19-7.4 mya). Selection in FNs for language emerged with an early hominin ancestor (7.4-1.7 mya) and was later accompanied by adaptive evolution in FNs for strategic thinking during recent (0.8 mya-present) speciation of AMHs. This pattern mirrors increasingly complex cognitive demands and suggests that co-selection for language alongside strategic thinking may have separated AMHs from their archaic Denisovan and Neanderthal relatives.

Gruffi: an algorithm for computational removal of stressed cells from brain organoid transcriptomic datasets.

Organoids enable in vitro modeling of complex developmental processes and disease pathologies. Like most 3D cultures, organoids lack sufficient oxygen supply and therefore experience cellular stress. These negative effects are particularly prominent in complex models, such as brain organoids, and can affect lineage commitment. Here, we analyze brain organoid and fetal single-cell RNA sequencing (scRNAseq) data from published and new datasets, totaling about 190,000 cells. We identify a unique stress signature in the data from all organoid samples, but not in fetal samples. We demonstrate that cell stress is limited to a defined subpopulation of cells that is unique to organoids and does not affect neuronal specification or maturation. We have developed a computational algorithm, Gruffi, which uses granular functional filtering to identify and remove stressed cells from any organoid scRNAseq dataset in an unbiased manner. We validated our method using six additional datasets from different organoid protocols and early brains, and show its usefulness to other organoid systems including retinal organoids. Our data show that the adverse effects of cell stress can be corrected by bioinformatic analysis for improved delineation of developmental trajectories and resemblance to in vivo data.

NMD is required for timely cell fate transitions by fine-tuning gene expression and regulating translation.

Cell fate transitions depend on balanced rewiring of transcription and translation programs to mediate ordered developmental progression. Components of the nonsense-mediated mRNA decay (NMD) pathway have been implicated in regulating embryonic stem cell (ESC) differentiation, but the exact mechanism is unclear. Here we show that NMD controls expression levels of the translation initiation factor Eif4a2 and its premature termination codon-encoding isoform (Eif4a2PTC ). NMD deficiency leads to translation of the truncated eIF4A2PTC protein. eIF4A2PTC elicits increased mTORC1 activity and translation rates and causes differentiation delays. This establishes a previously unknown feedback loop between NMD and translation initiation. Furthermore, our results show a clear hierarchy in the severity of target deregulation and differentiation phenotypes between NMD effector KOs (Smg5 KO > Smg6 KO > Smg7 KO), which highlights heterodimer-independent functions for SMG5 and SMG7. Together, our findings expose an intricate link between mRNA homeostasis and mTORC1 activity that must be maintained for normal dynamics of cell state transitions.

It Is Just a Matter of Time: Balancing Homologous Recombination and Non-homologous End Joining at the rDNA Locus During Meiosis.

Ribosomal RNA genes (rDNAs) are located in large domains of hundreds of rDNA units organized in a head-to-tail manner. The proper and stable inheritance of rDNA clusters is of paramount importance for survival. Yet, these highly repetitive elements pose a potential risk to the genome since they can undergo non-allelic exchanges. Here, we review the current knowledge of the organization of the rDNA clusters in Arabidopsis thaliana and their stability during meiosis. Recent findings suggest that during meiosis, all rDNA loci are embedded within the nucleolus favoring non-homologous end joining (NHEJ) as a repair mechanism, while DNA repair via homologous recombination (HR) appears to be a rare event. We propose a model where (1) frequent meiotic NHEJ events generate abundant single nucleotide polymorphisms and insertions/deletions within the rDNA, resulting in a heterogeneous population of rDNA units and (2) rare HR events dynamically change rDNA unit numbers, only to be observed in large populations over many generations. Based on the latest efforts to delineate the entire rDNA sequence in A. thaliana, we discuss evidence supporting this model. The results compiled so far draw a surprising picture of rDNA sequence heterogeneity between individual units. Furthermore, rDNA cluster sizes have been recognized as relatively stable when observing less than 10 generations, yet emerged as major determinant of genome size variation between different A. thaliana ecotypes. The sequencing efforts also revealed that transcripts from the diverse rDNA units yield heterogenous ribosome populations with potential functional implications. These findings strongly motivate further research to understand the mechanisms that maintain the metastable state of rDNA loci.

Predictive Antibiotic Susceptibility Testing by Next-Generation Sequencing for Periprosthetic Joint Infections: Potential and Limitations.

Joint replacement surgeries are one of the most frequent medical interventions globally. Infections of prosthetic joints are a major health challenge and typically require prolonged or even indefinite antibiotic treatment. As multidrug-resistant pathogens continue to rise globally, novel diagnostics are critical to ensure appropriate treatment and help with prosthetic joint infections (PJI) management. To this end, recent studies have shown the potential of molecular methods such as next-generation sequencing to complement established phenotypic, culture-based methods. Together with advanced bioinformatics approaches, next-generation sequencing can provide comprehensive information on pathogen identity as well as antimicrobial susceptibility, potentially enabling rapid diagnosis and targeted therapy of PJIs. In this review, we summarize current developments in next generation sequencing based predictive antibiotic susceptibility testing and discuss potential and limitations for common PJI pathogens.

Culture-Free Detection of Antibiotic Resistance Markers from Native Patient Samples by Hybridization Capture Sequencing.

The increasing incidence of antimicrobial resistance (AMR) is a major global challenge. Routine techniques for molecular AMR marker detection are largely based on low-plex PCR and detect dozens to hundreds of AMR markers. To allow for comprehensive and sensitive profiling of AMR markers, we developed a capture-based next generation sequencing (NGS) workflow featuring a novel AMR marker panel based on the curated AMR database ARESdb. Our primary objective was to compare the sensitivity of target enrichment-based AMR marker detection to metagenomics sequencing. Therefore, we determined the limit of detection (LOD) in synovial fluid and urine samples across four key pathogens. We further demonstrated proof-of-concept for AMR marker profiling from septic samples using a selection of urine samples with confirmed monoinfection. The results showed that the capture-based workflow is more sensitive and requires lower sequencing depth compared with metagenomics sequencing, allowing for comprehensive AMR marker detection with an LOD of 1000 CFU/mL. Combining the ARESdb AMR panel with 16S rRNA gene sequencing allowed for the culture-free detection of bacterial taxa and AMR markers directly from septic patient samples at an average sensitivity of 99%. Summarizing, the newly developed ARESdb AMR panel may serve as a valuable tool for comprehensive and sensitive AMR marker detection.

Sequencing of the Arabidopsis NOR2 reveals its distinct organization and tissue-specific rRNA ribosomal variants.

Despite vast differences between organisms, some characteristics of their genomes are conserved, such as the nucleolus organizing region (NOR). The NOR is constituted of multiple, highly repetitive rDNA genes, encoding the catalytic ribosomal core RNAs which are transcribed from 45S rDNA units. Their precise sequence information and organization remain uncharacterized. Here, using a combination of long- and short-read sequencing technologies we assemble contigs of the Arabidopsis NOR2 rDNA domain. We identify several expressed rRNA gene variants which are integrated into translating ribosomes in a tissue-specific manner. These findings support the concept of tissue specific ribosome subpopulations that differ in their rRNA composition and provide insights into the higher order organization of NOR2.

TMT-Opsins differentially modulate medaka brain function in a context-dependent manner.

Vertebrate behavior is strongly influenced by light. Light receptors, encoded by functional opsin proteins, are present inside the vertebrate brain and peripheral tissues. This expression feature is present from fishes to human and appears to be particularly prominent in diurnal vertebrates. Despite their conserved widespread occurrence, the nonvisual functions of opsins are still largely enigmatic. This is even more apparent when considering the high number of opsins. Teleosts possess around 40 opsin genes, present from young developmental stages to adulthood. Many of these opsins have been shown to function as light receptors. This raises the question of whether this large number might mainly reflect functional redundancy or rather maximally enables teleosts to optimally use the complex light information present under water. We focus on tmt-opsin1b and tmt-opsin2, c-opsins with ancestral-type sequence features, conserved across several vertebrate phyla, expressed with partly similar expression in non-rod, non-cone, non-retinal-ganglion-cell brain tissues and with a similar spectral sensitivity. The characterization of the single mutants revealed age- and light-dependent behavioral changes, as well as an impact on the levels of the preprohormone sst1b and the voltage-gated sodium channel subunit scn12aa. The amount of daytime rest is affected independently of the eyes, pineal organ, and circadian clock in tmt-opsin1b mutants. We further focused on daytime behavior and the molecular changes in tmt-opsin1b/2 double mutants, and found that-despite their similar expression and spectral features-these opsins interact in part nonadditively. Specifically, double mutants complement molecular and behavioral phenotypes observed in single mutants in a partly age-dependent fashion. Our work provides a starting point to disentangle the highly complex interactions of vertebrate nonvisual opsins, suggesting that tmt-opsin-expressing cells together with other visual and nonvisual opsins provide detailed light information to the organism for behavioral fine-tuning. This work also provides a stepping stone to unravel how vertebrate species with conserved opsins, but living in different ecological niches, respond to similar light cues and how human-generated artificial light might impact on behavioral processes in natural environments.

Timing strains of the marine insect Clunio marinus diverged and persist with gene flow.

Genetic divergence of populations in the presence of gene flow is a central theme in speciation research. Theory predicts that divergence can happen with full range overlap - in sympatry - driven by ecological factors, but there are few empirical examples of how ecologically divergent selection can overcome gene flow and lead to reproductive isolation. In the marine midge Clunio marinus (Diptera: Chironomidae) reproduction is ecologically restricted to the time of the lowest tides, which is ensured through accurate control of development and adult emergence by circalunar and circadian clocks. As tidal regimes differ along the coastline, locally adapted timing strains of C. marinus are found in different sites across Europe. At the same time, ecologically suitable low tides occur at both full and new moon and twice a day, providing C. marinus with four nonoverlapping temporal niches at every geographic location. Along the coast of Brittany, which is characterized by a steep gradient in timing of the tides, we found an unusually large number of differentially adapted timing strains, and the first known instances of sympatric C. marinus strains occupying divergent temporal niches. Analysis of mitochondrial genotypes suggests that these timing strains originated from a single recent colonization event. Nuclear genotypes show strong gene flow, sympatric timing strains being the least differentiated. Even when sympatric strains exist in nonoverlapping temporal niches, timing adaptations do not result in genome-wide genetic divergence, suggesting timing adaptations are maintained by permanent ecological selection. This constitutes a model case for incipient ecological divergence with gene flow.

Complex Evolution of Light-Dependent Protochlorophyllide Oxidoreductases in Aerobic Anoxygenic Phototrophs: Origin, Phylogeny, and Function.

Light-dependent protochlorophyllide oxidoreductase (LPOR) and dark-operative protochlorophyllide oxidoreductase are evolutionary and structurally distinct enzymes that are essential for the synthesis of (bacterio)chlorophyll, the primary pigment needed for both anoxygenic and oxygenic photosynthesis. In contrast to the long-held hypothesis that LPORs are only present in oxygenic phototrophs, we recently identified a functional LPOR in the aerobic anoxygenic phototrophic bacterium (AAPB) Dinoroseobacter shibae and attributed its presence to a single horizontal gene transfer event from cyanobacteria. Here, we provide evidence for the more widespread presence of genuine LPOR enzymes in AAPBs. An exhaustive bioinformatics search identified 36 putative LPORs outside of oxygenic phototrophic bacteria (cyanobacteria) with the majority being AAPBs. Using in vitro and in vivo assays, we show that the large majority of the tested AAPB enzymes are genuine LPORs. Solution structural analyses, performed for two of the AAPB LPORs, revealed a globally conserved structure when compared with a well-characterized cyanobacterial LPOR. Phylogenetic analyses suggest that LPORs were transferred not only from cyanobacteria but also subsequently between proteobacteria and from proteobacteria to Gemmatimonadetes. Our study thus provides another interesting example for the complex evolutionary processes that govern the evolution of bacteria, involving multiple horizontal gene transfer events that likely occurred at different time points and involved different donors.

A human tissue screen identifies a regulator of ER secretion as a brain-size determinant.

Loss-of-function (LOF) screens provide a powerful approach to identify regulators in biological processes. Pioneered in laboratory animals, LOF screens of human genes are currently restricted to two-dimensional cell cultures, which hinders the testing of gene functions requiring tissue context. Here, we present CRISPR-lineage tracing at cellular resolution in heterogeneous tissue (CRISPR-LICHT), which enables parallel LOF studies in human cerebral organoid tissue. We used CRISPR-LICHT to test 173 microcephaly candidate genes, revealing 25 to be involved in known and uncharacterized microcephaly-associated pathways. We characterized IER3IP1, which regulates the endoplasmic reticulum (ER) function and extracellular matrix protein secretion crucial for tissue integrity, the dysregulation of which results in microcephaly. Our human tissue screening technology identifies microcephaly genes and mechanisms involved in brain-size control.

Poly(ADP-ribose) glycohydrolase coordinates meiotic DNA double-strand break induction and repair independent of its catalytic activity.

Poly(ADP-ribosyl)ation is a reversible post-translational modification synthetized by ADP-ribose transferases and removed by poly(ADP-ribose) glycohydrolase (PARG), which plays important roles in DNA damage repair. While well-studied in somatic tissues, much less is known about poly(ADP-ribosyl)ation in the germline, where DNA double-strand breaks are introduced by a regulated program and repaired by crossover recombination to establish a tether between homologous chromosomes. The interaction between the parental chromosomes is facilitated by meiotic specific adaptation of the chromosome axes and cohesins, and reinforced by the synaptonemal complex. Here, we uncover an unexpected role for PARG in coordinating the induction of meiotic DNA breaks and their homologous recombination-mediated repair in Caenorhabditis elegans. PARG-1/PARG interacts with both axial and central elements of the synaptonemal complex, REC-8/Rec8 and the MRN/X complex. PARG-1 shapes the recombination landscape and reinforces the tightly regulated control of crossover numbers without requiring its catalytic activity. We unravel roles in regulating meiosis, beyond its enzymatic activity in poly(ADP-ribose) catabolism.

Structure of the space of taboo-free sequences.

Models of sequence evolution typically assume that all sequences are possible. However, restriction enzymes that cut DNA at specific recognition sites provide an example where carrying a recognition site can be lethal. Motivated by this observation, we studied the set of strings over a finite alphabet with taboos, that is, with prohibited substrings. The taboo-set is referred to as [Formula: see text] and any allowed string as a taboo-free string. We consider the so-called Hamming graph [Formula: see text], whose vertices are taboo-free strings of length n and whose edges connect two taboo-free strings if their Hamming distance equals one. Any (random) walk on this graph describes the evolution of a DNA sequence that avoids taboos. We describe the construction of the vertex set of [Formula: see text]. Then we state conditions under which [Formula: see text] and its suffix subgraphs are connected. Moreover, we provide an algorithm that determines if all these graphs are connected for an arbitrary [Formula: see text]. As an application of the algorithm, we show that about [Formula: see text] of bacteria listed in REBASE have a taboo-set that induces connected taboo-free Hamming graphs, because they have less than four type II restriction enzymes. On the other hand, four properly chosen taboos are enough to disconnect one suffix subgraph, and consequently connectivity of taboo-free Hamming graphs could change depending on the composition of restriction sites.

A compendium of genome-wide sequence reads from NBS (nucleotide binding site) domains of resistance genes in the common potato.

SolariX is a compendium of DNA sequence tags from the nucleotide binding site (NBS) domain of disease resistance genes of the common potato, Solanum tuberosum Group Tuberosum. The sequences, which we call NBS tags, for nearly all NBS domains from 91 genomes-representing a wide range of historical and contemporary potato cultivars, 24 breeding programs and 200 years-were generated using just 16 amplification primers and high-throughput sequencing. The NBS tags were mapped to 587 NBS domains on the draft potato genome DM, where we detected an average, over all the samples, of 26 nucleotide polymorphisms on each locus. The total number of NBS domains observed, differed between potato cultivars. However, both modern and old cultivars possessed comparable levels of variability, and neither the individual breeder or country nor the generation or time appeared to correlate with the NBS domain frequencies. Our attempts to detect haplotypes (i.e., sets of linked nucleotide polymorphisms) frequently yielded more than the possible 4 alleles per domain indicating potential locus intermixing during the mapping of NBS tags to the DM reference genome. Mapping inaccuracies were likely a consequence of the differences of each cultivar to the reference genome used, coupled with high levels of NBS domain sequence similarity. We illustrate that the SolariX database is useful to search for polymorphism linked with NBS-LRR R gene alleles conferring specific disease resistance and to develop molecular markers for selection.