Anciently duplicated genes continuously recruited to heart expression in vertebrate evolution are associated with heart chamber increase.

Although gene/genome duplications in the early stage of vertebrates have been thought to provide major resources of raw genetic materials for evolutionary innovations, it is unclear whether they continuously contribute to the evolution of morphological complexity during the course of vertebrate evolution, such as the evolution from two heart chambers (fishes) to four heart chambers (mammals and birds). We addressed this issue by our heart RNA-Seq experiments combined with published data, using 13 vertebrates and one invertebrate (sea squirt, as an outgroup). Our evolutionary transcriptome analysis showed that number of ancient paralogous genes expressed in heart tends to increase with the increase of heart chamber number along the vertebrate phylogeny, in spite that most of them were duplicated at the time near to the origin of vertebrates or even more ancient. Moreover, those paralogs expressed in heart exert considerably different functions from heart-expressed singletons: the former are functionally enriched in cardiac muscle and muscle contraction-related categories, whereas the latter play more basic functions of energy generation like aerobic respiration. These findings together support the notion that recruiting anciently paralogous genes that are expressed in heart is associated with the increase of chamber number in vertebrate evolution.

Identification of Germline Mutations in East-Asian Young Never-Smokers with Lung Adenocarcinoma by Whole-Exome Sequencing.

Recently, an increasing number of young never-smokers are diagnosed with lung cancer. The aim of this study is to investigate the genetic predisposition of lung cancer in these patients and discover candidate pathogenic variants for lung adenocarcinoma in young never-smokers. Peripheral blood was collected from 123 never-smoking east-Asian patients diagnosed with lung adenocarcinoma before the age of 40. Whole-exome sequencing (WES) was conducted on genomic DNA extracted from peripheral blood cells. As a result, 3,481 single nucleotide variants were identified. By bioinformatical tools and the published gene list associated with genetic predisposition of cancer, pathogenic variants were detected in ten germline genes: ATR, FANCD2, FANCE, GATA2, HFE, MSH2, PDGFRA, PMS2, SDHB, and WAS. Patients with pathogenic variants were more likely to occur in females (9/10, 90.0%) and have stage IV lung adenocarcinoma (4/10, 40%). Furthermore, germline mutations in 17 genes (ASB18, B3GALT5, CLEC4F, COL6A6, CYP4B1, C6orf132, EXO1, GATA4, HCK, KCP, NPHP4, PIGX, PPIL2, PPP1R3G, RRBP1, SALL4, and TTC28), which occurred in at least two patients, displayed potentially pathogenic effects. Gene ontology analysis further showed that these genes with germline mutations were mainly located in nucleoplasm and associated with DNA repair-related biological processes. The study provides spectrum of pathogenic variants and functional explanation for genetic predisposition of lung adenocarcinoma in young never-smokers, which sheds a light on prevention and early diagnosis of lung cancer. Supplementary Information: The online version contains supplementary material available at 10.1007/s43657-022-00062-1.

Draft Genome of White-blotched River Stingray Provides Novel Clues for Niche Adaptation and Skeleton Formation.

The white-blotched river stingray (Potamotrygon leopoldi) is a cartilaginous fish native to the Xingu River, a tributary of the Amazon River system. As a rare freshwater-dwelling cartilaginous fish in the Potamotrygonidae family in which no member has the genome sequencing information, P. leopoldi provides the evolutionary details in fish phylogeny, niche adaptation, and skeleton formation. In this study, we present its draft genome of 4.11 Gb comprised of 16,227 contigs and 13,238 scaffolds, with contig N50 of 3937 kb and scaffold N50 of 5675 kb in size. Our analysis shows that P. leopoldi is a slow-evolving fish that diverged from elephant sharks about 96 million years ago. Moreover, two gene families related to the immune system, immunoglobulin heavy constant delta genes and T-cell receptor alpha/delta variable genes, exhibit expantion in P. leopoldi only. We also identified the Hox gene clusters in P. leopoldi and discovered that seven Hox genes shared by five representative fish species are missing in P. leopoldi. The RNA sequencing data from P. leopoldi and other three fish species demonstrate that fishes have a more diversified tissue expression spectrum as compared with the corresponding mammalian data. Our functional studies suggest that the lack of the GC gene encoding vitamin D-binding protein in cartilaginous fishes (both P. leopoldi and Callorhinchus milii) could partly explain the absence of hard bone in their endoskeleton. Overall, this genome resource provides new insights into the niche adaptation, body plan, and skeleton formation of P. leopoldi as well as the genome evolution in cartilaginous fishes.

