Orphan Genes Research Articles

Genetic innovations: Transposable element recruitment and de novo formation lead to the birth of orphan genes in the rice genome

AbstractOrphan genes are genetic innovations that lack homologs in other lineages. Orphan genes can rapidly originate and become substantially functional, yet the mechanisms underlying their origins are still largely unknown in plants. Here, we investigated the origin of orphan genes in the Oryza sativa ssp. japonica “Nipponbare” genome using genome‐wide comparisons with 10 closely related Oryza species. We identified a total of 37 orphan genes in the Nipponbare genome that show short sequence lengths, elevated GC content, and absence of introns. Interestingly, half of the identified orphan genes originated by way of a distinctive mechanism that involved the generation of new coding sequences through independent and rapid divergence within the inserted transposable element. Our results provide valuable insight into genetic innovations in the model rice genome that formed on a very short timescale.

Radical Technology Meets Radical Application: An Interview with George Church.

Radical Technology Meets Radical Application: An Interview with George Church.

Fagin: synteny-based phylostratigraphy and finer classification of young genes

BackgroundWith every new genome that is sequenced, thousands of species-specific genes (orphans) are found, some originating from ultra-rapid mutations of existing genes, many others originating de novo from non-genic regions of the genome. If some of these genes survive across speciations, then extant organisms will contain a patchwork of genes whose ancestors first appeared at different times. Standard phylostratigraphy, the technique of partitioning genes by their age, is based solely on protein similarity algorithms. However, this approach relies on negative evidence ─ a failure to detect a homolog of a query gene. An alternative approach is to limit the search for homologs to syntenic regions. Then, genes can be positively identified as de novo orphans by tracing them to non-coding sequences in related species.ResultsWe have developed a synteny-based pipeline in the R framework. Fagin determines the genomic context of each query gene in a focal species compared to homologous sequence in target species. We tested the fagin pipeline on two focal species, Arabidopsis thaliana (plus four target species in Brassicaseae) and Saccharomyces cerevisiae (plus six target species in Saccharomyces). Using microsynteny maps, fagin classified the homology relationship of each query gene against each target genome into three main classes, and further subclasses: AAic (has a coding syntenic homolog), NTic (has a non-coding syntenic homolog), and Unknown (has no detected syntenic homolog). fagin inferred over half the “Unknown” A. thaliana query genes, and about 20% for S. cerevisiae, as lacking a syntenic homolog because of local indels or scrambled synteny.Conclusionsfagin augments standard phylostratigraphy, and extends synteny-based phylostratigraphy with an automated, customizable, and detailed contextual analysis. By comparing synteny-based phylostrata to standard phylostrata, fagin systematically identifies those orphans and lineage-specific genes that are well-supported to have originated de novo. Analyzing within-species genomes should distinguish orphan genes that may have originated through rapid divergence from de novo orphans. Fagin also delineates whether a gene has no syntenic homolog because of technical or biological reasons. These analyses indicate that some orphans may be associated with regions of high genomic perturbation.

Differences in gene expression in whitefly associated with CYSDV-infected and virus-free melon, and comparison with expression in whiteflies fed on ToCV- and TYLCV-infected tomato

