Jiaofang Shao

Analysis of banana transcriptome and global gene expression profiles in banana roots in response to infection by race 1 and tropical race 4 of Fusarium oxysporum f. sp. cubense

BackgroundCavendish, the most widely grown banana cultivar, is relatively resistant to Race 1 of Fusarium oxysporum f. sp. cubense (Foc1) which caused widespread Panama disease during the first half of the 20th century but is susceptible to Tropical Race 4 of Foc (Foc TR4) which is threatening world banana production. The genome of the diploid species Musa acuminata which is the ancestor of a majority of triploid banana cultivars has recently been sequenced. Availability of banana transcriptomes will be highly useful for improving banana genome annotation and for biological research. The knowledge of global gene expression patterns influenced by infection of different Foc races will help to understand the host responses to the infection.ResultsRNA samples from different organs of the Cavendish cultivar were pooled for deep sequencing using the Illumina technology. Analysis of the banana transcriptome led to identification of over 842 genes that were not annotated by the Musa genome project. A large number of simple nucleotide polymorphisms (SNPs) and short insertions and deletion (indels) were identified from the transcriptome data. GFP-expressing Foc1 and Foc TR4 were used to monitor the infection process. Both Foc1 and Foc TR4 were found to be able to invade banana roots and spread to root vascular tissues in the first two days following inoculation. Digital gene expression (DGE) profiling analysis reveal that the infection by Foc1 and Foc TR4 caused very similar changes in the global gene expression profiles in the banana roots during the first two days of infection. The Foc infection led to induction of many well-known defense-related genes. Two genes encoding the ethylene biosynthetic enzyme ACC oxidase and several ethylene-responsive transcription factors (ERF) were among the strongly induced genes by both Foc1 and Foc TR4.ConclusionsBoth Foc1 and Foc TR4 are able to spread into the vascular system of banana roots during the early infection process and their infection led to similar gene expression profiles in banana roots. The transcriptome profiling analysis indicates that the ethylene synthetic and signalling pathways were activated in response to the Foc infection.

A Genome-Wide Hybrid Incompatibility Landscape between Caenorhabditis briggsae and C. nigoni

Systematic characterization of ẖybrid incompatibility (HI) between related species remains the key to understanding speciation. The genetic basis of HI has been intensively studied in Drosophila species, but remains largely unknown in other species, including nematodes, which is mainly due to the lack of a sister species with which C. elegans can mate and produce viable progeny. The recent discovery of a C. briggsae sister species, C. nigoni, has opened up the possibility of dissecting the genetic basis of HI in nematode species. However, the paucity of dominant and visible marker prevents the efficient mapping of HI loci between the two species. To elucidate the genetic basis of speciation in nematode species, we first generated 96 chromosomally integrated GFP markers in the C. briggsae genome and mapped them into the defined locations by PCR and Next-Generation Sequencing (NGS). Aided by the marker, we backcrossed the GFP-associated C. briggsae genomic fragments into C. nigoni for at least 15 generations and produced 111 independent introgressions. The introgression fragments cover most of the C. briggsae genome. We finally dissected the patterns of HI by scoring the embryonic lethality, larval arrest, sex ratio and male sterility for each introgression line, through which we identified pervasive HI loci and produced a genome-wide landscape of HI between the two nematode species, the first of its type for any non-Drosophila species. The HI data not only provided insights into the genetic basis of speciation, but also established a framework for the possible cloning of HI loci between the two nematode species. Furthermore, the data on hybrids confirmed Haldane’s rule and suggested the presence of a large X effect in terms of fertility between the two species. Importantly, this work opens a new avenue for studying speciation genetics between nematode species and allows parallel comparison of the HI with that in Drosophila and other species.

Unraveling Regulatory Programs for NF-kappaB, p53 and MicroRNAs in Head and Neck Squamous Cell Carcinoma

