A perspective on information optimality in a neural circuit The nematode worm Caenorhabditis elegans has a relatively simple neural system for analysis of information transmission from sensory organ to muscle fiber. Therefore, an example of a neural circuit is analyzed that originates in the nematode worm, and a method is applied for measuring its information flow efficiency by use of a model of logic gates. This model-based approach is useful where the assumptions of a logic gate design are applicable. It is also an useful approach where there are competing mathematical models for explaining the role of a neural circuit since the logic gate model can estimate the computational complexity of a network, and distinguish which of the mathematical models require fewer computations. In addition, for generalization of the concept of information optimality in biological systems, there is an extensive discussion of its role in the genetic-based pathways of organisms.
Thalamic control of sensory enhancement and sleep spindle properties in a biophysical model of thalamoreticular microcircuitry Thalamoreticular circuitry is known to play a key role in attention, cognition and the generation of sleep spindles, and is implicated in numerous brain disorders, but the cellular and synaptic mechanisms remain intractable. Therefore, we developed the first detailed computational model of mouse thalamus and thalamic reticular nucleus microcircuitry that captures morphological and biophysical properties of [~]14,000 neurons connected via [~]6M synapses, and recreates biological synaptic and gap junction connectivity. Simulations recapitulate multiple independent network-level experimental findings across different brain states, providing a novel unifying cellular and synaptic account of spontaneous and evoked activity in both wakefulness and sleep. Furthermore, we found that: 1.) inhibitory rebound produces frequency-selective enhancement of thalamic responses during wakefulness, in addition to its role in spindle generation; 2.) thalamic interactions generate the characteristic waxing and waning of spindle oscillations; and 3.) changes in thalamic excitability (e.g. due to neuromodulation) control spindle frequency and occurrence. The model is openly available and provides a new tool to interpret spindle oscillations and test hypotheses of thalamoreticular circuit function and dysfunction across different network states in health and disease.
How do stomata evolve on the road to C4 photosynthesis? C4 photosynthesis optimizes plants carbon and water relations, allowing high photosynthetic rate with low stomatal conductance. Stomata have long been believed as a part of C4 syndrome. However, it remains unclear how stomata traits evolved along the path from C3 to C4. Stomatal patterning was examined in Flaveria genus, a model for studying C4 evolution. Comparative, transgenic and semi-in-vitro experiments were used to study molecular basis that underlies stomatal traits along C4 evolution. Novel results: the evolution from C3 to C4 species through intermediate species is accompanied by a stepwise rather than an abrupt change in the stomatal traits. The initial change occurs near Type II and dramatic change occurs at the C4-like species. On the road to C4, stomata become less in number but bigger in size and changes in stomatal density dominates changes in maximum stomatal conductance (gsmax). The reduction of FSTOMAGEN expression underlies altered gsmax between Flaveria species with different photosynthetic pathways and likely occurs in other C4 lineages. Our study provides insight into the pattern, mechanism and role of stomatal evolution along the road towards C4. This work highlights the stomatal traits in the current C4 evolutionary model and the co-evolution of photosynthetic pathway and stomata. The contributions from the progenitor genomes of the mesopolyploid Brassiceae are evolutionarily distinct but functionally compatible The members of the tribe Brassiceae share a whole genome triplication (WGT), and one proposed model for its formation is a "two-step" pair of hybridizations producing hexaploid descendants. However, evidence for this model is incomplete, and the evolutionary and functional constraints that drove evolution after the hexaploidy are even less understood. Here we report a new genome sequence of Crambe hispanica, a species sister to most sequenced Brassiceae. Using this new genome and three others that share the hexaploidy, we traced the history of gene loss after the WGT using POInT (the Polyploidy Orthology Inference Tool). We confirm the two-step formation model and infer that there was a significant temporal gap between those two allopolyploidizations, with about a third of the gene losses from the first two subgenomes occurring prior to the arrival of the third. We also, for the 90,000 individual genes in our study, make parental "subgenome" assignments, inferring, with measured uncertainty, which of the progenitor genomes of the allohexaploidy each gene derives from. We further show that each subgenome has a statistically distinguishable rate of homoeolog losses. There is little indication of functional distinction between the three subgenomes: the individual subgenomes show no patterns of functional enrichment, no excess of shared protein-protein or metabolic interactions between their members, and no biases in their likelihood of having experienced a recent selective sweep. We propose a "mix and match" model of allopolyploidy, where subgenome origin drives homoeolog loss propensities but where genes from different subgenomes function together without difficulty. Penalized Logistic Regression Analysis for Genetic Association Studies of Binary Phenotypes IntroductionIncreasingly, logistic regression methods for genetic association studies of binary phenotypes must be able to accommodate data sparsity, which arises from unbalanced case-control ratios and/or rare genetic variants. Sparseness leads to maximum likelihood estimators (MLEs) of log-OR parameters that are biased away from their null value of zero and tests with inflated type 1 errors. Different penalized-likelihood methods have been developed to mitigate sparse-data bias. We study penalized logistic regression using a class of log-F priors indexed by a shrinkage parameter m to shrink the biased MLE towards zero.

