Competing endogenous RNAs (ceRNAs) are RNA molecules that sequester shared microRNAs (miRNAs) thereby affecting the expression of other targets of the miRNAs. Whether genetic variants in ceRNA can affect its biological function and disease development is still an open question. Here we identified a large number of genetic variants that are associated with ceRNA's function using Geuvaids RNA-seq data for 462 individuals from the 1000 Genomes Project. We call these loci competing endogenous RNA expression quantitative trait loci or ‘cerQTL’, and found that a large number of them were unexplored in conventional eQTL mapping. We identified many cerQTLs that have undergone recent positive selection in different human populations, and showed that single nucleotide polymorphisms in gene 3΄UTRs at the miRNA seed binding regions can simultaneously regulate gene expression changes in both cis and trans by the ceRNA mechanism. We also discovered that cerQTLs are significantly enriched in traits/diseases associated variants reported from genome-wide association studies in the miRNA binding sites, suggesting that disease susceptibilities could be attributed to ceRNA regulation. Further in vitro functional experiments demonstrated that a cerQTL rs11540855 can regulate ceRNA function. These results provide a comprehensive catalog of functional non-coding regulatory variants that may be responsible for ceRNA crosstalk at the post-transcriptional level.

Infection after renal transplantation remains a major cause of morbidity and death, especially infection from the extensively drug-resistant bacteria, A. baumannii. A total of fourteen A. baumannii isolates were isolated from the donors’ preserved fluid from DCD (donation after cardiac death) renal transplantation and four isolates in the recipients’ draining liquid at the Kidney Disease Center, The First Affiliated Hospital, College of Medicine, Zhejiang University, from March 2013 to November 2014. An outbreak of A. baumannii emerging after DCD renal transplantation was tracked to understand the transmission of the pathogen. PFGE displayed similar DNA patterns between isolates from the same hospital. Antimicrobial susceptibility tests against thirteen antimicrobial agents were determined using the K-B diffusion method and eTest. Whole-genome sequencing was applied to investigate the genetic relationship of the isolates. With the clinical data and research results, we concluded that the A. baumannii isolates 3R1 and 3R2 was probably transmitted from the donor who acquired the bacteria during his stay in the ICU, while isolate 4R1 was transmitted from 3R1 and 3R2 via medical manipulation. This study demonstrated the value of integration of clinical profiles with molecular methods in outbreak investigation and their importance in controlling infection and preventing serious complications after DCD transplantation.

Whole genome sequencing (WGS) is a promising strategy to unravel variants or genes responsible for human diseases and traits. However, there is a lack of robust platforms for a comprehensive downstream analysis. In the present study, we first proposed three novel algorithms, sequence gap-filled gene feature annotation, bit-block encoded genotypes and sectional fast access to text lines to address three fundamental problems. The three algorithms then formed the infrastructure of a robust parallel computing framework, KGGSeq, for integrating downstream analysis functions for whole genome sequencing data. KGGSeq has been equipped with a comprehensive set of analysis functions for quality control, filtration, annotation, pathogenic prediction and statistical tests. In the tests with whole genome sequencing data from 1000 Genomes Project, KGGSeq annotated several thousand more reliable non-synonymous variants than other widely used tools (e.g. ANNOVAR and SNPEff). It took only around half an hour on a small server with 10 CPUs to access genotypes of ∼60 million variants of 2504 subjects, while a popular alternative tool required around one day. KGGSeq's bit-block genotype format used 1.5% or less space to flexibly represent phased or unphased genotypes with multiple alleles and achieved a speed of over 1000 times faster to calculate genotypic correlation.

Honey bees as other insects rely on the innate immune system for protection against diseases. The innate immune system includes the circulating hemocytes (immune cells) that clear pathogens from hemolymph (blood) by phagocytosis, nodulation or encapsulation. Honey bee hemocyte numbers have been linked to hemolymph levels of vitellogenin. Vitellogenin is a multifunctional protein with immune-supportive functions identified in a range of species, including the honey bee. Hemocyte numbers can increase via mitosis, and this recruitment process can be important for immune system function and maintenance. Here, we tested if hemocyte mediated phagocytosis differs among the physiologically different honey bee worker castes (nurses, foragers and winter bees), and study possible interactions with vitellogenin and hemocyte recruitment. To this end, we adapted phagocytosis assays, which—together with confocal microscopy and flow cytometry—allow qualitative and quantitative assessment of hemocyte performance. We found that nurses are more efficient in phagocytic uptake than both foragers and winter bees. We detected vitellogenin within the hemocytes, and found that winter bees have the highest numbers of vitellogenin-positive hemocytes. Connections between phagocytosis, hemocyte-vitellogenin and mitosis were worker caste dependent. Our results demonstrate that the phagocytic performance of immune cells differs significantly between honey bee worker castes, and support increased immune competence in nurses as compared to forager bees. Our data, moreover, provides support for roles of vitellogenin in hemocyte activity.