In head and neck squamous cell carcinoma (HNSCC), mutations of p53 usually coexist with aberrant activation of NF-kappaB (NF-κB), other transcription factors and microRNAs, which promote tumor pathogenesis. However, how these factors and microRNAs interact to globally modulate gene expression and mediate oncogenesis is not fully understood. We devised a novel bioinformatics method to uncover interactive relationships between transcription factors or microRNAs and genes. This approach is based on matrix decomposition modeling under the joint constraints of sparseness and regulator-target connectivity, and able to integrate gene expression profiling and binding data of regulators. We employed this method to infer the gene regulatory networks in HNSCC. We found that the majority of the predicted p53 targets overlapped with those for NF-κB, suggesting that the two transcription factors exert a concerted modulation on regulatory programs in tumor cells. We further investigated the interrelationships of p53 and NF-κB with five additional transcription factors, AP1, CEBPB, EGR1, SP1 and STAT3, and microRNAs mir21 and mir34ac. The resulting gene networks indicate that interactions among NF-κB, p53, and the two miRNAs likely regulate progression of HNSCC. We experimentally validated our findings by determining expression of the predicted NF-κB and p53 target genes by siRNA knock down, and by examining p53 binding activity on promoters of predicted target genes in the tumor cell lines. Our results elucidating the cross-regulations among NF-κB, p53, and microRNAs provide insights into the complex regulatory mechanisms underlying HNSCC, and shows an efficient approach to inferring gene regulatory programs in biological complex systems.

Systems‐level quantification of division timing reveals a common genetic architecture controlling asynchrony and fate asymmetry

Coordination of cell division timing is crucial for proper cell fate specification and tissue growth. However, the differential regulation of cell division timing across or within cell types during metazoan development remains poorly understood. To elucidate the systems‐level genetic architecture coordinating division timing, we performed a high‐content screening for genes whose depletion produced a significant reduction in the asynchrony of division between sister cells (ADS) compared to that of wild‐type during Caenorhabditis elegans embryogenesis. We quantified division timing using 3D time‐lapse imaging followed by computer‐aided lineage analysis. A total of 822 genes were selected for perturbation based on their conservation and known roles in development. Surprisingly, we find that cell fate determinants are not only essential for establishing fate asymmetry, but also are imperative for setting the ADS regardless of cellular context, indicating a common genetic architecture used by both cellular processes. The fate determinants demonstrate either coupled or separate regulation between the two processes. The temporal coordination appears to facilitate cell migration during fate specification or tissue growth. Our quantitative dataset with cellular resolution provides a resource for future analyses of the genetic control of spatial and temporal coordination during metazoan development.

Collaborative Regulation of Development but Independent Control of Metabolism by Two Epidermis-specific Transcription Factors in Caenorhabditis elegans

Background: NHR-25 and ELT-3 are required for development but not for initial specification of epidermis in C. elegans. Results: Genome-wide in vivo targets of NHR-25 are identified. Conclusion: NHR-25 and ELT-3 collaboratively regulate development but differentially control metabolism of epidermis. Significance: The results provide insight into how tissue-specific transcription factors enforce cell fate specification initiated by its master regulator. Cell fate specification is typically initiated by a master regulator, which is relayed by tissue-specific regulatory proteins (usually transcription factors) for further enforcement of cell identities, but how the factors are coordinated among each other to “finish up” the specification remains poorly understood. Caenorhabditis elegans epidermis specification is initiated by a master regulator, ELT-1, that activates its targets, NHR-25 and ELT-3, two epidermis-specific transcription factors that are important for development but not for initial specification of epidermis, thus providing a unique paradigm for illustrating how the tissue-specific regulatory proteins work together to enforce cell fate specification. Here we addressed the question through contrasting genome-wide in vivo binding targets between NHR-25 and ELT-3. We demonstrate that the two factors bind discrete but conserved DNA motifs, most of which remain in proximity, suggesting formation of a complex between the two. In agreement with this, gene ontology analysis of putative target genes suggested differential regulation of metabolism but coordinated control of epidermal development between the two factors, which is supported by quantitative analysis of expression of their specific or common targets in the presence or absence of either protein. Functional validation of a subset of the target genes showed both activating and inhibitory roles of NHR-25 and ELT-3 in regulating their targets. We further demonstrated differential control of specification of AB and C lineage-derived epidermis. The results allow us to assemble a comprehensive gene network underlying C. elegans epidermis development that is likely to be widely used across species and provides insights into how tissue-specific transcription factors coordinate with one another to enforce cell fate specification initiated by its master regulator.

The comprehensive transcriptional analysis in Caenorhabditis elegans by integrating ChIP-seq and gene expression data.

