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.     

Identification of two types of keratin polypeptides within the acidic cytokeratin subfamily I

Cytoskeletal filaments of the α-keratin type (cytokeratins) are a characteristic of epithelial cells. In diverse mammals (man, cow and rodents) these cytokeratins consist of a family of approximately 20 polypeptides, which may be divided into the more acidic (I) and the more basic (II) subfamilies. These two subfamilies show only limited amino acid sequence homology. In contrast, nucleic acid hybridization experiments and peptide maps have been interpreted to show that polypeptides of the same subfamily share extended sequence homology.We compare two polypeptides of the acidic cytokeratin subfamily, VIb (M_r 54,000) and VII (M_r 50,000), which are co-expressed in large amounts in bovine epidermal keratinocytes. These two epidermal keratins can be distinguished by specific antibodies and show different patterns of expression among several bovine tissues and cultured cells. In addition, they differ in the stability of their complexes with basic keratin polypeptides and in their tryptic peptide maps. The amino acid sequences deduced from the nucleotide sequences of complementary DNA clones containing the 3′ ends of the messenger RNAs for these keratins are compared with each other and with available amino acid sequences of human, murine and amphibian epidermal keratins. Bovine keratins VIb and VII share considerable sequence homology in the α-helical portion (68% residues identical) but lack significant homology in the extrahelical portion. Bovine keratin VIb shows, in its α-helical region, a pronounced sequence homology (88% identity) to the murine epidermal keratin of M_r 59,000. In addition, the non-helical carboxy-terminal regions of both proteins are glycinerich and contain a canonic sequence GGG^S_GYGG, which may be repeated several times. Moreover, their mRNAs present a highly conserved stretch of 236 nucleotides containing, in the murine sequence, the end of the coding and all of the non-coding region (81% identical nucleotides). Bovine keratin VII is considerably different from the murine M_r 59,000 keratin but is almost identical to the human cytokeratin number 14 of M_r 50,000, both in the α-helical and in the non-α-helical regions of the proteins, and the mRNAs of the human and the bovine keratins also display a high homology in their 3′ non-coding ends.The results show that in the same species keratins of the same subfamily can differ considerably, whereas equivalent keratin polypeptides of different species are readily identified by characteristic sequence homologies in the α-helical and the non-helical regions as well as in the 3′ non-coding portions of their mRNAs. Among the members of the acidic subfamily I of cytokeratin polypeptides that are co-expressed in bovine epidermis, at least two types can be distinguished by their carboxy-terminal sequences. One type is characterized by its abundance of glycine residues, a consensus GGG^S_GYGG heptapeptide sequence, which may be repeated several times, and an extended stretch of high RNA sequence homology in the 3′ non-coding part. The other type shows a predominance of serine and valine residues, a subterminal GGG^S_GYGG sequence (which has been maintained in Xenopus, cow and man) and also a high level of homology in the 3′ non-coding part of the mRNA. The data indicate that individual keratin type specificity overrides species diversity, both at the protein and the mRNA level. We discuss the evolutionary conservation and the tissue distribution of these two types of acidic keratin polypeptides as well as their possible biological functions.     

Dye-coupling visualizes networks of large-field motion-sensitive neurons in the fly

In the fly, visually guided course control is accomplished by a set of 60 large-field motion-sensitive neurons in each brain hemisphere. These neurons have been shown to receive retinotopic motion information from local motion detectors on their dendrites. In addition, recent experiments revealed extensive coupling between the large-field neurons through electrical synapses. These two processes together give rise to their broad and elaborate receptive fields significantly surpassing the extent of their dendritic fields. Here, we demonstrate that the electrical connections between different large-field neurons can be visualized using Neurobiotin dye injection into a single one of them. When combined with a fluorescent dye which does not cross electrical synapses, the injected cell can be identified unambiguously. The Neurobiotin staining corroborates the electrical coupling postulated amongst the cells of the vertical system (VS-cells) and between cells of the horizontal system (HS-cells and CH-cells). In addition, connections between some cells are revealed that have so far not been considered as electrically coupled.     

Comparison of the expression profiles induced by genotoxic and nongenotoxic carcinogens in rat liver

