The genome of the cereal pest Sitophilus oryzae: a transposable element haven BackgroundThe rice weevil Sitophilus oryzae is one of the most important agricultural pests, causing extensive damage to cereal in fields and to stored grains. S. oryzae has an intracellular symbiotic relationship (endosymbiosis) with the Gram-negative bacterium Sodalis pierantonius and is a valuable model to decipher host-symbiont molecular interactions.

ResultsWe sequenced the Sitophilus oryzae genome using a combination of short and long reads to produce the best assembly for a Curculionidae species to date. We show that S. oryzae has undergone successive bursts of transposable element (TE) amplification, representing 72% of the genome. In addition, we show that many TE families are transcriptionally active, and changes in their expression are associated with insect endosymbiotic state. S. oryzae has undergone a high gene expansion rate, when compared to other beetles. Reconstruction of host-symbiont metabolic networks revealed that, despite its recent association with cereal weevils (30 Kyear), S. pierantonius relies on the host for several amino acids and nucleotides to survive and to produce vitamins and essential amino-acids required for insect development and cuticle biosynthesis.

ConclusionsHere we present the genome of an agricultural pest beetle, which may act as a foundation for pest control. In addition, S. oryzae may be a useful model for endosymbiosis, and studying TE evolution and regulation, along with the impact of TEs on eukaryotic genomes.
Diversity of SARS-CoV-2 genome among various strain identified in Lucknow, Uttar Pradesh Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has emerged as a significant challenge worldwide. Rapid genome sequencing of SARS-CoV-2 is going on across the globe to detect mutations and genomic modifications in SARS-CoV-2. In this study, we have sequenced twenty-three SARS-CoV-2 positive samples collected during the first pandemic from the state of Uttar Pradesh, India. We observed a range of already reported mutations (2-22), including; D614G, L452R, Q613H, Q677H, T1027I in the S gene; S194L in the N gene; Q57H, L106F, T175I in the ORF3. Few unreported mutations such as P309S in the ORF1ab gene; T379I in the N gene; and L52F, V77I in the ORF3a gene were also detected. Phylogenetic genome analysis showed similarity with other SARS-CoV-2 viruses reported from Uttar Pradesh. The observed mutations may be associated with SARS-CoV-2 virus pathogenicity or disease severity.
Relationship between vestibular hair cell loss and deficits in two anti-gravity reflexes in the rat The tail-lift reflex and the air-righting reflex in rats are anti-gravity reflexes that depend on vestibular function. To begin identifying their cellular basis, this study examined the relationship between reflex loss and the graded lesions caused in the vestibular sensory epithelia by varying doses of an ototoxic compound. After ototoxic exposure, we recorded these reflexes using high speed video. The movies were used to obtain objective measures of the reflexes: the minimum angle formed by the nose, the back of the neck and the base of the tail during the tail-lift maneuver and the time to right in the air-righting test. The vestibular sensory epithelia were then collected from the rats and used to estimate the loss of type I (HCI), type II (HCII) and all hair cells (HC) in both central and peripheral parts of the crista, utricle, and saccule. As expected, tail-lift angles decreased, and air-righting times increased, while the numbers of HCs remaining in the epithelia decreased in a dose-dependent manner. The results demonstrated greater sensitivity of HCI compared to HCII to the IDPN ototoxicity, as well as a relative resiliency of the saccule compared to the crista and utricle. Comparing the functional measures with the cell counts, we observed that loss of the tail-lift reflex associates better with HCI than with HCII loss. In contrast, most HCI in the crista and utricle were lost before air-righting times increased. These data suggest that these reflexes depend on the function of non-identical populations of vestibular HCs. Feedforward regulatory logic underlies robustness of the specification-to-differentiation transition and fidelity of terminal cell fate during C. elegans endoderm development Development is driven by gene regulatory networks (GRNs) that progressively dictate specification and differentiation of cell fates. The architecture of GRNs directly determines the specificity and accuracy of developmental outcomes. We report here that the core regulatory circuitry for endoderm development in C. elegans is comprised of a recursive series of interlocked feedforward modules linking a cascade of six sequentially expressed GATA-type transcription factors. This structure results in a reiterated sequential redundancy, in which removal of a single factor or alternate factors in the cascade results in no, or a mild, effect on endoderm development and gut differentiation, while elimination of any two factors that are sequentially deployed in the cascade invariably results in a strong phenotype. The strength of the observed phenotypes is successfully predicted by a computational model based on the timing and levels of transcriptional states. The feedforward regulatory logic in the GRN appears to ensure timely onset of terminal differentiation genes and allows rapid and robust lockdown of cell fate during early embryogenesis. We further found that specification-to-differentiation transition is linked through a common regulator, the END-1 GATA factor that straddles the two processes. Finally, we revealed roles for key GATA factors in establishing spatial regulatory state domains by acting as transcriptional repressors that appear to define the boundaries of the digestive tract. Our findings support a comprehensive model of the core gene network that describes how robust endoderm development is achieved during C. elegans embryogenesis.

