Yonggui Fu

Adaptive Evolution of Newly Emerged Micro-RNA Genes in Drosophila

How often micro-RNA (miRNA) genes emerged and how fast they evolved soon after their emergence are some of the central questions in the evolution of miRNAs. Because most known miRNA genes are ancient and highly conserved, these questions can be best answered by identifying newly emerged miRNA genes. Among the 78 miRNA genes in Drosophila reported before 2007, only 5 are confirmed to be newly emerged in the genus (although many more can be found in the newly reported data set; e.g., Ruby et al. 2007; Stark et al. 2007; Lu et al. 2008). These new miRNA genes have undergone numerous changes, even in the normally invariant mature sequences. Four of them (the miR-310/311/312/313 cluster, denoted miR-310s) were duplicated from other conserved miRNA genes. The fifth one (miR-303) appears to be a very young gene, originating de novo from a non-miRNA sequence recently. We sequenced these 5 miRNA genes and their neighboring regions from a worldwide collection of Drosophila melanogaster lines. The levels of divergence and polymorphism in these miRNA genes, vis-à-vis those of the neighboring DNA sequences, suggest that these 5 genes are evolving adaptively. Furthermore, the polymorphism pattern of miR-310s in D. melanogaster is indicative of hitchhiking under positive selection. Thus, a large number of adaptive changes over a long period of time may be essential for the evolution of newly emerged miRNA genes.

Decelerated genome evolution in modern vertebrates revealed by analysis of multiple lancelet genomes

Vertebrates diverged from other chordates ~500 Myr ago and experienced successful innovations and adaptations, but the genomic basis underlying vertebrate origins are not fully understood. Here we suggest, through comparison with multiple lancelet (amphioxus) genomes, that ancient vertebrates experienced high rates of protein evolution, genome rearrangement and domain shuffling and that these rates greatly slowed down after the divergence of jawed and jawless vertebrates. Compared with lancelets, modern vertebrates retain, at least relatively, less protein diversity, fewer nucleotide polymorphisms, domain combinations and conserved non-coding elements (CNE). Modern vertebrates also lost substantial transposable element (TE) diversity, whereas lancelets preserve high TE diversity that includes even the long-sought RAG transposon. Lancelets also exhibit rapid gene turnover, pervasive transcription, fastest exon shuffling in metazoans and substantial TE methylation not observed in other invertebrates. These new lancelet genome sequences provide new insights into the chordate ancestral state and the vertebrate evolution.

HaploMerger: Reconstructing allelic relationships for polymorphic diploid genome assemblies

Whole-genome shotgun assembly has been a long-standing issue for highly polymorphic genomes, and the advent of next-generation sequencing technologies has made the issue more challenging than ever. Here we present an automated pipeline, HaploMerger, for reconstructing allelic relationships in a diploid assembly. HaploMerger combines a LASTZ-ChainNet alignment approach with a novel graph-based structure, which helps to untangle allelic relationships between two haplotypes and guides the subsequent creation of reference haploid assemblies. The pipeline provides flexible parameters and schemes to improve the contiguity, continuity, and completeness of the reference assemblies. We show that HaploMerger produces efficient and accurate results in simulations and has advantages over manual curation when applied to real polymorphic assemblies (e.g., 4%-5% heterozygosity). We also used HaploMerger to analyze the diploid assembly of a single Chinese amphioxus (Branchiostoma belcheri) and compared the resulting haploid assemblies with EST sequences, which revealed that the two haplotypes are not only divergent but also highly complementary to each other. Taken together, we have demonstrated that HaploMerger is an effective tool for analyzing and exploiting polymorphic genome assemblies.

Genome-wide alternative polyadenylation in animals: insights from high-throughput technologies

Alternative polyadenylation (APA) plays an important role in gene expression by affecting mRNA stability, translation, and translocation in cells. However, genome-wide APA events have only recently been subjected to more systematic analysis with newly developed high-throughput methods. In this review, we focus on the recent technological development of APA analyses on a genome-wide scale, as well as the impact of APA switches on a number of critical biological processes in animals, including cell proliferation, differentiation, and oncogenic transformation. With the highly enlarged scope of genome-wide APA analyses, the APA regulations of various biological processes have increasingly become a new paradigm for the regulation of gene transcription and translation.