BackgroundCucurbit yellow stunting disorder virus (CYSDV; genus Crinivirus, Closteroviridae) is transmitted in a semipersistent manner by the whitefly, Bemisia tabaci, and is efficiently transmitted by the widely prevalent B. tabaci cryptic species, MEAM1. In this study, we compared transcriptome profiles of B. tabaci MEAM1, after 24 h, 72 h and 7 days of acquisition feeding on melon plants infected with CYSDV (CYSDV-whiteflies) with those fed on virus-free melon, using RNA-Seq technology. We also compared transcriptome profiles with whiteflies fed on tomato plants separately infected with Tomato chlorosis virus (ToCV), a crinivirus closely related to CYSDV, and Tomato yellow leaf curl virus (TYLCV), a member of the genus Begomovirus, which has a distinctly different mode of transmission and their respective virus-free controls, to find common gene expression changes among viruliferous whiteflies feeding on different host plants infected with distinct (TYLCV) and related (CYSDV and ToCV) viruses.ResultsA total of 275 differentially expressed genes (DEGs) were identified in CYSDV-whiteflies, with 3 DEGs at 24 h, 221 DEGs at 72 h, and 51 DEGs at 7 days of virus acquisition. Changes in genes encoding orphan genes (54 genes), phosphatidylethanolamine-binding proteins (PEBP) (20 genes), and AAA-ATPase domain containing proteins (10 genes) were associated with the 72 h time point. Several more orphan genes (20 genes) were differentially expressed at 7 days. A total of 59 common DEGs were found between CYSDV-whiteflies and ToCV-whiteflies, which included 20 orphan genes and 6 lysosomal genes. A comparison of DEGs across the three different virus-host systems revealed 14 common DEGs, among which, eight showed similar and significant up-regulation in CYSDV-whiteflies at 72 h and TYLCV-whiteflies at 24 h, while down-regulation of the same genes was observed in ToCV-whiteflies at 72 h.ConclusionsDynamic gene expression changes occurred in CYSDV-whiteflies after 72 h feeding, with decreased gene expression changes associated with 7 days of CYSDV acquisition. Similarities in gene expression changes among CYSDV-whiteflies, ToCV-whiteflies and TYLCV-whiteflies suggest the possible involvement of common genes or pathways for virus acquisition and transmission by whiteflies, even for viruses with distinctly different modes of transmission.

The wheat SnRK1α family and its contribution to Fusarium toxin tolerance

Deoxynivalenol (DON) is a mycotoxin produced by phytopathogenic Fusarium fungi in cereal grain and plays a role as a disease virulence factor. TaFROG (Triticum aestivum Fusarium Resistance Orphan Gene) enhances wheat resistance to DON and it interacts with a sucrose non-fermenting-1 (SNF1)-related protein kinase 1 catalytic subunit α (SnRK1α). This protein kinase family is central integrator of stress and energy signalling, regulating plant metabolism and growth. Little is known regarding the role of SnRK1α in the biotic stress response, especially in wheat. In this study, 15 wheat (Triticum aestivum) SnRK1α genes (TaSnRK1αs) belonging to four homoeologous groups were identified in the wheat genome. TaSnRK1αs are expressed ubiquitously in all organs and developmental stages apart from two members predominantly detected in grain. While DON treatment had either no effect or downregulated the transcription of TaSnRK1αs, it increased both the kinase activity associated with SnRK1α and the level of active (phosphorylated) SnRK1α. Down-regulation of two TaSnRK1αs homoeolog groups using virus induced gene silencing (VIGS) increased the DON-induced damage of wheat spikelets. Thus, we demonstrate that TaSnRK1αs contribute positively to wheat tolerance of DON and conclude that this gene family may provide useful tools for the improvement of crop biotic stress resistance.

Genome-wide identification and evolutionary analysis of TGA transcription factors in soybean

The gain of function in genes and gene families is a continuous process and is a key factor in understanding gene and genome evolution in plants. TGACG-Binding (TGA) transcription factors (TFs) have long been known for their essential roles in plant defence in Arabidopsis, but their roles in legume symbiosis are yet to be explored. Here, we identified a total of 25 TGA (named GmTGA1-GmTGA25) genes in soybean. Through phylogenetic analysis, we discovered a clade of GmTGA proteins that appear to be legume-specific. Among them, two GmTGAs were unique by possessing the autophagy sequence in their proteins, while the third one was an orphan gene in soybean. GmTGAs were structurally different from AtTGAs, and their expression patterns also differed with the dominant expression of AtTGAs and GmTGAs in aerial and underground parts, respectively. Moreover, twenty-five GmTGAs showed a strong correlation among the gene expression in roots, nodules, and root hairs. The qRT-PCR analysis results revealed that among 15 tested GmTGAs, six were induced and four were suppressed by rhizobia inoculation, while 11 of these GmTGAs were induced by high nitrate. Our findings suggested the important roles of GmTGAs in symbiotic nodulation and in response to nitrogen availability in soybean.

Engineered Production of Hapalindole Alkaloids in the Cyanobacterium Synechococcus sp. UTEX 2973.