Graphic abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=65 SRC="FIGDIR/small/457588v2_ufig1.gif" ALT="Figure 1">
View larger version (23K):
org.highwire.dtl.DTLVardef@8515d6org.highwire.dtl.DTLVardef@15ca234org.highwire.dtl.DTLVardef@9d0fe7org.highwire.dtl.DTLVardef@17f28db_HPS_FORMAT_FIGEXP  M_FIG C_FIG Lifespan prolonging mechanisms and insulin upregulation without fat accumulation in long-lived reproductives of a higher termite Kings and queens of eusocial termites can live for decades, while queens sustain a nearly maximal fertility. To investigate the molecular mechanisms underlying their long lifespan, we carried out transcriptomics, lipidomics and metabolomics in Macrotermes natalensis on sterile short-lived workers, long-lived kings and five stages spanning twenty years of adult queen maturation. Reproductives share gene expression differences from workers in agreement with a reduction of several aging-related processes, involving upregulation of DNA damage repair and mitochondrial functions. Anti-oxidant gene expression is downregulated, while peroxidability of membranes in queens decreases. Against expectations, we observed an upregulated gene expression in fat bodies of reproductives of several components of the IIS pathway, including an insulin-like peptide, Ilp9. This pattern does not lead to deleterious fat storage in physogastric queens, while simple sugars dominate in their hemolymph and large amounts of resources are allocated towards oogenesis. Our findings support the notion that all processes causing aging need to be addressed simultaneously in order to prevent it. TCF11 Has a Potent Tumor-Repressing Effect than Its Prototypic Nrf1α by Definition of both Similar yet Different Regulatory Profiles, with a Striking Disparity from Nrf2 Nrf1 and Nrf2, as two principal CNC-bZIP transcription factors, regulate similar but different targets involved in a variety of biological functions for maintaining cell homeostasis and organ integrity. Of note, the unique topobiological behavior of Nrf1 makes its functions more complicated than Nrf2, because it is allowed for alternatively transcribing and selectively splicing to yield multiple isoforms (e.g., TCF11, Nrf1). In order to gain a better understanding of their similarities and differences in distinct regulatory profiles, all four distinct cell models for stably expressing TCF11, TCF11{Delta}N, Nrf1 or Nrf2 have been herein established by an Flp-In T-REx-293 system and then identified by transcriptomic sequencing. Further analysis revealed that Nrf1 and TCF11 have similar yet different regulatory profiles, although both contribute basically to positive regulation of their co-targets, which are disparate from those regulated by Nrf2. Such disparity in those gene regulation by Nrf1 and Nrf2 was further corroborated by scrutinizing comprehensive functional annotation of their specific and/or common target genes. Conversely, the mutant TCF11{Delta}N, resulting from a deletion of the N-terminal amino acids 2-156 from TCF11, resembles Nrf2 with the largely consistent structure and function. Interestingly, our further experimental evidence demonstrates that TCF11 acts as a potent tumor-repressor relative to Nrf1, albeit both isoforms possess a congruous capability to prevent malignant growth of tumor and upregulate those genes critical for improving the survival of patients with hepatocellular carcinoma. Macroscopic control of synchronous electrical signaling with chemically-excited gene expression Excitable cells can convert electrical signals into chemical outputs to facilitate the active transport of information across larger distances. This electrical-to-chemical conversion requires a tightly regulated expression of ion channels. Alterations of ion channel expression provide landmarks of numerous pathological diseases, such as cardiac arrhythmia, epilepsy, or cancer. Although the activity of ion channels can be locally regulated by external light or chemical stimulus, it remains challenging to coordinate the expression of ion channels on extended spatial-temporal scales in a non-invasive manner. Here, we have engineered yeast S. cerevisiae to read and convert local chemical concentrations into a dynamic electrical field distributed across cell populations. The core mechanism encodes a chemically-excitable dual-feedback gene circuit that precisely tunes the expression domain of potassium channels, globally coordinating cyclic firing of the plasma membrane potential (PMP). We demonstrate that this mechanism leverages an engineered constitutively open bacterial potassium channel KcsA to directly couple chemical stimuli with ion flux through gene expression and it can interface with the host ion channels through the pulsatile production of toxins. Our study provides a robust synthetic transcriptional toolbox underlying the conversion of local chemical environments into spatiotemporally organized electrical impulses for various cellular engineering, synthetic biology, and potential therapeutic applications. 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. Characterisation of protease activity during SARS-CoV-2 infection identifies novel viral cleavage sites and cellular targets with therapeutic potential SARS-CoV-2 is the causative agent behind the COVID-19 pandemic, and responsible for over 170 million infections, and over 3.7 million deaths worldwide. Efforts to test, treat and vaccinate against this pathogen all benefit from an improved understanding of the basic biology of SARS-CoV-2. Both viral and cellular proteases play a crucial role in SARS-CoV-2 replication, and inhibitors targeting proteases have already shown success at inhibiting SARS-CoV-2 in cell culture models. Here, we study proteolytic cleavage of viral and cellular proteins in two cell line models of SARS-CoV-2 replication using mass spectrometry to identify protein neo-N-termini generated through protease activity. We identify previously unknown cleavage sites in multiple viral proteins, including major antigenic proteins S and N, which are the main targets for vaccine and antibody testing efforts. We discovered significant increases in cellular cleavage events consistent with cleavage by SARS-CoV-2 main protease, and identify 14 potential high-confidence substrates of the main and papain-like proteases, validating a subset with in vitro assays. We showed that siRNA depletion of these cellular proteins inhibits SARS-CoV-2 replication, and that drugs targeting two of these proteins: the tyrosine kinase SRC and Ser/Thr kinase MYLK, showed a dose-dependent reduction in SARS-CoV-2 titres. Overall, our study provides a powerful resource to understand proteolysis in the context of viral infection, and to inform the development of targeted strategies to inhibit SARS-CoV-2 and treat COVID-19. Darting across space and time: Parametric modulators of sex-biased conditioned fear responding Pavlovian fear conditioning is a widely used behavioral paradigm for studying associative learning in rodents. Despite early recognition that subjects may engage in a variety of both conditioned and unconditioned responses, the last several decades have seen the field narrow its focus to measure freezing as the sole indicator of conditioned fear. We previously reported that female rats were more likely than males to engage in darting, an escape-like conditioned response that is associated with heightened shock reactivity. To determine how experimental parameters contribute to the frequency of darting in both males and females, we manipulated factors such as chamber size, shock intensity, and number of trials. To better capture fear-related behavioral repertoires in our animals, we developed ScaredyRat, an open-source custom Python tool that analyzes Noldus Ethovision-generated raw data files to identify Darters and quantify both conditioned and unconditioned responses. We find that like freezing, conditioned darting occurrences scale with experimental alterations. While most darting occurs in females, we find that with an extended training protocol, darting can emerge in males as well. Collectively, our data suggest that darting reflects a behavioral switch in conditioned responding that is a product of both an individual animals sex, shock reactivity, and experimental parameters, underscoring the need for careful consideration of sex as a biological variable in classic learning paradigms. Naturally enhanced neutralizing breadth to SARS-CoV-2 after one year Over one year after its inception, the coronavirus disease-2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) remains difficult to control despite the availability of several excellent vaccines. Progress in controlling the pandemic is slowed by the emergence of variants that appear to be more transmissible and more resistant to antibodies1,2. Here we report on a cohort of 63 COVID-19-convalescent individuals assessed at 1.3, 6.2 and 12 months after infection, 41% of whom also received mRNA vaccines3,4. In the absence of vaccination antibody reactivity to the receptor binding domain (RBD) of SARS-CoV-2, neutralizing activity and the number of RBD-specific memory B cells remain relatively stable from 6 to 12 months. Vaccination increases all components of the humoral response, and as expected, results in serum neutralizing activities against variants of concern that are comparable to or greater than neutralizing activity against the original Wuhan Hu-1 achieved by vaccination of naive individuals2,5-8. The mechanism underlying these broad-based responses involves ongoing antibody somatic mutation, memory B cell clonal turnover, and development of monoclonal antibodies that are exceptionally resistant to SARS-CoV-2 RBD mutations, including those found in variants of concern4,9. In addition, B cell clones expressing broad and potent antibodies are selectively retained in the repertoire over time and expand dramatically after vaccination. The data suggest that immunity in convalescent individuals will be very long lasting and that convalescent individuals who receive available mRNA vaccines will produce antibodies and memory B cells that should be protective against circulating SARS-CoV-2 variants. Embryo development in Carica papaya Linn. Papaya (Carica papaya Linn.) is a tropical plant whose draft genome has been sequenced. Papaya produces large fruits rich in vitamins A and C and is an important cash crop in developing countries. Nonetheless, little is known about how the female gametophyte develops, how it is fertilized and how it develops into a mature seed containing an embryo and an endosperm. The Papaya female gametophyte displays a Polygonum-type architecture consisting of two synergid cells, an egg cell, a central cell, and three antipodal cells. Reports are available of the presumed existence of varieties in which cross fertilization is bypassed and autonomous development of embryos occurs (e.g., apomixis). In this study, we analyzed the development of female gametophytes in a commercial Hawaiian parental line and in the presumed apomictic Costa Rican line L1. Samples were collected before and after anthesis to compare the overall structure, size and transcriptional patterns of several genes that may be involved in egg and endosperm cell fate and proliferation. These genes were the putative papaya homologs of ARGONAUTE9 (AGO9), MEDEA (MEA), RETINOBLASTOMA RELATED-1 (RBR1), and SLOW WALKER-1 (SWA1). Our results suggest that its feasibly to identify structural features of Polygonum-type development, and that in bagged female flowers of line L1 we might have observed autonomous development of embryo-like structures. Possible downregulation of papaya homologs for AGO9, MEA, RBR1 and SWA1 was observed in embryo sacs from line L1 before and after anthesis, which may suggest a link between suspected apomixis and transcriptional downregulation of genes for RNA-directed DNA methylation, histone remodelers, and rRNA processing. Most notably, the large size of the papaya embryo sac suggests that it could be a cytological alternative to Arabidopsis thaliana for study. Significant variation in embryo sac size was observed between the varieties under study, suggesting wide differences in the genetic regulation of anatomical features. Second Heart Field-derived Cells Contribute to Angiotensin II-mediated Ascending Aortopathies BackgroundThe ascending aorta is a common location for aneurysm and dissection. This aortic region is populated by a mosaic of medial and adventitial cells that are embryonically derived from either the second heart field (SHF) or the cardiac neural crest. SHF-derived cells populate areas that coincide with the spatial specificity of thoracic aortopathies. The purpose of this study was to determine whether and how SHF-derived cells contribute to ascending aortopathies.

