Hui Yu

Link-based quantitative methods to identify differentially coexpressed genes and gene Pairs

BackgroundDifferential coexpression analysis (DCEA) is increasingly used for investigating the global transcriptional mechanisms underlying phenotypic changes. Current DCEA methods mostly adopt a gene connectivity-based strategy to estimate differential coexpression, which is characterized by comparing the numbers of gene neighbors in different coexpression networks. Although it simplifies the calculation, this strategy mixes up the identities of different coexpression neighbors of a gene, and fails to differentiate significant differential coexpression changes from those trivial ones. Especially, the correlation-reversal is easily missed although it probably indicates remarkable biological significance.ResultsWe developed two link-based quantitative methods, DCp and DCe, to identify differentially coexpressed genes and gene pairs (links). Bearing the uniqueness of exploiting the quantitative coexpression change of each gene pair in the coexpression networks, both methods proved to be superior to currently popular methods in simulation studies. Re-mining of a publicly available type 2 diabetes (T2D) expression dataset from the perspective of differential coexpression analysis led to additional discoveries than those from differential expression analysis.ConclusionsThis work pointed out the critical weakness of current popular DCEA methods, and proposed two link-based DCEA algorithms that will make contribution to the development of DCEA and help extend it to a broader spectrum.

Transcript-level annotation of Affymetrix probesets improves the interpretation of gene expression data

BackgroundThe wide use of Affymetrix microarray in broadened fields of biological research has made the probeset annotation an important issue. Standard Affymetrix probeset annotation is at gene level, i.e. a probeset is precisely linked to a gene, and probeset intensity is interpreted as gene expression. The increased knowledge that one gene may have multiple transcript variants clearly brings up the necessity of updating this gene-level annotation to a refined transcript-level.ResultsThrough performing rigorous alignments of the Affymetrix probe sequences against a comprehensive pool of currently available transcript sequences, and further linking the probesets to the International Protein Index, we generated transcript-level or protein-level annotation tables for two popular Affymetrix expression arrays, Mouse Genome 430A 2.0 Array and Human Genome U133A Array. Application of our new annotations in re-examining existing expression data sets shows increased expression consistency among synonymous probesets and strengthened expression correlation between interacting proteins.ConclusionBy refining the standard Affymetrix annotation of microarray probesets from the gene level to the transcript level and protein level, one can achieve a more reliable interpretation of their experimental data, which may lead to discovery of more profound regulatory mechanism.

Slatkin’s Paradox: when direct observation and realized gene flow disagree. A case study in Ficus

In flying insects, there is frequently a lack of congruence between empirical estimates of local demographic parameters and the prediction that differentiation between populations should decrease with increasing dispersal, a puzzling phenomenon known as Slatkin’s Paradox. Here, we generalize Slatkin’s Paradox to other taxa, drawing from available information on dispersal to predict the relative importance of pollen vs. seed migration in structuring broad‐scale patterns of genetic variation in Ficus hirta, a dioecious fig whose pollen is dispersed by minute, species‐specific fig wasps and whose seeds are disseminated by a variety of vagile vertebrates (especially bats and birds). Local‐scale observational and genetic studies of dioecious understory figs suggest comparable rates of pollen and seed migration. In contrast, we found unusually low nuclear differentiation (FSTu2003=u20030.037, RSTu2003=u20030.074) and high chloroplast differentiation (GSTu2003=u20030.729, NSTu2003=u20030.798) among populations separated by up to 2850u2003km, leading us to reject the hypothesis of equal pollen and seed migration rates and to obtain an equilibrium estimate of the ratio of pollen to seed migration of ru2003=u200316.2–36.3. We reconcile this example of Slatkin’s Paradox with previously published data for dioecious figs and relate it to the picture of exceptionally long‐distance wasp‐mediated pollen dispersal that is emerging for large monoecious fig trees. More generally, we argue that Slatkin’s Paradox is a general phenomenon and suggest it may be common in plants and animals.

Nuclear and chloroplast DNA phylogeography of Ficus hirta: obligate pollination mutualism and constraints on range expansion in response to climate change.

