Identification of testis development-related genes by combining Iso-Seq and RNA-Seq in Zeugodacus tau.

Introduction: Zeugodacus tau (Walker) is an invasive pest. An effective method to control this pest is the sterile insect technique (SIT). To better apply this technique, it is necessary to understand testis development progression. Methods: Differentially expressed genes (DEGs) during testis development were analyzed by PacBio Iso-Seq and RNA-seq. Results: RNA-Seq library of Z. tau testes on day 1, 6, and 11 post eclosion were constructed. We identified 755 and 865 differentially expressed genes in the comparisons of T6 (testes on day 6) vs. T1 and T11 vs. T1, respectively. The KEGG pathway analysis showed that the DEGs were significantly enriched in retinol metabolism, vitamin B6 metabolism, and ascorbate and aldarate metabolism pathways. Knockdown of retinol dehydrogenase 12-like (rdh12-like), pyridoxal kinase (pdxk) and regucalcin (rgn), the representative gene in each of the above 3 pathways, reduced the hatching rate of Z. tau offspring. In addition, we identified 107 Drosophila spermatogenesis-related orthologous genes in Z. tau, of which innexin 2 (inx2) exhibited significantly up-regulated expression throughout testis development, and the knockdown of this gene reduced offspring hatching rate. Discussion: Our data indicated that rdh12-like, pdxk, rgn, and inx2 genes were related to testis development, and they were conserved in tephritid species. These results suggested that this gene might have the same function in tephritid. The findings provide an insight into testis development and spermatogenesis in tephritid species.

Identification and functional verification of Y-chromosome-specific gene typo-gyf in Bactrocera dorsalis.

Genes on the Y chromosome play important roles in male sex determination and development. The identification of Y-chromosome-specific genes not only provides a theoretical basis for the study of male reproductive development, but also offers genetic control targets for agricultural pests. However, Y-chromosome genes are rarely characterized due to their high repeatability and high heterochromatinization, especially in the oriental fruit fly. In this study, 1 011 Y-chromosome-specific candidate sequences were screened from 2 to 4 h Bactrocera dorsalis embryo datasets with the chromosome quotient method, 6 of which were identified as Y-chromosome-specific sequences by polymerase chain reaction, including typo-gyf, a 19 126-bp DNA sequence containing a 575-amino acid open reading frame. Testicular deformation and a significant reduction in sperm number were observed after typo-gyf knockdown with RNA interference in embryos. After typo-gyf knockout with clustered regularly interspaced palindromic repeats (CRISPR) / CRISPR-associated protein 9 in the embryonic stage, the sex ratio of the emergent adults was unbalanced, with far more females than males. A genotype analysis of these females with the Y-chromosome gene MoY revealed no sex reversal. Typo-gyf knockout led to the death of XY individuals in the embryonic stage. We conclude that typo-gyf is an essential gene for male survival, and is also involved in testicular development and spermatogenesis. The identification of typo-gyf and its functional verification provide insight into the roles of Y-chromosome genes in male development.

Comprehensive spatiotemporal mapping of single-cell lineages in developing mouse brain by CRISPR-based barcoding.

A fundamental interest in developmental neuroscience lies in the ability to map the complete single-cell lineages within the brain. To this end, we developed a CRISPR editing-based lineage-specific tracing (CREST) method for clonal tracing in Cre mice. We then used two complementary strategies based on CREST to map single-cell lineages in developing mouse ventral midbrain (vMB). By applying snapshotting CREST (snapCREST), we constructed a spatiotemporal lineage landscape of developing vMB and identified six progenitor archetypes that could represent the principal clonal fates of individual vMB progenitors and three distinct clonal lineages in the floor plate that specified glutamatergic, dopaminergic or both neurons. We further created pandaCREST (progenitor and derivative associating CREST) to associate the transcriptomes of progenitor cells in vivo with their differentiation potentials. We identified multiple origins of dopaminergic neurons and demonstrated that a transcriptome-defined progenitor type comprises heterogeneous progenitors, each with distinct clonal fates and molecular signatures. Therefore, the CREST method and strategies allow comprehensive single-cell lineage analysis that could offer new insights into the molecular programs underlying neural specification.