Biallelic mutations of TTC12 and TTC21B were identified in Chinese patients with multisystem ciliopathy syndromes.

BACKGROUND: Abnormalities in cilia ultrastructure and function lead to a range of human phenotypes termed ciliopathies. Many tetratricopeptide repeat domain (TTC) family members have been reported to play critical roles in cilium organization and function. RESULTS: Here, we describe five unrelated family trios with multisystem ciliopathy syndromes, including situs abnormality, complex congenital heart disease, nephronophthisis or neonatal cholestasis. Through whole-exome sequencing and Sanger sequencing confirmation, we identified compound heterozygous mutations of TTC12 and TTC21B in six affected individuals of Chinese origin. These nonsynonymous mutations affected highly conserved residues and were consistently predicted to be pathogenic. Furthermore, ex vivo cDNA amplification demonstrated that homozygous c.1464 + 2 T > C of TTC12 would cause a whole exon 16 skipping. Both mRNA and protein levels of TTC12 were significantly downregulated in the cells derived from the patient carrying TTC12 mutation c.1464 + 2 T > C by real-time qPCR and immunofluorescence assays when compared with two healthy controls. Transmission electron microscopy analysis further identified ultrastructural defects of the inner dynein arms in this patient. Finally, the effect of TTC12 deficiency on cardiac LR patterning was recapitulated by employing a morpholino-mediated knockdown of ttc12 in zebrafish. CONCLUSIONS: To the best of our knowledge, this is the first study reporting the association between TTC12 variants and ciliopathies in a Chinese population. In addition to nephronophthisis and laterality defects, our findings demonstrated that TTC21B should also be considered a candidate gene for biliary ciliopathy, such as TTC26, which further expands the phenotypic spectrum of TTC21B deficiency in humans.

Generation of vascularized brain organoids to study neurovascular interactions.

Brain organoids have been used to recapitulate the processes of brain development and related diseases. However, the lack of vasculatures, which regulate neurogenesis and brain disorders, limits the utility of brain organoids. In this study, we induced vessel and brain organoids, respectively, and then fused two types of organoids together to obtain vascularized brain organoids. The fused brain organoids were engrafted with robust vascular network-like structures and exhibited increased number of neural progenitors, in line with the possibility that vessels regulate neural development. Fusion organoids also contained functional blood-brain barrier-like structures, as well as microglial cells, a specific population of immune cells in the brain. The incorporated microglia responded actively to immune stimuli to the fused brain organoids and showed ability of engulfing synapses. Thus, the fusion organoids established in this study allow modeling interactions between the neuronal and non-neuronal components in vitro, particularly the vasculature and microglia niche.

Understanding how the organs form and how their cells behave is essential to finding the causes and treatment for developmental disorders, as well as understanding certain diseases. However, studying most organs in live animals or humans is technically difficult, expensive and invasive. To address this issue, scientists have developed models called 'organoids' that recapitulate the development of organs using stem cells in the lab. These models are easier to study and manipulate than the live organs. Brain organoids have been used to recapitulate brain formation as well as developmental, degenerative and psychiatric brain conditions such as microcephaly, autism and Alzheimer's disease. However, these brain organoids lack the vasculature (the network of blood vessels) that supplies a live brain with nutrients and regulates its development, and which has important roles in brain disorders. Partly due to this lack of blood vessels, brain organoids also do not develop a blood brain barrier, the structure that prevents certain contents of the blood, including pathogens, toxins and even certain drugs from entering the brain. These characteristics limit the utility of existing brain organoids. To overcome these limitations, Sun, Ju et al. developed brain organoids and blood vessel organoids independently, and then fused them together to obtain vascularized brain organoids. These fusion organoids developed a robust network of blood vessels that was well integrated with the brain cells, and produced more neural cell precursors than brain organoids that had not been fused. This result is consistent with the idea that blood vessels can regulate brain development. Analyzing the fusion organoids revealed that they contain structures similar to the blood-brain barrier, as well as microglial cells (immune cells specific to the brain). When exposed to lipopolysaccharide ­ a component of the cell wall of certain bacteria ­ these cells responded by initiating an immune response in the fusion organoids. Notably, the microglial cells were also able to engulf connections between brain cells, a process necessary for the brain to develop the correct structures and work normally. Sun, Ju et al. have developed a new organoid system that will be of broad interest to researchers studying interactions between the brain and the circulatory system. The development of brain-blood-barrier-like structures in the fusion organoids could also facilitate the development of drugs that can cross this barrier, making it easier to treat certain conditions that affect the brain. Refining this model to allow the fusion organoids to grow for longer times in the lab, and adding blood flow to the system will be the next steps to establish this system.