MethodsAscending aortic pathologies were examined in patients with sporadic thoracic aortopathies and angiotensin II (AngII)-infused mice. Ascending aortas without overt pathology from AngII-infused mice were subjected to mass spectrometry assisted proteomics, and molecular features of SHF-derived cells were determined by single cell transcriptomic analyses. Genetic deletion of either low-density lipoprotein receptor-related protein 1 (Lrp1) or transforming growth factor-{beta} receptor 2 (Tgfbr2) in SHF- derived cells was conducted to examine the impact of SHF-derived cells on vascular integrity.

ResultsPathologies in human ascending aortic aneurysmal tissues were predominant in outer medial layers and adventitia. This gradient was mimicked in mouse aortas following AngII infusion that was coincident with the distribution of SHF-derived cells. Proteomics indicated that brief AngII infusion, prior to overt pathology, evoked downregulation of SMC proteins and differential expression of extracellular matrix proteins, including several LRP1 ligands. LRP1 deletion in SHF-derived cells augmented AngII-induced ascending aortic aneurysm and rupture. Single cell transcriptomic analysis revealed that brief AngII infusion decreased Lrp1 and Tgfbr2 mRNA abundance in SHF-derived cells and induced a unique fibroblast population with low abundance of Tgfbr2 mRNA. SHF-specific Tgfbr2 deletion led to embryonic lethality at E12.5 with dilatation of the outflow tract and retroperitoneal hemorrhage. Integration of proteomic and single cell transcriptomics results identified plasminogen activator inhibitor 1 (PAI1) as the most increased protein in SHF-derived SMCs and fibroblasts during AngII infusion. Immunostaining revealed a transmural gradient of PAI1 in both ascending aortas of AngII-infused mice and human ascending aneurysmal aortas that mimicked the gradient of medial and adventitial pathologies.