Honey bees (Apis mellifera) provide a system for studying social and food-related behavior. A caste of workers performs age-related tasks: young bees (nurses) usually feed the brood and other adult bees inside the nest, while older bees (foragers) forage outside for pollen, a protein/lipid source, or nectar, a carbohydrate source. The workers' transition from nursing to foraging and their foraging preferences correlate with differences in gustatory perception, metabolic gene expression, and endocrine physiology including the endocrine factors vitellogenin (Vg) and juvenile hormone (JH). However, the understanding of connections among social behavior, energy metabolism, and endocrine factors is incomplete. We used RNA interference (RNAi) to perturb the gene network of Vg and JH to learn more about these connections through effects on gustation, gene transcripts, and physiology. The RNAi perturbation was achieved by single and double knockdown of the genes ultraspiracle (usp) and vg, which encode a putative JH receptor and Vg, respectively. The double knockdown enhanced gustatory perception and elevated hemolymph glucose, trehalose, and JH. We also observed transcriptional responses in insulin like peptide 1 (ilp1), the adipokinetic hormone receptor (AKHR), and cGMP-dependent protein kinase (PKG, or “foraging gene” Amfor). Our study demonstrates that the Vg–JH regulatory module controls changes in carbohydrate metabolism, but not lipid metabolism, when worker bees shift from nursing to foraging. The module is also placed upstream of ilp1, AKHR, and PKG for the first time. As insulin, adipokinetic hormone (AKH), and PKG pathways influence metabolism and gustation in many animals, we propose that honey bees have conserved pathways in carbohydrate metabolism and conserved connections between energy metabolism and gustatory perception. Thus, perhaps the bee can make general contributions to the understanding of food-related behavior and metabolic disorders.

Modeling of transcriptional regulatory networks (TRNs) has been increasingly used to dissect the nature of gene regulation. Inference of regulatory relationships among transcription factors (TFs) and genes, especially among multiple TFs, is still challenging. In this study, we introduced an integrative method, LogicTRN, to decode TF–TF interactions that form TF logics in regulating target genes. By combining cis-regulatory logics and transcriptional kinetics into one single model framework, LogicTRN can naturally integrate dynamic gene expression data and TF-DNA-binding signals in order to identify the TF logics and to reconstruct the underlying TRNs. We evaluated the newly developed methodology using simulation, comparison and application studies, and the results not only show their consistence with existing knowledge, but also demonstrate its ability to accurately reconstruct TRNs in biological complex systems.

Evolutionary games model a common type of interactions in a variety of complex, networked, natural systems and social systems. Given such a system, uncovering the interacting structure of the underlying network is key to understanding its collective dynamics. Based on compressive sensing, we develop an efficient approach to reconstructing complex networks under game-based interactions from small amounts of data. The method is validated by using a variety of model networks and by conducting an actual experiment to reconstruct a social network. While most existing methods in this area assume oscillator networks that generate continuous-time data, our work successfully demonstrates that the extremely challenging problem of reverse engineering of complex networks can also be addressed even when the underlying dynamical processes are governed by realistic, evolutionary-game type of interactions in discrete time.

Linnorm is a novel normalization and transformation method for the analysis of single cell RNA sequencing (scRNA-seq) data. Linnorm is developed to remove technical noises and simultaneously preserve biological variations in scRNA-seq data, such that existing statistical methods can be improved. Using real scRNA-seq data, we compared Linnorm with existing normalization methods, including NODES, SAMstrt, SCnorm, scran, DESeq and TMM. Linnorm shows advantages in speed, technical noise removal and preservation of cell heterogeneity, which can improve existing methods in the discovery of novel subtypes, pseudo-temporal ordering of cells, clustering analysis, etc. Linnorm also performs better than existing DEG analysis methods, including BASiCS, NODES, SAMstrt, Seurat and DESeq2, in false positive rate control and accuracy.

Insect immune systems can recognize specific pathogens and prime offspring immunity. High specificity of immune priming can be achieved when insect females transfer immune elicitors into developing oocytes. The molecular mechanism behind this transfer has been a mystery. Here, we establish that the egg-yolk protein vitellogenin is the carrier of immune elicitors. Using the honey bee, Apis mellifera, model system, we demonstrate with microscopy and western blotting that vitellogenin binds to bacteria, both Paenibacillus larvae – the gram-positive bacterium causing American foulbrood disease – and to Escherichia coli that represents gram-negative bacteria. Next, we verify that vitellogenin binds to pathogen-associated molecular patterns; lipopolysaccharide, peptidoglycan and zymosan, using surface plasmon resonance. We document that vitellogenin is required for transport of cell-wall pieces of E. coli into eggs by imaging tissue sections. These experiments identify vitellogenin, which is distributed widely in oviparous species, as the carrier of immune-priming signals. This work reveals a molecular explanation for trans-generational immunity in insects and a previously undescribed role for vitellogenin.

It remains challenging to predict regulatory variants in particular tissues or cell types due to highly context-specific gene regulation. By connecting large-scale epigenomic profiles to expression quantitative trait loci (eQTLs) in a wide range of human tissues/cell types, we identify critical chromatin features that predict variant regulatory potential. We present cepip, a joint likelihood framework, for estimating a variant’s regulatory probability in a context-dependent manner. Our method exhibits significant GWAS signal enrichment and is superior to existing cell type-specific methods. Furthermore, using phenotypically relevant epigenomes to weight the GWAS single-nucleotide polymorphisms, we improve the statistical power of the gene-based association test.