CRISPR/Cas9-induced Mutation of Sex Peptide Receptor Gene Bdspr Affects Ovary, Egg Laying, and Female Fecundity in Bactrocera dorsalis (Hendel) (Diptera: Tephritidae).

The oriental fruit fly, Bactrocera dorsalis (Hendel) (Diptera: Tephritidae), is an invasive and polyphagous pest of horticultural crops, and it can cause huge economic losses in agricultural production. The rapid development of CRISPR/Cas9 gene editing technology has provided new opportunities for the scientific control of agricultural pests. Here, we explore the applicability of the B. dorsalis sex peptide receptor (Bdspr) as a target gene for the CRISPR/Cas9-based sterile insect technique (SIT) in B. dorsalis. We screened two high-efficient single guide RNAs (sgRNAs) for gene editing. The results showed that both mutation efficiency and germline transmission rate were 100% in the surviving G0 females (8/8) from injected embryos, and that 75% of mosaically mutated G0 females (6/8) were sterile. The 50% of heterozygous G1 females (4/8) could not lay eggs; 100% of eggs laid by them could not survive; and 62.5% of individual females (5/8) had abnormal ovaries. These results indicate that Bdspr plays an important role in regulating fertility, egg viability, and ovary development in female B. dorsalis, suggesting that the spr gene can be used for CRISPR/Cas9-based SIT in B. dorsalis.

Integrated analysis of miRNA and mRNA expression profiles in response to gut microbiota depletion in the abdomens of female Bactrocera dorsalis.

Insect gut microbiota has been reported to participate in regulating host multiple biological processes including metabolism and reproduction. However, the corresponding molecular mechanisms remain largely unknown. Recent studies suggest that microRNAs (miRNAs) are involved in complex interactions between the gut microbiota and the host. Here, we used next-generation sequencing technology to characterize miRNA and mRNA expression profiles and construct the miRNA-gene regulatory network in response to gut microbiota depletion in the abdomens of female Bactrocera dorsalis. A total of 3016 differentially expressed genes (DEGs) and 18 differentially expressed miRNAs (DEMs) were identified. Based on the integrated analysis of miRNA and mRNA sequencing data, 229 negatively correlated miRNA-gene pairs were identified from the miRNA-mRNA network. Gene ontology enrichment analysis indicated that DEMs could target several genes involved in the metabolic process, oxidation-reduction process, oogenesis, and insulin signaling pathway. Finally, real-time quantitative polymerase chain reaction further verified the accuracy of RNA sequencing results. In conclusion, our study provides the profiles of miRNA and mRNA expressions under antibiotics treatment and provides an insight into the roles of miRNAs and their target genes in the interaction between the gut microbiota and its host.

miR-275/305 cluster is essential for maintaining energy metabolic homeostasis by the insulin signaling pathway in Bactrocera dorsalis.

Increasing evidence indicates that miRNAs play crucial regulatory roles in various physiological processes of insects, including systemic metabolism. However, the molecular mechanisms of how specific miRNAs regulate energy metabolic homeostasis remain largely unknown. In the present study, we found that an evolutionarily conserved miR-275/305 cluster was essential for maintaining energy metabolic homeostasis in response to dietary yeast stimulation in Bactrocera dorsalis. Depletion of miR-275 and miR-305 by the CRISPR/Cas9 system significantly reduced triglyceride and glycogen contents, elevated total sugar levels, and impaired flight capacity. Combined in vivo and in vitro experiments, we demonstrated that miR-275 and miR-305 can bind to the 3'UTR regions of SLC2A1 and GLIS2 to repress their expression, respectively. RNAi-mediated knockdown of these two genes partially rescued metabolic phenotypes caused by inhibiting miR-275 and miR-305. Furthermore, we further illustrated that the miR-275/305 cluster acting as a regulator of the metabolic axis was controlled by the insulin signaling pathway. In conclusion, our work combined genetic and physiological approaches to clarify the molecular mechanism of metabolic homeostasis in response to different dietary stimulations and provided a reference for deciphering the potential targets of physiologically important miRNAs in a non-model organism.

Human midbrain dopaminergic neuronal differentiation markers predict cell therapy outcomes in a Parkinson's disease model.