Application of recently developed gene expression techniques using microarrays in toxicological studies (toxicogenomics) facilitate the interpretation of a toxic compound's mode of action and may also allow the prediction of selected toxic effects based on gene expression changes. In order to test this hypothesis, we investigated whether carcinogens at doses known to induce liver tumors in the 2-year rat bioassay deregulate characteristic sets of genes in a short term in vivo study and whether these deregulated genes represent defined biological pathways. Male Wistar rats were dosed with the four nongenotoxic hepatocarcinogens methapyrilene (MPy, 60 mg/kg/day), diethylstilbestrol (DES, 10 mg/kg/day), Wy-14643 (Wy, 60 mg/kg/day), and piperonylbutoxide (PBO, 1200 mg/kg/day). After 1, 3, 7, and 14 days, the livers were taken for histopathological evaluation and for analysis of the gene expression profiles on Affymetrix RG_U34A arrays. The expression profile of the four nongenotoxic carcinogens were compared to the profiles of the four genotoxic carcinogens 2-nitrofluorene (2-NF), dimethylnitrosamine (DMN), 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK), and aflatoxin B1 (AB1) from a similar study reported previously. By using statistical and clustering tools characteristically deregulated genes were extracted and functionally classified. Distinct cellular pathways were affected by the nongenotoxic carcinogens compared to the genotoxic carcinogens which at least partly correlated with the two-stage model of carcinogenesis. Characteristic to genotoxic carcinogens were a DNA damage response and the activation of proliferative and survival signaling. Nongenotoxic carcinogens showed responses to oxidative DNA or protein damage, as well as cell cycle progression and signs of regeneration. Many of the gene alterations found with the nongenotoxic carcinogens imply compound-specific mechanisms. Although neither a single gene nor a single pathway will be sufficient to discriminate the two classes of carcinogens, it became evident that combinations of pathway-associated gene expression profiles may be used to predict a genotoxic or nongenotoxic carcinogenic potential of a compound in short-term studies.     

Some Xylaria species on termite nests

Xylaria arenicola, X. brasiliensis, X. escharoidea, X. furcata, X. nigripes, X. piperiformis and X. rhizomorpha represent ancient names of fungi known to inhabit abandoned termite nests. We attempt to redescribe them and to reduce the confusion among them. Xylaria tanganyikaensis and X. readeri, species that might be associated with termite nests, are described. We describe a new variety, X. furcata var. hirsuta, and discuss an unnamed fungus that probably represents a new species. Photographs and a key are presented to aid the identification of these taxa.     

Gene delivery with low molecular weight linear polyethylenimines

BACKGROUND: Linear polyethylenimine (LPEI) with a molecular weight (MW) of 22 kDa has been described as having a superior ability to induce gene transfer compared to its branched form. However, the transfection efficiency of the polymer cannot be enhanced beyond a certain limit due to cytotoxicity. We explored the potential of utilizing LPEIs with MWs ranging from 1.0 to 9.5 kDa to overcome this limitation. METHODS: Polyplexes of plasmid DNA encoding for the enhanced green fluorescent protein (EGFP) and various LPEIs were compared concerning their transfection efficiency and cytotoxicity in CHO-K1 and HeLa cells by flow cytometry. The involvement of endolysosomes in LPEI-mediated gene transfer was investigated by applying the proton pump inhibitor bafilomycin A1 and the lysosomotropic agent sucrose. Confocal laser scanning microscopy was applied to assess the size and shape of polyplexes under cell culture conditions, to detect their endolysosomal localization and to observe their translocation to the nucleus. RESULTS: The transfection efficiency could be altered by varying the MW and the amount of the polymer available for polyplex formation. The highest transfection efficiency (about 44%), i.e. the fraction of EGFP-positive cells, was obtained with LPEI 5.6 kDa, while the cytotoxicity remained low. The colocalization of polyplexes and endolysosomes was observed, and it appeared that the larger polyplexes escaped from the acidic organelles particularly quickly. For LPEI 5.0 and 9.0 kDa, the number of cells and nuclei that had taken up DNA after 6 hours was similar, as determined by flow cytometry. CONCLUSIONS: Our study suggests that LPEIs with low MWs are promising candidates for non-viral gene delivery, because they are more efficient and substantially less toxic than their higher MW counterparts.     

Genetic diversity and differentiation of central European freshwater pearl mussel (Margaritifera margaritifera L.) populations: implications for conservation and management