Biallelic DNAH9 mutations are identified in Chinese patients with defective left-right patterning and cilia-related complex congenital heart disease.

Defective left-right (LR) pattering results in a spectrum of laterality disorders including situs inversus totalis (SIT) and heterotaxy syndrome (Htx). Approximately, 50% of patients with primary ciliary dyskinesia (PCD) displayed SIT. Recessive variants in DNAH9 have recently been implicated in patients with situs inversus. Here, we describe six unrelated family trios and 2 sporadic patients with laterality defects and complex congenital heart disease (CHD). Through whole exome sequencing (WES), we identified compound heterozygous mutations in DNAH9 in the affected individuals of these family trios. Ex vivo cDNA amplification revealed that DNAH9 mRNA expression was significantly downregulated in these patients carrying biallelic DNAH9 mutations, which cause a premature stop codon or exon skipping. Transmission electron microscopy (TEM) analysis identified ultrastructural defects of the outer dynein arms in these affected individuals. dnah9 knockdown in zebrafish lead to the disturbance of cardiac left-right patterning without affecting ciliogenesis in Kupffer's vesicle (KV). By generating a Dnah9 knockout (KO) C57BL/6n mouse model, we found that Dnah9 loss leads to compromised cardiac function. In this study, we identified recessive DNAH9 mutations in Chinese patients with cardiac abnormalities and defective LR pattering.

Long noncoding RNA BS-DRL1 modulates the DNA damage response and genome stability by interacting with HMGB1 in neurons.

Long noncoding RNAs (lncRNAs) are known to regulate DNA damage response (DDR) and genome stability in proliferative cells. However, it remains unknown whether lncRNAs are involved in these vital biological processes in post-mitotic neurons. Here, we report and characterize a lncRNA, termed Brain Specific DNA-damage Related lncRNA1 (BS-DRL1), in the central nervous system. BS-DRL1 is a brain-specific lncRNA and depletion of BS-DRL1 in neurons leads to impaired DDR upon etoposide treatment in vitro. Mechanistically, BS-DRL1 interacts with HMGB1, a chromatin protein that is important for genome stability, and is essential for the assembly of HMGB1 on chromatin. BS-DRL1 mediated DDR exhibits cell-type specificity in the cortex and cerebellum in gamma-irradiated mice and BS-DRL1 knockout mice show impaired motor function and concomitant purkinje cell degeneration. Our study extends the understanding of lncRNAs in DDR and genome stability and implies a protective role of lncRNA against neurodegeneration.

Comprehensive Analysis of APA Events and Their Association With Tumor Microenvironment in Lung Adenocarcinoma.