MethodsWe propose a two-step approach to the analysis of a genetic association study: first, a set of variants that show evidence of association with the trait is used to estimate m; and second, the estimated m is used for log-F -penalized logistic regression analyses of all variants using data augmentation with standard software. Our estimate of m is the maximizer of a marginal likelihood obtained by integrating the latent log-ORs out of the joint distribution of the parameters and observed data. We consider two approximate approaches to maximizing the marginal likelihood: (i) a Monte Carlo EM algorithm (MCEM) and (ii) a Laplace approximation (LA) to each integral, followed by derivative-free optimization of the approximation.

ResultsWe evaluate the statistical properties of our proposed two-step method and compared its performance to other shrinkage methods by a simulation study. Our simulation studies suggest that the proposed log-F -penalized approach has lower bias and mean squared error than other methods considered. We also illustrate the approach on data from a study of genetic associations with "super senior" cases and middle aged controls.

Discussion/ConclusionWe have proposed a method for single rare variant analysis with binary phenotypes by logistic regression penalized by log-F priors. Our method has the advantage of being easily extended to correct for confounding due to population structure and genetic relatedness through a data augmentation approach. Ecological interactions in breast cancer: Cell facilitation promotes growth and survival under drug pressure The interplay of positive and negative interactions between drug sensitive and resistant cells determines the efficacy of treatment in heterogeneous cancer cell populations. We developed isogenic estrogen receptor positive (ER+) breast cancer cell lineages sensitive and resistant to ribociclib-induced CDK4/6 inhibition to investigate cell interactions under selective pressure. During treatment with ribociclib alone or in combination with endocrine therapies, sensitive cells survive and proliferate better when grown together with resistant cells than when grown in monoculture. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) assays detected estrone, a less active metabolite, in medium taken from each lineage but detected highly active estradiol only from resistant cells. Because sensitive cells are more dependent on estrogen signaling than resistant cells, they benefit from the increase in local estradiol from resistant cells, leading to facilitated growth when grown under treatment in coculture. We developed a suite of mathematical models to quantify the strength of competition and facilitation between cell lines and evaluate the effects of treatments that modify these interactions. These analyses support the conclusion that resistant cell lineages facilitate growth of sensitive lineages during therapy through the increased production and sharing of local estradiol and themselves suffer from increased competition and suppression as a result. These models show that blocking facilitation has the potential to control both resistant and sensitive cell populations and inhibit the emergence of a refractory population. This combination of experiments and models shows that interactions within heterogeneous cancer cell populations can combine positive and negative interactions to shape the dynamics of resistance. Decellularized Normal and Tumor Extracellular Matrix as Scaffold for Cancer Organoid Cultures of Colorectal Peritoneal Metastases Peritoneal metastases (PM) from colorectal cancer (CRC) are associated with poor survival. The extracellular matrix (ECM) plays a fundamental role in modulating the homing of CRC metastases to the peritoneum. The mechanisms underlying the interactions between metastatic cells and the ECM, however, remain poorly understood and the number of in vitro models available for the study of the peritoneal metastatic process is limited. Here, we show that decellularized ECM of the peritoneal cavity allows the growth of organoids obtained from PM, favoring the development of three-dimensional nodules that maintain the characteristics of in vivo PM. Organoids preferentially grow on scaffolds obtained from neoplastic peritoneum, which are characterized by greater stiffness than normal scaffolds. A gene expression analysis of organoids grown on different substrates reflected faithfully the clinical and biological characteristics of the organoids. An impact of the ECM on the response to standard chemotherapy treatment for PM was also observed.