Human pluripotent stem cell-based (hPSC-based) replacement therapy holds great promise for the treatment of Parkinson's disease (PD). However, the heterogeneity of hPSC-derived donor cells and the low yield of midbrain dopaminergic (mDA) neurons after transplantation hinder its broad clinical application. Here, we have characterized the single-cell molecular landscape during mDA neuron differentiation. We found that this process recapitulated the development of multiple but adjacent fetal brain regions including the ventral midbrain, the isthmus, and the ventral hindbrain, resulting in a heterogenous donor cell population. We reconstructed the differentiation trajectory of the mDA lineage and identified calsyntenin 2 (CLSTN2) and protein tyrosine phosphatase receptor type O (PTPRO) as specific surface markers of mDA progenitors, which were predictive of mDA neuron differentiation and could facilitate high enrichment of mDA neurons (up to 80%) following progenitor cell sorting and transplantation. Marker-sorted progenitors exhibited higher therapeutic potency in correcting motor deficits of PD mice. Different marker-sorted grafts had a strikingly consistent cellular composition, in which mDA neurons were enriched, while off-target neuron types were mostly depleted, suggesting stable graft outcomes. Our study provides a better understanding of cellular heterogeneity during mDA neuron differentiation and establishes a strategy to generate highly purified donor cells to achieve stable and predictable therapeutic outcomes, raising the prospect of hPSC-based PD cell replacement therapies.

Does COVID-19 make the firms' performance worse? Evidence from the Chinese listed companies.

The influence of pandemics is still a black box, and the mechanism is attracting the attention of policymakers and scholars to guide the policy design in the aftermath of Covid-19 pandemics. This paper takes an in-depth look at the performance impact of pandemics from the perspective of operation, which is essential in a comprehensive evaluation of the economic effects of pandemics. With the help of novel quarterly data of Chinese listed firms from 2019 Q1 to 2021 Q2, we find that the Covid-19 decreases the sale-related profitability. For the mechanism, this paper finds that the pandemics make the operation longer, increase the cost, and reduce the potential cash flows. In addition, the environmental tax can significantly weaken the adverse shocks. The policy implication is that the sale boosting or consumption stimulus is vital in economic recovery, and the governments should efficiently use the positive effect of environmental tax.

Early embryonic transcriptomes of Zeugodacus tau provide insight into sex determination and differentiation genes.

Zeugodacus tau (Walker) is an invasive pest. The sterile insect technique is an environment-friendly method for pest control. Understanding the mechanism of sex determination will contribute to improving efficiency of this technique. In this study, we identified the transformer (tra) gene in Z. tau. One female-specific and two male-specific isoforms of tra were found in Z. tau, and the male-specific splicing pattern of tra was found to occur 5 h after egg laying. We performed transcriptome sequencing at 1 h (E1), 5 h (E5), and 9 h (E9) after egg laying and obtained high-quality transcriptome libraries of early embryo development. We identified 13 297 and 11 713 differentially expressed genes (DEGs) from E5 versus E1 and E9 versus E1 comparisons, respectively. To explore the potential functions of the DEGs during embryonic development, Gene Ontology, Clusters of Orthologous Groups of proteins, and Kyoto Encyclopedia of Genes and Genomes analyses were performed. Twenty-six genes potentially involved in sex determination or differentiation, including Maleness-on-the-Y (MoY), were identified in Z. tau. To verify the transcriptome results, 15 genes were selected for quantitative real-time PCR validation. The results were consistent with the transcriptome sequencing results. Moreover, U2 small nuclear riboprotein auxiliary factor (U2AF-50), female lethal d (fl(2)d), and virilizer (vir) were highly expressed at E5, indicating that they may be related to the sex-specific splicing of tra. Further functional analysis is needed to confirm this speculation. Our data provide an insight into the mechanism underlying sex determination and differentiation in tephritid species.

In silico promoter analysis and functional validation identify CmZFH, the co-regulator of hypoxia-responsive genes CmScylla and CmLPCAT.