BACKGROUND: Alternative polyadenylation (APA) is a pervasive posttranscriptional mechanism regulating gene expression. However, the specific dysregulation of APA events and its potential biological or clinical significance in lung adenocarcinoma (LUAD) remain unclear. METHODS: Here, we collected RNA-Seq data from two independent datasets: GSE40419 (n = 146) and The Cancer Genome Atlas (TCGA) LUAD (n = 542). The DaPars algorithm was employed to characterize the APA profiles in tumor and normal samples. Spearman correlation was used to assess the effects of APA regulators on 3' UTR changes in tumors. The Cox proportional hazard model was used to identify clinically relevant APA events and regulators. We stratified 512 patients with LUAD in the TCGA cohort through consensus clustering based on the expression of APA factors. FINDINGS: We identified remarkably consistent alternative 3' UTR isoforms between the two cohorts, most of which were shortened in LUAD. Our analyses further suggested that aberrant usage of proximal polyA sites resulted in escape from miRNA binding, thus increasing gene expression. Notably, we found that the 3' UTR lengths of the mRNA transcriptome were correlated with the expression levels of APA factors. We further identified that CPSF2 and CPEB3 may serve as key regulators in both datasets. Finally, four LUAD subtypes according to different APA factor expression patterns displayed distinct clinical results and oncogenic features related to tumor microenvironment including immune, metabolic, and hypoxic status. INTERPRETATION: Our analyses characterize the APA profiles among patients with LUAD and identify two key regulators for APA events in LUAD, CPSF2 and CPEB3, which could serve as the potential prognostic genes in LUAD.

Convergent degeneration of olfactory receptor gene repertoires in marine mammals.

BACKGROUND: Olfactory receptors (ORs) can bind odor molecules and play a crucial role in odor sensation. Due to the frequent gains and losses of genes during evolution, the number of OR members varies greatly among different species. However, whether the extent of gene gains/losses varies between marine mammals and related terrestrial mammals has not been clarified, and the factors that might underlie these variations are unknown. RESULTS: To address these questions, we identified more than 10,000 members of the OR family in 23 mammals and classified them into 830 orthologous gene groups (OGGs) and 281 singletons. Significant differences occurred in the number of OR repertoires and OGGs among different species. We found that all marine mammals had fewer OR genes than their related terrestrial lineages, with the fewest OR genes found in cetaceans, which may be closely related to olfactory degradation. ORs with more gene duplications or loss events tended to be under weaker purifying selection. The average gain and loss rates of OR genes in terrestrial mammals were higher than those of mammalian gene families, while the average gain and loss rates of OR genes in marine mammals were significantly lower and much higher than those of mammalian gene families, respectively. Additionally, we failed to detect any one-to-one orthologous genes in the focal species, suggesting that OR genes are not well conserved among marine mammals. CONCLUSIONS: Marine mammals have experienced large numbers of OR gene losses compared with their related terrestrial lineages, which may result from the frequent birth-and-death evolution under varied functional constrains. Due to their independent degeneration, OR genes present in each lineage are not well conserved among marine mammals. Our study provides a basis for future research on the olfactory receptor function in mammals from the perspective of evolutionary trajectories.

Conversion of Astrocytes and Fibroblasts into Functional Noradrenergic Neurons.

Dysfunction of noradrenergic (NA) neurons is associated with a number of neuronal disorders. Diverse neuronal subtypes can be generated by direct reprogramming. However, it is still unknown how to convert non-neuronal cells into NA neurons. Here, we show that seven transcription factors (TFs) (Ascl1, Phox2b, AP-2α, Gata3, Hand2, Nurr1, and Phox2a) are able to convert astrocytes and fibroblasts into induced NA (iNA) neurons. These iNA neurons express the genes required for the biosynthesis, release, and re-uptake of noradrenaline. Moreover, iNA neurons fire action potentials, receive synaptic inputs, and control the beating rate of co-cultured ventricular myocytes. Furthermore, iNA neurons survive and integrate into neural circuits after transplantation. Last, human fibroblasts can be converted into functional iNA neurons as well. Together, iNA neurons are generated by direct reprogramming, and they could be potentially useful for disease modeling and cell-based therapies.

The prognostic signature of the somatic mutations in Ewing sarcoma: from a network view.