This study uses a phylogeographic approach to investigate how interspecific interactions in an obligate pollination mutualism enhance or constrain dispersal and the range distributions of species through time. Fifteen populations of Ficus hirta, a bird-dispersed fig pollinated by a species-specific fig wasp, were sampled from Thailand to the northern limits of the tropical forest in China. These populations were assayed for six nuclear microsatellite loci and two intergenic chloroplast DNA sequences. Analyses of range expansion and genetic clustering indicated a relatively slow rate of range expansion from two or more southern glacial refugia. Low nuclear differentiation, combined with high interpopulation differentiation, and phylogeographic structuring of chloroplast variation indicated that seed dispersal has had a greater constraint than obligate interactions with fig wasps on the rate of post-glacial range expansion. This study is the first to investigate the phylogeographic history of a widely distributed southeast Asian tropical plant whose distribution extends to the northern limits of tropical forest habitat in China. It is also the first study of Ficus utilizing molecular data to evaluate whether species-specific pollination is a limitation or an aid to range expansion in response to climate change.

Systems perspectives on erythromycin biosynthesis by comparative genomic and transcriptomic analyses of S. erythraea E3 and NRRL23338 strains.

BackgroundS. erythraea is a Gram-positive filamentous bacterium used for the industrial-scale production of erythromycin A which is of high clinical importance. In this work, we sequenced the whole genome of a high-producing strain (E3) obtained by random mutagenesis and screening from the wild-type strain NRRL23338, and examined time-series expression profiles of both E3 and NRRL23338. Based on the genomic data and transcriptpmic data of these two strains, we carried out comparative analysis of high-producing strain and wild-type strain at both the genomic level and the transcriptomic level.ResultsWe observed a large number of genetic variants including 60 insertions, 46 deletions and 584 single nucleotide variations (SNV) in E3 in comparison with NRRL23338, and the analysis of time series transcriptomic data indicated that the genes involved in erythromycin biosynthesis and feeder pathways were significantly up-regulated during the 60xa0hours time-course. According to our data, BldD, a previously identified ery cluster regulator, did not show any positive correlations with the expression of ery cluster, suggesting the existence of alternative regulation mechanisms of erythromycin synthesis in S. erythraea. Several potential regulators were then proposed by integration analysis of genomic and transcriptomic data.ConclusionThis is a demonstration of the functional comparative genomics between an industrial S. erythraea strain and the wild-type strain. These findings help to understand the global regulation mechanisms of erythromycin biosynthesis in S. erythraea, providing useful clues for genetic and metabolic engineering in the future.

Moving Your Sons to Safety: Galls Containing Male Fig Wasps Expand into the Centre of Figs, Away From Enemies

Figs are the inflorescences of fig trees (Ficus spp., Moraceae). They are shaped like a hollow ball, lined on their inner surface by numerous tiny female flowers. Pollination is carried out by host-specific fig wasps (Agaonidae). Female pollinators enter the figs through a narrow entrance gate and once inside can walk around on a platform generated by the stigmas of the flowers. They lay their eggs into the ovules, via the stigmas and styles, and also gall the flowers, causing the ovules to expand and their pedicels to elongate. A single pollinator larva develops in each galled ovule. Numerous species of non-pollinating fig wasps (NPFW, belonging to other families of Chalcidoidea) also make use of galled ovules in the figs. Some initiate galls, others make use of pollinator-generated galls, killing pollinator larvae. Most NPFW oviposit from the outside of figs, making peripherally-located pollinator larvae more prone to attack. Style length variation is high among monoecious Ficus spp. and pollinators mainly oviposit into more centrally-located ovules, with shorter styles. Style length variation is lower in male (wasp-producing) figs of dioecious Ficus spp., making ovules equally vulnerable to attack by NPFW at the time that pollinators oviposit. We recorded the spatial distributions of galled ovules in mature male figs of the dioecious Ficus hirta in Southern China. The galls contained pollinators and three NPFW that kill them. Pollinators were concentrated in galls located towards the centre of the figs, NPFW towards the periphery. Due to greater pedicel elongation by male galls, male pollinators became located in more central galls than their females, and so were less likely to be attacked. This helps ensure that sufficient males survive, despite strongly female-biased sex ratios, and may be a consequence of the pollinator females laying mostly male eggs at the start of oviposition sequences.

Combinatorial network of transcriptional regulation and microRNA regulation in human cancer