Oxygen (O2) plays an essential role in aerobic organisms including terrestrial insects. Under hypoxic stress, the cowpea bruchid (Callosobruchus maculatus) ceases feeding and growth. However, larvae, particularly 4th instar larvae exhibit very high tolerance to hypoxia and can recover normal growth once brought to normoxia. To better understand the molecular mechanism that enables insects to cope with low O2 stress, we performed RNA-seq to distinguish hypoxia-responsive genes in midguts and subsequently identified potential common cis-elements in promoters of hypoxia-induced and -repressed genes, respectively. Selected elements were subjected to gel-shift and transient transfection assays to confirm their cis-regulatory function. Of these putative common cis-elements, AREB6 appeared to regulate the expression of CmLPCAT and CmScylla, two hypoxia-induced genes. CmZFH, the putative AREB6-binding protein, was hypoxia-inducible. Transient expression of CmZFH in Drosophila S2 cells activated CmLPCAT and CmScylla, and their induction was likely through interaction of CmZFH with AREB6. Binding to AREB6 was further confirmed by bacterially expressed CmZFH recombinant protein. Deletion analyses indicated that the N-terminal zinc-finger cluster of CmZFH was the key AREB6-binding domain. Through in silico and experimental exploration, we discovered novel transcriptional regulatory components associated with gene expression dynamics under hypoxia that facilitated insect survival.

Suggestions of management on emergency responses to epidemic outbreaks in outpatient or emergency departments of tertiary teaching hospitals.

Image, graphical abstract.

Comparative genomics of Klebsiella michiganensis BD177 and related members of Klebsiella sp. reveal the symbiotic relationship with Bactrocera dorsalis.

BACKGROUND: Bactrocera dorsalis is a destructive polyphagous and highly invasive insect pest of tropical and subtropical species of fruit and vegetable crops. The sterile insect technique (SIT) has been used for decades to control insect pests of agricultural, veterinary, and human health importance. Irradiation of pupae in SIT can reduce the ecological fitness of the sterile insects. Our previous study has shown that a gut bacterial strain BD177 that could restore ecological fitness by promoting host food intake and metabolic activities. RESULTS: Using long-read sequence technologies, we assembled the complete genome of K. michiganensis BD177 strain. The complete genome of K. michiganensis BD177 comprises one circular chromosome and four plasmids with a GC content of 55.03%. The pan-genome analysis was performed on 119 genomes (strain BD177 genome and 118 out of 128 published Klebsiella sp. genomes since ten were discarded). The pan-genome includes a total of 49305 gene clusters, a small number of 858 core genes, and a high number of accessory (10566) genes. Pan-genome and average nucleotide identity (ANI) analysis showed that BD177 is more similar to the type strain K. michiganensis DSM2544, while away from the type strain K. oxytoca ATCC13182. Comparative genome analysis with 21 K. oxytoca and 12 K. michiganensis strains, identified 213 unique genes, several of them related to amino acid metabolism, metabolism of cofactors and vitamins, and xenobiotics biodegradation and metabolism in BD177 genome. CONCLUSIONS: Phylogenomics analysis reclassified strain BD177 as a member of the species K. michiganensis. Comparative genome analysis suggested that K. michiganensis BD177 has the strain-specific ability to provide three essential amino acids (phenylalanine, tryptophan and methionine) and two vitamins B (folate and riboflavin) to B. dorsalis. The clear classification status of BD177 strain and identification of unique genetic characteristics may contribute to expanding our understanding of the symbiotic relationship of gut microbiota and B. dorsalis.

Precise and Rapid Validation of Candidate Gene by Allele Specific Knockout With CRISPR/Cas9 in Wild Mice.

It is a tempting goal to identify causative genes underlying phenotypic differences among inbred strains of mice, which is a huge reservoir of genetic resources to understand mammalian pathophysiology. In particular, the wild-derived mouse strains harbor enormous genetic variations that have been acquired during evolutionary divergence over 100s of 1000s of years. However, validating the genetic variation in non-classical strains was extremely difficult, until the advent of CRISPR/Cas9 genome editing tools. In this study, we first describe a T cell phenotype in both wild-derived PWD/PhJ parental mice and F1 hybrids, from a cross to C57BL/6 (B6) mice, and we isolate a genetic locus on Chr2, using linkage mapping and chromosome substitution mice. Importantly, we validate the identification of the functional gene controlling this T cell phenotype, Cd44, by allele specific knockout of the PWD copy, leaving the B6 copy completely intact. Our experiments using F1 mice with a dominant phenotype, allowed rapid validation of candidate genes by designing sgRNA PAM sequences that only target the DNA of the PWD genome. We obtained 10 animals derived from B6 eggs fertilized with PWD sperm cells which were subjected to microinjection of CRISPR/Cas9 gene targeting machinery. In the newborns of F1 hybrids, 80% (n = 10) had allele specific knockout of the candidate gene Cd44 of PWD origin, and no mice showed mistargeting of the B6 copy. In the resultant allele-specific knockout F1 mice, we observe full recovery of T cell phenotype. Therefore, our study provided a precise and rapid approach to functionally validate genes that could facilitate gene discovery in classic mouse genetics. More importantly, as we succeeded in genetic manipulation of mice, allele specific knockout could provide the possibility to inactivate disease alleles while keeping the normal allele of the gene intact in human cells.