ConclusionSHF-derived cells exert a critical role in maintaining vascular integrity through LRP1 and TGF-{beta} signaling associated with increases of aortic PAI1.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=89 SRC="FIGDIR/small/930917v6_ufig1.gif" ALT="Figure 1">
View larger version (34K):
org.highwire.dtl.DTLVardef@185e3d9org.highwire.dtl.DTLVardef@131d702org.highwire.dtl.DTLVardef@11afa42org.highwire.dtl.DTLVardef@242786_HPS_FORMAT_FIGEXP  M_FIG C_FIG Clinical PerspectiveO_ST_ABSWhat is new?C_ST_ABS- SHF-derived SMCs and fibroblasts associate with AngII-induced aortic pathologies.
- AngII induces a distinct fibroblast sub-cluster that is less abundant for mRNAs related to major extracellular components and TGF{beta} ligands and receptors, but more abundant for proliferative genes.
- TGFBR2 deletion in SHF-derived cells are embryonic lethal with significant dilatation of the outflow tract in mice.
- SHF-specific deletion of LRP1 leads to aortic pathologies in mice, supporting the importance of SHF-derived cells in maintaining ascending aortic wall integrity.


What are the clinical implications?- Heterogeneity of the embryonic origins of SMCs and fibroblasts contributes to complex mechanisms of vasculopathy formation, which should be considered when investigating the pathogenesis of thoracic aortopathies. Holobiont Evolution: Population Theory for the Hologenome This article develops mathematical theory for the population dynamics of microbiomes with their hosts and theory for microbiome evolution caused by holobiont selection.