The fundamental step of learning transcriptional regulation mechanism is to identify the target genes regulated by transcription factors (TFs). Despite numerous target genes identified by chromatin immunoprecipitation followed by high-throughput sequencing technology (ChIP-seq) assays, it is not possible to infer function from binding alone in vivo. This is equally true in one of the best model systems, the nematode Caenorhabditis elegans (C. elegans), where regulation often occurs through diverse TF binding features of transcriptional networks identified in modENCODE. Here, we integrated ten ChIP-seq datasets with genome-wide expression data derived from tiling arrays, involved in six TFs (HLH-1, ELT-3, PQM-1, SKN-1, CEH-14 and LIN-11) with tissue-specific and four TFs (CEH-30, LIN-13, LIN-15B and MEP-1) with broad expression patterns. In common, TF bindings within 3 kb upstream of or within its target gene for these ten studies showed significantly elevated level of expression as opposed to that of non-target controls, indicated that these sites may be more likely to be functional through up-regulating its target genes. Intriguingly, expression of the target genes out of 5 kb upstream of their transcription start site also showed high levels, which was consistent with the results of following network component analysis. Our study has identified similar transcriptional regulation mechanisms of tissue-specific or broad expression TFs in C. elegans using ChIP-seq and gene expression data. It may also provide a novel insight into the mechanism of transcriptional regulation not only for simple organisms but also for more complex species.

Comparative proteome analysis between C . briggsae embryos and larvae reveals a role of chromatin modification proteins in embryonic cell division

Caenorhabditis briggsae has emerged as a model for comparative biology against model organism C. elegans. Most of its cell fate specifications are completed during embryogenesis whereas its cell growth is achieved mainly in larval stages. The molecular mechanism underlying the drastic developmental changes is poorly understood. To gain insights into the molecular changes between the two stages, we compared the proteomes between the two stages using iTRAQ. We identified a total of 2,791 proteins in the C. briggsae embryos and larvae, 247 of which undergo up- or down-regulation between the two stages. The proteins that are upregulated in the larval stages are enriched in the Gene Ontology categories of energy production, protein translation, and cytoskeleton; whereas those upregulated in the embryonic stage are enriched in the categories of chromatin dynamics and posttranslational modification, suggesting a more active chromatin modification in the embryos than in the larva. Perturbation of a subset of chromatin modifiers followed by cell lineage analysis suggests their roles in controlling cell division pace. Taken together, we demonstrate a general molecular switch from chromatin modification to metabolism during the transition from C. briggsae embryonic to its larval stages using iTRAQ approach. The switch might be conserved across metazoans.

Genome-Wide Introgression Revealed Pervasive Hybrid Incompatibilities (HI) between Caenorhabditis species

Systematic characterization of hybrid incompatibility (HI) between related species remains the key to understanding speciation. The genetic basis of HI has been intensively studied in Drosophila species, but remains largely unknown in other species, including nematodes. This is mainly due to the lack of a sister species with which C. elegans can mate and produce viable progeny. The recent discovery of a C. briggsae sister species, C. sp.9, opened up the possibility of dissecting the genetic basis of HI in nematode species. However, paucity of molecular and genetic tools has prevented the precise mapping of HI loci between the two species. To systematically isolate the HI loci between the nematode species pair, we first generated 96 chromosomally integrated, independent GFP insertions in the C. briggsae genome. We next mapped the GFP insertion site into defined locations using a method we had developed earlier. The dominant and visible markers facilitated the directional crossing of its linked genomic sequences into C. sp.9. We then backcrossed each individual marker into C. sp.9 for at least 15 generations and produced 111 independent introgression lines, which together represent most of the C. briggsae genome. We finally dissected the HI patterns by scoring embryonic lethality, larval arrest, sex ratio, fertility, male sterility and inviability in a subset of the introgression lines, and identified pervasive HIs between the two species. The study produced a genome-wide landscape of HI between nematode species for the first time. The initial crossing results confirmed the Haldane?s rule and the fertility data from homozygous introgressions supported the rule of large X effect. The large collection of introgression lines allows mapping of numerous HI loci into defined genomic regions between C. briggsae and C. sp.9, thus facilitating further characterization of their genetic and molecular mechanisms. Importantly, the study permits comparative analysis of speciation genetics between nematodes and other species.