Cyanobacteria produce numerous valuable bioactive secondary metabolites (natural products) including alkaloids, isoprenoids, nonribosomal peptides, and polyketides. However, the genomic organization of the biosynthetic gene clusters, complex gene expression patterns, and low compound yields synthesized by the native producers currently limits access to the vast majority of these valuable molecules for detailed studies. Molecular cloning and expression of such clusters in heterotrophic hosts is often precarious owing to genetic and biochemical incompatibilities. Production of such biomolecules in photoautotrophic hosts analogous to the native producers is an attractive alternative that has been under-explored. Here, we describe engineering of the fast-growing cyanobacterium Synechococcus elongatus UTEX 2973 to produce key compounds of the hapalindole family of indole-isonitrile alkaloids. Engineering of the 42-kbp "fam" hapalindole pathway from the cyanobacterium Fischerella ambigua UTEX 1903 into S2973 was accomplished by rationally reconstructing six to seven core biosynthetic genes into synthetic operons. The resulting Synechococcus strains afforded controllable production of indole-isonitrile biosynthetic intermediates and hapalindoles H and 12-epi-hapalindole U at a titer of 0.75-3 mg/L. Exchanging genes encoding fam cyclase enzymes in the synthetic operons was employed to control the stereochemistry of the resulting product. Establishing a robust expression system provides a facile route to scalable levels of similar natural and new forms of bioactive hapalindole derivatives and its structural relatives (e.g., fischerindoles, welwitindolinones). Moreover, this versatile expression system represents a promising tool for exploring other functional characteristics of orphan gene products that mediate the remarkable biosynthesis of this important family of natural products.

Toxins from scratch? Diverse, multimodal gene origins in the predatory robber fly Dasypogon diadema indicate a dynamic venom evolution in dipteran insects.

BackgroundVenoms and the toxins they contain represent molecular adaptations that have evolved on numerous occasions throughout the animal kingdom. However, the processes that shape venom protein evolution are poorly understood because of the scarcity of whole-genome data available for comparative analyses of venomous species.ResultsWe performed a broad comparative toxicogenomic analysis to gain insight into the genomic mechanisms of venom evolution in robber flies (Asilidae). We first sequenced a high-quality draft genome of the hymenopteran hunting robber fly Dasypogon diadema, analysed its venom by a combined proteotranscriptomic approach, and compared our results with recently described robber fly venoms to assess the general composition and major components of asilid venom. We then applied a comparative genomics approach, based on 1 additional asilid genome, 10 high-quality dipteran genomes, and 2 lepidopteran outgroup genomes, to reveal the evolutionary mechanisms and origins of identified venom proteins in robber flies.ConclusionsWhile homologues were identified for 15 of 30 predominant venom protein in the non-asilid genomes, the remaining 15 highly expressed venom proteins appear to be unique to robber flies. Our results reveal that the venom of D. diadema likely evolves in a multimodal fashion comprising (i) neofunctionalization after gene duplication, (ii) expression-dependent co-option of proteins, and (iii) asilid lineage-specific orphan genes with enigmatic origin. The role of such orphan genes is currently being disputed in evolutionary genomics but has not been discussed in the context of toxin evolution. Our results display an unexpected dynamic venom evolution in asilid insects, which contrasts the findings of the only other insect toxicogenomic evolutionary analysis, in parasitoid wasps (Hymenoptera), where toxin evolution is dominated by single gene co-option. These findings underpin the significance of further genomic studies to cover more neglected lineages of venomous taxa and to understand the importance of orphan genes as possible drivers for venom evolution.

De Novo, Divergence, and Mixed Origin Contribute to the Emergence of Orphan Genes in Pristionchus Nematodes.