BACKGROUND: Ewing sarcoma is a malignant bone tumor mainly affecting teenagers and young adults. Its main driver mutation, the EWS-FLI1 fusion gene, has been identified more than 20 years ago, whereas its other somatic mutations have been just recently reported. METHODS: In this study, we organized the somatic mutations from 216 Ewing sarcoma cases into 216 individual protein-protein interaction networks by using interactome information. These mutation networks were then classified into five different clusters based on their structural similarities. The prognostic effect of mutation genes was evaluated according to their network features. RESULTS: The cases in cluster two exhibited remarkably high metastasis and mortality rates, and STAG2, TP53 and TTN were the three most significantly mutated genes in this cluster. Microarray data demonstrate that the expression of STAG2, TP53 and TTN are down-regulated in the EWS-FLI1-knockdown Ewing sarcoma cells. However, the mutation effect analysis shows that the somatic mutations in TTN are less damaging than those in STAG2 and TP53. The analyses of functional network modules further revealed that STAG2, TP53 and their interacting gene partners participate in the oncogenic-related biological modules such as cell cycle and regulation of transcription from RNA polymerase II promoter while TTN, TP53 and their interacting gene partners constitute the modules less relevant to oncogenesis. The results of Gene Ontology analyses demonstrated that the uniquely mutated genes associated with poor prognosis in Clusters 1, 4 and 5 were distinctively enriched in epidermal growth factor-related functions and phosphoproteins. CONCLUSIONS: Our study identified the highly lethal mutation combination cases and characterized the possible prognostic genes in Ewing sarcoma from a network perceptive.

Integrative analysis identifies DNMTs against immune-infiltrating neutrophils and dendritic cells in colorectal cancer.

Molecular characterizations, including microsatellite instability (MSI) and the CpG island methylator phenotype (CIMP) showed strong associations in colorectal carcinoma (CRC) and provided a deeper understanding of the etiology of disease. However, the global relationship between epigenetic alternations and changes in mRNA expression in CRC remains largely undefined, especially regarding the roles of DNA methyltransferases (DNMTs). Here, we conducted a systematic network comparison to explore the global conservation between co-expressed and co-methylated modules. We successfully identified immune-related modules that were regulated by DNMTs and had strong associations with immune-infiltrating neutrophils and dendritic cells in CRC. Moreover, we found that genes in those modules were prognostic for CRC, with 97.1% (168/173) being significantly influenced by DNMTs. Thus, this study resolved an interaction between DNA methylation and mRNA expression through DNMTs. Additionally, we provided evidence that DNMTs control the global hypomethylation of oncogenes, including ALOX5AP and CSF3R that otherwise have high methylation in normal colons. Such genes were also more sensitive to DNMT changes, such as in CRC. Collectively, our analyzes provided a systems biology approach to investigate the association among different molecular phenotypes in diseases.

Single-cell transcriptomes reveal molecular specializations of neuronal cell types in the developing cerebellum.

The cerebellum is critical for controlling motor and non-motor functions via cerebellar circuit that is composed of defined cell types, which approximately account for more than half of neurons in mammals. The molecular mechanisms controlling developmental progression and maturation processes of various cerebellar cell types need systematic investigation. Here, we analyzed transcriptome profiles of 21119 single cells of the postnatal mouse cerebellum and identified eight main cell clusters. Functional annotation of differentially expressed genes revealed trajectory hierarchies of granule cells (GCs) at various states and implied roles of mitochondrion and ATPases in the maturation of Purkinje cells (PCs), the sole output cells of the cerebellar cortex. Furthermore, we analyzed gene expression patterns and co-expression networks of 28 ataxia risk genes, and found that most of them are related with biological process of mitochondrion and around half of them are enriched in PCs. Our results also suggested core transcription factors that are correlated with interneuron differentiation and characteristics for the expression of secretory proteins in glia cells, which may participate in neuronal modulation. Thus, this study presents a systematic landscape of cerebellar gene expression in defined cell types and a general gene expression framework for cerebellar development and dysfunction.

The evolution of shame and guilt.