SignificanceEvidence of the value of ex vivo 3D models obtained by combining patient-derived extracellular matrices depleted of cellular components and organoids to mimic the metastatic niche, to be used as a tool to develop new therapeutic strategies in a biologically relevant context, to personalize treatments and increase their efficacy. CryoEM of endogenous mammalian V-ATPase interacting with the TLDc protein mEAK-7 V-ATPases are rotary proton pumps that serve as signaling hubs with numerous proposed binding partners in cells. We used cryoEM combined with exhaustive focused classification to detect endogenous proteins that associate with V-ATPase from porcine kidney. A super-stoichiometric copy of subunit C was found in ~3% of complexes, while an additional ~1.6% of complexes bound mEAK-7, a protein with proposed roles in dauer formation in nematodes and mTOR signaling in mammals. High-resolution cryoEM of porcine kidney V-ATPase with recombinant mEAK-7 shows that mEAK-7s TLDc domain interacts with V-ATPases stator while its C-terminal  helix binds V-ATPases rotor. This crosslink would be expected to inhibit rotary catalysis. However, unlike inhibition of yeast V-ATPase by the TLDc protein Oxr1p, exogenous mEAK-7 does not inhibit V-ATPase and mEAK-7 overexpression in cells does not alter lysosomal or phagosomal pH. Instead, cryoEM suggests that interaction of mEAK-7 with V-ATPase is disrupted by ATP-induced rotation of the rotor. Comparison of Oxr1p and mEAK-7 binding explains this difference. Together, these results show that differences in V-ATPase binding by TLDc domain-containing proteins can lead to effects ranging from strong inhibition to formation of labile interactions that are sensitive to the enzymes activity. Characterising ontogeny of numerosity discrimination in zebrafish reveals the use of multiple, numerical and non- numerical mechanisms. A sense of non-symbolic numerical magnitudes is widespread in the animal kingdom and has been documented in adult zebrafish. Here we investigated the ontogeny of this ability using a group size preference task in juvenile zebrafish. Fish showed group size preference from 21 days post fertilization (dpf) and reliably chose the larger group when presented with discriminations of between 1 vs. 3, 2 vs. 5 and 2 vs. 3 conspecifics but not 2 vs. 4 conspecifics. When the ratio between the number of conspecifics in each group was maintained at 1:2, fish could discriminate between 1 vs. 2 individuals and 3 vs. 6, but again, not when given a choice between 2 vs. 4 individuals. These findings are in agreement with studies in other species suggesting the systems involved in quantity representation do not operate separately from other cognitive mechanisms. Rather they suggest quantity processing in fish may be the result of an interplay between attentional, cognitive and memory-related mechanisms as in humans and other animals. Our results emphasise the potential of the use of zebrafish to explore the genetic and neural processes underlying the ontogeny and function of number cognition. Characterization of methylation profiles in spontaneous preterm birth placental villous tissue Preterm birth is a global public health crisis which results in significant neonatal and maternal mortality. Yet little is known regarding the molecular mechanisms of idiopathic spontaneous preterm birth, and we have few diagnostic markers for adequate assessment of placental development and function. Previous studies of placental pathology and our transcriptomics studies suggest a role for placental maturity in idiopathic spontaneous preterm birth. It is known that placental DNA methylation changes over gestation. We hypothesized that if placental hypermaturity is present in our samples, we would observe a unique idiopathic spontaneous preterm birth DNA methylation profile potentially driving the gene expression differences we previously identified in our placental samples. Our results indicate the idiopathic spontaneous preterm birth DNA methylation pattern mimics the term birth methylation pattern suggesting hypermaturity. Only seven significant differentially methylated regions fitting the idiopathic spontaneous preterm birth specific (relative to the controls) profile were identified, indicating unusually high similarity in DNA methylation between idiopathic spontaneous preterm birth and term birth samples. We identified an additional 1,718 significantly methylated regions in our gestational age matched controls were the idiopathic spontaneous preterm birth DNA methylation pattern mimics the term birth methylation pattern, again indicating a striking level of similarity between the idiopathic spontaneous preterm birth and term birth samples. Pathway analysis of these regions revealed differences in genes within the WNT and Cadherin signaling pathways, both of which are essential in placental development and maturation. Taken together, these data demonstrate that the idiopathic spontaneous preterm birth samples are molecularly more mature than expected given their respective gestational age which likely impacts birth timing. Detection and functional resolution of soluble multimeric immune complexes by a comprehensive FcγR reporter cell panel Fc-gamma receptor (Fc{gamma}R) activation by soluble IgG immune complexes (sICs) represents a major mechanism of inflammation in certain autoimmune diseases such as systemic lupus erythematosus (SLE). A robust and scalable test system allowing for the detection and quantification of sIC bioactivity is missing. Previously described Fc{gamma}R interaction assays are limited to certain Fc{gamma}Rs, lack scalability and flexibility, are not indicative of receptor activation or lack sensitivity towards sIC size. We developed a comprehensive reporter cell panel detecting individual activation of Fc{gamma}Rs from humans and the mouse. The reporter cell lines were integrated into an assay format that provides flexible read-outs enabling the quantification of sIC reactivity via ELISA or a fast detection using flow cytometry. This identified Fc{gamma}RIIA(H) and Fc{gamma}RIIIA as the most sIC-sensitive Fc{gamma}Rs in our test system. Applying the assay we demonstrate that sICs versus immobilized ICs are fundamentally different Fc{gamma}R-ligands with regard to Fc{gamma}R preference and signal strength. Reaching a detection limit in the very low nanomolar range, the assay proved also to be sensitive to sIC stoichiometry and size enabling for the first time a complete reproduction of the Heidelberger-Kendall precipitation curve in terms of immune receptor activation. Analyzing sera from SLE patients and mouse models of lupus and arthritis proved that sIC-dependent Fc{gamma}R activation has predictive capabilities regarding severity of SLE disease. The new methodology provides a sensitive, scalable and comprehensive tool to evaluate the size, amount and bioactivity of sICs in all settings.