Complex genetic interactions underlying expression differences between Drosophila races: analysis of chromosome substitutions.

Regulation of gene expression is usually separated into cis and trans components. The separation may become artificial if much of the variation in expression is under multigenic and epistatic (e.g., cis-by-trans) control. There is hence a need to quantify the relative contribution of cis, trans, and cis-by-trans effects on expression divergence at different levels of evolution. To do so across the whole genome, we analyzed the full set of chromosome-substitution lines between the two behavioral races of Drosophila melanogaster. Our observations: (i) Only ≈3% of the genes with an expression difference are purely cis regulated. In fact, relatively few genes are governed by simple genetics because nearly 80% of expression differences are controlled by at least two chromosomes. (ii) For 14% of the genes, cis regulation does play a role but usually in conjunction with trans regulation. This joint action of cis and trans effects, either additive or epistatic, is referred to as inclusive cis effect. (iii) The percentage of genes with inclusive cis effect increases to 32% among genes that are strongly differentiated between the two races. (iv) We observed a nonrandom distribution of trans-acting factors, with a substantial deficit on the second chromosome. Between Drosophila racial groups, trans regulation of expression difference is extensive, and cis regulation often evolves in conjunction with trans effects.

APASdb: a database describing alternative poly(A) sites and selection of heterogeneous cleavage sites downstream of poly(A) signals