Homology is a fundamental concept in comparative biology. It is extensively used at the sequence level to make phylogenetic hypotheses and functional inferences. Nonetheless, the majority of eukaryotic genomes contain large numbers of orphan genes lacking homologs in other taxa. Generally, the fraction of orphan genes is higher in genomically undersampled clades, and in the absence of closely related genomes any hypothesis about their origin and evolution remains untestable. Previously, we sequenced ten genomes with an underlying ladder-like phylogeny to establish a phylogenomic framework for studying genome evolution in diplogastrid nematodes. Here, we use this deeply sampled data set to understand the processes that generate orphan genes in our focal species Pristionchus pacificus. Based on phylostratigraphic analysis and additional bioinformatic filters, we obtained 29 high-confidence candidate genes for which mechanisms of orphan origin were proposed based on manual inspection. This revealed diverse mechanisms including annotation artifacts, chimeric origin, alternative reading frame usage, and gene splitting with subsequent gain of de novo exons. In addition, we present two cases of complete de novo origination from non-coding regions, which represents one of the first reports of de novo genes in nematodes. Thus, we conclude that de novo emergence, divergence, and mixed mechanisms contribute to novel gene formation in Pristionchus nematodes.

A stemness screen reveals C3orf54/INKA1 as a promoter of human leukemia stem cell latency

AbstractThere is a growing body of evidence that the molecular properties of leukemia stem cells (LSCs) are associated with clinical outcomes in acute myeloid leukemia (AML), and LSCs have been linked to therapy failure and relapse. Thus, a better understanding of the molecular mechanisms that contribute to the persistence and regenerative potential of LSCs is expected to result in the development of more effective therapies. We therefore interrogated functionally validated data sets of LSC-specific genes together with their known protein interactors and selected 64 candidates for a competitive in vivo gain-of-function screen to identify genes that enhanced stemness in human cord blood hematopoietic stem and progenitor cells. A consistent effect observed for the top hits was the ability to restrain early repopulation kinetics while preserving regenerative potential. Overexpression (OE) of the most promising candidate, the orphan gene C3orf54/INKA1, in a patient-derived AML model (8227) promoted the retention of LSCs in a primitive state manifested by relative expansion of CD34+ cells, accumulation of cells in G0, and reduced output of differentiated progeny. Despite delayed early repopulation, at later times, INKA1-OE resulted in the expansion of self-renewing LSCs. In contrast, INKA1 silencing in primary AML reduced regenerative potential. Mechanistically, our multidimensional confocal analysis found that INKA1 regulates G0 exit by interfering with nuclear localization of its target PAK4, with concomitant reduction of global H4K16ac levels. These data identify INKA1 as a novel regulator of LSC latency and reveal a link between the regulation of stem cell kinetics and pool size during regeneration.

A wheat NAC interacts with an orphan protein and enhances resistance to Fusarium head blight disease.

SummaryTaxonomically‐restricted orphan genes play an important role in environmental adaptation, as recently demonstrated by the fact that the Pooideae‐specific orphan TaFROG (Triticum aestivum Fusarium Resistance Orphan Gene) enhanced wheat resistance to the economically devastating Fusarium head blight (FHB) disease. Like most orphan genes, little is known about the cellular function of the encoded protein TaFROG, other than it interacts with the central stress regulator TaSnRK1α. Here, we functionally characterized a wheat (T. aestivum) NAC–like transcription factor TaNACL‐D1 that interacts with TaFROG and investigated its’ role in FHB using studies to assess motif analyses, yeast transactivation, protein‐protein interaction, gene expression and the disease response of wheat lines overexpressing TaNACL‐D1. TaNACL‐D1 is a Poaceae‐divergent NAC transcription factor that encodes a Triticeae‐specific protein C‐terminal region with transcriptional activity and a nuclear localisation signal. The TaNACL‐D1/TaFROG interaction was detected in yeast and confirmed in planta, within the nucleus. Analysis of multi‐protein interactions indicated that TaFROG could form simultaneously distinct protein complexes with TaNACL‐D1 and TaSnRK1α in planta. TaNACL‐D1 and TaFROG are co‐expressed as an early response to both the causal fungal agent of FHB, Fusarium graminearum and its virulence factor deoxynivalenol (DON). Wheat lines overexpressing TaNACL‐D1 were more resistant to FHB disease than wild type plants. Thus, we conclude that the orphan protein TaFROG interacts with TaNACL‐D1, a NAC transcription factor that forms part of the disease response evolved within the Triticeae.

Orphan genes are involved in drought adaptations and ecoclimatic-oriented selections in domesticated cowpea.