A severe atherosclerosis mouse model on the resistant NOD background.

Atherosclerosis is a complex disease affecting arterial blood vessels and blood flow that could result in a variety of life-threatening consequences. Disease models with diverged genomes are necessary for understanding the genetic architecture of this complex disease. Non-obese diabetic (NOD) mice are highly polymorphic and widely used for studies of type 1 diabetes and autoimmunity. Understanding atherosclerosis development in the NOD strain is of particular interest as human atherosclerosis on the diabetic and autoimmune background has not been successfully modeled. In this study, we used CRISPR/Cas9 genome editing to genetically disrupt apolipoprotein E (ApoE) and low-density lipoprotein receptor (LDLR) expression on the pure NOD background, and compared phenotype between single-gene-deleted mice and double-knockout mutants with reference to ApoE-deficient C57BL/6 mice. We found that genetic ablation of Ldlr or Apoe in NOD mice was not sufficient to establish an atherosclerosis model, in contrast to ApoE-deficient C57BL/6 mice fed a high-fat diet (HFD) for over 12 weeks. We further obtained NOD mice deficient in both LDLR and ApoE, and assessed the severity of atherosclerosis and immune response to hyperlipidemia in comparison to ApoE-deficient C57BL/6 mice. Strikingly, the double-knockout NOD mice treated with a HFD developed severe atherosclerosis with aorta narrowed by over 60% by plaques, accompanied by destruction of pancreatic islets and an inflammatory response to hyperlipidemia. Therefore, we succeeded in obtaining a genetic model with severe atherosclerosis on the NOD background, which is highly resistant to the disease. This model is useful for the study of atherosclerosis in the setting of autoimmunity.

Speed genome editing by transient CRISPR/Cas9 targeting and large DNA fragment deletion.

Genetic engineering of cell lines and model organisms has been facilitated enormously by the CRISPR/Cas9 system. However, in cell lines it remains labor intensive and time consuming to obtain desirable mutant clones due to the difficulties in isolating the mutated clones and sophisticated genotyping. In this study, we have validated fluorescent protein reporter aided cell sorting which enables the isolation of maximal diversity in mutant cells. We further applied two spectrally distinct fluorescent proteins DsRed2 and ECFP as reporters for independent CRISPR/Cas9 mediated targeting, which allows for one-cell-one-well sorting of the mutant cells. Because of ultra-high efficiency of the CRISPR/Cas9 system with dual reporters and large DNA fragment deletion resulting from independent loci cleavage, monoclonal mutant cells could be easily identified by conventional PCR. In the speed genome editing method presented here, sophisticated genotyping methods are not necessary to identify loss of function mutations after CRISPR/Cas9 genome editing, and desirable loss of function mutant clones could be obtained in less than one month following transfection.

A genetically enhanced sterile insect technique against the fruit fly, Bactrocera dorsalis (Hendel) by feeding adult double-stranded RNAs.

RNAi based sterile insect technique (SIT) is an authentic insect management approach but requires proper target genes. During this study, spermless males were developed by interfering with germ cell differentiation and azoospermia related genes. Data demonstrates significant reductions in the target genes expressions (boul, zpg, dsx M , fzo and gas8) after oral dsRNAs administration. Knock down of target genes significantly affected the reproductive ability of males and reduced egg-hatching as compared to the control group. Furthermore, different combinations of selected gene dsRNAs (boul + zpg, boul + dsx M and zpg + dsx M ) were made, which resulted up to 85.40% of male sterility. The most effective combination was selected to prepare different concentrations of dsRNA, 250, 500, 750 and 1000 ng/µl, that caused 18.97%, 38.68%, 58.02% and 85.40% male sterility, respectively. Subsequently, 1000 ng/µl of the same combination of ds-RNAs was used against differently aged adult flies (1, 5, 7, 10 days) which lead to 85.40%, 31.42%, 21.76% and 9.90% male sterility, respectively. SIT developed in this study showed that, boul + zpg combination of dsRNA feeding for 6 hours significantly reduced the number of spermatozoa and viability of sperm in 1-day-old B. dorsalis flies. In short, this study provides an effective SIT technique for long-term B. dorsalis management.