A holobiont is the union of a host with its microbiome. Holobiont selection is natural selection on the combined host and microbiome as a unit. The hologenome is the union of genes in the microbiome with genes in the host nucleus. The deep integration within many holobionts of microbiomes with their hosts extends to development, metabolism, physiology and behavior. The kind of eco-evolutionary processes producing this integration is unclear.

For a microbe to join a hosts microbiome, its population-dynamic parameters must sufficiently align with the hosts population-dynamic parameters. Accordingly, a microbe might be excluded from the microbiome, might enter the microbiome in a way leading to its subsequent gradual shedding from the holobiont, might join the microbiome and coexist with the host, or might, for a deleterious microbe, enter the microbiome culminating in the extinction of both it and the host.

A horizontally transmitted microbiome constitutes a genetic system that features "collective inheritance" capable of underwriting Darwinian descent with modification. The microbial source pool in the environment is the counterpart for the microbiome of the gamete pool for the nuclear genome. Poisson sampling of the microbial source pool is the counterpart of binomial sampling of the gamete pool. However, holobiont selection on the microbiome does not lead to a counterpart of the Hardy-Weinberg Law, a concept of selective neutrality, or directional selection that produces fixation of the microbial gene with the highest fitness. The effectiveness of holobiont selection on microbial genes is limited compared to the effectiveness of natural selection on nuclear genes. Nonetheless, holobiont selection can affect the composition of the microbiome by permitting the coexistence of a competitively inferior microbial strain with a superior strain in a colonization-extinction equilibrium and by permitting a polymorphism between altruistic and selfish microbial strains.