BackgroundBoth transcriptional control and microRNA (miRNA) control are critical regulatory mechanisms for cells to direct their destinies. At present, the combinatorial regulatory network composed of transcriptional regulations and post-transcriptional regulations is often constructed through a forward engineering strategy that is based solely on searching of transcriptional factor binding sites or miRNA seed regions in the putative target sequences. If the reverse engineering strategy is integrated with the forward engineering strategy, a more accurate and more specific combinatorial regulatory network will be obtained.ResultsIn this work, utilizing both sequence-matching information and parallel expression datasets of miRNAs and mRNAs, we integrated forward engineering with reverse engineering strategies and as a result built a hypothetical combinatorial gene regulatory network in human cancer. The credibility of the regulatory relationships in the network was validated by random permutation procedures and supported by authoritative experimental evidence-based databases. The global and local architecture properties of the combinatorial regulatory network were explored, and the most important tumor-regulating miRNAs and TFs were highlighted from a topological point of view.ConclusionsBy integrating the forward engineering and reverse engineering strategies, we manage to sketch a genome-scale combinatorial gene regulatory network in human cancer, which includes transcriptional regulations and miRNA regulations, allowing systematic study of cancer gene regulation. Our work establishes a pipeline that can be extended to reveal conditional combinatorial regulatory landscapes correlating to specific cellular contexts.

Time-Resolved Transcriptome Analysis of Bacillus subtilis Responding to Valine, Glutamate, and Glutamine

Microorganisms can restructure their transcriptional output to adapt to environmental conditions by sensing endogenous metabolite pools. In this paper, an Agilent customized microarray representing 4,106 genes was used to study temporal transcript profiles of Bacillus subtilis in response to valine, glutamate and glutamine pulses over 24 h. A total of 673, 835, and 1135 amino-acid-regulated genes were identified having significantly changed expression at one or more time points in response to valine, glutamate, and glutamine, respectively, including genes involved in cell wall, cellular import, metabolism of amino-acids and nucleotides, transcriptional regulation, flagellar motility, chemotaxis, phage proteins, sporulation, and many genes of unknown function. Different amino acid treatments were compared in terms of both the global temporal profiles and the 5-minute quick regulations, and between-experiment differential genes were identified. The highlighted genes were analyzed based on diverse sources of gene functions using a variety of computational tools, including T-profiler analysis, and hierarchical clustering. The results revealed the common and distinct modes of action of these three amino acids, and should help to elucidate the specific signaling mechanism of each amino acid as an effector.

First record of an apparently rare fig wasp feeding strategy: obligate seed predation

1. Fig trees require host‐specific agaonid fig wasps for pollination, but their figs also support numerous non‐pollinating fig wasps (NPFW) that gall fig tissues or develop as parasitoids.

Lack of genetic isolation by distance, similar genetic structuring but different demographic histories in a fig- pollinating wasp mutualism

Historical abiotic factors such as climatic oscillations and extreme climatic events as well as biotic factors have shaped the structuring of species genetic diversity. In obligate species‐specific mutualisms, the biogeographic histories of the interacting species are tightly linked. This could be particularly true for nuclear genes in the Ficus‐pollinating wasp mutualistic association as the insects disperse pollen from their natal tree. In this study, we compare spatial genetic structure of plant and pollinator for the Ficus hirta–Valisia javana association throughout southeast China including Hainan Island, for both nuclear and cytoplasmic markers. We show that dispersal of the insect leads to plant and insect presenting similar signatures of lack of genetic isolation by distance for nuclear genes on the continent over a distance of 1000 km. But we also show that the demographic histories of plant and insect are strikingly different. This is in agreement with extreme climatic events leading to transient regional extinctions of the insects, associated with local survival of the plants. We also observe evidence of genetic differentiation for both wasps and fig‐tree between the continent and Hainan Island, although the Qiongzhou Strait is only on average 30 km wide, suggesting that geographic isolation by itself has not been sufficient to generate this differentiation. Hence, our results suggest that in highly dispersive mutualistic systems, isolation‐by‐dispersal limitation across a geographic barrier could be supplemented by isolation by adaptation, and maybe by coevolution, allowing further genetic divergence. In such systems, species may frequently be composed of a single population.