Shame and guilt seem to be two synonymous moral emotions but actually lead to contrasting human behaviors or behavioral tendencies. Shame drives people to hide or deny their wrongdoings while guilt drives people to amend their mistakes. How shame and guilt evolved in humans is still obscure. Here we present a computer model featured with reciprocal altruism and gregarious lifestyle for studying this question. We tested ten different strategies in our model and the pairwise contests show that shame-driven-hiding strategy can dominate the other strategies such as tit-for-tat and Pavlov in more than half of parameter combinations. The mathematical analysis of our model demonstrates that shame-driven-hiding strategy is an evolutionary stable strategy within a group as long as hiding can let an individual evade the retaliations to his wrongdoings. However, the problem of hiding is that it reduces an individual's social circle, i.e. living in a smaller group. Our analysis also shows that guilt-driven-amending strategy can outperform shame-driven-denying strategy at both individual and group levels if the cooperative behavior is sustainable within a group (b/(b-c) < T/n). Thus, we propose that shame is more adaptive at the individual level while guilt is more advantageous in the context of intergroup competition.

Double agents: genes with both oncogenic and tumor-suppressor functions.

The role of genetic components in cancer development is an area of interest for cancer biologists in general. Intriguingly, some genes have both oncogenic and tumor-suppressor functions. In this study, we systematically identified these genes through database search and text mining. We find that most of them are transcription factors or kinases and exhibit dual biological functions, e.g., that they both positively and negatively regulate transcription in cells. Some cancer types such as leukemia are over-represented by them, whereas some common cancer types such as lung cancer are under-represented by them. Across 12 major cancer types, while their genomic mutation patterns are similar to that of oncogenes, their expression patterns are more similar to that of tumor-suppressor genes. Their expression profile in six human organs propose that they mainly function as tumor suppressor in normal tissue. Our network analyses further show they have higher network degrees than both oncogenes and tumor-suppressor genes and thus tend to be the hub genes in the protein-protein interaction network. Our mutation, expression spectrum, and network analyses might help explain why some cancer types are specifically associated with them. Finally, our results suggest that the functionally altering mutations in "double-agent" genes and oncogenes are the main driving force in cancer development, because non-silent mutations are biasedly distributed toward these two gene sets across all 12 major cancer types.

Evolutionary Dynamics of MERS-CoV: Potential Recombination, Positive Selection and Transmission.

Middle East respiratory syndrome coronavirus (MERS-CoV) belongs to beta group of coronavirus and was first discovered in 2012. MERS-CoV can infect multiple host species and cause severe diseases in human. We conducted a series of phylogenetic and bioinformatic analyses to study the evolution dynamics of MERS-CoV among different host species with genomic data. Our analyses show: 1) 28 potential recombinant sequences were detected and they can be classified into seven potential recombinant types; 2) The spike (S) protein of MERS-CoV was under strong positive selection when MERS-CoV transmitted from their natural host to human; 3) Six out of nine positive selection sites detected in spike (S) protein are located in its receptor-binding domain which is in direct contact with host cells; 4) MERS-CoV frequently transmitted back and forth between human and camel after it had acquired the human-camel infection capability. Together, these results suggest that potential recombination events might have happened frequently during MERS-CoV's evolutionary history and the positive selection sites in MERS-CoV's S protein might enable it to infect human.

The Effects of Estradiol and Glucocorticoid on Human Osteosarcoma Cells: Similarities and Differences.

BACKGROUND: Human U2OS osteosarcoma cells were first derived from a tibial osteosarcoma of a teenage girl. They are a good cell model for osteosarcoma research in vitro. We compared the expression profiles of osteosarcoma cells after they were treated with estradiol and glucocorticoid, respectively. MATERIALS AND METHODS: We used published microarray data to compare the expression profiles of human osteosarcoma cells treated with estradiol and glucocorticoid. We investigated their differences and similarities with various bioinformatics tools. Oncogenes and tumor-suppressor genes differentially expressed after the hormone treatments were identified using online databases. The expression levels of cellular markers for proliferation and apoptosis were also compared before and after treatment with estradiol or glucocorticoid. RESULTS: Genes in human osteosarcoma cells differentially expressed after treatment with estradiol or glucocorticoid were detected. Our analysis showed that their similarity is more functionally prominent than their difference. Both estradiol and glucocorticoid can inhibit purine metabolic and biosynthetic pathway in human osteosarcoma cells. We also identified oncogenes and tumor-suppressor genes among the differentially expressed genes. The functional enrichment analyses of the identified cancer-associated genes suggests that estradiol has antagonistic effects on regulation of cell proliferation, while glucocorticoid can both arrest the cell cycle and prompt apoptosis. The effect of estradiol or glucocorticoid treatment on expression levels of cellular markers for proliferation and apoptosis support this argument. CONCLUSION: Our study suggests that glucocorticoid is more efficient in controlling osteosarcoma cell proliferation and apoptosis than estradiol.