An altruistic microbe might strike an optimal balance between lowering within-host fitness to increase holobiont fitness. However, altruistic microbes are replaced by otherwise identical microbes that contribute nothing to holobiont fitness. This replacement can be reversed by hosts that discriminate against colonization by selfish microbes--likely brought about by hosts making antibodies against selfish microbes. This discrimination leads to a process of microbial species sorting orchestrated by the host immune system. Holobiont assembly with host-orchestrated species sorting (HOSS) followed by microbial competition, rather than coevolution or multi-level selection, is the likely cause of host-microbiome integration. P2X2 receptor subunit interfaces are missense variant hotspots where mutations tend to increase apparent ATP affinity Background and PurposeP2X receptors (P2XRs) are trimeric ligand-gated ion channels (LGICs) that open a cation-selective pore in response to ATP binding to their large extracellular domain (ECD). The seven known P2XR subtypes can assemble as homo- or heterotrimeric complexes and contribute to numerous physiological functions, including nociception, inflammation and hearing. The overall structure of P2XRs is well established, but little is known about the spectrum and prevalence of human genetic variations and the functional implications in specific domains.

Experimental ApproachHere we examine the impact of P2X2 receptor (P2X2R) inter-subunit interface missense variants identified in the human population or through structural predictions. We test both single and double mutants through electrophysiological and biochemical approaches.

Key resultsWe demonstrate that predicted ECD inter-subunit interfaces display a higher-than-expected density of missense variations and that the majority of mutations that disrupt putative inter-subunit interactions result in channels with higher apparent ATP affinity. Lastly, we show that double mutants at the subunit interface show significant energetic coupling, especially if located in close proximity.

Conclusions and ImplicationsWe provide the first structural mapping of the mutational burden across the human population in a LGIC and show that the density of missense mutations is constrained between protein domains, indicating evolutionary selection at the domain level. Our data may indicate that, unlike other LGICs, P2X2Rs have evolved an intrinsically high threshold for activation, possibly to allow for additional modulation or as a cellular protection mechanism against overstimulation.

Bullet point summary What is already known: O_LIP2X2 receptors are ATP-activated ion channels implicated in hearing and nociceptice pathways
C_LI
 What this study adds: O_LIA structural mapping of missense variants observed in the human population
C_LIO_LIWe identify the intersubunit-interface as a variant hotspot and decipher functional impact of mutations
C_LI
 Clinical significance: O_LIThe development of both inhibitors and activators of P2X2 receptor function may be required