Increasing amounts of genes have been shown to utilize alternative polyadenylation (APA) 3′-processing sites depending on the cell and tissue type and/or physiological and pathological conditions at the time of processing, and the construction of genome-wide database regarding APA is urgently needed for better understanding poly(A) site selection and APA-directed gene expression regulation for a given biology. Here we present a web-accessible database, named APASdb (http://mosas.sysu.edu.cn/utr), which can visualize the precise map and usage quantification of different APA isoforms for all genes. The datasets are deeply profiled by the sequencing alternative polyadenylation sites (SAPAS) method capable of high-throughput sequencing 3′-ends of polyadenylated transcripts. Thus, APASdb details all the heterogeneous cleavage sites downstream of poly(A) signals, and maintains near complete coverage for APA sites, much better than the previous databases using conventional methods. Furthermore, APASdb provides the quantification of a given APA variant among transcripts with different APA sites by computing their corresponding normalized-reads, making our database more useful. In addition, APASdb supports URL-based retrieval, browsing and display of exon-intron structure, poly(A) signals, poly(A) sites location and usage reads, and 3′-untranslated regions (3′-UTRs). Currently, APASdb involves APA in various biological processes and diseases in human, mouse and zebrafish.

The Requirement of CD8+ T Cells To Initiate and Augment Acute Cardiac Inflammatory Response to High Blood Pressure

Macrophage infiltration and activation in myocardium are hallmarks of acute cardiac inflammatory response to high blood pressure. However, the underlying mechanisms remain elusive. In this article, we report that CD8+ T cells are required for cardiac recruitment and activation of macrophages. First, mice with CD8 gene-targeted (CD8 knockout) or CD8+ T cells depleted by Ab showed significantly reduced cardiac inflammatory response to the elevation of blood pressure after angiotensin II (Ang II) infusion, whereas CD8 knockout mice reconstituted with CD8+ T cells restored the sensitivity to Ang II. More importantly, CD8+ T cells were required for macrophage infiltration in myocardium and subsequent activation to express proinflammatory cytokines and chemokines. Furthermore, macrophage activation required direct contact with activated CD8+ T cells, but with TCR dispensable. TCR-independent activation of macrophages was further confirmed in MHC class I–restricted OVA-specific TCR transgenic mice, which showed a CD8+ T cell activation and cardiac proinflammatory response to Ang II similar to that of wild-type mice. Finally, only myocardium-infiltrated, but not peripheral, CD8+ T cells were specifically activated by Ang II, possibly by the cardiac IFN-γ that drove IFN-γR+ CD8+ T cell infiltration and activation. Thus, this work identified a TCR-independent innate nature of CD8+ T cells that was critical in initiating the sterile immune response to acute elevation of blood pressure.

An Ancient Balanced Polymorphism in a Regulatory Region of Human Major Histocompatibility Complex Is Retained in Chinese Minorities but Lost Worldwide

The coding regions of many of the major histocompatibility complex (MHC) (human leukocyte antigen [HLA] in humans) molecules are believed to be subject to balancing selection. But it is less certain whether the regulatory regions of such coding sequences are also subject to the same type of selection. Here, we studied the polymorphism of the regulatory regions of the HLA-DPA1 and HLA-DPB1 genes among ethnic minorities in southwestern China. Phylogenetic analysis revealed two deep clades >10 million years old. There is almost complete linkage disequilibrium between the regulatory and coding regions of HLA-DPA1, which hints at coadaptive balancing selection on the entire region. Thus, the molecular mechanism of balancing selection in MHC may involve expression modulation in addition to coding-region polymorphisms. Although the frequency of clade II is >30% in some ethnic minorities, it decreases to <5% among southern Han Chinese and vanishes among Europeans. As suspected, some ancient balanced polymorphisms, lost in major populations, still exist in isolated ethnicities. These isolated populations may thus contribute disproportionately to the total diversity of modern humans.

The role of alternative polyadenylation in the antiviral innate immune response

Alternative polyadenylation (APA) is an important regulatory mechanism of gene functions in many biological processes. However, the extent of 3′ UTR variation and the function of APA during the innate antiviral immune response are unclear. Here, we show genome-wide poly(A) sites switch and average 3′ UTR length shortens gradually in response to vesicular stomatitis virus (VSV) infection in macrophages. Genes with APA and mRNA abundance change are enriched in immune-related categories such as the Toll-like receptor, RIG-I-like receptor, JAK-STAT and apoptosis-related signalling pathways. The expression of 3′ processing factors is down-regulated upon VSV infection. When the core 3′ processing factors are knocked down, viral replication is affected. Thus, our study reports the annotation of genes with APA in antiviral immunity and highlights the roles of 3′ processing factors on 3′ UTR variation upon viral infection.

A global analysis of tandem 3'UTRs in eosinophilic chronic rhinosinusitis with nasal polyps.

Background Alternative polyadenylation (APA) is emerging as a widespread mechanism of gene regulation. The usage of APA sites allows a single gene to encode multiple mRNA transcripts with different 3′-untranslated region (3′UTR) lengths. Many disease processes reflect the importance of the regulation of APA site switching. The objective of this study was to explore the profiling of tandem APA sites in nasal polyps compared with nasal uncinate process mucosa. Methods Sequencing of APA sites (SAPAS) based on second-generation sequencing technology was undertaken to investigate the use of tandem APA sites and identify gene expression patterns in samples from the nasal polyps and nasal uncinate process mucosa of two patients with chronic rhinosinusitis with nasal polyps. The findings of the SAPAS analysis were validated via quantitative reverse-transcription polymerase chain reaction (qRT-PCR). Results First, the results showed a switching of 3′UTR lengths in nasal polyps compared with nasal uncinate process mucosa. From the two patients, 105 genes that were detected in both patients in the nasal polyps were switched to distal poly(A) sites, and 90 such genes were switched to proximal poly(A) sites. Several Gene Ontology terms were enriched in the list of genes with switched APA sites, including transcription regulation, cell cycle, apoptosis, and metabolism. Second, we detected genes that showed differential expression with at least a 3-fold difference between nasal polyp tissue and nasal uncinate process mucosa. Between the two sample types, 627 genes exhibited differential expression. The qRT-PCR results confirmed our SAPAS results. Conclusion APA site-switching events of 3′UTRs are prevalent in nasal polyp tissue, and the regulation of gene expression mediated by APA may play an important role in the formation and persistence of nasal polyps. Our results may provide new insights into the possible pathophysiologic processes involved in nasal polyps.