A Seed Expansion Graph Clustering Method for Protein Complexes Detection in Protein Interaction Networks.

Most proteins perform their biological functions while interacting as complexes. The detection of protein complexes is an important task not only for understanding the relationship between functions and structures of biological network, but also for predicting the function of unknown proteins. We present a new nodal metric by integrating its local topological information. The metric reflects its representability in a larger local neighborhood to a cluster of a protein interaction (PPI) network. Based on the metric, we propose a seed-expansion graph clustering algorithm (SEGC) for protein complexes detection in PPI networks. A roulette wheel strategy is used in the selection of the seed to enhance the diversity of clustering. For a candidate node u, we define its closeness to a cluster C, denoted as NC(u, C), by combing the density of a cluster C and the connection between a node u and C. In SEGC, a cluster which initially consists of only a seed node, is extended by adding nodes recursively from its neighbors according to the closeness, until all neighbors fail the process of expansion. We compare the F-measure and accuracy of the proposed SEGC algorithm with other algorithms on Saccharomyces cerevisiae protein interaction networks. The experimental results show that SEGC outperforms other algorithms under full coverage.

The role of the transformer gene in sex determination and reproduction in the tephritid fruit fly, Bactrocera dorsalis (Hendel).

Transformer (tra) is a switch gene in the somatic sex-determination hierarchy that regulates sexual dimorphism based on RNA splicing in many insects. In tephritids, a Y-linked male determining gene (M) controls sex in the sex-determination pathway. Here, homologues of Drosophila tra and transformer-2 (tra-2) genes were isolated and characterized in Bactrocera dorsalis (Hendel), one of the most destructive agricultural insect pests in many Asian countries. Two male-specific and one female-specific isoforms of B. dorsalis transformer (Bdtra) were identified. The presence of multiple TRA/TRA-2 binding sites in Bdtra suggests that the TRA/TRA-2 proteins are splicing regulators promoting and maintaining, epigenetically, female sex determination by a tra positive feedback loop in XX individuals during development. The expression patterns of female-specific Bdtra transcripts during early embryogenesis shows that a peak appears at 15 h after egg laying. Using dsRNA to knock-down Bdtra expression in the embryo and adult stages, we showed that sexual formation is determined early in the embryo stage and that parental RNAi does not lead to the production of all male progeny as in Tribolium castaneum. RNAi results from adult abdominal dsRNA injections show that Bdtra has a positive influence on female yolk protein gene (Bdyp1) expression and fecundity.

The effects of RNA interference targeting Bactrocera dorsalis ds-Bdrpl19 on the gene expression of rpl19 in non-target insects.

Double-stranded RNA (dsRNA) designed to target pest genes emerges as a promising strategy for improving pest control. Therefore, it is necessary to assess the effects of dsRNA on non-target insects, such as native enemies and beneficial insects, to determine the environmental safety of such treatments. In this paper, we investigated the effects of dsRNA targeting rpl19 from Bactrocera dorsalis on non-target insects in citrus ecological systems by feeding the dsRNA to Bactrocera minax, Apis mellifera and Diachasmimorpha longicaudata. The results showed that when B. dorsalis were fed rpl19 CDS dsRNA or 3'UTR dsRNA, the expression of rpl19 was dramatically decreased. Feeding the Bdrpl19 CDS dsRNA to adult B. minax and D. longicaudata caused their respective rpl19 genes to be knocked down over 50-70 and 40%, respectively, but it had no effect on the expression of the rpl19 gene in A. mellifera. The Bdrpl19 3'UTR dsRNA did not have any silencing effects on the expression levels of rpl19 in non-target insects. This study provides evidence that dsRNA can impact non-target organisms, but the 3'UTR dsRNA may not have effects in non-target organisms.