The evolutionary panorama of organ-specifically expressed or repressed orthologous genes in nine vertebrate species.

RNA sequencing (RNA-Seq) technology provides the detailed transcriptomic information for a biological sample. Using the RNA-Seq data of six organs from nine vertebrate species, we identified a number of organ-specifically expressed or repressed orthologous genes whose expression patterns are mostly conserved across nine species. Our analyses show the following results: (i) About 80% of these genes have a chordate or more ancient origin and more than half of them are the legacy of one or multiple rounds of large-scale gene duplication events. (ii) Their evolutionary rates are shaped by the organ in which they are expressed or repressed, e.g. the genes specially expressed in testis and liver generally evolve more than twice as fast as the ones specially expressed in brain and cerebellum. The organ-specific transcription factors were discriminated from these genes. The ChIP-seq data from the ENCODE project also revealed the transcription-related factors that might be involved in regulating human organ-specifically expressed or repressed genes. Some of them are shared by all six human organs. The comparison of ENCODE data with mouse/chicken ChIP-seq data proposes that organ-specifically expressed or repressed orthologous genes are regulated in various combinatorial fashions in different species, although their expression features are conserved among these species. We found that the duplication events in some gene families might help explain the quick organ/tissue divergence in vertebrate lineage. The phylogenetic analysis of testis-specifically expressed genes suggests that some of them are prone to develop new functions for other organs/tissues.

Function and evolution of two forms of SecDF homologs in Streptomyces coelicolor.

The general secretion (Sec) pathway plays a prominent role in bacterial protein export, and the accessory component SecDF has been shown to improve transportation efficiency. Inspection of Streptomyces coelicolor genome reveals the unexpected presence of two different forms of secDF homologous genes: one in fused form (secDF) and the other in separated form (secD and secF). However, the functional role of two SecDF homologs in S. coelicolor has not yet been determined. Transcriptional analysis of secDF homologs reveals that these genes are constitutively expressed. However, the transcript levels of secD and secF are much higher than that of secDF in S. coelicolor. Deletion of secDF or/and secD/secF in S. coelicolor did result in reduced secretion efficiency of Xylanase A and Amylase C, suggesting that they may have redundant functions for Sec-dependent translocation pathway. Moreover, our results also indicate that SecD/SecF plays a more prominent role than SecDF in protein translocation. Evolutionary analysis suggests that the fused and separated SecDF homologs in Streptomyces may have disparate evolutionary ancestries. SecD/SecF may be originated from vertical transmission of existing components from ancestor of Streptomyces species. However, SecDF may be derived from bacterial ancestors through horizontal gene transfer. Alternately, it is also plausible that SecDF may have arisen through additional gene duplication and fusion events. The acquisition of a second copy may confer a selective benefit to Streptomyces by enhancing protein transport capacity. Taken together, our results provide new insights into the potential biological function and evolutionary aspects of the prokaryotic SecDF complex.

Phylogenomic distance method for analyzing transcriptome evolution based on RNA-seq data.

Thanks to the microarray technology, our understanding of transcriptome evolution at the genome level has been considerably advanced in the past decade. Yet, further investigation was challenged by several technical limitations of this technology. Recent innovation of next-generation sequencing, particularly the invention of RNA-seq technology, has shed insightful lights on resolving this problem. Though a number of statistical and computational methods have been developed to analyze RNA-seq data, the analytical framework specifically designed for evolutionary genomics remains an open question. In this article we develop a new method for estimating the genome expression distance from the RNA-seq data, which has explicit interpretations under the model of gene expression evolution. Moreover, this distance measure takes the data overdispersion, gene length variation, and sequencing depth variation into account so that it can be applied to multiple genomes from different species. Using mammalian RNA-seq data as example, we demonstrated that this expression distance is useful in phylogenomic analysis.