One Sentence SummaryIn this study we established a comprehensive Fc{gamma}R reporter cell assay enabling the detection and quantification of soluble immune complexes generated in experimental and clinical settings. Mechanisms of A-type lamin targeting to nuclear ruptures are disrupted in LMNA- and BANF1-associated progerias Mutations in the genes LMNA and BANF1 can lead to accelerated aging syndromes called progeria. The protein products of these genes, A-type lamins and BAF, respectively, are nuclear envelope (NE) proteins that interact and participate in various cellular processes, including nuclear envelope rupture and repair. BAF localizes to sites of nuclear rupture and recruits NE-repair machinery, including the LEM-domain proteins, ESCRT-III complex, A-type lamins, and membranes. Here, we show that it is a mobile, nucleoplasmic population of A-type lamins that is rapidly recruited to ruptures in a BAF-dependent manner via BAFs association with the Ig-like {beta} fold domain of A-type lamins. These initially mobile lamins become progressively stabilized at the site of rupture. Farnesylated prelamin A and lamin B1 fail to localize to nuclear ruptures, unless that farnesylation is inhibited. Progeria-associated LMNA mutations inhibit the recruitment affected A-type lamin to nuclear ruptures, due to either permanent farnesylation or inhibition of BAF binding. A progeria-associated BAF mutant targets to nuclear ruptures but is unable to recruit A-type lamins. Together, these data reveal the mechanisms that determine how lamins respond to nuclear ruptures and how progeric mutations of LMNA and BANF1 impair recruitment of A-type lamins to nuclear ruptures. Molecular diagnosis of new isolate of tomato yellow leaf curl virus (TYLCV) in Iraq Tomato yield and quality in Iraq have been threatened by a variable range of infections caused by tomato yellow leaf curl virus. In previous studies, the TYLCV isolates were partially characterized using molecular tests for small fragments not the entire length of the virus. Sample of TYLCV-infected tomato has applied in this study to diagnose complete sequence of TYLCV isolate. Three sets of primers that belong to three well-identified strains in Iraq were used in a PCR technique and interestingly the results were negative. A new Iraqi isolate has been characterized as a first novel Iraqi isolate detected ever using next generation (NGS) and bioinformatics techniques. The NGS platform has produced about 78,232,062 paired reads of the TYLCV-infected tomato var. Oula F1. The complete raw reads of the infected variety have been analyzed using RepeatExplorer pipeline and Map to reference. The full sequence of TYLCV was reconstructed and extracted to consist of 2770 nt and then deposited in Genbank under accession number MT583814. Copy numbers and genome proportion of this sequence have been calculated that were 3523 and 0,086% respectively. Phylogenetic analyses of full nucleotide sequences confirmed close relationship to Iranian isolate (TYLCV-Kahnooj) than other published viruses. Additionally, a and {beta} DNA satellites have not discovered in the TYLCV-infected sample. Common coupling map advances GPCR-G protein selectivity Two-thirds of human hormones and one-third of clinical drugs act on membrane receptors that couple to G proteins to achieve appropriate functional responses. While G protein transducers from literature are annotated in the Guide to Pharmacology database, two recent large-scale datasets now expand the receptor-G protein  couplome. However, these three datasets differ in scope and reported G protein couplings giving different coverage and conclusions on GPCR-G protein signaling. Here, we report a common coupling map uncovering novel couplings supported by both large-scale studies, the selectivity/promiscuity of GPCRs and G proteins, and how the co-coupling and co-expression of G proteins compare to the families from phylogenetic relationships. The coupling map and insights on GPCR-G protein selectivity will catalyze advances in receptor research and cellular signaling towards the exploitation of G protein signaling bias in design of safer drugs. Fatal COVID-19 ARDS associated with incomplete AEC1 differentiation from the transitional state without senescence or fibrosis ARDS due to COVID-19 and other etiologies results from injury to the alveolar epithelial cell (AEC) barrier resulting in noncardiogenic pulmonary edema, which causes acute respiratory failure; clinical recovery requires epithelial regeneration. During physiologic regeneration in mice, AEC2s proliferate, exit the cell cycle, and transiently assume a transitional state before differentiating into AEC1s; persistence of the transitional state is associated with pulmonary fibrosis in humans. It is unknown whether transitional cells emerge and differentiate into AEC1s without fibrosis in human ARDS and why transitional cells differentiate into AEC1s during physiologic regeneration but persist in fibrosis. We hypothesized that incomplete but ongoing AEC1 differentiation from transitional cells without fibrosis may underlie persistent barrier permeability and fatal acute respiratory failure in ARDS. Immunostaining of postmortem ARDS lungs revealed abundant transitional cells in organized monolayers on alveolar septa without fibrosis. They were typically cuboidal or partially spread, sometimes flat, and occasionally expressed AEC1 markers. Immunostaining and/or interrogation of scRNAseq datasets revealed that transitional cells in mouse models of physiologic regeneration, ARDS, and fibrosis express markers of cell cycle exit but only in fibrosis express a specific senescence marker. Thus, in severe, fatal early ARDS, AEC1 differentiation from transitional cells is incomplete, underlying persistent barrier permeability and respiratory failure, but ongoing without fibrosis; senescence of transitional cells may be associated with pulmonary fibrosis. Depletion of Sun1/2 induces heterochromatin accrual in mesenchymal stem cells during adipogenesis Critical to the mechano-regulation of mesenchymal stem cells (MSC), Linker of the Nucleoskeleton and Cytoskeleton (LINC) complex transduces cytoskeletal forces to the nuclei. The LINC complex contains outer nuclear membrane Nesprin proteins that associate with the cytoskeleton and their inner nuclear membrane couplers, Sun proteins. In addition to coupling Nesprin-associated cytoskeletal elements to inner nuclear membrane, Sun proteins also function in regulating nuclear mechanics and chromatin tethering to inner nuclear membrane. This suggests that release of LINC-mediated cytoskeletal connections from cell nuclei may have different effects on chromatin organization and MSC differentiation than those due to ablation of intra-nuclear Sun proteins. To test this hypothesis we compared Sun1/2 depletion with expression of the dominant-negative KASH domain (dnKASH) that inhibits Nesprin-Sun association. In cells cultured under adipogenic conditions, disconnecting LINC from the cytoskeleton through dnKASH expression, increased adipogenic gene expression and fat droplet formation; heterochromatin H3K27 and H3K9 methylation was unaltered. In contrast, Sun1/2 depletion inhibited adipogenic gene expression and fat droplet formation; as well the anti-adipogenic effect of Sun1/2 depletion was accompanied by increased H3K9me3, which was enriched on Adipoq, silencing this fat locus. We conclude that releasing the nucleus from cytoskeletal constraints via dnKASH accelerates adipogenesis while depletion of Sun1/2 increases heterochromatin accrual on adipogenic genes in a fashion independent of LINC complex function. Therefore, while these two approaches both disable LINC functions, their divergent effects on the epigenetic landscape indicate they cannot be used interchangeably to study mechanical regulation of cell differentiation. TeamTree analysis: a new approach to evaluate scientific production Advances in science and technology depend on the work of research teams and the publication of results through peer-reviewed articles representing a growing socio-economic resource. Current methods to mine the scientific literature regarding a field of interest focus on content, but the workforce credited by authorship remains largely unexplored, and appropriate measures of scientific production are debated. Here, a new bibliometric approach named TeamTree analysis is introduced that visualizes the development and composition of the workforce driving a field. A new citation-independent measure that scales with the H index estimates impact based on publication record, genealogical ties and collaborative connections. This author-centered approach complements existing tools to mine the scientific literature and to evaluate research across disciplines. Shelterin components modulate nucleic acids condensation and phase separation in the context of telomeric DNA Telomeres are nucleoprotein complexes that protect the ends of chromosomes and are essential for chromosome stability in Eukaryotes. In cells, individual telomeres form distinct globules of finite size that appear to be smaller than expected for bare DNA. Moreover, upon changes in their protein composition, telomeres can cluster to form telomere-induced-foci (TIFs) or co-localize with promyelocytic leukemia (PML) nuclear bodies. The physical basis for collapse of individual telomeres and coalescence of multiple ones remains unclear, as does the relationship between these two phenomena. By combining single-molecule measurements, optical microscopy, turbidity assays, and simulations, we show that the telomere scaffolding protein TRF2 can condense individual DNA chains and drives coalescence of multiple DNA molecules, leading to phase separation and the formation of liquid-like droplets. Addition of the TRF2 binding protein hRap1 modulates phase boundaries and tunes the specificity of solution demixing while simultaneously altering the degree of DNA compaction. Our results suggest that the condensation of single telomeres and formation of biomolecular condensates containing multiple telomeres are two different outcomes driven by the same set of molecular interactions. Moreover, binding partners, such as other telomere components, can alter those interactions to promote single-chain DNA compaction over multiple-chain phase separation. Mapping the planet's critical natural assets Sustaining the organisms, ecosystems, and processes that underpin human well-being is necessary to achieve sustainable development. Here we identify critical natural assets, natural and semi-natural ecosystems that provide 90% of the total current magnitude of 14 types of natures contributions to people (NCP). Critical natural assets for maintaining local-scale NCP (12 of the 14 NCP mapped) comprise 30% of total global land area and 24% of national territorial waters, while 44% of land area is required for maintaining all NCP (including those that accrue at the global scale, carbon storage and moisture recycling). At least 87% of the worlds population lives in the areas benefiting from critical natural assets for local-scale NCP, while only 16% lives on the lands containing these assets. Critical natural assets also overlap substantially with areas important for biodiversity (covering area requirements for 73% of birds and 66% of mammals) and cultural diversity (representing 96% of global Indigenous and non-migrant languages). Many of the NCP mapped here are left out of international agreements focused on conserving species or mitigating climate change, yet this analysis shows that explicitly prioritizing critical natural assets for NCP could simultaneously advance development, climate, and conservation goals. Crafting policy and investment strategies that protect critical natural assets is essential for sustaining human well-being and securing Earths life support systems. GUIdEStaR (G-quadruplex, uORF, IRES, Epigenetics, Small RNA, Repeats), the integrated metadatabase in conjunction with neural network methods GUIdEStaR integrates existing databases of various types of G-quadruplex, upstream Open Reading Frame (uORF), Internal Ribosome Entry Site (IRES), methylation to RNA and histone protein, small RNA, and repeats. GUIdEStaR consists of approx. 40,000 genes and 320,000 transcripts. An mRNA transcript is divided into 5 regions (5UTR, 3UTR, exon, intron, and biological region) where each region contains presence-absence data of 169 different types of elements. Recently, artificial intelligence (AI) based analysis of sequencing data has been gaining popularity in the area of bioinformatics. GUIdEStaR generates datasets that can be used as inputs to AI methods. At the GUIdEStaR homepage, users submit gene symbols by clicking a "Send" button, and shortly result files in CSV format are available for download at the result website. Users have an option to send the result files to their email addresses. Additionally, the entire database and the example Java codes are also freely available for download. Here, we demonstrate the database usage with three neural network classification studies-1) small RNA study for classifying transcription factor (TF) genes into either one of TF mediated by small RNA originated from SARS-CoV-2 or by human microRNA (miRNA), 2) cell membrane receptor study for classifying receptor genes as either with virus interaction or without one, and 3) nonsense mediated mRNA decay (NMD) study for classifying cell membrane and nuclear receptors as either NMD target or non-target. GUIdEStaR is available for access to the easy-to-use web-based database at www.guidestar.kr and for download at https://sourceforge.net/projects/guidestar. Ubiquity and origins of structural maintenance of chromosomes (SMC) proteins in eukaryotes Structural maintenance of chromosomes (SMC) protein complexes are common in Bacteria, Archaea, and Eukaryota. SMC proteins, together with the proteins related to SMC (SMC-related proteins), constitute a superfamily of ATPases. Bacteria/Archaea and Eukaryotes are distinctive from one another in terms of the repertory of SMC proteins. A single type of SMC protein is dimerized in the bacterial and archaeal complexes, whereas eukaryotes possess six distinct SMC subfamilies (SMC1-6), constituting three heterodimeric complexes, namely cohesin, condensin, and SMC5/6 complex. Thus, to bridge the homodimeric SMC complexes in Bacteria and Archaea to the heterodimeric SMC complexes in Eukaryota, we need to invoke multiple duplications of an SMC gene followed by functional divergence. However, to our knowledge, the evolution of the SMC proteins in Eukaryota had not been examined for more than a decade. In this study, we reexamined the ubiquity of SMC1-6 in phylogenetically diverse eukaryotes that cover the major eukaryotic taxonomic groups recognized to date and provide two novel insights into the SMC evolution in eukaryotes. First, multiple secondary losses of SMC5 and SMC6 occurred in the eukaryotic evolution. Second, the SMC proteins constituting cohesin and condensin (i.e., SMC1-4), and SMC5 and SMC6 were derived from closely related but distinct ancestral proteins. Based on the above-mentioned findings, we discuss how SMC1-6 have diverged from the archaeal homologs. Cluster-specific gene markers enhance Shigella and Enteroinvasive Escherichia coli in silico serotyping Shigella and enteroinvasive Escherichia coli (EIEC) cause human bacillary dysentery with similar invasion mechanisms and share similar physiological, biochemical and genetic characteristics. The ability to differentiate Shigella and EIEC from each other is important for clinical diagnostic and epidemiologic investigations. The existing genetic signatures may not discriminate between Shigella and EIEC. However, phylogenetically, Shigella and EIEC strains are composed of multiple clusters and are different forms of E. coli. In this study, we identified 10 Shigella clusters, 7 EIEC clusters and 53 sporadic types of EIEC by examining over 17,000 publicly available Shigella/EIEC genomes. We compared Shigella and EIEC accessory genomes to identify the cluster-specific gene markers or marker sets for the 17 clusters and 53 sporadic types. The gene markers showed 99.63% accuracy and more than 97.02% specificity. In addition, we developed a freely available in silico serotyping pipeline named Shigella EIEC Cluster Enhanced Serotype Finder (ShigEiFinder) by incorporating the cluster-specific gene markers and established Shigella/EIEC serotype specific O antigen genes and modification genes into typing. ShigEiFinder can process either paired end Illumina sequencing reads or assembled genomes and almost perfectly differentiated Shigella from EIEC with 99.70% and 99.81% cluster assignment accuracy for the assembled genomes and mapped reads respectively. ShigEiFinder was able to serotype over 59 Shigella serotypes and 22 EIEC serotypes and provided a high specificity with 99.40% for assembled genomes and 99.38% for mapped reads for serotyping. The cluster markers and our new serotyping tool, ShigEiFinder (https://github.com/LanLab/ShigEiFinder), will be useful for epidemiologic and diagnostic investigations.