Orphan genes (OGs) are genes that are restricted to a single species or a particular taxonomic group. To date, little is known about the functions of OGs in domesticated crops. Here, we report our findings on the relationships between OGs and environmental adaptation in cowpea (Vigna unguiculata). We identified 578 expressed OGs, of which 73.2% were predicted to be non-coding. Transcriptomic analyses revealed a high rate of OGs that were drought inducible in roots when compared with conserved genes. Co-expression analysis further revealed the possible involvement of OGs in stress response pathways. Overexpression of UP12_8740, a drought-inducible OG, conferred enhanced tolerance to osmotic stresses and soil drought. By combining Capture-Seq and fluorescence-based Kompetitive allele-specific PCR (KASP), we efficiently genotyped single nucleotide polymorphisms (SNPs) on OGs across a 223 accession cowpea germplasm collection. Population genomic parameters, including polymorphism information content (PIC), expected heterozygosity (He), nucleotide diversity (π), and Tajima's D statistics, that were calculated based on these SNPs, showed distinct signatures between the grain- and vegetable-type subpopulations of cowpea. This study reinforces the idea that OGs are a valuable resource for identifying new genes related to species-specific environmental adaptations and fosters new insights that artificial selection on OGs might have contributed to balancing the adaptive and agronomic traits in domesticated crops in various ecoclimatic conditions.

Transcriptomic atlas of mushroom development reveals conserved genes behind complex multicellularity in fungi

The evolution of complex multicellularity has been one of the major transitions in the history of life. In contrast to simple multicellular aggregates of cells, it has evolved only in a handful of lineages, including animals, embryophytes, red and brown algae, and fungi. Despite being a key step toward the evolution of complex organisms, the evolutionary origins and the genetic underpinnings of complex multicellularity are incompletely known. The development of fungal fruiting bodies from a hyphal thallus represents a transition from simple to complex multicellularity that is inducible under laboratory conditions. We constructed a reference atlas of mushroom formation based on developmental transcriptome data of six species and comparisons of >200 whole genomes, to elucidate the core genetic program of complex multicellularity and fruiting body development in mushroom-forming fungi (Agaricomycetes). Nearly 300 conserved gene families and >70 functional groups contained developmentally regulated genes from five to six species, covering functions related to fungal cell wall remodeling, targeted protein degradation, signal transduction, adhesion, and small secreted proteins (including effector-like orphan genes). Several of these families, including F-box proteins, expansin-like proteins, protein kinases, and transcription factors, showed expansions in Agaricomycetes, many of which convergently expanded in multicellular plants and/or animals too, reflecting convergent solutions to genetic hurdles imposed by complex multicellularity among independently evolved lineages. This study provides an entry point to studying mushroom development and complex multicellularity in one of the largest clades of complex eukaryotic organisms.

Origination and evolution of orphan genes and de novo genes in the genome of Caenorhabditis elegans.

Orphan genes that lack detectable homologues in other lineages could contribute to a variety of biological functions. However, their origination and function mechanisms remain largely unknown. Herein, through a comprehensive and systematic computational pipeline, we identified 893 orphan genes in the lineage of C. elegans, of which only a low fraction (0.9%) were derived from transposon elements. Six new protein-coding genes that de novo originated from non-coding DNA sequences in the genome of C. elegans were also identified. The authenticity and functionality of these orphan genes and de novo genes are supported by three lines of evidences, consisting of transcriptional data, and in silico proteomic data, and the fixation status data in wild populations. Orphan genes and de novo genes exhibited simple gene structures, such as, short in protein length, of fewer exons, and are frequently X-linked. RNA-seq data analysis showed these orphan genes are enriched with expression in embryo development and gonad, and their potential function in early development was further supported by gene ontology enrichment analysis results. Meanwhile, de novo genes were found to be with significant expression in gonad, and functional enrichment analysis of the co-expression genes of these de novo genes suggested they may be functionally involved in signaling transduction pathway and metabolism process. Our results presented the first systematic evidence on the evolution of orphan genes and de novo origin of genes in nematodes and their impacts on the functional and phenotypic evolution, and thus could shed new light on our appreciation of the importance of these new genes.

Topological evolution of coexpression networks by new gene integration maintains the hierarchical and modular structures in human ancestors.