C_LI Morphology of immatures of the thelytokous ant, Monomorium triviale Wheeler (Formicidae: Myrmicinae: Solenopsidini) with descriptions of the extraordinary last-instar queen larvae The ant genus Monomorium is one of the most species-rich but taxonomically problematic groups in the hyperdiverse subfamily Myrmicinae. An East Asian species, M. triviale Wheeler, produces both reproductive queens and sterile workers via obligate thelytokous parthenogenesis. Here, we describe the immature forms of M. triviale based on light and scanning electron microscopy observations, with a note on the striking caste dimorphism in the last larval instar. The last-instar queen larvae were easily recognized by their large size, "aphaenogastroid" body shape, and rows of doorknob-like tubercles on the lateral and dorsal body surface. This type of queen-specific structure has not been found in ants in general, let alone congeneric species found in Japan. In stark contrast to the queen larvae, worker larvae showed a "pheidoloid" body shape and a body surface similar to other ants. The worker larvae were estimated to have three instars, consistent with previously described congeners. The pupae of both castes had no cocoon, a characteristic commonly described in other Myrmicinae species. In total, the developmental period from egg to adult worker averaged 59 days under 25{degrees}C. We discuss possible functions of the tubercles of queen larvae based on previous studies. The cortical bone metabolome of C57BL/6J mice is sexually dimorphic Cortical bone quality, which is sexually dimorphic, depends on bone turnover and therefore the activities of remodeling bone cells. However, sex differences in cortical bone metabolism are not yet defined. Adding to the uncertainty about cortical bone metabolism, the metabolomes of whole bone, isolated cortical bone without marrow, and bone marrow have not been compared. We hypothesized that the metabolome of isolated cortical bone would be distinct from that of bone marrow and would reveal sex differences. Metabolite profiles from LC-MS of whole bone, isolated cortical bone, and bone marrow were generated from humeri from 20-week-old female C57Bl/6J mice. The cortical bone metabolomes were then compared for 20-week-old female and male C57Bl/6J mice. Femurs from male and female mice were evaluated for flexural material properties and were then categorized into bone strength groups. The metabolome of isolated cortical bone was distinct from both whole bone and bone marrow. We also found sex differences in the isolated cortical bone metabolome. Based on metabolite pathway analysis, females had higher lipid metabolism, and males had higher amino acid metabolism. High-strength bones, regardless of sex, had greater tryptophan and purine metabolism. For males, high strength bones had upregulated nucleotide metabolism, whereas lower strength bones had greater pentose phosphate pathway metabolism. Since the higher strength groups (females compared with males, high strength males compared with lower strength males) had higher serum CTX1/P1NP, we estimate that the metabolomic signature of bone strength in our study at least partially reflects differences in bone turnover. These data provide novel insight into bone bioenergetics and the sexual dimorphic nature of bone material properties in C57Bl/6 mice. Succinate metabolism uncouples retinal pigment epithelium mitochondria Succinate is central in energy metabolism. It is thought to be made and consumed within mitochondria, and not to cross the plasma membrane. Contrary to the canonical role of succinate, it is produced in the retina and released into the extracellular milleau. It is unknown how this affects adjacent retinal pigment epithelium, within eyecup tissue. By measuring consumption of 13C4-succinate and O2 we probe metabolism of the eyecup, we find robust succinate consumption by eyecup but not retina tissue. Succinate-dependant stimulation of eyecup respiration is oligomycin resistant; indicating that in eyecup tissue succinate uncouples mitochondrial ATP synthesis from respiration. Minimizing spectral overlap in multicolor flow cytometry experiments Selecting fluorochromes for polychromatic panels for flow cytometry is complex and time-consuming. Poorly designed panels can result in overlap between the emission spectra of the different fluorochromes, making their signals difficult to separate. While assessing all possible panels is simple to do programmatically, the combinatorial complexity of most real world problems renders brute force computation impractical. Here we present a novel complexity optimization algorithm for fast design of minimal spectral overlap fluorochrome combinations. To aid researchers in designing fluorochrome panels we implemented the algorithm in a web server, Spectracular, that also considers instrument laser configuration and allows users to define various constraints related to antibody-fluorochrome availability and marker co-expression. 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. Identifying and correcting for misspecifications in GWAS summary statistics and polygenic scores Publicly available genome-wide association studies (GWAS) summary statistics exhibit uneven quality, which can impact the validity of follow-up analyses. First, we present an overview of possible misspecifications that come with GWAS summary statistics. Then, in both simulations and real data analyses, we show that additional information such as imputation INFO scores, allele frequencies, and per-variant sample sizes in GWAS summary statistics can be used to detect possible issues and correct for misspecifications in the GWAS summary statistics. One important motivation for us is to improve the predictive performance of polygenic scores built from these summary statistics. Unfortunately, due to the lack of reporting standards for GWAS summary statistics, this additional information is not systematically reported. We also show that using well-matched LD references can improve model fit and translate into more accurate prediction. Finally, we discuss how to make polygenic score methods such as lassosum and LDpred2 more robust to these misspecifications to improve their predictive power. The compound of baicalein, wogonin, and oroxylin-A inhibiting EMT in A549 cell line via PI3K/AKT-TWIST1-glycolysis pathway Non-small cell lung cancer (NSCLC) is a worldwide disease with high morbidity and mortality, which is most derived from its metastasis. Some researches show that the epithelial-mesenchymal transition (EMT) process promotes lung cancer cells migration and invasion, leading to NSCLC metastasis. Total Flavonoid Aglycones Extract (TFAE) isolated from Scutellaria baicalensis was reported to inhibit tumor growth and induce apoptosis. In this study, we found that baicalein, wogonin, and oroxylin-A were the active compounds of TFAE. After reconstructing with these three compounds (baicalein (65.8%), wogonin (21.2%), and oroxylin-A (13.0%)), the reconstructed TFAE (reTFAE) present inhibitory effect on the EMT process of A549 cells. Then, bioinformatic technology was employed to elucidate the potential pharmacodynamic mechanism network of reTFAE. We identified the relationship between reTFAE and PI3K/Akt signaling pathways, with TWIST1 as the key protein. LY294002, the inhibitor of PI3K/Akt signaling pathway, and knock-down TWIST1 could significantly enhance the efficacy of reTFAE, with increasing expression of epithelial markers and decreasing expression of mesenchymal markers in A549 cells at the same time. Furthermore, stable isotope dimethyl-labeled proteomics technology was conducted to complement the follow-up mechanism that the EMT-inhibition process may be realized through glycolysis pathway. In conclusion, we claim that TWIST1-targeted flavonoids could provide a new strategy to inhibit EMT progress for the treatment of NSCLC. An End-to-End Deep Learning Speech Coding and Denoising Strategy for Cochlear Implants Cochlear implant (CI) users struggle to understand speech in noisy conditions. To address this problem, we propose a deep learning speech denoising sound coding strategy that estimates the CI electric stimulation patterns out of the raw audio data captured by the micro-phone, performing end-to-end CI processing. To estimate the relative denoising performance differences between various approaches, we compared this technique to a classic Wiener filter and to a convTasNet. Speech enhancement performance was assessed by means of signal-to-noise-ratio improvement and the short-time objective speech intelligibility measure. Additionally, 5 CI users were evaluated for speech intelligibility in noise to assess the potential benefits of each algorithm. Our results show that the proposed method is capable of replacing a CI sound coding strategy while preserving its general use for every listener and performing speech enhancement in noisy environments, without sacrificing algorithmic latency. A tool for assessing alignment of biomedical data repositories with open, FAIR, citation and trustworthy principles Increasing attention is being paid to the operation of biomedical data repositories in light of efforts to improve how scientific data is handled and made available for the long term. Multiple groups have produced recommendations for functions that biomedical repositories should support, with many using requirements of the FAIR data principles as guidelines. However, FAIR is but one set of principles that has arisen out of the open science community. They are joined by principles governing open science, data citation and trustworthiness, all of which are important aspects for biomedical data repositories to support. Together, these define a framework for data repositories that we call OFCT: Open, FAIR, Citable and Trustworthy. Here we developed an instrument using the open source PolicyModels toolkit that attempts to operationalize key aspects of OFCT principles and piloted the instrument by evaluating eight biomedical community repositories listed by the NIDDK Information Network (dkNET.org). Repositories included both specialist repositories that focused on a particular data type or domain, in this case diabetes and metabolomics, and generalist repositories that accept all data types and domains. The goal of this work was both to obtain a sense of how much the design of current biomedical data repositories align with these principles and to augment the dkNET listing with additional information that may be important to investigators trying to choose a repository, e.g., does the repository fully support data citation? The evaluation was performed from March to November 2020 through inspection of documentation and interaction with the sites by the authors. Overall, although there was little explicit acknowledgement of any of the OFCT principles in our sample, the majority of repositories provided at least some support for their tenets. Extended Graphical Lasso for Multiple Interaction Networks for High Dimensional Omics Data There has been a spate of interest in association networks in biological and medical research, for example, genetic interaction networks. In this paper, we propose a novel method, the extended joint hub graphical lasso (EDOHA), to estimate multiple related interaction networks for high dimensional omics data across multiple distinct classes. To be specific, we construct a convex penalized log likelihood optimization problem and solve it with an alternating direction method of multipliers (ADMM) algorithm. The proposed method can also be adapted to estimate interaction networks for high dimensional compositional data such as microbial interaction networks. The performance of the proposed method in the simulated studies shows that EDOHA has remarkable advantages in recognizing class-specific hubs than the existing comparable methods. We also present three applications of real datasets. Biological interpretations of our results confirm those of previous studies and offer a more comprehensive understanding of the underlying mechanism in disease.

Author summaryReconstruction of multiple association networks from high dimensional omics data is an important topic, especially in biology. Previous studies focused on estimating different networks and detecting common hubs among all classes. Integration of information over different classes of data while allowing difference in the hub nodes is also biologically plausible. Therefore, we propose a method, EDOHA, to jointly construct multiple interaction networks with capacity in finding different hub networks for each class of data. Simulation studies show the better performance over conventional methods. The method has been demonstrated in three real world data. Limited selection towards lighter seeds in Danish grasslands over an eight-year period Time-series analysis of cover data from 236 Danish grassland sites demonstrated community selection on seed mass. Across four grassland habitat types and during an eight-year period there was significant community selection towards lighter seeds. In the analysis, the continuous seed size variable was used to group plant species into functional types. This method was chosen in order to account for the sampling process of the abundance data, and is in contrast to most other analyses of trait selection, where the community weighted mean of the traits is used as the dependent variable. Generally, little is known on the underlying causes of possible plant trait community selection in semi-natural grasslands that are not subjected to successional processes.