Impact statementThe differentiation of Shigella strains from enteroinvasive E. coli (EIEC) is important for clinical diagnosis and public health epidemiologic investigations. The similarities between Shigella and EIEC strains make this differentiation very difficult as both share common ancestries within E. coli. However, Shigella and EIEC are phylogenetically separated into multiple clusters, making high resolution separation using cluster specific genomic markers possible. In this study, we identified 17 Shigella or EIEC clusters including five that were newly identified through examination of over 17,000 publicly available Shigella and EIEC genomes. We further identified an individual or a set of cluster-specific gene markers for each cluster using comparative genomic analysis. These markers can then be used to classify isolates into clusters and were used to develop an in silico pipeline, ShigEiFinder (https://github.com/LanLab/ShigEiFinder) for accurate differentiation, cluster typing and serotyping of Shigella and EIEC from Illumina sequencing reads or assembled genomes. This study will have broad application from understanding the evolution of Shigella/EIEC to diagnosis and epidemiology.

Data summarySequencing data have been deposited at the National Center for Biotechnology Information under BioProject number PRJNA692536.

RepositoriesRaw sequence data are available from NCBI under the BioProject number PRJNA692536. First Report of Lactoria cornuta from Karnataka, West coast of India, Role of tropical cyclones in Range Extension India is known for its rich biodiversity and is fortunate to have several endemic species from different classes of vertebrates. India is home to 7.5% of the global fish diversity, with 91 endemic species of ray finned fishes, the actinopterygians. Some fish species have never been reported until recently, one such example is the long-horned Cowfish (Ostracidae), best known for inhabiting only certain areas of the Indo-Pacific oceans. However, this fish has been reported from Bay of Bengal, Costs of Goa and Lakshadweep. This is the first time this fish has been found on the Coasts of Karnataka, west coast of India. Given its morphology, the migration seems highly improbable since it can only endure lethargic swimming. In this article, we have explored the role of other external forces that could have contributed to its range extension in the east and west coast of India. In this paper we are exploring the role of Super Cyclones as a vital force in determining the expansion of its range. The fact of collection longhorn cowfish for the first time from the coast of Karwar coincides with the passage of Cyclone Ockhi on the west coast of India this suggests the possible role of super cyclone in its range extension. A year later, another specimen was recovered from a location very close to the first, which indicates the successful establishment of Lacturia cornuta in its new environment. Related morphometric and meristic studies of our specimens are consistent with previous studies reported in the Bay of Bengal. The role of Kabuki Syndrome genes KMT2D and KDM6A in development: Analysis in Human sequencing data and compared to mice and zebrafish KMT2D and KDM6A are epigenetic regulators that have been implicated in Kabuki Syndrome, a rare congenital birth defect with multiple tissue and organ abnormalities, including craniofacial and heart defects. Our previous study identified human families with mutations in the epigenetic modifiers KMT2D and KDM6A, which is implicated in 32% and 10% of Kabuki Syndrome patients respectively. To understand the connection to Kabuki syndrome patients, and the transcriptional targets of KMT2D and KDM6A in humans, we performed RNA sequencing (seq) of lymphoblastoid cells from Kabuki Syndrome patients carrying mutations in KMT2D and KDM6A. We identified 1995 significant changes in transcriptional targets for KMT2D and 1917 for KDM6A, as compared to control. When compared with RNA-seq datasets obtained from other mouse and zebrafish studies, our analysis revealed KMT2D mutations affect the expression of 76 orthologous genes across all three datasets. Similarliy, KDM6A afftects the expprssion of 7 orthologous genes across three datasets. Despite the differences in cell types, stages, and species in the comparison between the transcriptomic datasets, there are common gene expression changes associated with KMT2D and KDM6A mutations. qPCR on novel zebrafish mutants confirmed the differentially expression in KMT2D or KDM6A mutant backgrounds. Taken together, our results show that KMT2D and KDM6A regulate common and unique genes across humans, mice, and zebrafish for early craniofacial and cardiac development and this information contributes to the understanding of epigenetic dysregulation during development of Kabuki syndrome. Diverse stem-chondrichthyan oral structures and evidence for an independently acquired acanthodid dentition The teeth of sharks famously form a series of transversely-organised files with conveyor-belt replacement that are borne directly on the jaw cartilages, in contrast to the dermal plate-borne dentition of bony fishes that undergoes site-specific replacement. A major obstacle in understanding how this system evolved is the poorly understood relationships of the earliest chondrichthyans and the profusion of morphologically and terminologically diverse bones, cartilages, splints and whorls that they possess. Here we use tomographic methods to investigate mandibular structures in several early branching  acanthodian-grade stem-chondrichthyans. We show that the dentigerous jaw bones of disparate genera of ischnacanthids are united by a common construction, being growing bones with non-shedding dentition. Mandibular splints, which support the ventro-lateral edge of the Meckels cartilage in some taxa, are formed from dermal bone and may be an acanthodid synapomorphy. We demonstrate that the teeth of Acanthodopsis are borne directly on the mandibular cartilage and that this taxon is deeply nested within an edentulous radiation, representing an unexpected independent origin of teeth. Many or even all of the range of unusual oral structures may be apomorphic, but they should nonetheless be considered when building hypotheses of tooth and jaw evolution, both in chondrichthyans and more broadly. Design, Mutate, Screen: High-throughput creation of genetic clocks with different period-amplitude characteristics Directed evolution has become an invaluable tool in protein engineering and has also greatly influenced the construction of synthetic gene circuits. The ability to generate diversity at precise targets for directed evolution approaches has improved vastly, allowing researchers to create large, specific mutant libraries with relative ease. Screening approaches for large mutant libraries have similarly come a long way, especially when the desired behavior can easily be tested for with static, single time-point assays. For more complex gene circuits with dynamic phenotypes that change over time, directed evolution approaches to controlling and tuning circuit behavior have been hindered by the lack of sufficiently high-throughput screening methods to isolate variants with desired characteristics. Here we utilize directed mutagenesis and multiplexed microfluidics to develop a workflow for creating, screening and tuning dynamic gene circuits that operate at the population level. Specifically, we create a mutant library of an existing oscillator, the synchronized lysis circuit, and tune its dynamics while uncovering principles regarding its behavior. Lastly, we utilize this directed evolution workflow to construct a new synchronized genetic oscillator that exhibits robust dynamics over long time scales.