We analyze the global structure and evolution of human gene coexpression networks driven by new gene integration. When the Pearson correlation coefficient is greater than or equal to 0.5, we find that the coexpression network consists of 334 small components and one "giant" connected subnet comprising of 6317 interacting genes. This network shows the properties of power-law degree distribution and small-world. The average clustering coefficient of younger genes is larger than that of the elderly genes (0.6685 vs. 0.5762). Particularly, we find that the younger genes with a larger degree also show a property of hierarchical architecture. The younger genes play an important role in the overall pivotability of the network and this network contains few redundant duplicate genes. Moreover, we find that gene duplication and orphan genes are two dominant evolutionary forces in shaping this network. Both the duplicate genes and orphan genes develop new links through a "rich-gets-richer" mechanism. With the gradual integration of new genes into the ancestral network, most of the topological structure features of the network would gradually increase. However, the exponent of degree distribution and modularity coefficient of the whole network do not change significantly, which implies that the evolution of coexpression networks maintains the hierarchical and modular structures in human ancestors.

Comparative analysis of gene family size provides insight into the adaptive evolution of vertebrates

Copy numbers of homologous gene families vary greatly among different species, which is caused by the differences in the rates of gene gain and loss. It is well known that gene copy number variation can be responsible for the phenotypic novelties of particular species. In this study, 64 species that represent the main vertebrate groups spanning evolutionary period of about 600 million years were selected and the homology of gene families across these species were established, thereby revealing the evolutionary patterns of gene family size in vertebrates. The results show that among the 6857 gene families inferred to be present in the most recent common ancestor of the vertebrates, 6712 had changed their sizes in at least one lineage, and these gene families had contracted in most cases. Gene families in Choloepus hoffmanni and Danio rerio had undergone the greatest contraction and expansion, respectively. Based on the highly dynamic evolution of vertebrate gene family size, we sought to identify any genomic signals that might be related to the evolution of specific vertebrate populations from the perspective of the distinct gene family size changes. We observed a high proportion of gene family amplification occurred, probably due to genome-wide duplication in the recent common ancestral genome of teleosts, which was followed by contraction in the decedents due to the extensive gene fractionation. Furthermore, we found evidence that orphan genes in the bony fish might contribute to the adaptive evolution of fish in aquatic environment. For example, some orphan genes were involved in fin development, tail development and kidney physiology. Overall, our work provides novel insights into the evolution of vertebrate gene family size and provides several lines of evidence for understanding the relationship between the genome evolution and phenotypic diversity in vertebrates.

Orphan Genes Shared by Pathogenic Genomes Are More Associated with Bacterial Pathogenicity.

Orphan genes (also known as ORFans [i.e., orphan open reading frames]) are new genes that enable an organism to adapt to its specific living environment. Our focus in this study is to compare ORFans between pathogens (P) and nonpathogens (NP) of the same genus. Using the pangenome idea, we have identified 130,169 ORFans in nine bacterial genera (505 genomes) and classified these ORFans into four groups: (i) SS-ORFans (P), which are only found in a single pathogenic genome; (ii) SS-ORFans (NP), which are only found in a single nonpathogenic genome; (iii) PS-ORFans (P), which are found in multiple pathogenic genomes; and (iv) NS-ORFans (NP), which are found in multiple nonpathogenic genomes. Within the same genus, pathogens do not always have more genes, more ORFans, or more pathogenicity-related genes (PRGs)-including prophages, pathogenicity islands (PAIs), virulence factors (VFs), and horizontal gene transfers (HGTs)-than nonpathogens. Interestingly, in pathogens of the nine genera, the percentages of PS-ORFans are consistently higher than those of SS-ORFans, which is not true in nonpathogens. Similarly, in pathogens of the nine genera, the percentages of PS-ORFans matching the four types of PRGs are also always higher than those of SS-ORFans, but this is not true in nonpathogens. All of these findings suggest the greater importance of PS-ORFans for bacterial pathogenicity. IMPORTANCE Recent pangenome analyses of numerous bacterial species have suggested that each genome of a single species may have a significant fraction of its gene content unique or shared by a very few genomes (i.e., ORFans). We selected nine bacterial genera, each containing at least five pathogenic and five nonpathogenic genomes, to compare their ORFans in relation to pathogenicity-related genes. Pathogens in these genera are known to cause a number of common and devastating human diseases such as pneumonia, diphtheria, melioidosis, and tuberculosis. Thus, they are worthy of in-depth systems microbiology investigations, including the comparative study of ORFans between pathogens and nonpathogens. We provide direct evidence to suggest that ORFans shared by more pathogens are more associated with pathogenicity-related genes and thus are more important targets for development of new diagnostic markers or therapeutic drugs for bacterial infectious diseases.