Despite the fact that mollusc species play an important role in many aquatic ecosystems, little is known about their biodiversity and conservation genetics. Freshwater pearl mussel (Margaritifera margaritifera L.) populations are seriously declining all over Europe and a variety of conservation programs are being established to support the remaining endangered central European populations. In order to provide guidelines for conservation strategies and management programs, we investigated the genetic structure of 24 freshwater pearl mussel populations originating from five major central European drainages including Elbe, Danube, Rhine, Maas and Weser, representing the last and most important populations in this area. We present a nondestructive sampling method of haemolymph for DNA analyses, which is applicable for endangered bivalves. The analyses of nine microsatellite loci with different levels of polymorphism revealed a high degree of fragmented population structure and very different levels of genetic diversity within populations. These patterns can be explained by historical and demographic effects and have been enforced by anthropogenic activities. Even within drainages, distinct conservation units were detected, as revealed from high F(ST) values, private alleles and genetic distance measures. Populations sampled close to contact zones between main drainage systems showed lowest levels of correct assignment to present-day drainage systems. Populations with high priority for conservation should not only be selected by means of census population size and geographical distance to other populations. Instead, detailed genetic analyses are mandatory for revealing differentiation and diversity parameters, which should be combined with ecological criteria for sustainable conservation and recovery programs.     

Reproductive conflict in social insects: male production by workers in a slave-making ant

In insect societies, workers cooperate but may also pursue their individual interests, such as laying viable male eggs. The case of obligatory slave-making ants is of particular interest because workers do not engage in maintenance activities and foraging. Therefore, worker egg laying is expected to be less detrimental for colony efficiency than in related, nonparasitic species. Furthermore, as slave-making workers usually do not perform brood care and thus might have little power in manipulating sex allocation, they might be more strongly selected to increase their direct fitness by producing their own sons than workers in nonparasitic species. In this study we investigated worker reproduction in four natural colonies of the slave-making ant Polyergus rufescens, using highly variable microsatellite markers. Our results show that workers produce up to 100% of the males. This study thus presents the first direct evidence of an almost complete takeover of male reproduction by workers in ants.     

The identification of 1,3-oxazolidine-2-thiones and 1,3-thiazolidine-2-thiones from the reaction of glucose with benzyl isothiocyanate

The structure of interaction products resulting from the reaction of unmodified glucose with benzyl isothiocyanate is reported. Prior to their identification, the main products of this reaction were isolated using solid-phase extraction (SPE) as well as preparative HPLC. They were then identified by NMR and MS as 3-benzyl-4-hydroxy-5-(D-arabino-1,2,3,4-tetrahydroxybutyl)-1,3-oxazolidine-2-thione, 3-benzyl-4-hydroxy-4-hydroxymethyl-5-(D-erythro-1,2,3-trihydroxypropyl)-1,3-oxazolidine-2-thione, N-benzyl-(D-gluco-4,5-dihydroxy-6-hydroxymethyl-tetrahydropyrano)[2,3-b]oxazolidine-2-thione and 3-benzyl-4-(N-benzyl amino)-5-(D-arabino-1,2,3,4-tetrahydroxybutyl)-1,3-thiazolidine-2-thione. The identity of the last compound was secured by X-ray crystal structure data.     

Mutational and structural studies of the diisopropylfluorophosphatase from Loligo vulgaris shed new light on the catalytic mechanism of the enzyme

The active site, the substrate binding site, and the metal binding sites of the diisopropylfluorophosphatase (DFPase) from Loligo vulgaris have been modified by means of site-directed mutagenesis to improve our understanding of the reaction mechanism. Enzymatic characterization of mutants located in the major groove of the substrate binding pocket indicates that large hydrophobic side chains at these positions are favorable for substrate turnover. Moreover, the active site residue His287 proved to be beneficial, but not essential, for DFP hydrolysis. In most cases, hydrophobic side chains at position 287 led to significant catalytic activities although reduced relative to the wild-type enzyme. With respect to the Ca-1 binding site, where catalysis occurs, various mutants indicated that the net charge at this calcium-binding site as well as the relative positions of the charged calcium ligands is crucial for catalytic activity. The importance of the electrostatic potential at the active site was furthermore revealed by various mutations of residues lining the interior of the central water-filled tunnel, which traverses the entire protein structure. In this respect, the structural features of residue His181, which is located at the opposite end of the DFPase tunnel relative to the active site, were characterized extensively. It was concluded that a tunnel-spanning hydrogen bond network, which includes a large number of apparently slow exchanging water molecules, relays any modifications in the electrostatics of the system to the active site, thus affecting the catalytic reactivity of the enzyme.