Symbiotic organs shaped by distinct modes of genome evolution in cephalopods

Microbes have been critical drivers of evolutionary innovation in animals. To understand the processes that influence the origin of specialized symbiotic organs, we report the sequencing and analysis of the genome of Euprymna scolopes, a model cephalopod with richly characterized host-microbe interactions. We identified large-scale genomic reorganization shared between E. scolopes and Octopus bimaculoides and posit that this reorganization has contributed to the evolution of cephalopod complexity. To reveal genomic signatures of host-symbiont interactions, we focused on two specialized organs of E. scolopes: the light organ, which harbors a monoculture of Vibrio fischeri, and the accessory nidamental gland (ANG), a reproductive organ containing a bacterial consortium. Our findings suggest that the two symbiotic organs within E. scolopes originated by different evolutionary mechanisms. Transcripts expressed in these microbe-associated tissues displayed their own unique signatures in both coding sequences and the surrounding regulatory regions. Compared with other tissues, the light organ showed an abundance of genes associated with immunity and mediating light, whereas the ANG was enriched in orphan genes known only from E. scolopes Together, these analyses provide evidence for different patterns of genomic evolution of symbiotic organs within a single host.

Sequence Analysis of Alpha Crystalline Domain of Small Heat Shock Proteins (sHSPs) Family in Different Plant Species

Alpha crystalline domain (ACD/HSP20) proteins belong to a diverse family of molecular chaperones which play a central role in different plant species under stress condition. We have studied about 100 sequences having alpha crystalline domains from Arabidopsis thaliana, Zea mays, Oryza sativa, Sorghum bicolor, and Cicer arietinum. High residual conservancy in ACD domain and variation in the N-terminal region were observed in these plant species. We also observed the orphan genes like transcripts. Phylogenetic studies classified them into three well-categorized clades with eleven subgroups/families according to their amino acid sequences of a conserved domains. Duplication events and basic to the intricate pattern of conservancy within the conserved domain of subgroups existed. Exon/intron boundaries and cellular localization pattern further support the phylogenetic-based categorization of these genes. About 50% of them are cytoplasm-specific while the others are chloroplast-, mitochondria-, and extracellular matrix-specific tags. We also detected the conserved patterns of functional motifs within the coding sequences of heat shock protein genes of above-mentioned plant species. Overall expression pattern was found to be much higher in Arabidopsis, rice, and maize under heat, salt, and drought stress, respectively, than cold stress while homogenous expression was observed at tissue-specific level. These results would be useful to understand the further functional characterization of small heat shock proteins.

The genome of the jellyfish Aurelia and the evolution of animal complexity

We present the genome of the moon jellyfish Aurelia, a genome from a cnidarian with a medusa life stage. Our analyses suggest that gene gain and loss in Aurelia is comparable to what has been found in its morphologically simpler relatives—the anthozoan corals and sea anemones. RNA sequencing analysis does not support the hypothesis that taxonomically restricted (orphan) genes play an oversized role in the development of the medusa stage. Instead, genes broadly conserved across animals and eukaryotes play comparable roles throughout the life cycle. All life stages of Aurelia are significantly enriched in the expression of genes that are hypothesized to interact in protein networks found in bilaterian animals. Collectively, our results suggest that increased life cycle complexity in Aurelia does not correlate with an increased number of genes. This leads to two possible evolutionary scenarios: either medusozoans evolved their complex medusa life stage (with concomitant shifts into new ecological niches) primarily by re-working genetic pathways already present in the last common ancestor of cnidarians, or the earliest cnidarians had a medusa life stage, which was subsequently lost in the anthozoans. While we favour the earlier hypothesis, the latter is consistent with growing evidence that many of the earliest animals were more physically complex than previously hypothesized.